Reduced-order modeling of transonic flows around an airfoil submitted to small deformations

A reduced-order model (ROM) is developed for the prediction of unsteady transonic flows past an airfoil submitted to small deformations, at moderate Reynolds number. Considering a suitable state formulation as well as a consistent inner product, the Galerkin projection of the compressible flow Navier–Stokes equations, the high-fidelity (HF) model, onto a low-dimensional basis determined by Proper Orthogonal Decomposition (POD), leads to a polynomial quadratic ODE system relevant to the prediction of main flow features. A fictitious domain deformation technique is yielded by the Hadamard formulation of HF model and validated at HF level. This approach captures airfoil profile deformation by a modification of the boundary conditions whereas the spatial domain remains unchanged. A mixed POD gathering information from snapshot series associated with several airfoil profiles can be defined. The temporal coefficients in POD expansion are shape-dependent while spatial POD modes are not. In the ROM, airfoil deformation is introduced by a steady forcing term. ROM reliability towards airfoil deformation is demonstrated for the prediction of HF-resolved as well as unknown intermediate configurations.     

The Elastoplast Discontinuous Galerkin (EDG) method for the Navier–Stokes equations

The present work details the Elastoplast (this name is a translation from the French “sparadrap”, a concept first applied by Yves Morchoisne for Spectral methods [1]) Discontinuous Galerkin (EDG) method to solve the compressible Navier–Stokes equations. This method was first presented in 2009 at the ICOSAHOM congress with some Cartesian grid applications. We focus here on unstructured grid applications for which the EDG method seems very attractive. As in the Recovery method presented by van Leer and Nomura in 2005 for diffusion, jumps across element boundaries are locally eliminated by recovering the solution on an overlapping cell. In the case of Recovery, this cell is the union of the two neighboring cells and the Galerkin basis is twice as large as the basis used for one element. In our proposed method the solution is rebuilt through an L² projection of the discontinuous interface solution on a small rectangular overlapping interface element, named Elastoplast, with an orthogonal basis of the same order as the one in the neighboring cells. Comparisons on 1D and 2D scalar diffusion problems in terms of accuracy and stability with other viscous DG schemes are first given. Then, 2D results on acoustic problems, vortex problems and boundary layer problems both on Cartesian or unstructured triangular grids illustrate stability, precision and versatility of this method.     

Continuous and discontinuous Galerkin methods for a scalable three-dimensional nonhydrostatic atmospheric model: Limited-area mode

This paper describes a unified, element based Galerkin (EBG) framework for a three-dimensional, nonhydrostatic model for the atmosphere. In general, EBG methods possess high-order accuracy, geometric flexibility, excellent dispersion properties and good scalability. Our nonhydrostatic model, based on the compressible Euler equations, is appropriate for both limited-area and global atmospheric simulations. Both a continuous Galerkin (CG), or spectral element, and discontinuous Galerkin (DG) model are considered using hexahedral elements. The formulation is suitable for both global and limited-area atmospheric modeling, although we restrict our attention to 3D limited-area phenomena in this study; global atmospheric simulations will be presented in a follow-up paper. Domain decomposition and communication algorithms used by both our CG and DG models are presented. The communication volume and exchange algorithms for CG and DG are compared and contrasted. Numerical verification of the model was performed using two test cases: flow past a 3D mountain and buoyant convection of a bubble in a neutral atmosphere; these tests indicate that both CG and DG can simulate the necessary physics of dry atmospheric dynamics. Scalability of both methods is shown up to 8192 CPU cores, with near ideal scaling for DG up to 32,768 cores.     

Strong and weak constraint variational assimilations for reduced order fluid flow modeling

In this work we propose and evaluate two variational data assimilation techniques for the estimation of low order surrogate experimental dynamical models for fluid flows. Both methods are built from optimal control recipes and rely on proper orthogonal decomposition and a Galerkin projection of the Navier Stokes equation. The techniques proposed differ in the control variables they involve. The first one introduces a weak dynamical model defined only up to an additional uncertainty time-dependent function whereas the second one, handles a strong dynamical constraint in which the dynamical system’s coefficients constitute the control variables. Both choices correspond to different approximations of the relation between the reduced basis on which is expressed the motion field and the basis components that have been neglected in the reduced order model construction. The techniques have been assessed on numerical data and for real experimental conditions with noisy particle image velocimetry data.     

风力机翼型气动计算的降阶模型研究

