通量展开节块法求解六角形几何二维中子扩散方程

提出了一种在二维六角形几何节块内数值求解中子扩散方程的节块法，节块内的各群通量分布用解析基函数近似展开，节块之间采用面偏流0次矩和1次矩进行耦合．为了提高计算效率，将响应矩阵技术应用于迭代求解过程．该方法不仅避免了常规方法中偏通量积分出现的奇异项问题，而且具备重构六角形几何节块内精细功率分布的能力．基于提出的模型，发展了二维六角形组件中子扩散计算程序FEMHEX．通过对水水动力反应堆基准问题的校验计算表明，该方法能高效、准确地给出有效增值因数及节块功率．     

基于通量展开节块法的六角形时-空动力学方程数值解

提出了一种基于节块内瞬态中子通量展开的六角形几何时一空动力学方程数值解法．在该方法中，各群中子通量分布用解析基函数和二阶正交多项式近似展开，而包含各组缓发中子先驱核浓度的固定源项则利用多项式进行近似．将面平均偏流及其一次矩作为节块之间的耦合条件，不但明显改善了节块耦合关系，而且使得响应矩阵技术比较容易地应用于迭代求解过程．对二维、三维基缝问题计算表明，该方法能高教、准确地给出各时间寿内的堆芯总功率和节块功率分布。     

一种新型的数值计算方法——基函数法

我们提出一种新型的数值计算方法——基函数法．此方法直接在非结构网格上离散微分算子．采用基函数展开逼近真实函数,构造出了导数的中心格式和迎风格式．取多项式为基函数并采用通量分裂法及中心格式和迎风格式相结合的技术以消除激波附近的非物理波动,我们构造出数值求解无黏可压缩流动一阶多项式的基函数格式．通过一、二、三维多个无黏可压缩流动典型算例的数值计算表明本方法是一种高精度的,对激波具有高分辨率的无波动新型数值计算方法,与网格自适应技术相结合可得到十分满意的结果．     

通量展开节块法求解六角形几何三维多群中子扩散方程

提出了一种在三维六角形几何节块内数值求解多群中子扩散方程的节块法，该方法把节块内各群中子通量分布用解析基函数近似展开．为了改善节块耦合关系，采用了一种新的节块边界条件：面平均偏流0次矩和1次矩同时保持连续．将响应矩阵技术应用于迭代求解过程，使得该方法具有较高的计算效率．通过对基准问题的校验计算表明，该方法能准确地给出有效增值系数以及节块功率；对于二维多群问题，所有基准题的组件最大功率偏差均小于1％．     

基于非结构网格的中子输运方程计算方法研究

针对非结构网格下中子输运方程的求解方法问题，分别研究了球谐函数（Pn）有限元方法、离散纵标（Sn）有限元方法、三角形穿透概率方法和三角形节块Sn方法的计算模型．结果表明：有限元方法具有比较高的计算精度；Pn方法避免了Sn方法中的“射线效应”现象，而Sn有限元方法则具有容易处理真空边界条件的优点；三角形穿透概率方法的计算时间最快；三角形节块Sn方法则可以保证在计算较快的情况下同时具有较高的计算精度．根据上述4种方法分别研制了相应的计算程序，利用基准题进行了验证并分析了各自的优缺点，计算结果表明每种方法均能得到满意的计算精度．     

格林函数节块法求解中子伴随通量方程

将中子伴随通量方程变换处理以后，利用格林函数节块法的原理求解中子伴随通量，以此代替过去采用最多的细网有限差分法。对ＴＷＩＧＬ二维模型进行了校算，并与有限差分法ＣＩＴＡＴＩＯＮ程序结果进行了比较。两种方法的计算结果十分吻合，这表明采用粗网格的格林函数节块法求解中子伴随通量是可行的，而且这种粗网格的方法具有空间网格大，计算精度高，计算速度快的特点，因而具有较高的计算效率。     

求解具有复杂地形二维浅水方程的修正HLL格式

采用非结构网格的有限体积方法,对具有复杂地形的二维浅水方程进行数值模拟.采用HLL近似Riemann解计算界面数值通量,基于三角形网格,底坡源项采用简单的斜底模型离散,保证了地形的离散精度,摩阻源项采用全隐方式求解以保证格式的稳定性.采用多维重构及多维限制器的方法获得空间二阶精度的格式,时间离散采用三阶Runge-Kutta法以获得高阶时间精度.为保证格式的和谐性,对经典的HLL格式计算的数值通量中的静水压力项进行了修正.数值计算的结果验证此格式具有良好的高精度捕捉间断的能力,可以应用到地形复杂的二维浅水问题计算中去.     

二维无结构网格浅水波方程的高阶非振荡有限体积法

基于二维浅水波方程,对无结构网格给出了一种三阶精度非振荡有限体积方法,方法的主要思想是先对每一个三角形单元采用最小二乘的思想构造一个二次插值多项式,而在计算交界面的流通量时采用了两点高斯公式以保证格式的整体精度,时间离散采用三阶TVD Runge-Kutta方法。最后用该格式对圆柱形溃坝问题和倾斜水跃问题进行了数值试验,得到了满意的结果。     

