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.     

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.     

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.     

COHESIVE ZONE FINITE ELEMENT-BASED MODELING OF HYDRAULIC FRACTURES

Hydraulic fracturing is a powerful technology used to stimulate fluid production from reservoirs. The fully 3-D numerical simulation of the hydraulic fracturing process is of great importance to the efficient application of this technology, but is also a great challenge because of the strong nonlinear coupling between the viscous flow of fluid and fracture propagation. By taking advantage of a cohesive zone method to simulate the fracture process, a finite element model based on the existing pore pressure cohesive finite elements has been established to investigate the propagation of a penny-shaped hydraulic fracture in an infinite elastic medium. The effect of cohesive material parameters and fluid viscosity on the hydraulic fracture behaviour has been investigated. Excellent agreement between the finite element results and analytical solutions for the limiting case where the fracture process is dominated by rock fracture toughness demonstrates the ability of the cohesive zone finite element model in simulating the hydraulic fracture growth for this case.     

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.     

Generation of acoustic solitary waves in a lattice of Helmholtz resonators

This paper addresses the propagation of high amplitude acoustic pulses through a 1D lattice of Helmholtz resonators connected to a waveguide. Based on the model proposed by Sugimoto (1992), a new numerical method is developed to take into account both the nonlinear wave propagation and the different mechanisms of dissipation: the volume attenuation, the linear viscothermal losses at the walls, and the nonlinear absorption due to the acoustic jet formation in the resonator necks. Good agreement between numerical and experimental results is obtained, highlighting the crucial role of the nonlinear losses. Different kinds of solitary waves are observed experimentally with characteristics depending on the dispersion properties of the lattice.     

Computational aeroacoustics applications based on a discontinuous Galerkin method

caa simulation requires the calculation of the propagation of acoustic waves with low numerical dissipation and dispersion error, and to take into account complex geometries. To give, at the same time, an answer to both challenges, a Discontinuous Galerkin Method is developed for Computational AeroAcoustics. Euler's linearized equations are solved with the Discontinuous Galerkin Method using flux splitting technics. Boundary conditions are established for rigid wall, non-reflective boundary and imposed values. A first validation, for induct propagation is realized. Then, applications illustrate: the Chu and Kovasznay's decomposition of perturbation inside uniform flow in term of independent acoustic and rotational modes, Kelvin–Helmholtz instability and acoustic diffraction by an air wing. To cite this article: Ph. Delorme et al., C. R. Mecanique 333 (2005).     

The effect of virtual mass on the basic equations for unsteady one-dimensional heterogeneous flows

A separated two-component flow model is presented which includes virtual mass forces coupling the momentum equations of the two components. It is shown that for physically realistic situations four real roots of the characteristics determinant can exist. These are associated with the acoustic propagation velocities and the fiow velocities of the constitutive phases. Direct analytical solution of the full characteristic determinant is difficult. However, for low Mach number flows an acoustic propagation velocity is obtained which falls between the well-known true separated and homogeneous wave speeds, and compares favorably with experimental data for glass/water and air/water mixtures.     

Linear and nonlinear approaches and digital signal processing in aeroacoustics

The linear and nonlinear approaches to the calculation of small acoustic disturbance propagation and evolution in nonuniform flows are compared. In the conventional linear approach it is the linearized equations of time-dependent, ideal (inviscid and non-heat-conducting) or viscous gas flow that are integrated. In the nonlinear approach the original nonlinear equations governing the same time-dependent flow (Euler equations for an ideal gas) are integrated; these are the same equations that, together with time relaxation procedure, are used in the linear approach for calculating the stationary background. It is shown that the application of digital signal processing, widely used in acoustic experiments, makes it possible to isolate the harmonic acoustic waves from the results of integration of the nonlinear equations, though their intensity is smaller than that of the noise due to computational errors, including inadequate attainment of the stationary background.     

Nonlinear disturbance evolution in a two-dimensional boundary-layer along an elastic plate and induced radiated sound

The interaction between a boundary-layer flow and an elastic plate is addressed by direct numerical simulation, taking into account the full coupling between the fluid flow and the flexible wall. The convectively unstable flow state is harmonically forced and two-dimensional nonlinearly saturated wavelike disturbances are computed along archetype-plates with respect to stiffness and natural frequencies. In the aim of determining the low-Mach-number radiated sound for the system, the simulation data are used to solve the Lighthill's equation in terms of a Green's function in the wavenumber-frequency space. Different degrees of fluid-structure coupling are implemented in the radiated sound model and the resulting acoustic pressure levels are compared. The sound radiation levels are shown to be increased in the presence of flexible walls with however significant differences in the radiated pressure levels for different coupling assumptions.     

Linearized and non‐linear acoustic/viscous splitting techniques for low Mach number flows

