平面热传导方程Douglas交替方向隐格式的稳定性与收敛性

1 引言 解抛物方程的交替方向差分方法产生于本世纪五十年代,由于其稳定性好、易编制计算程序,几十年来在科学与工程计算中得到了广泛的应用,最早的交替方向差分格式由Peace- man,D.W.和Rachford,H.H.以及Douglas,J,Jr.和Rachford,H.H.给出,Douglas,J.Jr.于1962年提出了后人称之为Douglas格式的一类交替方向隐格式,其特点是有二阶截断误差而且适用于三维情形.下面简要给出其格式。 考虑平面方形区域上热传导方程 (1.1) (1.2) (1.3)其中Ω=(0,1)×(0,1),为其边界,t为时间变量,以及u_o为已知函数. 令x与y方向网格步长均为h,x_i=ih,y_j=jh,N=1/h为正整数;时间步长为△t,t~n=n△t,记w_(ij)~n=w(x_i,y_j,t~n)以及     

A parameter-uniform numerical method for a Sobolev problem with initial layer

The present study is concerned with the numerical solution, using finite difference method of a one-dimensional initial-boundary value problem for a linear Sobolev or pseudo-parabolic equation with initial jump. In order to obtain an efficient method, to provide good approximations with independence of the perturbation parameter, we have developed a numerical method which combines a finite difference spatial discretization on uniform mesh and the implicit rule on Shishkin mesh(S-mesh) for the time variable. The fully discrete scheme is shown to be convergent of order two in space and of order one expect for a logarithmic factor in time, uniformly in the singular perturbation parameter. Some numerical results confirming the expected behavior of the method are shown.      

A globally well-posed finite element algorithm for aerodynamics applications

A finite element CFD algorithm is developed for Euler and Navier-Stokes aerodynamic applications. For the linear basis, the resultant approximation is at least second-order-accurate in time and space for synergistic use of three procedures: (1) a Taylor weak statement, which provides for derivation of companion conservation law systems with embedded dispersion-error control mechanisms; (2) a stiffly stable second-order-accurate implicit Rosenbrock-Runge-Kutta temporal algorithm; and (3) a matrix tensor product factorization that permits efficient numerical linear algebra handling of the terminal large-matrix statement. Thorough analyses are presented regarding well-posed boundary conditions for inviscid and viscous flow specifications. Numerical solutions are generated and compared for critical evaluation of quasi-one- and two-dimensional Euler and Navier-Stokes benchmark test problems. Of critical importance, essentially non-oscillatory solutions are uniformly attained for a range of supercritical flow situations with shocks.     

EXISTENCE AND UNIQUENESS OF THE ENTROPY SOLUTION TO A NONLINEAR HYPERBOLIC EQUATION

ＥＸＩＳＴＥＮＣＥＡＮＤＵＮＩＱＵＥＮＥＳＳＯＦＴＨＥＥＮＴＲＯＰＹＳＯＬＵＴＩＯＮＴＯＡＮＯＮＬＩＮＥＡＲＨＹＰＥＲＢＯＬＩＣＥＱＵＡＴＩＯＮ￥Ｒ．ＥＹＭＡＲＤ；Ｔ．ＧＡＬＬＯＵＥＴ；Ｒ．ＨＥＲＢＩＮ（ＬａｂｏｒａｔｏｉｒｅＣｅｎｔｒａｌｄｅｓＰｏ...     

On streamline diffusion arising in Galerkin FEM with predictor/multi‐corrector time integration

In the present paper, the author shows that the predictor/multi‐corrector (PMC) time integration for the advection–diffusion equations induces numerical diffusivity acting only in the streamline direction, even though the equations are spatially discretized by the conventional Galerkin finite element method (GFEM). The transient 2‐D and 3‐D advection problems are solved with the PMC scheme using both the GFEM and the streamline upwind/Petrov Galerkin (SUPG) as the spatial discretization methods for comparison. The solutions of the SUPG‐PMC turned out to be overly diffusive due to the additional PMC streamline diffusion, while the solutions of the GFEM‐PMC were comparatively accurate without significant damping and phase error. A similar tendency was seen also in the quasi‐steady solutions to the incompressible viscous flow problems: 2‐D driven cavity flow and natural convection in a square cavity. Copyright © 2002 John Wiley & Sons, Ltd.     

Banach空间中的相补问题

本文在Ｂａｎａｃｈ空间中研究了三类相补问题解的存在性。所得结果是［４，５，６，９，１１－１４］中相应结果的深入和发展。     

A PARALLEL COMPUTATION SCHEME FOR IMPLICIT RUNGE-KUTTA METHODS AND THE ITERATIVELY B-CONVERGENCE OF ITS NEWTON ITERATIVE PROCESS

In this paper, based on the implicit Runge-Kutta(IRK) methods, we derive a class of parallel scheme that can be implemented on the parallel computers with Ns(N is a positive even number) processors efficiently, and discuss the iteratively B-convergence of the Newton iterative process for solving the algebraic equations of the scheme, secondly we present a strategy providing initial values parallelly for the iterative process. Finally, some numerical results show that our parallel scheme is higher efficient as N is not so large.     

三维热传导型半导体器件瞬态模拟问题的Crank-Nicolson差分-流线扩散有限元法及其数值分析

本文研究三维热传导型半导体器件瞬态模拟问题的数值方法.针对数学模型中各方程不同的特点,分别提出不同的有限元格式.特别针对浓度方程组是对流为主扩散问题的特点,使用Crank-Nicolson差分-流线扩散计算格式,提高了数值解的稳定性.得到的L2误差估计关于空间剖分步长是拟最优的,关于时间步长具有二阶精度.     

二维三温热传导方程组的分数步隐式差分格式

本文给出一个解二维三温热传导方程组的分数步隐式有限差分格式，利用离散变分形式及能量方法，给出差分格式的最优阶离散H^1范数先验误差及稳定性估计。