基于多层增量未知元方法的一类三维对流扩散方程的研究

对于一类一般形式的三维对流扩散方程, 运用有限差分方法, 在增量未知元方法(IU)下, 可以得到一个IU型正定但非对称的线性方程组.其系数矩阵条件数要远远优于不用IU方法的情形^[1]. 考虑到IU方法的这一优点, 作者在文中将IU方法与几种经典的迭代方法相结合, 来求解上述系统. 作者从理论上对该系统的IU型系数矩阵条件数进行了估计, 并通过数值试验验证了这几种IU型迭代方法的有效性.     

Incremental unknowns in three-dimensional stationary problem

The main purpose of this work is to set up the explicit matrix framework appropriate to three-dimensional partial differential equations by means of the incremental unknowns method. Multilevel schemes of the incremental unknowns are presented in the three space dimensions, and through numerical experiments, we confirm that the incremental unknowns method is efficient and the hierarchical preconditioning based on the incremental unknowns can be applied in a more general form.      

Incremental unknowns for solving partial differential equations

Summary Incremental unknowns have been proposed in [T] as a method to approximate fractal attractors by using finite difference approximations of evolution equations. In the case of linear elliptic problems, the utilization of incremental unknown methods provides a new way for solving such problems using several levels of discretization; the method is similar but different from the classical multigrid method.In this article we describe the application of incremental unknowns for solving Laplace equations in dimensions one and two. We provide theoretical results concerning two-level approximations and we report on numerical tests done with multi-level approximations.     

Bi-parameter incremental unknowns ADI iterative methods for elliptic problems

Bi-parameter incremental unknowns (IU) alternating directional implicit (ADI) iterative methods are proposed for solving elliptic problems. Condition numbers of the coefficient matrices for these iterative schemes are carefully estimated. Theoretical analysis shows that the condition numbers are reduced significantly by IU method, and the iterative sequences produced by the bi-parameter incremental unknowns ADI methods converge to the unique solution of the linear system if the two parameters belong to a given parameter region. Numerical examples are presented to illustrate the correctness of the theoretical analysis and the effectiveness of the bi-parameter incremental unknowns ADI methods.     

Tangential frequency filtering decompositions for unsymmetric matrices

Summary. In this paper, tangential frequency filtering decompositions (TFFD) for unsymmetric matrices are introduced. Different algorithms for the construction of unsymmetric tangential frequency filtering decompositions are presented. These algorithms yield for a specified class of matrices equivalent decompositions. The convergence rates of an iterative scheme, which uses a sequence of TFFDs as preconditioners, are independent of the number of unknowns for this class of matrices. Several numerical experiments verify the efficiency of these methods for the solution of linear systems of equations which arise from the discretisation of convection-diffusion differential equations. Received April 1, 1996 / Revised version received July 4, 1996     

The effect of block red-black ordering on blockILU preconditioner for sparse matrices

It is well known that the ordering of the unknowns can have a significant effect on the convergence of a preconditioned iterative method and on its implementation on a parallel computer. To do so, we introduce a block red-black coloring to increase the degree of parallelism in the application of the blockILU preconditioner for solving sparse matrices, arising from convection-diffusion equations discretized using the finite difference scheme (five-point operator). We study the preconditioned PGMRES iterative method for solving these linear systems.     

提高反应—扩散方程有限差分格式的稳定性问题

This paper deals with the special nonlinear reaction-diffusion equation.The finite difference scheme with incremental unknowns approximating to the differential equation (2.1) is set up by means of introducing incremental unknowns methods.Through the stability analyzing for the scheme,it was shown that the stability conditions of the finite difference schemes with the incremental unknowns are greatly improved when compared with the stability conditions of the corresponding classic difference scheme.     

Nonlinear Galerkin method in the finite difference case and wavelet-like incremental unknowns

Summary The IMG algorithm (Inertial Manifold-Multigrid algorithm) which uses the first-order incremental unknowns was introduced in [20]. The IMG algorithm is aimed at numerically implementing inertial manifolds (see e.g. [19]) when finite difference discretizations are used. For that purpose it is necessary to decompose the unknown function into its long wavelength and its short wavelength components; (first-order) Incremental Unknowns (IU) were proposed in [20] as a means to realize this decomposition. Our aim in the present article is to propose and study other forms of incremental unknowns, in particular the Wavelet-like Incremental Unknowns (WIU), so-called because of their oscillatory nature.In this report, we first extend the general convergence results in [20] by proving them under slightly weaker conditions. We then present three sets of incremental unknowns (i.e. the first-order as in [20], the second-order and wavelet-like incremental unknowns). We show that these incremental unknown can be used to construct convergent IMG algorithms. Special stress is put on the wavelet-like incremental unknowns since this set of unknowns has theL ² orthogonality property between different levels of unknowns and this should make them particularly appropriate for the approximation of evolution equations by inertial algorithms.     

Spacetime discontinuous Galerkin methods for convection-diffusion equations

We have developed two new methods for solving convection-diffusion systems, with particular focus on the compressible Navier-Stokes equations. Our methods are extensions of a spacetime discontinuous Galerkin method for solving systems of hyperbolic conservation laws [3]. Following the original scheme, we use entropy variables as degrees of freedom and entropy stable numerical fluxes for the nonlinear convection term. We examine two different approaches for incorporating the diffusion term: the interior penalty method and the local discontinuous Galerkin approach. For both extensions, we can show an entropy stability result for convection-diffusion systems. Although our schemes are designed for systems, we focus on scalar convectiondiffusion equations in this contribution. This allows us to highlight our main ideas behind the stability proofs, which are the same for scalar equations and systems, in a simplified setting.     

