Farassat, F. — Page 3 | Engineering Papers

NASA NTRS1999

A Simple Derivation of the Acoustic Boundary Condition in the Presence of Flow

In modeling the noise propagation in a duct of a ducted fan engine, the authors needed the boundary condition satisfied by the acoustic pressure on a lined duct wall in the presence of flow. The results, in the form used by the authors, was derived by M.K. Myers, who presented a formal derivation of a related results by K. Taylor. Here, a brief and simple derivation of the acoustic boundary condition of Myers is given. The main difference in the derivation from that by Myers is that the Gaussian coordinate system (q(exp 1, q(exp 2), q(exp 3)(q(exp 1, q(exp 2), t)) is used to specify the instantaneous position, i.e. Lagrangian variable, of the point on the mean position of the wall with curvilinear (Gaussian) coordinates (q(exp 1), q(exp 2)). Myers uses a locally orthogonal coordinate system which is somewhat less specific than what is used here.

Farassat, F.↗

NASA NTRS1999

A Simple Derivation of the Acoustic Boundary Condition in the Presence of Flow

In modelling the noise propagation in a duct of a ducted fan engine, the authors needed the boundary condition satisfied by the acoustic pressure on a lined duct wall in the presence of flow. The result in the form used by the authors was derived by M. K. Myers who presented a formal derivation of a related result by K. Taylor. A good review of the subject can be found in reference. Here, a brief and simple derivation of the acoustic boundary condition of Myers is given. The main difference in the derivation from that in reference is that the Gaussian co-ordinate system (q(sup 1), q(sup 2), q(sup 3)(q(sup 1), q(sup 2), t)) is used to specify the instantaneous position, i.e., Lagrangian variable, of the point on the mean position of the wall with curvilinear (Gaussian) co-ordinates (q(sup 1), q(sup 2)). Myers uses a locally orthogonal coordinate system which is somewhat less specific than what is used here.

Farassat, F.↗

NASA NTRS1999

The Application of a Boundary Integral Equation Method to the Prediction of Ducted Fan Engine Noise

The prediction of ducted fan engine noise using a boundary integral equation method (BIEM) is considered. Governing equations for the BIEM are based on linearized acoustics and describe the scattering of incident sound by a thin, finite-length cylindrical duct in the presence of a uniform axial inflow. A classical boundary value problem (BVP) is derived that includes an axisymmetric, locally reacting liner on the duct interior. Using potential theory, the BVP is recast as a system of hypersingular boundary integral equations with subsidiary conditions. We describe the integral equation derivation and solution procedure in detail. The development of the computationally efficient ducted fan noise prediction program TBIEM3D, which implements the BIEM, and its utility in conducting parametric noise reduction studies are discussed. Unlike prediction methods based on spinning mode eigenfunction expansions, the BIEM does not require the decomposition of the interior acoustic field into its radial and axial components which, for the liner case, avoids the solution of a difficult complex eigenvalue problem. Numerical spectral studies are presented to illustrate the nexus between the eigenfunction expansion representation and BIEM results. We demonstrate BIEM liner capability by examining radiation patterns for several cases of practical interest.

Dunn, M. H.↗

NASA NTRS1999

Acoustic Energy Estimates in Inhomogeneous Moving Media

In ducted fan engine noise research, there is a need for defining a simple and easy to use acoustic energy conservation law to help in quantification of noise control techniques. There is a well known conservation law relating acoustic energy and acoustic energy flux in the case of an isentropic irrotational flow. Several different approaches have been taken to generalize this conservation law. For example, Morfey finds an identity by separating out the irrotational part of the perturbed flow. Myers is able to find a series of indentities by observing an algebraic relationship between the basic conservation of energy equation for a background flow and the underlying equations of motion. In an approximate sense, this algebraic relationship is preserved under perturbation. A third approach which seems to have not been pursued in the literature is a result known as Noether's theorem. There is a Lagrangian formulation for the Euler equation of fluid mechanics. Noether's theorem says that any group action that leaves the Lagrangian action invariant leads to a conserved quantity. This presentation will include a survey of current results regarding acoustic energy and preliminary results on the symmetries of the Lagrangian.

Farassat, F.↗

NASA NTRS1999

Two-Dimensional Fourier Transform Analysis of Helicopter Flyover Noise

A method to separate main rotor and tail rotor noise from a helicopter in flight is explored. Being the sum of two periodic signals of disproportionate, or incommensurate frequencies, helicopter noise is neither periodic nor stationary. The single Fourier transform divides signal energy into frequency bins of equal size. Incommensurate frequencies are therefore not adequately represented by any one chosen data block size. A two-dimensional Fourier analysis method is used to separate main rotor and tail rotor noise. The two-dimensional spectral analysis method is first applied to simulated signals. This initial analysis gives an idea of the characteristics of the two-dimensional autocorrelations and spectra. Data from a helicopter flight test is analyzed in two dimensions. The test aircraft are a Boeing MD902 Explorer (no tail rotor) and a Sikorsky S-76 (4-bladed tail rotor). The results show that the main rotor and tail rotor signals can indeed be separated in the two-dimensional Fourier transform spectrum. The separation occurs along the diagonals associated with the frequencies of interest. These diagonals are individual spectra containing only information related to one particular frequency.

SantaMaria, Odilyn L.↗

NASA NTRS1999

A Mode Detection Method Using the Azimuthal Directivity of a Turbofan Model

The azimuthal, far field directivity of a scale fan model was measured in high resolution. The model is a 12 inch diameter rotor with 16 blades followed by 40 stator vanes. The tests were conducted at the nominal 100% speed corresponding to a tip speed of 905 ft/sec. Measurement of the radiated sound field, forward of the fan, was made in an anechoic chamber with an inflow control device and a baffle separating the aft and forward radiated interaction noise. The acoustic field was surveyed with a circular hoop array of 16 microphones which was moved to 14 axial stations. At each axial station the hoop was rotated in half-degree increments to take 736 points in the azimuthal angle. In addition to sound pressure level, the phase angle relative to a reference microphone was measured at each point. The sound pressure level is shown to vary in patterns by 10-15 dB especially for the fundamental tone but also for the first and second harmonic. A far field mode detection method has been developed and used with the data which determines the modes generated by the fan and which then interact to form the azimuthal directivity.

Thomas, R. H.↗

NASA NTRS1999

Verification and Analysis of Formulation 4 of Langley for the Study of Noise From High Speed Surfaces

There are several approaches to the prediction of the noise from sources on high speed surfaces. Two of these are the Kirchhoff and the Ffowcs williams-Hawkings methods. It can be shown that both of these methods depend on the solution of the wave equation with mathematically similar inhomogeneous source terms. Two subsonic solutions known as Formulation 1 and 1A of Langley are simple and efficient for noise prediction. The supersonic solution known as Formulation 3 is very complicated and difficult to code. Because of the complexity of the result, the computation time is longer than the subsonic formulas. Furthermore, it is difficult to assess the accuracy of noise prediction. We have been searching for a new and simpler supersonic formulation without these shortcomings. In the last AIAA Aeroacoustics Conference in Toulouse, Farassat, Dunn and Brentner presented a paper in which such a result was presented and called Formulation 4 of Langley. In this paper we will present two analytic tests of the validity this Formulation: 1) the noise from dipole distribution on the unit circle whose strength varies radially with the square of the distance from the center and 2) the noise from dipole distribution on the unit sphere whose strength varies with the cosine of the angle from the polar axis. We will discuss the question of singularities of Formulation 4.

