Next Article in Journal
Construction of Superhydrophobic Coating on Iron Surface with Enhanced Anti-Corrosion, Anti-Adhesive and Anti-Bacterial Properties
Previous Article in Journal
Effect of Basalt Powder on Hydration, Rheology, and Strength Development of Cement Paste
Previous Article in Special Issue
Ultrasound Study of Magnetic and Non-Magnetic Nanoparticle Agglomeration in High Viscous Media
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Finite Element Simulation of Ultrasonic Scattering by Rough Flaws with Multi-Scale Distortions

1
Beijing Institute of Aeronautical Materials, Beijing 100095, China
2
School of Traffic and Transportation Engineering, Central South University, Changsha 410075, China
*
Author to whom correspondence should be addressed.
Materials 2022, 15(23), 8633; https://doi.org/10.3390/ma15238633
Submission received: 18 October 2022 / Revised: 17 November 2022 / Accepted: 28 November 2022 / Published: 3 December 2022

Abstract

:
The roughness of a flaw’s surface significantly affects the scattering behavior of ultrasonic waves. It is vital to understand the impact of roughness on flaw echoes, especially when performing ultrasonic nondestructive inspection on safety-critical components. However, the current approach for creating rough flaw models fails to reconstruct complicated cracks with secondary cracks. Here, a multi-scale distortion method is developed to generate a rough flaw by using an optical microscope image of a real flaw. The finite element (FE) is then implemented to simulate the near-surface rough flaws in nickel-based bars, which are detected by an offsetting immersion transducer with mode-converted transverse waves. Numerical results show that the randomness and complexity of flaw echoes from rough flaws are exceptionally high. The gap between the maximum and minimum normalized amplitude values of flaw echoes from a rough crack with secondary cracks can reach 7.125 dB. Meanwhile, the maximum time of flight (TOF) is almost twice as large as the minimum TOF. Therefore, the present method can generate effective rough flaw models in terms of macroscopic rough geometry and microscopic rough surface. Moreover, the impact of the rough flaw surface on the flaw echoes goes beyond amplitude changes and may make flaw location challenging.

1. Introduction

Ultrasonic nondestructive inspection is dependent on the scattering of elastic waves by flaws. Hence, the scattering behavior of flaws must be investigated to improve the inspection scheme. Schmerr et al. [1] highlighted that the flaw scattering amplitude is strongly connected to the flaw echo signal and plays a crucial role in flaw sizing. However, real flaws are generally not standard reflectors with regular shapes (e.g., notch, void, and column), and changes in macroscopic rough geometry will alter the fundamental characteristics of flaw scattering. Meanwhile, real flaws do not have an ideally smooth surface. In addition, scattering caused by microscopic rough surfaces can significantly degrade signal amplitude [2,3,4] and should not be ignored. Thus, the impact of the surface roughness of flaws on ultrasonic inspection and quantitative evaluation should be analyzed at multiple scales: from macro to micro.
Most studies on elastic wave scattering from a rough surface rely on computer simulations, with rough surface reconstruction serving as the most important part. However, the development of rough surface reconstruction is currently primarily focused on microscopic rough surfaces. For example, Ogilvy [5] used the moving average process, the Gaussian height distribution, and the Gaussian surface correlation function (SCF) to reconstruct the rough surface. Monte Carlo simulations were then performed to obtain the statistical features of ultrasonic scattering. However, the point number of moving average cannot be excessively large or excessively small. Zhang et al. [6] and Pettit et al. [7] also employed the Gaussian SCF to reconstruct a one-dimensional (1D) rough surface, where the convolution operation for the Gaussian white noise was emphasized. Monte Carlo finite element (FE) simulations were conducted using the 1D rough surface. Their goal was to verify the analytical models using the Kirchhoff approximation. However, in reality, there is a significant difference between a 1D rough surface and three-dimensional (3D) rough surface. Jarvis and Cegla [8,9] also reconstructed a 1D rough surface with a SCF-based method, but the scattering events were simulated by the distributed point source method (DPSM).
Shi et al. [10,11] used the SCF-based method to reconstruct a two-dimensional (2D) rough surface. They utilized the 3D FE boundary integration approach to analyze the rough-surface scattering. Recently, Choi and Shi et al. [12] enhanced the SCF-based method by moving average and auto-regressive (AR) methods. They found that the AR surfaces showed the best agreement with the real rough surfaces. Even though the SCF-based approach has made significant progress in numerical simulations for microscopic surface roughness, it cannot simultaneously account for the macroscopic rough geometry in surface modeling. In addition, the Weierstrass-Mandelbrot (W-M) function [13] can be used to reconstruct a fractal-regular rough surface, but the W-M function cannot be used to create a complex rough surface with a primary profile. Obviously, the existing approaches are inadequate for simulating a primary crack with secondary cracks, because it is difficult to extract and reset the profile of a crack surface. Consequently, this work aims to develop a universal method for surface modeling at multiple scales, which will lead to a more realistic ultrasonic FE simulation.
In this work, we introduce an image distortion technique [14,15] into the ultrasonic FE simulation. This is a random perturbation approach that can produce rough-surfaced flaws with multi-scale distortions, according to an optical microscope (OM) image of a real flaw. The FE simulations were placed in a scenario where an offsetting immersion transducer used mode-converted transverse waves to detect near-surface flaws in nickel-based alloy bars. Finally, the simulation results of the flaw echoes from secondary cracks, sub-surface pores, and sub-surface inclusion were statistically analyzed and the effects of flaw scattering behavior with various rough surfaces were discovered.

