# Monte Carlo Simulations of Coupled Diffusion and Surface Reactions during the Aqueous Corrosion of Borosilicate Glasses

## Abstract

Borosilicate nuclear waste glasses develop complex altered layers as a result of coupled processes such as hydrolysis of network species, condensation of Si species, and diffusion. However, diffusion has often been overlooked in Monte Carlo models of the aqueous corrosion of borosilicate glasses. Therefore, three different models for dissolved Si diffusion in the altered layer were implemented in a Monte Carlo model and evaluated for glasses in the compositional range (75-x) mol% SiO2 (12.5+x/2) mol% B2O3 and (12.5+x/2) mol% Na2O, where 0 ≤ x ≤ 20%, and corroded in static conditions at a surface-to-volume ratio of 1000 m-1. The three models considered instantaneous homogenization (M1), linear concentration gradients (M2), and concentration profiles determined by solving Fick’s 2nd law using a finite difference method (M3). Model M3 revealed that concentration profiles in the altered layer are not linear and show changes in shape and magnitude as corrosion progresses, unlike those assumed in model M2. Furthermore, model M3 showed that, for borosilicate glasses with a high forward dissolution rate compared to the diffusion rate, the gradual polymerization and densification of the altered layer is significantly delayed compared to models M1 and M2. Models M1 and M2 were found to be appropriate modelsmore »

- Authors:

- Publication Date:

- Research Org.:
- Pacific Northwest National Laboratory (PNNL), Richland, WA (US), Environmental Molecular Sciences Laboratory (EMSL)

- Sponsoring Org.:
- USDOE

- OSTI Identifier:
- 1208740

- Report Number(s):
- PNNL-SA-102399

47580; AF5805010

- DOE Contract Number:
- AC05-76RL01830

- Resource Type:
- Journal Article

- Resource Relation:
- Journal Name: Journal of Non-crystalline Solids, 408:142-149

- Country of Publication:
- United States

- Language:
- English

- Subject:
- Monte Carlo; borosilicates; nuclear waste glasses; Environmental Molecular Sciences Laboratory

### Citation Formats

```
Kerisit, Sebastien N., Pierce, Eric M., and Ryan, Joseph V.
```*Monte Carlo Simulations of Coupled Diffusion and Surface Reactions during the Aqueous Corrosion of Borosilicate Glasses*. United States: N. p., 2015.
Web. doi:10.1016/j.jnoncrysol.2014.07.020.

```
Kerisit, Sebastien N., Pierce, Eric M., & Ryan, Joseph V.
```*Monte Carlo Simulations of Coupled Diffusion and Surface Reactions during the Aqueous Corrosion of Borosilicate Glasses*. United States. doi:10.1016/j.jnoncrysol.2014.07.020.

```
Kerisit, Sebastien N., Pierce, Eric M., and Ryan, Joseph V. Thu .
"Monte Carlo Simulations of Coupled Diffusion and Surface Reactions during the Aqueous Corrosion of Borosilicate Glasses". United States.
doi:10.1016/j.jnoncrysol.2014.07.020.
```

```
@article{osti_1208740,
```

title = {Monte Carlo Simulations of Coupled Diffusion and Surface Reactions during the Aqueous Corrosion of Borosilicate Glasses},

author = {Kerisit, Sebastien N. and Pierce, Eric M. and Ryan, Joseph V.},

abstractNote = {Borosilicate nuclear waste glasses develop complex altered layers as a result of coupled processes such as hydrolysis of network species, condensation of Si species, and diffusion. However, diffusion has often been overlooked in Monte Carlo models of the aqueous corrosion of borosilicate glasses. Therefore, three different models for dissolved Si diffusion in the altered layer were implemented in a Monte Carlo model and evaluated for glasses in the compositional range (75-x) mol% SiO2 (12.5+x/2) mol% B2O3 and (12.5+x/2) mol% Na2O, where 0 ≤ x ≤ 20%, and corroded in static conditions at a surface-to-volume ratio of 1000 m-1. The three models considered instantaneous homogenization (M1), linear concentration gradients (M2), and concentration profiles determined by solving Fick’s 2nd law using a finite difference method (M3). Model M3 revealed that concentration profiles in the altered layer are not linear and show changes in shape and magnitude as corrosion progresses, unlike those assumed in model M2. Furthermore, model M3 showed that, for borosilicate glasses with a high forward dissolution rate compared to the diffusion rate, the gradual polymerization and densification of the altered layer is significantly delayed compared to models M1 and M2. Models M1 and M2 were found to be appropriate models only for glasses with high release rates such as simple borosilicate glasses with low ZrO2 content.},

doi = {10.1016/j.jnoncrysol.2014.07.020},

journal = {Journal of Non-crystalline Solids, 408:142-149},

number = ,

volume = ,

place = {United States},

year = {Thu Jan 01 00:00:00 EST 2015},

month = {Thu Jan 01 00:00:00 EST 2015}

}