Nondestructive Evaluation of Special Defects Based on Ultrasound Metasurface

We demonstrate the nondestructive evaluation by means of directional ultrasound emitted from a planar metasurface. The ultrasound metasurface is designed to generate the collimated and directional ultrasound efficiently in a planar configuration, which is endowed with the full-2π-range phase manipulation ability and high transmittance up to 80%. We employ the directional emission based on the ultrasound metasurface to innovate the traditional nondestructive evaluation methods, benefited from the freely controlled directivity and the superior fitness to sample surface of the planar metasurface. Merits of this innovative application are evidenced by the remarkable accuracy (higher than 98%) in the thickness evaluation, and precise detection (accuracy higher than 96%) of the special defect inside the V-shaped workpiece which is intractable to be inspected conventionally. The implementation of the metasurface-based directional ultrasound emission in the nondestructive evaluation bears the advantages of high coupling efficiency, superior fitness, high accuracy, and applicability to special defect, providing new solutions to the challenges in conventional defect detection and promotes the development in the nondestructive evaluation applications.


INTRODUCTION
Measurements based on acoustic waves have been widely explored in many applications such as nondestructive evaluation and biomedical engineering. In industrial nondestructive evaluation, the ultrasound measurements (Chatillon et al., 2000;Donskoy et al., 2001;Babich et al., 2004;Dutton et al., 2011a;Edwards et al., 2011;Edwards et al., 2012;Yassin et al., 2018;Taheri and Hassen, 2019) have been proved to be a powerful means for measuring and characterizing the workpieces. By analyzing the incident and reflected/transmitted ultrasound signals from the targeted workpiece, one can effectively characterize and obtain information about the workpiece, such as thickness (Hsu et al., 1994), density, speed of sound (Hans et al., 1996), and position of the defect inside the workpiece (Blackshire and Sathish, 2002;Yu et al., 2018). Ultrasound measurements have also been used in some fields like gas pipeline leakage detection (Liang et al., 2013) and food safety testing (Elvira et al., 2009). For the detection of targets in different scenarios, the collimated and directional ultrasound emission (Ren et al., 2010;Layman et al., 2011;Ren, 2015;Tang and Ren, 2017a;Tong et al., 2020) with varied incident angles is desired, which is conventionally realized by the combination of a planar transducer and the wedges of different angles (Bermes et al., 2008;Pruell et al., 2009). However, there are some limitations of these traditional methods for the directional ultrasound emission in the nondestructive detection (Dutton et al., 2011b). On the one hand, the employment of the wedge would introduce the impedance mismatch and influence the fitness to the sample surface, which would bring about the reduction in the transmission efficiency. On the other hand, for some special sample measurements such as the V-shaped workpiece, it might meet difficulties in placing the equipment for inspecting of the buried defect and thus affect the defection accuracy.
The emergence of acoustic metamaterial and metasurface (Chen and Chan, 2007;Chen and Chan, 2010;Zhao et al., 2013;Al Jahdali and Wu, 2016;Zhu et al., 2017;Chu et al., 2018;Zhao et al., 2019;Zhu and Assouar, 2019) could provide a new avenue to meet these challenges in nondestructive evaluation. One of the uniqueness of the metasurface lies in its ability to freely manipulate the wavefront by freely controlling the phase and amplitude of acoustic wave (Assouar et al., 2018). Benefiting from the special properties and subwavelength geometry of the metasurface, a variety of exotic acoustic phenomena, including acoustic cloaking (He et al., 2020), acoustic illusion (Shen et al., 2019a), super-resolution focusing Qian et al., 2017;Chen et al., 2018;Shen et al., 2019b), selfaccelerating beams (Li and Assouar, 2015), acoustic vortex with twisted wavefront (Jiang et al., 2016;Ye et al., 2016), non-diffraction beams (Durnin et al., 1987), and asymmetric transmission (Li et al., 2017), have been successfully realized with acoustic metasurface. The efficient wave control ability of the metasurface also enables the directional and collimated emission in an ultrathin planar configuration, which could be employed to promote the applications in the traditional nondestructive evaluation and biomedical engineering. The directivity can be improved (Tang and Ren, 2017a;Tang and Ren, 2017b), and the emission direction can be controlled by adjusting the phase profile Xia et al., 2018) on the metasurface based on the generalized Snell's law (Yu et al., 2011).
In this work, we employ the acoustic metasurface as a powerful instrument to innovate the nondestructive evaluation, evidenced by the precise thickness evaluation and special defect inspection unavailable in the conventional methods. These are achieved by the collimated and directional ultrasound emission with the acoustic metasurface processing both the full-2π-range phase manipulation and the high transmittance (up to 80%). Based on the directional ultrasound emission, we numerically demonstrate the thickness evaluation and defect detection inside the special V-shaped workpiece, with remarkable accuracies higher than 98% and 96%, respectively. With the advantages of high coupling efficiency, superior fitness to sample surface, high accuracy, and applicability to special defect, the directional ultrasound emission based on the proposed planar metasurface would offer an alternative route to ultrasound measurements and boost the development in nondestructive evaluation.

