二维土壤溶质输运方程的有限体积元格式

本文导出二维的土壤溶质输运方程的有限体积元格式, 并分析其误差.通过数值例子说明, 有限体积元格式比有限元格式稳定.     

非定态中子迁移方程的离散纵标—多群逼近理论

本文对有界凸的非均匀介质中具各向异性散射和裂变的连续能量中子迁移的非定态方程,将方向和能量两个变量同时离散的所谓离散纵标——多群逼近方法建立起系统的数学理论,证明了: 1 非定态迁移方程的解,可由相应的非定态离散纵标——多群迁移系统的解逼近。 2 原迁移算子的占优本征值,可由离散纵标——多群迁移算子所确定的具非负本征函数且实部为最大的本征值逼近。 3 原迁移算子的占优本征值所相应的正本征函数,可由离散纵标——多群迁移算子的实部为最大的本征值所相应的非负本征函数逼近。 4 估计了各种逼近的阶。     

非线性中子输运方程全离散有限元解的误差估计

本文研究非线性中子输运方程的全离散有限元方法,我们用Brower不动点定理证明了全离散格式非线性方程组的可解性,并给出有限元解的误差估计。     

二维柱坐标系辐射输运方程保球对称的离散纵标格式

 离散纵标格式是计算辐射输运方程的常用格式之一. 但是, 传统的离散纵标格式求解二维柱坐标系辐射输运方程模拟一维球对称问题时, 会出现明显的非对称现象, 球对称性被破坏. 针对该问题, 本文分析了传统离散纵标格式不能够保持球对称性的原因, 提出了空间基于柱坐标系、方向基于球坐标系的辐射输运方程, 并对该方程设计了新的离散纵标格式, 从理论上证明了当空间网格取球对称剖分时该离散格式能够保持一维球对称性的充分必要条件. 通过对真空球区域辐射输运、与物质耦合辐射输运等球对称算例的数值模拟, 验证了该格式的保球对称性, 球对称误差能够达到机器精度. 非对称辐射驱动模型以及非对称网格剖分条件下的数值模拟等算例也取得了较好的结果.     

算子方程离散格式判稳的充分条件

萨马斯基曾给出判稳的充分条件,但它涉及估计算子的模或特征值,这在一般情况下是困难的,不容易检验.本文给出一种易于检验的充分条件,即把稳定性与一代数方程组的Jacobi迭代法的迭代阵的最大特征值联系起来,从而可利用迭代法收敛的某些已知结果来判别稳定性,其中特别方便的是利用矩阵的对角占优条件.本方法的特点是适用于一般的非均匀有限元剖分.     

非线性二维Volterra积分方程的一个高阶数值格式

对非线性二维Volterra积分方程构造了一个高阶数值格式.block-byblock方法对积分方程来说是一个非常常见的方法,借助经典block-by-block方法的思想,构造了一个所谓的修正block-by-block方法.该方法的优点在于除u(x_1,y),u(x_2,y),u(x,y_1)和u(x,y_2)外,其余的未知量不需要耦合求解,且保存了block-by-block方法好的收敛性.并对此格式的收敛性进行了严格的分析,证明了数值解逼近精确解的阶数是4阶。     

二维抛物型方程的一族高精度显式差分格式

构造了一族解二维抛物型方程的高精度显格式 ,其稳定性条件为r=Δt/Δx2 =Δt/Δy2 <1 /2 ,截断误差为O(Δt3 +Δx4)     

带有不连续系数的线性输运方程差分格式的收敛性

提出了一个基于三角形网格的显式差分格式逼近带有不连续系数的线性输运方程. 通过对数值解的有界性、TVD(total variation decreasing)和空间、时间方向的平移估计, 利用Kolmogorov紧性原理证明了数值解在L¹_loc模下收敛于初值问题的唯一弱解.从而得到了初值问题解的存在唯一性和关于初值的稳定性. 数值算例表明本文提出的格式计算方便而且比 Lax-Friedrichs格式更有效.       

二维抛物型方程的一族高精度分支稳定显格式

