- •Preface
- •Objectives of the Book
- •Style
- •Prerequisites
- •The Big Picture
- •Contents
- •1.1 Review of Complex Numbers
- •1.2 Complex Numbers in Polar Form
- •1.3 Four Equivalent Forms to Represent Harmonic Waves
- •1.4 Mathematical Identity
- •1.5 Derivation of Four Equivalent Forms
- •1.5.1 Obtain Form 2 from Form 1
- •1.5.2 Obtain Form 3 from Form 2
- •1.5.3 Obtain Form 4 from Form 3
- •1.6 Visualization and Numerical Validation of Form 1 and Form 2
- •1.8 Homework Exercises
- •1.9 References of Trigonometric Identities
- •1.9.1 Trigonometric Identities of a Single Angle
- •1.9.2 Trigonometric Identities of Two Angles
- •1.10 A MATLAB Code for Visualization of Form 1 and Form 2
- •2.2 Equation of Continuity
- •2.3 Equation of State
- •2.3.1 Energy Increase due to Work Done
- •2.3.2 Pressure due to Colliding of Gases
- •2.3.3 Derivation of Equation of State
- •2.4 Derivation of Acoustic Wave Equation
- •2.5 Formulas for the Speed of Sound
- •2.5.1 Formula Using Pressure
- •2.5.2 Formula Using Bulk Modulus
- •2.5.3 Formula Using Temperature
- •2.5.4 Formula Using Colliding Speed
- •2.6 Homework Exercises
- •3.1 Review of Partial Differential Equations
- •3.1.1 Complex Solutions of a Partial Differential Equation
- •3.1.2 Trigonometric Solutions of a Partial Differential Equation
- •3.2 Four Basic Complex Solutions
- •3.3 Four Basic Traveling Waves
- •3.4 Four Basic Standing Waves
- •3.5 Conversion Between Traveling and Standing Waves
- •3.6 Wavenumber, Angular Frequency, and Wave Speed
- •3.7 Visualization of Acoustic Waves
- •3.7.1 Plotting Traveling Wave
- •3.7.2 Plotting Standing Wave
- •3.8 Homework Exercises
- •4.2 RMS Pressure
- •4.2.1 RMS Pressure of BTW
- •4.2.2 RMS Pressure of BSW
- •4.3 Acoustic Intensity
- •4.3.1 Acoustic Intensity of BTW
- •4.3.2 Acoustic Intensity of BSW
- •4.5.1 Issues with Real Impedance
- •4.6 Computer Program
- •4.7 Homework Exercises
- •4.8 References
- •4.8.1 Derivatives of Trigonometric and Complex Exponential Functions
- •4.8.2 Trigonometric Integrals
- •5.1 Spherical Coordinate System
- •5.2 Wave Equation in Spherical Coordinate System
- •5.3 Pressure Solutions of Wave Equation in Spherical Coordinate System
- •5.4 Flow Velocity
- •5.4.1 Flow Velocity in Real Format
- •5.4.2 Flow Velocity in Complex Format
- •5.5 RMS Pressure and Acoustic Intensity
- •5.7 Homework Exercises
- •6.1 Review of Pressure and Velocity Formulas for Spherical Waves
- •6.2 Acoustic Waves from a Pulsating Sphere
- •6.3 Acoustic Waves from a Small Pulsating Sphere
- •6.4 Acoustic Waves from a Point Source
- •6.4.1 Point Sources Formulated with Source Strength
- •6.4.2 Flow Rate as Source Strength
- •6.5 Acoustic Intensity and Sound Power
- •6.6 Computer Program
- •6.7 Project
- •6.8 Objective
- •6.9 Homework Exercises
- •7.1 1D Standing Waves Between Two Walls
- •7.2 Natural Frequencies and Mode Shapes in a Pipe
- •7.3 2D Boundary Conditions Between Four Walls
- •7.3.1 2D Standing Wave Solutions of the Wave Equation
- •7.3.2 2D Nature Frequencies Between Four Walls
- •7.3.3 2D Mode Shapes Between Four Walls
- •7.4 3D Boundary Conditions of Rectangular Cavities
- •7.4.1 3D Standing Wave Solutions of the Wave Equation
- •7.4.2 3D Natural Frequencies and Mode Shapes
- •7.5 Homework Exercises
- •8.1 2D Traveling Wave Solutions
- •8.1.2 Wavenumber Vectors in 2D Traveling Wave Solutions
- •8.2 Wavenumber Vectors in Resonant Cavities
- •8.3 Traveling Waves in Resonant Cavities
- •8.4 Wavenumber Vectors in Acoustic Waveguides
- •8.5 Traveling Waves in Acoustic Waveguides
- •8.6 Homework Exercises
- •9.1 Decibel Scale
- •9.1.1 Review of Logarithm Rules
- •9.1.2 Levels and Decibel Scale
- •9.1.3 Decibel Arithmetic
- •9.2 Sound Pressure Levels
- •9.2.2 Sound Power Levels and Decibel Scale
- •9.2.3 Sound Pressure Levels and Decibel Scale
- •9.2.4 Sound Pressure Levels Calculated in Time Domain
- •9.2.5 Sound Pressure Level Calculated in Frequency Domain
- •9.3 Octave Bands
- •9.3.1 Center Frequencies and Upper and Lower Bounds of Octave Bands
- •9.3.2 Lower and Upper Bounds of Octave Band and 1/3 Octave Band
- •9.3.3 Preferred Speech Interference Level (PSIL)
- •9.4 Weighted Sound Pressure Level
- •9.4.1 Logarithm of Weighting
- •9.5 Homework Exercises
- •10.1 Sound Power, Acoustic Intensity, and Energy Density
- •10.2.2 Room Constant
- •10.2.3 Reverberation Time
- •10.3 Room Acoustics
- •10.3.1 Energy Density due to an Acoustic Source
- •10.3.2 Sound Pressure Level due to an Acoustic Source
- •10.4 Transmission Loss due to Acoustical Partitions
- •10.4.2 Transmission Loss (TL)
- •10.5 Noise Reduction due to Acoustical Partitions
- •10.5.1 Energy Density due to a Partition Wall
- •10.5.2 Sound Pressure Level due to a Partition Wall
- •10.5.3 Noise Reduction (NR)
- •10.6 Homework Exercises
- •11.1 Complex Amplitude of Pressure and Acoustic Impedance
- •11.1.3 Transfer Pressure
- •11.2 Complex Acoustic Impedance
- •11.3 Balancing Pressure and Conservation of Mass
- •11.4 Transformation of Pressures
- •11.5 Transformation of Acoustic Impedance
- •11.7 Numerical Method for Molding of Pipelines
- •11.8 Computer Program
- •11.9 Project
- •11.10 Homework Exercises
- •12.1.1 Equivalent Acoustic Impedance of a One-to-Two Pipe
- •12.2 Power Transmission of a One-to-Two Pipe
- •12.3 Low-Pass Filters
- •12.4 High-Pass Filters
- •12.5 Band-Stop Resonator
- •12.6 Numerical Method for Modeling of Pipelines with Side Branches
- •12.7 Project
- •12.8 Homework Exercises
- •Nomenclature
- •Appendices
- •Appendix 1: Discrete Fourier Transform
- •Discrete Fourier Transform
- •Fourier Series for Periodical Time Function
- •Formulas of Discrete Fourier Series
- •Appendix 2: Power Spectral Density
- •Power Spectral Density
- •Accumulated Sound Pressure Square
- •Sound Pressure Level in Each Band
- •References
- •Index
6.3 Acoustic Waves from a Small Pulsating Sphere |
141 |
|||
at r ¼ a : |
uða, tÞ ¼ |
A |
ka þ θ ϕaÞ ¼ Ua cos ðωt þ θoÞ |
|
|
cos ðωt |
|||
aρoc cos ðϕaÞ |
where the subscript a indicates the value of the associated constant at r ¼ a and then ϕa will be:
ϕa ¼ tan 1 1 ka
The unknown constant A can be obtained by comparing the magnitude of the velocity evaluated at r ¼ a to the magnitude of the given velocity at r ¼ a as:
A
aρoc cos ðϕaÞ ¼ Ua
! A ¼ aρoc cos ðϕaÞUa
This concludes the validation of formulas for velocity and sound pressure from a pulsating spherical source.
6.3Acoustic Waves from a Small Pulsating Sphere
This section will study velocity and sound pressure of spherical waves radiated from a small spherical source. Formulas for velocity and sound pressure can be simplified when the wave is radiating from a small sphere (a 1). Such a source is called a small spherical source. Since a is small, ka and ϕa become:
ω
ka ¼ c a 0
!ϕa ¼ tan 1 ka1 π2
ka
! cos ðϕaÞ ¼ q ka
1 þ ðkaÞ2
Substitute the above approximated values into the following equations:
a
pðr, tÞ ¼ Uaρoc cos ðϕaÞ r cos ðωt kðr aÞ þ θo þ ϕaÞ
uðr, tÞ ¼ Ua cos ðϕaÞ a cos ðωt kðr aÞ þ θo þ ϕa ϕÞ cos ðϕÞ r
The pressure and velocity of the small spherical source are simplified into:
142 |
|
|
|
|
|
6 Acoustic Waves from Spherical Sources |
|
|
a |
|
|
|
π |
|
|
pðr, tÞ Uaρocka |
|
|
cos ωt kr þ θo þ 2 |
|
|||
r |
|
||||||
|
ka |
|
|
a |
π |
|
|
uðr, tÞ Ua |
|
|
|
cos ωt kr þ θo þ 2 |
ϕ |
||
cos ðϕÞ |
r |
The error of this simplification for the small spherical source can be calculated by using the amplitude and phase. The errors of amplitude and phase are:
|
|
|
|
|
a |
|
|
q |
|
|
|
||||||
Eamplitude |
¼ |
ka cos ðϕaÞ |
¼ |
1 |
þ ð |
ka |
Þ |
2 |
|
1 |
|||||||
|
|||||||||||||||||
|
cos |
ϕ |
|
Þ |
|
|
|
|
|||||||||
|
π |
|
ð |
|
|
π |
|
|
|
|
|
|
1 |
|
|
|
|
Ephase ¼ |
ðϕa þ kaÞ ¼ |
|
tan 1 |
|
|
þ ka |
|||||||||||
2 |
2 |
|
|||||||||||||||
ka |
Define a total error of the approximation as:
q
E ¼ E2amplitude þ E2phase
Depending on the application, the threshold of the total error can be determined. For example, if ka ¼ 0.1, E ¼ 0.005; if ka ¼ 0.01, E ¼ 0.00005.
6.3.1Near-Field Solutions of a Small Spherical Source (kr 1)
It is also of interest to study the behavior of sound pressure and velocity at short distances (r 1) from the source or at very low frequencies (ω 1) of a wave.
When the solution (measurement point) is very close to the source, kr and ϕ are approximately given by:
kr ¼ ωc r 0
ϕ ¼ tan 1 kr1 π2
kr
cos ðϕÞ ¼ q kr
1 þ ðkrÞ2
Substituting the above approximated values back to the pressure and velocity equations induced by the pulsating small sphere yields:
6.3 Acoustic Waves from a Small Pulsating Sphere |
143 |
|||||
|
a |
|
π |
|||
pðr, tÞ Uaρoc ka |
|
|
|
cos ωt kr þ θo þ |
2 |
|
r |
||||||
|
a |
2 |
|
|||
uðr, tÞ Ua |
|
|
cos ðωt kr þ θoÞ |
|
||
r |
|
Hence, the pressure is nearly 90 out of phase with the velocity near the source.
6.3.2Far-Field Solutions of a Small Spherical Source (kr 1)
It is also of interest to study the behavior of sound pressure and velocity at large distances (r 1) from the source or at high frequencies (ω 1) of a wave. When the solution (measurement point) is far away from the source, kr and ϕ are approximated as:
kr ¼ ωc r 1
ϕ ¼ tan 1 1 0 kr
kr
cos ðϕÞ ¼ q 1
1 þ ðkrÞ2
Substituting the above approximated values back to the pressure and velocity relations induced by the pulsating small sphere yields:
|
|
a |
|
|
π |
|
pðr, tÞ Uaρoc ka |
|
cos ωt kr þ θo |
þ |
2 |
||
r |
||||||
a |
|
|
|
π |
|
|
uðr, tÞ Uaka |
|
|
cos ωt kr þ θo þ |
2 |
||
r |
|
A comparison of the above pressure and velocity shows that the pressure is in phase with the velocity at large distances from the source or at high frequencies of waves.
Example 6.1 (A Small Spherical Source)
An acoustic pressure p(r, t) is created by a surface vibration of a spherical source with a radius a ¼ 1001 ½m&. Given the surface velocity of the sphere as:
|
1 |
|
m |
1 |
|
m |
||
uðr ¼ a, tÞ ¼ Ua cos 2πft |
|
πh |
|
i |
¼ Re Uae jðωt 2 |
πÞh |
|
i |
2 |
s |
s |
144 |
6 Acoustic Waves from Spherical Sources |
where the surface velocity is Ua ¼ 20 [m/s] and frequency of radiation is f ¼ 680 [Hz].
Assuming that this spherical source can be treated as a small spherical source, calculate (a) the flow velocity, (b) acoustic pressure, (c) intensity, and (d) power radiating from the small spherical source.
Use 415 [rayls] for characteristic impedance (ρ0c) and 340 [m/s] for the speed of sound. Show units in the Meter-Kilogram-Second (MKS) system.
Example 6.1 Solution
The angular frequency can be calculated as:
ω ¼ 2πf ¼ 2π 680 h1si ¼ 1360π h1si
The wavenumber can be calculated as:
k c |
|
c |
|
|
|
|
m |
s |
|
|
4π m |
|||
¼ |
ω |
|
2πf |
2π |
|
680 |
1 |
|
h |
1 |
i |
|||
|
¼ |
|
¼ |
|
|
|
|
|
¼ |
|
||||
|
|
|
340 |
s |
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Part (a)
Compare the given the surface velocity of the sphere:
|
1 |
m |
||
uðr ¼ a, tÞ ¼ Ua cos 2πft |
|
π h |
|
i |
2 |
s |
to the boundary condition (BC) of a small pulsating spherical source in the summary table below:
Source types |
Pressure: p(r, t) |
|
Velocity: u(r, t) |
|
|
|
|
|
|
|
|
|
|
||||||||
Small pulsating |
rp(r, |
|
|
|
|
ruðr, tÞ |
|
A |
cos ðωt kr þ θ ϕÞ |
||||||||||||
|
|
|
|
ρo c cos ðϕÞ |
|||||||||||||||||
spherical source |
t) |
|
A cos (ωt |
|
kr + θ) |
||||||||||||||||
where ϕ |
¼ |
tan 1 |
1 |
|
! |
cos |
ϕ |
|
kr |
||||||||||||
|
|
|
|
|
|
kr |
|
|
2 |
|
|||||||||||
a 1 |
where |
|
2 |
|
|
|
|
|
ð Þ ¼ p1þðkrÞ |
||||||||||||
A |
|
Ua ρoc ka |
|
Near field: |
|
|
|
|
|
|
|
|
|
|
|||||||
Based on BC: |
¼ |
|
|
|
|
|
|
|
Far field: |
|
|||||||||||
u(a, t)¼ |
|
|
π |
|
|
cos(ϕ) kr |
|
|
|
|
|
cos(ϕ) 1 |
|
|
|||||||
θ ¼ θo þ |
2 |
|
|
|
|
|
|
|
|
|
|||||||||||
Ua cos (ωt + θo) |
|
|
|
|
|
|
Remarks: |
|
|
|
|
|
|
Remarks: |
|
|
|||||
!ka 0π |
|
|
|
|
|
|
•π p-u out of phase |
|
• p-u in phase, |
|
|
||||||||||
! ϕa 2 |
|
|
|
|
|
|
by 2 |
|
|
|
|
|
|
|
,ϕ ¼ 0 |
|
|
||||
! cos (ϕa) ka |
|
|
|
|
|
|
spherical wave |
|
|
|
|
spherical wave |
Comparing the given BC to the formulas in the summary table above gives:
a ¼ 1001 ½m&
6.3 Acoustic Waves from a Small Pulsating Sphere |
145 |
Ua ¼ 20½m=s&
π
θo ¼ 2
The formula of velocity from a small source as shown is shown in the summary table as:
ruðr, tÞ |
|
A |
|
|
||
|
|
cos ðωt kr þ θ ϕÞ |
||||
ρoc cos ðϕÞ |
||||||
where |
|
|
|
|
|
|
A ¼ Ua ρoc ka2 ¼ 20 ∙ 415 ∙ 4π ∙ 0:012 ¼ 10:4 |
||||||
|
|
π |
π |
π |
||
θ ¼ θo þ 2 ¼ |
2 |
þ 2 ¼ 0 |
||||
Therefore: |
|
|
|
|
|
|
|
10:4 |
|
|
|
|
|
ruðr, tÞ |
|
cos ðωt kr ϕÞ |
||||
415 cos ðϕÞ |
||||||
|
0:0251 |
|
|
|
m |
h i
! uðr, tÞ r cos ðϕÞ cos ð1360πt 4πr ϕÞ s
where:
ϕ ¼ tan 1 1 kr
Part (b)
The formula for acoustic pressure from a small source is shown in the summary table as:
pðr, tÞ |
|
A |
cos ðωt kr þ θÞ½Pa& |
|
|
r |
|||
|
10:4 |
cos ð1360πt 4πrÞ ½Pa& |
||
|
|
r |
Part (c)
The formula of the intensity of spherical waves is derived in Chap. 5 as: