On the stability of projection-based model order reduction for convection-dominated laminar and turbulent flows

Abstract In the literature on nonlinear projection-based model order reduction for computational fluid dynamics problems, it is often claimed that due to modal truncation, a projection-based reduced-order model (PROM) does not resolve the dissipative regime of the turbulent energy cascade and therefore is numerically unstable. Efforts at addressing this claim have ranged from attempting to model the effects of the truncated modes to enriching the classical subspace of approximation in order to account for the truncated phenomena. The objective of this paper is to challenge this claim. Exploring the relationship between projection-based model order reduction and semi-discretization and using numerical evidence from three relevant flow problems, this paper argues in an orderly manner that the real culprit behind most if not all reported numerical instabilities of PROMs for turbulence and convection-dominated turbulent flow problems is the Galerkin framework that has been used for constructing the PROMs. The paper also shows that alternatively, a Petrov-Galerkin framework can be used to construct numerically stable and accurate PROMs for convection-dominated laminar as well as turbulent flow problems, without resorting to additional closure models or tailoring of the subspace of approximation. It also shows that such alternative PROMs deliver significant speed-up factors.

[1]  Pierre Sagaut,et al.  Intermodal energy transfers in a proper orthogonal decomposition–Galerkin representation of a turbulent separated flow , 2003, Journal of Fluid Mechanics.

[2]  Benjamin Stamm,et al.  Model Order Reduction for Problems with Large Convection Effects , 2018, Computational Methods in Applied Sciences.

[3]  Charbel Farhat,et al.  Stabilization of projection‐based reduced‐order models , 2012 .

[4]  Charbel Farhat,et al.  Nonlinear model order reduction based on local reduced‐order bases , 2012 .

[5]  Mario Ohlberger,et al.  Nonlinear reduced basis approximation of parameterized evolution equations via the method of freezing , 2013 .

[6]  C. Farhat,et al.  Efficient non‐linear model reduction via a least‐squares Petrov–Galerkin projection and compressive tensor approximations , 2011 .

[7]  S. Osher,et al.  Efficient implementation of essentially non-oscillatory shock-capturing schemes,II , 1989 .

[8]  A. Kolmogorov,et al.  The local structure of turbulence in incompressible viscous fluid for very large Reynolds numbers , 1991, Proceedings of the Royal Society of London. Series A: Mathematical and Physical Sciences.

[9]  R. D. Richtmyer,et al.  A Method for the Numerical Calculation of Hydrodynamic Shocks , 1950 .

[10]  Pavel B. Bochev,et al.  Least-Squares Finite Element Methods , 2009, Applied mathematical sciences.

[11]  W. K. Anderson,et al.  Comparison of Finite Volume Flux Vector Splittings for the Euler Equations , 1985 .

[13]  Omer San,et al.  Learning-based robust stabilization for reduced-order models of 2D and 3D Boussinesq equations , 2017 .

[14]  Nadine Aubry,et al.  The dynamics of coherent structures in the wall region of a turbulent boundary layer , 1988, Journal of Fluid Mechanics.

[15]  R. Henderson.  Details of the drag curve near the onset of vortex shedding , 1995 .

[16]  Thomas J. R. Hughes,et al.  A new finite element formulation for computational fluid dynamics: IX. Fourier analysis of space-time Galerkin/least-squares algorithms , 1991 .

[17]  Traian Iliescu,et al.  Variational multiscale proper orthogonal decomposition: Navier‐stokes equations , 2012, 1210.7389.

[18]  Earl H. Dowell,et al.  Minimal subspace rotation on the Stiefel manifold for stabilization and enhancement of projection-based reduced order models for the compressible Navier-Stokes equations , 2015, J. Comput. Phys..

[19]  C. Farhat,et al.  Structure‐preserving, stability, and accuracy properties of the energy‐conserving sampling and weighting method for the hyper reduction of nonlinear finite element dynamic models , 2015 .

[20]  W. Ritz.  Über eine neue Methode zur Lösung gewisser Variationsprobleme der mathematischen Physik. , 1909 .

[21]  Earl H. Dowell,et al.  Low-dimensional modelling of high-Reynolds-number shear flows incorporating constraints from the Navier–Stokes equation , 2013, Journal of Fluid Mechanics.

[22]  B. V. Leer,et al.  Towards the ultimate conservative difference scheme V. A second-order sequel to Godunov's method , 1979 .

