电磁内爆X射线数值计算

叙述一维辐射磁流体力学方程组（三温）的数值模拟，用于电磁内爆X射线计算。动量方程、三温方程组采用分裂法求解，以克服强耦合非线性可能引起的数值不稳定性。采用隐式差分格式避免小时间步长限制。隐式格式采用Newton-Raphson法迭代求解。使用了分层网格、加权平均，迭代初值公式，自动调整时间步长，奇点处理等特殊处理方法。最后给出Pegasusl装置一个实验模型的计算结果，获得15J的X射线能量，峰值功率为16GW。