混合流体-PIC物理建模及其在激光聚变中的应用

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

计算物理 ›› 2023, Vol. 40 ›› Issue (2): 159-168.DOI: 10.19596/j.cnki.1001-246x.8602

所属专题：贺贤土院士从事科学研究工作60周年暨激光聚变相关研究进展专刊

• 贺贤土院士从事科学研究工作60周年暨激光聚变相关研究进展专刊 • 上一篇下一篇

混合流体-PIC物理建模及其在激光聚变中的应用

蔡洪波¹^,²(), 许育培³, 姚沛霖⁴, 张恩浩³, 黄翰霄⁵, 朱少平¹^,³, 贺贤土¹^,²

1. 北京应用物理与计算数学研究所, 北京 100094
2. 北京大学应用物理与技术研究中心, 工学院, 北京 100871
3. 中国工程物理研究院研究生院, 北京 100088
4. 清华大学工程物理系, 北京 100084
5. 北京大学应用物理与技术研究中心, 物理学院, 北京 100871

收稿日期:2022-07-29 出版日期:2023-03-25 发布日期:2023-07-05
作者简介:
蔡洪波, 男, 博士, 研究员, 博士生导师, 研究方向为激光聚变物理, E-mail: cai_hongbo@iapcm.ac.cn
基金资助:
国家自然科学基金(11975055)

Hybrid Fluid-PIC Modeling and Its Application in Laser Fusion

Hongbo CAI¹^,²(), Yupei XU³, Peilin YAO⁴, Enhao ZHANG³, Hanxiao HUANG⁵, Shaoping ZHU¹^,³, Xiantu HE¹^,²

1. Institute of Applied Physics and Computational Mathematics, Beijing 100094, China
2. Center for Applied Physics and Technology, and College of Engineering, Peking University, Beijing 100871, China
3. Graduate School, China Academy of Engineering Physics, Beijing 100088, China
4. Department of Engineering Physics, Tsinghua University, Beijing 100084, China
5. Center for Applied Physics and Technology, School of Physics and College of Engineering, Peking University, Beijing 100871, China

Received:2022-07-29 Online:2023-03-25 Published:2023-07-05

摘要/Abstract

摘要：

本文探索性地建立混合流体-PIC(particle-in-cell)模拟方法, 电子用无质量的流体描述, 多组分离子用PIC粒子描述; 流体运动通过求解电子磁流体方程组来获得, 电磁场通过求解Ohm定律和Faraday定律来获得。针对激光聚变等离子体条件, 利用混合流体-PIC模拟研究高能量密度条件下冲击波结构及其特性、流体力学界面不稳定性演化的影响等问题。混合流体-PIC物理建模为研究等离子体效应对冲击波和界面不稳定性的影响提供了新的研究手段。

关键词: 激光聚变, 混合流体-particle-in-cell物理建模, 冲击波, 流体力学界面不稳定性

Abstract:

The influence of plasma effect on the shock wave and hydrodynamic instabilities is a widely concerned problem in the current research of laser fusion. However, due to the limitations of the numerical simulation methods, there is still a lack of research tools on this issue. In this work, a hybrid fluid PIC(particle-in-cell) simulation method is established tentatively. Electrons are described by a massless fluid, and multi-component ions are described by PIC method; The fluid motion is obtained by solving the equations of electro-magnetohydrodynamic, and the electromagnetic fields are obtained by solving Ohm's law and Faraday's law. Aiming at the plasma condition of laser fusion, we use hybrid fluid-PIC simulation to study the shock wave structure and its characteristics in high energy density conditions, and the influence of plasma effect on the evolution of hydrodynamic instabilities. Hybrid fluid-PIC physical modeling provides a new research method for studying the effect of plasma effect on shock wave and hydrodynamic instabilities under high energy density.

Key words: laser fusion, hybrid fluid-particle-in-cell modeling, shock wave, hydrodynamic instability

蔡洪波, 许育培, 姚沛霖, 张恩浩, 黄翰霄, 朱少平, 贺贤土. 混合流体-PIC物理建模及其在激光聚变中的应用[J]. 计算物理, 2023, 40(2): 159-168.

Hongbo CAI, Yupei XU, Peilin YAO, Enhao ZHANG, Hanxiao HUANG, Shaoping ZHU, Xiantu HE. Hybrid Fluid-PIC Modeling and Its Application in Laser Fusion[J]. Chinese Journal of Computational Physics, 2023, 40(2): 159-168.

图/表 4

图1 混合流体-PIC程序Ascent-H的物理建模与流程图

Fig.1 Physical model of the Hybrid fluid-PIC code Ascent-H

图2 马赫数(a) M = 1.43和(c) M = 3.79时等离子体冲击波中的电子、离子的温度空间分布(黑色、红色、蓝色实线分别表示电子、D+、3He2+)；(b)、(d) 离子相空间分布(黑色实线为离子密度分布，红色实线为电场分布。)

