Linear System Modeling of Phased Arrays



Fig. 9.1
The elements of a phased array immersion system, where a pair of sending and receiving elements are shown





A314073_1_En_9_Fig2_HTML.gif


Fig. 9.2
The elements of a phased array contact testing system, where a pair of sending elements is shown

First consider the driving circuits in Figs. 9.1 and9.2. If we assume that these driving circuits can be represented by a linear active network, then Thévenin’s theorem says that in the frequency domain we can replace these driving circuits by an equivalent voltage source,
$${{V}_{i}}( \omega),$$
and electrical impedance,
$$Z_{i}^{\text{e}}( \omega),$$
as shown in Fig. 9.3 [Schmerr-Song]. These equivalent circuit elements will generate the same output voltage,
$${{V}_{1}}( \omega),$$
and current,
$${{I}_{1}}( \omega),$$
as the original driving circuit. Note that these equivalent parameters are affected by the choice of the phased array instrument settings that control the driving circuits.



A314073_1_En_9_Fig3_HTML.gif


Fig. 9.3
An equivalent voltage source and electrical impedance model of the driving circuits

The driving circuits are connected to the sending element through the cabling present in the sound generation process. The voltage and current
$$( {{V}_{1}},{{I}_{1}} )$$
of the driving circuits at the input side of the cable produce corresponding voltage and current,
$$( {{V}_{2}},{{I}_{2}} )$$
at the output side. At the typical MHz frequencies present in an ultrasonic experiment, the electrical characteristics of the cabling can be important if the cabling lengths are on the order of a meter or more in length since in these cases the cable does not just pass the inputs unchanged to the outputs. It is reasonable to assume the cable is a simple passive linear, and reciprocal device that can be treated as a two-port system [Schmerr-Song] where the inputs and outputs can be related through a 2 × 2 transfer matrix ,
$${{\mathbf{T}}^{\text{G}}}\text{,}$$
as shown in Fig. 9.4. More explicitly we can write



A314073_1_En_9_Fig4_HTML.gif


Fig. 9.4
A two-port representation of a cable as a 2 × 2 transfer matrix





$$ \left[ \begin{matrix} {{V}_{1}}\\ {{I}_{1}}\\\end{matrix} \right]=\left[ \begin{matrix} T_{11}^{\text{G}} & T_{12}^{\text{G}}\\ T_{21}^{\text{G}} & T_{22}^{\text{G}}\\\end{matrix} \right]\left[ \begin{matrix} {{V}_{2}}\\ {{I}_{2}}\\\end{matrix} \right], $$

(9.1)

where the elements of the transfer matrix are frequency dependent. Since the cable is assumed to be reciprocal, the determinant of the transfer matrix must be unity [Schmerr-Song], i.e.





$$ \det [ {{\mathbf{T}}^{\text{G}}} ]=T_{11}^{\text{G}}T_{22}^{\text{G}}-T_{12}^{\text{G}}T_{21}^{\text{G}}=1. $$

(9.2)

The array driving element is normally made of a piezoelectric material and is attached to both backing and facing materials. Its function is to transform the electrical inputs of the cable attached to the element,
$$( {{V}_{2}},{{I}_{2}} ),$$
into mechanical outputs at its face. The face is in contact with either a fluid (for the immersion case) or a solid (for the contact case). In the contact case a thin fluid layer between the face of the driving array and the underlying solid is normally present to ensure good coupling of sound into the solid . For both the immersion and contact cases the acoustic output properties of the element can be characterized by two “lumped” acoustic parameters—a compressive force,
$${{F}_{\text{t}}}( \omega),$$
and velocity,
$${{v}_{\text{t}}}( \omega),$$
both acting on the adjacent material. However, the detailed physics of the acoustic fields on the face of the ­element is modeled differently in the immersion and contact cases. In the ­immersion case , it is generally assumed that the face of the array at the element moves with the uniform velocity,
$${{v}_{\text{t}}}( \omega),$$
normal to the face of the transducer and generates a spatially varying pressure distribution,
$$p( \mathbf{x},\omega),$$
over its face. The force,
$${{F}_{\text{t}}},$$
is then just the integral of this pressure distribution over the active area,
$${{S}_{\text{A}}},$$
of the driving element, i.e.





