节块内嵌离散纵标(SN)方法的算法研究及程序开发

介绍节块内嵌离散纵标(S_N)方法求解三维堆芯中子输运/扩散方程的算法框架.在基于扩散理论的三维粗网节块展开方法(NEM)的算法体系中,用基于输运理论的径向二维细网节块离散纵标方法(NDOM)的内迭代过程,替代节块展开方法(NEM)内迭代的径向求解过程.该算法充分考虑了核电厂反应堆堆芯的三维结构特点,另一方面,也充分利用了已经成熟的三维粗网节块展开方法(NEM)和二维离散纵标方法(S_N)的研究成果,同时有效避免了利用离散纵标方法(S_N)求解三维中子输运方程所面临的计算内存和计算时间的瓶颈问题.编制开发二维多群节块离散纵标方法(NDOM)模块程序NS_NM和三维多群节块展开方法(NEM)模块程序MGNEM,并以此为基础编制开发节块内嵌S_N方法的模块程序HANWIND;其中,NS_NM为HANWIND求解两维问题的功能模块.针对OECD/NEA-2D C5G7MOX基准问题以及两环路核电厂三维堆芯的数值验算结果表明,节块内嵌S_N方法的算法开发及程序编制有效、切实可行.     

基于节块展开法的Jacobian-Free Newton Krylov联立求解物理-热工耦合问题

目前反应堆物理热工耦合程序通常采用固定点迭代思路, 这可能导致部分工况收敛速度慢, 甚至出现不收敛的现象, 严重影响了计算效率. 基于此, 本文将高效的粗网节块展开法(NEM)与Jacobian-Free Newton-Krylov (JFNK)方法结合, 成功地开发出了一套新方法NEM_JFNK, 实现了联立求解物理热工耦合问题. 首先将NEM推广到热工问题的求解, 之后使用NEM来离散物理-热工耦合问题的所有控制方程, 使得所有变量都能在粗网格下进行离散, 从而大大减小求解问题的规模; 其次将NEM离散后的方程经过某些特殊的处理, 成功地嵌入JFNK的计算框架, 最终开发出了基于线性预处理的NEM_JFNK, 即LP_NEM_JFNK. 此外, 为了充分利用原有的迭代程序, 避免JFNK残差方程的重新建立, 本文还开发了无需重构残差方程的NEM_JFNK, 即NRC_NEM_JFNK, 并实现“黑箱”耦合. 文中以一维中子-热工模型为例, 给出LP_NEM_JFNK和NRC_NEM_JFNK数学模型, 并对计算结果进行分析. 结果表明:新方法无论是收敛速度还是计算效率都具有明显优势.     

基于平面通量展开的三维Pin by Pin输运SP3计算

基于Pin by Pin输运SP₃计算是下一代反应堆物理计算方法中一个重要的研究方向。节块法一般采用高阶通量和中子源展开以适应较大的节块尺寸，因此对于Pin by Pin计算来说效率偏低。由于反应堆堆芯不均匀性更多发生在径向，本文提出一种径向基于二维平面通量展开结合轴向常规节块法的综合方法进行三维Pin by Pin输运SP₃计算，其三个方向的节块耦合迭代均采用基于中子净流耦合的格式。通过IAEA 2D/3D基准问题和典型输运基准问题计算，验证该计算方法对于堆芯扩散计算和输运计算的正确性和可行性。     

改进的节块展开方法求解圆柱几何对流扩散方程

为高效求解球床高温气冷堆物理-热工耦合问题,发展改进节块展开法求解圆柱几何下的对流扩散方程.针对圆柱几何和对流扩散方程的特殊性,采用三阶多项式和指数函数作为r向横向积分方程的展开函数,在节块展开法的框架下高效求解对流扩散方程.数值验证表明,改进的节块展开方法具有固有的迎风特性,在使用粗网节块时依然能保持稳定性和较高的计算精度.     

基于二维广义横向积分综合一维半解析的三维中子扩散节块法

目前, 轻水堆堆芯计算广泛采用基于横向积分技术的中子扩散方程节块展开法, 该方法需要对横向泄漏进行多项式近似而不严格, 而且堆芯设计还需要额外的精细功率重构模块用于获得组件内各栅元的功率分布。本文提出两维广义横向积分方法, 直接采用源展开以及表面流耦合方法, 可以避免上述两个不足。由于反应堆堆芯不均匀性更多发生在径向, 因此采用径向基于广义横向积分方法结合轴向常规节块法的综合方法进行三维中子扩散计算。通过基准问题的数值计算, 验证了该方法对于堆芯扩散计算的正确性和可行性。     

全堆芯Pin by Pin中子扩散节块法CMFD加速

采用节块方法中有效的半解析方法求解中子扩散方程，在径向采用pin by pin的计算方式且中子通量基于两阶展开，在轴向采用基于4阶展开中子通量的节块法，并且进行全堆粗网有限差分（CMFD）加速，以达到计算精度与计算效率的平衡.通过IAEA 2D/3D基准问题和自定义3D pin by pin参考问题的数值计算，验证所提出的方法能够达到预期的精度，同时CMFD能够取得很好的加速效果.     

改进准静态方法的研究

处理反应堆三维动力学问题时，在时间方面采用了改进准静态方法，并加入了综合法改进，它通过对中子通量的因式分解，达到放大时间步长的目的，减少进行空间计算的次数，在必要的空间计算中，该方法又利用了格林函数节块法的原理求解形状函数的空间分布，使空间网格距可以达到有限差分法的近二十倍。     

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

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

三维圆柱几何格林函数节块法中子扩散计算

发展了中子扩散计算三维圆柱几何格林函数节块法。首先通过横向积分将中子扩散方程化为三个互相耦合的一维偏通量方程。对于径向偏通量方程，将径向扩散微分算符分解为平板几何的扩散微分算符和一个修正项之和，将修正项移到方程右端作为修正源项。这样，三个方程都化为平板几何的一维方程莆式，再借助平板几何第二类边界条件格林函数，对主几相应体源作积分，建立偏通量积分方程，对于修正源项，通过分部积分方法将偏通量导数项转化     

变系数(2+1)维Nizhnik-Novikov-Vesselov方程的三孤子新解

本文为获得非线性发展方程的相互作用解,研究了辅助方程法,并扩展应用辅助方程法和(G'/G)展开法, 获得了变系数非线性(2+1)维Nizhnik-Novikov-Vesselov方程的由椭圆函数、双曲函数、 三角函数和有理函数混合构成的新相互作用解. 关键词： G'/G)展开法')" href="#">(G'/G)展开法 辅助方程法 三孤子解     

积分节点迎风偏移的节点积分无单元Galerkin方法

建立求解稳态对流-扩散方程的-种稳定、高效的无单元Galerkin方法.该方法计算积分时采用基于局部Taylor展开的节点积分,并根据对流占优的程度对积分节点进行自适应迎风偏移.与传统的使用稳定化的无单元Galerkin方法相比,该方法是-种不依赖于背景网格积分的纯无网格方法,具有更好的稳定性和较高的计算效率,其程序实施更为简便.     

Modified nodal point method for the measurement of the focal length of divergent lenses

The well-known nodal point method has been employed with some modification to measure the focal lengths of divergent lenses. In this modification, a telescope was used to locate precisely the nodal/principal point of the lens under test, and a travelling microscope was employed to locate precisely the focal point of the lens as well as to measure the focal distance. The accuracy of the measurement was found to be ± 0.5%, with the focal length of the order of about 10 cm.     

用动力学法测杨氏模量实验及其实验装置的研制

介绍了新研制的用动力学法杨氏模量的实验装置,它具有功耗低,灵敏度高的特点;教学中测量悬线不同位置时的共振频率,用作图法拟合出其精确值,以减少由于悬线不在节点引起的系统性误差,测量精密度比传统的静态拉伸法高一个量级以上。     

XAuBi (X=Eu, Ba)拓扑nodal chain半金属的第一性原理研究

由于在磁性材料体系中缺失时间反演对称性,导致nodal chain被破坏,所以nodal chain通常存在于非磁材料中。但是,磁性材料EuAuBi是与常规磁性材料不同。本工作以第一性原理计算为研究方法,预言了在不考虑自旋轨道相互作用时,磁性材料EuAuBi体系为新型拓扑nodal chain半金属;当考虑自旋轨道耦合时,EuAuBi会退化为外尔半金属。对于非磁材料BaAuBi来说,在不考虑自旋轨道相互作用时,它同样是一种拓扑nodal chain半金属;当考虑自旋轨道相互作用时,由于C3旋转对称性的存在,BaAuBi会退化为狄拉克半金属。在XAuBi (X=Eu, Ba)中发现nodal chain半金属,会促进对六角材料的拓扑性质研究以及开拓其新实际应用领域。