Farassat, F.↗

NASA NTRS1999

Two-Dimensional Fourier Transform Analysis of Helicopter Flyover Noise

A method to separate main rotor and tail rotor noise from a helicopter in flight is explored. Being the sum of two periodic signals of disproportionate, or incommensurate frequencies, helicopter noise is neither periodic nor stationary. The single Fourier transform divides signal energy into frequency bins of equal size. Incommensurate frequencies are therefore not adequately represented by any one chosen data block size. A two-dimensional Fourier analysis method is used to separate main rotor and tail rotor noise. The two-dimensional spectral analysis method is first applied to simulated signals. This initial analysis gives an idea of the characteristics of the two-dimensional autocorrelations and spectra. Data from a helicopter flight test is analyzed in two dimensions. The test aircraft are a Boeing MD902 Explorer (no tail rotor) and a Sikorsky S-76 (4-bladed tail rotor). The results show that the main rotor and tail rotor signals can indeed be separated in the two-dimensional Fourier transform spectrum. The separation occurs along the diagonals associated with the frequencies of interest. These diagonals are individual spectra containing only information related to one particular frequency.

SantaMaria, Odilyn L.↗

NASA NTRS1999

Two-Dimensional Fourier Transform Analysis of Helicopter Flyover Noise

A method to separate main rotor and tail rotor noise from a helicopter in flight is explored. Being the sum of two periodic signals of disproportionate, or incommensurate frequencies, helicopter noise is neither periodic nor stationary. The single Fourier transform divides signal energy into frequency bins of equal size. Incommensurate frequencies are therefore not adequately represented by any one chosen data block size. A two-dimensional Fourier analysis method is used to separate main rotor and tail rotor noise. The two-dimensional spectral analysis method is first applied to simulated signals. This initial analysis gives an idea of the characteristics of the two-dimensional autocorrelations and spectra. Data from a helicopter flight test is analyzed in two dimensions. The test aircraft are a Boeing MD902 Explorer (no tail rotor) and a Sikorsky S-76 (4-bladed tail rotor). The results show that the main rotor and tail rotor signals can indeed be separated in the two-dimensional Fourier transform spectrum. The separation occurs along the diagonals associated with the frequencies of interest. These diagonals are individual spectra containing only information related to one particular frequency.