Fig.2 (a), (c) Electron and ion temperature profile in plasma shock wave with Mach number (a) M = 1.43 and (c) M = 3.79 (Black, red, and blue solid lines denote the electron, D+ and 3He2+ respectively.); (b), (d) Ion phase space of the shock (The black solid line denotes the ion density profile and the red solid line denotes the electric field.)

图3 从离子热涨落增长起来的瑞利-泰勒不稳定性密度分布

Fig.3 Density snapshots of Rayleigh-Taylor instability growth from the ion thermal fluctuations

图4 不同时刻激光功率(a)~(c)、电子压力(d)~(f)、自生磁场(g)~(i)的空间分布图(a)、(d)、(g)中t = 2 500 T0；(b)、(e)、(h) t=10 000 T0；(c)、(f)、(i) t = 20 000 T0

Fig.4 Profiles of Laser intensity, electron pressure, and self-generated magnetic fields at different times (a), (d), (g) t=2 500 T0; (b), (e), (h) t=10 000 T0; (c), (f), (i) t=20 000 T0

参考文献 31

1	HURRICANE O A , CALLAHAN D A , CASEY D T , et al. Fuel gain exceeding unity in an inertially confined fusion implosion[J]. Nature, 2014, 506, 343- 348. DOI
2	HURRICANE O A , CALLAHAN D A , SPRINGER P T , et al. Beyond alpha-heating: Driving inertially confined fusion implosions toward a burning-plasma state on the National Ignition Facility[J]. Plasma Phys Control Fusion, 2019, 61, 014033. DOI
3	RINDERKNECHT H G , AMENDT P A , ROSENBERG M J , et al. Ion kinetic dynamics in strongly-shocked plasmas relevant to ICF[J]. Nuclear Fusion, 2017, 57 (6): 064014.
4	宋鹏, 翟传磊, 李双贵, 等. 激光间接驱动惯性约束聚变二维总体程序——LARED集成程序[J]. 强激光与粒子束, 2015, 27, 032007. DOI
5	MARINAK M M , KERBEL G D , GENTILE N A , et al. Three-dimensional HYDRA simulations of National Ignition Facility targets[J]. Phys Plasmas, 2001, 8, 2275- 2280. DOI
6	RINDERKNECHT H G , AMENDT P A , WILKS S C , et al. Kinetic physics in ICF: Present understanding and future directions[J]. Plasma Phys Control Fusion, 2018, 60, 064001. DOI
7	XU A G , SONG J H , CHEN F , et al. Modeling and analysis methods for complex fields based on phase space[J]. Chinese Journal of Computational Physics, 2021, 38 (6): 631- 660.
8	ZHANG H , WU S Z , ZHOU C T , et al. Numerical method of relativistic Fokker-planck equation for energy deposition of fast electrons[J]. Chinese Journal of Computational Physics, 2017, 34 (5): 555- 562. DOI
9	AMENDT P . Entropy generation from hydrodynamic mixing in inertial confinement fusion indirect-drive targets[J]. Phys Plasmas, 2021, 28, 072701. DOI
10	HOPKINS L F , PAPE S LE , DIVOL L , et al. Near-vacuum hohlraums for driving fusion implosions with high density carbon ablators[J]. Phys Plasmas, 2015, 22, 056318. DOI
11	RYGG J R , FRENJE J A , LI C K , et al. Time-dependent nuclear measurements of mix in inertial confinement fusion[J]. Phys Rev Lett, 2007, 98, 215002. DOI
12	GU J F , GE F J , DAI Z S , et al. Influence of penetrating defect in DD ice layer on ICF cryogenic capsule implosion performance[J]. Chinese Journal Of Computational Physics, 2022, 39 (3): 265- 276.
13	GLASSER A H , SOVINEC C R , NEBEL R A , et al. The NIMROD code: A new approach to numerical plasma physics[J]. Plasma Phys Control Fusion, 1999, 41, A747- A755. DOI
14	FU G Y , PARK W , STRAUSS H R , et al. Global hybrid simulations of energetic particle effects on the n = 1 mode in Tokamaks: Internal kink and fishbone instability[J]. Phys Plasmas, 2006, 13, 052517. DOI
15	COWEE M M , WINSKE D , GARY S P . Hybrid simulations of plasma transport by Kelvin-Helmholtz instability at the magnetopause: Density variations and magnetic shear[J]. J Geophy Res, 2010, 115, A06214.
16	XU H , ZHUO H B , YANG X H , et al. Hybrid particle-in-cell/fluid model for hot electron transport in dense plasmas[J]. Chinese Journal of Computational Physics, 2017, 34 (5): 505- 525. DOI
17	JONES M E , LEMONS D S , MASON R J , et al. A Grid-based Coulomb collision model for PIC codes[J]. J Comput Phys, 1996, 123, 169- 181. DOI
18	BRAGINSKII S I . Transport processes in a plasma: In reviews of plasma physics[M]. New York: Consultants Bureau, 1965: 205- 311.
19	JI J Y , HELD E D . Closure and transport theory for high-collisionality electron-ion plasmas[J]. Phys Plasmas, 2013, 20, 042114. DOI
20	LEE Y T , MORE R M . An electron conductivity model for dense plasma[J]. Phys Fluids, 1984, 27, 1273- 1286. DOI
21	CAI H B , YAN X X , YAO P L , et al. Hybrid fluid-particle modeling of shock-driven hydrodynamic instabilities in a plasma[J]. Matter Radiat Extremes, 2021, 6, 035901. DOI
22	LEMBEGE B , SIMONET F . Hybrid and particle simulations of an interface expansion and of collisionless shock: A comparative and quantitative study[J]. Phys Plasmas, 2001, 8, 3967- 3981. DOI
23	MEEZAN N B , WOODS D T , IZUMI N , et al. Evidence of restricted heat transport in National Ignition Facility hohlraums[J]. Phys Plasmas, 2020, 27, 102704. DOI
24	ZHANG E H , CAI H B , ZHANG W S , et al. Influence of the electron thermal conduction and ion kinetic effects on the structure of collisional plasma shocks[J]. Phys Plasmas, 2022, 29, 082110. DOI
25	SCHURTZ G P , NICOLAI PH D , BUSQUET M . A nonlocal electron conduction model for multidimensional radiation hydrodynamics codes[J]. Phys Plasmas, 2000, 7, 4238- 4249. DOI
26	张钧, 常铁强. 激光核聚变靶物理基础[M]. 北京: 国防工业出版社, 2004.
27	SENTOKU Y , KEMP A J . Numerical methods for particle simulations at extreme densities and temperatures: Weighted particles, relativistic collisions and reduced currents[J]. J Comput Phys, 2008, 227, 6846- 6861. DOI
28	NANBU K , YONEMURA S . Weighted particles in Coulomb collision simulations based on the theory of a cumulative scattering angle[J]. J Comput Phys, 1998, 145, 639- 654. DOI
29	ROSS J S , HIGGINSON D P , RYUTOV D , et al. Transition from collisional to collisionless regimes in interpenetrating plasma Flows on the National Ignition Facility[J]. Phys Rev Lett, 2017, 118, 185003. DOI
30	HIGGINSON D P , LINK A , SCHMIDT A . A pairwise nuclear fusion algorithm for weighted particle-in-cell plasma simulations[J]. J Comput Phys, 2019, 388, 439- 453. DOI
31	YAN X X , CAI H B , YAO P L , et al. Ion kinetic effects on the evolution of Richtmyer-Meshkov instability and interfacial mix[J]. New J Phys, 2021, 23, 053010. DOI

