P000115
基于JASMIN框架的辐射流体与粒子输运耦合计算
*任健 (北京应用物理与计算数学研究所)
魏军侠 (北京应用物理与计算数学研究所)
曹小林 (北京应用物理与计算数学研究所)
在应用研究领域中,许多复杂物理问题通常由多个不同种类的物理过程耦合构成。数值模拟此类复杂物理问题时,通常将其分解为多种不同类型的物理过程,然后分别针对单个物理过程提出数学模型方程,并且设计相应的数值计算方法与并行程序,最后将这些独立开发的应用程序耦合起来共同完成。例如,在惯性约束聚变(ICF)复杂物理问题的数值模拟中,需要耦合求解多个物理过程,主要包括流体力学运动、辐射传输(扩散近似或者输运描述)、核反应动力学、粒子输运和激光传输等。
在以往的多物理耦合集成程序研究过程中,采用上述的研究方式遇到了诸多困难。首先,不同时期、不同人员开发的应用程序集成难度大、周期长;其次,集成程序的耦合模式扩展性较差,不利于程序中物理建模与数学方法的改进;最后,“串行程序并行化”的研究模式,使得集成程序的并行化研究滞后。如何将求解不同类型物理过程的数值模拟程序耦合在一起,既要保证程序各自的独立开发,又要保证集成程序总体模拟的高效运行,并且实现集成程序的并行计算,成为多物理耦合集成程序研究的重要课题。采用JASMIN 框架的相关支撑技术,可以帮助集成程序克服以上困难,提高集成程序的数值模拟研究能力。
JASMIN 框架(J Adaptive Structured Mesh applications INfrastructure)的全称为并行自适应结构网格应用支撑软件框架,其中J代表研制单位北京应用物理与计算数学研究所。JASMIN 框架围绕科学与工程计算中普遍采用的结构网格或非结构拼接的结构网格,提供丰富的数据结构、算法构件与解法器,屏蔽高性能并行计算(数据通信、负载平衡)和网格自适应(细化、粗化)的实现,支撑多物理模型的耦合计算方法研究与高效并行应用程序研制。
本报告以惯性约束聚变(ICF)中辐射流体与粒子输运的多物理耦合过程为例,基于JASMIN 框架的“联邦计算”,联接辐射流体RH2D与粒子输运Sn2D两个串行程序,实现集成程序的多物理耦合并行数值模拟。集成程序RHSn2D中的邦元具有各自独立的网格划分与并行算法,同时借助框架技术,可以屏蔽邦元间的并行数据传递。算例表明,对于应用问题规模(90720个网格单元,辐射流体100个Patch,粒子输运2835个Patch,Sn方向48,16群),集成程序RHSn2D采用1024个处理器可以达到36%的并行效率。