结构网格上局部一维显式离散格式的并行自适应计算研究.pdfVIP

2008年伞围高性能计算学术年会 结构网格上局部一维显式离散格式的并行自适应计算 刘青凯 北京应用物理与计算数学研究所 高性能计算中心) accn) (1iuqk@iapcm 摘 要：局部一维显式离散格式是构造高维流体 力学计算格式的常用计算方法。，本文给出一个并行 p 2+尸 自适应显式时间积分算法，用于实现局部一维显式 pu pu 离散格式的并行自适应计算，、基于IASMIN框架． ∥ 0E Euo．“+ P 我们以局部一维FCT格式为例实现该算法，形成一 Meshkov 个用于模拟多介质流体力学Rlchtlnyer p为密度，R，v为速度，P为压力，E为总 不稳定性的并行自适应数值模拟程序。，数值实验袁 能量。从算子分裂的角度，基于均匀矩形结构网格 明，该程序具有很好的并行可扩展性。， 的显式离散格式可分为整体离散格式： 关键词：并行计算 自适应计算局部一维显式 离散格式TASMIN 氓。=饿，尝k。(卅F。，㈦一昙|G。(卅G，j(洲∽’ 和局部一维离散格式： 1 引 言 ㈦ WI，=wll，一去k¨，．∽’)一F，‘，【w“》 相列于串行计算，并行计算可充分利用当前的 高性能计算环境，扩大求解问题规模，提高求解精 w≯wil,，詈抚+!㈥G。÷w沩㈡ 度。，相对于单层网格计算，多层白适应网格计算采 从式(3)和(4)可知，各个坐标轴方向上的 用局部细化网格的方法，在物理敏感的局部计算区 时间积分之间具有数掂依赖关系(中问量Wv)， 域保持相同数值精度的情况下，可大量减少网格单 必须按坐标轴方向依次积分。 元，节省计算资源基于多层网格的并行自适应计 目前，基于均匀矩形结构网格的自适应计算方 算融合了这两种方法，已经成为高性能、大规模科 法已广泛应用于流体力学数值模拟。该方法在物理 学数值模拟的重要手段。 敏感的局部计算区域加密网格单元，形成多层均匀 住高维流体力学计算中，基于算子分裂的局部 矩形结构网格，如图1所示。其中，最粗网格层覆 一维显式离散格式通常用来构造高维计算格式，这 盖整个计算区域，绌网格层由相邻粗网格层的部分 种方法将高维问题分解为若干个一维问题，大大降 网格单元按一定比例细化得到，这个比例称为网格 低了理论分析和数值求解的难度，所以被广泛廊用 细化率、在计算过程中，细网格层动态变化，始终