并行计算方法.pptVIP

2006-12-10 抛物型方程并行计算方法袁光伟 姚彦忠北京应用物理与计算数学研究所 内容简介 1. 问题背景 2. 研究历史 3. 守恒型并行格式 4. 数值实验 5. 小结 1. 问题背景 天体物理，(等离子体物理……) 大气与海洋流场模拟 (工业，ICF，高新技术，……) 地质、油藏模拟…… 1. 问题背景 物质热传导 单温热传导 1. 问题背景 1. 问题背景：精细建模 三温模型 物质温度还可以细分为离子温度和电子温度，这两者在极端情况下也会分离； (多群扩散模型，多群输运模型) 辐射能量可以区分为与辐射频率相关的分群能量密度。 磁流体模型 磁约束聚变中还需要考虑磁场的影响。 1. 物理模型 1. 物理模型 1. 物理模型 1.非线性抛物型方程组 辐射、电子或离子热传导方程、以及磁扩散方程都具有下列形式 重要性： 直接影响系统能量的传输、分配和交换； 确保高置信度的关键之一。 1. 计算方法 2. 研究历史 无条件稳定 数据强相关，不能直接并行求解 大型稀疏线性代数方程组，迭代求解，计算量很大，全局通信 2. 研究历史 易于并行 条件稳定， 网格比限制 不实用 多介质，材料的特性相差大，需加密网格 模拟时间长，计算效率不高 2. 研究历史 并行格式设计目标 (1) (稳定性) 保持与隐格式同样稳定； (2) (高精度) 保持与隐格式同样精度； (3) (高并行度) 自由分割求解区域，局部通信； (4) (易于实施) 适用于现有隐式程序并行化，便于求解辐射流体力学耦合问题； (5) (经济性) 子区域上的问题的规模变小，条件数变小，迭代收敛加速； (6) (守恒性) 越过分区界面保持守恒。 两层并行格式 两层并行格式,界面三层 2. 研究历史 2. 研究历史 A. 交替分块显隐三层并行格式 无条件稳定，设计复杂 无预估校正 B. 界面显式处理的两层并行格式 条件稳定, 或：无条件稳定、一阶精度 (无)预估校正 2. 研究历史 C. 基于界面守恒修正的迭代并行格式 无条件稳定，二阶精度，易于实施 实现途径： 不同内边界条件耦合； 预估(守恒)校正技术与迭代方法耦合 3. 研究历史 3. 守恒型并行格式 3. 守恒型并行格式 全隐Picard迭代方法 其中 s 为非线性迭代次数 3. 守恒型并行格式 首先，内界面处采用Dirichlet边界条件 其次，在各个分区子区域独立进行求解 3. 守恒型并行格式 然后，进行非线性迭代直至收敛； 最后，在分区界面采用Neumann边界条件，再做一次非线性迭代 扩散格式 4.数值实验 算例 1 算例 2 三维可视化图 5. 小结 三温能量方程组，磁扩散方程 构造了守恒型并行格式 易于实施、保持隐格式精度、无条件稳定、局部通讯、高并行效率 实际应用于Z-箍缩问题求解。 5. 小结 “瓶颈”问题： 不适应大规模并行计算的数值方法； 全局通讯 稳定性，精度 程序实现 5. 小结 并行格式设计的基本原则 分而治之 界面协调 整体优化 分而治之 把一个问题分裂成若干个小规模的子问题，然后独立地并行求解这些子问题。 “之”是指所分解的子问题 “治”主要是指如何求解这些子问题。 高并行度 步长稳定性限制、收敛阶降低、 计算复杂性 界面协调 核心：在分区界面处，如何进行预估修正 保持较高并行度 改进界面值的精度，整体精度提高 适当恢复数据强相关、全局依赖和全局影响的本性 保持无条件稳定性。 九点格式 (李德元等) 迭代法 守恒型并行格式 内边界条件类型 （1） Dirichlet 边界条件 无条件稳定，不守恒 （2） Neumann边界条件（流条件） 条件稳定，精度高，守恒 Dirichlet 边界条件 + Neumann边界条件 无条件稳定 + 守恒 ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ?