Skip to main content

TECHNOLOGY AND CODE article

Front. Bioeng. Biotechnol., 13 May 2022
Sec. Biomaterials

A Computational Pitting Corrosion Model of Magnesium Alloys

Chia-Jung ChangChia-Jung Chang1Chih-Han ChangChih-Han Chang1Tin-Kan Hung
Tin-Kan Hung2*
  • 1Department of Biomedical Engineering, National Cheng Kung University, Tainan, Taiwan
  • 2Department of Bioengineering, University of Pittsburgh, Pittsburgh, PA, United States

Controlling the corrosion rate of implants to maintain mechanical properties during tissue healing is significant in developing magnesium alloy implants. In addition to surface treatment and material properties, the study of geometric alteration and mechanical strength are also vital for implant development. In this study, we developed a three-dimensional model for semi-autonomous computational pitting corrosion. It is based on the Monte Carlo method, modeling magnesium alloy implants toward clinical application. The corrosion probability is based on the number of exposed surfaces to saline and the oxidation characteristics of the elements. The computational results are well compared with the experimental measurement using micro-computed tomography (micro-CT) in 500 h. Subsequently, the computational analysis is extended to 3,000 h of corrosion analysis. The 3D model appears promising to assist the development of biodegradable implants.

Introduction

Magnesium (Mg) alloys are prospective biodegradable materials for medical applications in recent years. The alloys have low density, biocompatibility, and mechanical and biodegradable properties comparable to human bones (Soya et al., 2014; Živić et al., 2014; Bowen et al., 2018). Magnesium bone implants can minimize the effects of stress shielding in traditional metallic implants (Xu et al., 2007). Biodegradable magnesium alloy stents have a better bracing force and tissue compatibility than absorbable polymer stents (Li et al., 2019). However, there remains difficulty to control the corrosion rate of implants. Magnesium alloys corrode too fast for maintaining mechanical properties, and the implant function before tissue regeneration and healing are well developed (Song and Song, 2007; Kirkland et al., 2012; Kirkland and Birbilis, 2014).

Bioengineering strategies for controlling the corrosion rate of Mg alloys were discussed in many studies on surface modification (Li et al., 2013; Córdoba et al., 2016), element alloying (Hort et al., 2010; Jin et al., 2015), and microstructure modification (Wu et al., 2017). Improvements in corrosion resistance and implant surfaces remain the main tasks in research and development. The spatial characteristics are also important factors affecting the mechanical behaviors of the implant. Pitting corrosion could cause stress concentration and crack formation (Suresh, 1998; Gu et al., 2010), leading toward peeling off, fragmented fracture (Li et al., 2016), and the collapse of implants.

With advances in technology and computing power, computational models are becoming practical for analyzing complex corrosion phenomena. Quite a number of studies were focused on modeling pitting corrosion and were well compared in review (Jafarzadeh et al., 2019). Non-autonomous models employed the classical transport equation for numerical solutions to transport kinetics in the electrolyte (Sharland et al., 1989; Krouse et al., 2014). They provided a good framework, but tracking pit growth is addressed as a separate part that increases the modeling complexity (Jafarzadeh et al., 2019). Autonomous models such as cellular automata (CA) techniques (Cui et al., 2019), the finite volume method (FVM) (Scheiner and Hellmich, 2009), peridynamic (PD) formulations (Jafarzadeh et al., 2018), and phase-field (PF) models (Tsuyuki et al., 2018) were focused on calculating the motion of corrosion front directly. Nonetheless, there remained many limitations with various autonomous models (Jafarzadeh et al., 2019). The CA model had difficulty to calibrate, and its grid size and time are not physical parameters (Stafiej et al., 2013; Pérez-Brokate et al., 2016), while the FVM did not capture mechanical damage and diffuse the corrosion layer (Jafarzadeh et al., 2019); the PD model had nonlocal effects at boundaries (Le and Bobaru, 2018), and the damage evolution in PF models appeared limited because of its dependency in selecting the energy functional (Geelen et al., 2019) with the dense mathematical framework.