SantaMaria, Odilyn L.↗

NASA NTRS1998

Ducted-Fan Engine Acoustic Predictions using a Navier-Stokes Code

A Navier-Stokes computer code is used to predict one of the ducted-fan engine acoustic modes that results from rotor-wake/stator-blade interaction. A patched sliding-zone interface is employed to pass information between the moving rotor row and the stationary stator row. The code produces averaged aerodynamic results downstream of the rotor that agree well with a widely used average-passage code. The acoustic mode of interest is generated successfully by the code and is propagated well upstream of the rotor; temporal and spatial numerical resolution are fine enough such that attenuation of the signal is small. Two acoustic codes are used to find the far-field noise. Near-field propagation is computed by using Eversman's wave envelope code, which is based on a finite-element model. Propagation to the far field is accomplished by using the Kirchhoff formula for moving surfaces with the results of the wave envelope code as input data. Comparison of measured and computed far-field noise levels show fair agreement in the range of directivity angles where the peak radiation lobes from the inlet are observed. Although only a single acoustic mode is targeted in this study, the main conclusion is a proof-of-concept: Navier-Stokes codes can be used both to generate and propagate rotor/stator acoustic modes forward through an engine, where the results can be coupled to other far-field noise prediction codes.

Rumsey, C. L.↗

NASA NTRS1998

Ducted-Fan Engine Acoustic Predictions Using a Navier-Stokes Code

A Navier-Stokes computer code is used to predict one of the ducted-fan engine acoustic modes that results from rotor-wake/stator-blade interaction. A patched sliding-zone interface is employed to pass information between the moving rotor row and the stationary stator row. The code produces averaged aerodynamic results downstream of the rotor that agree well with a widely used average-passage code. The acoustic mode of interest is generated successfully by the code and is propagated well upstream of the rotor, temporal and spatial numerical resolution are fine enough such that attenuation of the signal is small. Two acoustic codes are used to find the far-field noise. Near-field propagation is computed by using Eversman's wave envelope code, which is based on a finite-element model. Propagation to the far field is accomplished by using the Kirchhoff formula for moving surfaces with the results of the wave envelope code as input data. Comparison of measured and computed far-field noise levels show fair agreement in the range of directivity angles where the peak radiation lobes from the inlet are observed. Although only a single acoustic mode is targeted in this study, the main conclusion is a proof-of-concept: Navier Stokes codes can be used both to generate and propagate rotor-stator acoustic modes forward through an engine, where the results can be coupled to other far-field noise prediction codes.

