A high‐fidelity low‐cost aerodynamic model using proper orthogonal decomposition

In this paper a high‐fidelity aerodynamic model is presented for use in parametric studies of weapon aerodynamics. The method employs a reduced‐order model obtained from the proper orthogonal decomposition (POD) of an ensemble of computational fluid dynamics (CFD) solutions with varying parameters. This decomposition produces an optimal linear set of orthogonal basis functions that best describe the ensemble of numerical solutions. These solutions are then projected onto this set of basis functions to provide a finite set of scalar coefficients that represent the solutions. A pseudo‐continuous representation of these projection coefficients is constructed, which allows predictions to be made of parameter combinations not in the original set of observations. The paper explores the performance of a few design‐of‐experiment approaches for the generation of the initial ensemble of computational experiments. Response surface construction methods based on parametric and non‐parametric models for the pseudo‐continuous representation of the projection coefficients are also evaluated. The model has been applied to two‐flow problems related to high‐speed weapon aerodynamics, inviscid flow around a flare‐stabilized hypersonic projectile and supersonic turbulent flow around a fin‐stabilized projectile with drooping nose control. Comparisons of model predictions with high‐fidelity CFD simulations suggest that the POD provides a reliable and robust approach to the construction of reduced‐order models. The practicality of the model is shown to be sensitive to the technique used to generate the ensemble of observations from which the model is constructed, while the accuracy of the approach depends on the pseudo‐continuous representation of the projection coefficients. Copyright © 2009 John Wiley & Sons, Ltd.     

An optimizing reduced order FDS for the tropical Pacific Ocean reduced gravity model

Proper orthogonal decomposition (POD) and singular value decomposition (SVD) methods are used to study a finite difference discretization scheme (FDS) for the tropical Pacific Ocean reduced gravity model. Ensembles of data are compiled from transient solutions computed from the discrete equation system derived by FDS for the tropical Pacific Ocean reduced gravity model. The optimal orthogonal bases are used to reconstruct the elements of the ensemble with POD and SVD. Combining the above approach with a Galerkin projection procedure yields a new optimizing FDS model of lower dimensions and high accuracy for the tropical Pacific Ocean reduced gravity model. An error estimate of the new reduced order optimizing FDS model is then derived. Numerical examples are presented illustrating that the error between the POD approximate solution and the full FDS solution is consistent with previously obtained theoretical results, thus validating the feasibility and efficiency of POD method. Copyright © 2007 John Wiley & Sons, Ltd.     

An optimizing implicit difference scheme based on proper orthogonal decomposition for the two‐dimensional unsaturated soil water flow equation

An optimizing reduced implicit difference scheme (IDS) based on singular value decomposition (SVD) and proper orthogonal decomposition (POD) for the two‐dimensional unsaturated soil water flow equation is presented. An ensemble of snapshots is compiled from the transient solutions derived from the usual IDS for a two‐dimensional unsaturated flow equation. Then, optimal orthogonal bases are reconstructed by implementing SVD and POD techniques for the ensemble of snapshots. Combining POD with a Galerkin projection approach, a new lower dimensional and highly accurate IDS for the two‐dimensional unsaturated flow equation is obtained. Error estimates between the true solution, the usual IDS solution, and the reduced IDS solution based on POD basis are derived. Finally, it is shown by means of a numerical example using the technology of local refined grids that the computational load is greatly diminished by using the reduced IDS. Also, the error between the POD approximate solution and the usual IDS solution is proved to be consistent with the derived theoretical results. Thus, both feasibility and efficiency of the POD method are validated. Copyright © 2011 John Wiley & Sons, Ltd.     

A dual‐weighted trust‐region adaptive POD 4D‐VAR applied to a finite‐element shallow‐water equations model

We consider a limited‐area finite‐element discretization of the shallow‐water equations model. Our purpose in this paper is to solve an inverse problem for the above model controlling its initial conditions in presence of observations being assimilated in a time interval (window of assimilation). We then attempt to obtain a reduced‐order model of the above inverse problem, based on proper orthogonal decomposition (POD), referred to as POD 4‐D VAR. Different approaches of POD implementation of the reduced inverse problem are compared, including a dual‐weighed method for snapshot selection coupled with a trust‐region POD approach. Numerical results obtained point to an improved accuracy in all metrics tested when dual‐weighing choice of snapshots is combined with POD adaptivity of the trust‐region type. Results of ad‐hoc adaptivity of the POD 4‐D VAR turn out to yield less accurate results than trust‐region POD when compared with high‐fidelity model. Directions of future research are finally outlined. Copyright © 2009 John Wiley & Sons, Ltd.     

Non‐intrusive reduced order modelling with least squares fitting on a sparse grid

