Non-linear model reduction for the Navier-Stokes equations using residual DEIM method

  • Authors:
  • D. Xiao;F. Fang;A. G. Buchan;C. C. Pain;I. M. Navon;J. Du;G. Hu

  • Affiliations:
  • Applied Modelling and Computation Group, Department of Earth Science and Engineering, Imperial College London, Prince Consort Road, London SW7 2BP, UK11http://amcg.ese.imperial.ac.uk. and State Ke ...;Applied Modelling and Computation Group, Department of Earth Science and Engineering, Imperial College London, Prince Consort Road, London SW7 2BP, UK11http://amcg.ese.imperial.ac.uk.;Applied Modelling and Computation Group, Department of Earth Science and Engineering, Imperial College London, Prince Consort Road, London SW7 2BP, UK11http://amcg.ese.imperial.ac.uk.;Applied Modelling and Computation Group, Department of Earth Science and Engineering, Imperial College London, Prince Consort Road, London SW7 2BP, UK11http://amcg.ese.imperial.ac.uk.;Department of Scientific Computing, Florida State University, Tallahassee, FL 32306-4120, USA;Institute of Atmospheric Physics, Chinese Academy of Sciences, Beijing 100029, China;State Key Laboratory of Geological Processes and Mineral Resources, China University of Geosciences, Wuhan 430074, China

  • Venue:
  • Journal of Computational Physics
  • Year:
  • 2014

Quantified Score

Hi-index 31.45

Visualization

Abstract

This article presents a new reduced order model based upon proper orthogonal decomposition (POD) for solving the Navier-Stokes equations. The novelty of the method lies in its treatment of the equation's non-linear operator, for which a new method is proposed that provides accurate simulations within an efficient framework. The method itself is a hybrid of two existing approaches, namely the quadratic expansion method and the Discrete Empirical Interpolation Method (DEIM), that have already been developed to treat non-linear operators within reduced order models. The method proposed applies the quadratic expansion to provide a first approximation of the non-linear operator, and DEIM is then used as a corrector to improve its representation. In addition to the treatment of the non-linear operator the POD model is stabilized using a Petrov-Galerkin method. This adds artificial dissipation to the solution of the reduced order model which is necessary to avoid spurious oscillations and unstable solutions. A demonstration of the capabilities of this new approach is provided by solving the incompressible Navier-Stokes equations for simulating a flow past a cylinder and gyre problems. Comparisons are made with other treatments of non-linear operators, and these show the new method to provide significant improvements in the solution's accuracy.