Simulation-Based Investigation of a Probability of Detection (POD) Model Using Phased Array Ultrasonic Testing (PAUT) Technique

Probability of detection (POD) as a metric for quantifying the capability of inspection procedures in nondestructive evaluation (NDE), has been applied and evolved in industries since 1970s. Progress had been noted when certain statistical functions were brought up to model POD behavior, including log-normal model (also referred as Probit model). This model had been concluded to be the best fit and therefore has been widely used in many studies, while the involved assumptions and conditions were not carefully addressed and explored. To make flexible POD datasets available for specific inspection procedures and reduce the number of expensive experiments needed, model-assisted POD (MAPOD) is an alternative. This paper addresses a pure simulation-based POD procedure of an inspection scenario involving phased array ultrasonic testing (PAUT) on lack-of-fusion defects in additive manufactured (AM) components. The mathematical simulations are performed by an ultrasonic testing (UT) simulation software, simSUNDT, developed at Chalmers University of Technology in Sweden. Resulted inspection datasets with the proposed data processing steps are evaluated in terms of the assumptions and conditions of log-normal POD model, with the purpose of discussing the POD model validity under different circumstances. Simulation-based POD curves are finally compared with several discrete POD values at some defect sizes, calculated through massive computations from physics-model based metamodel. Comparisons and observations confirm satisfactory application of log-normal POD model despite some violations in model hypotheses.

assets and assuring of structural integrity. However, consistent results are not always necessarily guaranteed through repeated inspections on a specific defect due to the probabilistic characterization of inspection capability. This makes it much essential to quantify the performance and reliability of a nondestructive evaluation (NDE) procedure, especially when risk-based inspection (RBI) methodologies and component life-cycle estimations are introduced for in-service property management. Probability of detection (POD), as a statistic metric initiated since 1970s from aero industry [1,2], has been a tool used for assessing the reliability of NDE procedures. It helps describe the accuracy of an inspection and reveal how well an inspection procedure can detect vital defects. Inspection data is in this method transformed into a relation between the probability of defects detection and a characteristic parameter of the defects [3], usually the defect size [4], as it is crucial for structure integrity and estimation of component life. A constant POD value was proposed at the beginning for all defect types of the same size, while an analysis of a study showed discrepancies in repeated inspections on a certain crack and on different cracks of the same size [5]. In fact, there are many influential factors that may introduce variations when determining the POD values, including test methods, equipment and material conditions, defect features as well as human factors. It is therefore important to understand the process of obtaining a POD value and to question the validity of its application and limitations. The experimental data required to determine the POD value can be collected from for example field service records, artificially created defects in components and generic test blocks [6]. Experiments are repeated in many runs following a controlled condition specified by procedure to record several results. These results can be categorized in hit/miss type if binary conclusions are drawn based on certain detection criteria, or in signal response type if result amplitudes are of interest. The interpretation of POD value for a defect size is then the proportion of times the defect being detected (hit) among all trials. A POD curve is approximated with best fit through all POD values among a range of defect sizes using statistical analysis models, e.g., log-odds and log-normal model [7], also referred as Probit model [6]. Obviously, the experiment processes can provide practical measures of the capability of inspection instruments and the actual environments. However, these practical processes are always expensive and resource consuming, and even more challenging when introducing representative features to the artificial defects. It might also end up in poor statistics and the results could show scatters [8]. Besides, the experiments tend to have bias because the trials are conducted with defective test pieces, which could be different from the reality [9]. Parametric studies for investigating influential parameters in procedures are also unfeasible by practical trials.
With the development of physics-based nondestructive testing (NDT) model in recent decades, mathematical simulations are widely applied in many related applications, among which assisting or replacing part of the experimental works in POD estimations is included [8,10]. There is also a European project PICASSO exploring the possibility of building simulation-supported POD curves [11]. This model-assisted probability of detection (MAPOD) approach, realized by either transfer function [12] or using full model [13], has the advantages of being resource-and time-saving. It also provides possibilities of parametric studies of performance and reliability. Models can also be utilized as training tools for inspectors and be part of the qualification and validation processes in NDT applications. Human factor, the most influential parameter within the manual experimental processes, can to a large extent be isolated, but can also be studied and approximated by modelling approach [14,15]. Along with the formation of MAPOD Working Group at Center for Nondestructive Evaluation of IOWA State University in 2005 [13,16], there are more groups of researchers devoted into MAPOD research and development. The National NDT Center of UK is an example, now as part of ESR Technology [17]. There has also been a group of researchers in Netherlands developing a so-called POD-generator [18], in which 3 individual modules are contained and combined to simulate the entire POD estimation process from defect initiation and growth to the prediction of failure probabilities. A commercial NDT simulation software, CIVA, has also been developed for MAPOD purposes with corresponding module [19,20]. Another ultrasonic testing (UT) simulation software, simSUNDT, developed at Chalmers University of Technology in Sweden, has undergone an investigation regarding POD estimation in terms of conventional UT inspection according to a testing procedure UT-01 for pipes [21]. Its mathematic kernel, UTDefect [22][23][24][25], which performs the actual UT modeling and calculation, can be further applied independently in various parametric studies, such as sound field optimizations for phased array ultrasonic testing (PAUT) technique in terms of maximizing the received echo amplitude [26].
In this paper, the above-mentioned mathematic kernel for UT simulation, UTDefect, is used for generating UT inspection signal response data for pure simulation-based POD estimation. The studied inspection scenario is based on PAUT on lack-of-fusion defect, which is plausible in components produced by some additive manufacturing (AM) processes (e.g., laser metal deposition, LMD). The results are transformed into corresponding POD curve using the widely accepted log-normal POD model (Probit). The POD curve is then compared with a set of discrete POD value points obtained through a large number of metamodel-based simulations, to verify the validity of the POD model. In addition, a modified distance amplitude correction and time varied gain (DAC/TVG-mod) framework is considered afterwards to contribute to simulation results and POD curves for further investigation and comparison.

