Non‐intrusive reduced‐order modelling of the Navier–Stokes equations based on RBF interpolation

We present a new non‐intrusive model reduction method for the Navier–Stokes equations. The method replaces the traditional approach of projecting the equations onto the reduced space with a radial basis function (RBF) multi‐dimensional interpolation. The main point of this method is to construct a number of multi‐dimensional interpolation functions using the RBF scatter multi‐dimensional interpolation method. The interpolation functions are used to calculate POD coefficients at each time step from POD coefficients at earlier time steps. The advantage of this method is that it does not require modifications to the source code (which would otherwise be very cumbersome), as it is independent of the governing equations of the system. Another advantage of this method is that it avoids the stability problem of POD/Galerkin. The novelty of this work lies in the application of RBF interpolation and POD to construct the reduced‐order model for the Navier–Stokes equations. Another novelty is the verification and validation of numerical examples (a lock exchange problem and a flow past a cylinder problem) using unstructured adaptive finite element ocean model. The results obtained show that CPU times are reduced by several orders of magnitude whilst the accuracy is maintained in comparison with the corresponding high‐fidelity models. Copyright © 2015 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.     

A POD-based reduced-order finite difference extrapolating model with fully second-order accuracy for non-stationary Stokes equations

A proper orthogonal decomposition (POD) reduced-order finite difference (FD) extrapolating model is established for the channel flow with local expansion denoted by non-stationary Stokes equations. The POD-based reduced-order numerical model to produce the solutions on the time span [T₀, T] (T₀ ? T) are obtained by extrapolation and iteration from the very short time span [0, T₀] information. The guides to choose the number of POD basis and renew POD basis are provided, and an implementation for solving the POD-based reduced-order FD extrapolating model is given. Some numerical experiments are used to show that the POD-based reduced-order FD extrapolating model is feasible and efficient for simulating the channel flow with local expansion.     

Discrete filter operators for large‐eddy simulation using high‐order spectral difference methods

The combination of a high‐order unstructured spectral difference (SD) spatial discretization scheme with sub‐grid scale (SGS) modeling for large‐eddy simulation is investigated with particular focus on the consistent implementation of a structural mixed model based on the scale similarity hypothesis. The difficult task of deriving a consistent formulation for the discrete filter within the SD element of arbitrary order led to the development of a new class of three‐dimensional constrained discrete filters. The discrete filters satisfy a set of selected criteria and are completely local within the SD element. Their weights can be automatically computed at run time from the number of solution points within each element and the expected filter cutoff length scale. The novel discrete filters can be applied to any SGS model involving explicit filtering and to a broad class of high‐order discontinuous finite element numerical schemes. The code is applied to the computation of turbulent channel flows at three Reynolds numbers, namely Re_τ = 180, 395, and 590 (based on the friction velocity u_τ and channel half‐width δ). Results from computations with and without the SGS model are compared against results from direct numerical simulation. The numerical experiments suggest that the results are sensitive to the use of the SGS model, even when a high‐order numerical scheme is used, especially when the grid resolution is kept relatively low and mostly in terms of resolved Reynolds stresses. Results obtained using existing filters based on the projection of the solution over lower‐order polynomial bases are also shown and demonstrate that these filters are inadequate for SGS modeling purposes, mostly because of their inability to enforce the selected cutoff length scale with sufficient accuracy. The use of the similarity mixed formulation proved to be particularly accurate in reproducing SGS interactions, confirming that its well‐known potential can be realized in conjunction with state‐of‐the‐art high‐order numerical schemes.Copyright © 2012 John Wiley & Sons, Ltd.     

A finite element variational multiscale method for incompressible flows based on the construction of the projection basis functions

In this article, we present a finite element variational multiscale (VMS) method for incompressible flows based on the construction of projection basis functions and compare it with common VMS method, which is defined by a low‐order finite element space L_h on the same grid as X_h for the velocity deformation tensor and a stabilization parameter α. The best algorithmic feature of our method is to construct the projection basis functions at the element level with minimal additional cost to replace the global projection operator. Finally, we give some numerical simulations of the nonlinear flow problems to show good stability and accuracy properties of the method. Copyright © 2011 John Wiley & Sons, Ltd.     

A reduced‐order approach for optimal control of fluids using proper orthogonal decomposition

