-
研究热斑点火物理过程对于实现聚变点火和高增益聚变放能具有重要意义, 但是, 迄今为止的大部分相关研究都是针对等压构型预压缩等离子体进行的, 对等容构型预压缩等离子体中的热斑点火过程研究尚不充分. 本文针对双锥对撞点火方案产生的等容预压缩高密度等离子体, 建立了描述热斑边界演化和核聚变燃烧的半解析模型. 该模型表明, 在等容预压缩高密度等离子体中的热斑边界, 可以用对热斑产生的
$ {\mathrm{\alpha }} $ 粒子具有强烈吸收作用的激波波峰的位置来定义, 且等容预压缩等离子体中的热斑点火过程也存在$ {\mathrm{\alpha }} $ 粒子射程主导的自调节现象. 通过考虑$ {\mathrm{\alpha }} $ 粒子的空间不均匀沉积效应, 可以利用该模型描述等容预压缩等离子体中热斑内部的温度和密度演化. 使用该模型分析热斑在点火初期时刻的劳森参数和平均热斑温度发现, 在快电子总能量相同情况下, 能量较低的快电子束更有利于实现点火. 辐射流体模拟程序O-SUKI-N的验证计算表明, 本文提出的半解析模型的计算结果具有较强的合理性.-
关键词:
- 等容预压缩等离子体 /
- 快点火 /
- 热斑演化 /
- $ {\mathrm{\alpha }} $粒子沉积
The formation and evolution of hotspots is important for achieving ignition and high energy gain in inertial fusion process. However, most of relevant studies are carried out on pre-compressed plasmas with an isobaric configuration, the evolution of the hotspot in a plasma with an isochoric configuration is rarely studied. In this paper, a semi-analytical model is developed to describe the evolution of the hotspot boundary and propagation of fusion burning waves for a high-density pre-compressed plasma with an isochoric configuration in the double-cone ignition scheme. For the shock wave, the strong shock wave approximation and the quasi-isobaric approximation are reasonable. The quasi-isobaric approximation shows that as the plasma density behind the shock wave increases, the plasma temperature decreases. Because of these, the range of$ {\mathrm{\alpha }} $ -particles decreases rapidly behind the shock wave, forming an$ {\mathrm{\alpha }} $ -particle absorption peak. Therefore, considering that the hotspot is the main region where$ {\mathrm{\alpha }} $ -particles are produced and deposited, the position of the shock peak can be used to identify the boundary of the hotspot in a high-density plasma with an isochoric configuration. It also shows that a “self-regulating burning process” exists in the burning process of the isochoric hotspot, most of$ {\mathrm{\alpha }} $ -particles are deposited in the stable region and behind the shock, and finally, transport through the shock peak and heat the cold fuel, resulting in the temperature rising. In the high-density hotspots of plasma with an isochoric configuration, the deposition of α-particles behaves as an obvious non-uniform distribution effect. By analyzing the non-uniform deposition of α-particles, the deposition rate of α-particles at the edge of spherical uniform hotspot is calculated, then the temperature and density evolution of the isochoric hotspot can be well described. The model can be used to estimate the Lawson parameter of the hotspots at the end of the early stage of ignition. It is found that a lower fast electron energy is more beneficial to ignition and high gain operation of fusion plasma. It is also shown that the high density of the hotspots in the isochoric plasma will lead to a higher fusion burning rate, which can offset the negative influence of the shock wave and even achieve higher energy gain. The semi-analytical model is verified by the hydrodynamic simulations of O-SUKI-N.-
Keywords:
- pre-compressed isochoric plasma /
- fast ignition /
- hotspot evolution /
- $ {\mathrm{\alpha }} $-particles deposition
[1] Nuckolls J, Wood L, Thiessen A, Zimmerman G 1972 Nature 239 139
Google Scholar
[2] Tabak M, Hammer J, Glinsky M E, Kruer W L, Wilks S C, Woodworth J, Campbell E M, Perry M D, Mason R J 1994 Phys. Plasmas 1 1626
Google Scholar
[3] Zhang J, Wang W M, Yang X H, Wu D, Ma Y Y, Jiao J L, Zhang Z, Wu F Y, Yuan X H, Li Y T, Zhu J Q 2020 Philos. Trans. R. Soc. London, Ser. A 378 20200015
[4] Lan K 2022 Matter Radiat. Extrem. 7 055701
Google Scholar
[5] Tollefson J, Gibney E 2022 Nature 612 597
Google Scholar
[6] Lawson J D 1957 Proc. Trans. R. Soc. London, Ser. B 70 6
Google Scholar
[7] Atzeni S, Meyer-ter-Vehn J 2004 The Physics of Inertial Fusion (New York : Oxford University Press) pp44, 48
[8] Olson R E, Leeper R J, Yi S A, et al. 2016 J. Phys. Conf. Ser. 717 012042
Google Scholar
[9] Patel P K, Springer P T, Weber C R, et al. 2020 Phys. Plasmas 27 050901
Google Scholar
[10] 杨孟奇, 吴福源, 陈致博, 张翼翔, 陈一, 张晋川, 陈致真, 方志凡, Rafael Ramis, 张杰 2022 物理学报 71 225202
Google Scholar
Yang M Q, Wu F Y, Chen Z B, Zhang Y X, Chen Y, Zhang J C, Chen Z Z, Fang Z F, Rafael R, Zhang J 2022 Acta Phys. Sin. 71 225202
Google Scholar
[11] 张喆, 远晓辉, 张翌航, 刘浩, 方可, 张成龙, 刘正东, 赵旭, 董全力, 刘高扬, 戴羽, 谷昊琛, 李玉同, 郑坚, 仲佳勇, 张杰 2022 物理学报 71 155201
Google Scholar
Zhang Z, Yuan X H, Zhang Y H, Liu H, Fang K, Zhang C L, Liu Z D, Zhao X, Dong Q L, Liu G Y, Dai Y, Gu H C, Li Y T, Zheng J, Zhong J Y, Zhang J 2022 Acta Phys. Sin. 71 155201
Google Scholar
[12] Wu F, Yang X, Ma Y, et al. 2022 High Power Laser Sci. Eng. 10 e12
Google Scholar
[13] Song H, Wu F, Sheng Z, Zhang J 2023 Phys. Plasmas 30 092707
Google Scholar
[14] Xu Z, Wu F, Jiang B, Kawata S, Zhang J 2023 Nucl. Fusion 63 126062
Google Scholar
[15] Nakamura H, Uchibori K, Kawata S, Karino T, Sato R, Ogoyski A I 2022 Comput. Phys. Commun. 272 108223
Google Scholar
[16] Ramis R, Meyer-ter-Vehn J, Ramírez J 2009 Comput. Phys. Commun. 180 977
Google Scholar
[17] Wu F, Chu Y, Ramis R, et al. 2018 Matter Radiat. Extrem. 3 248
Google Scholar
[18] Krokhin O N, Rozanov V B 1973 Sov. J. Quantum Electron. 2 393
Google Scholar
[19] Fraley G S 1974 Phys. Fluids 17 474
Google Scholar
[20] Frenje J A, Grabowski P E, Li C K, et al. 2015 Phys. Rev. Lett. 115 205001
Google Scholar
[21] Spitzer L 2006 Physics of Fully Ionized Gases (2nd rev. ed.) (Mineola, NY: Dover Publications) p170
[22] Zylstra A B, Hurricane O A 2019 Phys. Plasmas 26 062701
Google Scholar
[23] Solodov A A, Betti R 2008 Phys. Plasmas 15 042707
Google Scholar
[24] Betti R, Hurricane O A 2016 Nat. Phys. 12 435
Google Scholar
[25] Zylstra A B, Hurricane O A, Callahan D A, et al. 2022 Nature 601 542
Google Scholar
-
图 2 (a) 等容预拼装等离子体密度温度示意图; (b) O-SUKI-N程序给出的等容预拼装等离子体在燃烧过程中的温度、密度和压强分布, 靶丸被划分为中心燃烧区域、激波区域和冷燃料区域. 此时初始热斑温度为$ 9\;{\mathrm{k}}{\mathrm{e}}{\mathrm{V}} $, 图中所取基准值${T}_{0}=5\;{\mathrm{k}}{\mathrm{e}}{\mathrm{V}}, $ $ {\rho }_{0}=100\;{\mathrm{g}}/{\mathrm{c}}{\mathrm{c}} $, ${P}_{0}={10}^{3}\;{\mathrm{G}}{\mathrm{b}}{\mathrm{a}}{\mathrm{r}} $
Fig. 2. (a) Density and temperature of isochoric pre-assembled plasma. (b) The temperature, density, and pressure distribution of the isochoric pre-assembly plasma in the burning process given by O-SUKI-N. The target is divided into central burning zone, shock wave, and cold fuel by density. The initial hotspot temperature is $ 9\;{\mathrm{k}}{\mathrm{e}}{\mathrm{V}} $, and the reference value is $ {T}_{0}=5\;{\mathrm{k}}{\mathrm{e}}{\mathrm{V}}, ~ {\rho }_{0}= $$ 100\;{\mathrm{g}}/{\mathrm{c}}{\mathrm{c}}, {P}_{0}={10}^{3}\;{\mathrm{G}}{\mathrm{b}}{\mathrm{a}}{\mathrm{r}} $.
图 3 O-SUKI-N程序给出的燃烧过程中等离子体参数和$ {\mathrm{\alpha }} $粒子吸收峰, 此时初始热斑温度为$ 9\;{\mathrm{k}}{\mathrm{e}}{\mathrm{V}} $ (a) 密度曲线(蓝色)与$ {\mathrm{\alpha }} $粒子沉积$ {n}_{{\mathrm{\alpha }}}{\omega }_{{\mathrm{\alpha }}} $分布曲线(红色); (b) 等离子体温度T、密度$\rho $和$ {\mathrm{\alpha }} $粒子沉积$ {n}_{{\mathrm{\alpha }}}{\omega }_{{\mathrm{\alpha }}} $分布
Fig. 3. The plasma parameters and α-particle absorption peak during burning process given by O-SUKI-N, and the initial hotspot temperature is $ 9\;{\mathrm{k}}{\mathrm{e}}{\mathrm{V}} $: (a) Density distribution (blue), α-particle deposition $ {n}_{{\mathrm{\alpha }}}{\omega }_{{\mathrm{\alpha }}} $ distribution (red); (b) temperature T, density of plasma $\rho $ and $ {\mathrm{\alpha }} $-particle deposition $ {n}_{{\mathrm{\alpha }}}{\omega }_{{\mathrm{\alpha }}} $.
图 5 初始热斑温度为$ 9\;{\mathrm{ }}{\mathrm{k}}{\mathrm{e}}{\mathrm{V}} $条件下, O-SUKI-N程序给出的燃烧过程中$ -\dfrac{1}{T}\dfrac{\partial T}{\partial r} $(蓝色实线), $ \dfrac{1}{\rho }\dfrac{\partial \rho }{\partial r} $(蓝色虚线)和$ {\mathrm{\alpha }} $粒子沉积分布$ {n}_{{\mathrm{\alpha }}}{\omega }_{{\mathrm{\alpha }}} $(红色实线)分布, 图中黑色虚线对应激波波峰的位置, $ {{\overline T}_{h}} $是热斑平均温度 (a)$ {{\overline T}_{h}}=10.99\;{\mathrm{keV}} $; (b) $ {{\overline T}_{h}}=14.46\;{\mathrm{keV}} $; (c) $ {{\overline T}_{h}}=22.65\;{\mathrm{keV}} $
Fig. 5. Distribution of $ -\dfrac{1}{T}\dfrac{\partial T}{\partial r} $ (blue solid line), $ \dfrac{1}{\rho }\dfrac{\partial \rho }{\partial r} $ (blue dashed line) and α-particle deposition $ {n}_{{\mathrm{\alpha }}}{\omega }_{{\mathrm{\alpha }}} $ (red solid line) in the burning process when the initial hotspot temperature is $ 9\;{\mathrm{ }}{\mathrm{k}}{\mathrm{e}}{\mathrm{V}} $ given by O-SUKI-N. The black dashed line corresponds to the position of the shock peak, and $ {{\overline T}_{{\mathrm{h}}}} $ is the average temperature of the hotspot (a) ${{\overline T}_{h}}=10.99\;{\mathrm{keV}} $; (b) $ {{\overline T}_{h}}=14.46\;{\mathrm{keV}} $; (c) ${{\overline T}_{h}}=22.65\;{\mathrm{keV}} $.
图 7 不均匀沉积模型与O-SUKI-N程序稳定区域边缘不均匀沉积函数$ {f}_{{\mathrm{n}}{\mathrm{o}}{\mathrm{n}}}\left(R\right) $和$ {f}_{{\mathrm{O}}{\mathrm{S}}{\mathrm{U}}}\left(R\right) $对比
Fig. 7. Comparison of non-uniform deposition function $ {f}_{{\mathrm{n}}{\mathrm{o}}{\mathrm{n}}}\left(R\right) $ and $ {f}_{{\mathrm{O}}{\mathrm{S}}{\mathrm{U}}}\left(R\right) $ at the stable zone edge of the non-uniform deposition mode model with O-SUKI-N.
图 8 快电子总能量$ 3\;{\mathrm{ }}{\mathrm{k}}{\mathrm{J}} $条件下半解析模型与数值模拟程序O-SUKI-N对比 (a) 点火初期结束时刻的热斑劳森参数$ P\tau $(蓝色)和热斑平均温度$ \overline{T} $(红色)对比; (b) $ {\mathrm{\alpha }} $粒子加热倍率$ {Q}_{{\mathrm{\alpha }}} $
Fig. 8. Comparison between the semi-analytical model and the simulation program O-SUKI-N when the total fast electron energy is 3 kJ: (a) Comparison of the Lawson parameter $ P\tau $ (blue) and average temperature $ \overline{T} $ (red) of the hotspot at the end of the early stage of ignition; (b) comparison of $ {Q}_{\alpha } $.
表 1 均匀球形热斑中$ {\mathrm{\alpha }} $粒子的产生与沉积区域
Table 1. Production and deposition of $ {\mathrm{\alpha }} $-particles in uniform spherical hotspots.
$ {\mathrm{\alpha }} $粒子产生区域 $ {\mathrm{\alpha }} $粒子沉积区域 符号 中心球体 ($ 0\leqslant r < R-{\mathrm{d}}R $) 中心球体 ($ 0\leqslant r < R-{\mathrm{d}}R $) $ {E}_{{\mathrm{d}}{\mathrm{e}}{\mathrm{p}}}(R-{\mathrm{d}}R\to R-{\mathrm{d}}R) $ 中心球体 $ (0\leqslant r < R-{\mathrm{d}}R $) 球壳 ($ R-{\mathrm{d}}R\leqslant r\leqslant R $) $ {E}_{{\mathrm{d}}{\mathrm{e}}{\mathrm{p}}}(R-{\mathrm{d}}R\to {\mathrm{d}}R) $ 球壳$ (R-{\mathrm{d}}R\leqslant r\leqslant R) $ 中心球体 ($ 0\leqslant r < R-{\mathrm{d}}R $) $ {E}_{{\mathrm{d}}{\mathrm{e}}{\mathrm{p}}}({\mathrm{d}}R\to R-{\mathrm{d}}R) $ 球壳 ($ R-{\mathrm{d}}R\leqslant r\leqslant R $) 球壳 ($ R-{\mathrm{d}}R\leqslant r\leqslant R $) $ {E}_{{\mathrm{d}}{\mathrm{e}}{\mathrm{p}}}({\mathrm{d}}R\to {\mathrm{d}}R) $ -
[1] Nuckolls J, Wood L, Thiessen A, Zimmerman G 1972 Nature 239 139
Google Scholar
[2] Tabak M, Hammer J, Glinsky M E, Kruer W L, Wilks S C, Woodworth J, Campbell E M, Perry M D, Mason R J 1994 Phys. Plasmas 1 1626
Google Scholar
[3] Zhang J, Wang W M, Yang X H, Wu D, Ma Y Y, Jiao J L, Zhang Z, Wu F Y, Yuan X H, Li Y T, Zhu J Q 2020 Philos. Trans. R. Soc. London, Ser. A 378 20200015
[4] Lan K 2022 Matter Radiat. Extrem. 7 055701
Google Scholar
[5] Tollefson J, Gibney E 2022 Nature 612 597
Google Scholar
[6] Lawson J D 1957 Proc. Trans. R. Soc. London, Ser. B 70 6
Google Scholar
[7] Atzeni S, Meyer-ter-Vehn J 2004 The Physics of Inertial Fusion (New York : Oxford University Press) pp44, 48
[8] Olson R E, Leeper R J, Yi S A, et al. 2016 J. Phys. Conf. Ser. 717 012042
Google Scholar
[9] Patel P K, Springer P T, Weber C R, et al. 2020 Phys. Plasmas 27 050901
Google Scholar
[10] 杨孟奇, 吴福源, 陈致博, 张翼翔, 陈一, 张晋川, 陈致真, 方志凡, Rafael Ramis, 张杰 2022 物理学报 71 225202
Google Scholar
Yang M Q, Wu F Y, Chen Z B, Zhang Y X, Chen Y, Zhang J C, Chen Z Z, Fang Z F, Rafael R, Zhang J 2022 Acta Phys. Sin. 71 225202
Google Scholar
[11] 张喆, 远晓辉, 张翌航, 刘浩, 方可, 张成龙, 刘正东, 赵旭, 董全力, 刘高扬, 戴羽, 谷昊琛, 李玉同, 郑坚, 仲佳勇, 张杰 2022 物理学报 71 155201
Google Scholar
Zhang Z, Yuan X H, Zhang Y H, Liu H, Fang K, Zhang C L, Liu Z D, Zhao X, Dong Q L, Liu G Y, Dai Y, Gu H C, Li Y T, Zheng J, Zhong J Y, Zhang J 2022 Acta Phys. Sin. 71 155201
Google Scholar
[12] Wu F, Yang X, Ma Y, et al. 2022 High Power Laser Sci. Eng. 10 e12
Google Scholar
[13] Song H, Wu F, Sheng Z, Zhang J 2023 Phys. Plasmas 30 092707
Google Scholar
[14] Xu Z, Wu F, Jiang B, Kawata S, Zhang J 2023 Nucl. Fusion 63 126062
Google Scholar
[15] Nakamura H, Uchibori K, Kawata S, Karino T, Sato R, Ogoyski A I 2022 Comput. Phys. Commun. 272 108223
Google Scholar
[16] Ramis R, Meyer-ter-Vehn J, Ramírez J 2009 Comput. Phys. Commun. 180 977
Google Scholar
[17] Wu F, Chu Y, Ramis R, et al. 2018 Matter Radiat. Extrem. 3 248
Google Scholar
[18] Krokhin O N, Rozanov V B 1973 Sov. J. Quantum Electron. 2 393
Google Scholar
[19] Fraley G S 1974 Phys. Fluids 17 474
Google Scholar
[20] Frenje J A, Grabowski P E, Li C K, et al. 2015 Phys. Rev. Lett. 115 205001
Google Scholar
[21] Spitzer L 2006 Physics of Fully Ionized Gases (2nd rev. ed.) (Mineola, NY: Dover Publications) p170
[22] Zylstra A B, Hurricane O A 2019 Phys. Plasmas 26 062701
Google Scholar
[23] Solodov A A, Betti R 2008 Phys. Plasmas 15 042707
Google Scholar
[24] Betti R, Hurricane O A 2016 Nat. Phys. 12 435
Google Scholar
[25] Zylstra A B, Hurricane O A, Callahan D A, et al. 2022 Nature 601 542
Google Scholar
计量
- 文章访问数: 1488
- PDF下载量: 90
- 被引次数: 0