simSUNDT Software
The simSUNDT software consists of Windows®-based preand post-processor as well as a mathematic kernel, UTDefect, for the actual modeling and calculation, which has been experimentally validated [23][24][25]27] to some extent by comparing with available experimental data. This mathematic kernel takes specified text file with input parameter values and generates output information in plain text format, so that external programs can call this kernel for automatic operations with certain cases. To model the probes and the interaction with defects (scattering), a series of integral transforms and integral equations are employed. Together with the possibility of calibration, i.e., towards a reference reflector including side-drilled hole (SDH) and flat-bottomed hole (FBH), the software can simulate the entire testing procedure under different testing scenarios.
The model is fully three-dimensional, while the simulated component in the software at this stage is limited to be an infinite plate with finite or infinite thickness and bounded by the scanning surface. The probe is modelled as a boundary condition on the surface of an elastic half-space. This surface is traction free except the active part of the contact area beneath the probe. This enables the flexibility of probe simulation regarding its shape, wave type, element size, angles, etc. The receiver in the UT system is modelled by using a widely known and used reciprocity argument [28]. The available defect types are volumetric and crack-like defects, and the methods of solution towards the defect modeling is the T-matrix method [29], where the transition matrix includes all information regarding the defects.
For the phased array (PA) probe technique, each individual element follows the principle of probe modelling above. The individual boundary conditions are then translated into the main coordinate system, so that a PA wave front with certain beam angle is formulated by constructive phase interference. In order to enable larger beam angles, a model of a wedge is also included as an option. This PA probe model has been experimentally validated qualitatively [26] and quantitatively [30].

