Two-Grid Algorithm of $H^{1}$-Galerkin Mixed Finite Element Methods for Semilinear Parabolic Integro-Differential Equations

In this paper, we present a two-grid discretization scheme for semilinear parabolic integro-differential equations by $H^{1}$-Galerkin mixed finite element methods. We use the lowest order Raviart-Thomas mixed finite elements and continuous linear finite element for spatial discretization, and backward Euler scheme for temporal discretization. Firstly, a priori error estimates and some superclose properties are derived. Secondly, a two-grid scheme is presented and its convergence is discussed. In the proposed two-grid scheme, the solution of the nonlinear system on a fine grid is reduced to the solution of the nonlinear system on a much coarser grid and the solution of two symmetric and positive definite linear algebraic equations on the fine grid and the resulting solution still maintains optimal accuracy. Finally, a numerical experiment is implemented to verify theoretical results of the proposed scheme. The theoretical and numerical results show that the two-grid method achieves the same convergence property as the one-grid method with the choice $h=H^2$.     

Mixed Method for Compressible Miscible Displacement with Dispersion in Porous Media

Compressible miscible displacement of one fluid by another in porous media is modelled by a nonlinear parabolic system. A finite element procedure is introduced to approximate the concentration of one fluid and the pressure of the mixture. The concentration is treated by a Galerkin method while the pressure is treated by a parabolic mixed finite element method. The effect of dispersion, which is neglected in [1], is considered. Optimal order estimates in L2 are derived for the errors in the approximate solutions.     

多孔介质中两相不可压混溶驱替的差分流线扩散法及其误差分析

