interFoam.C

1、 interFoam 求解流程 初始化流场 修正流率 解相连续性方程 是否进行 子循环 子循环是 否结束 解相连续性方程 Y N 求解动量方程 Y N 求解压力方程 修正流率和速度 PISO循环是 否结束 End Y N 质量流率修正 2、 interFoam 控制方程组 基本控制方程组 2.1 相连续性方程： ∂∂αtg + ∇⋅ + ∇⋅ − = ( ) ( ) U U 1 α α α g g g ⎡ ⎣ r ⎤ ⎦ 0 （1） U n min C ,max r α S S ⎡ ⎤ φ φ ⎛ ⎞ = ⎢ ⎜ ⎟ f ⎢ ⎥ ⎣ ⎦ f f ⎜ ⎟ ⎝ ⎠⎥ （2） ( ) ( ) α α δ ∇ = ∇ + f f f n n 1/3 ε δ = ⎛ ⎞ ⎜ ⎟ ⎜ ⎟ ⎜ ⎟ ⎝ ⎠ n ∑ i N V N φ = ⋅ U S f ρ = + − α ρ α ρ g g g (1 ) l , Fm = ρφ ρ φ = ⋅ ， Cα 压缩因子，在计算气泡上