Simulation-Based Probability of Detection
POD data from physical experiments are usually retrieved through repeating many inspection trials for a set of defects under a controlled condition specified by corresponding procedure. The POD value for a defect size is interpreted as the proportion of times this defect is marked as detected among all trials on this defect. As said in Sect. 1, experimental trials can provide measures of the capability of used inspection instruments and reflect on inspection surroundings, but they are also expensive and challenging because large amounts of specimens containing a specific defect type are needed for computing valid statistical parameters of POD function. The number of defects recommended for hit/miss type of data is at least 60 and their sizes should be uniformly distributed between the minimum and maximum defect size of interest. At least 30 defects should be available for signal response type of data thanks to added information in response signal [7,31]. This type of data is sometimes referred as "a" vs. "â" data, where "a" stands for defect size and "â" for quantified signal response, recorded in terms of a parameter. Yet, the statistics might still be poor from this limited number of available defects and shows scatters in the results.
Physics-based mathematical simulation models compared to physical experiments of inspection, have the superiority of reducing the massive cost and extensive operations. The number of defects is not a limitation by simulations, which in turn facilitates satisfaction of appropriate statistical sampling requirement in design of experiments (DOE).
Since observed from physical experiments of probabilistic nature of NDE processes, repeated inspections of a certain defect will not necessarily produce the same signal responses due to minor variations in some inspection parameters, e.g., experimental setup and calibration. Variations could also come from minor differences of some defect-features related inspection parameters despite that the defect sizes, the characteristic parameter of the defects, against which the POD curve will be plotted, are the same. One way of representing these variations in simulation models is to use the uncertainty propagation method [32], where the inspection parameters are specified within individual uncertainty ranges (parameter space) and corresponding distributions. A series of simulations are then performed in the parameter space and the resulted signal responses are transformed into a POD curve using log-normal model, which was concluded the best fit among other assessed models [6]. This model is based on a linear relation between defect size (a) and signal response ( a): where β 0 and β 1 are the intercept and slope, and δ is a random error term addressing the difference between observed and estimated signal response. Berens [7] assumes that δ is normally distributed with mean value of zero and constant standard deviation σ δ , which is independent on defect size: δ ∼ N (0, σ δ ). The POD curve generated by log-normal model is finally expressed in a form: where denotes standard normal distribution function, and: a dec in Eq. (3) stands for a specified decision threshold for signal response. The defect is marked detected only if its signal response exceeds this threshold value. Upon the POD curve is obtained, a large number of additional simulations are continued in this work for a defect size to determine its detection percentage, as a specific estimate of POD for this defect size. This process is repeated for several defect sizes to acquire multiple discrete POD value points, aiming at comparing the modelled POD curve with specific In addition, this work also involves a modified distance amplitude correction and time varied gain (DAC/TVG-mod) framework to contribute to simulation results as a postprocessing step. The POD curve and values are also added with these contributions in evaluation.

