- •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
34 |
2 Derivation of Acoustic Wave Equation |
2.3Equation of State
The third equation needed for deriving the acoustic wave equation is called the equation of state which is the most difficult equation to derive and comprehend among the three equations. This is because it involves both the kinetic theory of gases and the principle of conservation of energy. Despite this difficulty, the physics involved in the equation of state is both interesting and important for understanding acoustics.
The equation of state describes the relationship between air pressure change and mass density change as shown in the figure below:
Pressure: |
|
|
|
Equation of State |
||||||
|
|
|
||||||||
Density: |
|
|
|
− |
= |
1 |
|
− |
||
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
: Absolute density
: Averaged density
: Absolute pressure
Pressure: |
↑ |
|
: Averaged pressure |
Density: |
↑ |
|
|
|
|
||
|
|
The equation of state can be expressed as:
P Po |
¼ |
γ |
ρ ρo |
Po |
|
ρo |
where:
Po is the average air pressure. P is the instantaneous pressure.
ρo is the averaged air mass density. ρ is the instantaneous mass density.
γis the ratio of specific heats and is defined as:
γ3α þ 2 3α
where α is a ratio between the total energy (rotational energy + translational energy) and the translational energy of molecules as:
2.3 Equation of State |
|
35 |
||
α |
|
21 Iaθca2 þ 21 mavca2 |
||
|
|
21 mavca2 |
|
or:
12 Iaθ2ca ¼ ðα 1Þ 12 mav2ca
where ma is the mass of a particle, Ica is the mass moment of inertia, vca is the collision speed, and θca is the spin speed of a particle.
The equation of state can be derived from (A) the law of conservation of energy and (B) pressure due to the kinetic theory of gases.
The following is a flowchart of the derivation of the equation of state:
Derivation of Equation of State
A 2 |
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
= |
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Energy |
|
|
|
Work |
|
|
|
|
|
|
|
|
|
|
|
|
|
Increase |
|
|
|
Done |
|
|
|
|
|
|
|
|||
|
|
|
|
|
= 0 |
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
||||||
= |
|
|
|
|
|
|
|
||||||||||
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
B
= |
1 |
= − |
3 |
Colliding Pressure |
the same |
|
= |
1 |
= − |
|
3 |
|
2.3.1Energy Increase due to Work Done
The energy of a particle with mass ma occupying a volume Va is the summation of translational and rotational energy:
Ka ¼ 12 mav2ca þ 12 Iaθ2ca
Assume the rotational kinetic energy is linearly related to the translational kinetic energy:
36 |
2 Derivation of Acoustic Wave Equation |
12 Iaθ2ca ¼ ðα 1Þ 12 mav2ca
where α is equal to or greater than 1. When the particle does not have rotational energy, α¼1; when the particle has rotational energy, α >1:
Kinetic Energy Considering |
|
|
|
|
|
|
|
Kinetic energy ( |
) |
|
|||||
Rotation of Air Particle |
|
||||||
|
1 |
|
1 |
|
̇ |
||
|
|
|
|
|
|
|
|
|
= 2 |
+ |
2 |
|
|||
|
|
|
=
̇
|
1 |
|
Colliding Pressure: |
||||
= |
|
|
1 |
|
|||
|
|
= |
|
||||
3 |
|||||||
|
|
|
|
||||
|
|
3 |
|
||||
|
|
|
|
|
Using this linear relationship between rotational energy and translational energy, the total energy of a particle is given as follows:
Ka ¼ 12 mav2ca þ 12 Iaθ2ca ¼ 12 mav2ca þ ðα 1Þ 12 mav2ca ¼ α2 mav2ca
According to the conservation of energy, in a closed system, the increase of energy is equal to the work done to the system:
∆ = |
|
|
|
= − |
2 |
|
|||
|
|
|
|
|
|
|
Energy |
|
Work |
|
|
Increase |
|
Done |
|
|
|
|
|
The total energy Ka includes both translational energy and rotational energy and is equal to α2 mav2ca , as shown before. And the work done by a constant external pressure Po due to a system of decreased volume dVa is PodVa; therefore:
Ka ¼ W
!α2 d mav2ca ¼ PodVa
2.3 Equation of State |
37 |
2.3.2Pressure due to Colliding of Gases
According to the kinetic theory, the pressure is directly related to the translational energy and is independent of the rotational energy. The pressure due to the particle collision is:
Po ¼ 1 ma v2
3 Va ca
The derivation of the collision pressure is based on the fundamentals of kinetic energy, as shown in the figure below:
Colliding Pressure (Kinetic theory of gases)
1 = 3
4 = 3
According to Newton’s law of motion, the force is equal to the change of momentum when the particle collides with the wall of the sphere as:
Fa |
|
ðmavcaÞ |
2mavca |
mavca2 |
||||
|
|
|
|
|
||||
¼ |
¼ vcaa |
¼ |
Ra |
|||||
|
t |
|||||||
|
|
|
|
2R |
|
|
|
The pressure on the surface of the sphere is the force divided by the total surface area of the sphere as:
2
Fa |
|
mavca |
|
1 ma 2 |
|||||||
|
|
Ra |
|
|
|||||||
Po ¼ |
|
¼ |
|
|
|
¼ |
|
|
|
|
vca |
Sa |
4πRa2 |
3 |
Va |
Note that the pressure is directly related to the translational energy only (not include the rotational energy), and rewrite the relation as: