of-v1706中的IsoAdvector



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



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



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



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



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



  • @李东岳of-v1706中的IsoAdvector 中说:

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

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




  • 管理员

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

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


  • 管理员

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


  • 管理员

    @mohui

    Fluent理论指南

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



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



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



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



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



  • @mohui 不是大神,不是大神,互相交流:laughing:

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


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


登录后回复
 

与 CFD中文网 的连接断开,我们正在尝试重连,请耐心等待