非线性迭代解析节块方法的数值稳定方法

针对非线性迭代解析节块方法在特定条件下出现的数值不稳定问题,提出了一种解决方法,并给出了数值结果。     

厚非线性介质瞬态热光非线性效应的研究

瞬态热光非线性效应是指由脉冲宽度可以和非线性效应的建立时间相比拟的纳秒脉冲激光入射产生的热光非线性效应.用Zscan法测量了不同厚度、不同束腰半径碘的酒精溶液的瞬态热光非线性效应曲线.从声波的传播方程和光波的非线性传播方程出发,用数值计算方法模拟了这一非线性过程.并对实验结果和数值模拟结果进行了比较 关键词： 瞬态热光非线性效应 非线性光学 数值计算     

六角形燃料组件内精细核密度分布计算

介绍一种可用于六角形燃料组件的精细计算方法，在理论模型上编制了精细计算程序ＪＸＪＳ：ＪＸＪＳ计算结果与用文（１）提供的方程计算出的结果进行了比较。并对结果进行了误差分析。     

非线性抛物型方程迭代加速计算方法

研究非线性抛物型方程隐式格式的迭代加速求解方法，包括三方面内容：一是构造具有二阶收敛性的非线性迭代方法，二是迭代初值的选取方法，三是证明迭代方法的保正性。     

液体燃料熔盐堆三维瞬态耦合分析

根据熔盐堆的流动特性,采用多群扩散理论,建立液体燃料熔盐堆的三维中子动力学模型和流动传热模型,开发针对液体燃料熔盐堆的三维核热耦合程序,并分析瞬态情况下MOSART堆的物理热工特性.结果表明：堆芯热工参数的瞬态变化直接影响物理参数的变化,液体熔盐堆具有较明显的温度负反馈特性.     

克尔介质中瞬态热光非线性效应的作用

研究了介质中瞬态热光非线性效应对于克尔效应Z扫描测量结果的影响，观察到了Z扫描峰谷特性的转变，相应的数值模拟结果与实验结果很好地吻合.实验也发现，克尔效应和瞬态热光非线性效应相互抵消时Z扫描曲线呈现出新奇的单峰或平台状曲线.对于用纳秒脉冲激光进行的Z扫描测量结果的正确分析具有重要意义. 关键词： 瞬态热光非线性效应 克尔效应 Z扫描     

基于JSNT程序的压水堆屏蔽计算

反应堆屏蔽计算是评估核电站安全性能的基础,是指导电站设计、运行的重要手段之一。JSNT程序是中物院高性能数值模拟软件中心研制的大规模并行离散纵标输运程序,具有较高的计算精度和计算效率。利用JSNT程序对某压水堆进行了建模计算,给出了中子通量密度的分布结果,并与实验测量值进行了对比。对比结果显示：无论是采用S8计算还是S16计算,计算结果都能满足工程要求;相比S8而言,采用S16计算可以显著提高计算精度,能够将某些测点处的相对误差降至1%以内。     

基于JSNT程序的压水堆屏蔽计算

反应堆屏蔽计算是评估核电站安全性能的基础,是指导电站设计、运行的重要手段之一。JSNT程序是中物院高性能数值模拟软件中心研制的大规模并行离散纵标输运程序,具有较高的计算精度和计算效率。利用JSNT程序对某压水堆进行了建模计算,给出了中子通量密度的分布结果,并与实验测量值进行了对比。对比结果显示:无论是采用S8计算还是S16计算,计算结果都能满足工程要求;相比S8而言,采用S16计算可以显著提高计算精度,能够将某些测点处的相对误差降至1%以内。     

基于FLUENT的核热耦合程序反应性反馈参数敏感性

将计算流体力学模型与中子动力学模型耦合来进行反应堆瞬态安全分析的方法,由于可以开展复杂几何结构的三维流动传热分析,因此受到很大的关注。基于FLUENT用户自定义功能（UDF）开发了一套可用于池式铅堆瞬态安全分析的核热耦合程序,程序耦合了临界/次临界点堆中子动力学模型和燃料棒模型。由于反应堆处于不同寿期时,随着燃料燃耗、可燃毒物积累等因素导致反应性反馈系数有较大变化,因此使用开发的核热耦合程序对中国科学技术大学提出的小型自然循环铅冷快堆进行不同关键反馈系数下无保护的瞬态超功率事故安全分析。调整点堆模块考虑到的四个反应性反馈系数,可以发现燃料多普勒系数对堆安全的影响最大,同时定量的分析结果表明超功率事故引入时间长短对事故演化有重要影响。     

基于FLUENT的核热耦合程序反应性反馈参数敏感性

