基于脉管制冷机的二维可压缩流动数值计算程序的开发与应用

本文在SIMPLEC计算程序的基础上,开发了适用于各种流动Mach的计算程序,既能用于计算不可压缩流体,又可以计算可压缩流动。利用此程序计算了一些典型算例,和有关文献结果符合较好。最后利用可压缩程序计算了小孔型脉管制冷机当小孔处于最佳开度时小孔处的气流分布。     

基于FLUKA计算压力容器快中子注量的可行性研究

基于知识产权的考虑,通过与蒙特卡罗程序MCNP计算结果比对,研究使用FLUKA程序替代MCNP程序进行反应堆压力容器快中子注量计算的可行性。通过修改和调用子程序对次级粒子堆栈进行操作,解决了关闭裂变中子这一关键问题,FLUKA程序的计算结果与MCNP程序的计算结果相对偏差在5%以内,符合得较好,证明使用FLUKA程序替代MCNP程序用于计算反应堆压力容器快中子注量在技术上是可行的。     

封装结构强电磁脉冲多物理效应并行计算程序研制

介绍了自主研发的强电磁脉冲多物理效应并行计算程序JEMS-CDS-System的情况,该程序采用时域有限元方法,基于JAUMIN并行自适应结构网格支撑框架研制,并行效能高,可扩展性强,且支持动态负载平衡。通过算例测试表明,该程序对于键合线的电-热-应力失效过程的最高温度与范式等效应力计算结果与COMSOL软件计算结果吻合较好;SiP功率放大模块的热-应力耦合天河2高性能计算平台并行计算结果表明,该程序在CPU1024核时,具有38.1%并行效率。     

用P TA量表评价学生计算题解题能力

在教学过程中, 如何客观准确地评价学生的计算题解题能力? 笔者分析构成计算题解题能力的基本 要素以及确定这些要素的评价指标, 分配指标的权重, 细化评判的标准, 编制了P TA量表. 综合多次P TA量表的得 分对学生的计算题解题能力做出客观的评价, 利用客观的评价结果与家长, 学生交流. 制定相应的策略去提高学生 的计算题解题能力, 取得了不错的效果     

基于认知规律的牵连速度进阶式教学设计

在教学过程中, 如何客观准确地评价学生的计算题解题能力? 笔者分析构成计算题解题能力的基本 要素以及确定这些要素的评价指标, 分配指标的权重, 细化评判的标准, 编制了P TA量表. 综合多次P TA量表的得 分对学生的计算题解题能力做出客观的评价, 利用客观的评价结果与家长, 学生交流. 制定相应的策略去提高学生 的计算题解题能力, 取得了不错的效果     

JMCT的高浓铀模型验证计算北大核心CSCD

为了对JMCT程序的准确性和适用性进行验证,从国际临界基准库中选取了56个模型,通过使用JLAMT程序进行建模,然后利用JMCT程序计算了模型的有效增殖因子,并与基准值、MCNP程序计算结果进行对比。计算结果表明,JMCT计算值基本都在基准值的3倍标准差以内,且JMCT计算结果和基准值的相对偏差基本都在2%以下。在一定程度上验证了JMCT程序的准确性和适用性,为后续JMCT程序的使用验证提供借鉴。     

涡轮叶栅气热耦合数值模拟

本文建立了透平机械叶栅流场与叶片温度场的气热耦合计算模型,编制了三维非正交曲线坐标系下的温度场求解程序,计算了实心及空心叶片的温度场分布。介绍了区域分解方法在耦合计算及并行计算中的应用,以矩形区域温度场求解程序为例研究了程序的并行化,最后进仃了叶栅流场计算与叶片导热计算程序的分区耦合计算。     

激光烧蚀瑞利－泰勒不稳定性模拟

 给出了激光烧蚀流体不稳定性计算程序EUL2D的物理方程,介绍了计算中使用的活动网格和一些技术问题处理。EUL2D程序的计算结果与Takabe公式、FAST2D程序和LASNEX程序的结果较好符合。数值计算日本大阪大学激光烧蚀瑞利－泰勒不稳定性实验,再现了实验结果。发现了横向电子热传导烧蚀在长波长扰动的非线性瑞利－泰勒不稳定性演变中起重要作用。     