In this study, a semi-autonomous computational model based on the Monte Carlo process is developed for pitting corrosion, with elements having different material characteristics and the surrounding condition. Effort is also made on the oxidized surface characteristics for the specimen with pitting corrosion. The computational results are compared with the experimental 3D images reconstructed by micro-computed tomography (micro-CT, μCT) for 500 h experiments, making the combined use of computation and experiment feasible to project long-term corrosion development.

Materials and Methods

Corrosion Experiment

Corrosion of Mg alloy, AZ61, in 0.90% saline was studied individually using 11 specimens of size 15 mm × 15 mm × 2 mm. Weight compositions (wt%) of the AZ61 alloy are listed in Table 1. With initial 570 mm2 surface area (SO), each specimen was immersed and suspended in 114 ml of solution at room temperature (25°C ± 2°C). Based on the ASTM31-12a protocol (NACE, 2012), it was dried every 25 h for micro-CT scanning and mass measurement. The saline was then refreshed for the next 25 h of corrosion measurement. Including the initial scan, 21 scans were obtained for each experiment over a period of 500 h. They were obtained using Bruker Micro-CT Sky Scan 1,076 with the following parameters: source voltage = 100 kV, source current = 100 μA, filter = Al 1.0 mm, exposure time = 316 ms, image pixel size = 34.44 μm, and rotation step = 1.000°. Sky Scan software and NRecon software were employed for volumetric reconstruction of each scan measurement, resulting in a total of 22 × 21 three-dimensional (3D) imaging.

TABLE 1
www.frontiersin.org

TABLE 1. Compositions (wt%) of AZ61 in this study.

The micro-CT scans the specimen mass at each scanning stage, which was calculated from the volume multiplied by its initial density of 1.7 g/cm3. The mass loss rate .ML per total surface (SO) per hour (mg/mm2 h.) is calculated by the following equation:

.ML=Mo– MNSo T (1)

where MO (about 760 mg) is the initial mass measured by a microbalance, MN (mg) is the calculated mass of the specimen at the Nth scan via the measured volume from micro-CT, and T is the submerged time in saline.

Corrosion Simulation

Modeling of the Mg alloy specimen begins with setting uniform hexahedral elements. The corrosion probability (CPi) for each micro element is assessed at each time step (Δt). Its value for the ith element is based on two factors: the expose attribute (EAi) and the oxide attribute (OAi). The EAi value indicates how many surfaces of an element are exposed to saline (Kelly et al., 2002). A hexahedral element has six surfaces; thus, six chances in contact with saline solution are possible. Initially, all the inner elements are not in contact with saline; EAi is zero in the computational algorithm. For elements on the top and bottom surfaces and four sides, EAi is set to one; on the edge, EAi = 2 and EAi = 3 for the corner element. Song and Atrens (1999) reported that the inhomogeneous property on the oxide layer and different corrosion potential between the oxide layer and the Mg alloy matrix inside can cause pitting corrosion. To simulate the phenomena, the corrosion speed of an element is indicated by its OAi value. The code of simulation model is written in Matlab, and the output data are stored as a txt file or an .inp file, which is suitable for several numerical analysis software.

Corrosion rates of the Mg oxide surface are much lower than those of the inner alloy; they are needed for assessing OAi values. Because the corrosion rate of the interior AZ61 Mg alloy element in saline is 13.7 times faster than that of the external oxide layer elements estimated from the experimental data of Wu et al. (2015) and Tkacz et al. (2016), the OAi value of the oxide layer is assigned to 0.07 when OAi = 1 is set for micro elements of intact magnesium. However, OAi values for the spotted microporous surface are not known a priori for a specimen. To identify the locations of initial porous elements on the oxide surface, the pit distribution of a specimen measured at the first 200 h is used for specifying OAi values on the surface. The final selection of OAi values between 0.4 and 0.6 was further assessed by trial and error via the first 200 h of simulation. This step leads to a proper comparison between simulation and experiment. The combined effect of OAi and EAi is expressed as the corrosion attribute (CAi):

CAi=EAiOAi(2)

The corrosion probability (CPi) of the ith element is determined as follows:

