Residual‐based stabilization of the finite element approximation to the acoustic perturbation equations for low Mach number aeroacoustics

The acoustic perturbation equations (APE) are suitable to predict aerodynamic noise in the presence of a non‐uniform mean flow. As for any hybrid computational aeroacoustics approach, a first computational fluid dynamics simulation is carried out from which the mean flow characteristics and acoustic sources are obtained. In a second step, the APE are solved to get the acoustic pressure and particle velocity fields. However, resorting to the finite element method (FEM) for that purpose is not straightforward. Whereas mixed finite elements satisfying an appropriate inf–sup compatibility condition can be built in the case of no mean flow, that is, for the standard wave equation in mixed form, these are difficult to implement and their good performance is yet to be checked for more complex wave operators. As a consequence, strong simplifying assumptions are usually considered when solving the APE with FEM. It is possible to avoid them by resorting to stabilized formulations. In this work, a residual‐based stabilized FEM is presented for the APE at low Mach numbers, which allows one to deal with the APE convective and reaction terms in its full extent. The key of the approach resides in the design of the matrix of stabilization parameters. The performance of the formulation and the contributions of the different terms in the equations are tested for an acoustic pulse propagating in sheared‐solenoidal mean flow, and for the aeolian tone generated by flow past a two‐dimensional cylinder. Copyright © 2016 John Wiley & Sons, Ltd.     

On the signal degradation induced by the interpolation and the sampling rate reduction in aeroacoustics hybrid methods

Taking place within more extensive work that focuses on hybrid methods in aeroacoustics, the present study is devoted to the data transfer operations that are to be performed between two stages of a hybrid calculation. More precisely, the article focuses on two typical operations that usually accompany such data transfer, which are (i) the sampling rate reduction and (ii) the interpolation of the unsteady perturbations to be transmitted from one stage to the other. First part of the paper analyzes the two main issues of such operations, which are the spuriousing and the aliasing phenomena. For doing so, the usual notions of the interpolation theory are revisited before they are synthesized within an original approach. The here proposed formalism allows to understand better both the spuriousing and the aliasing phenomena, as well as to accurately predict the impact of the latter on the data to be transmitted in terms of signal degradation. Second part of the paper provides an illustration and a validation of these theoretical developments via a direct application to a typical aerodynamic noise problem (aeroacoustic emission by a 2D cylinder cross flow). There, it is further shown how the here proposed formalism can help in improving aeroacoustics hybrid calculations by predicting (and thus possibly minimizing) the bias to be induced on the acoustic extrapolation stage because of the aliasing and/or spuriousing effects inherited from the sampling rate reduction and/or interpolation of CFD data—which is likely to occur in any hybrid scenario. Copyright © 2012 John Wiley & Sons, Ltd.     

Toward low‐noise synthetic turbulent inflow conditions for aeroacoustic calculations

The accuracy of boundary conditions for computational aeroacoustics is a well‐known challenge, due in part to the necessity of truncating the flow domain and replacing the analytical boundary conditions at infinity with numerical boundary conditions. In particular, the inflow boundary condition involving turbulent velocity or scalar fields is likely to introduce spurious waves into the domain, therefore degrading the flow behavior and deteriorating the physical acoustic waves. In this work, a method to generate low‐noise, divergence‐free, synthetic turbulence for inflow boundary conditions is proposed. It relies on the classical view of turbulence as a superposition of random eddies convected with the mean flow. Within the proposed model, the vector potential and the requirement that the individual eddies must satisfy the linearized momentum equations about the mean flow are used. The model is tested using isolated eddies convected through the inflow boundary and an experimental benchmark data for spatially decaying isotropic turbulence. Copyright © 2013 John Wiley & Sons, Ltd.     

Resolution of computational aeroacoustics problems on unstructured grids with a higher-order finite volume scheme

Computational fluid dynamics (CFD) has become increasingly used in the industry for the simulation of flows. Nevertheless, the complex configurations of real engineering problems make the application of very accurate methods that only work on structured grids difficult. From this point of view, the development of higher-order methods for unstructured grids is desirable. The finite volume method can be used with unstructured grids, but unfortunately it is difficult to achieve an order of accuracy higher than two, and the common approach is a simple extension of the one-dimensional case. The increase of the order of accuracy in finite volume methods on general unstructured grids has been limited due to the difficulty in the evaluation of field derivatives. This problem is overcome with the application of the Moving Least Squares (MLS) technique on a finite volume framework. In this work we present the application of this method (FV-MLS) to the solution of aeroacoustic problems.     

