解二维中子输运问题的三角网格配置法

本文研究了二维轴对称中子输运(包括各向同性散射和各向异性散射)问题的一种有限元数值方法。给出了一系列的数值计算结果,并与国内外的SN、间断有限元方法的数值结果及实验结果进行比较,取得了令人满意的结果。     

Towel型域到方形计算域全局同胚网格的数值生成方法

对任意维Towel型物理域 D到方形计算域D_c 的全局同胚映射数值实现方法进行研究,给出了任意维等参数正变换的数值迭代求解方法,并就迭代求解方法的误差进行分析,给出了数值网格点互不重合的一个充分性条件.最后,根据该方法,给出了钱塘江口的二维不重叠数值网格生成算例.     

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

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

各向异性扩散DLA的标度性质

本文利用坐标变换的方法讨论了各向异性扩散对DLA集团标度行为的影响。对二维正方形点阵,解析计算了各向异性扩散DLA集团的Hausdorff维数。结果表明,此时生长集困的Hausdorff维数D随各向异性扩散概率p连续变化,其最大值为5/3,最小值为3/2,还将这里所得到的理论结果与Jullien等人的数值结果进行了比较,发现它们符合得很好。最后,还解析讨论并计算了各向异性扩散DLA的广义维数D_q(p)。     

渤海潮波问题的有限元方法

本文对文[1]中提出的二维长波方程给出了数值求解方法。由于此方法无需处理边界上的积分,不仅使数值计算的数据准备工作量大大减少,缩短了计算时间,降低了计算误差,而且没出现文[2]提到的短波长噪音干扰。用本文建立的方法,计算了我国渤海的M₂分潮,计算结果与实测结果吻合较好。     

流函数Navier-Stokes方程的交替方向部分隐式和完全隐式解法

本文给出了数值求解以流函数表示的Navier-Stokes方程的交替方向部分隐式和交替方向完全隐式方法;作为算例,计算了低雷诺数的二维方腔流动问题,并与已知结果进行了比较。     

非线性薛定谔方程的高阶分裂改进光滑粒子动力学算法

为提高传统光滑粒子动力学(SPH)方法求解高维非线性薛定谔(nonlinear Schr?dinger/Gross-Pitaevskii equation, NLS/GP)方程的数值精度和计算效率,本文首先基于高阶时间分裂思想将非线性薛定谔方程分解成线性导数项和非线性项,其次拓展一阶对称SPH方法对复数域上线性导数部分进行显式求解,最后引入MPI并行技术,结合边界施加虚粒子方法给出一种能够准确、高效地求解高维NLS/GP方程的高阶分裂修正并行SPH方法.数值模拟中,首先对带有周期性和Dirichlet边界条件的NLS方程进行求解,并与解析解做对比,准确地得到了周期边界下孤立波的奇异性,且对提出方法的数值精度、收敛速度和计算效率进行了分析;随后,运用给出的高阶分裂粒子方法对复杂二维和三维NLS/GP问题进行了数值预测,并与其他数值结果进行比较,准确地展现了非线性孤立波传播中的奇异现象和玻色-爱因斯坦凝聚态中带外旋转项的量子涡旋变化过程.     

一种适用于各向异性材料的修正PUFF物态方程

给出一种适用于各向异性材料的修正PUFF物态方程,既能体现各向异性材料在压缩和膨胀过程中材料体积变化的非线性特征,又能反映材料的各向异性强度效应.以正交各向异性材料二维碰撞问题和X射线辐照下的二维热击波为例,利用有限元方法编写程序进行数值模拟,对根据传统PUFF物态方程和修正PUFF物态方程计算得到的平均正应力值作对比研究,结果表明,低压和拉伸状态下二者存在明显差异,在高压下二者趋于一致.     

量子自旋系统的块自旋方案及其簇团Monte Carlo模拟方法

本文提出了一种关于一维和二维量子自旋系统的Monte Carlo模拟新算法。具体给出了对于一维和二维Heisenberg模型的计算方案。在s=1/2的铁磁和反铁磁链情况下的计算结果证实了此方法有效地克服了“临界慢化(CSD)”效应和能够达到更低的温度区域进行计算。该方法可推广到各向异性和高自旋(s=1,3/2,…)量子系统、二维经典8-顶角格点模型和低维费密系统。 关键词：     

定常不可压Navier-Stokes型方程的非线性Galerkin有限元算法的收敛性

