惯性约束聚变内爆中基于多块结构网格的高效辐射扩散并行算法

doi:10.19596/j.cnki.1001-246x.8419

计算物理 ›› 2022, Vol. 39 ›› Issue (3): 277-285.DOI: 10.19596/j.cnki.1001-246x.8419

惯性约束聚变内爆中基于多块结构网格的高效辐射扩散并行算法

胡晓燕(), 范征锋*()

北京应用物理与计算数学研究所, 北京 100088

收稿日期:2021-07-07 出版日期:2022-05-25 发布日期:2022-09-02
通讯作者: 范征锋
作者简介:
胡晓燕, 副研究员, 研究方向为辐射流体数值算法, E-mail: hu_xiaoyan@iapcm.ac.cn
基金资助:
国家自然科学基金(11871114)

An Efficient Parallel Algorithm on Multi-block Structure Meshes for Radiation Diffusion in Inertial Confinement Fusion Implosion

Xiaoyan HU(), Zhengfeng FAN*()

Institute of Applied Physics and Computational Mathematics, Beijing 100088, China

Received:2021-07-07 Online:2022-05-25 Published:2022-09-02
Contact: Zhengfeng FAN

摘要/Abstract

摘要：

针对三维球形靶丸内爆高效模拟需求和传统笛卡尔正交网格上辐射加源困难的问题, 发展一种多块结构非正交网格生成方法, 并基于此种计算网格提出高效的三维扩散格式并行算法, 将其应用于辐射流体方程组的求解和三维内爆不对称性的数值模拟, 数值结果显示了算法的有效性。并行性能测试显示该算法可扩展到5400个核上, 并行效率达到69%。

关键词: 多块结构网格, 扩散方程, 惯性约束核聚变内爆, 并行算法

Abstract:

For high efficiency simulation of three-dimensional spherical implosion and the difficulty of adding radiation source on traditional Cartesian orthogonal grid, a multi-block non-orthogonal mesh generation method is developed. An efficient parallel algorithm of three-dimensional diffusion scheme on these meshes is proposed. It was applied to the simulation of radiation hydrodynamic equations and the asymmetry study in three-dimensional implosion problems. Effectiveness of the algorithm is shown. Scalability is demonstrated with a parallel efficiency above 69% on 5400 processors.

Key words: multi-block structure meshes, diffusion equation, inertial confinement fusion implosion, parallel algorithm

胡晓燕, 范征锋. 惯性约束聚变内爆中基于多块结构网格的高效辐射扩散并行算法[J]. 计算物理, 2022, 39(3): 277-285.

Xiaoyan HU, Zhengfeng FAN. An Efficient Parallel Algorithm on Multi-block Structure Meshes for Radiation Diffusion in Inertial Confinement Fusion Implosion[J]. Chinese Journal of Computational Physics, 2022, 39(3): 277-285.

图/表 10

图1 (a) 球几何网格生成; (b) 六块球面网格; (c) 计算控制体

Fig.1 (a) Spherical grid generation; (b) Six-block spherical grid; (c) Computational control volume

图2 (a) 单块球面上的9点图和(b)多块块间边界

Fig.2 (a) 9-point on a single block and (b) boundary between multiple blocks

图3 块内区域分解

Fig.3 Domain decomposition in a block

图4 单块半径方向l=0和l=1的进程

Fig.4 l=0 and l = 1 processes on radius

图5 空间收敛阶

Fig.5 Order of convergence in space

图6 减速阶段模型初始状态

Fig.6 Initial states in deceleration phase

图7 0.4 ns时的密度、辐射温度和速度分布

Fig.7 Density, radiation temperature and velocity profiles at 0.4 ns

图8 (a) 初始靶丸构型; (b) 辐射温度源; (c) 驱动不对称加源

Fig.8 (a) Initial target configuration; (b) Radiation temperature source; (c) Drive asymmetric source profile

图9 阻滞时刻电子温度及密度等值面

Fig.9 Isosurfaces of electron temperature and density at stagnation time

表1 并行性能测试表

Table 1 Parallel performance

CPU核数	计算时间/min	并行效率/%
675	437	100
1 350	237	92.2
2 700	133	82.1
5 400	79	69.1

参考文献 18

