輸送現象まとめファイル2017_01

Transcription

1 第 5 章流れの解析 5.1 平行平板内流れ 流体の連続の式と運動量方程式を連立して解き 流体内部の速 度分布を求める 第一の例は 平行平板内流れの定常状態速度分 布である 非圧縮性流体 ( 密度 ρ が一定 ) 定常状態 ( t の項 が0) 重力以外の外力が働かないとき デカルト座標系の -! 二次元速度分布 v, ( ) を決める方程式は次の様に書ける 連続の式 : u + v = 0 (5-1) Fig. 5-1 平行平板層流速度分布 軸方向運動量方程式 : u u + v u = 1 p ρ +ν u + u + g (5-) 軸方向運動量方程式 : u v + v v = 1 p ρ +ν v + v + g (5-3) 境界条件は 次の通りである =0 で u=u0 v=0 p=p0 (5-4a) =0 で u=v=0 (5-4b) 未知数は u, v, p 式の数は 3 つなので 境界条件を使って式を解くことができる 方向速度は 入口で一様流 入速度 u=u0 から次第に変化し Fig.5-1 の様にある区間 ( 速度助走区間 ) を経過すると u は 座標のみの関数と なる また 方向速度 v は v=0 となり 方向圧力損失 ( dp d = Δp Δ ) は一定になる 従って dp d 以外の を含む項と v を含む項を無視する h を中心から平行平板までの距離とし 方向に重力が働かないので 発達領 域で (5-) 式は次式のように簡略化される 0 = 1 dp ρ d +ν d u d (5-5) =0 で u=0 =h で du d = 0 (5-6) 発達速度領域で 平行平板間強制対流速度分布は次式となる u = u ma h h, u = h Δp ma µ Δ (5-7) 次の例は Fig.5- のように 重力のみが働き厚さ h の液膜が傾斜 した坂の上を流れているとき g = ρg sinθ となる また自由表面流 Fig.5- 落下液膜の流速分布 であり外からの圧力は均一一定である 発達速度状態の速度分布は 上と同様の手段で求められ 最大速度は次式 で得られる u ma = ρh g µ sinθ ( 液膜落下のとき ) (5-8) 6

2 5. 円管内速度分布 円管内速度分布は広く流体 熱 物質移動によく現れる (,,z) の直交座標系から (r,θ,z) の円筒座標系への変数変 換の詳細は省略し結果のみを書くと 密度 ρ と粘度 µ が一定のときの三次元連続の式と運動量方程式は次式になる 円筒座標連続の式 : 1 r r rv r ( ) + z + 1 r 円筒座標 r 方向と z 方向運動量方程式 : v θ θ = 0 (5-9) v r t + v v r r r + v v r z z + v v θ r r θ v θ r = 1 p ρ r + µ 1 ( ρ r r r rv r ) + v r z t + v r r + v z z + v v θ z r θ = 1 p ρ z + µ ρ 1 r r r r + z + 1 r θ + 1 r v r θ v r r + g θ r (5-10) + g z (5-11) r-z 二次元運動量方程式では θ の項と v θ を含む項が無視できる 定常状態で t の項も消える 円管内流動 で 流体が均一分布で入口から入り 入口からしばらく進むと r 方向速度分布 vr を含む項は消え 発達した速度分 布状態で z = 0 であり z の項も消える 最終的に (5-11) 式の z 方向方程式のみが残る 重力の項を無視 すると (5-11) 式は次式のように簡単になる 0 = dp dz + µ 1 r d dr r d dr (5-1) 境界条件は r=0 で d dr = 0 r=r0 で vz=0 である 従って円管内速度分布 : = 1 r v ma r v = r 0 ma 0 4µ を得る Δp Δz (5-13) 一般的な速度分布の場合 例えば密度一定のとき 連続の式 (4-4) 式 運動量方程式 (4-8), (4-9), (4-10) 式の4つの! 方程式を定められた境界条件で解き 三次元速度 v u,v, w ( ) 圧力 p を求めることになる しかし運動量方程式は非 線形であり 例え 次元でも解析的に解くことは不可能である 最近では 大型計算機を使った直接数値解法の手 段が発達しており結果が得られている しかし直接解法だけでは パラメーターの影響あるいは作用する支配的力 等に関する情報を得る事は なかなか難しい そこで 以後に示す無次元化と無次元数を用いて どの項が支配因 子なのか さらには 次元解析を用いて整理する手法が有効となる 5.3 速度境界層流体がある物体の回りを流れるとき 物体近辺の境界壁表面近くと 境界壁から遠く離れた位置とでは 流体に働く要素が異なってくる 例えば Fig.5-3 の翼周りの流れを考えると 翼から十分離れた領域では 粘性の効果が十分に小さいと考えられる すなわち 運動量方程式の右辺粘性項が無視できる 一般に非圧縮性の運動量方程式 (Navier-Stokes 方程式 ) で粘性項を無視した場合の流れをポテンシャル流れ ( 完全流体 ) と呼ぶ 完全流体は理論数 学手段を使いかなり解き明かされており 詳しい内容は流体 力学教科書 ( 例えば 今井功著 流体力学 岩波書店 ) に説 7 Fig. 5-3 物体の周りに発生する境界層

