・層流 : 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ν(1−ne)2gΦ2n3ed2 ν=μ/ρよりα=ρgk=180μ(1−ne)2Φ2n3ed2
2. 層流~乱流への遷移状態、乱流(粒径大、高速)
Forchheimer:
−∇h=au+b|u|u
Ergun(1952)
a=1k=150v(1−ne)2gn3ed2 b=1.75(1−ne)gn3ed
Kadlec and Knight(1996):角ばった粒子に適用。
a=1k=255v(1−ne)gn3.7ed2 b=2(1−ne)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|u−∇p=ρgau+ρgb|u|u=αu+β|u|u=μDu+12ρF|u|u
Darcy:
α=ρgk=μDD=ρgkμF=0
Darcy-Kozeny Carman:
α=ρgk=180μ(1−ne)2Φ2n3ed2=μDD=180(1−ne)2Φ2n3ed2F=0
Forchheimer-Ergun:
α=ρga=ρgk=150μ(1−ne)2n3ed2=μDD=150(1−ne)2n3ed2β=ρgb=1.75ρ(1−ne)n3ed=12ρFF=3.5(1−ne)n3ed
PersianSPHでは、式と係数を選択。ソース内で係数変更可。
OpenFOAMでは、D,F 入力で表現。
0 件のコメント:
コメントを投稿