Local POD Plus Galerkin Projection in the Unsteady Lid-Driven Cavity Problem

A local proper orthogonal decomposition (POD) plus Galerkin projection method is applied to the unsteady lid-driven cavity problem, namely the incompressible fluid flow in a two-dimensional box whose upper wall is moved back and forth at moderately large values of the Reynolds number. Such a method was recently developed for one-dimensional parabolic problems. Its extension to fluid dynamics problems is nontrivial (especially if rough CFD codes are used) and consists of using a computational fluid dynamics (CFD) code and a Galerkin system (GS) in a sequence of interspersed intervals $I_{CFD}$ and $I_{GS}$, respectively. The POD manifold is calculated retaining the most energetic POD modes resulting from the snapshots computed in the $I_{CFD}$ intervals; in fact, the POD manifold is completely calculated in the first $I_{CFD}$ interval but only updated in subsequent $I_{CFD}$ intervals. Intended to mimic industrial solvers, the CFD code contains unphysical terms that are introduced for purely numerical reasons to accelerate runs. The use of such a CFD code poses the question of which equations should be Galerkin projected onto the POD manifold, the exact governing equations or the CFD numerical scheme. Also, Galerkin projection can be made using either the standard $L_2$ inner product or a nonstandard one, based on a limited number of mesh points. After addressing these issues, a method is constructed that is able to accelerate the CFD code by a factor of the order of 5-15, depending on the Reynolds number and the nature (steady, periodic, or quasi-periodic) of the forcing velocity.

[1]  B. R. Noack,et al.  A hierarchy of low-dimensional models for the transient and post-transient cylinder wake , 2003, Journal of Fluid Mechanics.

[2]  I. Aranson,et al.  The world of the complex Ginzburg-Landau equation , 2001, cond-mat/0106115.

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

[4]  P. Sagaut,et al.  Towards an adaptive POD/SVD surrogate model for aeronautic design , 2011 .

[5]  Lucia Russo,et al.  On POD reduced models of tubular reactor with periodic regimes , 2008, Comput. Chem. Eng..

[6]  Arthur Veldman,et al.  Proper orthogonal decomposition and low-dimensional models for driven cavity flows , 1998 .

[7]  P. Holmes,et al.  The Proper Orthogonal Decomposition in the Analysis of Turbulent Flows , 1993 .

[8]  Jeffrey P. Thomas,et al.  Using Automatic Differentiation to Create a Nonlinear Reduced-Order-Model Aerodynamic Solver , 2010 .

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

[10]  Kelly Cohen,et al.  Low-dimensional modelling of a transient cylinder wake using double proper orthogonal decomposition , 2008, Journal of Fluid Mechanics.

[11]  Siep Weiland,et al.  Missing Point Estimation in Models Described by Proper Orthogonal Decomposition , 2004, IEEE Transactions on Automatic Control.

[12]  S. S. Ravindran,et al.  Reduced-Order Adaptive Controllers for Fluid Flows Using POD , 2000, J. Sci. Comput..

[13]  D. Rempfer.  LOW-DIMENSIONAL MODELING AND NUMERICAL SIMULATION OF TRANSITION IN SIMPLE SHEAR FLOWS , 2003 .

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

[15]  P. Sagaut,et al.  Calibrated reduced-order POD-Galerkin system for fluid flow modelling , 2005 .

[16]  J. Yorke,et al.  CHAOTIC ATTRACTORS IN CRISIS , 1982 .

[17]  O. G. Martynenko,et al.  Convective heat transfer , 1989 .

[18]  Juan J. Alonso,et al.  Investigation of non-linear projection for POD based reduced order models for Aerodynamics , 2001 .

[19]  José M. Vega,et al.  Reduced order models based on local POD plus Galerkin projection , 2010, J. Comput. Phys..

[20]  H. B. Keller,et al.  Driven cavity flows by efficient numerical techniques , 1983 .

[21]  Wei Shyy,et al.  Proper Orthogonal Decomposition for Time-Dependent Lid-Driven Cavity Flows , 2002 .

[22]  Stefan Siegel,et al.  FEEDBACK CONTROL OF A CIRCULAR CYLINDER WAKE IN EXPERIMENT AND SIMULATION (INVITED) , 2003 .

[23]  P. Duck.  Oscillatory flow inside a square cavity , 1982, Journal of Fluid Mechanics.

[24]  Haym Benaroya,et al.  Modeling Fluid Structure Interaction , 2000 .

[25]  G. Karniadakis,et al.  A spectral viscosity method for correcting the long-term behavior of POD models , 2004 .

[26]  A. Velazquez,et al.  A method to generate computationally efficient reduced order models , 2009 .

[27]  Gilead Tadmor,et al.  Mean field representation of the natural and actuated cylinder wake , 2010 .

[28]  Charbel Farhat,et al.  Reduced-order fluid/structure modeling of a complete aircraft configuration , 2006 .

[29]  Robert L. Street,et al.  An analysis and comparison of the time accuracy of fractional‐step methods for the Navier–Stokes equations on staggered grids , 2002 .

[30]  G. Karniadakis,et al.  Stability and accuracy of periodic flow solutions obtained by a POD-penalty method , 2005 .

[31]  U. Ghia,et al.  High-Re solutions for incompressible flow using the Navier-Stokes equations and a multigrid method , 1982 .

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

[33]  I. Kevrekidis,et al.  Equation-free/Galerkin-free POD-assisted computation of incompressible flows , 2005 .

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

[35]  D. Gottlieb,et al.  Numerical analysis of spectral methods : theory and applications , 1977 .