Toward accurate hybrid prediction techniques for cavity flow noise applications

A large variety of hybrid computational aeroacoustics (CAA) approaches exist differing from each other in the way the source region is modeled, in the way the equations are used to compute the propagation of acoustic waves in a non-quiescent medium, and in the way the coupling between source and acoustic propagation regions is made. This paper makes a comparison between some commonly used numerical methods for aeroacoustic applications. The aerodynamically generated tonal noise by a flow over a 2D rectangular cavity is investigated. Two different cavities are studied. In the first cavity (L/D=4, M=0.5), the sound field is dominated by the cavity wake mode and its higher harmonics, originating from a periodical vortex shedding at the cavity leading edge. In the second cavity (L/D=2, M=0.6), shear-layer modes, due to flow-acoustic interaction phenomena, generate the major components in the noise spectrum. Source domain modeling is carried out using a second-order finite-volume large eddy simulation. Propagation equations, taking into account convection and refraction effects, are solved using high-order finite-difference schemes for the linearized Euler equations and the acoustic perturbation equations. Both schemes are compared with each other for various coupling methods between source region and acoustic region. Conventional acoustic analogies and Kirchhoff methods are rewritten for the various propagation equations and used to obtain near-field acoustic results. The accuracy of the various coupling methods in identifying the noise-generating mechanisms is evaluated. In this way, this paper provides more insight into the practical use of various hybrid CAA techniques to predict the aerodynamically generated sound field by a flow over rectangular cavities. Copyright © 2009 John Wiley & Sons, Ltd.     

A hybrid approach for nonlinear computational aeroacoustics predictions

In many aeroacoustics applications involving nonlinear waves and obstructions in the far-field, approaches based on the classical acoustic analogy theory or the linearised Euler equations are unable to fully characterise the acoustic field. Therefore, computational aeroacoustics hybrid methods that incorporate nonlinear wave propagation have to be constructed. In this study, a hybrid approach coupling Navier–Stokes equations in the acoustic source region with nonlinear Euler equations in the acoustic propagation region is introduced and tested. The full Navier–Stokes equations are solved in the source region to identify the acoustic sources. The flow variables of interest are then transferred from the source region to the acoustic propagation region, where the full nonlinear Euler equations with source terms are solved. The transition between the two regions is made through a buffer zone where the flow variables are penalised via a source term added to the Euler equations. Tests were conducted on simple acoustic and vorticity disturbances, two-dimensional jets (Mach 0.9 and 2), and a three-dimensional jet (Mach 1.5), impinging on a wall. The method is proven to be effective and accurate in predicting sound pressure levels associated with the propagation of linear and nonlinear waves in the near- and far-field regions.     

Strictly stable high order difference approximations for computational aeroacoustics

High order finite difference approximations with improved accuracy and stability properties have been developed for computational aeroacoustics (CAA). One of our new difference operators corresponds to Tam and Webb's DRP scheme in the interior, but is modified near the boundaries to be strictly stable. A unified formulation of the nonlinear and linearized Euler equations is used, which can be extended to the Navier–Stokes equations. The approach has been verified for 1D, 2D and axisymmetric test problems. We have simulated the sound propagation from a rocket launch before lift-off. To cite this article: B. Müller, S. Johansson, C. R. Mecanique 333 (2005).     

A lattice Boltzmann model for two‐dimensional sound wave in the small perturbation compressible flows

A multi‐entropy‐level lattice Boltzmann model for two‐dimensional sound wave equation in the small perturbation flows is presented. In this model, we used higher‐order moment method, multi‐scale technique and the Chapman–Enskog expansion, and multi‐entropy‐level to obtain sound wave equation with isentropic equation. As numerical examples, the Doppler effects in the sound wave propagation, the sound scattering from circular cylinder are simulated. The numerical results show that this model can be used to simulate sound wave propagation. Copyright © 2010 John Wiley & Sons, Ltd.     

Hybrid methods for airframe noise numerical prediction