3 明されている 一方翼に近い部分では 固体壁上で速度 0 の境界条件が課されるので 速度が主流速度から小さくなる この領 域を速度境界層と呼び 流れを解析する上で 重要な概念である 速度境界層内では 流体の粘性が無視できない 逆に圧力項は 固体壁に近くなっても境界条件の影響を受けにくいので 境界層内では 圧力は境界層外部の値を 使用する 従って まず圧力分布を粘性の影響を考慮しないポテンシャル流で解析し 求めた圧力分布を既知とし て 粘性の影響を考慮した境界層方程式で解析する ( 境界層近似 ) 手法がよく用いられる 例えば 後で出てくる Fig. 6-4 の平板境界層から十分離れた主流流れ u は 一定で一様と考えられる 従って 重力を無視し非圧縮性 次元定常状態の平板境界層における流れの方程式は次式になる u + v = 0 (5-14) u u + v u = ν " u + u % $ ' (5-15) # & 通常 (5-15) の右辺第一項 ( 方向速度変化の 方向変化 ) は第二項に比べて十分に小さいので 無視される 主流速度が遅いとき 流れは層流であり u, v とも時間的に変化しない この条件で 流れ関数を使う解法があ る 流れ関数 φ を次の様に定義すると (5-14) 式は自動的に満足される u = φ v = φ (5-16) さらに次の相似変数 η と無次元流れ関数 ζ を使って整理する η = u ν φ = u ν ζ (5-17) (5-16) 式に (5-17) 式を代入して整理すると 次の速度式を得る u = u dζ dη v = u ν η dζ dη ζ (5-18) その結果 (5-15) 式は次式に変形される d 3 ζ dη + ζ d ζ 3 dη = 0 (5-19) 境界条件は次のように表せる η=0 で ζ =0 ζ =0 η= で ζ =1 (5-0) ζ を η の多項式で表し (5-19) 式 (5-0) 式に代入して多 項式の係数を順次求める 現在では表計算ソフトを使っ て パソコン上で数知的に求める事が出来る 解いた結 果は Blasius の解と呼ばれ Fig 5-4 に示されている Fig. 5-4 層流境界層内の速度分布 要するに 境界層内速度は η のみの関数として表わせる η 定義式中の Re = u ν ( ) は先の次元解析で現れた無次元数である 流速 粘度等が異なっても Re で整理する 事により 層流境界層内の速度分布がみごとに一つの関係で表す事が出来る これが境界層の相似と無次元数の意 味である 8

4 5.4 無次元化と無次元数 直角座標系における連続の式 運動量方程式 エネルギー式 拡散方程式の独立変数は,, z, t であり 従属 変数は u, v, w, p, T, CA である それぞれ以下のように 代表長さ L 代表速度 U 代表温度差 ΔT 代表濃度差 ΔCA を使って無次元化する 代表長さ L として 例えば内部流れのときダクトの直径を取り 外部流れのとき着目 している物体の長さ等をとる 代表速度 U は外部流れのとき主流速度を取り 内部流れのときは平均流速 um を取 る 代表温度差 ΔT は二つの離れた位置での温度差 (ΔT=T1-T0) を 同様に代表濃度差 ΔCA は二つの離れた位置で の成分 A の濃度差 (ΔCA=CA,1 CA,0) を取る 無次元変数を次のように定義する 時間 :!t = tu L 速度 :!u = u U!v = v U!w = w U 距離 :! = L! = L!z = z L 圧力 : p = 発熱項 :! Q H = Q H L ρc P U ΔT 濃度 :!c A = c c A A,0 発生項: S! Δc A = S L A A U Δc A 以上の無次元変数を使うと 連続の式 運動量方程式 熱収支式 拡散式は次のように変形できる 無次元連続の式 :!u! +!v! +!w!z = 0 (5-1) p ρu 温度 :! T = T T 0 ΔT 無次元 ( 方向 ) 運動量方程式 :!u!u!u!u +!u +!v +!w!t!!!z = 1!u Re! +!u! +!u!z!p! + Gr Re (5-) 非定常項対流項粘性項圧力項浮力項 無次元エネルギー式 :!T!t +!u!t! +!v!t! +!w!t!z = 1!T Re Pr! +!T! +!T!z + Q! H (5-3) 無次元拡散式 :!c A!t 非定常項対流項熱伝導項発熱項 +!u!c A! +!v!c A! +!w!c A!z = 1 Re Sc!c A! +!c A! +!c A!z 非定常項対流項拡散項発生項 ここで Re, Pr, Sc, Gr は次に定義する重要な無次元数である Re = ρul µ Pr = µc P µ Sc = Gr = gl3 k T ρd AB ν (5-5) 以上の無次元化により 次の事が分かる + S! A (5-4) (1) 事なった体系にある異なった流体 熱移動 物質輸送の方程式を無次元化するとすべて同じ式で表される そ れぞれの状況に応じた Re, Sc, Pr, Gr を代入する事で 各場合の移動量を定量化する事ができる () 運動量方程式 エネルギー式 拡散式ともそれぞれの独立変数 u, T, CA を入れ替えると 同じ形式で表される すなわち運動量輸送 熱輸送 物質輸送は類似性 ( アナロジー ) が成立し エネルギー式で現れる現象は 同 様に物質輸送でも現れる (3) Re 数は 流体の慣性力と粘性力の比を表し この比が大きくなると粘性によって流れの乱れを抑制させる働き が弱まる 9

