基于大规模并行计算的三维多群中子扩散方程有限差分方法

三维多群中子扩散方程的精确、高效求解是核动力堆芯设计及燃料管理的基础。应用有限差分方法求解该方程具有简便、精确、成熟的优点;然而,该方法的计算量和存储量均较大,极大地限制了它的计算规模和应用范围。本文基于大规模并行计算,研究三维多群中子扩散方程有限差分方法:采用中心有限差分格式离散中子扩散方程;基于MPI并行编程模型,采用空间区域分解的方式实现大规模并行计算;采用多群多区域耦合PGMRES算法进行并行加速。在集群服务器上开发了ParaFiDi程序,并采用IAEA3D,PHWR等多个基准题对该程序进行验证。数值结果表明,ParaFiDi程序具有较高的计算精度和计算效率。     

用SIP和SSIP方法求解三维反应堆中子扩散方程

本文首先简要介绍三维区域中的中子扩散方程及其差分方程的建立及求解的内、外迭代过程。然后导出内迭代过程中求解差分方程的SIP和SSIP公式,以及内、外迭代的加速方法和迭代参数的选取。最后给出一些数值结果。     

一类非线性微分方程的近似解析解

从理论上研究了一类广义扩散方程的求解问题. 利用相似变换和解析拆分技巧给出了求解该类非线性微分方程近似解的一种有效方法, 方程的解可以表示为一个收敛的幂级数. 近似解结果和数值结果非常符合，证明了所提出的方法的准确性和可靠性, 该方法可以用于解决其他科学和工程技术问题. 关键词： 广义扩散方程 非线性边界值问题 解析拆分 近似解析解     

解二维输运问题的三角网间断有限元法

本文研究了求解二维柱对称几何中子输运问题的三角网间断有限元法,证明了数值方法的稳定性,给出了一些数值计算结果,并与精确解和DSN等方法的计算结果进行了比较,取得了满意的结果。     

节块内嵌离散纵标(SN)方法的算法研究及程序开发

介绍节块内嵌离散纵标(S_N)方法求解三维堆芯中子输运/扩散方程的算法框架.在基于扩散理论的三维粗网节块展开方法(NEM)的算法体系中,用基于输运理论的径向二维细网节块离散纵标方法(NDOM)的内迭代过程,替代节块展开方法(NEM)内迭代的径向求解过程.该算法充分考虑了核电厂反应堆堆芯的三维结构特点,另一方面,也充分利用了已经成熟的三维粗网节块展开方法(NEM)和二维离散纵标方法(S_N)的研究成果,同时有效避免了利用离散纵标方法(S_N)求解三维中子输运方程所面临的计算内存和计算时间的瓶颈问题.编制开发二维多群节块离散纵标方法(NDOM)模块程序NS_NM和三维多群节块展开方法(NEM)模块程序MGNEM,并以此为基础编制开发节块内嵌S_N方法的模块程序HANWIND;其中,NS_NM为HANWIND求解两维问题的功能模块.针对OECD/NEA-2D C5G7MOX基准问题以及两环路核电厂三维堆芯的数值验算结果表明,节块内嵌S_N方法的算法开发及程序编制有效、切实可行.     

二维柱对称辐射输运基准模型及程序考核

对二维辐射扩散和输运模型问题,分别构造具有含外源项的精确解,解具有一定的物理特征且形式简单,并且辐射输运方程的解可以直接退化到辐射扩散的情形.给出LARED-R程序求解扩散模型数值解及分析.     

一维球几何菱形差分中子输运方程的扩散综合加速解法

本文导出了一维球几何定态中子输运方程菱形格式的扩散综合加速方程,并给出了差分公式。所给出的加速方法可以加速菱形格式的输运方程的迭代求解。并给出了部分模型的数值计算结果。     

二维柱几何中子输运方程的并行区域分解方法