2. Materials and Methods

2.1. Materials

Nickel-based alloy bars are widely used in the aerospace industry, due to their superior high-temperature strength, oxidation resistance, corrosion resistance, and fatigue resistance [16]. However, flaws are inevitably produced in nickel-based bars during the manufacturing process, which affects their performance and service safety. For instance, cracks, pores, and inclusions are frequently observed on the surface and subsurface of these bars [17,18,19,20]. Over the history of our extensive inspection experience, we have seldom found smooth flaws in nickel-based alloy bars, and the accurate ultrasonic estimation of flaw size is complicated by the presence of rough surfaces. Thus, nickel-based alloy bars are selected for the ultrasonic FE simulation in this work. Table 1 shows the input material properties that are associated with the nickel-based alloy, and one common inclusion of this alloy.

2.2. Generating Rough Flaws with Multi-Scale Distortions

Image distortion is defined as the aberrant or warped representation of an object by an optical imaging system. This phenomenon must be suppressed and corrected in the field of optical imaging. Figure 1 illustrates the pincushion, barrel, and ripple distortions of a 3D cube and a 2D checkerboard, where the pincushion distortion and barrel distortion phenomena originate from a lens’s misaligned edge-to-center magnification, and the ripple distortion phenomenon originates from the imaging difficulties caused by water-surface ripples. However, we suggest that the image distortion technique could be employed to artificially create rough flaws into nickel-based alloy bars. It is possible to generate distorted images at multiple scales by randomly perturbing the seed of a flaw image with different distortion parameters. This paves the way for an accurate simulation of the ultrasonic response characteristics of rough flaws, which will be described in the next section.
First, image distortion is essentially coordinate mapping, which can be achieved by converting the coordinates according to various mapping functions, and assigning the values under the original coordinates to the new coordinates, as follows:
x d = x + δ x x , y , z , y d = y + δ y x , y , z , z d = z + δ z x , y , z ,
where x , y and z are the original coordinates; x d , y d and z d are the new coordinates after distortion; and δ x , δ y and δ z are distortions in three directions. The mathematical definition of ripple distortion, pincushion distortion, and barrel distortion can usually be expressed as [14]:
ripple distortion:
δ x = K cos 2 π Τ y + φ cos 2 π Τ z + φ , δ y = K cos 2 π Τ z + φ cos 2 π Τ x + φ , δ z = K cos 2 π Τ x + φ cos 2 π Τ y + φ ,
pincushion distortion and barrel distortion:
δ x = sgn x K cos 2 π Τ y 2 + z 2 + φ , δ y = sgn y K cos 2 π Τ z 2 + x 2 + φ , δ z = sgn z K cos 2 π Τ x 2 + y 2 + φ ,
where K , T and φ are the coefficients of image distortion; K is the distortion amplitude; T is the period term; and φ is the phase term. It should be noted that a positive K value results in pincushion distortion, while a negative K value results in barrel distortion. When the z coordinate is unavailable, the 3D distortion degrades to a 2D distortion (let z = 0 ). It is difficult to count the flaw response under the independent identically distributed (IID) condition, because pincushion distortion and barrel distortion tend to cause significant changes in the flaw’s size. Therefore, we chose the ripple distortion technique for multi-scale perturbation of the seed of the flaw image.
As shown in the flowchart in Figure 2, the original flaw image is preprocessed and progressively distorted according to specified distortion parameters. For macroscopic distortion, the random distortion amplitude is K i U K 1 i , K 2 i , the random period is T i U T 1 i , T 2 i and the random phase is φ U 0 , 2 π , where U · denotes the uniform distribution; K 1 i < K 2 i and; K 1 , 2 i > K 1 , 2 i + 1 and T 1 , 2 i > T 1 , 2 i + 1 ; i = 0 , 1 , , N 1 ; and N is the number of macroscopic distortions. For microscopic distortion, the random distortion amplitude is K j U K 3 j , K 4 j , the random period is T j U T 3 j , T 4 j and the random phase is also φ U 0 , 2 π , where K 3 j < K 4 j and T 3 j < T 4 j ; K 3 , 4 j > K 3 , 4 j + 1 and T 3 , 4 j > T 3 , 4 j + 1 ; j = 0 , 1 , , M 1 ; and M is the number of microscopic distortions. K 1 , 2 , K 3 , 4 , T 1 , 2 , and T 3 , 4 are the hyper-parameters of the uniform distributions. Notice that K 3 , 4 K 1 , 2 and T 3 , 4 T 1 , 2 reflect the division between macroscopic and microscopic scales. Continuous distortions can be achieved by nested calculations of Equations (1) and (2). Considering Equation (2) as part of the Fourier series, we may use the idea of spatial Fourier transform to help understand our method. In the spatial Fourier transform, the convolution calculation and transfer function can change the input image. Here, the amplitude, period, and phase parameters of the ripple distortion can be related to the module, wave number, and direction of the wave vector of the transfer function, respectively. Thus, any image transformation might be achieved by combining the appropriate parameters of the ripple distortion. In this work, we have randomly chosen a few discrete parameter combinations as samples. However, it is not necessary to always use uniform distribution; one may obtain the desired rough surface reconstruction by carefully designing the parameters combination for ripple distortion based on experimental observations in the future. Finally, repeating this process Q times results in Q flaws with various degrees of roughness.
Figure 3 shows three samples of seeds and the resulting multi-scale ripple distortion after multiple random perturbations. Let Q = 100 , and three typical rough flaws generated within 100 repetitions are exhibited for each seed. Figure 3a depicts a primary crack with several secondary cracks, which has seldom been discussed in ultrasonic simulation. This crack is equivalently ~0.02 mm wide and ~3 mm deep in the original OM images. Figure 3b,c demonstrate the multi-scale distortion results of a pore and a non-metallic inclusion, whose equivalent diameters in the original OM images are 0.1 mm and 1 mm, respectively. As observed in the pseudo-color images of Figure 3, flaws distorted by random ripples are notably more different in macroscopic rough geometry and microscopic rough surfaces than the seeds of the original flaw images.
Additionally, based on ISO 21920-2:2021, the arithmetic mean roughness R a can be expressed as:
R a = 1 0 Z x d x
where Z and are the height and length of the roughness profile. To produce the roughness profile, the flaw’s primary profile or waviness profile must be eliminated. Thus, the edge extraction, moving average, and skeletonization methods were employed to calculate the primary profile, such that the roughness profile can be obtained by deviating from the primary profile in the normal direction. Consequently, the surface roughness R a of 100 cracks, pores, and inclusions are: 10.85 ± 2.85 μm, 0.60 ± 0.37 μm, and 9.05 ± 3.15 μm, respectively. The roughness of the pore is one order of magnitude smaller than that of the other two flaws because the size of the pore is also one order of magnitude smaller.

