openfoam温度场计算



  • 各位大神,我在做一个圆管的数值模拟一个入口,一个出口。用buoyantBoussinesqSimpleFoam这个求解器模拟温度场,现在遇到了瓶颈,关于各个参数的边界条件设置。请问我的设置是否正确???多谢各位
    ![alt text](

    • image url)![alt text](0_1490665908800_QQ图片20170328094458.png image url)0_1490665918438_QQ图片20170328094647.png 0_1490665933684_QQ图片20170328094739.png 0_1490665944643_QQ图片20170328094744.png 0_1490665953688_QQ图片20170328094748.png 0_1490665959312_QQ图片20170328094752.png 0_1490665966793_QQ图片20170328094755.png

  • 网格教授 OpenFOAM教授 管理员

    你好,估计问题已经解决了,你的设置看起来没有问题。



  • @李东岳 谢谢您!我模拟的是管道流动,您觉得我的速度壁面边界条件设置成zeroGradient可以吗???



  • @zhanghan 不行,管道壁面上速度应该设置为无滑移边界条件,在OpenFOAM 中就是

    type       fixedValue;
    value     uniform (0 0 0);
    


  • @xpqiu 感谢您的回复,我在fvsolution中添加参考压力设置,如下
    SIMPLE

    nNonOrthogonalCorrectors 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.
    我实在没弄明白是哪里出错。我参考压力设置不正确吗??还是哪里出了问题呢!望指点!!多谢多谢!!



  • @zhanghan 出口压力不要用 zeroGradient,用

    type      fixedValue;
    value    uniform 0; // 或者别的值,根据你的需要来设置
    


  • @xpqiu

    多谢您的回答!!!
    我添加了重力场,这样设置,出口流场和压力场都出现了问题(回流)



  • @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§),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
    

  • OpenFOAM讲师

    @haining-luo 参考压力应该是只有在没有进出口边界条件的情况下才发挥作用,例如封闭罐体。所以你的算例里设置pRef实际上并没有发挥作用。