Multi-cavity element on the metasurface
We realize the directional ultrasound emission by designing the metasurface consisting of six identical cascaded elements. The schematic of an individual element is shown in Figure 1A, which is the combination of a rectangular cavity located in the center of the element and two symmetric inclined channels on both sides of the rectangular cavity. The structure is made of stainless steel, whose mass density, sound speed, Young's modules, and Poisson's ratio are ρ s 7800 kg/m 3 , c s 6100 m/s, E 215 GPa, and σ 0.3, respectively, and the background medium is water (ρ w 1000 kg/m 3 , c w 1500 m/s). The acoustic impedance of stainless steel is approximately 31 times higher than that of water. The lengths (heights) of the element and cavity are t (h) and t 1 (h 1 ), respectively, and the height and inclined angle of the two channels are h 2 and θ. All these parameters contribute to the phase control of the transmitted wave. By analyzing the ultrasound response dependence on these parameters, we choose the parameter t 1 as a knob to manipulate the phase delay since the phase of the transmitted ultrasound is found to be relatively sensitive to t 1 . As a result, this single-parameterbased control makes it much easier in designing the ultrasound metasurface. The metasurface is designed for the underwater ultrasound at 300 kHz, with the wavelength λ 5 mm and other parameters fixed to be t 1.4 mm, h 1.25 mm, h 1 1.1 mm, h 2 0.1 mm, and θ 45°. The transmission performance of the metasurface is simulated by employing the acoustic-structure interface, frequency domain with the finite element method based on COMSOL Multiphysics software. Figure 1B shows the phase delay (red solid line) and the transmittance (blue solid line) of an element as functions of parameter m (m t 1 /t), where t 1 is varied in the range from 0.7 mm to 1.3 mm. It can be observed that the phase delay covers a 0.4π range in this parameter region where the transmittance remains higher than 82%, which guarantees the high efficiency in the following nondestructive evaluation applications.
For achieving the full-2π-phase control, we connect six identical elements in series in the x-direction to construct a single unit, with the periodicity of the elements being a 2.3 mm. The schematic of a unit which is constructed by six elements is illustrated in Figure 1C. The ultrasound incident from the left accumulates the phase delay after passing through the structure, while keeping the transmittance higher than 80%, as shown in Figure 1D.

