L-asymptotic behavior for a finite element approximation in parabolic quasi-variational inequalities related to impulse control problem

In this paper, the parabolic quasi-variational inequalities are transformed into a noncoercive elliptic quasi-variational inequalities. A new iterative discrete algorithm is proposed to show the existence and uniqueness, and a simple proof to asymptotic behavior in uniform norm is also given using the theta time scheme combined with a finite element spatial approximation. The proposed approach stands on a discrete L^∞-stability property with respect to the right-hand side and obstacle defined as an impulse control problem.     

New a priori analysis of first‐order system least‐squares finite element methods for parabolic problems

We provide new insights into the a priori theory for a time‐stepping scheme based on least‐squares finite element methods for parabolic first‐order systems. The elliptic part of the problem is of general reaction‐convection‐diffusion type. The new ingredient in the analysis is an elliptic projection operator defined via a nonsymmetric bilinear form, although the main bilinear form corresponding to the least‐squares functional is symmetric. This new operator allows to prove optimal error estimates in the natural norm associated to the problem and, under additional regularity assumptions, in the L² norm. Numerical experiments are presented which confirm our theoretical findings.     

Conjugate gradient-boundary element solution using multiple reciprocity method for distributed elliptic optimal control problems

This paper deals with the numerical solution of optimal control problems, where the state equations are given by the fourth order elliptic partial differential equations. An iterative algorithm for this class of problems is developed. This new proposal is obtained by combining the Conjugate Gradient Method (CGM) with the Boundary Element Method (BEM) and Multiple Reciprocity Method (MRM). The local error estimates based on the stability of this scheme in the H² norm, L² norm and L^∞ norm are obtained. Finally, the numerical results on a test case show that this method is correct and feasible.     

A posteriori error estimates for optimal distributed control governed by the evolution equations

We describe a technique for a posteriori error estimates suitable to the optimal control problem governed by the evolution equations solved by the method of lines. It is applied to the control problem governed by the parabolic equation, convection-diffusion equation and hyperbolic equation. The error is measured with the aid of the L²-norm in the space-time cylinder combined with a special time weighted energy norm.     

A new fourth‐order numerical algorithm for a class of three‐dimensional nonlinear evolution equations

In this article, a new compact alternating direction implicit finite difference scheme is derived for solving a class of 3‐D nonlinear evolution equations. By the discrete energy method, it is shown that the new difference scheme has good stability and can attain second‐order accuracy in time and fourth‐order accuracy in space with respect to the discrete H¹ ‐norm. A Richardson extrapolation algorithm is applied to achieve fourth‐order accuracy in temporal dimension. Numerical experiments illustrate the accuracy and efficiency of the extrapolation algorithm. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

A high order finite volume method for one dimensional nonlocal reactive flows of parabolic type

In this paper, we consider a high order finite volume approximation of one‐dimensional nonlocal reactive flows of parabolic type. The method is obtained by discretizing in space by arbitrary order vertex‐centered finite volumes, followed by a modified Simpson quadrature scheme for the time stepping. Compared to the existed finite volume methods, this new finite volume scheme could achieve the desired accuracy with less data storage by employing higher‐order trial spaces. The finite volume approximations are proved to possess optimal order convergence rates in the H¹‐norm and L²‐norm, which are also confirmed by numerical tests. Copyright © 2015 John Wiley & Sons, Ltd.     

A new robust C-type nonconforming triangular element for singular perturbation problems

In this paper, a new robust C⁰ triangular element is proposed for the fourth order elliptic singular perturbation problem with double set parameter method and bubble function technique, and a general convergence theorem for C⁰ nonconforming elements is presented. The convergence of the new element is proved in the energy norm uniformly with respect to the perturbation parameter. Numerical experiments are also carried out to demonstrate the efficiency of the new element.     

A posteriori error estimates for finite element approximations to the wave equation with discontinuous coefficients

We derive residual‐based a posteriori error estimates of finite element method for linear wave equation with discontinuous coefficients in a two‐dimensional convex polygonal domain. A posteriori error estimates for both the space‐discrete case and for implicit fully discrete scheme are discussed in L^∞(L²) norm. The main ingredients used in deriving a posteriori estimates are new Clément type interpolation estimates in conjunction with appropriate adaption of the elliptic reconstruction technique of continuous and discrete solutions. We use only an energy argument to establish a posteriori error estimates with optimal order convergence in the L^∞(L²) norm.     

A method of asymptotic constructions with improved accuracy for a quasilinear singularly perturbed parabolic convection-diffusion equation

The Dirichlet problem on an interval for quasilinear singularly perturbed parabolic convection-diffusion equation is considered. The higher order derivative of the equation is multiplied by a parameter ε that takes any values from the half-open interval (0, 1]. For this type of linear problems, the order of the ε-uniform convergence (with respect to x and t) for the well-known schemes is not higher than unity (in the maximum norm). For the boundary value problem under consideration, grid approximations are constructed that converge ε-uniformly at the rate of O(N ^?2ln² N + N ^?2 ₀), where N + 1 and N ₀ + 1 are the numbers of the mesh points with respect to x and t, respectively. On the x axis, piecewise uniform meshes that condense in the boundary layer are used. If the parameter value is small compared to the effective step of the spatial grid, the domain decomposition method is used, which is motivated by “asymptotic constructions.” Monotone approximations of “auxiliary” subproblems describing the main terms of the asymptotic expansion of the solution outside a neighborhood of the boundary layer neighborhood are used. In the neighborhood of the boundary layer (of the width O(ε ln N)) the first derivative with respect to x is approximated by the central difference derivative. These subproblems are successively solved in the subdomains on uniform grids. If the parameter values are not sufficiently small (compared to the effective step of the mesh with respect to x), the classical implicit difference schemes approximating the first derivative with respect to x by the central difference derivative are applied. To improve the accuracy in t, the defect correction technique is used. Notice that the calculation of the solution of the constructed difference scheme (the scheme based on the method of asymptotic constructions) can be considerably simplified for sufficiently small values of the parameter ε.     

Single and multi‐interval Legendre spectral methods in time for parabolic equations

In this article, we take the parabolic equation with Dirichlet boundary conditions as a model to present the Legendre spectral methods both in spatial and in time. Error analysis for the single/multi‐interval schemes in time is given. For the single interval spectral method in time, we obtain a better error estimate in L²‐norm. For the multi‐interval spectral method in time, we obtain the L²‐optimal error estimate in spatial. By choosing approximate trial and test functions, the methods result in algebraic systems with sparse forms. A parallel algorithm is constructed for the multi‐interval scheme in time. Numerical results show the efficiency of the methods. The methods are also applied to parabolic equations with Neumann boundary conditions, Robin boundary conditions and some nonlinear PDEs. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2006     

Crank-Nicolson split least-squares procedure for nonlocal reactive flow in porous media

In this paper, we introduce a Crank-Nicolson split least-squares Galerkin finite element procedure for parabolic integro-differential equations, arising in the modeling of nonlocal reactive flows in porous media. By selecting the least-squares functional properly, the procedure can be split into two independent sub-procedures, one of which is for the primitive unknown and the other is for the flux. By carefully choosing projections, we get optimal order H ¹(Ω) and L ²(Ω) norm error estimates for u and sub-optimal (L ²(Ω)) ^d norm error estimate for σ with second-order accuracy in time increment. The numerical examples are given to testify the efficiency of the introduced scheme.     

Finite element methods for semilinear elliptic problems with smooth interfaces

The purpose of this paper is to study the finite element method for second order semilinear elliptic interface problems in two dimensional convex polygonal domains. Due to low global regularity of the solution, it seems difficult to achieve optimal order of convergence with straight interface triangles [Numer. Math., 79 (1998), pp. 175–202]. For a finite element discretization based on a mesh which involve the approximation of the interface, optimal order error estimates in L ² and H ¹-norms are proved for linear elliptic interface problem under practical regularity assumptions of the true solution. Then an extension to the semilinear problem is also considered and optimal error estimate in H ¹ norm is achieved.     

On a parabolic logarithmic Sobolev inequality

In order to extend the blow-up criterion of solutions to the Euler equations, Kozono and Taniuchi [H. Kozono, Y. Taniuchi, Limiting case of the Sobolev inequality in BMO, with application to the Euler equations, Comm. Math. Phys. 214 (2000) 191-200] have proved a logarithmic Sobolev inequality by means of isotropic (elliptic) BMO norm. In this paper, we show a parabolic version of the Kozono-Taniuchi inequality by means of anisotropic (parabolic) BMO norm. More precisely we give an upper bound for the L^∞ norm of a function in terms of its parabolic BMO norm, up to a logarithmic correction involving its norm in some Sobolev space. As an application, we also explain how to apply this inequality in order to establish a long-time existence result for a class of nonlinear parabolic problems.     