2.3. Ultrasonic Finite Element Method

The ultrasonic inspection of nickel-based bars with a mode-converted transverse wave is simulated and analyzed with the discontinuous Galerkin time-domain (DGTD) method, which solves the sound field and pulse-echo signals from a focal transducer. Here, Q = 100 is the total number of simulation repetitions for each case. In each successive simulation, the rough defect model is altered by employing the multi-scale distortion technique. In addition, the attenuation coefficient of the nickel-based bar is considered to predict the flaw scattering amplitude with high accuracy. This simulation is built in COMSOL® and MATLAB® using an Intel CPU i9-10900K (10 cores, 5 GHz). The simulations will be verified by comparing them to the experimental results in Section 3.

2.3.1. FE Simulation Setup

Figure 4 shows the FE simulation setup for the ultrasonic inspection of nickel-based bars with a mode-converted transverse wave, which was simplified to a 2D model under the assumption of plane strain. The focal transducer was simulated using a concave circular arc, with an element radius of 3.175 mm, a focal length of 12.7 mm, and a water path of 8.7 mm. The radius of the nickel-based bar was set to 10 mm. According to Table 1, the corresponding longitudinal wave velocity of the nickel-based bar was 5640 m/s, its transverse wave velocity was 3067 m/s, and the longitudinal wave velocity in water was 1486 m/s. Based on Snell’s law, the first critical angle of the fluid/solid interface was 15.28°, and the second critical angle was 28.98°. To regulate the refraction angle to 45°, the incidence angle should be 21.77°. Hence, the transducer’s offset distance was adjusted to 3.7 mm. As a baseline for further simulations of rough flaws, an ideal open crack or notch was analyzed at first. The width of the notch was 0.2 mm and its depth was 3 mm, according to the real crack in Figure 3a. A perfectly matched layer of sound absorption is applied to reduce sidewall interference and boundary reflections. Considering that we must perform a large number of repeated simulations for a class of flaws, the center frequency of the transducer was set to 2.5 MHz to obtain simulation results within a certain computational time.
To make the simulation more realistic, material attenuation is added to COMSOL’s solid mechanics and acoustics modules to quantify the attenuation in solids and in water; their values are given in Figure 5. The frequency-dependent absorbing attenuation coefficient in water is α w a t e r dB / m = 0.22 f 2 MHz , which is 1.36 dB/m at 2.5 MHz. The longitudinal and transverse wave attenuation coefficients in solids can be obtained from Weaver’s scattering attenuation model [23]. If the average grain diameter of the nickel-based bar is 50 µm, then the scattering attenuation coefficients of the longitudinal and transverse waves are 1.33 and 6.10 dB/m, respectively.
Figure 6 gives the typical FEM results of flaw scattering from a smooth notch in the nickel-based bar. From the snapshots, we can see how the incident wave impacts the bar and converts into a transverse wave. Then, as the refracted wave travels to the smooth notch, the flaw scatters both the longitudinal and transverse wave components. Subsequently, the flaw scattering waves return to the solid/fluid interface. Finally, the flaw echo is received by the focal transducer. The imaging of snapshots is governed by the second stress invariant, with red representing the longitudinal wave component and blue representing the transverse wave component. The results that are presented in Figure 6 are consistent with the expectations.

