Loading [MathJax]/jax/output/CommonHTML/jax.js

2022年6月11日土曜日

地下水~地表流まで

・層流 : i=au
・層流から乱流への遷移状態、乱流 : i=au+bu2
・乱流 : i=bu2

1. 層流(粒径小、速度低)

Darcy:h=1ku=au h=p/(ρg)よりp=ρg1ku=ρgau=αu
Darcy-Kozeny Carman:
a=1k=180ν(1ne)2gΦ2n3ed2 ν=μ/ρよりα=ρgk=180μ(1ne)2Φ2n3ed2


2. 層流~乱流への遷移状態、乱流(粒径大、高速)

Forchheimer:
h=au+b|u|u
Ergun(1952)
a=1k=150v(1ne)2gn3ed2 b=1.75(1ne)gn3ed

Kadlec and Knight(1996):角ばった粒子に適用。
a=1k=255v(1ne)gn3.7ed2 b=2(1ne)gn3ed


3.乱流

Manning:
u=1nR23I12I=n2R43u2=bu2
h:水頭(m)
p:ポテンシャル(Pa)
ρ:密度(kg/m3)
k:透水係数(m/s)
u:速度(m/s)
α:係数(s/m)
β:係数(s2/m2)
Φ:球形度(0~1、球=1)
d:有効粒径(m)
μ:粘性係数(Pa・s)
ν:動粘性係数(m2/s)=μ/ρ
ne:間隙率
g:重力加速度(m/s2)
n:マニングの粗度係数(m−1/3・s)
R:径深(m)
I:勾配


OpenFOAM
https://openfoamwiki.net/index.php/DarcyForchheimer
https://phreeqc.blogspot.com/2022/02/darcy-forchheimer-model.html?m=0
h=au+b|u|up=ρgau+ρgb|u|u=αu+β|u|u=μDu+12ρF|u|u
Darcy:
α=ρgk=μDD=ρgkμF0
Darcy-Kozeny Carman:
α=ρgk=180μ(1ne)2Φ2n3ed2=μDD=180(1ne)2Φ2n3ed2F0
Forchheimer-Ergun:
α=ρga=ρgk=150μ(1ne)2n3ed2=μDD=150(1ne)2n3ed2β=ρgb=1.75ρ(1ne)n3ed=12ρFF=3.5(1ne)n3ed

PersianSPHでは、式と係数を選択。ソース内で係数変更可。
OpenFOAMでは、D,F 入力で表現。

0 件のコメント:

コメントを投稿