[1]	吴宗铎, 满庆洋, 严谨, 庞建华, 孙一方. 三介质Mie-Grüneisen混合模型在水下爆炸防护层问题中的应用[J]. 计算物理, 2023, 40(5): 556-569.
[2]	张文帅, 蔡洪波, 张恩浩, 杜报, 邹士阳, 朱少平. 马赫数对等离子体碰撞冲击波结构的影响———全动理学模拟研究[J]. 计算物理, 2023, 40(2): 199-209.
[3]	薛社生, 李守先. 正冲击波后固壁表面颗粒的上升运动[J]. 计算物理, 2021, 38(3): 280-288.
[4]	李双贵, 杭旭登, 杨容, 宋鹏, 翟传磊, 齐进. LARED集成程序辐射输运模拟的性能优化[J]. 计算物理, 2017, 34(3): 320-326.
[5]	冯士德, 任荣彩, 毛江玉. 一个完全守恒的格子Boltzmann模型[J]. 计算物理, 2001, 18(5): 407-411.
[6]	田宙, 郭永辉, 郝保田. 气液爆轰波及冲击波传播的数值研究[J]. 计算物理, 2000, 17(S1): 131-136.
[7]	张钧, 赖东显, 张维岩. 辐射烧蚀界面的单模Rayleigh-Taylor不稳定性分析[J]. 计算物理, 1998, 15(4): 483-488.
[8]	郝保田, 张海波. 三维多级扩散室对冲击波传播规律的影响[J]. 计算物理, 1997, 14(S1): 572-573,571.
[9]	曾先才, 李沄生, 常铁强. 神光Ⅰ强激光腔靶X光辐射烧蚀铝平面靶的理论研究[J]. 计算物理, 1996, 13(2): 217-222.
[10]	张钧, 陈光南, 常铁强, 裴文兵, 尤希文, 隋成之, 徐少则, 古培俊, 章实, 张兴宏. 激光产生X光的数值模拟研究[J]. 计算物理, 1994, 11(1): 1-8.

混合流体-PIC物理建模及其在激光聚变中的应用

Hybrid Fluid-PIC Modeling and Its Application in Laser Fusion

RichHTML

PDF

可视化

摘要/Abstract

引用本文

使用本文

图/表 4

参考文献 31

相关文章 10

编辑推荐

Metrics

作者中心

审稿中心

期刊浏览

期刊介绍