分析不同的区域分解方法及优先级插入算法对二维柱几何下中子输运方程S_n间断有限元方程并行效率的影响,给出基于最小面体比的正方形区域分解方法及沿径向的优先级插入算法,并通过将正方形区域分解方法与径向优先级插入算法进行组合,形成新的算法.新算法更适应于二维柱几何下输运方程S_n间断有限元方法的并行计算.数值试验表明,在通信延迟较高的大型国产并行机上,新算法用数百个CPU还可以取得较好的并行效果,比已有方法具有更良好的可扩展性.     

基于界面捕捉的三维多介质辐射流体力学方程MMALE计算方法

针对界面重构产生的混合介质网格,建立基于能量守恒的子网格封闭模型,给出基于混合介质网格上的三维扩散方程求解方法.在此基础上,结合界面重构和流体力学方程MMALE （multi-material arbitrary Lagrangian-Eulerian）算法,给出一种针对三维辐射流体力学的MMALE计算方法.采用解析解算例,验证算法的正确性和精度.通过对典型辐射驱动问题的模拟,展示方法的健壮性.     

用改进的准静态方法求解反应堆三维瞬态方程

本文介绍了三维中子瞬态方程的改进准静态求解方法。给出了两个瞬态问题计算结果的比较;TWIGL结果与基准直接数值结果比较,其功率最大偏差为1.57%;CANDU堆模型,在功率上升达到最大值时的热中子通量与参考解的偏差为0.88%。因此,用改进的准静态方法求解三维中子瞬变问题,其结果是满意的,大大节省了计算时间。     

严谨而精密的创新——瑞利进行研究的特点

着重介绍了瑞利在测量氮气密度过程中发现惰性气体氩的研究和创新历程,以及瑞利对光学、声学等物理学的各个领域的主要贡献,指出正是求实的科学态度、锲而不舍的钻研精神和以创新为乐趣的事业追求是瑞利一生中取得丰硕的研究成果的主要原因。     

基于二维广义横向积分综合一维半解析的三维中子扩散节块法

目前, 轻水堆堆芯计算广泛采用基于横向积分技术的中子扩散方程节块展开法, 该方法需要对横向泄漏进行多项式近似而不严格, 而且堆芯设计还需要额外的精细功率重构模块用于获得组件内各栅元的功率分布。本文提出两维广义横向积分方法, 直接采用源展开以及表面流耦合方法, 可以避免上述两个不足。由于反应堆堆芯不均匀性更多发生在径向, 因此采用径向基于广义横向积分方法结合轴向常规节块法的综合方法进行三维中子扩散计算。通过基准问题的数值计算, 验证了该方法对于堆芯扩散计算的正确性和可行性。     

Multi-scale Analysis of Compressible Viscous and Rotating Fluids

We study a singular limit for the compressible Navier-Stokes system when the Mach and Rossby numbers are proportional to certain powers of a small parameter . If the Rossby number dominates the Mach number, the limit problem is represented by the 2-D incompressible Navier-Stokes system describing the horizontal motion of vertical averages of the velocity field. If they are of the same order then the limit problem turns out to be a linear, 2-D equation with a unique radially symmetric solution. The effect of the centrifugal force is taken into account.     

The Bethe lattice spin glass revisited

So far the problem of a spin glass on a Bethe lattice has been solved only at the replica symmetric level, which is wrong in the spin glass phase. Because of some technical difficulties, attempts at deriving a replica symmetry breaking solution have been confined to some perturbative regimes, high connectivity lattices or temperature close to the critical temperature. Using the cavity method, we propose a general non perturbative solution of the Bethe lattice spin glass problem at a level of approximation which is equivalent to a one step replica symmetry breaking solution. The results compare well with numerical simulations. The method can be used for many finite connectivity problems appearing in combinatorial optimization. Received 27 September 2000     

