Microstructural Characterization of Additively Manufactured Metal Components Using Linear and Nonlinear Ultrasonic Techniques

Metal additive manufacturing (AM) is an innovative manufacturing technology that uses a high-power laser for the layer-by-layer production of metal components. Despite many achievements in the field of AM, few studies have focused on the nondestructive characterization of microstructures, such as grain size and porosity. In this study, various microstructures of additively manufactured metal components were characterized non-destructively using linear/nonlinear ultrasonic techniques. The contributions of this study are as follows: (1) presenting correlation analyses of various microstructures (grain size and texture, lack of fusion, and porosity) and ultrasonic properties (ultrasonic velocity, attenuation, and nonlinearity parameters), (2) development of nondestructive microstructural characterization techniques for additively manufactured components; and (3) exploring the potential for the online monitoring of AM processes owing to the nondestructive nature of the proposed technique. The performance of the proposed technique was validated using additively manufactured samples under varying laser beam speed conditions. The characteristics of the target microstructures characterized using the proposed technique were consistent with the results obtained using destructive optical microscopy and electron back-scattered diffraction methods.


Introduction
Metal additive manufacturing (AM), also referred to as metal 3D printing, is an emerging technology for producing metal components [1]. AM is used to produce 3D metal components via layer-by-layer deposition of metal wires, sheets, or powders, wherein a target region is melted using a heat source, such as a laser, arc, or electron beam [2]. AM facilitates innovative, flexible, and lightweight designs for 3D structures with minimal tooling, reduced material waste and lead times, and fast prototyping compared with conventional subtractive or casting manufacturing [3].
Despite the many advantages of AM, the quality assurance of AM components is a significant barrier to its adoption [4]. Typically, the presence of inhomogeneous microstructures in the interior of AM components owing to insufficient energy irradiation, material overheating, or prolonged exposure to varying temperatures in the deposited layer, results in uneven grain structures, gas pores, and a lack of fusion (LOF) [5] in the final products. To address these problems, several nondestructive testing (NDT) methods, such as vision-based testing [6], pyrometry [7], thermography [8], and radiography [9] have been proposed. Although the aforementioned NDT methods are effective for evaluating are affected by the porosity content and distribution [35], and the grain size and texture [30] and can be obtained as [36]: where x is the ultrasonic propagation distance and t is the arrival time of the ultrasonic wave. The parameter a, which is another linear ultrasonic property, quantifies the ultrasonic amplitude loss due to scattering, absorption, or friction when ultrasonic waves propagate through a material [16]. a is often dependent on microstructural features such as grain size, boundary [37], textures [38], pores [12], and inclusions [39]. For conventional metals, nondestructive evaluation of various microstructural features using ultrasonic forward or backscattering has been extensively reported [40]. This ultrasonic attenuation has been actively applied by several researchers in AM for melt-pool monitoring during the AM process [41] and mechanical property estimation [30,31]. Generally, the mechanisms of varying a differ according to the ratio of the ultrasonic wavelength and the existing scatterers in the material [10]. Because the range of ultrasonic frequency commonly used for NDT is 2-10 MHz, the ultrasonic wavelengths for common AM metals (stainless steel and titanium alloys) range from 0.6-3 mm, which is larger than the material microstructural features. Then, the Rayleigh scattering theorem becomes a dominant attenuation mechanism and can be obtained as follows [10]: where S 0 and S x are the initial amplitude of the ultrasonic wave and the amplitude after the propagation distance, respectively. Parameter β is an index related to the material second-order nonlinear elastic constants [17]. This parameter depends on the monochromatic fundamental ultrasonic wave and the higher harmonic components generated by the distortion of the fundamental wave. Because β is known to be sensitive to the material microstructural characteristics, such as dislocations [19], precipitates [42], microscale cracks [14], and grain boundaries [37], its application has recently been actively investigated. β can be derived as [17]: where k is the ultrasonic wavenumber and A 1 and A 2 are the displacement amplitudes of the fundamental and second-order harmonic frequency components of the propagating ultrasonic wave, respectively. Owing to the difficulty in precisely measuring ultrasonic displacements, several researchers have used a relative parameter β instead of β, which is obtained as follows [43]: where A 1 and A 2 indicate the electrical signal amplitudes for the respective frequency component measured by an ultrasonic detector and D is the ultrasonic attenuation compensation factor. Although β is the relative value of ultrasonic nonlinearity, it can be used to compare the before/after microstructural variations under constant material conditions [17].
Notably, ultrasonic attenuation also results in a decrease in A 1 and then results in β being overrated during the ultrasonic measurement [44]. To prevent this problem, D is included to distinguish the effects of nonlinear harmonic generation from those of attenuation, which is derived as follows [44]: where a 1 and a 2 are the attenuation parameters at the fundamental and second-order harmonic frequencies, respectively.