5 (4) 運動量方程式の左辺は非線形であり 流れが速くなると 乱れが大きくなり 層流から乱流に至る Re 数が円管内強制対流では 300 以上になると層流から乱流に遷移する (5) Gr 数は 自然対流の大きさを表す無次元数であり この数が大きくなると 自然対流の大きさが大きくなる 最終的は乱流になる (6) Pr=1 Sc=1 のとき 壁近傍の濃度変化と温度変化は同じ形で変化する 5.5 次元解析 いままで見てきたように 例えば流体中で運動する物体 は 連続の式と運動量方程式 (Navier-Stokes 方程式 ) で Fig.5-5 記述できる しかし実際の体系で この運動量方程式は非線形であり 解く事が非常に難しい そこで系統的に実験をおこない データーを整理し 未知の条件における値を既知条件の値から合理的に求める方法が次元解析である ここでの次元解析とは 現象に関係する物理量間にある関数関係を 次元を持たない無次元数で記述する方法である 無次元数は 現象に関係する物理量の組 み合わせからなり 次元解析を用いる事により 対象とする現象の特徴をうまく 表すように導く事が出来る すなわち 着目する物理現象が複雑でこれを記述す る理論解を求める事が不可能である場合がほとんどである この場合 現象に関 Fig. 5-6 静止流体中を動 く球に働く力 与する物理量 ( 速度 温度 圧力 物性値 系の大きさ等 ) をパラメーターにして 多くの実験をおこない 関係する物理量への依存性をうまく表すように実験結果を整理する 次元解析を用いて 実験によって得られた物理量間の複雑な関係を一般的な実験相関式の形で表そうとするものである (a) 次元解析の例 Fig.5-6 の様に 密度 ρ[kg/m 3 ] 粘性係数 μ[pa s] の静止流体中を速度 v[m/s] で運動する直径 d[m] の球に作用す る抗力 ( 抵抗 )Fd[N] を考える Fd は F d = f d,v,ρ,µ ( ) の関数として表される Fd への d, v, ρ, µ の各物理量の依存性 を調べるため いろいろな状況の Fd 値を実験で求めデーター整理する 次元解析では Fd の一般式が各物理量の塁乗積として表せると考える ( 指数法 ) 従って次式を得る F d = K ρ ( ) a 1 ( µ ) a ( v) a 3 ( d) a 4 (5-6) K は無次元の係数 a1 a4 は未知指数である (5-6) 式が成立するためには 基本次元が両辺で同じである必要が ある 基本次元として 長さ L 質量 M 時間 τ とすると (5-6) 式に含まれる各物理量の次元は 次式で表せる F d = ML τ ρ = M L µ 3 = M τ L v = L τ d = L (5-6) 式の両辺で 3 個の基本次元が一致するためには 次式が成立する必要がある 30