Rumsey, C. L.↗

NASA NTRS1998

Liner Optimization Studies Using the Ducted Fan Noise Prediction Code TBIEM3D

In this paper we demonstrate the usefulness of the ducted fan noise prediction code TBIEM3D as a liner optimization design tool. Boundary conditions on the interior duct wall allow for hard walls or a locally reacting liner with axially segmented, circumferentially uniform impedance. Two liner optimization studies are considered in which farfield noise attenuation due to the presence of a liner is maximized by adjusting the liner impedance. In the first example, the dependence of optimal liner impedance on frequency and liner length is examined. Results show that both the optimal impedance and attenuation levels are significantly influenced by liner length and frequency. In the second example, TBIEM3D is used to compare radiated sound pressure levels between optimal and non-optimal liner cases at conditions designed to simulate take-off. It is shown that significant noise reduction is achieved for most of the sound field by selecting the optimal or near optimal liner impedance. Our results also indicate that there is relatively large region of the impedance plane over which optimal or near optimal liner behavior is attainable. This is an important conclusion for the designer since there are variations in liner characteristics due to manufacturing imprecisions.

Dunn, M. H.↗

NASA NTRS1998

The Application of Nonstandard Analysis to the Study of Inviscid Shock Wave Jump Conditions

The use of conservation laws in nonconservative form for deriving shock jump conditions by Schwartz distribution theory leads to ambiguous products of generalized functions. Nonstandard analysis is used to define a class of Heaviside functions where the jump from zero to one occurs on an infinitesimal interval. These Heaviside functions differ by their microstructure near x = 0, i.e., by the nature of the rise within the infinitesimal interval it is shown that the conservation laws in nonconservative form can relate the different Heaviside functions used to define jumps in different flow parameters. There are no mathematical or logical ambiguities in the derivation of the jump conditions. An important result is that the microstructure of the Heaviside function of the jump in entropy has a positive peak greater than one within the infinitesimal interval where the jump occurs. This phenomena is known from more sophisticated studies of the structure of shock waves using viscous fluid assumption. However, the present analysis is simpler and more direct.

Farassat, F.↗

NASA NTRS1997

Azimuthal Patterns of the Radiated Sound Field from a Turbofan Model

The azimuthal directivity of a scale fan model was measured extensively. The model is a 12 inch diameter fan with 16 rotors and 40 stator vanes and tests were done at a tip speed of 905 ft/sec. Tests were conducted in an anechoic chamber with an inflow control device installed on the stationary fan model. The acoustic far field of the fan was surveyed with a circular hoop, with a diameter of six fan diameters, centered on the fan axis and was moved along the fan axis at polar angles from 20 to 110 degrees in increments of 10 degrees. The hoop, with 16 microphones evenly spaced at intervals of 22.5 degrees was rotated in 24 increments in the azimuthal direction for a total 384 points. From this extensive mapping of the directivity it is shown that the azimuthal directivity of the fundamental and first two harmonics is significant and can vary up to 15 dB. The broadband can also have an azimuthal directivity with as much as a 4 dB variation. A theory is proposed with relates the radiated modes with the generation of the far field patterns which produce the azimuthal directivity.

Thomas, R. H.↗

NASA NTRS1997

Inlet Shape Effects on the Far-Field Sound of a Model Fan

A wind tunnel test was conducted to determine the effects of inlet shape on fan radiated noise. Four inlet geometries, which included a long standard flight type inlet, a short, aggressive flight inlet a scarf inlet, and an elliptical inlet were investigated in the study. The fan model used in the study was a 0.1 scale of the Pratt and Whitney Advanced Ducted Propeller (ADP), an ultra high bypass ratio turbofan engine. Acoustic data are presented for a fan speed of 70% (12,000 rpm) and a tunnel speed of 0.10 Mach number, The fan was configured with a 16-bladed rotor and a 40 stator vane set that were separated by 2.0 chord lengths. The radiated noise was measured with 15 microphones on a boom that traversed the length of the tunnel test section. Data from these microphones are presented in the form of sideline angle directivity plots. Noise associated with the test inlets was also predicted using a ray acoustics code. Inlet shape has been found to have a significant effect on both tone and broadband noise, and the non-axisymmetric inlet shape can be used for a noise reduction method.

