A streamline derivative POD-ROM for advection-diffusion-reaction equations

We introduce a new streamline derivative projection-based closure modeling strategy for the numerical stabilization of Proper Orthogonal Decomposition-Reduced Order Models (PODROM). As a first preliminary step, the proposed model is analyzed and tested for advection-dominated advection-diffusion-reaction equations. In this framework, the numerical analysis for the Finite Element (FE) discretization of the proposed new POD-ROM is presented, by mainly deriving the corresponding error estimates. Numerical tests for advection-dominated regime show the efficiency of the proposed method, as well the increased accuracy over the standard POD-ROM that discovers its well-known limitations very soon in the numerical settings considered, i.e. for low diffusion coefficients.Nous introduisons une nouvelle stratégie de modélisation de type streamline derivative basée sur projection pour la stabilisation numérique de modèles d’ordre réduit de type POD (PODROM). Comme première étape préliminaire, le modèle proposé est analysé et testé pour les équations d’advection-diffusion-réaction dominées par l’advection. Dans ce cadre, l’analyse numérique de la discrétisation par éléments finis (FE) du nouveau POD-ROM proposé est présentée, en dérivant principalement les estimations d’erreur correspondantes. Des tests numériques pour le régime dominé par l’advection montrent l’efficacité de la méthode proposée, ainsi que la précision accrue par rapport à la méthode POD-ROM standard qui d´ecouvre très rapidement ses limites bien connues dans le cas des paramètres numériques considérés, c’est-à-dire pour de faibles coefficients de diffusion

Robust error estimates in weak norms for advection dominated transport problems with rough data

We consider mixing problems in the form of transient convection--diffusion equations with a velocity vector field with multiscale character and rough data. We assume that the velocity field has two scales, a coarse scale with slow spatial variation, which is responsible for advective transport and a fine scale with small amplitude that contributes to the mixing. For this problem we consider the estimation of filtered error quantities for solutions computed using a finite element method with symmetric stabilization. A posteriori error estimates and a priori error estimates are derived using the multiscale decomposition of the advective velocity to improve stability. All estimates are independent both of the P\'eclet number and of the regularity of the exact solution

An adaptive fixed-mesh ALE method for free surface flows

In this work we present a Fixed-Mesh ALE method for the numerical simulation of free surface flows capable of using an adaptive finite element mesh covering a background domain. This mesh is successively refined and unrefined at each time step in order to focus the computational effort on the spatial regions where it is required. Some of the main ingredients of the formulation are the use of an Arbitrary-Lagrangian–Eulerian formulation for computing temporal derivatives, the use of stabilization terms for stabilizing convection, stabilizing the lack of compatibility between velocity and pressure interpolation spaces, and stabilizing the ill-conditioning introduced by the cuts on the background finite element mesh, and the coupling of the algorithm with an adaptive mesh refinement procedure suitable for running on distributed memory environments. Algorithmic steps for the projection between meshes are presented together with the algebraic fractional step approach used for improving the condition number of the linear systems to be solved. The method is tested in several numerical examples. The expected convergence rates both in space and time are observed. Smooth solution fields for both velocity and pressure are obtained (as a result of the contribution of the stabilization terms). Finally, a good agreement between the numerical results and the reference experimental data is obtained.Postprint (published version

Variational Multiscale Stabilization and the Exponential Decay of Fine-scale Correctors

This paper addresses the variational multiscale stabilization of standard finite element methods for linear partial differential equations that exhibit multiscale features. The stabilization is of Petrov-Galerkin type with a standard finite element trial space and a problem-dependent test space based on pre-computed fine-scale correctors. The exponential decay of these correctors and their localisation to local cell problems is rigorously justified. The stabilization eliminates scale-dependent pre-asymptotic effects as they appear for standard finite element discretizations of highly oscillatory problems, e.g., the poor $L^2$ approximation in homogenization problems or the pollution effect in high-frequency acoustic scattering