多孔介质中两相不可压混熔驱替问题可描述为椭圆和抛物耦合的非线性偏微分方程组，对椭圆方程采用混合元方法，而对抛物方程采用差分流线扩散法，本文构造了求解该问题的差分流线扩散-混合元格式，最后，给出所构造格式按L^∞(L^2）模的拟最优误差阶估计。     

A two-grid finite element method for the Allen-Cahn equation with the logarithmic potential

In this paper, we present a two-grid finite element method for the Allen-Cahn equation with the logarithmic potential. This method consists of two steps. In the first step, based on a fully implicit finite element method, the Allen-Cahn equation is solved on a coarse grid with mesh size H. In the second step, a linearized system whose nonlinear term is replaced by the value of the first step is solved on a fine grid with mesh size h. We give the energy stabilities of the traditional finite element method and the two-grid finite element method. The optimal convergence order of the two-grid finite element method in H¹ norm is achieved when the mesh sizes satisfy h = O(H²). Numerical examples are given to demonstrate the validity of the proposed scheme. The results show that the two-grid method can save the CPU time while keeping the same convergence rate.     

多孔介质中可压缩可混溶驱动问题的有限体积元法

有界区域上多孔介质中可压缩可混溶驱动问题由两个非线性抛物型方程耦合而成：压力方程和饱和度方程均是抛物型方程．运用有限体积元法对两个方程进行数值分析,给出了全离散有限体积元格式,并通过详细的理论分析,得到了近似解与原问题真解的最优H^1模误差估计。     

可压溶混流的混合元逼近

二维和三维空间中，多孔介质里可压溶混流被非线性偏微分方程组所描述．浓度方程采用Galerkin方法逼近，而压力方程采用混合有限元逼近．我们导出了浓度、压力、速度及其时间导数的最优L2误差估计，同时得到了浓度和压力的拟最优L∞误差估计．本文处理了带分子弥散的非线性问题．     

IMPLICIT-EXPLICIT MULTISTEP FINITE ELEMENT-MIXED FINITE ELEMENT METHODS FOR THE TRANSIENT BEHAVIOR OF A SEMICONDUCTOR DEVICE

The transient behavior of a semiconductor device consists of a Poisson equa-tion for the electric potential and of two nonlinear parabolic equations for the electrondensity and hole density.The electric potential equation is discretized by a mixed finiteelement method. The electron and hole density equations are treated by implicit-explicitmultistep finite element methods. The schemes are very efficient. The optimal order errorestimates both in time and space are derived.     

Two-Grid Finite Element Method with Crank-Nicolson Fully Discrete Scheme for the Time-Dependent Schrödinger Equation

In this paper, we study the Crank-Nicolson Galerkin finite element method and construct a two-grid algorithm for the general two-dimensional time-dependent Schrödinger equation. Firstly, we analyze the superconvergence error estimate of the finite element solution in $H^1$ norm by use of the elliptic projection operator. Secondly, we propose a fully discrete two-grid finite element algorithm with Crank-Nicolson scheme in time. With this method, the solution of the Schrödinger equation on a fine grid is reduced to the solution of original problem on a much coarser grid together with the solution of two Poisson equations on the fine grid. Finally, we also derive error estimates of the two-grid finite element solution with the exact solution in $H^1$ norm. It is shown that the solution of two-grid algorithm can achieve asymptotically optimal accuracy as long as mesh sizes satisfy $H = \mathcal{O}(h^{\frac{1}{2}})$.     

TWO-GRID CHARACTERISTIC FINITE VOLUME METHODS FOR NONLINEAR PARABOLIC PROBLEMS＊

In this work, two-grid characteristic finite volume schemes for the nonlinear parabolic problem are considered. In our algorithms, the diffusion term is discretized by the finite volume method, while the temporal differentiation and advection terms are treated by the characteristic scheme. Under some conditions about the coefficients and exact solution, optimal error estimates for the numerical solution are obtained. Furthermore, the two- grid characteristic finite volume methods involve solving a nonlinear equation on coarse mesh with mesh size H, a large linear problem for the Oseen two-grid characteristic finite volume method on a fine mesh with mesh size h = O（H2） or a large linear problem for the Newton two-grid characteristic finite volume method on a fine mesh with mesh size h = 0（I log hll/2H3）. These methods we studied provide the same convergence rate as that of the characteristic finite volume method, which involves solving one large nonlinear problem on a fine mesh with mesh size h. Some numerical results are presented to demonstrate the efficiency of the proposed methods.     

A splitting positive definite mixed element method for miscible displacement of compressible flow in porous media

A miscible displacement of one compressible fluid by another in a porous medium is governed by a nonlinear parabolic system. A new mixed finite element method, in which the mixed element system is symmetric positive definite and the flux equation is separated from pressure equation, is introduced to solve the pressure equation of parabolic type, and a standard Galerkin method is used to treat the convection‐diffusion equation of concentration of one of the fluids. The convergence of the approximate solution with an optimal accuracy in L²‐norm is proved. © 2001 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 17: 229–249, 2001     

A PRIORI L_2 ERROR ESTIMATES FOR A NONLINEAR PARABOLIC SYSTEM BASED ON COMBINING THE METHOD OF CHARACTERISTICS WITH MIXED FINITE ELEMENT PROCEDURE

A nonlinear parabolic system is derived to describe incompressible nuclear waste-disposal contamination in porous media. A sequential implicit tirne-stepping is defined, in which the pressure and Darcy velocity of the mixture are approximated simultaneously by a mixed finite element method and the brine, radionuclid and heat are treated by a combination of a Galerkin finite element method and the method of characteristics. Optimal-order convergence in L2 is proved. Time-truncation errors of standard procedures are reduced by time stepping along the characteristics of the hyperbolic part of the brine, radionuclide and heal equalios, temporal and spatial error are lossened by direct compulation of the velocity in the mixed method, as opposed to differentiation of the pressure.     

一类四阶抛物型积分-微分方程的混合间断时空有限元法

构造四阶抛物型积分-微分方程的混合间断时空有限元格式,利用混合有限元方法将高阶方程降阶,利用空间连续而时间允许间断的时空有限元方法离散方程,证明离散解的稳定性,存在唯一性和收敛性.     

带广义边界条件的四阶抛物型方程的混合间断时空有限元法

构造具有广义边界条件的四阶线性抛物型方程的混合间断时空有限元格式,利用混合有限元方法将高阶方程降阶,利用空间连续而时间允许间断的时空有限元方法离散方程,证明了离散解的存在唯一性,稳定性和收敛性,并给出数值算例验证了方法的有效性.     

多孔介质驱动问题的混合元最小二乘特征有限元方法及其收敛分析

１．引言多孔介质二相驱动问题的数学模型是偶合的非线性偏微分方程组的初边值问题．该问题可转化为压力方程和浓度方程［１－４］．浓度方程一般是对流占优的对流扩散方程,它的对流速度依赖于比浓度方程的扩散系数大得多的Ｆａｒｃｙ速度．因此Ｄａｒｃｙ速度的求解精度直接影响着浓度的求解精度．为了提高速度的求解精度,７０年代Ｐ．Ａ．Ｒａｖｉａｔ和Ｊ．Ｍ．Ｔｈｏｍａｓ提出混合有限元方法［５］．Ｊ．ＤｏｕｇｌａｓＪｒ,Ｔ．Ｆ．Ｒｕｓｓｅｌｌ,Ｒ．Ｅ．Ｅｗｉｎｇ,Ｍ．Ｆ．Ｗｈｅｅｌｅｒ［１］－［４］,［９］,［１２］袁…     

A Two-Grid Finite Element Approximation for Nonlinear Time Fractional Two-Term Mixed Sub-Diffusion and Diffusion Wave Equations

In this paper, we develop a two-grid method (TGM) based on the FEM for 2D nonlinear time fractional two-term mixed sub-diffusion and diffusion wave equations. A two-grid algorithm is proposed for solving the nonlinear system, which consists of two steps: a nonlinear FE system is solved on a coarse grid, then the linearized FE system is solved on the fine grid by Newton iteration based on the coarse solution. The fully discrete numerical approximation is analyzed, where the Galerkin finite element method for the space derivatives and the finite difference scheme for the time Caputo derivative with order $\alpha\in(1,2)$ and $\alpha_{1}\in(0,1)$. Numerical stability and optimal error estimate $O(h^{r+1}+H^{2r+2}+\tau^{\min\{3-\alpha,2-\alpha_{1}\}})$ in $L^{2}$-norm are presented for two-grid scheme, where $t,$ $H$ and $h$ are the time step size, coarse grid mesh size and fine grid mesh size, respectively. Finally, numerical experiments are provided to confirm our theoretical results and effectiveness of the proposed algorithm.     

A Two-Grid Algorithm of Fully Discrete Galerkin Finite Element Methods for a Nonlinear Hyperbolic Equation

A two-grid finite element approximation is studied in the fully discrete scheme obtained by discretizing in both space and time for a nonlinear hyperbolic equation. The main idea of two-grid methods is to use a coarse-grid space ($S_H$) to produce a rough approximation for the solution of nonlinear hyperbolic problems and then use it as the initial guess on the fine-grid space ($S_h$). Error estimates are presented in $H^1$-norm, which show that two-grid methods can achieve the optimal convergence order as long as the two different girds satisfy $h$ = $\mathcal{O}$($H^2$). With the proposed techniques, we can obtain the same accuracy as standard finite element methods, and also save lots of time in calculation. Theoretical analyses and numerical examples are presented to confirm the methods.     

APPLICATION OF SUPERCONVERGENCE TO A MODEL FOR COMPRESSIBLE MISCIBLE DISPLACEMENT

A modification of a finite element method of Douglas and Roberts for approximating the solution of the equations describing compressible miscible displacement in a porous medium is proposed and analyzed. The pressure is treated by a parabolic mixed finite element method using a Raviart-Thomas space of index rover a quasiregular partition, An extension of the Darcy velocity along Gauss lines is used in the evaluation of the coefficients in the Galerkin procedure for the concentration. A simple computational procedure allows the superconvergence property of the fluid velocity to be retained in our total algorithm.     

三维热传导型半导体问题的交替方向特征有限元方法及理论分析

１．引言 三维热传导型半导体器件瞬态问题的数学模型由四个非线性偏微分方程描述 ［１,２］．工程研究中一般考虑绝流边条件,由于绝流条件可以看作一反射条件来处理、为了数值分析方便,我们在此考虑三维周期问题： 其中,　＝［０,１］３,未知函数是电子位势　;电子,空穴浓度ｅ,ｐ;温度函数Ｔ．方程（１,１）－（１.４）中出现的系数均有正的上下界,且是　周期的． ａ＝Ｑ／ε,Ｑ,ε分别表示电子负荷和介电系数,均为正常数．Ｎ（ｘ）是给定的函数．Ｄｓ（ｘ）为扩散系数,μｓ（ｘ）为迁移率,ｓ＝ｅ,Ｐ．Ｒ（ｅ,ｐ,Ｔ）…     

Partition of Unity for a Class of Nonlinear Parabolic Equation on Overlapping Non-Matching Grids

A class of nonlinear parabolic equation on a polygonal domain Ω  R2 is inves- tigated in this paper. We introduce a finite element method on overlapping non-matching grids for the nonlinear parabolic equation based on the partition of unity method. We give the construction and convergence analysis for the semi-discrete and the fully discrete finite element methods. Moreover, we prove that the error of the discrete variational problem has good approximation properties. Our results are valid for any spatial dimensions. A numerical example to illustrate the theoretical results is also given.     

基于二重网格的定常Navier-Stokes方程的局部和并行有限元算法

对二维定常的不可压缩的Navier-Stokes方程的局部和并行算法进行了研究．给出的算法是多重网格和区域分解相结合的算法,它是基于两个有限元空间:粗网格上的函数空间和子区域的细网格上的函数空间．局部算法是在粗网格上求一个非线性问题,然后在细网格上求一个线性问题,并舍掉内部边界附近的误差相对较大的解．最后,基于局部算法,通过有重叠的区域分解而构造了并行算法,并且做了算法的误差分析,得到了比标准有限元方法更好的误差估计,也对算法做了数值试验,数值结果通过比较验证了本算法的高效性和合理性．