Simulation of Real Defect Geometry and Its Detection Using Passive Magnetic Inspection (PMI) Method

: Reinforced concrete is the most commonly used material in urban, road, and industrial structures. Quantifying the condition of the reinforcing steel can help manage the human and ﬁnancial risks that arise from unexpected reinforced concrete structure functional failure. Also, a quantitative time history of reinforcing steel condition can be used to make decisions on rehabilitation, decommissioning, or replacement. The self-magnetic behavior of ferromagnetic materials is useful for quantitative condition assessment. In this study, a ferromagnetic rebar with artiﬁcial defects was scanned by a three-dimensional (3D) laser scanner. The obtained point cloud was imported as a real geometry to a ﬁnite element software platform; its self-magnetic behavior was then simulated under the inﬂuence of Earth’s magnetic ﬁeld. The various passive magnetic parameters that can be measured were reviewed for different conditions. Statistical studies showed that 0.76% of the simulation-obtained data of the rebar surface was related to the defect locations. Additionally, acceptable coincidences were conﬁrmed between the magnetic properties from numerical simulation and from experimental outputs, most noticeably at hole locations.


Introduction
Reinforced concrete as a composite infrastructure material is widely used in construction because of its excellent properties [1] and construction ease. Three factors control the behavioral responses of reinforced concrete: the reinforcing steel (generically referred to as rebar in this article), which has a noticeable ductile nature; the concrete itself, which has a noticeable brittle nature (low tensile strength but high compressive strength); and the condition of the rebar-concrete bonding (to achieve reliable stress transfer) [2].
Reinforced concrete is commonly used in infrastructure such as buildings, bridges, and highway construction [3]. The quality of a country's transportation system is mostly based on the conditions of its highway bridges, all of which contain steel. At the present time, apparently, approximately 28% of concrete bridge decks in the United States (US) and 33% of highway bridges in Canada can actually be considered operationally deficient or in a condition warranting the cessation of active service, mainly because of rebar corrosion [4].
Rebar corrosion is common in environmentally-exposed structures; it reduces the service life of these structures and impacts load-carrying capacity [5]. In the worst cases, structural failure occurs because corrosion reduces a stressed rebar's cross-sectional area [6] to the point of rupture. Rebar corrosion also degrades the bonding quality and can create cracks in the structure from Based on potential theory, volume charge density (λ V )-Equation (3)-and surface charge density (σ S )-equations (4) and (5)-are other parameters related to magnetization (M)-Equation (6)-and can be implemented for computing stray fields. Surface charge density is calculated by Equation (4) when there is just one magnetic medium; Equation (5) M(r) = J(r)/J s (6) According to Equation (7), the stray field energy at a position (r) can be also calculated through the negative gradient of the potential of the stray field energy at a place (Φ d (r)) [55], where Φ d (r)-Equation (8)-is a function of magnetization saturation (J s ), volume charge density (λ V ), surface charge density (σ S ), and the derivative of the position vector (r ). Next, the magnetic field energy is obtained from Equation (9) through the integration functions of surface charge density and volume charge density over the volume and surface, respectively. (9) For conducting this research article, we scanned 373.87 mm of the surface of a ferromagnetic rebar (low-carbon steel), with a diameter of 16 mm, and two artificial defects (Table 1) [37], using a high-resolution 3D laser scanner (Figure 1a) [56]. The shape of the rebar was created with cloud points (Figure 1b) that were modified and converted to a mesh by Mesh Lab V1.3.2 (http://meshlab. sourceforge.net/). Subsequently, the produced mesh was imported to COMSOL R software and converted to a discretized surface and solid, respectively ( Figure 1c). The solid rebar was simulated via COMSOL R software with regard to the magnetic field of the Earth, different components of magnetic flux density were investigated at different spacing, related simulation results were compared with our previous experimental results, and statistical approaches were introduced. According to Equation (7), the stray field energy at a position (r) can be also calculated through the negative gradient of the potential of the stray field energy at a place ( ( )) [55], where ( )-Equation (8)-is a function of magnetization saturation ( ), volume charge density ( ), surface charge density ( ), and the derivative of the position vector (ˊ). Next, the magnetic field energy is obtained from Equation (9) through the integration functions of surface charge density and volume charge density over the volume and surface, respectively.
For conducting this research article, we scanned 373.87 mm of the surface of a ferromagnetic rebar (low-carbon steel), with a diameter of 16 mm, and two artificial defects (Table 1) [37], using a high-resolution 3D laser scanner (Figure 1a) [56]. The shape of the rebar was created with cloud points (Figure 1b) that were modified and converted to a mesh by Mesh Lab V1.3.2 (http://meshlab.sourceforge.net/). Subsequently, the produced mesh was imported to COMSOL R software and converted to a discretized surface and solid, respectively ( Figure 1c). The solid rebar was simulated via COMSOL R software with regard to the magnetic field of the Earth, different components of magnetic flux density were investigated at different spacing, related simulation results were compared with our previous experimental results, and statistical approaches were introduced.

Simulations and Results
After converting the rebar mesh to solid in COMSOL R software, the magnetic behavior simulation was undertaken. Considering the variations of Earth's magnetic field in time and location, to obtain consistent and realistic results the average (within a year) of the different components of the magnetic field for the Waterloo, Ontario region (the location of the experiments) was adopted for the simulations (Table 2). Moreover, since the unitless relative magnetic permeability of low-carbon steels (ASTM 1020) range from 50 to 100 [57,58], a relative magnetic permeability of 75 was selected for this study. The duration of exposure to an external magnetic field will affect the magnetic behavior of ferromagnetic materials. In reality, ferromagnetic materials are affected by the magnetic field of the Earth from the beginning of their production process. There may also be some unknown external magnetic sources in the surrounding environment that affect the magnetic behavior of ferromagnetic objects [59]. However, as accurately as possible, we can apply the magnetic field of the Earth to the object and simulate its magnetic behavior, although some divergence will exist between the simulation and the experimental results.
To consider the Earth's magnetic field in the simulation, the rebar was located in a regular space ( Figure 2) with dimensions of 100 mm × 150 mm × 410 mm, which included the magnetic field presented in Table 2 and Figure 3. To have better control of simulation parameters, the box and rebar were meshed separately with tetrahedral meshes according to the specifications in Table 3 ( Figure 4a,b). Then, the rebar and box were jointly subjected to the simulation process as a single system (Figure 4c). The values of the different components (X, Y, and Z) of the magnetic flux densities were recorded for the Y direction of the rebar (i.e., the path parallel to the rebar's length). This path is at the surface of the rebar, and extends from one side (Edge A) to the other side of the box (Edge B) ( Figure 5).

