A Pseudo-3D Model for Electromagnetic Acoustic Transducers (EMATs)

: Previous methods for modelling Rayleigh waves produced by a meander-line-coil electromagnetic acoustic transducer (EMAT) consisted mostly of two-dimensional (2D) simulations that focussed on the vertical plane of the material. This paper presents a pseudo-three-dimensional (3D) model that extends the simulation space to both vertical and horizontal planes. For the vertical plane, we combines analytical and ﬁnite-difference time-domain (FDTD) methods to model Rayleigh waves’ propagation within an aluminium plate and their scattering behaviours by cracks. For the horizontal surface plane, we employ an analytical method to investigate the radiation pattern of Rayleigh waves at various depths. The experimental results suggest that the models and the modelling techniques are valid.


Introduction
A wide group of non-destructive testing (NDT) techniques are commonly used in biomedical industries, such as ultrasonic techniques, electromagnetic techniques, and laser testing [1][2][3][4]. Due to the non-contact nature, more and more attention has been paid to the NDT technique with electromagnetic acoustic transducers (EMATs), and EMATs have gradually been used in industrial applications, such as thickness gauging and defect detection [5][6][7][8].
A classic EMAT sensor is made of a meander-line-coil and a permanent magnet ( Figure 1). There are two major coupling principles for EMATs: magnetostriction is for ferromagnetic metallic materials, and the Lorentz force mechanism is for conductive and ferromagnetic materials [9]. This work focussed on only the Lorentz force mechanism performing on an aluminium plate. The Lorentz force mechanism is: the meander-line-coil placed above the sample generates eddy currents J within the sample. A permanent magnet placed above the coil generates a static magnetic field B to the sample. The interaction between J and B produces Lorentz force density F, as shown in Equation (1): Substantial works have been reported on EMATs modelling, which comprises an electromagnetic (EM) model and an ultrasonic (US) model [10][11][12]. The EM model was accomplished by the finite element method (FEM) and the analytical method, while the US model was accomplished with FEM, the finite-difference time-domain method (FDTD), and the analytical solutions. Some of the previous Most of EMATs' simulations were two-dimensional (2D), and can only focus on the specific plane. This article is attempting to build a pseudo-three-dimensional (3D) model in order to further study EMATs by combining a surface plane 2D model (the x-y cross-section shown in Figure 1) and a vertical plane 2D model (the y-z cross-section shown in Figure 1) together. More specifically, the Lorentz force density obtained from the vertical plane of the sample is imported to the surface plane of the sample as the driving source to generate Rayleigh waves. Previously, only the beam directivity at the surface of the sample (z = 0 in Figure 1) was investigated [18,19]. However, Rayleigh waves not only distribute along the surface (z = 0), they also distribute within a depth of the Rayleigh waves' wavelength. Some industrial defects are within the sample instead of on the surface of the sample; in order to use Rayleigh waves to detect such defects, the beam directivity of the Rayleigh waves at various depths is worth investigating. In this article, the equations to study the beam directivity are more complete compared to the approximate equations presented in [18] by Xie et al., and the beam directivity at various depths are investigated. This study lays a solid industrial foundation for nearsurface defects detection using Rayleigh waves, and can be a starting point to build an advanced 3D EMAT model in the future. Except for near-surface detection, the proposed strategy can be used to perform body detection, i.e., to model bulk waves, including longitudinal waves and shear waves. In addition, the proposed 3D EMAT model can be used to characterise other EMAT structures to generate surface waves, such as unidirectional Rayleigh waves EMATs and multiple directional Rayleigh waves EMATs. The 2D simulation on the vertical plane utilizes the analytical method and FDTD; the 2D simulation on the vertical plane and experimental validations are introduced in Section 2. The pseudo-3D model is presented in Section 3 to investigate the radiation pattern of Rayleigh waves at various depths by utilizing a wholly analytical solution. This work is an extension of the work published in [15,18] by Xie et al.

Vertical Plane Modelling
Previously, the authors have conducted EMAT modelling for Rayleigh waves focussing on the vertical plane of the material. This model combines the analytical method and FDTD to model EMATs [15]. The dimension and material of the test piece, the coil, and the permanent magnet are the same as the ones used in [15] by Xie et al. Based on such design, the working frequency used to form the interference of Rayleigh waves is 483 kHz. Reproduced with permission from [15], IEEE, 2016. Most of EMATs' simulations were two-dimensional (2D), and can only focus on the specific plane. This article is attempting to build a pseudo-three-dimensional (3D) model in order to further study EMATs by combining a surface plane 2D model (the x-y cross-section shown in Figure 1) and a vertical plane 2D model (the y-z cross-section shown in Figure 1) together. More specifically, the Lorentz force density obtained from the vertical plane of the sample is imported to the surface plane of the sample as the driving source to generate Rayleigh waves. Previously, only the beam directivity at the surface of the sample (z = 0 in Figure 1) was investigated [18,19]. However, Rayleigh waves not only distribute along the surface (z = 0), they also distribute within a depth of the Rayleigh waves' wavelength. Some industrial defects are within the sample instead of on the surface of the sample; in order to use Rayleigh waves to detect such defects, the beam directivity of the Rayleigh waves at various depths is worth investigating. In this article, the equations to study the beam directivity are more complete compared to the approximate equations presented in [18] by Xie et al., and the beam directivity at various depths are investigated. This study lays a solid industrial foundation for near-surface defects detection using Rayleigh waves, and can be a starting point to build an advanced 3D EMAT model in the future. Except for near-surface detection, the proposed strategy can be used to perform body detection, i.e., to model bulk waves, including longitudinal waves and shear waves. In addition, the proposed 3D EMAT model can be used to characterise other EMAT structures to generate surface waves, such as unidirectional Rayleigh waves EMATs and multiple directional Rayleigh waves EMATs. The 2D simulation on the vertical plane utilizes the analytical method and FDTD; the 2D simulation on the vertical plane and experimental validations are introduced in Section 2. The pseudo-3D model is presented in Section 3 to investigate the radiation pattern of Rayleigh waves at various depths by utilizing a wholly analytical solution. This work is an extension of the work published in [15,18] by Xie et al.
where A is the vector potential generated by I, ω and I are the angular frequency and the density of the applied alternating current (AC), respectively, ε, µ and σ are the permittivity, permeability and conductivity of the test piece respectively, and E and J are the induced electric field and eddy current density, respectively [20]. As described in [15] by Xie et al., for a small radius circular coil, the distribution of the vector potential A at z = 0 (surface of the sample) is not symmetrical with the radius due to the bent wire.
The coil used in this work consists of straight wires; thus, the analytical solution for a straight wire is needed. The adapted solution for a straight wire has been presented in [15] by Xie et al. Here is a brief introduction. A hypothesis is proposed: if the radius of the circular coil, compared with its width, is very large, a bent wire can be viewed as a straight wire, and the distribution of A should be symmetrical. To validate such a hypothesis, a model is built with a large-radius circular coil above the aluminium plate. The aluminium sample used has a dimension of 80 mm × 30 mm, and the inner radius and the outer radius of the circular coil are set to 5.0395 m and 5.0405 m, respectively. At kHz, the current density applied to the circular coil is 1 A/m 2 , and the lift-off of the coil is 1 mm. The permeability and the conductivity of the aluminium plate are 1.257 × 10 −6 H/m and 3.8 × 10 7 Siemens/m, respectively.
The distribution of the magnitude of A based on the adapted solution is shown in Figure 2. The red marker denotes the maximum magnitude of the vector potential. The distribution of the magnitude of A is symmetrical with the radius of 5.04 m, where the coil is located. The result verifies the hypothesis that, when the radius of the circular coil is very large, a bent wire serves as a straight wire.

