中子输运计算界面流方法的微扰计算

以六角形几何中子积分输运计算界面流算法及其相对应的数学共扼方程计算为基础,利用微扰原理计算了当反应堆六角形组件中栅元核参数发生微量变化时系统反应性的变化。计算结果表明,本文所开发的基于六角形几何中子积分输运算法的微扰计算方法是正确的。     

界面流方法二维六角形组件计算

讨论了用界面流方法计算二维六角形组件中子通量分布。从积分输运方程出发,导出了一种简便的数学模型,在子区内采用平源通量近似,并假设中子发射和散射为各向同性。在子区表面上,中子通量的空间分布为常数,中子角通量分布通过伴随勒让特多项式展开表示,采用DP_1近似。推导出界面流方程组,给出了泄漏、穿透几率矩阵及其矩阵元素的表达式及计算方法。根据提出的数学模型,编制了TPHEX程序,对二维六角形组件进行了计算,本程序可用于水堆六角形燃料组件计算。     

中子积分输运动态方程

通过理论分析给出了中子积分输运动态方程 ,发展了中子积分输运理论 ,使中子积分输运理论不仅可以用来分析反应堆栅格非均匀效应和计算反应堆参数等稳态问题 ,而且还可以处理反应堆动态问题。中子积分输运动态方程是一个多群多点 (一个空间分区为一点 )中子动态方程 ,在单群情况下就是多点反应堆动态方程。多点动态方程可以用来分析与空间有关的反应堆动态问题。介绍了中子积分输运动态方程的应用个例 ,通过中子积分输运动态方程分析了中国先进研究堆中子代时间的构成 (刚性和柔性中子代时间 )问题。     

用递推公式计算栅元有效共振积分

本文讨论数值求解反应堆栅元中子积分输运慢化方程,计算燃料及其它吸收体有效共振积分的方法。文中推导了适用范围较广的慢化源递推公式,这些慢化源递推公式的应用大大简化了慢化源的计算。程序计算表明方法是有效的,并用以分析计算了某些情况下慢化剂窄其振近似及完全窄共振近似产生的误差。     

基于三棱柱网格的穿透概率方法研究

为了求解三维非结构几何中子输运方程,研究了基于三棱柱的穿透概率方法.三棱柱子区内部中子源采用平源分布近似,子区界面中子角注量率空间分布采用均匀分布近似,角度分布采用简化P1近似.在此基础上编制了适应于三维几何的中子输运问题求解的程序TPMTD,并对三维Takeda输运基准问题进行了校算,计算结果与参考值吻合良好.     

基于非结构网格的三维中子输运方程求解方法

本文使用离散纵标-间断有限元方法求解了三维中子输运方程,它对能量变量采用多群近似离散,对方向变量采用离散纵标法离散,对空间变量采用间断有限元离散;并研究了每个S_N离散方向的有限元网格的排序以及中子输运DFEM方程中几个矩阵的矩阵元的精确求解方法,并据此开发了基于非结构网格的三维输运计算程序TetTran1.0。基准例题校核结果表明,该程序具有很高的计算精度。     

非结构网格中子输运方程的球谐函数解法研究

从新的二阶自共扼角通量密度(SAAF：Self-Adjoint Angular Flux)中子输运方程出发．利用球谐函数对角度变量进行展开,导出了一组关于空间变量的偏微分方程组,中子通量密度的各个分量相互耦合,应用一定的迭代策略进行迭代求解。针对每一个方程,应用有限元方法对非结构网格进行离散求解。据此编写了二维球谐函数方法输运计算程序,对一系列基准题进行校算的数值结果表明,该方法具有较高的计算精度,克服了射线效应,并能用于非结构网格。     

基于三角形网格的穿透概率方法研究

研究了基于三角形网格的中子输运方程穿透概率方法.网格子区内部中子源采用平源分布,子区界面中子角注量密度空间分布设为均匀,角度分布采用简化4Pi近似,并采用一种按照非结构网格编码顺序进行网格扫描的新方法.编制了适应于二维任意形状几何中子输运特征值问题求解的程序TPTRI,对一些二维输运基准问题作了计算.与MG-MCNP3B、SURCU、TEPFEM等程序结果进行了比较分析,结果吻合很好.     

刚性限制法在瞬态中子输运计算中的应用