Clark, L. R.↗

NASA NTRS1997

An Analytical Comparison of the Acoustic Analogy and Kirchhoff Formulation for Moving Surfaces

The Lighthill acoustic analogy, as embodied in the Ffowcs Williams-Hawkings (FW-H) equation, is compared with the Kirchhoff formulation for moving surfaces. A comparison of the two governing equations reveals that the main Kirchhoff advantage (namely nonlinear flow effects are included in the surface integration) is also available to the FW-H method if the integration surface used in the FW-H equation is not assumed impenetrable. The FW-H equation is analytically superior for aeroacoustics because it is based upon the conservation laws of fluid mechanics rather than the wave equation. This means that the FW-H equation is valid even if the integration surface is in the nonlinear region. This is demonstrated numerically in the paper. The Kirchhoff approach can lead to substantial errors if the integration surface is not positioned in the linear region. These errors may be hard to identify. Finally, new metrics based on the Sobolev norm are introduced which may be used to compare input data for both quadrupole noise calculations and Kirchhoff noise predictions.

Brentner, Kenneth S.↗

NASA NTRS1997

Supersonic Quadrupole Noise Theory for High-Speed Helicopter Rotors

High-speed helicopter rotor impulsive noise prediction is an important problem of aeroacoustics. The deterministic quadrupoles have been shown to contribute significantly to high-speed impulsive (HSI) noise of rotors, particularly when the phenomenon of delocalization occurs. At high rotor-tip speeds, some of the quadrupole sources lie outside the sonic circle and move at supersonic speed. Brentner has given a formulation suitable for efficient prediction of quadrupole noise inside the sonic circle. In this paper, we give a simple formulation based on the acoustic analogy that is valid for both subsonic and supersonic quadrupole noise prediction. Like the formulation of Brentner, the model is exact for an observer in the far field and in the rotor plane and is approximate elsewhere. We give the full analytic derivation of this formulation in the paper. We present the method of implementation on a computer for supersonic quadrupoles using marching cubes for constructing the influence surface (Sigma surface) of an observer space- time variable (x; t). We then present several examples of noise prediction for both subsonic and supersonic quadrupoles. It is shown that in the case of transonic flow over rotor blades, the inclusion of the supersonic quadrupoles improves the prediction of the acoustic pressure signature. We show the equivalence of the new formulation to that of Brentner for subsonic quadrupoles. It is shown that the regions of high quadrupole source strength are primarily produced by the shock surface and the flow over the leading edge of the rotor. The primary role of the supersonic quadrupoles is to increase the width of a strong acoustic signal.

Farassat, F.↗

NASA NTRS1996

The Kirchhoff Formulas for Moving Surfaces in Aeroacoustics - The Subsonic and Supersonic Cases

One of the active areas of computational aeroacoustics is the application of the Kirchhoff formulas to the problems of the rotating machinery noise predictions. The original Kirchhoff formula was derived for a stationary surface. In 1988, Farassat and Myers derived a Kirchhoff Formula obtained originally by Morgans using modem mathematics. These authors gave a formula particularly useful for applications in aeroacoustics. This formula is for a surface moving at subsonic speed. Later in 1995 these authors derived the Kirchhoff formula for a super-sonically moving surface. This technical memorandum presents the viewgraphs of a day long workshop by the author on the derivation of the Kirchhoff formulas. All necessary background mathematics such as differential geometry and multidimensional generalized function theory are discussed in these viewgraphs. Abstraction is kept at minimum level here. These viewgraphs are also suitable for understanding the derivation and obtaining the solutions of the Ffowcs Williams-Hawkings equation. In the first part of this memorandum, some introductory remarks are made on generalized functions, the derivation of the Kirchhoff formulas and the development and validation of Kirchhoff codes. Separate lists of references by Lyrintzis, Long, Strawn and their co-workers are given in this memorandum. This publication is aimed at graduate students, physicists and engineers who are in need of the understanding and applications of the Kirchhoff formulas in acoustics and electromagnetics.

Farassat, F.↗