EMAT-EM Model
This section introduces the EMAT-EM model to analyse the distribution of (Lorentz force density). Firstly, the classic Dodd and Deeds solution [20] to the vector potential is described, and the strategy of adapting the circular analytical solutions for a straight wire is introduced (Section 2.1.1). The FEM is employed to validate the adapted solution (Section 2.1.2). Finally, based on the adapted analytical solutions, the distribution of Lorentz force density is presented (Section 2.1.3).

Adapted Analytical Solutions to the Vector Potential for a Straight Wire
The governing equations for calculating eddy currents are described in Equations (2)-(4): where A is the vector potential generated by I, ω and I are the angular frequency and the density of the applied alternating current (AC), respectively, ε, µ and σ are the permittivity, permeability and conductivity of the test piece respectively, and E and J are the induced electric field and eddy current density, respectively [20]. As described in [15] by Xie et al., for a small radius circular coil, the distribution of the vector potential A at z = 0 (surface of the sample) is not symmetrical with the radius due to the bent wire.
The coil used in this work consists of straight wires; thus, the analytical solution for a straight wire is needed. The adapted solution for a straight wire has been presented in [15] by Xie et al. Here is a brief introduction. A hypothesis is proposed: if the radius of the circular coil, compared with its width, is very large, a bent wire can be viewed as a straight wire, and the distribution of A should be symmetrical. To validate such a hypothesis, a model is built with a large-radius circular coil above the aluminium plate. The aluminium sample used has a dimension of 80 mm × 30 mm, and the inner radius and the outer radius of the circular coil are set to 5.0395 m and 5.0405 m, respectively. At 1 kHz, the current density applied to the circular coil is 1 A/m 2 , and the lift-off of the coil is 1 mm. The permeability and the conductivity of the aluminium plate are 1.257 × 10 −6 H/m and 3.8 × 10 7 Siemens/m, respectively.
The distribution of the magnitude of A based on the adapted solution is shown in Figure 2. The red marker denotes the maximum magnitude of the vector potential. The distribution of the magnitude of A is symmetrical with the radius of 5.04 m, where the coil is located. The result verifies the hypothesis that, when the radius of the circular coil is very large, a bent wire serves as a straight wire.

Comparison between the Adapted Solution and FEM
In order to compare the adapted solutions to FEM, Maxwell Ansoft, which is a FE solver, is utilised. The FEM model has a rectangular cross-sectional coil located above the cross-sectional aluminium plate, and is surrounded by a vacuum region that is four times larger than the sample. The FEM subdivides the large model to smaller elements, and this FEM solves the calculation by minimising the energy error. In this work, when the elements number is beyond 20,000, the energy error is as low as 0.05%, which is sufficiently accurate for the FEM computation. In this work, the mesh number used is 20,395, and the boundary used is a balloon boundary to simulate an infinite space. The distribution of A at the sample's surface (z = 0) is presented in Figure 3. The analytical solution and FEM present a good agreement at an operational frequency of 1 kHz. However, at a working frequency of 1 MHz, the distribution of A from the FEM is not smooth compared to that from the analytical solution; the reason is that the FEM is affected by the elements density and numerical approximation is unavoidable, etc. Therefore, the adapted analytical method presents a more accurate result compared to FEM, especially for a high working frequency. In order to compare the adapted solutions to FEM, Maxwell Ansoft, which is a FE solver, is utilised. The FEM model has a rectangular cross-sectional coil located above the cross-sectional aluminium plate, and is surrounded by a vacuum region that is four times larger than the sample. The FEM subdivides the large model to smaller elements, and this FEM solves the calculation by minimising the energy error. In this work, when the elements number is beyond 20,000, the energy error is as low as 0.05%, which is sufficiently accurate for the FEM computation. In this work, the mesh number used is 20,395, and the boundary used is a balloon boundary to simulate an infinite space. The distribution of A at the sample's surface (z = 0) is presented in Figure 3. The analytical solution and FEM present a good agreement at an operational frequency of 1 kHz. However, at a working frequency of 1 MHz, the distribution of A from the FEM is not smooth compared to that from the analytical solution; the reason is that the FEM is affected by the elements density and numerical approximation is unavoidable, etc. Therefore, the adapted analytical method presents a more accurate result compared to FEM, especially for a high working frequency. Reproduced with permission from [15], IEEE, 2016.

EMAT-Lorentz Force Calculation
The analytical solution to a straight wire has been described in Section 2.1.2. A (vector potential) generated by a meander-line-coil is the addition of A generated by every single straight wire. The zone where the meander-line-coil mainly operates on is selected to model EM simulation to increase The red curve is the result from the FEM, and the blue curve is the result from the analytical solutions. Reproduced with permission from [15], IEEE, 2016.