Samples
Eight test samples with different microstructural characteristics were fabricated using a direct-metal-printing (DMP) 3D-printer (ProX DMP 320; 3D Systems, South Caroline, SC 45, USA) with LaserForm TM 316 L stainless steel powder. The microstructures considered were the grain size, texture, LOF porosity, and microporosity, which are dominantly attributable to the mechanical properties of the AM components [45]. To induce these microstructures in the test samples, we used eight LBS conditions of the 3D-printer: 0.6 m/s, 0.68 m/s, 0.92 m/s, 1 m/s, 1.08 m/s, 1.16 m/s, 1.4 m/s, and 1.48 m/s. They were divided into three groups according to the LBS conditions: low (0.6 m/s and 0.68 m/s), medium (0.92 m/s, 1 m/s, 1.08 m/s, and 1.16 m/s), and high (1.4 m/s and 1.48 m/s) LBS conditions. Each sample was 20 mm long, 20 mm wide, and 8 mm high and labeled as #1 through #8, as shown in Figure 1, according to the LBS conditions used. A previous study reported that variations in the LBS induce various microstructures in the AM components targeted in our study [45]. The other AM process parameters, including the laser power, hatch space, and layer thickness, were fixed at 200 W, 0.1 mm, and 0.03 mm, respectively. Notably, ultrasonic measurement is inconsistent under the rough surface conditions of the test samples because the signal-to-noise ratio is lower in the case of a rough surface [46]. To avoid this problem, the tested samples were treated via wire electrical discharge machining, wherein the arithmetic mean roughness (Ra) of the surfaces was smoothed to approximately 0.33 µm (initial Ra: approximately 5.0 µm). where 1 and 2 are the attenuation parameters at the fundamental and second-order harmonic frequencies, respectively.

Samples
Eight test samples with different microstructural characteristics were fabricated using a direct-metal-printing (DMP) 3D-printer (ProX DMP 320; 3D Systems, South Caroline, SC 45, USA) with LaserForm TM 316 L stainless steel powder. The microstructures considered were the grain size, texture, LOF porosity, and microporosity, which are dominantly attributable to the mechanical properties of the AM components [45]. To induce these microstructures in the test samples, we used eight LBS conditions of the 3D-printer: 0.6 m/s, 0.68 m/s, 0.92 m/s, 1 m/s, 1.08 m/s, 1.16 m/s, 1.4 m/s, and 1.48 m/s. They were divided into three groups according to the LBS conditions: low (0.6 m/s and 0.68 m/s), medium (0.92 m/s, 1 m/s, 1.08 m/s, and 1.16 m/s), and high (1.4 m/s and 1.48 m/s) LBS conditions. Each sample was 20 mm long, 20 mm wide, and 8 mm high and labeled as #1 through #8, as shown in Figure 1, according to the LBS conditions used. A previous study reported that variations in the LBS induce various microstructures in the AM components targeted in our study [45]. The other AM process parameters, including the laser power, hatch space, and layer thickness, were fixed at 200 W, 0.1 mm, and 0.03 mm, respectively. Notably, ultrasonic measurement is inconsistent under the rough surface conditions of the test samples because the signal-to-noise ratio is lower in the case of a rough surface [46]. To avoid this problem, the tested samples were treated via wire electrical discharge machining, wherein the arithmetic mean roughness (Ra) of the surfaces was smoothed to approximately 0.33 μm (initial Ra: approximately 5.0 μm).

Ultrasonic Measurements
In this study, pulse-echo (PE) [47] and through-transmission (TT) [48] modes were used, where the former measures parameters and and the latter measures parameter ′. For the PE mode, a single ultrasonic transducer was used, and two consecutive backwall echoes were measured. Here, the respective (ultrasonic propagation distance) and (arrival-time) correspond to twice the sample thickness and time-of-flight (TOF) between two consecutive echoes. Figure 2 shows the experimental setup for the PE and TT modes. A commercial pulser/receiver generated an electrical pulsed signal and excited a commercial piezoelectric contact transducer with a frequency of 5 MHz. A longitudinal ultrasonic wave excited by the transducer was transmitted to the test sample, propagated through it, and reflected on the other side. The ultrasonic signals of the first and second back-wall echoes reflected were received by the same transducer and pulser/receiver and then digitalized on a commercial oscilloscope with a 0.1 ns time resolution and 300 times averaging. The measurements were repeated 10 times for each sample. Because the imperfect contact conditions between the test sample and transducer cause a measurement error owing to the multireflection and attenuation in the contact gap [49], a specially designed pneumatic device was used, in which the contact condition in each measurement was constant at 0.4 MPa, as shown in Figure 2b. Notably, when the ultrasonic wave propagates over a certain distance (the so-called Fresnel near-field zone [43]), the ultrasonic wave in a material may diffract and result in diffraction attenuation, which is independent of the microstructural characteristics. Here, the wave propagation distance is 16 mm

