“聚龙一号”装置上磁驱动铝飞片实验的数值模拟

采用二维磁驱动数值模拟程序MDSC2,对大电流脉冲装置"聚龙一号"上的151发次磁驱动370μm厚铝飞片实验、164发次磁驱动330μm厚铝飞片实验进行了数值模拟和分析。数值模拟表明:370μm厚铝飞片和330μm厚铝飞片的磁驱动过程中,VISAR测量的速度不是飞片自由面的速度,而是飞片中邻近自由面最近的固体反射面的速度。这是由于磁驱动飞片发射过程中,飞片自由面部分被烧蚀,密度低于固体密度状态,而飞片自由面和加载面中间的飞片还保持固体密度状态。VISAR测量的激光将穿过自由面的低于固体密度状态的飞片部分,到达飞片自由面最近的固体密度位置再反射回去,获得这一位置的飞片速度。数值模拟的飞片固体反射面速度历史和VISAR测量的速度历史相吻合。     

准等熵压缩实验中金属铝的热力学耗散分析

为获取高压下材料的纯热力学压力-比容参考线和完全物态方程,减去应力-应变曲线中的其它信息,对准等熵压缩实验中由加载应变率引起的黏性耗散和热传导引起的热耗散做了分析讨论。基于反积分计算和流体动力学积分计算相结合的方法,根据激光加载(约108 s-1)和磁驱动准等熵压缩(约105 s-1)的实验数据,对材料声速、应力-应变曲线、温度和熵增等物理量进行计算,分析了不同应变率与该物理量的关系;还对热传导和SCG本构模型进行了计算,分析了热传导引起的温度变化对材料屈服强度、剪切模量和拉格朗日声速的影响。结果表明:激光加载实验中,应变率引起的温升差异约为180K,熵增差异约为250J/(kg·K),热传导导致温度下降40K;磁驱动准等熵压缩应变率较低,引起的熵增变化小于8J/(kg·K)。     

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

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

磁驱动高速飞片实验样品设计

利用理论与数值模拟相结合的方法对磁驱动飞片实验样品的尺寸设计进行了分析探讨,给出了横 向尺寸设计的估算公式、长度设计原则、影响厚度设计的主要因素等。     

衬套内爆ALE方法二维MHD数值模拟

早期关于二维磁流体力学方程组的描述都是在Lagrangian或Eulerian坐标内进行的。在Lagrangian描述中，计算网格固定在物体上随物体一起运动，网格点与物质点在物体的变形过程中始终保持一致，因此能准确描述物体的移动界面。但对于大变形问题，物质的扭曲将导致计算网格的畸形而使计算无法进行下去。在Eulerian描述中，网格固定在空间中，因而计算网格在物体的变形过程中保持不变，因此很容易处理物质的扭曲。但对于运动界面，需要引入非常复杂的数学映射，将可能导致较大的误差。     

电磁驱动45钢准等熵压缩的实验研究

在已有的众多准等熵加载技术中,磁驱动准等熵加载技术具有准等熵程度高、压力范围大、实验材料种类多、效费比高等特点。利用中物院流体物理研究所建成的磁驱动准等熵压缩和高速飞片实验装置CQ-1.5(最高加载压力为50GPa),成功开展了45钢的准等熵压缩实验,对装置的主要参数进行了介绍;利用激光干涉测试系统DPS获得了45钢飞片的自由面速度历史,通过反积分处理给出了材料准等熵压缩的p-V关系。通过分析实验数据,获得了45钢3种形式的等熵方程的参数。实验获得的最高等熵压力为47.5GPa。     

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

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

微放电等离子体多负辉区结构融合过程数值模拟研究

高气压下的微型电热推进器(MPT)中的放电等离子体存在多负辉区结构，其负辉区有融合趋势。对矩形微放电等离子体推进器(RMPT)的负辉区融合过程进行了二维模拟分析，在方法上采用了非平衡态的自洽流体模型，并考虑了离子电流加热和三体碰撞过程。结果显示:矩形微放电等离子体推进器(RMPT)在低电流条件下存在两个稳定的负辉区，当超过某一电流阈值条件后，两个负辉区会在腔体中心重合。分析了这一过程的成因，认为其融合过程本质上是空心阴极的导通过程，其融合与否与鞘层电压有关。     

自由面被烧蚀磁驱动飞片的数值模拟

为了理解磁驱动飞片自由面烧蚀的机理、自由面被烧蚀后飞片的物质状态、自由面被烧蚀后激光速度干涉仪(VISAR)测量的机理等,采用磁驱动数值模拟程序MDSC2对聚龙一号装置上PTS-151发次磁驱动飞片实验中370 m厚飞片进行了模拟和分析。数值模拟表明,334 ns之前飞片自由面部分保持固体状态,334 ns之后飞片自由面部分已经熔化,到340 ns后整个飞片都被熔化。飞片自由面烧蚀的主要机制是电流焦耳加热,热扩散和压缩做功的贡献很小。数值模拟的固体密度反射面速度历史和VISAR测量的速度历史一致。飞片自由面熔化后,VISAR测量的速度是距离自由面最近的固体密度反射面的速度。     

拉氏方法模拟二维多介质可压缩流体的运动

在用拉格朗日方法模拟二维多介质可压缩流体的运动时,网格发生大变形往往是模拟不能正常进行下去的重要原因。网格动态局域重分可以有效改善网格的畸变程度,使计算得以持续。针对三角形计算网格提出了一种新的动态局域重分方法,包含"对角线交换""长边劈裂""短边融合"和"帽子戏法"4种基本操作,其中前3种操作不仅作用于同种介质内部,还可将其拓展到多介质界面处,与"帽子戏法"一起处理界面附近的大变形网格。在网格动态局域重分后,将旧网格上的物理量映射到新网格上,先计算出新三角形的质量和内能,再根据动量守恒和能量守恒对新三角形的格点速度及内能进行修正。利用该方法对冲击波与气泡相互作用和R-T不稳定性问题进行了数值模拟,取得了良好的效果。在R-T不稳定性算例中,采用同种介质和不同介质两种模型进行对比,模拟结果验证了该方法的有效性。