二维磁驱动数值模拟程序MDSC2的验证与确认

为了对磁驱动实验提供高置信度的数值模拟,需要开展磁流体力学程序的验证与确认。采用人为解比较法、网格收敛性研究和与成熟程序比较等方法,对二维磁驱动数值模拟程序MDSC2进行了程序验证。数值模拟表明：MDSC2程序正确地表示了磁流体力学模型,其中热扩散、磁扩散的离散格式具有二阶收敛精度。采用与磁驱动实验相比较的方法,进行了MDSC2程序的确认。对聚龙一号装置上的PTS-061发次磁驱动单侧飞片发射和PTS-122发次磁驱动双侧飞片发射实验进行了模拟,模拟的飞片自由面速度与实验测量的飞片自由面速度相一致;对FP-1装置上的固体套筒实验进行了模拟,模拟的套筒内外半径与实验测量结果相一致。MDSC2程序能正确模拟磁驱动单侧飞片发射、磁驱动双侧飞片发射和磁驱动固体套筒等磁驱动实验。     

电磁线圈炮计算程序的收敛性与实验验证

根据电磁线圈炮中推力线圈与电枢之间的感应耦合方程和电枢动态响应方程,考虑发射过程中推力线圈和电枢的欧姆加热,开发了电磁线圈炮的计算程序,给出了程序的算法和流程。分析了程序的收敛性,认为最大时间步长和电枢最大网格均依赖于驱动脉冲波形的上升时间。通过模拟结果与单级线圈炮实验测量结果的比较,验证了程序的有效性。模拟的线圈电流波形与实验波形吻合较好,模拟的出口速度比测量速度大约6.5%。程序收敛性和有效性的验证表明该程序可用于电磁线圈炮系统的初步设计。     

组件计算程序PANDA研发及初步验证

PANDA是上海核工程研究设计院研发的压水堆组件计算程序。该程序采用基于特征线方法（MOC）的一步化计算流程,即在不引入能群压缩和栅元均匀化的情况下直接进行组件层面的两维非均匀输运计算。多群数据库采用基于ENDF/B-VI制作的70群结构中子数据库,基于ENDF/B-VII的新版数据库也正在开发中。共振自屏计算采用了空间相关丹可夫方法（SDDM）,既具备燃料芯块分区计算的能力,又保留了传统Stammler方法的计算效率。多群非均匀输运计算采用二维模块化MOC方法,并辅以双重粗网有限差分(CMFD)加速技术,具有良好的计算精度和效率。对传统线性子链解析（TTA）方法以及多种矩阵指数方法进行了研究,选取了适合PANDA程序燃耗链的燃耗方程求解技术。基于以上基本模型开发了PANDA程序,并从程序模块、总体集成和核设计程序系统确认等三个层面,初步验证了PANDA程序的计算性能,表明了PANDA程序的工程设计计算能力。     

组件计算程序PANDA研发及初步验证

PANDA是上海核工程研究设计院研发的压水堆组件计算程序。该程序采用基于特征线方法（MOC）的一步化计算流程,即在不引入能群压缩和栅元均匀化的情况下直接进行组件层面的两维非均匀输运计算。多群数据库采用基于ENDF/B-VI制作的70群结构中子数据库,基于ENDF/B-VII的新版数据库也正在开发中。共振自屏计算采用了空间相关丹可夫方法（SDDM）,既具备燃料芯块分区计算的能力,又保留了传统Stammler方法的计算效率。多群非均匀输运计算采用二维模块化MOC方法,并辅以双重粗网有限差分(CMFD)加速技术,具有良好的计算精度和效率。对传统线性子链解析（TTA）方法以及多种矩阵指数方法进行了研究,选取了适合PANDA程序燃耗链的燃耗方程求解技术。基于以上基本模型开发了PANDA程序,并从程序模块、总体集成和核设计程序系统确认等三个层面,初步验证了PANDA程序的计算性能,表明了PANDA程序的工程设计计算能力。     

JMCT在CPR1000工程设计中的验证与确认

CPR1000系列反应堆是目前国内广泛应用的第二代压水堆型号之一,蒙特卡罗程序在CPR1000系列反应堆的验证与确认是该程序实现反应堆工程设计应用的关键环节。基于某CPR1000机组实际参数,使用由国内单位研发的蒙特卡罗程序JMCT在该机组开展了粒子输运建模计算,分别进行了临界计算和固定源计算,并进行了验证与确认。对于临界计算,采用JMCT建立了全堆芯pin-by-pin模型,计算了堆芯有效增殖因子和功率分布。对于固定源计算,建立适用于屏蔽分析的反应堆模型和辐照监督管精细结构模型,计算了两个核电机组多个循环的辐照监督管探测器位置累积快中子注量。通过将JMCT的计算结果与参考程序的计算结果、反应堆实际测量值进行了对比,验证了JMCT程序在CPR1000反应堆工程设计中的实际使用效果,证明了JMCT程序具备工程级的计算精度。     

