Multi-species diffusion analysis in porous media under various dry density and temperature conditions: molecular dynamics simulation, homogenization analysis, and finite element method

• Published in: International journal for numerical and analytical methods in geomechanics 2014-11, Vol.38 (16), p.1744-1760
• Main Authors: Choi, Jung Hae, Kawamura, Katsuyuki, Kimoto, Kazushi, Ichikawa, Yasuaki, Chae, Byung-Gon
• Format: Article
• Language: English
• Subjects: Applied sciences; Buildings. Public works; Computation methods. Tables. Charts; Computer simulation; Density; Drying; Exact sciences and technology; finite element method; Geotechnics; homogenization analysis; Homogenizing; Hydroxyapatite; Mathematical analysis; Mathematical models; Media; molecular dynamics; multi-species diffusion; Soil investigations. Testing; Structural analysis. Stresses
• ISSN: 0363-9061; 1096-9853
• DOI: 10.1002/nag.2303

SUMMARY In porous media, chemical species that dissolve in pore water can be transported via diffusion mechanisms or advective fluxes, close to or far away from where precipitation occurs. In the case of a high‐level radioactive waste disposal system, compacted bentonite is used in a buffer material in an engineering barrier system to minimize the amount of specific nuclides that breach into the surrounding host rock. To minimize breaching, it is very important to understand the transport mechanism of multiple chemical species in porous media. In the following research, we introduced FEM analysis methods using the results of the molecular dynamics simulation and homogenization analysis (MD/HA) method. First, the diffusion coefficients of ions (Cl−, I−, and Na+) in different water layers of Na‐beidellite were calculated using the MD/HA procedure under various dry density (1.2, 1.6, and 2.0 Mg/m3) and temperature (293, 323, and 363 K) conditions. Next, using FEM analysis that used the MD/HA results as input parameters, the diffusion behaviors of ions in porous media were calculated. The results indicate that the diffusion coefficients of the interlayer water in Na‐beidellite are different from the diffusion coefficients under dry density conditions. Further, the concentration profiles (Ct/C0) of iodine and chloride are proportional to temperature but inversely proportional to dry density. Copyright © 2014 John Wiley & Sons, Ltd.