In this article, a reduced‐order modeling approach, suitable for active control of fluid dynamical systems, based on proper orthogonal decomposition (POD) is presented. The rationale behind the reduced‐order modeling is that numerical simulation of Navier–Stokes equations is still too costly for the purpose of optimization and control of unsteady flows. The possibility of obtaining reduced‐order models that reduce the computational complexity associated with the Navier–Stokes equations is examined while capturing the essential dynamics by using the POD. The POD allows the extraction of a reduced set of basis functions, perhaps just a few, from a computational or experimental database through an eigenvalue analysis. The solution is then obtained as a linear combination of this reduced set of basis functions by means of Galerkin projection. This makes it attractive for optimal control and estimation of systems governed by partial differential equations (PDEs). It is used here in active control of fluid flows governed by the Navier–Stokes equations. In particular, flow over a backward‐facing step is considered. Reduced‐order models/low‐dimensional dynamical models for this system are obtained using POD basis functions (global) from the finite element discretizations of the Navier–Stokes equations. Their effectiveness in flow control applications is shown on a recirculation control problem using blowing on the channel boundary. Implementational issues are discussed and numerical experiments are presented. Copyright © 2000 John Wiley & Sons, Ltd.     

A new stable space–time formulation for two‐dimensional and three‐dimensional incompressible viscous flow

A space–time finite element method for the incompressible Navier–Stokes equations in a bounded domain in ?^d (with d=2 or 3) is presented. The method is based on the time‐discontinuous Galerkin method with the use of simplex‐type meshes together with the requirement that the space–time finite element discretization for the velocity and the pressure satisfy the inf–sup stability condition of Brezzi and Babu?ka. The finite element discretization for the pressure consists of piecewise linear functions, while piecewise linear functions enriched with a bubble function are used for the velocity. The stability proof and numerical results for some two‐dimensional problems are presented. Copyright © 2001 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.     

Orthogonality of modal bases in hp finite element models

In this paper, we exploit orthogonality of modal bases (SIAM J. Sci. Comput. 1999; 20 :1671–1695) used in hp finite element models. We calculate entries of coefficient matrix analytically without using any numerical integration, which can be computationally very expensive. We use properties of Jacobi polynomials and recast the entries of the coefficient matrix so that they can be evaluated analytically. We implement this in the context of the least‐squares finite element model although this procedure can be used in other finite element formulations. In this paper, we only develop analytical expressions for rectangular elements. Spectral convergence of the L₂ least‐squares functional is verified using exact solution of Kovasznay flow. Numerical results for transient flow over a backward‐facing step are also presented. We also solve steady flow past a circular cylinder and show the reduction in computational cost using expressions developed herein. Copyright © 2007 John Wiley & Sons, Ltd.     

A dual‐weighted trust‐region adaptive POD 4‐D Var applied to a finite‐volume shallow water equations model on the sphere

In this paper we study solutions of an inverse problem for a global shallow water model controlling its initial conditions specified from the 40‐yr ECMWF Re‐analysis (ERA‐40) data sets, in the presence of full or incomplete observations being assimilated in a time interval (window of assimilation) with or without background error covariance terms. As an extension of the work by Chen et al. (Int. J. Numer. Meth. Fluids 2009), we attempt to obtain a reduced order model of the above inverse problem, based on proper orthogonal decomposition (POD), referred to as POD 4D‐Var for a finite volume global shallow water equation model based on the Lin–Rood flux‐form semi‐Lagrangian semi‐implicit time integration scheme. Different approaches of POD implementation for the reduced inverse problem are compared, including a dual‐weighted method for snapshot selection coupled with a trust‐region POD adaptivity approach. Numerical results with various observational densities and background error covariance operator are also presented. The POD 4‐D Var model results combined with the trust‐region adaptivity exhibit similarity in terms of various error metrics to the full 4D Var results, but are obtained using a significantly lesser number of minimization iterations and require lesser CPU time. Based on our previous and current work, we conclude that POD 4‐D Var certainly warrants further studies, with promising potential of its extension to operational 3‐D numerical weather prediction models. Copyright © 2011 John Wiley & Sons, Ltd.     

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.     

A high‐order Petrov–Galerkin finite element method for the classical Boussinesq wave model

A high‐order Petrov–Galerkin finite element scheme is presented to solve the one‐dimensional depth‐integrated classical Boussinesq equations for weakly non‐linear and weakly dispersive waves. Finite elements are used both in the space and the time domains. The shape functions are bilinear in space–time, whereas the weighting functions are linear in space and quadratic in time, with C⁰‐continuity. Dispersion correction and a highly selective dissipation mechanism are introduced through additional streamline upwind terms in the weighting functions. An implicit, conditionally stable, one‐step predictor–corrector time integration scheme results. The accuracy and stability of the non‐linear discrete equations are investigated by means of a local Taylor series expansion. A linear spectral analysis is used for the full characterization of the predictor–corrector inner iterations. Based on the order of the analytical terms of the Boussinesq model and on the order of the numerical discretization, it is concluded that the scheme is fourth‐order accurate in terms of phase velocity. The dissipation term is third order only affecting the shortest wavelengths. A numerical convergence analysis showed a second‐order convergence rate in terms of both element size and time step. Four numerical experiments are addressed and their results are compared with analytical solutions or experimental data available in the literature: the propagation of a solitary wave, the oscillation of a flat bottom closed basin, the oscillation of a non‐flat bottom closed basin, and the propagation of a periodic wave over a submerged bar. Copyright © 2008 John Wiley & Sons, Ltd.     

