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...

This paper is devoted to the study of a posteriori error estimates for the scalar nonlinear convection-diffusion-reaction equation $c_t+\nabla ·\left(𝐮f\left(c\right)\right)-\nabla ·\left(D\nabla c\right)+\lambda c=0$. The estimates for the error between the exact solution and an upwind finite volume approximation to the solution are derived in the $L^1$-norm, independent of the diffusion parameter $D$. The resulting a posteriori error estimate is used to define an grid adaptive solution algorithm for the finite volume scheme. Finally numerical experiments underline the applicability...

This paper is devoted to the study of a posteriori error estimates for the scalar nonlinear convection-diffusion-reaction equation $c_t+\nabla ·\left(𝐮f\left(c\right)\right)-\nabla ·\left(D\nabla c\right)+\lambda c=0$. The estimates for the error between the exact solution and an upwind finite volume approximation to the solution are derived in the L1-norm, independent of the diffusion parameter D. The resulting a posteriori error estimate is used to define an grid adaptive solution algorithm for the finite volume scheme. Finally numerical experiments underline the applicability...

We deal with the numerical solution of the nonstationary heat conduction equation with mixed Dirichlet/Neumann boundary conditions. The backward Euler method is employed for the time discretization and the interior penalty discontinuous Galerkin method for the space discretization. Assuming shape regularity, local quasi-uniformity, and transition conditions, we derive both a posteriori upper and lower error bounds. The analysis is based on the Helmholtz decomposition, the averaging interpolation...

We build a non-dissipative second order algorithm for the approximate resolution of the one-dimensional Euler system of compressible gas dynamics with two components. The considered model was proposed in [1]. The algorithm is based on [8] which deals with a non-dissipative first order resolution in Lagrange-remap formalism. In the present paper we describe, in the same framework, an algorithm that is second order accurate in time and space, and that...

We consider highly accurate schemes for nonlinear time fractional Schrödinger equations (NTFSEs). While an $L1$ strategy is employed for approximating the Caputo fractional derivative in the temporal direction, compact CCD finite difference approaches are incorporated in the space. A highly effective hybrid $L1$-CCD method is implemented successfully. The accuracy of this linearized scheme is order six in space, and order $2-\gamma$ in time, where $0<\gamma <1$ is the order of the Caputo fractional derivative involved. It...

We propose a simple numerical method for capturing the steady state solution of hyperbolic systems with geometrical source terms. We use the interface value, rather than the cell-averages, for the source terms that balance the nonlinear convection at the cell interface, allowing the numerical capturing of the steady state with a formal high order accuracy. This method applies to Godunov or Roe type upwind methods but requires no modification of the Riemann solver. Numerical experiments on scalar...

We propose a simple numerical method for capturing the steady state solution of hyperbolic systems with geometrical source terms. We use the interface value, rather than the cell-averages, for the source terms that balance the nonlinear convection at the cell interface, allowing the numerical capturing of the steady state with a formal high order accuracy. This method applies to Godunov or Roe type upwind methods but requires no modification of the Riemann solver. Numerical experiments on scalar...

This paper studies the exact controllability of a finite dimensional system obtained by discretizing in space and time the linear 1-D wave system with a boundary control at one extreme. It is known that usual schemes obtained with finite difference or finite element methods are not uniformly controllable with respect to the discretization parameters $h$ and $\Delta t$. We introduce an implicit finite difference scheme which differs from the usual centered one by additional terms of order $h^2$ and $\Delta t^2$. Using a discrete...

This paper studies the exact controllability of a finite dimensional system obtained by discretizing in space and time the linear 1-D wave system with a boundary control at one extreme. It is known that usual schemes obtained with finite difference or finite element methods are not uniformly controllable with respect to the discretization parameters h and Δt. We introduce an implicit finite difference scheme which differs from the usual centered one by additional terms of order h2 and Δt2. Using...

We present an alternative framework for designing efficient numerical schemes for non-conservative hyperbolic systems. This approach is based on the design of entropy conservative discretizations and suitable numerical diffusion operators that mimic the effect of underlying viscous mechanisms. This approach is illustrated by considering two model non-conservative systems: Lagrangian gas dynamics in non-conservative form and a form of isothermal Euler equations. Numerical experiments demonstrating...

We present an alternative framework for designing efficient numerical schemes for non-conservative hyperbolic systems. This approach is based on the design of entropy conservative discretizations and suitable numerical diffusion operators that mimic the effect of underlying viscous mechanisms. This approach is illustrated by considering two model non-conservative systems: Lagrangian gas dynamics in non-conservative form and a form of isothermal Euler equations. Numerical experiments demonstrating...

These lecture notes present adaptive multiresolution schemes for evolutionary PDEs in Cartesian geometries. The discretization schemes are based either on finite volume or finite difference schemes. The concept of multiresolution analyses, including Harten’s approach for point and cell averages, is described in some detail. Then the sparse point representation method is discussed. Different strategies for adaptive time-stepping, like local scale dependent time stepping and time step control, are...

We propose an adaptive finite element method for the solution of a coefficient inverse problem of simultaneous reconstruction of the dielectric permittivity and magnetic permeability functions in the Maxwell's system using limited boundary observations of the electric field in 3D. We derive a posteriori error estimates in the Tikhonov functional to be minimized and in the regularized solution of this functional, as well as formulate the corresponding adaptive algorithm. Our numerical experiments...

We discuss the occurrence of oscillations when using central schemes of the Lax-Friedrichs type (LFt), Rusanov’s method and the staggered and non-staggered second order Nessyahu-Tadmor (NT) schemes. Although these schemes are monotone or TVD, respectively, oscillations may be introduced at local data extrema. The dependence of oscillatory properties on the numerical viscosity coefficient is investigated rigorously for the LFt schemes, illuminating also the properties of Rusanov’s method. It turns...