- •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
124 5 Solutions of Spherical Wave Equation
rpðr, tÞ ¼ ψðr, tÞ ¼ A cos ðωt kr þ θÞ |
ðREPÞ |
||
¼ |
1 |
hAe jðωt krþθÞ þ Ae jðωt krþθÞi |
ðCEPÞ |
2 |
5.4Flow Velocity
Similar to the plane wave, the flow velocity can be obtained from sound pressure using Euler’s force equation.
The one-dimensional Euler’s force equation in the spherical coordinates can be obtained by replacing the vector gradient operator with the spherical coordinate gradient operator:
p ¼ ρo |
∂ |
|
|
ðVector formatÞ |
|||
|
u |
|
|
||||
∂t |
|
|
|||||
|
∂ |
∂ |
|
||||
! |
|
pðr, tÞ ¼ ρo |
|
uðr, tÞ |
ðSpherical coordÞ |
||
∂r |
∂t |
5.4.1Flow Velocity in Real Format
Sound pressure p(r, t) in the spherical coordinate system was derived in the previous section as:
rpðr, tÞ ¼ A cos ðωt kr þ θÞ
Flow velocity u(r, t) in spherical coordinates can be derived from Euler’s force equation and expressed in real explicit phase [REP] format as:
|
|
A |
|
||
ruðr, tÞ ¼ |
|
|
cos ðωt kr þ θ ϕÞ |
ðREPÞ |
|
ρoc cos ϕ |
|||||
or: |
|
|
|
|
|
uðr, tÞ ¼ |
|
A |
|
||
|
cos ðωt kr þ θ ϕÞ |
ðREPÞ |
|||
rρoc cos ϕ |
The following will derive the above equation in real explicit phase [REP] format:
5.4 Flow Velocity |
125 |
Flow velocity expressed in real explicit phase [REP] format can be derived independently using Euler’s force equation:
∂∂r pðr, tÞber ¼ ρo ∂∂t uðr, tÞ
For a given pressure p(x, t), flow velocity u(x, t) can be calculated as:
uðr, tÞ ¼ ρo |
Z ∂r pðr, tÞ dt |
1 |
∂ |
Flow velocity in real number format will be derived in this section:
ruðr, tÞ ¼ |
A |
ðREPÞ |
ρoc cos ϕ cos ðωt kr þ θ ϕÞ |
Flow velocity formulated in real explicit phase [REP] format is derived from pressure p(r, t) as:
rpðr, tÞ ¼ A cos ðωt kr þ θÞ
Substituting the above pressure p(r, t) into the formula for flow velocity yields:
ð |
|
Þ ¼ |
ρo |
|
Z |
∂r |
ð |
|
Þ |
|
¼ |
|
ρo |
Z |
|
∂r hr |
|
|
|
ð |
|
|
|
|
|
þ |
|
Þi |
|||||||||||||||||||||||||
u |
r, t |
|
|
1 |
|
|
|
|
∂ |
p |
r, t |
|
|
|
|
dt |
|
|
|
A |
|
|
|
|
∂ |
1 |
|
cos |
|
ωt |
|
|
kr |
|
|
|
θ |
|
|
dt |
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||
Carrying out the derivative with respect to r yields: |
|
|
|
|
|
|
|
|
|
|
|
|
Þ |
|
|
||||||||||||||||||||||||||||||||||||||
|
ð |
|
Þ ¼ |
|
ρo |
Z r |
|
|
ð |
|
|
|
|
|
þ |
|
|
Þ r2 |
|
|
|
ð |
|
|
|
|
|
|
þ |
|
|
|
|||||||||||||||||||||
|
u |
r, t |
|
|
|
A |
|
|
|
k |
sin |
|
|
ωt |
|
|
kr |
|
|
θ |
|
|
|
1 |
|
cos |
|
ωt |
|
|
|
kr |
|
|
θ |
|
|
dt |
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||
Carrying out the integral for t yields: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
ð |
|
Þ ¼ |
|
r |
|
ω |
|
|
|
|
ð |
|
|
|
|
|
þ |
|
|
Þ |
|
r2 |
|
ω |
|
|
|
|
ð |
|
|
|
|
|
|
þ |
|
||||||||||||||||
u |
r, t |
|
|
A |
|
k |
|
1 |
|
cos |
|
|
ωt |
|
|
kr |
|
|
θ |
|
|
1 |
|
|
1 |
|
|
sin |
|
|
ωt |
|
|
|
kr |
|
|
θ |
|
||||||||||||||
|
|
|
|
¼ rρo |
ω h cos ðωt kr þ θÞ þ kr sin ðωt kr þ θÞi |
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
A |
k |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Simplify the above equation and replace |
ω with c to arrive at: |
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
k |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
uðr, tÞ ¼ rρoc h cos ðωt kr þ θÞ þ kr sin ðωt kr þ θÞi |
|
|
|
||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
A |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Combining the cosine and sine functions into one cosine function yields:
126 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
5 Solutions of Spherical Wave Equation |
|||||||
|
|
|
uðr,tÞ¼ |
|
A |
|
hcosðωt krþθÞþ |
|
1 |
sinðωt kr þ |
θÞi |
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
|
|
rρoc |
kr |
|
|||||||||||||||||||||||
|
|
|
|
¼ |
A |
|
1 |
½krcosðωt kr þθÞþsinðωt kr þθÞ& |
|
||||||||||||||||||
|
|
|
|
|
|
|
|||||||||||||||||||||
|
|
|
|
rρockr |
3 |
||||||||||||||||||||||
|
|
q |
|
|
|
|
|
|
|
|
|
cos ωt |
kr |
θ |
|
|
|
sin ωt kr θ |
|||||||||
|
A |
1þðkrÞ2 2 |
|
|
kr |
|
|
ð þ Þþ |
1 |
|
|
||||||||||||||||
¼rρoc |
|
kr |
6 |
1þðkrÞ2 |
|
1þðkrÞ2 |
|
ð þ Þ7 |
|||||||||||||||||||
|
|
|
|
|
|
4q |
|
|
|
|
|
|
|
|
q |
|
|
5 |
|||||||||
Define a new variable ϕ as: |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
cos ϕ |
¼ |
|
|
|
|
kr |
|
; |
sin |
ϕ ¼ |
1 |
|
|
|
||||||||
|
|
|
|
|
q |
q |
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
1 þ ðkrÞ2 |
1 þ ðkrÞ2 |
|
|
|||||||||||||||
Then, u(r, t) can be simplified as: |
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
uðr, tÞ ¼ |
A |
1 |
|
½ cos ðϕÞ cos ðωt kr |
þ θÞ þ sin ðϕÞ sin ðωt kr þ θÞ& |
||||||||||||||||||||||
|
|
|
|||||||||||||||||||||||||
rρoc |
|
cos ðϕÞ |
|||||||||||||||||||||||||
|
|
! ruðr, tÞ ¼ |
|
A |
|
|
|
|
1 |
|
½ cos |
ðωt kr þ θ ϕÞ& |
|
|
ð5:4Þ |
||||||||||||
|
|
|
ρoc |
|
cos ϕ |
|
|
where:
ϕ ¼ tan 1 1 kr
Eq. (5.4) is the flow velocity of spherically symmetric acoustic wave, for a given pressure function, as described in Eq.(5.3).
5.4.2Flow Velocity in Complex Format
Sound pressure p(r, t) in the spherical coordinate system was derived in the previous ion as:
rpðr, tÞ ¼ A cos ðωt kr þ θÞ
It can be formulated in a complex explicit phase [CEP] as shown below using Euler’s formula:
5.4 Flow Velocity |
127 |
rpðr, tÞ ¼ A |
1 |
ðCEPÞ |
2 he jðωt krþθÞ þ e jðωt krþθÞi |
The corresponding flow velocity in a complex exponential format as shown below can be derived using Euler’s force equation:
ruðr, tÞ ¼ |
A |
1 |
he jðωt krþθ ϕÞ þ e jðωt krþθ ϕÞi |
ðCEPÞ |
|
ρoc cos ϕ |
|
2 |
where:
ϕ ¼ tan 1 1 kr
The derivation of the flow velocity formulated in complex explicit phase [CEP] format is not derived here. The derivation is part of the homework in this chapter.
Example 5.2
What is the phase difference between the pressure and the flow velocity of a
spherical wave at 3 [m] from the origin at 1000 [Hz]?
Use 340 ms for the speed of sound (c) in air.
Example 5.2 Solution The general form of pressure of a spherical wave is:
rpðr, tÞ ¼ A cos ðωt kr þ θÞ
The formula of flow velocity is derived from pressure using Euler’s force equation:
A
ruðr, tÞ ¼ ρoc cos ϕ cos ðωt kr þ θ ϕÞ
where:
ϕ ¼ tan 1 1 kr
Hence, by inspecting the exponents of pressure and flow velocity, the phase difference between pressure and flow velocity of a spherical wave is equal to ϕ, with flow velocity lagging behind the pressure.
For the given conditions:
f ¼ 1000 ½Hz&