CPi=MLMeCAii=1nCAi(3)

where the ratio of mass loss (ML=Me˙SΔt) to element mass (Me) indicates the number of elements being corroded at the time step. The ratio of CAi to the summation of CAi values of all the n-elements stands for the affecting corrosion rate of each element. Notice that for intact magnesium micro elements, CPi remains nil as OAi = 0. They will not be corroded until they are in contact with saline through pitting surfaces. The increase in the pitting surface with saline could increase the corrosion speed with time.

A Monte Carlo (Sharland et al., 1989) algorithm is employed to randomly select ML/Me (number of elements) elements to be corroded and those elements would become void in the computational model at the next time step. The algorithm involves generating a random number for each element to compare CPi values to determine if this element is corroded and removed. The procedure can be repeated N time steps to relieve elements (N is the number of cycles of corrosion, which can be defined according to user needs, each cycle is default to 25 h). The results of the first 500 h are obtained and compared with the experimental results, and the total simulation time was extended to 3,000 h in this study.

Image Processing and Quantization

To identify the corrosion characteristics, image processing procedures are employed using Matlab to assess the pitting characteristics. Both the micro-CT images and simulation results are converted to 3D models in STL files. To obtain the pitting maximum depth, the maximum depth of each corrosion pit is calculated by the difference between the maximum and minimum of z as shown in Supplementary Figure S1 in the Appendix. To obtain the number of corrosion pits and the projected corrosion area on the surface, the cross sections from the 3D model are exported for morphological feature recognition. After using the Canny method (Canny, 1986) to detect the edge of the pit pattern, the enclosed region from the pits is filled for distinguishing the corroded and non-corroded elements as shown in Supplementary Figure S2.

Results

The experimental results of mass loss with the corrosion time for AZ61 specimen in saline are presented in Figure 1A. Each I-bar represents micro-CT measurements of 11 specimens with a mean value indicated by a circular point. The mass loss per surface area (2 mm × 15 mm × 15 mm) can be fitted by MassLossPerArea(MLPA)=0.003T in 500 h. The squares in the figure are the results of simulation; they are in good agreement with the mean values of experimental measurements. 3D graphics in Figure 1B are via Autodesk Fusion 360 for making magnesium alloy rendering effects.

FIGURE 1
www.frontiersin.org

FIGURE 1. Comparison of overall AZ61 corrosion in simulation and experiment. (A) Mass loss per area of experimental specimens and simulation model at different corrosion times. (B) 3D graphics.

Computational results are obtained after identifying the effects of mesh sizes and time steps. The MLPA calculated from simulation is based on four sizes (0.250, 0.175, 0.135, and 0.125 mm) of hexahedral elements. As shown in Supplementary Figure S5, the difference between element sizes of 0.135 and 0.125 mm is under 1.5%, and the micro element of 0.125 mm is used for this study.

Figure 2 compares the corrosion simulation with experimental results of the mid-section (y = 7.5 mm) of a specimen for corrosion time at 50, 250, and 500 h. The experimental data are also assessed by time variations of the maximum corrosion depth. As shown in Figure 2C, the increase in the mean value with time (Dm indicated by a circular point) for 11 specimens can be fitted by an exponential function as follows:

Dm=1.3(1-e0.0044T)(4)

FIGURE 2
www.frontiersin.org

FIGURE 2. Comparison of experimental and simulated perpendicular cross sections. (A,B) Comparison of experimental and simulated corrosion patterns across three parallel cross sections at y = 7.5 mm at T = 50, 250, and 500 h. (C) Maximum depth of specimen at different corrosion times.

They are in general agreement with computational simulation depicted by square points in Figure 2C.

Comparison of corrosion patterns between simulation and experimental results is presented in Figure 3A for these three planes Z = 0, 0.125, and 0.25 mm at three times T = 100, 300, and 500 h, respectively. The results are further analyzed by measuring the total corroded surface AC on these three planes over 500 h. Figure 3B presents the mass loss per surface area (MLPA) obtained from the computational analysis and experimental measurement. Three features can be noticed in this figure:

