The Sinc–Legendre collocation method for a class of fractional convection–diffusion equations with variable coefficients

This paper deals with the numerical solution of classes of fractional convection–diffusion equations with variable coefficients. The fractional derivatives are described based on the Caputo sense. Our approach is based on the collocation techniques. The method consists of reducing the problem to the solution of linear algebraic equations by expanding the required approximate solution as the elements of shifted Legendre polynomials in time and the Sinc functions in space with unknown coefficients. The properties of Sinc functions and shifted Legendre polynomials are then utilized to evaluate the unknown coefficients. Several examples are given and the numerical results are shown to demonstrate the efficiency of the newly proposed method.     

Superconvergence of finite element method for singularly perturbed convection–diffusion equations in 1D

In this article, a singularly perturbed convection–diffusion equation is solved by a linear finite element method on a Shishkin mesh. By means of an analysis exploiting symmetries in the convective term of the bilinear form, a new superconvergence rate, which improves the existing result, is obtained.     

Numerical solutions for convection–diffusion equation using El-Gendi method

The method of El-Gendi [El-Gendi SE. Chebyshev solution of differential integral and integro-differential equations. J Comput 1969;12:282–7; Mihaila B, Mihaila I. Numerical approximation using Chebyshev polynomial expansions: El-gendi’s method revisited. J Phys A Math Gen 2002;35:731–46] is presented with interface points to deal with linear and non-linear convection–diffusion equations.The linear problem is reduced to two systems of ordinary differential equations. And, then, each system is solved using three-level time scheme.The non-linear problem is reduced to three systems of ordinary differential. Each one of these systems is, then, solved using three-level time scheme. Numerical results for Burgers’ equation and modified Burgers’ equation are shown and compared with other methods. The numerical results are found to be in good agreement with the exact solutions.     

Exp-function method for solving Kuramoto–Sivashinsky and Boussinesq equations

In this paper, we use the Exp-function method to construct the generalized solitary and periodic solution of the Kuramoto–Sivashinsky and Boussinesq equations. These equations play very important role in mathematical physics and engineering sciences. The suggested algorithm is quite efficient and is practically well suited for use in these problems. The results show the reliability and efficiency of the proposed method.     

A study of isogeometric analysis for scalar convection–diffusion equations

Isogeometric analysis (IGA), in combination with the streamline upwind Petrov–Galerkin (SUPG) stabilization, is studied for the discretization of steady-state convection–diffusion equations. Numerical results obtained for the Hemker problem are compared with results computed with the SUPG finite element method of the same order. Using an appropriate parameterization for IGA, the computed solutions are much more accurate than those obtained with the finite element method, both in terms of the size of spurious oscillations and of the sharpness of layers.     

Application of the Laplace decomposition method for solving linear and nonlinear fractional diffusion–wave equations

In this paper, the Laplace decomposition method is employed to obtain approximate analytical solutions of the linear and nonlinear fractional diffusion–wave equations. This method is a combined form of the Laplace transform method and the Adomian decomposition method. The proposed scheme finds the solutions without any discretization or restrictive assumptions and is free from round-off errors and therefore, reduces the numerical computations to a great extent. The fractional derivative described here is in the Caputo sense. Some illustrative examples are presented and the results show that the solutions obtained by using this technique have close agreement with series solutions obtained with the help of the Adomian decomposition method.     

A skew-symmetric iterative method for solving a steady convection–diffusion–reaction equation with an indefinite reaction coefficient

An iterative product-type triangular skew-symmetric method (PTSM) is used to solve systems of linear algebraic equations (SLAEs) obtained by approximation with a central-difference scheme of a first-type boundary value problem for convection–diffusion–reaction and standard grid ordering. Sufficient conditions for non-negative definiteness of the SLAE matrix resulting from this approximation are obtained for the indefinite reaction coefficient. This property provides convergence of a wide class of iterative methods (in particular, the PTSM). In test problems, agreement of the theory with computational experiments is shown, and a comparison of the PTSM and SSOR is done.     

A generalized Laguerre–Legendre spectral collocation method for solving initial-boundary value problems

