一种有效的非线性规划算法的应用

一、引言目标函数为非线性可微函数时,我们采用一种变尺度算法,结合比例因子法的一维寻查作为非线性规划方法。在微型计算机上处理了用TDPAC方法研究辐射损伤实验数据的非线性多参数曲线拟合问题,获得了理想结果。     

基于非线性预处理JFNK的中子-热工联立求解

利用Jacobian-free Newton-Krylov(JFNK)方法联立求解中子-热工耦合问题,采用非线性预处理方式,以避免求解非线性残差,使得JFNK具有可以充分利用原有的中子-热工计算程序,易于实现"黑箱"耦合的特点。对非线性预处理的相关性质进行分析,同时对非线性预处理与线性预处理的区别与联系以及计算效率进行理论分析。以二维简化中子-热工耦合模型作为算例,对比非线性预处理/线性预处理JFNK方法、传统耦合求解方法的计算效率。结果表明:非线性预处理/线性预处理JFNK方法的计算效率比传统方法具有明显优势,线性预处理的计算效率高于非线性预处理。     

计算α本征值的蒙特卡罗期望估计方法

从齐次中子输运方程出发,推导了计算α本征值的数学表达式,并给出了三种蒙特卡罗期望估计式,它们分别是:径迹长度期望估计;碰撞期望估计;逃逸期望估计。不同模型的计算结果表明,期望估计方法的计算效率可比MCNP中的(α,k)回归迭代方法有量级的提高。尤其对含慢化材料的次临界模型,期望估计方法有着特殊的优势。     

反应堆节块计算中一种数值求解边界反照率的方法

提出一种数值近似求解反照率的计算方法。在节块扩散计算中，用该方法求得的反照率边界条件，替代压水堆中的反射层。该方法数学模型简单，且能够处理各种几何边界条件，包括“Ｌ”型和“”型拐角。一些基准题的数值结果表明，以新的边界条件进行的堆芯功率和临界计算所得结果与精确表示反射层的节块计算结果符合良好。     

隐式重启的Arnoldi方法及其在高阶谐波求解中的应用

Krylov子空间方法的出现是近年来大型线性方程组和特征值问题求解领域的重大进展,介绍其中一类适用于求解反应堆k-本征值问题的隐式重启的Arnoldi方法(IRAM),以及该方法在高阶谐波求解中的应用。研究结果表明,IRAM方法求解高阶谐波具有和源修正法同样的精度,但计算速度更快,尤其是当所求的谐波阶次较高时,IRAM方法可获得10倍以上的速度优势。同时,IRAM方法还具备较好的处理重特征值问题的能力。     

多种影响因素下的设备易损性计算

抗震裕度评估是核电厂地震安全评估的方法之一，通过地震易损性分析计算高置信度低失效概率的抗震能力值是抗震裕度评估中很重要的一步。本文对于同时受到多种失效模式影响的设备易损性计算进行了研究，讨论了蒙特卡罗抽样方法和拉丁超立方分布抽样方法在设备易损性计算中的应用，对两种抽样方法的计算效率和准确度进行了评价。结果表明，在小样本抽样计算时拉丁超立方抽样方法有更好的计算效率和收敛速度，在1 000次样本数量时，两种抽样方法计算结果均可达到收敛。     

基于因子函数的α-IHS抽样方法

核电厂概率安全分析中,由于数据源存在不确定性,导致无法进行准确评估,因此需开展不确定性分析。样本在空间分布的均匀特性是不确定性分析的关键因素,不同的样本分布导致不确定分析结果差异较大。传统的拉丁超立方抽样方法在样本空间分布均匀性方面未进行优化,改进分布式超立方抽样(IHS)方法通过保持样本点之间的最优距离来实现空间均匀分布,但其最优距离只能在理想分布下达到最优。为改进IHS设计上的缺陷,提出了基于因子函数的α-IHS方法,利用修正因子α来优化IHS方法中的最优距离。结果表明,该方法较IHS方法具有更均匀的空间分布,提高了抽样效率。     

隐伏铀矿勘查中非放射性方法最佳组合的应用

本文阐述了铀矿勘查方法发生的历史性变化,并以实例介绍电磁、重力及磁测等非放射性方法的最佳组合在隐伏铀矿勘查中的应用。