将基于POD的降阶模型应用于风力机翼型的气动研究。首先应用CFD数值模拟得到一系列快照结果;应用基于本征正交分解(POD)的方法得到流场的一组基模态,认为对于所研究的问题,任一流场可以由这些基模态通过线性表达得到;对控制方程进行Galerkin投影,得到降阶模型,将离散求解N-S方程的问题转化为一组只有十几个自由度的常微分方程,从而减少计算时间,提高计算效率。并对二维翼型的绕流的定常和非定常问题进行了分析,计算结果表明,降阶模型可以较好地捕捉流动的特征,与直接CFD模拟相比计算精度相当,但大幅有效地提高了计算效率。     

A stable high-order Spectral Difference method for hyperbolic conservation laws on triangular elements

Numerical schemes using piecewise polynomial approximation are very popular for high order discretization of conservation laws. While the most widely used numerical scheme under this paradigm appears to be the Discontinuous Galerkin method, the Spectral Difference scheme has often been found attractive as well, because of its simplicity of formulation and implementation. However, recently it has been shown that the scheme is not linearly stable on triangles. In this paper we present an alternate formulation of the scheme, featuring a new flux interpolation technique using Raviart–Thomas spaces, which proves stable under a similar linear analysis in which the standard scheme failed. We demonstrate viability of the concept by showing linear stability both in the semi-discrete sense and for time stepping schemes of the SSP Runge–Kutta type. Furthermore, we present convergence studies, as well as case studies in compressible flow simulation using the Euler equations.     

Error estimation and anisotropic mesh refinement for 3d laminar aerodynamic flow simulations

This article considers a posteriori error estimation and anisotropic mesh refinement for three-dimensional laminar aerodynamic flow simulations. The optimal order symmetric interior penalty discontinuous Galerkin discretization which has previously been developed for the compressible Navier–Stokes equations in two dimensions is extended to three dimensions. Symmetry boundary conditions are given which allow to discretize and compute symmetric flows on the half model resulting in exactly the same flow solutions as if computed on the full model. Using duality arguments, an error estimation is derived for estimating the discretization error with respect to the aerodynamic force coefficients. Furthermore, residual-based indicators as well as adjoint-based indicators for goal-oriented refinement are derived. These refinement indicators are combined with anisotropy indicators which are particularly suited to the discontinuous Galerkin (DG) discretization. Two different approaches based on either a heuristic criterion or an anisotropic extension of the adjoint-based error estimation are presented. The performance of the proposed discretization, error estimation and adaptive mesh refinement algorithms is demonstrated for 3d aerodynamic flows.     

Shear-thinning-induced chaos in taylor-couette flow

The effect of weak shear thinning on the stability of the Taylor-Couette flow is explored for a Carreau-Bird fluid in the narrow-gap limit. The Galerkin projection method is used to derive a low-order dynamical system from the conservation of mass and momentum equations. In comparison with the Newtonian system, the present equations include additional nonlinear coupling in the velocity components through the viscosity. It is found that the critical Taylor number, corresponding to the loss of stability of the base (Couette) flow, becomes lower as the shear-thinning effect increases. That is, shear thinning tends to precipitate the onset of Taylor vortex flow. Similar to Newtonian fluids, there is an exchange of stability between the Couette and Taylor vortex flows, which coincides with the onset of a supercritical bifurcation. However, unlike the Newtonian model, the Taylor vortex cellular structure loses its stability in turn as the Taylor number reaches a critical value. At this point, a Hopf bifurcation emerges, which exists only for shear-thinning fluids.     

A Hermite WENO reconstruction-based discontinuous Galerkin method for the Euler equations on tetrahedral grids

A Hermite WENO reconstruction-based discontinuous Galerkin method RDG(P₁P₂), designed not only to enhance the accuracy of discontinuous Galerkin method but also to ensure linear stability of the RDG method, is presented for solving the compressible Euler equations on tetrahedral grids. In this RDG(P₁P₂) method, a quadratic polynomial solution (P₂) is first reconstructed using a least-squares method from the underlying linear polynomial (P₁) discontinuous Galerkin solution. By taking advantage of handily available and yet invaluable information, namely the derivatives in the DG formulation, the stencils used in the reconstruction involve only von Neumann neighborhood (adjacent face-neighboring cells) and thus are compact and consistent with the underlying DG method. The final quadratic polynomial solution is then obtained using a WENO reconstruction, which is necessary to ensure linear stability of the RDG method. The developed RDG method is used to compute a variety of flow problems on tetrahedral meshes to demonstrate its accuracy, efficiency, robustness, and versatility. The numerical experiments demonstrate that the developed RDG(P₁P₂) method is able to maintain the linear stability, achieve the designed third-order of accuracy: one order accuracy higher than the underlying DG method without significant increase in computing costs and storage requirements.     

