Search NASA⌕ Search

SEARCH · Search NASA

Results for “Computational Aeroacoustics”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 145 records · Page 8

CFD methods development considerations for unsteady aerodynamic analysis

The development of computational fluid dynamics (CFD) methods for unsteady aerodynamic analysis is described. Special emphasis is placed on considerations that are required for application of the methods to unsteady aerodynamic flow problems. Two broad categories of topics are presented to illustrate the major points. Although primary application of these CFD methods is to relatively low frequency oscillatory phenomena such as flutter, the ideas that are presented may be of value to developers of computational aeroacoustic methods for predicting high frequency acoustics.

Batina, John T.↗

Numerical predictions in acoustics

Computational Aeroacoustics (CAA) involves the calculation of the sound produced by a flow as well as the underlying flowfield itself from first principles. This paper describes the numerical challenges of CAA and recent research efforts to overcome these challenges. In addition, it includes the benefits of CAA in removing restrictions of linearity, single frequency, constant parameters, low Mach numbers, etc. found in standard acoustic analyses as well as means for evaluating the validity of these numerical approaches. Finally, numerous applications of CAA to both classical as well as modern problems of concern to the aerospace industry are presented.

Hardin, Jay C.↗

The method of space-time conservation element and solution element-applications to one-dimensional and two-dimensional time-marching flow problems

A nontraditional numerical method for solving conservation laws is being developed. The new method is designed from a physicist's perspective, i.e., its development is based more on physics than numerics. Even though it uses only the simplest approximation techniques, a 2D time-marching Euler solver developed recently using the new method is capable of generating nearly perfect solutions for a 2D shock reflection problem used by Helen Yee and others. Moreover, a recent application of this solver to computational aeroacoustics (CAA) problems reveals that: (1) accuracy of its results is comparable to that of a 6th order compact difference scheme even though nominally the current solver is only of 2nd-order accuracy; (2) generally, the non-reflecting boundary condition can be implemented in a simple way without involving characteristic variables; and (3) most importantly, the current solver is capable of handling both continuous and discontinuous flows very well and thus provides a unique numerical tool for solving those flow problems where the interactions between sound waves and shocks are important, such as the noise field around a supersonic over- or under-expansion jet.

Chang, Sin-Chung↗

An extension of the Lighthill theory of jet noise to encompass refraction and shielding

A formalism for jet noise prediction is derived that includes the refractive 'cone of silence' and other effects; outside the cone it approximates the simple Lighthill format. A key step is deferral of the simplifying assumption of uniform density in the dominant 'source' term. The result is conversion to a convected wave equation retaining the basic Lighthill source term. The main effect is to amend the Lighthill solution to allow for refraction by mean flow gradients, achieved via a frequency-dependent directional factor. A general formula for power spectral density emitted from unit volume is developed as the Lighthill-based value multiplied by a squared 'normalized' Green's function (the directional factor), referred to a stationary point source. The convective motion of the sources, with its powerful amplifying effect, also directional, is already accounted for in the Lighthill format: wave convection and source convection are decoupled. The normalized Green's function appears to be near unity outside the refraction dominated 'cone of silence', this validates our long term practice of using Lighthill-based approaches outside the cone, with extension inside via the Green's function. The function is obtained either experimentally (injected 'point' source) or numerically (computational aeroacoustics). Approximation by unity seems adequate except near the cone and except when there are shrouding jets: in that case the difference from unity quantifies the shielding effect. Further extension yields dipole and monopole source terms (cf. Morfey, Mani, and others) when the mean flow possesses density gradients (e.g., hot jets).

Ribner, Herbert S.↗

Benchmark problems and solutions

The scientific committee, after careful consideration, adopted six categories of benchmark problems for the workshop. These problems do not cover all the important computational issues relevant to Computational Aeroacoustics (CAA). The deciding factor to limit the number of categories to six was the amount of effort needed to solve these problems. For reference purpose, the benchmark problems are provided here. They are followed by the exact or approximate analytical solutions. At present, an exact solution for the Category 6 problem is not available.

Tam, Christopher K. W.↗

Accurate Finite Difference Algorithms

Two families of finite difference algorithms for computational aeroacoustics are presented and compared. All of the algorithms are single step explicit methods, they have the same order of accuracy in both space and time, with examples up to eleventh order, and they have multidimensional extensions. One of the algorithm families has spectral like high resolution. Propagation with high order and high resolution algorithms can produce accurate results after O(10(exp 6)) periods of propagation with eight grid points per wavelength.

Goodrich, John W.↗

Application of the Discontinuous Galerkin Method to Acoustic Scatter Problems

