Login| Sign Up| Help| Contact|

Patent Searching and Data


Title:
CONTINUOUS WAVE ULTRASOUND FOR ANALYSIS OF A SURFACE
Document Type and Number:
WIPO Patent Application WO/2017/109490
Kind Code:
A1
Abstract:
A continuous wave ultrasound analysis device for analysing a surface of a component includes an input to receive ultrasound data relating to continuous wave, swept frequency ultrasound propagated through the component; a processing section to derive a standing wave amplitude spectrum from the received ultrasound data, and determine a property of an unconstrained layer on the surface based on the standing wave amplitude spectrum; and an output to output the determined property of the layer. A method of continuous wave ultrasound for analysing a surface is also provided.

Inventors:
MILLS ROBIN (GB)
DWYER-JOYCE ROBERT S (GB)
MARSHALL MATTHEW (GB)
Application Number:
PCT/GB2016/054024
Publication Date:
June 29, 2017
Filing Date:
December 21, 2016
Export Citation:
Click for automatic bibliography generation   Help
Assignee:
UNIV SHEFFIELD (GB)
International Classes:
G01N29/04; G01N29/12; G01N29/34; G01N29/44
Foreign References:
JPH09250919A1997-09-22
US6286370B12001-09-11
JP2009288045A2009-12-10
JP2008116406A2008-05-22
JP2007010543A2007-01-18
Attorney, Agent or Firm:
JENNINGS, Vincent (GB)
Download PDF:
Claims:
CLAIMS

1. A continuous wave ultrasound analysis device for analysing a surface of a component, comprising:

an input to receive ultrasound data relating to continuous wave, swept frequency ultrasound propagated through the component;

a processing section to:

derive a standing wave amplitude spectrum from the received ultrasound data, and

determine a property of an unconstrained layer on the surface based on the standing wave amplitude spectrum; and

an output to output the determined property of the layer.

2. The device of claim 1 , wherein the processing section is to approximate the amplitude of the ultrasound using

^o ei^(^-t)-n0-(n-l)0 ] e2nLaffnff'("-l)

where:

Ao is an amplitude of a driving signal of the ultrasound,

n is a number of passages of the wave through the component,

L is a thickness of the component,

t is time,

ø is a phase shift at the surface of the component,

ø is a phase shift at a face of the component at which the ultrasound is input, ω is an angular frequency of the input ultrasound,

a is an attenuation coefficient of component material,

R is a reflection coefficient at the surface of the component,

R' is a reflection coefficient at the face of the component at which the ultrasound is input, and

c is a speed of sound in the component.

3. The device of claim 1 or claim 2, wherein the processing section is to determine the property based on a comparison between the received ultrasound data and ultrasound data associated with a reference state of the component.

4. The device of claim 3, wherein the reference state is a state in which transmission of ultrasound at the surface of the component is negligible.

5. The device of claim 3, wherein the device is arranged to assume that the reference state is a state in which transmission of ultrasound at the surface of the component is negligible.

6. The device of any one of claims 1 to 5, wherein the processing section is to determining the property of the unconstrained layer based on a determination of a superimposed standing wave reflection coefficient.

7. The device of any one of claims 1 , or 6 to 8, wherein the processing section is to evaluate S, where S is given by

where:

L is a thickness of the component,

ø is a phase shift at the surface of the component,

ø is a phase shift at a face of the component at which the ultrasound is input, ω is an angular frequency of the input ultrasound,

a is an attenuation coefficient of component material,

R is a reflection coefficient at the surface of the component,

R' is a reflection coefficient at the face of the component at which the ultrasound is input, and

c is a speed of sound in the component.

8. The device of claim 7, wherein the property of the unconstrained layer is determined based on S. 9. The device of claim 8, wherein the property of the unconstrained layer is determined by: determining R based on S, and

determining the property of the unconstrained layer based on R. 10. The device of claim 9, wherein R is taken to be:

(1 - ζ12) cos J + ι(ζ3213) sin ( -^— J

c3 /

R =

ζί2) cos (^r^) + i(z3/z2 +z1/z3) sin (^r^)

1 1. The device of claim 9, wherein R is taken to be: '

12. The device of any one of claims 7 to 1 1 , wherein S is a root mean square superimposed standing wave reflection coefficient.

The device of any one of claims 2 to 12, wherein the property includes at least one of: a stiffness of the layer,

a thickness of the layer,

a bulk modulus of the layer,

a speed of sound in the layer,

a density of the layer,

a viscosity of the layer, or

a presence of the layer.

14. The device of claim 1 , wherein the property includes a thickness of the layer, and wherein the processing section is to determine the thickness based on a shift in a standing wave peak frequency relative to a reference state of the component.

15. The device of claim 14, wherein the determination of the thickness is based on a standing wave peak corresponding to a harmonic having wave mode, n, greater or equal to 2.

16. The device of claim 14 or claim 15, wherein the determination of the thickness uses

Lc + lct + hci (Lc + lct)