Ultrasonic Measurements
In this study, pulse-echo (PE) [47] and through-transmission (TT) [48] modes were used, where the former measures parameters v and a and the latter measures parameter β . For the PE mode, a single ultrasonic transducer was used, and two consecutive backwall echoes were measured. Here, the respective x (ultrasonic propagation distance) and t (arrival-time) correspond to twice the sample thickness and time-of-flight (TOF) between two consecutive echoes. Figure 2 shows the experimental setup for the PE and TT modes. A commercial pulser/receiver generated an electrical pulsed signal and excited a commercial piezoelectric contact transducer with a frequency of 5 MHz. A longitudinal ultrasonic wave excited by the transducer was transmitted to the test sample, propagated through it, and reflected on the other side. The ultrasonic signals of the first and second back-wall echoes reflected were received by the same transducer and pulser/receiver and then digitalized on a commercial oscilloscope with a 0.1 ns time resolution and 300 times averaging. The measurements were repeated 10 times for each sample. Because the imperfect contact conditions between the test sample and transducer cause a measurement error owing to the multireflection and attenuation in the contact gap [49], a specially designed pneumatic device was used, in which the contact condition in each measurement was constant at 0.4 MPa, as shown in Figure 2b. Notably, when the ultrasonic wave propagates over a certain distance (the so-called Fresnel near-field zone [43]), the ultrasonic wave in a material may diffract and result in diffraction attenuation, which is independent of the microstructural characteristics. Here, the wave propagation distance is 16 mm (twice the sample thickness) within the Fresnel zone (approximately 19 mm [43]); therefore, this effect was neglected.
The β parameter was measured using the TT mode, wherein a 5 MHz transducer was used as the transmitter and a 10 MHz transducer was used as the detector to measure the second harmonic component more suitably. In the TT mode, an ultrasonic wave with a long pulse length is used to obtain a narrower frequency spectrum than that of the PE mode because the narrow-band spectrum is beneficial for dividing the first and second harmonic components in the frequency domain [43]. Therefore, a voltage sinusoidal signal comprising nine cycles was used in the TT mode, generated by a pulser (RAM-5000; RITEC, Warwick, 44-74300, UK) comprising a synthesizer and high-power gated amplifier. This signal passed through a 12 dB attenuator, a high-power 50 Ω termination, and a 7 MHz low-pass filter to suppress the initial harmonic frequency components generated in the pulser and drive the 5 MHz transmitting transducer. Suppressing the initial harmonic components is important for measuring just the harmonic components that correspond to the microstructural characteristics of the test samples [17]. The wave propagated through the sample was received on the other side of the sample by the detector (10 MHz transducer). The measurement was repeated by increasing the input power to obtain a linear fitting plot of A 1 2 and A 2 (proportional to β ). A specially designed pneumatic control system was also used to obtain accurate measurements [50]. Here, the D value (attenuation compensation factor (D)) was obtained from the aforementioned PE mode. A schematic and a photograph of the TT mode are shown in Figure 2c,d.
(twice the sample thickness) within the Fresnel zone (approximately 19 mm [43]); therefore, this effect was neglected. The ′ parameter was measured using the TT mode, wherein a 5 MHz transducer was used as the transmitter and a 10 MHz transducer was used as the detector to measure the second harmonic component more suitably. In the TT mode, an ultrasonic wave with a long pulse length is used to obtain a narrower frequency spectrum than that of the PE mode because the narrow-band spectrum is beneficial for dividing the first and second harmonic components in the frequency domain [43]. Therefore, a voltage sinusoidal signal comprising nine cycles was used in the TT mode, generated by a pulser (RAM-5000; RITEC, Warwick, 44-74300, UK) comprising a synthesizer and high-power gated amplifier. This signal passed through a 12 dB attenuator, a high-power 50 Ω termination, and a 7 MHz low-pass filter to suppress the initial harmonic frequency components generated in the pulser and drive the 5 MHz transmitting transducer. Suppressing the initial harmonic components is important for measuring just the harmonic components that correspond to the microstructural characteristics of the test samples [17]. The wave propagated through the sample was received on the other side of the sample by the detector (10 MHz transducer). The measurement was repeated by increasing the input power to obtain a linear fitting plot of ( 1 ′ ) 2 and A 2 ′ (proportional to ′). A specially designed pneumatic control system was also used to obtain accurate measurements [50]. Here, the value (attenuation compensation factor (D)) was obtained from the aforementioned PE mode. A schematic and a photograph of the TT mode are shown in Figure 2c,d.

Metallography
After the ultrasonic measurements, X-ray diffraction (XRD), OM, and EBSD were performed on small samples, which were detached from the test samples for a comparison with the ultrasonic results. Phase identification was performed using a commercial X-ray diffractometer with a step size of 0.01° and 2θ ranging from 20 to 80°. Porosity, grain size,

Metallography
After the ultrasonic measurements, X-ray diffraction (XRD), OM, and EBSD were performed on small samples, which were detached from the test samples for a comparison with the ultrasonic results. Phase identification was performed using a commercial X-ray diffractometer with a step size of 0.01 • and 2θ ranging from 20 to 80 • . Porosity, grain size, and texture measurements were performed using a commercial optical microscope and EBSD instrument, respectively, at the cross-section parallel to the building direction of the additional samples. The grain sizes and porosity contents were measured and averaged using "Quatro S" and "ImageJ", respectively, which are commercial software programs for the size measurement. In these procedures, additional mechanical polishing with diamond suspensions was performed on surfaces, down to 0.04 µm.

Metallography Results
An XRD result for the 1.16 m/s LBS condition is presented in Figure 3, which indicates that the major phase with the highest peak intensity of the sample was cubic FCC-austenite γ, in accordance with previous studies [51]. For the other specimens, similar XRD results were observed. Based on the XRD results, the grain size and texture were measured for all the test samples via EBSD.  Figure 4b,d,f show the crystallographic textures for the representative 0.6 m/s, 1.0, and 1.4 m/s LBS conditions. A high texture intensity of up to 4.9, along the 101 direction, was observed at a 0.6 m/s LBS, and the textures were more random with the increased LBS conditions. An increase in LBS has been reported to reduce grain size and texture intensity in a specific direction due to the decreased laser energy input in the melt pool and interaction time between the powders and laser, consequently restricting grain growth [52]. and texture measurements were performed using a commercial optical microscope and EBSD instrument, respectively, at the cross-section parallel to the building direction of the additional samples. The grain sizes and porosity contents were measured and averaged using "Quatro S" and "ImageJ", respectively, which are commercial software programs for the size measurement. In these procedures, additional mechanical polishing with diamond suspensions was performed on surfaces, down to 0.04 μm.