2.3.2. Model Verification

To demonstrate the reliability of the rough flaw simulations, ultrasonic experiments with a standard smooth notch were conducted to verify the simulations. A roller mechanism and a 2.5 MHz line focal transducer were utilized to perform a circular B-scan on a segment of nickel-based bar with a 10 mm radius and a 20 mm height, using a 1° angular scanning step. To be consistent with the simulation shown earlier, the offset distance was adjusted to 3.7 mm. This indicates that the incident angle is 21.77° and the corresponding refracted wave is 50°. In the radial direction, an open notch with a width of 0.2 mm and a depth of 3 mm was made into the nickel-based bar using low-speed wire cutting, and the roughness of the notch was 0.67 ± 0.18 μm.
The B-scan images from the simulation and the experiment are shown in Figure 7a and Figure 7b, respectively. There is substantial agreement between these two results. Because the refracted transverse wave is not always perpendicular to the notch, the time of flight (TOF) of the flaw echo is not a linear function of the scanning angle, as demonstrated in Figure 7 by the blue short-dot line. Here, TOF is the time difference between the front-wall echo and flaw echo on the time axis. To compare the simulated and measured TOFs of the notch, we recorded TOF values at a 45° scanning angle, which are marked by the upper triangles in Figure 7. In Figure 7a,b, the TOFs measure 2.52 and 2.64 μs, which indicates that their relative error is just 2.7%. Furthermore, Figure 8 shows the normalized flaw echoes corresponding to the markers. There is no obvious discrepancy between the simulated and measured flaw echoes. The discrepancy in the front-wall echoes might be due to the focus line not being on the surface, and the 2D simulation and 3D experiment at the non-focal point are different. Therefore, in subsequent simulations, a smooth notch will be replaced with various types of rough flaws to examine the influence of rough surfaces on flaw echoes.

