Chapter - 7c Directivity Index
Chapter - 7c Directivity Index
Chapter - 7c Directivity Index
The axial pressure amplitude for the far field (Fraunhofer zone) is given by (for θ = 0)
p ( r, 0) kr a2 p ( r ,0,0) 1a
= sin 1 + 2 − 1 and = ka
2 ρo cUo 2 r 2 ρocU o 4 r
where the solid line represents a λ = 3 , the dashed line represents the a λ = 7 . Did one reduce to the
other? Why not?
So,
k
k
− j (R 1 +r )
p ( r , t ) = 2 j ρocU oe jω t e 2 sin ( R1 − r ) where R1 = r 2 + a 2
2
and
jωρo a 2Uo j (ω t− kr ) 2 J1 ( ka sinθ )
p ( r, θ , t ) = e
2r ka sinθ
Homework EXERCISE: Show that these two expressions are the same, on-axis, in the Fraunhofer
Zone.
ANSWER: The equations for the line source and circular transducer are listed below:
kL sin θ
H L (θ ) = sinc (Line source)
2
2 J ( ka sinθ )
HC (θ ) = 1 (Circular transducer)
ka sinθ
The plots for the two cases are shown below. Note that the continuous-line source has the x-axis
1
“expanded out” by a factor of two because of the factor of 0.5 in the expression kL sin (θ ) . Thus,
2
one can immediately see that the circular transducer (i.e. the circular piston source) has greater
directivity, but at the same time has more side lobes than the continuous-line source. Also note that
L for the line source is the “full length” whereas L (a) for the piston source is “1/2 the length.”
Log Comparison of Directivity
10
Continuous-line Source
0 Circular Piston Source
-10
-20
20 log |H(theta)|
-30
-40
-50
-60
-70
-80
0 5 10 15
kL*sin(theta)
0.80
0.60
|H(theta)|
0.40
0.20
0.00
0 5 10 15
kL*sin(theta)
How would the plots be different if L = D , the diameter of the circular transducer?
******************************************************************
Oelze ECE/TAM 373 Notes - Chapter 7 pg 28
The far field pattern can be modified by changing U0 over the face of the source. For example, if U0
decreases toward the edges of the circular piston then the side lobe levels will decrease and the main
lobe will be broader.
Ex 1 – High frequencies in a loudspeaker move the center more than outer parts.
Transducer dimension (thickness) changes in response to electric field. “Shading” lowers the field
in the outer portions of the disk.
Beamwidth
Defined as angular width where the intensity is down a specified amount from that on the axis.
Ex 0.5 → 3 dB down
0.25 → 6 dB down
0.1 → 10 dB down
Thus if θ = 13° for the acoustic amplitude down by 10 dB then the –10 dB beam width is 2θ = 26°.
For –10 dB the intensity directivity function will have a value of 0.1
2
Alternatively, at a specific distance from the source the beam width might be specified in terms
of a linear distance as opposed to an angular width.
I
I0
0.5 I0
Beam patterns for circular sources with uniform velocity amplitude distributions.
6 z0
0 1/2 z0 z0 3 z0
2θ1
p
Spherical
Waves-Field
comp. ⊥ surface
→
dfs df s
zr = ∫ r
u
The mechanical impedance of the transducer is analogous to discussion from Chapter 1 for the mass
loaded, damped driven harmonic spring system.
Piston → mass m
mechanical resistance Rm
stiffness s
jωt r
Driven by external force f = Fe to give u0 = U 0 e jωt = jωξ 0 for simple harmonic time dep.
We can plug all of this into a second order differential equation (just like the damped, driven
harmonic oscillator)
r r dξ 0 d 2 ξ0
f – f s – Rm – sξ 0 = m 2
dt dt
or
du0
f = f s + Rmu0 + s ∫ u0dt + m
dt
We can assume a simple harmonic time dependence like the driving force so that
s
f = f s + R m + jω m + u0
jω
= ( z%r + z%m ) u% 0
where
f%
u%0 =
z%m + z%r
jωt
dividing by e gives
F
U0 =
z% m + z% r
where
N–s
z% r = ( zr ) e jθ = Rr + jX r
m
↑ ↑
radiation radiation
resistance reactance
1
Power radiated Π r = U 20 Rr
2
Let’s look at the Simple Source
ρ0 c( kS ) 2
Rr = unbaffled, where S is the surface area
4π
ρ c( kS ) 2 baffled simple source
Rr = 0
2π ( twice Rr and Π r )
_________________
Looking at Ztotal = Zm + Z r
s
= Rm + j ωm – + Rr + jX r
ω
a positive Xr means that we are loading the system with mass from the surrounding fluid so that the
resonance frequency is decreased for the oscillator.
Xr
Thus, we define a radiation mass m r = .
ω
Since the resonance frequency of the system is
The radiation mass reduces the resonance frequency. mr is related to mass of material moved in
front of source. Less mass at higher frequency.
What does this mean for a transducer operating in air compared to a transducer operating in water?
Rr = ρ0cSR1 ( 2ka )
X r = ρ0cSX 1 ( 2ka )
where
2 J 1 ( 2ka )
R1 ( 2 ka ) = 1 −
2ka
2 H1 ( 2 ka )
X1 ( 2ka ) = (To be shown in homework)
2 ka
and S = π a (the area of the piston face).
2
1.0
R1(x)
0.5 X1(x)
Figure 7.5.2 Radiation resistance and
reactance for a plane circular piston of
radius a radiating sound of wave
0 number k(x = 2ka).
0 4 8 12
x
For ka << 1 (a small compared to λ)
Physically, this means that the piston appears to be loaded with a cylindrical volume of fluid mass
with cross sectional area S (the same as the piston surface) and a depth of the fluid mass cylinder
8a
reaching to z = ; 0.85a .We see how the fluid changes the resonant frequency.
3π
Oelze ECE/TAM 373 Notes - Chapter 7 pg 33
For ka >> 1 (a large compared to λ)
Again we use the series expansion for J and H at large ka.
Rr ≅ ρ0cS (1)
2
Xr ≅ ρ0cS ⇒ 0 for ka very large
π ka
2ρ a
mr ≅ 02 ⇒ 0 for ka very large
k
Thus, at large ka or the radiation impedance is almost real with no mass loading
a
λ
So, the power is given by
1 1
Π r = U 20 Rr ≅ ρ0cSU 20 .
2 2
This is the same power as a plane wave propagating through area S with characteristic impedance
ρ0 c .
Πr = 1
U
2
Rr
0.61 1,316 (W)
2 0
X1 0.1464 0.693
Xr 0.172 2,861
mr 2.73 × 10–6 4.56 × 10–2 (kg)
Physically it takes a lot more work done against Rr in water than in air. Hence, the more power
expended in water than in air by the transducer.
******************************************************************
The general expression for beam pattern, b(θ , φ ) , represents the variation of intensity level with
angle. The expression is given by
I ( r ,θ ,φ )
b(θ ,φ ) = 10 log
axI (r ) Note: is generally a function of both θ and φ
P( r ,θ ,φ ) since it may not have circular symmetry.
= 20 log
pax (r )
= 20 log ( H (θ ,φ ) )
(d) Directivity
nonspherical
I ax ( r ) directional
source P ax2 (r )
directivity = D = =
I S ( r )(spherical) P S2 (r )
for a simple source having the same acoustic power as the directive source. For the simple
source
1
Πs = 4π r 2 P S2 ( r )
2 ρ0 c
and there is no preferred direction. For the directional source
1
2 ρ0c ∫4π
Π directional = P 2 ( r ,θ , φ ) r 2 d Ω where dΩ = sin(θ) dθ dφ and the integration is
1 2 2
Π directional = r P ax ( r ) ∫ H 2 (θ ,φ )d Ω
2 ρ0 c 4π
Since the acoustic powers are the same, we can solve for P ax2 (r ) and P 2S (r ) in terms of the
acoustic power giving
2Π ρ0c
r2 ∫ H 2 (θ ,φ )d Ω
2
P ax
( r)
4π
∴D = = 4π
=
2Π ρ0c
∫ H (θ ,φ ) d Ω
2
P S2 4π
4π r 2
The larger the ratio a/λ the more directional the beam since ∫
4π
H 2 (θ ,φ )d Ω is smaller.
4π 4π
D= =
∫ H (θ ,φ ) d Ω
2
2 J1 (ka sin θ )
2
π 2
4π
∫
0 ka sin θ 2π sinθ dθ
(ka ) 2
D= again Bessel functions are tabulated in appendix A6
J (2 ka)
1– 1
ka
for ka >> 1
2
D ≅ (ka)
which is greater for greater a/λ, or the beam is more highly directive when the frequency is
increased
10
Take a λ = > 1 ka = 10 D ≅ 100
2π
15
aλ= > 2 ka = 15 D ≅ 225
2π
Thus the source starts to become quite directive even for a/λ not much greater than 1.
******************************************************************
Reversible Transducers
The extent of the near field for a source is defined by the two distances rmin (the distance from the
field point to the closest element of the source) and rmax (the distance from the field point to the
furthest element of the source).
When the field point is moved towards the source such that r max – rmin ≈ λ 2 then the axial
pressure is sufficiently shifted from the far field axial response to alter the far field axial pressure
amplitude. The above relation can be shown through geometry to imply the extent of the near field
response is quantified by
2
L
rmin 1 L L2 2
≅ or rmin = =
L 4λ 4λ λ
L1 L2 > L1
L2 L 22
rmin =
4λ
λ
sin θ1 ≅
L θ
L
λ
θ1 decreases with L 2
Estimate of Directivity
4π
D≈ ∴ determine solid angle for main lobe–for side lobes low
Ω effective
1 2
D≅ k L1L2 L2
4
increases with dimensions and
L1
frequency