刚体与SPH耦合的三维数值方法研究

 研究了刚性体与SPH（Smoothed Particle Hydrodynamics）方法相耦合的计算方法和技术,以此技术为核心,在原有三维综合拉氏程序CL3D的基础上,研制了适用于斜侵彻问题数值研究的刚性弹体与SPH靶体耦合计算的拉氏三维数值模拟程序CRSPH3D。针对金属材料编入了Johnson-Cook本构模型计算模块,针对混凝土材料编入了Johnson-Holmquist本构模型及与之配套的损伤模型和状态方程等的计算模块。用研制的计算程序对弹体分别正碰和斜碰穿透中厚铝靶和侵彻厚混凝土靶的过程进行了数值模拟计算和分析。结果表明,对于铝靶穿透问题,模拟计算得到的弹体剩余速度、弹体动能损失和贯穿图像与实验结果符合得较好。对于混凝土靶的侵彻问题也得到了与实验现象一致的、定性合理的图像。     

量子算法核磁共振实现脉冲序列设计程序问题研究

从两量子位核磁共振量子处理器物理模型出发,利用Raedt小组提出的自旋-1/2代数理论,根据量子控制非门的定义及Grover量子算法原理,介绍了量子控制非门的4种不同脉冲序列及两量子位Grover量子算法的两种不同脉冲序列的设计过程,通过数值求解含时薛定谔方程模拟量子控制非门和两量子位Grover量子算法,等价于执行量...     

BEAVRS基准模型的cosRMC精细建模与临界测试验证

BEAVRS基准题是来自麻省理工学院(MIT)的基于商用压水堆全堆芯精细建模的国际基准模型,具有详细的实测数据,是由多种富集度燃料及控制棒组件组成的大型压水堆分析验证基准题。使用我国自主研发的蒙卡程序cosRMC对BEAVRS基准题进行精细建模,主要计算了热零功率(HZP)状态下的临界本征值、全堆功率分布和控制棒价值,并与实测值以及国际知名蒙卡程序MCNP,OpenMC,MC21等结果进行对比。HZP状态下,cosRMC临界本征值结果与MCNP计算结果相差仅7.1 pcm,符合较好;不同控制棒组件插入情况下的临界本征值与理论值1.000的偏差小于0.74%,控制棒价值结果与实测值误差小于100 pcm,计算精度与同类软件相当;此外还对比分析了全堆功率分布与实测值结果的吻合程度及误差产生原因。初步验证了cosRMC程序对复杂堆芯精细建模计算的可行性和准确性,为程序以后的应用及完善打下基础。     

BEAVRS基准模型的cosRMC精细建模与临界测试验证

BEAVRS基准题是来自麻省理工学院(MIT)的基于商用压水堆全堆芯精细建模的国际基准模型,具有详细的实测数据,是由多种富集度燃料及控制棒组件组成的大型压水堆分析验证基准题。使用我国自主研发的蒙卡程序cosRMC对BEAVRS基准题进行精细建模,主要计算了热零功率(HZP)状态下的临界本征值、全堆功率分布和控制棒价值,并与实测值以及国际知名蒙卡程序MCNP,OpenMC,MC21等结果进行对比。HZP状态下,cosRMC临界本征值结果与MCNP计算结果相差仅7.1 pcm,符合较好;不同控制棒组件插入情况下的临界本征值与理论值1.000的偏差小于0.74%,控制棒价值结果与实测值误差小于100 pcm,计算精度与同类软件相当;此外还对比分析了全堆功率分布与实测值结果的吻合程度及误差产生原因。初步验证了cosRMC程序对复杂堆芯精细建模计算的可行性和准确性,为程序以后的应用及完善打下基础。     

涡轮叶栅气热耦合数值方法与验证

应用三阶精度TVD格式和自由型曲面阿格技术,对涡轮内冷叶片进行气热耦合数值模拟.结合商用软件采用多种湍流模型计算方法的结果进行实验验证和对比分析,结果表明,商用软件采用传统湍流模型对温度边界层的模拟是不合理的,采用考虑转捩的湍流模型,调整转捩起始雷诺数,可以较好地模拟温度边界层的热传导.本文对非转捩区边界层温度的计算结果与实验结果吻合较好.     