3. Results and Discussion

In this section, ultrasonic inspections for various types of rough defects will be simulated based on the validated FE model (which was described earlier). The impact of the rough flaws on the flaw echoes will be explored using the crack, pore, and inclusion after multi-scale distortion, as shown in Figure 3. The sizes of all three flaws have previously been mentioned. However, the starting point of each rough crack lies on the bar’s surface, and the line between the starting point and centroid point of each rough crack passes through the center of the bar, while the pores and inclusions were positioned 1 mm beneath the surface. The other FE simulation setup is identical.
Figure 9 shows the statistical findings of 100 groups of different flaw echoes in the time and frequency domains for a rough crack. By use of the Hilbert transform, the signal envelope is obtained from the time-domain waveform signal, which is then normalized by the maximum front-wall echo. Meanwhile, the corresponding frequency-domain signal for the flaw echo is acquired in 14 to 16 µs through fast Fourier transform. Figure 9a shows that the randomness and complexity of flaw echoes are exceptionally high. Figure 9b shows that the peak frequencies occur randomly in the range of 2.5 ± 1.25 MHz. However, there may have been nonlinear flaw scattering processes, such as the generation of subharmonics. In Figure 9c,d, the dark red line represents the ensemble average, whereas the red region represents the range of one standard deviation. The amplitude of the ensemble average of flaw echoes from the rough crack drops by −7.13 dB when compared with the flaw waves of the same-sized smooth notch (0.02 mm width and 3 mm depth). The flaw echo of rough flaws is less than that of the standard reflector, which is in accordance with previous experience and experimental results. In addition, the variation in the rough defect wave at this frequency is more noticeable because the center frequency of the transducer is 2.5 MHz.
The joint distribution of the amplitude and TOF of the peak of flaw echoes are shown in Figure 10. The maximum and minimum values of the normalized amplitude are 0.91 and 0.40, and their difference is greater than double. However, the maximum amplitude of the rough crack is still less than that of the smooth notch. The maximum and minimum values of the TOF are 3.48 and 1.73 µs, and the gap between them is about double. The peak of the TOF distribution is close to the TOF of the smooth notch. These results show that after multi-scale distortion, the flaw scattering of rough cracks becomes unexpected. This makes it challenging to locate and quantify them.
Figure 11 shows the statistical results for the rough pore and rough inclusion. From the range of one standard deviation, the fluctuations of rough pore and rough inclusion are 1.22 dB and 1.93 dB, respectively. This may happen because their equivalent sizes are all smaller than the wavelength of the transverse wave, with the pore size being 0.08 times the wavelength and the inclusion size being 0.81 times the wavelength, whereas the wavelength of the transverse wave at 2.5 MHz is 1.23 mm in the nickel-based bar. Therefore, a 2.5 MHz transducer will be unsuitable if small flaws must be detected during the actual inspection. We anticipate that the roughness-induced fluctuation will increase when increasing the inspection frequency or shortening the wavelength of the transverse wave. Additionally, it is interesting to observe that after multi-scale distortion, both the flaw echoes of rough pore and rough inclusion are higher than those of the equivalent-sized standard reflector. More specifically, their leading-edge responses are noticeably stronger than those caused by circular flaws. This may be caused by the macroscopic rough geometry of the distorted flaw. This requires a more detailed study in the future.

4. Conclusions

In this work, the image distortion technique from the optical imaging field is first introduced for the generation of a rough flaw for the FE simulation of ultrasonic flaw scattering. The image distortion technique was briefly reviewed, and a multi-scale ripple distortion method for operating the seed of the OM image of a real flaw was developed, which resulted in the generation of rough flaws that have a feeling of realism. The proposed method has demonstrated its superiority in generating rough cracks with secondary cracks, rough pores, and rough inclusions. We then selected the ultrasonic inspection of surface or sub-surface flaws in nickel-based bars, using a mode-converted transverse wave as an example. To perform the 2D FE simulations, the rough flaw was positioned on the surface or sub-surface of the bar. The numerical results show that the macroscopic rough geometry and the microscopic rough surface of the flaw may be to blame for the randomness of amplitude and TOF in the flaw echoes. Furthermore, our simulations revealed that it would be extremely challenging to quantify the size of a rough crack with secondary cracks.
Finally, it should be noted that the present work has observed several limitations. First, the parameter combination settings of ripple distortion that are presented here are for illustrative purposes only. The use of uniform distribution has not been experimentally verified and is not mandatory. Future industrial applications require research into the explainable parameter combination of ripple distortion, which will be useful for controlling the correlation length and the arithmetic mean roughness. Second, limited by computing resources, our simulation frequency is lower than the frequency of real inspection. Third, only 2D simulations are performed, but 3D simulations for rough crack, rough pore, and rough inclusion are necessary to compare the experiments. Finally, experimental verification for rough flaw inspection is unavailable. These limitations also provide us insight into where to focus our efforts in the future. In addition, the practical applications may extend beyond the field of ultrasonic nondestructive evaluation to include medical ultrasound diagnosis, quantitative seismology, and seabed sediment acoustics.