Metallography Results
An XRD result for the 1.16 m/s LBS condition is presented in Figure 3, which indicates that the major phase with the highest peak intensity of the sample was cubic FCC-austenite γ, in accordance with previous studies [51]. For the other specimens, similar XRD results were observed. Based on the XRD results, the grain size and texture were measured for all the test samples via EBSD.  Figure 4b,d,f show the crystallographic textures for the representative 0.6 m/s, 1.0, and 1.4 m/s LBS conditions. A high texture intensity of up to 4.9, along the 101 direction, was observed at a 0.6 m/s LBS, and the textures were more random with the increased LBS conditions. An increase in LBS has been reported to reduce grain size and texture intensity in a specific direction due to the decreased laser energy input in the melt pool and interaction time between the powders and laser, consequently restricting grain growth [52].   EBSD instrument, respectively, at the cross-section parallel to the building direction of additional samples. The grain sizes and porosity contents were measured and avera using "Quatro S" and "ImageJ", respectively, which are commercial software progra for the size measurement. In these procedures, additional mechanical polishing with d mond suspensions was performed on surfaces, down to 0.04 μm.

Metallography Results
An XRD result for the 1.16 m/s LBS condition is presented in Figure 3, which indica that the major phase with the highest peak intensity of the sample was cubic FCC-aus ite γ, in accordance with previous studies [51]. For the other specimens, similar XRD sults were observed. Based on the XRD results, the grain size and texture were measu for all the test samples via EBSD.  Figure 4b,d,f show the crystallographic textures for the representa 0.6 m/s, 1.0, and 1.4 m/s LBS conditions. A high texture intensity of up to 4.9, along 101 direction, was observed at a 0.6 m/s LBS, and the textures were more random with increased LBS conditions. An increase in LBS has been reported to reduce grain size texture intensity in a specific direction due to the decreased laser energy input in the m pool and interaction time between the powders and laser, consequently restricting gr growth [52].   The porosity was evaluated using OM and its content was calculated from 2D images as the ratio of the pore area to the total area [53]. The typical OM images for 0.6 m/s, 1.0 m/s, and 1.4 m/s conditions are presented in Figure 5a,b and Figure 5c, resp tively. The occurrence of various types of pores, including gas pores and small pits m/s LBS), and LOF with unmelted powders (1.4 m/s LBS), indicated by black spots, w observed, whereas almost no pores were observed in the case of the 1.0 m/s LBS. The creased LBS at 0.6 m/s results in excessive laser energy input in the melt pool, genera various discontinuities, such as wavy surfaces and welded particles on the deposited la [54]. In general, these discontinuities can be eliminated when the next AM layer is dep ited in the following layer (the so-called re-melting process). However, the lower la energy input according to the decreased LBS often results in a melt pool with a shal shape, thus resulting in insufficient penetration of the previously deposited layer [ Consequently, pores with small pits were observed between the interlayers. Furtherm excessive energy input increases gas solubility, resulting in increased gas porosity [ The observation of higher porosity contents for the 1.4 m/s LBS condition is attribute the lack of laser energy input in the melt pool during the AM [54], wherein the insuffic energy causes LOF with the unfused powders between layers or within layers. The porosity was evaluated using OM and its content was calculated from 2D OM images as the ratio of the pore area to the total area [53]. The typical OM images for the 0.  [54]. In general, these discontinuities can be eliminated when the next AM layer is deposited in the following layer (the so-called re-melting process). However, the lower laser-energy input according to the decreased LBS often results in a melt pool with a shallow shape, thus resulting in insufficient penetration of the previously deposited layer [55]. Consequently, pores with small pits were observed between the interlayers. Furthermore, excessive energy input increases gas solubility, resulting in increased gas porosity [56]. The observation of higher porosity contents for the 1.4 m/s LBS condition is attributed to the lack of laser energy input in the melt pool during the AM [54], wherein the insufficient energy causes LOF with the unfused powders between layers or within layers. The microstructural characteristics of each LBS condition, based on the metallographic results, are summarized in Table 1. Under LBS conditions lower than 0.68 m/s (low LBS condition), approximately 100 μm grain size, approximately 5.2 texture intensity along the 101 direction, and submicrometer microporosity (porosity content: 1.8%) were observed. Under LBS conditions between 0.92 m/s and 1.16 m/s (medium LBS condition), the average grain size and texture intensity along the 101 direction were 68 μm and 2.6 μm, respectively, and almost no pores were observed (porosity content: 1.1%). As the LBS further increased (high LBS condition), the average grain size (53 μm) and texture intensity (1.9 along the 101 direction) decreased, i.e., textures were more random, and submillimeter-scale LOF porosity (porosity content: 7.0%) occurred. Accordingly, each microstructural characteristic was labeled as "State I," "State II," and "State III," respectively, as denoted in Table 1.  The microstructural characteristics of each LBS condition, based on the metallographic results, are summarized in Table 1. Under LBS conditions lower than 0.68 m/s (low LBS condition), approximately 100 µm grain size, approximately 5.2 texture intensity along the 101 direction, and submicrometer microporosity (porosity content: 1.8%) were observed. Under LBS conditions between 0.92 m/s and 1.16 m/s (medium LBS condition), the average grain size and texture intensity along the 101 direction were 68 µm and 2.6 µm, respectively, and almost no pores were observed (porosity content: 1.1%). As the LBS further increased (high LBS condition), the average grain size (53 µm) and texture intensity (1.9 along the 101 direction) decreased, i.e., textures were more random, and submillimeter-scale LOF porosity (porosity content: 7.0%) occurred. Accordingly, each microstructural characteristic was labeled as "State I," "State II," and "State III," respectively, as denoted in Table 1.  Figure 6 presents an ultrasonic signal acquired from the PE mode, where the parameters v and a were obtained using Equations (1) and (2), respectively. Here, v was calculated from x (ultrasonic propagation distance) and TOF between the first and second peaks of the back-wall echoes, as shown in Figure 6, wherein TOF was measured using a cross-correlation function [24]. When calculating a, the peak amplitudes of the first and second back-wall echoes were used as the S 0 and S x values, as shown in Figure 6. Figure 7a presents the ultrasonic signal comprising a sinusoidal wave with nine cycles in the time domain measured from the TT mode. Figure 7b presents the signal in the frequency domain that was Fourier transformed from the time-domain signal using the fast Fourier transform (FFT). Here, a Hanning window was used to obtain a good frequency resolution, thus resulting in high amplitude accuracy in the frequency domain [49]. In the frequency domain, the voltage signal amplitudes for the 5 MHz fundamental and 10 MHz second-harmonic frequency components were extracted and presented as A 1 and A 2 , respectively, as shown in Figure 7b. Figure 7c presents the relationship between A 1 2 and A 2 with increasing voltage input to extract the β parameter. A linear relationship with a high correlation coefficient (0.999) was observed. Finally, β was calculated using Equation (4). As mentioned above, the value of D was obtained using the PE mode. Here, a 1 and a are identical, and a 2 is obtained using an additional 10 MHz transducer [44].

Ultrasonic Signal Acquisition Results
culated from (ultrasonic propagation distance) and TOF b peaks of the back-wall echoes, as shown in Figure 6, wherein cross-correlation function [24]. When calculating , the peak second back-wall echoes were used as the 0 ′ and ′ value ure 7a presents the ultrasonic signal comprising a sinusoidal time domain measured from the TT mode. Figure 7b present domain that was Fourier transformed from the time-domain transform (FFT). Here, a Hanning window was used to obta tion, thus resulting in high amplitude accuracy in the freque quency domain, the voltage signal amplitudes for the 5 MH second-harmonic frequency components were extracted an respectively, as shown in Figure 7b. Figure 7c presents the and 2 ′ with increasing voltage input to extract the ′ para with a high correlation coefficient (0.999) was observed. Fina Equation (4). As mentioned above, the value of was obtain 1 and are identical, and 2 is obtained using an addition    Figure 6 presents an ultrasonic signal acquired from the PE mode, where the param eters and were obtained using Equations (1) and (2), respectively. Here, was cal culated from (ultrasonic propagation distance) and TOF between the first and second peaks of the back-wall echoes, as shown in Figure 6, wherein TOF was measured using cross-correlation function [24]. When calculating , the peak amplitudes of the first and second back-wall echoes were used as the 0 ′ and ′ values, as shown in Figure 6. Fig  ure 7a presents the ultrasonic signal comprising a sinusoidal wave with nine cycles in th time domain measured from the TT mode. Figure 7b presents the signal in the frequency domain that was Fourier transformed from the time-domain signal using the fast Fourie transform (FFT). Here, a Hanning window was used to obtain a good frequency resolu tion, thus resulting in high amplitude accuracy in the frequency domain [49]. In the fre quency domain, the voltage signal amplitudes for the 5 MHz fundamental and 10 MHz second-harmonic frequency components were extracted and presented as 1 ′ and 2 ′ respectively, as shown in Figure 7b. Figure 7c presents the relationship between ( 1 ′ ) and 2 ′ with increasing voltage input to extract the ′ parameter. A linear relationship with a high correlation coefficient (0.999) was observed. Finally, ′ was calculated using Equation (4). As mentioned above, the value of was obtained using the PE mode. Here 1 and are identical, and 2 is obtained using an additional 10 MHz transducer [44].       Figure 8a, the parameter increases as the LBS increases from 0.6 to 1.1 in State I and II), before sharply decreasing as LBS increases further (i.e., in Stat hough increased by 0.42% in State II, it sharply decreased by 0.65% in the zone of State II and III. At the highest LBS condition of 1.48 m/s in State III, the n parameter was the lowest, i.e., 0.74% less than that at the 1.16 m/s LBS (= 0 c When State I transitions to State II, almost no slope variation is observed. In Fig  decreases from State I to the first LBS condition of State III, corresponding to 0.6-1.4 m/s; however, a reversed trend is observed as the LBS further increases to In Figure 8c, the ′ value increases as the LBS value increases from 0.92 m/s t (in State II and III); however, a notable slope variation is observed at the LBS c In Figure 8a, the parameter v increases as the LBS increases from 0.6 to 1.16 m/s (i.e., in State I and II), before sharply decreasing as LBS increases further (i.e., in State III). Although v increased by 0.42% in State II, it sharply decreased by 0.65% in the transition zone of State II and III. At the highest LBS condition of 1.48 m/s in State III, the normalized parameter v was the lowest, i.e., 0.74% less than that at the 1.16 m/s LBS (=v 0 condition). When State I transitions to State II, almost no slope variation is observed. In Figure 8b, a decreases from State I to the first LBS condition of State III, corresponding to an LBS of 0.6-1.4 m/s; however, a reversed trend is observed as the LBS further increases to 1.48 m/s. In Figure 8c, the β value increases as the LBS value increases from 0.92 m/s to 1.48 m/s (in State II and III); however, a notable slope variation is observed at the LBS conditions of 0.6 m/s and 0.68 m/s (in State I) as compared to those of State II and III. Although a 3.8% change in the normalized β was observed in State II, only a 0.16% change occurred in the transition zone of State I and II. At the minimum LBS of 0.6 m/s, the slope was reversed, and the normalized β increased by 0.23% compared to that for 0.68 m/s LBS. Almost no slope variation was observed between states II and III.

Microstructural Characterization Using Ultrasonic NDT
From the ultrasonic measurement results, ultrasonic parameters that are effective in characterizing the microstructural evolution at the respective states I, II, and III were selected. State II (0.92-1.16 m/s LBSs), where only the grain size and crystallographic texture vary with LBS, was analyzed first. Figure 9 presents a comparison of the normalized v, a, and β parameters correlated with varying LBS in State II. The three parameters show a consistent variation, wherein each normalized v, a, and β parameter decreased, slightly increased, and increased, respectively. Comparing the average variations for each measured parameter, the most sensitive parameter depending on LBS (i.e., grain size and crystallographic texture) in State II was β . The absolute variation in β with the 0.24 m/s LBS increment was 3.8%, which is approximately nine and two times higher than those of v (0.42%) and a (2.0%), respectively. Because only the grain size and texture varied in State II, it can be concluded that β is the most effective parameter for evaluating grain characteristics. Concurrently, despite being less than β , a also exhibited a good correlation. aterials 2022, 15, x FOR PEER REVIEW , and ′ parameters correlated with varying LBS in State II. The three par a consistent variation, wherein each normalized , , and ′ parame slightly increased, and increased, respectively. Comparing the average vari measured parameter, the most sensitive parameter depending on LBS (i.e., crystallographic texture) in State II was ′. The absolute variation in ′ wi LBS increment was 3.8%, which is approximately nine and two times highe (0.42%) and (2.0%), respectively. Because only the grain size and tex State II, it can be concluded that ′ is the most effective parameter for ev characteristics. Concurrently, despite being less than ′, also exhibited a tion. The high sensitivity of ′ in State II (0.92 m/s-1.16 m/s LBSs) is attrib lationship between the AM microstructural features, including grain size an the produced ultrasonic nonlinearity. An increase in LBS leads to a fast cooli the AM process, which results in randomly distributed textures, as show These randomly distributed textures have been reported to produce high u linearity [33]. Furthermore, these random textures owing to the fast cooling reported to contain a high dislocation density, which can also contribute to ′ [33]. In addition, a previous study [37] reported that a decrease in the creases the number of grain boundaries in a unit volume, in which the gra often act as interfacial discontinuities and generate higher harmonics for th ultrasonic wave. To further substantiate this observation, the high-angle g maps for the 0.92 m/s (grain size: 75 μm) and 1.16 m/s (grain size: 59 μm) L were obtained using EBSD and compared, as shown in Figure 10, where gra The high sensitivity of β in State II (0.92 m/s-1.16 m/s LBSs) is attributed to the relationship between the AM microstructural features, including grain size and texture, and the produced ultrasonic nonlinearity. An increase in LBS leads to a fast cooling rate during the AM process, which results in randomly distributed textures, as shown in Figure 4. These randomly distributed textures have been reported to produce high ultrasonic nonlinearity [33]. Furthermore, these random textures owing to the fast cooling rate have been reported to contain a high dislocation density, which can also contribute to the increased β [33]. In addition, a previous study [37] reported that a decrease in the grain size increases the number of grain boundaries in a unit volume, in which the grain boundaries often act as interfacial discontinuities and generate higher harmonics for the propagating ultrasonic wave. To further substantiate this observation, the highangle grain-boundary maps for the 0.92 m/s (grain size: 75 µm) and 1.16 m/s (grain size: 59 µm) LBS conditions were obtained using EBSD and compared, as shown in Figure 10, where grain boundaries with a misorientation angle exceeding 15 • are illustrated by gray lines [57]. More grain boundaries were observed at 1.16 m/s LBS than 0.92 m/s LBS. Although linear ultrasonic properties based on the scattering attenuation mechanism are reportedly sensitive to grain effects in conventional metals [40], β was shown to be the most sensitive factor in State II in our experimental results. This may be attributed to the unique thermomechanical manufacturing process of AM, which causes multiple microstructural variations, such as varying grain size and randomly distributed textures; however, further studies are needed to elucidate this observation in detail. A similar observation in AM samples has been reported in the literature [33], which also includes the observation of varying β in the AM samples according to the different heat-treatment times. However, the high correlation between the grain characteristics and β in the AM samples fabricated under various processing conditions was uniquely observed and investigated in this study. The second-highest correlation of ′ in State II is shown in Figure 9. This observati is attributed to the ultrasonic attenuation at the grain boundaries, which has been exte sively investigated [10,40]. An ultrasonic wave propagating through the grain boundar scatters and its amplitude decreases according to the Rayleigh scattering theory, as f lows [10]: is the grain size, is the ultrasonic frequency, and superscript b is a consta value determined by the Rayleigh-stochastic zone transition. This theory can explain w decreased with reducing grain size. This result indicates that the measurement is a go candidate for replacing the ′ measurement for characterizing State II because of its me urement simplicity [58]. The decreased grain size has been reported to result in a slig increase in owing to multiple scatterings and the mode conversion of the ultrason wave propagating through the grain boundary in a material [59]. Therefore, a correlati between the grain size and in State II was also demonstrated; however, its sensitiv was less than that of the other parameters.
When the LBS conditions deviated from State II, the combined effects of the gra size, texture, and porosity on the ultrasonic parameters were observed. Figure 11 presen the trend variations for each measured ultrasonic parameter as State II transitions to Sta I or III. Here, a label along the y-axis represents the change rate of the respective , , a ′ parameters, measured according to the 0.08 m/s LBS increment. The second-highest correlation of a in State II is shown in Figure 9. This observation is attributed to the ultrasonic attenuation at the grain boundaries, which has been extensively investigated [10,40]. An ultrasonic wave propagating through the grain boundaries scatters and its amplitude decreases according to the Rayleigh scattering theory, as follows [10]: where d is the grain size, f is the ultrasonic frequency, and superscript b is a constant value determined by the Rayleigh-stochastic zone transition. This theory can explain why a decreased with reducing grain size. This result indicates that the a measurement is a good candidate for replacing the β measurement for characterizing State II because of its measurement simplicity [58]. The decreased grain size has been reported to result in a slight increase in v owing to multiple scatterings and the mode conversion of the ultrasonic wave propagating through the grain boundary in a material [59]. Therefore, a correlation between the grain size and v in State II was also demonstrated; however, its sensitivity was less than that of the other parameters. When the LBS conditions deviated from State II, the combined effects of the grain size, texture, and porosity on the ultrasonic parameters were observed. Figure 11 presents the trend variations for each measured ultrasonic parameter as State II transitions to State I or III. Here, a label along the y-axis represents the change rate of the respective v, a, and β parameters, measured according to the 0.08 m/s LBS increment.
was less than that of the other parameters.
When the LBS conditions deviated from State II, the combined effects of the grain size, texture, and porosity on the ultrasonic parameters were observed. Figure 11 presents the trend variations for each measured ultrasonic parameter as State II transitions to State I or III. Here, a label along the y-axis represents the change rate of the respective , , and ′ parameters, measured according to the 0.08 m/s LBS increment. When State II transitions to State I (below 0.68 m/s LBS), the grain intensity along the 101 direction increased, and microporosity occurred. rate of ′ according to the 0.08 m/s LBS was the most prominent among eters, as shown in Figure 11c. The change rate of the ′ increment decre mately 1.2% and 1.5% at the transition of State I to II and in State I, resp pared to that of State II. This observation may be attributed to the overgeneration, the size of which is at the sub-micrometer level, as shown i material nonlinear elastic behaviors caused by the imperfect contact grains near the micropores can generate higher harmonics of the propa wave [18]. Furthermore, the closed interfaces of the micropores also gen higher harmonics [60]. This type of acoustic nonlinearity is not related to earity, such as the grain boundary effect, and results in a much stronger n material nonlinearity. Notably, although the porosity size is very small ( the micrometer level), it can degrade the mechanical properties of the co the and parameter results, no significant trend variations are ob State I and II because these parameters are less sensitive than ′ to m material of several micrometers; thus, these parameters have a stronger size variations [58]. If laser energy input greater than that for State I is AM process, the fabricated AM components may be overheated [55]. Un tions, not only microporosity but also small pits and bubbles correspon porosity in size (approximately above 100 μm) can be produced. Then, properties, including or rather than ′, can be effective in damage In State III (above 1.4 m/s LBS), the grain size and texture intensi direction decreased compared to that in State II, and macroscale LOF p Therefore, the ultrasonic parameters were affected to a greater extent by When State II transitions to State I (below 0.68 m/s LBS), the grain size and texture intensity along the 101 direction increased, and microporosity occurred. Here, the change rate of β according to the 0.08 m/s LBS was the most prominent among the three parameters, as shown in Figure 11c. The change rate of the β increment decreased by approximately 1.2% and 1.5% at the transition of State I to II and in State I, respectively, as compared to that of State II. This observation may be attributed to the over-melting porosity generation, the size of which is at the sub-micrometer level, as shown in Figure 5a. The material nonlinear elastic behaviors caused by the imperfect contact between several grains near the micropores can generate higher harmonics of the propagating ultrasonic wave [18]. Furthermore, the closed interfaces of the micropores also generate additional higher harmonics [60]. This type of acoustic nonlinearity is not related to material nonlinearity, such as the grain boundary effect, and results in a much stronger nonlinearity than material nonlinearity. Notably, although the porosity size is very small (at approximately the micrometer level), it can degrade the mechanical properties of the components [5]. In the v and a parameter results, no significant trend variations are observed between State I and II because these parameters are less sensitive than β to microvariations in material of several micrometers; thus, these parameters have a stronger relation to grain size variations [58]. If laser energy input greater than that for State I is used during the AM process, the fabricated AM components may be overheated [55]. Under these conditions, not only microporosity but also small pits and bubbles corresponding to the LOF porosity in size (approximately above 100 µm) can be produced. Then, linear ultrasonic properties, including v or a rather than β , can be effective in damage diagnosis.
In State III (above 1.4 m/s LBS), the grain size and texture intensity along the 101 direction decreased compared to that in State II, and macroscale LOF porosity occurred. Therefore, the ultrasonic parameters were affected to a greater extent by the LOF porosity than by the grain characteristics. In Figure 11a, a drastic decrease in the rate of change of v is observed, i.e., approximately 0.30% and 0.43% at the transition of State II to III and at State III, respectively, as compared to that in State II. This observation is attributed to the LOF porosity, which is in the approximate range of 100 to 200 µm and is much larger than the over-melting porosity (i.e., micropores), as shown in Figure 5c. A large number of LOF pores cause the air-filled cavities inside a component lower the elastic modulus of the material and result in a decreased v according to the following relationship: where E i is the elastic modulus obtained from the optimally manufactured component (i.e., porosity-free component), h is a constant value, V p is the porosity content, and ρ is the material density. This relationship indicates that the reduced v results obtained in our experiments are reasonable. In addition, a previous study has also reported that the macroscale porosity is related to an increase in a, which agrees with the observed increase in the increased a result (0.52% increase in the change rate) in State III, as shown in Figure 11b. This is attributed to the additional scattering at the boundary of the LOF porosity, the principle of which is similar to grain scattering (Equation (6)). The relationships between the linear ultrasonic properties, including v and a, and macroscale porosity have been reported in the literature [12]. Owing to the combined effects of the LOF porosity and grain size, no significant increase in a in the transition from State II to III may be observed.
Notably, although β is the best indicator in State II, monitoring only β may result in a false diagnosis of LOF porosity occurring in the transition zone between State II and III. This is because no notable trend fluctuation for β was observed in this transition zone, as shown in Figure 11c. This issue may be addressed via the simultaneous monitoring of β and v, which is the best indicator for this transition zone, as mentioned above.
In addition to the microstructural features mentioned above, the different precipitate and residual stress conditions for the tested samples can affect the varying ultrasonic parameters. However, in our experiments, no precipitates were observed in the 316 L stainless steel samples. Generally, precipitates in 316 L stainless steel can be generated in the form of (1) chromium carbides owing to the sensitization heat treatment or (2) delta ferrite under certain manufacturing conditions [22]. However, the samples used did not meet either of the above conditions. Previous studies [61,62] have reported that the residual stress during the AM process primarily depends on different 3D-printing strategies or varying building heights. In our experiments, the above factors were applied equally to all tested samples. Varying AM processing parameters such as laser power and LBS have been reported to contribute to the variation in residual stress conditions [63], which may cause certain experimental errors. To analyze the correlation among AM processing parameters, residual stress, and ultrasonic properties in detail, more detailed experimental designs are required in future studies. Note that when the LBSs are outside the conditions used in this study (from 0.6 m/s to 1.48 m/s), wavy surfaces and large pits were formed in the specimens where ultrasounds could not be detected well.

Conclusions
In this study, the microstructural characterization of additively manufactured 316 L stainless steel components was investigated using linear and nonlinear ultrasonic techniques. The contributions of this study are the presentation of (1) correlation analyses between various microstructural features (grain size and texture, LOF porosity, and microporosity) and linear/nonlinear ultrasonic properties (ultrasonic velocity, attenuation, and nonlinearity parameters), (2) the development of NDT for the microstructural characterization of additively manufactured components, and (3) the potential for the online monitoring of AM processes owing to its nondestructive nature. The key research outcomes of this study are as follows in Table 2. The various microstructural features with varying LBSs were characterized using three ultrasonic parameters: v, a, and β . State I (low LBS condition: below 0.68 m/s) could be effectively characterized using β because of the higher harmonics generated owing to the microporosity at levels of several micrometers. The β parameter was also effective in characterizing State II (medium LBS condition: 0.92-1.16 m/s) because of the variations in texture, grain boundary, and dislocation density effects. The ultrasonic scattering attenuation at the grain boundary also allowed parameter a to be used for characterizing State II. The generation of the macroscale LOF porosity (approximately greater than 100 µm) in State III (high LBS condition: above 1.4 m/s) showed a good correlation with the v parameter, which was directly related to the material elastic modulus and density.
Further research is required to develop a novel method for the simultaneous monitoring of β and v to estimate various microstructural features and use laser ultrasonics for microstructural characterization, which can be conducted in a contactless and nondestructive manner. In addition to the LBS conditions, microstructural variations caused by other process parameters including laser power and hatch space will also be considered in future studies.