On positivity-preserving high order discontinuous Galerkin schemes for compressible Euler equations on rectangular meshes

We construct uniformly high order accurate discontinuous Galerkin (DG) schemes which preserve positivity of density and pressure for Euler equations of compressible gas dynamics. The same framework also applies to high order accurate finite volume (e.g. essentially non-oscillatory (ENO) or weighted ENO (WENO)) schemes. Motivated by Perthame and Shu (1996) [20] and Zhang and Shu (2010) [26], a general framework, for arbitrary order of accuracy, is established to construct a positivity preserving limiter for the finite volume and DG methods with first order Euler forward time discretization solving one-dimensional compressible Euler equations. The limiter can be proven to maintain high order accuracy and is easy to implement. Strong stability preserving (SSP) high order time discretizations will keep the positivity property. Following the idea in Zhang and Shu (2010) [26], we extend this framework to higher dimensions on rectangular meshes in a straightforward way. Numerical tests for the third order DG method are reported to demonstrate the effectiveness of the methods.     

A study of spectral element and discontinuous Galerkin methods for the Navier–Stokes equations in nonhydrostatic mesoscale atmospheric modeling: Equation sets and test cases

We present spectral element (SE) and discontinuous Galerkin (DG) solutions of the Euler and compressible Navier–Stokes (NS) equations for stratified fluid flow which are of importance in nonhydrostatic mesoscale atmospheric modeling. We study three different forms of the governing equations using seven test cases. Three test cases involve flow over mountains which require the implementation of non-reflecting boundary conditions, while one test requires viscous terms (density current). Including viscous stresses into finite difference, finite element, or spectral element models poses no additional challenges; however, including these terms to either finite volume or discontinuous Galerkin models requires the introduction of additional machinery because these methods were originally designed for first-order operators. We use the local discontinuous Galerkin method to overcome this obstacle. The seven test cases show that all of our models yield good results. The main conclusion is that equation set 1 (non-conservation form) does not perform as well as sets 2 and 3 (conservation forms). For the density current (viscous), the SE and DG models using set 3 (mass and total energy) give less dissipative results than the other equation sets; based on these results we recommend set 3 for the development of future multiscale research codes. In addition, the fact that set 3 conserves both mass and energy up to machine precision motives us to pursue this equation set for the development of future mesoscale models. For the bubble and mountain tests, the DG models performed better. Based on these results and due to its conservation properties we recommend the DG method. In the worst case scenario, the DG models are 50% slower than the non-conservative SE models. In the best case scenario, the DG models are just as efficient as the conservative SE models.     

h-Multigrid for space-time discontinuous Galerkin discretizations of the compressible Navier–Stokes equations

Being implicit in time, the space-time discontinuous Galerkin discretization of the compressible Navier–Stokes equations requires the solution of a non-linear system of algebraic equations at each time-step. The overall performance, therefore, highly depends on the efficiency of the solver. In this article, we solve the system of algebraic equations with a h-multigrid method using explicit Runge–Kutta relaxation. Two-level Fourier analysis of this method for the scalar advection–diffusion equation shows convergence factors between 0.5 and 0.75. This motivates its application to the 3D compressible Navier–Stokes equations where numerical experiments show that the computational effort is significantly reduced, up to a factor 10 w.r.t. single-grid iterations.     

h-Multigrid for space-time discontinuous Galerkin discretizations of the compressible Navier–Stokes equations

Being implicit in time, the space-time discontinuous Galerkin discretization of the compressible Navier–Stokes equations requires the solution of a non-linear system of algebraic equations at each time-step. The overall performance, therefore, highly depends on the efficiency of the solver. In this article, we solve the system of algebraic equations with a h-multigrid method using explicit Runge–Kutta relaxation. Two-level Fourier analysis of this method for the scalar advection–diffusion equation shows convergence factors between 0.5 and 0.75. This motivates its application to the 3D compressible Navier–Stokes equations where numerical experiments show that the computational effort is significantly reduced, up to a factor 10 w.r.t. single-grid iterations.     

A reconstructed discontinuous Galerkin method for the compressible Navier–Stokes equations on arbitrary grids

