管道内可压缩空气流动出口压力边界条件设置问题 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】

边界条件:

ParameterInletOutletWall
VelocityzeroGradientzeroGradient0
Pressure600000Pa101325PazeroGradient
Temperature300KzeroGradientzeroGradient

输运参数: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:

 

  • 0
    点赞
  • 5
    收藏
    觉得还不错? 一键收藏
  • 1
    评论
OpenFOAM 是一个开源的计算流体力学(CFD)软件包,它提供了广泛的功能和工具来模拟流体流动和传热问题。在 OpenFOAM 中,二次边界条件是指在边界上定义的非线性边界条件。 要定义一个二次边界条件,你需要在边界条件文件中指定一个函数,该函数描述了该边界上的非线性关系。下面是一个示例,展示了如何在 OpenFOAM 中定义一个二次边界条件: ```cpp // 在边界条件文件中定义二次边界条件 #include "fvCFD.H" // 定义边界条件 void fixedValueParabolicVelocityFvPatchVectorField::updateCoeffs() { const volScalarField& alpha = mesh().alphaEff(); const surfaceScalarField& phi = boundaryField()[patchI()].phi(); const scalarField& U = this->patchInternalField(); const vectorField& gradU = fvc::grad(U); const scalarField A = 1.0 - alpha; const vectorField gradA = fvc::grad(A); const scalarField C = A*phi - 0.5*alpha*magSqr(gradU) - dot(gradA, gradU); const vectorField D = U - dot(U, gradA)*gradA; const surfaceScalarField& alphaSf = alpha.boundaryField()[patchI()]; forAll(faceFlux_, faceI) { faceFlux_[faceI] = alphaSf[faceI]*C[faceI] + dot(D[faceI], Sf()[faceI]); } } ``` 上述示例中的函数 `fixedValueParabolicVelocityFvPatchVectorField::updateCoeffs()` 定义了一个二次边界条件,其中 `alpha` 是边界上的非线性系数,`phi` 是边界上的标量场,`U` 是速度场。这个函数根据定义的关系计算边界上的二次边界条件。 需要注意的是,定义二次边界条件需要对 OpenFOAM 的源代码进行修改,并重新编译。因此,在使用二次边界条件之前,请确保你已经熟悉了 OpenFOAM 的基本操作和编译流程。 希望这个回答能帮到你!如有任何问题,请随时提问。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值