Studied Inspection Scenario for Evaluation
The inspection scenario studied in this paper is PAUT on lack-of-fusion defects in AM components, produced by e.g., laser metal deposition (LMD). The component is made of material with acoustic properties of 5573 m/s and 3150 m/s for longitudinal and transverse wave speed, respectively. To simplify the result obtaining process, factors such as material attenuations, noises etc. are excluded from this study. Within the scope and availability of the UT simulation software simSUNDT, penny-shaped circular crack [33] is chosen as representative of the actual defect in simulations. A general appearance and parameter definitions of a circular crack is shown in Fig. 1, where: -the rotation axis of the circular crack passes through the center of the crack θ is the crack tilt angle (the angle between the crack rotation axis and the vertical z-axis perpendicular to scanning surface) -φ is the crack skew angle (the angle between the probe scan line along x-axis and the projection of the crack rotation axis onto the scanning surface) Relative bandwidth (%) 74 Figure 1 also shows the inspection setup considered in this work, i.e., contact test using PA probe under pulse-echo (one probe only) situation. Simulated UT instruments are the ones used in previous validations [26,30]. This includes a PAUT device named Topaz64 [34] from Zetec as data acquisition hardware, and a longitudinal wave linear PA probe labelled LM-5MHz [35] from Zetec. Key information of this PA probe is summarized in Table 1.
To account for probabilistic characteristics of inspection capability, the inspection parameters in this studied case are formulated by uncertainty propagation as follows.
Using the PA probe with proper delay law, the generated sound beam inside the test component should have a nominal beam angle of 0°(α in Fig. 1) with uncertainty of ± 2°. The nominal probe skew angle of 0°is set, which is along the bead direction of AM component, with uncertainty of ± 5°due to the play in physical situation. Test volume addressed includes defects with nominal depth of 25 mm with uncertainty of ± 5 mm. Focusing effect of the PA probe is used so that most wave energies could be reflected from the test volume. Each defect has a nominal tilt and skew angle of 0°but is specified with uncertainty of ± 10°and ± 90°, respectively. The diameters of circular cracks, taken as the defect size (characteristic defect parameter), are limited between 0.5 and 5 mm. Table 2 summarizes these preliminary inspection parameters as well as their distributed ranges. The maximum echo amplitudes from all simulations are taken as the signal response. To assess if the defect is detected, a decision threshold ( a dec ) of − 6 dB drop from a calibration defect is set. The calibration defect in this work is selected as a SDH with 0.5 mm in diameter and 25 mm in depth.

Sensitivity Analysis of Inspection Parameters
Practically, a number of inspection parameters are evaluated in terms of their impact on the resulted signal response amplitude, aiming at reducing the number of involved factors in analysis [21]. Identified preliminary inspection parameters in Table 2 thus undergo a sensitivity analysis. UT simulation kernel, UTDefect, is incorporated and launched through a simulation process automation software, mod-eFrontier, to generate and analyze a series of simulation results. For the purpose of sensitivity analysis, multi-level full factorial sequence of DOE is adopted in modeFrontier for all inspection parameters according to Table 2, except for the characteristic parameter (defect size). This so-called multi-level full factorial sequence of DOE helps sample the involved inspection parameters within each parameter range at multiple levels (values), and plans the experiments or simulations at all possible combinations of these values. It contributes to a comprehensive study of the effect of each parameter on the final response. When all simulations are completed, sensitivity analysis is performed in modeFrontier using Smoothing Spline ANOVA (SS-ANOVA) proprietary algorithm [36], which ends up in relative contribution index of each selected inspection parameter to the simulated maximum echo amplitude. Analysis shows that the inspection parameters can be ordered according to their relative contribution indexes, where larger index indicates higher impact of inspection parameter to the echo amplitude. Note these contribution indexes in parenthesis indicate only relative impact of each parameter to the result under this analyzed situation: Contribution index for beam angle and probe skew angle is almost 0 from the analysis. These two inspection parameters are therefore set constant at their nominal value of 0°i n simulation parameter space of the study. Other inspection parameters are then treated as essential parameters with proper uncertainties, summarized in Table 3. Though as noted on focusing depth with the highest relative contribution index, it is however treated as constant and set at the nominal depth of test volume, i.e., 25 mm. This could be accomplished and assured in practice by performing corresponding calibration beforehand to exclude the impact of beam focusing effect. Defect sizes are evenly distributed within the range between 0.5 and 5 mm.

