请问这是OF中的假扩散现象吗?



  • 之前一直听说OF中的VOF方法存在假扩散现象,后来还有了IsoAdvector方法。最近做计算的时候遇到了类似的现象,因此来和大家讨论一下:

    我要模拟的是一个波在斜坡上传播的过程,水面和斜坡后方一样高,采取自己写的动边界造波,大概是这样:
    0_1504686684499_all09.jpg

    可是在造波的过程中遇到了类似于假扩散的情况,初始时刻:
    0_1504686724073_problems0.jpg
    造波刚刚开始,直接在交界面处出现了一个速度:
    1_1504686724073_problems1.jpg
    后面明显看到,波浪生成并传播的过程中,坡顶存在一个沿着波浪传播方向的速度,并且在传播,但是此时自由面并未移动,可以观察到一些相体积分数的“扩散”现象:
    2_1504686724074_problems2.jpg
    3_1504686724074_problems3.jpg
    4_1504686724074_problems4.jpg
    5_1504686724074_problems5.jpg
    6_1504686724074_problems6.jpg
    最后导致坡顶附近的速度测点得到如下结果:
    0_1504687272814_QQ截图20170906164043.jpg
    而实际上,无论是根据常识还是模型试验,那个扰动是不应该存在的,坡顶初始干潜礁的速度应该是从0开始的。

    因此特来问问各位大牛,这个现象是怎么回事,是否就是OF中VOF方法的问题,计算中是否有人遇到过,该如何解决呢?


  • 网格教授 OpenFOAM教授 管理员

    假扩散是类似一种扩散的行为,你这个看起来不是假扩散,边界造波不应该影响远处的内部场,你的代码除了边界改动之后,内部场也改动了?你的内部网格动了么?



  • @李东岳 东岳老师你好,我的内部改动了,但是理论上应该不会引起这个问题。

    主要是可以观察到相体积分数的扩散(波还没传过来)。可能和我这个地方的网格有关系?我再去看看