[23]  Alain Dervieux,et al.  A low-diffusion MUSCL scheme for LES on unstructured grids , 2004 .

[24]  C. Farhat,et al.  Dimensional reduction of nonlinear finite element dynamic models with finite rotations and energy‐based mesh sampling and weighting for computational efficiency , 2014 .

[25]  T. Hughes,et al.  The Galerkin/least-squares method for advective-diffusive equations , 1988 .

[26]  Stéphane Lanteri,et al.  Two-dimensional viscous flow computations on the Connection Machine: unstructured meshes, upwind schemes and massively parallel computations , 1993 .

[27]  Karen Willcox,et al.  Model Reduction for Large-Scale Systems with High-Dimensional Parametric Input Space , 2008, SIAM J. Sci. Comput..

[28]  T. Hughes,et al.  Streamline upwind/Petrov-Galerkin formulations for convection dominated flows with particular emphasis on the incompressible Navier-Stokes equations , 1990 .

[29]  Sebastian Grimberg,et al.  Hyperreduction of CFD Models of Turbulent Flows using a Machine Learning Approach , 2020 .

[30]  L. W.,et al.  The Theory of Sound , 1898, Nature.

[31]  Nirmal J. Nair,et al.  Transported snapshot model order reduction approach for parametric, steady‐state fluid flows containing parameter‐dependent shocks , 2018, International Journal for Numerical Methods in Engineering.

[32]  Fabien Casenave,et al.  Time Stable Reduced Order Modeling by an Enhanced Reduced Order Basis of the Turbulent and Incompressible 3D Navier–Stokes Equations , 2019, Mathematical and Computational Applications.

[33]  ShakibFarzin,et al.  A new finite element formulation for computational fluid dynamics , 1991 .

[34]  J. Debonis.  Solutions of the Taylor-Green Vortex Problem Using High-Resolution Explicit Finite Difference Methods , 2013 .

[35]  R. Murray,et al.  Model reduction for compressible flows using POD and Galerkin projection , 2004 .

[36]  J. Boris,et al.  Flux-Corrected Transport , 1997 .

[37]  L. Sirovich.  Turbulence and the dynamics of coherent structures. II. Symmetries and transformations , 1987 .

[38]  Nirmal J. Nair,et al.  Transported snapshot model order reduction approach for parametric, steady‐state fluid flows containing parameter‐dependent shocks , 2017, International Journal for Numerical Methods in Engineering.

[39]  Jean-Frédéric Gerbeau,et al.  Approximated Lax pairs for the reduced order integration of nonlinear evolution equations , 2014, J. Comput. Phys..

[40]  Gianluigi Rozza,et al.  A reduced order variational multiscale approach for turbulent flows , 2018, Advances in Computational Mathematics.

[41]  Earl H. Dowell,et al.  Stabilization of projection-based reduced order models of the Navier–Stokes , 2012 .

[42]  A. Jameson,et al.  Numerical solution of the Euler equations by finite volume methods using Runge Kutta time stepping schemes , 1981 .

[43]  Gianluigi Rozza,et al.  On a Certified Smagorinsky Reduced Basis Turbulence Model , 2017, SIAM J. Numer. Anal..

[44]  Gianluigi Rozza,et al.  Reduced Basis Approximation of Parametrized Advection-Diffusion PDEs with High Péclet Number , 2013, ENUMATH.

[45]  Jay P. Boris,et al.  Flux-corrected transport. I. SHASTA, a fluid transport algorithm that works , 1973 .

[46]  C. Wieselsberger.  New Data on the Laws of Fluid Resistance , 1922 .

[47]  Alain Dervieux,et al.  Mixed-element-volume MUSCL methods with weak viscosity for steady and unsteady flow calculations , 2000 .

[48]  Harbir Antil,et al.  Galerkin v. least-squares Petrov-Galerkin projection in nonlinear model reduction , 2015, J. Comput. Phys..

[49]  I. Daubechies,et al.  Iteratively reweighted least squares minimization for sparse recovery , 2008, 0807.0575.

[50]  Karen Willcox,et al.  Nonlinear Model Order Reduction via Lifting Transformations and Proper Orthogonal Decomposition , 2018, AIAA Journal.

[51]  Lawrence Sirovich,et al.  Karhunen–Loève procedure for gappy data , 1995 .

[52]  A. Pinkus.  n-Widths in Approximation Theory , 1985 .