Simulation Scheme
As in physical experiments, operators tend to have a continuous exploration over the scanning surface to find where the potential signal response from the defect of interest could appear, then a local search proceeds around this limited region to find the ultimate maximum signal response. It is especially the case to search in a larger region when the defect features are not witting beforehand because these unknown features can change the position of the ultimate maximum signal response. For mathematical simulations, this overall idea and strategy can to a large extent help avoid inefficient and unintended simulation steps within a complete mesh grid with discrete computation points, with the purpose of capturing the potential maximum signal response, as operators behave in physical experiments. To realize this strategy, a so-called macro-micro scanning pattern is developed and implemented in simulation workflow. This pattern enables that a simulation firstly is performed with a coarse mesh setup with mesh size of s x1 and s y1 in x-and y-direction, respectively, over the complete scanning area to find a global maximum of signal response, i.e., echo amplitude a 1 in this work. Based on the location (x 1 , y 1 ) of this global maximum, the next scanning region shrinks to x ∈ [x 1 − s x1 , x 1 + s x1 ] and y ∈ [y 1 − s y1 , y 1 + s y1 ], and new mesh size is decreased 5 times from previous size, i.e., s x2 s x1 /5 and s y2 s y1 /5. This process is repeated until the mesh size is decreased below 0.1 mm. The last mesh size is restricted to 0.1 mm because it already gives convergent results from UTDefect kernel. An illustration of the process is shown in Fig. 2.  It is in the end implemented as an automation workflow in modeFrontier as shown in Fig. 3 for the specific study case in this paper, where s x1 s y1 s 1 1mm (corresponds to "inc1" in the figure) and s x2 s y2 s 1 /5 s 2 0.2mm ("inc0.2" in the figure). After a UTDefect simulation, a MATLAB script is used in the workflow for generating the new mesh sizes and updating the scanning area based on the resulted C-scan file. This helps save enormous simulation runtime. As an example of a defect with diameter of 5 mm, simulation with mesh size of 0.1 mm within the entire scanning region takes around 50 h for just one case, but is surprisingly reduced down to 45 min using the macro-micro scanning pattern. This makes heavy computations more feasible.

Simulation Base for Metamodel Training
To get statistically sufficient number of inspection results for estimating a POD value of a certain defect size, a method of response surface model (RSM), in other words metamodel, is to be used. The metamodel is an approximate mathematical representation of complex model responses, constructed or trained based on a selection of time-consuming simulation results from original mathematical model. Metamodel could then provide efficient and fast estimations of the model response instead of calculating through the original mathematical model.
To construct a valid metamodel for the current study, a multi-level full factorial DOE sequence is applied to gener-ate a total of 625 inspection cases in modeFrontier, within the parameter space specified in Table 3. Simulation of these DOEs using UTDefect are sequentially and automatically launched by modeFrontier with application of the developed macro-micro scanning pattern. A metamodel with defect size, depth, tilt and skew angle as input parameters and maximum echo amplitude as output value is finally built using Duchon's Polyharmonic Spline based radial basis function [37]. This metamodel correlates well with results from UTDefect simulations in verification step, with mean absolute error of 0.5 dB and mean relative error of 0.6%. The output from this metamodel is called virtual simulations in the following context.

Modified Distance Amplitude Correction and Time Varied Gain (DAC/TVG) Function
Practically, DAC/TVG function is applied in inspections to compensate for material attenuation and beam spreading etc., so that the differences of signal amplitudes due to sound wave travel distance could be reduced to some extent. This function is additionally considered and evaluated in this work in a certain way, modified as followed, to examine its influence on final POD results. A series of reference reflectors, usually FBHs, of the same size at different depth are commonly used in practice to generate DAC/TVG curve. A DAC curve plots variations of signal amplitudes of these reflectors as a function of their depth, and a corresponding TVG presents compensating gains for these reflectors such that all signal amplitudes are brought to the same height, e.g., normally 80% screen height of the instrument, see an illustration in Fig. 4.
Considering the current study case and simulation results from UTDefect, a modified DAC/TVG function framework is proposed, denoted as DAC/TVG-mod. Two reference blocks with the same material properties as the current study case are assumed. Reference block 1 contains 11 FBHs of size 0.5 mm as reference reflectors, situating at depth from 20 to 30 mm with increment of 1 mm, covering the entire test volume of the study case. Block 2 as an extra complement, is the same as Block 1 but the FBH sizes are 5 mm instead. For each reference block, simulations using the same inspection setup as the study case, i.e., 0 degree beam angle and probe skew angle with a fixed beam focusing depth of 25 mm, are performed on all contained FBHs. By the availability of defect types in simSUNDT, the FBHs are well-represented by open circular cracks in simulations, which are experimentally verified. Their maximum echo amplitudes are recorded as signal responses, where the largest response value from each block is then treated as the "standard level", A s . Amplitude differences from other reflectors on the same block would be the corresponding gains at individual depth and the gain value