This paper describes some significant steps made towards the numerical simulation of the noise radiated by the high-lift devices of a plane. Since the full numerical simulation of such configuration is still out of reach for present supercomputers, some hybrid strategies have been developed to reduce the overall cost of such simulations. The proposed strategy relies on the coupling of an unsteady nearfield CFD with an acoustic propagation solver based on the resolution of the Euler equations for midfield propagation in an inhomogeneous field, and the use of an integral solver for farfield acoustic predictions.In the first part of this paper, this CFD/CAA coupling strategy is presented. In particular, the numerical method used in the propagation solver is detailed, and two applications of this coupling method to the numerical prediction of the aerodynamic noise of an airfoil are presented.Then, a hybrid RANS/LES method is proposed in order to perform some unsteady simulations of complex noise sources. This method allows for significant reduction of the cost of such a simulation by considerably reducing the extent of the LES zone. This method is described and some results of the numerical simulation of the three-dimensional unsteady flow in the slat cove of a high-lift profile are presented. While these results remain very difficult to validate with experiments on similar configurations, they represent up to now the first 3D computations of this kind of flow.     

Analysis of the anisotropy of group velocity error due to spatial finite difference schemes from the solution of the 2D linear Euler equations

Numerical differencing schemes are subject to dispersive and dissipative errors, which in one dimension, are functions of a wavenumber. When these schemes are applied in two or three dimensions, the errors become functions of both wavenumber and the direction of the wave. For the Euler equations, the direction of flow and flow velocity are also important. Spectral analysis was used to predict the error in magnitude and direction of the group velocity of vorticity–entropy and acoustic waves in the solution of the linearised Euler equations in a two‐dimensional Cartesian space. The anisotropy in these errors, for three schemes, were studied as a function of the wavenumber, wave direction, mean flow direction and mean flow Mach number. Numerical experiments were run to provide confirmation of the developed theory. Copyright © 2012 John Wiley & Sons, Ltd.     

Finite volume solution to integrated shallow surface–saturated groundwater flow

This paper describes development of an integrated shallow surface and saturated groundwater model (GSHAW5). The surface flow motion is described by the 2‐D shallow water equations and groundwater movement is described by the 2‐D groundwater equations. The numerical solution of these equations is based on the finite volume method where the surface water fluxes are estimated using the Roe shock‐capturing scheme, and the groundwater fluxes are computed by application of Darcy's law. Use of a shock‐capturing scheme ensures ability to simulate steady and unsteady, continuous and discontinuous, subcritical and supercritical surface water flow conditions. Ground and surface water interaction is achieved by the introduction of source‐sink terms into the continuity equations. Two solutions are tightly coupled in a single code. The numerical solutions and coupling algorithms are explained. The model has been applied to 1‐D and 2‐D test scenarios. The results have shown that the model can produce very accurate results and can be used for simulation of situations involving interaction between shallow surface and saturated groundwater flows. Copyright © 2005 John Wiley & Sons, Ltd.     

An implicit meshless method for application in computational fluid dynamics

An implicit meshless scheme is developed for solving the Euler equations, as well as the laminar and Reynolds‐averaged Navier–Stokes equations. Spatial derivatives are approximated using a least squares method on clouds of points. The system of equations is linearised, and solved implicitly using approximate, analytical Jacobian matrices and a preconditioned Krylov subspace iterative method. The details of the spatial discretisation, linear solver and construction of the Jacobian matrix are discussed; and results that demonstrate the performance of the scheme are presented for steady and unsteady two dimensional fluid flows. Copyright © 2012 John Wiley & Sons, Ltd.     

A preconditioning technique for steady Euler solutions

This paper presents a preconditioning technique for solving a two‐dimensional system of hyperbolic equations. The main attractive feature of this approach is that, unlike a technique based on simply extending the solver for a one‐dimensional hyperbolic system, convergence and stability analysis can be investigated. This method represents a genuine numerical algorithm for multi‐dimensional hyperbolic systems. In order to demonstrate the effectiveness of this approach, applications to solving a two‐dimensional system of Euler equations in supersonic flows are reported. It is shown that the Lax–Friedrichs scheme diverges when applied to the original Euler equations. However, convergence is achieved when the same numerical scheme is employed using the same CFL number to solve the equivalent preconditioned Euler system. Copyright © 1999 John Wiley & Sons, Ltd.     

Application of the method of lines to unsteady compressible Euler equations