A Crank‐Nicolson Petrov‐Galerkin method with quadrature for semi‐linear parabolic problems

We propose and analyze an application of a fully discrete C² spline quadrature Petrov‐Galerkin method for spatial discretization of semi‐linear parabolic initial‐boundary value problems on rectangular domains. We prove second order in time and optimal order H¹ norm convergence in space for the extrapolated Crank‐Nicolson quadrature Petrov‐Galerkin scheme. We demonstrate numerically both L² and H¹ norm optimal order convergence of the scheme even if the nonlinear source term is not smooth. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2005.     

Biquadratic finite volume element methods based on optimal stress points for parabolic problems

In this paper, the semi-discrete and full discrete biquadratic finite volume element schemes based on optimal stress points for a class of parabolic problems are presented. Optimal order error estimates in H¹ and L² norms are derived. In addition, the superconvergences of numerical gradients at optimal stress points are also discussed. A numerical experiment confirms some results of theoretical analysis.     

Finite element methods with numerical quadrature for elliptic problems with smooth interfaces

The purpose of this paper is to study the effect of the numerical quadrature on the finite element approximation to the exact solution of elliptic equations with discontinuous coefficients. Due to low global regularity of the solution, it seems difficult to achieve optimal order of convergence with classical finite element methods [Z. Chen, J. Zou, Finite element methods and their convergence for elliptic and parabolic interface problems, Numer. Math. 79 (1998) 175-202]. We derive error estimates in finite element method with quadrature for elliptic interface problems in a two-dimensional convex polygonal domain. Optimal order error estimates in L² and H¹ norms are shown to hold even if the regularity of the solution is low on the whole domain. Finally, numerical experiment for two dimensional test problem is presented in support of our theoretical findings.     

Weak Galerkin finite element methods for a fourth order parabolic equation

This paper is devoted to a newly developed weak Galerkin finite element method with the stabilization term for a linear fourth order parabolic equation, where weakly defined Laplacian operator over discontinuous functions is introduced. Priori estimates are developed and analyzed in L² and an H² type norm for both semi‐discrete and fully discrete schemes. And finally, numerical examples are provided to confirm the theoretical results.     

Approximate controllability and homogenization of a semilinear elliptic problem

The L²- and H¹-approximate controllability and homogenization of a semilinear elliptic boundary-value problem is studied in this paper. The principal term of the state equation has rapidly oscillating coefficients and the control region is locally distributed. The observation region is a subset of codimension 1 in the case of L²-approximate controllability or is locally distributed in the case of H¹-approximate controllability. By using the classical Fenchel-Rockafellar's duality theory, the existence of an approximate control of minimal norm is established by means of a fixed point argument. We consider its asymptotic behavior as the rapidly oscillating coefficients H-converge. We prove its convergence to an approximate control of minimal norm for the homogenized problem.     

A low order anisotropic nonconforming characteristic finite element method for a convection-dominated transport problem

A low order anisotropic nonconforming rectangular finite element method for the convection-diffusion problem with a modified characteristic finite element scheme is studied in this paper. The O(h²) order error estimate in L²-norm with respect to the space, one order higher than the expanded characteristic-mixed finite element scheme with order O(h), and the same as the conforming case for a modified characteristic finite element scheme under regular meshes, is obtained by use of some distinct properties of the interpolation operator and the mean value technique, instead of the so-called elliptic projection, which is an indispensable tool in the convergence analysis of the previous literature. Lastly, some numerical results of the element are provided to verify our theoretical analysis.     

A second order scheme for a time-dependent,singularly perturbed convection–diffusion equation

We consider a numerical scheme for a one-dimensional, time-dependent, singularly perturbed convection–diffusion problem. The problem is discretized in space by a standard finite element method on a Bakhvalov–Shishkin type mesh. The space error is measured in an L² norm. For the time integration, the implicit midpoint rule is used. The fully discrete scheme is shown to be convergent of order 2 in space and time, uniformly in the singular perturbation parameter.