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

[CFDngu]

大家好，最近我在连续性方程中添加了质量源，但是求解时发现相体积分数最大值存在越界的情况，用的interFoam求解器。
我的质量源是正弦形式：

我修改了连续性方程，并打开了mules修正，然而得到的相体积分数最大值是这样的：

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

希望做过类似（添加过源项的）的童鞋老师指点一下：）

[李东岳]

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

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

[CFDngu]

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

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

[李东岳]

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

[CFDngu]

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

我写一下我的连续性方程和相方程吧，在源项区域内：

其中S为我添加的质量源。相方程用mules求解的时候就会有源项区域网格的相体积分数最大值大于1的情况。

[李东岳]

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

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

[CFDngu]

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

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

[mohui]

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

[CFDngu]

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

[mohui]

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

[CFDngu]

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

[Joann]

您好，我也在学习用OpenFOAM实现源造波。
请问质量源造波是只能通过在interFoam的pEqu.H中的压力方程的等号右边加入源项吗？
我了解到为了方便源项有fvOptions，就是在case里增加一个fvOptions文件，在这个文件里说明源项等相关设置，请问楼主在方面是否有尝试？
我目前在做源项造jonswap不规则波，但是在引入的源项方面有些不确定，原理是Lin和Liu的1999年的文献中的这一段：

所以想了解楼主加入源项的源项是怎么表达的？加入源项后是通过在求解器里wmake编译就可以在case里使用源造波功能了吗？
希望能得到老师和各位前辈们的指导，感激~