In the sense of method of lines, numerical solution of the unsteady compressible Euler equations in 1D, 2D and 3D is split into three steps: First, space discretization is performed by the first‐order finite volume method using several approximate Riemann solvers. Second, smoothness and Lipschitz continuity of RHS of the arising system of ordinary dimensional equations (ODEs) is analysed and its solvability is discussed. Finally, the system of ODEs is integrated in time by means of implicit and explicit higher‐order adaptive schemes offered by ODE packages ODEPACK and DDASPK, by a backward Euler scheme based on the linearization of the RHS and by higher‐order explicit Runge–Kutta methods. Time integrators are compared from several points of view, their applicability to various types of problems is discussed, and 1D, 2D and 3D numerical examples are presented. Copyright © 2003 John Wiley & Sons, Ltd.     

Assessment of WENO schemes for multi‐dimensional Euler equations using GPU

This paper presents a detailed study on the implementation of Weighted Essentially Non‐Oscillatory (WENO) schemes on GPU. GPU implementation of up to ninth‐order accurate WENO schemes for the multi‐dimensional Euler equations of gas dynamics is presented. The implementation detail is discussed in the paper. The computational times of different schemes are obtained and the speedups are reported for different number of grid points. Furthermore, the execution times for the main kernels of the code are given and compared with each other. The numerical experiments show the speedups for the WENO schemes are very promising especially for fine grids. Copyright © 2014 John Wiley & Sons, Ltd.     

The non‐reflective interface: an innovative forcing technique for computational acoustic hybrid methods

The present article concerns a commonly used methodology for the numerical simulation of acoustic emission and propagation phenomena. We consider the so‐called multi‐stage hybrid acoustic approach, in which a given noise problem is simulated via a sequence of weakly coupled computations of noise generation and acoustic propagation stages, wherein the simulation of the propagation stage is based on advanced Computational AeroAcoustics (CAA) techniques. The paper introduces an original forcing technique, namely, the Non‐Reflective Interface (NRI), to enable the transfer of an acoustic signal from an a priori noise generation stage into a CAA‐based acoustic propagation phase. Unlike most existing forcing techniques, the NRI is non‐reflective (or anechoic) in nature and, therefore, can properly handle the backscattering effects arising during the noise propagation stage. This attribute makes the NRI‐based weak‐coupling procedure and the associated CAA‐based hybrid approach compatible with a larger variety of realistic noise problems (such as those involving installed configurations in wind tunnel experiments, for instance). The NRI technique is first validated via several test cases of increasing complexity and is then applied to two aerodynamic noise problems. Copyright © 2015 John Wiley & Sons, Ltd. Copyright © 2015 John Wiley & Sons, Ltd.     

Verification of Euler/Navier–Stokes codes using the method of manufactured solutions

The method of manufactured solutions is used to verify the order of accuracy of two finite‐volume Euler and Navier–Stokes codes. The Premo code employs a node‐centred approach using unstructured meshes, while the Wind code employs a similar scheme on structured meshes. Both codes use Roe's upwind method with MUSCL extrapolation for the convective terms and central differences for the diffusion terms, thus yielding a numerical scheme that is formally second‐order accurate. The method of manufactured solutions is employed to generate exact solutions to the governing Euler and Navier–Stokes equations in two dimensions along with additional source terms. These exact solutions are then used to accurately evaluate the discretization error in the numerical solutions. Through global discretization error analyses, the spatial order of accuracy is observed to be second order for both codes, thus giving a high degree of confidence that the two codes are free from coding mistakes in the options exercised. Examples of coding mistakes discovered using the method are also given. Copyright © 2004 John Wiley & Sons, Ltd.     

Effects of the Jacobian evaluation on Newton's solution of the Euler equations

Newton's method is developed for solving the 2‐D Euler equations. The Euler equations are discretized using a finite‐volume method with upwind flux splitting schemes. Both analytical and numerical methods are used for Jacobian calculations. Although the numerical method has the advantage of keeping the Jacobian consistent with the numerical residual vector and avoiding extremely complex analytical differentiations, it may have accuracy problems and need longer execution time. In order to improve the accuracy of numerical Jacobians, detailed error analyses are performed. Results show that the finite‐difference perturbation magnitude and computer precision are the most important parameters that affect the accuracy of numerical Jacobians. A method is developed for calculating an optimal perturbation magnitude that can minimize the error in numerical Jacobians. The accuracy of the numerical Jacobians is improved significantly by using the optimal perturbation magnitude. The effects of the accuracy of numerical Jacobians on the convergence of the flow solver are also investigated. In order to reduce the execution time for numerical Jacobian evaluation, flux vectors with perturbed flow variables are calculated only for neighbouring cells. A sparse matrix solver that is based on LU factorization is used. Effects of different flux splitting methods and higher‐order discretizations on the performance of the solver are analysed. Copyright © 2005 John Wiley & Sons, Ltd.     

