高速CT阻滞产生X射线辐射的数值模拟研究-核聚变与等离子体物理.PDF

第 28 卷 第 1 期 核 聚 变 与 等 离 子 体 物 理 Vol.28, No.1 2 0 0 8 年 3 月 Nuclear Fusion and Plasma Physics March 2008 文章编号：0254−6086(2008)01−0034−06 高速 CT 阻滞产生 X 射线辐射的数值模拟研究 王震涛，王尚武，肖亚斌 (国防科学技术大学光电科学与工程学院，长沙 410073) 摘 要：以 Ar 等离子体为例，应用一维辐射流体力学模型，对高速度、高密度的紧凑等离子体环(CT)撞击 静态靶减速产生软 X 射线辐射的物理过程进行了数值模拟研究。利用分裂格式方法计算得到了描述碰撞物理过程 的各力学量随时间的演化关系，以及固定 CT 动能时使停滞 CT 成为最有效辐射源的初始速度和质量的最优搭配。 对计算过程中相关物理量的跟踪分析表明，在 1%的误差范围内计算格式满足能量守恒方程。 关键词：紧凑等离子体环；软 X 射线辐射；辐射流体力学；辐射产额；辐射塌缩 中图分类号：O536 文献标识码：A 1 引言 紧凑等离子体环(CT) 由其内部电流产生的环 向和极向磁场约束，处在 Woltjer-Taylor(W-T)态， 使得 CT 可以被稳定地压缩和加速而具有高能量密 度。因此，CT 可以作为聚变的驱动源、软 X 射线 源、高功率微波源和等离子体开关等[1~5]。高速度、 高密度的 CT 与静态靶碰撞减速产生高功率的软 X 射线辐射，可以作为武器效应试验、惯性约束核聚 变试验以及武器物理试验的能量源。 本文对高速度、高密度的 CT 与静态靶碰撞产生 软 X 射线辐射的物理过程进行了数值模拟。国外特 图 1 一维数值模拟的理想几何模型 别是美国在这方面取得了一系列成果，其中研究比 较成熟的主要有劳伦斯利弗莫尔国家实验室和美 在一维情况下，不考虑磁场影响，并假定静态 国空军的菲利浦实验室[6~9] ，日本在上世纪 90 年代 靶是理想的，没有电子和辐射热传导，辐射输运采 也开展了相关研究[10] ，而国内研究尚处于起步阶段。 用限流扩散近似。假设物质处于局域热力学平衡状 态，采用拉格朗日方法数值求解 RHD 方程组，对 2 基本原理 CT 与静态靶碰撞减速产生软 X 射线辐射的物理过 一维数值模拟的理想几何模型如图 1 所示，用 程进行数值模拟。 高速度、高密度的 CT 撞击静态靶，在撞击之前 CT 在一维平面对称坐标系下，定义某一流体质团 动能主要由离子携带，温度相对较低；碰撞压缩使 的拉格朗日坐标(质量坐标) z ( ) CT 受热，动能转化为离子热能，而离子热能通过 m ∫0 ρξ, t 0 dξ ，用 与电子相互作用转化为电子内能。此时，CT 温度 t=0 时刻从欧拉坐标 z 到 z=0 之间流体质团的初始 很高且高度电离，高能电子跃迁产生光子以辐射能 总质量来定义 m 。下面的方程构成了闭合的 RHD 形式释放出来。 方程组。 收稿日期：2007−05−16；修订日期：2007−12−17 基金项目：国防科学技术大学预研基金资助项目 作者简介：王震涛(1982−) ，男，山东淄博人，硕士研究生，研究方向：粒子束物理与技