给出数值求解二维定常不可压Navier-Stokes型方程的非线性Galerkin有限元算法,并分析了数值解的正则性和收敛性,当粗网格参数H和细网格参数h满足关系式H=O(h^1/2)时,该算法具有和Galerkin有限元算法同阶的收敛精度,然而在计算上比Galerkin有限元算法更为简单,可以节省可观的计算量.最后给出了数值试验,验证了上述结果。     

The discrete ordinates method for anisotropic scattering in the neutron transport theory: The critical sphere problem

The spherical harmonics method for anisotropic scattering in the neutron transport theory related to the critical sphere problem was investigated by Yildiz [The spherical harmonics method for anisotropic scattering in neutron transport theory: the critical sphere problem. JQSRT 2001;71:25-37]. Some numerical results and figures that they provided are incorrect. The correct numerical results for the critical radius are obtained and tabulated for different scattering parameters by using the discrete ordinates method.     

解各向异性散射输运问题的双向间断有限元方法

本文以球几何为例研究了解各向异性散射中子输运问题的双向间断有限元方法、介绍了算法,给出了逼近解的范数估计,证明了解的唯一性和算法的稳定性,并给出了数值结果。     

球几何中子输运保正线性间断有限元格式

针对球几何中子输运方程线性间断有限元方法计算的负中子通量问题,构造了保正线性间断有限元格式,该格式保持中子角通量0阶矩和1阶矩。现有方法计算中子角通量非负时,采用传统的线性间断有限元方法,求解线性方程组;原方法计算出现负通量,则采用构造的保正格式,求解非线性方程组。编制了球几何中子输运问题保正格式程序模块,并集成到应用程序。数值算例表明构造的保正格式计算的中子通量非负,有效降低数值误差,提高数值计算的精度。     

球几何中子输运保正线性间断有限元格式

针对球几何中子输运方程线性间断有限元方法计算的负中子通量问题,构造了保正线性间断有限元格式,该格式保持中子角通量0阶矩和1阶矩。现有方法计算中子角通量非负时,采用传统的线性间断有限元方法,求解线性方程组;原方法计算出现负通量,则采用构造的保正格式,求解非线性方程组。编制了球几何中子输运问题保正格式程序模块,并集成到应用程序。数值算例表明构造的保正格式计算的中子通量非负,有效降低数值误差,提高数值计算的精度。     

The solution of the neutron transport equation in a slab with anisotropic scattering

The neutron transport equation for a slab geometry with the extremely anisotropic scattering kernel is considered. The albedo and transmission factors are calculated using the variation method. The effect of the extremely anisotropic parameter on the variation of the slab albedo and transmission factor is calculated. The obtained results are compared with the published data.     

On pseudo-continuities in the eigenvalue continuum of one-speed neutron transport theory

We report on a new closed-form approximant to the singular eigenfunction transient solution of the one-speed, one-dimensional neutron transport equation with an anisotropic scattering kernel. It is proved that the associated eigenvalue continuum must be pointwise perforated.     

Space-angle distribution of the reflected charged particles adjusted for spin calculations

The efficiency of the proposed method of solution for the spin-considered highly anisotropic multiple scattering problem is based, following optics, upon the admission that the complete solution of the vectorial transport equation boundary problem is the superposition of the anisotropic and smooth parts, calculated separately. We use a vectorial small angle modification of the spherical harmonics method to evaluate the anisotropic part and the vectorial discrete ordinates method to obtain a smooth one. Some calculation examples for reflected radiation are given.     

Milne's problem for anisotropic scattering medium with specular reflecting boundary

The Milne problem is investigated subject to reflecting boundary conditions. The original version of the problem with vacuum boundary condition is generalized assigning, to the surface x=0, a specular reflection coefficient . Linearly anisotropic case is studied. The integral version of the transport equation solved using trial functions based on Case's eigenvalues and exponential integral function. Solution of the Milne problem is formulated in terms of characteristic parameters such as extrapolated end point, emergent angular distribution and total neutron density. Numerical results for the analytically evaluated parameters are then present. Some of our numerical results are compared with the available published results.     

Chapman-Enskog-maximum entropy method on time-dependent neutron transport equation

The time-dependent neutron transport equation in semi and infinite medium with linear anisotropic and Rayleigh scattering is proposed. The problem is solved by means of the flux-limited, Chapman-Enskog-maximum entropy for obtaining the solution of the time-dependent neutron transport. The solution gives the neutron distribution density function which is used to compute numerically the radiant energy density E(x,t), net flux F(x,t) and reflectivity R_f. The behaviour of the approximate flux-limited maximum entropy neutron density function are compared with those found by other theories. Numerical calculations for the radiant energy, net flux and reflectivity of the proposed medium are calculated at different time and space.