We consider the numerical approximation of a first order stationary hyperbolic equation by the method of characteristics with pseudo time step k using discontinuous finite elements on a mesh ${𝒯}_h$. For this method, we exhibit a “natural” norm || ||h,k for which we show that the discrete variational problem $P_h^k$ is well posed and we obtain an error estimate. We show that when k goes to zero problem $\left(P_h^k\right)$ (resp. the || ||h,k norm) has as a limit problem (Ph) (resp. the || ||h norm) associated to the...

Finite element methods with piecewise polynomial spaces in space for solving the nonstationary heat equation, as a model for parabolic equations are considered. The discretization in time is performed using the Crank-Nicolson method. A new a priori estimate is proved. Thanks to this new a priori estimate, a new error estimate in the discrete norm of ${𝒲}^{1,\infty }\left({ℒ}^2\right)$ is proved. An ${ℒ}^{\infty }\left({\mathscr{H}}^1\right)$-error estimate is also shown. These error estimates are useful since they allow us to get second order time accurate approximations...

This work is concerned with the introduction of a new numerical scheme based on the discontinuous Galerkin (DG) method. We propose to follow the methodology of higher order finite volume schemes and introduce a reconstruction operator into the DG scheme. This operator constructs higher order piecewise polynomial reconstructions from the lower order DG scheme. Such a procedure was proposed already in [2] based on heuristic arguments, however we provide a rigorous derivation, which justifies the increased...

The motion of an incompressible fluid confined to a shallow basin with a slightly varying bottom topography is considered. Coriolis force, surface wind and pressure stresses, together with bottom and lateral friction stresses are taken into account. We introduce appropriate scalings into a three-dimensional anisotropic eddy viscosity model; after averaging on the vertical direction and considering some asymptotic assumptions, we obtain a two-dimensional model, which approximates the three-dimensional...

The motion of an incompressible fluid confined to a shallow basin with a slightly varying bottom topography is considered. Coriolis force, surface wind and pressure stresses, together with bottom and lateral friction stresses are taken into account. We introduce appropriate scalings into a three-dimensional anisotropic eddy viscosity model; after averaging on the vertical direction and considering some asymptotic assumptions, we obtain a two-dimensional model, which approximates the three-dimensional...

This paper is devoted to Hermite interpolation with Chebyshev-Lobatto and Chebyshev-Radau nodal points. The aim of this piece of work is to establish the rate of convergence for some types of smooth functions. Although the rate of convergence is similar to that of Lagrange interpolation, taking into account the asymptotic constants that we obtain, the use of this method is justified and it is very suitable when we dispose of the appropriate information.

We discuss the formulation of a simulator in three spatial dimensions for a multicomponent, two phase (air, water) system of groundwater flow and transport with biodegradation kinetics and wells with multiple screens. The simulator has been developed for parallel, distributed memory, message passing machines. The numerical procedures employed are a fully implicit expanded mixed finite element method for flow and either a characteristics-mixed method or a Godunov method for transport and reactions...

This paper presents a postprocessing technique for estimating the local regularity of numerical solutions in high-resolution finite element schemes. A derivative of degree p ≥ 0 is considered to be smooth if a discontinuous linear reconstruction does not create new maxima or minima. The intended use of this criterion is the identification of smooth cells in the context of p-adaptation or selective flux limiting. As a model problem, we consider a 2D convection equation discretized with bilinear finite...

We analyze Euler-Galerkin approximations (conforming finite elements in space and implicit Euler in time) to coupled PDE systems in which one dependent variable, say $u$, is governed by an elliptic equation and the other, say $p$, by a parabolic-like equation. The underlying application is the poroelasticity system within the quasi-static assumption. Different polynomial orders are used for the $u$- and $p$-components to obtain optimally convergent a priori bounds for all the terms in the error energy norm....

We analyze Euler-Galerkin approximations (conforming finite elements in space and implicit Euler in time) to coupled PDE systems in which one dependent variable, say u, is governed by an elliptic equation and the other, say p, by a parabolic-like equation. The underlying application is the poroelasticity system within the quasi-static assumption. Different polynomial orders are used for the u- and p-components to obtain optimally convergent a priori bounds for all the terms in the error energy...

In this paper, we extend the reduced-basis methods and associated a posteriori error estimators developed earlier for elliptic partial differential equations to parabolic problems with affine parameter dependence. The essential new ingredient is the presence of time in the formulation and solution of the problem – we shall “simply” treat time as an additional, albeit special, parameter. First, we introduce the reduced-basis recipe – Galerkin projection onto a space $W_N$ spanned by solutions of the...

In this paper, we extend the reduced-basis methods and associated a posteriori error estimators developed earlier for elliptic partial differential equations to parabolic problems with affine parameter dependence. The essential new ingredient is the presence of time in the formulation and solution of the problem – we shall “simply” treat time as an additional, albeit special, parameter. First, we introduce the reduced-basis recipe – Galerkin projection onto a space WN spanned by solutions...

Phase-field models, the simplest of which is Allen–Cahn’s problem, are characterized by a small parameter $\epsilon$ that dictates the interface thickness. These models naturally call for mesh adaptation techniques, which rely on a posteriori error control. However, their error analysis usually deals with the underlying non-monotone nonlinearity via a Gronwall argument which leads to an exponential dependence on ${\epsilon }^{-2}$. Using an energy argument combined with a topological continuation argument and a spectral...

Phase-field models, the simplest of which is Allen–Cahn's problem, are characterized by a small parameter ε that dictates the interface thickness. These models naturally call for mesh adaptation techniques, which rely on a posteriori error control. However, their error analysis usually deals with the underlying non-monotone nonlinearity via a Gronwall argument which leads to an exponential dependence on ε-2. Using an energy argument combined with a topological continuation argument and...

The paper presents an a posteriori error estimator for a (piecewise linear) nonconforming finite element approximation of the heat equation in ${ℝ}^d$, $d=2$ or 3, using backward Euler’s scheme. For this discretization, we derive a residual indicator, which use a spatial residual indicator based on the jumps of normal and tangential derivatives of the nonconforming approximation and a time residual indicator based on the jump of broken gradients at each time step. Lower and upper bounds form the main results...

The paper presents an a posteriori error estimator for a (piecewise linear) nonconforming finite element approximation of the heat equation in ${ℝ}^d$, d=2 or 3, using backward Euler's scheme. For this discretization, we derive a residual indicator, which use a spatial residual indicator based on the jumps of normal and tangential derivatives of the nonconforming approximation and a time residual indicator based on the jump of broken gradients at each time step. Lower and upper bounds form the main...

Systems of parabolic differential equations are studied in the paper. Two a posteriori error estimates for the approximate solution obtained by the finite element method of lines are presented. A statement on the rate of convergence of the approximation of error by estimator to the error is proved.