EMAT-Lorentz Force Calculation
The analytical solution to a straight wire has been described in Section 2.1.2. A (vector potential) generated by a meander-line-coil is the addition of A generated by every single straight wire. The zone where the meander-line-coil mainly operates on is selected to model EM simulation to increase the modelling effectiveness. The distribution of A and F at z = 0 are shown in Figures 4 and 5. The generated periodic fields have different directions for any neighbouring wires, since their applied ACs are opposite, and therefore, the periodic distribution of A and F has six positive values and six negative values. The outermost A and outermost F are largest, because A is under the outermost wires, and thus is only determined by the fields on one side.

Elastodynamic Equations
Elastodynamic equations (Equations (5) and (6)) are a group of partial differential equations to model the wave propagation: where and are the density and the fourth stiffness tensor of the material, and and are the force source and strain tensor rate source, respectively.

Combination of EMAT-EM and EMAT-US Models
As described in Section 2.1.3, F is obtained from the EMAT-EM calculation. In this section, F,

Elastodynamic Equations
Elastodynamic equations (Equations (5) and (6)) are a group of partial differential equations to model the wave propagation: where and are the density and the fourth stiffness tensor of the material, and and are the force source and strain tensor rate source, respectively.

Combination of EMAT-EM and EMAT-US Models
As described in Section 2.1.3, F is obtained from the EMAT-EM calculation. In this section, F,

Elastodynamic Equations
Elastodynamic equations (Equations (5) and (6)) are a group of partial differential equations to model the wave propagation: where ρ and c ijkl are the density and the fourth stiffness tensor of the material, and f i and θ ij are the force source and strain tensor rate source, respectively.

Combination of EMAT-EM and EMAT-US Models
As described in Section 2.1.3, F is obtained from the EMAT-EM calculation. In this section, F, which is used as the force source, is imported to the EMAT-US model to produce ultrasound ( Figure 6). Since F is calculated in the frequency domain and FDTD is a time-domain solver, the excitation signal for the EMAT-US model is a time sequence signal with the peak equalling the peak values of F. The excitation signal used is a Gaussian-modulated sinusoidal with a fractional width of 0.18. A crack and a receiver R are located within the sample, as shown in Figure 6. Regarding the FDTD setup in the ultrasonic model, there are two main parameters: the spatial step, and the time step. The spatial step used is 0.2 mm, which approximately equals to 1/30th of the wavelength. The time step is 0.0222 µs, which is calculated based on the Courant-Friedrichs-Lewy (CFL) condition. Free surface conditions are utilised on the surface of the sample. Perfectly-matched layers (PML) with a thickness of 16 mm are utilised to absorb ultrasound.

Wave Propagations
Based on FDTD calculations, the velocity fields of ultrasound waves are obtained. The velocity fields at 27 µs and 83 µs are shown in Figure 7, which describes the Rayleigh waves' propagation and their scattering behaviours, respectively. The white arrows denote the propagation path.

EMAT-Reception Simulation
The EMAT reception process has been reported quite a lot [21]. The received signal from simulations is shown in Figure 8, where DRW and RRW denote the directly transmitted Rayleigh waves and the reflected Rayleigh waves, respectively. The propagation distance of DRW and RRW is 100 mm and 300 mm, respectively, and the velocity of Rayleigh waves is 2.93 mm/µs; hence, the theoretically arrival time of DRW and RRW is 34 µs and 102.4 µs, respectively. Figure 8 shows the numerically arrival time of DRW and RRW; these numerical and experimental results present a good

Wave Propagations
Based on FDTD calculations, the velocity fields of ultrasound waves are obtained. The velocity fields at 27 µs and 83 µs are shown in Figure 7, which describes the Rayleigh waves' propagation and their scattering behaviours, respectively. The white arrows denote the propagation path.

Wave Propagations
Based on FDTD calculations, the velocity fields of ultrasound waves are obtained. The velocity fields at 27 µs and 83 µs are shown in Figure 7, which describes the Rayleigh waves' propagation and their scattering behaviours, respectively. The white arrows denote the propagation path.

EMAT-Reception Simulation
The EMAT reception process has been reported quite a lot [21]. The received signal from simulations is shown in Figure 8, where DRW and RRW denote the directly transmitted Rayleigh waves and the reflected Rayleigh waves, respectively. The propagation distance of DRW and RRW is 100 mm and 300 mm, respectively, and the velocity of Rayleigh waves is 2.93 mm/µs; hence, the theoretically arrival time of DRW and RRW is 34 µs and 102.4 µs, respectively. Figure 8 shows the numerically arrival time of DRW and RRW; these numerical and experimental results present a good

EMAT-Reception Simulation
The EMAT reception process has been reported quite a lot [21]. The received signal from simulations is shown in Figure 8, where DRW and RRW denote the directly transmitted Rayleigh waves and the reflected Rayleigh waves, respectively. The propagation distance of DRW and RRW is 100 mm and 300 mm, respectively, and the velocity of Rayleigh waves is 2.93 mm/µs; hence, the theoretically arrival time of DRW and RRW is 34 µs and 102.4 µs, respectively. Figure 8 shows the numerically arrival time of DRW and RRW; these numerical and experimental results present a good agreement.

Experimental Validations
Experiments were carried out to validate the proposed modelling method. The setup of the experiments is the same as that of our previous work [15]. The received signals from experiments are shown in Figure 9, where three signals are observed. The "Main bang" is the interference signal due to a high power excitation, arriving before DRW and RRW. The red curve and the blue curve are the envelope and the time series signal, respectively.

Experimental Validations
Experiments were carried out to validate the proposed modelling method. The setup of the experiments is the same as that of our previous work [15]. The received signals from experiments are shown in Figure 9, where three signals are observed. The "Main bang" is the interference signal due to a high power excitation, arriving before DRW and RRW. The red curve and the blue curve are the envelope and the time series signal, respectively.

Experimental Validations
Experiments were carried out to validate the proposed modelling method. The setup of the experiments is the same as that of our previous work [15]. The received signals from experiments are shown in Figure 9, where three signals are observed. The "Main bang" is the interference signal due to a high power excitation, arriving before DRW and RRW. The red curve and the blue curve are the envelope and the time series signal, respectively.    permission from [15], IEEE, 2016. Figure 10 shows the envelope of DRW and RRW from experiments and simulations. The experimental arrival times of DRW is 34 µs, which is consistent with the numerical results. However, the experimental arrival time of RRW is slightly different from the simulations. This is because of the approximated model used and the inevitable experimental noise.

