- •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
134 |
5 Solutions of Spherical Wave Equation |
imaginary number; when kr 1, ϕ 0, z ¼ ρoc is a real number and is the same as a plane wave.
5.7Homework Exercises
Exercise 5.1
Acoustic pressure magnitude P0 of a spherical wave is given as 4.149 [Pa] at 1 [m] from the center of a source. If the wavelength is 0.1 [m], answer the following questions, and show units in the Meter-Kilogram-Second (MKS) system:
a)What is the particle velocity magnitude at any point in space?
b)What is the acoustic pressure magnitude at any point in space?
c)What is the acoustic intensity at 1 and 25 [m] from the center of the source?
d)What is the phase difference between the pressure and flow velocity of the sphere wave at 1 and 25 [m]?
Use 415 [rayls] for the characteristic impedance (ρ0c) of air. (Answers):
p
a) The velocity magnitude U |
0:0001591 |
1þ |
3948r2 |
[m/s] |
|
|
2 |
||||
r ¼ |
|
r |
|
|
b)The pressure magnitude at r: Pr ¼ 4:149r [Pa]
c)The acoustic intensity at r ¼ 1 [m]: 0.02074 [ w/m2]
The acoustic intensity at r ¼ 25 [m]: 0.00003318 [w/m2]
d) The phase difference at r ¼ 1 [m]: ϕ1 ¼ 0.01591 rad ¼ 0.9118 [Deg] The phase difference at r ¼ 25 [m]: ϕ25 ¼ 0.00064 rad ¼ 0.0365 [Deg]
Exercise 5.2
The acoustic pressure magnitude from a spherical source, measured at 10 m from the center of the source, is 0.8 Pa. If the frequency of radiation is 250 Hz:
a)What is the surface velocity magnitude of the source if the source radius is 0.05 [m]?
b)What is the acoustic pressure magnitude on the surface of the source?
Use 415 [rayls] for the characteristic impedance (ρ0c) of air and 340 ms for the speed of sound (c) in air.
(Answers): (a) 1:713 ms ; (b) 160[Pa]
Exercise 5.3
Flow velocity magnitude U0 of a spherical wave is given as 0.02 [m/s] at 2 [m] from the center of a source. If the frequency is 1700 [Hz].
5.7 Homework Exercises |
135 |
a)What is the flow velocity magnitude at any point in space?
b)What is the acoustic pressure magnitude at any point in space?
c)What is the acoustic intensity at 5 [m] from the center of the source?
Use 415 [rayls] for the characteristic impedance (ρ |
|
|
m |
for the |
|||||||
speed of sound (c) in air. |
|
|
r |
|
s |
0c) of air and 340 |
s |
|
m |
||
|
|
r |
½ & |
|
|
||||||
|
|
3 p |
100π2r2 |
|
|
|
|
|
|||
(Answers): (a) 1:2731 10 |
|
1þ 2 |
|
m ; (b) 16:598 |
Pa ; (c) 0:0133 |
|
w |
|
|||
|
|
2 |
Exercise 5.4
Derive the complex flow velocity u(r, t) from the following complex sound pressure p(r, t) using Euler’s force equation:
|
pðr, tÞ ¼ |
A |
ðCEPÞ |
Pressure : |
r e jðωt krþθÞ |
Hint: The flow velocity u(r, t) can be obtained by finding the integral of the derivative of the pressure p(r, t) in the Euler’s force equation shown below:
Euler0s force equation : |
|
∂ |
|
∂ |
||||
|
|
pðr, tÞ ¼ ρo |
|
uðr, tÞ |
||||
∂r |
∂t |
|||||||
(Answers): Flow velocity : u r, t |
|
|
|
A |
e jðωt krþθ ϕÞ |
|||
Þ ¼ rρoc cos ϕ |
||||||||
ð |
|
|
|
Exercise 5.5
Derive the complex flow velocity u(r, t) from the following complex sound pressure p(r, t) using Euler’s force equation:
|
|
A |
|
|
|
||
Pressure : |
pðr, tÞ ¼ |
|
e jðωt krþθÞ |
ðCEPÞ |
|||
r |
|||||||
(Answers): Flow velocity : u r, t |
|
A |
e jðωt krþθ ϕÞ |
||||
Þ ¼ rρoc cos ϕ |
|||||||
|
|
ð |
|
Chapter 6
Acoustic Waves from Spherical Sources
In the previous chapter, formulas for sound pressure, flow velocity, acoustic intensity, and specific acoustic impedance were formulated in spherical coordinates. Note that in the previous chapter, the formulas are just general solutions and are not based on any physical source. In this chapter, you will learn how to calculate sound pressure radiated from a spherical source of different radii. Note that this spherical source has physical quantities such as radius and surface vibration velocity.
The end goal of this chapter is to obtain the formula for calculating sound pressure radiating from a point source. Even though this formula is derived from a spherical source with physical quantities of radius and velocity, the radius is actually eliminated in this formula which makes it possible to model any shape for vibration surfaces. The formula for acoustic waves from a point source is summarized below and in the table on the next page:
|
|
|
ρ ck |
|
|
π |
|
pðr, tÞ ¼ |
|
o |
Qs cos ωt kr þ θo þ |
2 |
|
||
|
4πr |
|
|||||
|
|
k |
|
|
π |
|
|
uðr, tÞ ¼ |
|
Qs cos ωt kr þ θo þ 2 |
ϕ |
||||
4πr cos ðϕÞ |
where Qs is the source strength and is defined as:
Z
Qs Unds ¼ 4πa2Ua
s
This formula can be used for numerical calculations of sound pressure radiating from an arbitrary vibrating surface. The vibrating surface is numerically modeled as a large number of singular point sources. Each point source has a source strength representing the vibration power of its corresponding vibration surface. This technique is used in the project in this chapter.
© The Author(s), under exclusive license to Springer Nature Switzerland AG 2021 |
137 |
H. Lin et al., Lecture Notes on Acoustics and Noise Control, https://doi.org/10.1007/978-3-030-88213-6_6
138 |
6 Acoustic Waves from Spherical Sources |
To obtain the formulas for point sources, we will follow the three formulas as described below:
Formula 1: Formula for a pulsating spherical source with an arbitrary radius a Formula 2: Formula for a pulsating spherical source with a small radius (a 1) Formula 3: Formula for a point source with a zero radius (a ¼ 0)
Each formula is developed based on the previous formula. For example, Formula 3 is based on Formula 2; Formula 2 is based on Formula 1; Formula 1 is based on the general spherical wave derived in the previous chapter.
The following is a summary of formulations of acoustic waves from a spherical source that will be derived in this chapter:
Source types |
|
|
|
Pressure: p(r, t) |
|
|
Velocity: u(r, t) |
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
General spherical |
|
|
rp(r, t) |
|
|
|
|
|
|
|
ruðr, tÞ ¼ |
|
|
|
A |
|
|
cos ðωt kr þ θ ϕÞ |
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
ρo c cos ðϕÞ |
||||||||||||||||||||||||||||
wave |
|
|
|
|
|
|
¼ |
A cos (ωt |
|
kr + θ) |
|
||||||||||||||||||||||||||
|
|
|
|
|
|
|
where |
ϕ |
¼ |
tan 1 |
|
1 |
! |
cos |
ϕ |
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
kr |
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
where r is the distance |
|
|
kr |
|
|
|
|
|
|
|
|
|
ð Þ ¼ |
||||||||||||||
|
|
|
|
|
|
|
|
|
|
from the origin of the |
|
p1þðkrÞ2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
coordinate |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
A and θ are determined by |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
BC |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Pulsating spherical |
|
|
rp(r, t) |
|
|
|
|
|
|
|
ruðr, tÞ ¼ |
|
|
|
A |
|
|
cos ðωt kr þ θ ϕÞ |
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
ρo c cos ðϕÞ |
||||||||||||||||||||||||||||
source of radius a |
|
|
¼ |
A cos (ωt |
|
kr + θ) |
|
||||||||||||||||||||||||||||||
|
|
|
where |
ϕ |
¼ |
tan |
1 |
|
1 |
|
|
|
|
|
|||||||||||||||||||||||
Based on BC: |
|
|
|
|
|
|
|
|
|
|
|
|
kr |
|
|
|
|
||||||||||||||||||||
|
|
|
where |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
u(a, |
|
a |
|
|
o |
) |
A ¼ Uaρoc cos (ϕa)a |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
t) |
¼ |
|
|
|
|
¼ |
ka + |
1 |
1 |
ϕa |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
U |
cos (ωt + θ |
θ |
|
θo |
+ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
ϕa ¼ tan |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
ka |
ka |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
! cos ðϕaÞ ¼ |
p1þðkaÞ2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
Small pulsating |
|
|
rp(r, t) |
|
|
|
|
|
|
ruðr, tÞ ¼ |
|
|
|
A |
|
|
cos ðωt kr þ θ ϕÞ |
||||||||||||||||||||
|
|
|
|
|
|
|
ρo c cos ðϕÞ |
||||||||||||||||||||||||||||||
spherical source if |
|
|
¼A cos (ωt kr + θ) |
|
|||||||||||||||||||||||||||||||||
|
|
|
where |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
ka |
|
|
1 |
|
|
|
|
|
where |
|
|
|
|
2 |
|
|
ϕ |
¼ |
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|||
Based on BC: |
|
|
|
A |
|
Ua ρoc ka |
|
|
fi |
|
|
|
|
|
|
|
|
|
fi |
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
tan 1 |
|
kr |
|
|
|
|
|
|
|
|
|
|
||||||
u(a, t)¼ Ua cos (ωt + θo) |
|
|
|
π |
|
|
|
|
|
Near |
|
eld: |
|
|
|
|
|
|
|
|
Far eld: |
|
|||||||||||||||
θ θo þ |
2 |
|
|
|
|
|
kr 0 |
|
|
|
|
|
|
|
|
|
|
kr 1 |
|
||||||||||||||||||
! |
ka |
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
kr |
|
|
|
|
|
|
|
1 |
||||||||||
|
|
π |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ϕ) |
|
|
|
|
|
|
|
|
cos(ϕ) |
|
|||||||||
! ϕa 2 |
|
ka |
|
|
|
|
|
|
|
|
|
|
|
|
cos( |
π |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
! |
cos (ϕa) |
|
|
|
|
|
|
|
|
|
|
|
|
|
ϕ |
2 |
(out of phase) |
|
ϕ 0 (in phase) |
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
Point source |
|
|
|
rp(r, t) |
|
|
|
|
|
|
|
ruðr, tÞ ¼ |
|
|
|
A |
|
|
cos ðωt kr þ θ ϕÞ |
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
ρo c cos ðϕÞ |
|||||||||||||||||||||||||||
Based on BC: |
|
|
|
|
A cos (ωt |
|
|
kr + θ) |
|
||||||||||||||||||||||||||||
|
|
|
|
|
|
|
where |
ϕ |
¼ |
tan 1 |
|
1 |
|
|
|
|
|
||||||||||||||||||||
u(a, t) Ua cos (ωt + θo) |
where A |
|
Qs |
o4π |
|
|
|
|
|
|
|
||||||||||||||||||||||||||
|
|
|
|
¼ |
|
|
|
|
|
¼ |
|
|
¼π |
|
ρ c k |
|
|
|
|
kr |
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
2 |
Ua |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
Qs R sUnds ¼ 4πa |
|
θ ¼ θo þ 2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|