Simulations and Results
After converting the rebar mesh to solid in COMSOL R software, the magnetic behavior simulation was undertaken. Considering the variations of Earth's magnetic field in time and location, to obtain consistent and realistic results the average (within a year) of the different components of the magnetic field for the Waterloo, Ontario region (the location of the experiments) was adopted for the simulations (Table 2). Moreover, since the unitless relative magnetic permeability of low-carbon steels (ASTM 1020) range from 50 to 100 [57,58], a relative magnetic permeability of 75 was selected for this study. The duration of exposure to an external magnetic field will affect the magnetic behavior of ferromagnetic materials. In reality, ferromagnetic materials are affected by the magnetic field of the Earth from the beginning of their production process. There may also be some unknown external magnetic sources in the surrounding environment that affect the magnetic behavior of ferromagnetic objects [59]. However, as accurately as possible, we can apply the magnetic field of the Earth to the object and simulate its magnetic behavior, although some divergence will exist between the simulation and the experimental results.
To consider the Earth's magnetic field in the simulation, the rebar was located in a regular space ( Figure 2) with dimensions of 100 mm × 150 mm × 410 mm, which included the magnetic field presented in Table 2 and Figure 3. To have better control of simulation parameters, the box and rebar were meshed separately with tetrahedral meshes according to the specifications in Table 3 (Figure 4a,  b). Then, the rebar and box were jointly subjected to the simulation process as a single system ( Figure  4c). The values of the different components (X, Y, and Z) of the magnetic flux densities were recorded for the Y direction of the rebar (i.e., the path parallel to the rebar's length). This path is at the surface of the rebar, and extends from one side (Edge A) to the other side of the box (Edge B) ( Figure 5).            As observed in Figure 6, at first, the values of all of the components of magnetic flux densities are equal to the background magnetic flux (the magnetic field of the Earth). When the Y distance reaches about 18.065 mm, at the end of the rebar, the values of all of the components begin to reflect the impact of the magnetic properties of the ferromagnetic rebar on the magnetic fluxes. The values of all of the components have a harmonic variation because of the corrugated rebar shape. When the Y distance reaches the end of the rebar, all of the components of magnetic flux densities revert to the magnitudes of the background magnetic field. However, there is a distinguishable irregularity in the direction and values of all of the components at the location of Hole 2 (~301 mm from Edge A of the box). This irregularity is in the form of a minimum peak in the values of the Z and X magnetic flux densities, and in the form of a sudden change in the gradient of the Y-component of the magnetic flux density (a spike above the zero line, followed by a sudden dip below the zero line, then a sharp jump back to the zero line).
There are some outlier values in the different components of magnetic flux densities, which are related to the specifications of the elements used in this simulation. In order to have mesh element independent results, more accurate element specifications are implemented (Table 4). Then, the As observed in Figure 6, at first, the values of all of the components of magnetic flux densities are equal to the background magnetic flux (the magnetic field of the Earth). When the Y distance reaches about 18.065 mm, at the end of the rebar, the values of all of the components begin to reflect the impact of the magnetic properties of the ferromagnetic rebar on the magnetic fluxes. The values of all of the components have a harmonic variation because of the corrugated rebar shape. When the Y distance reaches the end of the rebar, all of the components of magnetic flux densities revert to the magnitudes of the background magnetic field. However, there is a distinguishable irregularity in the direction and values of all of the components at the location of Hole 2 (~301 mm from Edge A of the box). This irregularity is in the form of a minimum peak in the values of the Z and X magnetic flux densities, and in the form of a sudden change in the gradient of the Y-component of the magnetic flux density (a spike above the zero line, followed by a sudden dip below the zero line, then a sharp jump back to the zero line).
There are some outlier values in the different components of magnetic flux densities, which are related to the specifications of the elements used in this simulation. In order to have mesh element independent results, more accurate element specifications are implemented (Table 4). Then, the The values of all of the components have a harmonic variation because of the corrugated rebar shape. When the Y distance reaches the end of the rebar, all of the components of magnetic flux densities revert to the magnitudes of the background magnetic field. However, there is a distinguishable irregularity in the direction and values of all of the components at the location of Hole 2 (~301 mm from Edge A of the box). This irregularity is in the form of a minimum peak in the values of the Z and X magnetic flux densities, and in the form of a sudden change in the gradient of the Y-component of the magnetic flux density (a spike above the zero line, followed by a sudden dip below the zero line, then a sharp jump back to the zero line).
There are some outlier values in the different components of magnetic flux densities, which are related to the specifications of the elements used in this simulation. In order to have mesh element independent results, more accurate element specifications are implemented (Table 4). Then, the maximum and minimum values of the Z-component magnetic flux density (as a representative metric) from 295.0592 mm to 307.0592 mm (values symmetric about the fixed extent of Hole 2) are extracted. The location of Hole 2 was chosen because of its importance in our investigation. The differences between the maximum and minimum of these values are also used to verify the convergence of the simulation outcomes.
The values of the maximum and minimum magnetic flux densities become stable at mesh numbers 4 and 8 (Table 4), respectively (Figures 7 and 8). The difference between the maximum and minimum magnetic flux densities in the Z-component stabilize at rebar mesh #5 ( Figure 9). Hence, the result of mesh #8 is used for continuing the simulation. The magnetic flux density values for mesh #8 have no out-of-range or disorder trend, compared with the trend of rebar mesh #1, which was the initial simulation ( Figure 10). The location of Hole 2 was chosen because of its importance in our investigation. The differences between the maximum and minimum of these values are also used to verify the convergence of the simulation outcomes. The values of the maximum and minimum magnetic flux densities become stable at mesh numbers 4 and 8 (Table 4), respectively (figures 7 and 8). The difference between the maximum and minimum magnetic flux densities in the Z-component stabilize at rebar mesh #5 ( Figure 9). Hence, the result of mesh #8 is used for continuing the simulation. The magnetic flux density values for mesh #8 have no out-of-range or disorder trend, compared with the trend of rebar mesh #1, which was the initial simulation ( Figure 10).   The location of Hole 2 was chosen because of its importance in our investigation. The differences between the maximum and minimum of these values are also used to verify the convergence of the simulation outcomes. The values of the maximum and minimum magnetic flux densities become stable at mesh numbers 4 and 8 (Table 4), respectively (figures 7 and 8). The difference between the maximum and minimum magnetic flux densities in the Z-component stabilize at rebar mesh #5 (Figure 9). Hence, the result of mesh #8 is used for continuing the simulation. The magnetic flux density values for mesh #8 have no out-of-range or disorder trend, compared with the trend of rebar mesh #1, which was the initial simulation ( Figure 10).    To determine the effect of spacing, values of the magnetic flux density (rebar mesh #8) with different spacing were investigated. It was understood that increasing the spacing between the rebar and the recording point would result in some outliers in the trend of the Z-component magnetic flux density, related to the specifications of the tetrahedral elements used in the box. To make the results of the simulation independent of the element specifications, more accurate element specifications were applied to the box ( Table 5). As a representative result, the magnetic flux densities for the Zcomponent at a distance of 16 mm were extracted ( Figure 11). The maximum and minimum values from 295.0592 mm to 307.0592 mm (values related to Hole 2) were reviewed. Subsequently, the difference between the maximum and minimum values was investigated, and we found that the values became stable with box mesh #5 (figures [12][13][14].  To determine the effect of spacing, values of the magnetic flux density (rebar mesh #8) with different spacing were investigated. It was understood that increasing the spacing between the rebar and the recording point would result in some outliers in the trend of the Z-component magnetic flux density, related to the specifications of the tetrahedral elements used in the box. To make the results of the simulation independent of the element specifications, more accurate element specifications were applied to the box ( Table 5). As a representative result, the magnetic flux densities for the Zcomponent at a distance of 16 mm were extracted ( Figure 11). The maximum and minimum values from 295.0592 mm to 307.0592 mm (values related to Hole 2) were reviewed. Subsequently, the difference between the maximum and minimum values was investigated, and we found that the values became stable with box mesh #5 (figures [12][13][14]. To determine the effect of spacing, values of the magnetic flux density (rebar mesh #8) with different spacing were investigated. It was understood that increasing the spacing between the rebar and the recording point would result in some outliers in the trend of the Z-component magnetic flux density, related to the specifications of the tetrahedral elements used in the box. To make the results of the simulation independent of the element specifications, more accurate element specifications were applied to the box ( Table 5). As a representative result, the magnetic flux densities for the Z-component at a distance of 16 mm were extracted ( Figure 11). The maximum and minimum values from 295.0592 mm to 307.0592 mm (values related to Hole 2) were reviewed. Subsequently, the difference between the maximum and minimum values was investigated, and we found that the values became stable with box mesh #5 (Figures 12-14).     (Table 4).  (Table 4).     (Table 4).  (Table 4).      (Table 4).  (Table 4).  (Table 4).
According to figures 7-9 and 12-14, the outcomes from the simulation of the rebar with mesh #8 and box mesh #5, the optimum mesh specifications, were chosen for the rest of the investigations. After the optimum mesh specifications, achieved by increasing the meshing accuracy, the biggest difference in the maximum and minimum values of magnetic flux densities (of the same place) are respectively less than 0.01 µT and 0.04 µT, which are considered negligible. Additionally, the same values are observed for the maximum and minimum values of magnetic data when the element accuracy is increased. The same magnetic values mean that the results converge and are independent of the mesh specifications.
Carrying out simulations with optimum mesh specifications led to a graphical representation (Figure 15), which shows the behavior of the Z-component magnetic flux density at the location of Hole 2. Also, a planar slice of the magnetic field under the rebar in Figure 15 shows the conditions of the stray magnetic field around the rebar. As the distance from the rebar increases, the stray magnetic field around the rebar decreases relatively uniformly and symmetrically.  Figure 16 shows the values of magnetic flux densities of rebar with optimum mesh specifications at different spacings from the center of the rebar, ranging from 8 mm (surface level) to 72 mm (maximum distance from the surface). The behavior of the Z-component magnetic flux density is distinguishable at Hole 2 at a maximum of 16 mm from the rebar's center (Figure 16), which is a distance equal to 8 mm from the rebar's surface (this distance represents how thick the concrete above the bar can be for detection of the buried defects). According to the simulation results, it seems that the technique can be used only for very thin concrete layers with a maximum thickness of 8 mm. It should be mentioned that the simulations were performed under the Earth's present magnetic field,  (Table 4).
According to Figures 7-9 and 12-14, the outcomes from the simulation of the rebar with mesh #8 and box mesh #5, the optimum mesh specifications, were chosen for the rest of the investigations. After the optimum mesh specifications, achieved by increasing the meshing accuracy, the biggest difference in the maximum and minimum values of magnetic flux densities (of the same place) are respectively less than 0.01 µT and 0.04 µT, which are considered negligible. Additionally, the same values are observed for the maximum and minimum values of magnetic data when the element accuracy is increased. The same magnetic values mean that the results converge and are independent of the mesh specifications.
Carrying out simulations with optimum mesh specifications led to a graphical representation (Figure 15), which shows the behavior of the Z-component magnetic flux density at the location of Hole 2. Also, a planar slice of the magnetic field under the rebar in Figure 15 shows the conditions of the stray magnetic field around the rebar. As the distance from the rebar increases, the stray magnetic field around the rebar decreases relatively uniformly and symmetrically.  (Table 4).
According to figures 7-9 and 12-14, the outcomes from the simulation of the rebar with mesh #8 and box mesh #5, the optimum mesh specifications, were chosen for the rest of the investigations. After the optimum mesh specifications, achieved by increasing the meshing accuracy, the biggest difference in the maximum and minimum values of magnetic flux densities (of the same place) are respectively less than 0.01 µT and 0.04 µT, which are considered negligible. Additionally, the same values are observed for the maximum and minimum values of magnetic data when the element accuracy is increased. The same magnetic values mean that the results converge and are independent of the mesh specifications.
Carrying out simulations with optimum mesh specifications led to a graphical representation (Figure 15), which shows the behavior of the Z-component magnetic flux density at the location of Hole 2. Also, a planar slice of the magnetic field under the rebar in Figure 15 shows the conditions of the stray magnetic field around the rebar. As the distance from the rebar increases, the stray magnetic field around the rebar decreases relatively uniformly and symmetrically.  Figure 16 shows the values of magnetic flux densities of rebar with optimum mesh specifications at different spacings from the center of the rebar, ranging from 8 mm (surface level) to 72 mm (maximum distance from the surface). The behavior of the Z-component magnetic flux density is distinguishable at Hole 2 at a maximum of 16 mm from the rebar's center (Figure 16), which is a distance equal to 8 mm from the rebar's surface (this distance represents how thick the concrete above the bar can be for detection of the buried defects). According to the simulation results, it seems that the technique can be used only for very thin concrete layers with a maximum thickness of 8 mm. It should be mentioned that the simulations were performed under the Earth's present magnetic field,  Figure 16 shows the values of magnetic flux densities of rebar with optimum mesh specifications at different spacings from the center of the rebar, ranging from 8 mm (surface level) to 72 mm (maximum distance from the surface). The behavior of the Z-component magnetic flux density is distinguishable at Hole 2 at a maximum of 16 mm from the rebar's center (Figure 16), which is a distance equal to 8 mm from the rebar's surface (this distance represents how thick the concrete above the bar can be for detection of the buried defects). According to the simulation results, it seems that the technique can be used only for very thin concrete layers with a maximum thickness of 8 mm.
It should be mentioned that the simulations were performed under the Earth's present magnetic field, but ferromagnetic materials are considered saturated by the natural magnetic field, and may show stronger magnetic behavior.
For further investigation, the data-recording distance was increased to the maximum possible distance from the rebar, aligning with the inside edge of the box. At larger distances, the magnetic flux density trend becomes smoother and straighter, and approaches the background magnetic field.       We reviewed the different components of magnetic flux densities at the surface of the rebar, which were extracted from the optimum mesh specifications (Figure 18). The noise and out-of-range values at their minimum and results correlate well with the experimental results reported previously ( Figure 19) [37]. The laboratory flux magnetic density measurements were conducted using the PMI device that was specifically developed for this work in our lab, and is on its way to being commercialized. The PMI device works through scanning the SMFL arising from ferromagnetic structures [37].
values at their minimum and results correlate well with the experimental results reported previously ( Figure 19) [37]. The laboratory flux magnetic density measurements were conducted using the PMI device that was specifically developed for this work in our lab, and is on its way to being commercialized. The PMI device works through scanning the SMFL arising from ferromagnetic structures [37].
The patterns of laboratory and simulated outputs at the holes' locations follow the same trend; both curves generally have an up and down trend due to the corrugated shape of the rebar, along with a minimum value at the center of the Hole 2. As seen in Figure 18, the top hole that is ~301 mm from Edge A of the box (equal to ~282 mm from the rebar's start point (Figure 19)) is substantially easier to detect than the hole on the side of the rebar. Finding the irregularity in the magnetic data related to Hole 2 is easier because of the difference between the magnetic property of its surface (filled with air) and the rest of the rebar's surface (which has a different magnetic property). No detectable irregularity can be sensed in the surface magnetic flux densities for Hole 1, because more metal lies between it and the scanning line.   values at their minimum and results correlate well with the experimental results reported previously ( Figure 19) [37]. The laboratory flux magnetic density measurements were conducted using the PMI device that was specifically developed for this work in our lab, and is on its way to being commercialized. The PMI device works through scanning the SMFL arising from ferromagnetic structures [37]. The patterns of laboratory and simulated outputs at the holes' locations follow the same trend; both curves generally have an up and down trend due to the corrugated shape of the rebar, along with a minimum value at the center of the Hole 2. As seen in Figure 18, the top hole that is ~301 mm from Edge A of the box (equal to ~282 mm from the rebar's start point (Figure 19)) is substantially easier to detect than the hole on the side of the rebar. Finding the irregularity in the magnetic data related to Hole 2 is easier because of the difference between the magnetic property of its surface (filled with air) and the rest of the rebar's surface (which has a different magnetic property). No detectable irregularity can be sensed in the surface magnetic flux densities for Hole 1, because more metal lies between it and the scanning line.   The patterns of laboratory and simulated outputs at the holes' locations follow the same trend; both curves generally have an up and down trend due to the corrugated shape of the rebar, along with a minimum value at the center of the Hole 2. As seen in Figure 18, the top hole that is~301 mm from Edge A of the box (equal to~282 mm from the rebar's start point (Figure 19)) is substantially easier to detect than the hole on the side of the rebar. Finding the irregularity in the magnetic data related to Hole 2 is easier because of the difference between the magnetic property of its surface (filled with air) and the rest of the rebar's surface (which has a different magnetic property). No detectable irregularity can be sensed in the surface magnetic flux densities for Hole 1, because more metal lies between it and the scanning line.