Author Contributions

Data curation, Z.W. and Z.Z.; Software, Z.W. and Z.Z.; Validation, Z.Z.; Investigation, Z.Z.; Writing—original draft preparation, Z.Z.; Resources, Z.W.; Methodology, Y.S.; conceptualization, Y.S.; Writing—review and editing, Y.S.; Writing—review and editing, X.L.; Supervision, X.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China, grant number 92060111; the National Equipment Pre-research Foundation of China, grant number 6141C56; and the collaborative research project between Central South University and AECC Beijing Institute of Aeronautical Materials, contract number 2021-07-2227-WX.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Due to the nature of this research, since the participants of this study did not agree for their data to be shared publicly, the supporting data are not available.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Schmerr, L.W. Fundamentals of Ultrasonic Nondestructive Evaluation, 2nd ed.; Springer: Cham, Switzerland, 2016; pp. 419–523. [Google Scholar]
  2. Ogilvy, J.A. Theoretical Comparison of Ultrasonic Signal Amplitudes from Smooth and Rough Defects. NDT Int. 1986, 19, 371–385. [Google Scholar] [CrossRef]
  3. Ogilvy, J.A.; Culverwell, I.D. Elastic Model for Simulating Ultrasonic Inspection of Smooth and Rough Defects. Ultrasonics 1991, 29, 490–496. [Google Scholar] [CrossRef]
  4. Ogilvy, J.A. Model for the Ultrasonic Inspection of Rough Defects. Ultrasonics 1989, 27, 69–79. [Google Scholar] [CrossRef]
  5. Ogilvy, J.A. Computer Simulation of Acoustic Wave Scattering from Rough Surfaces. J. Phys. Appl. Phys. 1988, 21, 260–277. [Google Scholar] [CrossRef]
  6. Zhang, J.; Drinkwater, B.W.; Wilcox, P.D. Longitudinal Wave Scattering from Rough Crack-like Defects. IEEE Trans. Ultrason. Ferroelectr. Freq. Control 2011, 58, 2171–2180. [Google Scholar] [CrossRef]
  7. Pettit, J.R.; Walker, A.E.; Lowe, M.J.S. Improved Detection of Rough Defects for Ultrasonic Nondestructive Evaluation Inspections Based on Finite Element Modeling of Elastic Wave Scattering. IEEE Trans. Ultrason. Ferroelectr. Freq. Control 2015, 62, 1797–1808. [Google Scholar] [CrossRef] [Green Version]
  8. Jarvis, A.J.C.; Cegla, F.B. Application of the Distributed Point Source Method to Rough Surface Scattering and Ultrasonic Wall Thickness Measurement. J. Acoust. Soc. Am. 2012, 132, 1325. [Google Scholar] [CrossRef]
  9. Jarvis, A.J.C.; Cegla, F.B. Scattering of near Normal Incidence SH Waves by Sinusoidal and Rough Surfaces in 3-D: Comparison to the Scalar Wave Approximation. IEEE Trans. Ultrason. Ferroelectr. Freq. Control 2014, 61, 1179–1190. [Google Scholar] [CrossRef]
  10. Shi, F.; Lowe, M.J.S.; Xi, X.; Craster, R.V. Diffuse Scattered Field of Elastic Waves from Randomly Rough Surfaces Using an Analytical Kirchhoff Theory. J. Mech. Phys. Solids 2016, 92, 260–277. [Google Scholar] [CrossRef] [Green Version]
  11. Shi, F.; Choi, W.; Lowe, M.J.S.; Skelton, E.A.; Craster, R.V. The Validity of Kirchhoff Theory for Scattering of Elastic Waves from Rough Surfaces. Proc. R. Soc. Math. Phys. Eng. Sci. 2015, 471, 20140977. [Google Scholar] [CrossRef]
  12. Choi, W.; Shi, F.; Lowe, M.J.; Skelton, E.A.; Craster, R.V.; Daniels, W.L. Rough Surface Reconstruction of Real Surfaces for Numerical Simulations of Ultrasonic Wave Scattering. NDT E Int. 2018, 98, 27–36. [Google Scholar] [CrossRef] [Green Version]
  13. Wang, S.; Shen, J.; Kin Chan, W. Determination of the Fractal Scaling Parameter from Simulated Fractal-Regular Surface Profiles Based on the Weierstrass-Mandelbrot Function. J. Tribol. 2007, 129, 952–956. [Google Scholar] [CrossRef]
  14. Weng, J.; Cohen, P.; Herniou, M. Camera Calibration with Distortion Models and Accuracy Evaluation. IEEE Trans. Pattern Anal. Mach. Intell. 1992, 14, 965–980. [Google Scholar] [CrossRef] [Green Version]
  15. Wang, J.; Shi, F.; Zhang, J.; Liu, Y. A New Calibration Model and Method of Camera Lens Distortion. In Proceedings of the IEEE/RSJ International Conference on Intelligent Robots and Systems, Beijing, China, 9 October 2006. [Google Scholar]
  16. Akca, E.; Gürsel, A. A Review on Superalloys and IN718 Nickel-Based INCONEL Superalloy. Period. Eng. Nat. Sci. PEN 2015, 3, 15–27. [Google Scholar] [CrossRef]
  17. Sawafuji, Y. Automatic Ultrasonic Testing of Non-Metallic Inclusions Detectable with Size of Several Tens of Micrometers Using a Double Probe Technique along the Longitudinal Axis of a Small-Diameter Bar. ISIJ Int. 2021, 61, 248–257. [Google Scholar] [CrossRef]
  18. Deutsch, W.A.K.; Joswig, M.; Kattwinkel, R.; Roye, W.; Maxam, K.; Ranzeng, M. Automated Ultrasonic Testing Systems for Bars and Tubes, Examples with Mono-Element and Phased Array Probes. In Proceedings of the ECNDT 2014, Prague, Czech Republic, 6 October 2014. [Google Scholar]
  19. Fan, Y.; Sinclair, A.N.; Honarvar, F. Ultrasonic Characterization of Continuously Cast Rod by Resonance Acoustic Spectroscopy. Nondestruct. Test. Eval. 2003, 19, 15–28. [Google Scholar] [CrossRef]
  20. Yi, L. The Development of the Automatic Ultrasonic Testing System for Bar Materials. Master’s Thesis, Chongqing University, Chongqing, China, 2006. [Google Scholar]
  21. Haldipur, P.; Margetan, F.J.; Thompson, R.B. Estimation of Single-Crystal Elastic Constants from Ultrasonic Measurements on Polycrystalline Specimens. AIP Conf. Proc. 2004, 700, 1061–1068. [Google Scholar] [CrossRef]
  22. Yang, S.; Yang, S.; Qu, J.; Du, J.; Gu, Y.; Zhao, P.; Wang, N. Inclusions in Wrought Superalloys: A Review. J. Iron Steel Res. Int. 2021, 28, 921–937. [Google Scholar] [CrossRef]
  23. Haldipur, P.; Margetan, F.J.; Thompson, R.B. Estimation of Single-Crystal Elastic Constants of Polycrystalline Materials from Back-Scattered Grain Noise. AIP Conf. Proc. 2006, 820, 1133–1140. [Google Scholar] [CrossRef]
