个人认为如果你在某一版本上开发了，就继续下去吧，不用移植了。。。。。

@李东岳 对 有

额 好像发现个问题，质量守恒方程应该是这个

其中𝑆(𝜌)与 𝑆(𝛾)形式一样，只不过自变量由相体积分数变成了密度，不知道这样对不对？

感觉有点弄乱了，求大佬指点，不胜感激

各位老师同学大家好，最近遇到了点小问题，真心求帮助。

假设流体不可压缩。MULES算法中，添加了人工压缩项以保证界面尖锐：

按我的理解，其实就相当于在VOF方程中引入了一个源项S

虽然界面压缩项只在自由表面附近起作用，但是严格来讲，连续性方程应该修改成如下形式才对吧？(因为对于不可压缩流体，相体积分数的控制方程与连续性方程是一致的)

但是在interfoam里面并未发现类似的代码，请问这是怎么回事呢？是对于这个问题，这个源项可以忽略吗？

推广一下，如果想在VOF方程中添加其他的源项，是不是也要修改质量守恒方程？

此外，对于相体积分数为0的网格，上式中除以相体积分数这个操作是不是需要特殊处理一下？我试了一下，假设这样

程序运行倒是没什么问题，但是这个意义我没搞清楚，这么一除不是得到无限大的速度散度了吗？目前，我倾向于这么做：对于相体积分数为0的网格，我直接将源项S的值设为0，即源项不起作用。

说的可能有点乱，但主要的意思就是这样。。。

谢谢大家！

@李东岳 在 interfoam控制方程的问题 中说：

@cfdngu http://dyfluid.com/interFoam.html 不是方程26么 这里面的推导跟nonNewtonianIcoFoam是一样的 那个更简单 你可以看看

好的，非常感谢东岳老师

各位老师好，下图是我推的动量方程， 但是我发现很多文献里面写的并不是这样。

我推的（忽略表面张力项）：

文献中常见的：

被我圈出来的这一项为什么不一样呢？是我链式法则用错了吗？

另外，有哪位老师能提供些intercompressiblefoam的学习资料？或者针对低速可压缩流体的教材的名称？之前一直做不可压缩流体，现在想做一些可压缩流体的工作，想尽快入门，谢谢！

@blookcfd
你好，感谢回复。

我发现我的问题可能是由于snappyHexMesh选定的区域与重叠区域重合导致的，现在已经解决了。

但是也有没解决的部分，我自己在原来overset网格的基础上加了一点东西，自己编译了另一个lib，这个lib与SHM还是有点冲突，运行的时候会报错，但实际上SHM执行了，网格也能用，所以我下一步好好看看我自己编的lib哪里有问题。

当用toposet + createPatch时，提示找不到stl文件，但是我的stl文件就在对应文件夹下放着，所以这个问题我还没找到原因。

我还在找寻替代方法，如果您知道其他方法能把stl放入计算域中，请不吝赐教，谢谢。

大家好，今天发现overset网格和snappyHexMesh好像有点冲突，导致用stl文件在计算域中创建边界不成功。请问有其他办法可以用stl文件(在计算域中)创建边界吗？我想把地形放到计算域里。

我尝试了topoSet+createMesh，但是使用topoSet的时候提示找不到stl文件，但我的stl文件没有变过，就在constant/triSurface下面，SHM能找到。

其他的方法也行，主要是要在overset 网格下中实现这个功能。

谢谢

@jlx 在 关于overset网格的问题 中说：

@cfdngu 您好，我最近遇到了一个差不多的问题，想请教一下。
在OFv1912中使用重叠网格模拟波浪和浮体的相互作用，在监测重叠网格区域内的波高时，发现监测的是背景网格不变的波高而不是重叠网格区域的波高（如下图），因此监测的波高值一直保持不变是0，这和您使用probe提取浮体边界某点流场信息遇到的问题相似，因此想请教您该怎样得到重叠区域的波高呢？谢谢！

你好，

很遗憾，我没解决这个问题，目前我认为需要自己调整sample的代码。我的问题通选择patchProbes暂时解决了，但是你那个可能有难度。

思路大概就两种：
1、将重叠区网格细化，因此重叠区域得以变小，让浪高仪在非重叠区测波面；
2、改代码，重点关注如何根据位置坐标找到对应网格编号的，过程中可能要用到cellZone或者cellType场来加以区分结构物网格与背景网格

@bestucan 好的，非常感谢，我再研究一下~

@bestucan 你好，非常感谢你的回复，我其实周日晚上就看到了，受到很大启发，也还有些不懂的地方，所以思考了两天才回复。

第一个问题，我认为应该如你所说，是显式和隐式的区别；
第二个问题，我将ALE框架下的积分形式的方程手动离散了一下，发现和MULES或isoAdvector中对应的代码一致，用的是相对通量求解的，而rhoPhi也是用相对通量组建的，所以不需要再次进行makeRelative操作；
第三个问题，我还是没太明白，按照我的理解，fvc::intepolate(U)就应该直接是新的面上速度了。我认为最后一句是不是类似非正交修正的操作？(尤其是相邻网格体心连线没有经过面心的时候)

大家好，最近重新研究interFoam，关于动网格的处理有些困惑，希望大家指点。

若网格进行运动，首先interFoam.C中要运行如下内容：

if (correctPhi)
{
        // 计算绝对通量
         phi = mesh.Sf() & Uf();

        // 求解pcorr，利用其修正phi，保证质量守恒
         #include "correctPhi.H"

        // 绝对通量转换为相对通量
         fvc::makeRelative(phi, U);

         mixture.correct();
}

这段代码所做的工作我在注释中标注出来了，不知道是否准确。我的第一个问题是，计算绝对通量为什么不是如下代码：

phi = fvc::interpolate(U) & mesh.Sf()

按我的理解，直接用Uf()访问的是上一迭代得出的面上速度，网格运动后其不一定与新网格面在一个平面了，而fvc::interpolate(U)是按照更新的网格插值得出的，这样貌似更合理些？

随后是求解相方程的内容，采用相对通量求解没有问题。

#include "alphaControls.H"
#include "alphaEqnSubCycle.H"

mixture.correct();

相方程求解后，我没有发现对rhoPhi采用类似makeRelative()的操作，即程序直接组建速度矩阵然后求解压力泊松方程，并在求解后才执行makeRelative()，位置在pEqn.H:70。这一点我没理解，为什么在动量方程的求解中未使用相对通量，而在求解后才显式地将phi转换为相对通量？

最后还有个小问题，在求解完压力方程后，将phi转化为相对值之前，即pEqn.H:67，程序还对面上速度进行了更新：

Uf() = fvc::interpolate(U);
surfaceVectorField n(mesh.Sf()/mesh.magSf());
Uf() += n*(phi/mesh.magSf() - (n & Uf()));

后两句存在的意义是什么呢？尤其是最后一句，实在没想明白。

劳烦大家指点下，我哪里理解有误还是哪个地方知识储备不够，这三个问题困扰我好几天了，查了不少资料也没找到原因，谢谢！