Design of the metasurfaces for directional ultrasound emission
The units consisting of six cascaded elements are arranged in the y-direction according to the generalized Snell's law to design the ultrasound metasurface for directional emission. Considering the ultrasound wave normally incident from the x-direction, the transmitted ultrasound with the refraction angle θ t (y) can be realized by the metasurface with the phase distribution: where φ(y) represents the phase distribution on the metasurface, k 2πf/c is the wave vector in water, and f 300 kHz is the frequency of the incident ultrasound. The continuous distribution of the phase profile is discretized to facilitate the implementation with the metasurface, and the discretization resolution is determined by the unit height (h 1.25 mm) on the metasurface. Two sets of directional ultrasound emission with the refraction angles of 30°and 45°are realized by constructing two different planar metasurfaces, both of which has 80 units. We numerically simulate the directional ultrasound emission based on COMSOL Multiphysics software. The distributions of the directional ultrasound fields with the refraction angles of 30°F  Frontiers in Materials | www.frontiersin.org January 2022 | Volume 8 | Article 802001 3 and 45°are illustrated in Figures 2A,B. The normally incident plane wave propagates along the x-direction and impinges on the metasurface, whose phase distribution is tuned by the metasurface. As a result, the transmitted ultrasound wave is endowed with the additional transverse momentum along the y-direction and passes through the metasurface with the targeted refraction angle. In Figures 2A,B, it can be observed that the incident ultrasound wave is deflected by 30°and 45°after passing through the metasurfaces, respectively, agreeing well with the targeted directions. In addition, the amplitude of the transmitted wave remains relatively high compared with the incident wave, which proves the high transmission of the metasurface and contributes to its application in nondestructive evaluation.

RESULTS
The directional and collimated ultrasound emission based on the metasurface is employed in the thickness evaluation and defect detection inside the special V-shaped workpiece and provides an alternative route for nondestructive evaluation by taking advantage of the unique advantages of the high coupling efficiency, superior fitness, high accuracy, and applicability to special defect, which will be shown as follows.

Thickness evaluation
We first demonstrate the thickness evaluation of a plexiglass block underwater with the directional ultrasound emission of the refraction angle θ t 45°. The density and sound speed of plexiglass block are ρ s 1180 kg/m 3 and c s 2730 m/s, respectively. A planar metasurface consisting of 100 elements are constructed according to the phase profile, with the total height of H 125 mm.
Owing to the planar configuration, the designed metasurface can be directly attached to the surface of the detected workpiece which helps to reduce the mismatch caused by the wedge in conventional methods and improve the efficiency of the evaluation. A plane ultrasound at 300 kHz is normally incident from the left and transmits through the metasurface. As shown in the ultrasound intensity field in Figure 3A, ultrasound wave can almost completely penetrate the plexiglass block (indicated by the red solid box), and the transmitted beam is deflected to the direction of 45°by the metasurface (an enlarged view of it is shown in the green dashed box). The ultrasound beam transmitted from the metasurface is directionally emitted to the plexiglass block from the left surface and then further reflected by the right surface of the plexiglass block and detected on the left surface. The position of maximum ultrasound intensity on the left surface of the plexiglass block reveals the geometrical relationship between the emitted and reflected ultrasound trajectories and can be used to evaluate the thickness of the target. The trajectories of the directional ultrasound beam emitted from the metasurface and reflected from the right surface of the plexiglass block are illustrated in Figure 3B. The distance from the bottom of the metasurface on the left surface of the plexiglass block is denoted as l 1 . According to the geometric relationship ( Figure 3B), the thickness of the block can be calculated by: where l 1_max is the position of the maximum ultrasound intensity.
An ultrasound intensity profile as a function of l 1 is plotted in Figure 3C, which shows that the intensity reaches the maximum value at position l 1 max 343 mm. The evaluated thickness of the block is then calculated to be d 1 202.75 mm by Eq. 2, which has a good agreement with the actual value d a1 200 mm, with the relative error ε 1 (d 1 − d a1 )/d a1 1.375% smaller than 2%.
To compare the thickness evaluation performance based on ultrasound metasurface for the workpiece with different thicknesses, we conduct the evaluation of the plexiglass blocks with the thickness d a1 ranging from 100 to 250 mm. The evaluation accuracy as a function of the thickness d a1 is plotted in Figure 3D, where it is consistently higher than 98% over this thickness range. The consistency between the evaluated and actual thicknesses strongly proves the high accuracy and effectiveness of the implementation of the directional ultrasound emission based on metasurface in thickness evaluation, which provides broad interests on the applications in industrial engineering.

