Nagoya Institute of Technology

Nagoya Institute of Technology
中部CAE懇話会
流体伝熱基礎講座
第3回 午後
名古屋工業大学大学院
創成シミュレーション工学専攻
後藤俊幸
Nagoya Institute of Technology
粘性流体
y
U
H
F
A
u(y,t)
-F
x
単位面積当たりのせん断応力
Nagoya Institute of Technology
Newton 流体
t
線形関係
Nagoya Institute of Technology
応力テンソル
t ij
力
力の方向 面(法線)
z
n=(0,0,1)
t zz t
y
n=(0,1,0)
t zy
t yy
yz
t xz
t xy
t zxt
yx
t xx
n=(1,0,0)
x
Nagoya Institute of Technology
一般の線形関係は
しかし
流体は等方的と仮定
変形速度テンソル
(rate of strain tensor)
回転は寄与しない
非圧縮性流体
Nagoya Institute of Technology
圧縮性も考慮した粘性流体の応力テンソル
非圧縮性粘性流体の応力テンソル
Nagoya Institute of Technology
t xy
流体に働く粘性力
dy
t xx
o
2点間の力の差が働く
非圧縮性条件より
dx
Nagoya Institute of Technology
粘性流体の運動方程式 (Navier Stokes方程式) と連続の式
3次元で書き下すと
Nagoya Institute of Technology
初期条件
境界条件
境界上で
Nagoya Institute of Technology
重要で簡単な厳密解
A.定常な平行流
A1.クエット流れ
A2.平面ポアズイユ流れ
A3.円管ポアズイユ流れ
B.非定常な平行流
B1.レイリーの流れ
Nagoya Institute of Technology
定常な平行流
L
y
H
P0+Dp
U
u(y)
P0
x
流れに平行な方向に
x 軸をとる
Nagoya Institute of Technology
0
0
0
0
0
0
Nagoya Institute of Technology
U
U=0
Dp=0
x
x
Couette flow
Poiseuille flow
U
x
流量
Nagoya Institute of Technology
Hagen Poiseuille流れ
r
z
R
流量
Nagoya Institute of Technology
非定常な平行流
レイリーの問題
y
u(y,t)
x
U0
Nagoya Institute of Technology
Nagoya Institute of Technology
無次元化と相似則
流体力学の実験はどうやって行う?
r U
L
Nagoya Institute of Technology
無次元化された Navier Stokes 方程式
レイノルズ数
Nonliner term
= Viscous term
非線形性の強さ
レイノルズの相似法則
境界の幾何学的形状が相似でかつレイノルズ数が等しい流れは
力学的に同等である
Nagoya Institute of Technology
水
空気
例1 時速4km/hで歩く人のまわりの流れ
例2 1mm/secで泳ぐプランクトンのまわりの流れ
Nagoya Institute of Technology
Multi scale
Reynolds No.
10 12
台風(カトリーナ)
10 8
潜水艦
10 6
人
laminar
1
turbulence
ゾウリムシ
10 -4
バクテリア鞭毛
(2μm)
DNA
10 -9
scale
10 -6
10 -3
1
10 3
10 6
10 9
m
Nagoya Institute of Technology
レイノルズの相似法則を用いて流れの実験を行うことができる
模型 (幾何学形状、 長さ L )
流速 (U )
粘性率(n )
Re real =
UL
n
= Re model
Nagoya Institute of Technology
Wright 兄弟 の使った風洞
(スミソニアン博物館)
Nagoya Institute of Technology
NASA Ames Research Center 80 x 120 Wind Tunnel
Fans
http://www.windtunnels.arc.nasa.gov/
Nagoya Institute of Technology
F18
Nagoya Institute of Technology
F18 parts
MD11
Nagoya Institute of Technology
泳ぐ、飛ぶ生物と相似則
本川達雄 「ゾウの時間ネズミの時間」 中公新書
Nagoya Institute of Technology
Wikipedia
本川達雄 「ゾウの時間ネズミの時間」 中公新書
Nagoya Institute of Technology
本川達雄 「ゾウの時間ネズミの時間」 中公新書
Nagoya Institute of Technology
流れの安定性
初期条件
境界条件
境界上
Re が与えられたら 流れはただ一つに定まる?
y
例
ポアズイユ流れ
1
x
Nagoya Institute of Technology
かく乱の方程式
流れ=主流+かく乱
u’について2次の項は微小量として無視
Nagoya Institute of Technology
かく乱の方程式
u’が成長するかどうかで、流れの安定性を調べる
を代入
固有値問題に帰着
流れは不安定
流れは中立安定
流れは安定
Nagoya Institute of Technology
例
2次元ポアズイユ流れ
Orr-Sommerfeld の方程式
W. Heisenberg
k
中立曲線
Recr
Re
Copyright © 2012 Shin-ya Murakami [murashin _at_ gfd-dennou.org]
Nagoya Institute of Technology
Taylor-Couette 流れ
𝑅1 = 5
𝑅2 = 6
𝐻=8
𝑑 = 𝑅2 − 𝑅1 = 1
𝜂 = 𝑟𝑜 /𝑟𝑖 = 0.716
Ω1 𝑅1 𝑅2 − 𝑅1
𝑅𝑒 =
𝜈
Nagoya Institute of Technology
内円筒近傍での渦度の絶対値(Re=1000)
r-Z断面
q-Z断面
10pp
m
50ppm
Nagoya Institute of Technology
乱流数値解析講座
予告編
Provided by Nagoya Institute of Technology and Chubu CAE Forum
Nagoya Institute of Technology
乱流の特徴
1. 微小撹乱に対する不安定性
2. 強い非線形性
Nonlinear term
=
Viscous term
3. 巨大な輸送能力
Turbulent viscosity
=
Molecular viscosity
4. 無限大自由度
5. 散逸系
をもつランダムな流れ場
Nagoya Institute of Technology
Kolmogorov の乱流理論 (1941, K41)
Prandtl, Onsager, von Weizeker, Heisenberg
乱流における各スケールでの統計を考える
乱流理論の一里塚
Velocity increments
v
u
r
乱流の解析的統計理論
乱流の計算科学
Nagoya Institute of Technology
乱流数値解析講座 (アドバンストコース):乱流のk-ε モデル
乱流数値解析を行なうのに必要な基礎知識を講義
レイノルズ平均化された方程式のレイノルズ応力を渦動粘性係数でモデル化
渦動粘性係数を構成する乱流エネルギーkの輸送方程式の導出を解説
 k-εモデルの計算で解く散逸方程式の取扱に
ついて解説
 発達したチャネル内乱流熱伝達場をk-εモデ
ルで計算する実習付
 有限体積法による乱流熱伝達場解析を体験
k-εモデルで解析したチャネル内乱流熱伝達場
の平均速度分布(実習で算出を体験)
Nagoya Institute of Technology
12号館実験装置
自然対流(鉛直平板)
加熱に伴う浮力により引き起こされる対流熱伝達
強制対流(境界層風洞)
減速流れにおける熱伝達
高さ4m,幅1mの銅製加熱壁
強制対流(成層風洞)
都市環境(ビルや建物)における熱・物質移動