刚性限制法（SCM）可有效缓解中子动力学方程中的刚性问题,可采用较大时间步长获得同等计算精度,提高计算效率。现有SCM主要用于求解两群瞬态中子扩散方程。本文将SCM应用于求解多群瞬态中子输运方程,在原有中子输运方程特征线方法求解程序PEACH的基础上,增添了瞬态求解功能,开发了PEACH-K程序。采用OECD/NEA最新发布的基准题C5G7-TD对PEACH-K程序进行数值验证,结果表明,PEACH-K程序在大时间步长下仍具有很高的计算精度,且具有良好的数值稳定性。     

反应性平衡方程及其应用

１引言稳态堆物理计算的基本出发点是中子守恒。由中子守恒推导出的中子输运方程及其扩散近似是描述堆物理现象的基本方程。反应堆核计算的绝大部分内容是用扩散方程完成的。但对某些堆物理现象的研究来说，扩散方程应用起来并不方便，需要寻找中子守恒的另一种表达形式。...     

穿透概率法求解二维六角形轻水堆燃料组件中子通量密度分布

研究利用穿透概率法求解二维六角形轻水堆燃料组件内中子通量密度分布。子区内中子源采用线性分布,子区表面通量密度在方向上采用简化6P1近似。提出了六角形组件周边水隙的处理方法。根据提出的模型,编制了TPHEX－C程序,并对六角形组件进行了计算,结果与蒙特卡罗方法计算的结果符合良好。     

压水堆六角形燃料组件均匀化 计算软件包TPFAP-HEX

介绍了所研制的具有工程实用价值的压水堆六角形燃料组件均匀化计算软件包。该组件中子空间能谱的计算采用穿透概率法与响应矩阵法相结合的方法，在六角形几何内求解中子积分输运方程。在此方法中，栅元内中子源采用空间线性或二次近似，栅元表面中子通量密度角分布采用简化6P     

离散节块输运方法在二维柱坐标系中的应用

讨论粗网离散节块法在二维柱坐标系内中子输运方程数值求解中的应用.给出方法的数学模型,在节块内中子通量采用二次近似,表面泄漏采用常数近似.根据所提出的理论模型,编制了考虑各向异性散射的二维节块输运程序DNSN/2D,并对一系列检验和基准问题作了计算.计算结果表明,与传统S_N方法相比,它具有很高的计算效率,在很粗的节块分割下具有非常高的精度.     

一种求解反应堆中子扩散方程和输运方程的方法

结合反应堆中子物理方程的特点和非结构网格技术,提出了一种在任意网格下求解反应堆中子物理方程的方法。以反应堆中子扩散方程和一阶离散纵标的输运方程为例,从空间离散、方程离散和边界条件实施等方面介绍了该方法的实施过程。利用该方法编制的程序计算了BN-600基准题、CFR1000概念堆和BWR栅元3种情况的有效增殖因数,并与其他程序的计算结果进行了对比,初步验证了该方法的合理性。     

基于OpenMP的中子输运方程特征线法并行计算研究

特征线法是目前求解反应堆中子输运方程的主要计算方法之一。本文开发了基于OpenMP的中子输运方程特征线法并行计算程序,以提高特征线法的计算效率。OpenMP是共享存储体系结构上的一个并行编程模型,采用Fork-Join并行执行方式,适合于SMP共享内存多处理系统和多核处理器体系结构。通过相关基准题测试验证,表明所开发的程序在有效增殖因数以及相对中子通量(归一化栅元功率)分布等参数上都能取得良好的精度,且使用OpenMP能取得良好的加速效果,使计算时间显著减少。     

六角形轻水堆组件中子通量密度分布的计算

介绍利用穿透概率法求解二维六角形轻水堆燃料组件中子通量密度分布。子区内中子源及通量密度在空间上采用二次分布 ,子区表面通量密度在空间上采用平通量密度分布 ,在方向上采用简化 6P1近似。根据提出的模型 ,编制了TPHEX D程序 ,并对一些轻水堆六角形组件问题作了计算 ,计算结果与MC结果进行了比较 ,符合良好。本程序可用于六角形轻水堆燃料组件计算。     

二维中子输运方程的非结构网格离散纵标数值解法

从一阶中子输运方程出发,对方向变量采用离散纵标法展开,得到一系列关于空间变量的偏微分方程,对这些方程采用最小二乘有限元方法进行离散.编制了二维中子输运方程的非结构网格离散纵标计算程序.对一系列基准问题做了验算,计算结果表明,该方法能用于非结构网格,并具有较高的计算精度.