The application of the quadrature-free form of the discontinuous Galerkin method to two problems from Category 1 of the Second Computational Aeroacoustics Workshop on Benchmark problems is presented. The method and boundary conditions relevant to this work are described followed by two test problems, both of which involve the scattering of an acoustic wave off a cylinder. The numerical test performed to evaluate mesh-resolution requirements and boundary-condition effectiveness are also described.

Atkins, H. L.↗

Numerical Simulation of the Generation of Axisymmetric Mode Jet Screech Tones

An imperfectly expanded supersonic jet, invariably, radiates both broadband noise and discrete frequency sound called screech tones. Screech tones are known to be generated by a feedback loop driven by the large scale instability waves of the jet flow. Inside the jet plume is a quasi-periodic shock cell structure. The interaction of the instability waves and the shock cell structure, as the former propagates through the latter, is responsible for the generation of the tones. Presently, there are formulas that can predict the tone frequency fairly accurately. However, there is no known way to predict the screech tone intensity. In this work, the screech phenomenon of an axisymmetric jet at low supersonic Mach number is reproduced by numerical simulation. The computed mean velocity profiles and the shock cell pressure distribution of the jet are found to be in good agreement with experimental measurements. The same is true with the simulated screech frequency. Calculated screech tone intensity and directivity at selected jet Mach number are reported in this paper. The present results demonstrate that numerical simulation using computational aeroacoustics methods offers not only a reliable way to determine the screech tone intensity and directivity but also an opportunity to study the physics and detailed mechanisms of the phenomenon by an entirely new approach.

Shen, Hao↗

Prediction of Sound Waves Propagating Through a Nozzle Without/With a Shock Wave Using the Space-Time CE/SE Method

The benchmark problems in Category 1 (Internal Propagation) of the third Computational Aeroacoustics (CAA) Work-shop sponsored by NASA Glenn Research Center are solved using the space-time conservation element and solution element (CE/SE) method. The first problem addresses the propagation of sound waves through a nearly choked transonic nozzle. The second one concerns shock-sound interaction in a supersonic nozzle. A quasi one-dimension CE/SE Euler solver for a nonuniform mesh is developed and employed to solve both problems. Numerical solutions are compared with the analytical solution for both problems. It is demonstrated that the CE/SE method is capable of solving aeroacoustic problems with/without shock waves in a simple way. Furthermore, the simple nonreflecting boundary condition used in the CE/SE method which is not based on the characteristic theory works very well.

Wang, Xiao-Yen↗

Sound Emission of Rotor Induced Deformations of Generator Casings

The casing of large electrical generators can be deformed slightly by the rotor's magnetic field. The sound emission produced by these periodic deformations, which could possibly exceed guaranteed noise emission limits, is analysed analytically and numerically. From the deformation of the casing, the normal velocity of the generator's surface is computed. Taking into account the corresponding symmetry, an analytical solution for the acoustic pressure outside the generator is round in terms of the Hankel function of second order. The normal velocity or the generator surface provides the required boundary condition for the acoustic pressure and determines the magnitude of pressure oscillations. For the numerical simulation, the nonlinear 2D Euler equations are formulated In a perturbation form for low Mach number Computational Aeroacoustics (CAA). The spatial derivatives are discretized by the classical sixth-order central interior scheme and a third-order boundary scheme. Spurious high frequency oscillations are damped by a characteristic-based artificial compression method (ACM) filter. The time derivatives are approximated by the classical 4th-order Runge-Kutta method. The numerical results are In excellent agreement with the analytical solution.

Polifke, W.↗

Entropy Splitting for High Order Numerical Simulation of Vortex Sound at Low Mach Numbers

A method of minimizing numerical errors, and improving nonlinear stability and accuracy associated with low Mach number computational aeroacoustics (CAA) is proposed. The method consists of two levels. From the governing equation level, we condition the Euler equations in two steps. The first step is to split the inviscid flux derivatives into a conservative and a non-conservative portion that satisfies a so called generalized energy estimate. This involves the symmetrization of the Euler equations via a transformation of variables that are functions of the physical entropy. Owing to the large disparity of acoustic and stagnation quantities in low Mach number aeroacoustics, the second step is to reformulate the split Euler equations in perturbation form with the new unknowns as the small changes of the conservative variables with respect to their large stagnation values. From the numerical scheme level, a stable sixth-order central interior scheme with a third-order boundary schemes that satisfies the discrete analogue of the integration-by-parts procedure used in the continuous energy estimate (summation-by-parts property) is employed.

Mueller, B.↗