Discussion
Assuming that the magnetic flux densities of different locations on the rebar are independent of one another, the probability graph method was used for fitting the magnetic flux values to a probability distribution. This method involves equating the empirical distribution of magnetic data (P i ) with the chosen cumulative distribution function (CDF). Next, the CDF function is written in linear form, which can be expressed as Equation (10) for Gamma CDF (GAMMADIS). The linearity is then used as a basis for determining whether the data can be modeled by a particular distribution. Additionally, the goodness-of-fit of the model is given by the coefficient of determination R 2 .
The magnetic flux density data were plotted against various probability distributions (normal, log-normal, Weibull, and gamma distributions); a gamma distribution was chosen based on the minimum least-squared error ( Figure 20). This distribution is based on a function of two parameters: α and β (Equation (11)); these were calculated by the mean and standard deviations (SD), which are 87.8 µT and 25.6 µT, respectively. As observed in Figure 21, the gamma function correlates well with the histogram frequency of data, and this approximation may be useful for estimation in practical cases.

Discussion
Assuming that the magnetic flux densities of different locations on the rebar are independent of one another, the probability graph method was used for fitting the magnetic flux values to a probability distribution. This method involves equating the empirical distribution of magnetic data ( ) with the chosen cumulative distribution function (CDF). Next, the CDF function is written in linear form, which can be expressed as Equation (10) for Gamma CDF (GAMMADIS). The linearity is then used as a basis for determining whether the data can be modeled by a particular distribution. Additionally, the goodness-of-fit of the model is given by the coefficient of determination 2 .
The magnetic flux density data were plotted against various probability distributions (normal, log-normal, Weibull, and gamma distributions); a gamma distribution was chosen based on the minimum least-squared error ( Figure 20). This distribution is based on a function of two parameters: and (Equation (11)); these were calculated by the mean and standard deviations (SD), which are 87.8 µT and 25.6 µT, respectively. As observed in Figure 21, the gamma function correlates well with the histogram frequency of data, and this approximation may be useful for estimation in practical cases. According to Figure 18, a Z-component magnetic flux density of less than 76 µT (without considering the edge effect and background magnetic field) corresponds to the location of Hole 2. Importing this value into the obtained CDF shows that 0.76% of the data is related to the defective locations. In other words, 0.76% of the rebar surface (at the scanned section) can be considered imperfect. This result can be verified by the Monte Carlo simulation method (based on inverse values of the obtained gamma distribution function). Figure 22 presents the probability of defects considering the mean, SD, and limit state, showing that the probability of defectiveness fluctuates until the first 300 trails are completed, and then stabilizes at the value of~0.75%. According to Figure 18, a Z-component magnetic flux density of less than 76 µT (without considering the edge effect and background magnetic field) corresponds to the location of Hole 2. Importing this value into the obtained CDF shows that 0.76% of the data is related to the defective locations. In other words, 0.76% of the rebar surface (at the scanned section) can be considered imperfect. This result can be verified by the Monte Carlo simulation method (based on inverse values of the obtained gamma distribution function). Figure 22 presents the probability of defects considering the mean, SD, and limit state, showing that the probability of defectiveness fluctuates until the first 300 trails are completed, and then stabilizes at the value of ~0.75%.   According to Figure 18, a Z-component magnetic flux density of less than 76 µT (without considering the edge effect and background magnetic field) corresponds to the location of Hole 2. Importing this value into the obtained CDF shows that 0.76% of the data is related to the defective locations. In other words, 0.76% of the rebar surface (at the scanned section) can be considered imperfect. This result can be verified by the Monte Carlo simulation method (based on inverse values of the obtained gamma distribution function). Figure 22 presents the probability of defects considering the mean, SD, and limit state, showing that the probability of defectiveness fluctuates until the first 300 trails are completed, and then stabilizes at the value of ~0.75%. For our statistical investigations, we considered the magnetic data as independent variables. Those independent variables were described by the chosen probability distribution with its particular distribution parameters, knowing that distribution allowed us to estimate an interval over which the  For our statistical investigations, we considered the magnetic data as independent variables. Those independent variables were described by the chosen probability distribution with its particular distribution parameters, knowing that distribution allowed us to estimate an interval over which the unknown future values may lie (with a stipulated level of confidence). Using the CDF of the gamma distribution, about 98% of all of the data are from 76 µT to 100 µT (Equation (12)). Hence, regarding the recorded magnetic data of the rebar, it can be predicted with 98% confidence that if the rebar was longer (by how much is irrelevant), the next values indicating a flawless rebar would be somewhere between 76-100 µT. Values outside this range should be reviewed as suspected defect locations.  Figure 23 indicates that the values of the X and Y-components of magnetic flux density still include some noise that is attributable to the mesh specifications. This issue can be investigated by increasing the mesh density in both the rebar and the box containing the rebar. Additionally, increasing the quality of clouds points defining the bar geometry can help achieve more accurate outcomes. For instance, Hole 1 in the solid part produced from the captured cloud points was not as deep as the real depth (measured directly), but the mesh was not corrected, as the intent was in part to test the laser scan cloud point data without additional data input management. unknown future values may lie (with a stipulated level of confidence). Using the CDF of the gamma distribution, about 98% of all of the data are from 76 µT to 100 µT (Equation (12)). Hence, regarding the recorded magnetic data of the rebar, it can be predicted with 98% confidence that if the rebar was longer (by how much is irrelevant), the next values indicating a flawless rebar would be somewhere between 76-100 µT. Values outside this range should be reviewed as suspected defect locations.
(100 µT) − (76 µT) = 0.99 − 0.01 = 0.98 (12) Figure 23 shows the values of different components of magnetic flux densities, 16 mm from the center of the rebar. For better observation of the irregularities related to Hole 2, all of the graphs are presented from 200 mm to 370 mm at the appropriate scale. The anomaly in the magnetic values at the location of Hole 2 can be detected in all of the components of the magnetic flux densities. Figure  23 indicates that the values of the X and Y-components of magnetic flux density still include some noise that is attributable to the mesh specifications. This issue can be investigated by increasing the mesh density in both the rebar and the box containing the rebar. Additionally, increasing the quality of clouds points defining the bar geometry can help achieve more accurate outcomes. For instance, Hole 1 in the solid part produced from the captured cloud points was not as deep as the real depth (measured directly), but the mesh was not corrected, as the intent was in part to test the laser scan cloud point data without additional data input management.

