Analysis of a BDF–DGFE scheme for nonlinear convection–diffusion problems

We deal with the numerical solution of a scalar nonstationary nonlinear convection–diffusion equation. We employ a combination of the discontinuous Galerkin finite element method for the space semi-discretization and the k-step backward difference formula for the time discretization. The diffusive and stabilization terms are treated implicitly whereas the nonlinear convective term is treated by a higher order explicit extrapolation method, which leads to the necessity to solve only a linear algebraic problem at each time step. We analyse this scheme and derive a priori asymptotic error estimates in the discrete L ^∞(L ²)-norm and the L ²(H ¹)-seminorm with respect to the mesh size h and time step τ for k = 2,3. Numerical examples verifying the theoretical results are presented. This work is a part of the research project MSM 0021620839 financed by the Ministry of Education of the Czech Republic and was partly supported by the Grant No. 316/2006/B-MAT/MFF of the Grant Agency of the Charles University Prague. The research of M. Vlasák was supported by the project LC06052 of the Ministry of Education of the Czech Republic (Jindřich Nečas Center for Mathematical Modelling).     

A priori error estimation for the dual mixed finite element method of the elastodynamic problem in a polygonal domain, II

In this paper we analyze a new dual mixed formulation of the elastodynamic system in polygonal domains by using an implicit scheme for the time discretization. After the analysis of stability of the fully discrete scheme, L^∞ in time, L² in space a priori error estimates for the approximation of the displacement, the strain, the pressure and the rotational are derived. Numerical tests are presented which confirm our theoretical results.     

A priori error estimation for the dual mixed finite element method of the elastodynamic problem in a polygonal domain, I

In this paper we analyze a new dual mixed formulation of the elastodynamic system in polygonal domains. In this formulation the symmetry of the strain tensor is relaxed by the rotation of the displacement. For the time discretization of this new dual mixed formulation, we use an explicit scheme. After the analysis of stability of the fully discrete scheme, L^∞ in time, L² in space a priori error estimates are derived for the approximation of the displacement, the strain, the pressure and the rotation. Numerical experiments confirm our theoretical predictions.     

Optimal error estimate and superconvergence of the DG method for first-order hyperbolic problems

We consider the original discontinuous Galerkin method for the first-order hyperbolic problems in d-dimensional space. We show that, when the method uses polynomials of degree k, the L₂-error estimate is of order k+1 provided the triangulation is made of rectangular elements satisfying certain conditions. Further, we show the O(h^2k+1)-order superconvergence for the error on average on some suitably chosen subdomains (including the whole domain) and their outflow faces. Moreover, we also establish a derivative recovery formula for the approximation of the convection directional derivative which is superconvergent with order k+1.     

A new high accuracy locally one-dimensional scheme for the wave equation

In this paper, a new locally one-dimensional (LOD) scheme with error of O(Δt⁴+h⁴) for the two-dimensional wave equation is presented. The new scheme is four layer in time and three layer in space. One main advantage of the new method is that only tridiagonal systems of linear algebraic equations have to be solved at each time step. The stability and dispersion analysis of the new scheme are given. The computations of the initial and boundary conditions for the two intermediate time layers are explicitly constructed, which makes the scheme suitable for performing practical simulation in wave propagation modeling. Furthermore, a comparison of our new scheme and the traditional finite difference scheme is given, which shows the superiority of our new method.     

Domain decomposition and splitting methods for Mortar mixed finite element approximations to parabolic equations

Summary We introduce in this article a new domain decomposition algorithm for parabolic problems that combines Mortar Mixed Finite Element methods for the space discretization with operator splitting schemes for the time discretization. The main advantage of this method is to be fully parallel. The algorithm is proven to be unconditionally stable and a convergence result in (Δt/h ^1/2) is presented.     

Spectral element methods for parabolic problems

A spectral element method for solving parabolic initial boundary value problems on smooth domains using parallel computers is presented in this paper. The space domain is divided into a number of shape regular quadrilaterals of size h and the time step k   is proportional to h²$h^2$. At each time step we minimize a functional which is the sum of the squares of the residuals in the partial differential equation, initial condition and boundary condition in different Sobolev norms and a term which measures the jump in the function and its derivatives across inter-element boundaries in certain Sobolev norms. The Sobolev spaces used are of different orders in space and time. We can define a preconditioner for the minimization problem which allows the problem to decouple. Error estimates are obtained for both the h and p versions of this method.     