移动网格上解二维Euler方程的WENO型有限体积法

对二维Euler方程在基于弹簧技术的移动非结构三角形网格上给出了一种WENO型有限体积格式,方法的主要思想是先对每一个三角形单元构造插值多项式,而在计算交界面的流通量采用两点高斯积分公式以保证格式的整体精度,时间离散采用三阶TVD Runge-Kutta方法.最后用该格式对一些典型算例进行了数值模拟,结果表明该方法计算速度快,对间断有很好的分辨能力.     

二阶迎风有限体积法方腔流数值模拟

介绍了一种基于非结构网格的二阶迎风有限体积离散格式，在对流项的离散过程中，为了达到二阶精度，在界面上对物理量咖作Taylor展开，在处理展开项中的梯度时经入NND格式的优点，克服了中心差分格式不稳定的缺点，并对方腔流动进行了系统的数值模拟．计算网格采用三角形网格，节点数为12960，单元数为25600．压强修正基于SIMPLEC方法．给出了雷诺数Re达到5000时的定常流动结果，与以往方腔流计算的标准解非常吻合，但所用的网格数要少．计算结果表明该二阶迎风有限体积算法有很好的收敛性和稳定性，且松弛因子对计算结果影响很小．     

六角形几何多群中子扩散的解析节块法

通过寻找二维亥姆霍兹方程解析试验解的方式来求解六角形几何多群中子扩散方程,并利用对称性原理和群论确定节块内的中子通量分布．求解过程不采用横向积分技术,且得到的解在节块内任意点上都满足扩散方程．对ＩＡＥＡ基准题的校核计算表明,该方法具有很高的精度,最大组件功率误差小于２％．     

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

从一阶三维中子输运方程出发，对方向变量采用离散纵标方法展开，得到一系列关于空间变量的偏微分方程，从而避免了二阶方程由于分母上存在截面，不能准确描述内含真空介质的问题．对这些关于空间变量的方程采用最小二乘有限元方法进行离散，形成的刚度矩阵是对称的，因此可以采用快速迭代方法求解．据此编制了三维中子输运方程的非结构网格离散纵标计算程序，并采用三棱柱元素和四面体元素剖分对一系列基准问题做了验算．计算结果表明，该方法能用于非结构网格，并具有较高的计算精度，对多数问题，有效增值系数的误差都小于0．3％，通量误差都小于3．0％．     

桁梁结构空间计算的三节点截面桁段有限元法

利用二次多项式位移函数导出一个桁架结构的三节点截面桁段单元,提出了桁梁空间计算的三节点截面桁段有限元法,该方法适当地考虑了主桁弦杆节点刚性,桥面系横梁局部弯曲,斜桥门架及横联横梁局部弯曲的影响,算例表明,计算结果可靠、自由度较少,比较适用于大跨度桁架桥的空间结构计算分析。     

基于Burton-Miller边界元法的不同类型单元计算精度对比

为了提高边界元法的计算精度和对具有复杂边界形状实际问题的应用能力,发展并应用非连续线性和二次边界单元进行数值计算.使用传统边界积分方程计算外声场,通过带有解析解数值算例,对比不同类型单元的计算精度,得到最有效的单元类型.然而使用传统边界积分法,在某些虚假特征频率处会产生解的非唯一性问题,Burton-Miller方法可以有效地克服这一问题.基于Burton-Miller法得到的非连续线性和二次单元的优化节点位置并不在勒让德多项式零点位置上,虽然表现得不像传统边界元法那样规律和统一,但是合适的经验值仍然被给出.     

三维离散纵标中子输运稳态及瞬态动力学计算

基于著名输运程序DOT4．2，开发了三维离散纵标中子输运稳态及瞬态动力学计算程序．对时间变量采用直接的无条件稳定且具有高精度的全隐式向后差分格式进行离散处理，将瞬态输运方程转化为各个离散时间步上的固定源型稳态中子输运方程进行求解．通过对TAKEDA基准题及输运瞬态基准题的校验计算表明：该三维稳态程序能准确地给出有效增殖系数以及各区的区域平均中子通量密度，其相对误差分别小于0．1％和4．0%；对于缓发超临界和瞬发超临界问题，三维瞬态程序的计算结果均与参考值吻合良好．本三维程序可以为复杂反应堆堆芯的临界计算和瞬态动力学特性分析提供一个有效而准确的分析工具和手段．     

第二类边界条件先进格林函数节块法

目前广泛应用的节块格林函数法是基于第三类边界条件,这类格林函数比较复杂,并且不便直接使用不连续因子,必须经过转换后才能使用和功率重构。该文发展了第二类边界条件三维几何先进格林函数节块法。该方法的优点是格林函数计算较简单,便于在交界面处引入通量不连续因子,并通过功率重构得出堆芯内精细通量分布。基准计算表明,该程序精度高、速度快,成为水堆堆芯物理设计和燃料管理计算的核心程序。