A sensitivity analysis on the parameter of the GLS method for a second‐gradient theory of incompressible flow

Using a non‐conforming C⁰‐interior penalty method and the Galerkin least‐square approach, we develop a continuous–discontinuous Galerkin finite element method for discretizing fourth‐order incompressible flow problems. The formulation is weakly coercive for spaces that fail to satisfy the inf‐sup condition and consider discontinuous basis functions for the pressure field. We consider the results of a stability analysis through a lemma which indicates that there exists an optimal or quasi‐optimal least‐square stability parameter that depends on the polynomial degree used to interpolate the velocity and pressure fields, and on the geometry of the finite element in the mesh. We provide several numerical experiments illustrating such dependence, as well as the robustness of the method to deal with arbitrary basis functions for velocity and pressure, and the ability to stabilize large pressure gradients. We believe the results provided in this paper contribute for establishing a paradigm for future studies of the parameter of the Galerkin least square method for second‐gradient theory of incompressible flow problems. Copyright © 2015 John Wiley & Sons, Ltd.     

Reduced-order finite element method based on POD for fractional Tricomi-type equation

The reduced-order finite element method (FEM) based on a proper orthogonal decomposition (POD) theory is applied to the time fractional Tricomi-type equation. The present method is an improvement on the general FEM. It can significantly save memory space and effectively relieve the computing load due to its reconstruction of POD basis functions. Furthermore, the reduced-order finite element (FE) scheme is shown to be unconditionally stable, and error estimation is derived in detail. Two numerical examples are presented to show the feasibility and effectiveness of the method for time fractional differential equations (FDEs).     

2D Burgers equation with large Reynolds number using POD/DEIM and calibration

Model order reduction of the two‐dimensional Burgers equation is investigated. The mathematical formulation of POD/discrete empirical interpolation method (DEIM)‐reduced order model (ROM) is derived based on the Galerkin projection and DEIM from the existing high fidelity‐implicit finite‐difference full model. For validation, we numerically compared the POD ROM, POD/DEIM, and the full model in two cases of Re = 100 and Re = 1000, respectively. We found that the POD/DEIM ROM leads to a speed‐up of CPU time by a factor of O(10). The computational stability of POD/DEIM ROM is maintained by means of a careful selection of POD modes and the DEIM interpolation points. The solution of POD/DEIM in the case of Re = 1000 has an accuracy with error O(10^?3) versus O(10^?4) in the case of Re = 100 when compared with the high fidelity model. For this turbulent flow, a closure model consisting of a Tikhonov regularization is carried out in order to recover the missing information and is developed to account for the small‐scale dissipation effect of the truncated POD modes. It is shown that the computational results of this calibrated ROM exhibit considerable agreement with the high fidelity model, which implies the efficiency of the closure model used. Copyright © 2016 John Wiley & Sons, Ltd.     

An efficient high order direct ALE ADER finite volume scheme with a posteriori limiting for hydrodynamics and magnetohydrodynamics

In this paper, we present a new family of direct arbitrary–Lagrangian–Eulerian (ALE) finite volume schemes for the solution of hyperbolic balance laws on unstructured meshes in multiple space dimensions. The scheme is designed to be high‐order accurate both in space and time, and the mesh motion, which provides the new mesh configuration at the next time step, is taken into account in the final finite volume scheme that is based directly on a space‐time conservation formulation of the governing PDE system. To improve the computational efficiency of the algorithm, high order of accuracy in space is achieved using the a posteriori MOOD limiting strategy that allows the reconstruction procedure to be carried out with only one reconstruction stencil for any order of accuracy. We rely on an element‐local space‐time Galerkin finite element predictor on moving curved meshes to obtain a high‐order accurate one‐step time discretization, while the mesh velocity is computed by means of a suitable nodal solver algorithm that might also be supplemented with a local rezoning procedure to improve the mesh quality. Next, the old mesh configuration at time level tⁿ is connected to the new one at t^n + 1 by straight edges, hence providing unstructured space‐time control volumes, on the boundary of which the numerical flux has to be integrated. Here, we adopt a quadrature‐free integration, in which the space‐time boundaries of the control volumes are split into simplex sub‐elements that yield constant space‐time normal vectors and Jacobian matrices. In this way, the integrals over the simplex sub‐elements can be evaluated once and for all analytically during a preprocessing step. We apply the new high‐order direct ALE algorithm to the Euler equations of compressible gas dynamics (also referred to as hydrodynamics equations) as well as to the magnetohydrodynamics equations and we solve a set of classical test problems in two and three space dimensions. Numerical convergence rates are provided up to fifth order of accuracy in 2D and 3D for both hyperbolic systems considered in this paper. Finally, the efficiency of the new method is measured and carefully compared against the original formulation of the algorithm that makes use of a WENO reconstruction technique and Gaussian quadrature formulae for the flux integration: depending on the test problem, the new class of very efficient direct ALE schemes proposed in this paper can run up to ≈12 times faster in the 3D case. Copyright © 2016 John Wiley & Sons, Ltd.     