Numerical Solutions to the Third CAA Workshop Benchmark Problems

This paper presents numerical solutions to the problems of propagation of sound waves through a transonic nozzle, shock-sound interactions, and automobile noise involving feedback of the third NASA Computational Aeroacoustics (CAA) Workshop on benchmark problems. The numerical algorithm is based on a dual time scheme for temporal discretization and a third-order finite volume scheme for spatial discretization. The aims of this study are to apply a dual time stepping scheme to treat aeroacoustic problems of sound propagation and to validate our CAA solver with the benchmark problems for developing a numerical tool for noise analysis and control.

Loh, Roy H.↗

Comparison of Analytical and CAA Solution(s): Unsteady Response of a Rectilinear Swept Cascade to an Incident Gust

This problem is aimed at testing the ability of computational aeroacoustic methods to calculate the acoustic pressure field generated as a result of the interaction of a convected harmonic gust with a rectilinear cascade of swept flat plates. All relevant geometric information is the same as the unsteady response of an isolated finite span swept airfoil to an incident gust problem with the exception of the sweep angle alpha, which is now taken to be a variable. Assume a cascade stagger angle of zero (with the x-axis aligned with the chord) and a gap-to-chord ratio of h/c = 1. The mean flow Mach number M(sub 0), gust frequency omega, gust amplitude A, and chordwise wavenumber k(sub x) are the same. But for the cascade problem take k(sub y) = pi and k(sub z) = 0. The appropriate physical scales are also the same, as is the governing equation and boundary conditions. Naturally, the impermeability condition now applies to the entire cascade. For this problem, determine the amplitude of upstream-radiated rms acoustic pressure as a function of the sweep angle at the specified frequency. Specifically, show the variations of rms acoustic pressure amplitude at the upstream location (-5c, 0, l/2) for sweep angle alpha in the range (0.0 deg, 30.0 deg). Use sweep angle increments no larger than 2.5 deg. Express the results in dB using the rms pressure value for alpha = 0.0 deg as the reference level. The appropriate FORTRAN output statement should read "WRITE(IUNIT,*) alpha, 20-log10 (p(sub rms)(alpha)/p(sub rms)(0))".

Envia, Edmane↗

High Order Numerical Simulation of Sound Generated by the Kirchhoff Vortex

An improved high order finite difference method for low Mach number computational aeroacoustics (CAA) is described. The improvements involve the conditioning of the Euler equations in perturbation form to minimize numerical cancellation error, and the use of a stable non-dissipative sixth-order central spatial differencing for the interior points and third-order at the boundary points. The spatial difference operator satisfies the summation-by-parts property to guarantee strict stability for linear hyperbolic systems. Spurious high frequency oscillations are damped by a third-order characteristic-based filter. The objective of this paper is to apply these improvements in the simulation of sound generated by the Kirchhoff vortex.

Mueller, Bernhard↗

Numerical Simulation of Screech Tones from Supersonic Jets: Physics and Prediction

The objectives of this project are to: (1) perform a numerical simulation of the jet screech phenomenon; and (2) use the data of the simulations to obtain a better understanding of the physics of jet screech. The original grant period was for three years. This was extended at no cost for an extra year to allow the principal investigator time to publish the results. We would like to report that our research work and results (supported by this grant) have fulfilled both objectives of the grant. The following is a summary of the important accomplishments: (1) We have now demonstrated that it is possible to perform accurate numerical simulations of the jet screech phenomenon. Both the axisymmetric case and the fully three-dimensional case were carried out successfully. It is worthwhile to note that this is the first time the screech tone phenomenon has been successfully simulated numerically; (2) All four screech modes were reproduced in the simulation. The computed screech frequencies and intensities were in good agreement with the NASA Langley Research Center data; (3) The staging phenomenon was reproduced in the simulation; (4) The effects of nozzle lip thickness and jet temperature were studied. Simulated tone frequencies at various nozzle lip thickness and jet temperature were found to agree well with experiments; (5) The simulated data were used to explain, for the first time, why there are two axisymmetric screech modes and two helical/flapping screech modes; (6) The simulated data were used to show that when two tones are observed, they co-exist rather than switching from one mode to the other, back and forth, as some previous investigators have suggested; and (7) Some resources of the grant were used to support the development of new computational aeroacoustics (CAA) methodology. (Our screech tone simulations have benefited because of the availability of these improved methods.)

Tam, Christopher K. W.↗

Multiresolution Wavelet Based Adaptive Numerical Dissipation Control for Shock-Turbulence Computations