闪光X射线照相蒙特卡罗程序的并行化

研究了闪光X射线辐射照像蒙特卡罗程序(FXRMC)在MPI平台下的并行计算实现,给出了实现过程中并行随机数的产生方法。对并行程序的测试结果表明:并行程序与串行程序结果一致,加速比比较理想,呈线性增长,并行效率在16个处理器上可达80%以上。算例的结果说明了并行化可有效地解决程序计算散射技术性能时的耗时问题,从而有效化解FXRMC耗时和大规模计算的困难,提高了FXRMC程序的计算规模和计算速度,达到了研究要求。 (Institute of Fluid Physics, CAEP, P. O. Box 919-105, Mianyang 621900, China)     

有限元方法在二维六角形组件堆芯临界问题计算中的应用

应用有限元方法对二维六角组件维芯的少群中子扩散临界计算进行了研究，并编制了二维少群中子扩散方程有限元计算程序２ＤＦＭＥ。以快堆３００ＭＷ－ＬＭＦＢＲ临界问题为例，作了数值计算。且与国内外有限元程序和有限差分程序ＣＩＴＡＴＩＯＮ的计算结果进行了分析比较。     

泄压口大小对约束空间爆炸准静态超压载荷的影响规律

约束空间爆炸载荷主要包含瞬态冲击波和持续时间较长的准静态超压,为了研究泄压口大小对约束空间爆炸准静态超压载荷的影响规律,采用三阶WENO有限差分格式,编写了约束空间内炸药爆炸冲击波三维数值计算程序。利用Sod激波管、双爆轰波碰撞、激波与熵波相互作用等经典算例验证了数值程序的可靠性。基于验证的程序开展了带有泄压口的约束空间内爆炸载荷数值计算。在分析数值计算结果基础上,提出了描述准静态超压载荷简化理论公式,该理论公式与数值计算结果吻合较好。研究结果可为工程抗爆结构设计提供一定的参考和指导。     

JMCT程序临界安全基准校验计算与分析

JMCT程序是北京计算数学与应用物理研究所自主开发的蒙特卡罗方法粒子输运计算程序。从国际临界安全实验数据库中选取棒栅类、U系、Pu系、233U系实验方案对JMCT程序进行验证,验证涵盖了常见裂变元素（如U,Pu等）从低富集度到高富集度的快、热能区,同时也包含了常见的中子毒物以及反射层材料。将JMCT程序的计算结果与基准实验值进行对比,并且与MCNP,MONK程序的计算结果进行了比较。结果表明,在检验的范围内,JMCT程序具有与国际通用蒙特卡罗方法粒子输运程序一样良好的计算精度。     

短波自由电子激光放大器数值模拟

介绍在理论上研究自由电子激光(FEL)的MICFEL程序。MICEFL是以单粒子理论为基础的二维放大器程序。利用MICFEL计算的Paladin实验结果与FRED程序计算结果以及25米Paladin实验结果相比均基本符合。理论计算的结果显示了FEL增益光导的特性。     

特殊炸药的爆轰参数计算

 针对BKW（Becker-Kistiakowsky-Wilson）状态方程,编写了爆炸参数计算程序。该程序除了可以计算常规含C、H、O、N元素炸药的爆轰参数,还能计算含有其它金属元素炸药的爆轰参数。通过对含Fe、Mn元素的乳化炸药以及含Al元素的水胶炸药的爆轰参数计算表明,该计算程序得到的计算结果与实验结果基本一致,相对误差不超过1%。采用该程序计算爆轰参数时,只需输入炸药的化学式、常温常压密度和生成焓即可。     

基于非线性迭代的压水堆瞬态计算程序开发