Riemann solvers with evolved initial conditions

The scope of this paper is three fold. We first formulate upwind and symmetric schemes for hyperbolic equations with non‐conservative terms. Then we propose upwind numerical schemes for conservative and non‐conservative systems, based on a Riemann solver, the initial conditions of which are evolved non‐linearly in time, prior to a simple linearization that leads to closed‐form solutions. The Riemann solver is easily applied to complicated hyperbolic systems. Finally, as an example, we formulate conservative schemes for the three‐dimensional Euler equations for general compressible materials and give numerical results for a variety of test problems for ideal gases in one and two space dimensions. Copyright © 2006 John Wiley & Sons, Ltd.     

An efficient discontinuous Galerkin method for aeroacoustic propagation

An efficient discontinuous Galerkin formulation is applied to the solution of the linearized Euler equations and the acoustic perturbation equations for the simulation of aeroacoustic propagation in two‐dimensional and axisymmetric problems, with triangular and quadrilateral elements. To improve computational efficiency, a new strategy of variable interpolation order is proposed in addition to a quadrature‐free approach and parallel implementation. Moreover, an accurate wall boundary condition is formulated on the basis of the solution of the Riemann problem for a reflective wall. Time discretization is based on a low dissipation formulation of a fourth‐order, low storage Runge–Kutta scheme. Along the far‐field boundaries a perfectly matched layer boundary condition is used. For the far‐field computations, the integral formulation of Ffowcs Williams and Hawkings is coupled with the near‐field solver. The efficiency and accuracy of the proposed variable order formulation is assessed for realistic geometries, namely sound propagation around a high‐lift airfoil and the Munt problem. Copyright © 2011 John Wiley & Sons, Ltd.     

A Three-Dimensional Hybrid LES-Acoustic Analogy Method for Predicting Open-Cavity Noise

A three-dimensional (3D) hybrid LES-acoustic analogy method for computational aeroacoustics (CAA) is presented for the prediction of open-cavity noise. The method uses large-eddy simulation (LES) to compute the acoustic source while the Ffowcs Williams-Hawkings (FW-H) acoustic analogy is employed for the prediction of the far-field sound. As a comparison, a two-dimensional (2D) FW-H analogy is also included. The hybrid method has been assessed in an open-cavity flow at a Mach number of 0.85 and a Reynolds number of Re=1.36×10⁶, where some experimental data are available for comparison. The study has identified some important technical issues in the application of the FW-H acoustic analogy to cavity noise prediction and CAA in general, including the proper selection of the integration period and the modes of sound sources in the frequency domain. The different nature of 2D and 3D wave propagation is also highlighted, which calls for a matching acoustic solver for each problem. The developed hybrid method has shown promise to be a feasible, accurate and computationally affordable approach for CAA.     

Sound wave propagation through incompressible flows

This paper derives a variable coefficient, multidimensional Burgers equation which models the propagation of a nonlinear sound wave through an incompressible background flow. The equation is derived from the compressible Euler equations by the combination of a weakly nonlinear acoustics expansion for the sound wave and an incompressible expansion for the background flow. The main effect of the incompressible flow on the sound wave is the advection of the sound wave by the transverse velocity component of the flow.     

Numerical simulation of vortical ideal fluid flow through curved channel

A numerical algorithm to study the boundary‐value problem in which the governing equations are the steady Euler equations and the vorticity is given on the inflow parts of the domain boundary is developed. The Euler equations are implemented in terms of the stream function and vorticity. An irregular physical domain is transformed into a rectangle in the computational domain and the Euler equations are rewritten with respect to a curvilinear co‐ordinate system. The convergence of the finite‐difference equations to the exact solution is shown experimentally for the test problems by comparing the computational results with the exact solutions on the sequence of grids. To find the pressure from the known vorticity and stream function, the Euler equations are utilized in the Gromeka–Lamb form. The numerical algorithm is illustrated with several examples of steady flow through a two‐dimensional channel with curved walls. The analysis of calculations shows strong dependence of the pressure field on the vorticity given at the inflow parts of the boundary. Plots of the flow structure and isobars, for different geometries of channel and for different values of vorticity on entrance, are also presented. Copyright © 2003 John Wiley & Sons, Ltd.