$$ {{F}_{\text{t}}}( \omega)=\int_{{{S}_{\text{A}}}}{p( \mathbf{x},\omega)\text{d}S( \mathbf{x} )}. $$

(9.3)

In contrast, in the contact case, since the face of the array is in contact with a stiff solid material through a thin fluid couplant, it is more reasonable to assume that piezoelectric element generates a constant pressure,
$$p( \omega),$$
at its face and the total force is just
$${{F}_{\text{t}}}( \omega)=p( \omega){{S}_{\text{A}}},$$
while the velocity field normal to the face has the distribution
$$v( \mathbf{x},\omega).$$
In this case the “lumped” parameter,
$${{v}_{\text{t}}}( \omega),$$
can be taken to be the average velocity at the face:





$$ {{v}_{\text{t}}}( \omega)=\frac{1}{{{S}_{\text{A}}}}\int_{{{S}_{\text{A}}}}{v( \mathbf{x},\omega)\text{d}S( \mathbf{x} )}. $$

(9.4)

Thus, in either the immersion or contact testing setups, if we assume that the driving array elementA acts as a linear reciprocal system that transforms the electrical inputs
$$( {{V}_{2}},{{I}_{2}} )$$
to the mechanical outputs
$$( {{F}_{\text{t}}},{{v}_{\text{t}}} )$$
we can model this element also as a two-port system as shown in Fig. 9.5.



A314073_1_En_9_Fig5_HTML.gif


Fig. 9.5
A two-port model of an array element as a 2 × 2 transfer matrix

Explicitly, we then have





$$ \left[ \begin{matrix} {{V}_{2}}\\ {{I}_{2}}\\\end{matrix} \right]=\left[ \begin{matrix} T_{11}^{A} & T_{12}^{A}\\ T_{21}^{A} & T_{22}^{A}\\\end{matrix} \right]\left[ \begin{matrix} {{F}_{\text{t}}}\\ {{v}_{\text{t}}}\\\end{matrix} \right]. $$

(9.5)

Since the array element is assumed to satisfy reciprocity , we also have





$$ \det ( {{\mathbf{T}}^{A}} )=T_{11}^{A}T_{22}^{A}-T_{12}^{A}T_{21}^{A}=1. $$

(9.6)

At the acoustic output port of an array element the compressive force and velocity are not independent, reflecting the fact that if either the velocity on the face of the element or the pressure on the face of the element is specified then the other corresponding output parameter is also determined. This relationship is expressed as





$$ {{F}_{\text{t}}}( \omega)=Z_{\text{r}}^{A;a}( \omega){{v}_{\text{t}}}( \omega), $$

(9.7)

where
$$Z_{\text{r}}^{A;a}( \omega)$$
is theacoustic radiation impedance of the element. For ultrasonic immersion measurement systems that use large, single element transducers of circular cross section one can derive an explicit expression for this radiation impedance as [3])





$$ {Z_{\text{r}}^{A;a}}/{{{\rho }_{1}}{{c}_{p1}}{{S}_{\text{A}}}}\;=1-{[ {{J}_{1}}( 2ka )-\text{i}{{S}_{1}}( 2ka ) ]}/{ka}\; $$

(9.8)

where
$${{J}_{1}}$$
is a Bessel function of order one and
$${{S}_{1}}$$
is a Struve function. Here
$$( {{\rho }_{1}},{{c}_{p1}} )$$
are the density and compressional wave speed of the fluid,a is the radius of the transducer,
$$k=\omega /{{c}_{p1}}$$
is the wave number, and
$${{S}_{\text{A}}}=\text{ }\!\!\pi\!\!\text{ }{{a}^{2}}$$
. For most single element transducers
$$ka$$
can be a value of 100 or larger and Eq. (9.8) gives approximately
$$Z_{\text{r}}^{A;a}={{\rho }_{1}}{{c}_{p1}}{{S}_{\text{A}}},$$
which is just the acoustic impedance of a plane wave. While this same plane wave limit is also true for rectangular elements that are many wavelengths long in both dimensions, the rectangular elements often used in linear and 2-D arrays have at least one dimension which is not many wavelengths in size so that the acoustic radiation impedance will inherently be a complex function of frequency.