6 L の指数について : 1 = 3a 1 a + a 3 + a 4 M の指数について : 1= a 1 + a τ の指数について : = a a 3 4 個の未知数に対して 方程式が 3 個であるので a を使って残りの未知の指数を表すと a1=1 - a a3= - a a4= - a を得る この関係を再度 (5-19) 式に代入し 整理すると 次式を得る F d ρv d = K ρvd µ a (5-7) 球の投影面積 A(=πd /4) 流れの動圧 ρv / を使って ( 無次元 ) 抗力係数 Cd を次式で定義する F d C d = A ρv (5-8) (5-7) 式の右辺は粘性力 µv d に対する慣性力 ρv の比を表す無次元数で 前の Renolds 数である Re = ρvd µ = vd ν (5-9) (5-7) 式に (5-8), (5-9) 式を考慮すると 抗力係数は 次の様な Renolds 数のみの関数で表せる事が予測できる C D = f Re ( ) (5-30) 過去多くの異なった流体 流速 形状の条件で集められた球の抗力係数データーをまとめたものが Fig.5-7 である データーが1 本のカーブ f(re) で相関でき 次元解析を用いて得られた無次元数の関数関係から 異なる条件での抗力係数を求める事が出来る Re 数の小さい領域は Stokes 流領域であり CD=4/Re の関係で表せる 続いて Re=1 付近から球背後に渦を発生し始め Re=100 付近で周期的な Karman 渦が発生する そのため 抗力係数は Stokes 流の仮定からはなれ CD は大きくなる その後 Re>1500 になると球背後の流れは乱流になり 大きな渦が発生し CD=0.6 付近でほぼ一定の状況が続く さらに Re>310 5 以上になると 球表面での流れ も乱流に遷移し CD は急に小さくなる Fig. 5-7 球の抗力係数の無次元整理 流れはいくら複雑になっても 当初の基 本次元と関係する物理量に変化がなく CD が Re 数で整理できることに変わりはない (b) バッキンガムの π 定理 次元解析では 現象に関係する因子 ( 物理量 ) を組み合わせた無次元数の種類を知る事が出来るが 一般には 31

7 次のバッキンガムのπ 定理によりその個数が規定される すなわちある現象に関与する物理量が φ1, φ, φ3, - - φn の n 個あり これらの間に次式のような関係があるものとする ( ) = 0 (5-31) f 1 φ 1,φ,φ 3,,φ n このとき 物理量 φ1, φ, φ3, - - φn を表すのに必要な k 個の基本次元が必要であるとすると この物理量は m=n-k 個の無次元数 π1, π, - -,πm を用いて次のように表す事が出来る これがバッキンガムのπ 定理である f ( ) = 0 (5-3) π 1,π,,π m 静止流体中に動く球に対する抗力係数の次元解析では 抗力を含む物理量の数は n=5 基本次元の数は m=3 であり 独立な無次元量の数は k= である 基本次元は 基本単位と呼ぶこともあるが ここでは次元と単位を区別して扱う 次元が ある物理的な挙動や性質の特性を表す 長さ 質量 時間等の物理特性であるのに対して 単位は [m], [kg], [s] 等の次元の大きさを表す尺度の基準である 次元解析で用いる基本次元には 力学現象において現れる物理量に対して 長さ L 質量 M 時間 τの次元を 伝熱現象には熱量 H 温度 θ 物質移動現象には 濃度 C を独立次元として加えて用いる しかしながら 熱移動における熱量は H=ML /τ であるから L, M, τ およびθを基本次元として次元解析をおこなう事が出来る ( 補講 5) 乱流応力 乱流において 流れは不規則に時間的に変化しながら流れる 任意の瞬間における,,z 方向速度 v, v, vz と圧 力 p を 時間平均値 v, v,, p と変動値 v,, p に分け次式で考える v = v v = v = + p = p + p (5-33) 微小時間 dt 内に 軸に直角な微小面積 da を通過する流体の質量は ρv dadt であり 方向運動量流束は J = ρv dadt である JX の時間平均を取ったものを dj で表すと 簡単のため密度一定の仮定のもとで 次式のよ うに表す事が出来る dj = ρv da = ρ( v ) da (5-34) 同様に, 方向と z 方向の運動量流束の時間平均 dj,dj z も次式で表せる dj = ρ( v v )da = ρ( v v v )da (5-35) dj z = ρ( v )da = ρ( v )da (5-36) 各時間変動分の時間平均について v = v = = 0 であるが (5-34)-(5-36) 式右辺第二項中に現れる二つの方向の変 動成分を乗じたものは 時間平均をとっても 0 にはならない これらは 軸方向に垂直な面の単位面積に働く力に なる 方向の垂直応力が ρ v ( ) あるいは z 方向の剪断応力が ρ v v + ( v v ) ρ v ( ) である 夫々の第一項は 3.1 節で示した粘性応力に相当する 第二項は乱れがあるときに現れる応力で 乱流応力あるいは Renolds 応力と呼ばれる 結果を 3 次元に拡張すると 3 次元の乱流応力は次式のように書ける 3

8 σ,τ,τ z v v τ,σ,τ z = ρ v v τ z,τ z,σ z v, v, (5-37) Prandtle は 乱流においても分子の平均自由行程に対応する混合距離 l が存在し 混合距離離れた位置間で 運動 量が輸送されると考えた その理論によれば 例えば 乱流応力成分 τ は次式で表せる τ = ρ v v = ρε v M (5-38) 導入された混合距離は場所や対象としている流れによって変化するが 平均速度に依らないとされている 式中の ε M は渦動粘性係数であり これに適切なモデルを導入し N-S 方程式中で計算される場合がある ε M = v (5-39) 33