采用目前工业上成熟的非线性迭代计算策略,基于两群粗网有限差分方法和多群UNM节块方法,开发了针对压水堆工况的三维瞬态扩散计算程序。UNM方法采用解析基函数作为基函数,通过方程变换解决了解析节块法在临界节块计算不稳定的问题,提高了计算精度。热工计算采用单通道模型和燃料棒一维导热模型,相比目前堆芯在用的经验关系式方法,该模型可以更加准确地计算燃料棒温度分布。采用基于横向积分方程的三节块方法,可以有效减轻控制棒尖齿效应对瞬态计算的影响。为测试程序性能,采用NEACRP等基准算例对程序进行了校验。数值结果表明,开发的程序计算结果正确,适用于压水堆堆芯瞬态过程的模拟。     

三维多物质欧拉界面处理的并行算法研究

针对现有的算法无法对爆炸冲击进行有效的研究这一问题,在改进了Youngs界面技术的基础上设计了三维多物质Euler界面处理的并行算法。然后采用MPI标准进行了算法的程序设计并对程序进行了并行性能测试和不同分区的算例测试,测试结果表明并行算法对加速比有较大的提高,而不同分区对并行算法的计算结果并没有影响。最后应用所编写的程序对空中爆炸和聚能射流进行了仿真,结果表明程序的模拟结果与前人的研究数据是符合的,说明了程序的有效性。     

DFVLR压气机转子中跨声速流动的两类流面准三元迭代解

本文依据三元流动的通用理论,完成了一个跨声速压气机转子流场的准三元迭代计算。计算表明,使用S_1跨声速正问题程序和S_2反问题程序的迭代方法是十分有效的。计算中一系列合理措施的采用,解决了迭代的收敛问题,收到了颇为满意的效果。文中的算例是西德具有详细内部流动测量数据的DFVLR跨声速转子的设计点工况,计算结果与实验值进行了比较。     

中国实验快堆燃料组件精细功率分布计算

在热工中,需要中子学计算给出燃料组件内各元件棒功率的相对分布。利用蒙特卡罗程序对中国实验快堆（CEFR）燃料组件内元件棒功率分布进行了理论计算分析,并保证计算结果相对统计误差小于0.8%。使用另一个基于六角形节块扩散理论的钠冷快堆中子学设计软件NAS计算得到的结果对蒙特卡罗程序计算结果进行了对比计算。结果表明,蒙特卡罗程序与NAS计算得到的元件棒相对功率分布结果的最大相对偏差小于3%。使用蒙特卡罗程序对CEFR燃料组件内精细功率分布的计算是可靠的,可用于设计计算当中。     

改进的线性子链方法在中子活化计算中的应用

在传统线性子链方法（TTA）采用的Bateman表达式基础上,利用极限运算推导出允许重复特征值出现的广义解析表达式,克服了重复核素对线性子链方法计算稳定性的限制,并通过回溯算法的应用实现了对核反应链的自动搜索。通过以上对线性子链方法的改进,显著提高了计算效率,最终形成了用于中子活化计算的程序ITACT。基于ITACT程序和EAF-2007数据库,计算了压水堆燃料包壳材料和聚变堆第一壁材料活化等例题,并将计算结果与欧洲活化程序FISPACT进行对比。结果表明:对于长寿命核素,两者的计算结果较为接近,对于短寿命核素,ITACT程序的计算精度更高,初步验证了ITACT程序的可用性和准确性。     

J-TEXT 装置等离子体平衡重建与集成运算

利用Matlab 语言将EFIT 程序反演计算过程进行优化整合,简化了程序操作步骤,并增加了对计算结果的分析与存储功能。利用整合后的EFIT 程序,计算分析了J-TXET 装置放电过程中等离子体磁场位形和其它参数的演变,绘制了等离子体压强剖面图、等离子体电流密度和安全因子分布等磁面信息的重建结果。此外,还将反演结果与软X 射线测算数据进行比对。结果表明,改善后的EFIT 程序能够初步满足对实验放电的反演计算分析要求,为J-TXET 装置的放电控制和运行提供重要的参考数据。