1 引言 在渗流、扩散、热传导等领域中经常会遇到求解二维抛物型方程的初边值问题 {(6)u/(6)=a((6)2u/(6)x2+(6)2u/(6)y2), 0＜x,y＜L,t＞0,a＞0u(x, y, 0) =φ(x, y), 0 ≤ x, y ≤ L (1)u(0,y,t) =f1(y,t),u(L,y,t) =f2...     

二维变系数扩散方程交替分组显格式的稳定性分析

1.引言 由于高性能并行计算机的出现和并行计算的推动,十多年来,抛物型方程有限差分并行算法设计与分析一直受到关注. D.J.Evalns和A.R.B.Abdullah(1983,[1,2]利用Saul’yev非对称格式对常系数抛物方程设计了AGE(交替分组显格式)算法,并用矩阵分析的方法证明了该算法的无条件稳定性.该算法有明显的并行性,倍受推崇,且计算的实践([8],[9])表明它对变系数的抛物方程也是可行的,但稳定性的分析成为一个难点.张宝琳([3])在一维情     

Convergence of the modified discrete ordinates method for the anisotropic scattering transport equation

Criticality problem of nuclear tractors generally refers to an eigenvalue problem for the transport equations. In this paper, we deal with the eigenvalue of the anisotropic scattering transport equation in slab geometry. We propose a new discrete method which was called modified discrete ordinates method. It is constructed by redeveloping and improving discrete ordinates method in the space of L¹(X). Different from traditional methods, norm convergence of operator approximation is proved theoretically. Furthermore, convergence of eigenvalue approximation and the corresponding error estimation are obtained by analytical tools.     

守恒型扩散方程非线性离散格式的性质分析和快速求解

对于守恒型扩散方程,研究其二阶时间精度非线性全隐有限差分离散格式的性质,证明了其解的存在唯一性.研究了二阶时间精度的Picard-Newton迭代格式,证明了迭代解对原问题真解的二阶时间和空间收敛性,以及对非线性离散解的二次收敛速度,实现了非线性问题的快速求解.本文中方法也适用于一阶时间精度格式的分析,并可推广至对流扩散问题.数值实验验证了二阶时间精度Picard-Newton迭代格式的高精度和高效率.     

一类偏积分微分方程一阶差分全离散格式

本文给出了数值求解一类偏积分微分方程的一阶差分全离散格式。时间方向采用了一阶向后差分格式，空间方向采用二阶差分格式，给出了稳定性的证明，误差估计及收敛性的结果，并给出了数值例子。     

An economical difference scheme for heat transport equation at the microscale

Heat transport at the microscale is of vital importance in microtechnology applications. In this article, we proposed a new ADI difference scheme of the Crank‐Nicholson type for heat transport equation at the microscale. It is shown that the scheme is second order accurate in time and in space in the H¹ norm. Numerical result implies that the theoretical analysis is correct and the scheme is effective. © 2004 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2004     

Convergence of discrete schemes for the Perona-Malik equation

We prove the convergence, up to a subsequence, of the spatial semidiscrete scheme for the one-dimensional Perona-Malik equation ut=(?^′x(ux)), , when the initial datum is 1-Lipschitz out of a finite number of jump points, and we characterize the problem satisfied by the limit solution. In the more difficult case when has a whole interval where is negative, we construct a solution by a careful inspection of the behaviour of the approximating solutions in a space-time neighbourhood of the jump points. The limit solution u we obtain is the same as the one obtained by replacing ?(⋅) with the truncated function min(?(⋅),1), and it turns out that u solves a free boundary problem. The free boundary consists of the points dividing the region where |ux|>1 from the region where |ux|?1. Finally, we consider the full space-time discretization (implicit in time) of the Perona-Malik equation, and we show that, if the time step is small with respect to the spatial grid h, then the limit is the same as the one obtained with the spatial semidiscrete scheme. On the other hand, if the time step is large with respect to h, then the limit solution equals , i.e., the standing solution of the convexified problem.     

Monotone difference schemes stabilized by discrete mollification for strongly degenerate parabolic equations

The discrete mollification method is a convolution‐based filtering procedure suitable for the regularization of ill‐posed problems and for the stabilization of explicit schemes for the solution of PDEs. This method is applied to the discretization of the diffusive terms of a known first‐order monotone finite difference scheme [Evje and Karlsen, SIAM J Numer Anal 37 (2000) 1838–1860] for initial value problems of strongly degenerate parabolic equations in one space dimension. It is proved that the mollified scheme is monotone and converges to the unique entropy solution of the initial value problem, under a CFL stability condition which permits to use time steps that are larger than with the unmollified (basic) scheme. Several numerical experiments illustrate the performance and gains in CPU time for the mollified scheme. Applications to initial‐boundary value problems are included. © 2010 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 28: 38–62, 2012     

Convergence rates to the discrete travelling wave for relaxation schemes

This paper is concerned with the asymptotic convergence of numerical solutions toward discrete travelling waves for a class of relaxation numerical schemes, approximating the scalar conservation law. It is shown that if the initial perturbations possess some algebraic decay in space, then the numerical solutions converge to the discrete travelling wave at a corresponding algebraic rate in time, provided the sums of the initial perturbations for the -component equal zero. A polynomially weighted norm on the perturbation of the discrete travelling wave and a technical energy method are applied to obtain the asymptotic convergence rate.

     

The global decay to discrete shocks for scalar monotone schemes

Given a family of discrete shocks of a monotone scheme, we prove that the discrete shock profile with rational shock speed is asymptotically stable with respect to the topology : if , then as under no restriction conditions of the initial data to the interval . The asymptotic wave profile is uniquely identified from the above family by a mass parameter.

     

Concatenating construction of the multisymplectic schemes for 2+1-dimensional sine-Gordon equation

In this paper, taking the 2+1-dimensional sine-Gordon equation as an example, we present the concatenating method to construct the multisymplectic schemes. The-method is to discretizee independently the PDEs in different directions with symplectic schemes, so that the multisymplectic schemes can be constructed by concatenating those symplectic schemes. By this method, we can construct multisymplectic schemes, including some widely used schemes with an accuracy of any order. The numerical simulation on the collisions of solitons are also proposed to illustrate the efficiency of the multisymplectic schemes.