Non-linear Petrov–Galerkin methods for reduced order modelling of the Navier–Stokes equations using a mixed finite element pair

Authors:

D. Xiao ^a,b, F. Fang ^b,J.Du ^d, C.C. Pain ^b, I.M. Navon ^c,* , A.G. Buchan ^b, A.H. ElSheikh ^b,G.Hu ^a
a State Key Laboratory of Geological Processes and Mineral Resources, China University of Geosciences,Wuhan 430074, China
^b Applied Modelling and Computation Group, Department of Earth Science and Engineering, Imperial College London, Prince Consort Road, London SW7 2AZ,UK
^c Department of Scienti？c Computing, Florida State University, Tallahassee, FL 32306-4120, USA
^d Institute of Atmospheric Physics, Chinese Academy of Sciences, Beijing 100029, China

Abstract：

A new nonlinear Petrov–Galerkin approach has been developed for proper orthogonal decomposition (POD) reduced order modelling (ROM) of the Navier–Stokes equations. The new method is based on the use of the cosine rule between the advection direction in Cartesian space–time and the direction of the gradient of the solution. A finite element pair, P1DGP2, which has good balance preserving properties is used here, consisting of a mix of discontinuous (for velocity components) and continuous (for pressure) basis functions. The contribution of the present paper lies in applying this new non-linear Petrov–Galerkin method to the reduced order Navier–Stokes equations, and thus improving the stability of ROM results without tuning parameters. The results of numerical tests are presented for a wind driven 2D gyre and the flow past a cylinder, which are simulated using the unstructured mesh finite element CFD model in order to illustrate the numerical performance of the method. The numerical results obtained show that the newly proposed POD Petrov–Galerkin method can provide more accurate and stable results than the POD Bubnov–Galerkin method.

Key words：

Finite element; Petrov–Galerkin; Navier–Stokes; Proper orthogonal decomposition; Discontinuous-Galerkin

Citation:

D. Xiao et al. Non-linear Petrov–Galerkin methods for reduced order modelling of the Navier–Stokes equations using a mixed finite element pair. Comput. Methods Appl. Mech. Engrg. 255 (2013) 147–157http://dx.doi.org/10.1016/j.cma.2012.11.002