- •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
312 11 Power Transmission in Pipelines
ZiPipe(nPipe)); |
|
% calculate Tw of the Main Pipe with a side-pipe-line |
|
ZiMain=loC/SiMain; |
|
ZoMain=loC/SoMain; |
|
ZoPara=1.0/(1.0/ZcPipe(1)+1.0/ZoMain); |
|
PoPipe_PiMain=2*ZoPara/(ZoPara+ZiMain); |
|
PoMain_PiMain=PoPipe_PiMain; |
|
PiMain=1; |
|
Tw(iw,2)=get_Tw_PipeLine(PiMain,ZiMain,PoMain_PiMain,ZoMain); |
|
end |
|
%--------------------------------------------------------------- |
– |
% |
|
% Section 4: Plotting |
|
%--------------------------------------------------------------- |
– |
% |
|
%plot the Tw of the pipe line figure(1); plot(freq,Tw(:,1)); xlabel('frequency [Hz]') ylabel('Tw')
title('Tw of a Pipe line') saveas(gcf,'Figure_Tw_of_a_Pipe_line','emf')
%plot the Tw of main pipe with a side branch of the pipe line figure(2);
plot(freq,Tw(:,2)); xlabel('frequency [Hz]') ylabel('Tw')
title('Tw of a Main Pipe with a Side Branch of a Pipe Line')
%saveas(gcf,'Figure_Tw_of_a_Main_Pipe_with_a_Side_Branch','emf') end
11.9Project
Any pipeline can be constructed by connecting n pipes as shown below. Acoustic waves enter the left-hand side of the pipeline and exit the right-hand side:
Pipe 1 |
|
Pipe 2 |
|
Pipe 3 |
|
Pipe |
|
|
|
|
|
|
|
= =
=
=
= 0 |
= |
= + |
11.10 Homework Exercises |
313 |
For given dimensions of the length L and the cross-sectional area S of each pipe, write a computer program using the provided MATLAB function (in Sect. 11.8) to calculate and plot power transmission coefficient Tw of the pipeline as a function of frequency.
Hint: Doing this example by hand can be difficult because (1) formulas of acoustic impedance are in complex number format and (2) power transmission coefficients need to be calculated at multiple frequencies. You can use the MATLAB functions provided in Computer Code Section (Sect. 11.7) or any suitable programming language for this project.
11.10Homework Exercises
Exercise 11.1
In Pipe 2 as shown in the figure below, a pressure pf2(x, t) of a forward wave and a pressure pb2(x, t) of a backward ware are formulated in terms of the complex amplitudes of the pressures Pf2R and Pb2R at x ¼ x2 (RHS of Pipe 2) as:
p f 2ðx, tÞ ¼ 12 nP f 2Re jkðx x2Þejωt þ cco
pb2ðx, tÞ ¼ 12 nPb2Rejkðx x2Þejωt þ cco
Pipe 2
Show that:
Part a (Backward Plane Waves)
The complex amplitude of the pressure Pb2L of the forward wave pb2(x, t) at x ¼ x2 (RHS of Pipe 2) can be formulated in terms of the complex amplitude of the pressure Pb2R and length L2 as:
Pb2L ¼ Pb2Re jkL2 |
ð(Formula 2BÞ |
314 |
11 Power Transmission in Pipelines |
Part b (LHS of a Pipe)
The acoustic impedance Zc2L at x ¼ x1 (LHS of Pipe 2) can be formulated in terms of the complex amplitudes of the pressures Pf2R, Pb2R, the length L2, and the crosssectional area S2 as:
Z |
c2L ¼ |
Z |
f 2 |
P f 2RejkL2 |
þ Pb2Re jkL2 |
ð |
Formula 3A |
Þ |
|
P f 2RejkL2 |
Pb2Re jkL2 |
||||||||
|
|
|
where Z f 2 ¼ ρoc.
S2
Exercise 11.2
The cross-sectional areas of Pipe 1 and Pipe 2 are S1 and S2, respectively, as shown in the figure below. The pressures of the forward wave and backward wave in Pipe 1 at x ¼ x2 (RHS of Pipe 1) are Pf1R and Pb1R, respectively. The pressures of the forward wave and backward wave in Pipe 2 at x ¼ x3 (RHS of Pipe 2) are Pf2R and Pb2R, respectively:
Pipe 1 |
|
Pipe 2 |
|
|
|
|
, |
Show that the pressures Pf1R and Pb1R (RHS of Pipe 1) can be formulated in terms of the pressure Pf2R and Pb2R (RHS in Pipe 2) as:
P f 1R ¼ 2 |
1 þ Z f 2 |
P f 2RejkL2 |
þ 2 |
1 Z f 2 |
Pb2Re jkL2 |
|||
1 |
|
Z f 1 |
|
1 |
Z f 1 |
|
||
Pb1R ¼ 2 |
1 Z f 2 |
P f 2RejkL2 |
þ 2 |
1 þ Z f 2 |
Pb2Re jkL2 |
|||
1 |
|
Z f 1 |
|
1 |
|
Z f 1 |
|
Exercise 11.3 (Closed-Open Pipe)
A pipe has a length L and a cross-section area Si. The far end is open (Zo ¼ 0) and the near end is closed by a rigid cap (Z1 = 1):
11.10 Homework Exercises |
315 |
Near end |
Pipe 1 |
Far end |
Pipe 0 |
|||||
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
= ∞ |
|
|
|
= 0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(a)Calculate the eigenfrequencies of the pipe in terms of the pipe length L.
(b)Calculate the eigenmodes of the pressure and plot the first three modes of the pipe using the equation below:
Z |
|
= Z |
Zc0L þ jZ f 1 tan ðkLÞ |
ð |
Formula 5 |
Þ |
||||||
c1L |
f 1 jZc0L tan ðkLÞ þ Z f 1 |
|||||||||||
|
|
|
|
|
|
|||||||
where: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ρoc |
|
|
|
|
|||
|
|
|
|
|
Z f 1 ¼ |
|
|
|
|
|
|
|
|
|
|
|
|
S1 |
|
|
xi, n ¼ |
||||
(Answers): (a) ωn ¼ c |
ð |
2L |
, n ¼ 0, 1, 2. . .; (b) pn ¼ 2jP sin h |
2L |
||||||||
|
2nþ1Þπ |
|
|
|
|
ð2nþ1Þπ |
|
|
||||
0, 1, 2 . . . |
|
|
|
|
|
|
|
|
|
|
|
Chapter 12
Filters and Resonators
In the previous chapter, we defined and derived formulas of complex amplitudes, acoustic impedance, and power transmission in the pipeline. In this chapter, we will extend the analysis of pipelines (pipes without side branches) to pipeline with side branches. Most of the formulas derived in the previous chapter will be reused in this chapter. In fact, a pipeline with a side branch can be treated as a one-to-two pipeline. The acoustic analysis of a one-to-two pipeline can be treated as a one-to-one pipeline by an equivalent acoustic impedance for the two pipes.
Three basic types of filters will be introduced in the second half of this chapter: low-pass filters (Sect. 12.3), high-pass filters (Sect. 12.4), and band-stop filters (Sect. 12.5) based on the principle of the Helmholtz resonator. This chapter will introduce the fundamental principles, the formulas, and the power transmission coefficient of each type of filter.
A numerical method for modeling pipelines with side branches is included in this chapter. A numerical project to model filters is provided without a solution. This chapter is organized into the following sections:
Section 12.1 introduces one-to-two pipes. Formulas for calculating the equivalent acoustic impedance of a one-to-two pipe are derived (Formulas 8A–8B) in this section.
Section 12.2 introduces the power transmission of one-to-two pipes. Formulas for calculating power reflection and transmission will be derived (Formulas 13A–13B) in this section.
Section 12.3 introduces low-pass filters. A formula for the power transmission coefficient of an open side branch as a low-pass filter will be derived in this section.
Section 12.4 introduces high-pass filters. A formula for the power transmission coefficient of a three pipeline as a high-pass filter will be derived in this section.
Section 12.5 introduces band-stop filters. A formula for the power transmission coefficient of the Helmholtz resonator as a band-stop filter is derived in this section.
Section 12.6 introduces a numerical method for calculating the power transmission coefficients for filter designs. This method can be used to model all three basic types of filters: high-pass filters, low-pass filters, and band-stop filters.
© The Author(s), under exclusive license to Springer Nature Switzerland AG 2021 |
317 |
H. Lin et al., Lecture Notes on Acoustics and Noise Control, https://doi.org/10.1007/978-3-030-88213-6_12
318 |
12 Filters and Resonators |
12.1Pressure in a One-to-Two Pipe
The power transmission coefficients of a one-to-two pipe as shown below will be calculated and analyzed in this section.
The dimensions and acoustic impedance of a one-to-two pipe are given as:
• Cross-section areas of the pipes are Si, S1, and S2.
• Acoustic impedances of the pipes are Zi, Zr,Z1, and Z2:
1 = 1 + 1
=
= =
+
where:
Si: cross-section area of the inlet pipe
S1: cross-section area of outlet Pipe 1 (side branch) S2: cross-section area of outlet Pipe 2 (main pipe) pi: pressure of the incident wave in the inlet pipe
p1: pressure of the transmitted wave in outlet Pipe 1 (side branch) p2: pressure of transmitted wave to outlet Pipe 2 (main branch)
Pi ¼ PiR: complex pressure of the incident wave at the RHS of the inlet pipe
P1: complex pressure of the transmitted wave in the side branch pipe at the intersection
P2 ¼ P2L: complex pressure of the transmitted wave at the LHS of the outlet pipe Zi: acoustic impedance of the incident wave at the RHS of the inlet pipe
Z1: acoustic impedance of the transmitted wave inside the branch pipe at the intersection
Z2: acoustic impedance of the transmitted wave at the LHS of the outlet pipe
Pressures and Velocities in a One-to-Two Pipe
The incident and transmission pressures are forward plane waves (Pi, P1, and P2) and have the same wave function format. The pressure and velocity functions of an incident wave are shown below for reference: