2009, Bartels, Sören, Müller, Rüdiger

Optimal a posteriori error estimates in L∞(L2) are derived for the finite element approximation of Allen-Cahn equations. The estimates depend on the inverse of a small parameter only in a low order polynomial and are valid past topological changes of the evolving interface. The error analysis employs an elliptic reconstruction of the approximate solution and applies to a large class of conforming, nonconforming, mixed, and discontinuous Galerkin methods. Numerical experiments illustrate the theoretical results.

2014, Eigel, Martin, Samrowski, Tatiana

A functional type a posteriori error estimator for the finite element discretisation of the stationary reaction-convection-diffusion equation is derived. In case of dominant convection, the solution for this class of problems typically exhibits boundary layers and shock-front like areas with steep gradients. This renders the accurate numerical solution very demanding and appropriate techniques for the adaptive resolution of regions with large approximation errors are crucial. Functional error estimators as derived here contain no mesh-dependent constants and provide guaranteed error bounds for any conforming approximation. To evaluate the error estimator, a minimisation problem is solved which does not require any Galerkin orthogonality or any specific properties of the employed approximation space. Based on a set of numerical examples, we assess the performance of the new estimator. It is observed that it exhibits a good efficiency also with convection-dominated problem settings.

2010, Ganesan, Sashikumaar

We present a heterogeneous finite element approximation of the solution of a population balance equation, which depends both the physical and internal property coordinates. We employ the operator-splitting method to split the high-dimensional population balance equation into two low-dimensional equations, and discretize the low-dimensional equations separately. In particular, we discretize the physical and internal spaces with the standard Galerkin and Streamline Upwind Petrov Galerkin (SUPG) finite elements, respectively. It is demonstrated that the variational form of the operator-split population balance equation is equivalent to the variational form of the standard equation up to a perturbation term of order $tau^2$ in the backward Euler scheme, where $tau$ is a time step. Further, the stability and error estimates have been derived for the heterogeneous finite element discretization scheme applied to the population balance equation. It is shown that a slightly more regularity, $i.e,$ the mixed partial derivatives of the solution has to be bounded, is necessary for the solution of the population balance equation with the operator-splitting finite element method. Numerical results are presented to demonstrate the accuracy of the numerical scheme.

2021, John, Volker, Moreau, Baptiste, Novo, Julia

A reduced order model (ROM) method based on proper orthogonal decomposition (POD) is analyzed for convection-diffusion-reaction equations. The streamline-upwind Petrov--Galerkin (SUPG) stabilization is used in the practically interesting case of dominant convection, both for the full order method (FOM) and the ROM simulations. The asymptotic choice of the stabilization parameter for the SUPG-ROM is done as proposed in the literature. This paper presents a finite element convergence analysis of the SUPG-ROM method for errors in different norms. The constants in the error bounds are uniform with respect to small diffusion coefficients. Numerical studies illustrate the performance of the SUPG-ROM method.

2010, Bartels, Sören, Müller, Rüdiger

A fully computable upper bound for the finite element approximation error of Allen-Cahn and Cahn-Hilliard equations with logarithmic potentials is derived. Numerical experiments show that for the sharp interface limit this bound is robust past topological changes. Modifications of the abstract results to derive quasi-optimal error estimates in different norms for lowest order finite element methods are discussed and lead to weaker conditions on the residuals under which the conditional error estimates hold.

2016, Carstensen, Carsten, Eigel, Martin

A hierarchical a posteriori error estimator for the first-order finite element method (FEM) on a red-refined triangular mesh is presented for the 2D Poisson model problem. Reliability and efficiency with some explicit constant is proved for triangulations with inner angles smaller than or equal to π/2 . The error estimator does not rely on any saturation assumption and is valid even in the pre-asymptotic regime on arbitrarily coarse meshes. The evaluation of the estimator is a simple post-processing of the piecewise linear FEM without any extra solve plus a higher-order approximation term. The results also allows the striking observation that arbitrary local averaging of the primal variable leads to a reliable and efficient error estimation. Several numerical experiments illustrate the performance of the proposed a posteriori error estimator for computational benchmarks.

2020, Eigel, Martin, Trunschke, Philipp, Schneider, Reinhold

We consider best approximation problems in a nonlinear subset of a Banach space of functions. The norm is assumed to be a generalization of the L2 norm for which only a weighted Monte Carlo estimate can be computed. The objective is to obtain an approximation of an unknown target function by minimizing the empirical norm. In the case of linear subspaces it is well-known that such least squares approximations can become inaccurate and unstable when the number of samples is too close to the number of parameters. We review this statement for general nonlinear subsets and establish error bounds for the empirical best approximation error. Our results are based on a restricted isometry property (RIP) which holds in probability and we show sufficient conditions for the RIP to be satisfied with high probability. Several model classes are examined where analytical statements can be made about the RIP. Numerical experiments illustrate some of the obtained stability bounds.

2018, Barrenechea, Gabriel R., John, Volker, Knobloch, Petr, Rankin, Richard

Recent results on the numerical analysis of Algebraic Flux Correction (AFC) finite element schemes for scalar convection-diffusion equations are reviewed and presented in a unified way. A general form of the method is presented using a link between AFC schemes and nonlinear edge-based diffusion scheme. Then, specific versions of the method, this is, different definitions for the flux limiters, are reviewed and their main results stated. Numerical studies compare the different versions of the scheme.

2006, Daoud, Daoud, Geiser, Jürgen

In this paper we consider the first order fractional splitting method to solve decomposed complex equations with multi-physical processes for applications in porous media and phase-transitions. The first order fractional splitting method is also considered as basic solution for the overlapping Schwarz-Waveform-Relaxation method for an overlapped subdomains. The accuracy and the efficiency of the methods are investigated through the solution of different model problems of scalar, coupling and decoupling systems of convection reaction diffusion equation.

2019, John, Volker, Knobloch, Petr, Wilbrandt, Ulrich

Discretizations of incompressible flow problems with pairs of finite element spaces that do not satisfy a discrete inf-sup condition require a so-called pressure stabilization. This paper gives an overview and systematic assessment of stabilized methods, including the respective error analysis.