Horizontal Surface Plane Modelling-Directivity Analysis of Rayleigh Waves
The main mode of propagation on the horizontal surface of the material is the Rayleigh wave. The radiation pattern and the beam directivity of Rayleigh waves at z = 0 ( Figure 1) were reported in [18] by Xie et al. Rayleigh waves not only distribute along the surface (z = 0), they also distribute within a depth of the Rayleigh waves' wavelengths. However, Rayleigh waves not only concentrate on the surface, they also concentrate within a depth of one wavelength. In this section, the beam directivity of Rayleigh waves at various depths are investigated utilising an analytical method.

The Analytical Solution to the Displacement of Rayleigh Waves
N. A. Haskell proposed an analytical solution to Rayleigh waves [22,23]. Ref. [24] by Love introduced these solutions in detail, and investigated the beam directivity of Rayleigh waves at z = 0 with approximated equations. However, Rayleigh waves not only concentrate on the surface; they also concentrate within a depth of one wavelength. The complete equations of the Rayleigh waves' displacement at various depths are: 0. The envelope of the received signals. The blue curve and the red curve are the envelope of ved signal from simulations and experiments, respectively.

l Surface Plane Modelling-Directivity Analysis of Rayleigh Waves
in mode of propagation on the horizontal surface of the material is the Rayleigh wave. pattern and the beam directivity of Rayleigh waves at z = 0 ( Figure 1) were reported in t al. Rayleigh waves not only distribute along the surface (z = 0), they also distribute th of the Rayleigh waves' wavelengths. However, Rayleigh waves not only concentrate ce, they also concentrate within a depth of one wavelength. In this section, the beam Rayleigh waves at various depths are investigated utilising an analytical method.
ytical Solution to the Displacement of Rayleigh Waves askell proposed an analytical solution to Rayleigh waves [22,23]. Ref. [24] by Love hese solutions in detail, and investigated the beam directivity of Rayleigh waves at z = 0 imated equations. However, Rayleigh waves not only concentrate on the surface; they rate within a depth of one wavelength. The complete equations of the Rayleigh waves' t at various depths are: d are the in-plane and the out-of-plane displacement to be calculated, respectively; nce between the source point and the field point ( Figure 11); is the excitation source; ity of the material; is the angle between the force vector and the in-plane displacement the operational angular frequency; z is the depth; and , , , and ĸ are the velocity dinal waves, shear waves, Rayleigh waves, and the wave number, respectively.

. Horizontal Surface Plane Modelling-Directivity Analysis of Rayleigh Waves
The main mode of propagation on the horizontal surface of the material is the Rayleigh wave. he radiation pattern and the beam directivity of Rayleigh waves at z = 0 ( Figure 1) were reported in 18] by Xie et al. Rayleigh waves not only distribute along the surface (z = 0), they also distribute ithin a depth of the Rayleigh waves' wavelengths. However, Rayleigh waves not only concentrate n the surface, they also concentrate within a depth of one wavelength. In this section, the beam irectivity of Rayleigh waves at various depths are investigated utilising an analytical method.
= ĸ( − 1) here: here and are the in-plane and the out-of-plane displacement to be calculated, respectively; is the distance between the source point and the field point ( Figure 11); is the excitation source; is the density of the material; is the angle between the force vector and the in-plane displacement ector; is the operational angular frequency; z is the depth; and , , , and ĸ are the velocity f the longitudinal waves, shear waves, Rayleigh waves, and the wave number, respectively.

Horizontal Surface Plane Modelling-Directivity Analysis of Rayleigh Waves
The main mode of propagation on the horizontal surface of the material is the Rayleigh wave. The radiation pattern and the beam directivity of Rayleigh waves at z = 0 ( Figure 1) were reported in [18] by Xie et al. Rayleigh waves not only distribute along the surface (z = 0), they also distribute within a depth of the Rayleigh waves' wavelengths. However, Rayleigh waves not only concentrate on the surface, they also concentrate within a depth of one wavelength. In this section, the beam directivity of Rayleigh waves at various depths are investigated utilising an analytical method.

The Analytical Solution to the Displacement of Rayleigh Waves
N. A. Haskell proposed an analytical solution to Rayleigh waves [22,23]. Ref. [24] by Love introduced these solutions in detail, and investigated the beam directivity of Rayleigh waves at z = 0 with approximated equations. However, Rayleigh waves not only concentrate on the surface; they also concentrate within a depth of one wavelength. The complete equations of the Rayleigh waves' displacement at various depths are: = ĸ( − 1) where: (ĸ, , ) = ĸ 2 4 ( 2 2 ĸ 3 ) = cos ( ) where and are the in-plane and the out-of-plane displacement to be calculated, respectively; is the distance between the source point and the field point ( Figure 11); is the excitation source; is the density of the material; is the angle between the force vector and the in-plane displacement vector; is the operational angular frequency; z is the depth; and , , , and ĸ are the velocity of the longitudinal waves, shear waves, Rayleigh waves, and the wave number, respectively. where: A( . Haskell proposed an analytical solution to Rayleigh waves [22,23]. Ref. [24] by Love these solutions in detail, and investigated the beam directivity of Rayleigh waves at z = 0 oximated equations. However, Rayleigh waves not only concentrate on the surface; they ntrate within a depth of one wavelength. The complete equations of the Rayleigh waves' ent at various depths are: = ĸ( − 1) = cos ( ) and are the in-plane and the out-of-plane displacement to be calculated, respectively; stance between the source point and the field point ( Figure 11); is the excitation source; nsity of the material; is the angle between the force vector and the in-plane displacement is the operational angular frequency; z is the depth; and , , , and ĸ are the velocity itudinal waves, shear waves, Rayleigh waves, and the wave number, respectively.

