一维磁化套筒惯性聚变模拟程序的设计与校验

在调研跟踪国外磁化套筒惯性聚变(MagLIF)研究领域最新进展的基础上,以一维磁流体力学方程为基础,结合氘氚燃料能量交换方程,建立一维磁化套筒惯性聚变物理模型,并使用Fortran语言编写完成一维数值模拟程序MagLIF-1D;通过与圣地亚实验室LASNEX,HYDRA等程序计算结果对比展示,完成程序校验工作,讨论后认为程序计算结果存在差异的主要原因可能来自不同程序对于材料状态方程库的选择;通过计算,MagLIF-1D程序可以直接获得内爆速度、燃料压力、燃料密度、聚变产额等关键物理量,这为后续更好地开展磁化套筒惯性聚变实验设计提供了有力工具。     

磁化套筒惯性聚变一维集成化数值模拟

磁化套筒惯性聚变(magnetized liner inertial fusion,MagLIF)结合了传统磁约束聚变和惯性约束聚变的优点,理论上可以显著地降低聚变实现的难度,具有极大的应用潜力.以研究MagLIF中的关键问题为目标,建立能够综合考虑磁化、预加热、套筒内爆、聚变反应、端面效应、磁通压缩等多种复杂机制在内的集成化物理模型,特别是通过引入流体喷射模型,使得可以在一维计算条件下考虑具有二维特性的端面损失情况,并额外考虑Nernst扩散项对磁通损失的影响.在此基础上编写实现一维集成化MagLIF数值模拟程序MIST(magnetic implosion simulation tools),与FP-1装置(2 MA,7.2μs)上铝套筒内爆实验结果的对比验证了程序磁流体模块的正确性;将聚变模块纳入后与国外同类程序LASNEX和HYDRA计算结果进行整体比较,所得数值结果总体接近,主要差异体现在燃料温度的计算上,对可能影响的原因进行了简要分析.所建立的集成化模型与程序将为未来开展MagLIF聚变实验研究提供坚实的理论基础和重要工具.     

磁化套筒惯性聚变研究进展

磁化套筒惯性聚变（MagLIF）是一种新的聚变构型,它结合了传统惯性约束聚变和磁约束聚变的优点,理论上可以显著地降低聚变实现的难度,未来必将朝着点火的目标进一步发展,具备极大的应用潜力。针对这一特殊构型,分别从理论、实验和工程三个部分介绍了国际上该领域主要的研究进展,内容覆盖理论研究、数值模拟、实验加载、测量与诊断、负载设计与加工、分解实验、构型改进等多个方面,通过该文能够对该领域的研究现状有相对完善的了解,对未来发展趋势也有一定的认知。

     

磁化套筒惯性聚变典型物理过程及特征参量

磁化套筒惯性聚变(MagLIF)结合了传统磁约束聚变(MCF)与惯性约束聚变(ICF)的优势,理论上在有限的驱动能力下可以有效降低聚变实现的难度,具有极大的应用潜力。基于一维集成化物理模型编写了数值模拟程序,以ZR装置典型驱动能力27 MA为出发点,以时间演化为顺序,通过数值模拟系统性地总结分析了典型负载参数下MagLIF构型初始化、加速内爆及迟滞3个关键过程中重要特征参量的分布及演化情况。数值模拟结果有助于理解MagLIF构型从预加热经由燃料压缩到最终发生聚变这一快速而复杂的过程,从而为建立相应的物理图像和认知提供了重要支撑,与传统ICF典型参数的对比也体现了该构型的优势所在,为后续研究奠定了基础。     

磁化套筒惯性聚变初步探索研究与实验可行性分析

磁化套筒惯性聚变(MagLIF)构型可充分利用现有大型脉冲功率驱动装置,如聚龙一号等。基于磁流体力学方程组和1∶1比例氘氚(DT)混合燃料聚变模型,开发了零维MagLIF数值模拟程序并进行了初步探索研究。计算结果表明初始负载参数(如轴向磁场强度,预加热温度、时刻,负载半径等)与聚变产额之间有着密切的联系,在给定条件下,可依据计算给出的定性关系进行负载优化设计。值得注意的是,根据计算结果,即使在理想条件下,氘氚燃料要实现能量收支平衡,则驱动器的电流必须大于21.2 MA。这意味着聚龙一号装置(10 MA)无法开展集成化的MagLIF实验,进一步的校验计算验证了上述观点,并在此基础上提出铝套筒分解实验的建议和负载设计参数。所取得的计算结果有利于加深对MagLIF套筒压缩阶段物理过程的认知和理解。     

