添加源项后相方程求解越界,请问这种情况有人遇到过吗?



  • 大家好,最近我在连续性方程中添加了质量源,但是求解时发现相体积分数最大值存在越界的情况,用的interFoam求解器。
    我的质量源是正弦形式:
    0_1500434513347_1.jpg
    我修改了连续性方程,并打开了mules修正,然而得到的相体积分数最大值是这样的:
    0_1500434569171_2.jpg

    源项最大值对应相体积分数最大值,反之也是,计算没有报错,然而就是有着5%的平均误差,这是咋回事?

    希望做过类似(添加过源项的)的童鞋老师指点一下:)



  • 源项最大值对应相体积分数最大值,反之也是,计算没有报错,然而就是有着5%的平均误差,

    第二个是你的相分数?你的相分数一直都是大于1?不是很理解你的只有5%的平均误差的意思。相分数应该在0>1之间。



  • @李东岳 东岳老师你好,我这个是相体积分数的最大值。

    目前计算时,质量源区域附近的相体积分数最大值超过1.



  • 检查你的源项设置,从数学的角度,如果不考虑对流,相方程为:
    \begin{equation}
    \frac{\partial \rho\alpha}{\partial t}=m
    \end{equation}
    最大的m会使得alpha变为1,如果你的m过大,这不符合物理,并且导致alpha要超过1了。



  • @李东岳 不好意思东岳老师你的公式乱码了,我没有看懂您的回复。

    我写一下我的连续性方程和相方程吧,在源项区域内:
    0_1500598315908_QQ截图20170721085114.jpg
    其中S为我添加的质量源。相方程用mules求解的时候就会有源项区域网格的相体积分数最大值大于1的情况。



  • 公式乱码目前由于网页技术问题没有解决,可以刷新一下网页试试。

    我确认了下你的方程这两个是一致的。你的压力方程做了对应的改动么?



  • @李东岳
    我改了,interFoam中求解器中,连续性方程的修改不就是通过改压力方程实现的嘛。。。。。

    我在这句加的源项:
    fvm::laplacian(rAU, p) == fvc::div(phiHbyA)



  • @CFDngu 你的这个问题和空化是一个道理,话说,你的第二相是通过1-alpha得到吗?



  • @mohui 你好,我的第二相是通过1-alpha1得到的,我没做过空化不太了解,难道问题在这里?



  • @CFDngu 这个我也不是很清楚,可能算的结果就很难保证有界吧,有时候算空化也是会出现越界的现象。这个看你所能接受的精度吧。



  • @mohui
    目前来讲精度不怎么好。。因为相体积分数大于1的网格对其他的网格也会有影响,这样整体的求解多少(比如速度场)会有些问题。


登录后回复
 

与 CFD 中国 的连接断开,我们正在尝试重连,请耐心等待