乱流構造に基づくSGSモデルを用いた 微気象場の数値計算手法の開発

JpGU 2015
乱流構造に基づくSGSモデルを用いた
微気象場の数値計算手法の開発
Development of computational method for
micrometeorological field using SGS model
based on turbulence structures
May 25, 2015
京都大学防災研究所
気象水文リスク情報(日本気象協会)研究分野
井上 実
DPRI, Kyoto Univ.
INOUE, Minoru
1
Purpose
非等温場における熱や物質拡散は局所的な集中豪雨や大気汚染に係わ
る重要な現象であり、このような現象を理解することは集中豪雨の早期予測
や大気汚染の緩和に役立つものと考えられる。
時空間スケールの細かい複雑な現象
LESが有効な手段の一つ
・本研究では、非等温場における微細な乱流構造を捉えることができる
LESモデルの開発を目的とする。
Development of a computational method for the LES
・微気象場を扱うことができるように、流体の運動方程式の他に、熱や
水蒸気の輸送方程式および凝結過程を導入する。
Investigation of the turbulent transport of heat or vapor in the
micrometeorological field
・SGSの乱流モデルには、Kobayashi(2005)が提案した乱流構造に基づく
新しい手法を適用する。
Coherent-structure Smagorinsky Model
ここでは、熱や水蒸気輸送に対する本手法の妥当性確認の結果を示す。
Demonstration of the validity of the present method for the
turbulent transport of heat or vapor
2
Computational Method
Characteristics
・Conservation law of mass, momentum, heat and total water
・Generalized curvilinear coordinate system
・Finite volume method
Governing Equations(filtered)
Mass
Momentum
 n  u dS  0
S

udV   n  TdS    g  v   0 dV
S
V
t V

T   pI  uu  (   t ) u  (u)Tn
Liquid water potential temperature

: coefficientofthermalexpantion
: acceralationofgravity
: virtualpotentialtemperature
1 0.61
: referencepotentialtemperature
: turbulentviscosity
and total water specific humidity

dV   n   u     t  dS
S
t V
,
:turbulentPrandtlnumber
3
Condensation Process
Condensation Process
タイムステップ毎の水蒸気量(比湿)が飽和に達しているか否かを
判断し、飽和に達している場合は多い部分の水蒸気を凝結させ、未
飽和の場合は飽和比湿に満たない雲水を蒸発させる。その際、潜
熱による気温の上昇、下降を考慮する。
Condensation only occurs if the mean state of the grid box becomes
over saturated.
: Watervaporspecifichumidity,
condensation
liquid water
temperature rise
saturation
: Liquidwater
evaporation
temperature fall
liquid
water
unsaturation
Image of condensation process
4
Turbulence Model
Coherent-structure Smagorinsky Model (CSM, Kobayashi, 2005)
∆
:Smagorinskycoefficient
Δ: filterwidth
2
Where,
,
0.05
: secondinvariantnormalizedbythe
magnitudeofavelocitygradienttensor
1
2
rotation
1
2
strain
1
2
1
2
0isinterpretedasavortextube
5
Validation of Turbulent Transport of Heat and Vapor
LES of Shallow Cumulus Convection
Siebesmaら(2003)はカリブ海沖で発生した境界層積雲を対象とし、10の研究機
関で開発されたLESモデルの比較実験を行った。ここでは彼らと同様のLESを
実行し、本手法の熱や水蒸気輸送に対する妥当性を確認する。
Comparison between the computational results of shallow cumulus convection using this
method and the intercomparison study with 10 LES models by Siebesma et al. (2003)
Computational Condition
Initial conditions
Computational domain
・6400×6400×3000m (cyclic)
Forcing terms
・Coriolis force
・large-scale subsidence
・moisture tendency
・radiative cooling
Boundary conditions
・cyclic boundary condition
・momentum flux
・sensible heat flux
fixed values
・latent heat flux
6
(Siebesma et al., 2003)
Iso-surface of Liquid Water ql
Trade wind
Iso‐surface of liquid water ql(1.0×10‐6) and
contour map of vertical velocity component w (t=6h)
7
Vertical Profiles of ql and Cloud Cover (t=5~6h)
2500
liquid water
ql
Siebesma
et al.
present
height(m)
2000
1500
1000
500
0
0
0.002
0.004 0.006
ql(g/kg)
0.008
0.01
2500
cloud cover
cloud
core
height(m)
2000
1500
1000
500
0
0
0.02
0.04
fraction
0.06
0.08
8
Turbulent Flux Profiles of
and
(t=3~6h)
2500
total water
Siebesma
et al.
present
height(m)
2000
1500
1000
500
0
0
25
50
75 100 125
w'qt'(W/m 2)
150
175
2500
Liquid water
potential
temperature
height(m)
2000
1500
1000
500
0
-40
-30
-20
-10
2
w' θ'l(W/m )
0
10
9
Cumulus Convection over Mountain Top
Computational domain
・3000×3000×3000m (cyclic boundary)
・mountain height = 300m
Initial condition
・u, v & w = 0.0 m/s
・θ & qv : Siebesma et al.
・surface potential temperature = +5K
Boundary condition
・Louis equation (bulk method)
・roughness height = 0.05m
・evaporation efficiency rate = 0.3
Forcing term
・Coriolis force
Computational domain
10
Time Sequence of Cumulus Convection
t=9000
~12000s
気塊の上昇とともに積雲
が発達し、逆転層で雲水
が蒸発し下降気流が生
じる様子が捉えられた。
Velocity vector field and contour map of , line : cloud shape
11
Turbulence structure around Cumulus
t=9000
~12000s
Q>0 : vortex tube
積雲の発達に伴い先端
付近に渦輪が形成され
る。また、積雲が消滅す
る際にも渦輪が外側に
広がりながら崩壊する様
子が捉えられた。
Iso‐surface of Q (=0.0001, yellow) and ql (=0.00001, white)
12
Conclusions
・非等温場の乱流による熱や物質拡散の解析を目的として、
熱や水蒸気輸送を伴う微気象場のLESモデルを開発した。
SGS乱流モデルには乱流構造に基づく新しい手法を適用した。
・基本的な乱流場や境界層積雲の問題に本モデルを適用し、
他の実験値や計算結果と概ね一致する結果が得られた。
また、非定常な積雲対流中の微細な乱流構造を捉えること
ができた。
・以上のことから、本計算手法を用いた微気象場の乱流解析
に対する可能性が示された。
A computational method for the LES in the micrometeorological field was developed,
and some numerical experiments were carried out to demonstrate the validity of this
method. The computational results suggested the possibility of this method for the
analysis of the turbulent flow in the micrometeorological field.
13