管道內可壓縮空氣流動出口壓力邊界條件設置問題 in openfoam

       在這個案例中我們來討論管道內亞音速/超音速出口壓力邊界條件:1)fixedValue  uniform 101325。2)waveTransimissive。對仿真結果的影響。

1、問題描述

計算域:

【直徑:0.03mm】【長度:40mm】【長徑比:1】

2、控制方程

質量守恆:\frac{\partial \rho }{\partial t}+\bigtriangledown \cdot (\rho {\color{Green} U})=0

動量守恆:\frac{\partial \rho {\color{Green} U}}{\partial t}+\bigtriangledown \cdot (\rho {\color{Green} UU})-\bigtriangledown \cdot \tau =-\bigtriangledown p

能量守恆:\frac{\partial \rho h}{\partial t}+\bigtriangledown \cdot (\rho {\color{Green} U}h)+\frac{\partial \rho K}{\partial t}+\bigtriangledown \cdot (\rho {\color{Green} U}K)-\bigtriangledown \cdot (\alpha_{eff}\bigtriangledown h)=\frac{\partial p}{\partial t}

狀態方程:p=\rho RT

初始條件:【U=0m/s】【p=101325pa】【T=300K】

邊界條件:

Parameter Inlet Outlet Wall
Velocity zeroGradient zeroGradient 0
Pressure 600000Pa 101325Pa zeroGradient
Temperature 300K zeroGradient zeroGradient

輸運參數:dynamic viscosity of air u=1.79e-05kg/m-s

熱力學參數:【摩爾質量molweight:28.9】【Cp:1005】【Prandtl:1】【Rspecific:240.29j/mol*k】【conductivity:k=32.3uw/mk】

設計聲速:【 c=\sqrt{\gamma RT}=307m/s】【C_{p}-C_{v}=R】【\gamma =\frac{C_{p}}{C_{v}}=\frac{C_{p}}{C_{p}-R}=1.31

3、網格生成

網格採用ICEM來劃分,此處只截取入口處的網格示意。

4、openfoam運算

fixedValue 101325 boundary condition:

waveTransimissive boundary condition:

 

發表評論
所有評論
還沒有人評論,想成為第一個評論的人麼? 請在上方評論欄輸入並且點擊發布.
相關文章