The method of fundamental solutions with eigenfunction expansion method for nonhomogeneous diffusion equation

In this article we describe a numerical method to solve a nonhomogeneous diffusion equation with arbitrary geometry by combining the method of fundamental solutions (MFS), the method of particular solutions (MPS), and the eigenfunction expansion method (EEM). This forms a meshless numerical scheme of the MFS‐MPS‐EEM model to solve nonhomogeneous diffusion equations with time‐independent source terms and boundary conditions for any time and any shape. Nonhomogeneous diffusion equation with complex domain can be separated into a Poisson equation and a homogeneous diffusion equation using this model. The Poisson equation is solved by the MFS‐MPS model, in which the compactly supported radial basis functions are adopted for the MPS. On the other hand, utilizing the EEM the diffusion equation is first translated to a Helmholtz equation, which is then solved by the MFS together with the technique of the singular value decomposition (SVD). Since the present meshless method does not need mesh generation, nodal connectivity, or numerical integration, the computational effort and memory storage required are minimal as compared with other numerical schemes. Test results for two 2D diffusion problems show good comparability with the analytical solutions. The proposed algorithm is then extended to solve a problem with irregular domain and the results compare very well with solutions of a finite element scheme. Therefore, the present scheme has been proved to be very promising as a meshfree numerical method to solve nonhomogeneous diffusion equations with time‐independent source terms of any time frame, and for any arbitrary geometry. © 2006 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2006     

Numerical approximations for the nonlinear time fractional reaction–diffusion equation

In this article, we propose an implicit pseudospectral scheme for nonlinear time fractional reaction–diffusion equations with Neumann boundary conditions, which is based upon Gauss–Lobatto–Legendre–Birkhoff pseudospectral method in space and finite difference method in time. A priori estimate of numerical solution is given firstly. Then the existence of numerical solution is proved by Brouwer fixed point theorem and the uniqueness is obtained. It is proved rigorously that the fully discrete scheme is unconditionally stable and convergent. Furthermore, we develop a modified scheme by adding correction terms for the problem with nonsmooth solutions. Numerical examples are given to verify the theoretical analysis.     

Convergence of Three-Step Taylor Galerkin Finite Element Scheme Based Monotone Schwarz Iterative Method for Singularly Perturbed Differential-Difference Equation

Monotone Schwarz iterative methods for parabolic partial differential equations are well known for their advantage of eliminating the search for an initial solution. In this article, we propose a monotone Schwarz iterative method for singularly perturbed parabolic retarded differential-difference equations based on a three-step Taylor Galerkin finite element scheme. The stability and ε-uniform convergence of the three-step Taylor Galerkin finite element method have been discussed. Further, by using maximum principle and induction hypothesis, the convergence of the proposed monotone Schwarz iterative method has been established.     

A high order finite difference method with Richardsonextrapolation for 3D convection diffusion equation

In this paper, we extend the Sun and Zhang’s [24] work on high order finite difference method, which is based on the Richardson extrapolation technique and an operator interpolation scheme for the one and two dimensional steady convection diffusion equations to the three dimensional case. Firstly, we employ a fourth order compact difference scheme to get the fourth order accurate solution on the fine and the coarse grids. Then, we use the Richardson extrapolation technique by combining the two approximate solutions to get a sixth order accurate solution on coarse grid. Finally, we apply an operator interpolation scheme to achieve the sixth order accurate solution on the fine grid. During this process, we use alternating direction implicit (ADI) method to solve the resulting linear systems. Numerical experiments are conducted to verify the accuracy and effectiveness of the present method.     

Tikhonov regularization method for a backward problem for the inhomogeneous time-fractional diffusion equation

Fractional (nonlocal) diffusion equations replace the integer-order derivatives in space and time by their fractional-order analogs and they are used to model anomalous diffusion, especially in physics. In this paper, we study a backward problem for an inhomogeneous time-fractional diffusion equation with variable coefficients in a general bounded domain. Such a backward problem is of practically great importance because we often do not know the initial density of substance, but we can observe the density at a positive moment. The backward problem is ill-posed and we propose a regularizing scheme by using Tikhonov regularization method. We also prove the convergence rate for the regularized solution by using an a priori regularization parameter choice rule. Numerical examples illustrate applicability and high accuracy of the proposed method.     

A combined finite volume–finite element scheme for the discretization of strongly nonlinear convection–diffusion–reaction problems on nonmatching grids

