再次求教不可压流体计算时,边界处压力震荡的问题



  • 上次提了问题,就是发现利用不可压二相动网格求解器计算波浪冲击直墙时,压力出现震荡,地址如下:
    http://www.cfd-china.com/topic/784/求教关于计算开始时压力震荡的问题/5

    多谢@李东岳 和 @wwzhao 两位版主的回复,我回去又研究了一下,发现以前有不少算例都出现震荡情况,只不过后来的压力变化相对于之前的震荡值量级比较大,因此掩盖了这个问题,但是问题依旧存在,所以我加了一个直墙,设置的边界条件与右边边界相同:p_rgh中为 fixedFluxPressure ,U中为fixedValue,大小都为0.

    再交代一下,我采用动网格造波,即通过边界的运动模拟实现类似于推板造波机的功能。之前造波效果都不错,直到最近发现计算压力存在震荡现象,这次我又算了几个例子,首先我的计算区域是这样的:
    算例长21m,设置一个直墙,位于x=20m。右端设置为直墙,左端动边界,距离动边界越远,网格运动的幅度越小。
    0_1492567693235_0s-chushi.jpg
    初始速度场:
    0_1492567614566_u-0s.jpg
    计算开始后,我发现速度场在之直墙两端有突变:
    t=0.5s,整体速度场:
    0_1492567862778_u-0.5s.jpg
    直墙附近压力震荡:
    0_1492567891321_u-0.5s-jubu.jpg
    t=1s,整体速度场:
    0_1492567920059_u-1s.jpg
    直墙附近压力震荡:
    0_1492567927302_u-1s-jubu.jpg

    基本上计算开始几个时间步后,就会在直墙的水、气交界面开始出现这种问题,最后就会在得到震荡的压力。
    我还特意试了一下wave2foam速度入口造波,发现没有这个问题。虽然我知道不可压缩流体的速度和压力瞬间更新全场,但是这个地方怎么这么诡异?还是动网格导致的呢?如何能够优化或者避免这种情况?

    请大牛们不吝赐教,多谢了:)