Superconvergence of the direct discontinuous Galerkin method for convection‐diffusion equations

This paper is concerned with superconvergence properties of the direct discontinuous Galerkin (DDG) method for one‐dimensional linear convection‐diffusion equations. We prove, under some suitable choice of numerical fluxes and initial discretization, a 2k‐th and ‐th order superconvergence rate of the DDG approximation at nodes and Lobatto points, respectively, and a ‐th order of the derivative approximation at Gauss points, where k is the polynomial degree. Moreover, we also prove that the DDG solution is superconvergent with an order k + 2 to a particular projection of the exact solution. Numerical experiments are presented to validate the theoretical results. © 2016 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 290–317, 2017     

The Vlasov–Maxwell–Fokker–Planck system in two space dimensions

We consider the Cauchy problem for the Vlasov–Maxwell–Fokker–Planck system in the plane. It is shown that for smooth initial data, as long as the electromagnetic fields remain bounded, then their derivatives do also. Glassey and Strauss have shown this to hold for the relativistic Vlasov–Maxwell system in three dimensions, but the method here is totally different. In the work of Glassey and Strauss, the relativistic nature of the particle transport played an essential role. In this work, the transport is nonrelativistic, and smoothing from the Fokker–Planck operator is exploited. Copyright © 2015 John Wiley & Sons, Ltd.     

Local discontinuous Galerkin methods with explicit-implicit-null time discretizations for solving nonlinear diffusion problems

In this paper,we discuss the local discontinuous Galerkin methods coupled with two specific explicitimplicit-null time discretizations for solving one-dimensional nonlinear diffusion problems Ut=(a(U)Ux)x.The basic idea is to add and subtract two equal terms a0 Uxx the right-hand side of the partial differential equation,then to treat the term a0 Uxx implicitly and the other terms(a(U)Ux)x-a0 Uxx explicitly.We give stability analysis for the method on a simplified model by the aid of energy analysis,which gives a guidance for the choice of a0,i.e.,a0≥max{a(u)}/2 to ensure the unconditional stability of the first order and second order schemes.The optimal error estimate is also derived for the simplified model,and numerical experiments are given to demonstrate the stability,accuracy and performance of the schemes for nonlinear diffusion equations.     

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     

Numerical solution of Fokker‐Planck equation using the cubic B‐spline scaling functions

In this article a numerical technique is presented for the solution of Fokker‐Planck equation. This method uses the cubic B‐spline scaling functions. The method consists of expanding the required approximate solution as the elements of cubic B‐spline scaling function. Using the operational matrix of derivative, the problem will be reduced to a set of algebraic equations. Some numerical examples are included to demonstrate the validity and applicability of the technique. The method is easy to implement and produces very accurate results. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2009     

An extension of the Gegenbauer pseudospectral method for the time fractional Fokker‐Planck equation

The time fractional Fokker‐Planck equation has been used in many physical transport problems which take place under the influence of an external force field. In this paper we examine pseudospectral method based on Gegenbauer polynomials and Chebyshev spectral differentiation matrix to solve numerically a class of initial‐boundary value problems of the time fractional Fokker‐Planck equation on a finite domain. The presented method reduces the main problem to a generalized Sylvester matrix equation, which can be solved by the global generalized minimal residual method. Some numerical experiments are considered to demonstrate the accuracy and the efficiency of the proposed computational procedure.     

Local discontinuous Galerkin method for solving an N‐carrier system

In this article, we present local discontinuous Galerkin (LDG) method for solving a model of energy exchanges in an N ‐carrier system with Neumann boundary conditions. This model extends the concept of the well‐known parabolic two‐step model for microheat transfer to the energy exchanges in a generalized N ‐carrier system with heat sources. The energy norm stability and error estimate of the LDG method is proved for solving N ‐carrier system. Some numerical examples are given. The numerical results when compared with the exact solution and other numerical results indicate that the present method is seen to be a very good alternative to some existing techniques for realistic problems. © 2011 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2011     

Local discontinuous Galerkin approximation of non-Fickian diffusion model in viscoelastic polymers

In this paper, we investigate a fully discrete local discontinuous Galerkin approximation of a non-linear non-Fickian diffusion model in viscoelastic polymers. For the spatial discretization, we adopt local discontinuous Galerkin finite element method and for the time discretization we use backward Euler method. We derive the stability estimate and a priori error estimate for the discrete scheme. Numerical examples are given to verify the theoretical findings.     

Convergence order estimates of the local discontinuous Galerkin method for instationary Darcy flow

We present a new a priori stability and convergence analysis for the local discontinuous Galerkin method applied to the instationary Darcy problem formulated on a d‐dimensional polytope with nonhomogeneous Neumann and Dirichlet boundary conditions. In addition to including a spatially and temporally varying permeability tensor into all estimates, the utilized analysis technique produces a convergence order result for the primary and the flux variables. The only stabilization in the proposed scheme is represented by a penalty term in the primary unknown, and our analysis provides some insights into the role played by this particular choice of stabilization. © 2016 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 1374–1394, 2017     

Interior penalty discontinuous Galerkin methods with implicit time‐integration techniques for nonlinear parabolic equations

We prove existence and numerical stability of numerical solutions of three fully discrete interior penalty discontinuous Galerkin methods for solving nonlinear parabolic equations. Under some appropriate regularity conditions, we give the l²(H¹) and l^∞(L²) error estimates of the fully discrete symmetric interior penalty discontinuous Galerkin–scheme with the implicit θ ‐schemes in time, which include backward Euler and Crank–Nicolson finite difference approximations. Our estimates are optimal with respect to the mesh size h. The theoretical results are confirmed by some numerical experiments. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