Afs = 2 (L + l + h)2 ~ 2 (L + 2

where:

nAfs is a shift of the standing wave peak corresponding to wave mode n,

L is a path length of the wave in the component,

c is a speed of sound in the component,

/ is a path length of the wave in an emitter for emitting the ultrasound,

Ct is a speed of sound in the emitter,

h is a path length in the layer, and

Ci is a speed of sound in the layer.

17. The device of any one of claims 14 to 16, wherein the ultrasound data relates to a measurement state of the component, and the determination of thickness is to include:

superimposing the standing wave amplitude spectrum of the measurement state with a standing wave amplitude spectrum of the reference state and

determining the shift in a standing wave peak frequency based on the superimposition.

18. The device of any one of claims 14 to 17, wherein the reference state is a state in which transmission of ultrasound at the surface of the component is negligible or assumed to be negligible.

19. The device of any preceding claim, wherein the ultrasound data is to relate to ultrasound transmitted into the component at a first face of the component, and to ultrasound data measured at the first face of the component. 20. The device of any preceding claim, wherein the ultrasound data is to relate to ultrasound swept initially over a first range, and subsequently over a second range, the second range being narrower than, and within, the first range.

21. A method of continuous wave ultrasound analysis for analysing a surface of a component, comprising:

receiving ultrasound data relating to continuous wave, swept frequency ultrasound propagated through the component;

deriving a standing wave amplitude spectrum from the received ultrasound data, and determining a property of an unconstrained layer on the surface based on the standing wave amplitude spectrum.

22. The method of claim 21 , wherein the deriving includes approximating the amplitude of the ultrasound using

where:

Ao is an amplitude of a driving signal of the ultrasound,

n is a number of passages of the wave through the component,

L is a thickness of the component,

t is time,

ø is a phase shift at the surface of the component,

ø is a phase shift at a face of the component at which the ultrasound is input, ω is an angular frequency of the input ultrasound,

a is an attenuation coefficient of component material,

R is a reflection coefficient at the surface of the component,

R' is a reflection coefficient at the face of the component at which the ultrasound is input, and c is a speed of sound in the component.

23. The method of claim 21 , wherein the property includes a thickness of the layer, and wherein the processing section is to determine the thickness based on a shift in a standing wave peak frequency relative to a reference state of the component.

24. A computer program that, when executed by an electronic device, causes the electronic device to operate as the device of any one of claims 1 to 20, or to perform the method of any one of claims 21 to 23.

25. A computer readable medium storing the program of claim 16, or storing instructions thereon that when executed by an electronic device cause the electronic device to operate as the device of any one of claims 1 to 20, or to perform the method of any one of claims 21 to 23. 26. A continuous wave ultrasound analysis device substantially as described herein with reference to the drawings.

27. A method of continuous wave ultrasound analysis for analysing a surface of a component as described herein with reference to the drawings.

Description:
Continuous Wave Ultrasound for Analysis of a Surface

[0001] This disclosure relates to the use of ultrasound for analysis of physical systems, such as mechanical components.

BACKGROUND

[0002] Ultrasound may be used for analysis of physical systems, such as for non-destructive testing of mechanical components. Ultrasound techniques typically derive information by measuring the time-of-flight or amplitude of an ultrasound pulse in an object.

BRIEF DESCRIPTION OF THE DRAWINGS

[0003] Examples are further described hereinafter with reference to the accompanying drawings, in which:

Figure 1 illustrates attenuation of an ultrasonic wave in a material.

Figure 2a illustrates reflection of an ultrasonic wave at a solid-air boundary.

Figure 2b illustrates reflection of an ultrasonic wave at a boundary between two materials. Figure 2c illustrates reflection of an ultrasonic wave at a thin layer between two materials. Figure 3 illustrates the superposition of a series of reflected waves to form a superimposed standing wave.

Figure 4 illustrates the relationship between the driving frequency and the amplitude of the superimposed wave. At f=nf s maximum coherence occurs and the standing wave frequency is maximum.

Figure 5 illustrates locations of a wave as it is emitted from a transducer, propagates, attenuates, and reflects through a component.

Figure 6a illustrates reflection amplitude and phase represented on a polar plot for a reference measurement in which the surface of interest is a boundary between the component and air.

Figure 6b illustrates reflection amplitude and phase represented on a polar plot for a measurement where there is a medium at the surface of interest.

Figure 7 illustrates an ideal frequency sweep to be applied as the driving frequency.

Figure 8 illustrates a driving frequency sweep transmitted into the substrate, in which the ideal signal of Figure 7 is attenuated by the real transducer response function.

Figure 9 is a schematic diagram of the response of a component to the input function of figure 8. The response consists of the superposition of all the internally reflected waves within that response window. Figure 10 is a schematic diagram of the superimposed standing wave root mean square (RMS) amplitude spectrum.

Figures 11 a and 11 b schematically illustrate changes in measured response due to the addition of a surface layer.

Figure 12 schematically illustrates the change in resonant frequency due to an increase in the wave path length.

Figure 13 schematically illustrates the components of an exemplary measurement apparatus.

Figure 14 schematically illustrates a device according to some examples.

Figure 15 illustrates a method according to some examples.

Figure 16 schematically illustrates paint spraying and instrumented sample.

Figure 17 illustrates a measured standing wave amplitude spectrum

Figure 18 illustrates acquisition response to subsequent layers of paint after drying.

Figure 19 is a Schematic diagram of an apparatus used to measure the viscosity of a range of Newtonian oils.

Figure 20 is a plot of the standing wave reflection coefficient variation with oil viscosity. The model (equations (10) and (22)) are compared with experimental data.

DETAILED DESCRIPTION

[0004] According to some examples, a novel superimposed standing wave method is provided to measure the presence and properties of a medium present on a surface. According to some examples, ultrasound transmitting/receiving equipment may be provided on a "back face" of the component, such that the component is between the ultrasound equipment and the surface of interest (the surface on which the medium is present). This may facilitate a non-invasive measurement that may be carried out during a coating or layer generation process, for example. According to some examples various properties and phenomena may be determined, such as:

• The thickness of a film or coating (such as paint, ice, or residual oil)

• The thickness of an oil film between components

• The viscosity of a layer of liquid, e.g. within a sealed vessel

[0005] According to some examples, a continuous ultrasonic wave is created inside the component. This may involve the use of a piezoelectric transducer driven with a continuous sinusoidal varying voltage, for example. The multiple wave reflections at the component boundaries interfere with the incident wave and superimpose to generate a standing wave. The amplitude and frequency of the standing wave may be recorded with time. This amplitude depends on the proportion of wave energy reflected at the boundary with the surface medium, known as the reflection coefficient. By extracting this reflection coefficient it is possible to deduce the properties of the medium at the surface. The frequency of the standing wave depends on the properties and geometry of the component and surface medium assembly. The thickness of the surface medium may be deduced by monitoring the frequency of the standing wave.

[0006] Examples may have one of more of the following features:

1 . A high measurement sensitivity may be achieved in some examples. The wave may be reflected many times from the interface/surface of interest, and so the effect of the interface (e.g. changes of the interface) on the signal may be magnified. This may allow the measurement of previously un-measureable properties (e.g. thin surface films, or viscosity in thin acoustically unmatched layers).

2. A very fast response time may be achieved in some examples. In conventional pulse-echo ultrasonics, the time between successive measurements depends on the pulse repetition rate. When continuous wave measurements are used each wave peak may be a measurement point. Accordingly, according to some examples time resolution of the technique is the wave frequency.

3. Examples may be implanted using very low cost instrumentation. The use of continuous waves avoids the need for a conventional ultrasonic pulsing and receiving system. This may allow, for example, the use of low cost signal generator and capture integrated circuits (ICs).

Ultrasonic Propagation, Attenuation, Reflection, and Phase Shift

[0007] Ultrasound is a mechanical vibration at frequencies above that of human hearing. The vibration consists of tiny mechanical deflections that propagate through a gas, liquid, or solid. As an ultrasonic wave passes through a material it will be attenuated as a function of distance, eventually decaying to zero amplitude. Figure 1 schematically shows an ultrasonic wave decaying in a material. There is a simple relationship between the initial amplitude of a wave, Ao, and the amplitude, A of the wave after it has traveled through a path length of L:

A=A 0 e La (1 )

Where a is known as the attenuation coefficient, which is a property of the material in which the ultrasound propagates.

[0008] If the wave strikes a boundary or interface, a reflection will occur, the magnitude of which will depend upon the conditions at the boundary/interface. The refection coefficient, R is defined as the amplitude of the reflected wave (AR) divided by the amplitude of the incident wave (A): R = ½ (2)

A, '

[0009] The magnitude and phase of R depends on the nature of the interface. Figure 2a schematically shows the reflection of an ultrasonic wave at a boundary with air, Figure 2b schematically shows the reflection of an ultrasonic wave at a boundary another material, and Figure 2c schematically shows the reflection of an ultrasonic wave at the boundary of a thin layer between two solid components. Each of these cases will reflect a different amount of the incident wave amplitude. If the interface consists of two materials perfectly bonded together (e.g. in the arrangement of figure 2b), then the reflection depends on the acoustic impedance mismatch between the two materials according to: R = - z ^ (3)

z 2 +z 1 '

where z? and Z2 are the acoustic impedance of the materials either side of the interface. For longitudinally polarised ultrasound the acoustic impedance is the product of the material density and wave speed:

z x = p 1 c 1 and z 2 = p 2 c 2 (4) For the case of a reflection with a solid-air boundary (e.g. figure 2a) then z∑«zi and practically all the wave is reflected back, R→1 and AR=AI.

[0010] Figure 2 also illustrates that the phase of the reflected signal may change at a boundary. If the wave is reflected at a boundary between an acoustically less dense to more dense medium, the phase changes by Φ=π, if dense to less dense then Φ=0.

Reflection Coefficient and Surface Medium Properties

A Three Layer System

[0011] For the case when the wave is reflected back from a three layer system (e.g. figure 2c), the reflection coefficient, R depends on the acoustic impedance of the three materials (Kinsler et al., Fundamentals of Acoustics, Wiley & Sons):

where ω is the angular frequency of the wave (=2nf) and subscript 3 corresponds to the thin layer)