将计算流体力学模型与中子动力学模型耦合来进行反应堆瞬态安全分析的方法,由于可以开展复杂几何结构的三维流动传热分析,因此受到很大的关注。基于FLUENT用户自定义功能（UDF）开发了一套可用于池式铅堆瞬态安全分析的核热耦合程序,程序耦合了临界/次临界点堆中子动力学模型和燃料棒模型。由于反应堆处于不同寿期时,随着燃料燃耗、可燃毒物积累等因素导致反应性反馈系数有较大变化,因此使用开发的核热耦合程序对中国科学技术大学提出的小型自然循环铅冷快堆进行不同关键反馈系数下无保护的瞬态超功率事故安全分析。调整点堆模块考虑到的四个反应性反馈系数,可以发现燃料多普勒系数对堆安全的影响最大,同时定量的分析结果表明超功率事故引入时间长短对事故演化有重要影响。     

基于线性子链法的压水堆裂变产物源项一体化计算方法

核电厂设计基准源项计算可为核电厂安全评审提供依据,同时也是辐射屏蔽计算的基础。基于压水堆堆芯、一回路和气载源项的研究基础,类比衰变常数引入了迁移常数和核反应常数的概念,进而总结了一体化计算上述源项中裂变产物源项的源项方程。针对源项方程变系数、大型、稀疏和刚性的特点,在时间离散近似的基础上,基于线性子链算法编写程序求解了上述方程。通过与典型压水堆工程文件对比,证明了程序的正确性和必要性。     

基于线性子链法的压水堆裂变产物源项一体化计算方法

核电厂设计基准源项计算可为核电厂安全评审提供依据,同时也是辐射屏蔽计算的基础。基于压水堆堆芯、一回路和气载源项的研究基础,类比衰变常数引入了迁移常数和核反应常数的概念,进而总结了一体化计算上述源项中裂变产物源项的源项方程。针对源项方程变系数、大型、稀疏和刚性的特点,在时间离散近似的基础上,基于线性子链算法编写程序求解了上述方程。通过与典型压水堆工程文件对比,证明了程序的正确性和必要性。     

大型先进压水堆堆芯机械补偿控制模式

研究了AP1000堆芯物理计算程序的计算方法,分析了该程序尤其是机械补偿（MSHIM）控制模式计算功能的理论模型、计算方法及应用范围,对比了大型先进压水堆堆芯设计与AP1000的差异,评估了AP1000反应堆MSHIM计算功能在大型先进压水反应堆堆芯MSHIM计算功能的适用性。基于大型先进压水堆堆芯物理建模,针对大型先进压水反应堆首循环及平衡循环的典型燃耗点进行典型的100%-70%-100%和100%-50%-100%负荷跟踪运行模式计算分析,并依据计算结果对大型先进压水堆的MSHIM运行模式进行了分析。针对大型先进压水反应堆首循环及平衡循环,开展机械补偿控制初始启动运行、再启动运行计算分析,研究机械补偿控制模式的堆芯初始启动和再启动运行能力。计算结果表明：采用MSHIM运行模式的大型先进压水堆,不调节硼浓度的情况下,在首循环、平衡循环典型燃耗下具有一定的负荷跟踪能力;启动、再启动运行模式则需要配合调节硼浓度才能完成。     

大型先进压水堆堆芯机械补偿控制模式

研究了AP1000堆芯物理计算程序的计算方法,分析了该程序尤其是机械补偿（MSHIM）控制模式计算功能的理论模型、计算方法及应用范围,对比了大型先进压水堆堆芯设计与AP1000的差异,评估了AP1000反应堆MSHIM计算功能在大型先进压水反应堆堆芯MSHIM计算功能的适用性。基于大型先进压水堆堆芯物理建模,针对大型先进压水反应堆首循环及平衡循环的典型燃耗点进行典型的100%-70%-100%和100%-50%-100%负荷跟踪运行模式计算分析,并依据计算结果对大型先进压水堆的MSHIM运行模式进行了分析。针对大型先进压水反应堆首循环及平衡循环,开展机械补偿控制初始启动运行、再启动运行计算分析,研究机械补偿控制模式的堆芯初始启动和再启动运行能力。计算结果表明:采用MSHIM运行模式的大型先进压水堆,不调节硼浓度的情况下,在首循环、平衡循环典型燃耗下具有一定的负荷跟踪能力;启动、再启动运行模式则需要配合调节硼浓度才能完成。     

基于运行反馈的压水堆氚排放量研究北大核心CSCD

以压水堆核电厂中氚的产生机理和氚源项计算模型为基础,结合对国内外大量压水堆核电厂的氚排放运行数据的系统性分析,识别出冷却剂硼酸活化和次级中子源活化是压水堆氚排放量的主要来源,其中对中国广核集团运行机组,锑铍中子活化后的产氚量对氚年排放量的贡献可达到40%,而氚从完整的锆合金包壳的燃料棒中的释放是可以忽略不计的。由于优化次级中子源是降低压水堆氚排放量的唯一有效措施,通过分析建议压水堆核电厂采用双层不锈钢包壳的次级中子源或者取消次级中子源以降低压水堆氚排放。     

基于不连续因子的沸水堆芯三维瞬态数值模拟

基于不连续因子校正的粗网格有限差分法是实现堆芯瞬态三维数值模拟的高效方法之一, 粗网节块的界面不连续因子与边界反照率的计算方法决定了实时数值模拟过程中的精度。在计算不连续因子的过程中, 省去了细网节块计算与粗网均匀化过程, 直接在粗网格划分情况下, 基于节块展开法和非线性迭代策略, 推导了粗网格界面不连续因子比率与边界反照率的计算公式, 并编制了相应的计算程序。沸水堆典型算例的三维瞬态模拟证实该方法可在空间域和时间域两方面, 使静态、瞬态精度均达到与先进节块法相等同的程度, 并且计算效率优于先进节块法, 为核电站全范围模拟机三维堆芯的实时仿真模型开发提供了一种切实可行的选择。     

基于Christopherson迭代的超精密加工流场分析方法

随着特种超精密加工技术的发展,复杂流体被越来越多地用于超精密加工工艺中。超精密加工流场分析具有几何特征复杂、流体本构特性多样、流体边界为自有边界等特点,传统流体数值分析方法难以实现可靠分析。从流体的一般特性出发,将D. G. Christopherson提出的非负二阶偏微分系统的超松弛迭代方法用于超精密加工流场分析,建立了适应性与可靠性兼顾的流场分析方法。以磁流变抛光为例,开展了抛光区域压力场数值计算,结果表明所得压力分布形态正确,且分布从x轴正半轴延伸到负半轴,与郑立功等人的实验测定结果一致。另外,基于Kistler力传感器对磁流变抛光过程的法向压力在0.1~0.3 mm浸深段进行了在位测量,发现计算与实验结果偏差均小于20%。表明了该方法的有效性与准确性。     

Conformable variational iteration method,conformable fractional reduced differential transform method and conformable homotopy analysis method for non-linear fractional partial differential equations

In this paper, we introduce conformable variational iteration method (C-VIM), conformable fractional reduced differential transform method (CFRDTM) and conformable homotopy analysis method (C-HAM). Between these methods, the C-VIM is introduced for the first time for fractional partial differential equations (FPDEs). These methods are new versions of well-known VIM, RDTM and HAM. In addition, above-mentioned techniques are based on new defined conformable fractional derivative to solve linear and non-linear conformable FPDEs. Firstly, we present some basic definitions and general algorithm for proposal methods to solve linear and non-linear FPDEs. Secondly, to understand better, the presented new methods are supported by some examples. Finally, the obtained results are illustrated by the aid of graphics and the tables. The applications show that these new techniques C-VIM, CFRDTM and C-HAM are extremely reliable and highly accurate and it provides a significant improvement in solving linear and non-linear FPDEs.     

蒙特卡罗与离散纵标耦合方法在压水堆堆腔漏束计算中的应用

为克服蒙特卡罗(MC)方法计算时间长和离散纵标(SN)方法复杂几何描述不精确的困难,采用SN-MC耦合计算流程,研究了基于蒙特卡罗方法和离散纵标方法的耦合方法。耦合方法的主要思想是根据离散纵标程序提供的中子角通量,利用接口程序计算出面源的分布概率,然后由修改后的源抽样子程序生成包括上下圆面源和圆柱面源的组合源,提供给蒙特卡罗程序进行计算分析。初步计算结果表明,该耦合方法是正确的,可用于压水堆堆腔漏束的计算分析。     

基于压缩感知理论的非线性γ谱分析方法

γ谱分析是一种重要的放射性核素定量分析方法. 弱峰的检测和重叠峰的分解是γ 谱分析中的难点. 为了解决这两个问题, 基于压缩感知理论, 提出了一种新的γ 谱分析方法. 这一方法从谱仪对γ 谱调制的物理机理出发, 通过数学建模, 将γ 谱分析转化为一个以真实γ 谱为解的求逆问题, 并在压缩感知理论框架下, 运用γ 谱特征峰的稀疏性, 进行逆问题的求解, 直接获得γ 谱的估计结果. 数值模拟结果和蒙特卡洛模拟结果表明: 该方法能在降低统计涨落的同时, 有效减小谱仪调制带来的能谱展宽, 从而提高γ 谱分析精度. 关键词： γ谱分析')" href="#">γ谱分析 压缩感知 非线性 逆问题