1) The slope of AC/AO on the surface is practically the same for experimental and computational data.

2) The slope of corrosion for the inner layer (Z = 0.125 mm) is much smaller than that on the surface for T < 250 h; they are higher than those on the mid-section (Z = 0.25 mm).

3) The slopes for these two inner layers continue to increase and approach those of the surface as T reaches 500 h.

FIGURE 3
www.frontiersin.org

FIGURE 3. Comparison of experimental and simulated horizontal cross sections. (A) Comparison of experimental and simulated corrosion patterns across three parallel cross sections at z = 0, 0.125, and 0.25 mm at T = 100, 300, and 500 h, respectively. (B) Projected corrosion area AC on three planes at z = 0, 0.125, and 0.25 mm.

Figure 4A presents all the results of MLPA obtained by the computational model. The corrosion rate becomes higher when T approaches 1,000 h. The increase in mass loss exceeds the straight line, which is extended from the 500 h experimental measurement shown in Figure 1. The salient results in Figure 4B indicate that the corrosion in the inner layer (Z = 0.125 mm) exceeds that on the surface when T > 500 h. The phenomena appear due to the oxide surface effect on pitting corrosion. In other words, the inner layer is more prone to pitting corrosion.

FIGURE 4
www.frontiersin.org

FIGURE 4. Long-term predictions. (A) Mass loss per area of the simulation model compares to the initial experimental slope at different corrosion times in the long-term. (B) 3D model and corrosion patterns across two parallel cross-sections at z = 0, and 0.125 mm at T = 500, 1,000, 1,500, and 2,000 h, respectively.

Discussion

Developing degradable implants requires controlling the corrosion rate with changes in geometry, environment, and mechanical characteristics for clinical application. The 3D corrosion simulation model appears promising to complement the experimental development of implants. It could facilitate geometric changes of the model and couple with 3D modeling of tissue growth in implant voids. The altered geometry is needed for stress and strain analysis. All the biomechanical parameters could affect corrosion patterns. Also, electrochemical factors might interact with corrosion in remodeling processes. Modeling an implant weakness can be made by simulation for improvements before or during the experimental study. The effect on the oxide layer can be changed by coatings. The challenge is to determine if an implant can maintain its functional capability in the human body.

Although the non-autonomous models (mostly based on FEM analysis) are in a good framework and commercialized, they still have difficulties in boundary tracking and domain updating (Jafarzadeh et al., 2019), and several autonomous models still have many limitations for developing implants. In our semi-autonomous computational model, the effect of oxide layers on corrosion is considered and modeled by the oxide attribute (OAi). The computational results are in agreement with the micro-CT measurement for 20 time steps in 500 h. The computation appears effective to pair with experimental measurements, especially to project long-term corrosion patterns. Also, the 3D asymmetrical corroded morphology is practical in comparison with 2D models (Sharland, 1987; Pérez-Brokate et al., 2016).

In the overall mass loss evaluation as shown in Figure 1A, the data at 200–300 h have a slight tendency to reduce the degradation rate as the past report shows that corrosion products of magnesium decreased the corrosion rate (Lindström et al., 2002; Liu et al., 2019). However, due to the measurement and recombination errors of micro-CT, there are some fluctuations in the experimental data. In the study, we employed simple linear fitting for small data, and it seems to have a minor effect on the comparison of simulation results with experiments in the different indexes in our study. The increase in the corrosion rate in inner layers shown in Figure 3A indicates the surface effect.

Evaluation of corrosion pits is good for assessing function in medical applications (Apostolopoulos et al., 2013). In this study, the positions of the porous spots for simulation are based on experimental results after the specimen has been submerged in saline for 200 h. The fragile and loose places on the oxide layer could be identified by a more advanced method. A case-specific simulation can be accomplished when more accurate porous spots are prescribed as the initial oxide surface for the first 200 h simulation to specify OAi values.

The area and the type of corrosion pit patterns are important indicators in the evaluation of medical supplies. In terms of the consistency between the simulation and experiment, the geometric features of the corrosion could have been highly related to the oxide layer. The property and structure of the oxide layer may even be the dominant factor in corrosion patterns.

