- •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
150 |
|
|
|
|
6 Acoustic Waves from Spherical Sources |
|
|
|
ρ ck |
|
π |
|
|
pðr, tÞ ¼ |
|
o |
QH cos ωt kr þ θo þ 2 |
|
||
|
2πr |
|
||||
|
|
k |
|
π |
|
|
uðr, tÞ ¼ |
|
QH cos ωt kr þ θo þ 2 |
ϕ |
|||
2πr cos ðϕÞ |
Therefore, it requires only half of the source strength to produce the same pressure field as a full spherical pressure field. In other words, if the same source strength is put on the half space, half of the source strength will be reflected to the half space and produce twice the pressure field in the full spherical space created by the same source strength.
6.5Acoustic Intensity and Sound Power
Sound intensity I is the acoustic intensity, formulated in the previous section as:
|
1 |
T |
|
1 A2 |
ρoc k |
2 Qs2 |
||||||||
IðrÞ |
|
|
Z0 |
ðpuÞdt ¼ |
|
|
|
¼ |
|
|
|
|
|
|
T |
2ρoc |
r2 |
2 |
4π |
|
r2 |
where A was calculated for a point source as:
A ¼ Qs ρoc k
4π
For a spherical source with radius r ¼ a , acoustic intensity I is constant on the surface of the sphere. Therefore, sound power can be formulated without integration:
Z
w IðrÞ ds
s
¼ 4πa2Iðr ¼ aÞ
¼ 4πa |
2 1 |
|
A2 |
2πA2 |
ρock2 2 |
½w& |
||||
|
|
|
|
¼ |
|
¼ |
|
Qs |
||
|
2 |
r2ρoc |
ρoc |
8π |
Even though the above formula for sound power is formulated at r ¼ a, this formula is valid for any other choice of surface for integration. It makes sense because the sound power is the total radiation energy of the source and will not vary with the choices of the surface for integration. This can be checked by comparing the above sound power, at r ¼ a, to the sound power at r ¼ 2a as shown below:
6.5 Acoustic Intensity and Sound Power |
151 |
|
|
1 |
|
|
2 |
|
2πA2 |
|
w Z |
I ds ¼ 4πð2aÞ2Iðr ¼ 2aÞ ¼ 4πð2aÞ2 |
|
|
A |
|
¼ |
ρoc |
|
2 |
2 |
ρoc |
||||||
|
s |
|
ð2aÞ |
|
|
|
Example 6.2 (A Point 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 |
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 point source, calculate
(a) the source strength, (b) flow velocity, (c) acoustic pressure, (d) intensity, and
(e) power radiating from the point source.
Use 415 [rayls] for characteristic impedance (ρ0c) of air and 340 [m/s] for the speed of sound in air. Show units in the Meter-Kilogram-Second (MKS) system.
Example 6.2 Solution
Similar to Example 6.1, the angular frequency can be calculated as:
ω ¼ 2πf ¼ 2π 680 h1si ¼ 1360π h1si
And the wave number can be calculated as:
k c |
|
c |
|
|
|
|
m |
s |
|
|
4π m |
|||
¼ |
ω |
|
2πf |
2π |
|
680 |
1 |
|
h |
1 |
i |
|||
|
¼ |
|
¼ |
|
|
|
|
|
¼ |
|
||||
|
|
|
340 |
s |
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Compare the given surface velocity of the sphere:
|
1 |
m |
||
uðr ¼ a, tÞ ¼ Ua cos 2πft |
|
π h |
|
i |
2 |
s |
to the boundary condition of a small pulsating spherical source in the summary table below:
Source types |
Pressure: p(r, t) |
|
|
|
Velocity: u(r, t) |
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Point source |
rp(r, |
|
|
|
|
|
|
ruðr, tÞ ¼ |
|
|
A |
|
cos ðωt kr þ θ ϕÞ |
||||||||
|
|
|
|
|
|
ρo c cos ðϕÞ |
|||||||||||||||
Based on BC: |
t) ¼ A cos (ωt kr + θ) |
|
|||||||||||||||||||
|
where |
ϕ |
¼ |
tan 1 |
|
1 |
; |
cos |
ϕ |
||||||||||||
|
|
|
kr |
||||||||||||||||||
u(a, t) |
¼ |
|
where A |
¼ |
Qs ρo4cπ k; |
θ |
¼ |
|
kr |
|
|
|
|
|
|
|
ð Þ ¼ |
||||
|
|
|
π |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
p1 kr |
|
|
|
|
|
|
|
|
|
|
|
||||
Ua cos (ωt + θo) |
θo 2 |
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
||||
Qs ¼ |
2 |
Ua |
þ |
|
|
|
|
|
|
þð |
Þ |
|
|
|
|
|
|
|
|
|
|
4πa |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
152 |
6 Acoustic Waves from Spherical Sources |
Comparing the given boundary condition to the formulas in the summary table above gives:
a ¼ 1001 ½m&
Ua ¼ 20[m/s]
π
θo ¼ 2
Part (a)
When this source is modeled as a point source, the source strength is given by:
Qs ¼ 4πa2Ua ¼ 4π ∙ 0:012 ∙ 20 ¼ 0:0251 m3 s
Part (b)
The formula of the velocity from a point source is shown in the summary table as:
uðr, tÞ |
|
A |
|
|
|
cos ðωt kr þ θ ϕÞ |
|||||||||||
r ρoc cos ðϕÞ |
|||||||||||||||||
where: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
A |
¼ |
Q |
ρoc k |
|
¼ |
0:0251 |
415 4π |
¼ |
10:4 |
||||||||
4π |
|||||||||||||||||
|
|
s |
|
|
|
|
|
4π |
|
||||||||
|
|
|
θ ¼ θo þ |
π |
π |
π |
|
|
|||||||||
|
|
|
2 ¼ |
2 |
þ 2 ¼ 0 |
|
|||||||||||
Therefore: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
uðr, tÞ |
|
|
A |
|
|
|
|
cos ðωt kr þ θ ϕÞ |
|||||||||
|
|
|
|||||||||||||||
r ρoc cos ðϕÞ |
|||||||||||||||||
|
¼ |
|
|
10:4 |
|
|
|
cos |
ð1360πt 4πr ϕÞ |
||||||||
|
415 r cos ðϕÞ |
||||||||||||||||
|
|
0:0251 |
cos ð1360πt 4πr ϕÞ |
||||||||||||||
|
¼ |
|
|
||||||||||||||
|
r |
cos ðϕÞ |
where:
6.5 Acoustic Intensity and Sound Power |
153 |
ϕ ¼ tan 1 1 kr
Part (c)
The formula for the acoustic pressure from a point source is shown in the summary table as:
pðr, tÞ ¼ Ar cos ðωt kr þ θÞ
¼ 10r:4 cos ð1360πt 4πrÞ
Part (d)
Acoustic intensity:
|
|
|
|
|
ρoc k |
|
|
2 |
Qs2 |
|
|
w |
|
|
|
||||||
|
|
IðrÞ ¼ |
|
|
|
|
|
|
|
h |
|
|
i |
|
|
|
|||||
|
|
|
2 |
4π |
|
|
r2 |
m2 |
|
|
|
||||||||||
|
415 |
4π 2 |
0 |
02512 |
|
w |
|
|
0 131 |
|
w |
||||||||||
¼ |
|
|
|
|
|
: |
hm2i ¼ |
|
: |
h |
|
i |
|||||||||
2 |
4π |
|
r2 |
r2 |
m2 |
Part (e)
Sound power:
w |
¼ |
|
ρock2 |
Q2 |
w |
415 ð4πÞ2 |
0:02512 |
w |
1:65 w |
|
8π |
||||||||
|
|
s |
½ & ¼ |
8π |
|
½ & ¼ |
½ & |
Remarks
The results show that the acoustic quantities calculated from a small spherical formulation (Example 6.1) and a point source formulation (Example 6.2) are identical, as expected.
For a small sphere, the constant A is calculated by comparing the magnitudes of the velocity:
|
A |
|
k |
|
|
¼ |
|
Qs |
|
|
rρoc cos ðϕÞ |
4πr cos ðϕÞ |
to arrive at:
154 |
|
|
|
6 Acoustic Waves from Spherical Sources |
A ¼ |
ρock |
Qs ¼ |
415 ∙ 4π |
∙ 0:0251 ¼ 10:4 |
4π |
4π |
This shows that the constant A calculated from a point source is the same or close to the value calculated from a small spherical source. Therefore, we can conclude that this spherical source can be modeled as a point source.
Example 6.3
The rectangular plate has the width Lx and the height Ly with the following given dimensions:
Lx ¼ 0:8 ½m&
Ly ¼ 0:6 ½m&
= cos 2 +
Use m × n point sources
m=2 (x-direction) n=2 (y-direction)
at (0.1, 0, 0.1) [m]
= 0.8 [m] |
̂ |
|
= 0.6 [m] |
||
Drawing not to scale |
||
|
||
|
Assume that this plate vibrates at frequency f and the amplitude U of surface velocity is:
f ¼ 340 ½Hz&
hmi
U ¼ 20 s
θo ¼ 0 ½rad&
Therefore, the surface velocity is:
hmi uðtÞ ¼ U cos ð2πftÞ ¼ 20 cos ð2π 340 tÞ s
6.5 Acoustic Intensity and Sound Power |
155 |
We will equally divide the rectangular plate into m sections in the x-direction and n sections in the y-direction. The values of m and n are:
m ¼ 2; |
n ¼ 2 |
Calculate (a) sound pressure and (b) flow velocity (three components) generated by surface vibration of the plate at Point A (0.1, 0, 0.1) [m].
Use 415 [rayls] for characteristic impedance (ρ0c) of air and 340 [m/s] for the speed of sound in air. Show units in the Meter-Kilogram-Second (MKS) system.
You can use the function POINTSOURCE (in Sect. 6.6) for calculating sound pressure and flow velocity of a vibration plate. Or you can implement this function using your preferred tool such as Python or Excel spreadsheet.
Example 6.3 Solution
We will equally divide the rectangular plate into m sections in the x-direction and n sections in the y-direction. The values of m and n are:
m ¼ 2; |
n ¼ 2 |
Each divided piece has the same dimension as:
x ¼ |
Lx |
¼ |
0:8 |
½m& ¼ 0:4 ½m& |
|
|
m |
2 |
|
|
|||
y ¼ |
Ly |
¼ |
0:6 |
½m& ¼ 0:3 ½m& |
|
|
n |
2 |
s |
||||
Qs ¼ x y U ¼ 0:4 0:3 20 ¼ 2:4 |
||||||
|
|
|
|
|
|
m3 |
Due to the hemispherical point source, the source strength is doubled as:
QH ¼ Qs 2
Use a time resolution of t ¼ 0.0001 [s] for the time history plot.
The locations of the point sources are user-defined arguments in the main program as:
src(1,:)=[-0.2, 0.15, 0, 2.4]; % sx,sy,sz,Qs src(2,:)=[ 0.2, 0.15, 0, 2.4]; % sx,sy,sz,Qs src(3,:)=[-0.2,-0.15, 0, 2.4]; % sx,sy,sz,Qs src(4,:)=[ 0.2,-0.15, 0, 2.4]; % sx,sy,sz,Qs
The location of Point A is also a user-defined argument in the main program as:
x=0.1; y=0 ; z=0.1; % [m] Point A
156 |
6 Acoustic Waves from Spherical Sources |
The function POINTSOURCE.m and the main program VibrationPlate.m of this example can be found on Moodle.
The outputs of the first 11 time steps are:
time [s] |
pressure [Pa] |
velocity x [m/s] |
velocity y [m/s] |
velocity z [m/s] |
|
|
|
|
|
0.0000 |
13902.74 |
–6.87 |
0.00 |
15.36 |
|
|
|
|
|
0.0001 |
13738.71 |
–5.46 |
0.00 |
16.93 |
|
|
|
|
|
0.0002 |
12950.06 |
–3.80 |
0.00 |
17.74 |
0.0003 |
11572.66 |
–1.97 |
0.00 |
17.73 |
0.0004 |
9669.12 |
–0.05 |
0.00 |
16.92 |
|
|
|
|
|
0.0005 |
7325.99 |
1.87 |
0.00 |
15.34 |
|
|
|
|
|
0.0006 |
4649.79 |
3.71 |
0.00 |
13.07 |
0.0007 |
1762.19 |
5.38 |
0.00 |
10.19 |
0.0008 |
–1205.52 |
6.81 |
0.00 |
6.86 |
|
|
|
|
|
0.0009 |
–4118.43 |
7.92 |
0.00 |
3.21 |
|
|
|
|
|
0.0010 |
–6844.10 |
8.68 |
0.00 |
-0.58 |
|
|
|
|
|
The output of the first 51 time steps can be plotted as: