Large-eddysimulationsofthewind-inducedturbulentEkmanlayer.pptVIP

Direct numerical simulation ofturbulent channel flowunder stable density stratification Iida, O., N. Kasagi and Y. Nagano Int. J. Heat and Mass Trans. (2002) 45:1693-1703 * 水圏地球物理学ゼミナール 論文紹介 2006.11.1 1 先行研究 ○乱流に関する安定成層の効果 　1.乱流を抑制　→再層流化の場合も　2.内部重力波(IGW)の励起 (Tritton 1988, Gossard and Hooke 1975,Hooke and Jones 1986)　　　IGWは運動量を輸送する　　　しかし、線形の範囲では、鉛直流速(v)と擾乱水温(θ)の位相差がπ/2　　　　→熱は輸送しない (Stull 1988, Webster 1964) ○成層の下での境界層乱流に関する多くの水槽実験と観測 　勾配リチャードソン数Rigによって流れの状態が決定　　0.25以下の領域：Kelvin-Helmholz不安定による間欠的乱流 (Tritton 1988, Turner 1973)　　0.25以上の領域：IGWが支配的 　2つの領域間で相互作用　圧力相関項を通して、IGWが乱流境界層からエネルギーを運び出すと示唆(Stull 1988)　　ただし圧力擾乱の測定は困難　→　数値実験の必要性 1 目的 Channel(平行平板間)乱流について一連の直接数値実験(DNS)を行い、安定密度成層が境界層乱流に与える基本的な効果を調べる 　特に、乱流とIGWの相互作用、　　　　　再層流化に注目 目次：　2　数値モデル　3.1　乱流熱?運動量輸送　3.2　乱流構造　3.3　乱流運動エネルギーと乱流熱フラックスの収支 2 Computational conditions of DNS X (cyclic) y z (cyclic) 一定の圧力勾配 実験ケース：　バルクのリチャードソン数：0～0.54 粘着壁、温度T0 粘着壁、温度T1 支配方程式：　ブシネスク近似の下での通常の流体力学方程式 3.1 抵抗係数、Nusselt数 中立成層実験に対する比 3.1節：乱流熱?運動量輸送に対する安定成層の効果 バルクリチャードソン数 解像度は十分水槽実験とも一致 Rib大　→　Cf,Nu低下 3.1 相関係数(uとv,vとθ) 実験(×):ばらつきがあるが傾向は一致 成層強(Rig大)　→　相関低下　(乱流運動量?熱輸送低下) 先行研究とよく一致Rig0.1でばらつくのは粘性の影響 3.1 平均流速U、温度T分布(対数軸) ( )+:摩擦速度、摩擦温度、粘性係数による無次元化 y+ 10: 成層強→平均流速大　　　　　　　　　　　　対数領域小 変化は平均流速と似ている　y+ 10: 成層強→温度勾配大 channel中央で変化最大 3.1 レイノルズ応力、乱流熱フラックス 壁近傍を除く、全領域で低下 channel中央で大きく減少　→Tの大きな変化 3.1 等温度面(瞬間場) 乱流混合が抑制 Rib=0 Rib=0.35 大きく波うつ　→混合大 3.2 流速ベクトル、圧力(瞬間場) 3.2節：乱流構造に対する安定成層の効果 壁 壁 壁 壁 黒：低圧 青：低圧 streamwise vorticeに対応 (Kasagi et al. 1995) Rib=0.35 低圧部：水平に広がる 3.2 流速、温度の分散 u v,w u,v,w (成層強) θ 壁近傍では乱流維持 IGW 成層による抑制→v小 dU/dz大→u大 Rib=0.54で再層流化 変化最大←浮力の影響大 3.2 擾乱流速のスペクトル 成層：中央を除いて、低波数帯で低下　→大スケール擾乱の弱化 　←v+2を低下させる浮力拡散項vθのエネルギーは低波数に集中 中央ではk=1~4の範囲でエネルギーが増加　←IGWの発生 3.2 vとθの位相差 vとθの位相差：π/2　→　IGWと言える y+=10,100では位相差～0　→IGWではなく乱流が支配 3.2 勾配リチャードソン数 流れの構造(乱流かIGW)はRigで判断可能 乱流 IGW 3.3 TKE収支 3.3節：TKEと乱流熱フラックスの収支に対する安定成層の効果 浮力の直接の影響 Viscous diffusion以外 主に生成と散逸でバランス　浮力散逸は小さい 逆では？ 生成項減少→散逸項とバランス 生成項増大 (dU/dｙ大)→乱流拡散とバランス 3.3 TKEz成分の収支 浮力散逸項 →圧力歪み相関項の