Exploring the function of a medical material not only involves assessing the overall amount of corrosion but also includes a discussion of changes in geometric, local conditions and mechanical properties over time. The advantage of this simulation method is that it established a corrosion simulation for the 3D model with pit features that have the potential to be used in combination with the finite element analysis (FEA) with remodeling processes. The FE method can be used to obtain the local stress or other characteristics of the material. These stress values can be applied to the corrosion model to determine the probability of corrosion. The geometric changes due to corrosion are useful for re-evaluating the stress distribution by the FE method. FEA software is good for analyzing the mechanical (Bian et al., 2021) and electrochemical factors. These data can be implemented in a corrosion model with tissue in the remodeling process. The weaknesses of medical materials and implants such as bone implants, dental implants, and stents can be studied for optimization and improvement. Also, the oxide layer on this model can be changed to a coating layer to simulate a case of improved medical devices. This model also has the potential to simulate different composite materials such as biodegradable polymer implant made in PGA, PLA, and PLGA, and its feasibility needs further verification.

The ultimate goal is to simulate conditions that mimic the biodegradable medical material in the human body in order to better assess whether medical materials and implants in the body can still maintain their functional capability during the process of degradation to achieve the goal of developing degradable materials and implants.

Conclusion

A simulation model for the pitting corrosion of the Mg alloy is established in this study. Pitting development is modeled by specifying the non-homogeneous oxide layer based on experimental data of the first 200 h corrosion of the specimen. Pitting corrosion is quantified by measuring the surface area developed across three sections of 11 specimens along with the maximum pitting depths. The results are well compared with those obtained by computational simulation for 500 h of the corrosion pattern. The computational mass loss is further studied from the results at 500–3,000 h. Combining a novel 3D simulation model with the experimental measurement using micro-CT for implant development is demonstrated. Proper setting of characteristics for the oxide layer is vital to computational modeling of pitting corrosion. The 3D model is useful for predicting long-term corrosion patterns and implant function. It could complement experimental efforts toward developing biomedical degradable implants.

Data Availability Statement

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

Author Contributions

C-JC conceived and proposed the study followed by developing the simulation code, conducting experimental measurement, and preparing the manuscript. T-KH and C-HC designed and supervised the project and worked with Chia-Jung for this report.

Funding