Multistep characteristic method for incompressible flow in porous media

In this paper we present a multistep difference scheme for the problem of miscible displacement of incompressible fluid flow in porous media. The discretization involves a three-level time scheme based on the characteristic method and a five-point finite difference scheme for space discretization. We prove that the convergence is of order O(h²+(Δt)²), which is in contrast to the convergence of order O(h+Δt) proved for a singlestep characteristic with the same space discretization. Numerical experiments demonstrate the stability and second-order convergence of the scheme.     

A note onC o Galerkin methods for two-point boundary problems

Summary As is known [4]. theC ^o Galerkin solution of a two-point boundary problem using piecewise polynomial functions, hasO(h ^2k ) convergence at the knots, wherek is the degree of the finite element space. Also, it can be proved [5] that at specific interior points, the Gauss-Legendre points the gradient hasO(h ^k+1) convergence, instead ofO(h ^k ). In this note, it is proved that on any segment there arek–1 interior points where the Galerkin solution is ofO(h ^k+2), one order better than the global order of convergence. These points are the Lobatto points.     

Stability of a family of multi-time-level difference schemes for the advection equation

Summary For the linear advection equation we consider explicit multi-time-level schemes of highest order which are one step in space direction only. If a stencil involvesk time steps we show that it is stable in theL ₂-sense for Courant numbers in the interval (0, 1/k). Since the order is 2k–1 one can use these schemes for high order discretization of the boundary conditions in hyperbolic initial value problems.Part of this work has been performed in the project Mehrschritt-Differenzenschemata of the Schwerpunktprogramm Finite Approximationen in der Strömungsmechanik which has been supported by the DFG     

Error estimate for finite volume scheme

We study the convergence of a finite volume scheme for the linear advection equation with a Lipschitz divergence-free speed in R ^d . We prove a h ^1/2-error estimate in the L ^∞(0,t;L ¹)-norm for BV data. This result was expected from numerical experiments and is optimal.     

Non-polynomials spline methods for solving a fourth-order parabolic equation

In this paper a fourth-order variable coefficient parabolic partial differential equation, that governs the behaviour of a vibrating beam, is solved by using a three level method based on non-polynomial quintic spline in space and finite difference discretization in time. We also obtain two new high accuracy schemes of O (k⁴, h⁶) and O (k⁴, h⁸) and two new schemes which are analogues of Jain's formula for the non-homogeneous case. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Correction of finite element eigenvalues for problems with natural or periodic boundary conditions

Recent results of Andrew and Paine for a regular Sturm-Liouville problem with essential boundary conditions are extended to problems with natural or periodic boundary conditions. These results show that a simple asymptotic correction technique of Paine, de Hoog and Anderssen reduces the error in the estimate of thekth eigenvalue obtained by the finite element method, with linear hat functions and mesh lengthh, fromO(k ⁴ h ²) toO(kh ²). Numerical results show the correction to be useful even for low values ofk.     

Error estimate of a first-order time discretization scheme for the geodynamo equations

In this paper, we analyze a first-order time discretization scheme for a nonlinear geodynamo model and carry out the convergence analysis of this numerical scheme. It is concluded that our numerical scheme converges with first-order accuracy in the sense of L²$L^2$-norm with respect to the velocity field u$\mathit{u}$ and the magnetic field B$\mathit{B}$ and with half-order accuracy in time for the total kinematic pressure P.     

Convergence analysis of a discontinuous Galerkin method for a sub-diffusion equation

We employ a piecewise-constant, discontinuous Galerkin method for the time discretization of a sub-diffusion equation. Denoting the maximum time step by k, we prove an a priori error bound of order k under realistic assumptions on the regularity of the solution. We also show that a spatial discretization using continuous, piecewise-linear finite elements leads to an additional error term of order h ² max (1,logk ^− 1). Some simple numerical examples illustrate this convergence behaviour in practice. We thank the University of New South Wales for financial support provided by a Faculty Research Grant.     

Correction of finite element estimates for Sturm-Liouville eigenvalues

Summary It is shown that a simple asymptotic correction technique of Paine, de Hoog and Anderssen reduces the error in the estimate of thekth eigenvalue of a regular Sturm-Liouville problem obtained by the finite element method, with linear hat functions and mesh lengthh, fromO(k ⁴ h ²) toO(k h ²). The result still holds when the matrix elements are evaluated by Simpson's rule, but if the trapezoidal rule is used the error isO(k ² h ²). Numerical results demonstrate the usefulness of the correction even for low values ofk.     

Error estimates of a linear decoupled Euler–FEM scheme for a mass diffusion model

We present error estimates of a linear fully discrete scheme for a three-dimensional mass diffusion model for incompressible fluids (also called Kazhikhov–Smagulov model). All unknowns of the model (velocity, pressure and density) are approximated in space by C ⁰-finite elements and in time an Euler type scheme is used decoupling the density from the velocity–pressure pair. If we assume that the velocity and pressure finite-element spaces satisfy the inf–sup condition and the density finite-element space contains the products of any two discrete velocities, we first obtain point-wise stability estimates for the density, under the constraint lim_(h,k)→0 h/k = 0 (h and k being the space and time discrete parameters, respectively), and error estimates for the velocity and density in energy type norms, at the same time. Afterwards, error estimates for the density in stronger norms are deduced. All these error estimates will be optimal (of order O(h+k){\mathcal{O}(h+k)}) for regular enough solutions without imposing nonlocal compatibility conditions at the initial time. Finally, we also study two convergent iterative methods for the two problems to solve at each time step, which hold constant matrices (independent of iterations).     

Boundary controllability of a linear semi-discrete 1-D wave equation derived from a mixed finite element method

In this article one discusses the controllability of a semi-discrete system obtained by discretizing in space the linear 1-D wave equation with a boundary control at one extremity. It is known that the semi-discrete models obtained with finite difference or the classical finite element method are not uniformly controllable as the discretization parameter h goes to zero (see [8]). Here we introduce a new semi-discrete model based on a mixed finite element method with two different basis functions for the position and velocity. We show that the controls obtained with these semi-discrete systems can be chosen uniformly bounded in L²(0,T) and in such a way that they converge to the HUM control of the continuous wave equation, i.e. the minimal L²-norm control. We illustrate the mathematical results with several numerical experiments. Supported by Grant BFM 2002-03345 of MCYT (Spain) and the TMR projects of the EU ``Homogenization and Multiple Scales" and ``New materials, adaptive systems and their nonlinearities: modelling, control and numerical simulations". Partially Supported by Grant BFM 2002-03345 of MCYT (Spain), Grant 17 of Egide-Brancusi Program and Grant 80/2005 of CNCSIS (Romania).     

Convergence of a finite volume scheme for nonlinear degenerate parabolic equations

Summary. One approximates the entropy weak solution u of a nonlinear parabolic degenerate equation by a piecewise constant function using a discretization in space and time and a finite volume scheme. The convergence of to u is shown as the size of the space and time steps tend to zero. In a first step, estimates on are used to prove the convergence, up to a subsequence, of to a measure valued entropy solution (called here an entropy process solution). A result of uniqueness of the entropy process solution is proved, yielding the strong convergence of to{\it u}. Some on a model equation are shown. Received September 27, 2000 / Published online October 17, 2001     

A Priori Error Estimates of Crank–Nicolson Finite Volume Element Method for a Hyperbolic Optimal Control Problem

In this article, a Crank–Nicolson linear finite volume element scheme is developed to solve a hyperbolic optimal control problem. We use the variational discretization technique for the approximation of the control variable. The optimal convergent order O(h² + k²) is proved for the numerical solution of the control, state and adjoint‐state in a discrete L²‐norm. To derive this result, we also get the error estimate (convergent order O(h² + k²)) of Crank–Nicolson finite volume element approximation for the second‐order hyperbolic initial boundary value problem. Numerical experiments are presented to verify the theoretical results.© 2016 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 1331–1356, 2016