[53]  Traian Iliescu,et al.  Proper orthogonal decomposition closure models for turbulent flows: A numerical comparison , 2011, 1106.3585.

[54]  Sebastian Grimberg,et al.  Mesh sampling and weighting for the hyperreduction of nonlinear Petrov–Galerkin reduced‐order models with local reduced‐order bases , 2020, International Journal for Numerical Methods in Engineering.

[55]  Charbel Farhat,et al.  Accelerated mesh sampling for the hyper reduction of nonlinear computational models , 2017 .

[56]  D. Rempfer,et al.  On Low-Dimensional Galerkin Models for Fluid Flow , 2000 .

[57]  A. W. Vreman.  An eddy-viscosity subgrid-scale model for turbulent shear flow: Algebraic theory and applications , 2004 .

[58]  Charbel Farhat,et al.  Progressive construction of a parametric reduced‐order model for PDE‐constrained optimization , 2014, ArXiv.

[59]  Gianluigi Rozza,et al.  Supremizer stabilization of POD–Galerkin approximation of parametrized steady incompressible Navier–Stokes equations , 2015 .

[60]  Thomas J. R. Hughes,et al.  Recent progress in the development and understanding of SUPG methods with special reference to the compressible Euler and Navier-Stokes equations†‡ , 1987 .

[61]  L. Sirovich.  Turbulence and the dynamics of coherent structures. I. Coherent structures , 1987 .

[62]  Benjamin Peherstorfer,et al.  Localized Discrete Empirical Interpolation Method , 2014, SIAM J. Sci. Comput..

[63]  Charbel Farhat,et al.  A Minimum Overlap Restricted Additive Schwarz Preconditioner and Applications in 3D Flow Simulations , 1998 .

[64]  Rémi Abgrall,et al.  Model reduction usingL1-norm minimization as an application to nonlinear hyperbolic problems: Model reduction for hyperbolic problems , 2018 .

[65]  D. Rovas,et al.  Reliable Real-Time Solution of Parametrized Partial Differential Equations: Reduced-Basis Output Bound Methods , 2002 .

[66]  Charbel Farhat,et al.  The GNAT method for nonlinear model reduction: Effective implementation and application to computational fluid dynamics and turbulent flows , 2012, J. Comput. Phys..

[67]  Charles-Henri Bruneau,et al.  Enablers for robust POD models , 2009, J. Comput. Phys..

[68]  Gianluigi Rozza,et al.  Reduced basis approximation and a posteriori error estimation for Stokes flows in parametrized geometries: roles of the inf-sup stability constants , 2013, Numerische Mathematik.

[69]  B. R. Noack,et al.  On the need for a nonlinear subscale turbulence term in POD models as exemplified for a high-Reynolds-number flow over an Ahmed body , 2013, Journal of Fluid Mechanics.

[70]  Gianluigi Rozza,et al.  Stabilized Weighted Reduced Basis Methods for Parametrized Advection Dominated Problems with Random Inputs , 2017, SIAM/ASA J. Uncertain. Quantification.

[71]  Kookjin Lee,et al.  Model reduction of dynamical systems on nonlinear manifolds using deep convolutional autoencoders , 2018, J. Comput. Phys..

[72]  Per-Olof Persson,et al.  Stage-parallel fully implicit Runge-Kutta solvers for discontinuous Galerkin fluid simulations , 2017, J. Comput. Phys..

[73]  Hans Petter Langtangen,et al.  High performance Python for direct numerical simulations of turbulent flows , 2016, Comput. Phys. Commun..

[74]  Bernard Haasdonk,et al.  Model reduction of parametrized evolution problems using the reduced basis method with adaptive time partitioning , 2011 .

[75]  Alfio Quarteroni,et al.  An algebraic least squares reduced basis method for the solution of nonaffinely parametrized Stokes equations , 2019, Computer Methods in Applied Mechanics and Engineering.

[76]  Danny C. Sorensen,et al.  Nonlinear Model Reduction via Discrete Empirical Interpolation , 2010, SIAM J. Sci. Comput..

[77]  Rémi Abgrall,et al.  High‐order CFD methods: current status and perspective , 2013 .

[78]  Charbel Farhat,et al.  On the Use of Discrete Nonlinear Reduced-Order Models for the Prediction of Steady-State Flows Past Parametrically Deformed Complex Geometries , 2016 .

[79]  Matthew F. Barone,et al.  Stable Galerkin reduced order models for linearized compressible flow , 2009, J. Comput. Phys..