This study was supported by the Department of Biomedical Engineering, the Medical Device Innovation Center, and the Instrument Development Center of the National Cheng Kung University. Grant support: Ministry of Science and Technology (MOST) (#110-2221-E-006 -012 -MY2).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Acknowledgments

C-JC was awarded a visiting scholarship from the Ministry of Education, Taiwan, to visit T-KH. The academic guidance and discussion are enriched via the NCKU-Pitt Initiative at the University of Pittsburgh, Pennsylvania.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fbioe.2022.887444/full#supplementary-material

Supplementary Figure S1 | STL files of magnesium alloy specimens.

Supplementary Figure S2 | Image processing procedures for quantitating the pitting characteristics.

Supplementary Figure S3 | Flowchart of simulation.

Supplementary Figure S4 | Experimental data of AZ61 in saline. (A) Number of pits. (B) Total projected area of pits on the surface.

Supplementary Figure S5 | Comparison of results for four element sizes. (A) Convergence analysis by different grid sizes. (B) Reproducibility analysis by the same grid size.

Supplementary Figure S6 | Ratio of the surface area in the experiment.

References

Apostolopoulos, C. A., Demis, S., and Papadakis, V. G. (2013). Chloride-induced Corrosion of Steel Reinforcement - Mechanical Performance and Pit Depth Analysis. Construction Building Mater. 38, 139–146. doi:10.1016/j.conbuildmat.2012.07.087

CrossRef Full Text | Google Scholar

Bian, D., Zhou, X., Liu, J., Li, W., Shen, D., Zheng, Y., et al. (2021). Degradation Behaviors and In-Vivo Biocompatibility of a Rare Earth- and Aluminum-free Magnesium-Based Stent. Acta Biomater. 124, 382–397. doi:10.1016/j.actbio.2021.01.031

PubMed Abstract | CrossRef Full Text | Google Scholar

Bowen, P. K., Seitz, J. M., Guillory, R. J., Braykovich, J. P., Zhao, S., Goldman, J., et al. (2018). Evaluation of Wrought Zn-Al Alloys (1, 3, and 5 Wt % Al) through Mechanical and In Vivo Testing for Stent Applications. J. Biomed. Mater. Res. B Appl. Biomater. 106 (1), 245–258. doi:10.1002/jbm.b.33850

PubMed Abstract | CrossRef Full Text | Google Scholar

Canny, J. (1986). A Computational Approach to Edge Detection. IEEE Trans. Pattern Anal. Mach. Intell. PAMI-8 (6), 679–698. doi:10.1109/tpami.1986.4767851

CrossRef Full Text | Google Scholar

Córdoba, L., Montemor, M., and Coradin, T. (2016). Silane/TiO 2 Coating to Control the Corrosion Rate of Magnesium Alloys in Simulated Body Fluid. Corrosion Sci. 104, 152–161. doi:10.1016/j.corsci.2015.12.006

CrossRef Full Text | Google Scholar

Cui, C., Ma, R., Chen, A., Pan, Z., and Tian, H. (2019). Experimental Study and 3D Cellular Automata Simulation of Corrosion Pits on Q345 Steel Surface under Salt-spray Environment. Corrosion Sci. 154, 80–89. doi:10.1016/j.corsci.2019.03.011

CrossRef Full Text | Google Scholar

Geelen, R. J. M., Liu, Y., Hu, T., Tupek, M. R., and Dolbow, J. E. (2019). A Phase-Field Formulation for Dynamic Cohesive Fracture. Comput. Methods Appl. Mech. Eng. 348, 680–711. doi:10.1016/j.cma.2019.01.026

CrossRef Full Text | Google Scholar

Gu, X. N., Zhou, W. R., Zheng, Y. F., Cheng, Y., Wei, S. C., Zhong, S. P., et al. (2010). Corrosion Fatigue Behaviors of Two Biomedical Mg Alloys - AZ91D and WE43 - in Simulated Body Fluid. Acta Biomater. 6 (12), 4605–4613. doi:10.1016/j.actbio.2010.07.026

PubMed Abstract | CrossRef Full Text | Google Scholar

Hort, N., Huang, Y., Fechner, D., Störmer, M., Blawert, C., Witte, F., et al. (2010). Magnesium Alloys as Implant Materials - Principles of Property Design for Mg-RE Alloys☆. Acta Biomater. 6 (5), 1714–1725. doi:10.1016/j.actbio.2009.09.010

PubMed Abstract | CrossRef Full Text | Google Scholar

Jafarzadeh, S., Chen, Z., and Bobaru, F. (2019). Computational Modeling of Pitting Corrosion. Corrosion Rev. 37 (5), 419–439. doi:10.1515/corrrev-2019-0049

CrossRef Full Text | Google Scholar

Jafarzadeh, S., Chen, Z., and Bobaru, F. (2018). Peridynamic Modeling of Repassivation in Pitting Corrosion of Stainless Steel. Corrosion 74 (4), 393–414. doi:10.5006/2615

CrossRef Full Text | Google Scholar

Jin, W., Wu, G., Feng, H., Wang, W., Zhang, X., and Chu, P. K. (2015). Improvement of Corrosion Resistance and Biocompatibility of Rare-Earth WE43 Magnesium alloy by Neodymium Self-Ion Implantation. Corrosion Sci. 94, 142–155. doi:10.1016/j.corsci.2015.01.049

CrossRef Full Text | Google Scholar

Kelly, R. G., Scully, J. R., Shoesmith, D. W., and Buchheit, R. G. (2002). Electrochemical Techniques in Corrosion Science and Engineering. CRC Press.

Google Scholar

Kirkland, N. T., and Birbilis, N. (2014). Magnesium Biomaterials: Design, Testing, and Best Practice. Berlin, Germany: Springer.

Google Scholar

Kirkland, N. T., Birbilis, N., and Staiger, M. P. (2012). Assessing the Corrosion of Biodegradable Magnesium Implants: a Critical Review of Current Methodologies and Their Limitations. Acta Biomater. 8 (3), 925–936. doi:10.1016/j.actbio.2011.11.014

PubMed Abstract | CrossRef Full Text | Google Scholar

Krouse, D., Laycock, N., and Padovani, C. (2014). Modelling Pitting Corrosion of Stainless Steel in Atmospheric Exposures to Chloride Containing Environments. Corrosion Eng. Sci. Technol. 49 (6), 521–528. doi:10.1179/1743278214y.0000000221

CrossRef Full Text | Google Scholar

Le, Q. V., and Bobaru, F. (2018). Surface Corrections for Peridynamic Models in Elasticity and Fracture. Comput. Mech. 61 (4), 499–518. doi:10.1007/s00466-017-1469-1

CrossRef Full Text | Google Scholar

Li, H., Liu, Y., Pang, S., Liaw, P. K., and Zhang, T. (2016). Corrosion Fatigue Behavior of a Mg-Based Bulk Metallic Glass in a Simulated Physiological Environment. Intermetallics 73, 31–39. doi:10.1016/j.intermet.2016.03.003

CrossRef Full Text | Google Scholar

Li, M., Cheng, Y., Zheng, Y. F., Zhang, X., Xi, T. F., and Wei, S. C. (2013). Plasma Enhanced Chemical Vapor Deposited Silicon Coatings on Mg alloy for Biomedical Application. Surf. Coat. Technol. 228, S262–S265. doi:10.1016/j.surfcoat.2012.05.128

CrossRef Full Text | Google Scholar

Li, Y., Wang, L., Chen, S., Yu, D., Sun, W., and Xin, S. (2019). Biodegradable Magnesium alloy Stents as a Treatment for Vein Graft Restenosis. Yonsei Med. J. 60 (5), 429. doi:10.3349/ymj.2019.60.5.429

PubMed Abstract | CrossRef Full Text | Google Scholar

Lindström, R., Svensson, J.-E., and Johansson, L.-G. (2002). The Influence of Carbon Dioxide on the Atmospheric Corrosion of Some Magnesium Alloys in the Presence of NaCl. J. Electrochem. Soc. 149 (4), B103.

Google Scholar

Liu, H., Cao, F., Song, G.-L., Zheng, D., Shi, Z., Dargusch, M. S., et al. (2019). Review of the Atmospheric Corrosion of Magnesium Alloys. J. Mater. Sci. Technol. 35 (9), 2003–2016. doi:10.1016/j.jmst.2019.05.001

CrossRef Full Text | Google Scholar

Nace, A. (2012). ASTMG31–12a Standard Guide for Laboratory Immersion Corrosion Testing of Metals. West Conshohocken, PA: ASTM International.

Google Scholar

Pérez-Brokate, C. F., Di Caprio, D., Feron, D., De Lamare, J., and Chausse, A. (2016). Three Dimensional Discrete Stochastic Model of Occluded Corrosion Cell. Corrosion Sci. 111, 230–241. doi:10.1016/j.corsci.2016.04.009

CrossRef Full Text | Google Scholar

Scheiner, S., and Hellmich, C. (2009). Finite Volume Model for Diffusion-And Activation-Controlled Pitting Corrosion of Stainless Steel. Comput. Methods Appl. Mech. Eng. 198 (37-40), 2898–2910. doi:10.1016/j.cma.2009.04.012

CrossRef Full Text | Google Scholar

Sharland, S. M. (1987). A Review of the Theoretical Modelling of Crevice and Pitting Corrosion. Corrosion Sci. 27 (3), 289–323. doi:10.1016/0010-938x(87)90024-2

CrossRef Full Text | Google Scholar

Sharland, S. M., Jackson, C. P., and Diver, A. J. (1989). A Finite-Element Model of the Propagation of Corrosion Crevices and Pits. Corrosion Sci. 29 (9), 1149–1166. doi:10.1016/0010-938x(89)90051-6

CrossRef Full Text | Google Scholar

Song, G. L., and Atrens, A. (1999). Corrosion Mechanisms of Magnesium Alloys. Adv. Eng. Mater. 1 (1), 11–33. doi:10.1002/(sici)1527-2648(199909)1:1<11::aid-adem11>3.0.co;2-n

CrossRef Full Text | Google Scholar

Song, G., and Song, S. (2007). A Possible Biodegradable Magnesium Implant Material. Adv. Eng. Mater. 9 (4), 298–302. doi:10.1002/adem.200600252

CrossRef Full Text | Google Scholar

Soya, Y., Yoshihara, S., Ohmura, Y., Donald, B. J. M., and Galvin, E. (2014). “Corrosion Behavior of Engineering Materials in Flow Field,” in Advanced Materials Research (Trans Tech Publ). doi:10.4028/www.scientific.net/amr.922.722

CrossRef Full Text | Google Scholar

Stafiej, J., Di Caprio, D., and Bartosik, Ł. (2013). Corrosion-passivation Processes in a Cellular Automata Based Simulation Study. J. Supercomput 65 (2), 697–709. doi:10.1007/s11227-013-0933-8

CrossRef Full Text | Google Scholar

Suresh, S. (1998). Fatigue of Materials. Cambridge, United Kingdom: Cambridge University Press.

Google Scholar

Tkacz, J., Slouková, K., Minda, J., Drábiková, J., Fintová, S., Doležal, P., et al. (2016). Corrosion Behavior of Wrought Magnesium Alloys AZ31 and AZ61 in Hank's Solution. Koroze a Ochrana Materiál 60 (4), 101–106. doi:10.1515/kom-2016-0016

CrossRef Full Text | Google Scholar

Tsuyuki, C., Yamanaka, A., and Ogimoto, Y. (2018). Phase-field Modeling for pH-dependent General and Pitting Corrosion of Iron. Sci. Rep. 8 (1), 12777–12814. doi:10.1038/s41598-018-31145-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, F., Liang, J., and Li, W. (2015). Electrochemical Deposition of Mg(OH)2/GO Composite Films for Corrosion protection of Magnesium Alloys. J. Magnesium Alloys 3 (3), 231–236. doi:10.1016/j.jma.2015.08.004

CrossRef Full Text | Google Scholar

Wu, T.-C., Ho, Y.-H., Joshi, S. S., Rajamure, R. S., and Dahotre, N. B. (2017). Microstructure and Corrosion Behavior of Laser Surface-Treated AZ31B Mg Bio-Implant Material. Lasers Med. Sci. 32 (4), 797–803. doi:10.1007/s10103-017-2174-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu, L., Yu, G., Zhang, E., Pan, F., and Yang, K. (2007). In Vivo corrosion Behavior of Mg-Mn-Zn alloy for Bone Implant Application. J. Biomed. Mater. Res. 83A (3), 703–711. doi:10.1002/jbm.a.31273

PubMed Abstract | CrossRef Full Text | Google Scholar

Živić, F., Grujović, N., Manivasagam, G., Richard, C., Landoulsi, J., and Petrović, V. (2014). The Potential of Magnesium Alloys as Bioabsorbable/Biodegradable Implants for Biomedical Applications. Tribology in Industry 36 (1), 67–73.

Google Scholar

Keywords: computational simulation, magnesium alloy, micro-CT, Monte Carlo method, pitting corrosion

Citation: Chang C-J, Chang C-H and Hung T-K (2022) A Computational Pitting Corrosion Model of Magnesium Alloys. Front. Bioeng. Biotechnol. 10:887444. doi: 10.3389/fbioe.2022.887444

Received: 01 March 2022; Accepted: 22 March 2022;
Published: 13 May 2022.

Edited by:

Yin Fang, Nanyang Technological University, Singapore

Reviewed by:

Junyu Wang, Brown University, United States
Chunguang Yang, Institute of Metal Research (CAS), China

Copyright © 2022 Chang, Chang and Hung. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Tin-Kan Hung, dGtodW5nQHBpdHQuZWR1

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.