Horizontal Surface Plane Modelling-Directivity Analysis of Rayleigh Waves
The main mode of propagation on the horizontal surface of the material is the Rayleigh wave. The radiation pattern and the beam directivity of Rayleigh waves at z = 0 ( Figure 1) were reported in [18] by Xie et al. Rayleigh waves not only distribute along the surface (z = 0), they also distribute within a depth of the Rayleigh waves' wavelengths. However, Rayleigh waves not only concentrate on the surface, they also concentrate within a depth of one wavelength. In this section, the beam directivity of Rayleigh waves at various depths are investigated utilising an analytical method.

The Analytical Solution to the Displacement of Rayleigh Waves
N. A. Haskell proposed an analytical solution to Rayleigh waves [22,23]. Ref. [24] by Love introduced these solutions in detail, and investigated the beam directivity of Rayleigh waves at z = 0 with approximated equations. However, Rayleigh waves not only concentrate on the surface; they also concentrate within a depth of one wavelength. The complete equations of the Rayleigh waves' displacement at various depths are: = ĸ( − 1) where: (ĸ, , ) = ĸ 2 4 ( 2 2 ĸ 3 ) = cos ( ) where and are the in-plane and the out-of-plane displacement to be calculated, respectively; is the distance between the source point and the field point ( Figure 11); is the excitation source; is the density of the material; is the angle between the force vector and the in-plane displacement vector; is the operational angular frequency; z is the depth; and , , , and ĸ are the velocity of the longitudinal waves, shear waves, Rayleigh waves, and the wave number, respectively.

Horizontal Surface Plane Modelling-Directivity Analysis of Rayleigh Waves
The main mode of propagation on the horizontal surface of the material is the Rayleigh wave. The radiation pattern and the beam directivity of Rayleigh waves at z = 0 ( Figure 1) were reported in [18] by Xie et al. Rayleigh waves not only distribute along the surface (z = 0), they also distribute within a depth of the Rayleigh waves' wavelengths. However, Rayleigh waves not only concentrate on the surface, they also concentrate within a depth of one wavelength. In this section, the beam directivity of Rayleigh waves at various depths are investigated utilising an analytical method. [22,23]. Ref. [24] by Love introduced these solutions in detail, and investigated the beam directivity of Rayleigh waves at z = 0 with approximated equations. However, Rayleigh waves not only concentrate on the surface; they also concentrate within a depth of one wavelength. The complete equations of the Rayleigh waves' displacement at various depths are:

N. A. Haskell proposed an analytical solution to Rayleigh waves
= ĸ( − 1) where: (ĸ, , ) = ĸ 2 4 ( 2 2 ĸ 3 ) = cos ( ) where and are the in-plane and the out-of-plane displacement to be calculated, respectively; is the distance between the source point and the field point ( Figure 11); is the excitation source; is the density of the material; is the angle between the force vector and the in-plane displacement vector; is the operational angular frequency; z is the depth; and , , , and ĸ are the velocity of the longitudinal waves, shear waves, Rayleigh waves, and the wave number, respectively.

Horizontal Surface Plane Modelling-Directivity Analysis of Rayleigh Waves
The main mode of propagation on the horizontal surface of the material is the Rayleigh wave. The radiation pattern and the beam directivity of Rayleigh waves at z = 0 ( Figure 1) were reported in [18] by Xie et al. Rayleigh waves not only distribute along the surface (z = 0), they also distribute within a depth of the Rayleigh waves' wavelengths. However, Rayleigh waves not only concentrate on the surface, they also concentrate within a depth of one wavelength. In this section, the beam directivity of Rayleigh waves at various depths are investigated utilising an analytical method.

The Analytical Solution to the Displacement of Rayleigh Waves
= cos ( ) where and are the in-plane and the out-of-plane displacement to be calculated, respectively; is the distance between the source point and the field point ( Figure 11); is the excitation source; is the density of the material; is the angle between the force vector and the in-plane displacement vector; is the operational angular frequency; z is the depth; and , , , and ĸ are the velocity of the longitudinal waves, shear waves, Rayleigh waves, and the wave number, respectively.

Horizontal Surface Plane Modelling-Directivity Analysis of Rayleigh Waves
The main mode of propagation on the horizontal surface of the material is the Rayleigh wave. The radiation pattern and the beam directivity of Rayleigh waves at z = 0 ( Figure 1) were reported in [18] by Xie et al. Rayleigh waves not only distribute along the surface (z = 0), they also distribute within a depth of the Rayleigh waves' wavelengths. However, Rayleigh waves not only concentrate on the surface, they also concentrate within a depth of one wavelength. In this section, the beam directivity of Rayleigh waves at various depths are investigated utilising an analytical method. [22,23]. Ref. [24] by Love introduced these solutions in detail, and investigated the beam directivity of Rayleigh waves at z = 0 with approximated equations. However, Rayleigh waves not only concentrate on the surface; they also concentrate within a depth of one wavelength. The complete equations of the Rayleigh waves' displacement at various depths are:

