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

基于离散节块方法（DNTM）的思想，研究了一种新的适用于二维三角形节块的离散坐标数值解法．使用坐标变换将任意三角形变换为正三角形，源项和通量的空间分布使用二次正交多项式展开，横向泄漏项的空间二次分布通过使用二元二次多项式近似正三角形节块的通量分布来获得、针对小尺寸易出现不收敛的情况，提出了一种改进的节块平衡有限差分方法．数值计算结果表明，该方法能满足非结构几何区域问题的求解需求，与细网SN方法相比，该方法在较粗的网格下具有高的精度，同时还适用于较细网格下的中子输运计算．     

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

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

基于Monte Carlo方法的中子标识输运计算

为了考察中子与介质相互作用的细致物理过程,得到不同种类中子对所关心物理量的贡献,利用Monte Carlo方法的模拟能力,由中子标识记录中子及其次级粒子的输运历史,在程序中赋予中子一个整型标识变量,记录中子的出身、碰撞过程等历史,在记数时根据该变量的值进行分门别类的统计以达到求解众多物理量的目的。结果表明,该方法能得到其它方法不易得到的物理量,如不同种类中子出壳流等。     

非线性中子输运方程的流线扩散有限元法

0 引言在核反应系统中,若考虑中子间的相互作用,则导致非线性中子输运方程.文[1]讨论了一类非线性中子输运方程解的存在性和唯一性.文[2～4]分别研究了非线性中子输运方程的半离散和全离散Galerkin有限元法,得到了解的存在性和相应误差估计.本文的目的在于讨论非线性中子输运方程的流线扩散法. 流线扩散法(SD),又名SUPG方法,最早被Hughes,Johnson及其合作者们研究.对于非定常问题,其数学分析基于空间-时间有限元法.到目前,流线扩散法已成为解决对流占优扩     

非等距网格高精度差分方法用于气动声学问题计算

研究了非等距网格下高阶精度有限差分方法用于气动声学问题的可行性．通过Taylor级数展开法构造了不等距网格下的七点六阶精度空间离散格式，分析了格式适用的波数范围，时间积分采用显式四阶精度推进格式，可直接计算非定常欧拉方程用于气动声学问题．算例采用随机变化的网格间距，对一维单波方程和球形波方程进行了计算，验证了非等距网格格式模拟波动问题的能力．对二维情况，计算了亚音速均匀流中初始声、涡和熵脉冲波问题，得到了很好的结果．     

基于非结构Cartesian网格的电磁散射场计算

计算流体力学中的网格技术和数值方法对于计算电磁学的应用具有借鉴意义。该文结合计算流体力学中的数值方法,建立了基于非结构Cartesian网格上的时域有限体积法,用于求解电磁散射问题。为保证这种方法在时间和空间上具有二阶精度,在时间离散上采用二步Runge-Kutta方法,空间离散采用非结构的二阶NND格式。通过求解时域中的Maxwell方程组,计算了二维翼型和三维乘波体外形等典型完全导电散射体的雷达截面。计算结果表明:该方法计算精度高,适用于复杂外形的电磁散射问题的求解。     

湄洲湾潮流和污染物输运的高精度非结构网格模拟

针对湄洲湾不规整岸边界和复杂的潮流变化等特征,模型采用高精度的二维非结构网格浅水流模型对湾内潮流水动力进行模拟。基本方程离散和有限体积法思想,并利用改进的QUICK格式提高对流项计算精度,在扩散通量估算上引入具有良好数值性能的超松弛校正方法。不合理的压力波由Rhie-Chow的动量插值方法加以抑制,水位校正方程结合SIMPLEC程式导出。通过对湄洲湾小潮的水动力模拟,结果表明,潮位、流速、流向的实测值和计算值吻合良好。另外,模型还对海湾内肖厝峰尾和鲤鱼尾两处主要污水排放口的COD排放进行了计算,模拟结果较好地反应了湄洲湾内污染物输运特征,为更好地进行近海水污染控制提供有力的手段。     

三维Euler方程的非结构网格非振荡有限体积法