Spectral/hp penalty least‐squares finite element formulation for unsteady incompressible flows

In this paper, we present spectral/hp penalty least‐squares finite element formulation for the numerical solution of unsteady incompressible Navier–Stokes equations. Pressure is eliminated from Navier–Stokes equations using penalty method, and finite element model is developed in terms of velocity, vorticity and dilatation. High‐order element expansions are used to construct discrete form. Unlike other penalty finite element formulations, equal‐order Gauss integration is used for both viscous and penalty terms of the coefficient matrix. For time integration, space–time decoupled schemes are implemented. Second‐order accuracy of the time integration scheme is established using the method of manufactured solution. Numerical results are presented for impulsively started lid‐driven cavity flow at Reynolds number of 5000 and transient flow over a backward‐facing step. The effect of penalty parameter on the accuracy is investigated thoroughly in this paper and results are presented for a range of penalty parameter. Present formulation produces very accurate results for even very low penalty parameters (10–50). Copyright © 2008 John Wiley & Sons, Ltd.     

BOUNDARY CONDITIONS IN SHALLOW WATER MODELS—AN ALTERNATIVE IMPLEMENTATION FOR FINITE ELEMENT CODES

Finite element solution of the shallow water wave equations has found increasing use by researchers and practitioners in the modelling of oceans and coastal areas. Wave equation models, most of which use equal-orderC⁰ interpolants for both the velocity and the surface elevation, do not introduce spurious oscillation modes, hence avoiding the need for artificial or numerical damping. An important question for both primitive equation and wave equation models is the interpretation of boundary conditions. Analysis of the characteristics of the governing equations shows that for most geophysical flows a single condition at each boundary is sufficient, yet there is not a consensus in the literature as to what that boundary condition must be or how it should be implemented in a finite element code. Traditionally (partly because of limited data), surface elevation is specified at open ocean boundaries while the normal flux is specified as zero at land boundaries. In most finite element wave equation models both of these boundary conditions are implemented as essential conditions. Our recent work focuses on alternative ways to numerically implement normal flow boundary conditions with an eye towards improving the mass-conserving properties of wave equation models. A unique finite element formulation using generalized functions demonstrates that boundary conditions should be implemented by treating normal fluxes as natural conditions with the flux interpreted as external to the computational domain. Results from extensive numerical experiments show that the scheme does conserve mass for all parameter values. Furthermore, convergence studies demonstrate that the algorithm is consistent, as residual errors at the boundary diminish as the grid is refined.     

Experimental Study on the Pore Structure Fractals and Seepage Characteristics of a Coal Sample Around a Borehole in Coal Seam Water Infusion

In this study, a two-dimensional fully coupled computational model is developed for simulation of proppant settlement in hydro-fractures with the use of the extended finite element framework. The porous domain is governed by the well-known \((\mathbf{u}-p)\) formulation, which consists of the momentum balance equation of the bulk, in conjunction with the momentum balance and continuity equations of the pore fluid. The hydro-fracture inflow is modeled as a 1D flow on the basis of the Darcy law, in which fracture permeability is incorporated by means of the cubic law. Contact constraints are elaborated to eliminate the overlap of fracture edges and the leak-off flow. Proppant settlement is conducted on the basis of Stokes’ law for particle terminal velocity, in which the effects of fracture walls, concentration, viscosity and bridging are incorporated into the model. A fixed-point algorithm is introduced to achieve the optimum combination for the proppant injection. Using the extended finite element method, the strong discontinuity in the displacement field due to crack body, as well as the weak discontinuity in the pressure field due to leakage, is included in the model with the use of the Heaviside and modified level set enrichment functions, respectively. The robustness and versatility of the proposed numerical algorithm in determining the optimum proppant injection is examined through several numerical simulations.