An Eulerian‐Lagrangian discontinuous Galerkin method for transient advection‐diffusion equations

We develop an Eulerian‐Lagrangian discontinuous Galerkin method for time‐dependent advection‐diffusion equations. The derived scheme has combined advantages of Eulerian‐Lagrangian methods and discontinuous Galerkin methods. The scheme does not contain any undetermined problem‐dependent parameter. An optimal‐order error estimate and superconvergence estimate is derived. Numerical experiments are presented, which verify the theoretical estimates.© 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2007     

A uniformly optimal‐order error estimate for a bilinear finite element method for degenerate convection‐diffusion equations

We prove an optimal‐order error estimate in a degenerate‐diffusion weighted energy norm for bilinear Galerkin finite element methods for two‐dimensional time‐dependent convection‐diffusion equations with degenerate diffusion. In the estimate, the generic constants depend only on certain Sobolev norms of the true solution but not the lower bound of the diffusion. This estimate, combined with a known stability estimate of the true solution of the governing partial differential equations, yields an optimal‐order estimate of the Galerkin finite element method, in which the generic constants depend only on the Sobolev norms of the initial and right side data. Preliminary numerical experiments were conducted to verify these estimates numerically. © 2011 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2011     

A uniformly optimal‐order estimate for finite volume method for advection‐diffusion equations

We prove an optimal‐order error estimate in a weighted energy norm for finite volume method for two‐dimensional time‐dependent advection–diffusion equations on a uniform space‐time partition of the domain. The generic constants in the estimates depend only on certain norms of the true solution but not on the scaling parameter. These estimates, combined with a priori stability estimates of the governing partial differential equations with full regularity, yield a uniform estimate of the finite volume method, in which the generic constants depend only on the Sobolev norms of the initial and right side data but not on the scaling parameter. We use the interpolation of spaces and stability estimates to derive a uniform estimate for problems with minimal or intermediate regularity, where the convergence rates are proportional to certain Besov norms of the initial and right‐hand side data. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 17‐43, 2014     

A robust algebraic multilevel preconditioner for non‐symmetric M‐matrices

Stable finite difference approximations of convection‐diffusion equations lead to large sparse linear systems of equations whose coefficient matrix is an M‐matrix, which is highly non‐symmetric when the convection dominates. For an efficient iterative solution of such systems, it is proposed to consider in the non‐symmetric case an algebraic multilevel preconditioning method formerly proposed for pure diffusion problems, and for which theoretical results prove grid independent convergence in this context. These results are supplemented here by a Fourier analysis that applies to constant coefficient problems with periodic boundary conditions whenever using an ‘idealized’ version of the two‐level preconditioner. Within this setting, it is proved that any eigenvalue λ of the preconditioned system satisfies for some real constant c such that . This result holds independently of the grid size and uniformly with respect to the ratio between convection and diffusion. Extensive numerical experiments are conducted to assess the convergence of practical two‐ and multi‐level schemes. These experiments, which include problems with highly variable and rotating convective flow, indicate that the convergence is grid independent. It deteriorates moderately as the convection becomes increasingly dominating, but the convergence factor remains uniformly bounded. This conclusion is supported for both uniform and some non‐uniform (stretched) grids. Copyright © 2000 John Wiley & Sons, Ltd.     

An Eulerian‐Lagrangian control volume scheme for two‐dimensional unsteady advection‐diffusion problems

We develop a mass conservative Eulerian‐Lagrangian control volume scheme (ELCVS) for the solution of the transient advection‐diffusion equations in two space dimensions. This method uses finite volume test functions over the space‐time domain defined by the characteristics within the framework of the class of Eulerian‐Lagrangian localized adjoint characteristic methods (ELLAM). It, therefore, maintains the advantages of characteristic methods in general, and of this class in particular, which include global mass conservation as well as a natural treatment of all types of boundary conditions. However, it differs from other methods in that class in the treatment of the mass storage integrals at the previous time step defined on deformed Lagrangian regions. This treatment is especially attractive for orthogonal rectangular Eulerian grids composed of block elements. In the algorithm, each deformed region is approximated by an eight‐node region with sides drawn parallel to the Eulerian grid, which significantly simplifies the integration compared with the approach used in finite volume ELLAM methods, based on backward tracking, while retaining local mass conservation at no additional expenses in terms of accuracy or CPU consumption. This is verified by numerical tests which show that ELCVS performs as well as standard finite volume ELLAM methods, which have previously been shown to outperform many other well‐received classes of numerical methods for the equations considered. © 2011 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2012     

Taylor‐Galerkin B‐spline finite element method for the one‐dimensional advection‐diffusion equation

The advection‐diffusion equation has a long history as a benchmark for numerical methods. Taylor‐Galerkin methods are used together with the type of splines known as B‐splines to construct the approximation functions over the finite elements for the solution of time‐dependent advection‐diffusion problems. If advection dominates over diffusion, the numerical solution is difficult especially if boundary layers are to be resolved. Known test problems have been studied to demonstrate the accuracy of the method. Numerical results show the behavior of the method with emphasis on treatment of boundary conditions. Taylor‐Galerkin methods have been constructed by using both linear and quadratic B‐spline shape functions. Results shown by the method are found to be in good agreement with the exact solution. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2010     

A uniform estimate for the MMOC for two‐dimensional advection‐diffusion equations

We prove an optimal‐order error estimate in a weighted energy norm for the modified method of characteristics (MMOC) and the modified method of characteristics with adjusted advection (MMOCAA) for two‐dimensional time‐dependent advection‐diffusion equations, in the sense that the generic constants in the estimates depend on certain Sobolev norms of the true solution but not on the scaling diffusion parameter ε. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2010     

Two‐grid methods for mixed finite‐element solution of coupled reaction‐diffusion systems

We develop 2‐grid schemes for solving nonlinear reaction‐diffusion systems: where p = (p, q) is an unknown vector‐valued function. The schemes use discretizations based on a mixed finite‐element method. The 2‐grid approach yields iterative procedures for solving the nonlinear discrete equations. The idea is to relegate all the Newton‐like iterations to grids much coarser than the final one, with no loss in order of accuracy. The iterative algorithms examined here extend a method developed earlier for single reaction‐diffusion equations. An application to prepattern formation in mathematical biology illustrates the method's effectiveness. © 1999 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 15: 589–604, 1999     

Coexistence states for a Lotka‐Volterra symbiotic system with cross‐diffusion

In this work, we study coexistence states for a Lotka‐Volterra symbiotic system with cross‐diffusion under homogeneous Dirichlet boundary conditions. By using topological degree theory and bifurcation theory, we prove the existence and multiplicity of positive solutions under certain conditions on the parameters. Asymptotic behaviors of positive solutions are respectively studied as the cross‐diffusion coefficient tends to infinity and the interaction rate tends to zero. Finally, we compare our results with those of the Lotka‐Volterra predator and competition systems.     

An Eulerian‐Lagrangian substructuring domain decomposition method for unsteady‐state advection‐diffusion equations

We develop an Eulerian‐Lagrangian substructuring domain decomposition method for the solution of unsteady‐state advection‐diffusion transport equations. This method reduces to an Eulerian‐Lagrangian scheme within each subdomain and to a type of Dirichlet‐Neumann algorithm at subdomain interfaces. The method generates accurate and stable solutions that are free of artifacts even if large time‐steps are used in the simulation. Numerical experiments are presented to show the strong potential of the method. © 2001 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 17:565–583, 2001     

Diffusion in poro‐plastic media

A model is developed for the flow of a slightly compressible fluid through a saturated inelastic porous medium. The initial‐boundary‐value problem is a system that consists of the diffusion equation for the fluid coupled to the momentum equation for the porous solid together with a constitutive law which includes a possibly hysteretic relation of elasto‐visco‐plastic type. The variational form of this problem in Hilbert space is a non‐linear evolution equation for which the existence and uniqueness of a global strong solution is proved by means of monotonicity methods. Various degenerate situations are permitted, such as incompressible fluid, negligible porosity, or a quasi‐static momentum equation. The essential sufficient conditions for the well‐posedness of the system consist of an ellipticity condition on the term for diffusion of fluid and either a viscous or a hardening assumption in the constitutive relation for the porous solid. Copyright © 2004 John Wiley & Sons, Ltd.     

Gauss‐Lobatto‐Legendre‐Birkhoff pseudospectral approximations for the multi‐term time fractional diffusion‐wave equation with Neumann boundaryconditions

A second‐order finite difference/pseudospectral scheme is proposed for numerical approximation of multi‐term time fractional diffusion‐wave equation with Neumann boundary conditions. The scheme is based upon the weighted and shifted Grünwald difference operators approximation of the time fractional calculus and Gauss‐Lobatto‐Legendre‐Birkhoff (GLLB) pseudospectral method for spatial discretization. The unconditionally stability and convergence of the scheme are rigorously proved. Numerical examples are carried out to verify theoretical results.     

Nonexistence of coexisting steady‐state solutions of Dirichlet problem for a cross‐diffusion model

We study a special case of Shigesada–Kawasaki–Teramoto (SKT) model for two competing species with the Dirichlet boundary conditions. In our case, one of the species is not influenced by self‐diffusion or cross‐diffusion. We specify the explicit range of parameters by contradiction such that there are no coexisting steady‐state solutions to the model.     

A single‐node characteristic collocation method for unsteady‐state convection‐diffusion equations in three‐dimensional spaces

We develop a nonconventional single‐node characteristic collocation method with piecewise‐cubic Hermite polynomials for the numerical simulation to unsteady‐state advection‐diffusion transport partial differential equations. This method greatly reduces the number of unknowns in the conventional collocation method, and generates accurate numerical solutions even if very large time steps are taken. The reduction of number of nodes has great potential for problems defined on high space dimensions, which appears in such problems as quantification of uncertainties in subsurface porous media. The method developed here is easy to formulate. Numerical experiments are presented to show the strong potential of the method. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 27: 786–802, 2011     

Extension of high‐order compact schemes to time‐dependent problems

In this article, we present an extension of our previous approaches for steady‐state higher‐order compact (HOC) difference methods to time‐dependent problems. The formulation also provides a framework for similar treatment of other HOC spatial schemes. A stability analysis is provided for transient convection‐diffusion in 1D and transient diffusion in 2D. Supporting numerical experiments are included to illustrate stability and accuracy as well as oscillatory and dissipative behavior. © 2001 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 17: 657–672, 2001     

A high‐order B‐spline collocation scheme for solving a nonhomogeneous time‐fractional diffusion equation

A high‐accuracy numerical approach for a nonhomogeneous time‐fractional diffusion equation with Neumann and Dirichlet boundary conditions is described in this paper. The time‐fractional derivative is described in the sense of Riemann‐Liouville and discretized by the backward Euler scheme. A fourth‐order optimal cubic B‐spline collocation (OCBSC) method is used to discretize the space variable. The stability analysis with respect to time discretization is carried out, and it is shown that the method is unconditionally stable. Convergence analysis of the method is performed. Two numerical examples are considered to demonstrate the performance of the method and validate the theoretical results. It is shown that the proposed method is of order O(Δx⁴ + Δt^{2 ? α}) convergence, where α ∈ (0,1) . Moreover, the impact of fractional‐order derivative on the solution profile is investigated. Numerical results obtained by the present method are compared with those obtained by the method based on standard cubic B‐spline collocation method. The CPU time for present numerical method and the method based on cubic B‐spline collocation method are provided.     

A CFL‐free explicit characteristic interior penalty scheme for linear advection‐reaction equations

We develop a CFL‐free, explicit characteristic interior penalty scheme (CHIPS) for one‐dimensional first‐order advection‐reaction equations by combining a Eulerian‐Lagrangian approach with a discontinuous Galerkin framework. The CHIPS method retains the numerical advantages of the discontinuous Galerkin methods as well as characteristic methods. An optimal‐order error estimate in the L² norm for the CHIPS method is derived and numerical experiments are presented to confirm the theoretical estimates. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2010     

Relaxation‐time limit of the three‐dimensional hydrodynamic model with boundary effects

In this paper, we study three‐dimensional (3D) unipolar and bipolar hydrodynamic models and corresponding drift‐diffusion models from semiconductor devices on bounded domain. Based on the asymptotic behavior of the solutions to the initial boundary value problems with slip boundary condition, we investigate the relation between the 3D hydrodynamic semiconductor models and the corresponding drift‐diffusion models. That is, we discuss the relation‐time limit from the 3D hydrodynamic semiconductor models to the corresponding drift‐diffusion models by comparing the large‐time behavior of these two models. These results can be showed by energy arguments. Copyrightcopyright 2011 John Wiley & Sons, Ltd.     

Graded Galerkin methods for the high‐order convection‐diffusion problem

We develop a Galerkin method using the Hermite spline on an admissible graded mesh for solving the high‐order singular perturbation problem of the convection‐diffusion type. We identify a special function class to which the solution of the convection‐diffusion problem belongs and characterize the approximation order of the Hermite spline for such a function class. The approximation order is then used to establish the optimal order of uniform convergence for the Galerkin method. Numerical results are presented to confirm the theoretical estimate.© 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2009     

Blow‐up solutions for localized reaction‐diffusion equations with variable exponents

This paper deals with radial solutions to localized reaction‐diffusion equations with variable exponents, subject to homogeneous Dirichlet boundary conditions. The global existence versus blow‐up criteria are studied in terms of the variable exponents. We proposed that the maximums of variable exponents are the key clue to determine blow‐up classifications and describe blow‐up rates for positive solutions. Copyright © 2011 John Wiley & Sons, Ltd.