用θ方法数值求解六角形节块时—空动力学方程

提出了一种求解六角形节块时－空动力学方程的数值方法，对空间变量采用节块展开法处理，用θ方法处理时间变量。为检验该算法的正确性，对几个瞬态算例进行了验算分析，计算结果和有关文献符合良好。最后对不同的时间步长及不同的θ因子对计算格式稳定性的影响进行了初步的讨论。     

节块积分方法求解圆柱几何对流扩散方程

为改进高温气冷堆中热工场方程的计算方法,研究了求解圆柱几何对流扩散方程的节块积分方法。针对圆柱几何下的r向横向积分方程的特殊性,提出了两种可行的近似方法--移项处理和常数近似,并进行相应的误差分析。数值计算结果表明：节块积分方法求解圆柱几何对流扩散方程的数值解具有迎风特性,一维和多维问题的计算结果均与解析解符合得很好;当节块在r向靠近零点时,常数近似带来的误差较移项处理带来的误差小,但当节块远离零点时,二者误差基本相当。     

求解中子扩散方程的半解析节块方法

为快速、精确实现反应堆堆芯多群中子扩散计算,采用基于横向积分技术的半解析节块方法（SANM）求解中子扩散方程,并结合解析粗网有限差分（ACMFD）方法,导出了基于半解析节块方法的粗网有限差分方程（CMFD）耦合系数。在半解析节块方法中,散射源和裂变源采用勒让德多项式,并在此假设下解析求解中子扩散横向积分方程。分别采用了零次、二次和四次勒让德多项式展开,以适应粗网和细网的计算。数值计算结果表明,所提出的方法具有很高的计算精度和计算效率。     

一种节块展开法

采用分离热中子方程的技巧 ,节块展开法解中子扩散方程 ,可提高计算速度 ,节省计算时间约 40 %。     

用于两群瞬态中子扩散计算的节块方法

采用非线性迭代节块方法作为空间近似方法 ,同时使用隐式差分和时间积分方法作为时间域近似方法 ,求解两群多维瞬态中子扩散方程 ,并给出基准问题的计算结果。     

节块展开法求解对流扩散方程的稳定性和数值耗散特性分析

本文研究了节块展开法求解对流扩散方程的稳定性和数值耗散特性。通过离散方程精确解和数值实验方法分析不同阶节块展开法的稳定性和数值耗散特性,并将其与有限体积法中的中心差分和一阶迎风格式进行对比。结果表明:偶数阶节块展开法的稳定是有条件的,即Peclet数(Pe)小于限值,且Pe限值会随展开阶数的增大而增大,其稳定性范围和精度均优于中心差分格式;奇数阶节块展开法是无条件稳定的,但随Pe的增大,数值耗散增大、计算误差增大,且当Pe大于一定值后,产生的数值耗散大于一阶迎风格式。     

解析函数展开节块法的数值稳定性分析

解析函数展开节块法是求解六角形几何堆芯中子注量率的常用方法,其核心思想是利用严格满足中子扩散方程的解析基函数,直接把节块内的各群中子注量率近似展开。研究发现,此方法在求解一些特殊问题时,出现计算发散的情况。针对此类问题,提出了两种解决方法：截断近似方法和泰勒展开方法。通过对改造后的VVER440基准题验证表明：这两种方法在保证计算精度的基础上很好地解决了解析函数展开节块法的数值不稳定性问题。     

非线性迭代半解析节块方法在压水堆燃料管理计算中的应用

利用半解析节块法的非线性迭代计算方法求解多维多群中子扩散方程。数值结果表明 ,这种非线性迭代半解析节块法具有较高的精度 ,并且在计算效率上比传统方法有很大的提高。将这种方法应用于SIMULATE E中 ,研制了压水堆燃料管理计算程序CSIM 3B ,对秦山核电厂第一循环进行了跟踪校核计算 ,取得了较好的计算结果     

基于非线性迭代的节块方法

使用粗网有限差分方法 ,同时利用解析方法通过非线性迭代对差分近似进行修正 ,解两群扩散方程。数值计算结果证明 ,这种非线性方法具有较高的精度 ,并且较传统方法计算效率成倍提高     

用于轻水堆扩散计算的格林函数节块展开法

提出一种用于轻水堆多维扩散计算的新方法─—格林函数节块展开法（ＧＮＥＭ），它是对节块展开法（ＮＥＭ）和节快格林函数方法（ＮＧＦＭ）的改进。参照ＮＥＭ把节块内偏中于通量用高阶多项式展开，利用格林函数求出节块表面偏中子通量与中子流的耦合关系，并用净中子流作为中间变量以简化计算。根据提出的理论编制了二维计算程序ＧＮＥＭ，对ＬＷＲ基准问题的数值计算表明，ＧＮＥＭ的计算精度与ＮＧＦＭ相当，而计算速度较ＮＧＦＭ和ＮＥＭ分别提高了１倍与１／３倍。应用ＧＮＥＭ取代ＰＳＵＩ－ＬＥＯＰＡＲＤ／ＡＤＭＡＲＣ中的扩散程序，对三里岛核电站（ＴＭＩ－１）的第１、第６循环进行了燃耗计算，其计算速度提高了３倍多，且具有更高的精度。     

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

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

三维圆柱几何格林函数节块法

发展了中子扩散计算三维圆柱几何格林函数节块法。通过横向积分将中子扩散方程化为３个互相耦合的一维偏通量方程。对于径向偏通量方程，将径向扩散微分算符分解为平板几何的扩散微分算符和１个附加项之和，将附加项移到方程右端作为１个附加源项，这样，３个方程都化为平板几何一维方程的形式。再借助平板几何第二类边界条件格林函数，建立偏通量积分方程。方程推导中，对圆柱形曲面几何的线积分和横向积分均需对相应的广义线元作积分，对于修正源项，通过分部积分方法将偏通量导数项转化为对格林函数的求导，通过源迭代法求解方程。基准计算表明，该计算精度高、速度快、成为三维圆柱几何堆芯计算的有效方法。     

热中子分离模型节块法

利用两群中子扩散方程的解析解，导出快、热中子通量之间的关系，去掉快群与热群方程之间的耦合，将原来的解两群方程的问题，理论上严格地转化成解单独的快群方程问题。快群方程用格林函数节块法进行求解。依照上述理论，研制了三维节块法程序ＮＧＭＤ。对一系列轻水堆基准问题的验算表明，ＮＧＭＤ的计算精度与知名的格林函数节块法程序ＮＧＦＭ相当，而计算时间约是ＮＧＦＭ的１３。