The recently developed essentially fourth-order or higher low dissipative shock-capturing scheme of Yee, Sandham and Djomehri (1999) aimed at minimizing nu- merical dissipations for high speed compressible viscous flows containing shocks, shears and turbulence. To detect non smooth behavior and control the amount of numerical dissipation to be added, Yee et al. employed an artificial compression method (ACM) of Harten (1978) but utilize it in an entirely different context than Harten originally intended. The ACM sensor consists of two tuning parameters and is highly physical problem dependent. To minimize the tuning of parameters and physical problem dependence, new sensors with improved detection properties are proposed. The new sensors are derived from utilizing appropriate non-orthogonal wavelet basis functions and they can be used to completely switch to the extra numerical dissipation outside shock layers. The non-dissipative spatial base scheme of arbitrarily high order of accuracy can be maintained without compromising its stability at all parts of the domain where the solution is smooth. Two types of redundant non-orthogonal wavelet basis functions are considered. One is the B-spline wavelet (Mallat & Zhong 1992) used by Gerritsen and Olsson (1996) in an adaptive mesh refinement method, to determine regions where re nement should be done. The other is the modification of the multiresolution method of Harten (1995) by converting it to a new, redundant, non-orthogonal wavelet. The wavelet sensor is then obtained by computing the estimated Lipschitz exponent of a chosen physical quantity (or vector) to be sensed on a chosen wavelet basis function. Both wavelet sensors can be viewed as dual purpose adaptive methods leading to dynamic numerical dissipation control and improved grid adaptation indicators. Consequently, they are useful not only for shock-turbulence computations but also for computational aeroacoustics and numerical combustion. In addition, these sensors are scheme independent and can be stand alone options for numerical algorithm other than the Yee et al. scheme.

Sjoegreen, B.↗

A Model for Shear Layer Effects on Engine Noise Radiation

Prediction of aircraft engine noise is an important aspect of addressing the issues of community noise and cabin noise control. The development of physics based methodologies for performing such predictions has been a focus of Computational Aeroacoustics (CAA). A recent example of code development in this area is the ducted fan noise propagation and radiation code CDUCT-LaRC. Included within the code is a duct radiation model that is based on the solution of FfowcsWilliams-Hawkings (FW-H) equation with a penetrable data surface. Testing of this equation for many acoustic problems has shown it to provide generally better results than the Kirchhoff formula for moving surfaces. Currently, the data surface is taken to be the inlet or exhaust plane for inlet or aft-fan cases, respectively. While this provides reasonable results in many situations, these choices of data surface location lead to a few limitations. For example, the shear layer between the bypass ow and external stream can refract the sound waves radiated to the far field. Radiation results can be improved by including this effect, as well as the rejection of the sound in the bypass region from the solid surface external to the bypass duct surrounding the core ow. This work describes the implementation, and possible approximation, of a shear layer boundary condition within CDUCT-LaRC. An example application also illustrates the improvements that this extension offers for predicting noise radiation from complex inlet and bypass duct geometries, thereby providing a means to evaluate external treatments in the vicinity of the bypass duct exhaust plane.

Nark, Douglas M.↗

Technologies for Turbofan Noise Reduction

An overview presentation of NASA's engine noise research since 1992 is given for subsonic commercial aircraft applications. Highlights are included from the Advanced Subsonic Technology (AST) Noise Reduction Program and the Quiet Aircraft Technology (QAT) project with emphasis on engine source noise reduction. Noise reduction goals for 10 EPNdB by 207 and 20 EPNdB by 2022 are reviewed. Fan and jet noise technologies are highlighted from the AST program including higher bypass ratio propulsion, scarf inlets, forward-swept fans, swept/leaned stators, chevron nozzles, noise prediction methods, and active noise control for fans. Source diagnostic tests for fans and jets that have been completed over the past few years are presented showing how new flow measurement methods such as Particle Image Velocimetry (PIV) have played a key role in understanding turbulence, the noise generation process, and how to improve noise prediction methods. Tests focused on source decomposition have helped identify which engine components need further noise reduction. The role of Computational AeroAcoustics (CAA) for fan noise prediction is presented. Advanced noise reduction methods such as Hershel-Quincke tubes and trailing edge blowing for fan noise that are currently being pursued n the QAT program are also presented. Highlights are shown form engine validation and flight demonstrations that were done in the late 1990's with Pratt & Whitney on their PW4098 engine and Honeywell on their TFE-731-60 engine. Finally, future propulsion configurations currently being studied that show promise towards meeting NASA's long term goal of 20 dB noise reduction are shown including a Dual Fan Engine concept on a Blended Wing Body aircraft.

Huff, Dennis↗