Modified DAC/TVG Function
As mentioned in Sect. 5.4, two reference blocks are proposed in this work with FBHs of different sizes. Simulations using UTDefect provide echo amplitudes for DAC curves, see Table 4 and Fig. 5. Note that the resulted amplitude from UTDefect is expressed in log scale in terms of decibel (dB), which has no practical meaning unless it is calibrated. The values with superscript "a" in Table 4 indicate the "standard level", A s , for corresponding reference block.
Corresponding TVG are retrieved through amplitude differences to the "standard level", A s , at each reflector depth, see Table 5 and Fig. 6. Required gain in this work for a defect depth is thus available through these results by interpolation.

Simulation Results Processing
With the help of the developed metamodel in this work, 20 inspection cases for each defect size are conducted, among which the essential parameters of the inspection cases are  Table 3. In total, 30 defect sizes are evenly distributed within the formulated defect size range, i.e., from 0.5 to 5 mm, which means 600 virtual simulations in total are computed. Figure 7 plots the distribution of these calibrated echo amplitudes against the defect sizes. Based on report [38], extra attention should be focused on evaluating the performance of the default log-log transformation of defect size (a) and signal response ( a) in log-normal POD model. This is to ensure the linear relation between these two quantities, assumed by Berens [7] in Eq. (1). The report points out that log-log scale is not always the best transformation of this linear relation, but either linear or log scale of these quantities is to be considered for a combination that best fits the current dataset, to satisfy the basis of log-normal POD model. In this study, the calibrated echo amplitude results from simulations are already in unit of decibel (dB) as in log scale. It is for this reason that only defect sizes are taken in both log and linear scales in Fig. 7 for assessment. Additionally, to correctly apply the log-normal POD model, assumptions of normal distribution and constant standard deviation for the random It is observed from Fig. 7 (left) that echo amplitude distribution under log scale of defect size can be better approximated using a straight line, and the assumptions of modelling error are also satisfied. Thus, using log-scale of defect size could ensure the validity of applying log-normal POD model for this dataset. It is also noticed that a data point from a large defect size (with circle) seems to stay out of the normal distributed region of other points. This will be addressed and discussed later.
Taking the DAC/TVG-mod function into consideration using the proposed two reference blocks, resulted echo amplitudes are hereby called TVG-compensated echo amplitudes in this paper. An example addressing 5 different defect sizes (0.5 mm, 1.6 mm, 2.8 mm, 3.9 mm and 5 mm) and 20 inspection cases per size is shown in Fig. 8 using TVG from two blocks. Note that these data from different blocks are plotted with a little shift in abscissa for clear visualization purpose only.
It is noted in Fig. 8 that reference Block 1 (with FBHs of size 0.5 mm) helps bring the calibrated echo amplitudes from defects of size 0.5 mm (Def. 1) to a very similar level (the Table 5 Corresponding TVG (in dB) based on DAC and "standard level" for two proposed reference blocks Reflector depth (mm) 20   circle signs of Def. 1 "shrink" very much from corresponding plus signs), as the purpose of reducing the impact of sound wave travel distance (i.e., defect depth in this work) using DAC/TVG-mod function. However, Block 1 helps little when the defect sizes deviate from 0.5 mm (the circle signs of Def. 5 do not "shrink" that much from corresponding plus signs). Reference block 2 (with FBHs of size 5 mm) reveals similar behavior as Block 1 on the same defect size (Def. 5), but it also reduces the impact of defect depth on much different defect sizes (Def. 1), which differ much to FBH sizes of 5 mm. This observation is further confirmed by sensitivity analysis, investigating the impact of defect depth to resulted echo amplitudes with or without using Block 1 and 2. The relative contribution index of defect depth from each scenario can be ordered in parenthesis as below, which indicates that DAC/TVG-mod function by Block 2 helps better in reducing the impact of defect depth to echo amplitudes: No DAC/TVG-mod function (0.144) > With DAC/TVG-mod using Block 1 (0.052) > With DAC/TVG-mod using Block 2 (0.003) Applying these DAC/TVG-mod functions to original dataset in Fig. 7 gives distribution of TVG-compensated echo amplitudes, shown in Figs. 9 and 10 by using Block 1 and 2, respectively. Log scale of defect size (figures to the left) could still present better modelling behavior of linearity comparing to linear scale (figures to the right). This is consistent as concluded from Fig. 7, thus is employed in corresponding modelling process of POD curves. It is however noted from these log-scaled figures when it comes to other hypotheses of POD log-normal model, that the standard deviation of resulted data is not uniform about the estimated straight line, which represents the linear relation in Eq. (1). In other words, the standard deviation of random errors is dependent on defect sizes when DAC/TVG-mod function is involved. There are also some data points from large defect sizes fall outside of the normal distributed region along the regression line. In convenience of discussion for some of them, we number the data points No.1-5 in Fig. 10 (left), where data point No.1-4 fall above the decision threshold and point No.5 falls below the threshold. This point No.5 is also the circled data point indicated in Fig. 7. The combination of the inspection parameters of these resulted data points are listed in Table  6. It can be noticed that these data points come from the defects with large tilt and skew angle. However, it is inappropriate to conclude that a combination of large tilt and skew angle of defect will give a weak echo amplitude if this DAC/TVG-mod function is applied. In other words, if the DAC/TVG-mod function is applied, then the large tilt and skew angle of defect is the necessary and insufficient condition of having a weak echo amplitude. This is because a weak echo amplitude could be compensated with a high TVG depending on the defect depth. It is hereby a joint effect of these inspection parameters. These specific cases of No.1-5, especially the No.5 data point could however question the resulted POD curve when the log-normal POD model is fitted in upcoming work.