N. A. Haskell proposed an analytical solution to Rayleigh waves
where: (ĸ, , ) = ĸ 2 4 ( = cos ( ) where and are the in-plane and the out-of-plane displacement to be calculated, respectively; is the distance between the source point and the field point ( Figure 11); is the excitation source; is the density of the material; is the angle between the force vector and the in-plane displacement vector; is the operational angular frequency; z is the depth; and , , , and ĸ are the velocity of the longitudinal waves, shear waves, Rayleigh waves, and the wave number, respectively.
here and are the in-plane and the out-of-plane displacement to be calculated, respectively; is the distance between the source point and the field point ( Figure 11); is the excitation source; is the density of the material; is the angle between the force vector and the in-plane displacement ector; is the operational angular frequency; z is the depth; and , , , and ĸ are the velocity f the longitudinal waves, shear waves, Rayleigh waves, and the wave number, respectively.

Horizontal Surface Plane Modelling-Directivity Analysis of Rayleigh Waves
The main mode of propagation on the horizontal surface of the material is the Rayleigh wave. The radiation pattern and the beam directivity of Rayleigh waves at z = 0 ( Figure 1) were reported in [18] by Xie et al. Rayleigh waves not only distribute along the surface (z = 0), they also distribute within a depth of the Rayleigh waves' wavelengths. However, Rayleigh waves not only concentrate on the surface, they also concentrate within a depth of one wavelength. In this section, the beam directivity of Rayleigh waves at various depths are investigated utilising an analytical method. [22,23]. Ref. [24] by Love introduced these solutions in detail, and investigated the beam directivity of Rayleigh waves at z = 0 with approximated equations. However, Rayleigh waves not only concentrate on the surface; they also concentrate within a depth of one wavelength. The complete equations of the Rayleigh waves' displacement at various depths are:

N. A. Haskell proposed an analytical solution to Rayleigh waves
where: = cos ( ) where and are the in-plane and the out-of-plane displacement to be calculated, respectively; is the distance between the source point and the field point ( Figure 11); is the excitation source; is the density of the material; is the angle between the force vector and the in-plane displacement vector; is the operational angular frequency; z is the depth; and , , , and ĸ are the velocity of the longitudinal waves, shear waves, Rayleigh waves, and the wave number, respectively.

Horizontal Surface Plane Modelling-Directivity Analysis of Rayleigh Waves
The main mode of propagation on the horizontal surface of the material is the Rayleigh wave. The radiation pattern and the beam directivity of Rayleigh waves at z = 0 ( Figure 1) were reported in [18] by Xie et al. Rayleigh waves not only distribute along the surface (z = 0), they also distribute within a depth of the Rayleigh waves' wavelengths. However, Rayleigh waves not only concentrate on the surface, they also concentrate within a depth of one wavelength. In this section, the beam directivity of Rayleigh waves at various depths are investigated utilising an analytical method. [22,23]. Ref. [24] by Love introduced these solutions in detail, and investigated the beam directivity of Rayleigh waves at z = 0 with approximated equations. However, Rayleigh waves not only concentrate on the surface; they also concentrate within a depth of one wavelength. The complete equations of the Rayleigh waves' displacement at various depths are:

N. A. Haskell proposed an analytical solution to Rayleigh waves
where: where and are the in-plane and the out-of-plane displacement to be calculated, respectively; is the distance between the source point and the field point ( Figure 11); is the excitation source; is the density of the material; is the angle between the force vector and the in-plane displacement vector; is the operational angular frequency; z is the depth; and , , , and ĸ are the velocity of the longitudinal waves, shear waves, Rayleigh waves, and the wave number, respectively.

Horizontal Surface Plane Modelling-Directivity Analysis of Rayleigh Waves
The main mode of propagation on the horizontal surface of the material is the Rayleigh wave. The radiation pattern and the beam directivity of Rayleigh waves at z = 0 ( Figure 1) were reported in [18] by Xie et al. Rayleigh waves not only distribute along the surface (z = 0), they also distribute within a depth of the Rayleigh waves' wavelengths. However, Rayleigh waves not only concentrate on the surface, they also concentrate within a depth of one wavelength. In this section, the beam directivity of Rayleigh waves at various depths are investigated utilising an analytical method. [22,23]. Ref. [24] by Love introduced these solutions in detail, and investigated the beam directivity of Rayleigh waves at z = 0 with approximated equations. However, Rayleigh waves not only concentrate on the surface; they also concentrate within a depth of one wavelength. The complete equations of the Rayleigh waves' displacement at various depths are:

N. A. Haskell proposed an analytical solution to Rayleigh waves
where: where and are the in-plane and the out-of-plane displacement to be calculated, respectively; is the distance between the source point and the field point ( Figure 11); is the excitation source; is the density of the material; is the angle between the force vector and the in-plane displacement vector; is the operational angular frequency; z is the depth; and , , , and ĸ are the velocity of the longitudinal waves, shear waves, Rayleigh waves, and the wave number, respectively.

Horizontal Surface Plane Modelling-Directivity Analysis of Rayleigh Waves
The main mode of propagation on the horizontal surface of the material is the Rayleigh wave. The radiation pattern and the beam directivity of Rayleigh waves at z = 0 ( Figure 1) were reported in [18] by Xie et al. Rayleigh waves not only distribute along the surface (z = 0), they also distribute within a depth of the Rayleigh waves' wavelengths. However, Rayleigh waves not only concentrate on the surface, they also concentrate within a depth of one wavelength. In this section, the beam directivity of Rayleigh waves at various depths are investigated utilising an analytical method.

The Analytical Solution to the Displacement of Rayleigh Waves
N. A. Haskell proposed an analytical solution to Rayleigh waves [22,23]. Ref. [24] by Love introduced these solutions in detail, and investigated the beam directivity of Rayleigh waves at z = 0 with approximated equations. However, Rayleigh waves not only concentrate on the surface; they also concentrate within a depth of one wavelength. The complete equations of the Rayleigh waves' displacement at various depths are: where: = cos ( ) where and are the in-plane and the out-of-plane displacement to be calculated, respectively; is the distance between the source point and the field point ( Figure 11); is the excitation source; is the density of the material; is the angle between the force vector and the in-plane displacement vector; is the operational angular frequency; z is the depth; and , , , and ĸ are the velocity

Horizontal Surface Plane Modelling-Directivity Analysis of Rayleigh Waves
The main mode of propagation on the horizontal surface of the material is the Rayleigh wave. The radiation pattern and the beam directivity of Rayleigh waves at z = 0 ( Figure 1) were reported in [18] by Xie et al. Rayleigh waves not only distribute along the surface (z = 0), they also distribute within a depth of the Rayleigh waves' wavelengths. However, Rayleigh waves not only concentrate on the surface, they also concentrate within a depth of one wavelength. In this section, the beam directivity of Rayleigh waves at various depths are investigated utilising an analytical method.

The Analytical Solution to the Displacement of Rayleigh Waves
N. A. Haskell proposed an analytical solution to Rayleigh waves [22,23]. Ref. [24] by Love introduced these solutions in detail, and investigated the beam directivity of Rayleigh waves at z = 0 with approximated equations. However, Rayleigh waves not only concentrate on the surface; they also concentrate within a depth of one wavelength. The complete equations of the Rayleigh waves' displacement at various depths are: where: = cos ( ) where and are the in-plane and the out-of-plane displacement to be calculated, respectively; is the distance between the source point and the field point ( Figure 11); is the excitation source; is the density of the material; is the angle between the force vector and the in-plane displacement

Horizontal Surface Plane Modelling-Directivity Analysis of Rayleigh Waves
The main mode of propagation on the horizontal surface of the material is the Rayleigh wave. The radiation pattern and the beam directivity of Rayleigh waves at z = 0 ( Figure 1) were reported in [18] by Xie et al. Rayleigh waves not only distribute along the surface (z = 0), they also distribute within a depth of the Rayleigh waves' wavelengths. However, Rayleigh waves not only concentrate on the surface, they also concentrate within a depth of one wavelength. In this section, the beam directivity of Rayleigh waves at various depths are investigated utilising an analytical method.

The Analytical Solution to the Displacement of Rayleigh Waves
N. A. Haskell proposed an analytical solution to Rayleigh waves [22,23]. Ref. [24] by Love introduced these solutions in detail, and investigated the beam directivity of Rayleigh waves at z = 0 with approximated equations. However, Rayleigh waves not only concentrate on the surface; they also concentrate within a depth of one wavelength. The complete equations of the Rayleigh waves' displacement at various depths are: where: where and are the in-plane and the out-of-plane displacement to be calculated, respectively; is the distance between the source point and the field point ( Figure 11); is the excitation source;

Horizontal Surface Plane Modelling-Directivity Analysis of Rayleigh Waves
The main mode of propagation on the horizontal surface of the material is the Rayleigh wave. The radiation pattern and the beam directivity of Rayleigh waves at z = 0 ( Figure 1) were reported in [18] by Xie et al. Rayleigh waves not only distribute along the surface (z = 0), they also distribute within a depth of the Rayleigh waves' wavelengths. However, Rayleigh waves not only concentrate on the surface, they also concentrate within a depth of one wavelength. In this section, the beam directivity of Rayleigh waves at various depths are investigated utilising an analytical method.

The Analytical Solution to the Displacement of Rayleigh Waves
N. A. Haskell proposed an analytical solution to Rayleigh waves [22,23]. Ref. [24] by Love introduced these solutions in detail, and investigated the beam directivity of Rayleigh waves at z = 0 with approximated equations. However, Rayleigh waves not only concentrate on the surface; they also concentrate within a depth of one wavelength. The complete equations of the Rayleigh waves' displacement at various depths are: where: where and are the in-plane and the out-of-plane displacement to be calculated, respectively; is the distance between the source point and the field point ( Figure 11); is the excitation source; < ω/c s

Horizontal Surface Plane Modelling-Directivity Analysis of Rayleigh Waves
The main mode of propagation on the horizontal surface of the material is the Rayleigh wave. The radiation pattern and the beam directivity of Rayleigh waves at z = 0 ( Figure 1) were reported in [18] by Xie et al. Rayleigh waves not only distribute along the surface (z = 0), they also distribute within a depth of the Rayleigh waves' wavelengths. However, Rayleigh waves not only concentrate on the surface, they also concentrate within a depth of one wavelength. In this section, the beam directivity of Rayleigh waves at various depths are investigated utilising an analytical method.

The Analytical Solution to the Displacement of Rayleigh Waves
N. A. Haskell proposed an analytical solution to Rayleigh waves [22,23]. Ref. [24] by Love introduced these solutions in detail, and investigated the beam directivity of Rayleigh waves at z = 0 with approximated equations. However, Rayleigh waves not only concentrate on the surface; they also concentrate within a depth of one wavelength. The complete equations of the Rayleigh waves' displacement at various depths are: where: where u r and u z are the in-plane and the out-of-plane displacement to be calculated, respectively; r is the distance between the source point and the field point ( Figure 11); F is the excitation source; ρ is the density of the material; θ is the angle between the force vector and the in-plane displacement vector; ω is the operational angular frequency; z is the depth; and c L , c S , c R , and where: (ĸ, , ) = ĸ 2 4 ( 2 2 ĸ 3 ) = cos ( ) (10) where and are the in-plane and the out-of-plane displacement to be calculated, respectively; is the distance between the source point and the field point ( Figure 11); is the excitation source; is the density of the material; is the angle between the force vector and the in-plane displacement vector; is the operational angular frequency; z is the depth; and , , , and ĸ are the velocity of the longitudinal waves, shear waves, Rayleigh waves, and the wave number, respectively.
are the velocity of the longitudinal waves, shear waves, Rayleigh waves, and the wave number, respectively.

Linking EMAT-EM and EMAT-US Models
On the surface plane of the material, the combination between the EMAT-EM model and the EMAT-US model is described in Figure 12. The calculated F values, acting as the excitation source, are imported to each surface layer at various depths. The Rayleigh waves' distribution is the addition of the Rayleigh waves generated by each point source. Table 1 illustrates the parameters used for the EMAT-US model.   Figure 13 shows the calculated Rayleigh waves' radiation pattern at z = 0, which is symmetrical with the center of the coil. Rayleigh waves are mainly generated along the y direction (main lobe) and some undeniable directions (side lobe). The characteristics of Rayleigh waves are quantitatively

Linking EMAT-EM and EMAT-US Models
On the surface plane of the material, the combination between the EMAT-EM model and the EMAT-US model is described in Figure 12. The calculated F values, acting as the excitation source, are imported to each surface layer at various depths. The Rayleigh waves' distribution is the addition of the Rayleigh waves generated by each point source. Table 1 illustrates the parameters used for the EMAT-US model.

Linking EMAT-EM and EMAT-US Models
On the surface plane of the material, the combination between the EMAT-EM model and the EMAT-US model is described in Figure 12. The calculated F values, acting as the excitation source, are imported to each surface layer at various depths. The Rayleigh waves' distribution is the addition of the Rayleigh waves generated by each point source. Table 1 illustrates the parameters used for the EMAT-US model.   Figure 13 shows the calculated Rayleigh waves' radiation pattern at z = 0, which is symmetrical with the center of the coil. Rayleigh waves are mainly generated along the y direction (main lobe) and some undeniable directions (side lobe). The characteristics of Rayleigh waves are quantitatively   Figure 13 shows the calculated Rayleigh waves' radiation pattern at z = 0, which is symmetrical with the center of the coil. Rayleigh waves are mainly generated along the y direction (main lobe) and some undeniable directions (side lobe). The characteristics of Rayleigh waves are quantitatively investigated by means of beam directivity, as shown in the red arc in Figure 13   The beam directivity of Rayleigh waves at z = 0 is shown in Figure 14. A main lobe contains a larger magnitude compared to the side lobes, which contain a smaller magnitude. The main lobe is centered at 0°, and the side lobes are roughly centered at −25.5°, −18°, −10°, 10°, 18°, and 25.5°, respectively. The largest magnitude of the side lobes is 25.87% that of the main lobe. In most applications, side lobes are usually undesirable. The beam directivity of the Rayleigh waves at various depths is shown in Figure 15. The magnitude of the beam directivity is normalised. From this image, the magnitude of the Rayleigh waves decays with the depth, especially for the depth larger than one Rayleigh wavelength. At a depth equalling to one Rayleigh wave's wavelength, z = 6 mm, the magnitude of the Rayleigh waves is 34.9% of that at z = 0. At a depth of 7 mm, the magnitude of the Rayleigh waves decays to 22.8% of that at the surface of the test piece (z = 0). This observation confirms that Rayleigh waves mainly distribute within a depth equal to one Rayleigh wavelength. The beam directivity of Rayleigh waves at z = 0 is shown in Figure 14. A main lobe contains a larger magnitude compared to the side lobes, which contain a smaller magnitude. The main lobe is centered at 0 • , and the side lobes are roughly centered at −25.5 • , −18 • , −10 • , 10 • , 18 • , and 25.5 • , respectively. The largest magnitude of the side lobes is 25.87% that of the main lobe. In most applications, side lobes are usually undesirable. investigated by means of beam directivity, as shown in the red arc in Figure 13 (r = 250 mm, θ1 = −70°, and θ2 = 70°). Figure 13. The radiation pattern of Rayleigh waves generated by the meander-line-coil EMAT. Reproduced with permission from [18], Elsevier, 2017.

Analysis of the Beam Directivity of Rayleigh Waves
The beam directivity of Rayleigh waves at z = 0 is shown in Figure 14. A main lobe contains a larger magnitude compared to the side lobes, which contain a smaller magnitude. The main lobe is centered at 0°, and the side lobes are roughly centered at −25.5°, −18°, −10°, 10°, 18°, and 25.5°, respectively. The largest magnitude of the side lobes is 25.87% that of the main lobe. In most applications, side lobes are usually undesirable. The beam directivity of the Rayleigh waves at various depths is shown in Figure 15. The magnitude of the beam directivity is normalised. From this image, the magnitude of the Rayleigh waves decays with the depth, especially for the depth larger than one Rayleigh wavelength. At a depth equalling to one Rayleigh wave's wavelength, z = 6 mm, the magnitude of the Rayleigh waves is 34.9% of that at z = 0. At a depth of 7 mm, the magnitude of the Rayleigh waves decays to 22.8% of that at the surface of the test piece (z = 0). This observation confirms that Rayleigh waves mainly distribute within a depth equal to one Rayleigh wavelength. The beam directivity of the Rayleigh waves at various depths is shown in Figure 15. The magnitude of the beam directivity is normalised. From this image, the magnitude of the Rayleigh waves decays with the depth, especially for the depth larger than one Rayleigh wavelength. At a depth equalling to one Rayleigh wave's wavelength, z = 6 mm, the magnitude of the Rayleigh waves is 34.9% of that at z = 0. At a depth of 7 mm, the magnitude of the Rayleigh waves decays to 22.8% of that at the surface of the test piece (z = 0). This observation confirms that Rayleigh waves mainly distribute within a depth equal to one Rayleigh wavelength.

Experimental Validations
In Figure 15, based on the analytical simulations, the distribution of the Rayleigh waves at various depths is presented. In this part, the measured results at z = 0 are picked to compare with the simulation results. The experimental setup was the same as the one in Section 2.4. The measured beam directivity at z = 0 is obtained by placing the receiver along the scanning path, as shown in Figure 16. Figure 17 shows the beam directivity results at z = 0 from simulations and experiments. Thirty-three measuring points on the scanning path were picked with a moving step of 2.5°. The measured beam directivity at z = 0 agrees well with the simulated beam directivity, which validates the proposed method. There are some non-overlapping points between these two curves due to several factors, which include the experimental noise and the tolerance of the receiver's position, etc.

Experimental Validations
In Figure 15, based on the analytical simulations, the distribution of the Rayleigh waves at various depths is presented. In this part, the measured results at z = 0 are picked to compare with the simulation results. The experimental setup was the same as the one in Section 2.4. The measured beam directivity at z = 0 is obtained by placing the receiver along the scanning path, as shown in Figure 16. Figure 17 shows the beam directivity results at z = 0 from simulations and experiments. Thirty-three measuring points on the scanning path were picked with a moving step of 2.5 • . The measured beam directivity at z = 0 agrees well with the simulated beam directivity, which validates the proposed method. There are some non-overlapping points between these two curves due to several factors, which include the experimental noise and the tolerance of the receiver's position, etc.

Experimental Validations
In Figure 15, based on the analytical simulations, the distribution of the Rayleigh waves at various depths is presented. In this part, the measured results at z = 0 are picked to compare with the simulation results. The experimental setup was the same as the one in Section 2.4. The measured beam directivity at z = 0 is obtained by placing the receiver along the scanning path, as shown in Figure 16. Figure 17 shows the beam directivity results at z = 0 from simulations and experiments. Thirty-three measuring points on the scanning path were picked with a moving step of 2.5°. The measured beam directivity at z = 0 agrees well with the simulated beam directivity, which validates the proposed method. There are some non-overlapping points between these two curves due to several factors, which include the experimental noise and the tolerance of the receiver's position, etc.

Experimental Validations
In Figure 15, based on the analytical simulations, the distribution of the Rayleigh waves at various depths is presented. In this part, the measured results at z = 0 are picked to compare with the simulation results. The experimental setup was the same as the one in Section 2.4. The measured beam directivity at z = 0 is obtained by placing the receiver along the scanning path, as shown in Figure 16. Figure 17 shows the beam directivity results at z = 0 from simulations and experiments. Thirty-three measuring points on the scanning path were picked with a moving step of 2.5°. The measured beam directivity at z = 0 agrees well with the simulated beam directivity, which validates the proposed method. There are some non-overlapping points between these two curves due to several factors, which include the experimental noise and the tolerance of the receiver's position, etc.