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.

2016, Ahmed, Naveed, John, Volker, Matthies, Gunar, Novo, Julia

The local projection stabilization (LPS) method in space is considered to approximate the evolutionary Oseen equations. Optimal error bounds independent of the viscosity parameter are obtained in the continuous-in-time case for the approximations of both velocity and pressure. In addition, the fully discrete case in combination with higher order continuous Galerkin-Petrov (cGP) methods is studied. Error estimates of order k + 1 are proved, where k denotes the polynomial degree in time, assuming that the convective term is time-independent. Numerical results show that the predicted order is also achieved in the general case of time-dependent convective terms.

2015, John, Volker, Kaiser, Kristine, Novo, Julia

Finite element error estimates are derived for the incompressible Stokes equations with variable viscosity. The ratio of the supremum and the infimum of the viscosity appears in the error bounds. Numerical studies show that this ratio can be observed sometimes. However, often the numerical results show a weaker dependency on the viscosity.

2015, Frutos, Javier de, John, Volker, Novo, Julia

Weighted essentially non-oscillatory (WENO) finite difference schemes have been recommended in a competitive study of discretizations for scalar evolutionary convectiondiffusion equations [?]. This paper explores the applicability of these schemes for the simulation of incompressible flows. To this end, WENO schemes are used in several nonincremental and incremental projection methods for the incompressible Navier-Stokes equations. Velocity and pressure are discretized on the same grid. A pressure stabilization Petrov-Galerkin (PSPG) type of stabilization is introduced in the incremental schemes to account for the violation of the discrete inf-sup condition. Algorithmic aspects of the proposed schemes are discussed. The schemes are studied on several examples with different features. It is shown that the WENO finite difference idea can be transferred to the simulation of incompressible flows. Some shortcomings of the methods, which are due to the splitting in projection schemes, become also obvious.

2017, John, Volker, Knobloch, Petr, Novo, Julia

The contents of this paper is twofold. First, important recent results concerning finite element methods for convection-dominated problems and incompressible flow problems are described that illustrate the activities in these topics. Second, a number of, in our opinion, important problems in these fields are discussed.

2013, de Frutos, Javier, García-Archilla, Bosco, John, Volker, Novo, Julia

An adaptive algorithm for the numerical simulation of time-dependent convectiondiffusion-reaction equations will be proposed and studied. The algorithm allows the use of the natural extension of any error estimator for the steady-state problem for controlling local refinement and coarsening. The main idea consists in considering the SUPG solution of the evolutionary problem as the SUPG solution of a particular steady-state convectiondiffusion problem with data depending on the computed solution. The application of the error estimator is based on a heuristic argument by considering a certain term to be of higher order. This argument is supported in the one-dimensional case by numerical analysis. In the numerical studies, particularly the residual-based error estimator from [18] will be applied, which has proved to be robust in the SUPG norm. The effectivity of this error estimator will be studied and the numerical results (accuracy of the solution, fineness of the meshes) will be compared with results obtained by utilizing the adaptive algorithm proposed in [6]

2011, John, Volker, Novo, Julia

Finite element and finite difference discretizations for evolutionary convection-diffusion-reaction equations in two and three dimensions are studied which give solutions without or with small under- and overshoots. The studied methods include a linear and a nonlinear FEM-FCT scheme, simple upwinding, an ENO scheme of order 3, and a fifth order WENO scheme. Both finite element methods are combined with the Crank--Nicolson scheme and the finite difference discretizations are coupled with explicit total variation diminishing Runge--Kutta methods. An assessment of the methods with respect to accuracy, size of under- and overshoots, and efficiency is presented, in the situation of a domain which is a tensor product of intervals and of uniform grids in time and space. Some comments to the aspects of adaptivity and more complicated domains are given. The obtained results lead to recommendations concerning the use of the methods.

2010, John, Volker, Novo, Julia

Conditions on the stabilization parameters are explored for different approaches in deriving error estimates for the SUPG finite element stabilization of time-dependent convection-diffusion-reaction equations that is combined with the backward Euler method. Standard energy arguments lead to estimates for stabilization parameters that depend on the length of the time step. The stabilization vanishes in the time-continuous limit. However, based on numerical experiences, this seems not to be the correct behavior. For this reason, the time-continuous case is analyzed under certain conditions on the coefficients of the equation and the finite element method. An error estimate with the standard order of convergence is derived for stabilization parameters of the same form that is optimal for the steady-state problem. Numerical studies support the analytical results.

2013, John, Volker, Novo, Julia

Optimal error estimates for the pressure stabilized Petrov-Galerkin (PSPG) method for the continuous-in-time discretization of the evolutionary Stokes equations are proved in the case of regular solutions. The main result is applicable to higher order finite elements. The error bounds for the pressure depend on the error of the pressure at the initial time. An approach is suggested for choosing the discrete initial velocity in such a way that this error is bounded. The "instability of the discrete pressure for small time steps", which is reported in the literature, is discussed on the basis of the analytical results. Numerical studies confirm the theoretical results, showing in particular that this instability does not occur for the proposed initial condition.

2015, Frutos, Javier de, García-Archilla, Bosco, John, Volker, Novo, Julia

The approximation of the time-dependent Oseen problem using inf-sup stable mixed finite elements in a Galerkin method with grad-div stabilization is studied. The main goal is to prove that adding a grad-div stabilization term to the Galerkin approximation has a stabilizing effect for small viscosity. Both the continuous-in-time and the fully discrete case (backward Euler method, the two-step BDF, and CrankNicolson schemes) are analyzed. In fact, error bounds are obtained that do not depend on the inverse of the viscosity in the case where the solution is sufficiently smooth. The bounds for the divergence of the velocity as well as for the pressure are optimal. The analysis is based on the use of a specific Stokes projection. Numerical studies support the analytical results