1	MIHALAS D , MIHALAS B W , et al. Foundations of radiation hydrodynamics[M]. New York: Oxford University Press, Oxford, 1984: 1- 10.
2	阿蔡塞. 惯性约束聚变物理[M]. 沈百飞, 译. 北京: 科学出版社, 2008: 37-61.
3	袁光伟, 盛志强, 杭旭登, 等. 扩散方程计算方法[M]. 北京: 科学出版社, 2015: 1- 20.
4	YUAN G W , HANG X D , et al. Progress in numerical methods for radiation diffusion equations[J]. Chinese Journal of Computational Physics, 2009, 26 (4): 475- 500.
5	李德元, 陈广南. 抛物型方程差分方法引论[M]. 北京: 科学出版社, 1995: 94- 177.
6	CHEN G N , LI D Y , WANG Z S . Discrete scheme based on variation principle for three-dimensional diffusion equation[J]. Chinese Journal of Computational Physics, 2003, 20 (4): 291- 297.
7	CHEN G N , LI D Y , WANG Z S . Difference scheme by integral interpolation method for three-dimensional diffusion equation[J]. Chinese Journal of Computational Physics, 2003, 20 (3): 205- 209.
8	郭红, 莫则尧, 张爱清. JASMIN框架中多块结构网格拼接并行计算及其应用[J]. 计算机工程与科学, 2012, 34 (8): 69- 74. DOI
9	HU X Y , NI G X , FAN Z F , et al. Algorithm of radiation hydrodynamics with nonorthogonal mesh for 3D implosion problem[J]. Journal of Computational Physics, 2021, 437 (15): 11309.
10	PEI W B . The construction of simulation algorithms for laser fusion[J]. Commun Comput Phys, 2007, 2 (2): 255- 270.
11	YU C X , XUE C , LIU J , et al. Multiple eigenmodes of the Rayleigh-Taylor instability observed for a fluid interface with smoothly varying density[J]. Phys Rev E, 2018, 97 (1): 013102. DOI
12	WU J F , MIAO W Y , WANG L F , et al. Indirect-drive ablative Rayleigh-Taylor growth experiments on the Shenguang-II laser facility[J]. Physics of Plasmas, 2014, 21 (4): 042707. DOI
13	FAN Z F , HE X T , LIU J , et al. A wedged-peak-pulse design with medium fuel adiabat for indirect-drive fusion[J]. Physics of Plasmas, 2014, 21 (10): 100705. DOI
14	LI Y S , GU J F , WU C S , et al. Effects of the P2 M-band flux asymmetry of laser-driven gold hohlraums on the implosion of ICF ignition capsule[J]. Physics of Plasmas, 2016, 23 (7): 072705. DOI
15	FAN Z F , LIU Y Y , LIU B , et al. Non-equilibrium between ions and electrons inside hot spots from National Ignition Facility experiments[J]. Matter and Radiation at Extremes, 2017, 2 (1): 3- 8. DOI
16	ZHANG C B , YU C X , YANG C , et al. A theoretical model for low-mode asymmetries in ICF implosions[J]. Physics of Plasmas, 2019, 26 (2): 022707. DOI
17	SERDUKE F J D , MINGUEZ E , DAVIDSON S J , et al. WorkOp-IV summary: Lessons from iron opacities[J]. Journal of Quantitative Spectroscopy, 2000, 65 (1): 527.
18	徐锡申, 张万箱. 实用物态方程理论导引[M]. 北京: 科学出版社, 1986.

[1]	冯舒婷, 戴厚平, 宋通政. 二维空间分数阶反应扩散方程组的格子Boltzmann方法[J]. 计算物理, 2022, 39(6): 666-676.
[2]	朱家莉, 尚月强. 不可压缩流的并行两水平稳定有限元算法[J]. 计算物理, 2022, 39(3): 309-317.
[3]	蔡颖, 张存波, 刘旭, 范征锋, 刘元元, 徐小文, 张爱清. 二维球坐标系中子输运方程的一种并行SN算法[J]. 计算物理, 2022, 39(2): 143-152.
[4]	张淑娴, 杭旭登. 二维扩散方程的局部二阶线性保极值节点计算方法及应用[J]. 计算物理, 2022, 39(1): 17-32.
[5]	魏雪丹, 戴厚平, 李梦军, 郑洲顺. 一维空间Riesz分数阶对流扩散方程的格子Boltzmann方法[J]. 计算物理, 2021, 38(6): 683-692.
[6]	徐金景, 袁光伟. 多流管网格上九点格式的节点插值方法[J]. 计算物理, 2021, 38(2): 153-164.
[7]	赵菲, 盛志强, 袁光伟. 基于二阶格式的有限体积保正格式[J]. 计算物理, 2020, 37(4): 379-392.
[8]	丁琪, 尚月强. 非定常Navier-Stokes方程基于两重网格离散的有限元并行算法[J]. 计算物理, 2020, 37(1): 10-18.
[9]	张守慧, 梁栋. 二维抛物型问题的Strang型交替分段区域分裂格式[J]. 计算物理, 2018, 35(4): 413-428.
[10]	朱晓钢, 聂玉峰, 王俊刚, 袁占斌. 分数阶对流扩散方程的特征有限元方法[J]. 计算物理, 2017, 34(4): 417-424.
[11]	祁美玲, 杨琼, 王苍龙, 田园, 杨磊. 结构材料辐照损伤的分子动力学程序GPU并行化及优化[J]. 计算物理, 2017, 34(4): 461-467.
[12]	郭少冬, 章明宇, 周海兵, 熊俊, 张树道. 三维非匹配网格上的扩散方程数值求解[J]. 计算物理, 2017, 34(1): 19-28.
[13]	吕桂霞, 孙顺凯. 扩散方程的有限方向差分方法[J]. 计算物理, 2015, 32(6): 649-661.
[14]	舒适, 岳孝强, 周志阳, 徐小文. SAMR网格上扩散方程有限体格式的逼近性与两层网格算法[J]. 计算物理, 2014, 31(4): 390-402.
[15]	邓志红, 孙玉良, 李富, Rizwan-uddin. 改进的节块展开方法求解圆柱几何对流扩散方程[J]. 计算物理, 2013, 30(4): 475-482.

惯性约束聚变内爆中基于多块结构网格的高效辐射扩散并行算法

An Efficient Parallel Algorithm on Multi-block Structure Meshes for Radiation Diffusion in Inertial Confinement Fusion Implosion

RichHTML

PDF

可视化

摘要/Abstract

引用本文

使用本文

图/表 10

参考文献 18

相关文章 15

编辑推荐

Metrics

作者中心

审稿中心

期刊浏览

期刊介绍