Conclusions
Being able to detect defects in steel infrastructure would substantially improve risk management and condition evaluation over time. To this end, mathematical simulations were carried out on a pre-flawed specimen that was laser-scanned to generate a point cloud surface map. This map was then used as a basis to develop a model. The intent was to establish detectability limits for very small flaws in order to reduce Type I and Type II errors in anomaly detection.
The magnetic behavior of the ferromagnetic rebar specimen was simulated with a finite element-based software considering the background magnetic field. Different components of magnetic flux densities on the surface showed consistent harmonic trends because of the corrugated shape of the rebar. However, there were specific irregularities in the direction and values for the different components of magnetic flux densities at the location of Hole 2. Simulated patterns could be correlated with the experimental data at the holes' locations, so the top hole (Hole 2) was easily located, but Hole 1 was not, because of its orientation in the magnetic field and because the point cloud model did not replicate its true depth. The gamma probability distribution was chosen to statistically assess the magnetic flux density behavior of the rebar. Two main outcomes were extracted: 0.76% of the scanned section of the rebar was considered defective, and if the rebar specimen were longer, the Z-component magnetic flux density values indicating flawless rebar would be predicted to lie between 76-100 µT with 98% confidence.
The values of the different components of magnetic flux densities at different distances from the rebar were reviewed. Increasing the vertical distance of the data recording line led to a logarithmic reduction of magnetic flux density values. As this distance was increased, the magnetic flux density values became approximately constant and close to the background magnetic field. In conclusion:

•
The pattern of the simulation results at defect locations were similar to the outputs of previous physical experiments; • The background magnetic field had a significant effect on the trend and values of different components of the magnetic flux density; • All of the magnetic flux density components displayed correctly located anomalies corresponding to the defect on the top surface of the rebar; • Increasing the distance from the rebar changed the trend and values of the magnetic flux densities such that at some distance, the anomaly became undetectable; • To detect various shapes and sizes of defects at different places along a rebar specimen, additional magnetic parameters should be considered. For instance, the Z-component of the magnetic flux density was totally constant on the sides of the rebar, and could not detect the anomaly arising from Hole 1; • The stray magnetic field around the rebar decreased relatively symmetrically by increasing the distance from the rebar; and Funding: This research received no direct funding.