丝阵Z箍缩X射线辐射机制数值分析.pdfVIP

第 24卷第 4期 强 激 光 与 粒 子 束 Vol_24，No．4 2012年 4月 H1GH POW ER LASER AND PARTICLE BEAMS Apr．，2012 文章编号： 1001—4322(2012)04—0776—04 丝阵Z箍缩 X射线辐射机制数值分析 段书超， 黄显宾， 杨礼兵， 李 晶，但加坤， 章征伟， 任晓东 (中国工程物理研究院 流体物理研究所 ，四川 绵 阳 621900) 摘 要 ： 基于简化模型编写了结构简单、运行高效的Z箍缩计算程序 ，并 已应用于 “阳”加速器上开展 的 z箍缩物理研究。在 “阳”加速器和 PTS的参数范围内进行了数值计算 ，结果表 明：对于质量较轻 的负载，激波 加热足产生辐射的主要原因；对于质量较重 的负载 ，辐射南激波加热和绝热压缩共 同主导 。激波加热功率与绝 热疗{缩功率之 比随丝阵质量增大而减小。烧蚀产生 的先驱质世减弱激波加热，对峰值功率的提高不利 ，但这部 分质量能抑制磁瑞利一泰勒不稳定性 。因此综合激波加热干¨磁瑞利一泰勒不稳定性两方面的考虑 ，应该存在一 个最优的烧蚀率 。 关键词 ： Z箍缩 ； 数值模拟 ； X射线辐射 ； 能量转换 中图分类 号： O532； O411 文 献标 志码 ： A doi：10．3788／HPLP0776 基于脉冲功率技术的 z箍缩 ，既可以作为一种聚变等离子体约束方式 ，也可 以作为一种强 X射线辐射源 ， 有着重要的科学研究价值和广泛的技术应用价值 。近二十年来 ，z箍缩研究在实验 、理论 、数值模拟三方面都 取得了快速进展l1]，鉴于此，国内多家单位也对 Z箍缩开展了多方面多种形式 的研究 ]。Z箍缩的多种应用 都是基于其 内爆末期的X射线辐射 ，因此深人理解 z箍缩的x射线辐射机制是 z箍缩研究的重要内容之一。 数值模拟是 z箍缩理论研究和负载设计及诊断分析的有力工具 。本文基于简化模型编写了结构简单 、运行高 效 的Z箍缩计算程序，并数值研究 了1～9MA百纳秒脉冲电流驱动的丝阵 Z箍缩中的X射线辐射机制 。 1 数值模型 z箍缩过程是个 3维 、非线性、多物理 、多尺度的复杂过程，其物理涉及物态相变 、磁流体力学、辐射输运 、 原子物理过程 ，其空间尺度从 m到 cm，其时间尺度从 fs到百 ns。为 了简单 ，将这一复杂过程简化近似为 1 维辐射磁流体力学过程 ，包含的物理过程有理想气体状态方程 、绝热近似、激波加热、辐射输运、单流体磁流体 力学。 数字离散方案为：流体力学部分采 VailLeer保单调格式 ]，磁扩散采用 Crank—Nicolson格式 ，辐射输 运采用 RTVD格式 。 ，时间步进为二阶 Runge—Kutta。按照这种方案编写的FOI—la结构简单 、运行高效 、守 恒 、保单调。 2 边界条件和数值诊断 大量实验证据表明，存在一个延长的烧蚀期决定其后 的内爆动力学[】H 。由于是单相流体近似 ，不能考 虑金属丝阵从固态到等离子态的初始化相变过程 ，丝阵质量在计算一开始就处在等离子态。由于是 1维近似 ， 丝阵烧蚀形成芯冕结构并进一步形成内爆壳结构的动力学过程也无法考虑，但丝烧蚀形成质量再分布效应能 在一定程度上予以考虑。将丝烧蚀作为边界上一个 向内流入的等离子体质量源，忽略质量源 的厚度和烧蚀边 界层的径向厚度 ，并忽略角向和轴向尺度 ，就得到一个 1维径向的质量分布，这个 1维分布对应于 2维或 3维 烧蚀形成的质量分布的角向和轴 向积分。这种质量人流边界条件 的表式为[1¨]：(pv)1： 一一fn(f)， ()一 丁广 ^]、 ( ) g·Cm～·ns～，其中定标系数a=O．2～0．4，定标指数志=1～2均为唯象系数，Dl为等离子体密 、 、L … ， 度，为等离子体流体速度。总质量按电流的平方定标，这是为了保持内爆动力学的相似性，使各算例问具有 可 比性 。 * 收稿 日期 ：2011-10—31； 修订 日期 ：201卜12—25 基金项且：中国工程物理研究院流体物理研究所发展基金项 目(SF 作者简介 ：段书超(198O～)，男，主要从事 z箍缩实验诊