チャネル乱流についての質問
- チャネル乱流の条件設定について質問です。
- 現在、DNSやLESでチャネル乱流を行っているが、流れ場の乱れと平均速度分布が一致しない。
- 条件（無次元化、平均圧力勾配など）が誤っている可能性が高いため、質問したい。
- ベストアンサー
チャネル乱流について
現在、DNSやLESでチャネル乱流をやっているのですが、流れ場は乱れるものの、平均速度分布などがまったく合いません。 条件（無次元化、平均圧力勾配など）が間違っている可能性が高いため、それについて質問いたします。 記号（変数）は Re_tau:壁指標のレイノルズ数 U_tau:壁面摩擦速度 delta:チャンネル半幅 v:動粘性係数 dp/dx:平均圧力勾配 rho:密度 です。 非圧縮のNavier-Stokes方程式を考え、主流方向に外力（平均圧力勾配）を与えています。 Re_tau=180として、U_tau=1とするため、v=1/360としています。（delta=0.5） また、平均圧力勾配は dp/dx = -rho * U_tau^2/delta = -rho * v^2 * Re_tau^2/delta^3 から決めています。 解析結果の流速は乱流の分布と同じようなものになっていますが、値は250程度になります。 U_tau=1としているため、結果の流速u=u+（壁指標の流速）と考えております。 上の条件の設定は間違っているのでしょうか？ どなたかご存知の方がいましたら、教えてください。 どうかよろしくお願いいたします。
- fluct
- お礼率100% (1/1)
- 物理学
- 回答数1
- ありがとう数2
- みんなの回答 （1）
- 専門家の回答
質問者が選んだベストアンサー
私もチャネルの計算をやっています． この計算のポイントは圧力勾配です． レイノルズ応力の分布を見ると圧力勾配が正しいか どうかのチェックが出来ます． で，壁面摩擦速度，チャンネル幅で‘無次元化’する と，密度や動粘性係数は陽にはでてきません． 実際には圧力勾配はチャンネル幅を2δとして代表長さに チャンネル半値とした場合，dp/dx=1と決定されます．
関連するQ&A
- 流体力学についての質問なのですが・・・
流体力学についての質問なのですが・・・ 水平に設置された円管内を非圧縮性流体が一定の流量で流れていて、円管の途中で直径が半分になっていたとき 管内平均流速・最大流速・壁面せん断応力・レイノルズ数は何倍になっているか、という問題なのですが 円管直径をdとすると、ハーゲンポアズイユの式から 平均流速 Um(R=d)=-d^2/32*dp/dx um(R=d/2)=-d^2/32*1/4*dp/dx=1/4*Um と解答を導いて1/4倍という答えを求めたのですが、よく考えてみると円管を細くすると流速は上がると思うのですが計算で導けなくて困っています。 この解答では間違っているのでしょうか？ せん断応力Τ(R=d)=d/4*ΔP/L τ(R=d/2)=d/8*ΔP/L レイノルズ数 Re(R=d)=Um・L/v Re(R=2/d)=1/4*Um・L/v （自分の答え） 平均流速 1/4倍 最大流速 1/4倍 せん断応力 1/2倍 レイノルズ数 1/4倍
- 締切済み
- 物理学
- 乱流における管断面平均流速について
乱流の場合、指数法則によると、 U/Umax = (1-r/R)^（1/n） U：速度分布、Umax：最大速度、r：仮想管半径、R：円管半径 なる実験式が得られますが、 n=7 のときの、管断面平均速度を求めたいため、 V = 1/(πR^2)∫2πrUdr （積分範囲は、0からR）に、上記の実験式を用いて計算をしたいのですが、 7乗根をうまく積分する方法がわかりません。 テキストでは中間の計算式が省いてあるため、困っております。 ご教授をお願いいたします。
- ベストアンサー
- 物理学
- 流体
次の問題を解いてみたのですが途中からわかりません。また自分で解いた答えも間違えているかもしれないのでわかる方がいましたら解法を教えてください。お願いします。 [問題] 下図は静止したシリンダの中で、シリンダと同軸の半径Rのドラムが一定角速度Ωで回転している様子を表したものである。シリンダには溝1と溝２がつけられている。流れは図に示された平面内の二次元流れとなっている。経路CDの隙間間隔は経路AEBのみを通って溝２より流出すると考えてよい。点Bでの圧力は点Aでの圧力よりも△Pだけ高くなっている。点Aを原点としてシリンダ表面に沿ってx座標、これと垂直にｙ座標をとる。経路AEBのすきま間隔はｈでRよりも十分小さく、すきま内の流れは十分発達した層流として、以下の問いに答えよ。 (１)すきま内の流れは以下の式で近似できる。 d^2u/dy^2=(1/μ）(dP/dx) ただしuは流速、Pは圧力、μは粘性係数である。この式を適切な条件下で解き、流速を求めよ。(できたら速度分布を図示してください） (２)紙面に垂直な単位幅あたりの体積流量Qを求めなさい。 (３)(２)の結果を用いて、経路AEBのAB間の距離lとして圧力差△Pを求めよ。 (４)ドラム表面に働くせん断応力τを求めなさい。 (５）ドラム表面から毎秒与えられる仕事のうち、Q△Pだけが流体に与えられるとして、流体輸送の効率をもとめなさい。ただし、ドラム表面からは経路AEBのすきまを通して仕事が与えられるものとする。 ---------------------------------------------------------------------------------- [自分の答え] (1) 二回ｙで積分して u=(1/2μ)(dP/dx)(y^2-hy) （ ただし境界条件y=0のときu=0,y=hのときu=0より） (２) Q=∫(0→h)udh=-(h^3/12)(dP/dx) (3) Q=-(h^3/12)(△P/l) △P＝－12μl/h^3 (4) ニュートンの粘性法則より τ＝μdu/dy=(1/2)(dp/dx)(2y-h) (5) ?
- ベストアンサー
- 物理学
- 流体力学の問題について
以下の流体力学の問題がわからないので、どなたか教えていただけませんか。途中式もよろしくお願いします。 2-1）粘性流体の方程式は次式で与えられる。流速ベクトルｖが3成分v=（u,v,w)を持つとして、運動方程式のｘ成分を書き下しなさい。 ∂ｖ/∂ｔ+（ｖ・∇）ｖ＝－１/ｐ+ｖ∇²ｖ 2-2）いま、ｘ軸方向の長さが１の2枚の板(紙面に直交するｚ軸方向には無限に長い）がｙ軸方向に距離Lを隔てて平行に置かれている(下の画像参照）。この2枚の板の間に粘性係数および密度が一様な非圧縮性流体が存在しているとする。X軸方向およびZ軸方向には流速は一様であると仮定する。いま、両端で定常な圧力差が加えられており、流体中ではｄｐ/ｄｘは一様であると仮定する。このような条件下では流速はX軸方向の流れのみとなりu=u(y)の形をとると考えてよい。 この場合の、運動方程式のｘ成分を書き下しなさい。さらに、u=u(y)を求め、流れが2次曲線となることを示しなさい。 以上のような問題です。誠に勝手ながら、どうかよろしくお願いいたします。
- ベストアンサー
- 物理学
- 壁にぶつかる噴流による圧力について。
流速u(m/s)、断面積S(m2)で一様な噴流が壁にぶつかって放射状に散らばっているとき、壁に一秒間にぶつかる噴流の質量はρ•u•S(kg)で、ぶつかった後の運動量は合計ゼロと考えることができるので、一秒間に壁に与える力積は壁にぶつかる運動量と等しく、Ft = ρ•u•S•u = ρ•(u^2)•S (N•s)となり、壁はF = ρ•(u^2)•S (N)で押されていることになります。従って壁の直近の圧力はρ•(u^2)•S/S = ρ•(u^2) (Pa)になると思います。 一方、この噴流は単位体積あたりρ•(u^2)/2の運動エネルギー(動圧)を持っているので、流線上でロス無く流速をゼロにした場合はベルヌーイの定理によればρ•(u^2)/2の圧力(静圧)が生じる筈で、先の運動量から考えた計算の半分の圧力となります。ロスがある場合は圧力はさらに小さくなって、先の運動量から考えた計算との開きはさらに大きくなります。 ベルヌーイの定理の場合に何故ρ•(u^2)の半分になるのかという点について考えてみました。流速がuから0に瞬時に変化することは不可能だから、流速がuから0に変化する区間で圧力は0からρ•(u^2)まで変化する。だからこの区間の平均圧力を取ればρ•(u^2)/2になる。区間を極限まで狭くしたとしても平均がρ•(u^2)/2になることに変わりはないから圧力はρ•(u^2)/2になるのだろうと思います。 後者のベルヌーイの定理の動圧から考えた圧力と、前者の運動量から考えた圧力は、どう考えれば矛盾なく理解できるでしょうか。
- ベストアンサー
- 物理学
- 円管内の層流について
円管内の層流についてなんですが、式を証明していくと、μ/r・d/dr(r・du/dr)=dp/dxという式が出てくるんですが、私が計算していきますと、μ/r・d/dr(u+r・du/dr)=dp/dxとなるんですが、u=0として良いのでしょうか？として良いのなら、なぜ、なんでしょうか？どなたかお願いします。条件は、速度は流れ方向に変化しないため、v=0,かつ、uはx方向に変化しないためu=u(r),そして、定常流れ、かつ、体積力がないと仮定しています。お願いします。
- ベストアンサー
- 物理学
- 超音波流量計の測定原理
伝播時間差法に基づく超音波流量計の測定原理を調べているのですが、1点分からないところがあり、ご教示いただければ幸いです。 超音波流量計で実際に測定されるのは超音波の伝播路に沿った線平均流速Vですが、流量の計算に用いるのは管の横断面での面平均流速Vaです。線平均流速と面平均流速は異なるため、これを補正する補正係数Kを用いて流量Qは次のように計算されます。（rは管の内半径です） Q＝π･r^2･Va＝π･r^2･K･V 流れが層流の場合は補正係数Kは理論的に一定ですが、流れが乱流の場合は補正係数Kはレイノルズ数Reの関数になります。 流れが乱流の場合に、流量計の校正により、レイノルズ数Reと補正係数Kの関係を求めることはできると思います。 分からないのはこここからです。 実際に超音波流量計で流量を測定するときに実測できるのは伝播時間であり、装着条件により決まる入射角，管内半径を用いて線平均速度が計算されます。しかし、この段階ではまだ流量は得られていないのでレイノルズ数Reは未知であり、したがって、補正係数Kも未知です。どのような計算で流量を計算するのか分かりません。 ご教示よろしくお願いします。
- 締切済み
- その他（品質管理）
- ハーゲン・ポアズイユ流れの式について
こんにちは。 さっそく質問させていただきます。 半径aのまっすぐな円管を考える。このときの層流の流速分布は u(r)=g*I/4ν(a^2-r^2)・・・・・(1) である。この式を流量Qの形にすると、 Q=（πgI /8ν）*a^4・・・・・(2) となる。(1)の式から(2)の式を導けという、問題です。 半径a 層流分布u(r) 重力加速度ｇ 動粘性係数ν 動水勾配I です。 いくら考えても、全く解くことができません。 どなたか、導くための式を教えてください。よろしくお願いします。
- ベストアンサー
- 物理学
- 翼を通り過ぎる流体の2次元定常流れの解析
下の図で示した2次元定常流れで翼を通り過ぎたあとの流体の密度ρ1，流速u1，圧力p1がわかりません． 領域Vにある流体は矢印の方向へ流れており十分な上流では密度ρ0，流速u0，圧力p0となっています．図にある翼は等間隔Hで並べられており，翼が受ける抗力はFx，揚力はFyと設定されています．流体は十分な下流では水平であるx軸から角度θ(<90°)下方へ向いた一様流となっています．この時の密度ρ1，流速u1，圧力p1を運動量保存則から求められないかと考えたのですが，x方向の流速u，y軸方向の流速vとした場合で x方向の運動量保存則より ∂(ρudxdy)/∂t=(ρu^2dy+ρuvdx)ー{(ρu^2dy+∂(ρu^2)dx/∂x)dy+(ρuvdx+∂(ρuv)dy/∂y)dx}ーFx ⇔∂(ρu^2)/∂x+∂(ρuv)/∂y=ーFx/dxdy y軸方向の運動量保存則より (中略)⇔∂(ρuv)/∂x+∂(ρv^2)/∂y=ーFy/dxdy からどのようにしてρ1やu1を求めればよいのかわかりません．p1については前の2つの値がわかればベルヌーイの定理から求められそうだということは予測できるのですが…． 問題の問では「領域Vに流入・流出する運動量のy成分とV内の流体に働くちからのy成分の釣り合いを考えることで下流の流速u1，密度ρ1を求めよ」とあります． 長い式が出てきましたが回答のほどよろしくお願い致します.
- 締切済み
- 物理学
お礼
ご回答いただきまして、ありがとうございます。 返答が遅れまして、本当に申し訳ありませんでした。 現在は、無次元化したN-S式においてRe_tauと圧力勾配のみを与えればよいと考えております。 チャンネル幅と圧力勾配の関係は、 文献などでは、仰るとおり、2δで1、δでは2としているようなのですが、 これはどのようにして決まるものなのでしょうか？ 非常に初歩的な質問で申し訳ないのですが、レイノルズ応力を求める際に、 <u'u'>(<>:空間平均,u:x方向の変動速度)などを計算する必要があると思うのですが、このときのu'は u'=u-<u> と単純に計算してもよろしいのでしょうか？ 身近に乱流の解析について知っている人がいないのと、私のCFDの経験が浅いため、ごく簡単なこともわからない状態です。 申し訳ありませんが、時間のおありになるときでかまいませんので、ご回答いただけますようよろしくお願いいたします。