二维磁流体力学程序MFCG Ⅱ与MFCG Ⅰ计算结果的对比分析

根据磁场、电流密度和感应电流的解析公式研制了爆炸磁通压缩发生器充电和爆磁压缩过程的计算程序MFCG Ⅱ。在相同条件下用该程序和二维磁流体力学程序MFCG Ⅰ分别计算相同模型所得结果符合较好。因此用MFCG Ⅱ可分析、解释MFCG实验中有关的物理问题。     

我们心目中的楷模

于敏先生是我国核科技事业的卓越指挥者，是为“两弹一星”研制做出突出贡献的科学家。我们认识他并在他指导下工作是在氢弹突破之后的20世纪60年代末，我跟着他参加过几个大项目的理论设计研究工作。在工作中，需要用计算机对设计方案进行物理验证，一个方案需经反复验证，多次修改才能方案定型．通常初步计算结果出来之后，搞物理的找出问题所在，迅速做出修改，并交由计算工作者修改程序，再上机，在理论设计比较成熟的基础上，通过冷实验加以验证，得出结果后再做物理分析，再上机计算，最终由热试验加以检验。这一过程有点像在涌浪波涛间行船，能够最终抵达彼岸需要经历许多困难和曲折。     

二维磁流体力学程序MFCG Ⅱ与MFCG Ⅰ计算结果的对比分析

 根据磁场、电流密度和感应电流的解析公式研制了爆炸磁通压缩发生器充电和爆磁压缩过程的计算程序MFCG Ⅱ。在相同条件下用该程序和二维磁流体力学程序MFCG Ⅰ分别计算相同模型所得结果符合较好。因此用MFCG Ⅱ可分析、解释MFCG实验中有关的物理问题。     

一个高压比离心叶轮的CFD结果确认

本文采用NUMECA的FINE/Turbo软件包数值分析了一个高压比离心叶轮的三维定常粘性流场。选用二阶精度的中心格式和Spalart-Allmaras的一方程湍流模型。通过与现有实验结果(包括性能曲线、外环壁静压和截面马赫数分布等)的比较对计算结果就行了详细的确认,同时讨论了网格密度和间隙尺度、以及不同的数据处理方法对计算结果确认的影响。     

LARED-H程序的激光黑腔靶耦合数值模拟

LARED-H程序是一个可用于激光黑腔靶耦合数值模拟研究的二维辐射流体力学程序。黑腔等离子体所形成的复杂流场使单纯的拉格朗日网格在计算中产生严重扭曲，影响计算精度，并导致计算中断。拉氏加网格重分是计算激光黑腔靶耦合常用的算法。对LARED-H程序的积分网格重分方法在网格重构和物理量重映方面作了较大改进，并利用改进后的LARED-H程序模拟了“神光”-Ⅱ和“神光”-Ⅲ条件下的激光空腔靶耦合物理全过程。     

金属圆管内爆轰波相互作用效应的数值模拟

炸药两端同时起爆时金属圆管运动规律的研究，粗略的实验不少，但在数值模拟方面，特别是在对碰部位，或是不能计算出鼓包现象，或是鼓包的计算结果与实验差别较大。文中采用动力学有限元程序，针对大变形带来的网格扭曲造成计算的锁死现象进行了ALE技术处理，对金属圆管内爆轰波的相互作用效应进行数值模拟，给出了爆轰波对碰作用下金属圆管运动规律的计算结果。     

高功率微波土壤击穿的数值验证研究

高功率微波在土壤中传播时, 会引起土壤击穿电离而导致土壤电阻率的非线性变化, 土壤电阻率的变化又将反作用于传播过程, 加剧高功率微波衰减, 影响其能量传输效率. 通过对土壤动态电离过程的分析, 结合Maxwell方程组构建了高功率微波土壤传播模型, 采用时域有限差分法对该模型进行数值验证. 数值结果显示了高功率微波在土壤中传播、衰减等物理图像, 以及土壤电阻率的非线性变化过程. 理论分析验证了这些数值结果.     

二维磁流体力学程序MFCGⅢ计算结果与测试结果的对比分析

 根据磁场、流函数和感应电流的解析公式，研制了间接馈电爆磁压缩发生器充电和爆磁压缩过程计算程序MFCGⅢ。用MFCGⅢ计算了实验模型，其结果与测试结果符合较好。因此MFCGⅢ可对实验提供优化设计，并且可分析、解释MFCG实验中有关的物理问题。     

Euler-Lagrange耦合计算的GEL方法

研究了一种Euler-Lagrange耦合数值方法——GEL（Ghost-fluid Euler-Lagrange）方法，编写了GEL二维计算程序。