Numerical study of the incremental unknowns method

Through numerical experiments we explore the incremental unknowns method; linear stationary as well as evolutionary problems are considered. For linear stationary problems, the method appears as a nearly optimal two-step preconditioning technique. At each step, we observe that the convergence behavior of the iterative methods employed is dramatically different, depending upon whether or not preconditioning is used. For linear evolutionary problems, successful and sharply accurate long-term integration is observed when the incremental unknowns (other than that of the coarsest level) are, effectively, small quantities. Otherwise, systematical aliasing arises. © 1994 John Wiley & Sons, Inc.     

Spectral methods for some singularly perturbed third order ordinary differential equations

Spectral methods with interface point are presented to deal with some singularly perturbed third order boundary value problems of reaction-diffusion and convection-diffusion types. First, linear equations are considered and then non-linear equations. To solve non-linear equations, Newton’s method of quasi-linearization is applied. The problem is reduced to two systems of ordinary differential equations. And, then, each system is solved using spectral collocation methods. Our numerical experiments show that the proposed methods are produce highly accurate solutions in little computer time when compared with the other methods available in the literature.      

A characteristic difference method for the transient fractional convection-diffusion equations

A new characteristic finite difference method for solving the two-sided space-fractional convection-diffusion equations is presented, by combining characteristic methods and fractional finite difference methods. Stability, consistency and (therefore) convergence of the new method are discussed in this paper. An error estimate is given. Numerical experiments of this method are carried out and compared with other known methods.     

一类多维非线性反应扩散方程有限差分格式的稳定性问题(英文)

本文研究了一类多维线性反应扩散方程差分格式的稳定性.利用量未知元方法,建立了具有增量未知元的有限差分格式;然后利用非线性Galerkin方法,得到该差分格式的稳定性条件.通过对该格式的稳定性分析,说明和经典的差分格式的稳定性相比较,带有增量未知元的有限差分格式的稳定性得到了提高.     

Finite Volume Multilevel Approximation of the Shallow Water Equations

The authors consider a simple transport equation in one-dimensional space and the linearized shallow water equations in two-dimensional space, and describe and implement a multilevel finite-volume discretization in the context of the utilization of the incremental unknowns. The numerical stability of the method is proved in both cases.     

计算有弥散和吸附的径向渗流问题的局部化间断Galerkin方法

在Corkburn和Shu新近发展的求解对流扩散方程的局部化间断Galerkin方法的基础上,针对有弥散和吸附的径向渗流问题中出现的推广的对流扩散方程的形式,构造了一种计算有弥散和吸附的径向渗流问题的局部化间断Galerkin有限元方法,为径向渗流问题的求解提供了一个高阶的新方法．对对流-弥散和对流-弥散-吸附两种情况进行了数值实验, 所得结果的相应部分与已知的一些精确解结果和数值结果是一致的,表明方法是可靠的．从计算速度上看,方法也是可行的．     

Preconditioning analysis of the one dimensional incremental unknowns method on nonuniform meshes

The condition number of the incremental unknowns matrix on nonuniform meshes associated to the elliptic problem is analyzed. Comparing to the usual nodal unknowns matrix, the condition number of the incremental unknowns matrix is reduced significantly even if the meshes are nonuniform. Furthermore, if a diagonal scaling is used, the condition number of the preconditioned incremental unknowns matrix comes out to be O(1). Numerical experiments are performed respectively on Shishkin mesh and Chebyshev mesh. Computational results with respect to the two particular nonuniform meshes confirm our theoretical analysis.     

非线性对流扩散方程的守恒律

利用直接方法研究了非线性对流扩散方程的守恒律,得到了关于非线性对流扩散方程的守恒律乘子性质的一个定理.利用这个定理,可以简化守恒律乘子的确定方程.随后通过对确定方程中的变量函数进行分析,发现在四种情况下乘子的确定方程是可解的.最后解出这些守恒律乘子,利用积分公式法分别得到了四种情况下对应于各个守恒律乘子的守恒律.     

A semi-circulant preconditioner for the convection-diffusion equation

In this paper we define and analyze a semi-circulant preconditioner for the convection-diffusion equation. We derive analytical formulas for the eigenvalues and the eigenvectors of the preconditioned system of equations. We show that for mesh Péclet numbers less than 2, the rate of convergence depends only on the mesh Péclet number and the direction of the convective field and not on the spatial grid ratio or the number of unknowns. Received February 20, 1997 / Revised version received November 19, 1997     

Analysis of the cell boundary element methods for convection dominated convection-diffusion equations

We introduce two kinds of the cell boundary element (CBE) methods for convection dominated convection-diffusion equations: one is the CBE method with the exact bubble function and the other with inexact bubble functions. The main focus of this paper is on inexact bubble CBE methods. For inexact bubble CBE methods we introduce a family of numerical methods depending on two parameters, one for control of interior layers and the other for outflow boundary layers. Stability and convergence analysis are provided and numerical tests for inexact bubble CBEs with various choices of parameters are presented.     

对流占优扩散方程的一种特征差分算法

A new kind of characteristic-difference scheme for convection-diffusion equations is constructed by characteristic method and bilinear interpolation method. The convergence of the scheme is proved. The advantages of this scheme are to obtain the solutions of the convection diffusion equations with variable coefficient expediently and to reduce the numerical oscillations of the convectiondominanted diffusion equations effectively.