用多极理论分析圆柱对称微波谐振腔

 用多极理论计算具有复杂几何形状、圆柱对称微波腔的谐振频率,推导出用多极理论计算圆柱对称微波腔谐振频率的本征值方程。三个工程实例的计算结果表明,用多极理论计算圆柱对称微波腔谐振频率,不仅具有较高的计算精度, 而且可以很方便地应用于复杂几何形状圆柱对称微波腔工程问题的设计与计算,多极理论是计算圆柱对称微波腔谐振频率的一种有效方法。     

Method of solving diffusion boundary-value problems for a region with a boundary moving in accordance with an arbitrary law

A general method is proposed for solving the boundary-value problem of the diffusion equation in a limited region with a boundary that moves in accordance with an arbitrary law. The method is used to solve the first linear diffusion problem. Other boundary-value problems can be solved in similar fashion.Translated from Izvestiya Vysshikh Uchebnykh Zavedenii, Fizika, No. 12, pp. 97–101, December, 1970.     

Millimeter Wave Scattering of 2-D Frequency Selective Surface by Efficient Method of Lines with Preconditioned CG Technique

In this paper, both banded and symmetric successive overrelaxation (SSOR) preconditioned conjugate gradient (PCG) techniques are combined with method of lines (MOL) to further enhance the computational efficiency of this semi-analytic method. The electromagnetic wave scattering of 2-D frequency-selective surface is used as the examples to describe its implementation, whose analysis usually needs fast algorithms because of electrically large dimension. For arbitrary incident wave, helmholtz equation and boundary condition are used to calculate the impedance matrix and then to obtain reduced current-voltage linear matrix equation in spatial domain. Both banded and effective symmetric successive overrelaxation preconditioned conjugate gradient iterative method are chosen to solve this matrix equation. Our numerical results show that PCG methods can converge to accurate solution in much fewer iteration steps for analysis of the electromagnetic wave scattering from 2-D frequency-selective surface.     

Perturbed Lie Symmetry and Systems of Non-Linear Diffusion Equations

Abstract

The method of one parameter, point symmetric, approximate Lie group invariants is applied to the problem of determining solutions of systems of pure one-dimensional, diffusion equations. The equations are taken to be non-linear in the dependent variables but otherwise homogeneous. Moreover, the matrix of diffusion coefficients is taken to differ from a constant matrix by a linear perturbation with respect to an infinitesimal parameter. The conditions for approximate Lie invariance are developed and are applied to the coupled system. The corresponding prolongation operator is derived and it is shown that this places a power law and logarithmic constraints on the nature of the perturbed diffusion matrix. The method is used to derive an approximate solution of the perturbed diffusion equation corresponding to impulsive boundary conditions.     

A Boundary Meshless Method for Solving Heat Transfer Problems Using the Fourier Transform

Fourier transform is applied to remove the time-dependent variable in the diffusion equation. Under non-harmonic initial conditions this gives rise to a non-homogeneous Helmholtz equation, which is solved by the method of fundamental solutions and the method of particular solutions. The particular solution of Helmholtz equation is available as shown in [4, 15]. The approximate solution in frequency domain is then inverted numerically using the inverse Fourier transform algorithm. Complex frequencies are used in order to avoid aliasing phenomena and to allow the computation of the static response. Two numerical examples are given to illustrate the effectiveness of the proposed approach for solving 2-D diffusion equations.     

基于平面通量展开的三维Pin by Pin输运SP3计算

基于Pin by Pin输运SP₃计算是下一代反应堆物理计算方法中一个重要的研究方向。节块法一般采用高阶通量和中子源展开以适应较大的节块尺寸,因此对于Pin by Pin计算来说效率偏低。由于反应堆堆芯不均匀性更多发生在径向,本文提出一种径向基于二维平面通量展开结合轴向常规节块法的综合方法进行三维Pin by Pin输运SP₃计算,其三个方向的节块耦合迭代均采用基于中子净流耦合的格式。通过IAEA 2D/3D基准问题和典型输运基准问题计算,验证该计算方法对于堆芯扩散计算和输运计算的正确性和可行性。