This paper presents a non‐intrusive reduced order model for general, dynamic partial differential equations. Based upon proper orthogonal decomposition (POD) and Smolyak sparse grid collocation, the method first projects the unknowns with full space and time coordinates onto a reduced POD basis. Then we introduce a new least squares fitting procedure to approximate the dynamical transition of the POD coefficients between subsequent time steps, taking only a set of full model solution snapshots as the training data during the construction. Thus, neither the physical details nor further numerical simulations of the original PDE model are required by this methodology, and the level of non‐intrusiveness is improved compared with existing reduced order models. Furthermore, we take adaptive measures to address the instability issue arising from reduced order iterations of the POD coefficients. This model can be applied to a wide range of physical and engineering scenarios, and we test it on a couple of problems in fluid dynamics. It is demonstrated that this reduced order approach captures the dominant features of the high‐fidelity models with reasonable accuracy while the computation complexity is reduced by several orders of magnitude. Copyright © 2016 John Wiley & Sons, Ltd.     

Reduced‐order modeling based on POD of a parabolized Navier–Stokes equation model I: forward model

A proper orthogonal decomposition (POD)‐based reduced‐order model of the parabolized Navier–Stokes (PNS) equations is derived in this article. A space‐marching finite difference method with time relaxation is used to obtain the solution of this problem, from which snapshots are obtained to generate the POD basis functions used to construct the reduced‐order model. In order to improve the accuracy and the stability of the reduced‐order model in the presence of a high Reynolds number, we applied a Sobolev H₁ norm calibration to the POD construction process. Finally, some numerical tests with a high‐fidelity model as well as the POD reduced‐order model were carried out to demonstrate the efficiency and the accuracy of the reduced‐order model for solving the PNS equations compared with the full PNS model. Different inflow conditions and different selections of snapshots were experimented to test the POD reduction technique. The efficiency of the H₁ norm POD calibration is illustrated for the PNS model with increasingly higher Reynolds numbers, along with the optimal dissipation coefficient derivation, yielding the best root mean square error and correlation coefficient between the full and reduced‐order PNS models. Copyright © 2011 John Wiley & Sons, Ltd.     

Interpolation-based reduced-order modelling for steady transonic flows via manifold learning

This paper presents a parametric reduced-order model (ROM) based on manifold learning (ML) for use in steady transonic aerodynamic applications. The main objective of this work is to derive an efficient ROM that exploits the low-dimensional nonlinear solution manifold to ensure an improved treatment of the nonlinearities involved in varying the inflow conditions to obtain an accurate prediction of shocks. The reduced-order representation of the data is derived using the Isomap ML method, which is applied to a set of sampled computational fluid dynamics (CFD) data. In order to develop a ROM that has the ability to predict approximate CFD solutions at untried parameter combinations, Isomap is coupled with an interpolation method to capture the variations in parameters like the angle of attack or the Mach number. Furthermore, an approximate local inverse mapping from the reduced-order representation to the full CFD solution space is introduced. The proposed ROM, called Isomap+I, is applied to the two-dimensional NACA 64A010 airfoil and to the 3D LANN wing. The results are compared to those obtained by proper orthogonal decomposition plus interpolation (POD+I) and to the full-order CFD model.     

Principal interval decomposition framework for POD reduced‐order modeling of convective Boussinesq flows

A principal interval decomposition (PID) approach is presented for the reduced‐order modeling of unsteady Boussinesq equations. The PID method optimizes the lengths of the time windows over which proper orthogonal decomposition (POD) is performed and can be highly effective in building reduced‐order models for convective problems. The performance of these POD models with and without using the PID approach is investigated by applying these methods to the unsteady lock‐exchange flow problem. This benchmark problem exhibits a strong shear flow induced by a temperature jump and results in the Kelvin–Helmholtz instability. This problem is considered a challenging benchmark problem for the development of reduced‐order models. The reference solutions are obtained by direct numerical simulations of the vorticity and temperature transport equations using a compact fourth‐order‐accurate scheme. We compare the accuracy of reduced‐order models developed with different numbers of POD basis functions and different numbers of principal intervals. A linear interpolation model is constructed to obtain basis functions when varying physical parameters. The predictive performance of our models is then analyzed over a wide range of Reynolds numbers. It is shown that the PID approach provides a significant improvement in accuracy over the standard Galerkin POD reduced‐order model. This numerical assessment of the PID shows that it may represent a reliable model reduction tool for convection‐dominated, unsteady‐flow problems. Copyright © 2015 John Wiley & Sons, Ltd.     

Reduced‐order modelling of an adaptive mesh ocean model

A novel proper orthogonal decomposition (POD) model has been developed for use with an advanced unstructured mesh finite‐element ocean model, the Imperial College Ocean Model (ICOM, described in detail below), which includes many recent developments in ocean modelling and numerical analysis. The advantages of the POD model developed here over existing POD approaches are the ability:

• 1. To increase accuracy when representing geostrophic balance (the balance between the Coriolis terms and the pressure gradient). This is achieved through the use of two sets of geostrophic basis functions where each one is calculated by basis functions for velocities u and v.
• 2. To speed up the POD simulation. To achieve this a new numerical technique is introduced, whereby a time‐dependent matrix in the discretized equation is rapidly constructed from a series of time‐independent matrices. This development imparts considerable efficiency gains over the often‐used alternative of calculating each finite element over the computational domain at each time level.
• 3. To use dynamically adaptive meshes in the above POD model.

Copyright © 2008 John Wiley & Sons, Ltd.