Your browser does not seem to support JavaScript. As a result, your viewing experience will be diminished, and you have been placed in read-only mode.
Please download a browser that supports JavaScript, or enable it if it's disabled (i.e. NoScript).
大家好,最近我在连续性方程中添加了质量源,但是求解时发现相体积分数最大值存在越界的情况,用的interFoam求解器。 我的质量源是正弦形式: 我修改了连续性方程,并打开了mules修正,然而得到的相体积分数最大值是这样的:
源项最大值对应相体积分数最大值,反之也是,计算没有报错,然而就是有着5%的平均误差,这是咋回事?
希望做过类似(添加过源项的)的童鞋老师指点一下:)
源项最大值对应相体积分数最大值,反之也是,计算没有报错,然而就是有着5%的平均误差,
第二个是你的相分数?你的相分数一直都是大于1?不是很理解你的只有5%的平均误差的意思。相分数应该在0>1之间。
@李东岳 东岳老师你好,我这个是相体积分数的最大值。
目前计算时,质量源区域附近的相体积分数最大值超过1.
检查你的源项设置,从数学的角度,如果不考虑对流,相方程为: \begin{equation} \frac{\partial \rho\alpha}{\partial t}=m \end{equation} 最大的m会使得alpha变为1,如果你的m过大,这不符合物理,并且导致alpha要超过1了。
@李东岳 不好意思东岳老师你的公式乱码了,我没有看懂您的回复。
我写一下我的连续性方程和相方程吧,在源项区域内: 其中S为我添加的质量源。相方程用mules求解的时候就会有源项区域网格的相体积分数最大值大于1的情况。
公式乱码目前由于网页技术问题没有解决,可以刷新一下网页试试。
我确认了下你的方程这两个是一致的。你的压力方程做了对应的改动么?
@李东岳 我改了,interFoam中求解器中,连续性方程的修改不就是通过改压力方程实现的嘛。。。。。
我在这句加的源项: fvm::laplacian(rAU, p) == fvc::div(phiHbyA)
@CFDngu 你的这个问题和空化是一个道理,话说,你的第二相是通过1-alpha得到吗?
@mohui 你好,我的第二相是通过1-alpha1得到的,我没做过空化不太了解,难道问题在这里?
@CFDngu 这个我也不是很清楚,可能算的结果就很难保证有界吧,有时候算空化也是会出现越界的现象。这个看你所能接受的精度吧。
@mohui 目前来讲精度不怎么好。。因为相体积分数大于1的网格对其他的网格也会有影响,这样整体的求解多少(比如速度场)会有些问题。
您好,我也在学习用OpenFOAM实现源造波。 请问质量源造波是只能通过在interFoam的pEqu.H中的压力方程的等号右边加入源项吗? 我了解到为了方便源项有fvOptions,就是在case里增加一个fvOptions文件,在这个文件里说明源项等相关设置,请问楼主在方面是否有尝试? 我目前在做源项造jonswap不规则波,但是在引入的源项方面有些不确定,原理是Lin和Liu的1999年的文献中的这一段: 所以想了解楼主加入源项的源项是怎么表达的?加入源项后是通过在求解器里wmake编译就可以在case里使用源造波功能了吗? 希望能得到老师和各位前辈们的指导,感激~