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

@李东岳 对 有

@bestucan 在 请问如何理解这一项？ 中说：

我猜是任意方向的加速度。

速度本身是矢量，做梯度就是个二阶张量了，前后乘以n和n的转置，应该是类似于对角化操作（基于某方向的对角化）。把各方向的加速度解耦出来，使之不在其他方向有分量。

这个“任意方向”就是右边那项分母上n的方向，咦，n就是法向哈

那应该就是了

但是要注意右边那项分子，速度有下标的，应该是n方向上的速度。
那这一项的意义应该是：“n方向上的速度，在n方向的加速度”

这也和左边对角化操作印证，嗯，就这么猜

理解了，这项应该是一个标量，即“沿自由表面法向的加速度”。

我水平有限，还是需要多研究。

非常感谢您的不吝赐教！

@bestucan 在 请问如何理解这一项？ 中说：

这，速度对坐标的偏微分，不就是加速度么。

三维上的，所以可以用左边的速度梯度（最大加速度方向的加速度）和坐标换算各个方向的加速度。

大概这个意思

感谢老师指点，我可以这么理解吗：这一项是沿自由表面法向的加速度？

各位老师好，最近我在文献中看到这一项：
\begin{equation}
\textbf{n}^T \cdot \nabla\textbf{u} \ \cdot \textbf{n} = \frac{\partial u_n}{\partial n}
\end{equation}
其中$\textbf{n}$是自由表面法向量，$\textbf{u}$为速度，但是右边那一项我看的不太明白，请问代表什么？文献中没有说明，好像大家都知道就我不知道

出处：
Haghshenas M, Wilson J A, Kumar R. Finite volume Ghost Fluid Method implementation of interfacial forces in PISO loop. Journal of Computational Physics, 2019, 376:20-27.

求大神解答，谢谢！

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

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

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

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

假设流体不可压缩。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场来加以区分结构物网格与背景网格