A Trapped Thin Film

[0012] For the case when the intermediate layer between two solids is thin compared to the wavelength, equation (5a) can be simplified, and R depends on the stiffness, K, of the intermediate layer (see Schoenberg, 1980, J. Acoust. Soc. Am. 68, 1516). This complex reflection coefficient has both magnitude and phase, Φ, where both depend on the stiffness of the layer/film between the solids (see Reddyhoff et al, 2005, Proc IMechE J, 219, 387):

[0013] For a liquid layer, the stiffness is a function of its bulk modulus, B and the layer thickness, h according to:

K = - (7)

h '

And so:

RHzz+z^-jzz-z 2

h =

1-R 2 (8)

[0014] Thus, the film thickness can be determined from measurements of either R or Φ, or vice versa, if the bulk modulus is known (see Dwyer-Joyce et al, Proc Roy Soc A459, p.957). Alternatively if the film thickness is known the bulk modulus of the liquid may be determined. The bulk modulus may then be used to determine the density and speed of sound in the liquid. Shear Waves and Viscosity

[0015] Equations (1), (2), and (3) are true for either longitudinal or shear polarised ultrasound. However for shear polarised ultrasound the acoustic impedance of a solid and liquid are given by: z s = p s c s and z ( = ^ρ^ωη (9) [0016] Note that equation (9) is true for a Newtonian fluid. Relationships are available for non- Newtonian fluids that involve the relaxation time of the liquid, and such relationships could be used in place of equation (9), where necessary/appropriate. Combining (9) with (3), where subscript s refers to the solid and subscript / refers to the liquid gives:

[0017] In this way the viscosity of the liquid may be determined from the shear wave reflection at the solid liquid boundary.

Formation of a Standing Wave

[0018] An ultrasonic wave propagates through the material until it reaches a boundary, at which point the wave is then reflected along the path of the incident wave (assuming the wave is normally incident on the boundary), returning to the source of the original wave where it will be reflected again. This process of successive reflection continues until the wave decays to zero amplitude. Figure 3 shows a schematic diagram of a series of ultrasonic wave reflections. At each boundary the wave amplitude is reduced and (in the example illustrated) the phase changes; attenuation occurs inside the material, also reducing the wave amplitude. The successively reflecting waves will superimpose and set up a local standing wave within the material, also shown in Figure 3.

The Effect of the Driving Wave Frequency

[0019] The amplitude spectrum, A(f) of the superimposed standing wave will be a function of the material properties of the component through which it propagates, its geometry and the driving frequency of the source. At a particular driving frequency, f, the wavelength will be such that the superposed waves will constructively interfere and the combined wave will have maximum amplitude. This occurs at a frequency referred to herein as the standing wave superposition frequency f = f s . The value of f s is not unique of any particular system, since there will be multiple frequencies that will each produce a (local) maximum. The lowest value of f s , (i.e. at n= , where n is the mode number) is, assuming a simple geometry, when half the wavelength equals the length of the component (A=2L), and so:

[0020] If the frequency of the wave moves away from this frequency, f≠f s , the waves will be less coherent and the amplitude of the superimposed wave will reduce relative to the amplitude at f=f s . At a frequency of f=2f s the waves will again be coherent and the amplitude will be at a maximum. Figure 4 shows how the frequency of the ultrasonic wave and the amplitude of the superimposed reflections are related. The peaks of maximum amplitude occur at intervals of f s . Mathematical Description of Standing Wave Superposition

[0021] Figure 5 schematically shows numbered locations of a wave as it is emitted from the transducer 510, propagates, attenuates, and reflects through component 520. The waves reflect from surface 530 (surface 530 could alternatively be a boundary or interface, for example). Surface 530 which represents a free surface, may be a surface bearing a

constrained or unconstrained layer, may be in contact with another (solid) element, or may represent a region in which the properties (e.g. acoustic impedance) of component 520 change. In any case, the ultrasound waves are reflected from surface/boundary 530 and also reflected from the surface 540 in contact with the transducer 510. The vertical displacement of the wave at each reflection in Figure 5 is for clarity, so that subsequent passages of the wave through the material are not illustrated as overlapping the previous passages through component 520. Numbers in circles denote the progress of a particular portion of the wave through component 520, as described in more detail below.

[0022] At location 1 the wave has been emitted into the component 520 by the transducer 510. For the present example, the emitter function may be denoted by: 1 . A 0 e iMt (12)

[0023] At location 2, just before the first reflection, the wave has been attenuated according to equation (1 ), where L and c are the distance travelled and speed of sound in the component 520, respectively. Thus, as location 2, the amplitude of the wave is given by:

2. _4 0 e i6 ( ½ _t) e La (1 3) [0024] Immediately after the first reflection the amplitude is given by:

3. _4 o e i6 ( ½ _t) e La Ke - i0 (14)

[0025] where R' and Φ' are the reflection coefficient and phase change respectively at the surface/boundary 530. This corresponds with the component of the incident wave that is reflected, where the incident wave is given by equation (13). The wave then passes back through the component 520 and immediately before the reflection at the transducer 510 location (surface 540) has an amplitude of:

4. A o e iM ( -^e 2La Re - i0 (1 5)

[0026] This corresponds with the wave of equation (14) with attenuation for a passage through the component 520. The amplitude immediately after the reflection at the interface between the component 520 and transducer 510 interface is given by:

5. A o e iM ( -^e 2La Re - i0 R'e - i0 ' (16) where R' and Φ' are the reflection coefficient and phase change respectively at the transducer boundary 540.

[0027] This corresponds with the component of the incident wave that is reflected from surface 540, where the incident wave is given by equation (15). After another passage of the wave through the material to surface 530 and back to surface 540, inclusive of reflection at surface 430 and attenuation, the amplitude at location 8 is:

[0028] After n passages through the component (where a passage consists of travelling from the transducer 510 to the surface 530 and back again), the amplitude of the wave is: Π. AQ e i[^-t)-n0-(n-l)0'] e 2nLa R n R I (n-l) ( 1 8)

[0029] The influence of reflection and phase change can be visualised using a spiral plot (figure 6) where successive reflection vectors are plotted on a polar grid with a cumulative phase shift. A greater a, smaller R and ft' will lead to a rapidly decaying spiral. The

measurement system is at its most sensitive (to measure R) when the transducer reflection, R' is close to unity and a is small.

[0030] Figure 5 illustrates an arrangement in which there is a direct line of sight between the transducer 510 and the surface 530. Further, in the arrangement of Figure 5, the transducer 510/surface 540 in contact with the transducer is parallel to the surface 530. This arrangement (or arrangements in which sufficiently approximate these conditions) allow for reduced complexity in molding the system, simplifying the analysis.

[0031] Figure 6a illustrates reflection amplitude and phase represented on a polar plot for a reference measurement in which the surface of interest 530 is a boundary between the component and air. Figure 6b illustrates reflection amplitude and phase represented on a polar plot for a measurement where there is a medium at the surface of interest 530, in this case the signal decays more rapidly (than that in Figure 6a) as energy is transmitted into the medium. The medium may be a layer, coating, film, etc. and these terms are used interchangeably herein.

[0032] Note that equations (12) to (16) merely describe a particular component of the wave, and do not take into account the effect of interference by portions of the wave emitted by the transducer 510 at different times and travelling through the component 520. All the reflected waves located at the transducer of Figure 5 superimpose, so the wave amplitude detected at the transducer is then the sum of (15), (17), and (18):

A R = Ao ^ e'^-^-^-^^R'^e 2 ^ (19) [0033] In the case where the reflection is between a component and air, then R=^ and Φ=0 and the amplitude is given by:

[0034] This arrangement (in which the surface of interest 530 is a boundary between the component 520 and air) may be used as a reference arrangement. Other reference

arrangements are possible (such as a known starting condition of a system, for example).

[0035] If we denote the superimposed standing wave reflection coefficient, S as the proportion of the superimposed incident wave that is reflected back to the transducer relative to the reference arrangement of equation (20), then S(t) is given by a ratio of equations (19) and (20): e i[w(-^ ~ t)-n0-(n-l)0'] R R ,(n- l ) e 2nLa

S(t) = (21 )

y∞_ β ί[ω(— -t)-(n- i )0'] R ,("- i ) e 2nLa

[0036] S is periodic with time. Taking the signal root mean square (RMS) provides the DC response of the system. This may be implemented practically by a rectification and smoothing process, for example.

[0037] If S is measured, e.g. by measuring the amplitude at the transducer in a reference state and a state of interest (also referred to herein as a measurement state) and taking a ratio, and knowing the parameters, R', Φ', L, and a, then R may be found using equation (22). Depending on the nature of the surface medium, R may be substituted from equation (5a/b), (8) or (10) to find the surface medium properties.

Practical Generation of a Superimposed Standing Wave

Driving Frequency Sweep

[0038] For any given component and layer system the driving frequency that causes the maximum amplitude superimposed standing wave is not necessarily known in advance.

According to some examples, the transducer is driven so as to sweep through a series of frequencies and record the response for all those frequencies. Accordingly, driving frequencies that produce a maximum measured amplitude may be determined in a convenient manner. Figure 7 schematically shows an ideal input frequency sweep. Once completed, the sweep may repeated continuously (three such sweeps are shown in the figure). Other sweeping schemes are possible, such as a scheme in which the frequency is increased in one sweep and decreased in the following sweep.

[0039] In some examples, each frequency sweep is identical. However, in some examples successive frequency sweeps may differ from each other, for example by being over different ranges, having different frequency increments, and/or a different rate of change of frequency, for example. For example, a driving frequency sweep is applied to determine the standing wave maximum frequencies (i.e. frequencies at which the amplitude of the standing wave is locally maximal in the frequency domain); however, once it is known that those frequencies occur within a particular range or bandwidth, a more closely defined sweep (or series of driving frequencies) might be used.

[0040] In some examples the frequency sweep may include discontinuous subranges of frequencies. For example, where the frequencies of a plurality of standing wave maxima are known, the frequency sweep may include frequencies only at/around each of these maxima, omitting frequencies at/around minima of the superimposed standing wave.

Transducer Response

[0041] Real piezoelectric transducers have a finite bandwidth and so do not generate or receive the full range of frequencies in the manner shown in figure 7. Figure 8 schematically shows a typical transducer response function or envelope for the three frequency sweeps of figure 7. The lower and higher frequencies of each sweep are generated at substantially reduced amplitudes.

Formation of the Superimposed Standing Wave

[0042] The driving frequency sweep of figure 8 superimposes with all its reflections and sub- reflections, as described above with reference to Figure 5, to create the superimposed standing wave shown schematically in figure 9. At some points in time (i.e. corresponding to a specific driving frequency) the received signal is at a maximum as constructive interference occurs. At other times, the amplitude is minimum, corresponding to destructive interference. Essentially, the component itself has a modulating effect on the input frequencies.

[0043] In Figure 9 the constituent wave amplitudes are shown only in the first sweep, for clarity. Only the envelope of the peaks is shown in the second and third sweeps. In practice, rather than using the peak envelope it is more convenient to use the RMS of the signal. And since every time corresponds to a driving frequency (since the driving frequency is swept in time), the RMS of the signal may be plotted against the driving frequency. This is shown schematically in Figure 10.

[0044] The RMS amplitude of this superimposed standing wave (with respect to frequency) may be used to derive properties of the surface of interest. The pitch of the maxima in the RMS signal is the standing wave superposition frequency, f s .

Effect of a Medium at the Surface

[0045] The addition of a surface layer to the surface of interest of the component causes the reflection coefficient at the interface, R to reduce and the phase Φ to change relative to a reference state with air at the surface of interest. If the layer is thick then it will also have resonant frequencies, f s i that superimpose on those in the component. The state of the component that is to be measured (in this case, the state in which the surface layer has been added) is referred to herein as the state of interest. This change in R and Φ changes the coherent frequency at which a superimposed standing wave forms. Accordingly, the superposed standing wave in the state of interest is changed, relative to the reference state, in amplitude and in the frequency, f s , at which the maximum amplitude occurs. Figure 1 1a shows the RMS standing wave pattern before and after the layer is added (i.e. in the reference state and state of interest). The layer causes the amplitude and frequency to reduce and also a further resonance is superimposed. A Figure 11 b shows the same figure zoomed in around the low frequency regime. For very thin layers the resonant frequency of the layer f may be outside the transducer bandwidth and so not observed in the spectrum. The spectrum would then resemble Figure 1 1 b throughout the bandwidth. Figure 11 b illustrates an example in which the amplitude of the standing wave is reduced by AA and the location of the maximum amplitude frequency is reduced by an amount Af s .

Reduction in Amplitude

[0046] Equation (22) describes the ratio of the amplitude of the superimposed standing wave for the measurement case divided by the reference case. As such, a ratio of the two RMS signals shown in figure 11 gives the standing wave reflection coefficient, S. Once S has been obtained, equation (22) may be used obtain R for the solid-liquid interface. In principle, S may be used at any of the measurement frequencies (and that frequency may be used in equation (22)). However, in practice the best signal to noise is obtained at frequencies around the maximum amplitude points, nf s . Once R is obtained this can be used in the relationships (5a/b), (8) or (10) to determine the properties of the liquid or thin liquid layer.

Change in Frequency

[0047] The presence of the layer, provided it is thick enough to cause resonance within the bandwidth frequencies, will superimpose those resonances onto the signal (f on figure 1 1a). The thickness of the layer can be obtained directly from this resonance frequency and the speed of sound in the layer material.

[0048] The component, the surface medium, and the bonded transducer (if present) behave as a composite system (i.e. response of the layer combines to give total system response and hence, peak positions). When there is no surface medium present the equivalent speed of sound through the component (path length, L) and the transducer (path length, I) combined is: (23)

L

where c is the wave speed in the component, and c t is the wave speed in the transducer material. The assembly first resonates (n=1) when the wavelength A=2(L+i) and so the location of the peak frequency is at:

_ (Lc+lc t ) (24)

J s - 2(L+0 2 K '

[0049] When there is a surface medium present the equivalent speed of sound through the component, transducer, and surface layer combined is:

Lc <25 \ where h is the path length in the surface layer and G is the wave speed in the surface medium material. The peak frequency is then increased above f s by an amount Af s :

= Lc + lc t+ nci

JS Js 2(L+l+h) 2 '

[0050] Af s may have a negative value, as in Figure 11 , for example.

[0051] Combining equations (24) and (26) gives the change in the peak frequency as:

_ Lc+lc t +hci (Lc+lc t )

h ~ 2{L+l+h) 2 2(L+0 2 ^ '

[0052] The layer thickness, h, may be determined by solving the quadratic equation (27). Higher Order Resonances

[0053] If the wave amplitudes from the reference (without a surface layer/medium) and coated component (with a surface layer/medium) are compared then an interference pattern is created (as shown in figure 12). The first order peaks are separated by Af s , the second order by 2Af s , and the nth order by nAf s . As the resonant order, n increases the separation increases until the peaks are back in phase when

[0054] Measuring the change in resonant frequency at a higher order harmonic, nAf s provides a very accurate way of measuring Af s and hence for deriving the layer thickness using equation (27). The higher the harmonic order, n the larger the value of nAfs and hence the more accurately Af s can be determined. Eventually the separation is nAf s =f s and the peaks are back in phase again.

[0055] If the standing wave amplitude spectrum for the measurement state is superimposed over the spectrum for the reference state an interference pattern is produced (as shown in Figure 12c). Determination of the peaks and valleys in this interference pattern is another accurate way to determine the change in resonant frequency f s .

Apparatus

[0056] An exemplary apparatus 300 suitable for use in some examples may include four sub- systems, as shown in Figure 13. The four subsystems are:

(i) A pair of transmitting 310 and receiving 315 transducers may be coupled to a back face of the component 520, where the back face is a face other than the surface of interest.

According to some examples the piezoelectric transducers 310, 315 may be

designed/arranged to operate in the 1-50 MHz frequency range.

(ii) A signal generator 320 may be arranged to provide a sinusoidal varying voltage at a tuneable range of frequencies. Typically driving voltages in the order of ±5V are used. The varying voltage is provided to the transmitting transducer 310 (emitter) to drive the transmitting transducer 310. The varying voltage may also be provided to a digitizer 330. The signal generator may be controlled by a PC 340. In some examples, the PC 340 may receive information from the signal generator indicative of the varying voltage provided to the transmitting transducer 310.

(iii) Digitizer 330 (illustrated in Figure 13 as a digital oscilloscope) may be arranged to receive a signal from the receiving transducer 315, and may record the waveform and/or pass it to PC 340 for processing.

(iv) PC 340 may be used for data storage and/or processing. In some examples the PC 340 may be used to control the signal generator 320 and/or digitizer 330. Any suitable computing device may be used in place of the PC.

[0057] The transmitting transducer 310 may be excited with a sinusoidal varying voltage that can be varied on command. The transducer converts this voltage to a continuous mechanical wave that propagates and reflects through the component 520. The wave arriving at the receiving transducer 315 is converted to an analogue voltage signal. This signal is captured, digitized by the digitizer 330, and passed to the PC 340 for processing.

[0058] According to some examples, the apparatus of Figure 13 may be realised using standard off the shelf instrumentation, such as an arbitrary function generator, a digital oscilloscope, and a conventional PC. Alternatively, in some examples, off-the-shelf integrated circuits (IC's) may be used to perform the key functions.

[0059] The type of transducer 310 used to generate a mechanical wave in the component 520 is not particularly limited. The type of transducer (its topology, polarisation, frequency, and geometry) may be selected to match the interface effect to be measured. By way of example, the following is a non-exclusive list of possible transducer types:

1. A longitudinally polarised piezoelectric element. A frequency rage of 1-50MHz may be used to measure the thickness of a surface layer, or the thickness of a trapped liquid film.

2. A shear polarised piezoelectric element. A frequency range 1-50MHz may be used to measure the viscosity of an adjacent liquid or of a liquid surface layer.

3. Non-contact driving may be achieved using an Electromagnetic Acoustic Transducer (EMAT).

4. A sputtered piezoelectric film may be used as the transmitter and receiver.

[0060] According to some examples the emitter and receiver may be separate elements (e.g. different transducers or a split transducer (i.e. the emitter and receiver are discrete

components). In some examples the emitter and receiver may be implemented using a single transducer. However, when the emitter and receiver are separate, the measured data is more easily extracted from the driving signal for the emitter.

[0061] Figure 14 illustrates a device 100 according to some examples. The device includes a processing section 110, suitable for carrying out processing operations.

[0062] The processing section 1 10 may be arranged to receive, via an input 120, ultrasound data relating to continuous wave, swept frequency ultrasound propagated through the component. The processing section 1 10 may be arranged to derive a standing wave amplitude spectrum from the received ultrasound data, and determine a property of an unconstrained layer on the surface based on the standing wave amplitude spectrum. The processing section 1 10 may be further arranged to output the determined property of the layer via output 130.

[0063] The device of Figure 14 may be, or may be included in, the PC of Figure 13.

[0064] The processing section may be arranged to derive the standing wave amplitude spectrum and determine a property of the unconstrained layer in accordance with the techniques described herein.

[0065] Herein an unconstrained layer is a layer having a free surface (e.g. in contact with a gas). Accordingly, the layer illustrated in Figure 2c is not an unconstrained layer, since both of its (major) surfaces are in contact with a solid component.

[0066] The ultrasound data may relate to a measurement state of the component 520, i.e. a state of interest of the component.

[0067] The input 120 and output 130 may be, for example physical ports, wiring for carrying data (e.g. electrical leads, data busses, etc.), registers or memory for storing data, etc.

Method

[0068] The transmitting element 310 may be excited with a frequency sweep of 100ms duration, for example. The frequency sweep may be a sequence of increasing frequencies that are generated and then repeated every 100ms, for example. The RMS of the signal at the receiving element 315 is monitored and the superimposed standing wave pattern is recorded.

[0069] Initially, when there is no surface layer or other element in contact with the component a reference measurement is recorded; in this arrangement there is negligible transmission into the air. This is the reference measurement as given by equation (20). Alternatively, other arrangements could be used as a reference state. For example, any arrangement having negligible transmission at the surface of interest, e.g. where the surface of interest is in contact with vacuum or a gas other than air. Arrangements in which there is non-negligible

transmission at the surface of interest may also be used if the parameters (such as the reflection coefficient and phase shift at the surface of interest in the reference state) are known. Use of an arrangement having negligible transmission at the surface of interest may reduce complexity relative to an arrangement having non-negligible transmission at the surface of interest.

[0070] When the component 520 is in the state of interest (e.g. when a surface layer has been applied) a superimposed standing wave is formed that is different from the superimposed standing wave in the reference state (assuming the state of interest is different form the reference state). The superimposed standing wave in the state of interest is given by equation (19). The RMS of this new signal may be recorded.

[0071] According to some examples, the measurement signal from the state of interest (e.g. having a layer on the surface of interest) is divided by that from the reference measurement to evaluate S. According to some examples S may be evaluated at a particular frequency and equation (22) may be used to evaluate R. In principle, any frequency can be used to evaluate R in this manner, and essentially the same layer properties should be derived. However, where the frequency used is at or around the standing wave maxima nf s the signal to noise ratio is expected to be at its highest, and so more accurate results are to expected compared with a frequency away from the standing wave maxima. In some examples, S may be derived at a plurality of frequencies and may be used together to derive a value of R; for example by determining respective values of R for each of the frequencies and obtaining a single value of R based on the values determined for the different frequencies (e.g. by averaging).

[0072] In some examples, S may be evaluated at a plurality of frequencies (e.g. at each frequency within the sweep, within the accuracy of the transmitting and/or measuring

equipment) in order to evaluate S(f). Equation (22) may then be used to find R, for example by performing a fit of equation (22) to the evaluated S(f) using a plot fitting or optimization technique. For example, one or more of the parameters of equation (22) may be allowed to vary and an optimization process may be performed to minimize the RMS error between S(f) based on the measured data and S(f) according to equation (22).

[0073] When R has been determined, it may be used to find the properties of the surface medium using equations (8) or (10).

[0074] According to some examples, Af s may be determined based on the change in the frequency of the standing wave maxima compared to the reference signal. This value of Af s may be used in equation (27) to determine the thickness of the surface layer.

[0075] The spacing between any two corresponding superimposed standing wave maxima (one corresponding to the reference state and one corresponding to the state of interest) may be used to deduce f s , provided that the order n is known. However, the higher the wave number, n, the larger that spacing nAf s is.

[0076] Figure 15 illustrates a method 200 according to some examples. The method begins at 210, and at 220 ultrasound data is received, the date relating to continuous wave, swept frequency ultrasound that has propagated through the component. At 230 a standing wave amplitude spectrum is derived from the received ultrasound data. This may be performed in accordance with the preceding description, for example. At 240 a property of an unconstrained layer on the surface is determined based on the standing wave amplitude spectrum. This may be performed in a manner consistent with the previous description. For example, one or more of equations (8), (10), (22) and/or (27) may be used.

Conclusion

[0077] According to some examples, a continuous wave is generated in a component 520. The wave propagates, reflects, and attenuates inside the component 520. The reflections superimpose to create a 'superimposed standing wave'.

[0078] The amplitude and frequency spectra of the superimposed standing wave are highly sensitive to the presence and properties of the surface (e.g. surface of interest) of the component 520, such as the properties of a medium at the surface. The ratio of the amplitude reflected from the surface (e.g. reflected by a medium) to that in a reference state (e.g. when there is no medium present) is referred to herein as the standing wave reflection coefficient, S. A mathematical description of the superimposed standing wave can be used along with the value of the superimposed standing wave derived from measurement to find the reflection coefficient at the surface of interest (e.g. the surface/medium interface.)

[0079] Exemplary applications of some examples are measurement of: surface layers (such as paint thickness), viscosity of a thin layer; thickness and properties of a thin surface film. Some examples may have one or more of increased sensitivity, lower response time, and/or low cost instrumentation when compared with conventional pulse-echo type ultrasonic techniques.

[0080] References herein to particular implementations of hardware, such as the use of transducers, are not particularly limiting and any suitable hardware may be used to implement the various examples.

[0081] According to some examples, devices described herein, or elements thereof may include, for example, one or more processors, a personal computer, an Application Specific Integrated Circuit (ASIC), and/or memory. The devices or elements may be arranged to execute software and/or firmware instructions. For example, processing section 1 10 may be embodied in this manner.

[0082] Throughout the description and claims of this specification, the words "comprise" and "contain" and variations of them mean "including but not limited to", and they are not intended to (and do not) exclude other moieties, additives, components, integers or steps. Throughout the description and claims of this specification, the singular encompasses the plural unless the context otherwise requires. In particular, where the indefinite article is used, the specification is to be understood as contemplating plurality as well as singularity, unless the context requires otherwise.

[0083] Features, integers or characteristics described in conjunction with a particular aspect or example of the invention are to be understood to be applicable to any other aspect or example described herein unless incompatible therewith. All of the features disclosed in this

specification (including any accompanying claims, abstract and drawings), and/or all of the steps of any method or process so disclosed, may be combined in any combination, except combinations where at least some of such features and/or steps are mutually exclusive. The invention is not restricted to the details of any foregoing examples. The invention extends to any novel one, or any novel combination, of the features disclosed in this specification

(including any accompanying claims, abstract and drawings), or to any novel one, or any novel combination, of the steps of any method or process so disclosed.

[0084] The reader's attention is directed to all papers and documents which are filed concurrently with or previous to this specification in connection with this application and which are open to public inspection with this specification, and the contents of all such papers and documents are incorporated herein by reference.

Nomenclature

A Amplitude of the recorded signal

Ao Amplitude of driving signal emitted by the transducer

A, Amplitude of a wave incident at a boundary

A R Amplitude of a wave reflected at a boundary

c Speed of sound in the component 520

Cequiv Equivalent speed of sound for a layered body

Cpzt Speed of sound in the piezoelectric transducer

f Frequency of the ultrasonic wave

fs Frequency at which the superimposed standing wave amplitude is maximum

h Surface medium layer thickness

K Stiffness of the interface

L Component thickness

1 Transducer thickness

n Number of passages of the wave through the component

n Wave mode number

R Reflection coefficient at the interface of interest

R' Reflection coefficient at transducer face

S Ratio of the superimposed standing wave amplitude from a component with surface medium present, divided by that in the reference state without a surface medum. t Time

z Acoustic impedance of the material (product of the wave speed, c and the density p) a Attenuation coefficient of component material

Φ Phase shift at interface of interest

Φ' Phase shift at the transducer face

ω Angular frequency of the wave {=2πή

λ Wavelength

Case Study - Paint film thickness measurements

[0085] According to some examples, the thickness of multiple layers of paint sprayed onto a metallic surface may be accurately measured. The measurement of paint may be performed non-invasively from the counter face (back face) of the component as the paint is deposited. Metals used for cases of engineering components, for example, may consist of thin metallic shells of aluminium alloys, titanium alloys and steels, e.g. between 1.5-0.9 mm thick. Paints and coatings used to protect metallic surfaces from environmental affects may consist of an organic base with pigments and fillers.

[0086] Piezoelectric transducers were bonded to the back face of a piece of aluminium

(0.96mm thick). Acrylic paint was sprayed onto the sample in approximately 10μηι layers and the change in the standing wave maximum frequency after each layer had cured was used to determine the paint thickness. Each paint layer thickness was measured using a profilometer for verification of the ultrasound results. Figure 16 schematically illustrates the arrangement.

[0087] Figure 17 shows the standing wave amplitude spectrum measured from this system; this is the experimentally measured equivalent of the schematic Figure 10. A clear spacing of standing wave peaks is observed.

[0088] Figure 18 shows one such peak extracted from this data set, with measurements acquired as layers of paint are applied. With each additional layer the shift of peak frequency occurs as described above in relation to Figure 11. After the application of a paint layer, the location of the standing wave peaks move to the left as a result of the systems additional thickness.

[0089] Values for Af s were determined from the data in Figure 18 for each of the paint layer thicknesses and compared with the physically measured results (using the profilometer after the addition of 4 approximately equal layers; the profilometer is accurate to 10 μηι). Table 1 shows the results, and good agreement is observed.

[0090] Profilometer Ultrasonic

Paint layer number Change in f s (MHz)

measurement (μηι) measurement (μηι)

1 2.61 10 10.95

2 2.56 20 22.18

3 2.54 30 30.80

4 2.51 40 39.60

TABLE 1 : Paint thickness resull ts

[0091] According to some examples thickness measurement accuracies may yield precision to 100's of nm.

Case Study - Viscosity measurement of liquid layer

[0092] This case study describes measurement of the viscosity of Newtonian oil at the surface of a metal component from the back face. The method is able to measure the viscosity of fluid non-invasively in real time, and so can be used in process applications such as pipes, pressure vessels, and engine bearings.

[0093] Figure 19 shows a transducer pair (illustrated as a single element) bonded to an aluminium block with a layer of viscous oil on the free surface. The transducer was pulsed continuously using the method as described above and the superimposed standing wave recorded. Several Newtonian viscosity standard oils ranging from 38cP to 10008cP were used and the standing wave reflection coefficient, S determined for each.

[0094] Figure 20 shows a plot of the variation in S derived from the measurements with the oil viscosity as points. S is sensitive over a very wide range of viscosities. S was also calculated based on viscosities measured with a viscometer and equation (22), and this is plotted over the data as a line. Good agreement is observed.

[0095] Table 2 shows the measured standing wave reflection coefficient for the range of oils. This was used in equation (22) to find the reflection coefficient, R, this R was used in equation (10) to find the viscosity. The viscosity measured by a conventional viscometer compares well with that measured by the ultrasonic standing wave method.

[0096] Viscosity Measured by Standing Wave Reflection Ultrasonic Measured Viscometer, Coefficient (S) Viscosity mPa.s (cP) mPa.s (cP)

29.1 0.9904 21.1

102 0.9794 104

200 0.9751 156

754 0.9582 466

1329 0.9433 1160

1499 0.9456 1277

2508 0.9258 2348

3156 0.9234 2931

4096 0.9145 3934

5798 0.9052 5258

TABLE 2: A range of viscosity standard oils were used to evaluate the relationship between standing wave amplitude and viscosity.

[0097] The results demonstrate that the viscosity of a liquid can be determined in situ from the outer surface of a vessel. A clear correlation was seen between the viscosity measured from a conventional viscometer and the standing wave reflection coefficient S.