Passive scalar convection and mixing in turbulence

Nagoya Institute of Technology
中部CAE懇話会 乱流数値解析講座
乱流現象と基礎理論
名古屋工業大学 大学院
創成シミュレーション工学専攻
後藤俊幸
Nagoya Institute of Technology
流れの方程式の特徴
運動量保存則
(粘性流体)
(Navier Stokes 方程式)
流れる
圧力
粘性
1.非線形
レイノルズ数 =
Nonliner term
Viscous term
=
2.連続体 (無限大自由度)
3.粘性


数学的解析が著しく困難 (限られた数の厳密解があるだけ)
同じ初期、境界条件のもとでも解は1通りでない(解の安定性)
Nagoya Institute of Technology
乱流の特徴
1. 微小撹乱に対する不安定性
2. 強い非線形性
Nonlinear term
=
Viscous term
3. 巨大な輸送能力
Turbulent viscosity
=
Molecular viscosity
4. 無限大自由度
5. 散逸系
をもつランダムな流れ場
Nagoya Institute of Technology
十分発達した乱流の基礎的理論
Nagoya Institute of Technology
乱流研究の巨人達 1
A. N. Kolmogorov
O.(1903-1987)
Reynolds
(1842-1912)
L. Prandtl
(1875-1953)
T. Von Karman
(1881-1963)
G. I. Taylor
(1886-1975)
Nagoya Institute of Technology
乱流研究の巨人達 2
A. N. Kolmogorov
(1903-1987)
L. Onsager
(1903-1976)
L. D. Landau
(1908-1968)
S. Chandrasekhar
(1910-1995)
Nagoya Institute of Technology
乱流研究の巨人達 3
J. von Neumann
(1903-1957)
R. P. Feynman
(1918-1988)
G. K. Batchelor
(1920–2000)
R. H. Kraichnan
(1928-2008)
Nagoya Institute of Technology
Nagoya Institute of Technology
Kolmogorov の乱流理論 (1941, K41)
Onsager, von Weizeker, Heisenberg
乱流における各スケールでの統計を考える
乱流理論の一里塚
Velocity increments
v
u
r
Nagoya Institute of Technology
2つの仮説と次元解析
1. R >> 1 でかつ r << L のとき
U と V の確率密度関数は n
n
0
のとき
と
,
乱流場は局所的に等方的で、
e で一意に決まる
e
は有限
単位質量あたりの
平均エネルギー散逸率
動粘性率
h= (n / e )
3
1/4
Kolmogorov 長さ
最も小さいスケール
2. h << r
のとき, 確率密度関数は粘性によらない
Nagoya Institute of Technology
V
n   のとき
抵抗
e 有限
Cd : 抵抗係数
S
パワー
Fluid dynamics for physicists
Faber (Cambridge 1995)
Nagoya Institute of Technology
Sketch of Turbulence
Energy input
ein
L
P
eout
h
Scale
Heat
Nagoya Institute of Technology
慣性領域
h << r << L
PDF
(Probability Density Function)
Moment
Normal scaling
Kolmogorov’s Spectrum
Kolmogorov’s 4/5 Law
(Onsager 1945)
(asymptotically exact)
Characteristic time
Nagoya Institute of Technology
十分発達した乱流の解析的統計理論
Nagoya Institute of Technology
フーリエ空間でのNS 方程式
圧力は非圧縮条件より消去
外力
Lbox >> L
Lbox
(乱流の巨視的な長さ)
Nagoya Institute of Technology
速度揺らぎのエネルギースペクトル
k3
k
dk
k2
k1
Nagoya Institute of Technology
乱流の解析的統計理論
完結の問題 (Closure problem) 移流項, 圧力項, 粘性項
NS eq.
Moment eq.
PDF eq.
Conditional average of
the Lagrangian acceleration
Random forcing
Nagoya Institute of Technology
乱流のスペクトル理論 (たとえばE(k,t) を計算すること)
Eulerian theory
•Quasi Normal
(Milionshtchikov 1941, Proudman and Reid 1954,Tatsumi 1957)
•Direct Interaction Approximation (Kraichnan 1959)
•Modified Quasi Normal Markovian (Tatsumi, Kida and Mizushima 1978)
Phenomenologcal Theory
•Eddy Damped Quasi Normal Markovian (Orszag 1966)
•Test Field Model (Kraichnan 1971)
Lagrangian Theory
•Lagrangian History Direct Interaction Approximation (Kraichnan 1965)
•Lagrangian Renormalized Approximation (Kaneda 1981)
• ラグランジュ的物理量を用いて方程式を閉じている
• 任意定数を含まない
• システマティック
Nagoya Institute of Technology
乱流の解析的統計理論
Quasi Normal Theory (準正規理論)
(Milionshtchikov 1941, Proudman and Reid 1954,Tatsumi 1957)
Moment eq.
Nagoya Institute of Technology
Key players
Average rate of the energy dissipation per unit mass
Equation of the energy spectrum
P (k)
Energy transfer flux
k
Nagoya Institute of Technology
Kolmogorov スペクトル
K : Kolmogorov constant
=1.62 (Exp. Sreenivasan 1995)
=1.64 (DNS, Gotoh et al. 2002)
=1.72 (LRA, Kaneda 1986)
スペクトル理論の成果
•
•
•
•
•
エネルギーカスケードの存在
乱流粘性の一般的表現
乱流拡散
ラグランジュ的およびオイラー的速度相関
スカラー分散のスペクトル
辻(2002)
Tsuji 2002
Nagoya Institute of Technology
log E(k)
乱流のエネルギースペクトル
連続スペクトル
Inertial range
Energy
containing
range
ein
in a steady state
~
ein = P =eout
~
P
Dissipation
range
eout
log k
kL=1/L
Energy input
kd=1/h
Heat
Nagoya Institute of Technology
エネルギースペクトルには切れ目がない
スケールの分離がない
あるスケール kc 以下(以上)の自由度の運動のモデル化の困難
kc に依存したモデとなって普遍性が弱くなる
予測と制御の不確かさ
E(r)
スケールの分離がある例
ミクロスケールの現象
粘性、拡散など
巨視的な現象
r=kc-1
r
Nagoya Institute of Technology
十分発達した乱流と計算科学
Nagoya Institute of Technology
問題点
解析の困難さ
理論に導入された仮定
誤差評価
適用限界
強い非線形性
散逸系
集団平均測度が不明
シミュレーション
スケールの分離がない!
すべての乱流スケールをまるごとシミュレートすることが必要
一切のモデル化の排除
乱流の直接数値計算 (Direct Numerical Simulation, DNS)
乱流のDNS は解像度が命
できるだけシンプルな乱流を選ぶ
より簡単な問題
バーガース乱流, Kraichnan Model (Passive scalar )
Nagoya Institute of Technology
Canonical Turbulence
乱流の基礎物理の解明に向けて、
1.できるだけ単純な境界条件
計算機資源のほとんどを時間・空間解像度に投入
2.研究の基準となるさまざまなデータを供給
3.理論の構築や検証に対する規範
一様等方性乱流
一様せん断乱流
平行平板間流れ
Nagoya Institute of Technology
一様等方性乱流
(Homogeneous Isotropic Turbulence, HIT )
• 無限に広がった空間内において、統計が空間並進と回転に対して不変な乱流
• 統計物理量の独立変数、従属変数が少ない
• 普遍性への期待(境界条件や初期条件によらない)
continuum, convection, pressure (incompressibility )
Simple Model
Real life flow
Lorenz
Shell model
Cascade model
(Multi Fractal,
Log Poisson)
…….
Poiseullie flow
Heat convection
Airplane
Atmosphere
Ocean
HIT
…….
Complexity
Nagoya Institute of Technology
Taylor Micro Scale
l
L
r
L : 乱流の巨視的スケール
l : Taylor microscale
Taylor microscale Reynolds number
Nagoya Institute of Technology
乱流のDNS は解像度が命
最小の乱流渦まで正しく解像することが必須
乱流の自由度
L : 乱流の巨視的スケール
ジェット旅客機
Re=10 8 N~10 18
台風
Re>10 12 N>10 27
Nagoya Institute of Technology
乱流の計算科学
一様等方乱流のDNS
スペクトル法
周期的境界条件(統計的一様性)のもとで NS 方程式を数値的に解く
フーリエ級数に展開
外力
Lbox >> L
Lbox
(乱流の巨視的な長さ)
Nagoya Institute of Technology
スペクトル法
長所
微分量の高精度評価
ポアッソン方程式の高精度解
(非圧縮条件あるいは質量保存則を高精度で満たす)
非線形項のエネルギー保存に優れている( k < kmax )
短所
の計算負荷が大きい
Nagoya Institute of Technology
直接和
計算量膨大
スペクトル法
u(k)
F [u(x)u(x)]
FFT (Fast Fourier Transform)
F- -1
F
高効率のFFTが重要
u(x)
u(x)u(x)
Nagoya Institute of Technology
乱流の DNS の発展
10000
4000
Peta Machine
460
170
1000
vpp5000/56 NUCC
Caltech Delta machine
1200
Rl
Cray-2 CCVR
Earth Simulator
100
vpp500 NUCC
Cray-1S NCAR
40
IBM Model 360-95
Cray-1 NCAR
280
2/3
R l ∝ Nx
10
10
100
1000
10000
100000
Nx
「乱流の謎にせまる計算科学」 後藤、石原 (パリティ 2002年10月号)
Nagoya Institute of Technology
一様等方性乱流
Rl =732
|w |> |w| + 4s
Kaneda and Ishihara ( JOT, 2005)
乱流の大規模DNSで
見えてきたこと
階層構造をもつ
渦管の大集団
乱流の統計力学
l
10l
L
10h
Nagoya Institute of Technology
乱流による温度や物質の輸送
Watanabe and Gotoh 2006
Nagoya Institute of Technology
T=<T>+q = Gz+q
G
Rl= 468, Sc=1, N=1024 3
Nagoya Institute of Technology
Scalar dissipation
Rl= 468, Sc=1, N=1024 3
Nagoya Institute of Technology
Great wall in the universe
Sketch of Pleiades cluster
by Galileo (1610)
Galileo’s telescope (1610)
de Lapparent et al. (1986)
http://www.nhk.or.jp/school/junior/yougo26.html#010
Subaru Telescope (1999)
Nagoya Institute of Technology
規格化されたエネルギー散逸率
Doering and Foias (2002)
eL
b=
3
urms
Donzis et al. (2005)
Doering and Foias (2002)
Nagoya Institute of Technology
Compensated Energy Spectra
N=4096 3
N=1024 3
(Kaneda, Ishihara, Yokokawa, Itakura, and Ueno, 2002)
Nagoya Institute of Technology
スケールごとのエネルギー
KHK equation
4/5 law
(Gotoh et al. 2002)
•KHK 方程式(エネルギー収支)は満たされている
•Rl の増大と共に 4/5 法則へゆっくりと漸近
Nagoya Institute of Technology
間欠性
Kolmogorov 理論のやぶれ
Nagoya Institute of Technology
Intermittency
deviation from Kolmogorov (K41) theory
Structure functions of velocity increments
Kolmogorov (K41)
DNS (Gotoh et al. 2002
Nagoya Institute of Technology
Kolmogorov 理論によるスケーリング指数 zp=p/3
(Gotoh et al. 2002)
(Chen et al 2005)
Kolmogorov 理論 とは異なるスケーリング
?
Nagoya Institute of Technology
間欠性 (Intermittency)
u(t)
PDF
t
du(t)
dt
PDF
t
大気接地層(地上7m, Rl=3300 )
Praskovski (1997)
Nagoya Institute of Technology
PDF of longitudinal velocity increment
N=2048 3
Rl =585
r decrease
r decrease
In Kolmogorov Theory
PDF の曲線は重ならず、スケールに依存して変化
分布関数はスケールの減少とともにどのように変化するか?
Navier Stokes方程式から定量的に計算できるか?
普遍性はどこにあるのか(中心部分それとも裾野)?
h << r << L
Nagoya Institute of Technology
Particle trajectory and acceleration in turbulence
Low
High
1,500 times the acceleration of gravity (equivalent to 40 s )
La Porta et. al. Nature 409, 1017-1019 (2001)
Nagoya Institute of Technology
PDF of acceleration
Experiment Rl = 690
Mordant, Crawford and Bordenschatz
Physica D (2004)
DNS Rl = 381
Gotoh and Kraichnan Physica D (2004)
Nagoya Institute of Technology
経済活動における特異なゆらぎ
Log (P/<P>)
Nagoya Institute of Technology
Energy flux
Forward
Average flux
Reverse
Wavenumber
k=1/r
Pr
: エネルギー輸送流束の大きな揺らぎと小さな平均値
Navier-Stokes の動力学(慣性力)が直接関わっている
Cf. Volume averaged dissipation
er
Nagoya Institute of Technology
乱流の計算科学
大きな特徴

実験では得にくい物理量が計算可能
圧力場、圧力勾配場(加速度)、 .....

可視化に優れている
渦度場、散逸場、圧力、濃度、 …..

仮想実験が可能
場の外科手術、方程式の各項の変形、次元の影響、…..
いろいろな乱流現象にも
パッシブスカラー(ベクター)乱流、対流、超流動乱流、…..
Nagoya Institute of Technology
乱流理論と計算科学 のこれから
統計力学、熱力学: 常にNA~10 23 の極限状態
もしNA~10 3 だったら?
乱流: R→∞ の漸近状態はめったに実現しない
大気、海洋、宇宙など
Kolmogorov 理論、これまでの統計理論
その多くがR→∞の極限での理論
乱流の HITのDNS はようやく漸近状態が
見えつつあるところまで来た しかしまだ不十分
Rl∝N 2/3 でしか増大しない
Nagoya Institute of Technology
もしのパワーをもつ計算機があったら?

本当に漸近状態はあるのか? あるとすれば

乱流の漸近状態の現象論
スペクトル、分布関数、応答、非等方性、エネルギー輸送、…

どれほど普遍的なのか?
記述のための物理理論と要求精度との兼ね合い
巨大なデータに埋もれる?
埋もれないためには
個々のスケールやパターン
(たまねぎの1枚1枚の皮) を見る眼と
全体を(たまねぎの全体)見渡す理論の構築
乱流の統計力学と巨視的理論
漸近状態での乱流理論 から 有限状態での乱流理論へ
Nagoya Institute of Technology
まとめ

Canonical Turbulence DNS を用いて乱流の基本的問題にアプローチ

Kolmogorov 理論の検証が進んでいる
スペクトル、4/5 法則はほぼOK, 高次モーメントはNO

スペクトル理論はある程度成功しているものの十分ではない
詳細な検証が必要

小さいスケールでの乱流場の統計は Kolmogorov 理論からずれる
間欠性の問題はいまだ解かれていない
普遍性、漸近状態(n  )

乱流のDNSにより乱流場の詳しい情報が得られるようになった
近似的理論の構成に役立てることができる
より大きなレイノルズ数の乱流場の実現が求められる

場の特異性や空間構造や配位と統計法則との関連が重要

仮想実験と乱流 DNS の有効性
Nagoya Institute of Technology
工学への応用と実用に向けての
乱流の数値シミュレーション
Nagoya Institute of Technology
数値シミュレーション
工学や物理現象解明のための重要な手法
詳細で信頼性の高い天気予報
環境・エネルギー問題
経済的な製品開発
危機予測と回避
………
高解像度とより忠実な物理プロセスの導入
信頼性のあるモデル化
Nagoya Institute of Technology
流れの方程式の数値シミュレーション
時空間をサイズ Dx, Dt のメッシュに分割して微分を表現
y
D
x
Nagoya Institute of Technology
流体力学における数値計算
Richardsonの夢、Neumannの挑戦
L. F. Richardson (1922)
「数値的手法による天気予報」
流体力学の式を使えば天気用法が出来る(はず)
手回し計算機
1920年5月20日の天気図
6時間後のミュンヘンでの気圧変化
145hPa!
失敗の理由
大きなスケールの高低気圧による変化と乱流(小さな渦)による変化が
区別できなかった (乱流による変化分>高低気圧による変化分)
大正12年3月18日 虎印計算機械 第一号
大本寅治郎(1887-1961)
http://www.tiger-inc.co.jp/temawashi/temawashi.html
Nagoya Institute of Technology
J. G. Charney,E. Fjortoft,J. von. Neumann (1946~1956)
Princeton, IAS(Institute Advanced Study)
ENIAC
24時間後の500hPaの高度面の予測に成功 (15x18 格子点)
成功の理由
大規模な流体運動のみの方程式を流体力学より導出した
(準水平運動、準地衡風近似によるノイズ(乱流成分)の除去)
気象現象(流体運動)におけるスケールの概念の導入
ENIAC
真空管17468本
メモリ 20個の変数 と300個の定数
幅24m、高さ2.5m、奥行き0.9m、
総重量30トン,消費電力は150kW
Wikipedia より
Nagoya Institute of Technology
乱流現象の謎を大規模コンピューターシミュレーションにより解明
たいていの流れは乱流 しかし、 いまだ解けていない大問題
地球: 大気と海洋
地球シミュレータによる海洋のシミュレーション
Nagoya Institute of Technology
Earth Simulator による太平洋の50年積分
http://www.jamstec.go.jp/esc/gallery/index.html#simulation_movie
Nagoya Institute of Technology
log E(k)
Inertial range
~
P
Energy
containing
range
grid scale
(GS)
in a steady state
~
ein = P =eout
Dissipation
range
subgrid scale
(SGS)
ein
eout
log k
kL=1/L
Energy input
Kc=p/dx
kd=1/h
Heat
Nagoya Institute of Technology
乱流のモデリング
乱流のすべてのスケールを解像する
コスト高い
dx = kc-1
Large Eddy Simulation (LES)
巨視的な乱流運動のみをシミュレート
格子サイズ以下の乱流運動をモデル化
grid scale
(GS)
subgrid scale
(SGS)
Nagoya Institute of Technology
Filtering
GD
x
GS
SGS
r
Nagoya Institute of Technology
乱流のモデリング
どのように LES のよしあしを判定するのか?
?
DNS はLES モデルの開発に必要不可欠
Nagoya Institute of Technology
理論的整合性のある乱流モデルの開発
乱流理論、乱流の大規模DNS
乱流モデルの構築
• 定性的、定量的評価
• 誤差の由来と適用限界
Nagoya Institute of Technology
どのようにLESのよしあしを判定するのか
アプリオリテスト
モデルの直接の比較
?
アポステリオリ テスト
乱流場の特性についてDNS と LESを比較
1.エネルギースペクトル、レイノルズストレス、分布関数、
時間相関
…….
2.速度場のパターン
3.誤差成長
|uDNS - uLES|
DNS はLES モデルの開発に必要不可欠
Nagoya Institute of Technology
参考文献
「乱流理論の基礎」 後藤俊幸
朝倉書店 1998
特集 「乱流」 パリティ 丸善 2002年10月号