Acoustics
CAA: Computer Aided Acoustics
Aeroacoustics deals with production, propagation (transmission: airborne or solidborne), reflection (as sound moves from one medium to another), absorption, scattering, attenuation and reception (by human ears) of sound waves. The sound is primarily generated by fluid fluctuations forced due to unsteady motion of fluids. Sound is a small perturbation of pressure over a mean (steady state) pressure, p’/p_{0} which propagates as a wave. Aeroacoustics is different from the study on generation of sound such as organ pipe and loudspeakers. The later is referred as "classical acoustics".
Excerpts from COMSOL user manual: "Acoustics is the physics of sound. Sound is the sensation, as detected by the ear, of very small rapid changes in the acoustic pressure p above and below a static value. This static value is the atmospheric pressure (about 100,000 Pa). The acoustic pressure variations are typically described as pressure waves propagating in space and time." During their propagation, sound waves can be reflected, refracted or attentuated by the medium and can have interference with sounds of different frequencies.
Concepts in Acoustics
A detailed mathematical approach to acoustics can be found in this ebook available free for personal and educational purposes.
 Threshold of Hearing: This it the lowest level of sound intensity that can be sensed (heard) by most humans at a specified frequency of 1000 [Hz]. In air it is 20 [μPa]. The threshold in terms of frequency is the range of 20  20,000 [Hz] (threshold of hearing: 0.02 mPa, loud noise: 20 Pa in rockbands, pain is experienced as sound pressure ~ 60 Pa). The peak of the human hearing range is between 1000 ~ 3000 [Hz] (reference STARCCM+ tutorial document). Analogous to light wave, sound waves with frequency < 20 [Hz] are called infrasonic (similar to infrared) and that of frequency > 20,000 [Hz] are called ultrasonic (similar to ultraviolet).
 Description of waves: frequency  ν is number of vibrations per seconds, wavelength  λ is the distance after which wave pattern repeats itself, speed of sound c  the velocity of propagation of acoustic waves = ν.λ, wave number  k is the number of waves required to cover a specified distance = 2π/λ. Acoustic waves are characterized by a much lower amplitudes and higher frequency than the pressure waves. Frequency is also known as 'pitch'. For example, a man's voice is lowpitched (fundamental frequency in the band 65 Hz ~ 250 Hz with mean around 120 Hz) and a female's voice is termed as highpitched (fundamental frequency band 100 ~ 500 Hz with mean around 210 Hz). Approximate sound speeds for 3 media: air, water and steel are 340 [m/s], 1500 [m/s] and 5000 [m/s] respectively.
A linear, timeharmonic wave equation in 3D is given by ∇^{2}p' + k^{2}p' = 0 where p' is the acoustic pressure. This is also known as Helmholtz equation for sound pressure.
 Particle displacement of a wave traveling in xdirection: y = y_{0} sin(kx  ωt)
 Particle velocity: u' = dt/dt = ωy_{0} cos(kx  ωt)
 Acoustic pressure: p' = B dy/dx where B = bulk or volumetric modulus of elasticity. Thus, p' = B.y_{0}k cos(kx  ωt)
 In most of the textbooks, for the sake of brevity wave equations or wave functions are expressed as complex numbers. For example, the general form of an onedimensional wave function is ψ = A cos(ωt  kx  φ) ≡ A cos(kx  ωt + φ). The complex wave function is ψ = ψ_{0} e^{i[kx  ωt]} with complext constant ψ_{0} = e^{iφ}.
 All these representations complex wave function are based on fundamental mathematical identity, known as Euler's theorem expressed as e^{iφ} ≡ cos(φ) + i sin(φ).
 The main advantage of the complex representation over direct real representation is single complex amplitude for for both real amplitude and phase angle. This is just a way to present 'real' phenomena using what is called imaginary numbers.
 Air Particle Velocity  Acoustics waves are called longitudinal waves which refers to the phenomena that the particles of medium osciallates along the direction of the motion of the wave. Thus, air particles will vibrate about its mean position along the direction of wave propagation. It is usually expressed as RMS (Root Mean Square) velocity [m/s]. Other approach to express it is "volume velocity" for flow inside ducts which equals the product of the rms velocity of the air particles and the crosssectional area of the tube and has units of [m^{3}/s].
For a given pressure field p'(x, t), the particle velocity is given by u'(x, t) = p'(x, t)/[ρ·c]
 Noise vs. music: Noise  unwanted sound. Music  pleasant sound. Thus, human perception is important aspect of acoustic studies. Rattle is another type of noise attributed to gears.
 Near Field  This is one of the ways to characterize the qualities of the sound field created by a source in its vicinity. The near field refers to zones closer to the source (typically involving spherical waves and the consequent acoustical behavior). Cabin noise level prediction due to HVAC louvres and flow over outside rear view mirrors (ORVM) are examples of near filed where receiver is closer to the source. The equation of a spherical waves is given by