Figure 1. Schematic diagram of the distortions of a 3D cube and a 2D checkerboard: (a1) original image; (b1) pincushion distortion; (c1) barrel distortion; (d1) ripple distortion; (a2) original image; (b2) pincushion distortion; (c2) barrel distortion; (d2) ripple distortion.
Figure 1. Schematic diagram of the distortions of a 3D cube and a 2D checkerboard: (a1) original image; (b1) pincushion distortion; (c1) barrel distortion; (d1) ripple distortion; (a2) original image; (b2) pincushion distortion; (c2) barrel distortion; (d2) ripple distortion.
Materials 15 08633 g001
Figure 2. Flowchart of multi-scale distortion.
Figure 2. Flowchart of multi-scale distortion.
Materials 15 08633 g002
Figure 3. Distorted image: (a) a crack; (b) a pore; (c) an inclusion. The similarities and differences between three typical repetitions are highlighted using pseudo-color imaging, where the shared areas are in white.
Figure 3. Distorted image: (a) a crack; (b) a pore; (c) an inclusion. The similarities and differences between three typical repetitions are highlighted using pseudo-color imaging, where the shared areas are in white.
Materials 15 08633 g003
Figure 4. Setup for ultrasonic inspection of nickel-based bars with a mode-converted transverse wave.
Figure 4. Setup for ultrasonic inspection of nickel-based bars with a mode-converted transverse wave.
Materials 15 08633 g004
Figure 5. Attenuation coefficients in the frequency domain.
Figure 5. Attenuation coefficients in the frequency domain.
Materials 15 08633 g005
Figure 6. Snapshots of wave field for a smooth-surface notch.
Figure 6. Snapshots of wave field for a smooth-surface notch.
Materials 15 08633 g006
Figure 7. B scan image: (a) simulation; (b) experiment. The width of the notch was 0.2 mm.
Figure 7. B scan image: (a) simulation; (b) experiment. The width of the notch was 0.2 mm.
Materials 15 08633 g007
Figure 8. Pulse-echo signal when the notch is at 45° (normalized with the amplitude of each flaw echo).
Figure 8. Pulse-echo signal when the notch is at 45° (normalized with the amplitude of each flaw echo).
Materials 15 08633 g008
Figure 9. The flaw echoes from the rough crack in time and frequency domain. (a) Multiple repetitions in time domain; (b) multiple repetitions in frequency domain; (c) statistical results in time domain; (d) statistical results in frequency domain. The width of the notch was 0.02 mm.
Figure 9. The flaw echoes from the rough crack in time and frequency domain. (a) Multiple repetitions in time domain; (b) multiple repetitions in frequency domain; (c) statistical results in time domain; (d) statistical results in frequency domain. The width of the notch was 0.02 mm.
Materials 15 08633 g009
Figure 10. Scatter plot of the amplitude and TOF.
Figure 10. Scatter plot of the amplitude and TOF.
Materials 15 08633 g010
Figure 11. Statistical results for the rough pore and rough inclusion. (a) Rough pore, time domain; (b) rough pore, frequency domain; (c) rough inclusion, time domain; (d) rough inclusion, frequency domain.
Figure 11. Statistical results for the rough pore and rough inclusion. (a) Rough pore, time domain; (b) rough pore, frequency domain; (c) rough inclusion, time domain; (d) rough inclusion, frequency domain.
Materials 15 08633 g011
Table 1. Properties of reference nickel-based alloy and non-metallic inclusion.
Table 1. Properties of reference nickel-based alloy and non-metallic inclusion.
MaterialsElastic Modulus
(GPa)
Density
(g/cm3)
Poisson Ratio
Nickel-based alloy [21]2008.240.29
Al2O3-MgO-TiN [22]2556.580.25
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wang, Z.; Zeng, Z.; Song, Y.; Li, X. Finite Element Simulation of Ultrasonic Scattering by Rough Flaws with Multi-Scale Distortions. Materials 2022, 15, 8633. https://doi.org/10.3390/ma15238633

AMA Style

Wang Z, Zeng Z, Song Y, Li X. Finite Element Simulation of Ultrasonic Scattering by Rough Flaws with Multi-Scale Distortions. Materials. 2022; 15(23):8633. https://doi.org/10.3390/ma15238633

Chicago/Turabian Style

Wang, Zheng, Zhanhong Zeng, Yongfeng Song, and Xiongbing Li. 2022. "Finite Element Simulation of Ultrasonic Scattering by Rough Flaws with Multi-Scale Distortions" Materials 15, no. 23: 8633. https://doi.org/10.3390/ma15238633

APA Style

Wang, Z., Zeng, Z., Song, Y., & Li, X. (2022). Finite Element Simulation of Ultrasonic Scattering by Rough Flaws with Multi-Scale Distortions. Materials, 15(23), 8633. https://doi.org/10.3390/ma15238633

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop