随机扰动三维流体界面不稳定性的并行计算

对三维流体界面不稳定性的数值模拟引进了新的数值计算方法,并在 MPI 并行计算环境下进行了数值模拟。利用 Level Set 方法确定界面位置,零水平集对应界面位置。对应离散 Level Set 方程和界面两侧的两套 Euler 方程,借助于 Ghost 网格方法完成离散。如何确定临近界面的网格点的状态量,涉及到跨界面点的两套 Euler方程,借助于 Ghost 网格点思想实现在全流场用二阶 TVD 格式的统一计算。对最后网格点上的两套状态量的辨认依赖于该点的 LevelSet 值的符号。数值计算结果表明,TVD 格式结合 Ghost 方法较好地保持了跨界面密度的 Heaviside 剖面。采用区域分解方法对多处理机分配计算负载,每一个处理机有自己的内存和 CPU。讨论了缓冲区设定和数据交换方式。实际计算网格数为 102.4 万。得到Rayleigh-Taylor 不稳定性随机扰动的图象。

Parallel Computation of Fluid Interface Instability in Three Dimentions with Random Perturbation (In Chinese)

Computation of three dimensional fluid interface instability with random perturbation is performed on a 8 CPU parallel personal computer system. The parallel computation platform is message passage interface (MPI). A second order TVD scheme with Ghost method is applied with a fully parallel algorithm to the 3D Euler equations. Level Set function is used to track the motion of a fluid interface in an Eulerian framework. Buffer zone and data communication is discussed. The number of computation meshes is about 1 million. Bubber evolution with random perturbation in Rayleigh-Taylor instability is obtained by numerical simulation.