openfoam温度场计算
-
@xpqiu 感谢您的回复,我在fvsolution中添加参考压力设置,如下
SIMPLEnNonOrthogonalCorrectors 0; consistent yes; pRefPoint (2 0.0025 0.0025); pRefValue 0;
并将出口压力设为zeroGradient,系统提示如下,
continuity error cannot be removed by adjusting the outflow.
please check the velocity boundary condition and/or potentialFoam to initial list the outflow.
我实在没弄明白是哪里出错。我参考压力设置不正确吗??还是哪里出了问题呢!望指点!!多谢多谢!! -
@xpqiu
cavity 里面是因为没有进出口所以才可以设置成 p zeroGradient ? cavity里面对p的BC都是Neuman所以必须设定pRef。(U全部是fixedValue)我的模拟里面有入口和出口,出口可能离上游复杂流动不够远但觉得 p fixedValue uniform有点太强,换成zeroGradient,(所以p的BC全部变成了Neuman) 即使我设定了pRef =>
continuity error cannot be removed by adjusting the outflow.
please check the velocity boundary condition and/or potentialFoam to initial list the outflow.
数学上有什么解释吗?有没有其他的边界条件可选?跟cavity里面不同的地方在我的算例里面有进口1,2,U都设定fixedValue类的BC,出口U zeroGradient。但求解U的过程没有问题(求解U用的是fvc::grad(p)),icoFoam走到压力方程才报错。按道理p方程都给了Neumann,给了pRef,为啥还是adjuctPhi error?
参考贴:http://www.cfd-china.com/topic/305/adjustphi的作用是检查边界条件
这里是第一个时间步:
Time = 1e-05 Courant Number mean: 1.07669e-05 max: 0.00554049 U_tau ajust for bulk velocity 0.3 max of yPlus = 84.6308 min of yPlus = 0.733485 U_tau ajust for bulk velocity 0.3 max of yPlus = 84.6308 min of yPlus = 0.733485 DILUPBiCCCG: Solving for Ux, Initial residual = 1, Final residual = 2.61292e-07, No Iterations 1 DILUPBiCCCG: Solving for Uy, Initial residual = 1, Final residual = 2.73053e-07, No Iterations 1 DILUPBiCCCG: Solving for Uz, Initial residual = 1, Final residual = 6.39136e-07, No Iterations 1 --> FOAM FATAL ERROR: Continuity error cannot be removed by adjusting the outflow. Please check the velocity boundary conditions and/or run potentialFoam to initialise the outflow. Total flux : 1.58773e-07 Specified mass inflow : 2.54763e-05 Specified mass outflow : 0 Adjustable mass outflow : 0
-
@haining-luo 参考压力应该是只有在没有进出口边界条件的情况下才发挥作用,例如封闭罐体。所以你的算例里设置pRef实际上并没有发挥作用。