Since the output force and velocity on the face of an element are related through the acoustic radiation impedance, we can model the driving circuits, cabling, and sending element as a series of two-port systems terminated at both ports, as shown in Fig. 9.6a, and we can replace this series of systems by a single input, single output linear time-shift invariant (LTI) system whose input is the Thevenin equivalent voltage,
$${{V}_{\text{i}}}( \omega),$$
of the driving circuits. On the output side, we can take either the velocity,
$${{v}_{\text{t}}}( \omega),$$
or the compressive force on the face of the array element,
$${{F}_{\text{t}}}( \omega),$$
as the quantity to describe the acoustic fields. When modeling large, single element immersion transducers, it has been customary to use the force as the output [Schmerr-Song] where, as discussed above, the acoustic radiation impedance simply is a constant, i.e.
$${{F}_{\text{t}}}( \omega)={{\rho }_{1}}{{c}_{p1}}{{S}_{\text{A}}}{{v}_{\text{t}}}( \omega).$$
For array elements in immersion tests we will continue to take the output as
$${{\rho }_{1}}{{c}_{p1}}{{S}_{\text{A}}}{{v}_{\text{t}}}( \omega),$$
as shown in Fig. 9.6b but one should realize that this is really just an output proportional to the velocity having the dimensions of a force and not the actual compressive force on the face of the array element,
$${{F}_{\text{t}}}$$
which is given as
$${{F}_{\text{t}}}( \omega)=Z_{\text{r}}^{A;a}( \omega){{v}_{\text{t}}}( \omega).$$
We have made this choice because, as we will see, we can then develop complete models of ultrasonic phased arrays inspections without having to explicitly evaluate the radiation impedances of the elements. For array elements this evaluation must typically be done numerically, although it is not a difficult calculation for simple element shapes such as rectangles [4]. The transfer function for an immersion setup that describes the LTI system of Fig. 9.6b is denoted as
$${{t}_{\text{G}}}( \omega).$$
For a contact setup, since the pressure is modeled as given at the transducer face, we will instead take the output of our LTI model to be the actual force,
$${{F}_{\text{t}}}( \omega),$$
given by Eq. (9.3) and let
$${{{t}'}_{\text{G}}}( \omega)$$
be the transfer function in the contact case, as shown in Fig. 9.6c. Note that these transfer functions are explicit functions of all the system components shown in Fig. 9.6a but we will not need to write those explicit relationships here since
$${{t}_{\text{G}}}$$
and
$${{{t}'}_{\text{G}}}$$
will not need to be independently known in a complete description of the ultrasonic system, as will be discussed in Sect. 9.4.



A314073_1_En_9_Fig6_HTML.gif


Fig. 9.6
a The sound generation process components, andb their replacement for an immersion inspection by a single input, single output linear time shift invariant system characterized by the transfer function,
${{t}_{\text{G}}}( \omega),$
andc their replacement for a contact inspection by the transfer function
${{{t}'}_{\text{G}}}( \omega)$



9.2 Linear System Modeling and Sound Reception


The LTI models developed in the previous section can be used on sound generation for either immersion or contact phased array systems. In developing similar models for the sound reception processes it will be necessary to treat the immersion and contact cases separately. Consider first an immersion setup. Figure 9.7

Only gold members can continue reading. Log In or Register to continue

Stay updated, free articles. Join our Telegram channel

May 6, 2016 | Posted by in ULTRASONOGRAPHY | Comments Off on Linear System Modeling of Phased Arrays

Full access? Get Clinical Tree

Get Clinical Tree app for offline access