High Temperature Ultrasonic Thickness
High Temperature Ultrasonic Thickness
High Temperature Ultrasonic Thickness
taper
angle. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 32
21 Amplitude of the rst trailing echo relative to the rst arrival versus
taper angle. Each data point represents a separate waveform acquisition. 33
22 Waveforms resulting from nite dierence simulations for dierent ta-
per angles on a waveguide of length 152.4 mm (6) and diameter 9.54
mm (3/8) (a) straight rod. (b) 0.25
taper
angle. (d) 0.75
taper angle. . . . . . . . . . . . . 34
23 Impact of the taper angle on the trailing echo attenuation with the
nite dierence model. . . . . . . . . . . . . . . . . . . . . . . . . . . 35
24 Waveforms resulting from ray tracing simulations for dierent taper
angles on a waveguide of length 152.4 mm (6) and diameter 9.54mm
(3/8) (a) straight rod. (b) 0.25
taper angle.
(d) 0.75
taper angle. . . . . . . . . . . . . . . . . 36
25 Same as Figure 24 after zooming on the trailing echoes (a) straight
rod. (b) 0.25
taper
angle. (e) 1
taper angle. . . . . . . . . . . . . . . . . . . . . . . . . . 37
26 Impact of the taper angle on the trailing echo attenuation with the ray
tracing model. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 37
27 Comparison of trailing echo amplitude ratios from experiments and
both simulations. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 38
28 Change in arrival times of trailing echo rays after introduction of a
taper angle. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 39
29 Variability of the widening eect with the trailing echo order. (a)
First order trailing echo order. (b) Second order trailing echo. . . . . 41
x
30 Example of imperfect machining. . . . . . . . . . . . . . . . . . . . . 43
31 Comparison of two tapered rods with dierent quality of machining.
(a) Bad machining. (b) Good machining. . . . . . . . . . . . . . . . . 44
32 Comparison of the wave paths corresponding to the backwall echo for
the pulse-echo and pitch-catch congurations. . . . . . . . . . . . . . 47
33 Replacement of the calibration echo through notch digging. (a) Waveg-
uide with a notch. (b) Impact of the notch on the pulse-echo waveform. 49
34 Pulse-echo waveforms in cylindrical waveguides of length 152.4 mm
(6) (a) Rod of diameter 9.54 mm (3/8) with a notch of reective
area of 56% (d
n
=3.21 mm). (b) Rod of diameter 12.7 mm (1/2) with
a notch of reective area of 33% (d
n
=1.73 mm). . . . . . . . . . . . . 50
35 Pulse-echo measurements with a tapered rod. . . . . . . . . . . . . . 51
36 Design of a cylindrical waveguide enabling pulse-echo thickness mea-
surements. (a) Random waveguide and its pulse-echo waveform. (b)
Optimization of the rod dimensions for the rst arrival to be more im-
portant than it trailing echoes. (c) Notch dimensioning to minimize the
trailing echoes associated with the the calibration echo while increasing
the amplitude of the backwall echo. . . . . . . . . . . . . . . . . . . . 52
37 Thickness measurement in a pitch-catch conguration with waveguides. 54
38 Experimental setup for concept verication of the pitch-catch measure-
ments (a) Picture of the actual setup. (b) Parameters of the experiment. 56
39 Experimental waveforms for pitch-catch measurements. (a) Pulse-echo
waveform. (b) Through-transmission waveform. . . . . . . . . . . . . 56
40 Generation of trailing echoes in a cylindrical waveguide. . . . . . . . . 58
41 Thickness measurements at room temperature. . . . . . . . . . . . . . 58
42 Results of thickness measurements at elevated temperatures. The ini-
tial calibration was performed at 22.5
C. . . . . . . . . . . . . . . . . 59
xi
SUMMARY
The time required for inspection and maintenance of piping systems and vessels
needs to be reduced to both minimize down time and decrease operational costs of
petrochemical plants. Current ultrasonic inspection systems are not suited for on-
line monitoring, with the main issues being the resistance of transducers and their
coupling to high temperatures and the removal of insulation to access structures.
The use of welded cylindrical rods is thus proposed to transfer transducers and
their coupling to a lower temperature environment. However, when traveling through
waveguides, waves are subjected to numerous mode conversions that result in the
creation of new delayed echoes called trailing echoes. These echoes can impede
measurements and their attenuation needs to be investigated to enable use of wave-
guides. Among other solutions proposed in the literature, tapered rods are chosen
for their good results and their independence from material properties, which makes
them good candidates for welding. The impact of the taper angle on trailing echoes
is studied through both simulations and experiments. While clean signals can always
be obtained by increasing the taper angle in long rods, that is not always the case for
short rods, which are considered here.
In addition, temperature variations have a non-negligible impact on the arrival
time of the backwall echo when performing measurements with a waveguide, and on-
line compensation is essential. Since the interface echo between the rod and the pipe
wall may be suppressed after the welding operation, a notch is machined at the end
of the rod to create a reected echo, which can be used for on-line compensation.
Finally, the implementation of waveguides is considered for both pulse-echo and
xii
pitch-catch modes. In the pitch-catch mode, the backwall echo and the notch echo
are received by dierent transducers and signals of interst are both rst arrivals. As
a result, trailing echoes do not impede measurements and their attenuation becomes
unnecessary. In contrast, pulse-echo measurements are sensitive to trailing echoes
and the waveguides design plays an essential role in the feasibility of measurements.
However, the environment also imposes a set of constraints on waveguide dimensions
that complicates the implementation of pulse-echo measurements. Being more exi-
ble, the pitch-catch conguration is chosen for nal implementation. Experiments are
performed to verify the concept feasibility, and the accuracy of measurements with
thickness and temperature changes is also conrmed.
The major contributions of this thesis are the development of a 2D ray tracing
model to enable a better understanding of wave propagation in short tapered rods,
the understanding of the taper angle inuence on trailing echoes in short rods, and
the development and comparison of potential solutions for high temperature thickness
monitoring.
xiii
CHAPTER I
INTRODUCTION
The American Petroleum Institute has reported that almost every unplanned renery
shutdown is linked to failures of piping systems [1]. The negative eects of such
incidents are amplied by the diculty of maintenance personnel to anticipate needed
repairs. The development of on-line monitoring systems for pressure vessels and piping
systems could thus play a critical role in the minimization of both plant down time
and maintenance expenses.
Existing nondestructive inspection methods based on wave propagation are in
common usage for measuring the wall thickness of pipes and pressure vessels [2, 3].
These methods are suitable for periodic inspection, either manually or with auto-
mated scanning equipment. More recently, the trend is towards the implementation
of permanently mounted monitoring systems able to provide near real-time data [4].
The failure of piping systems is caused by corrosion and erosion that result in wall
thinning. As of today, the integrity of pipelines is veried through a time and cost
intensive inspection process. Production is stopped, the insulation is removed, and
transducers are mounted directly on pipe walls for thickness measurements [5].
Each time an inspection is performed, the productivity of the plant is reduced,
resulting in a signicant loss of money. The use of an on-line monitoring system could
minimize these losses by allowing measurements with minimal eect on production.
Additionally, it would reduce the testing personnel since the insulation would not
need to be removed, and it would give faster results. However, the average operating
temperature of pipe walls is around 500
1 2
Purely longitudinal
paths
Longitudinal Wave
Shear Wave
t
Paths with one shear
segment
Figure 4: Generation of trailing echoes in a cylindrical waveguide.
of mode conversions they are subjected to. This results in the creation of new delayed
echoes corresponding to dierent numbers of mode conversions, called trailing echoes.
The time delay between the rst arrival and the rst trailing echo is the same
as the time delay between two consecutive trailing echoes, and can be calculated
theoretically by the limiting case when equals zero [13],
t =
d
c
T
cos(90 )
d
tan(90 )
c
L
=
d
tan()
c
2
L
c
2
T
c
L
c
2
T
. (4)
Here d is the thickness of the plate, t is the time delay between two trailing echoes,
c
L
and c
T
are the respective longitudinal and shear wave velocities, and is the
complement of the refracted angle of the shear wave shown in Figure 4 and obtained
via,
= 90 sin
1
(
c
T
c
L
). (5)
8
Note that in the case when an echo arriving after the rst arrival is needed for a
thickness measurement, it can be obscured by the trailing echoes, making thickness
measurements inaccurate or impossible. In addition, the time delay between each
echo is proportional to the diameter of the rod. As a result, the smaller the diameter
of the rod, the greater the aect on thickness measurements.
2.3 Thickness Measurements with a Waveguide
Even though several methods exist for thickness measurements, they are all based on
the same idea. The time-of-ight through the inspected specimen is measured and
the wave velocity is then used to calculate the thickness. When a waveguide is used,
waves propagate through it before entering the specimen and every echo is delayed.
In the context of high temperature measurements, a temperature gradient is ex-
pected in the rods. This gradient may vary because the temperatures at its ends can
change. If there is no interface echo between the rod and the part, such as might be
the case when there is excellent coupling and no impedance mismatch, only a mode
1 measurement is possible. However, the temperature-induced velocity changes are
generally too large for this approach to have the needed accuracy. The common ap-
proach consists in requiring that an echo exists from the interface between the rod and
the specimen so that a mode 2 measurement can be made. Specically, the arrival
time of the reected echo is subtracted from the arrival time of the backwall echo to
deduce the time-of-ight through the specimen, and the assumption is made that the
velocity in the specimen is constant.
2.3.1 Attenuation of Trailing Echoes
The interface echo and the backwall echo are both followed by their trailing echoes. As
a result, the trailing echoes associated with the interface echo can hide the backwall
echo, making any thickness measurement impossible. To avoid this situation, several
methods have been proposed to attenuate the trailing echoes. These are illustrated
9
1 2
1 2
t
1
t
2
t=t
2
-t
1
Figure 5: Thickness measurements with a waveguide.
in Figure 6.
(a) (b) (c) (d)
Figure 6: Waveguide modications for reducing trailing echoes. (a) Grooved rod.
(b) Tapered Rod. (c) Clad Rod. (d) Bundle Waveguide.
Most of the solutions lie in the modication of the boundaries since trailing echoes
are generated via mode conversions at the rod boundaries. The rst one consists in
using grooved rods and was rst considered in 1959 by McSkimin [14]. This proved
to attenuate trailing echoes, but it did not suppress them completely. The use of
tapered rods was rst considered in 1990 [15]. They were shown to attenuate trailing
10
echoes, but reduced the energy transmitted to the specimen. In parallel, clad rods
were investigated and proved to provide very good attenutation of trailing echoes
[16]. They are widely used today for measurements and sometimes combined with
rod tapering for optimal results [17] [18].
Besides boundary modications, the use of bundle waveguides was considered [19].
The underlying idea is that trailing echoes are not generated in waveguides of diam-
eter smaller than a wavelength. However, tiny waveguides do not provide a sucient
signal-to-noise ratio for thickness measurements. It is thus necessary to use several of
them packed tightly together to perform measurements. This solution is potentially
interesting, but the attachment of the bundle waveguides to both the transducer and
part is problematic, which reduces its practicality.
Table 1: Comparison of the dierent solutions proposed for trailing echoes attenu-
ation.
Solution Strengths Weaknesses
Grooved Rod Inexpensive and easy to ma-
chine
Limited signal improvements
Tapered Rod Good results in any material Reduction in signal-to-noise ra-
tio
Clad Rod Best signal improvement Expensive and material depen-
dant
Bundle Waveguide No trailing echoes Dicult attachment
11
CHAPTER III
STUDY OF WAVE PROPAGATION IN RODS
This chapter proposes an approach for the study of high temperature measurements
and presents the tools used to accomplish it.
3.1 Overview
Though robust methods exist for ultrasonic thickness measurements, they are not
designed for continuous monitoring in high temperature environments. The survival
of the transducer and its coupling become major issues when the temperature in-
creases. As a result, methods have to be adapted to support the high temperature
environment.
The transducer can be protected from high-temperatures in dierent ways, but
the solution chosen has to allow for reliable coupling. The rst option consists in
using high temperature transducers, which have a resistance to temperature up to
482
C [20]. Another solution consists in cooling the transducer with a cooling system
[21]. Finally, the background section suggested the use of waveguides as a solution to
keep transducers at a reasonable temperature.
Experiments performed by Mechanical Integrity, Inc. showed that a wide variety
of coupling methods did not last longer than a few months at high temperatures [22].
Lowering the operational temperature of the couplant is thus necessary to enable
long-term measurements. Given the thermal inertia of pipes, it is not possible to cool
their walls locally and the only solution consists in transferring the coupling area to a
cooler environment. In this perspective, the use of waveguides appears to be the best
solution since they can be welded on pipe walls. Then, the transducer and its coupling
can be transferred through the insulation to a lower temperature environment (i.e.
12
the opposite end of the buer rod) and both issues are addressed simultaneously.
As a result, the decision was made to study the implementation of welded wave-
guides as a solution to the problem of high temperature measurements. The issue then
becomes the management of trailing echoes, which can interfere with echoes of inter-
est and make thickness measurements impossible, as mentioned in the background
chapter.
It is thus necessary to address the generation of trailing echoes and their man-
agement to allow thickness measurements with buer rods. In this perspective, both
experiments and simulations are performed.
3.2 Experimental Methods
For experiments, plate sections are used instead of pipes, which are a good approxima-
tion since the radius of pipes is much bigger than the dimensions of the measurement
area. Three steel plates of respective thickness 9.54 mm (3/8), 12.7 mm (1/2) and
19.05 mm (3/4) were used for the measurements. The waveguides were machined
from steel rods of diameter 9.54 mm (3/8), 12.7 mm (1/2) and 19.05 mm (3/4).
Rods of dierent lengths and taper angle were machined for experiments in this study.
Testing itself was performed with 10 MHz Panametrics transducers, model V552-
SM and diameter 9.54mm (3/8), excited by a Panametrics pulser-receiver, model
5072 PR. Transducers were coupled to the tested parts with Panametrics couplant A
or C. Waveforms were acquired by a Tektronix TDS5034B digital oscilloscope with a
sampling frequency of either 125 MHz or 250MHz, and were saved for further analysis.
The pulser-receiver has two measurement modes. In pulse-echo mode, only one
transducer is used and it serves as both transmitter and receiver. In through-
transmission mode, the transmitter and the receiver are distinct transducers. Figure
7 presents signals obtained in each of these two modes. For pulse-echo mode, a single
transducer was mounted on one end of the rod, and for through-transmission mode,
13
the two transducers were mounted on opposite ends of the rod.
0 10 20 30 40 50 60 70 80 90 100
-0.2
-0.1
0
0.1
0.2
0.3
Time (s)
0 10 20 30 40 50 60 70 80 90 100
-0.2
-0.1
0
0.1
0.2
0.3
Time (s)
(a)
(b)
Figure 7: Waveform obtained with a cylindrical waveguide of length 152.4 mm
(6) and diameter 9.54 mm (3/8) (a) Pulse-Echo Mode. (b) Through-Transmission
Mode.
For better coupling, threaded holders were fabricated to hold the transducers to
rods of diameter 9.54 mm (3/8) and 12.7 mm (1/2). Plastic holders were also
fabricated to hold rods at a xed angle to the plate surface.
3.3 Simulation Tools
Simulation tools enable a deeper exploration of wave propagation phenomena while
providing more repeatable results than experiments. In addition they reduce the
cost of experiments by enabling pre-selection of relevant designs. However, it is
important to keep in mind that simulations are approximations and they need to be
conrmed by experiments. The simulations reported here are limited to 2D models
for simplication and reduction of the computational time. This approximation has
to be kept in mind when interpreting the simulation results.
14
3.3.1 Finite Dierence Model
The rst model considered for simulations is based on the nite dierence method.
A commercial software package, Wave2000 [23], was purchased. To permit use of this
software, the geometry of the waveguide was represented using a two-dimensional
shape. An example is shown in Figure 8, which is the radial-axial cross-section of
a tapered waveguide. The source for all Wave2000 simulations was modeled as a
10 MHz sinusoid over two cycles with the built-in source modeling tool. Although
it is slightly dierent from a real tone burst, the main goal of this tool is not to
produce a perfect match with the experimental waveforms, but rather to understand
the propagation in waveguides to better design them.
(a)
(b)
Figure 8: (a) Geometry used for nite dierence simulation of a 0.5 degree tapered
rod. (b) Transducer source function.
Once the source and the geometry are dened, the software performs nite dier-
ence calculations from the transmitter to the receiver to generate the output wave-
form. The waveform obtained after propagation in a cylindrical waveguide in a
through-transmission conguration is provided in Figure 9 as an example. The echo
amplitudes do not match the experimental waveform, but the model does capture
the trailing echoes. Also the eect of a nite grid appears to introduce noise in the
15
simulated waveform, as is evident in Figure 9(a).
25 30 35 40 45 50 55 60
-1
0
1
Time(s)
25 30 35 40 45 50 55 60
-1
-0.5
0
0.5
1
Time (s)
(a)
(b)
Figure 9: Waveform obtained for through-transmission propagation in a cylindrical
waveguide of length 152.4 mm (6) and diameter 9.53 mm (3/8). (a) Simulated with
Wave2000. (b) Experimental.
3.3.2 Ray Tracing Model
Ray tracing is a geometrical method based on paths of rays emanating from the source
[24]. Specically, the transmitter is modeled as an array of point sources distributed
uniformly across its surface, with each point source emitting energy distributed uni-
formly in angle. Then the path of each ray is computed between the source and
the receiver, and the output waveform is generated by adding the contribution of all
the rays reaching the receiver. In this approach, there is grid for computation; the
points of reection on the boundaries are determined mathematically without any
grid eect at the edges. Thus, there is no numerical noise as was the case for nite
dierence simulations. The mode conversions on boundaries are taken into account
to reproduce the trailing echoes.
The same 2D cross section as the one used for nite dierences was used for the 2D
representation of the rod. A 3D correction was made by weighting the contribution
16
from each point source by the annular area corresponding to its radius.
Each point source radiates isotropically in all directions within the 2D cross-
sectional plane, and rays were generated every 0.05 degree for the subsequent cal-
culations. Note that rod curvature eects are not taken into account in this model
because point sources only emit in the 2D cross-section considered.
Since all the rays emit simultaneously in time at the source, and because material
dispersion and attenuation are neglected, the source waveform can be scaled, shifted
and added at the receiver based upon the computed arrival time and amplitude of
each ray. The arrival time is determined from the path calculation by using the
appropriate wave velocity for each propagating mode. The amplitude is updated
at each mode conversion using theoretical coecients for plane wave reection at a
planar boundary [25]. Once this information is gathered for each ray, the received
waveform can be generated. For the model to be consistent with experiments, the
source waveform was chosen to be a 10 MHz sinusoid multiplied by a Hanning window
over three cycles, and is illustrated in Figure 10.
0 0.05 0.1 0.15 0.2 0.25 0.3
-1
-0.5
0
0.5
1
Time (s)
Figure 10: Wavelet used to model the source emission in the ray tracing model.
As an example, the evolution of a ray from its emission to its arrival at the receiver
is detailed here and illustrated in Figure 11. Ray 1 is subjected to a mode conversion
and splits into rays 2 and 3. Then, ray 2 directly reaches the receiver while ray 3
17
is subjected to another mode conversion, resulting in rays 4 and 5, both hitting the
receiver.
1
3
4
2
Longitudinal Wave
Shear Wave
5
Transmitter
Receiver
Figure 11: Propagation of a ray in a straight rod
As a result, ray 1 produces three arrivals at the receiver. However, ray 5 reaches
the receiver as a shear wave and is thus discarded; the assumption is made that the
receiver is only sensitive to longitudinal waves. To create the waveform at the receiver,
it is then needed to determine the arrival times and amplitudes of rays 2 and 4.
According to Snells law, rays 1 and 4 propagate with the same angle , ray 2
propagates at an angle of ( ), ray 3 propagates at an angle of , and ray 5
propagates at an angle of ( ). All angles are measured counterclockwise from a
vertical line, and and are related by,
= sin
1
(
c
T
c
L
sin), (6)
where c
L
and c
T
are the longitudinal and shear wave respective velocities. The arrival
18
times are calculated from the ray paths. The rst arrival corresponds to the combi-
nation of ray 1 and ray 2. The propagation is longitudinal along the entire path, and
the arrival time is given by,
t
2
=
l
1
+ l
2
c
L
, (7)
where l
1
is the length of ray 1, l
2
is the length of ray 2, and c
L
is the longitudinal
speed of waves in steel. Concerning its amplitude, ray 2 results from the reection of
ray 1 on the boundary. The reection coecient is,
A
2
= A
1
sin(2)sin(2)
2
cos
2
(2)
sin(2)sin(2) +
2
cos
2
(2)
(8)
where is Poissons ratio, A
1
is the amplitude of ray 1, and A
2
the amplitude of
ray 2. The second arrival corresponds to the combination of rays 1, 3 and 4. The
propagation is longitudinal for rays 1 and 4, and shear for ray 3. Consequently, the
arrival time of ray 4 is
t
4
=
l
1
+ l
4
c
L
+
l
3
c
T
, (9)
where l
1
is the length of ray 1, l
3
is the length of ray 3, l
4
is the length of ray 4, c
L
is the longitudinal wave velocity in steel, and c
T
is the shear wave velocity in steel.
Finally, the amplitude of ray 4 is determined by considering the amplitude partition
resulting from the two mode conversions along its path. The relation between the
amplitude of ray 1 and the amplitude of ray 4 is given by,
A
4
= A
1
2
2
sin(4)sin(2)cos(2)
(sin(2)sin(2) +
2
cos
2
(2))
2
, (10)
where is Poissons ratio, A
1
is the amplitude of ray 1, and A
4
the amplitude of ray 4.
The contributions of ray 1 to the nal waveform can then be plotted by adding
19
scaled wavelets at the corresponding arrival times. Figure 12 presents the contribution
of ray 1, while Figure 13 gives the nal waveform obtained via the entire ray tracing
simulation. Both gures are for the through transmission in a cylindrical waveguide
of length 152.4 mm (6) and diameter 9.54 mm (3/8).
As it was the case for the Wave 2000 simulations, the amplitude partition between
the initial and trailing echoes for ray tracing does not match the experiments. How-
ever, ray tracing does not have simulation noise resulting from the nite dierence
grid.
Table 2: Numerical values associated with the propagation of ray 1.
Rod Length Rod Diameter A
1
A
2
A
4
2
56.54
D
N
Figure 16: Model of transducer emission.
The velocity of waves depends upon the material, which has to be the same as the
pipe walls to enable the welding operation. However, the diameter and the frequency
of the transducer can be tuned to change the angle of the beam. Their impacts are
similar and thus there is no need to study both. It is thus decided to limit the study
to the inuence of the transducer diameter on trailing echoes.
The formula being valid in the far eld only, the average near eld distance needs
to be estimated. It is 40 mm for the transducer used here, which is not negligible in
the rods considered for experiments. In the near eld, waves travel in every direction
and many mode conversions can occur, thus contributing to trailing echo generation.
Theoretical studies based on the formula 11 are not sucient to study the im-
pact of the transducer diameter on trailing echoes. Besides, the limited number of
transducers in the diameter range studied does not allow precise experimental stud-
ies. As a result, it is decided to study the impact of the transducer diameter on
trailing echoes through simulations. Experimental measurements are then made for
verication purposes.
27
Through-transmission waveforms were simulated through ray tracing for cylindri-
cal waveguides of length 152.4 mm (6) and diameter 9.54 mm (3/8), with trans-
ducer diameters ranging between 0.5 mm and 9.5 mm. The distance between two
point sources used to model the transducers is set constant and equal to 0.2 mm.
Diameter variations were then implemented by adding or removing point sources. To
compensate for the variation in transmitted energy from one transducer to the other,
waveforms were normalized before comparison. In a through-transmission congura-
tion with the waveguide dimensions chosen here, most of the energy is split between
the rst arrival and the rst trailing echo. Thus, the attenuation of trailing echoes can
reasonably be quantied through the amplitude ratio between the rst trailing echo
and the rst arrival. The other trailing echoes experimentally proved to be smaller
than the rst one for every transducer diameter.
The results of simulations are presented in Figures 17 and 18. It appears that there
exists an optimal transmitter diameter such that the cleanest waveform is obtained
for through transmission in the waveguide considered. However, the improvement
brought by the choice of the optimal diameter is not that signicant and denitely
not sucient to signicantly improve measurements with waveguides. Experiments
were made with three transducers of diameters 3.175 mm (0.125), 6.4 mm (0.25),
and 9.5 mm (0.375), but the variation of the attenuation was too small to be detected.
It is thus concluded that the inuence of the transducer diameter on trailing echoes
is negligible.
4.3 Inuence of Boundary Modications on Trailing Echoes
The background section discusses several solutions for the attenuation of trailing
echoes based on waveguide design. Although giving the best results, clad rods are
not considered here because they are judged to be incompatible with the welding
operation. In addition, their manufacturing is expensive, which reduces their appeal
28
25 26 27 28 29 30 31
-1
0
1
Time (s)
0.5 mm
-1
0
1
9.5 mm
-1
0
1
8.5 mm
-1
0
1
7.5 mm
-1
0
1
6.5 mm
-1
0
1
5.5 mm
-1
0
1
4.5 mm
-1
0
1
3.5 mm
-1
0
1
2.5 mm
-1
0
1
1.5 mm
Figure 17: Through-transmission waveforms in a cylindrical rod of length 152.4 mm
(6) and diameter 9.54 mm (3/8) for several transducer diameters.
29
0 1 2 3 4 5 6 7 8 9 10
-16
-15.5
-15
-14.5
-14
-13.5
-13
-12.5
-12
-11.5
-11
Transducer Diameter (mm)
A
t
t
e
n
u
a
t
i
o
n
R
a
t
i
o
(
d
B
)
Figure 18: Relative attenuation of the trailing echoes with transducers diameter.
for large scale implementation. For the specic needs of the application considered
here, tapered rods provide a reasonable compromise. They can be machined from
any material, are inexpensive, and give reasonably good results for the attenuation
of trailing echoes.
According to the literature, cleaner signals are obtained when the taper angle
increases [15]. However, the link between the taper angle and the attenuation of the
trailing echoes has not been explicitly shown. In order to get the best design available,
it was thus decided to study the inuence of the taper angle on trailing echoes.
Preliminary measurements were performed by acquiring the signals resulting from
pulse-echo and through-transmission measurements for a 0.75 degree taper angle.
Figure 19 presents the corresponding results. As expected by reciprocity [27], the two
through-transmission signals are almost identical. In contrast, the two pulse-echo
signals are totally dierent. The waveform obtained with the source at the small end
is relatively clean, although lower in amplitude, while the signal obtained with the
source at the large end shows improvements on the rst trailing echo only. The other
trailing echoes are amplied. According to table 3, two trailing echoes fall in the
range of arrival times of the backwall echo for the rod of diameter 9.54 mm (3/8)
chosen. Additionally, the trailing echoes tend to have earlier arrival times when the
30
taper angle increases. Thickness measurements are not possible in this conguration
and they are not considered.
20 30 40 50 60 70 80 90 100
-2
0
2
Time (s)
20 30 40 50 60 70 80 90 100
-2
0
2
Time (s)
20 30 40 50 60 70 80 90 100
-1
0
1
Time (s)
20 30 40 50 60 70 80 90 100
-1
0
1
Time (s)
(a)
(b)
(d)
(c)
Figure 19: Waveforms recorded in dierent measurement congurations on a 0.75
tapered rod of length 152.4 mm (6) and larger diameter 9.54 mm (3/8). (a)
Through-transmission conguration with source at the small end. (b) Through-
transmission conguration with source at the large end. (c) Pulse-echo conguration
with source at the small end. (d) Pulse-echo conguration with source at the large
end.
Finally, the trends are the same for pulse-echo measurements from the small
end and for through-transmission. As a result, the study is carried out in through-
transmission mode for simplication of the simulations.
4.3.1 Experimental Observations
In order to capture the major trends linked with the introduction of a taper angle,
through-transmission signals were acquired for ve rods of respective taper angles
31
of 0
, 0.25
, 0.5
, 0.75
, and 1
taper angle.
(c) 0.50
taper angle.
As expected, the relative amplitude of the rst trailing echo decreases when a
32
taper angle is introduced. Once again, the attenuation was quantied through the
amplitude ratio between the rst trailing echo and the rst arrival. The resulting
data are presented in Figure 21.
The limited data shown here indicates that for larger taper angles, the increase
of the taper angle does not necessarily reduce the relative amplitude of the trailing
echoes. A close look at the waveforms shows that the trailing echoes tend to split and
spread out in time as the taper angle increases, which contributes to noise or clutter
behind the rst arrival. As a result, after a certain amount of taper, any increase
in the taper angle results in a deterioration of the waveform as the predominant
eect becomes echo spreading without a reduction in amplitude. For the specic rods
considered here, the optimal angle marking the limit between signal improvement and
signal deterioration is between 0.25
and 0.75
.
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
-25
-20
-15
-10
-5
0
Taper Angle (degrees)
A
m
p
l
i
t
u
d
e
R
a
t
i
o
(
d
B
)
Figure 21: Amplitude of the rst trailing echo relative to the rst arrival versus
taper angle. Each data point represents a separate waveform acquisition.
Experimental results are not precise enough for a ne study of the impact of the
taper angle. In addition, the cost of machining rods can be avoided by going to sim-
ulations. The 2D models introduced in chapter 3 are thus used to obtain a better
understanding of the taper angle inuence on trailing echoes and attempt to deter-
mine the bounds for the optimal angle. Note that models always have imperfections,
especially here where a 2D model is used to represent a 3D environment. Results
33
thus need to be analyzed keeping in mind the approximate nature of the simulations.
4.3.2 Results of Simulations
Simulations were performed for 21 taper angles from 0
to 1
with an increment of
0.05
. The attenuation of trailing echoes was quantied by the amplitude ratio be-
tween the rst trailing echo and the rst arrival. Results obtained for nite dierences
are presented in Figures 22 and 23.
-1
0
1
-1
0
1
-1
0
1
-1
0
1
25 26 27 28 29 30
-1
0
1
Time (s)
(a)
(b)
(c)
(e)
(d)
Figure 22: Waveforms resulting from nite dierence simulations for dierent taper
angles on a waveguide of length 152.4 mm (6) and diameter 9.54 mm (3/8) (a)
straight rod. (b) 0.25
taper angle.
Though attenuation of trailing echoes is observed when the taper angle increases,
this model has several imperfections. First, it does not reect the amplitude decrease
of the rst arrival when the taper angle increases. At the same time, it underestimates
34
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
-35
-30
-25
-20
-15
-10
-5
Taper Angle (degrees)
A
m
p
l
i
t
u
d
e
R
a
t
i
o
(
d
B
)
Figure 23: Impact of the taper angle on the trailing echo attenuation with the nite
dierence model.
the amplitude of trailing echoes. Finally, the nite dierence grid creates simulation
noise that precludes trailing echo analysis for taper angles larger than about 0.4
degree. Any determination of an optimal angle is thus impossible with this model.
However, the results obtained with nite dierences can be compared to the results
of ray tracing to determine if imperfections come from the lack of 3D considerations.
Ray tracing results are presented in Figures 24, 25 and 26. Figure 25 shows the
same signals as Figure 24 but with dierent vertical scales for each signal to better see
the trailing echoes. As expected, this model allows the exploration of angles larger
than 0.4 degrees since there is no simulation noise. In addition, it captures the echoes
shifting and spreading that could not be observed with nite dierences. However,
the respective amplitudes of each echo are still not in agreement with experiments.
To better understand the imperfection of models and thus learn about wave prop-
agation in tapered rods, all data are presented in Figure 27. Interestingly, the results
obtained for ray tracing and nite dierence are in reasonably close agreement for
taper angles smaller than 0.4 degree. After this amount of taper, results are not
comparable anymore due to the simulation noise introduced in the nite dierence
model. This observation indicates that the misrepresentation of amplitudes in simu-
lated waveforms is likely linked with the lack of 3D considerations.
35
-2
0
2
-2
0
2
25 26 27 28 29 30 31
-2
0
2
Time (s)
-2
0
2
-2
0
2
(a)
(b)
(c)
(e)
(d)
Figure 24: Waveforms resulting from ray tracing simulations for dierent taper
angles on a waveguide of length 152.4 mm (6) and diameter 9.54mm (3/8) (a)
straight rod. (b) 0.25
taper angle.
36
-0.5
0.5
0
-1.1
0
1.1
-0.2
0
0.2
-0.3
0.3
0
25 26 27 28 29 30 31
-0.1
0
0.1
Time (s)
(a)
(b)
(c)
(e)
(d)
Figure 25: Same as Figure 24 after zooming on the trailing echoes (a) straight rod.
(b) 0.25
taper
angle.
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
-45
-40
-35
-30
-25
-20
-15
-10
-5
Taper Angle (degrees)
A
m
p
l
i
t
u
d
e
R
a
t
i
o
(
d
B
)
Figure 26: Impact of the taper angle on the trailing echo attenuation with the ray
tracing model.
37
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
-50
-40
-30
-20
-10
0
10
Taper Angle (degrees)
A
m
p
l
i
t
u
d
e
R
a
t
i
o
(
d
B
)
Experiments
Ray Tracing
Finite Difference
Figure 27: Comparison of trailing echo amplitude ratios from experiments and both
simulations.
The representation of rods by their cross-sections limits the study to radial con-
tributions. However, in reality rays are emitted in every direction and o-axis rays
should be represented. These rays undergo more reections as they spiral down
the rod and they are subjected to curvature eects. As a result, their amplitudes are
primarily transferred to trailing echoes. The consideration of these rays in the ray
tracing model would be expected to give a more realistic amplitude partition between
the rst arrival and the trailing echoes.
To conclude, better results could be obtained with 3D models, but at the ex-
pense of a large increase in the computational time. The 2D models presented here,
particularly the ray tracing model, already allow a good understanding of the wave
propagation in tapered rods and are sucient to understand the impact of the taper
angle on trailing echoes.
4.3.3 Qualitative Analysis of the Results
Though precise quantication of results is not possible with the data reported here,
the simulations gave a better understanding of the taper angle impact on trailing
echoes. Specically, the modeling of waves through propagating rays greatly helps in
studying the evolution of waveforms. In this context, each echo is the result of the
38
constructive interference of multiple ray contributions. The study of how the echoes
evolve can then be understood by studying their components. The n
th
trailing echo
is subsequently referred as the trailing echo of order n to simplify the explanation.
In a straight rod, trailing echoes are composed of ray contributions with similar
arrival times that constructively interfere. When a taper angle is introduced, the
rod diameter is not constant and its variations result in a wider range of arrival
times for each contribution, which reduces the constructive interference eect and
thus reduces the amplitude of the trailing echoes. Figure 28 illustrates the inuence
of the diameter variations on the arrival times. The rst arrival is less impacted by
diameter variations since most of its contributions are provided by rays that never
hit boundaries.
Same arrival
time
Different
arrival times
Figure 28: Change in arrival times of trailing echo rays after introduction of a taper
angle.
Another eect of the taper angle is a variation in the number of mode conversions.
39
When traveling from the large end to the small end of the rod, the number of mode
conversions increases with respect to the straight rod case, while it decreases when
traveling in the opposite direction. Consequently, the amplitude of each contribution
as well as the energy partition between the echoes is modied. In addition, the
incidence angle on boundaries is modied by the taper angle, which also has an
impact on the energy partition.
The taper angle also interacts with the transducer beam angle. When traveling
from the small end to the large end, the beam tends to focus, while it spreads in the
opposite direction. As a result, when the transducer is at the larger end, waves are
subjected to more mode conversions on their round-trip than in the case when the
transducer is at the smaller end, and how the energy is partitioned is also aected.
To summarize, the introduction of a taper angle results in two primary eects.
First, the arrival times of the ray contributions to a single trailing echo split and
spread out, which reduces the amplitude of that echo. Second, the energy partition is
modied. Energy tends to be transferred to higher order trailing echoes when waves
propagate from the large end to the small end, while the opposite eect occurs in the
opposite direction. At the same time, more energy is transferred towards the trailing
echoes if the transducer is at the large end.
Signal improvements mainly result from the widening of the range of arrival times
for the ray contributions of each trailing echo. This eect is not as strong on every
trailing echo. The higher the trailing echo order, the smaller the impact of the taper
angle. In long rods, the impact is always important even on high order trailing echoes,
but in shorter rods, the widening eect only applies to low order trailing echoes.
This behaviour can be explained with geometrical considerations. Figure 29 presents
the paths of contributions to the two rst trailing echoes in a short tapered rod.
When considering the ray contributions to the rst trailing echo, there is a signicant
number of paths that are very dierent, and the introduction of the taper angle has
40
a real impact on the echo amplitude. Concerning the contributions to the second
trailing echo, the number of possible paths is reduced and the diameter variations
have similar impacts on multiple contributions. As a result, several contributions still
have very similar arrival times and they continue to interfere constructively.
1 3 4 2
Longitudinal Wave
Shear Wave
(a) (b)
Figure 29: Variability of the widening eect with the trailing echo order. (a) First
order trailing echo order. (b) Second order trailing echo.
With these considerations, it becomes possible to explain the lack of improvement
on the pulse-echo signal of Figure 19(d). The transducer was coupled to the large
end of a short tapered rod of angle 0.75
C [28].
47
Table 4: Error introduced by temperature changes in the measured thickness of a
12.7 mm (1/2) plate
Method T(
C) Thickness error (%)
contact transducer 10 0.11
contact transducer 100 1.08
contact transducer 500 5.63
waveguide 10 1.43
waveguide 100 14.04
waveguide 500 73.19
For the purpose of monitoring thickness changes caused by corrosion and erosion,
an error of 1% is acceptable. Thus, there is no need for calibration with contact
transducers for temperature changes up to 100
C. Then,
it was placed at room temperature, and measurements were performed as it cooled
down. The value of the highest temperature for which measurements were performed
was kept low, because the ecacy of waveguides in reducing the temperature still
needs to be veried and cooling solutions such as the addition of heat sinks might be
needed to avoid transducer damage.
Figure 41 compares the results of measurements with and without pulse-echo
57
12 14 16 18 20 22 24 26
5
6
7
8
9
Thickness (mm)
T
i
m
e
o
f
f
l
i
g
h
t
(
s
)
Linear Regression curve
: Experimental data
Figure 40: Generation of trailing echoes in a cylindrical waveguide.
compensation at room temperature (i.e. mode 2 vs mode 1). As expected, accurate
results are obtained through both methods since the initial calibration was performed
at room temperature.
16
18
20
22
24
26
T
h
i
c
k
n
e
s
s
(
m
m
)
5 10 15 20 25
12
14
16
Measurement Number
Actual Thickness
Measured Thickness without Online Calibration
Measured Thickness with Online Calibration
Actual Thickness
Thickness without Pulse-Echo Compensation
Thickness with Pulse-Echo Compensation
Figure 41: Thickness measurements at room temperature.
Figure 42 presents the evolution of the measured thickness at elevated tempera-
tures. With pulse-echo compensation, the measurement error is small. Some errors
58
are caused by the hand coupling of the waveguides. However, as expected, a larger
error is introduced without the pulse-echo compensation, which demonstrates the
usefulness and ecacy of the proposed method.
16
18
20
22
24
26
T
h
i
c
k
n
e
s
s
(
m
m
)
40C 35.9C 28.2C
32.1C
26 28 30 32 34 36 38
12
14
16
Measurement Number
Actual Thickness
Measured Thickness without Online Calibration
Measured Thickness with Online Calibration
42C
Actual Thickness
Thickness without Pulse-Echo Compensation
Thickness with Pulse-Echo Compensation
Figure 42: Results of thickness measurements at elevated temperatures. The initial
calibration was performed at 22.5
C.
59
CHAPTER VI
CONCLUSION
6.1 Summary and Conclusions
The objective of this thesis has been the development and verication of methods for
on-line thickness monitoring of piping systems and pressure vessels at high tempera-
ture. The decision to use welded waveguides had already been made based upon prior
work where multiple coupling methods were investigated [22], and thus the focus of
this thesis has been the determination of the precise conguration to employ.
A study was performed to understand both pulse-echo and through-transmission
wave propagation in straight and tapered buer rods. The primary issue was gener-
ation of trailing echoes in such rods and their impact on thickness measurements. A
prior literature search indicated that increasing the taper angle would reduce trailing
echoes in long buer rods. Experiments were performed, nite dierence simulations
were made, and a ray tracing model was developed to enable understanding of wave
propagation in tapered buer rods. It was found that the results for long buer rods
do not apply to shorter rods. Unlike reported results for longer rods, signals obtained
with the transducer mounted on the larger end of a short tapered rod are signicantly
noisier in terms of trailing echoes than signals obtained with the transducer mounted
on the other end. Although a modest taper angle does reduce the trailing echoes,
further increases in taper angle do not provide additional improvement because the
shorter length of the rod limits the number of mode conversions that can occur and
thus limits the variability in the arrival times of dierent beam paths. The conclusion
reached was that reasonable combinations of the diameter, length and taper angle of
the buer rod would not be adequate to reduce trailing echoes enough to permit
60
pulse-echo thickness measurements of pipe walls at high temperature.
A pitch-catch conguration was then considered using two adjacent buer rods
mounted at a small angle to the exposed surface. For this conguration, the backwall
echo is the rst arrival and is thus not aected by trailing echoes. The problem is
that temperature-induced velocity changes cause time shifts in this arrival that could
be mistaken for thickness changes. The proposed solution is to machine a notch
in the rods close to their attachment points on the surface and to simultaneously
record the pulse-echo signal. The arrival time of the echo reected from the notch
can then be used to compensate the error introduced by waveguides in the pitch-
catch arrival time, thereby reducing the sensitivity of the thickness measurements to
velocity perturbations. Measurements were performed to verify the ecacy of this
proposed solution.
6.2 Future Work
It is critical that the proposed pitch-catch method be veried using welded rods at
high-temperatures, and future work must include such experiments. The verication
should include the accuracy and robustness of the ultrasonic method as well as thermal
considerations. If the temperature at the transducer is too high to ensure survivability,
then the waveguide will need to be cooled. The preferred method would be installation
of passive heat sinks on the rods.
Once point thickness measurements are implemented, the ability to use the same
rods and transducers to generate guided waves in the pipe walls should be investigated.
The objective would be to enable large area monitoring between locations of point
measurements.
Finally, further investigation of a pulse-echo solution should be considered. If
the rod is made longer and the taper angle is increased, it may be possible to make
thickness measurements. The primary issue would be whether the signal-to-noise
61
ratio of the backwall echo is suciently large.
62
REFERENCES
[1] American Petroleum Institute. Technology Vision 2020: A Technology Vision
for the U.S Petroleum Rening Industry, 1999.
[2] D. E. Bray and R. K. Stanley. Nondestructive Evaluation : A Tool in Design,
Manufacturing, and Service. CRC Press, 1997.
[3] J. Krautkramer. Ultrasonic Testing of Materials. Springer, 1990.
[4] C. R. Farrar and K. Worden. An introduction to structural health monitoring.
Phil. Trans. R. Soc. A, 365:303315, 2006.
[5] M. Fields. Inservice pressure vessel inspections. Proceedings of SPIE, 3398:6674,
1998.
[6] K. A. Fowler, G. M. Elfbaum, K. A. Smith, and T. J. Nelligan. Theory and
application of precision ultrasonic thickness gauging. Insight - Non-Destructive
Testing and Condition Monitoring, 38(8):582587, 1996.
[7] Panametrics. Panametrics-NDT Ultrasonic Transducers, 2008. pp. 8-9.
[8] A. E. H. Love. The Mathematical Theory of Elasticity. Cambridge University
Press, London, fourth edition, 1934. p289.
[9] H. J. McSkimin. Propagation of longitudinal waves and shear waves in cylindrical
rods at high frequencies. J. Acoust. Soc. Am., 28(3):484493, 1956.
[10] W. P. Mason and H. J. McSkimin. Attenuation and scattering of high-frequency
sound waves in metals and glasses. J. Acoust. Soc. Am., 19(3):464473, 1947.
[11] E. P. Papadakis. Buer-rod system for ultrasonic attenuation measurements. J.
Acoust. Soc. Am., 44(5):14371441, 1968.
[12] V. S. K. Prasad. Viscosity measurements of melts at high-temperatures using
ultrasonic guided waves. J. Material Processing Technology, 207(1-3):315320,
2008.
[13] R. N. Thurston. Elastic waves in rods and clad rods. J. Acoust. Soc. Am.,
64(1):137, 1978.
[14] H. J. McSkimin. Measurement of ultrasonic wave velocities and elastic moduli for
small solid specimens at high-temperatures. J. Acoust. Soc. Am., 31(3):287295,
1959.
63
[15] C.-K. Jen, L. Piche, and J. F. Bussiere. Long isotropic buer rods. J. Acoust.
Soc. Am., 88(1):2325, 1990.
[16] C.-K. Jen, J.-G. Legoux, and L. Parent. Experimental evaluation of clad metal-
lic buer rods for high temperature ultrasonic measurements. NDT&E Interna-
tional, 33(3):145153, 2000.
[17] I. Ihara, D. Burhan, and H. Aso. Ultrasonic in-line sensors for inclusion detection
in liquid metals. IEEE Ultrasonics Symposium, 1:811814, 2002.
[18] Y. Ono, J.-F. Moisan, and C.-K. Jen. Ultrasonic techniques for imaging and
measurements in molten aluminum. IEEE Trans. Ultrason., Ferroelect. Freq.
Contr., 50(12):17111721, 2003.
[19] L. C. Lynnworth, Y. Liu, and J. A. Umina. Extensional bundle waveguide
techniques for measuring ow of hot uids. IEEE Trans. Ultrason., Ferroelect.
Freq. Contr., 52(4):538543, 2005.
[20] Piezotechnologies. High temperature transducers. http://www.
piezotechnologies.com/temp900.htm (Last visit 12 November 2008).
[21] S. P. Kelly, I. Atkinson, G. Chris, and K. J. Kirk. On-line ultrasonic inspection
at elevated temperatures. IEEE Ultrason. Symp., pages 904908, 2007.
[22] Mechanical Integrity, Inc. Private communication with Donald McNicol.
[23] Cyberlogic,Inc. Wave2000 Software for Computational Ultrasonics. www.
cyberlogic.org/wave2000.html (Last visit 12 november 2008).
[24] J. A. Johnson, N. M. Carlson, and L. A. Lott. Ultrasonic wave propagation in
temperature gradients. J. NDT&E, 6(3):147157, 1987.
[25] J. D. Achenbach. Wave Propagation in Elastic Solids. North-Holland, New York,
1999.
[26] Panametrics. Ultrasonic Transducers Technical Notes, 2006. p 43.
[27] B. A. Auld. Acoustic Fields and Waves in Solids. Robert E. Krieger Publishing
Company, Inc., 2nd edition, 1990.
[28] G. Mott. Temperature dependence of ultrasonic parameters. Review of progress
in Quantitative Nondestructive Evaluation, 3B:11371148, 1984.
[29] C. E. Jaske, J. A. Beavers, and N. G. Thompson. Improving plant reliability
through corrosion monitoring. In Proceedings of the Fourth International Con-
ference on Process Plant Reliability, pages 318, Houston, TX, 1995.
64