磁化套筒惯性聚变中端面损失效应的一维唯象模型与影响分析

得益于激光预加热和轴向磁场的作用,磁化套筒惯性聚变(magnetized liner inertial fusion,MagLIF)构型理论上能有效降低聚变实现的难度,具有极大的应用潜力.本文选择MagLIF过程中伴随激光预加热所必然存在端面损失效应作为研究目标,搭建了能够描述几何参数与腊肠不稳定性等高维效应的一维唯象物理模型,并分别通过与二维流体动力学程序和国外同类程序的计算对比完成参数拟合校验;在此基础上,获得端面损失效应对MagLIF内爆过程及预加热效果的影响规律.计算结果表明:不同喷射半径下MagLIF负载在内爆过程的绝大多数时间内保持了相近的流体动力学演化过程,并在迟滞阶段经历了相同的质量损失比例,且考虑端面效应后得到的预加热和内爆产额相对变差,但却不改变规律性的趋势.所建立的模型与结论有助于加深对MagLIF预加热和端面损失过程中物理图像的认知和理解.     

磁化套筒惯性聚变中轴向磁场演化特征与Nernst效应影响

轴向磁场是磁化套筒惯性聚变(magnetized liner inertial fusion,MagLIF)有别于其他惯性约束聚变构型的主要标志之一.本文在建立集成化物理模型并编写一维模拟程序的基础上,通过对ZR装置驱动能力下典型MagLIF负载参数的模拟,系统研究并获得MagLIF各个阶段轴向磁场演化与分布特征,发现...     

磁化套筒惯性聚变初步探索研究与实验可行性分析北大核心CSCD

磁化套筒惯性聚变(MagLIF)构型可充分利用现有大型脉冲功率驱动装置,如聚龙一号等。基于磁流体力学方程组和1∶1比例氘氚(DT)混合燃料聚变模型,开发了零维MagLIF数值模拟程序并进行了初步探索研究。计算结果表明初始负载参数(如轴向磁场强度,预加热温度、时刻,负载半径等)与聚变产额之间有着密切的联系,在给定条件下,可依据计算给出的定性关系进行负载优化设计。值得注意的是,根据计算结果,即使在理想条件下,氘氚燃料要实现能量收支平衡,则驱动器的电流必须大于21.2 MA。这意味着聚龙一号装置(10 MA)无法开展集成化的MagLIF实验,进一步的校验计算验证了上述观点,并在此基础上提出铝套筒分解实验的建议和负载设计参数。所取得的计算结果有利于加深对MagLIF套筒压缩阶段物理过程的认知和理解。     

预加热效果对磁化套筒惯性聚变放能影响的模拟研究

得益于预加热和轴向磁场的作用,理论上磁化套筒惯性聚变(MagLIF)构型能有效降低聚变实现的难度,具有极大的应用潜力。然而,当前研究过于重视对激光能量沉积效率的提升,而忽略了预加热自身参数对MagLIF过程和内爆结果的影响。为此,采用一维集成化数值模拟程序MIST,开展了MagLIF过程中预加热效果对聚变放能影响的模拟研究,基于参数扫描方法,从简单的模型着手,逐步深入探讨相关参数对内爆结果的影响。模拟结果表明:预加热是MagLIF构型能够成功的必要条件,最佳时间是套筒即将开始向内压缩燃料的时刻;燃料预加热的设计原则是让燃料获得尽可能平缓分布的高温,而中心局部加热方式对于未能达到点火条件的负载更有优势;激光预加热模式下,脉宽越短越好,对于以ZR装置驱动能力为目标的算例而言,最佳套筒高度为1.0 cm。研究结果有助于加深对MagLIF过程中预加热机制和效果的认知和理解,对于具体的负载参数设计也有较强的指导意义。     

重离子惯性聚变的研究

本文描述了重离子惯性聚变的研究现状和用于重离子惯性聚变的重离子加速器． The present researching status and the heavy ion accelerators used for inertial fusion of heavy ions are decribed in the paper.     

High-gain magnetized inertial fusion

Magnetized inertial fusion (MIF) could substantially ease the difficulty of reaching plasma conditions required for significant fusion yields, but it has been widely accepted that the gain is not sufficient for fusion energy. Numerical simulations are presented showing that high-gain MIF is possible in cylindrical liner implosions based on the MagLIF concept [S. A. Slutz et al Phys. Plasmas 17, 056303 (2010)] with the addition of a cryogenic layer of deuterium-tritium (DT). These simulations show that a burn wave propagates radially from the magnetized hot spot into the surrounding much denser cold DT given sufficient hot-spot areal density. For a drive current of 60 MA the simulated gain exceeds 100, which is more than adequate for fusion energy applications. The simulated gain exceeds 1000 for a drive current of 70 MA.     

激光间接驱动惯性约束聚变二维总体程序——LARED集成程序

介绍了辐射流体力学程序LARED集成程序的物理背景、模型方程、数值方法和数值算例。该程序主要应用于激光间接驱动惯性约束聚变的二维整体模拟,兼顾激光直接驱动、辐射驱动靶丸内爆过程和流体不稳定性等物理过程的数值模拟。通过与实验数据、一维辐射流体力学程序进行比对,验证了程序的可靠性。该程序实现了多群输运建模下NIF点火靶的全过程数值模拟,并已应用于惯性约束聚变的物理研究。     

激光间接驱动惯性约束聚变二维总体程序——LARED集成程序

介绍了辐射流体力学程序LARED集成程序的物理背景、模型方程、数值方法和数值算例。该程序主要应用于激光间接驱动惯性约束聚变的二维整体模拟,兼顾激光直接驱动、辐射驱动靶丸内爆过程和流体不稳定性等物理过程的数值模拟。通过与实验数据、一维辐射流体力学程序进行比对,验证了程序的可靠性。该程序实现了多群输运建模下NIF点火靶的全过程数值模拟,并已应用于惯性约束聚变的物理研究。     

ICF内爆靶丸能流分布三维数值仿真

针对惯性约束聚变(ICF)两区均匀模型中靶丸表面能流分布问题,采用三维离散视角因子方法进行模拟仿真。对辐照过程中出现腔壁自身遮挡的面元,利用四叉树结构进行快速剔除;实现了GPU上的靶丸能流分布的并行计算。基于该方法开发的三维仿真软件IRAD3D 1.0,其有效性通过与解析解的对比得到了验证。不同算例模型的仿真结果表明了软件具有通用性、高效性,给出了靶丸均匀性分析的应用实例。     

ICF内爆靶丸能流分布三维数值仿真

针对惯性约束聚变(ICF)两区均匀模型中靶丸表面能流分布问题,采用三维离散视角因子方法进行模拟仿真。对辐照过程中出现腔壁自身遮挡的面元,利用四叉树结构进行快速剔除;实现了GPU上的靶丸能流分布的并行计算。基于该方法开发的三维仿真软件IRAD3D 1.0,其有效性通过与解析解的对比得到了验证。不同算例模型的仿真结果表明了软件具有通用性、高效性,给出了靶丸均匀性分析的应用实例。     

Numerical simulation of Ge solar cells using D-AMPS-1D code

A solar cell is a solid state device that converts the energy of sunlight directly into electricity by the photovoltaic effect. When light with photon energies greater than the band gap is absorbed by a semiconductor material, free electrons and free holes are generated by optical excitation in the material. The main characteristic of a photovoltaic device is the presence of internal electric field able to separate the free electrons and holes so they can pass out of the material to the external circuit before they recombine. Numerical simulation of photovoltaic devices plays a crucial role in their design, performance prediction, and comprehension of the fundamental phenomena ruling their operation. The electrical transport and the optical behavior of the solar cells discussed in this work were studied with the simulation code D-AMPS-1D. This software is an updated version of the one-dimensional (1D) simulation program Analysis of Microelectronic and Photonic Devices (AMPS) that was initially developed at The Penn State University, USA. Structures such as homojunctions, heterojunctions, multijunctions, etc., resulting from stacking layers of different materials can be studied by appropriately selecting characteristic parameters. In this work, examples of cells simulation made with D-AMPS-1D are shown. Particularly, results of Ge photovoltaic devices are presented. The role of the InGaP buffer on the device was studied. Moreover, a comparison of the simulated electrical parameters with experimental results was performed.     

中子-中子碰撞非线性输运的确定论模拟

从几个方面着手提高确定论方法的计算精度：首先,中子输运计算的相空间离散采用间断有限元处理方法,并使用较大的角度离散数和散射阶数;其次,使用蒙特卡罗直接统计方法得到高精度多群截面;最后,引入收敛于真解的中子-中子碰撞源迭代。数值算例验证表明,经过改进的确定论方法具有良好的稳定性和精度,能以可靠的精度求解考虑中子-中子碰撞过程的非线性问题。