文章链接:https://doi.org/10.1016/j.cpc.2021.107849
算例部分在第四章。
4.1基础测试
4.1.1泊肃叶流动
设置一个8*8*2网格,模拟二维正方形水槽内的泊肃叶流动。参数[latex]||\nabla p||=\mu=\rho=1[/latex].在水槽左端右端设置压强边界条件。模拟至流场以1e-3精度稳定,约在1000步后,时刻1达成。使用显式viscous diffusion.上下边界处使用粘性边界条件u=0.
4.1.2二维斯托克斯流
和泊肃叶流动相同,但计算区域中心有一直径为0.5的圆盘。不进行advection以保证稳定。使用显式粘性力。如果使用隐式粘性,则收敛速度较快,但有一个量级为[latex]\tau[/latex]的误差项。除此之外还有两个测试。第一个是圆盘在重力作用下下落,第二个是设置流场入口/出口的速度边界条件而非用压差驱动。
4.1.3水滴对流
似乎是设置一个水滴让它在计算区域内移动,以证明VOF方法是正确的,能保证水滴的形状不变。
4.1.4圆柱对流
暂缺
4.1.5计算速度
大致每秒1.6e6个格子(单处理器)。
4.2表面张力波
二维计算域,两种密度不同的流体,在中间生成一个表面张力波。
4.3水滴和泡泡震荡
水滴是内部水外部空气,泡泡是外部水内部空气。一开始初始化为椭球,让它震荡变成球。用两个不同的分辨率,对比误差。
4.4雨滴下落
暂缺。特地选择了一个比较低的分辨率,因为低分辨率容易崩。