Computation of the acoustic disturbances generated by unsteady low‐speed flow fields including vortices and shear layers is considered. The equations governing the generation and propagation of acoustic fluctuations are derived from a two‐step acoustic/viscous splitting technique. An optimized high order dispersion–relation–preserving scheme is used for the solution of the acoustic field. The acoustic field generated by a corotating vortex pair is obtained using the above technique. The computed sound field is compared with the existing analytic solution. Results are in good agreement with the analytic solution except near the centre of the vortices where the acoustic pressure becomes singular. The governing equations for acoustic fluctuations are then linearized and solved for the same model problem. The difference between non‐linear and linearized solutions falls below the numerical error of the simulation. However, a considerable saving in CPU time usage is achieved in solving the linearized equations. The results indicate that the linearized acoustic/viscous splitting technique for the simulation of acoustic fluctuations generation and propagation by low Mach number flow fields seems to be very promising for three‐dimensional problems involving complex geometries. Copyright © 2003 John Wiley & Sons, Ltd.     

Effect of a Standing Acoustic Wave on the Development of Large-Scale Convection in a Horizontal Layer

The effect of a standing acoustic wave on the development of long-wave convective perturbations in a horizontal layer with thermally insulated boundaries is investigated. The main two-dimensional flow is determined. A nonlinear amplitude equation with spatially-periodic coefficients is derived for investigating the stability of the main flow and secondary convection flows in the neighborhood of the stability threshold. The intensity of the acoustic field is assumed to be low. It is shown that the acoustic action leads to destabilization of the layer. Plane and three-dimensional perturbations are critical at large and small Prandtl numbers, respectively. Nonlinear one-dimensional steady-state solutions of the amplitude equation are obtained and their stability is investigated.     

水下声波传播过程的时域间断Galerkin有限元方法

本文构建了声压波动方程的改进时域间断Galerkin有限元方法.传统时域连续有限元方法在计算高梯度、强间断特征水中声波传播问题时往往会出现虚假数值振荡现象,这些数值振荡会影响正常波动的计算精度.为了解决这一问题,本文通过引入人工阻尼的方式构建了改进的时域间断Galerkin有限元方法,并针对具有高梯度、强间断特征的多障...     

Finite-amplitude acoustic pulses in a waveguide layer with longitudinal inhomogeneity

This paper considers the propagation of a weakly nonlinear acoustic pulse in a slightly curved waveguide layer which is strongly inhomogeneous in the transverse direction and weakly inhomogeneous in the longitudinal direction. The basic system of hydrodynamic equations reduces to a nonlinear wave equation, whose coefficients are determined using the equations of state of the medium. It is established that as the adiabatic exponent passes through the value γ = 3/2, the nature of the pulse propagation changes: for large values of γ, the medium is focusing, and for smaller values, it is defocusing. It is shown that the pulse propagation process is characterized by three scales: the high-frequency filling is modulated by the envelope, whose evolution, in turn, is determined by the moderate-rate evolution of the envelope phase and slow amplitude variation. A generalized nonlinear Schrödinger equation with the coefficients dependent on the longitudinal coordinate is derived for the pulse envelope. An explicit soliton solution of this equation is constructed for some types of longitudinal inhomogeneity.     

Review of Computational Aeroacoustics Algorithms

This paper presents a review of recent advancements in computational methodology for aeroacoustics problems. High-order finite difference methods for computation of linear and nonlinear acoustic waves are the primary focus of the review. Schemes for numerical simulation of linear waves include explicit optimized and DRP finite-difference operators, compact schemes, wavenumber extended upwind schemes and leapfrog-like algorithms. Both spatial approximations and time-integration techniques, which include low-dissipation low-dispersion Adams-Bashforth and Runge-Kutta (RK) methods, are examined. Wave propagation properties are analysed in the wavenumber and frequency space. Different approaches to eliminate short-wave spurious numerical waves are also reviewed. Methods for simulating nonlinear acoustic phenomena include essentially non-oscillatory (ENO) schemes, numerical adaptive filtering for high-order explicit and compact finite-difference operators, MacCormack and adaptive compact nonlinear algorithms. A literature survey of other CAA methods is provided in the introductory part.     

Numerical simulation of flow‐induced noise using LES/SAS and Lighthill's acoustic analogy

We present a finite element (FE) formulation of Lighthill's acoustic analogy for the hybrid computation of noise generated by turbulent flows. In the present approach, the flow field is computed using large eddy simulation and scale adaptive simulation turbulence models. The acoustic propagation is obtained by solving the variational formulation of Lighthill's acoustic analogy with the FE method. In order to preserve the acoustic energy, we compute the inhomogeneous part of Lighthill's wave equation by applying the FE formulation on the fine flow grid. The resulting acoustic nodal loads are then conservatively interpolated to the coarser acoustic grid. Subsequently, the radiated acoustic field can be solved in both time and frequency domains. In the latter case, an enhanced perfectly matched layer technique is employed, allowing one to truncate the computational domain in the acoustic near field, without compromising the numerical solution. Our hybrid approach is validated by comparing the numerical results of the acoustic field induced by a corotating vortex pair with the corresponding analytical solution. To demonstrate the applicability of our scheme, we present full 3D numerical results for the computed acoustic field generated by the turbulent flow around square cylinder geometries. The sound pressure levels obtained compare well with measured values. Copyright © 2009 John Wiley & Sons, Ltd.     

