PARIS模拟器算例介绍

文章链接:https://doi.org/10.1016/j.cpc.2021.107849

算例部分在第四章。

4.1基础测试

4.1.1泊肃叶流动

设置一个8*8*2网格,模拟二维正方形水槽内的泊肃叶流动。参数\(||\nabla p||=\mu=\rho=1\).在水槽左端右端设置压强边界条件。模拟至流场以1e-3精度稳定,约在1000步后,时刻1达成。使用显式viscous diffusion.上下边界处使用粘性边界条件u=0.

4.1.2二维斯托克斯流

和泊肃叶流动相同,但计算区域中心有一直径为0.5的圆盘。不进行advection以保证稳定。使用显式粘性力。如果使用隐式粘性,则收敛速度较快,但有一个量级为\(\tau\)的误差项。除此之外还有两个测试。第一个是圆盘在重力作用下下落,第二个是设置流场入口/出口的速度边界条件而非用压差驱动。

4.1.3水滴对流

似乎是设置一个水滴让它在计算区域内移动,以证明VOF方法是正确的,能保证水滴的形状不变。

4.1.4圆柱对流

暂缺

4.1.5计算速度

大致每秒1.6e6个格子(单处理器)。

4.2表面张力波

二维计算域,两种密度不同的流体,在中间生成一个表面张力波。

4.3水滴和泡泡震荡

水滴是内部水外部空气,泡泡是外部水内部空气。一开始初始化为椭球,让它震荡变成球。用两个不同的分辨率,对比误差。

4.4雨滴下落

暂缺。特地选择了一个比较低的分辨率,因为低分辨率容易崩。

发表评论

您的电子邮箱地址不会被公开。 必填项已用*标注