A reconstruction-based discontinuous Galerkin (RDG) method is presented for the solution of the compressible Navier–Stokes equations on arbitrary grids. The RDG method, originally developed for the compressible Euler equations, is extended to discretize viscous and heat fluxes in the Navier–Stokes equations using a so-called inter-cell reconstruction, where a smooth solution is locally reconstructed using a least-squares method from the underlying discontinuous DG solution. Similar to the recovery-based DG (rDG) methods, this reconstructed DG method eliminates the introduction of ad hoc penalty or coupling terms commonly found in traditional DG methods. Unlike rDG methods, this RDG method does not need to judiciously choose a proper form of a recovered polynomial, thus is simple, flexible, and robust, and can be used on arbitrary grids. The developed RDG method is used to compute a variety of flow problems on arbitrary meshes to demonstrate its accuracy, efficiency, robustness, and versatility. The numerical results indicate that this RDG method is able to deliver the same accuracy as the well-known Bassi–Rebay II scheme, at a half of its computing costs for the discretization of the viscous fluxes in the Navier–Stokes equations, clearly demonstrating its superior performance over the existing DG methods for solving the compressible Navier–Stokes equations.     

The stability of simply supported tubes conveying a compressible fluid

This paper is concerned with establishing the conditions of static stability of a simply supported tube conveying a compressible fluid by application of Euler's method of equilibrium. Timoshenko beam theory is used to describe the tube motion whiler Euler's equations of motion govern the compressible flow through the tube. The eigenvalue problem associated with the linearized equations of motion first derived by Niordson is solved by using Muller's method. The effects on critical velocity of fluid sound speed, tube shear, and tube aspect ratio are parametrically studied. When the flow is subsonic, the aspect ratio increases the critical velocity predicted by the theory while increased aspect ratio decreases the critical velocity when the flow is supersonic. Reduced sound speed and tube shear modulus always reflect a reduced critical velocity for the onset of tube buckling or divergence instability.     

A new high-order discontinuous Galerkin spectral finite element method for Lagrangian gas dynamics in two-dimensions

This paper presents a new high-order cell-centered Lagrangian scheme for two-dimensional compressible flow. The scheme uses a fully Lagrangian form of the gas dynamics equations, which is a weakly hyperbolic system of conservation laws. The system of equations is discretized in the Lagrangian space by discontinuous Galerkin method using a spectral basis. The vertex velocities and the numerical fluxes through the cell interfaces are computed consistently in the Eulerian space by virtue of an improved nodal solver. The nodal solver uses the HLLC approximate Riemann solver to compute the velocities of the vertex. The time marching is implemented by a class of TVD Runge–Kutta type methods. A new HWENO (Hermite WENO) reconstruction algorithm is developed and used as limiters for RKDG methods to maintain compactness of RKDG methods. The scheme is conservative for the mass, momentum and total energy. It can maintain high-order accuracy both in space and time, obey the geometrical conservation law, and achieve at least second order accuracy on quadrilateral meshes. Results of some numerical tests are presented to demonstrate the accuracy and the robustness of the scheme.     

Data-Driven Model Reduction for Stochastic Burgers Equations

We present a class of efficient parametric closure models for 1D stochastic Burgers equations. Casting it as statistical learning of the flow map, we derive the parametric form by representing the unresolved high wavenumber Fourier modes as functionals of the resolved variable’s trajectory. The reduced models are nonlinear autoregression (NAR) time series models, with coefficients estimated from data by least squares. The NAR models can accurately reproduce the energy spectrum, the invariant densities, and the autocorrelations. Taking advantage of the simplicity of the NAR models, we investigate maximal space-time reduction. Reduction in space dimension is unlimited, and NAR models with two Fourier modes can perform well. The NAR model’s stability limits time reduction, with a maximal time step smaller than that of the K-mode Galerkin system. We report a potential criterion for optimal space-time reduction: the NAR models achieve minimal relative error in the energy spectrum at the time step, where the K-mode Galerkin system’s mean Courant–Friedrichs–Lewy (CFL) number agrees with that of the full model.     

A symmetry-preserving discretisation and regularisation model for compressible flow with application to turbulent channel flow

Most simulation methods for compressible flow attain numerical stability at the cost of swamping the fine turbulent flow structures by artificial dissipation. This article demonstrates that numerical stability can also be attained by preserving conservation laws at the discrete level. A new mathematical explanation of conservation in compressible flow reveals that many conservation properties of convection are due to the skew-symmetry of the convection operator. By preserving this skew-symmetry at the discrete level, a fourth-order accurate collocated symmetry-preserving discretisation with excellent conservation properties is obtained. Also a new symmetry-preserving regularisation subgrid-scale model is proposed. The proposed techniques are assessed in simulations of compressible turbulent channel flow. The symmetry-preserving discretisation for compressible flow has good stability without artificial dissipation and yields acceptable results already on coarse grids. Regularisation does not consistently improve upon no-model results, but often compares favourably with eddy-viscosity models.