Analysis of a discontinuous galerkin method for heterogeneous diffusion problems with low‐regularity solutions

We study the convergence of the Symmetric Weighted Interior Penalty discontinuous Galerkin method for heterogeneous diffusion problems with low‐regularity solutions only belonging to W^{2, p} with p ∈ (1,2]. In 2d, we infer an optimal algebraic convergence rate. In any space dimension d, we achieve the same result for p > 2 d/( d + 2). We also prove convergence without algebraic rates for exact solutions only belonging to the energy space. © 2011 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2012     

Exponentially fitted discontinuous Galerkin schemes for singularly perturbed problems

New discontinuous Galerkin schemes in mixed form are introduced for symmetric elliptic problems of second order. They exhibit reduced connectivity with respect to the standard ones. The modifications in the choice of the approximation spaces and in the stabilization term do not spoil the error estimates. These methods are then used for designing new exponentially fitted schemes for advection dominated equations. The presented numerical tests show the good performances of the proposed schemes. © 2011 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2011     

一种推广的对流扩散方程的局部化间断Galerkin方法

研究求解一种产生于径向渗流问题的推广的对流扩散方程的局部化间断Galerkin方法，对一般非线性情形证明了方法的L^2稳定性；对线性情形证明了，当方法取有限元空间为κ次多项式空间时，数值解逼近的L^∞（0，T；L^2)模的误差阶为κ。     

Composite generalized Laguerre spectral method for nonlinear Fokker–Planck equations on the whole line

In this paper, we propose a composite Laguerre spectral method for the nonlinear Fokker–Planck equations modelling the relaxation of fermion and boson gases. A composite Laguerre spectral scheme is constructed. Its convergence is proved. Numerical results show the efficiency of this approach and coincide well with theoretical analysis. Some results on the Laguerre approximation and techniques used in this paper are also applicable to other nonlinear problems on the whole line. Copyright © 2016 John Wiley & Sons, Ltd.     

Discontinuous Galerkin and multiscale variational schemes for a coupled damped nonlinear system of Schrödinger equations

In this article, we study a streamline diffusion‐based discontinuous Galerkin approximation for the numerical solution of a coupled nonlinear system of Schrödinger equations and extend the resulting method to a multiscale variational scheme. We prove stability estimates and derive optimal convergence rates due to the maximal available regularity of the exact solution. In the weak formulation, to make the underlying bilinear form coercive, it was necessary to supply the equation system with an artificial viscosity term with a small coefficient of order proportional to a power of mesh size. We justify the theory by implementing an example of an application of the time‐dependent Schrödinger equation in the coupled ultrafast laser. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

Global weak solutions to the Cauchy problem of compressible Navier–Stokes–Vlasov–Fokker–Planck equations

A fluid–particles system of the compressible Navier‐Stokes equations and Vlasov‐Fokker‐Planck equation (including the case of Vlasov equation) in three‐dimensional space is considered in this paper. The coupling arises from a drag force exerted by the fluid onto the particles. We study a Cauchy problem with large data, and establish the existence of global weak solutions through an approximation scheme, energy estimates, and weak convergence. Copyright © 2016 John Wiley & Sons, Ltd.     

Convergence of Galerkin approximations for the Kuramoto‐Tsuzuki equation

Standard Galerkin approximations, using smooth splines to solutions of the Kuramoto‐Tsuzuki equation are analyzed. The existence, uniqueness, and convergence of the fully discrete Crank‐Nicolson scheme are discussed. Furthermore, a second‐order convergent linearized Galerkin approximation are derived. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 21, 2005     

Convergence analysis of a coupled method for time‐dependent convection‐diffusion equations

The coupling of two locally mass conservative methods is formulated and analyzed for the time‐dependent convection‐diffusion problem. Finite volume method is used in some subdomains and interior penalty discontinuous Galerkin method is used in other subdomains. Numerical examples show the advantages of the proposed hybrid method, namely an accurate approximation obtained at a reduced computational cost. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 133–157, 2014     

Analysis of the parareal approach based on discontinuous Galerkin method for time-dependent Stokes equations

This paper analyzes a parareal approach based on discontinuous Galerkin (DG) method for the time-dependent Stokes equations. A class of primal discontinuous Galerkin methods, namely variations of interior penalty methods, are adopted for the spatial discretization in the parareal algorithm (we call it parareal DG algorithm). We study three discontinuous Galerkin methods for the time-dependent Stokes equations, and the optimal continuous in time error estimates for the velocities and pressure are derived. Based on these error estimates, the proposed parareal DG algorithm is proved to be unconditionally stable and bounded by the error of discontinuous Galerkin discretization after a finite number of iterations. Finally, some numerical experiments are conducted which confirm our theoretical results, meanwhile, the efficiency of the parareal DG algorithm can be seen through a parallel experiment.     

A backward euler orthogonal spline collocation method for the time‐fractional Fokker–Planck equation

We formulate and analyze a novel numerical method for solving a time‐fractional Fokker–Planck equation which models an anomalous subdiffusion process. In this method, orthogonal spline collocation is used for the spatial discretization and the time‐stepping is done using a backward Euler method based on the L1 approximation to the Caputo derivative. The stability and convergence of the method are considered, and the theoretical results are supported by numerical examples, which also exhibit superconvergence. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 1534–1550, 2015