Estimation of Model Parameters and POD Curves
As seen from Eq. (2), log-normal POD model is controlled by parameters β 0 , β 1 and σ δ . For the dataset of calibrated Fig. 9 Distribution of defect echo amplitudes against defect sizes of 600 simulation cases (30 defect sizes with 20 inspections per size) in log (left) and linear (right) scale of defect size (abscissa). DAC/TVG-mod function is applied using reference block 1 and decision threshold ( a dec ) is set at -6 dB after calibrating with a SDH Fig. 10 Distribution of defect echo amplitudes against defect sizes of 600 simulation cases (30 defect sizes with 20 inspections per size) in log (left) and linear (right) scale of defect size (abscissa). DAC/TVG-mod function is applied using reference block 2 and decision threshold ( a dec ) is set at -6 dB after calibrating with a SDH Table 6 Combination of inspection parameters for the five resulted data points in Fig. 10 No.
Defect echo amplitudes without and with TVG-compensation using reference block 1 and 2 (figures to the left in Figs. 7, 9 and 10, respectively), the estimated model parameters are summarized in Table 7. POD functions for these datasets using parameters in Table 7 are then expressed as below in form of Eq.
(2) and plotted in Fig. 11. Defect size of 90% POD with 95% confidence, a 90/95 , is 3.6 mm, 1.6 mm and 1.8 mm for three respective cases, which clearly indicates from Fig. 11 that POD is improved when DAC/TVG-mod function is in use (case 2 and case 3) under the same decision threshold of − 6 dB, i.e., smaller defect sizes could have a higher POD after TVG-compensation. This is because echo amplitudes from all defects are compensated with certain gains depending on their depth. The difference of resulted POD curves between case 2 and case 3 in Fig. 11 comes from the level of compensated gains by two reference blocks, as seen in Fig. 6. These gains help some echo amplitudes reach the detection threshold.
As concerned previously in Sect. 6.2, the resulted POD curves in Fig. 11 based on the estimated parameters in Table  7 could be doubtful due to non-uniform standard deviation about the estimated straight line and some data points outside of normal distributed region of most data points. This is to be further assessed and compared by estimating discrete POD value points at some defect sizes. A discrete POD value point for a defect size, according to POD definition, is obtained through performing a series of inspections on defects of this size and counting the proportion of times this defect size being detected among all trials. Taking the advantage of constructed metamodel, 5000 virtual simulations for a defect size can be rapidly accomplished. To account for a certain defect size range, 24 defect sizes are included in this investigation. The same detection criteria of − 6 dB after calibration as in POD curve model is used for these virtual simulations. These POD value points for three investigated cases are plotted with original POD curves in Fig. 12.
In general, Fig. 12 shows good correlations in trend between POD curve and discrete POD value points. The POD curves based on log-normal POD model underestimate the POD for defect sizes larger than about 1.3 mm for all cases. Defect sizes smaller than this have limited POD as seen from discrete values, while POD curve could prescribe a higher probability.
Specifically, it should be pointed out that not all discrete POD value points in Fig. 12 converge to 1 above a certain value of defect size. For example, the points of defect size from 3 to 5 mm have values between 0.95 and 0.98 for case 1, between 0.9992 and 1 for the points of defect size above 1.5 mm for case 2, and between 0.9994 and 1 for the points Fig. 12 Comparison between POD curve and corresponding discrete POD value points at some defect sizes of defect size above 1.7 mm for case 3. Though these POD values are very close to and could be treated as 1, they however indicate that there are still a small number of cases with combination of inspection parameters that can give resulted echo amplitudes below the decision threshold level. Part of these small number of echo amplitudes could also fall outside of the normal distributed region of most other results, as concerned in Sect. 6.2. These small number of cases, as being examined, come from the combination of large tilt and skew angles of defect, similar to earlier discussions. It should yet be emphasized again that the large tilt and skew angle of defect is only the necessary but insufficient condition of having a weak echo amplitude, if the DAC/TVG-mod function is applied. Now even if there are some echo amplitudes fall outside of the normal distributed region, which violate the corresponding assumption of log-normal POD model, we should still note that the number of these cases are very limited among 5000 cases. The comparisons in Fig. 12 between the discrete POD value points and corresponding POD curve show no need for this concern of violation of model assumptions.

Conclusion
This paper addresses a process of generating datasets for quantifying NDE detection capacity in terms of POD, of certain ultrasonic inspection scenario based on numerical simulation model. A set of inspection parameters for this specific inspection scenario undergo a sensitivity analysis to determine essential parameters in parameter space. These essential parameters, which have greater impact on the inspection results, are specified with reasonable uncertainties and distributions based on uncertainty propagation method. A series of well-organized UT simulations are performed in the parameter space using a proposed efficient macro-micro simulation scheme, which saves numerous runtimes of simulations. By these simulated results, a metamodel is constructed with the aim of providing fast computations, instead of running the original model. This can help perform thousands of virtual simulations rapidly for estimating statistically valid POD values. In the end, 30 defect sizes with 20 inspection cases each have been computed within the parameter space and they are transformed into POD curve using log-normal POD model. A modified DAC/TVG function is also applied on this dataset and improvement of POD results are seen from the new curves. Despite the fact that some of Berens hypotheses of log-normal POD model are not perfectly satisfied in datasets with DAC/TVG-mod compensation, discrete POD values from some defect sizes could still prove the log-normal function as a satisfactory model to signal response type of inspection datasets.