考虑三维Euler方程,对非结构四面体网格给出了一类时空分离的非震荡有限体积方法。空间方面采用Roe通量进行离散,同时基于物理量守恒运用最小二乘和极值原理在每个四面体单元构造非振荡线性重构多项式,时间离散采用两步Runge-Kutta TVD方法。文末对三维球形爆炸问题进行数值模拟,得到了满意的结果。     

Level Set方法在非结构网格上的应用

采用基于非结构网格的Level Set方法对多介质流动问题进行数值模拟.采用有限体积法求解Euler方程,控制面通量的计算采用HLLC(Hartern,Lax,van Leer,Contact)近似Riemann解方法.Level Set方程也采用有限体积法求解,采用Lax-Friedchs方法计算通量,通过窄带(Narrow Band)方法来减少计算量,重新初始化采用Fast Marching Method.界面的处理采用Ghost Fluid方法.Runge-Kutta显式时间推进,时间、空间都是二阶精度.对SOD激波管问题进行数值模拟,其结果和精确解相吻合;对空气/氦气泡相互作用问题进行模拟,取得较好的结果.     

非结构网格RKDG方法求解多介质界面问题

在多介质流体力学方程数值模拟中,界面的计算是一个非常重要的问题.采用非结构网格上RKDG(Runge-Kutta Discontinuous Galerkin)有限元方法,给出了捕捉多介质界面的Lev-el Set方程和重新初始化方程的计算方法.数值实验表明多介质界面在非结构网格计算上,RKDG有限元方法具有很大的优点.     

中子输运问题源项反演的反幂法

对于包括裂变反应在内的中子输运源项反演问题，研究关于源项有效倍增因子的惫一特征值问题的求解．基于球谐函数展开和有限差分离散，给出了中子输运方程的源项反演逼近的反幂算法，该方法的优势是在适当的初值条件下可以显著提高计算速度．计算结果表明，在对有效倍增因子有较好的预先估计值的情况下，反幂法迭代3步，误差就为0．04545％，而乘幂法迭代20步，误差为0．109％，由此可以看出反幂法计算速度更快，计算结果更精确。     

板对称中子迁移方程解的稳定性

本文对非均匀介质中具各向异性散射和裂变的连续能量的板对称中子迁移方程，在吸收不小于裂变的条件下，证明了其解随时间ｔ趋于无穷而指数稳定于零。     

穿透概率法求解二维六角形几何多群中子积分输运方程

讨论了用穿透概率法求解二维六角形几何轻水堆燃料组件中子通量的分布,提出了一种简易计算模型。子区内中子源采用各向同性平源分布,子区量在方向上采用简化６Ｐ１近似。根据提出的模型,编制了ＴＰＨＥＸ＿Ａ程序,并对一些轻水堆六角形组件问题作了计算,计算结果与ＤＲＡＧＯＮ程序计算结果以及ＭＣ结果进行了比较,符合较好,本程序可用于轻水堆六角形燃料组件计算。     

平板中子迁移方程多群逼近的一致收敛性

本文在连续函数空间中讨论了具连续能量的平板中子迁移方程.运用有界线性算子的积分半群理论证明了该方程的解可由相应的多群中子迁移方程的解一致逼近.     

Boltzmann积分—微分方程解的渐近行为

本文利用[1]中的占优本征值理论,改进[2][3]中的有关结果,获得与时间有关的中子迁移理论中动态解的更具体和更精确的性质估计。     

具广义边界条件的平板中子迁移方程解的渐近性质

引进上下解的概念，并利用其性质，讨论了平板几何上有一般反射边界条件的中子迁移方程解的渐近性质。     

输运方程数值解法分析

输运方程的数值解法有很多种,该文将介绍其中的两种数值方法：其一为向前差分法,其二为向后差分法,通过这两种方法求解同样输运方程会得到一些有趣的结果,并将计算结果与解析解的结果进行比较,最后对数值结果加以分析。     

与时间有关的中子迁移问题的经典解

本文研究了与时间有关的中子迁移方程经典解整体存在唯一性以及中子分布密度随时间的延长而衰减的渐近性质。