Nonreflecting boundary conditions based on nonlinear multidimensional characteristics

Nonlinear characteristic boundary conditions based on nonlinear multidimensional characteristics are proposed for 2‐ and 3‐D compressible Navier–Stokes equations with/without scalar transport equations. This approach is consistent with the flow physics and transport properties. Based on the theory of characteristics, which is a rigorous mathematical technique, multidimensional flows can be decomposed into acoustic, entropy, and vorticity waves. Nonreflecting boundary conditions are derived by setting corresponding characteristic variables of incoming waves to zero and by partially damping the source terms of the incoming acoustic waves. In order to obtain the resulting optimal damping coefficient, analysis is performed for problems of pure acoustic plane wave propagation and arbitrary flows. The proposed boundary conditions are tested on two benchmark problems: cylindrical acoustic wave propagation and the wake flow behind a cylinder with strong periodic vortex convected out of the computational domain. This new approach substantially minimizes the spurious wave reflections of pressure, density, temperature, and velocity as well as vorticity from the artificial boundaries, where strong multidimensional flow effects exist. The numerical simulations yield accurate results, confirm the optimal damping coefficient obtained from analysis, and verify that the method substantially improves the 1‐D characteristics‐based nonreflecting boundary conditions for complex multidimensional flows. Copyright © 2009 John Wiley & Sons, Ltd.     

Aeroacoustical coupling in a ducted shallow cavity and fluid/structure effects on a steam line

A pure tone phenomenon has been observed at 460 Hz in a piping steam line. The acoustical energy has been identified to be generated in an open gate valve and to be of cavity noise type. This energy is then transmitted to the main pipe by fluid/structure coupling. The objectives here are to display the mechanism of the flow acoustic coupling in the cavity and in the duct through an aeroacoustical analysis and to understand the way of energy transfer from the fluid to the main pipe through a vibroacoustical analysis. Concerning the first objective, an experimental study by means of 2/7 scale models in air is analysed by means of numerical flow simulation. The flow acoustic phenomena are modelled by computing the Euler equations. Two different computations are carried out: in the first one, a pure Euler modelling is used, in the second one, a boundary layer obtained from experimental data is introduced in the computation in order to have a realistic flow profile upstream the cavity. The boundary layer flow profile appears to be essential to recover the experimentally observed coupling between the shear-layer instability and the acoustical transverse mode of the pipe. The numerical results confirm that the second aerodynamic mode is responsible for the oscillation. While the predicted frequency agrees about 1% with the scale model experiments, the predicted amplitude is approximately 15 dB too low. For the second objective, fluid/structure coupling in the main pipe is studied using two fully coupled methods. The first method consists in a modal analysis of the line using a fluid–structure finite element model. The second one is based on the analysis of dispersion diagrams derived from the local equations of cylindrical shells filled with fluid. The way of energy transfer in transverse acoustical waves coupled with flexion-ovalization deformations of the pipe is highlighted using both methods. The dispersion diagrams allow a fast and accurate analysis. The modal analysis using a finite-element model may complete the first one with quantitative data. The link between the fluid/acoustic and the fluid/structure analysis is then the excitation of the transverse acoustical mode of the duct.     

Energy consideration of the nonlinear effects in a Rijke tube

The goal of this work is to characterize the excited states of a thermoacoustic system with mean flow. The properties of excited regimes are determined by the balance between thermoacoustic energy transformation and acoustic losses. In many systems, the sound intensity is not sufficient for nonlinear acoustic losses to be a major factor in defining nonlinear saturation of thermoacoustic instability. It is the nonlinearity of the heat transfer process that is responsible for limit-cycle stabilization of linearly unstable acoustic modes and for the appearance of higher harmonics. In the present study, both a nonlinear theory based on energy consideration and a model for the nonlinear convective heat transfer in unsteady flow are developed. Experimental data are obtained for the excited regimes of operation of an electric Rijke tube. Model results for hysteresis in the transition between stable and excited states and for limit-cycle parameters are compared with test data.     

Acoustic diode: Wave non-reciprocity in nonlinearly coupled waveguides

The paper describes a passive time-independent setting for non-reciprocal wave transmission in mechanical and acoustic systems with strong nonlinearities. In the proposed system vibro-impact elements with pre-defined clearances are used to couple two non-dispersive waveguides. The asymmetry necessary for the non-reciprocal behavior is realized through unequal grounding springs of the vibro-impact elements. We show that, for appropriate selection of the parameters, the proposed system acts as a mechanical diode, allowing the transmission of acoustic waves in one direction and completely preventing reverse transmission. Two different designs of the coupling elements are suggested, with the possibility of single-sided or double-sided impacts. A unique feature of the proposed non-reciprocal acoustic system is that minimal distortion of the harmonic content of the transmitted wave occurs, in contrast to current designs where nonlinear non-reciprocity is achieved at the expense of a rather strong distortion of the transmitted signals. For both designs, we derive exact solutions for propagation and reflection of the harmonic waves, and demonstrate the possibility for strong non-reciprocity. Stability properties of the observed solutions in the space of parameters are also explored.