of-v1706中的IsoAdvector

[mohui]

@赵一铭 质量守恒有待推敲考证吧

[mengweilm425]

@mohui 这个isoAdvector我大概看过，每个有自由液面的网格(0<alpha<1)的等值面不一定一样。根据每个网格中心的alpha值先插值得到网格节点的alpha值，再利用节点的值构建等值面，使等值面以下的水的几何体积等于该网格中alpha*V的值。得到等值面的信息，包括面心位置，法相和速度，可以得到在dt时间内，等值面以下网格面浸水面积随时间的变化，再由已知的phi值计算dt时间内水的通量。也就是在计算这个通量的时候，考虑了等值面的运动。为了确保这个算法的守恒还用了bounding，大概就是如果某个网格接收水以后，其alpha值大于1，多余的部分不是直接砍掉，而是转移的别的网格。

[mohui]

@mengweil抱歉怎么晚了才看到，首先感谢大神的分享，仔细看了几遍你所说的内容，还是有两个疑惑望大神解答.
一：一个网格中心的值可以插成周边网格节点的值吗？（不知道你是否了解FTM，也是一种界面捕捉的方法，好像也是一个中心插值出很多个点的值）
二：“利用节点的值构建等值面，使等值面以下的水的几何体积等于该网格中alphaV的值”这句话是指，构建出来的等值面后，液体的几何体积一定会等于alphaV的值，还是以此条件调整等值面来保证水的几何体积等于该网格中alpha*V的值？:kiss:

[mengweilm425]

@mohui 不是大神，不是大神，互相交流:laughing: 。

1. FTM没了解过，这个isoAdevtor里面，我的理解是，比如二维的正交网格，每个网格节点都有四个网格相连，这个网格节点的值是由周围四个网格中心的值差值而来。
2. 是以alphaV为条件，找到对应的等值面。这个构建等值面的过程是根据每个网格节点处的alpha值，比如二维问题，我们就有四个alpha值，这样可以构首先建四个等值面。这四个等值面对应着四个alphaV，判断一下哪两个alphaV中间包含着我们的目标alphaV，然后就在这个范围内构建一个等值面以下的体积关于等值面的值的函数V=func(isovalue)，这样就可以找到目标等值面。这里面都是纯几何操作。

[mohui]

@mengweilm425 谢谢您的分享，我结合您所说的再消化消化。:big_mouth:

[yhdthu]

请问前辈对isoAdvector怎么看？～

[李东岳]

这个图里面没有包含界面捕获类，因此没有VOF等。

没细看isoAdvector不好评价，但是从他提供的算例结果来看简直吊炸天，基金会的MULES基本可以扔了。但是据我接触的算法都有优点和缺点，我感觉isoAdvector应该也有缺点，缺点如果不去细玩这个的话，外人一般看不出来。:confused:

[yhdthu]

@李东岳 前辈说的极是～这个isoAdvector吸引我的地方主要是两点，一个是任意网格下可保持界面sharp（理论上的），另外就是要比PLIC重构吃资源少些。这个方法目前只是用于无相变的solver，我是打算在这个基础上拓展一下，植入相变source的影响，我看了代码，有些地方我发现还是有trick存在的，虽然我也不是很懂目的是啥，但是觉得有搞头:lol:

[李东岳]

我觉得有意思，可以写一写文章。

你用VOF捕获相变？ZHC用Euler-Euler模拟高速空化？我也有兴趣用Euler-Euler+PBM模拟空泡。

[yhdthu]

@李东岳 前辈，关于空话初生的话，目前不是我的研究重点，关注点在空化发展和破碎，虽然也包括空泡的动力学特性，但是被模化为群泡效应了，不过高精度界面捕捉方法是都适用的，我打算先解决这个共性问题，再深入其他方向:lol:

[李东岳]

你要用多少万网格捕获气泡群的界面？

[yhdthu]

@李东岳 目前准备用LES进行模拟，网格数量应该不会太少
我不知道前辈所说的气泡群界面指的是需要网格解析度足够，要将每个单独的气泡都捕捉么？
目前主流的空化模型的类型为FCM(fully cavitation model)，这种模型只能将泡群的外轮廓捕捉，但是每个气泡的界面无法分辨出来。
我想这其中需要加入一些动量传递的源项之类的进行模化吧～:laughing:

[yuan_neu]

@东岳 ZHC用Euler-Euler模拟高速空化？
请问ZHC是哪位呀

[李东岳]

我也忘了..