- •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
3.1 Review of Partial Differential Equations |
51 |
3.1Review of Partial Differential Equations
3.1.1Complex Solutions of a Partial Differential Equation
The acoustic wave equation is a second-order partial differential equation. The acoustic wave equation for a one-dimensional plane wave in Cartesian coordinates is:
∂2 pðx, tÞ ¼ 1 ∂2 pðx, tÞ ∂x2 c2 ∂t2
where p is sound pressure, x is a dimension, t is time, and c is the propagation speed of sound.
It can be directly shown (see Example 3.1) that when ωk ¼ c, the following four basic complex solutions:
pðx, tÞ ¼ Ae jðωtþkxþϕÞ
pðx, tÞ ¼ Ae jðωtþkxþϕÞ
pðx, tÞ ¼ Ae jðωt kxþϕÞ
pðx, tÞ ¼ Ae jðωt kxþϕÞ
satisfy the one-dimensional acoustic wave equation of p(x, t):
∂2 pðx, tÞ ¼ 1 ∂2 pðx, tÞ
∂x2 c2 ∂t2
where A is an amplitude (real number) and ϕ is a phase shift. The detailed derivation will be shown in Sect. 3.2. The following example is proof of the solutions:
Example 3.1 A complex exponential function p(x, t) is given as:
pðx, tÞ ¼ Ae jðωtþkxþϕÞ
where A is the amplitude (real number) and ϕ is the phase (angle). Show that when ωk ¼ c, the above complex exponential function p(x, t) satisfies the one-dimensional acoustic wave equation:
∂2 pðx, tÞ ¼ 1 ∂2 pðx, tÞ
∂x2 c2 ∂t2
52 3 Solutions of Acoustic Wave Equation
Example 3.1 Solution Step 1: The left-hand side of the given acoustic wave equation is ∂∂x22 pðx, tÞ.
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|||||||||
Calculate |
|
∂ |
p by substituting the given function p(x, t) into |
∂ |
pðx, tÞ to arrive at: |
||||||||||||||||||||||||
∂x2 |
∂x2 |
||||||||||||||||||||||||||||
|
|
|
|
|
∂2 |
|
∂ |
|
|
|
∂ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
∂x2 p ¼ |
|
∂x |
∂x Ae jðωtþkxþϕÞ |
|
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
h |
|
|
|
|
i |
|
|
|
||||||
|
|
|
|
|
|
|
|
∂ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
¼ |
|
|
|
|
njkhAe jðωtþkxþϕÞio |
|
|
|
|||||||||||||||||||||
|
∂x |
|
|
|
|||||||||||||||||||||||||
|
|
|
|
|
|
¼ k2hAe jðωtþkxþϕÞi |
|
∂2 |
|
||||||||||||||||||||
Step 2: The2 right-hand side of the given acoustic wave equation2 |
is |
pðx, tÞ. |
|||||||||||||||||||||||||||
∂t2 |
|||||||||||||||||||||||||||||
Calculate |
∂ |
p by substituting the given function p(x, t) into |
∂ |
pðx, tÞ to get: |
|||||||||||||||||||||||||
∂t2 |
∂t2 |
||||||||||||||||||||||||||||
|
|
|
|
|
∂2 |
|
∂ |
|
|
|
∂ |
Ae jðωtþkxþϕÞ |
|
|
|
||||||||||||||
|
|
|
|
|
∂t2 p ¼ |
|
∂t |
∂t |
|
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
h |
|
|
|
|
i |
|
|
|
||||||
|
|
|
|
|
|
|
|
∂ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
¼ |
|
|
|
hjωAe jðωtþkxþϕÞi |
|
|
|
||||||||||||||||||||||
|
∂t |
|
|
|
|||||||||||||||||||||||||
|
|
|
|
|
|
¼ ω2hAe jðωtþkxþϕÞi |
|
|
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
2 |
|
|
|
|
|
|
|
|
|
||||||
Step 3: Substituting the calculated |
∂ |
p and |
∂ |
p into the acoustic wave equation |
|||||||||||||||||||||||||
2 |
2 |
||||||||||||||||||||||||||||
will give: |
|
|
|
|
∂x |
|
|
∂t |
|
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
∂2 |
1 |
|
|
∂2 |
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
p ¼ c2 |
|
|
p |
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
∂x2 |
∂t2 |
|
|
|
|||||||||||||||||
|
|
|
|
! k2hAe jðωtþkxþϕÞi |
1 |
ω2hAe jðωtþkxþϕÞi |
|
|
|
||||||||||||||||||||
|
|
|
|
¼ |
|
|
|
|
|||||||||||||||||||||
|
|
|
|
c2 |
|
|
|
||||||||||||||||||||||
|
|
|
|
|
|
|
! |
|
|
|
k2 |
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
¼ c2 |
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
ω2 |
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
! |
|
|
|
|
ω |
¼ c |
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
k |
|
|
|
|
|
|
|
|
|
|
|
|
It shows that when ωk ¼ c, the complex exponential function p(x, t) satisfies the one-dimensional acoustic wave equation.
3.1 Review of Partial Differential Equations |
53 |
3.1.2Trigonometric Solutions of a Partial Differential Equation
Similar to the complex solutions of the partial differential equation covered in the previous section, it can be directly shown that when ωk ¼ c, the following four basic complex solutions:
pðx, tÞ ¼ Aþc cos
pðx, tÞ ¼ Aþs sin
pðx, tÞ ¼ A c cos
pðx, tÞ ¼ A s sin
ωt kx þ ϕþωt kx þ ϕþ
ðωt þ kx þ ϕ Þ ðωt þ kx þ ϕ Þ
satisfy the one-dimensional acoustic wave equation of p(x, t):
∂2 pðx, tÞ ¼ 1 ∂2 pðx, tÞ
∂x2 c2 ∂t2
where A is an amplitude (real number) and ϕ is a phase shift. The detailed derivation will be shown in Sect. 3.2.
Example 3.2 A backward traveling wave p (x, t) (Form 1:RIP) is given as:
p ðx, tÞ ¼ A s sin ðωt þ kxÞ
where A s is an amplitude. Show that if ωk ¼ c, the above trigonometric function p(x, t) satisfies the one-dimensional acoustic wave equation:
∂2 pðx, tÞ ¼ 1 ∂2 pðx, tÞ ∂x2 c2 ∂t2
Example 3.2 Solution Step 1: The left-hand side of the given acoustic wave equation is:
∂2
∂x2 pðx, tÞ
Calculate ∂2 p by substituting the given function p(x, t) into ∂2 pðx, tÞ to get:
∂x2 ∂x2