of-v1706中的IsoAdvector

[mohui]

问下大家对于新出的0f-v1706中的IsoAdvector有什么更加通俗的理解吗？看了文章还是不太明白，虽然我对原来的mules也不是很懂，但是仅从对比的结果来看新的算法确实很牛逼，望大神普及点知识，以及建议。

[李东岳]

CFD界将要推一篇文章，正在写，敬请期待。

[mohui]

@李东岳 特别期待，先感谢一波东岳君

[mohui]

@李东岳 东岳老师，我大致看了作者文章的介绍，这里貌似用到了后处理软件中构建等值面的内容，不知道我理解对不对？另外这个所谓的几何VOF和代数VOF具体的区别在哪里呢，或者说几何VOF如何保证质量守恒呢？还有就是看了源码之后发现主程序中alpha方程特别简单，您觉得可以是否类似于interphasechangfoam那样添加sp，su？（问题有点多，先感谢一波:kiss: ）

[李东岳]

不过不能只通过代码判断算法的简单，顶层代码都非常简答。代码我还没看，暂时还不能回答你的问题。那个文章CFD界要推，在计划中，不过最近8月份课程马上要开始要备课，任务很紧，其他的事情要延后，嘿嘿 :confused:

[李东岳]

@李东岳 在 of-v1706中的IsoAdvector 中说：

CFD界将要推一篇文章，正在写，敬请期待。

应该是正在研究，写了个标题，看了看isoadvector模拟的算例原理。CFD界文章，尤其这种算法类都要断断续续写几星期。其他文章写的比较快。

[mohui]

@李东岳 :cheeky:

[赵一铭]

恩，基本上新闻类，叙事类，只算写文章的大约3-4小时。排版我俩一起大约需要1小时，现在排版稍微快了点，几十分钟就可以。
算法类多长时间就不一定了。算法类文章更长。其实CFD界的文章，看起来没多少，但是写起来可不容易。现在做原创的东西，比转载麻烦的不是100倍啊！！

关键是找素材。。。。找起来太费事了...

[赵一铭]

我又看了一下isoAdvector那个动画，简直太神器了，秒杀MULES啊

[赵一铭]

@mohui

Fluent理论指南

这里面有你说的几何重构，代数重构

[mohui]

@赵一铭 看了你所给的链接的资料，从图中可以明显看出来几何重构确实界面比代数重构要好的多。但是我还是有点没明白这样子重构出来的界面是否满足质量守恒呢？

[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模拟空泡。