Defect detection inside the V-shaped workpiece
We further consider the defect detection inside the object with the metasurface-based directional emission. As an example, we demonstrate the testing of the defect inside a V-shaped workpiece, which is widely used in industry while quite challenging to detect for a conventional wedge-based method due to the special shape. Figure 4A shows the schematic of a V-shaped workpiece made of plexiglass with a strip of air defect inside it. We use the constructed metasurface which can emit the directional ultrasound beam to detect the strip defects inside the workpiece. Figure 4B shows the two-dimensional cross section of defect detection inside the V-shaped workpiece. The metasurface (indicated by the green dashed box) and the V-shaped workpiece (indicated by the red solid box) are immersed in water, and the planar metasurface is closely adherent to the surface of the workpiece. A plane ultrasound wave impinges on the upper surface of the metasurface and is transmitted through the metasurface with the refraction angle of 45°. The transmitted ultrasound beam is then emitted into the V-shaped workpiece and interacts with the defect inside, which leads to the reflection by the defect due to the acoustic impedance mismatch between the plexiglass and air stripe, and the reflected signal is detected on the top-right surface of the workpiece. Similarly, the position of the defect can be found by analyzing the trajectories of the emitted and reflected beams, and the geometric relationship of the V-shaped workpiece and the strip defect ( Figure 4C).
Specifically, the position of the maximum ultrasound intensity on the top-right surface l 2_max indicates the position of the defect inside the workpiece. The V-shaped workpiece tested here has the included angle of α 140°, and the distance between the workpiece concave corner and center of the metasurface is r 73.14 mm. The distance from the concave corner of the V-shaped workpiece along the top-right surface is denoted as l 2 . The ultrasound intensity profile along the top-right surface as a function of l 2 is plotted in Figure 4D, where the maximum ultrasound intensity is obtained at the position l 2 max 75 mm. Therefore, the height of the strip defect from the concave corner of the V-shaped workpiece can be calculated as which is calculated to be d 2 123.93 mm by substituting the l 2_max value. The actual position of strip defects is d a2 120 mm, and the relative error of the detection is ε 2 (d 2 − d a2 )/d a2 3.3%, which remains smaller than 4%. These results reveal that the directional ultrasound emission based on the metasurface can be expediently employed to the special defect detection with high precision. The defect with the sub-wavelength size could also be detected by integrating the directional emission with other technologies, such as the ultrasound focusing, beam forming, and artificial intelligent algorithm.

DISCUSSION
In conclusion, we have theoretically proposed and numerically demonstrated the nondestructive evaluation including the thickness evaluation and special defect detection by means of directional ultrasound waves emitted from a planar metasurface. The metasurface composed of parallel multi-cavity cascaded elements has been designed, which ensures a flexible 2π-range phase control while maintaining the transmittance up to 80%. The thickness evaluation with the accuracy higher than 98% has been achieved, and the comparison of the performance for evaluating different thicknesses has been discussed. Furthermore, the precise inspection of the defect located inside the V-shaped workpiece unavailable with conventional methods has been efficiently realized with the accuracy higher than 96%. The nondestructive evaluation with the metasurface-based directional ultrasound emission, with the merits of high coupling efficiency, superior fitness to sample surface, high accuracy, and applicability to special defect, would provide additional solutions to the challenges in conventional defect detection and promote the development in nondestructive evaluation applications.

DATA AVAILABILITY STATEMENT
The raw data supporting the conclusion of this article will be made available by the authors, without undue reservation.

AUTHOR CONTRIBUTIONS
HZ and XJ conceived the research. XJ and DT supervised the project. HZ performed the theoretical calculation and analysis. XJ, JH, and CZ contributed to perform the analysis with constructive discussions. All authors contributed to manuscript revision and read and approved the submitted version.