p' = S·k·c·ρ·/4πr·sin[k(ct  r)], u' = S·k/4πr [1/k/r·cos[k(ct  r)  sin[k(ct  r)], S = maximum rate of fluid emission from source, ρ is the density of the medium, r is radial distance from source. For spherical waves, the intensity decreases with the square of the distance but the sound pressure level decreases linearly with distance.
 Far Field  This is one of the way to characterize qualities of a sound field created by a source far away from it. The far field involves planar waves and the consequent acoustical behavior. Noise sevaral hundreds of chord length away from a wing of aircraft is an example of farfield situation where flow field may not be required to be modeled all the way till receiver. The dividing line betwen the two conditions is wavelength (frequency) dependent, in other words the location which signifies transition from nearfield to farfield will be depenent of frequency of noise generated. The equation for plane waves is given by
p' = a·k·c^{2}·ρ·sin[k(ct  x)], u' = a·k·c·sin[k(ct  x)], a = particle displacement, k = wave number.
Low frequency sounds have long wavelengths, 340 [m/s] / 50 [Hz] = 6.8 [m] at 50 Hz while the wavelength at 10 kHz is only 34 [mm]. Thus, high frequency waves can be considered planar only after a small fraction of a meter from the source while low frequency sounds require few meters of separation from the source to be safely considered planar. Real sounds or noise generated from a source may contain many component frequencies leading to a mixture of planar and spherical sources. This makes the measurements and calculation of acoustic field complicated. In theory, plane wave pressures do not decrease with distance as they do for spherical waves.
 Free Field and Diffuse Field The terms free field (also called direct field) and diffuse field (also known as reverberant field) refers to sounds measured in a closed zone such as a room. Near the source, the direct sound is dominant while further away the reverberant sound (sound reflected by the walls) prevails.
In direct field (like near field situation), sound pressure decreases linearly with distance while there is no decrease in sound pressure with distance in the diffuse field (as observed in plane waves). The distance at which the direct and reverberant sounds are at equal pressure is called the critical distance.
 SPL  Sound Pressure Level: This refers to the intensity of sound waves. This is represented as logarithm (because sound pressures have large ranges) of a ratio of magnitude of pressure fluctuations to the reference pressure at threshold of hearing and denoted as deciBel [dB]. SPL = 20 × log_{10}(p / p_{REF}) [dB] where p_{REF} = 20 [μPa]. Thus, "threshold of hearing" corresponds of SPL = 0 [dB]. In other form: SPL [dB] = 20 × log_{10}[p_{RMS}] + 94 [dB] where p_{RMS} is in [Pa]. Bel is another measure of SPL and is equal to 10 deciBels (dB). It is dimensionless since it is a ratio of two quantities. Threshold of pain is ~ 130 dB (such as gunshots, jet engines) and higher levels of the order of 200 dB may occur in explosions, sonic booms, rocket liftoff and similar extraordinarily loud events. Whisper has SPL of 30 dB. Disco and rock concerts has noise level in the range 100 ~ 110 dB.
 For a point source, the sound pressure ∝ 1/r from the measuring point to the sound source. Thus, doubling of the distance, sound pressure will decrease by half.
 log_{10}(2) = 0.301, thus doubling the distance, SPL will drop by 20 × 0.3 = 6 [dB]
 Due to the fact that SPL is based on a log scale, they cannot be added like a simple algebraic summation. As expressed below, adding two sound pressure levels of equal value will always result in an increase of 3 [dB].
 Thus, the maximum increase by combination two random noise sources is always < 3 [dB] above the noisier source.
 When the difference in noise level of two source is ≥ 10 [dB], the contribution of weaker noise source is negligible over stronger noise source.
 Example: if a plant have 5 exhaust fans each radiating sound pressure level of 60 [dB]. The sound pressure levels expected with various number of running fans (assuming the noise to be uncorrelated) are:
 2 fans: 60.0 [dB] + 60.0 [dB] = 63.0 [dB]
 3 fans: 63.0 [dB] + 60.0 [dB] = 64.8 [dB]
 4 fans: 63.0 [dB] + 63.0 [dB] = 66.0 [dB]
 5 fans: 66.0 [dB] + 60.0 [dB] = 67.0 [dB]
 If SPL_{1} < SPL_{2} and SPL_{1}  SPL_{2} = δ following table can be used to find addition of decibel levels, SPL_{3} = SPL_{2} + Δ:
δ  0.0  1.0  2.0  3.0  4.0  5.0  6.0  7.0  8.0  9.0  10.0 
Δ  3.0  2.5  2.1  1.8  1.5  1.2  1.0  0.8  0.6  0.5  0.4 
 Sound Intensity  I: Also known as soundenergy flux density, it is the timeaveraged value of "[instantaneous sound pressure × particle velocity]" at a point in a temporally stationary sound field (sound levels not varying with time). The particle velocity is the speed at which air molecules vibrate back and forth while transmitting the sound waves. Particle velocity is a vector quantity, while sound pressure is a scalar. Thus, sound intensity is a vector quantity having unit of [W/m^{2}]. Sound pressure level can be measured with a single microphone while sound intensity measurement requires two or more microphones in a specific arrangement.
In a freefield that is where no reflected sound waves are present and far away from a noise source, I = p^{2}_{RMS}/ρc where p_{RMS} is Root Mean Square acoustic pressure.
A similar concept is sound source strength which is the maximum rate of volume displacement [m^{3}/s] produced by a source when emitting a harmonic sound wave.
 Normal sound intensity level  L_{IN}: It equals 10 times the logarithm (base 10) of the ratio of the unsigned value of the normal component of the sound intensity (which is radiated by the sound source under test and determined in a direction perpendicular to the measurement surface) to the reference sound intensity. Like other logarithim definitions, it is expressed in decibels and the reference sound intensity is 10^{–12} [W/m^{2}] or 1 [pW/m^{2}].
 L_{IN} [dB] = 10 × log_{10}[I] + 120 [dB] where I is in [W/m^{2}].
 The sound intensity decreases inversely proportional to the squared distance, that is I ∝ 1/r^{2} from the measuring point to the sound source. Thus, doubling of the distance deceases the sound intensity by 4times of its initial value.
 Sound pressure (sound field) is not same as sound intensity (sound energy quantity). I ∝ p^{2}.
 Analogy between heat and sound  the sound pressure level, sound power level and sound intensity are analogous to temperature, wattage and heat flux respectively around an eletrical bulb.
 Band and Bandwidth  Band refers to a continuous range of frequencies between two limiting frequencies. Bandwidth is defined as range of frequencies, usually of standard size in acoustics. For example octave or onethird octave bands where octave bands are groups of frequencies named by the center frequency having upper limit always twice the lower limit of the range, f_{c} = [f_{L}.f_{U}]^{0.5}. The lower and upper frequencies are also known as the 3 dB or halfpower points. Based on the perception of sound by a human ear, the reference middle frequency in acoustics is 1000 [Hz] and octave bands are defined with respect to this value such as (15.6, 31.3, 62.5, 125, 250, 500, 1000, 2000, 4000, 8000, 16000) Hz.
 Centre or middle frequency value f_{0} of octave band is obtained by successive division and multipliation of 1000 [Hz] by 2. The upper band = f_{0} × 2^{1/2} and the lower band = f_{0} / 2^{1/2}
 A norder octave band (such as 1/3 octave) is obtained by successive division and multipliation of 1000 [Hz] by 2^{1/n} (e.g. 2^{1/3} for 1/3 octave). The upper band = f_{0} × 2^{1/2n} and the lower band = f_{0} / 2^{1/2n}.
 Loudness: Phons: Loudness is a psychological measure of physical strength (amplitude) of sound dependent collectively on human ear's perception along with sound pressure, frequency, bandwidth and duration. Phon is sound intensity level as a reference frequency of 1000 [Hz]. Two sounds or noise with equal intensity but at different frequencies will not have same loudness. A sound/noise of 75 dB at 1 kHz = 75 phons.
 Designated as Lw, Sound Power Level is the total acoustic energy output of a noise source independent of medium of propagation of the sound (the environment). Lw = 10 × log_{10}(P/P_{0}) where P_{0} = 1 [pW] = 10^{12} [W]. While sound power level is a constant value for each noise source, SPL are dependent on factors such as the distance from the noise source, presence of any reflective surfaces present near the source. Thus, sound pressure levels will always be higher than sound power levels. Sometimes a frequency weighting as per IEC 61672–1 and/or specific frequency bands are applied, this is denoted the Aweighted sound power level  L_{WA} and its unit is denoted by dB_{A} or dB(A).
L_{W} [dB] = 10 × log_{10}W + 120 [dB] where W is in [Watts]. Comparing noise source as a charge, as per Gauss's law: the total electric flux through a closed surface is zero if no charge is enclosed by the surface. Similarly, the total acoustic power through a closed surface is zero if no noise source is enclosed by the surface. Thus, sound power level measurement will not be affected by background noise if the measurement surface encompassing the noise source does not enclose source of background noise.
AWeighting and Equal Loudness Curve: SPL measurements in terms of dB take no account of how sound pressure / power is distributed with respect to sound frequency as it represents a comparison or ratio between the quantity being measured (such as sound pressure or power) and a reference value (reference pressure or reference power). As the sound pressure level (SPL) reduces, human ear's ability o detect low or high frequency noise also diminishes. On the other hand, almost all noises are generally depedent of frequency and consist of a mix of frequencies. A complete represenation is a frequency spectrum (dB vs. f curve) where noise levels are plotted as a function of noise frequencies. For purposes like comparing two noise sources, specifying noise reduction level ... a single value is typically required. Aweighted sound pressure levels in terms of dB(A) are widely used for this purpose.
The offset values mentioned in formula below (2.0, 0.17 and 0.06) are to ensure normalization to 0 [dB] at 1000 [Hz] and calculated respectively as 20·log_{10}[Ra(1000 Hz)], 20·log_{10}[Rb(1000 Hz)] and 20·log_{10}[Rc(1000 Hz)]. The Aweighting curve describes decibels of 'attenuation' or 'gain' at every frequency over which sound pressure level is measured. Few key features of Aweighting curve are: At 1 [kHz], the gain or attenuation is '0' and this freqeuncy value is reference for all other cases.
 Below 1 [kHz], the sound levels are attenuated (diminished) with respect to input levels.
 There is attenuation of 20 [dB] at frequencies 100 [Hz] and 40 [kHz].
 In the range of 1 ~ 6 [kHz], there is a gain of few decibels. The maximum gain is 1.27 [dB] at 2500 [Hz].
 For frequencies > 6 kHz the sound levels are attenuated.
 As gain (that too small relative to attenuation) is limited to a narrow range of frequency, dB(A) value will tend to be < [dB] value in most of the cases but this is not universally true.
Weighting comes from FletcherMunson Curves where 'Aweighting' is 40 Phon equal loudness contour 'Bweighting' is 70 Phon equal loudness contour and 'Cweighting' is 100 Phon equal loudness contour.
Steps to calculate Aweighting:
 Adjust the given sound level by the Aweighted adjustments as per graph or formula described above. Note that the subtraction or addition of two decibel values for this particular purposes need to be algebraic.
 Add the sound levels as per the actual addition formula to calculate the overall sound level.
 Example calculation: noise levels below 63 [Hz] and above 8 [kHz] are ignored.
f [Hz]  63  125  250  500  1000  2000  4000  8000 
Lp: measured [dB]  69  71  73  77  79  78  75  72 
Lp: total  84.48 [dB] 
Aweighting value: δ_{A}  26.2  16.1  8.6  3.2  0.0  1.2  1.0  1.1 
Lp' [dB] = Lp + δ_{A}  42.8  54.9  64.4  73.8  79.0  79.2  76.0  70.9 
Lp(A): Aweighted  83.83 [dB] 
Noise Type, Acoustic Impedance and Directivity
 Broadband Noise  This is also called wideband noise  a type of noise whose energy is distributed over a wide range of frequencies of pressure fluctuations (or audible range).
 Tonal Noise  Tone refers to a specified frequency. Tonal noise refers to sound wave forms that occur at a particular frequency such as RPM of a rotating wheel  also known as narrowband noise. Tonal noise are discrete and occur only at certain frequencies that is the frequency of a tonal noise remains constant for a range of mass flow rates or input flow velocity scale. Sometimes it is also called humming or whining or ticking noise.
In the case of tonal stimulation at varied stimulus frequency (f), the ratio of the sound pressure p(f) to the volume velocity Q(f) defines the acoustic behaviour of a system. In a linear system, a sinusoidal excitation produces sinusoidal response which may lead, lag or be in the phase of the excitation. Similarly, the responses produced by a tonal stimulus can differ from the input in terms of the magnitude which can be quantify by a rootmeansquare (RMS) value, and phase (which describes the relative timing of the sinusoidal responses  leading or lagging).
 Blade passing frequency: BPF is one of the method to estimated expected tonal frequency component for rotating machines like fas, turbines, blowers, turbochargers, camshaft drive for IV/EV... It is calculated as N/60×n_{B} where n_{B} is the number of blades and N is RPM. This is frequency at which a longitudinal plane is cut through the rotating blades.
 Combustion noise: This is applicable to internal combustion engines where the combustion is not a continuous process and valves open/close at frequency related to operating RPM of the engine (crankshaft rotation). ν = engine crankshaft rotation speed in Hz = N/60 where N = RPM, combustion frequency for a 4 stroke engine = ν/2.
 Tappet noise: A tappet or camfollower is the part that is made to move vertically by the action of the rotating cam which subsequently causes valves to open. The tappet noise refers to a characteristic hammering or "taktak" noise due to clearance between inlet/exhaust valves and corresponding driving or rocker arms.
 Acoustic Impedance  Analogous to the flow of electric current, acoustic impedance of a material it the resistance of the passage of sound waves. Impedenace is a type of resistance (as we learnt for capacitors and inductors whose resistance is function of frequency of electric current following through them and complicated by the fact that current and voltage may not be in phase). Similarly, acoustic impedance is defined as ratio of acoustic pressure p to acoustic volume flow rate Q. In particular the impedance magnitude Z(f) is defined by the ratio of the magnitudes of the sound pressure p(f) and the volume velocity Q(f).
 Z = p / Q [Pa.s/m^{3}] ≡ [kg/m^{4}/s] or acoustic Ω. Similar to electric impedance, the flow and pressure may not be in phase and hence like electrical systems, complex numbers method needs to be used to handle such impedances where the real part represents the inphase component and the imaginary part the outofphase component. Imaginary part of the acoustic impedance is called acoustic reactance and the real part of the acoustic impedance is called acoustic resistance.
 The acoustic impedance describes the sound pressure [Pa] needed to produce a unit measure of volume velocity [m^{3}/s] in the tube. For an infinitely long pipe of cross sectional area A which is filled with a medium of density ρ and at temperature T at which speed of sound c = (γ×R×T)^{0.5}, the acoustic impedance is [ρ×c/A].
 The acoustic impedance of a material determines fraction of sound power that will be transmitted and reflected when the wave encounters an interface created by two different materials. The larger the difference in acoustic impedance between two materials, the smaller the amount of transmitted energy will be. Similarly, the acoustic impedance of a sound attenuation device determines fraction of sound power that will be transmitted and reflected when the wave passes through it.
 Characteristic acoustic impedance [Pa.s/m] is defined as product of density of the medium and speed of sound in it.
 Specific acoustic impedance [Pa.s/m] in a sound field is the ratio of sound pressure and particle velocity [p' / u'] which equals [ρ·c]. This is analogous to Ohm's law where electrical resistance R = ΔV/I, I = current and ΔV = voltage difference.
 Analogous to capacitative and inductive resistance in circuits carrying alternating current, acoustic inertance (where pressure leads the velocity) andacoustic compliance (where pressure lags the velocity) are defined.
 Acoustic Admittance  The acoustic impedance and admittance are related quantities that depend on the ratio of sound pressure and sound flow (the volume velocity of sound). As electrical conductance is reciprocal to electrical resistance (impedance), the acoustic admittance is reciprocal to acoustic impedance: Y(f) = 1/Z(f). The units of acoustic admittance are [acoustic siemens], where 1 [acoustic siemen] = 1 [m^{3}/Pa.s]. Admittance and impedance are often collectively refereed as 'immittance'.
 Transmission loss  It is defined as the difference between the sound power incident at the entry of the sound attenuation device such as filter or muffler and the sound power transmitted after the device. Thus: TL = 10 × log(W_{IN}/W_{OUT}) = 20 × log(p_{IN}/p_{OUT}). Transmission loss is measured by mounting a test wall or panel between two reverberation rooms and measuring the sound pressure levels on each side as described below.
Note that the pressure fields for TL measurements or simulation case an be generated without a background mean flow  just by using a noise source at one end also known as time harmonic acoustic perturbations.
 Directivity Factor, Q  This is a measure of variation of angular dependence of sound intensity far away from source (that is in a region free of nearfield effects). For a spherical waves, the average sound intensity at a distance r is given by I_{AVG}(r) = W/4*πr^{2}. Directivity factor Q is defined as Q_{θ} = I_{θ}(θ, r) / I_{AVG}(r). Directivity index is defined as DI = 10 × log_{10}[Q_{θ]}. This evaluation parameter is important for speakers where sound source need to cover a large angle without significant drop in SPL.
 Power Spectral Density  PSD is a representation of the dependence of energy on frequency (or wavenumber). Most often PSD is not a smooth function of time or space and respresented in terms of the Fourier Coefficient. Hence, it is also viewed as a distribution of squared Fourier coefficients as a function of wavenumber (or frequency).
AeroAcoustic Calculations
Excerpts from "Surface integral methods in computational aeroacoustics — from the (CFD) nearfield to the (Acoustic) farfield" by Anastasios S. Lyrintzis, School of Aeronautics and Astronautics, Purdue University:
"Aerodynamically generated sound is governed by a nonlinear process. One class of problems is turbulence generated noise (e.g. jet noise). An accurate turbulence model is usually needed in this case. A second class of problems involves impulsive noise due to moving surfaces (e.g. helicopter rotor noise, propeller noise, fan noise etc.). In these cases an Euler/NavierStokes model or even a full potential model is adequate, because turbulence is not important. Furthermore, because the acoustic fluctuations are usually quite small (about three orders of magnitude less than the flow fluctuations), the use of nonlinear equations (whether NavierStokes or Euler) could result in errors. One usually has no choice but to separate the computation into two domains, one describing the nonlinear generation of sound, the other describing the propagation of sound."
Similarly, excerpts from "Implementation of Acoustical Analogies in Openfoam and CALFEM" by Johan Nilsson: There are several challenges to overcome when attempting to compute aeroacoustical sound. Two of these are the large energy and length scale disparity. The energy in the flow is much greater than the energy in the acoustical waves. To quote Crighton, "If all the acoustic energy radiated in the 45second takeoff of a large jet transport were recovered, that energy would be about enough to fry one egg!". Compare that to the energy in the flow which is used to lift the whole large jet aircraft to get a sense of the energy disparity.
Acoustic sources
James Lighthill in the 1950s using free flow assumptions (no reflections in solid surfaces) developed a wave equation including a source term caused by the instantaneous velocity and pressure fluctuations in the flow. He used the the same set of governing equations which describe flow field for fluids. The way to solve the acoustic wave equation is by integration of the sources using Green functions. This analogy holds true for cases like noise generated by jet engines in aeroplanes but may not be appropriate for closed volume like automotive passenger cabins.
Turbulent flow is considered in the source region and used as a excitation to the acoustic propagation. The sources are categorised as quadrupoles (within the field such as in turbulent wakes  derived based on time history of velocity field), dipoles (on stationary surfaces with flowinduced fluctuating pressure  derived based on time history of the static pressure on the surfaces), monopoles (massflow fluctuation across a boundary such as wall or when mass is added at an unsteady rate: in turbomachines it is called thickness noise due displacement of fluid caused by motion blades having finite thickness) and rotating dipoles (fluctuating pressure on moving surfaces such as blades of a fan and turbines, also known as loading noise  derived based on time history of the static pressure on the surfaces). A monopole source [1/s^{2}] can be thought as point source which has a equal intensity in all directions  a radially oscillating sphere. A dipole source [N/m^{3}] is stronger in two opposite directions and can be visualized as translational oscillating sphere.
Direct Method: Any acoustic calculation approach requires at least 2 steps: generation and propagation of sound. In direct method, both generation and propagation of sound waves are directly computed by solving the fluid dynamics equations modified to capture sound waves. This method requires very fine meshes encompassing receivers as well as acoustically nonreflecting boundary conditions.
Hybrid (CFD/CAA) Approach: This is called so because noise predictions are a twostep approach due to the large difference in requirements on the flow field and acoustic propagation. Ffowcs WilliamsHawkings formulation (bounded flow) and Lighthill's analogy (unbounded field) are examples of hybrid approach.
 Step  1 : Calculate the source of noise  this is typically achieved by a 3D simulation in a CFD program, either using RANS or LES approach. The Acoustic Analogy approach introduced by James Lighthill is used to divide computational domain into a non–linear source region and a wave propagation region.
CFD simulation can be either of the followings. ANSYS FLUENT recommends LES turbulence model because it resolves all eddies with scales larger than grid scale and hence wideband aeroacoustic noise can be predicted correctly using LES model. Similarly, to calculate presssure fluctuations on the wall, PRESTO! is recommended for spatial discretization of pressure as it is more accurate scheme for interpolating face pressure values from cell pressures. FLUENT also recommends "Bounded Central Differencing" for momentum in case of LES on unstructured mesh.
 Incompressible, steady state
 Incompressible, transient
 Compressible, transient
 Step  2 : Extract the noise sources using Lighthill's analogy for low Mach number flows and Mohring analogy for high Mach number flows. The fluid flow domain uses a finer mesh than the acoustical mesh. The later may include any objects that is required to see the effect of their shape and size in the acoustical domain. Also, the fluid and the acoustical domains do not need to encomapss the same space and acoustic domain is usually always bigger than the fluid domain.
 Step  3 : Use CAA tool such as SYSNOISE, ACTRAN, NUMECA to solve for noise propagation and SPL results.
A process flow based on user manual and tutorial examples in ANSYS FLUENT can be summarized as per steps described below:
Broadband Noise Source Model: In turbulent flows with no rotating components like fans and blowers, sound energy is a continuous distribution over a broad range of frequencies. In other words, no distinct spikes in sound levels are present at few discrete frequencies only. In such cases turbulence quantities [velocity field, turbulent kinetic energy  k and the dissipation rate  ε] computable from RANS equations can be utilized to estimate source of broadband noise. However being a nearfield approach, the model does not predict the sound at receivers far away from the flow field. This model can be employed to extract diagnostic information on the noise source to check portion of the flow primarily responsible for the noise generation.
This method available in ANSYS FLUENT quantifies the local contribution (per unit surface area or volume) to the total acoustic power generated by the flow and doed not require transient solutions to the flow field.
The FWH method requires a transient simulation with timesteps lower than highest frequency of sound to be resolved. For time step δt, the highest frequency computed in acoustic analysis is f = 0.5/δt .
Ways of Attenuation of Sound Levels
Sound Attenuation Devices
There are two ways to attenuate the level of sound: reflective systems  here the incident sound is scattered and canceled by destructive interference and dissipative system  here the incident sound energy is absorbed and hence has to be converted into the heat. Expansioncontraction chambers, resonators and HerschelQuincke tube fall under the category of reflective systems.
Automotive exhaust mufflers use perforated tubes. The STEP file and FreeCAD file of this geometry can be downloaded for mesh generation and simulation.
There are two basic types of mufflers or silencers:
 Reflective or reactive mufflers  the operating principle is to reflect acoustic waves by abrupt area expansions or changes of impedance. Such mufflers are recommended for the low frequency range where only plane waves can propagate in the system.
 Dissipative mufflers  these type of mufflers dissipate acoustic energy into heat through viscous losses either in fibrous materials and/or flowrelated hydraulic losses such perforated pipes shown above. This type of muffler is recommended for higher frequency ranges.
The plug flow muffler shown above works at all frequencies by dissipating sound energy when waves formed due vortices at perforation interfere with the waves of the sound field. Thus, the porosity (fraction of void space) and flow velocity through the holes will decide the performance of such mufflers. The perforates become more effective in automotive exhaust systems where flow is pulsating (as opposed to a constant mean flow).
 Most of the systems have noise components varying from low to high frequency range and hence no single type of muffler serves the purpose. Hence, dissipative mufflers can be made of hybrid type where flow losses can be used to work also at low frequencies. For example, automotive exhaust mufflers are of hybrid type using a combination of reflective and dissipative muffler elements where reflective part removes lowfrequency engine noise generated by combustion and valve motion while the dissipative part is designed to take care of higherfrequency noise.
Transmission and Insertion Loss for a Muffler: As described in the schematic below, the gas pressure in the duct is the result of superposition of waves coming from the source and reflections from the muffler cavity, it is impossible to directly measure the incident components of the pressure pulsation into the muffler. At the same time, all acoustic frequencies of interest (typically in the range 20 ~ 2000 Hz) must be excited which requires a need of largeband acoustic source during the test. Anechoic termination is also known as reflection free termination. TL is defined as difference in sound power between the incident wave entering and the transmitted wave exiting the muffler when termination is anechoic (no reflecting waves present). However, inline with other definition, it is calculated as logarithm of the ratio of sound powers. Attenuation is the difference in incident and transmitted sound powers through the muffler but the termination need not be anechoic. Thus, attenuation is depedent on outlet pipe length. TL is independent of the inlet and outlet pipe length and depends only on the geometry of muffler.
Insertion loss = SPL_{1}  SPL_{2} which is the difference between the radiated sound power with a reference exhaust system (without any muffler) and the radiated sound power with a muffler system. IL is dependent on both design the muffler as well as the impedances of source and the radiation. Similarly, pressure loss (PL) is an important performance parameter which depends only on the shape and size of the muffler and is defined as difference between the (static) pressure drop across the muffler. The sound (or acoustic) pressure level difference across the muffler is termed as Noise Reduction (NR).
Design parameters  Insertion loss  Transmission loss 
Dependent on geometry?  Y  Y 
Influence by source impedance?  Y  N 
Dependent on radiation impedance?  Y  N 
Anechoic termination required?  N  Y 
Calculation  Difficult  Easy 
Measurement  Easy  Difficult 
Noise Control Methods
Analogous to HVAC (where cooling effectiveness is achieved by reducing heat load and maintaining uniform cooling), most of the noise control techniques use a combination of 'absorption', 'reflection' and "transmission or diversion of sound" energy through airborne and stuctureborne paths. Such arrangements and walls are evaluated in terms of "transmission coefficients" and "sound transmission loss" values.  Transmission coefficient is a frequencydependent physical property of noiseattenuating materials define as τ = intensity of transmitted sound / intensity of incident sound energy. It also denotes ability to reflect or block sound energy.
 An ideal reflecting wall or material has transmission coefficient = 0 that is STL = ∞, the transmission coefficient of an opening is 1.0 that is STL =0.
 Sound transmission loss is logarithm of inverse of transmission coefficient: STL = 10 × log_{10}1/τ
 Sound absorbing materials allow sound waves to enter into its skin and dissipate the energy through viscous dissipation. Such materials are generally porous such as foam and laminates.
Steps to control noise: Determine noise sources using intensity meter or sound camera (an array of microphones) and identify airborne, liquidborne (cavitation) or structurborne (vibration) sources. Note that the measurement of noise radiating from a particular machine can be accurately measured only when the background noise is lower by more than 3 [dB].
 Noise control through air flow
 Reduce pressure drop or flow rates at outlet grilles / louvres
 Avoid barriers or obstacles in the flow
 Improve geometry for smoother flow: elbow bend instead of sharp knee bend, use of screens and separators in bends
 Noise control for structureborne noise
 Check and correct imbalances in rotating parts
 Use sealing rubbers near mounting of protective covers such as belt or chain drives
 Use flexible hoses to reduce trasnmission of vibrations
 For solid particle handling devices (sitteres, mixers, blenders): increase impact time (reduce speed) and decrease impact/collision speed, minimize impact mass by using smaller particles and increase mass of housing of a mixing tank.
 Select appropriate bearing keeping in ming vibration levels, proper lubrication also dampens the vibration.
 Use acoustic (perforated) screens  single or multilayered with a layer of sound energy absorbing (damping) material attached towards sound energy source.
 Stiffen the panels by crossribbing
Application of AeroAcoustics  CAA
 Automotive Exhaust Silencer (UK English) or Mufflers (US English): intake and exhaust noise can be due to pipe resonance, shell vibration and turbulent fluctuations.
 Cabin Noise from HVAC module: tonal component due to blowers, broadband component due to turbulent flows and shell vibrations
 Instrument panel vibration, ORVM (Outside Rear View Mirror) and APillar noise
 Rotating devices: Compressor and Fan noise, high frequency noise in turbochargers  whistling, structureborne noise
 Vacuum cleaners
 Building airconditioning Systems: duct noise, sound at exits from grilles (louvres), design of concert halls and conference rooms where an even distribution of sound pressure is desired
Numerical Tools for Aero and Vibroacoustics
 LMS SYSNOISE (now part of Siemens): Based on Boundary Element method (BEM) which is a well established numerical method for solving acoustic radiation problems in unbounded domains. It can be coupled with many CFD softwares such as FLUENT.
 ACTRAN by MSC Software Corporation: For example: process flow chart for a vibroacoustic simulatio in ACTRAN is: Multibodydynamics (MBD  ADMAS) → NASTRAN → ACTRAN for structureborne noise. An example screenshot for aeroacoustic simulation for ORVM noise in cabin of a car is
 ANSYS FLUENT, Mechanical and Multiphysics, Acoustics ACT Extension: FLUENT is used primarily for acoustic field generation using broadband noise and FWH analogy. Mechanical and Multiphysics modules can be used for harmonic response analysis and aerovibroacoustics. FLUID29 and FLUID30 are used to model the fluid for 2D and 3D cases respectively. FLUID129 and FLUID130 are companion elements to FLUID29/30 for modeling farfield or infinite envelope (absorbs incoming pressure waves) around the FLUID29/30 elements.
Broadband noise source calculated in ANSYS FLUENT for flow over a bluff body is shown below.
ANSYS FLUENT also has "Jet Noise Source Model" based on work of Goldstein to account for anisotropy of turbulence in axisymmetric turbulent jets.
 OpenFOAM: The official release has postprocessing utility called 'noise' for noise analysis of pressure data. The nose postprocessing utility requires pressure data in a file which can be generated either by probe function at few select locations or using 'sample' utility. The later can be used to write or export pressure fluctuation data over a patch or surface using sampleDict utility. This is required to calculated noise from a surface source.
OpenFOAMdev has a library called libAcoustics farfield noise computation which can be accessed from here.
 GTPower from Gamma Technologies and Wave from Ricardo: Primarily focused on automotive systems for noise calculations in engines, HVAC systems and gear chattering in gearboxes.
 NUMECA: General purpose CFD software with acoustics capabilities named FINE^{TM}/Acoustics. It can calculate tonal noise, broadband noise sources and its propagation with Boudnary Element Method (BEM) / Finite Element Method (FEM). It also has Ffowcs WilliamsHawkings (FWH) solver for radiation analysis with capability to import unsteady flow field data from most commercial CFD tool. It computes vibroacoustic noise through fluidstructure coupling by importing structural modes from NASTRAN, ABAQUS and ANSYS.
 COMSOL: It is a general purpose multiphysics solver which includes acoustics capabilities. COMSOL has two interfaces, transient interface and frequency domain interface for modeling the sound waves in the time domain and frequency domain respectively. Transient interface computes the pressure variation by modeling the propagation of acoustic waves in fluids at quiescent background conditions. Frequency Domain interface computes the pressure variation for propagation of acoustic waves in fluids at quiescent background conditions. It is suited for all frequencydomain simulations with harmonic variations of the pressure field.
Both ANSYS Mechanical and COMSOL uses PML (perfectly matched layer) boundary type which emulates a nonreflecting boundary condition independent of the shape and frequency of the incident wave front.
It also has a FSI module for (vibroacoustics) acousticstructure interaction to model multiphysics phenomenon where the acoustic pressure causes a fluid pressure on the solid surface and the structural acceleration (surface velocity) acts on the fluid domain as a normal acceleration across the fluidsolid boundary.
 STARCCM+ The CFD program has 3 models for broadband noise calculations: Curle surface acoustic power model (dipole source), Proudman surface acoustic power (acoustic quadrupole sources such as on symmetry planes) and Goldstein broadband noise source model (accounts for anisotropy of turbulence such as axisymmetric turbulent jets). Proudman model calculates acoustic power per unit volume as generated by turbulence structures such as wake and flow sepration regions around a cylinder. As per the tutorial document: "The Broadband Noise Source models are intended to be used as tools for qualifying a volume mesh before a transient acoustic analysis. Curle acoustic power model assumes an isotropic turbulent flow field and computes the sound that dipole sources generates. This model therefore represents the fluctuating surface pressure with which the solid boundary acts on the flow  noise that the turbulent boundary layer emits over a surface at low Mach number."
STARCCM+ has FWH model that can be used with transient flow field and requires FWH surfaces and locations of the FWH receivers (points at specific distances from the noise sources). Then FFT (Fast Fourier Transforms) is used to determine sound pressure as a function of frequency at each of the receiver(s). This is described as "Point Time Fourier Transform (G[p])". Similar conversion from time domain to frequency domain can be performed for a line and surface using "Line Time Fourier Transform (G[h[l]])" and "Surface Time Fourier Transform (G[h[s]])" respectively.
VibroAeroacoustics
A generic approach is described for an automotive engine. Note that this method is based on a weak fluidstructure coupling that is one way coupling in structural vibrition and acoustic field. Thus, this is not a typical FSI (FluidStrucure Interaction) approach.
Calculation step  Method, inputs 
Source: Multibody simulation for estimating the forces acting on the engine during operation   Combustion pressure
 Piston slap force (side thrust load)
 Bearing loads
 Valvetrain impact loads

Vibroacoustics transfer: FEA powertrain system model for cylinder block / head, oil pan, cam or tappet cover: Structural response of the engine as a function RPM  FEA model with loads calculated in previous step: to compute the structural normal velocity on a vibrating surface, structural FEA is used for modal and mode superposition: linear combination of the mode shapes of the body to get deformation shapes of the structure at different excitation frequencies for all RPM loads. 
Acoustic modeling: Calculation of noise radiated from engie  ACTRAN, SYSNOISE  The vibroacoustic relationship between engine vibrations and the acoustic pressure field is performed in a separate program using BEM or ATV (Acoustic Transfer Vectors) concept. 
A process flow chart to carry out aerovibroacoustic simulation using ANSYS FLUENT, Mechanical and Acoustic ACT extension is described here:
Special Boundary Conditions for CAA
 Nonreflecting boundary at outlet in ANSYS FLUENT and waveTransmissive boundary condition in OpenFOAM: It does not let the wave reflect at the outlet. In OpenFOAM, it provides a wave transmissive outflow condition by solving DDt(psi, U) = 0 at the boundary where psi is compressibility. waveTransmisive b.c. not perfectly non reflecting and small acoustic waves are reflected to the inner domain.
 Anechoic chamber: it is a representation of walls which absorbs sound energy incident on it.
 Radiation boundary: A face or zone which reflects sound energy as per specified coefficient.
TMM (Transfer Matrix Method)
Also known as transmission matrix or fourpole parameter presentation, TMM is a generalpurpose method used to solve problems quickly in electromagentics, mechanical vibration as well as aeroacoustics. This method is used when a system can be condensed into a set or network of subsystems that interact only with adjacent subsystems (applicable to linear and passive components). E.g intake flow path of an automotive engine where entry, air cleaner, duct, intake plenum and manifold runner are subsystems of overall intake system. These subsystems are mathematically connected using transfer matrices which can be obtained either analytically or using numerical methods. This is analogous to a flow or thermal network where a system is broken into a network of nodes and elements with "lumped properties".
One dimensional, plane wave propagation is assumed and hence method is accurate only below cuton frequency (defined as the frequency below which only plane waves propagate inside a duct). For a duct of diamter D, f_{CUTON} = 1.8412 × c/πD [ref: A. D. Pierce: Acoustics – An Introduction to its Physical Principles and Applications]. Similarly, cutoff frequency can be defined as the frequency above which no plane waves propagate inside a duct. For elliptical ducts, f_{CUTOFF} = β × c/2πa where 2a = major diameter, 2b = minor diameter, β = f(e), e = sqrt(1b^{2}/a^{2})
 Step1: Condense a complex systems into simpler subsystems arranged into series, espcially closer to standard shapes for which fourpole parameters are readily available.
 Step2: Calculate transmission loss of each subsystem: example for a expansion chamber is demonstrated below. Here, p is acoustic pressure and q is volume velocity. Note that the transfer matrix is a frequencydependent property of a system, is usually symmetric and reciprocity rule requires the determinant of matrix, T_{11}·T_{22}  T_{12}·T_{21} = 1. T_{ij} are also called fourpole constants.
 Step3: Overall TL of the system is product of TL of each subsystem in series.
Noise Level Regulations
NR: The permissible exposure limit (PEL) is a combination of sound pressure level and time of exposure. As per OSHA Hearing Conservation Amendment (HCA) 1981–1983, sound pressure level of 80 dBA has time limit of 32 hours whereas a level of 110 dBA has time limit of only 0.5 hours.
 Determination of sound power levels of noise sources — Guidelines for the use of basic standards: ISO 1680:2013
 Test code for the measurement of airborne noise emitted by rotating electrical machines: ISO 1680:2013
 Measurement of sound pressure level from service equipment in buildings — Engineering method: ISO 16032:2004
 Requirements for the performance and calibration of reference sound sources used for the determination of sound power levels: ISO 6926:2016
 Determination of sound power levels and sound energy levels of noise sources using sound pressure — Precision methods for reverberation test rooms: ISO 3741:2010
 Determination of sound power levels and sound energy levels of noise sources using sound pressure — Engineering methods for an essentially free field over a reflecting plane: ISO 3744:2010
 Determination of sound power levels and sound energy levels of noise sources using sound pressure — Precision methods for anechoic rooms and hemianechoic rooms: ISO 3745:2012
 Attenuation of sound during propagation outdoors — Part 1: Calculation of the absorption of sound by the atmosphere: ISO 96131:1993
 Determination of sound power levels of noise sources using sound intensity — Part 3: Precision method for measurement by scanning: ISO 96143
 Passby Noise: ISO 7216:1992
 Noise for agriculture and forestry wheeled tractors in motion fitted with elastic tyres or rubber tracks: governed by standard ISO 7216:2015. Test method specified by this standard require an acoustical environment which can only be obtained in an wide open space. Noise limits:
 Tractor with no load condition < 1.5 ton: 85 dB
 Tractor with no load condition ≥ 1.5 ton: 89 dB
 Construction equipments
 Noise level at operator's position at stationary condition is governed by ISO 6394:2008. For exterior locations, ISO 6393:2008 applieds.
 Noise level at operator's position under dynamic condition is governed by ISO 6396:2008. For exterior locations, ISO 6395:2008 applies.
NC Curve
Noise Criteria (NC) rates the indoor background noise levels in buildings and rooms. It is used to judge the appropriateness of the acoustic environment for various activities. The American Society of Heating, Refrigeration and Air Conditioning Engineers (ASHRAE) makes recommendations for the NC levels for various spaces such as Concert Halls, Conference room, Classrooms ...
Speech Interference Level (SIL)
Excessive background noise interferes with speech which causes annoyance and irritation. When verbal communication is disrupted, worker's efficiency can suffer and the likelihood of error due to miscommunication is increased. Aweighted sound level is the simplest method used to predict speech intelligibility. It is a plot to specify talkedtolistener distance for justreliable communications. It is common for a person to raise his voice when the background noise is high. Typically, a person raises the voice by 3 ~ 6 [dB] for every 10 [dB] increase in background noise level above 50 dBA.
Frequency bands for octave and 1/3octaves
B and CWighting Functions: While dBA is used to emphaise the fact that low frequency and high frequency noise pose less risk to humans, dBC is used for selection of hearing protection devices.
The frequency weighting correction value in tabulate form is as follows:
Centre  Frequency Correction [dB] 
Frequency [hz]  Aweighting  Cweighting  Dweighting 
10  70.4  14.3   
12.5  63.4  11.2   
16  56.7  8.5   
20  50.5  6.2   
25  44.7  4.4   
31.5  39.4  3.0   
40  34.6  2.0   
50  30.2  1.3  12.8 
63  26.2  0.8  10.9 
80  22.5  0.5  9.0 
100  19.1  0.3  7.2 
125  16.1  0.2  5.5 
160  13.4  0.1  4.0 
200  10.9  0.0  2.6 
250  8.6  0.0  1.6 
315  6.6  0.0  0.8 
400  4.8  0.0  0.4 
500  3.2  0.0  0.3 
630  1.9  0.0  0.5 
800  0.8  0.0  0.6 
1000  0.0  0.0  0.0 
1250  0.6  0.0  2.0 
1600  1.0  0.1  4.9 
2000  1.2  0.2  7.9 
2500  1.3  0.3  10.6 
3150  1.2  0.5  11.5 
4000  1.0  0.8  11.1 
5000  0.5  1.3  9.6 
6300  0.1  2.0  7.6 
8000  1.1  3.0  5.5 
10000  2.5  4.4  3.4 
12500  4.3  6.2  1.4 
16000  6.6  8.5   
20000  9.3  11.2   
Lighthill's Equation
Most of the derivations available in textbooks and literature on acoustics use Einstein's shorthand notation. Hence, a clear understanding of the convention is a strong prerequisite to understand each steps described below. The shorthand notation is based on concepts of "free index" and "dummy index".
 A free index appears once and only once within each additive term. The same letter must be used for the free index in every additive term.
 The first free index in a term corresponds to the row, and the second corresponds to the column. Thus, a vector (which has only one free index) is written as a column of three rows and a second order tensor is represented as 3 × 3 size matrix.
 A dummy index defines an index which does not appear in the final expression any more and implies a summation over the range of the index.
 A dummy index appears twice and exactly twice within an additive term of an expression.
 Number of freeindex on left hand and right hand sides of an equation should be equal.
The continuity or mass conservation equation is given by:
NavierStokes equation in vector form is expressed as:
The same equation in scalar or x, y and zcomponent are written as set of 3 equations:
XMomentum equation:
YMomentum equation:
ZMomentum equation:
In order to get a single scalar equation, the momentum equation also needs to be converted into a scalar equation. This is achieved by taking divergence of equation [ii]. Thus:
Introducing pressure and density fluctuations:
If the right hand side of Lighthill's equation is assumed to be known or estimated independent of the left hand side, this equation can be viewed as an inhomogeneous wave equation in an isotropic medium at rest. This also means that the sound propagation phenomena is not affected by sound generation. In flows with low Mach number, this assumption is useful and holds true and most aeroacoustic simulations are based on this simplification. Note that Lighthill's analogy does not account for presence of boundaries. This analogy is also classified as "densitybased analogy" as the primary variable of the governing equation is the density of the medium. Other variants of this equations have evolved such as FfowcsWilliams and Hawkings equation that includes surfaces in motion. Thus, FWH analogy is used in solid surface interactions that are directly involved in the generation of flow sound, for example those due to rotors of a helicopter, propellers driving aeroplanes and marine vessels and prime movers like fans, compressors and turbines.
Similarly, a special version of FWH analogy is the Curle’s analogy that takes into account presence of solid rigid boundaries where only the dipole source terms takes care of the sound scattering caused by the stationary surface.
^{[6]}In addition to analogies presented above, numerious other analogies have been developed: Phillips’ analogy and Lilley’s analogy are based on scaled logarithmic pressure Π = 1/γ·ln[p'/p_{0}]. These are also known as phibased analogy and have been used as a starting point for aeroacoustics calculations of jet engines.
Share on
References
 https://www.ncbi.nlm.nih.gov/pmc/articles/PMC4906304/
 COMSOL: Acoustic Module User Guide
 ANSYS FLUENT User Manual and Theory Guide
 Computational Aeroacoustics Methods with OpenFOAM v. 4.1: V. Korchagova, M. Kraposhin, S. Strizhak
 OpenRotor Noise Assessment with CFD/CAA Chaining: A. Chelius, T. Le Garrec, D.C. Mincu, Onera  The French Aerospace Lab, F92322, Chatillon, France
 Foundations of acoustic analogies: Seppo Uosukainen, VTT Publications 757
 "AWeighting": Is it the metric you think it is?: McMinn, Terrance  Curtin University of Technology
 https://uiowa.edu/voiceacademy/malefemalevoices
 The frequency range of the voice fundamental in the speech of male and female adults: Hartmut Traunmüller and Anders Eriksson, Institutionen för lingvistik, Stockholms universitet, S106 91 Stockholm, Sweden.
 Helmholtz resonator with extended neck: Ahmet Selamet and Iljae Lee, Department of Mechanical Engineering and The Center for Automotive Research, The Ohio State University, Columbus, Ohio 432101107
 Noise: Patrick N. Breysse and Peter S.J. Lees, Johns Hopkins University
 http://www.who.int/occupational_health/publications/noise1.pdf
 Predicting wind noise and other flowinduced noise sources with FLUENT/SYSNOISE coupling: Z. El Hachemi, M. Tournour, F. Kelecy, Y. DAI