On the discontinuous Galerkin computation of N‐waves focusing

Sonic boom focusing phenomenon can be predicted using the solution to the nonlinear Tricomi equation which is a hybrid (hyperbolic‐elliptic) second‐order partial differential equation. In this paper, the hyperbolic conservation law form is derived, which is valid in the entire domain. In this manner, the presence of two regions where the equation behaves differently (hyperbolic in the upper and elliptic in the lower half‐plane) is avoided. On the upper boundary, a new mixed boundary condition for the acoustic pressure is employed. The discretization is carried out using a discontinuous Galerkin (DG) method combined with a Runge–Kutta total‐variation diminishing scheme. The results show the accuracy of DG methods to solve problems involving sharp gradients and discontinuities. Comparisons with analytical results for the linear case, and other numerical results using classical explicit and compact finite difference schemes and weighted essentially non‐oscillatory schemes are included. Copyright © 2010 John Wiley & Sons, Ltd.     

Aeroacoustic Simulations Using Compressible Lattice Boltzmann Method

This paper presents a lattice Boltzmann (LB) method based study aimed at numerical simulation of aeroacoustic phenomenon in flows around a symmetric obstacle. To simulate the compressible flow accurately, a potential energy double-distribution-function (DDF) lattice Boltzmann method is used over the entire computational domain from the near to far fields. The buffer zone and absorbing boundary condition is employed to eliminate the non-physical reflecting. Through the direct numerical simulation, the flow around a circular cylinder at $Re$=150, $M$=0.2 and the flow around a NACA0012 airfoil at $Re$=10000, $M$=0.8, $α$=$0^◦$ are investigated. The generation and propagation of the sound produced by the vortex shedding are reappeared clearly. The obtained results increase our understanding of the characteristic features of the aeroacoustic sound.     

Mathematical model of micropolar fluid in two-phase immiscible fluid flow through porous channel

This paper is concerned with the flow of two immiscible fluids through a porous horizontal channel. The fluid in the upper region is the micropolar fluid/the Eringen fluid, and the fluid in the lower region is the Newtonian viscous fluid. The flow is driven by a constant pressure gradient. The presence of micropolar fluids introduces additional rotational parameters. Also, the porous material considered in both regions has two different permeabilities. A direct method is used to obtain the analytical solution of the concerned problem. In the present problem, the effects of the couple stress, the micropolarity parameter, the viscosity ratio, and the permeability on the velocity profile and the microrotational velocity are discussed. It is found that all the physical parameters play an important role in controlling the translational velocity profile and the microrotational velocity. In addition, numerical values of the different flow parameters are computed. The effects of the different flow parameters on the flow rate and the wall shear stress are also discussed graphically.     

Contributions of Computational Aeroacoustics to Jet Noise Research and Prediction

A brief review of recent progress in the field of computational aeroacoustics (CAA) is proposed. This paper is complementary to the previous reviews of Tam [(1995a) “Computational aeroacoustics: issues and methods”, AIAA J. 33(10), 1788–1796], Lele [(1997) “Computational Aeroacoustics: a review”, AIAA Paper 97–0018, 35th Aerospace Sciences Meeting and Exhibit, Reno, Nevada] and Glegg [(1999) “Recent advances aeroacoustics: the influence of computational fluid dynamics”, 6th International Congress on Sound and Vibration, Copenhagen, Danemark, 5–8 July, 43–58] on advances in CAA. After a short introduction concerning the current motivations of jet noise studies, connections between computational fluid dynamics (CFD) and CAA using hybrid approaches are discussed in the first part. The most spectacular advances are probably provided by the direct computation of jet noise, and some recent results are shown in the second part.     

Radiation and Wall Boundary Conditions for Computational Aeroacoustics: A Review

A review of unsteady computational boundary conditions for computational aeroacoustics (CAA) problems is presented. This review is meant to serve as a general overview of previous work on solid wall, radiation and outflow boundary conditions that have been proposed and used in CAA calculations. Both the physical nature of the boundary condition problem as well as the numerical considerations affecting their implementation are discussed.     

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.