We propose and analyze in this paper a numerical scheme for nonlinear degenerate parabolic convection–diffusion–reaction equations in two or three space dimensions. We discretize the time evolution, convection, reaction, and source terms on a given grid, which can be nonmatching and can contain nonconvex elements, by means of the cell‐centered finite volume method. To discretize the diffusion term, we construct a conforming simplicial mesh with the vertices given by the original grid and use the conforming piecewise linear finite element method. In this way, the scheme is fully consistent and the discrete solution is naturally continuous across the interfaces between the subdomains with nonmatching grids, without introducing any supplementary equations and unknowns or using any interpolation at the interfaces. We allow for general inhomogeneous and anisotropic diffusion–dispersion tensors, propose two variants corresponding respectively to arithmetic and harmonic averaging, and use the local Péclet upstream weighting in order to only add the minimal numerical diffusion necessary to avoid spurious oscillations in the convection‐dominated case. The scheme is robust, efficient since it leads to positive definite matrices and one unknown per element, locally conservative, and satisfies the discrete maximum principle under the conditions on the simplicial mesh and the diffusion tensor usual in the finite element method. We prove its convergence using a priori estimates and the Kolmogorov relative compactness theorem and illustrate its behavior on a numerical experiment. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2010     

NONLINEAR GALERKIN METHOD AND CRANK-NICOLSONMETHOD FOR VISCOUS INCOMPRESSIBLE FLOW

1.Introduction'NonlinearGalerkinmethodisnumericalmethodfordissipativeevolutionpartialdifferentialequationswherethespatialdiscretizationreliesonanonlinearmanifoldinsteadofalinearspaceasintheclassicalGalerkinmethod.Morepreciselygoneconsidersafinitedimension…     

A Variational Analysis for the Moving Finite Element Method for Gradient Flows

By using the Onsager principle as an approximation tool, we give a novel derivation for the moving finite element method for gradient flow equations. We show that the discretized problem has the same energy dissipation structure as the continuous one. This enables us to do numerical analysis for the stationary solution of a nonlinear reaction diffusion equation using the approximation theory of free-knot piecewise polynomials. We show that under certain conditions the solution obtained by the moving finite element method converges to a local minimizer of the total energy when time goes to infinity. The global minimizer, once it is detected by the discrete scheme, approximates the continuous stationary solution in optimal order. Numerical examples for a linear diffusion equation and a nonlinear Allen-Cahn equation are given to verify the analytical results.     

An upwind-mixed method on changing meshes for two-phase miscible flow in porous media

Two-phase miscible flow in porous media is governed by a system of nonlinear partial differential equations. In this paper, the upwind-mixed method on dynamically changing meshes is presented for the problem in two dimensions. The pressure is approximated by a mixed finite element method and the concentration by a method which upwinds the convection and incorporates diffusion using an expanded mixed finite element method. The method developed is shown to obtain almost optimal rate error estimate. When the method is modified we can obtain the optimal rate error estimate that is well known for static meshes. The modification of the scheme is the construction of a linear approximation to the solution, which is used in projecting the solution from one mesh to another. Finally, numerical experiments are given.     

非定常不可压Navier-Stokes方程两重网格方法的稳定性和L2误差估计

讨论了二维非定常不可压Navier-Stokes方程的两重网格方法．此方法包括在粗网格上求解一个非线性问题,在细网格上求解一个Stokes问题．采用一种新的全离散(时间离散用Crank-Nicolson格式,空间离散用混合有限元方法)格式数值求解N-S方程．证明了该全离散格式的稳定性．给出了L2误差估计．对比标准有限元方法,在保持同样精度的前提下,TGM能节省大量的计算量．     

对流占优扩散方程的特征线楔形基无网格法

本文研究了一维对流占优扩散方程的初边值问题.利用特征线法与楔形基无网格法,获得了特征线楔形基无网格显格式与隐格式算法.数值实验表明算法具有精度高、计算简单等优点.     

A uniformly convergent scheme for a system of reaction–diffusion equations

In this work a system of two parabolic singularly perturbed equations of reaction–diffusion type is considered. The asymptotic behaviour of the solution and its partial derivatives is given. A decomposition of the solution in its regular and singular parts has been used for the asymptotic analysis of the spatial derivatives. To approximate the solution we consider the implicit Euler method for time stepping and the central difference scheme for spatial discretization on a special piecewise uniform Shishkin mesh. We prove that this scheme is uniformly convergent, with respect to the diffusion parameters, having first-order convergence in time and almost second-order convergence in space, in the discrete maximum norm. Numerical experiments illustrate the order of convergence proved theoretically.     

Numerical solutions for some coupled systems of nonlinear boundary value problems

Summary In the well-known Volterra-Lotka model concerning two competing species with diffusion, the densities of the species are governed by a coupled system of reaction diffusion equations. The aim of this paper is to present an iterative scheme for the steady state solutions of a finite difference system which corresponds to the coupled nonlinear boundary value problems. This iterative scheme is based on the method of upper-lower solutions which leads to two monotone sequences from some uncoupled linear systems. It is shown that each of the two sequences converges to a nontrivial solution of the discrete equations. The model under consideration may have one, two or three nonzero solutions and each of these solutions can be computed by a suitable choice of initial iteration. Numerical results are given for these solutions under both the Dirichlet boundary condition and the mixed type boundary condition.     

Accelerated monotone scheme for finite difference equations concerning steady-state prey–predator interactions

In previous work, by adapting a suitable finite difference method to a particular monotone scheme, the authors and A. Lazer have studied the numerical solution of a system of semilinear elliptic partial differential equations which determines the equilibria of the Volterra–Lotka equations describing prey–predator interactions with diffusion. In this paper, in order to improve the efficiency of the method, we show how Newton's method can be successfully combined with the previous scheme to greatly accelerate the convergence. In some particularly ‘difficult’ problems, the new method reduces the average number of iterations necessary to generate each element of the monotone sequences from 15 to about 3.     

On the maximum principle and its application to diffusion equations

In this article, an analog of the maximum principle has been established for an ordinary differential operator associated with a semi‐discrete approximation of parabolic equations. In applications, the maximum principle is used to prove O(h²) and O(h⁴) uniform convergence of the method of lines for the diffusion Equation (1). The system of ordinary differential equations obtained by the method of lines is solved by an implicit predictor corrector method. The method is tested by examples with the use of the enclosed Mathematica module solveDiffusion. The module solveDiffusion gives the solution by O(h²) uniformly convergent discrete scheme or by O(h⁴) uniformly convergent discrete scheme. © 2006 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2007     

A numerical scheme based on a solution of nonlinear advection–diffusion equations

A mathematical formulation of the two-dimensional Cole–Hopf transformation is investigated in detail. By making use of the Cole–Hopf transformation, a nonlinear two-dimensional unsteady advection–diffusion equation is transformed into a linear equation, and the transformed equation is solved by the spectral method previously proposed by one of the authors. Thus a solution to initial value problems of nonlinear two-dimensional unsteady advection–diffusion equations is derived. On the base of the solution, a numerical scheme explicit with respect to time is presented for nonlinear advection–diffusion equations. Numerical experiments show that the present scheme possesses the total variation diminishing properties and gives solutions with good quality.     

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

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.     

Convergence analysis of a DDFV scheme for a system describing miscible fluid flows in porous media

In this article, we prove the convergence of a discrete duality finite volume scheme for a system of partial differential equations describing miscible displacement in porous media. This system is made of two coupled equations: an anisotropic diffusion equation on the pressure and a convection‐diffusion‐dispersion equation on the concentration. We first establish some a priori estimates satisfied by the sequences of approximate solutions. Then, it yields the compactness of these sequences. Passing to the limit in the numerical scheme, we finally obtain that the limit of the sequence of approximate solutions is a weak solution to the problem under study. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 723–760, 2015     

Wavelet-Taylor Galerkin Method for the Burgers Equation

In this paper, we propose a wavelet-Taylor Galerkin method for the numerical solution of the Burgers equation. In deriving the computational scheme, Taylor-generalized Euler time discretization is performed prior to wavelet-based Galerkin spatial approximation. The linear system of equations obtained in the process are solved by approximate-factorization-based simple explicit schemes, and the resulting solution is compared with that from regular methods. To deal with transient advection-diffusion situations that evolve toward a convective steady state, a splitting-up strategy is known to be very effective. So the Burgers equation is also solved by a splitting-up method using a wavelet-Taylor Galerkin approach. Here, the advection and diffusion terms in the Burgers equation are separated, and the solution is computed in two phases by appropriate wavelet-Taylor Galerkin schemes. Asymptotic stability of all the proposed schemes is verified, and the L_∞ errors relative to the analytical solution together with the numerical solution are reported. AMS subject classification (2000) 65M70     

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

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