- •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
7.3 2D Boundary Conditions Between Four Walls |
177 |
Mode 2: (l ¼ 2)
λl ¼ 2π ¼ 2L ! λ2 ¼ L kl l 2 2
Mode 3: (l ¼ 3)
λl ¼ 2π ¼ 2L ! λ3 ¼ L kl l 2 3
Note that mode shapes in a pipe are standing waves which are equivalent to two traveling waves with the same amplitude and traveling in opposite directions.
7.32D Boundary Conditions Between Four Walls
7.3.12D Standing Wave Solutions of the Wave Equation
The two-dimensional wave equation in Cartesian coordinate is:
2p ¼ 1 ∂2p c2 ∂t2
!∂2p þ ∂2p ¼ 1 ∂2p
∂x2 ∂y2 c2 ∂t2
178 |
7 Resonant Cavities |
The differential equation is separable if a solution can be cast in the following form that satisfies the acoustic wave equation:
psðx, y, tÞ ¼ TðtÞXðxÞ YðyÞ
Substituting the above solution into the acoustic wave equation above yields the three separate ordinary differential equations (ODEs):
1 |
T}ðtÞ |
¼ |
X}ðxÞ |
þ |
Y}ðyÞ |
¼ |
k2 |
|
k2 |
k2 |
|
||||||||||
c2 TðtÞ |
XðxÞ |
YðyÞ |
x |
y |
¼ |
where k is a new variable that relates the spatial and temporal parts of the differential equations. The new variable k is separated into two independent variables kx and ky for two independent differential equations in two independent directions:
1 T}ðtÞ ¼ k2 c2 TðtÞ
X}ðxÞ ¼ k2
XðxÞ x
Y}ðyÞ ¼ k2
YðyÞ y
where the new variable k can be considered a combined wavenumber that related to the wavenumbers kx and ky as:
k2 ¼ k2x þ k2y
With the new variables k, kx, and ky, the solutions of the temporal and the spatial differential equations are:
TðxÞ ¼ At cos ðωt þ θtÞ
XðxÞ ¼ Ax cos ðkxx þ θxÞ
YðxÞ ¼ Ay cos kyy þ θy
Therefore, the standing wave solution of the 2D wave equation is:
psðx, y, tÞ ¼ TðtÞXðxÞYðyÞ
¼ AtAxAy cos ðωt þ θtÞ cos ðkxx þ θxÞ cos kyy þ θy
7.3 2D Boundary Conditions Between Four Walls |
179 |
According to Euler’s force equation, the velocity of the 2D standing wave is:
where:
uxðx,
uyðx,
!usðx, y, tÞ ¼ |
1 |
|
Z ½ psðx, y, tÞ&dt |
|
|||||||||||||||
|
ρo |
|
|
||||||||||||||||
|
|
¼ |
1 |
|
Z |
|
ex |
|
∂ |
þ ey |
∂ |
pðx, y, tÞ dt |
|||||||
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
ρo |
|
∂x |
∂y |
|||||||||||||
|
|
|
|
|
ux |
|
|
x, |
y, |
t e |
xb |
u |
y |
x, y, t e |
y |
||||
|
|
¼ |
ð |
|
|
|
|
b |
|
||||||||||
|
|
|
|
|
|
Þb |
þ |
|
|
ð |
Þb |
||||||||
|
1 |
|
k |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
y, tÞ ¼ |
|
|
x |
AtAxAy sin ðωt þ θtÞ sin ðkxx þ θxÞ cos kyy þ θy |
|||||||||||||||
ρ0c |
k |
||||||||||||||||||
|
1 |
|
k |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
y, tÞ ¼ |
|
|
y |
AtAxAy sin ðωt þ θtÞ cos ðkxx þ θxÞ sin kyy þ θy |
|||||||||||||||
ρ0c |
|
k |
Note that the pressure is a scaler, and the velocity is a vector as indicated in the equation above.
The detailed derivation of the 2D standing wave velocity above is shown in Appendix 7.7.1.
7.3.22D Nature Frequencies Between Four Walls
Assume these two walls are oriented parallel to each other with the surface normal to the x-axis (y-axis). Also, assume the left (bottom) wall is located at x ¼ 0 (y ¼ 0) and the right (top) wall is located at ¼Lx (y ¼ Ly) as shown below:
= 0 |
= 0 |
|
|
= 0
̂
= 0
180 |
7 Resonant Cavities |
The general 2D standing wave pressure and velocity are shown below:
psðx, y, tÞ ¼ AtAxAy cos ðωt þ θtÞ cos ðkxx þ θxÞ cos kyy þ θy
!
usðx, y, tÞ ¼ uxðx, y, tÞbex þ uyðx, y, tÞbey
where:
uxðx, y, tÞ ¼ ρ1 kx AtAxAy sin ðωt þ θtÞ sin ðkxx þ θxÞ cos kyy þ θy 0c k
uyðx, y, tÞ ¼ ρ1 ky AtAxAy sin ðωt þ θtÞ cos ðkxx þ θxÞ sin kyy þ θy 0c k
The boundary condition at the left wall (CLOSED at x = 0):
uxðx, y, tÞjx¼0 ! θx ¼ 0
The boundary condition at the right wall (CLOSED at x = Lx):
uxðx, y, tÞjx¼Lx ¼ sin ðkxlLxÞ ¼ 0, kxl ¼ |
π |
l, |
l ¼ 1, 2, . . . |
|
|
|
||||||||||||||||
Lx |
|
|
|
|||||||||||||||||||
The boundary condition at the bottom wall (CLOSED at y = 0): |
||||||||||||||||||||||
|
|
|
|
|
uyðx, y, tÞ y¼0 ! θy ¼ 0 |
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
The boundary condition at the top wall (CLOSED at y = Ly): |
|
|||||||||||||||||||||
uyðx, y, tÞ y¼Ly |
|
¼ sin kymLy ¼ 0, |
|
|
|
π |
|
|
|
|
m ¼ 1, 2, . . . |
|||||||||||
|
kym ¼ Ly m, |
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Therefore, the discretized flow velocity becomes: |
|
|
|
|
|
|
|
|
|
|||||||||||||
! |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
usðx, y, tÞ ¼ uxðx, y, tÞex þ uyðx, y, tÞey |
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
uxðx, y, tÞ ¼ |
|
1 |
|
kx |
A |
A A sin |
ð |
ωt |
θ |
|
sin |
ð |
k |
|
x |
Þ |
cos |
k |
|
y |
||
|
ρ0c k |
|
|
|
||||||||||||||||||
|
t |
bx y |
|
þb tÞ |
|
|
xl |
|
|
|
ym |
|
||||||||||
|
1 |
|
k |
|
|
|
|
|
þ θtÞ cos ðkxlxÞ sin kymy |
|||||||||||||
uyðx, y, tÞ ¼ |
|
|
|
|
y |
AtAxAy sin |
ðωt |
|||||||||||||||
|
ρ0c |
k |
||||||||||||||||||||
where: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
π |
|
|
|
π |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
kxl ¼ |
|
l; |
|
kym ¼ |
|
m; |
|
|
l, m ¼ 1, 2, . . . |
|
|
|||||||||||
Lx |
|
Ly |
|
|
|
|
And the discretized sound pressure becomes: