Yung Boon Chong PWE
Yung Boon Chong PWE
Yung Boon Chong PWE
by
Yung Boon Chong, BEng
in the
Department of Design, Development,
Environment and Materials
Faculty of Mathematics, Computing and Technology
The Open University
Milton Keynes
United Kingdom
November 2012
Certificate of Originality
This is to certify that I am responsible for the work
submitted in this thesis, that the original work is my
own except as specified in acknowledgements, in text
or in bibliography, and that neither the thesis nor the
original work contained therein has been submitted to
this or any other institution for a degree.
_____________________________ (Signature)
(Y.B. Chong)
_____________________________ (Date)
Acknowledgements
This dissertation would not have been possible without the guidance and the
help of several individuals who in one way or another contributed and
extended their valuable assistance in the preparation and completion of this
study.
This research has been undertaken with the support of Engineering and
Physical Science Research Council, United Kingdom (Grant number
EP/E062806/1) which is gratefully acknowledged. I should also mention with
gratitude that my graduate studies were supported by the University as part of
the employee subsidy scheme. I could not fulfil my dream to pursue this
higher degree without these sponsors. Sincere appreciations are given to our
project collaborators Dr. Olga Umnova and Dr. Anton Krynkin from the
University of Salford for providing their expertise and the many productive and
challenging discussions. I also wish to acknowledge the help from Dr. Anton
Krynkin for providing and guiding the use of the MST algorithms. Special
thanks are due to Dr. Juan V. Sanchez-Perez and Dr. Vincent Romero Garcia
from the Polytechnic University of Valencia for collaboration in several
aspects of the work related to this Thesis, notably in respect of the influence
of ground effect on sonic crystal performance.
The assistance given by the staff and fellow research student in the Acoustic
Research Group throughout the research periods are acknowledged and
appreciated. I am indebted to Dr. Roland Kruse, Dr. Toby Hill, Dr. Ho-chul
Shin and Dr. Adrien Mamou-mani for their constructive advice on my
research. Special thanks to Peter Seabrook (Project Officer) who contributed
much time and effort to source and fabricated specimens to be used for the
sonic crystal both in laboratory and outdoor measurements. I should mention
with gratitude to Roger Frith and Mikki Thomas (Project Officers) for their help
rendered during the laboratory measurements. Credit also goes to Stan Hiller
(Project Officer) from the Material’s Department for his expertise in performing
tensile strength testing for the latex specimen used in our design. I would also
like to thank Imran Bashir (PhD student) for his insightful suggestions,
discussions and help rendered.
Last, but not least, I am truly grateful to my family whose care and support
have been the important factors which let me seek this higher degree without
worries. This dissertation is dedicated to them.
I
Preface
This PhD Thesis entitled “Sonic Crystal Noise Barriers” contains the results of
Taherzadeh and Dr. David Brian Sharp who are all members of academic
staff at the Open University. This research was funded by the Engineering
collaboration with Dr. Olga Umnova and Dr. Anton Krynkin from University of
(MST) and computational model (rigid, elastic shell, composite scatterers and
sonic crystal above a ground surface) have been used extensively in this
Method (FEM) and experimental works are developed and performed at the
Open University and Diglis Weir, Worcester. The sonic crystals work on Diglis
Weir arose from an opportunity due to the “Organ of Corti” project, (see
details in Section 5.6). Based on the outcome of the research findings, there
II
List of Publications
International journals
Conference proceedings
III
Contents
Acknowledgments ------------------------------------------------ I
Preface ------------------------------------------------------------------ II
Contents ---------------------------------------------------------------- IV
Chapter 1: Introduction
1.1) Fundamentals of acoustics ------------------------------- 1
1.2) Complex number notation ------------------------------- 6
1.3) The conventional road traffic noise barrier ---------- 7
1.4) Sonic Crystals ----------------------------------------------- 9
1.4.1) History ----------------------------------------------- 9
1.4.2) Crystallography -------------------------------------- 13
1.5) Aims and thesis organisation ----------------------------- 23
IV
Chapter 3: Multiple Scattering Theory (MST)
3.1) Introduction ---------------------------------------------------- 43
3.2) Multipole method for circular scatterers in 2- ---------- 47
dimensional system
3.3) Plane wave scattering model ----------------------------- 51
3.4) Cylindrical wave scattering model ----------------------- 56
3.5) Simulation and results -------------------------------------- 58
V
5.6) Outdoor in situ measurements of a sonic -------------- 123
crystal at Diglis Weir, Worcester
5.6.1) Measurement arrangement --------------------- 125
5.6.2) Measurement results ------------------------------ 127
5.6.3) Modelling with Multiple Scattering Theory -- 128
(MST) assuming a single point source
5.6.4) Modelling with multiple point sources -------- 132
on a line
VI
7.6) Coupled resonating elements ----------------------------- 202
7.6.1) Semi-analytical formulation of single -------- 204
composite scatterer
7.6.2) FEM modelling of single composite ---------- 212
scatterer
7.6.3) Results for single composite scatterer ------ 213
7.6.4) Semi-analytical formulation of an array ----- 220
of composite cylinders
7.6.5) Results for array of composite scatterers -- 222
7.7) Influence of adjustable parameters on the ------------- 229
acoustical performance of composite scatterer
arrays outdoors
7.7.1) The predicted influence of the elastic ------- 230
shell outer diameter
7.7.2) The predicted influence of elastic shell ----- 231
wall thickness
7.7.3) Predicted influence of the slit widths --------- 232
in the outer PVC cylinder
VII
8.6) Laboratories studies on rigid and elastic cylinder --------- 271
arrays with their axes parallel to a finite
impedance ground and normal to a rigid ground
8.6.1) IL spectra for horizontal rigid cylinders ------- 271
array over a finite impedance ground
8.6.2) IL spectra for horizontal elastic shell ---------- 273
cylinders array over a finite impedance
ground
8.6.3) Laboratories studies on a vertical axis -------- 274
cylinder array over a horizontal acoustically
rigid ground (MDF 1)
8.7) Laboratories studies on vertical axis array over ------ 276
acoustically rigid ground (glass plate)
Appendix C
i) Experimental study of the use of Polyethylene ------ XXXV
scatterers in a sonic crystal noise barrier
ii) Experimental study of the use of periodic ------------- XXXVII
stiffening along the scatterer length
VIII
List of Figures
IX
2.2 The left hand panel shows the plane wave expansion ------- 34
computed band structure for a SC consisting of
cylindrical PVC scatterers embedded in air in a
square array with lattice constant L of 0.069 m and
filling fraction of 50%. The right hand panel shows
the laboratory measured IL spectrum for 7x3 array
for the [1, 0] direction. Inset: schematic of a
Brillouin zone. refers to [1,0] direction, refers
to [1,1] direction and refers to the wave vector
varying from [1,0] to [1,1]. The (red) dotted box
indicates a full bandgap and the (blue) shaded
regions indicate partial bandgap.
2.3 Plane wave expansion predictions of the effect of ------- 35
varying the lattice constant in a SC comprising of
0.055 m diameter PVC cylinders embedded in air.
Lattice constant L of (a) 0.08 m and (b) 0.135 m.
Left hand panels show the band structure of the
corresponding SC and right hand panels show the
measured IL for the different lattice parameters.
2.4 Plane Wave Expansion predictions of the band ------- 38
structure with Lattice constant fixed at 0.135 m and
varying scatterer diameters (a) 0.055, (b) 0.09, (c)
0.11 and (d) 0.13 m.
2.5 Plane Wave Expansion predictions of the influence of ------- 41
material parameters on the band structure. The
lattice constant is fixed at 0.069 m and scatterers are
made from different materials (a) PVC, (b) silicone
rubber, (c) steel (d) PMMA and (e) wood. (f) PVC
scatterers embedded in water.
3.1 Plan view of two identical cylinders and ------- 48
corresponding cartesian and polar coordinates.
3.2 Example scenario including an array of rigid ------- 59
scatterers.
3.3 MST predictions and measured Insertion Loss ------- 60
spectra for square lattice arrays of rigid circular
scatterers of 0.055 m diameter with lattice constants
of (a) 0.069, (b) 0.08 and (c) 0.135 m respectively.
Both plane and cylindrical waves are compared for
all three cases.
3.4 MST predictions and measured Insertion Loss ------- 61
spectra for square lattice arrays of rigid circular
scatterers of 0.055 m diameter with lattice constants
of 0.069. Different truncation number of M 1, 2, ... , 5
is used for the MST predictions for (a) plane wave
and (b) cylindrical wave.
4.1 Definition of domain discretized using triangular ------- 66
shape finite elements.
4.2 Definition of geometry for a rectangular domain ------- 69
modelled as air and an incoming cylindrical wave
source.
X
4.3 COMSOL® computed pressure map for a rectangular ------- 70
domain modelled as air and an incoming plane (a) or
cylindrical wave (b) source.
4.4 COMSOL® computed pressure maps for 7x3 array of ------- 71
sonic crystal (acoustically hard scatterer) modelled in
rectangular air domain. Cylindrical waves is
performed and pressure maps at 200 Hz, 1.2, 3, 4,
4.5 and 5 kHz are shown for figure (a), (b), (c), (d)
and (e) respectively.
4.5 FEM predictions and laboratory measured Insertion ------- 72
Loss spectra for 7x3 square lattice arrays of rigid
scatterers of radius 0.0275 m with lattice constant
0.135 m. FEM predictions are performed with plane
and cylindrical wave sources for a rectangular
domain modelled as air.
4.6 FEM predicted pressure plots for cylindrical waves ------- 75
with different mesh element sizes.
4.7 Sound pressure level spectra at single point position ------- 76
for FEM computations of the field due to a cylindrical
wave computed using different mesh element sizes.
4.8 FEM (COMSOL®) predicted pressure maps ------- 78
comparing three different arrangements (top to
bottom) and 2 frequencies pressure maps for each
arrangement (left 1.3 kHz and right 2.3 kHz) of 7x3
triangular rigid scatterer arrays (square lattice) with
lattice constant of 0.135 m.
4.9 FEM predictions of Insertion Loss spectra for 7x3 ------- 79
square lattice arrays of triangular rigid scatterers with
lattice constant 0.135 m. A cylindrical wave source is
assumed and a rectangular domain is modelled as
air.
4.10 FEM (COMSOL®) predicted pressure maps ------- 80
comparing two 7x3 square scatterer arrangements
(square lattice with lattice constant of 0.135 m) at 4.3
kHz.
4.11 FEM predictions of Insertion Loss spectra for two ------- 80
arrangements of 7x3 square lattice arrays of square
rigid scatterers and a reference array of cylindrical
rigid scatterers with lattice constant 0.135 m. A
cylindrical wave source is assumed and a
rectangular domain is modelled as air.
4.12 FEM (COMSOL®) predicted pressure maps for three ------- 82
7x3 arrangements of elliptical rigid scatterers (square
lattice) with lattice constant of 0.135 m at 1.5 kHz.
4.13 FEM predicted Insertion Loss spectra for 7x3 square ------- 82
lattice arrays of elliptical rigid scatterers with lattice
constant 0.135 m compared with that predicted for an
equivalent cylindrical scatterer array. For the FEM
predictions a cylindrical wave source is assumed and
the rectangular domain is modelled as air.
XI
4.14 FEM predictions of Insertion Loss spectra for 7x3 ------- 83
square lattice arrays of the best performing
triangular, square and elliptical rigid scatterer arrays
compared to that predicted for the reference
cylindrical scatterer array.
4.15 FEM (COMSOL®) predicted pressure maps at 1.5 ------- 84
kHz for 7x3 arrays of cylindrical rigid scatterers with
lattice constant of 0.135 m (a) hexagonal lattice and
(b) square lattice.
4.16 FEM predictions of Insertion Loss spectra for 7x3 ------- 85
hexagonal and square arrays of cylindrical rigid
scatterers with lattice constant 0.135 m. FEM
prediction are performed with cylindrical wave source
and a rectangular domain modelled as air. The
cylindrical wave source is located at 1.5 m away from
the front face of the array and the receiver is
positioned 0.05 m from the nearest face of the array
but on the opposite side of the source.
4.17 Implementation of PMLs around a 7x3 square lattice ------- 87
array of circular scatterers (a) location of PMLs and
(b) FEM-computed pressure map at 1.2 kHz.
4.18 FEM predictions (with and without PML) compared to ------- 87
laboratory measurements of Insertion Loss spectra
for a 7x3 square lattice array of circular rigid
scatterers with a lattice constant 0.135 m.
5.1 Flow diagram for a typical impulse response ------- 93
measurement.
5.2 Linear feedback shift register for generation of a MLS ------- 94
(of length 24 - 1= 15) signal.
5.3 MLSSA system with audio patch cable connecting ------- 95
analogue output to input.
5.4 System loopback test - MLS time signal. ------- 96
5.5 Transfer function for system loopback test - MLS ------- 96
power spectrum.
5.6 Laboratory measuring system incorporating MLSSA. ------- 99
5.7 (a) Supporting base plate for 7x3 arrays of PVC or ------- 100
elastic cylinders. 55 mm diameter holes for PVC
cylinders. 5 mm diameter holes for elastic shell using
bolt and nut securing system. (b) Lower end of an
elastic shell showing plastic pipe for air inlet and
mounting bolt.
5.8 (a) Plan view of the source, receiver and array in the ------- 101
laboratory measurements at normal incidence (b) the
corresponding side view. Refer to chapters 2, 3, 4, 6
and 7 for the outer diameters and lattice constants
used.
5.9 MLSSA time signals for direct and total transmitted ------- 102
fields.
5.10 MLSSA Frequency spectra for direct and total ------- 103
transmitted fields.
XII
5.11 Insertion Loss (IL) spectrum using MLSSA of 7x3 ------- 103
square lattice arrays of rigid PVC scatterers of outer
diameter 0.055 m and lattice constant 0.069 m.
Frequency resolution at 48 Hz.
5.12 Power spectra for different window size taken in ------- 104
calculating the FFT for total transmitted fields.
5.13 Insertion loss spectra for 7x3 array of PVC scatterers ------- 106
of outer diameter 0.055 m and lattice constant of
0.069 m using different FFT filters.
5.14 Insertion loss spectra for single latex scatterer. ------- 107
5.15 Insertion loss spectra for single PVC scatterer of ------- 108
outer diameter 0.11 m, wall thickness of 0.003 m with
4 symmetrical slits size of 0.012 m.
5.16 Arrangement for measuring Excess attenuation ------- 109
spectra above a Medium Density Fibreboard (MDF)
in the laboratory.
5.17 Drilled Metal panels used for securing the scatterers. ------- 111
5.18 The outdoor noise barrier site, (a) source side of ------- 112
barrier showing SC and conventional noise barrier
sections and (b) rear side of barrier.
5.19 Plan view of the outdoor measurement arrangement ------- 112
using simultaneous multiple channels.
5.20 Graphical User Interface (GUI) Traffic Noise Analyzer ------- 113
implemented in Matlab®.
5.21 Example time domain signals measured outdoors (a) ------- 115
during SC barrier measurements using B&K
microphones and (b) during fence measurements
using ACO microphones.
5.22 Example power spectra measured outdoors (a) ------- 116
during SC barrier measurements using B&K
microphones and (b) during fence measurements
using ACO microphones.
5.23 Example background noise power spectra measured ------- 117
at the barrier test site (a) during SC barrier
measurements using B&K microphones and (b)
during fence measurements using ACO
microphones.
5.24 Setup for B&K microphone calibration. ------- 118
5.25 Calibration results for B&K microphones. ------- 119
5.26 Sonic anemometer mounted above the ------- 121
instrumentation cabin.
5.27 Data from sonic anemometer. ------- 121
5.28 Example of a deduced wind vector diagram. ------- 122
5.29 Swept sine signal generated electrically. ------- 123
5.30 Direct and total transmitted fields measured using the ------- 123
swept sine method in the laboratory.
5.31 Comparison of IL spectra due to a single PVC ------- 124
scatterer of OD 0.11 m, WT 0.003 m with 4
symmetrical slits of 0.012 m measured using the
MLS and swept sine methods.
XIII
5.32 (a) Plan view schematic of the sonic crystal ------- 126
arrangement and microphone locations at Diglis weir,
Worcester. (b) Aerial map of the site showing where
the sonic crystal was situated (picture taken from
Google map).
5.33 Mean IL spectrum (average of 3 measurements) with ------- 127
error bars for receiver position Pt 1.
5.34 The original receiver position Pt 1 and 4 other ------- 128
perturbed positions made for MST modelling.
5.35 Predicted IL spectra of all the 5 individual positions ------- 129
(Pt 1 original and 4 perturbed positions) and the
averaged IL spectrum. Source is located at
coordinates (0, 0).
5.36 Measured (averaged of 3 measurements) and ------- 130
predicted (averaged of 5 receiver points) IL spectra
for receiver position Pt 1 with source located at
coordinates (0, 0).
5.37 Predicted IL spectra at all the 5 individual positions ------- 131
(Pt 1 original and 4 perturbed positions) and the
averaged predicted IL spectrum. The point source is
assumed to be located at coordinates (0, 10).
5.38 Measured (averaged of 3 measurements) and ------- 132
predicted (averaged of 5 points) IL spectra for
receiver position Pt 1 with source located at
coordinates (0, 10).
5.39 Picture illustration for the location of the 10 sources ------- 134
(S1 to S10) and the position of sonic crystal.
5.40 Measured (average of 3 measurements) and MST ------- 135
predicted with multiple point sources IL spectra for
receiver position Pt 1 (averaged over 5 receiver
points for each individual source).
6.1 Example predictions of the two basic mode shapes of ------- 141
a freely supported cylinder – (a) radial (b) beam
modes up to fourth order (c) the combination of the
first radial mode (mode 0) and up to the fourth order
of bending modes (d) the combination of the first
radial mode (mode 0) and up to the third order of
axial modes [103].
6.2 Cross-section of an elastic shell in the primary cell of ------- 142
doubly periodic array.
6.3 Measured (solid black line) and MST predicted ------- 154
(broken blue line) IL spectra for single latex scatterer
of outer diameter (OD) 0.055 m, wall thickness (WT)
0.00025 m and length of 2 m. For comparison the IL
spectrum measured for an acoustically “rigid” Poly
Vinyl Chloride (PVC) pipe with similar diameter and
length is shown also (broken red line).
6.4 2-Dimensional modal analysis showing static ------- 157
deformation plots in sequence (a-f) on the “Breathing
mode” shape of elastic element in vacuum.
XIV
6.5 2-Dimensional modal analysis showing static ------- 159
deformation plots in sequence (a-d) on the
“Breathing mode” shape of elastic element coupled
with air surrounded by a perfectly matched layers for
the domain boundaries.
6.6 IL spectra of measurements and MST predictions for ------- 161
different outer diameter 0.0275, 0.043 and 0.055 m in
figure (a), (b) and (c) respectively.
6.7 Measured and predicted IL spectra for outer diameter ------- 163
of 0.043 m (left panel) and 0.055 m (right panel) with
two wall thicknesses of 0.00025 ((a) and (b)) and
0.0005 m ((c) and (d)) in each case.
6.8 IL spectra of measurement and MST predictions of ------- 164
latex geometry of 0.043 m and wall thickness of
0.0005 m. MST predictions are made for the 3
different stiffness of latex (refer to table 6.2).
6.9 Measured IL spectra of single latex of OD 0.0275 m ------- 166
and WT 0.00025 m with different tensioning states
(tightly, normal and loosely tension).
6.10 Measured and MST predicted IL spectra for an array ------- 168
of latex shell scatterers of outer diameter (OD) 0.055
m and wall thickness (WT) 0.00025 m with lattice
constant 0.08 m.
6.11 Measured and MST predicted IL spectra for an array ------- 169
with lattice constant 0.1 m of latex shell scatterers
having outer diameter (OD) 0.055 m and wall
thickness (WT) 0.00025 m.
6.12 Measured and MST predicted IL spectra for an array ------- 169
with lattice constant 0.15 m of latex shell scatterers
having outer diameter (OD) 0.055 m and wall
thickness (WT) 0.00025 m.
6.13 Figure 6.13: (a) Plan view of the source, receiver and ------- xxx
array in the laboratory measurements at normal
incidence and microphone laterally displace (see red
arrow) up to 14 cm (b) the corresponding side view
with microphone vertically displace (see red arrow)
up to 45 from normal incidence.
6.14 Lateral displacement, measurement IL spectra for ------- 170
array with lattice constant 0.08 m of latex scatterers
having outer diameter (OD) 0.055 m and wall
thickness (WT) 0.00025 m.
6.15 Lateral displacement, measurement IL spectra for ------- 171
array with lattice constant 0.15 m of latex scatterers
having outer diameter (OD) 0.055 m and wall
thickness (WT) 0.00025 m.
6.16 Meshing of single latex with OD 0.055 m and wall ------- 173
thickness (WT) 0.00025 m (a). FEM predicted
pressure plot at 1050 Hz for the exterior environment
of the model (b).
XV
6.17 Measured and predicted (FEM and MST) IL spectra ------- 173
for a single latex scatterer having 0.055 m (OD) and
wall thickness 0.00025 m (WT).
6.18 Measurement setup for array of cylinders of 7x3 ------- 175
array of latex (a) and a single latex showing the
overlapping strip, mounting bolt and pipe for air inlet
(b).
6.19 Averaged of 21 measurements of IL with error bars ------- 177
(indicating the range of data at each frequency) and
MST predicted IL spectra for a single latex scatterer
of (outer diameter (OD) 0.055 m and wall thickness
(WT) 0.00025 m) (blue broken line) and (OD 0.057 m
and WT 0.00025 m) (red solid line).
6.20 2-Dimensional FE modelling of the effect on the ------- 178
glued overlapping section (section blown up for
clarity).
6.21 Measured and FEM-prediction IL spectra due to a ------- 179
single latex cylinder (with and without rigid inclusion)
of outer diameter (OD) 0.055 m and wall thickness
(WT) 0.00025 m.
6.22 Measured and FEM-prediction IL spectra due to a ------- 179
single latex cylinder (with and without rigid inclusion)
of outer diameter (OD) 0.0275 m and wall thickness
(WT) 0.00025 m.
6.23 Industrially-manufactured latex cylinders joined by a ------- 181
Perspex disc (“Rigid”) and mounted between
supports in the anechoic chamber. The insert on the
left shows the shape of an “Extreme flat region” with
averaged dimensions and the insert on the right
shows an “Inflated region” with averaged dimensions.
6.24 IL spectra of single industrial latex (Sample 1, S1) at ------- 181
all 5 positions (Extreme flat face/edge, rigid and
inflated face/edge).
6.25 IL spectra of consistency test for 2-single industrial ------- 182
latex (Sample 1 (S1) and 2 (S2)) at all 5 positions
(Extreme flat face (a) or edge (b), inflated face (c) or
edge (d) and rigid (e)).
6.26 Measured IL spectra due to 7x2 and 7x3 arrays of ------- 183
industrial latex cylinders and an MST predicted IL
spectrum for a 7x3 latex array. The insert on the right
shows a photograph of the industrial latex array.
7.1 Illustration of Helmholtz resonator and vibration ------- 187
absorber.
7.2 COMSOL® pressure plot at 410 Hz of a 2D Helmholtz ------- 189
resonator.
7.3 Example slotted pipe SRR. ------- 190
7.4 Example of laboratory set up for array of SRR. ------- 191
7.5 Measurements and FEM predictions (of IL spectra for ------- 192
the single cylinders of type (a) NRR1, (b) SRR1, (c)
NRR2 and (d) SRR2 (see table 7.1).
XVI
7.6 Pressure maps predicted using COMSOL® for ------- 193
incident cylindrical waves on (a) a single scatterer
SRR1 at 700 Hz and (b) a single scatterer SRR2 at
300 Hz.
7.7 Single SRR orientation. ------- 194
7.8 Measured Insertion Loss spectra of the radiation ------- 195
pattern of single SRR (refer to figure 7.7 for the
definition of the orientated angle).
7.9 FEM predictions of IL spectra for the single SRR1 ------- 195
with increasing slit width.
7.10 Predicted and measured Insertion Loss spectra due ------- 197
to square arrays of SRR1 (refer to table 7.1) with
lattice constants of (a) 135, (b) 80 and (c) 70 mm.
Also shown in the left hand panels are pressure
maps at the first Bragg diffraction frequency for the
corresponding array geometry.
7.11 Example outdoor measurement arrangement (a) ------- 199
lateral displacements of the microphone along the
length of barrier. (b) Vertical displacement of the
microphone along the height of barrier.
7.12 Predicted and measured Insertion Loss spectra at ------- 199
outdoor for square arrays of reference no slit
cylinders (NRR2) and SRR scatterers (SRR2) with
lattice constants, L, of 160 mm (refer to table 7.1). (a)
54x3 square lattice array of scatterer NRR2. (b) 54x3
square array of scatterer SRR2. Corresponding
pressure maps at 500 Hz are shown in the left-hand
panels.
7.13 Measured Insertion Loss spectra of various lateral ------- 200
angles at outdoor for square lattice arrays of SRR2
(see table 7.1) with lattice constant of 0.16 m.
7.14 Measured Insertion Loss spectra due to outdoor ------- 201
square arrays of split ring resonators, SRR2 (see
table 7.1) with lattice constant of 0.16 m at various
inclinations of the source-receiver axis.
7.15 (a) Cross section of composite element consisting of ------- 204
a concentric arrangement of an outer 4-slit rigid
cylinder and an inner elastic cylindrical shell. (b)
Geometry of the slit.
7.16 (a) Measured and FEM predicted IL spectra due to a ------- 213
single PVC cylinder with 4 symmetrical slits
(component of SRR3, see table 7.2). (b) FEM
(COMSOL®) predictions of pressure maps at 2 kHz
for a similar scatterer showing the interior and
exterior environments.
7.17 Measured and predicted (MST and FEM) IL spectra ------- 213
for a single composite scatterer, SRR3 (refer to table
7.2).
XVII
7.18 FEM (COMSOL®) predicted pressure maps at 1 kHz ------- 215
for SRR3 showing the exterior and annular cavity
environments (a), interior environment of the latex
cavity (b) and static deformation plot on the
“Breathing mode” shape of elastic element in air.
7.19 Measured and predicted (MST and FEM) IL spectra ------- 216
for a single composite scatterer, SRR4 (refer to table
7.2).
7.20 Measured IL spectra for a single composite scatterer, ------- 216
SRR3 and it’s own component (refer to table 7.2).
7.21 Measured (in laboratory) IL spectra for a single ------- 217
concentric cylinder (SRR4) and its components (i.e.
single 4 slit PVC and latex, see table 7.2).
7.22 (a) Measured and MST predicted IL spectra for 7x3 ------- 220
PVC array of lattice constant 0.08 m (PVC with 4
symmetrical slits, component of SRR3, see table
7.2). (b) FEM (COMSOL®) predicted pressure maps
at 2 kHz for a similar scatterer array showing the
interior and exterior environments.
7.23 Measured and FEM predicted IL spectra for 7x3 ------- 221
array of composite scatterers SRR3 (see table 7.2)
arranged in square lattice of lattice constant 0.08 m.
7.24 Measured IL spectra due to 7x3 arrays of composite ------- 222
scatterer, SRR3 and the two components separately
(refer to table 7.2).
7.25 Average of 8-day measurements (outdoor) and MST ------- 223
IL spectra for 54x3 arrays of composite scatterers
SRR4 (see table 7.2) arranged in square lattice of
lattice constant 0.16 m. Error bars are plotted for the
measured spectrum.
7.26 Measured IL spectral of 36x4 and 54x3 arrays of ------- 224
concentric cylinders (averaged of all 8
measurements).
7.27 Measured IL spectral of 54x3 arrays of concentric ------- 224
cylinders (averaged of 8 measurements) and rigid no
slit cylinders (average of 3 measurements).
7.28 Measured IL spectral of 54x3 arrays of concentric ------- 225
cylinders (averaged of 8 measurements) and fence
(average of 8 measurements).
7.29 IL spectra for the MST predicted effect on changing ------- 226
latex outer diameter for 54x3 square lattice array of
SRR4 scatterer with lattice constant of 0.16 m.
7.30 IL spectra for the MST predicted effect on changing ------- 227
latex wall thickness for 54x3 square lattice array of
SRR4 scatterer with lattice constant of 0.16 m.
7.31 MST predicted IL spectra for 54x3 square lattice ------- 228
array with lattice constant of 0.16 m as the slit widths
in the outer PVC cylinders are varied from 3 to 12
mm.
XVIII
8.1 Measurement configuration to determine the surface ------- 242
impedance. (a) Single microphone method
(laboratory) and (b) Transfer function method
(outdoor in situ).
8.2 Dimensions for the range and heights of source and ------- 243
receiver for each Excess Attenuation (EA)
measurement in laboratory.
8.3 Measured and predicted EA spectra for ground (MDF ------- 244
1) at source and receiver heights at (a) 0.10 m (b)
0.20 m and (c) 0.30 m (see Tables 8.1 and 8.2).
8.4 Measured and predicted EA spectra over MDF 2 ------- 246
ground with source height fixed at 0.23 m and
receiver heights of (a) 0.11 m (b) 0.23 m and (c) 0.35
m (see tables 8.1 and 8.2).
8.5 Measured and predicted EA spectra over ground ------- 248
(glass) with source and receiver heights at (a) 0.10 m
(b) 0.15 m and (c) 0.20 m respectively (see tables
8.1 and 8.2).
8.6 Measured and predicted EA spectra over ground ------- 250
(polyurethane foam) with source height fixed at 0.23
m and receiver heights of (a) 0.11 m (b) 0.23 m and
(c) 0.35 m (see tables 8.1 and 8.2).
8.7 Measured and predicted LD spectra over MDF 2 with ------- 251
source height fixed at 0.07 m and receiver heights of
0.07 m and 0.14 m (see tables 8.1 and 8.2).
8.8 Measured and predicted LD spectra over asphalt ------- 253
surface with measurement geometry 2 (table 8.3).
8.9 Measured and predicted LD spectra over grass ------- 254
surface with measurement geometry B (table 8.3).
8.10 Square lattice array above a perfectly reflecting ------- 255
plane. (a) Set of vectors used in equation (8.6). (b)
Set of vectors employed in equation (8.14).
8.11 (a) Schematic of the experimental configuration ------- 261
showing the source location, the sonic crystal array
and three receiver locations. (b) Photograph of
experimental arrangement with rigid cylinders above
rigid (MDF 2 board) ground. (c) Photograph of
experimental arrangement with latex shell array near
to rigid (MDF 1 board) ground (see table 8.1).
8.12 Measured (solid line) and predicted (broken line) IL ------- 263
spectra due to a square lattice array of 7x3 rigid
cylinders of diameter 0.055 m over acoustically-rigid
ground with source coordinates (0, 0.23) m and
receiver coordinates (a) (1.20, 0.11) m, (b) (1.20,
0.23) m and (c) (1.20, 0.35) m.
XIX
8.13 Measured insertion loss spectra due to a 7x3 rigid ------- 265
cylinder array over acoustically rigid ground (solid
line) compared with that of an identical array
measured in free field (broken line) as well as
measured inverted excess attenuation (dash dot line)
The source is at coordinates (0.0.235) m and the
receiver coordinates are (a) (1.20, 0.11) m, and (c)
(1.20, 0.35) m.
8.14 Measured (solid line) and predicted (broken line) IL ------- 267
spectra due to a square lattice array of 7x3 latex
shell cylinders of diameter 0.055 m over acoustically-
rigid ground with source coordinates (0, 0.23) m and
receiver coordinates (a) (0.8, 0.11) m, and (c) (0.8,
0.35) m.
8.15 Measured insertion loss spectra due to a 7x3 latex ------- 269
shell cylinders array over acoustically rigid ground
(solid line) compared with that of an identical array
measured in free field (broken line) as well as
measured inverted excess attenuation (dash dot line)
The source is at coordinates (0.0.235) m and the
receiver coordinates are (a) (1.20, 0.11) m, and (c)
(1.20, 0.35) m.
8.16 Measured insertion loss spectra due to a 7x3 rigid ------- 271
cylinders array over finite impedance (Polyurethane
foam) ground (solid line) compared with that due to
an identical array measured in free field (broken line)
as well as measured inverted excess attenuation
spectra of the ground (dash dot line) taken at same
source-receiver distances. The source is at
coordinates (0.0.235) m and the receiver coordinates
are (a) (1.20, 0.11) m, and (c) (1.20, 0.35) m.
8.17 Measured insertion loss spectra due to a 7x3 latex ------- 273
shell cylinder array over finite impedance
(Polyurethane foam) ground (solid line) compared
with that due to an identical array measured in free
field (broken line) as well as measured inverted
excess attenuation spectra due to the ground alone
(dash dot line) The source is at coordinates (0.0.235)
m and the receiver coordinates are (a) (1.20, 0.11)
m, (b) (1.20, 0.23) m and (c) (1.20, 0.35) m.
8.18 7x3 PVC cylinders in vertical array orientation above ------- 274
a MDF ground (MDF 1).
8.19 Measured insertion loss spectra of vertical array with ------- 275
source and receiver at 0.3m height with and without
the MDF 1 ground plane. Measured inverted excess
attenuation spectra of the ground (dash dot line).
XX
8.20 Measured insertion loss spectra of sonic crystal in ------- 276
vertical orientation with (solid line) and without
(broken line) a ground plane (glass plate). The
source and receiver heights are at (a) 0.1, (b) 0.15
and (c) 0.2 m respectively. Measured excess
attenuation (inverted) spectra of the ground is also
shown (dash dot line).
9.1 Proposed modification for current composite design. ------- 291
B1 Geometry descriptions for two-dimensional addition ------- XXXIII
theorems.
C1 (a) Single LFT and (b) 7x3 square lattice array of LFT ------- XXXV
with lattice constant 0.15 m.
C2 IL spectra for single LFT cylinders of diameter 0.055 ------- XXXVI
m and wall thicknesses (WT) of 0.0001 m and
0.00018 m.
C3 IL spectra for 7x3 square lattice array of LFT with ------- XXXVII
diameter 0.055 m, wall thickness (WT) of 0.0001 m
and lattice constant (L) of 0.15 m.
C4 Single elastic cylinder with periodic stiffening ------- XXXVIII
mounted inside the anechoic chamber.
C5 IL spectra of single elastic cylinder with (solid line) ------- XXXVIII
and without (broken line) periodic stiffening.
D1 Tensile test: stress-strain curve of typical material. ------- XXXIX
D2 Tensile test experiment arrangement for the latex ------- XL
specimen using MTS Universal Testing System.
D3 Results of tests to obtain the Young’s Modulus of ------- XLI
industrial rubber Sample size: length 60 mm, width
25 mm and thickness 0.50 mm.
D4 Young’s Modulus of Polyethylene tubing. Sample ------- XLII
size: length 60 mm, width 25 mm and thickness 0.1
mm.
XXI
List of Tables
XXII
List of Symbols
Symbol Descriptions
f Frequency (Hz)
Wavelength (m)
Angular frequency
Compressibility of air
t Time (s)
pi (3.1416)
i 1
XXIII
p Pressure (Pa)
F (G ) Structure factor
k G ' and eigen Eigenvector
ff Filling fraction
L Lattice constant
XXIV
x, y, z Cartesian coordinates
r ,
j j Polar coordinates centred at the j-th cylinder
2 2
2 Laplacian
x12 x22
Be Unbounded exterior
of n-order)
Theory model
XXV
LMLS Samples number in one period of an m order MLS signal
Acoustically-deduced temperature in ( C )
Diglis weir
Diglis weir
Relative impedance
XXVI
p0,d r Pressure from the direct point source with respect of
vector r (space)
vector r (space)
XXVII
Abstract
XXVIII
Chapter 1
Introduction
expressed by,
c
f , (1.1)
where,
(m).
compressible fluid (i.e. air) will be introduced. Such acoustic waves produce
the aural sensation of sound we encounter in our everyday life. There are also
ultrasonic and infrasonic waves whose frequencies are beyond the audible
Chapter 1 Page 1
limits of human which span from 20 Hz to 20 kHz. The molecules in the
(also known as compressional) waves, where the motion of the air molecules
At any instance, some of the particles move closer together and some are
further apart from their equilibrium positions. This type of wave mode can
travel through solid, liquid and gases. Another type of wave mode to consider
are periodically arranged in space and for simplicity the atoms of the solid
longitudinal wave, for transverse wave the particle move perpendicular to the
Chapter 1 Page 2
(a)
Particle at equilibrium (air/solid)
(b)
Direction of sound propagation Movement of air particle
Figure 1.1: Air molecule patterns during propagation of plane acoustic wave through
infinite space: (a) the arrangement of air molecules at equilibrium positions without
any external force excitation, (b) during plane acoustic waves propagating through the
air medium and (c) during transverse wave propagation through a homogeneous solid
material.
[1], the speed of sound for the longitudinal signal, cL , can be determined from
Chapter 1 Page 3
E 1
cL , (1.2)
1 1 2
where,
E Young’s Modulus,
whereas the speed of sound for the transverse wave can be calculated by,
E
cs
, (1.3)
1 2
wave model, in which sound waves are assumed to have the same direction
direction, the wave equation in terms of the pressure can be expressed as [2],
2 p 1 2 p
, (1.4)
x 2 c 2 t 2
1
where c is the acoustic wave and t is the time. Note that where
c2
and are the equilibrium density and the compressibility of air respectively. It
Chapter 1 Page 4
is worth to note that equation 1.4 is restricted to homogeneous, isotropic fluid
and that small wave amplitude is assumed. For this reason, equation 1.4 is
2 p k02 p 0 , (1.5)
2f
where p is the complex valued function and k is the wave number in
c
air.
The first part of equation 1.6 refers to wave travelling in the x direction and
the second part of the equation refers to wave travelling in the x direction
respectively.
Chapter 1 Page 5
1.2) Complex number notation
wave equation and its solution (see equation 1.5) due to the frequent interest
in simple harmonic waves. If the wave is not simple harmonic, the waveform
quantity. The complex quantity may also be written in polar form as,
p p ei , (1.8)
Chapter 1 Page 6
1.3) The conventional road traffic noise barrier
noise especially during the night time in urban areas and has been identified
vehicular traffic and the large number of people exposed to it, disturbance of
concern. Many studies have shown that exposure to road traffic noise may
and 6]. Road traffic noise control can be achieved through better engineering
road or tyre surfaces and controlling the flow of the vehicles in a particular
area. When these at-source strategies are insufficient to reduce noise, the
level at the receiver. Usually these noise barriers are airtight and sufficiently
dense to shield the noise from the source to the receiver. Most of the sound
energy reaches the receiver only as a result of diffraction around the barrier
edges [7]. In the UK, £5 million are spent annually on highway noise barrier
reduction within the protected areas [8]. The drawbacks of using such barrier
are the aesthetic impact such as restriction of view and natural lighting caused
Chapter 1 Page 7
constraints. The presence of barriers alters the wind profile and turbulence in
their vicinity and this can act to reduce the barrier effect in outdoors.
degradation in the spectral values of insertion loss for frequencies above 800
– 1000 Hz [9].
temporal effects from moving traffic and by vehicle composition and speed.
According to the British Standards (see figure 1.3), the normalised traffic
to 5 KHz, with the main noise energy centred at 1 kHz [10]. It should be noted
that this normalised traffic noise spectrum does not take accounts into the
Decibel (dB) shown in equation 1.9. pdirect and ptotal _ transmitted denote the
pdirect
IL 20 log10 , (1.9)
ptotal _ transmitted
Chapter 1 Page 8
-6
-8
-10
-12
Noise level (dBA)
-14
-16
-18
-20
-22 2 3
10 10
Frequency (Hz)
1.4.1) History
written by Brillouin 1946 [11]. The book covers the mathematical background
wave inside periodic systems [12] in the late 80s generated much research on
Chapter 1 Page 9
velocity of light [13], superlensing effect [14], designing highly efficient
cloaking devices [17], optical computer chips [18] and enhancing surface
structure - in the form of an outdoor modern art sculpture - was made in 1995
was used in their experiment (figure 1.4). In such case of the acoustic audible
known as Sonic Crystal (SC). The cylinders were fixed on a circular platform
which can be rotated around the vertical axis. Sound attenuation was
0.066 the experiment results showed several maxima (sound attenuation) and
gap in the band structure which has a centre frequency at 1.7 kHz could be
attributed to the geometry of the structure as shown in figure 1.4. Ever since
Chapter 1 Page 10
Figure 1.4: Minimalistic sculpture by Eusebio Sempere.
vegetation have been conducted also. It was found that there are two
greater than 2 kHz, noise was attenuated by foliage sound absorption. In both
Chapter 1 Page 11
trees arranged in periodic array in bid to achieve attenuation at low
less width, making them more effective noise screen than typically more
randomly spaced tree belts. More recent work on sonic crystal noise barriers
exploits the use of localised sound absorption properties (i.e. rigid perforated
cylindrical shells filled with recycled rubber crumb material). Both numerical
transmittance spectra of such a sonic crystal [25]. Such design offers the
et al [59]). It is also shown in this work that having three rows of scatterers are
factor with respect to economics and land take. The subject of wind generated
resonance bandgaps below the first Bragg bandgap (i.e. due to the periodicity
of the SC) between 400 and 1600 Hz which is important for traffic noise
Chapter 1 Page 12
reduction. Essentially this design is based on the use of concentric split ring
resonators (SRR). The basis for split ring resonator designs is investigated in
Chapter 7.
1.4.2) Crystallography
in space [27]. When this basic unit is repeated infinitely and joined up in space
by specific lattice points (used to define the structure), the result is the crystal
crystal structure looks the same when viewed from the point r as when
viewed from every point r ' translated by an integral multiple of the translation
r' r T , (1.10)
Chapter 1 Page 13
Thus it can be said that the crystal structure is invariant under translations
and, sometimes, under rotations. The crystal structures are classified into
three categories, that is, one-dimension (1D), two-dimensions (2D), and three-
dimensions (3D) crystal by means of the group theory (see figure 1.7). A
Figure 1.7: Schematic illustrations of crystal structure (a) one-dimensional, (b) two-
dimensional and (c) three-dimensional.
By means of the use of group theory, it has been proved that there are
only one unique 1D periodic system, five 2D and fourteen 3D different lattices.
The majority of this work considers the sonic crystal as 2D arrangement, for
which the five distinct Bravais lattices with their principle lattice vectors and
Chapter 1 Page 14
a1 a1 a1
a2 a2 a2
(d) (e)
Figure 1.8: The five fundamental 2-Dimensionals Bravais lattices (a) square, (b)
oblique, (c) rectangular, (d) centered rectangular and (e) hexagonal.
The parallelepiped defined by the vectors ai forms the well known
primitive cell which defines the smallest possible area of a cell that can be
used as a basis for the crystal structure. The translation of the unit cell
following the vectors ai in space produces the lattice of the periodic system.
As the periodic replication is done in the direct space, such lattice is defined
as the direct space, or direct lattice. Another important feature associated with
the direct lattice is the reciprocal lattice. The reciprocal lattice may be used for
particle momentum is occurring in this reciprocal space (see figure 1.9). The
vectors of the 2D primitive cell in the reciprocal lattice are defined from the
Chapter 1 Page 15
a2
b1 2 , (1.11)
a1 a2
a1
b2 2 , (1.12)
a1 a2
Where both the vectors of the direct ( a1 and a2 ) and the reciprocal ( b1 and b2
b j ai 2 ij , (1.13)
where ij 1 if i j and ij 0 if i j .
Figure 1.9: Illustration of direct lattice points (black dot), reciprocal lattice points (red
dots) and a shade region indicate the reciprocal space.
fraction, ff , and can be determined once the lattice constant and the size of
the scatterers are known. The filling fraction is defined as the ratio between
the volume occupied by the scatterers and the total volume occupied by the
unit cell. If we consider cylindrical scatterers with radius, r0 , the ff for both the
Chapter 1 Page 16
r02
ff square , (1.14)
a2
2r02
ff hexagonal , (1.15)
3a 2
an infinite system may not be attainable but under certain assumptions a finite
conditions if the size of the crystal is much smaller than the wavelength of the
from others. The energy span between the valence band and the conduction
band is called an energy bandgap. The properties of this energy bandgap are
Bloch’s Theorem [28]. The same concept of wave reflection in Bragg’s Law
which refers mainly to light diffraction can be analogously used for sound
waves. It was found that artificial crystals can produce strong peaks over a
range of frequencies that depends on the lattice constant of the array as given
by Bragg’s Law (see equation 1.16). When a wave passes through a crystal,
part of it will be reflected by the first layer while the rest of it will continue to
Chapter 1 Page 17
move to the second layer, some portion of it will again be reflected and the
remaining part will be reflected once more, hence the process continues.
Incident
A A
1
L
θ θ'
B
B
1
Consider two parallel planes of cylinders A-A1 and B-B1 (figure 1.10) having
sound waves W11 and W22 occurs at an angle, and ' to the planes if the
shift is a multiple of 2 .
2L sin n , (1.16)
Chapter 1 Page 18
where, L is the spacing between planes (lattice constant), is the
n is an integer.
the matrix) material. The geometry of the SC takes the form of a single unit
Bragg diffraction (this term and the equivalent term Bragg resonance are used
as it passes through the crystal. Sound waves are therefore forbidden to pass
through in such band gaps which physically can be well described by the
classical wave theory for the multiple scattering of waves within the periodic
wave propagation are also depicted in the figure. This denotes the first
irreducible Brillouin Zone (BZ). For example X refers to the [100] direction
and M refers to the [110] direction, while XM refers to the wave vector
varying from [100] and [110] on the side of the first BZ. The geometrical
Chapter 1 Page 19
parameters, L , which is the lattice constant of the scatterers, and r0 , which is
the radius of each cylindrical scatterer, together define the ff of the square
lattice array in the matrix (see equation 1.14). The Bragg resonance
conditions are essential for opening a wide spectral bandgap utilizing a sonic
crystal. The fundamental Bragg resonance frequencies are found where the
c
f Bragg _ X , (1.17)
2L
c
f Bragg _ M , (1.18)
2 ( 2 L)
where c is the speed of sound in the host medium (i.e. 344 m/s for sound
resonance condition. The width and depth of each Bragg resonance are
scatterer as well as the filling fraction (function of lattice constant). The filling
material around the scatterers. On the other hand, if the filling fraction
Chapter 1 Page 20
becomes too high, constructive wave interference between the scatterers
selection of materials with both the mass densities and modulus to yield the
desired acoustic impedance and velocity mismatch between the matrix and
scatterer.
Z V E , (1.19)
where and K are the density and bulk modulus of the material and the
velocity V , is given by
K
V , (1.20)
Chapter 1 Page 21
r (Zs, Vs)
(ZM, VM)
Г X
Figure 1.11: 2D schematic of square lattice array showing the Brillouin zone (triangle
depicted by points , X and M ). Z M ,Vm and Z s ,Vs denote the acoustic
impedance and velocity for the medium and scatterer respectively.
properties such as density and speed of sound between the scatterer and the
passband filters. The nature of the host materials can be used to differentiate
between sonic and phononic crystals. If the host material (matrix material) is
solid then the term ‘Phononic Crystal’ is used for the artificial crystal. In a
waves) may exist and couple with one another which will add to the
typically made of solid materials and the host matrix is air/fluid to give high
Chapter 1 Page 22
1.5) Aims and thesis organisation
This Thesis offers a study of the interaction of sound with periodic structures.
(ii) To find ways of improving the size and magnitude of the band gaps
(iv) Test and evaluate a full scale SC traffic noise barrier (in-situ).
Among these studies, the most important contributions in this thesis are:
noise barriers.
describe the general context of road traffic noise control. The basic theories
Chapter 1 Page 23
on acoustic wave propagation and concepts of crystallography were reviewed
Chapter 2 reviews the Plane Wave Expansion (PWE) prediction method and
is intended to serve as a starting point for the creation of acoustic band gaps
Theory (MST) on which the subsequent research work has been largely
for 2D finite sonic crystals and the results obtained can be directly compared
sonic crystal using the Finite Element Method (FEM). Such a numerical
sonic crystal design. It is used also for analysis of elastic shells and concentric
Chapter 1 Page 24
Chapter 6 studies numerically and experimentally the feasibility of replacing
the conventional rigid scatterer with thin elastic shell whereby a local
resonance (i.e. breathing mode) effect is exploited to enhance the band gap
design.
Chapter 7 first re-visits of the widely known Split Ring Resonator (SRR) based
(chapter 6) are made. A novel design combines the elastic shell with an outer
performance of the sonic crystal with and without acoustically “hard” or “soft”
Chapter 9 states conclusions for each chapter and offers some ideas for
future work.
Chapter 1 Page 25
Chapter 2
2.1) Introduction
structures. To design sonic crystals and determine how sound will propagate
[12], and the study of wave dispersion bands for electronic waves in solids
[30, 31, 32 and 33] provides the basis of understanding the properties of
waves, the concepts of dispersion relations and Brillouin Zones (BZ) can be
prospect of achieving a complete band gap; this is defined to be the stop band
in which, sound waves are prohibited for all values of the wave-vector (or
Bloch vector). There has been considerable subsequent research into the
analogous phononic or sonic crystal band gaps [34, 35 and 36]. There are
several techniques, but one of the most studied method is the Plane Wave
Expansion (PWE) method [37 and 38] as it can be applied to an infinite arrays
of any scatterer shape. It was used in some of the earliest studies of photonic
crystals and can be readily applied in acoustics. The PWE method uses the
periodicity of the given structure and the Bloch theorem to solve the wave
Chapter 2 Page 26
equation. This allows the computation of the allowable frequencies also
crystal and can be used also to determine the field distributions in the crystal
section presents the theory of acoustic band structures for periodic composite
structures.
one direction (i.e. the z axis). This means that the material parameters of
mass density and velocity of sound in this case depend on the coordinates x
and y . The propagation of the incidental wave is also limited to the xy plane
1 2 p(r )
. p(r ) 0, (2.1)
(r ) (r )cl (r )
2
Chapter 2 Page 27
At every point r the medium is characterized by two material parameters as
follows:
(r ) = mass density
cl (r ) = velocity of sound
Fourier series:
1
(G) ei G r , (2.2a)
(r ) G
1
(G) ei G r , (2.2b)
(r ) cl (r ) G
2
direct lattice vectors. Since (r ) and cl (r ) are both known parameters, the
equations 2.2 is made for all possible reciprocal lattices that correspond to the
Bravais lattice of the system. The solution of the sound pressure field p can
Chapter 2 Page 28
where is called the 2D Bloch wave vector and 2 f is the angular
For a 2D sonic crystal system with square lattice geometry we have to define
some variable in the PWE method. For such a system the cylinder material
where F (G ) is the structure factor, subscripts a and b denotes the host and
scatterer media respectively. Since the system has circular scatterers, the
F (G)
1
Auc
Acyl
d 2 re i Gr 2 fJ 1 Gr0 / Gr0 , (2.5)
such that Auc is defined as the area of the unit cell, Acyl is the area of the
considered cylinder and J 1 is the Bessel function of the first kind of order 1.
After substituting the expanded wave equations (2.2) and the Bloch form
Chapter 2 Page 29
Equation 2.6 is based on a simple modification which excludes the transverse
term in reference [39, equation 9]. For G taking all the possible values,
the eigenvector k G ' and the eigenfrequencies ( ) . The ( ) is introduced
be solved for equation 2.6. The following vector notations are considered for
p(G1 )
p .
(2.7)
p(G N N )
(G1 G1 ) (G1 GN N )
,
(2.8)
(GN N G1 ) (GN N GN N )
(G1 G1 ) (G1 GN N )
.
(2.9)
(GN N G1 ) (GN N GN N )
Chapter 2 Page 30
Having i as the diagonal matrix whose terms are generated by first-order
spatial derivatives:
i Gi 0 0
0 i Gi 0
i , (2.10)
0 i Gi
where i is equal to 1, 2, 3.
p
i 1
i i
2
p . (2.11)
G 2 / a N x N y . (2.12)
Where a is defined as the lattice constant, N x and N y assume all the integer
values. With the substitutions of equations (2.5) and (2.12) into equation (2.6)
and then by solving the system equation given in (2.11) for each Bloch
vectors in the irreducible area of the first BZ for N ,, N , 2 , are obtained
Chapter 2 Page 31
2.2) Plane Wave Expansion Results
used for the initial case. This geometry gives a filling fraction ( ff ) of 50%.
These parameters will describe the contrast parameter and shape function for
the composite structure in the PWE method. Table 2.1 shows the material
compute the band structure in the three principal symmetry directions, letting
k scan the periphery of the irreducible triangle of the first BZ. The dispersion
relation is obtained by plotting frequency against the Bloch wave vector in the
first BZ related to the incident direction of the wave. This is also known as the
reduced wave vector as it is the vector confined to only the first BZ. The PWE
method tends to have a large convergence problem due to the large number
of plane waves needed to calculate the band structure [41]. For reasonable
convergence in our numerical results. Figure 2.1 show the first 10 bands for
MATLAB® (see Appendix A) This indicates any regions where the band gaps
Chapter 2 Page 32
Material Density kgm-3 Speed of sound ms-1
Air 1.2 344
PVC 1380 2400
Table 2.1: Material parameters for sonic crystal system comprising of PVC cylinders
embedded in air.
Assuming linear dispersion, the lowest Bragg band gap central frequency is
approximately given by
cl
Bragg . (2.13)
2L
homogeneous material. The PWE method solved for the first ten bands with
material parameters equal to that of air, a b = 1.2 kg/m3 and cla clb = 344
6000
M
5000
Г X
4000
Frequency (kHz)
3000
2000
1000
0
M X M
Reduced Wave Vector
Figure 2.1: Plane wave expansion computed band structure for a homogeneous
medium with equivalent material properties of air. Inset: Brillouin zone. refers to
[1,0] direction, M refers to [1,1] direction and refers to the wave vector
varying from [1,0] to [1,1] on the extreme side of the Brillouin zone.
Chapter 2 Page 33
Clearly in figure 2.1, it is evident that such a homogeneous system has
a dispersion relation with an almost linear trend with phase velocity and group
velocity being equal. Such structure does not exhibit any band gap properties.
9000
8000 M
7000
Г X
Frequency (kHz)
6000
5000
4000
3000
2000
1000
0
M X M 0 10 20 30 40 50
Insertion Loss (dB)
Reduced Wave Vector
Figure 2.2: The left hand panel shows the plane wave expansion computed band
structure for a SC consisting of cylindrical PVC scatterers embedded in air in a square
array with lattice constant L of 0.069 m and filling fraction of 50%. The right hand
panel shows the laboratory measured IL spectrum for 7x3 array for the [1, 0]
direction. Inset: schematic of a Brillouin zone. refers to [1,0] direction, refers
to [1,1] direction and refers to the wave vector varying from [1,0] to [1,1]. The
(red) dotted box indicates a full bandgap and the (blue) shaded regions indicate partial
bandgaps.
The left hand panel in figure 2.2 shows the band structure due to a SC
shading indicates the location of the first band gap ( direction) which
spans from 1.6 kHz to 3.2 kHz with centre frequency at 2.4 kHz. The lowest
Bragg band gap central frequency, calculated using equation 2.13 is 2.5 kHz.
This calculated central frequency falls well within the PWE method
computation for the mid frequency in the first band observed in the band
Chapter 2 Page 34
bandgaps) in the direction. It is also worth noting that a complete band
gap can be achieved for such SC configuration marked by the red dotted box.
A study of how the lattice parameter influences the location of the band
gap has been performed. The diameter of the scatterer has kept constant at
10000
9000
8000
7000
Frequency (kHz)
6000
5000
4000
3000
2000
1000
(a)
0
M X M 0 10 20 30 40
Insertion Loss (dB)
Reduced Wave Vector
10000
9000
8000
7000
Frequency (kHz)
6000
5000
4000
3000
2000
1000
(b)
0
M X M 0 10 20
Insertion Loss (dB)
Reduced Wave Vector
Figure 2.3: Plane wave expansion predictions of the effect of varying the lattice
constant in a SC comprising of 0.055 m diameter PVC cylinders embedded in air.
Lattice constant L of (a) 0.08 m and (b) 0.135 m. Left hand panels show the band
structure of the corresponding SC and right hand panels show the measured IL for the
different lattice parameters.
Chapter 2 Page 35
For a lattice constant of 0.08 m, the plane wave expansion method is
solved for the first 10 bands whereas the first 40 bands would be required to
extend the range of frequencies to 9 kHz. Focusing on the lowest band gap at
direction for SC with lattice constant of 0.08 m (see figure 2.3a – area
marked blue), the band gap lies between 1.4 kHz to 2.6 kHz, centred at 2
kHz. These band gaps are shifted to lower frequency as compared with the
shift of the first Bragg band gap to an even lower frequency when the lattice
constant is further increased to 0.135 m (see figure 2.3b). The first band gap (
direction) predicted by PWE method spans from 1 kHz to 1.4 kHz, and is
centred at 1.2 kHz (area marked blue). The lowest Bragg band gap central
with the PWE method for lattice constants of 0.08 and 0.135 m. The
corresponding measured IL spectral for both lattice constant are also found to
agree well on the first Bragg band gap. From this investigation it is clear that
other hand, for the same number of array elements this increasing the lattice
Chapter 2 Page 36
2.2.3) Predicted influence of filling fraction
have been studied. It should be noted that in the previous section 2.2.2, the
current study the lattice constant is kept constant but the scatterer size is
Chapter 2 Page 37
10000 10000
(a) 9000
(b) 9000
8000 8000
7000 7000
Frequency (kHz)
Frequency (kHz)
6000 6000
5000 5000
4000 4000
3000 3000
2000 2000
1000 1000
0 0
M X M M X M
10,000 10000
(c) 9,000
(d) 9000
8,000 8000
7,000 7000
Frequency (kHz)
6000
Frequency (kHz)
6,000
5,000 5000
4,000 4000
3,000 3000
2000
2,000
1000
1,000
0
0 M X M
M X M
Reduced Wave Vector
Reduced Wave Vector
Figure 2.4: Plane Wave Expansion predictions of the band structure with Lattice constant fixed at 0.135 m and varying scatterer diameters (a)
0.055, (b) 0.09, (c) 0.11 and (d) 0.13 m.
Chapter 2 Page 38
According to the PWE predictions, the wave propagation in different
fractions. An indication of the complete band gap can be observed when the
filling fraction exceeds certain value. For systems with low filling fractions of
13% and 35%, narrow pseudogaps are observed for the direction in the
band structure (see figure 2.4a and 2.4b). Although there have the same
number of band gaps predicted for these two filling fractions at similar
frequencies, careful examination reveals that they are not identical. The
lowest band gap extends from 1030 to 1390 Hz for the case of 13% filling
fraction whereas the lowest band gap for 35% filling fraction extends from 850
to 1530 Hz. Also the mid frequency for the first band in both filling fractions
occurs at about 1.2 kHz. It is evident that as the filling fraction increases
further to 52% and 73%, the band gaps are widened and additional band gaps
are formed. Complete acoustic band gaps are allowed to form for such filling
overlapped. Two complete band gaps are observed for the range of
Chapter 2 Page 39
2.2.4) Predicted influence of material parameters
This section looks at how changes to the density and speed of sound in
the scatterer material and the embedding medium will influence the band gap
speed of sound are listed in table 2.3 [42]. The lattice constant and diameter
in air and the results are shown in figure 2.4. Secondly, water is used to
Chapter 2 Page 40
10000 10000 10000
Frequency (kHz)
Frequency (kHz)
Frequency (kHz)
6000 6000 6000
1000
(a) 1000
(b) 1000
(c)
0 0 0
M X M M X M M X M
10000 10000
10000
9000 9000
9000
8000 8000
8000
7000 7000
7000
Frequency (kHz)
Frequency (kHz)
Frequency (kHz)
1000
(d) 1000
(e) 1000 (f)
0 0 0
M X M M X M M X M
Figure 2.5: Plane Wave Expansion predictions of the influence of material parameters on the band structure. The lattice constant is fixed at 0.069
m and scatterers are made from different materials (a) PVC, (b) silicone rubber, (c) steel (d) PMMA, (e) wood and (f) PVC scatterers embedded
in water
Chapter 2 Page 41
Figures 2.5(a), (b), (c), (d) and (e) show the predicted band structures
for scatterers made from PVC, silicon rubber, steel, PMMA and wood
respectively. The width of the first Bragg band gaps are identical for all of
these scatterers and the first Bragg band gaps lie between 1490 to 3170 Hz.
This demonstrates that so long as the scatterers are acoustically hard they
have no effect on the band structure of the SC. On the other hand, when PVC
2.5(f)). In such system, fewer band gaps are predicted for the direction
since the impedance contrast between PVC and water is smaller than that
Chapter 2 Page 42
Chapter 3
3.1) Introduction
remain undisturbed while others are scattered and spread out from the
scatterer in all directions distorting and interfering with the incident wave. If
wave spreads uniformly in all directions from the scatterer and the rest is
destructively with the unchanged incident wave. On the other hand, if the
scatterer is small compared to the wavelength (which is often the case with
sound wave) then, the entire scattered wave is propagated in all directions
and there is no interference behind the scatterer [2]. For a single obstacle,
calculation of the total field includes the sum of the scattered field and the
incident field. This assumption is valid only if the scatterer placed far enough
from other scatterers and the incident wavelength is small compared to the
distance between the scatterers. Some authors have suggested that the
four times the radius of the scatterer for independent scattering to be valid
scatterers are placed within close proximity to each other: a situation often
Chapter 3 Page 43
waves dates back to 1913 [45] when a method for predicting scattering of
acoustic waves by two dimensional finite arrays was described. Ever since,
known to have taken place are electromagnetic waves [46 and 47],
condensed matter physics [48], water waves [49], optics [50 and 51] and
acoustics [52].
fields due to two or more scatterers [53]. As explained previously for single
total scattered field is just the sum of the incident field and the field scattered
scatterers within close proximity) takes account of the field scattered from one
scatterer which will induce further scattered fields from all the other scatterers,
which will again induce further scattered fields from all the other obstacles and
is valid for both periodically and randomly placed scatterers. In this section,
Chapter 3 Page 44
problems e.g. Mal and Bose [55] in 1974 and Linton and Evans in 1990 [49].
Martin’s book [53] explains each technique developed during these periods of
time in great detail. In this section, periodically arranged rigid cylinders will be
considered.
N
pinc p scj , (3.1)
j 1
where pinc is the incident wave and p scj is the wave field scattered by the j-th
scatterer.
Then the effective field over the k-th scatterer which is the radiation incident
on the k-th scatterer in the presence of all other scatterers (N-1) can be
defined as,
N
pk pinc pscj . (3.2)
j 1
j k
pscj j pinc
j
, (3.3)
where j is an operator relating the field incident on the j-th scatterer, pinc
j
, to
the field scattered by the j-th scatterer, pscj . Thus substituting equation 3.3
Chapter 3 Page 45
N
p k pinc j pinc
j
, (3.4)
j 1
j k
or equivalently,
N
j
p k pinc p sc ,
k
sc (3.5)
j 1
j k
by solving equation 3.4 for pk or equation 3.5 for psck where k = 1, 2…, N, the
N
p pinc j p j , (3.6)
j 1
This over simplified equation hides the complexity through the operator
j which has not been clearly defined at this stage and also where equation
3.5 and 3.6 is required to hold in space. The next section will show how the
multipoles method.
Chapter 3 Page 46
3.2) Multipole method for circular scatterers in 2-dimensional
system
matrix [56], iterative technique [57] and the multipole methods [45]. The
and spherical scatterers. This method was developed by Zaviska, 1913 for
the multipole method using both plane and cylindrical incident waves
coordinate system, some distance away from the array. We will introduce N+1
coordinate systems in the x-y plane normal to the cylinder axes r , centred
at the origin and r j , j centred at the j-th cylinder. The vertical circular
cylinders are fixed along the z-axis and it is assumed that waves propagate in
the plane perpendicular to the main axis of the cylinders. Since the boundary
Chapter 3 Page 47
conditions and the geometry do not change with z, the problem can be
reduced to a two-dimensional one where the scatterers are taken as its cross
constant) between the j-th and q-th scatterers. The distance of the centre of
y
j-th scatterer at (xj, yj)
Receiver (Px, Py)
rj
j j
rq
Rj
Rjq q
r
Arrow showing
incident wave q-th scatterer at (x , y )
q q
direction
j
β
x
Source (0,0)
Figure 3.1: Plan view of two identical cylinders and corresponding cartesian and polar
coordinates [49, figure 1].
1 2 p
p 2 2 ,
2
(3.7)
c t
2 2
where 2 is the Laplacian, c is the constant speed of sound t is
x12 x 22
Chapter 3 Page 48
Assuming that an external wave pinc with temporal dependence e it
can write,
pinc Re p e it . (3.8)
equation as shown,
2 p k02 p 0 in Be , (3.9)
2f
Where k 0 , is the wave number in air and Be is the unbounded exterior
c0
to the cylinders region occupied by the air medium. Thus, solution can be
Each individual cylinder scatters the waves which are incident upon it
and the total field around the j-th cylinder is a superposition of the external
field and the radiation scattered by the rest of cylinders. Solutions of the
[53] which can be expressed in term for the scattering of waves by cylindrical
ensures that the Sommerfeld radiation condition (in two dimensions) at infinity
Chapter 3 Page 49
is satisfied. Such a condition implies that there are no incoming scattered
r cos , r sin .
The general solution for solving the total wave field at any point exterior to the
cylinders (position that is within the host medium which is air in this case) will
px, y pinc Anj Z nj H n(1) k 0 r j e
N
in j
, (3.12)
j 1 n
where pinc is the incidence wave over a two dimensions system of scatterers
in the form of either a plane wave from infinity or a cylindrical wave from a
point. It is also worth noting that pinc can also be modelled with a spherical
wave which will extends the problem to three dimensions but will necessarily
of this equation will be discussed further in the next section where incident
wave pinc , coefficient Anj , and impedance factor Z nj will be defined. We
Chapter 3 Page 50
will assume in our model that all motion is time harmonic with angular
propagating from the negative values of x axis. k and r are defined as the
wave vector (with k 2f as the angular frequency) and location respectively.
wave with angle propagating to the j-th cylinder as shown in figure 3.1
Chapter 3 Page 51
pinc eikr cos . (3.16)
The plane wave impinging the j-th scatterer is then scattered to other
positions (e.g. the receiver position) thus, the scalar product rj k is now
replaced with R j rj k . Following a similar approach in equation 3.14 gives,
where,
Therefore,
ikr j cos j
pinc I j e , (3.19)
ik x j cos y j sin
e can be written as series of Bessel function [58], equation 3.19
n
i J kR e
in j
pinc n
n j , (3.20)
n
2f
where k , and J n is the Bessel function of n-th order of the first kind.
c
Chapter 3 Page 52
n
J kr e
in 2 j
pinc I j n j , (3.21)
n
emanating from the j-th cylinder can be expressed with respect to the origin of
coordinates of the j-th cylinder using n-th order Hankel functions of the first
kind,
n
A H kr e
in j
pscj n
j
n
1
j , (3.22)
n
where
rj x x y y
j
2
j
2
, (3.23)
and
y yj
j sin 1 , (3.24)
r
j
3.21) that includes a phase factor associated with the corresponding j-th
scatterer.
Substituting equation 3.21 and 3.22 into equation 3.2, we are able to
solve for the total wave over the j-th cylinder. However, in order to perform the
computation, all the terms must be expressed in terms of the same origin
Chapter 3 Page 53
coordinates rq , q . This can be done using Graf’s addition theorem (See
Appendix B). When using the addition theorem for Bessel functions it is
coordinates r , ,
q q we have to ensure that rq R jq for all j q . This
geometry restriction implies that the expression obtained is only valid if the
point rq , q is closer to the centre of cylinder q rather than to the centre of
any other cylinders or the source and thus we can write the expression
obtained for rq , q as,
I J k r e
n
Ank Z nk H n1 k0 rq e
in 2 q p in q
prq , q q n 0 q
n
N n m e i nm j , (3.25)
J k r H k R
im q
Anj Z nj 1 e
m 0 q nm 0 jq
j 1 n m
jq
for some set of unknown complex numbers Anj . Both Anj and Z nj are related by
particle velocity between the scatterer and the surrounding air medium. Thus,
for the case of the q-th scatterer, the boundary condition can be expressed
by,
Chapter 3 Page 54
Leading to,
1 ext 1 int
rq aq rq aq , (3.27)
rq q rq
where is the density of the surrounding medium (air) and q is the density
condition which specifies the values that the derivative of a solution is to take
ext
pinc
q
pscj 0, (3.28)
rq aq rq aq
rqp rq
where pinc the total incident pressure over the q-th scatterer (refer to equation
3.21 and pscj is the scattered wave by the j-th scatterer (refer to equation
3.22).
This boundary condition leads to the definition of impedance for the rigid
scatterers as follow,
J n' k 0 aq
Z j
H n1' k 0 aq
n , (3.29)
Applying this boundary conditions for pressure and velocity on the outer
surface of q-th scatterer and after using the orthogonality of the functions eimq
Chapter 3 Page 55
N M
H n1M k 0 R jq I q e im
i n m 2 j
Amq A n
j
Z nj e 2
. (3.30)
j 1 M
j q
The obtained algebraic system can then be truncated to the finite number
it was found that taking M = 6 gave results accurate to four significant figures,
The previous scattering model was for incident plane waves. The
source may not be placed sufficiently far enough from the scatterers for such
acoustic source as a line source located at origin r 0 . The wave in free field
0
2 k 2 p r 4 2 r in B ,
e (3.31)
2f
where k 0 , 2 is the 2-dimension delta-function and Be is the
c
where H 01 is the zero-th order of Hankel function and of the first kind.
Chapter 3 Page 56
Following similar procedure in equation 3.21 where we considered the
radiating wave emanating from the j-th cylinder can be expressed with respect
to the origin of coordinates of the j-th cylinder using n-th order Hankel
n
A H kr e
in j
pscj n
j
n
1
j . (3.33)
n
Similarly by applying Graf’s addition theorems to solve for the total wave
I
n
J n k 0 rq e Ank Z nk H n1 k 0 rq e
in 2 q p in q
prq , q q (3.34)
n
N n m ei n m 2 j .
J m k 0 rq H n1m k 0 R jq
im q
Anj Z nj
e
j 1 n m
q
J n' k 0 aq
Z nj
H n1' k 0 aq
, (3.35)
Using the same approach as for equation 3.30, i.e. applying boundary
conditions for pressure and velocity on the outer surface on the q-th scatterer
coefficients Amq with cylindrical incident waves [59]. The obtained algebraic
Chapter 3 Page 57
system can then be truncated to the finite number N 2M 1 of equations.
N n
H n1M k 0 R jq H m1 k 0 Rq e
i n m 2 j
Amq A j im p
n Z nj e . (3.36)
j 1 n
j q
in a 7x3 square lattice array. The radius of the scatterers is 0.0275 m. Three
lattice constants (L) of 0.069, 0.08 and 0.135 m were used and thus give filling
fractions of 50, 37 and 13% respectively. The source is located at 1.5 m away
from the front face of the array and the receiver is placed 0.05 m from the
nearest face of the array but on the opposite side of the source. Such
normal to the array orientation (see figure 3.2). The results are computed in
terms of Insertion Loss (IL) which has been defined in equation 1.9. However
without loss of generality and for clarity in term of the notation used in this
pinc
IL 20 log 10 , (3.37)
pr
incident wave (equation 3.31) and the total field pr defined in equation 3.12.
Chapter 3 Page 58
The calculated results using MST are validated with the laboratory
numerical solution of the MST is also done since there is an infinite linear
system of equations.
1.5 m 0.05 m
Source Receiver
0.0275 m
Chapter 3 Page 59
50
(a)
40
7x3 array, scatterer radius 0.0275 m and lattice constant 0.069 m
Measurement
MST prediction, cylindrical wave
30
MST prediction, plane wave
Insertion Loss (dB)
20
10
-10 2 3
10 10
Frequency (Hz)
50
(b)
40
7x3 array, scatterer radius 0.0275 m and lattice constant 0.08 m
Measurement
MST prediction, cylindrical wave
30
Insertion Loss (dB)
20
10
-10 2 3
10 10
Frequency (Hz)
50
(c)
40
7x3 array, scatterer radius 0.0275 m and lattice constant 0.135 m
Measurement
20
10
-10 2 3
10 10
Frequency (Hz)
Figure 3.3: MST predictions and measured Insertion Loss spectra for square lattice
arrays of rigid circular scatterers of 0.055 m diameter with lattice constants of (a)
0.069, (b) 0.08 and (c) 0.135 m respectively. Both plane and cylindrical waves are
compared for all three cases.
Chapter 3 Page 60
According to figure 3.3, the theoretical (MST) predictions for both plane
and cylindrical incident waves show relatively good agreement with the
laboratory measurement data for 7x3 square lattice array of rigid circular
scatterers. The shifting of the band gaps due to the effect of changing the
periodicity by varying the lattice constant for the structure can be observed.
The first Bragg diffraction frequency can be predicted by using the formulation
finite sonic crystal array. The peaks near 2.5, 2 and 1.3 kHz for 50, 37 and 13
% filling fractions respectively correspond to the first Bragg frequencies for the
relatively minor peaks occurring before the first Bragg frequencies can be
Chapter 3 Page 61
50
(a)
40
7x3 array, scatterer radius 0.0275 m and
lattice constant 0.069 m, plane wave
Measurement
30 MST prediction, M = 1
Insertion Loss (dB)
MST prediction, M = 2
MST prediction, M = 3
MST prediction, M = 4
20 MST prediction, M = 5
10
-10 2 3
10 10
Frequency (Hz)
50
(b)
40
7x3 array, scatterer radius 0.0275 m and
lattice constant 0.069 m, cylindrical wave
Measurement
30 MST prediction, M = 1
Insertion Loss (dB)
MST prediction, M = 2
MST prediction, M = 3
MST prediction, M = 4
20
MST prediction, M = 5
10
-10 2 3
10 10
Frequency (Hz)
Figure 3.4: MST predictions and measured Insertion Loss spectra for square lattice
arrays of rigid circular scatterers of 0.055 m diameter with lattice constants of 0.069.
Different truncation number of M 1, 2, ... , 5 is used for the MST predictions for (a)
plane wave and (b) cylindrical wave.
M, for both plane and cylindrical incident waves. This test is to determine the
0.069 m for both types of incident wave and increment M from 1 to 5. It can be
observed that for our geometry, the accuracy of the numerical results
converge as long as M 4 for the case of plane waves and M 3 in the case
Chapter 3 Page 62
of cylindrical waves. Using M 4 and M 3 for plane and cylindrical incident
standard office PC (Intel core 2 Dual CPU 2.4 GHz, 2GB RAM with 32 bit
Operating System).
Chapter 3 Page 63
Chapter 4
Numerical method using Finite Element Method
(FEM)
4.1) Introduction
The Finite Element Method (FEM), also known as Finite Element Analysis
and total fields in a wide range of physical and engineering problems. The
The history of FEM dates back to 1943 when Courant determined the
torsional rigidity of a hollow shaft by dividing the cross section into triangles
and interpolating values of the stress function at each intersecting point [60].
Subsequently, the first publication using FEM for determining the frequencies
and mode shapes of acoustic cavities was published by Gladwell [61]. Beside
of having the capability to solve for a structure with an irregular shape, the
Chapter 4 Page 63
decisions. In the early 70’s, the use of FEM was limited to the expensive
technique has become much more widely available. Driven by research and
and 63], optics [64 and 65], biomechanics [66], thermodynamics [67],
about 3800 papers on FEM were being published annually and that the
generally yields five standard problems. Firstly, the radiation problem where
we have a source (i.e. speaker) radiating sound into the surrounding space.
scattered wave. Thirdly, the transmission problem where the incident sound
Chapter 4 Page 64
interactions (also known as structural acoustics) where the scattering
(LISA) [72]. In this chapter, COMSOL® Multiphysics (version 3.5a) has been
crystals. This software is relatively easy to use due to its user friendly
graphical user interface simulation environment for all of the steps of the
figure 4.1). The particular arrangement of elements is called a mesh. The set
Chapter 4 Page 65
of equations for each element can be assembled to compute the behaviour of
the system over the entire domain. These elements are connected at specific
points, called nodes, and the assembly process requires that the solution be
Domain
Element
Boundary
node
Figure 4.1: Definition of domain discretized using triangular shape finite elements.
fluid medium previously shown as chapter 3, equation 3.9 is repeated here for
convenience,
2 p k02 p 0 , (4.1)
2 2 2f
where 2 is the Laplacian, k 0 is the wave number in air,
x1 x2
2 2
c0
Chapter 4 Page 66
The expressions for the pressure emitted by:
pin _ c
i
J 0 k 0 rc iY0 k 0 rc , (4.3)
4
location to the travelling wave position. There have been more advanced
studies allowing for oblique incidence of the waves [73] but, in our model, we
consider that the direction of waves is normal to the symmetry axis of the
scatterer.
the Neumann boundary condition) to model the rigid surface of the scatterers
1
n. p 0 , (4.4)
0
surrounding boundaries of the modelling domain (see figure 4.2). For the
Chapter 4 Page 67
outward travelling wave, this boundary condition ensures minimal or no
reflections from the model boundary. The radiation boundary conditions can
be expressed as [74],
1 p
n p ik 0 r ik 0 r i n 0 e i r ,
p
(4.5)
0 0 0
plane wave: r 0 , or
cylindrical wave: r
1
,
2r
in which r is define as the shortest distance from the point on the boundary to
the source.
The term on the right hand side represents an incoming pressure wave with
normal vector, n .
least 6 to obtain the desired accuracy. The accuracy of FEM will be discussed
in section 4.4. Equation 4.1 can now be solved using a parametric solver to
obtain the pressure field for the rectangular domain modelled as air. This is
step of 100 Hz to produce a frequency spectrum for both direct field and total
Chapter 4 Page 68
transmitted field (without and with sonic crystal respectively). The Insertion
Loss (IL) for the sonic crystal noise barrier can then be obtained using chapter
pdirect
IL 20 log10 , (4.6)
ptotal_ transmitted
The hardware configuration used to run the simulations was a 64-bit windows
platform with 3.07 GHz Intel core i7 processor and 6 GB of RAM. This gives a
maximum solution time of 12 minutes for the case with the highest mesh
Receiver point
Figure 4.2: Definition of geometry for a rectangular domain modelled as air and an
incoming cylindrical wave source.
Chapter 4 Page 69
4.3) FEM computed results
5 kHz
(a)
5 kHz
(b)
Figure 4.3: COMSOL® computed pressure map for a rectangular domain modelled as
air and an incoming plane (a) or cylindrical wave (b) source.
In figure 4.3, the FEM-predicted pressure maps are shown for incident
300 Hz 1.2 kHz
5 kHz (a) plane and (b) cylindrical waves for direct field computed in a
from the left boundary to the right and the radiation boundary conditions are
applied to both simulations. Clearly, in the case for plane wave, the
radiation boundary conditions applied may not be ideal for such geometry. On
the other hand, this radiation boundary conditions work very well for cylindrical
Chapter 4 Page 70
(a) (b)
3 kHz 4 kHz
(c) (d)
4.5 kHz 5 kHz
(e) (f)
Figure 4.4: COMSOL computed pressure maps for 7x3 array of sonic crystal
(acoustically hard scatterer) modelled in rectangular air domain. Cylindrical waves is
performed and pressure maps at 200 Hz, 1.2, 3, 4, 4.5 and 5 kHz are shown for figure
(a), (b), (c), (d) and (e) respectively.
Figure 4.4 shows the FEM-predicted pressure maps for a square lattice
array of 7x3 rigid circular scatterers with lattice constant of 0.135 m modelled
with similar boundary conditions to those used for the direct field
The cylindrical wave source is located at 1.5 m away from the front face of the
array and the receiver is positioned 0.05 m from the nearest face of the array
but on the opposite side of the source. The number of triangular elements for
this model is 157392. The pressure maps are shown for various frequencies
ranging from 300 Hz to 5 kHz. Clearly, these results support the introductory
Chapter 4 Page 71
pass through relatively unaffected when the wavelengths are large compared
to the scatterers (see figure 4.4(a)) but as the incident wavelengths decrease,
25
20
7x3 array, rigid scatterers, L = 135 mm
Measurement
15 FEM prediction, cylindrical wave 9.361 dB
FEM prediction, plane wave
Insertion Loss (dB)
10
6.8 dB
-5 6.358 dB
-10 2 3
10 10
Frequency (Hz)
Figure 4.5: FEM predictions and laboratory measured Insertion Loss spectra for 7x3
square lattice arrays of rigid scatterers of radius 0.0275 m with lattice constant 0.135
m. FEM predictions are performed with plane and cylindrical wave sources for a
rectangular domain modelled as air.
Figure 4.5 shows the Insertion loss (IL) spectra comparing laboratory
measurements and FEM predictions using plane and cylindrical wave source
with the boundary conditions similar to figure 4.4. The details of laboratory
4.6 and the calculated Bragg frequency (refer to PWE chapter, equation 2.13)
for this particular set up is 1.3 kHz. Results show that FEM predictions agree
fairly well with the laboratory measurements up to 5 kHz using plane waves
and up to 6 kHz in the case of cylindrical waves. The IL values at 1.3 kHz for
all the 3 spectra are shown in the figure. The discrepency between the
prediction and measurement at the Bragg frequcny is about 3 dB. The Bragg
Chapter 4 Page 72
agreement could be improved using higher mesh densities which inevitably
the case for plane wave. Peaks occurring in the measurements before the first
Bragg frequencies could be due to the finite height of the array (finite length of
the cylinders).
problems with complex geometries. But also the method lends itself easily to
using FEM, the mesh must offer sufficient resolution. Generally for FEM, this
requires mesh grids with around ten nodal points per wavelength [75]. The
edges are present in the geometry. This section will discuss how frequency
domain. The FEM computation is solved for the different mesh densities listed
in table 4.1 using a parametric solver (frequency sweep over 100 Hz to 8 kHz
Chapter 4 Page 73
with frequency step of 100 Hz) to obtain the pressure field. The computation
time taken to completely solve the 80 iterations for each mesh density is also
shown in the table. Notably, the computation time increases with each
First, the different set of mesh densities will be solved and the pressure maps
stimulating graphics for the eyes but it can be tedious to analyse for all
spectrum for a fixed point in all the sets. The Sound Pressure Level, SPL
p
SPL 20 log10 , (4.7)
p0
5
pressure 2 10 Pa . Note that the SPL is measured in root mean square
(rms) of the pressure wave and is representing the energy level of the sound
Chapter 4 Page 74
Mesh element size 200 Hz 500 Hz 1 kHz 3 kHz 5 kHz 8 kHz
0.5 m
0.2 m
0.1 m
0.04 m
0.007 m
Figure 4.6: FEM predicted pressure plots for cylindrical waves with different mesh element sizes.
Chapter 4 Page 75
Inaccuracy of the results will arise when there are an insufficient
(short wavelength) solutions when using a mesh with large element size
which leads to the non-convergence of the solution and so the pressure map
illustrate that solutions are unresolved at high frequencies ranges for low
mesh density and that the resolution of the computed pressure map at higher
80
60
40
Sound Pressure Level (dB)
-20
-40
-60
300 1000 8000
Frequency (Hz)
Figure 4.7: Sound pressure level spectra at single point position for FEM
computations of the field due to a cylindrical wave computed using different mesh
element sizes.
Figure 4.7 shows the SPL in dB for the range of frequencies. Evidently
insufficient mesh density (i.e. Roll off from the nominal SPL value of near 60
dB. The frequency spectrum for the highest mesh density demonstrates that
Chapter 4 Page 76
Bearing this convergence study in mind, the maximum element size of
This set the convergence criterion to about 3%. This allows good
maximum element size must be at least 6 times smaller than the shortest
wavelength.
The relative ease of modelling different shapes in FEM allows the study
scatterers and different lattice configurations. The first part of this section
looks at predictions of the insertion loss spectra due to 7x3 arrays with a filling
triangular, square and elliptical shapes, and compares them with predictions
for cylindrical scatterers in an array of similar size. These studies are also
performed with the similar source and receiver distances mentioned in section
scatterer considered previously. The array size and lattice constant are
Chapter 4 Page 77
identical with the reference sonic crystal with cylindrical scatterers. Three
different orientations are studied and the pressure maps at 1.3 and 2.3 kHz
are compared (see figure 4.8). Comparisons are made also with an equivalent
(a) (b)
Arrangement 2 (1.3 kHz) Arrangement 2 (2.3 kHz)
(c) (d)
Arrangement 3 (1.3 kHz) Arrangement 3 (2.3 kHz)
(e) (f)
®
Figure 4.8: FEM (COMSOL ) predicted pressure maps comparing three different
arrangements (top to bottom) and 2 frequencies pressure maps for each arrangement
(left 1.3 kHz and right 2.3 kHz) of 7x3 triangular rigid scatterer arrays (square lattice)
with lattice constant of 0.135 m.
Figure 4.8 shows the FEM computed pressure maps for the 3 different
For all the 3 arrangements, similar scattering patterns are observed at 1.3 kHz
which is the first (Bragg) band gap frequency due to the periodicity. At the
higher frequency of 2.3 kHz, the scattering patterns vary from one
Chapter 4 Page 78
arrangement to the other. Notably at this frequency, higher sound energy is
35
30
Triangle scatterers
Arrangement 1
25 Arrangement 2
Arrangement 3
20 Circular scatterers
Insertion Loss (dB)
15
10
-5
-10
500 1000 8000
Frequency (Hz)
Figure 4.9: FEM predictions of Insertion Loss spectra for 7x3 square lattice arrays of
triangular rigid scatterers with lattice constant 0.135 m. A cylindrical wave source is
assumed and a rectangular domain is modelled as air.
Figure 4.9 shows the FEM predictions for all three triangular scatterer
For all triangular scatterer arrangements, the Bragg band gaps near 1.3 kHz
are observed to result in nearly 3 dB higher IL than that for the cylindrical
an additional band gap is observed around 1.8 kHz and the third band gap
(2.3 kHz) is also higher in amplitude than those for the other designs. In
Chapter 4 Page 79
4.5.2) Square scatterers
studied in a similar fashion to that used for the triangular scatterer arrays. The
length of the square sides is modelled as 0.049 m to obtain the same filling
of square scatterer arrays. Cleary the scattering patterns at 4.3 kHz are
Chapter 4 Page 80
35
30
Square scatterers
25 Arrangement 1
Arrangement 2
Circular scatterers
20
Insertion Loss (dB)
15
10
-5
-10
500 1000 8000
Frequency (Hz)
Figure 4.11: FEM predictions of Insertion Loss spectra for two arrangements of 7x3
square lattice arrays of square rigid scatterers and a reference array of cylindrical rigid
scatterers with lattice constant 0.135 m. A cylindrical wave source is assumed and a
rectangular domain is modelled as air.
Figure 4.11 shows the FEM predicted IL spectra for the 2 square
Bragg band gaps (1.3 kHz) are predicted to be of similar size and amplitude
for all the designs. The results show no significant advantage in using square
scatterers.
reported that certain configurations can vary the direction of propagation due
to the geometry anisotropy of the elliptical scatterer [76]. In our model, the
0.0218 m. This geometry provides the same filling fraction as that assumed
Chapter 4 Page 81
for different scatterer shapes previously. The elliptical scatterer arrays are
(a)
(b)
(c)
Figure 4.12: FEM (COMSOL®) predicted pressure maps for three 7x3 arrangements
of elliptical rigid scatterers (square lattice) with lattice constant of 0.135 m at 1.5 kHz
Chapter 4 Page 82
Figure 4.12 shows FEM predicted pressure maps for 3 arrangements
of elliptical scatterer arrays. The scattering patterns for all arrangements are
25
20
Elliptical scatterers
Arrangement 1
15 Arrangement 2
Arrangement 3
Insertion Loss (dB)
Circular scatterers
10
-5
-10
500 1000 8000
Frequency (Hz)
Figure 4.13: FEM predicted Insertion Loss spectra for 7x3 square lattice arrays of
elliptical rigid scatterers with lattice constant 0.135 m compared with that predicted
for an equivalent cylindrical scatterer array. For the FEM predictions a cylindrical
wave source is assumed and the rectangular domain is modelled as air.
receiver position are compared with that measured for the equivalent
cylindrical scatterer array. Arrangement 1 which the elliptical scatterer has its
predicted to produce the highest amplitude Bragg band gap at 1.3 kHz. This
suggests that much of the incoming wave front is blocked by the larger
surface of the scatterer. Also the Bragg band gap is predicted to contain two
peaks instead of one for this particular arrangement. The lowest amplitude
Bragg band gap is predicted for arrangement 2 where the semi-minor axis is
Chapter 4 Page 83
incident waves. Arrangement 3 with scatterers having 45° orientation is
frequencies.
35
30
Triangle scatterers, arrangement 3
Elliptical scatterers, arrangement 1
25
Square scatterers, arrangement 1
Circular scatterers
20
Insertion Loss (dB)
15
10
-5
-10
500 1000 8000
Frequency (Hz)
Figure 4.14: FEM predictions of Insertion Loss spectra for 7x3 square lattice arrays of
the best performing triangular, square and elliptical rigid scatterer arrays compared to
that predicted for the reference cylindrical scatterer array.
Figure 4.14 compares the predicted IL spectra for the best performing
predicted to reduce the focusing effect (negative IL) and that with elliptical
scatterers in arrangement 1 offers the highest amplitude first Bragg band gap.
equivalent cylindrical arrays, the manufacturing cost for elliptical shapes will
be much higher. For example, it costs more than £1000 for 21 3 m long
Chapter 4 Page 84
4.6) Investigation of arrays with different lattice arrangements
This section will look at using a hexagonal lattice instead of the square
(1.5 kHz)
(a)
(1.5 kHz)
(b)
Figure 4.15: FEM (COMSOL®) predicted pressure maps at 1.5 kHz for 7x3 arrays of
cylindrical rigid scatterers with lattice constant of 0.135 m (a) hexagonal lattice and
(b) square lattice.
very different pressure map compared to the square lattice shape design as
shown in figures 4.15(a) and (b) respectively. This suggests that the
Chapter 4 Page 85
30
25
Hexagonal lattice
Square lattice
Insertion Loss (dB) 20
15
10
-5
500 1000 8000
Frequency (Hz)
Figure 4.16: FEM predictions of Insertion Loss spectra for 7x3 hexagonal and square
arrays of cylindrical rigid scatterers with lattice constant 0.135 m. FEM prediction are
performed with cylindrical wave source and a rectangular domain modelled as air.
The cylindrical wave source is located at 1.5 m away from the front face of the array
and the receiver is positioned 0.05 m from the nearest face of the array but on the
opposite side of the source.
lattice SC designs (lattice constant 0.135 m) are compared. For the hexagonal
lattice array, the first Bragg band gap is shifted up to a higher frequency (near
1.5 kHz) compared to that for the square lattice array. For a hexagonal lattice
array, the first Bragg band gap formula should be modified from that for a
cl
f Bragg _ Hex , (4.8)
3L
For these configurations taking account of the scatterer size and lattice
Chapter 4 Page 86
and 7 kHz (based on a single receiver position). However as a traffic noise
requires that outgoing wave propagate out towards infinity and therefore the
modelling enviroment), is critical for any numerical code concerned with the
condition and this has been used to treat these boundaries in the earlier
absorbing layer in which the wave equation has been modified with an
perpendicular to the interface with the physical domain (see figure 4.17). The
result is that waves entering the PML are absorbed only in the outgoing
direction, while the wave components tangential to the interface between the
PML and physical domain remain unaffected. The key property of a PML that
that waves incident upon the PML from a non-PML medium do not reflect at
the interface. The main advantages of the PML over radiation boundary
Chapter 4 Page 87
conditions are the relative ease of implementation using a complex
coordinates scaling [77]. Thus making it possible to minimise the size of the
7x3 array of circular scatterers with similar setup to that investigated in section
5.3 has been studied. The incident waves are modelled as plane waves due
(PML) (PML)
(PML)
(a) (PML)
1.2 kHz
(b)
Figure 4.17: Implementation of PMLs around a 7x3 square lattice array of circular
scatterers (a) location of PMLs and (b) FEM-computed pressure map at 1.2 kHz.
Chapter 4 Page 88
25
20
Measurement
FEM prediction, PML
15 FEM prediction, without PML
Insertion Loss (dB)
10
-5
-10
500 1000 8000
Frequency (Hz)
Figure 4.18: FEM predictions (with and without PML) compared to laboratory
measurements of Insertion Loss spectra for a 7x3 square lattice array of circular rigid
scatterers with a lattice constant 0.135 m.
As shown in figure 4.18, the FEM predicted IL spectra with and without PML
implementation are in fair agreement with the laboratory data. The time taken
to solve the problem using PML is 402.41 s which is 0.37 times faster than
that taken when using radiation boundary conditions with the maximum
Chapter 4 Page 89
Chapter 5
Measurement Techniques
5.1) Introduction
outdoors at the Open University’s noise barrier test site and at Diglis weir near
Worcester.
the acoustic space in which they are placed can be considered as linear and
Chapter 5 Page 90
time invariant and characterised by a system response ht [83]. A common
by a Digital to Analogue Converter (DAC) and fed to the system under test
with the unknown impulse response ht . The resulting output y t is sampled
together with the known signal to determine the sampled impulse response
ht . It is well known that the Signal-to-Noise Ratio (SNR) can be improved by
taking multiple averages of the measured output signal before the start of
DAC
ADC
PC
Chapter 5 Page 91
5.2) Laboratory measurments using Maximum-Length Sequence
System Analyzer (MLSSA).
[84]. MLSSA can measure and analyze many types of linear systems and is
widely used for audio and acoustics measurements [85, 86 and 87]. A MLS
characteristics of white noise (i.e. broadband and flat response). MLS may be
MLSSA is single channel but, given the periodic nature of the MLS, a single
channel analyzer with the advantage of doubling the useful bandwidth and
XOR
Figure 5.2: Linear feedback shift register for generation of a MLS (of length 2 4 - 1=
15) signal.
Chapter 5 Page 92
The number of samples of one period of an m order MLS signal is
analogue input using an audio patch cable as shown in figure 5.3. A single
and computes the first 16384 points of the impulse response. The time signal
reasonably flat up to about 25 kHz and then sharply rolls-off (see figure 5.5).
Analogue input
Patch cable
Analogue output
Figure 5.3: MLSSA system with audio patch cable connecting analogue output to
input.
Chapter 5 Page 93
0.8
System loopback test: MLS impulse response
0.7
0.6
0.5
0.4
Voltage (v)
0.3
0.2
0.1
-0.1
-0.2
10
-10
-20
Magnitude (dB)
-30
-40
-50
System loopback test: Power spectrum of MLS
-60
-70 2 3 4
10 10 10
Frequency (Hz)
Figure 5.5: Transfer function for system loopback test - MLS power spectrum.
Some ripples occur in the power spectrum (see figure 5.5) which
indicates that the Chebyshev filter transfer function is not perfectly flat.
Considering the small amplitudes of the ripple (about 0.3 dB), such deviation
Chapter 5 Page 94
5.2.2) MLSSA setup, data acquisition and analysis in
laboratory
principle) and thereby prevent reflections of sound waves [83]. They are also
insulated from external sources of noise which could influence the laboratory
results. Thus, the combination of these two factors can simulate a quiet open
largely on its size, the frequency of the signal used and also the size of the
object being tested. Scaled models are sometimes used to represent large
results shown in this work have been carried out in an anechoic chamber and
far away from any reflecting surfaces. The size of the anechoic chamber used
reflection from the human body or large measuring equipments, a small hole
is drilled on one side of the wall to allow cables to be routed through from the
remotely outside the chamber. Rock wool and polyurethane foam are used to
seal up the hole after laying the cable in order not to jeopardise the anechoic
Having tested the MLSSA system electronically using the loopback test
perform the actual measurement (see figure 5.6). The system analogue
Chapter 5 Page 95
output of the MLSSA is connected to a Cambridge Audio stereo A1 amplifier
(power output of 25 Watts) to drive the Brüel & Kjær (B&K) type 4295 point
size of the orifice and the shape of the loud speaker have been carefully
to a battery operated B&K type 5935 dual channel microphone power supply
(which also has amplifying function) for the measured signal to feedback to
the MLSSA system through the analogue input. Normally the microphone can
sound wave and therefore to increase the signal to a useable level, a high-
precisely, to make sure they’re very low in comparison with the wanted signal.
the loud speaker and microphone in place inside the chamber. Strictly
speaking, it is crucial that any reflecting object should be placed as far from
the set-up as possible. Care was taken to avoid unwanted reflection from
orientating them far away from the acoustic test space or by covering them
with sound absorbing material (see figure 6.17(a)). Although the anechoic
chamber was potentially fully anechoic (having a mesh floor grille over
absorbent material) these were found to interfere with the measurements. The
influence of the metal mesh floor grille is discussed in section 5.2.5. The
geometries used with specific scatterers (i.e. rigid, elastic shell, Split Ring
Chapter 5 Page 96
stated in the associated chapters. Supports for the 2 m long rigid and SRR
cylinder arrays were provided by holed wooden boards at the top and base of
the array. Figure 5.7(a) shows the bottom plate. A similar plate was used as
the top support. For the elastic shells, a bolt and nut securing systems was
used to mount the cylinders in place (see figure 5.7(b)). The 2 m long latex
sheets have been formed into cylinders with various outer diameters and wall
elastic shells were created by overlapping the edges of the latex sheets by a
few millimetres and gluing them together. To maintain their shape and vertical
Sonic crystal
Anechoic chamber
Chapter 5 Page 97
Elastic shell
5 mm diameter holes for elastic shells overlapped
Bolt
5 mm thick MDF
The direct sound field measurement method that has been used is
closely related with that described in BS EN 1793-6 [89]. During the array
nearest face of the array but on the opposite side of the source. Both source
and receiver were 1.2 m above the floor of the chamber. The loud speaker
was placed between 1.5 m and 1.63 m away from the array, such that the
source-receiver axis was normal to the array orientation (see figure 5.8).
Chapter 5 Page 98
y z
x x
d
1.5 m 1.5 m 50 mm
50 mm
L
(Refer to relevant chapters for d and L values). Ground
Figure 5.8: (a) Plan view of the source, receiver and array in the laboratory
measurements at normal incidence (b) the corresponding side view. Refer to chapters
2, 3, 4, 6 and 7 for the outer diameters and lattice constants used.
The input signals are first passed through the fixed analogue
frequency component beyond the selected range which is 10 kHz in our case.
The sampling rate was chosen to be greater than twice the maximum
using the MLSSA system, 16 time averages were carried out for each
measurement. The time required for each measurement was very short, only
a few seconds were necessary to obtain the impulse response for the
acoustic test space (see figure 5.8). Example plots of the time signals
corresponding to the direct field and total transmitted field through a sonic
crystal are shown in figure 5.9. In contrast with figure 5.4, the direct field
signal looks different to that in figure 5.4 since the complete system response
filter was used to frame the time signal and window out the unwanted
reflections. An FFT was applied to decompose the time domain signal to yield
Chapter 5 Page 99
its complex frequency components for frequency analysis with a frequency
higher number of FFT points. In our case, we are interested in the magnitude
of the measured signal (see figure 5.10). Insertion loss spectra have been
without (direct field) the cylinder array present as depicted in figure 5.11. The
spectra of signals extracted from impulse responses sampled with the same
calibration of the measurement chain with regard to the sound pressure level
-3
x 10
3
Direct field (MLS signal)
Total transmitted field (MLS signal)
2
1
Voltage (v)
-1
-2
-3
0 0.005 0.01 0.015 0.02 0.025 0.03 0.035 0.04
Time (s)
Figure 5.9: MLSSA time signals for direct and total transmitted fields.
-80
-90
-100
Direct field
Total transmitted field
-110
-120 2 3 4
10 10 10
Frequency (Hz)
Figure 5.10: MLSSA Frequency spectra for direct and total transmitted fields.
50
40
30
Insertion Loss (dB)
20
10
-10 2 3
10 10
Frequency (Hz)
Figure 5.11: Insertion Loss (IL) spectrum using MLSSA of 7x3 square lattice arrays
of rigid PVC scatterers of outer diameter 0.055 m and lattice constant 0.069 m.
Frequency resolution at 48 Hz.
chosen during the post-processing. A test was carried out to determine how
truncation of the time domain signal by the size of the window could affect the
power spectrum. The time domain signal of the total transmitted signal (see
and the end signal was varied from 10, 15, 20 and 30 ms corresponding to
window sizes of 3.8, 8.8, 13.8, and 23.8 ms respectively. The resulting power
spectra look different for each signal (see figure 5.12) in the low frequency
region, this is especially so for the smallest window size in which a dip
resulting from sampling discontinuities (i.e. integer number of periods are not
sampled by the FFT. This causes the energy contained in the signal to ‘leak’
from the signal frequency points into adjacent frequency points [90] (see
figure 5.12). This indicates that spectral leakage can occur if the length is not
properly chosen. Also worth noting is that the true frequency resolution or
is the duration of the measured signal used in the FFT calculation. This gives
window sizes of 3.8, 8.8, 13.8 and 23.8 ms. For the laboratory measurements
-60
-70
Magnitude (dB)
-80
-110
-120 2 3 4
10 10 10
Frequency (Hz)
Figure 5.12: Power spectra for different window size taken in calculating the FFT for
total transmitted fields.
the input to which is an integral extending over time. The implicit assumption
is made that the measured signal contained in the sampled block is repeated
to assign a weighting coefficient to each of the input signal and thus reducing
those samples that cause spectral leakage. In this section, two types of
window functions, namely the rectangular and half Blackman Harris were
tested. Half Blackman Harris is known to exhibit less spectral leakage than
rectangle window function [88]. A 4 term half Blackman Harris give excellent
side band rejection [83]. The comparison is made with both the windows
having the same size of 20 ms. In figure 5.13, the IL spectra shown are
calculated from the same direct field and total transmitted field signal (see
figure 5.9) for the two different window functions. The IL spectra for both
windows function does not differ much except for the spike at 2.5 kHz for the
windowing method with half Blackman Harris. This peak coincides with the
Bragg centre frequency for such lattice constant and since the spectrum does
not differ much for other frequency components, half Blackman Harris window
40
Half Blackman Harris window
Rectangular window
30
Insertion Loss (dB)
20
10
-10 2 3 4
10 10 10
Frequency (Hz)
Figure 5.13: Insertion loss spectra for 7x3 array of PVC scatterers of outer diameter
0.055 m and lattice constant of 0.069 m using different FFT filters.
As explained in section 5.2.2, the presence of the floor grille made from
steel in the anechoic chamber could affect the laboratory measurements. The
thickness of 0.25 mm has been measured with and without sections of the
floor grille. During the measurement with the floor sections removed, the
remaining grille edges and the steel beam (a permanent fixture to support the
mesh floor grille) were covered with sound absorbing materials as shown in
figure 6.17(a). The measured insertion loss spectra for the single latex
scatterer are shown in figure 5.14. Clearly the presence of the grille floor
6
Measurement, (no ground)
5 Measurement, (with ground)
4
Insertion Loss (dB)
-1
-2
-3 2 3
10 10
Frequency (Hz)
this section. Results of measurements made with a quarter-inch B&K free field
microphone type 4190 are compared with those obtained previously with the
single rigid (PVC) scatterer of outer diameter 0.11 m and wall thickness of
0.003 m with 4 symmetrical slits of size 0.012 m and the results are shown in
figure 5.15 According to these data the influence of the microphone size is
negligible.
14
10
Insertion Loss (dB)
-2 2 3 4
10 10 10
Frequency (Hz)
Figure 5.15: Insertion loss spectra for single PVC scatterer of outer diameter 0.11 m,
wall thickness of 0.003 m with 4 symmetrical slits size of 0.012 m.
processing for Excess Attenuation (EA) spectra are basically the same as for
IL spectra except that the EA is the reciprocal of Insertion Loss (see equation
elevated by using a table instead of using on the floor grille in the anechoic
chamber (see figure 5.16). In this way, the direct field measurement can be
easily taken by removing both the ground plane and table altogether without
changing the source and receiver locations. The measured EA spectra are not
influenced as much by the mesh floor ground (see section 5.2.5) since the
source to receiver distance for such measurement is shorter than the distance
Outdoor measurements have been made at the noise barrier test site
situated at the Open University Milton Keynes campus. This section will give
detailed descriptions of how the measurements (in situ) were performed using
composite scatterers (SRR4) as specified in table 7.2. The bottom ends of the
PVC pipes were capped and bolted to a metal panel which was flush with the
ground surface (see figure 5.17). The latex scatterers were inserted into the
PVC pipes from the other end and bolted to the top cover. Perspex clips for
each individual pipe were attached to the adjacent pipes at a height of around
2 m (see figures 5.18(a) and (b)). This reinforcement acted to maintain the
passing underneath the barrier. Conventional road side noise barriers panels
both sides of the SC. The measurement method used is closely related to the
conventional barrier panel using different sources and receivers (see figure.
5.19). Both the speakers and receiving microphones were placed 5 m away
from each face of the respective barriers (see figure 5.19) at a height of 1.6 m
above the ground. Another pair of microphones was used for reference sound
level measurements and they were placed at 1 m away from each source.
The speakers were located above non-porous asphalt (a car park) and the
invariance. Although this works well in laboratory (see section 5.2), the
wind speed and humidity [92]. Consequently a more robust modified white
noise measurement technique was used for the outdoor measurements. This
applies a traffic filter (represented by the road traffic noise spectrum specific in
(approximately 10 m away from the test site) from the noise barrier to reduce
data acquisition toolbox. The TN signal was fed into a high powered Mordaunt
during the measurement: one pair were B&K microphones and the other pair
(see figure 5.19). As a result of using a higher sampling period (20 s) in the
data. A Graphic User Interface (GUI) has been implemented in MATLAB ® for
Figure 5.17: Drilled Metal panels used for securing the scatterers.
clips
MDF
(a)
Measured
point mics
MDF
(b)
Figure 5.18: The outdoor noise barrier site, (a) source side of barrier showing SC and
conventional noise barrier sections and (b) rear side of barrier.
Fence
Speaker
1m 4m 5m
SC
1m 4m 5m
made from foam material was used. This gives a low resistance to the airflow
Figure 5.20: Graphical User Interface (GUI) Traffic Noise Analyzer implemented in
Matlab®.
Once the start button was pressed in the GUI, the measuring system
performed 12 repetitions using the TN signal and displayed the averaged time
spectrum for all of the 4 input channels. The buttons Imp_fig, Save_Imp,
Imp_fig2, and Save_Imp2 were then activated to save the time domain figures
and raw data from the input channels (see figure 5.21(a) and (b)). A DC off-
set was observed in some of the time signals (see for example figure 5.21(b))
FFT or FFT 2 button was pressed, transfer function estimates were performed
to yield the corresponding frequency component for the time domain signal by
comparing the measured signal with the input signal (see figure 5.22(a) and
(b)). The data analysis was somewhat different from the laboratory
measured free-field sound level, pRe ference _ field (reference microphone) was
adjusted to predict sound levels that would occur at the receiver positions in
the absence of the barrier based only on distance corrections. For example
the insertion loss with source distance correction ILDC for the SC based on
pTotal _ field
ILDC 20 log 20 log 5 Lsc 5 , (5.1)
p 1
Re ference _ field
Due to the high levels of ambient noise (the test site location is near to
a busy road), insertion loss results had to be corrected to take account of the
noise. This was done by recording the ambient noise (muting the source) and
performing the same procedure to obtain the power spectra of the ambient
noise (see figure 5.23(a) and (b)). In fact, this noise level is an approximation
as it did not represent the actual noise level during the barrier measurements.
Although the ambient noise levels were high, according to [10], the
measurement should have been immune to the noise as long as the signals
1.1
0.9
0.8
0 1 2 3 4 5 6 7
1
Voltage (v)
-1
-2
0 1 2 3 4 5 6 7
Time (s)
0.8
0.7
0.6
0.5
0 1 2 3 4 5 6 7
Time (s)
Time signal (Channel 4)
6
Fence, reference field (ACO mic)
4
Voltage (v)
-2
0 1 2 3 4 5 6 7
Time (s)
Figure 5.21: Example time domain signals measured outdoors (a) during SC barrier
measurements using B&K microphones and (b) during fence measurements using
ACO microphones.
Magnitude (dB)
100
80
60
40
20 2 3
10 10
Frequency (Hz)
FFT (Channel 2)
140
SC, reference field (B&K mic)
120
Magnitude (dB)
100
80
60
40
20 2 3
10 10
Frequency (Hz)
FFT (Channel 3)
(b) 140
120 Fence, total transmitted (B&K mic)
Magnitude (dB)
100
80
60
40
20 2 3
10 10
Frequency (Hz)
FFT (Channel 4)
140
Fence, reference field (B&K mic)
120
Magnitude (dB)
100
80
60
40
20 2 3
10 10
Frequency (Hz)
Figure 5.22 Example power spectra measured outdoors (a) during SC barrier
measurements using B&K microphones and (b) during fence measurements using
ACO microphones.
Magnitude (dB)
60
40
20
0 2 3
10 10
Frequency (Hz)
FFT (Channel 2)
100
Background noise at channel 2
80
Magnitude (dB)
60
40
20
0 2 3
10 10
Frequency (Hz)
FFT (Channel 3)
(b) 100
60
40
20
0 2 3
10 10
Frequency (Hz)
FFT (Channel 4)
100
80
Magnitude (dB)
40
20
0 2 3
10 10
Frequency (Hz)
Figure 5.23: Example background noise power spectra measured at the barrier test site
(a) during SC barrier measurements using B&K microphones and (b) during fence
measurements using ACO microphones.
Given that the required precision is high and that outdoor conditions
arises. One approach that has been used for measuring ground impedance is
to use more than two microphones. In that case, since more than one transfer
method [93].
subsequently correct the field data. This is a cheaper option but the question
of long term stability arises (i.e. whether the laboratory data is still valid under
example of how the calibration is done for the B&K microphones is briefly
discussed. To perform the calibration, the source was isolated and replaced
by a B&K sound calibrator type 4231 using the same NI acquisition system.
The microphones (mic reference no. 306 and 307) were individually tested by
inserting them into the calibrator as shown in figure 5.24. The calibrator emits
94 dB. Since the B&K dual-channel microphone power supply - type 5935
comes with selector switch to vary the gain, the calibration test was performed
gain levels are shown in figure 5.25. All the frequency spectra show the peak
at 1 kHz and for the corresponding gains level, the amplitudes of the sound
and 124 dB. Undesirable features in these spectra are the spikes that occur at
various frequencies. Since the spikes for both microphones occur at similar
However these spikes will cancel out after performing the correction
procedure. From these calibration results and assuming linearity for all other
frequencies, the differences between the microphones can be known and will
be taken into account in our outdoor measurement. Hence the final calculated
reference microphone signal, the background noise and the system response.
80
Magnitude (dB)
60
40
20
-20
2 3 4
10 10 10
Frequency (Hz)
Instruments type R3A-100) mounted about 2 m above the roof top of the
instrumentation cabin (figure 5.26). It is important to note that the U-axis of the
sonic anemometer is aligned to the North in order to determine the true wind
components along the three axes and the speed of sound (SOS)) which are
fed through the analogue inputs of the NI Data acquisition box. Wind vectors
were calculated from the U (x-axis) and V (y-axis) of the sonic anemometer.
Values for the W (z-axis) representing the vertical components were typically
from,
When the anemometer button on the GUI was pressed (see figure
5.20), the reading from the sonic anemometer was displayed as in figure 5.27.
Sonic
anemometer
Vector U
Vector U
0
-0.5
0 1 2 3 4 5 6 7
Vector V 0
Vector V
-0.5
-1
0 1 2 3 4 5 6 7
0.2
Vector W
0
Vector W
-0.2
0 1 2 3 4 5 6 7
Temp. ( C) Speed (m/s)
341
Speed of sound (m/s)
340.5
340
0 1 2 3 4 5 6 7
16
15
Temperature ( C)
14
0 1 2 3 4 5 6 7
Time (s)
Figure 5.28 shows an example of the wind vector plot (north is at 0°),
the black line indicates the barrier position and the red dotted line denotes the
direction of sound propagation. Blue solid lines calculated from vector U and
V show the wind speed according to this set of measurements which is about
0.8 m/s. The average speeds of sound and temperature for the entire 7
seconds (multiply by 4 repetitions) are also shown at the bottom of the wind
vector diagram.
0
30 330
60 300
Barrier
1 0.8 0.6
0.4 0.2
90
270
Direction of sound
propagation
120 240
150 210
which depending on the circumstances may prove better than others [94]. We
have considered another input source which uses the swept sine technique
(also known as chirp) developed by Farina [92]. The swept sine method uses
discussed in section 5.2. The MLS relies on the assumption of linear time
deconvolved impulse response when this condition is not met. It has been
found that the swept sine can overcome this limitation. The sweep duration of
sweep sine signal to avoid time aliasing. The electrically generated swept sine
1.5
Swept sine input signal
0.5
Voltage (v)
-0.5
-1
-1.5
0 1 2 3 4 5 6
Time (s)
single PVC scatterer of outer diameter 0.11 m and wall thickness of 0.003 m
measurements, 16 repetitions were made for the direct field and total
0.02
0.01
-0.01
Voltage (v)
-0.02
-0.03
-0.04
-0.05
Direct field, swept-sine signal
Total transmitted field, swept sine signal
-0.06
0 1 2 3 4 5 6
Time (s)
Figure 5.30: Direct and total transmitted fields measured using the swept sine method
in the laboratory.
into its frequency component and the Insertion Loss (IL) is calculated. The IL
spectra obtained using MLS and Swept sine methods are compared in figure
5.31 and show little difference. This swept sine method has also been
8.4.5.
12
MLS
Swept sine
Insertion Loss (dB) 10
-2 2 3 4
10 10 10
Frequency
New Music Award 2010 [78]. The project entitled “Organ of Corti” involved
sounds present in road traffic and falling water (i.e. weir). This device was
95, 96 and 97]. There has been increasing interest in application structures
with a low number of defects with respect to the total number of scatterers in
the structure, so that the periodicity is locally broken, to high precision wave
with lattice constant, L= 0.3 m (see figure 5.32(a)). This structure was placed
at Diglis Weir, river Severn, Worcester as part of the Worcester Music Festival
in August 2011 (see figure 5.32(b)). The cylinders were arranged in a quasi-
form a passageway. This allowed people to walk through the structure and
structure as well as while walking around it. For an SC with regular hexagonal
lattice arrangement of lattice constant 0.3 m, the calculated Bragg band gap
344 344
f Bragg _ Hex 660 Hz , (5.3)
3 L 3 0.3
This section will also provide some modelling results using Multiple Scattering
acting as the source. Time signal was recorded for 7 seconds with data
were made and for each set of measurement there were 16 repetitions (thus a
was used to obtain the frequency spectrum from the recorded time signal. The
figure 5.32(a) and table 5.1. During the experiment, there were some errors
made in the data acquisition system (due to crosstalk of signal) and thus, only
future measurements. The insertion loss (IL) was deduced by using the
the distance correction for the reference microphone, ambient noise and
Pt 2
1
Pt 1
-1
-3
0 2 4 6 8 10 12
(b)
Sonic crystal
weir
Figure 5.32: (a) Plan view schematic of the sonic crystal arrangement and microphone
locations at Diglis weir Worcester. (b) Aerial map of the site showing where the sonic
crystal was situated (picture taken from Google map).
Microphone locations x y
(coordinates)
Reference point (Ref pt) 9.7 1.69
Point 1 (Pt 1) 10.9 -1
Point 2 (Pt 2) 10.7 1.09
Point 3 (Pt 3) 11.6 1.69
Table 5.1: Coordinates of source and receivers positions during ‘Organ of Corti’
measurements.
30
Measurement Pt 1 (average of 3 measurements)
25
20
15
Insertion Loss (dB)
10
-5
-10
-15
-20 2 3
10 10
Frequency (Hz)
Figure 5.33: Mean IL spectrum (average of 3 measurements) with error bars for
receiver position Pt 1.
(see figure 5.33) shows a “modified” Bragg band gap at a frequency of 480
Hz. One contributing factor to the shifting of Bragg band gap could be the
effect of having a quasi ordered array rather than having a regular (full) lattice
(5.6.5). The spectrum shows additional attenuation peaks around 250 and
320 Hz which possibly suggest evanescent (or localised) modes [98] in the
between successive measurements (large error bars). The most likely cause
is other noise sources present in the vicinity. Nevertheless, the variation at the
peak of the “modified” Bragg band gap (480 Hz) is relatively small (< 6 dB).
(assumed line) source located at coordinate (0, 0) and (0, 10), see figure
(±0.05 m perturbation in both axes) from the original position (see figure 5.34).
1.
-0.96
-0.98
-1.02
-1.04
Perturbed position,
-1.06 (10.9, -1.05)
-1.08
10.8 10.82 10.84 10.86 10.88 10.9 10.92 10.94 10.96 10.98
Figure 5.34: The original receiver position Pt 1 and 4 other perturbed positions
made for MST modelling.
20
10
-10
-20 2 3
10 10
Frequency (Hz)
Figure 5.35: Predicted IL spectra of all the 5 individual positions (Pt 1 original and 4
perturbed positions) and the averaged IL spectrum. Source is located at coordinates
(0, 0).
Figure 5.35 shows the predicted IL results for all the 5 positions (Pt 1
original and 4 perturbed positions) and the averaged IL spectrum. The source
noticeable when there is a slight change in the receiver position but the
0.95) where almost 20 dB difference is observed near 450 Hz. This narrow
frequency point.
20
10
-10
-20 2 3
10 10
Frequency (Hz)
MST predicted (average of 5 receiver points) IL spectra with the point source
modelled at location (0, 0). The predictions agree fairly well with the measured
data at low frequency. The first peak (200 Hz) in the predictions appears to be
shifted in frequency. The second and thirds peaks are similar to those
measured. The third peak corresponding to the ‘modified’ Bragg band gap is
predicted well. One explanation for the negative IL shown in the measured
data (low frequencies) could be the effect of the ground (see chapter 8).
20
10
-10
-20 2 3
10 10
Frequency (Hz)
Figure 5.37: Predicted IL spectra at all the 5 individual positions (Pt 1 original and 4
perturbed positions) and the averaged predicted IL spectrum. The point source is
assumed to be located at coordinates (0, 10).
Figure 5.37 shows the predicted IL results for all the 5 positions (Pt 1
original and 4 perturbed positions) and the averaged IL spectrum. The source
spectra are observed when there is a slight change in the receiver position as
in figure 5.35.
20
10
-10
-20 2 3
10 10
Frequency (Hz)
point at location (0, 10). The prediction agrees slightly better at high
coordinates (0, 0) see figure 5.32(a). Three peaks are again predicted at 290,
380 and 500 Hz. The Bragg band gap is well predicted in terms of the centre
The fact that predictions using a single fixed point source at (0,0) are in
better agreement with the measured data than those assuming a source at
(0,10 suggests that using a line of point sources could improve the accuracy
of the sources are shown in figure 5.39 and the corresponding coordinates
are shown in table 5.2. The MST prediction was run for each source location
and receiver position, Pt 1 (see table 5.1). Similarly, averaging of all the 5
points (receiver Pt 1 and 4 perturbed positions) was carried out (see figure
5.34). It is reasonable to assume that all the point sources are coherent since
the weir noise is caused by more or less the same flow and drop over the
length of the weir. The direct pressure (without SC) in Decibel, p Direct (dB) was
n pDirect ( dB)
p Direct _ total( dB) 10 log 10 10 10 . (5.4)
0
Likewise for the total transmitted pressure (with SC) in Decibel, Ptotal_ tranmitted( dB)
Figure 5.39: Picture illustration for the location of the 10 sources (S1 to S10) and the
position of sonic crystal (photo taken from Google map).
Source no. x y
S1 -7.82 10.92
S2 -16.73 6.38
S3 -25.58 1.72
S4 -34.11 -3.55
S5 -42.83 -8.45
S6 -50.92 -14.47
S7 -58.93 -20.6
S8 -67.15 -26.35
S9 -75.63 -31.65
S10 -83.29 -38.29
Table 5.2: Coordinates of individual sources in a multiple point source model of the
weir.
40
Measurement, average of 3 measurements, point 1
MST Prediction, worcester original design, 10 sources
30
Insertion Loss (dB)
20
10
-10
-20 2 3
10 10
Frequency (Hz)
Figure 5.40: Measured (average of 3 measurements) and MST predicted with multiple
point sources IL spectra for receiver position Pt 1 (averaged over 5 receiver points for
each individual source).
receiver position 1. The predictions resulting from the 10 source model of the
weir show better agreement with measurements than those assuming a single
point source at either location (0,0) or (0,10) (see figures 5.36 and figure
5.38). The first IL peak around 250 - 350 Hz is well predicted and this could
be due to the evanescent mode suggested in figure 5.33. Also, the prediction
for the modified Bragg band gap (possibly shifted due to the quasi periodic
ordered array effect) near 450 Hz agrees fairly well with the measured data.
6.1) Introduction
transmission is to increase the thickness or the mass per unit area of the
results in large size and mass of the insulator. On the other hand the stop and
pass bands of a Sonic Crystal (SC) composed of solid rigid elements and their
has been previous numerical and theoretical work in respect to periodic arrays
cylindrical balloons in air. In these systems, the authors analysed the band
the latex wall thickness, the filling fraction and the pressure inside the
balloons. It’s reported that the pressure inside the balloon is kept at 1.1 atm
which is slightly above the ambient pressure and the latex wall thickness is
such matrix. It has been found that resonance attenuation peaks can be
to predict the measured effects due to the pre-stress condition i.e. to take
account of any static external load and thereby model scattering by elastic
variation of the radius through each cycle. Vibration analysis for the natural
complex than that for a beam or plate because the generality of the shell
characters. The modal shapes of the cylinder model are basically composed
of two kinds of mode patterns: the breathing mode and the freely supported
beam mode [103]. The combination of the two basic mode patterns support
three modal shapes – axial, bending and torsional modes (see figure 6.1). It is
possible to realise such a mode provided that the thickness of the cylinder
(Young’s Modulus) is of the order of 1 MPa [104]. In this chapter, the theory of
It is a semi-analytical method and in this work the terms ‘MST’ and ‘semi-
mode resonance is also verified by performing modal analysis with the Finite
used to predict the transmission problem for a single elastic shell. These
available non-vulcanized rubber (latex) and data confirm the presence and
identity of the resonances for the chosen material and geometry of the elastic
(b)
(c)
(d)
Figure 6.1: Example predictions of the two basic mode shapes of a freely supported
cylinder – (a) radial (b) beam modes up to fourth order (c) the combination of the first
radial mode (mode 0) and up to the fourth order of bending modes (d) the
combination of the first radial mode (mode 1) and up to the third order of axial modes
[103].
by our collaborators [111] will predict the results of the use of thin “soft” elastic
are the compressional and shear wave speeds of the elastic material
array. These stop bands are independent from those related to the periodic
results, the scattering problem is solved for a single scatterer in this section
where j takes positive integer values, each shell can be replaced by its cross-
section (i.e. an elastic ring), see figure 6.2. The acoustic environment
o c
outside/inside of each shell is described by density and sound speed o ,
i ci
Figure 6.2: Cross-section of an elastic shell in the primary cell of doubly periodic
array.
Without loss of generality, the origins of the Cartesian x, y and polar
coordinates r , used in this section coincide with the centre of the scatterer
C0 in the primary cell of size L, also known as the lattice constant. Centre O j
of each scatterer in the infinite periodic array is defined by the position vector
R j n1a1 n2 a2 ; where n1 and n2 take integer values and a1 and a2 are the
fluid medium previously shown in chapter 4, equation 4.1 is repeated here for
2 p k2 p 0 , (6.1)
i
pr R j e pr .
T
Rj
(6.2)
2 2 2f
Where 2 is the Laplacian, k is the wave number in air
x1 x2
2 2
c
defined as the ratio between angular frequency and sound speed of the
density to one of the acoustic media (i.e., it equals to either “o” or “i” for
2 p , (6.3)
i p , (6.4)
2 k12 0 , (6.5)
2 k 22 0 , (6.6)
where k1 and k 2 .
c1 c2
r ao and r ai :
rr 2 p , (6.7)
p 1
. (6.9)
r r r
2 c22
2
rr r
r 2 ,
2 (6.10)
r2 r r
2 c22
2
r r
r 2 , (6.11)
r2 r r
2
orthogonal series of Bessel and trigonometric functions. In the outer and inner
po r , A J k r B Y k r e
n
n n o n n o
in
, r ao , (6.12)
pi r , E J k r e
n n i
in
, r ai , (6.13)
n
r , D
n
1, n J n k 2 r D2, nYn k 2 r e in , ai r ao . (6.15)
2
(6.15) into equations (6.7) - (6.9) and taking inner product e im d will
0
An
a J k a , inJ
o
'
n o o n ko ao , aoYn' ko ao , inYn ko ao n Yn' ko ao B ,
(6.16)
ao J n' ko ao J n' ko ao
n
T
2 o o k o2
Cn n 0, 0, , 0 Bn , (6.17)
aoJ n k o ao
'
where ao outer boundary of the shell, J n' = Bessel function of n-th order of the
2f c dZ k r
(the wave number in air), , , Z n' k r n is the
c
ko
c c2 c2 dr
form,
' 2 k2 Z n k a
2Z n k 2 a
in
f 2, , Z Z k a , (6.20)
Z n k a
' n 2
a a
1 ' 2n 2 k 22 a2
g 2, , Z 2Z n k 2 a Z n k 2 a . (6.22)
a a
In the case of the rigid scatterer equation 6.16 is reduced to the well-known
The contrast between elastic material of the shell and air medium considered
impedance,
o co
1 . (6.23)
c2
radius R
ao ai , then h 1 , where h is half Wall Thickness (WT) of the
2 R
shell.
In order to have better interaction between elastic and acoustic media, the
and mid-surface radius of the elastic shell. This assumption can be written as
.
positioned in the outer region of the shell. For convenience, but without loss of
generality, the origin of the Cartesian x, y and polar coordinates r,
coincides with the source location. Potential po r in the outer acoustic
po
r
ikpo 0 r 1/ 2 , as r , (6.24)
po, s r , An Z n H n1 k o rˆ e in , rˆ ao ,
ˆ
(6.25)
n
so that the total field in the outer acoustic medium takes the following form,
J n' k o R
Zn . (6.27)
H n1 k o R iUˆ 1, n
Note that when Uˆ 1, n 0 factors Z n transform to those for the case of a rigid
cylinder.
Solution for the field in the inner acoustic medium can be described by the
formulation in equation 6.7. The following will describe the formulation of the
approximation for the wave equations of the thin-walled elastic shell based on
the previous assumption regarding the thickness and relative impedance can
be formulated.
The problem stated by equation 5.5 and 5.6 and boundary conditions in
equation 5.7 to 5.9 can be reduced to the asymptotic equations given by:
1 2u2 u1 2
k3 u 2 0 , (6.29)
R 2 2
co
components in the elastic shell, k3 and ko . The p-wave speed c3 for
c3 c2
E
c3
, (6.30)
1 2
loading applied to the walls of the shell which in our case is a membrane
Equation 6.28 and 6.29 are solved in conjunction with the conditions imposed
po r , pi r ,
, (6.31)
r rR r r R
po r ,
u1 . (6.32)
r r R
The displacements of the elastic shell can be expanded into the Fourier
series:
u2 u2,n e in , (6.34)
(Appendix B) so that all the terms can be expressed in terms of the same
J k rˆH k Q e
po rˆ,ˆ in
An Z n H n1 ko rˆ ein , rˆ a0 ,
1 ˆ
n o n o (6.35)
n
within which vector Q Qcos , sin is the radius vector to the centre of the
shell.
and 6.29, boundary conditions 6.31 and 6.32, and the modified solution 6.35.
After solving the appropriate algebraic system, coefficients An are derived as,
Et Eˆ 1 f t , (6.37)
t
f t Eˆ j 1 e j ,
(6.38)
j
relaxation time j .
Applying the Laplace transform
0
e st dt , one can derived the Young’s
E s Eˆ 1 sF s , (6.39)
Eˆ j
F s
s1 s
. (6.40)
j j
Eˆ j
E Eˆ 1 ,
1 i j
(6.41)
relationship by taking the Fourier transform
0
e i t dt .
We use viscoelasticity data measured for the material with properties similar
to latex used in our experiment [108]. Instead of equation 6.41, this requires
use of:
N
E Eˆ Eˆ j and E Eˆ Eˆ j . (6.43)
j 1
where the Young’s modulus corresponds to the equilibrium state and values
of and are taken for an elastomer with properties similar to latex [108] (see
table 6.1).
Relaxation time j (s) Dynamic Young’s Moduli, E j N / m
2
0 1.745361577x106
0 2.13x106
4.32x10-9 9.00x106
5.84x10-8 4.20x106
3.51x10-7 2.94x106
2.28x10-6 2.41x106
1.68x10-5 1.87x106
2.82x10-4 1.31x106
7.96x10-3 7.02x105
9.50x10-3 4.45x105
Table 6.1: Relaxation time and dynamic Young’s Modulus.
the same ( o = i = 1.25 kg/m3 and co = ci = 344 m/s). Also the latex
material properties used for prediction are kept constant (except in section
0.4997. Tensile testing to yield the E of latex was performed with the kind
assistance of the Materials Department. The test procedure and results are
provided in Appendix D and data confirmed that the E value is consistent with
the approximation used. Initially, the prediction is made for a latex cylinder
cylinders as well. The results shown in figure 6.3 are computed in terms of
H o1 ko r
IL 20 log10 . (6.44)
p0
3
Insertion Loss (dB)
-1
-2
-3 2 3
10 10
Frequency (Hz)
Figure 6.3: Measured (solid black line) and MST predicted (broken blue line) IL
spectra for single latex scatterer of outer diameter (OD) 0.055 m, wall thickness (WT)
0.00025 m and length of 2 m. For comparison the IL spectrum measured for an
acoustically “rigid” Poly Vinyl Chloride (PVC) pipe with similar diameter and length
is shown also (broken red line).
In figure 6.3, the measured and MST predicted IL spectra for a single
latex scatterer of Outer Diameter (OD) 0.055 m and Wall Thickness (WT)
for a single Poly Vinyl Chloride (PVC) pipe of similar length (2 m), outer
diameter (0.055 m) and WT (0.002 m) is shown for comparison with the latex
scatterer. Clearly, the latex scatterer produces a resonance peak around 1.1
kHz which is not observed in the data for the PVC scatterer. The MST
resonance peaks at around 1.1, 4 and 6.1 kHz. A negative IL peak between
500 Hz and 1 kHz is observed just below the latex resonance. The narrow
band higher order resonance mode that appears in the predicted IL spectrum
at around 6 kHz is not visible in the data. This might suggest that such narrow
are determined by the material properties (i.e. mass, damping and stiffness)
working with eigenfrequencies, f eigen and working with eigenvalues, eigen . The
meigen
f eigen . (6.45)
2
The 2-dimensional simulation is first performed for an elastic shell with similar
geometry and material properties stated in section 6.2 (figure 6.3). This is
performed under vacuum for the inside and outside environments of the
(c) (d)
(e) (f)
deformation) looks similar to mode 0 in figure 6.1 (a). The breathing mode
and 0.0025 m for the volumetric expansion and contraction respectively. This
coupling and this must be taken into account. The important characteristic is
the interaction between the two media, that is, the elastic structure will
influence the air flow resulting from changes in the distribution and the size of
the air load. It is imperative that the solution of the single elastic shell takes
account of the presence of air and thus the additional complexity of modelling
this model, the boundary conditions are done in a more rigorous manner
sensitive to the shape of the wave fronts. The PML formulation introduces a
(c) (d)
plots of the latex breathing mode results after using perfectly matched layers
5924
eigenfrequencies using equation 6.45. For example gives an
2
eigenfrequency of 942 Hz. Clearly with the introduction of air damping, the
resonance frequency also agrees fairly well with the semi-analytical and
results show that the resonance is indeed a breathing mode type resonance
studied both in laboratory and using MST. Three different types of diameter
(0.0275, 0.043 and 0.055 m) with similar wall thickness of 0.00025 m are
investigated.
-1
-2 2 3 4
10 10 10
Frequency (Hz)
6
(b)
5
Measurement, 0.043mm (OD), 0.00025 m (WT)
MST prediction
4
Insertion Loss (dB)
-1
-2 2 3 4
10 10 10
Frequency (Hz)
6
(c)
5
Measurement, 0.055 m (OD), 0.00025 m (WT)
MST Prediction
4
Insertion Loss (dB)
-1
-2 2 3 4
10 10 10
Frequency (Hz)
Figure 6.6: IL spectra of measurements and MST predictions for different outer
diameter 0.0275, 0.043 and 0.055 m in figure (a), (b) and (c) respectively.
Figure 6.6 shows the measured and MST predicted IL spectra for
0.055 m are observed to be near 1.8, 1.5 and 1.1 kHz respectively. This
clearly shows that the “breathing mode” resonance can be tuned to lower
frequency by increasing the diameter of the elastic shell scatterer. For all
three cases, the MST predictions agree well with the experimental results. It is
also noted that the accuracy of the MST prediction for the IL amplitude
using MST are performed for 2 different diameter sizes of elastic shell 0.043
and 0.055 m. For each of the different diameters, wall thicknesses of 0.00025
(figure 6.7(a) and (c)) and 0.0005 m (figure 6.7(b) and (d)) are used.
(a) 5
(b) 5
Measurement, 0.043mm (OD), 0.00025 m (WT) Measurement, 0.055 m (OD), 0.00025 m (WT)
MST prediction MST Prediction
4 4
Insertion Loss (dB)
2 2
1 1
0 0
-1 -1
-2 2 3 4
-2 2 3 4
10 10 10 10 10 10
Frequency (Hz) Frequency (Hz)
6 6
(c) 5
(d) 5
Measurement, 0.043 m (OD), 0.0005 m (WT) Measurement, 0.055mm (OD), 0.00075 m (WT)
MST prediction MST prediction
4 4
Insertion Loss (dB)
2 2
1 1
0 0
-1 -1
-2 2 3 4 -2 2 3 4
10 10 10 10 10 10
Frequency (Hz) Frequency (Hz)
Figure 6.7: Measured and predicted IL spectra for outer diameter of 0.043 m (left
panel) and 0.055 m (right panel) with two wall thicknesses of 0.00025 ((a) and (b))
and 0.0005 m ((c) and (d)) in each case.
single elastic shell scatterer for 2 different diameters of elastic shell 0.043 m
(left panel) and 0.055 m (right panel). In the case of an elastic scatterer with
resonance frequency for both experiment and MST prediction is shifted from
near 1.5 kHz down to 900 Hz as the wall thickness is increase from 0.00025
figures 6.7(b) and 6.7(d)). The accuracy of the MST prediction is again
observed to improve with increase in either the diameter or the wall thickness
increasing the diameter or the wall thickness of the elastic shell might be due
reported using the material parameters listed in table 6.2. The table specifies
the longitudinal and transverse wave speed, density and Poisson’s ratio for
the elastic material. The upper, lower limit and the median for elastic shell are
calculated.
5
Measurement, 0.043 m (OD), 0.0005 m (WT)
MST prediction, rho 920, c1 40, c2 11, poisson ratio 0.48
4
MST prediction, rho 1010, c1 95, c2 22, poisson ratio 0.485
MST prediction, rho 1100, c1 150, c2 29, poisson ratio 0.49
Insertion Loss (dB)
-1
-2 2 3 4
10 10 10
Frequency (Hz)
frequencies for minimum, median and maximum stiffness values for latex are
near 930, 880 and 840 Hz respectively. Increasing the material density is
result indicates the robustness of the breathing mode resonance for such
material where only a small deviation (about 100 Hz) of the resonance
frequency results from the change in latex density between lowest and
influence of the 3 dimensional aspect of the latex cylinder (i.e. the finite height
the method of end support has been investigated empirically by varying the
tension at the ends of the latex cylinders between ‘loose’ and ‘tight’. To
quantify the tension, the numbers of turns on the securing nuts are recorded.
First, the single latex cylinder is mounted with ‘normal’ tension. Then a further
15 turns of the securing nuts (top and bottom of supporting plate) has been
used to give ‘tight’ tensions and 15 less turns than required for ‘normal’
2.5
Measurements, 0.0275 m (OD), 0.00025 m (WT)
2 Tightly tension
Normal tension
1.5 Loosely tension
Insertion Loss (dB)
0.5
-0.5
-1
-1.5
-2 3 4
10 10
Frequency (Hz)
(‘tight’, ‘normal’ and ‘loose’). The tension is controlled by the 2 mounting nuts
at both top and bottom of the latex cylinder. The breathing mode resonance at
around 1.7 kHz is clear for the ‘tight’ and ‘normal’ tensioning states. However
lattice with square cells of size L (figure 6.2). Assuming that a cylindrical
where M is the number of scatterers in the array, the first sum is taken over all
numbers of shells in the array, the factors Z nm take the previously derived
from equation 6.27 specified to the scatterer with index m , and variables
rˆm r , and ˆm r , are the polar coordinates with origin in the center of shell
of index m . Note that the solutions for the shells and the inner acoustic
n
p p 1
(6.48)
p 1, p m q
, m 1,..., M ,
Use of the modified solution 6.48 in equations 6.28, 6.29, 6.31 and 6.32
together with the solutions of inner acoustic medium and elastic shell makes it
N
H q1n koQmp e
M
H n1 koQm ein m ,
i q n mp
Anm A Z q
p
q
p
p 1, p m q N (6.49)
m 1,..., M ,
M(2N+1) of equations, and then can be used to obtain the coefficients Anm ,
35
30
Measurement, 7x3 elastic shell and L = 0.08 m
MST prediction
25
Insertion Loss (dB)
20
15
10
-5
100 1000 5000
Frequency (Hz)
Figure 6.10: Measured and MST predicted IL spectra for an array of latex shell
scatterers of outer diameter (OD) 0.055 m and wall thickness (WT) 0.00025 m with
lattice constant 0.08 m.
30
Measurement, 7x3 elastic shell and L = 0.1 m
MST prediction
Insertion Loss (dB) 25
20
15
10
-5
100 1000 5000
Frequency (Hz)
Figure 6.11: Measured and MST predicted IL spectra for an array with lattice constant
0.1 m of latex shell scatterers having outer diameter (OD) 0.055 m and wall thickness
(WT) 0.00025 m.
30
25
Measurement, 7x3 elastic shell and L = 0.15 m
20 MST prediction
Insertion Loss (dB)
15
10
-5
-10
100 1000 3000
Frequency (Hz)
Figure 6.12: Measured and MST predicted IL spectra for an array with lattice constant
0.15 m of latex shell scatterers having outer diameter (OD) 0.055 m and wall
thickness (WT) 0.00025 m.
In figure 6.10 – 6.12 show the measured and MST predicted IL spectra
for a 7x3 square lattice array with lattice constants of 0.08 m, 0.1 m and 0.15
m respectively. Each latex shell used in all three array configurations has OD
0.055 m and WT 0.00025 m (similar to the single scatterer use in figure 6.3).
(chapter 2) for lattice constants of 0.08 m, 0.1 m and 0.15 m are 2.1 kHz, 1.7
kHz and 1.1 kHz respectively. The predictions show relatively good
agreement with measured data. The peaks near 1 kHz in figures 6.10 and
6.11 can be identified as the result of the resonant behaviour of the elastic
amplitude between theoretical and experiment results for the resonant peak
may be explained by the differences between assumed and actual source and
between data and predictions. The peak at about 750 Hz could be due to the
y z
(a) (b)
x x
d
1.5 m 1.5 m 50 mm
50 mm
1.2 m
L
L
(Lateral displacement up to 14 cm) Ground (Vertical displacement up to 45 )
Figure 6.13: (a) Plan view of the source, receiver and array in the laboratory
measurements at normal incidence and microphone laterally displace (see red arrow)
up to 14 cm (b) the corresponding side view with microphone vertically displace (see
red arrow) up to 45 from normal incidence.
25
Lat. displ = 0
1 cm
Insertion Loss (dB)
20 2 cm
3 cm
5 cm
15 14 cm
10
-5
300 1000 5000
Frequency (Hz)
Figure 6.14: Lateral displacement, measurement IL spectra for array with lattice
constant 0.08 m of latex scatterers having outer diameter (OD) 0.055 m and wall
thickness (WT) 0.00025 m.
Figure 6.14 shows the IL spectra measured for 7x3 square lattice array
figure 6.13(a). The constancy of the low frequency peaks around 1 kHz
receiver angle. The Bragg band gap central frequency for this lattice
configuration is around 2 kHz, similar to that for a rigid scatterer array with the
same dimensions. It seems that the classical band gap is not angle
was not possible without introducing diffraction from the array edges. The
peak at around 700 Hz attributed to the finite height of the array (refer to FEM
35
Measurement, 7x3 array,
30 elastic shell, L = 0.1 m
Vert. Displ = 0
25 13
27
Insertion Loss (dB)
20 45
15
10
-5
-10
300 1000 5000
Frequency (Hz)
Figure 6.15: Lateral displacement, measurement IL spectra for array with lattice
constant 0.15 m of latex scatterers having outer diameter (OD) 0.055 m and wall
thickness (WT) 0.00025 m.
Figure 6.15 shows the measured IL spectra due to a 7x3 square lattice
array with lattice constant of 0.1 m of latex scatterers having OD 0.055 m and
between 1 and 2 kHz as the elevation angle is changed from 0 to 45. One
possible explanation is that the breathing mode for the latex shell is supported
only by waves incident in the plane normal to the cylinders. A second possible
reason could be the influence of reflections from the end plates of the
microphone was not possible to avoid diffraction from the base of the
supporting structure. These results also support the argument that the peak at
results from the use of FEM. The transmission problem for a single latex
interaction model for an elastic shell and air both inside and outside of the
represents the elastic shell (figure 6.16 (a)). The domain boundaries are
parallel to the Ox axis of the elastic shell. The boundary condition at the
has been used and inside the domain the element size is at least 6 times
solver is run as explained in section 4.2 (chapter 4). The predicted pressure
plot is plotted at 1050 Hz for the exterior environment (figure 6.16 (b)). Clearly
particular frequency.
Figure 6.16: Meshing of single latex with OD 0.055 m and wall thickness (WT
0.00025 m (a). FEM predicted pressure plot at 1050 Hz for the exterior environment
of the model (b).
5 3.8 dB
Measurement 0.055 m (OD) and 0.00025 m (WT)
MST Prediction
4 FEM prediction 3.2 dB
3
3.6 dB
Insertion Loss (dB)
-1
-2
-3 2 3
10 10
Frequency (Hz)
Figure 6.17: Measured and predicted (FEM and MST) IL spectra for a single latex
scatterer having 0.055 m (OD) and wall thickness 0.00025 m (WT).
MST predictions and the measurement data for the single latex scatterer
(0.055 m (OD) and wall thickness 0.00025 m (WT)). For all measurements
and predictions, the receiver is positioned 0.05 m away from the scatterer.
The IL values at 1 kHz (breathing mode resonance frequency) for all the 3
spectra are shown in the figure. The discrepency between the prediction and
measurement at the Bragg frequency is about 0.6 dB. Higher order modes
are observed by FEM as well as MST but due to the narrow band
0.00025 m thick) have been formed into cylinders with outer diameter of
0.0275 and 0.055 m, respectively. This has been achieved by folding each
millimetres (< 0.01 m) and gluing them together using cyanoacrylate. Perspex
circular end caps with metal bolts were made for both ends of the latex
cylinder and used to attach them to holed wooden boards at the top and base
of the cylinders. To maintain their shape and vertical orientation, the latex
common pipe driven by a small pump. Nevertheless some deviations from the
pipe
bolt
(a) (b)
Figure 6.18: Measurement setup for array of cylinders of 7x3 arrays of latex (a) and a
single latex showing the overlapping strip, mounting bolt and pipe for air inlet (b).
table 6.4). The resulting mean diameter and standard deviation are 0.05732 m
section 6.11. This simultaneously tested the influence of the metal beams at
the side of the supporting structure (see figure 6.18a). During the test the
overlapping and glued section of each cylinder was positioned to face the
6
MST Prediction (OD = 0.057 m)
5 MST Prediction (OD = 0.055 m)
Measurement with error bar (averaged of 21 cylinders)
4
Insertion Loss (dB)
-1
-2 2 3 4
10 10 10
Frequency (Hz)
Figure 6.19: Averaged of 21 measurements of IL with error bars (indicating the range
of data at each frequency) and MST predicted IL spectra for a single latex scatterer of
(outer diameter (OD) 0.055 m and wall thickness (WT) 0.00025 m) (blue broken line)
and (OD 0.057 m and WT 0.00025 m) (red solid line).
comparison with the MST prediction are shown in figure 6.19. The agreement
between the measured and predicted IL spectra indicates that the latex
Given the small variation in the diameter (OD) of latex cylinder in the
The elastic shell with glued section has a non-uniform cross section
both in respect of thickness and stiffness which may have a significant effect
the shell has been modelled, using COMSOL® Multiphysics 3.5, as an elastic
boundary conditions imposed at the acoustic interface and fixed faces (zero
replaces the radiation condition and as such the domain size can be reduced.
elements and the active domain is meshed by at least six elements per
64-bits windows platform with 3.07 GHz Intel core i7 processor and 6 GB of
RAM. The estimated computational time for solving the full problem in the
simulation has been performed by our collaborators [111]. Also note that in
contrast with their simulation, the simulation done here is performed with a
cylindrical wave rather than a plane wave. The results are found to have
scatterers,
3
Insertion Loss (dB)
-1
-2
-3 2 3
10 10
Frequency (Hz)
Figure 6.21: Measured and FEM-prediction IL spectra due to a single latex cylinder
(with and without rigid inclusion) of outer diameter (OD) 0.055 m and wall thickness
(WT) 0.00025 m
5
Measurement, 0.0275 m (OD), 0.00025 m (WT)
FEM prediction without rigid inclusion
4 FEM prediction with rigid inclusion of width 1 cm
3
Insertion Loss (dB)
-1
-2
-3 2 3
10 10
Frequency (Hz)
Figure 6.22: Measured and FEM-prediction IL spectra due to a single latex cylinder
(with and without rigid inclusion) of outer diameter (OD) 0.0275 m and wall
thickness (WT) 0.00025 m
without rigid inclusion) are due to a single latex cylinder (OD 0.055 m and WT
very small. The first peak which corresponds to the breathing mode is
also agree well with the measured data except for over prediction of the IL
level in the case of the smaller latex cylinder. Nevertheless, in the frequency
regions where the breathing mode there is good agreement between data and
predictions. This study supports the hypothesis that increasing the material
experiments were delivered in 0.7 m lengths due to the limitating size of the
cylinder, a Perspex disc of slightly larger diameter was used to join the
individual latex cylinders together (see figure 6.23). This improvised approach
also served to restore and maintain their cylindrical shape as the latex
along the scatterers were observed to be more elliptical than cylindrical. The
difference is more prominent at the “extreme flat region” which has the largest
semi-major axis. Another anomalous region is the “inflated region”, about 0.05
m deviation from the joining disc, where the cylinder starts to deviate from the
cylindrical shape maintained by the disc. The average dimensions taken from
5 cylinders for these two regions are shown in the insert in figure 6.23.
the breathing mode of the introduction of rigid disc and the changes in shape
“extreme flat and inflated” regions are denoted as ‘edge’ (see figure 6.23).
Perspex
disc,
“Rigid”
4
S1
Extreme flat face
Extreme flat edge
3
Insertion Loss (dB)
Rigid (joint)
Inflated flat edge
Inflated flat face
2
-1 2 3
10 10
Frequency (Hz)
Figure 6.24: IL spectra of single industrial latex (Sample 1, S1) at all 5 positions
(Extreme flat face/edge, rigid and inflated face/edge).
latex cylinder at all 5 positions illustrated in figure 6.23. Clearly, the breathing
mode resonance near 900 Hz is preserved in all the positions despite the non-
uniformity of the scatterer. This indicates that the breathing mode resonance
5 3
Extreme flat face Extreme flat edge
2.5
4 S1 S1
S2 S2
2
Insertion Loss (dB)
0
0
-1 -0.5
-2 2
(a) -1 2
(b)
3 4 3 4
10 10 10 10 10 10
Frequency (Hz) Frequency (Hz)
6 3.5
5 3
Inflated edge
Inflated face
2.5 S1
4 S1
S2 S2
Insertion Loss (dB)
1.5
2
1
1
0.5
0 0
(c) (d)
-1 2 3 4
-0.5 2 3 4
10 10 10 10 10 10
Frequency (Hz) Frequency (Hz)
3.5
Rigid
3
S1
2.5 S2
Insertion Loss (dB)
1.5
0.5
-0.5
-1 2
(e) 3 4
10 10 10
Frequency (Hz)
Figure 6.25: IL spectra of consistency test for 2-single industrial latex (Sample 1 (S1)
and 2 (S2)) at all 5 positions (Extreme flat face (a) or edge (b), inflated face (c) or
edge (d) and rigid (e)).
Figure 6.25 shows the IL spectra comparing two single industrial latex
extremely flat face (b), an extremely flat edge (c) an inflated face (d) an
inflated edge and (e) rigid. The results show that the breathing mode
The deviation between the IL spectra at higher frequency could be due to the
arrays of elastic cylinders have been constructed and their performance has
35
30
Extreme flat face
Measurement, 7x2 array
25
Measurement, 7x3 array
MST prediction, 7x3 array
Insertion Loss (dB)
20
15
10
-5 2 3
10 10
Frequency (Hz)
Figure 6.26: Measured IL spectra due to 7x2 and 7x3 arrays of industrial latex
cylinders and an MST predicted IL spectrum for a 7x3 latex array. The insert on the
right shows a photograph of the industrial latex array.
Figure 6.26 shows measured IL spectra due to 7x2 and 7x3 arrays of
industrial latex cylinders at normal incidence. The receiver was located at 0.05
m from the centre of the nearest row of the array. The lattice constant for this
set up is 0.069 m and the calculated Bragg band gap central frequency is 2.5
kHz (refer to PWE chapter equation 2.13). Indeed the peak due to the Bragg
band gap is near 2.5 kHz in both measured IL spectra and this band gap is
enhanced by the breathing mode resonance of the latex which is at 900 Hz.
distinct band gap peaks and the other features in the predicted IL spectrum
agrees with those in the measured spectrum. These results indicate that
to the array but one has to consider the larger array area and costs that would
be incurred.
7.1) Introduction
In the previous chapter, it was shown that the insertion loss spectra
elements other than thin elastic shells. So in this chapter, the possibilities for
possibility is the use of the Split Ring Resonator (SRR). The concept of SRR
in which waves cannot propagate through the system [113 and 114]. This
resonant frequencies of the scatterer lie below the first stop band associated
studied intensively [115 and 116]. It has also been shown numerically that
air, with a small opening, or neck, resulting in a coupling between the air in
the neck and the external air in the room [118]. There are also similar air
dimensions of the resonator are small compared with the wavelength of sound
rigid cavity with a protruding neck that connects to the cavity (see figure 7.1).
absorber in that the volume of air in the neck of the Helmholtz resonator
behaves much like a vibration absorber mass and the volume of air acts like a
fluctuation acting over the opening of the neck, resulting in oscillations of the
volume of air in the neck. Damping is the result of radiation losses at the neck
ends and viscous losses due to friction of the oscillating air in the neck. When
air is forced into a cavity, the pressure inside increases. When the external
force pushing the air into the cavity is removed, the higher-pressure air inside
will flow out. The cavity will be left at a pressure slightly lower than the
Air in neck
of bottle
R
Air at
atmospheric Inside air
pressure, . compressed
approximately,
m S Ln , (7.1)
where is the density of the air, S is the cross section area of the neck and
Ln is the effective neck length which includes a correction factor for mass-
loading due to air entrainment near the neck extremities [120]. If this “plug” of
air descends a small distance x into the bottle, it compresses the air in the
container so that the air that previously occupied volume V now has V Sx .
higher value p A p .
pressure in respect of the decreasing volume are fast and temperature rises
process and the resulting equation involves a constant , the ratio of specific
heats, which is about 1.4 for air. The pressure change p produced by a small
p V Sx
. (7.2)
pA V V
F ma ,
d 2x F (7.3)
or 2 ,
dt m
d 2x pS Sp A
x , (7.4)
dt 2
SLn VLn
1
equation for Simple Harmonic Motion (SHM) at a frequency which is
2
The speed c of sound in air is determined by the density, the pressure and
c S
f , (7.6)
2 VLn
resonator at its resonant frequency has been performed by the Finite Element
Method (FEM) using COMSOL® and is shown in figure 7.2. The parameters
are Ln = 0.01 m, outer radius = 0.0275 m, inner radius = 0.0255 m and the
opening = 0.004 m. The pressure map at 410 Hz indicates that the maximum
values of the acoustic pressure occur inside the cavity. This result agrees with
gouging out 0.004 m width slits along the majority of their lengths (Figure 7.3).
However, the difference between the opening in a bottle and a slit on a pipe
account of both the length (i.e. the thickness of the walls) and the aperture of
the neck of the resonator in order to consider the slit pipe as a Helmholtz
used to create band gaps in the range of frequencies below the bandgap
associated with Bragg scattering. The dimensions of the various SRRs used
in our laboratory and outdoor in situ experiments are listed in table 7.1. To
(rigid PVC) and resonating (elastic shell breathing mode) reference elements
are performed as described for the rigid PVC and elastic shells discussed in
x α
d
1.5 m
50 mm
has been used to simulate the scattering of acoustic waves in 2D first with a
single slotted pipe SRR and then for 3 rows of 7 SRRs arranged in a square
lattice in air. For the FEM modelling, the dimensions of the various SRR
scatterers and the experimental setup used are described in table 7.1 and
figure 7.4. The geometry of the SRR can be easily modified from the complete
using the standard 2D draw tool menu. Following the descriptions on chapter
Layers (PML) is used for the domain boundaries. To obtain the level of
noted that the slits are all facing the source at the Ox plane as illustrated in
figure 7.4. The FEM results are then compared with data.
(i.e. cylinder made from PVC without any slit) of the same geometry. For
convenience, the non-resonating references for type SRR1 and SRR2 will be
5
(a) (c)
10
Measurement (NRR1) Measurement (NRR2)
FEM prediction FEM prediction
4
8
Insertion Loss (dB)
2
0
-1 0
-2 2 3
-2 2 3
10 10 10 10
Frequency (Hz) Frequency (Hz)
6 12
5
(b) (d)
10
Measurement (SRR1) Measurement (SRR2)
FEM prediction FEM prediction
4
8
3
6
2
4
1
2
0
-1 0
-2 2 3 -2 2 3
10 10 10 10
Frequency (Hz) Frequency (Hz)
Figure 7.5: Measurements and FEM predictions (of IL spectra for the single cylinders
of type (a) NRR1, (b) SRR1, (c) NRR2 and (d) SRR2 (see table 7.1)
cylinders of types NRR1, SRR1, NRR2 and SRR2 are plotted against FEM
predictions. The IL spectra for single SRR scatterers have low frequency
peaks near 700 Hz for SRR1 and 450 Hz for SRR2 (see figures 7.5(b) and
since they are not present in the IL spectra for identical PVC cylinders (i.e.
with the same diameter and wall thickness) without slits (see figures 7.5(a)
and 7.5(c)). Notably, the measured resonance peak amplitudes are smaller
than predicted for both types of SRR scatterers which could be the result of
damping that is not taken into account in the FEM model. Narrow band gaps
are predicted at higher frequencies for both SRR1 and SRR2. However these
higher order band gaps are predicted to have small amplitudes and are not
Figure 7.6: Pressure maps predicted using COMSOL® for incident cylindrical waves
on (a) a single scatterer SRR1 at 700 Hz and (b) a single scatterer SRR2 at 300 Hz .
the acoustic pressure are inside the cavity. Similar effects are observed
Intuitively, the configuration will work best for the SRR when the slit
opens towards the source since this will trap the majority of the incoming
Insertion Loss (IL) of the single SRR in laboratory while rotating the slit to
different angles (see figure 7.7). For this experiment, SRR1 is used (refer to
table 7.1).
Source
Microphone
3.5
0
3
45
90
2.5
135
Insertion Loss (dB)
180
2
1.5
0.5
-0.5 3
10
twice as much for the slit facing the source (i.e. at 180°) than any other
angles. This will lay the foundation of future sonic crystal designs using SRR.
The effects of the width of the slit opening have been studied using
5
FEM Prediction, 4mm slit width (SRR1)
4 FEM Prediction, 8mm slit width
-1
-2 2 3
10 10
Frequency (Hz)
Figure 7.9: FEM predictions of IL spectra for the single SRR1 with increasing slit
width.
FEM predicted IL spectra for a single SRR1 with slit widths of 4, 8 and
12 mm are shown in figure 7.9. The 750 Hz cavity resonance for the original
SRR1 is shifted to 850 and 950 Hz for slit widths of 8 and 12 mm respectively.
This suggests that increasing the slit width will shift the cavity resonance to
higher frequencies.
square lattice arrays of SRR1 scatterers with three different lattice constants
50
(a)
40
7x3 array, scatterer SRR1, L = 135 mm
FEM Prediction
Measurement
Insertion Loss (dB)
30
20
10
-10 2 3
10 10
Frequency (Hz)
50
(b)
40
7x3 array, scatterer SRR1, L = 80 mm
Measurement
FEM Prediction
Insertion Loss (dB)
30
20
10
-10 2 3
10 10
Frequency (Hz)
50
(c)
40
7x3 array, scatterer SRR1, L = 70 mm
Measurement
FEM Prediction
Insertion Loss (dB)
30
20
10
-10 2 3
10 10
Frequency (Hz)
Figure 7.10: Predicted and measured Insertion Loss spectra due to square arrays of
SRR1 (refer to table 7.1) with lattice constants of (a) 135, (b) 80 and (c) 70 mm. Also
shown in the left hand panels are pressure maps at the first Bragg diffraction
frequency for the corresponding array geometry.
7.6(b)). These additional peaks are not observed for reference arrays of pipes
of the same size but no slits (FEM chapter, figure 4.5). The band gaps due to
the periodicity of the arrays are preserved in these IL spectra for SRR1
suggest that the maximum values of the acoustic pressure are localised in the
cavities of some of the pipes with lengthwise slits. This is consistent with
making larger scale outdoor trials at the OU test barrier site. The outdoor
7.1). The SRR2 scatterers were arranged in square lattice arrays with a lattice
covered in chapter 5 but for the reader’s convenience, the geometry used for
the outdoor experiments is briefly discussed here. The plan view of the array
is shown in figure 7.11(a). The source and receiver are placed 5 m apart from
the array with the source positioned at the asphalt ground and the receiver
test the angle dependency of the SRR barrier performance, the source was
3 x 53
x = 22° x array
= 11°
5m 3 x 54 5m 5m 5m
array = 0° = 0°
1.6 m
= -11° 1.6 m = 6°
= -22° = 8°
Asphalt ground Panel Grass area
(a) (b)
15
10
-5
-10 2 3
10 10
Frequency (Hz)
35
30 (b)
3x54 array, scatterer SRR2,
L = 160 mm
25 Measurement
FEM Prediction
Insertion Loss (dB)
20
15
10
-5
-10 2 3
10 10
Frequency (Hz)
Figure 7.12: Predicted and measured Insertion Loss spectra at outdoor for square
arrays of reference no slit cylinders (NRR2) and SRR scatterers (SRR2) with lattice
constants, L, of 160 mm (refer to table 7.1). (a) 54x3 square lattice array of scatterer
NRR2. (b) 54x3 square array of scatterer SRR2. Corresponding pressure maps at 500
Hz are shown in the left-hand panels.
in figure 7.12(a). Both data and FEM predictions show a Bragg band gap
centred around 1 kHz and higher order modes. In this case it has been found
demands in making predictions for the exact outdoor geometry (array size as
well as distances for source and receiver). The IL peaks below the Bragg
band gap could be due to scattering from ground or other objects near the
occurs at around 450 Hz corresponding to the cavity resonance for the SRR2
scatterer (see figure 7.5(d)). This effect is predicted also. The measured data
also show higher order modes in comparison with arrays of cylinders with no
slits for the same geometry (figure 7.12(a)). We also note the anomalous
noise created by such a barrier with the presence of high downwind (source to
receiver) speeds.
40
35
54x3 array, scatterer SRR2,
L = 160 mm
30 0 deg
-11.2 deg
-21.6 deg
25
+11.2 deg
Insertion Loss (dB)
+21.6 deg
20
15
10
-5 2 3
10 10
Frequency (Hz)
Figure 7.13: Measured Insertion Loss spectra of various lateral angles at outdoor for
square lattice arrays of SRR2 (see table 7.1) with lattice constant of 0.16 m.
Figure 7.13 shows the IL spectra measured for the 54x3 array of SRR2
±22 degree from normal incidence (see figure 7.11(a)). The constancy of the
low frequency peaks near 450 Hz offers experimental support for the
35
30
3x54 array, scatterer SRR2,
25
L = 160 mm
0 degree
6 degree
20 8 degree
Insertion Loss (dB)
15
10
-5
-10 2 3
10 10
Frequency (Hz)
Figure 7.14: Measured Insertion Loss spectra due to outdoor square arrays of split
ring resonators, SRR2 (see table 7.1) with lattice constant of 0.16 m at various
inclinations of the source-receiver axis.
Figure 7.14 shows the IL spectra measured for the 54x3 array of SRR2
degrees from normal incidence (see figure 7.11(b)). Again the low frequency
kHz is due to the Bragg type scattering mechanism and higher order modes
[117 and 122] are shown to support band gaps generated by the multiple
resonances of each scatterer in the array. Since an array of SRR and thin
resonances including Helmholtz and annular cavity resonances and also the
outdoor noise barrier, protecting the shells using concentric outer PVC
by Cui et al suggest the possibility of using hollow tubes with multiple narrow
slits [123], such that all the slits occupy the quadrants facing the source, to
create a broader band gap. In contrast with the slit locations used by Cui et al,
our experiment. To solve the scattering problem for N-slit cylinders we use
boundary conditions dependent on polar angle [124]. In the present model the
finite thickness of the wall of the N-slit cylinder is taken into consideration.
This avoids the use of adjustable parameters introduced for numerical stability
[124 and 125] and provides a more accurate description of the real structure
conditions [126] that describe the slit interface as a moving piston. It should
be noted that use of jump conditions to replace the slits make it possible to
accurately predict the total wave field in the low- and mid-frequency range.
The resulting approximations are of particular interest for studying the low-
consecutive slits in the Oxy plane are denoted by d n , n 1,..., N and they are
infinitely long in the direction of the cylinder main axis Oz . The external length
of the rigid arc is referred to as Dm (see figure 7.15(a)). It is also assumed that
the thickness of the rigid shell is much smaller than its radius so that the
(i) Internal and external arc lengths of the rigid strip have the same
length.
(iii) The angle n subtended by the arc of n-th slit and its width d n are
related by dn n ro , n 1,..., N .
(iv) The angle m subtended by the arc between slit n 1 and slit n
(a) (b)
Figure 7.15: (a) Cross section of composite element consisting of a concentric
arrangement of an outer 4-slit rigid cylinder and an inner elastic cylindrical shell. (b)
Geometry of the slit.
through the middle of the first slit makes an angle of either zero or with the
Oz axis. The acoustic environments outside and inside the N-slit cylinder and
within the slits are assumed to be identical and are described by density
fluid medium previously shown in chapter 4, equation 4.1 is repeated here for
2 p k2 p 0 , (7.7)
2 2 2f
where 2 is the Laplacian, k is the wave number in air
x1 x2
2 2
c
defined as the ratio between angular frequency and sound speed of the
the properties to one of the acoustic media (i.e., it is “o”, “i” and “s” to denote
po 1
ikpo o r 2 , as r , (7.8)
r
where r x 2 y 2 .
cylinder we first solve an auxiliary problem for the slits. The geometry of this
problem is illustrated in figure 7.15(b). Waves propagating inside the slit are
Equation 7.9 is subject to the continuity conditions imposed on the slit faces
h
y which are
2
p ps , (7.10)
p
where q is the normal derivative on the slit faces and index = i or o.
n
Using solution 7.8 and its unknown coefficient, equation can be transformed
pi po hqo , (7.12)
qi k 2 hpo qo . (7.13)
This relates the wave field inside the N-slit cylinder to that of the outer region
7.12 and 7.13) as well as to the Neumann conditions imposed on the rigid
surface of the N-slit cylinder. The former and the latter can be combined into
pi po
f k 2 hpo , (7.14)
r r
po f
po pi 0 , (7.15)
r h
distribution of slits along the N-slit cylinder surface over the interval 0 0,2
so that
1 1 1 1
1 if , 2 N N ,2 N
f 2 2 2 2 , (7.16)
0 otherwise
For the non-zero value of this function equations 7.14 and 7.15 are
transformed to the jump conditions 7.12 and 7.13 whereas the Neumann
zero.
alternative form of f ,
N
2 n 1 n 2 n 1 n
f H H ,
n2 N 2 N 2
(7.17)
1
1 1
H 2 j H 2 j
j 0 2 2
po H 01 kr A H krˆ e
n n
1 in
, (7.18)
n
placed at the centre of the scatterer and An are unknown coefficients. The
pi B J krˆ C Y krˆ e
n
n n n n
inˆ
, (7.19)
derived in terms of Bn or set to zero according to the type of the core so that:
Cn BnCn , n , (7.20)
where
J n' ka1
Cn , rigid core, (7.22)
Yn' ka1
Cn J n' kR 1 k R
2
s
2
n2
, Elastic
J kRY kR
1 k R n n
(7.23)
k s R sRh
2 2 1
'
n n
'
s
2 2 shell core
equation 7.23, the elastic shell is described by wave number k s , density
cs
a1 a2 a a2
R and its half thickness hs 1 . The dilatational wave speed cs
2 2
predictions the size of each elastic shell is always defined by the outer radius
express r , in terms of r̂ , ˆ . Then, substituting (equation 7.18) and (7.19)
into the boundary conditions (equations 7.14 and 7.15) and taking the inner
variables as follows:
k 2 h 1
n m ,n
A 2hH 1'
kr Fn m H 1
kr H 1'
kr H kr
o F j m Fn j j
2
n o n o n o n n
n j
(7.2
H kQe
n
1
n
in
m, n 2hJ kro Fn m J n kro J kro n
'
n
'
n 5)
J n kro F j m Fn j j
2
k h
, m ,
2 j
where m,n is Kronecker delta, vector Q Qcos , sin is the radius vector to
N
for n 0
1
Fn N 1 , (7.27)
2 sin nl e 2in / N
for n 0
n 0 2
Note that factor Cn and the geometrical parameters of the slits only appear in
invariant with respect to scatterer core type and the arrangement and size of
the slits.
To find the numerical solution the infinite system (7.25) has to be truncated at
of slits and their angles (i.e. their length) as well as on the frequency.
dimensions of the various scatterer can be found in table 7.2. The acoustic
structural interaction model is exactly similar to that for the single latex
cylinder described in chapter 6.11 except for the addition of PVC outer
boundary conditions in chapter 6 (table 6.3) to include the PVC cylinder (see
table 7.3). The modified model is then solved in a similar way to that used for
insertion loss spectra for composite scatterers and the results are compared
with data from measurements (see chapter 5 for detailed descriptions). First,
considered. The geometry of the study of 4-slits PVC will use the component
straightforward manner. The FEM predictions are compared with data from
laboratory experiments.
7 (a)
6
Measurement
FEM Prediction
(b)
5
Insertion Loss (dB)
-1
-2 3 4
10 10
Frequency (Hz)
Figure 7.16: (a) Measured and FEM predicted IL spectra due to a single PVC
cylinder with 4 symmetrical slits (component of SRR3, see table 7.2). (b) FEM
(COMSOL®) predictions of pressure maps at 2 kHz for a similar scatterer showing
the interior and exterior environments.
measurement data for the PVC cylinder with 4 symmetrical slits. The same
size of scatterer but with only one slit was found to have a resonance at 700
near 2 kHz.
3
Measurement
FEM Prediction
2
MST Prediction
Insertion Loss (dB)
-1
-2
-3
-4 3
10
Frequency (Hz)
Figure 7.17: Measured and predicted (MST and FEM) IL spectra for a single
composite scatterer, SRR3 (refer to table 7.2)
MST and FEM) shows a strong resonance at around 1 to 2 kHz for composite
frequencies are also observed. For all measurements and predictions, the
predictions disagree with measurements after 2.5 kHz. This could be the
paper.
(b)
(c)
Figure 7.18: FEM (COMSOL®) predicted pressure maps at 1 kHz for SRR3 showing
the exterior and annular cavity environments (a), interior environment of the latex
cavity (b) and static deformation plot on the “Breathing mode” shape of elastic
element in air.
In figure 7.18, FEM predicted pressure maps at 1 kHz for SRR3 (refer
to table 7.2) are shown. According to these pressure maps, higher energies
are observed to be trapped inside annular cavity (see a) and interior of the
4
Measurement
FEM Prediction
3 MST Prediction
Insertion loss (dB)
-1
-2
-3 3
10
Frequency (Hz)
Figure 7.19: Measured and predicted (MST and FEM) IL spectra for a single
composite scatterer, SRR4 (refer to table 7.2)
In figure 7.19, the IL spectrum predicted by FEM agrees with both MST
predictions and our laboratory measurement data for the single composite
scatterer SRR4 (refer to table 7.2). Although the same truncation number (6)
is used in the MST computation as for SRR3, doubling the size of the
scatterer could result in better accuracy for the range of frequencies shown.
5
Measurements
SRR3
4
PVC component
Elastic shell component
Insertion Loss (dB)
-1
-2 3 4
10 10
Frequency (Hz)
Figure 7.20: Measured IL spectra for a single composite scatterer, SRR3 and it’s own
component (refer to table 7.2).
The measured IL spectrum for the latex cylinder in Fig. 7.20 shows the
breathing mode resonance near 1.5 kHz. The resonance for the 4 slits PVC
kHz. This indicates that the coupling between two resonators leads to the shift
5
Measurements
SRR4
4 PVC component
Elastic shell component
Insertion loss (dB)
-1
-2 3
10
Frequency (Hz)
Figure 7.21: Measured (in laboratory) IL spectra for a single composite cylinder
(SRR4) and its components (i.e. single 4 slit PVC and latex, see table 7.2).
mode resonance for a single latex cylinder of the type used in the outdoor
measurements is near 900 Hz. That measured for a single 4 slits PVC
cylinder is near 1.2 kHz. When the cylinders are combined to form a
concentric composite cylinder, the resulting IL spectrum shows two peaks with
reduced magnitude near 750 Hz and 900 Hz and negative IL (dips) before
and after the resonance. These phenomena are similarly observed in the
Following the work of our collaborators [128], this section has been
each element of the array [129] and boundary conditions (equations 7.14 and
7.15) imposed on the surface of all scatterers in the array. As a result the
N
po r , H 01 kr A H n1 krˆm e in m ,
m ˆ
n (7.29)
m 1 n
where N is the number of scatterers in the array, variables rˆm r , , ˆm r , are
the polar coordinates with origin in the centre of scatterer of index m and Anm
are unknown coefficients. The form of the solution inside of the m-th slit
cylinder is the same as in the case of single scatterer problem that is equation
7.19.
(equation 7.29) in terms of polar coordinates rˆm , ˆm of m-th scatterer that
gives
,
m 1,, N
where vector Qm Qm cos m , sin m is the radius vector to the centre of m-th
scatterer, vector Qmp Qmp cos mp , sin mp defines the position of p-th
scatterer with respect to m-th scatterer and Qmp rˆm that is the requirement of
the addition theorem. The expansion (equation 7.30), along with the inner
A
n
p
m, n 2hH n1' kro, p Fnp m H n1 kro, p H n1' kro, p I np
n
H n kro, p F jp m Fnp j I jp
2
k h 1
2 j
Ans H n1v kQps e
N
i n ps
s 1, s p n
m,
2hJ n' kro, p F p m J n kro, p J n' kro, p I np
(7.31)
J kro, p Fjp m Fvp j I jp
k 2h
2 j
Fn m J n kro, p J n' kro, p I np
J n kro, p F jp m Fnp j I jp , m , p 1,, N ,
2
k h
2 j
where factors I mp and Fmp are given by equations 7.26 and 7.27 respectively.
The superscript p relates these factors to p-th scatterer. To find the unknown
Note that for the slit angles N 2M 1 , equation 7.31 is transformed to the
the total wave field (equation 7.29) which can yield the insertion loss. Again
the results for an array of PVC cylinders with 4 symmetrical slits will be shown
before the composite scatterers. This simply expands the FEM model
cylinders due to the high computational resources required and therefore only
predictions of the MST theory described above have been compared with
experimental data.
25
(b)
(a)
20
Measurement
15 FEM Prediction
Insertion Loss (dB)
10
-5
-10 2 3
10 10
Frequency (Hz)
Figure 7.22: (a) Measured and FEM predicted IL spectra for 7x3 PVC array of lattice
constant 0.08 m (PVC with 4 symmetrical slits, component of SRR3, see table 7.2).
(b) FEM (COMSOL®) predicted pressure maps at 2 kHz for a similar scatterer array
showing the interior and exterior environments.
predictions of IL spectra for a 7x3 4-slits PVC array arranged in square lattice
with lattice constant 0.08 m. The cylinder is the PVC (4 slits) component of
SRR3 defined in table 7.2. The band gap effect related to the cavity
resonance (see figure 7.16) and associated with a high insertion loss peak is
predicted amplitudes of the band gap could be due to the fact that the slits on
the PVC cylinders are narrower at their mid points due to the change in their
30
25
Measurement (SRR3)
MST Prediction
20
Insertion loss (dB)
15
10
-5 2 3
10 10
Frequency (Hz)
Figure 7.23: Measured and MST predicted IL spectra for 7x3 array of composite
scatterers SRR3 (see table 7.2) arranged in square lattice of lattice constant 0.08 m.
Figure 7.23 compares the measured and MST predicted IL spectra due
to 7x3 arrays of composite scatterers of type SRR3. This shows the best
agreement between predictions and data. However the predictions use 0.3
mm instead of 0.25 mm for the thickness of the latex shell. This suggests that
the thickness of the latex may not be precisely 0.25 mm as stated by the
prediction. The MST predictions agree fairly well with measurements except
generate several band gaps that result in high insertion loss peaks. One of
these peaks is observed around 1 kHz and is related to the shifted breathing
mode resonance of latex elastic shell. The insertion loss peak associated with
the shifted Helmholtz resonance in the slit cylinder is observed around 2.6
kHz and it follows the first Bragg band gap peak observed and expected
Position Measured
number thickness (mm)
1 0.32
2 0.28
3 0.31
4 0.33
5 0.29
6 0.3
7 0.34
8 0.29
9 0.31
10 0.28
Average thickness: 0.305 mm
Standard deviation: 0.02 mm
Table 7.4: Measured thickness on latex strip with average and standard deviation
values.
table 7.4). According to the results, this yields an average thickness of 0.3 mm
and a standard deviation of 0.02 mm. In view of this, all MST predictions
35
Measurements
30 SRR3
PVC component
Elastic shell component
25
Insertion Loss (dB)
20
15
10
-5 2 3 4
10 10 10
Frequency (Hz)
Figure 7.24: Measured IL spectra due to 7x3 arrays of composite scatterer, SRR3 and
the two components separately (refer to table 7.2).
The results in figure 7.24 for 7x3 arrays of composite cylinders, show
that the measured Bragg band gap near 2 kHz matches with the calculated
breathing mode resonance (See figure 7.20). The enhanced peak near 750
Hz may be associated with the finite height of the array since it also appears
with reduced magnitude in the IL for the PVC pipe array. Similar
measurements for an array of SRR4 composites were not possible due to the
25
Measurements (average of 8 days) with error bar
MST Prediction
20
15
Insertion Loss (dB)
10
-5
-10
-15 2 3
10 10
Frequency (Hz)
Figure 7.25: Average of 8-day measurements (outdoor) and MST IL spectra for 54x3
arrays of composite scatterers SRR4 (see table 7.2) arranged in square lattice of lattice
constant 0.16 m. Error bars are plotted for the measured spectrum.
0.16m) at an outdoor test site on the Open University Milton Keynes campus
have been obtained and are shown together with corresponding MST
predictions in figure 7.25. The error bars represent the spread of data at each
frequency. The setup shown in figure 7.11(a) has been used with the receiver
positioned at θ = 0°. The dimensions of the inner latex tube are the largest
a strong insertion loss peak at 1 kHz which corresponds to both the breathing
mode of the inner elastic shell and the Bragg band gap. The proximity in
frequency of these effects is not ideal for a noise barrier design. For such an
up to about 2 kHz. We also note that there are negative insertion losses at
various frequencies.
20
Measurements
36x4 array
15 54x3 array
Insertion Loss (dB)
10
-5
-10
-15 2 3
10 10
Frequency (Hz)
Figure 7.26: Measured IL spectral of 36x4 and 54x3 arrays of concentric cylinders
(averaged of all 8 measurements).
Figure 7.26 compares measured IL spectra for 36x4 and 54x3 arrays of
average wind speed of 0.7 m/s and 0.9 m/s for the 36x4 and 54x3 concentric
array measurements respectively (see Appendix E(I and ii)). The IL shows
more or less similar performance around the band gap at 1 kHz despite
25
Measurements
54x3 concentric array
20
54x3 rigid array
Insertion Loss (dB)
15
10
-5
-10 2 3
10 10
Frequency (Hz)
concentric cylinders and rigid no slit cylinders. Although for the concentric
near 1.2 – 2 kHz, there is a worse performance between 500 Hz and 800 Hz.
10
-10
-20 2 3
10 10
Frequency (Hz)
crystal except for the distinctive improvement at 1.2 kHz. The IL dip for the
fence at this frequency could be attributed by factor such as holes along the
The elastic shell outer diameter, elastic shell thickness and the slit
The influence of the outer diameter of the inner elastic shell has been
for this study are 54x3 arrays of SRR4 scatterers arranged in a square lattice
with a lattice constant of 0.16 m. The shell outer diameter has been increased
front and back faces of the array with the source to receiver line of sight
centred at the gap between the 2 central scatterers (see figure 7.11).
35
30
MST Prediction (latex OD)
25 43 mm
53 mm
20 63 mm
Insertion Loss (dB)
83 mm
15
10
-5
-10 2 3
10 10
Frequency (Hz)
Figure 7.29: IL spectra for the MST predicted effect on changing latex outer diameter
for 54x3 square lattice array of SRR4 scatterer with lattice constant of 0.16 m.
increasing separation of the modified latex resonance from the 1 kHz Bragg
band gap. Also there is a relationship between the predicted amplitude of the
modified latex resonance and that of the IL peak associated with the Bragg
in the predicted Bragg band gap amplitude and vice-versa. However again
the inner elastic shell. Following the predictions in figure 7.29 which show that
frequency band gap. The wall thickness for such scatterer has been varied
25 MST Predictions
(latex thickness)
0.5 mm
20
1 mm
1.5 mm
2 mm
Insertion Loss (dB)
15
10
-5
-10 2 3
10 10
Frequency (Hz)
Figure 7.30: IL spectra for the MST predicted effect on changing latex wall thickness
for 54x3 square lattice array of SRR4 scatterer with lattice constant of 0.16 m.
thickness) as the latex wall thickness is increased. The Bragg band gap
The effects of varying the slit widths in the outer 4-slits PVC cylinder
have been predicted. The assumed slit width has been varied from 3 to 12
mm.
40
MST Prediction (slit width)
35
12 mm
30 6 mm
3 mm
Insertion loss (dB)
25
20
15
10
-5 2 3
10 10
Frequency (Hz)
Figure 7.31: MST predicted IL spectra for 54x3 square lattice array with lattice
constant of 0.16 m as the slit widths in the outer PVC cylinders are varied from 3 to
12 mm.
The resulting IL spectra predictions are shown in Figure 7.31 for 54x3
4-slits PVC arrays with lattice constant 0.16 m and different slit widths. The IL
peak associated with the Bragg band gap slowly broadens near 1 kHz as the
slits width decreases. This may explain the modification of the 1st Bragg band
gap in data for our current outdoor setup (see figure 7.25). The predicted
effects of varying the parameters of the composite elements can lead to better
array design for outdoor noise reduction. This will be pursued in chapter 9.
8.1) Introduction
absorption and extra attenuation (i.e. barrier and ground effect) [131].
at a rate that depends on the type of source (point source or line source).
propagation. The wind profile can also have pronounced effect. For example a
down wind condition (wind blowing towards the receiver from the source) will
the barrier. At distances less than 100 m, meteorological effects play a minor
role and the acoustical properties of the ground – typically given by the
different physical properties are known to produce acoustic band gaps, also
If the scatterers are solid and the host medium is air then these are called
bands within which the amplitudes of transmitted waves are relatively high.
Pass bands and band gaps are a consequence of multiple scattering of waves
barriers, recent research into sonic crystal barriers for highway noise has
ground plane. Clearly this will be unrealistic if SCs are to be used as noise
barriers since a ground will always be present. On the other hand, the
been investigated [132 and 133] and the sound diffraction effect of a
(see chapter 9.2), for a single noise barrier has been studied using a 2
chapter is to study the total effect of SC and a reflecting ground. First there is
on the superposition of the solution for a single scatterer. Although the most
axes parallel to the ground. The semi-analytical method is used to analyse the
arrays in the presence of rigid ground is compared with that of the same
arrays in the free field. Insertion Loss (IL) data from laboratory experiments
are compared with predictions. Laboratory studies have been performed for
described.
The prime objective of this chapter is to study the total effect of sonic
means they are perfectly reflecting, or ‘acoustically soft’, which implies that
they are perfectly absorbing. According to ISO 9613-2, any ground surface of
above a surface can be done through measuring the sound pressure (easiest
parameter than can be measured for sound field) at one or more locations. A
determined from the spectrum of the ratio of the direct source to receiver total
sound level without pFree _ field and with pTotal_ field a ground surface (equation
8.1), keeping the distance between source to receiver and the source strength
constant. This is called the Excess Attenuation (EA) due to the ground surface
properties. If expressed in decibels, a positive value will imply that the sound
pressure is lower than the sound pressure in free field. We note that this
[136].
pTotal _ field
EA 20 log . (8.1)
p Free _ field
single parameter model developed by Delany and Bazley [137] in 1970. They
material where the flow resistivity represents the ratio of the applied pressure
gradient to the induced volume flow rate per unit thickness of material due to
the porosity. This single parameter model has been widely used in deducing
kc 1 0.0978 X 0.700 i 0.189 X 0.595
.
c0 (8.2)
0
where X , 0 , is the density of air, , is an effective flow resistivity
dimension to the thickness of the viscous boundary layer and this ratio is
the model has been used with tolerable success to characterize outdoor
located above the surface and a receiver point also located above the surface
as shown in figure 8.1(a). This allows the total field pressure pTotal_ field to be
measured. The free field spectrum level is obtained either by removing the
ground plane or by raising the source and receiver sufficiently far above the
surface that ground effect may be ignored. The latter method proved to be
more prone to error due to the required precision for the source-receiver
distance (see equation 8.1). Thus in our laboratory measurement, a flat table
(figure 5.16) is used as the ground plane support so that to establish the free
field spectrum, both the ground plane and the supporting table were removed
The above measuring technique serves well for deducing the excess
typical measurement setup for the transfer function method is shown in figure
symmetry is assumed the transfer function between the two locations will
characterize the sound field which depends on the surface impedance and
geometry. The user has to measure the Level Difference (LD) using only
magnitude between the two microphone locations and swop around both
The level difference spectrum is then calculated as the difference between the
excess attenuation spectra at the two receivers. Since this method eliminates
hs Rreflected hr
(b) Rdirect
Source Receivers
Upper mic.
Rreflected
hru
hs Lower mic.
hrl
θ
Figure 8.1: Measurement configuration to determine the surface impedance. (a) Single
microphone method (laboratory) and (b) Transfer function method (outdoor in situ).
Two Medium Density Fibre (MDF) boards and a glass plate have been
of polyurethane form layer has been used as ‘acoustically soft’ ground. The
surfaces and their dimensions are listed in table 8.1. Various ranges (source
to receiver distance) and heights of source and receiver have been used
according to the different types of surface materials (see table 8.2). During
measurements, care was taken to ensure that the ground planes were
horizontal by testing with a tubular sprit level gauge (see figure 8.2).
Ground (Material) Range (m) Source height (m) Receiver height (m)
MDF 1 1 0.10 0.10
1 0.20 0.20
1 0.30 0.30
MDF
Figure 8.2: Levelling test for MDF board using tubular spirit level gauge.
10
0
Excess Attenuation (dB)
-5
-10
-15
Prediction, = 10 GPa s m -2
-25
(a)
-30 3 4
10 10
Frequency (Hz)
10
0
Excess Attenuation (dB)
-5
-10
-15
Prediction, = 10 GPa s m -2
-25
-30
(b)
3 4
10 10
Frequency (Hz)
10
0
Excess Attenuation (dB)
-5
-10
-15
Prediction, = 10 GPa s m -2
-25
(c)
-30 3 4
10 10
Frequency (Hz)
Figure 8.3 – Measured and predicted EA spectra for ground (MDF 1) at source and
receiver heights at (a) 0.10 m (b) 0.20 m and (c) 0.30 m (see Tables 8.1 and 8.2)
Figure 8.3 shows the excess attenuation spectra obtained over the
surface of MDF 1 with various source and receiver heights (see tables 8.1 and
effect are at 8.2 kHz (figure 8.3(a)). Higher order dips are observed for source
and receiver heights of 0.20 and 0.30 m (figures 8.3(b) and 8.3(c)) when the
source and receiver are elevated. The excess attenuation predicted using the
Delany and Bazley one parameter model with an effective flow resistivity of 10
GPa s m-2 gives reasonable agreement with the data in the frequency range
the angle of incident sound [138]. The formula used in the EA prediction is
10
0
Excess attenuation (dB)
-5
-10
-15
Prediction, = 10 GPa s m -2
-25
(a)
-30 2 3 4
10 10 10
Frequency (Hz)
10
0
Excess Attenuation (dB)
-5
-10
-15
Prediction, = 10 GPa s m -2
-25
(b)
-30 2 3 4
10 10 10
Frequency (Hz)
10
0
Excess Attenuation (dB)
-5
-10
-15
Prediction, = 10 GPa s m -2
-25
(c)
-30 2 3 4
10 10 10
Frequency (Hz)
Figure 8.4 – Measured and predicted EA spectra over MDF 2 ground with source
height fixed at 0.23 m and receiver heights of (a) 0.11 m (b) 0.23 m and (c) 0.35 m
(see tables 8.1 and 8.2)
Figure 8.4 shows the excess attenuation spectra obtained over the
MDF 2 board (see table 8.1) using a fixed source height of 0.23 m and three
receiver heights (see table 8.2). Despite the difference in dimensions of MDF
Chapter 8 Page 245
2 from MDF 1, the ground effect minima agree reasonably well with
predictions using the one parameter model for ground impedance with an
could be due to the effect of diffraction by the edges of the board as the
15
10
0
Excess attenuation (dB)
-5
-10
-15
-20
-35 2
(a)
3 4
10 10 10
Frequency (Hz)
15
10
0
Excess attenuation (dB)
-5
-10
-15
-20
Measurement, R = 1 m, h = 15 cm, h = 15 cm
-25 direct s r
-2
Prediction, = 30 GPa s m
-30
-35 2
(b)
3 4
10 10 10
Frequency (Hz)
15
10
0
Excess attenuation (dB)
-5
-10
-15
-20
-35 2
(c)
3 4
10 10 10
Frequency (Hz)
Figure 8.5 – Measured and predicted EA spectra over ground (glass) with source and
receiver heights at (a) 0.10 m (b) 0.15 m and (c) 0.20 m respectively (see tables 8.1
and 8.2).
Figure 8.5 shows the excess attenuation spectra obtained over the
glass plate (see table 8.1). The source and receiver heights were 0.10, 0.15
effect minima agree with predictions using the one parameter model for
This higher fitted flow resistivity value suggests that the glass plate is
10
-5
-10
Measurement, R = 0.8 m, h = 23 cm, h = 11 cm
direct s r
-2
Prediction, = 20 kPa s m
-15
(a)
-20 2 3
10 10
Frequency (Hz)
10
5
Excess Attenuation (dB)
-5
-10
Measurement, Rdirect = 0.8 m, hs = 23 cm, hr = 23 cm
Prediction, = 20 kPa s m -2
-15
(b)
-20 2 3
10 10
Frequency (Hz)
10
5
Excess Attenuation (dB)
-5
-10
Measurement, Rdirect = 0.8 m, hs = 23 cm, hr = 35 cm
Prediction, = 20 kPa s m -2
-15
(c)
-20 2 3
10 10
Frequency (Hz)
Figure 8.6 – Measured and predicted EA spectra over ground (polyurethane foam)
with source height fixed at 0.23 m and receiver heights of (a) 0.11 m (b) 0.23 m and
(c) 0.35 m (see tables 8.1 and 8.2)
Figure 8.6 shows the excess attenuation spectra obtained over the
polyurethane foam layer (see table 8.1). The source height is fixed at 0.23 m
respectively (see table 8.2). The predicted ground effect minima using the one
resistivity of 20 kPa s m-2 are in reasonable agreement with the data in the
The transfer function method described in section 8.3.2 was tested first
different from that used for excess attenuation measurement i.e. the one
channel MLSSA system. A similar omnidirectional source was used for both
provided in chapter 5.
25
10
Level difference (dB)
-5
-10
-15
-20 2 3 4
10 10 10
Frequency (Hz)
Figure 8.7 – Measured and predicted LD spectra over MDF 2 with source height fixed
at 0.07 m and receiver heights of 0.07 m and 0.14 m (see tables 8.1 and 8.2)
obtained over the MDF 2 board (see table 8.1) using a source height of 0.07
m and receiver heights of 0.07 and 0.14 m (see table 8.2). The measured
ground effects agree reasonably well with predictions using the one parameter
model for ground impedance with an effective flow resistivity of 10 GPa s m-2
as was the case for the data obtained with the one microphone method (see
figure 8.4). This confirms the reliability of the transfer function method for
Given the consistency of the flow resistivity deduced for MDF 2 board
by both methods in the laboratory (section 8.3.1 and 8.3.2); the transfer
30
10
-10
-40 2 3 4
10 10 10
Frequency (Hz)
Figure 8.8 – Measured and predicted LD spectra over asphalt surface with
measurement geometry 2 (table 8.3).
Figure 8.8 shows the level difference spectra obtained over the asphalt
surface at the vicinity of the outdoor barrier test site. The measurement
follows the geometry (B) suggested by ANSI (see table 8.3) with a source
height of 0.2 m and receiver heights of 0.05 and 0.2 m. The measured ground
effects agree reasonably well with predictions using the one parameter model
for ground impedance with an effective flow resistivity of 50 GPa s m-2, except
for the maxima and minima where the amplitudes are over predicted. This
suggests that asphalt ground is acoustically harder than glass plate and MDF
board.
The ground impedance for grass covered ground has been deduced
10
Level difference (dB)
-10 2 3 4
10 10 10
Frequency (Hz)
Figure 8.9 – Measured and predicted LD spectra over grass surface with measurement
geometry B (table 8.3).
Figure 8.9 shows the level difference spectra obtained over the grass
surface at the vicinity of the outdoor barrier test site. Predictions (using the
Delany and Bazley model give a better fit with the grass ground than obtained
with asphalt. The effective flow resistivity is predicted to be 200 kPa s m-2
which is in good agreement with typical flow resistivity values (150 – 300 kPa
scatterers suspended with their axes parallel to a rigid ground. In this section,
with that of the crystals placed in the free field. Insertion Loss (IL) spectra
space characterised by the sound speed in air c = 344 m/s and density = 1.2
kg/m3. The source is assumed to be a point source. Figure 8.10 illustrates the
geometry of the array and its image. The position of each scatterer C m ,
lattice configuration.
O O
x x
(a) (b)
Figure 8.10: Square lattice array above a perfectly reflecting plane. (a) Set of vectors
used in equation (8.6). (b) Set of vectors employed in equation (8.14).
fluid medium previously shown in chapter 4, equation 4.1 is repeated here for
2 2 2f
where 2 is the Laplacian, k is the wave number in air
x1 x 2
2 2
c
defined as the ratio between angular frequency and sound speed in air, c .
rigid ground) and on the surface of the scatterers (this condition has to be
chapter 6) that is
p
0. (8.5)
n
images [141] the general solution of the formulated problem can be written as
[142]
pr p0 r ps r . (8.6)
whereby contributions from the point source and its image are collected in p0
p0,r r H 01 kr0' , (8.9)
p0 r ps ,d r ps , r r , (8.10)
M
ps ,r r Anm Z nm H n1 krm' e in m .
'
(8.12)
m 1 n
The vector r0 r0 cos 0 , sin 0 connects the point source and the receiver
point (i.e. point P in figure 8.11(a)). The vector rm rm cos m , sin m connects
the centre of scatterer Cm and the receiver. The prime ' is used to indicate
the geometrical parameters for the image source and the image scatterers
coefficients. The solution for the unbounded acoustic space can be retrieved
from equation 8.6 by putting to zero in equations (8.7) and (8.10) all terms
related to the constructed images that are p0, r and ps ,r respectively. One can
also deduce from equation (8.6) and vector definitions in figure 8.10(a) that for
the source and receiver both on the ground, the acoustic pressure in a half-
The factors Z nm describe the type of conditions imposed on the surface of the
r J n kam
Z nm , (8.13)
r H n1 kam
polar coordinate r.
system of equations can be derived to find the unknown coefficients Anm . This
system is given by
A Z
H q1n kRmp e
M
i q n mp
Anm q
p
q
p
q p 1, p m
M
A
p 1
q
p
Z qp H q1n kRmp
'
e i q n mp
'
iq
(8.14)
H n1 kR0 p e in 0 m H n1 kR0' m e in 0 m , n Z , m 1M ,
'
with respect to point source and vector R mp Rmp cos mp , sin mp defines the
position of scatterer C p with respect to scatterer Cm . Again the system of
equations 8.14 can be transformed to that for the case of unbounded acoustic
image source and scatterers (see equation 3.36). To solve numerically for the
significant figures [49]. It is also noted that for the considered configurations
and frequency range the computation time required to solve system (equation
Also note that throughout this chapter the insertion loss (equation 3.37) is
calculated. In free field the cylinders locations in the lower half of the 7x3 array
are defined by the coordinates identical to those of the image cylinders in the
half-space problem.
insertion loss spectrum due to an array of elastic shells with their axes parallel
to a rigid ground. The identical elastic shells are characterised by their density
value of these parameters, the first elastic shell resonance associated with the
breathing mode can be observed below the first Bragg band gap associated
with the lattice constant of the array in the unbounded acoustic space (see
The asymptotic theory of thin elastic shells [107] has been used to derive the
factors
r J n kS
Z nm , (8.15)
r H n1 kS iU n
where
n 2 k32 S 2
Un
Sh1 n 2 k32 S 2 r J n kS
. (8.16)
c c 1 2
is the relative impedance, and 3 . If the relative
s c2 c2 E
spaced parallel rigid cylinders and elastic shells without and with ground
planes have been carried out in an anechoic chamber. The rigid cylinders
elastic shells were made from 0.25 mm thick sheets of latex by overlapping
the edges and gluing them together. The sound source was a Bruel & Kjaer
the receiver. Figure 8.12 (a), (b) and (c) shows example measurement
MDF boards at the top and base of each array. To maintain their shape and
hs hr1
Ground
d1 d2 d3
d
Figure 8.11: (a) Schematic of the experimental configuration showing the source
location, the sonic crystal array and three receiver locations. (b) Photograph of
experimental arrangement with rigid cylinder above rigid (MDF 2 board) ground. (c)
Photograph of experimental arrangement with latex shell array near to rigid (MDF 1
board) ground (see table 8.1).
An MDF board large enough to avoid the diffraction at the edges was
above a horizontal MDF board (figure 8.11 (b)). However since the latex
cylinders had to be arranged vertically to preserve their shape, the MDF board
was also supported vertically (figure 8.11 (c)). In both cases the cylinder axes
were parallel to the board. For simplicity all distances to the MDF board in
both setups are referred to as heights. As shown in figure 8.11(a), the loud
speaker point source was positioned d 1 = 0.75 m from the array of rigid
cylinders and d1 = 0.35 m from the array of latex shells at the height of the
horizontal mid-plane of the array ( hs = 0.23 m above the ground). The height
and it was placed in a vertical plane d 3 = 0.25 m from the back of the array.
The receiver heights were chosen to be below, at and above, the horizontal
mid-plane of the array. In all cases, the distance between the microphone and
sound levels recorded in the x direction (0°) at the same point behind the
50
Measurement, 7x3 array over rigid ground, Rdirect = 1.2 m, hs = 23 cm, hr1 = 11 cm
40 Prediction, sigma = 10 GPa s m-2
Excess Attenuation/Insertion Loss (dB)
30
20
10
-10
-20
-30
(a)
-40 2 3 4
10 10 10
Frequency (Hz)
50
40 Measurement, 7x3 array over rigid ground, Rdirect = 1.2 m, hs = 23 cm, hr2 = 23 cm
Prediction, = 10 GPa s m-2
Excess Attenuation//Insertion Loss (dB)
30
20
10
-10
-20
-30
-40 2
(b)
3 4
10 10 10
Frequency (Hz)
50
Measurement, 7x3 array over rigid ground, Rdirect = 1.2 m, hs = 23 cm, hr3 = 35 cm
40
Prediction, = 10 GPa s m-2
30
20
Insertion Loss (dB)
10
-10
-20
-30
-40 2
(c)
3 4
10 10 10
Frequency (Hz)
Figure 8.12 – Measured (solid line) and predicted (broken line) IL spectra due to a
square lattice array of 7x3 rigid cylinders of diameter 0.055 m over acoustically-rigid
ground with source coordinates (0, 0.23) m and receiver coordinates (a) (1.20, 0.11)
m, (b) (1.20, 0.23) m and (c) (1.20, 0.35) m.
for a 7x3 rigid cylinder array over rigid ground for three receiver heights. The
experiments described in section 8.5.3. Up to 1.5 kHz the predictions and data
are in close agreement. Above 1.5 kHz there are some discrepancies which
may be due to unwanted reflections and departures from the assumed ideal
geometry. Both data and predictions in figures 8.12(a) and 8.12(b) show IL
maxima near 2.5 kHz which are associated with the Bragg band gaps
expected in the unbounded domain. Both data and predictions for the
elevated receiver height (0.23 m and 0.35 m) show the adverse influence of
the rigid ground effect on the IL spectra near 2 kHz and 1.25 kHz. This point
will be further elaborated in figure 8.13 with the IL spectra plotted against the
excess attenuation spectra. Indeed in figure 8.12(b), for the receiver at 0.235
m height, it is clear that the destructive interference in the ground effect near
30
(a)
R = 1.2 m, h = 23 cm, h = 11 cm
direct s r1
25
Measurement, 7x3 array over rigid ground
Measurement, 7x3 array in the free field
20
Measurement, EA (inverted)
10
-5
-10
-15 2 3 4
10 10 10
Frequency (Hz)
35
30
(b)
R = 1.2 m, h = 23 cm, h = 23 cm
direct s r2
25 Measurement, 7x3 array over rigid ground
Measurement, 7x3 array in the free field
20 Measurement, EA (inverted)
Insertion Loss (dB)
15
10
-5
-10
-15 2 3 4
10 10 10
Frequency (Hz)
35
30
(c)
R = 1.2 m, h = 23 cm, h = 35 cm
direct s r3
25
Measurement, 7x3 array over rigid ground
Measurement, 7x3 array in the free field
20
Measurement, EA (inverted)
Insertion Loss (dB)
15
10
-5
-10
-15 2 3 4
10 10 10
Frequency (Hz)
Figure 8.13 – Measured insertion loss spectra due to a 7x3 rigid cylinders array over
acoustically rigid ground (solid line) compares with identical array measured in free
field (broken line) as well as measured inverted excess attenuation (dash dot line) The
source is at coordinates (0.0.235) m and the receiver coordinates are (a) (1.20, 0.11)
m, and (c) (1.20, 0.35) m.
pipe array with and without a rigid (MDF 2) ground at source height of 0.23 m
and three receiver heights of 0.11, 0.23 and 0.35 m as shown in (a), (b) and
(MDF 2) ground for a similar source-receiver geometry to that used for the IL
measurements is superimposed into the IL plot but inverted (see figure 8.4).
The Bragg band gap in this particular arrangement (lattice constant of 0.69 m)
is near 2.5 kHz. The array band gaps are impaired, resulting in negative
insertion loss, where they coincide with the minima in the ground.
50
Measurement, 7x3 array over rigid ground, R direct = 0.8 m, hs = 23 cm, hr1 = 11 cm
40 Prediction, = 10 GPa s m -2
30
Insertion Loss (dB)
20
10
-10
-20
(a)
-30 2 3 4
10 10 10
Frequency (Hz)
50
Measurement, 7x3 array over rigid ground, R = 0.8 m, h = 23 cm, h = 23 cm
direct s r2
40 Prediction, = 10 GPa s m -2
30
Insertion Loss (dB)
20
10
-10
-20
(b)
-30 2 3 4
10 10 10
Frequency (Hz)
50
Measurement, 7x3 array over rigid ground, R direct = 0.8 m, hs = 23 cm, hr2 = 23 cm
40 Prediction, = 10 GPa s m -2
30
Insertion Loss (dB)
20
10
-10
-20
-30 2
(c)
3 4
10 10 10
Frequency (Hz)
Figure 8.14 – Measured (solid line) and predicted (broken line) IL spectra due to a
square lattice array of 7x3 latex shell cylinders of diameter 0.055 m over acoustically-
rigid ground with source coordinates (0, 0.23) m and receiver coordinates (a) (0.8,
0.11) m, and (c) (0.8, 0.35) m.
.
over the rigid ground shown in figure 8.14 are similar to that described in
section 8.5.2. In particular, figure 8.14(c) shows that when the first ground
effect dip is in the vicinity of the breathing mode resonance of the shell (near
900 Hz) the corresponding positive IL peak is no longer present. Again, this
point will be further elaborated in figure 8.15 with the excess attenuation
-10
-20
(a)
-30 2 3 4
10 10 10
Frequency (Hz)
40
R = 0.8 m, h = 23 cm, h = 23 cm
direct s r2
30
Measurement, 7x3 array over rigid ground
Measurement, 7x3 array in the free field
Measurement, EA (inverted)
20
Insertion Loss (dB)
10
-10
-20
(b)
-30 2 3 4
10 10 10
Frequency (Hz)
40
Rdirect = 0.8 m, hs = 23 cm, hr3 = 35 cm
Measurement, 7x3 array over rigid ground
30 Measurement, 7x3 array in the free field
Measurement, EA (inverted)
20
Insertion Loss (dB)
10
-10
-20
(c)
-30 2 3 4
10 10 10
Frequency (Hz)
Figure 8.15 – Measured insertion loss spectra due to a 7x3 latex shell cylinders array
over acoustically rigid ground (solid line) compares with identical array measured in
free field (broken line) as well as measured inverted excess attenuation (dash dot line)
The source is at coordinates (0.0.235) m and the receiver coordinates are (a) (1.20,
0.11) m, and (c) (1.20, 0.35) m.
Figure 8.15 compares the measured IL spectra (see equation 1.9) for a
horizontal elastic shell cylinders array with and without a rigid (MDF 2) ground
at source height of 0.23 m and three receiver heights of 0.11, 0.23 and 0.35 m
as shown in (a), (b) and (c) respectively. The measured excess attenuation
used for the IL measurements is superimposed into the IL plot but inverted.
The Bragg band gap in this particular arrangement (lattice constant of 0.069
m) is near 2.5 kHz. The array band gaps are impaired, resulting in negative
insertion loss, where they coincide with the minima in the ground.
40
35
R = 1,2 m, h = 23 cm, h = 11 cm
direct s r1
30
Measurement, 7x3 array over soft ground
Measurement, 7x3 array in the free field
25
Measurement, EA (inverted)
Insertion Loss (dB)
20
15
10
-5
-10 2
(a)
3 4
10 10 10
Frequency (Hz)
40
35
R = 1,2 m, h = 23 cm, h = 23 cm
direct s r2
30 Measurement, 7x3 array over soft ground
Measurement, 7x3 array in the free field
25 Measurement, EA (inverted)
Insertion Loss (dB)
20
15
10
-5
-10 2
(b) 3 4
10 10 10
Frequency (Hz)
40
35
R = 1,2 m, h = 23 cm, h = 35 cm
direct s r3
Excess Attenuation/Insertion Loss (dB)
20
15
10
-5
-10 2
(c)
3 4
10 10 10
Frequency (Hz)
Figure 8.16: Measured insertion loss spectra due to a 7x3 rigid cylinders array over
finite impedance (Polyurethane foam) ground (solid line) comparing with identical
array measured in free field (broken line) as well as measured inverted excess
attenuation spectra of the ground (dash dot line) taken at same source-receiver
distances. The source is at coordinates (0.0.235) m and the receiver coordinates are (a)
(1.20, 0.11) m, and (c) (1.20, 0.35) m.
rigid cylinder arrays over finite impedance ground (Polyurethane foam) for
three receiver heights using the source location described in section 8.5.3.
The measure IL spectra due to identical arrays without the presence of the
ground and excess attenuation (inverted) due to the Polyurethane foam layer
Compared to the results for the rigid ground plotted in figure 8.13, the minima
associated with the finite impedance ground effect seem to enhance the IL
peaks and also we note that the ground effects are also shifted towards lower
barriers.
50
(a)
40
Rdirect = 0.8 m, hs = 23 cm, hr1 = 11 cm
Measurement, 7x3 array over soft ground
30 Measurement, 7x3 array in the free field
Insertion Loss (dB)
Measurement, EA (inverted)
20
10
-10 2 3 4
10 10 10
Frequency (Hz)
50
(b)
40
Rdirect = 0.8 m, hs = 23 cm, hr2 = 23 cm
Measurement, 7x3 array over soft ground
30
Measurement, 7x3 array in the free field
Insertion Loss (dB)
Measurement, EA (inverted)
20
10
-10
-20 2 3 4
10 10 10
Frequency (Hz)
50
(c)
40
Rdirect = 0.8 m, hs = 23 cm, hr3 = 35 cm
Measurement, 7x3 array over soft ground
30 Measurement, 7x3 array in the free field
Measurement, EA (inverted)
Insertion Loss (dB)
20
10
-10
-20 2 3 4
10 10 10
Frequency (Hz)
Figure 8.17 – Measured insertion loss spectra due to a 7x3 latex shell cylinders array
over finite impedance (Polyurethane foam) ground (solid line) compared with that due
to an identical array measured in free field (broken line) as well as measured inverted
excess attenuation spectra due to the ground alone (dash dot line) The source is at
coordinates (0.0.235) m and the receiver coordinates are (a) (1.20, 0.11) m, (b) (1.20,
0.23) m and (c) (1.20, 0.35) m.
latex cylinder arrays over finite impedance ground (Polyurethane foam) for
three receiver heights using the source location described in section (a), (b)
and (c) respectively. The measure IL spectra for identical arrays without the
superimposed on the plots. According to figure 8.17, the ground effects are
contrary to what has been observed for the PVC cylinders array over the finite
impedance ground (see figure 8.16). Again, the band gaps are impaired,
resulting in negative insertion loss, where they coincide with the first minima in
the ground. Subsequent dips are too small in amplitude and have negligible
Figure 8.18: 7x3 PVC cylinders in vertical array orientation above a MDF ground
(MDF 1).
40
R = 1.3 m, h = 30 cm, h = 30 cm
direct s r
Measurement, 7x3 vertical array over rigid ground
30 Measurement, 7x3 array in the free field
Measurement, EA (inverted)
10
-10
-20 2 3 4
10 10 10
Frequency (Hz)
Figure 8.19: Measured insertion loss spectra of vertical array with source and receiver
at 0.3m height with (solid line) and without (broken line) the MDF 1 ground plane.
Measured inverted excess attenuation (inverted) spectra of the ground (dash dot line).
Figure 8.19 compares IL spectra due a vertical axis array without and
with the MDF 1 surface and source and receiver at 0.3 m height. Without the
ground surface, there is a band gap centred at about 2.5 kHz, due to multiple
scattering by the (acoustically hard) cylinders. This corresponds with the first
the reflecting surface, the first and higher order band gaps are impaired when
they coincide with minima in the EA spectra. On the other hand, for this
geometry, the first ground attenuation maximum due to the MDF 1 board is
60
R = 1 m, h = 10 cm, h = 10 cm
direct s r
Measurement, 7x3 vertical array over glass ground
40 Measurement, 7x3 array in the free field
Measurement, EA (inverted)
Insertion Loss (dB)
20
-20
-40
(a)
2 3 4
10 10 10
Frequency (Hz)
60
R = 1 m, h = 15 cm, h = 15 cm
direct s r
Measurement, 7x3 vertical array over glass ground
40
Measurement, 7x3 array in the free field
Measurement, EA (inverted)
Insertion loss (dB)
20
-20
-40
(b)
2 3 4
10 10 10
Frequency (Hz)
60
R = 1 m, h = 20 cm, h = 20 cm
direct s r
Measurement, 7x3 vertical array over glass ground
40 Measurement, 7x3 array in the free field
Measurement, EA (inverted)
Insertion Loss (dB)
20
-20
-40
(c)
2 3 4
10 10 10
Frequency (Hz)
Figure 8.20 – Measured insertion loss spectra of sonic crystal in vertical orientation
with (solid line) and without (broken line) a ground plane (glass plate). The source
and receiver heights are at (a) 0.1, (b) 0.15 and (c) 0.2 m respectively. Measured
excess attenuation (inverted) spectra of the ground is also shown (dash dot line).
cylinder arrays over acoustically solid ground (glass plate) for three receiver
heights using the source and receiver location described in table 8.2. The
destruction of the band gap at around 2.3 kHz due to its coincidence with a
ground effect minimum for this geometry can be observed again (see figure
8.20 (c)).
derivations of some formulations for both fluid and elastic waves. Only linear
wave phenomena were considered. The issue of road traffic noise impact was
barrier. A review of past work on the use of sonic crystals as noise barrier was
that a periodic square array of rigid cylinders in air can give rise to complete
and partial stopbands in which sound cannot propagate through at all or only
at certain angles respectively. It was seen also that by tuning the lattice
predict the centre frequency of the Bragg band gap for a specific square
lattice array. This formula must be modified if the configuration of the array is
changed to another Bravais lattice (i.e. hexagonal, see equation 4.8). The
the laboratory results (finite size and presented as Insertion Loss (IL)) for a
Bragg band gap. An investigation of the effect of varying the size of the
scatterer (which inherently varies the filling fraction) was performed using the
chapter, enabled us to infer that it is, in fact, the density contrast that strongly
influences the occurrence of the band gaps. The ability to determine the band
structure through PWE could be a cost effective way for making the proper
transmission loss for a two dimensional (2D) periodic array of identical rigid
scatterer. The basis of this theory uses the multipole methods to combine
describe the scattering of waves by circular scatterers were involved. The first
the first kind to express an outgoing wave radiating from each cylinder.
infinite arrays, MST calculates the total wave field (pressure) for a finite size
array at any point exterior to the cylinders. This method thereby offers a
means of direct comparison with laboratory results. The incident wave can be
modelled using a plane or cylindrical waves where the latter can take account
of the near field effects (i.e. when source is placed close to the array). The
MST method has proved to be powerful, and self-consistent with the ability to
solve in principle all multiple-body (i.e. the configuration of the scatterers can
individual scatterer are already known. This method was used to calculate the
insertion loss for a 7x3 square lattice arrays consisting of rigid circular
and the results were compared with laboratory data. The MST predicted
results (both plane and cylindrical incident waves) were found to agree well
with the laboratory data and these results complimented the hypothesis that
the Bragg band gap could be tuned to the desired frequency by varying the
lattice constant of the rigid scatterers. The data indicated minor peaks at
constant equal to 0.069 and 0.08 m that could be the effect of finite width or
height of the sonic crystal array. This extra low frequency band gap could be
exploited in SC noise barrier design. The accuracy of the MST has been
investigated and for a 7x3 square lattice array of rigid circular scatterers with
and M 3 respectively are necessary to obtain reliable results for plane and
that the time for computation rapidly increases as the number of scatterers
increases.
and when several physical problems are coupled (see chapter 6). Both plane
and cylindrical incident waves were simulated and the Insertion Loss (IL)
calculated for each incident wave was compared with data for an equivalent
found to agree well with the laboratory data in term of predicting the bandgaps
(see figure 4.5). The computation time for FEM is much longer than that for
MST for the same SC configuration. Nevertheless, the computation time has
been found to be reasonable period even for the case with the largest number
results have shown that increasing the number of elements in the FEM model
shapes scatterer (triangle, square and elliptical) with similar cross section area
performances were compared with that of the circular shape scatterers and,
according to the FEM results, the triangular shape scatterer with arrangement
3 (see figure 4.8) offers a better IL performance than the other shapes.
Comparisons made for hexagonal and square lattice arrangement for circular
scatterers, have shown that the square lattice is more suitable as for road
traffic noise barrier design. An efficient FEM for calculating the acoustical
performed also for SC with vacancies which potentially could be used to guide
systems and signal processing. A robust indoor measuring system has been
set up to measure the insertion loss of the sonic crystal in the anechoic
(MLSSA). Studies on the effect of transducer size and post processing of the
measured signal (i.e. different window sizes and functions) were reported. To
and Split Ring Resonator) and the supporting structure in the chamber were
taken at the Open University barrier test site. For these outdoor
of MLSSA. Two contributing factors for this change are; the question of the
the source, a modified white noise (termed as Traffic Noise signal in this
Thesis) was used instead of an MLS signal. A Graphical User Interface (GUI)
was developed using MATLAB® for enabling the user to control the data
screen and the associate time or frequency domain signal spectrum taken
during the measurement will then be plotted on the display. Saving the data is
measurement based on the outputs from the sonic anemometer and has
chapter, an alternative swept sine source was also considered. The outdoor
measuring system (except that the source is replaced with a natural made
out.
rubber (latex). The additional bandgap due to the breathing mode resonance
of the shell is predicted to be the widest among those generated by the shell
resonances. Proof of concept for the breathing mode was performed using
found to offer a fast and reliable approach to solving the problem of scattering
material as well as radius and thickness of the shell so that the first resonance
frequencies could be positioned below the first Bragg bandgap of the array.
The computational time required by using the MST is less than one tenth of
that needed for traditional numerical methods such as FEM for modelling a
single and arrays of elastic shells in the low-frequency range was verified by
barriers suitable for traffic noise. Laboratory experiments using periodic arrays
Resonator (SRR) which is made from single slit cylindrical cylinder was
first Bragg band gap. The re-visiting of this work also served to provide a
better understanding of the possible influences (i.e. slit width, area of cavity
and slit orientation) on the performance of the SRR which were studied both
performed on the full scale SC barriers made with SRR and their performance
elastic shells arrange inside the rigid 4-slit pipes) and suggested that it is
possible both to protect the elastic cylinders and, effectively, to broaden the
band gaps. The IL of a 7x3 array of 55 mm outer diameter PVC cylinders with
4 mm wide slits and an inner latex cylinder with 43 mm outer diameter and
0.25 mm thick walls was measured and demonstrated the efficacy of this
with both circular and annular cavities. With an inner elastic shell, an
the presence of the outer 4-slit pipe. This is due to the coupling between the
represent the solution inside the slits. The method can be applied to various
arrays the concentric arrangements with inner elastic shells and outer 4-slit
below the first Bragg band gap while still preserving the Bragg band gaps.
elastic shells, hence is a potentially useful basis for a sonic crystal barrier
design. Full scale field trials of such a sonic crystal barrier (using 3 rows each
conventional noise fence were performed on a test site that has been
the numerical modelling results (see figure 9.1) indicated that the outdoor
sonic crystal band gaps produced by the effect of multiple scattering from
function method) were used to yield the Excess Attenuation (EA) spectra for
Those of asphalt and grass covered ground were measured in-situ. Fits to
the ground planes used were either acoustically-hard (MDF and glass) or –
has been derived for predicting multiple scattering effects of a finite arrays of
cylinders with axes parallel to either rigid or impedance ground. Results of the
semi-analytical method have been compared with laboratory data. 7x3 square
lattice arrays consisting of both rigid and elastic scatterers have been placed
with their axes parallel to the grounds and the insertion losses at various
have shown also that it is possible to augment the barrier effect of a sonic
crystal array by means of the ground effect i.e. the interference between
sound travelling directly from source to receiver and that reflected from the
ground beneath the barrier. The data show that sonic crystal band gaps are
much reduced when their frequency coincides with excess attenuation minima
destruction of the positive IL peak associated with the first Bragg band gap if
effect minima to lower frequencies so that the Bragg band gap is maintained.
This work revealed that it should be possible for ground effect and SC stop
axis SC arrays over an acoustically ‘hard’ ground were found to show the
frequency bands, called band gaps related to the spacing and size of the
SCs and explores means of enhancing the stop bands notably in the low
Multiple Scattering Theory (MST) and Finite Element Method (FEM) have
scattering that can lead to band gaps, strategies for improving the band gaps
for example (i) making each cylinder in the form of a Split Ring Resonator
(SRR) and (ii) substituting the acoustically rigid scatterer with thin elastic shell.
Such designs are shown to increase Insertion loss (IL) in the low-frequency
below the Bragg resonance stop band. With the better knowledge of the
factors affecting the frequency response for both the SRR and thin elastic
band gaps for SC barriers and the significant improvements when locally
resonant scatterers are used. Such resonant arrays are found to give rise to
acoustically “hard” or “soft” ground plane, have been made using predictions
and laboratory experiments. Some data have also been obtained also for
barriers.
Possible areas of research and further research along the lines of the
structure in the present study are for a either a plane or a fixed single
explored.
interesting to know how the diffraction at the ends of the finite cylinders
which has been widely researched recently [143 and 144] to increase
(which are protected from the rigid cylinder) work in principle they are
interesting to study the performance along bends which will cause the
5) The parameter values chosen for the full scale prototype are far from
the outer cylinders (grouted by hand in the prototype) are too wide and
the 0.043 m diameter of the latex cylinders, although the largest that
40
MST Prediction (current design SRR4)
MST Prediction (proposed modifications)
30
20
10
-10 2 3
10 10
Frequency (Hz)
to alter the shape of the outer rigid scatterer from circular to triangular
Plane Wave Expansion MATLAB® script modified from codes obtained from
tic;
num_Eigenvals = 10;
nmax = 10;
n1=(-nmax:nmax);
n2=(-floor((2*nmax+1)^2/2):floor((2*nmax+1)^2/2));
lat_a=0.135;
ro=0.0275;
f=pi*ro^2/lat_a^2;
pa=1380;
pb=1.2;
delp=(pb/pa-1)/(f*pb/pa+1-f);
va=2400;
vb=344;
pcla=(va^2)*pa;
pclb=(vb^2)*pb;
delt=(pclb/pcla-1)/(f*pclb/pcla+1-f);
nx_arr=repmat(n1,1,2*nmax+1);
ny_arr=round(n2/(2*nmax+1));
[nx,nxp]=meshgrid(nx_arr, nx_arr);
[ny,nyp]=meshgrid(ny_arr, ny_arr);
fill =ones(1,length(nx_arr));
fillt =diag(fill);
G_Gp=sqrt(4*pi*f*((nx-nxp).^2+(ny-nyp).^2))+fillt;
fG_Gp=2*f*besselj(1,G_Gp)./G_Gp;
mode=1;
if mode==1
for kyind=1:3
onoff=1;
if kyind == 1
kxt=kx;
kx=.55-kx;
ky=kx;
XXIX
elseif kyind==2
kxt=kx;
ky=0;
elseif kyind==3
kxt=kx;
ky=.5;
end
mag_k_p_g=((kx+nx).^2+(ky+ny).^2);
k_p_g_dot_kp_p_g=(kx+nx).*(kx+nxp)+(ky+ny).*(ky+nyp);
kro_del_ggp=(nxp==nx)&(nyp==ny);
M=mag_k_p_g.*kro_del_ggp+delp.*fG_Gp.*k_p_g_dot_kp_p_g.*(1-kro_del_ggp);
N=kro_del_ggp+fG_Gp.*delt.*(1-kro_del_ggp);
A=N^-1*M;
[a,b]=eig(A);
g=find(b~=0);
eigs=(sort(b(g)));
eigs(1:num_Eigenvals);
k=repmat(kxt,1,num_Eigenvals);
min(b(g));
if onoff == 1
omega=[((340)/lat_a)*real(sqrt(eigs(1:num_Eigenvals)))'];
karr=[k];
if kyind == 1
omegafull=[((340)/lat_a)*real(sqrt(eigs(1:num_Eigenvals)))'];
karrfull=[k];
end
onoff = 0;
else
omega=[omega,((340)/lat_a)*real(sqrt(eigs(1:num_Eigenvals)))'];
karr=[karr,k];
omegafull=cat(1,omegafull,((340)/lat_a)*real(sqrt(eigs(1:num_Eigenvals)))');
karrfull=cat(1,karrfull,k);
end
end
if kyind==1
myfig = figure;
ymax=10000;
subplot(1,3,1)
b1=plot(karr,omega,'*');
ax1=gca;
po=get(gca, 'position');
set(gca, 'YLim', [0,ymax]);
set(gca, 'XLim', [0,.5]);
set(gca, 'nextplot', 'add');
left_side=po(1);
set(gca, 'xtick', []);
width=po(3);
right_side=left_side+width;
set(gca, 'xticklabel', {});
ylabel('Frequency (kHz)')
elseif kyind==2
subplot(1,3,2)
b2=plot(karr,omega,'*');
po=get(gca, 'position');
po(1)=right_side;
XXX
ax2=gca;
set(gca, 'nextplot', 'add');
set(gca, 'YLim', [0,ymax]);
set(gca, 'XLim', [0,.5]);
set(gca, 'ytick', []);
left_side=po(1);
width=po(3);
right_side=left_side+width;
set(gca, 'xticklabel', {});
elseif kyind==3
subplot(1,3,3)
b3=plot(karr,omega,'*');
po=get(gca, 'position');
po(1)=right_side;
ax3=gca;
set(gca, 'nextplot', 'add');
set(gca, 'ytick', []);
set(gca, 'YLim', [0,ymax]);
set(gca, 'XLim', [0,.5]);
left_side=po(1);
width=po(3);
right_side=left_side+width;
set(gca, 'xticklabel', {});
text(-1,-.05*ymax,'M');
text(-.5,-.05*ymax,'\Gamma');
text(0,-.05*ymax,'X');
text(.5,-.05*ymax,'M');
text(-.3,-.1*ymax,'Reduced Wave Vector', 'HorizontalAlignment','Center');
text(-.3,1.075*ymax,['Disperaion Relation for Square Lattice of PVC Cylinder in
Air'],'HorizontalAlignment', 'Center');
text(-.3,1.025*ymax,['f=',num2str(f),'
a=',num2str(lat_a),'mr_{o}=',num2str(ro)],'HorizontalAlignment', 'Center');
end
end
end
to
XXXI
Appendix B
concepts for the multiple scattering theories. Generally speaking, they are
used to transform one expansion about some point in space into a similar
expansion about a different point. Since our multiple scattering model is two-
For clarity, the Graf’s additional theorem will be explained based on our
origins at the j-th and q-th scatterers. Let rk be the position vector of a
XXXII
y
j-th scatterer at (xj, yj)
Receiver (Px, Py)
rj
j j
rq
Rj
Rjq q
n
J m krj e J kR e J m n krq e
im j
in 2 j i m n q
n jq
n
a1 , (B1.1)
n
J kR e J n krq e
i m n 2 j in q
mn jq
n
XXXIII
n
H m1 krj e H kR e J n krq e
im j
i m n 2 j in q
1
mn jq a1 , (B1.2)
n
For rq R jq , and
n
H m1 krj e J kR e H n1 krq e
im j
i m n 2 j in q
mn jq a1 , (B1.3)
n
For rq rjq .
Where H n1 is the n-th order Hankel function of the first kind.
XXXIV
Appendix C
Lay Flat Tubing (LFT). Perspex end caps with attaching bolts (for mounting
purpose) were fitted to the LFT as shown in figure C1. A small hole is made
through one end of the end cap and a small pipe is glued to the opening.
Hoses are then attached to all the cylinders in order to inflate then with air
using a small pump. Measurements of single and 7x3 square lattice array of
(a) (b)
End cap
Air hose
Figure C1: (a) Single LFT and (b) 7x3 square lattice array of LFT with lattice constant
0.15 m.
XXXV
Results:
0.055 m but different wall thicknesses of 0.0001 and 0.00018 m. Clearly, very
scatterers.
2
Insertion Loss (dB)
1.5
0.5
0 2 3
10 10
Frequency (Hz)
Figure C2: IL spectra for single LFT cylinders of diameter 0.055 m and wall
thicknesses (WT) of 0.0001 m and 0.00018 m.
0.55 m, WT 0.0001 m and lattice constant 0.15 m (see figure C3). Only the
classical band gap (Bragg) at 1.1 kHz is observed in the plot and we conclude
XXXVI
12
10
8
Insertion Loss (dB)
0 2 3
10 10
Frequency (Hz)
Figure C3: IL spectra for 7x3 square lattice array of LFT with diameter 0.055 m, wall
thickness (WT) of 0.0001 m and lattice constant (L) of 0.15 m.
elastic thin-walled cylinders by periodic stiffening along their lengths has been
and gluing together multiple vertical sections using “gaffer” tape. The width of
the tape sections is 0.05 m and they are periodically spaced at 0.05 m
throughout the length of the elastic cylinder as shown in figure C4. The IL
spectrum for a single elastic cylinder with periodic stiffening has been
XXXVII
Figure C4: Single elastic cylinder with periodic stiffening mounted inside the
anechoic chamber.
resonance frequency also reduced from around 1.2 kHz to 800 Hz by the
increased stiffness. This laboratory test indicates that such a design is not
5
Measurement, latex 0.055 m (OD) and 0.00025 m (WT) with periodic stiffening
Measurement, latex 0.055 m (OD) and 0.00025 m (WT)
4
Insertion Loss (dB)
-1
-2 2 3 4
10 10 10
Frequency (Hz)
Figure C5: IL spectra of single elastic cylinder with (solid line) and without (broken
line) periodic stiffening.
XXXVIII
Appendix D
Tensile test
cylinders made from latex and polyethylene. Specifically the Young’s Modulus
which is defined as the ratio of the uniaxial stress over the uniaxial strain in
the range of stress in which Hooke’s Law holds has been obtained. For most
materials, during the initial portion of the test, the relationship between the
applied force and the elongation (change in length) is linear (see figure D1).
The region is known as the elastic region since here the material would return
to its original shape if the load is removed. Beyond this elastic limit, the
material starts to deform plastically and therefore the region is call the plastic
force will lead to a decrease of load with increasing elongation and eventually
XXXIX
A rectangular strip of non-vulcanized latex or polyethylene strip was cut to
mounted on an MTS universal test system. The tensile test involved applying
a continually increasing tensile force at one end of the specimen. The other
end of the specimen was fixed to ensure that the sample is slowly pulled apart
and the corresponding deformation (elongation) using the MTS Test Star II
data acquisition software. The inverse problem (relaxation) was performed for
Specimen
Figure D2: Tensile test experiment arrangement for the latex specimen using MTS
Universal Testing System.
of the latex specimen with wall thickness of 0.5 mm. The dotted line shows
either the data for either tension or relaxation and the best fit is obtained by a
linear regression method. As we can see the curve contains a straight line
where the load is proportional to the displacement. The gradient of the line
XL
gives the Young Modulus which is 0.89 and 0.81 MPa for tension and
relaxation respectively. The slight deviations of the data between tension and
mounting, precision of the stepper motor during the testing or the sensitivity
0.35
0.3
y = 0.8126x + 0.0308
2
R = 0.9823
0.25
E = 0.81 MPa
Tensile stress (MPa)
0.2
0.15
y = 0.8873x - 0.0013
0.1 R2 = 0.9937
E = 0.89 MPa
0.05
0
0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4
-0.05
Tensile strain
Axial Stress (relaxation) Axial Stress (tension)
Linear (Axial Stress (relaxation)) Linear (Axial Stress (tension))
Figure D3: Results of tests to obtain the Young’s Modulus of industrial rubber
Sample size: length 60 mm, width 25 mm and thickness 0.50 mm.
thickness of 0.15 mm. The dotted line shows the data for tension and the best
fit is again obtained by a linear regression method. As we can see the data
displacement. The gradient of the line gives the Young Modulus which is 178
MPa for tension. This result suggests that the specimen is a low density
XLI
Polyethylene where, according to a reference from public source, the Young’s
Modulus is around 200 MPa. This test supports the observed lack of breathing
1.8
y = 177.83x - 0.0419
1.6
R2 = 0.9977
1.4
Tensile stress (MPa)
1.2
0.8
0.6
0.4
0.2
0
0 0.002 0.004 0.006 0.008 0.01 0.012
-0.2
Tensile strain
Axial Stress (tension)
Linear (Axial Stress (tension))
Figure D4: Young’s Modulus of Polyethylene tubing. Sample size: length 60 mm,
width 25 mm and thickness 0.1 mm.
XLII
Appendix E(i): Meteorological conditions - Compass plot for outdoor measurement on 12/5/11
SC and Fence, 12/5/11 (1)
0 SC and Fence, 12/5/11 (2) SC and Fence, 12/5/11 (3)
0 SC and Fence, 12/5/11 (4)
30 330 0 0
30 330 30 30 330
330
Average wind speed = 0.9m/s
60 300 60
300 60
300 60 300
Average temp = 17.1°C
1
Barrier
0.8
Barrier
Barrier Barrier
Average SOS = 341.4m/s
0.8 0.6 0.6 1
0.4 0.4 0.8 1.5
0.2 0.2 0.6 0.4 1
90 270 90 0.2 0.5
270 90
270 90 270
60 300 60 300 60 60
300 300
XLIII
Appendix E(ii): Meteorological conditions - Compass plot for outdoor measurement on 27/5/11
SC and Fence, 270511 (1) SC and Fence, 27/5/11 (3) SC and Fence, 27/5/11 (4)
0 Fence and SC, 27/5/11 (2) 0
0
0
30 330 30 330
30 330 30 330
Average wind speed = 0.7m/s
60 300 60
300
60 300 60 300
Average temp = 15.6°C
Barrier Barrier Barrier Barrier
Average SOS = 340.5m/s
0.8 0.8 1 1 0.8
0.6 0.6 0.8 0.6 0.6
0.4 0.4 0.4 0.4
0.2 0.2 0.2 0.2
90 270 90 90 270 90 270
270
XLIV
References
[2] P.M. Moses and K.U. Ingard, “Theoretical acoustic”, McGraw Hill
(1968).
[4] Y. de Kluizenaar, S.A. Janssen, F.J. van Lenthe, H.M.E. Miedema and
J.P. Mackenbach, “Long-term road traffic noise exposure is associated
with an increase in morning tiredness”. J. Acoust. Soc. Am., 126 (2),
pp. 626-633 (2009).
[6] W.F. Babisch “Road traffic noise and cardiovascular risk”, Noise &
Health, 10, pp. 27-33 (2008).
[7] B.J. Smith, R.J. Peter and S. Owen, “Acoustics and noise control”,
Longman (1985).
XLV
[14] X. Hu and C. T. Chan, “Photonic crystals with silver nanowires as a
near-infrared superlens”. Applied Physics Letters, 85 (9), 1520 (2004).
[16] A. Mekis, J.C. Chen, I. Kurland, S. Fan, P.R. Villeneuve and L.D.
Joannopoulos, “High transmission through sharp bends in photonic
crystal waveguides”. Physical Review letters. 77, 3787–90 (1996).
[17] D. Schurig, J.J. Mock, B.J. Justice, S.A. Cummer, J.B. Pendry,
A.F. Starr, D.R. Smith, “Metamaterial Electromagnetic Cloak at
Microwave Frequencies”. Amer. Assoc Advancement Science, 314
(5801), 977-980 (2006).
[23] A. Gupta, C.H. Chew and K.M. Lim, “Effect of periodic structure on
sound propagation”. 4th International Conference on experimental
mechanics, 7522 - 75223D (2010).
XLVI
[26] D.P. Elford, L. Chalmers, F.V. Kusmartsev and G.M. Swallowe,
“Matryoshka locally resonant sonic crystal”. J. Acoust. Soc. Am., 130
(5), 2746-2755 (2011).
[27] C. Kittel, “Introduction to solid state physic”, John Wiley & Sons Inc, 8
editions (2005).
[32] R.D. Meade, K.D. Brommer, A.M. Rappe, and J.D. Joannopoulos,
“Existence of a photonic band gap in two dimensions”, Applied Physics
Letters, 61, 495-497 (1992).
[33] K.M. Ho, C.T. Chan, and C.M. Soukoulis, “Existence of a photonic gap
in periodic dielectric structures”. Physical Review Letters, 65, 3152-
3155 (1990).
[34] Y.M. Soliman, M.F. Su, Z.C. Leseman, C.M. Reinke, I. El-kady and
R.H.III Olsson, “Phononic crystals operating in the gigahertz range with
extremely wide band gaps”. Applied Physics Letters, 97, 193502
(2010).
XLVII
[39] M.S. Kushwaha, P. Halevi , G. Martínez, L. Dobrzynski and B. Djafari-
Rouhani, “Theory of acoustic band structure of periodic elastic
composites”. Physical Review B, 49, 2313-2322 (1994).
[44] M.I. Mishchenko, D.W. Mackowski and L.D. Travis, “Scattering of light
by bispheres with touching and separated components”. Applied
Optics, 34, 21, 4589 (1995).
[48] L.L. Araujo, G.F. Foran and M.C. Ridgway, “Multiple scattering effects
on the EXAFS of Ge nanocrystal”. Journal of Physics: Condense
Matter, 20, 165210 (2008).
[49] C.M. Linton and D.V. Evans, “The interaction of waves with arrays of
vertical circular cylinders”. J. Fluid Mech., 215, pp549-569 (1990).
[51] J.M. Elson, “Infrared light scattering from surfaces covered with
multiple dielectric overlayers”. Applied Optics 16 (11), 2872–2881
(1977).
XLVIII
[52] A. Aubry, A. Derode, “Multiple scattering of ultrasound in weakly
inhomogeneous media: Application to human soft tissues”. J. Acoust.
Soc. Am., 129 (1), 225-233 (2011).
[55] A.K. Mal and S.K. Bose, “Dynamic elastic moduli of a suspension of
imperfectly bonded spheres”. Proceeding of the Cambridge
Philosophical Society, 76, 587-600 (1974).
[58] I.S. Gradshteyn and I.M. Ryzhik, “Table of integrals, series and
products”, 6 editions, 8.51-8.52, MO 27, 926 (2000).
[60] R.D. Cook, D.S. Malkus, M.E. Plesha and R.J. Witt, “Concepts and
applications of finite element analysis”. John Wiley & Sons, 4 editions,
Chapter 1, pp 10 (2002).
[62] B.R. Mace, D. Duhamel and M.J. Brennan, “Finite element prediction of
wave motion in structural waveguides”. J. Acoust. Soc. Am., 117(5),
2835-2843 (2005).
[64] T. Delort, “Finite-element method for gratings”. J. Opt. Soc. Am., A, 10,
12, 2592-2601 (1993).
[65] X. Wei, A.J. Wachters and H.P. Urbach, “Finite-element model for
three-dimensional optical scattering problems”, J. Opt. Soc. Am., A, 24,
3, 866-881 (2007).
XLIX
[66] C. Lally, F. Dolan and P.J. Prendergast, “Cardiovascular stent design
and vessel stresses: a finite element analysis”. Journal of
Biomechanics, 38 (8), 1574-1581 (2005).
L
[79] T. Miyashita, W. Sato and Y. Nakaso, “Experimental studies on two-
dimensional defect-mode waveguides in a sonic/phononic crystal”,
Japanese Journal of Applied Physic Part 1-Regular papers brief
communications & review papers, 46 (7B), 4684-4687 (2007).
[83] J.D. Turner and A.J. Pretlove, “Acoustic for engineers”, Macmillan
(1991).
[89] BS EN 1793-6, “Road traffic noise reducing devices. Test method for
determining the acoustic performance. Part 6. Intrinsic characteristics.
In situ values of airborne sound insulation under direct sound field
conditions” (2010).
[90] F.J. Harris, “On the use of windows for Harmonic Analysis with the
Discrete Fourier Transform”. Proceeding of the IEEE, 66, 1 (1987).
LI
[92] A. Farina, “Simultaneous measurement of impulse response and
distortion with a Swept-sine technique”. 108th Audio Engineering
Society Convention, 5093 (D-4) (2000).
[97] F.G. Wu, Z.L. Hou, Z.Y. Liu, and Y.Y. Liu. “Point defect states in two-
dimensional phononic crystals”. Physics Letter A, 292:198 (2001).
[98] M. Sigalas, “Eastic wave band gaps and defect states in two-
dimensional composites”. J. Acoust. Soc. Am., 101, 1257 (1997).
[100] M. Hirsekorn, P.P. Delsanto, N.K. Batra and P. Matic, “Modelling and
simulation of acoustic wave propagation in locally resonant sonic
materials”. Ultrasonics, 42, 231-235 (2004).
[101] M.S. Kushwaha and P. Halevi, “Stop bands for cubic arrays of
spherical balloons”. J. Acoust. Soc. Am., 101, 619-622 (1997).
LII
[105] N.A. Nicorovici, R.C. Mcphedran, “Photonic band gaps for arrays of
perfectly conducting cylinders”. Physical Review E, 52, 1135-1145
(1995).
[107] J.D. Kaplunov, L. Yu. Kossovich and E.V. Nolde, “Dynamic of thin
walled elastic bodies”. Academic Press London, 110 (1998).
[112] J.B. Pendry, A.J. Holden, W.J. Stewart, and I. Youngs, “Magnetism
from conductors, and enhanced non-linear phenomena.” IEEE Trans.
Microw. Theory Tech., 47, 2975 (1999).
[114] L.Y. Wu, L.W. Chen, M.L. Wu, “The nondiffractive wave propagation in
the sonic crystal consisting of rectangular rods with a slit”, Journal of
physics – Condensed Matter, 20, 295229 (2008).
LIII
[118] W. Elmore and M. Heald, “Physics of Waves”, New York: McGraw-Hill,
p. 148 (1969).
[119] C.H. Sohn and J.H. Park, “A comparative study on acoustic damping
induced by half-wave, quarter-wave and Helmholtz resonators”.
Aerospace Science and Technology, 15, 606 (2011).
[121] F.P. Mechel, “Formulas of acoustics”, Springer Verlag, New York Inc
(2008).
[122] S.G.L. Smith and A.M.J. Davis, "The split ring resonator". Proceedings
of the Royal Society A, 466, 3117-3134 (2010).
[123] Z.Y. Cui, T.N. Chen and Y.P. Su, “Experimental and calculated
research on a large band gap constituting of tubes with periodic narrow
slits”. Applied Acoustics 70, 1087 (2009).
LIV
[132] W. K. Lui and K.M. Li, “The scattering of sound by a long cylinder
above an impedance boundary”. J. Acoust. Soc. Am., 127 (2), 664-674
(2010).
[134] G. Watts, “Harmonise prediction model for road traffic noise”, Transport
Research Laboratory report, PPP034 (2005).
[137] M.E. Delany and E.N. Bazley, “Acoustic properties of fibrous absorbent
materials”. Applied Acoustics, 3, pp 105-116 (1970).
[145] M.R. Monazzam and Y.W. Lam, “Performance of Profiled Single Noise
Barriers Covered with Quadratic Residue Diffusers”, Applied Acoustics
66, pp.709-730 (2005)
LV
[146] J. Berenger "A perfectly matched layer for the absorption of
electromagnetic waves". Journal of Computational Physics 114 (2),
185–200 (1994).
LVI