In this research, a mixed spectral collocation method based on Kronecker product is proposed for solving initial-boundary value problems. New implementation is suggested to achieve more accurate approximation at longer times. Test problems are also studied to demonstrate how this method is implemented. Numerical experiments reveal that the new method is very effective and convenient.     

Asymptotic profiles of solutions to convection–diffusion equations

The large time behavior of zero-mass solutions to the Cauchy problem for the convection–diffusion equation $\text{u}{}_{\mathrm{t}}\text{?u}{}_{\mathrm{xx}}\text{+(|u|}{}^{\mathrm{q}}\text{)}{}_{\mathrm{x}}\text{=0,}\text{u(x,0)=u}{}_0\text{(x)}$ is studied when q>1 and the initial datum u₀ belongs to $\text{L}{}^1\text{(}\text{R}\text{,(1+|x|)}\text{d}\text{x)}$ and satisfies $\text{∫}{}_{\mathrm{<mtext>R</mtext>}}\text{u}{}_0\text{(x)}\text{d}\text{x=0}$. We provide conditions on the size and shape of the initial datum u₀ as well as on the exponent q>1 such that the large time asymptotics of solutions is given either by the derivative of the Gauss–Weierstrass kernel, or by a self-similar solution of the equation, or by hyperbolic N-waves. To cite this article: S. Benachour et al., C. R. Acad. Sci. Paris, Ser. I 338 (2004).     

A non-overlapping DDM combined with the characteristic method for optimal control problems governed by convection–diffusion equations

In this paper, we consider a non-overlapping domain decomposition method combined with the characteristic method for solving optimal control problems governed by linear convection–diffusion equations. The whole domain is divided into non-overlapping subdomains, and the global optimal control problem is decomposed into the local problems in these subdomains. The integral mean method is utilized for the diffusion term to present an explicit flux calculation on the inter-domain boundary in order to communicate the local problems on the interfaces between subdomains. The convection term is discretized along the characteristic direction. We establish the fully parallel and discrete schemes for solving these local problems. A priori error estimates in \(L^2\)-norm are derived for the state, co-state and control variables. Finally, we present numerical experiments to show the validity of the schemes and verify the derived theoretical results.     

A composite collocation method for the nonlinear mixed Volterra–Fredholm–Hammerstein integral equations

This paper presents a computational technique for the solution of the nonlinear mixed Volterra–Fredholm–Hammerstein integral equations. The method is based on the composite collocation method. The properties of hybrid of block-pulse functions and Lagrange polynomials are discussed and utilized to define the composite interpolation operator. The estimates for the errors are given. The composite interpolation operator together with the Gaussian integration formula are then used to transform the nonlinear mixed Volterra–Fredholm–Hammerstein integral equations into a system of nonlinear equations. The efficiency and accuracy of the proposed method is illustrated by four numerical examples.     

On numerical methods and error estimates for degenerate fractional convection–diffusion equations

First we introduce and analyze a convergent numerical method for a large class of nonlinear nonlocal possibly degenerate convection diffusion equations. Secondly we develop a new Kuznetsov type theory and obtain general and possibly optimal error estimates for our numerical methods—even when the principal derivatives have any fractional order between 1 and 2! The class of equations we consider includes equations with nonlinear and possibly degenerate fractional or general Levy diffusion. Special cases are conservation laws, fractional conservation laws, certain fractional porous medium equations, and new strongly degenerate equations.     

IIM-based ADI finite difference scheme for nonlinear convection–diffusion equations with interfaces

Based on Li’s immersed interface method (IIM), an ADI-type finite difference scheme is proposed for solving two-dimensional nonlinear convection–diffusion interface problems on a fixed cartesian grid, which is unconditionally stable and converges with two-order accuracy in both time and space in maximum norm. Correction terms are added to the right-hand side of standard ADI scheme at irregular points. The nonlinear convection terms are treated by Adams–Bashforth method, without affecting the stability of difference schemes. A new method for computing the correction terms is developed, in which the Adams–Bashforth method is employed. Thus we can get an explicit approximation for the computation of corrections, when the jump condition is solution-dependent. Three numerical experiments are displayed and analyzed. The numerical results show good agreement with the exact solutions and confirm the convergence order.