19–22 May 2025
US/Mountain timezone

Upscale rarefied volatile diffusion in porous media: a probability-based pore network modeling approach

20 May 2025, 12:05
15m
Oral Presentation (MS09) Pore-scale modelling MS09

Speaker

sunpeng zhou (College of Engineering, Peking University)

Description

In certain regions of airless bodies, such as the Moon and cometary nuclei, volatile components like water, carbon dioxide, and methanol can persist and migrate within the regolith due to extremely low temperatures[1-5]. Studying the diffusion behavior of volatiles in porous media is crucial for the exploration and exploitation of these extraterrestrial resources.

Under such extremely high vacuum conditions[6], gas molecules undergo Knudsen diffusion, where the average free path is more than 10 orders of magnitude larger than the size of regolith particles[7] (Knudsen number Kn > 1010). At this almost infinitely large Knudsen number, gas molecules rarely collide with each other, and the diffusion trajectories resemble chords (free paths) between solid surfaces, determined solely by porous structure. Additionally, the number of molecules within a single pore is too small to satisfy statistical continuity. Previous studies on Knudsen diffusion have been limited to either simple tube-like structures[8-11] or Monte Carlo simulation[12, 13] that only provide chord length distributions and effective diffusion coefficients. A direct model linking pore structure to dynamic diffusion behavior is still lacking.

We first conduct numerical experiments by employing test particle Monte Carlo (TPMC) simulations, to obtain molecular trajectories and the diffusion coefficients in complex porous media. We find that pore-bodies and pore-throat behaves very differently. Pore-bodies exhibit isotropy (the exit direction of molecules in pores is uniformly distributed) and satisfy the Markov property (the exit probability of molecules in pores is independent of their residence time), while pore-throats are anisotropic and do not satisfy Markov property.

Based on this duality of local diffusive property and the rarefied feature, we develop a probability-based pore network modeling (pb-PNM) approach to upscale Knudsen diffusion from pore-scale to representative elementary volume (REV) scale. The model uses the probability of a molecule residing in a specific structure (i.e., pore or throat) at a given time to represent the spatial distribution and transfer probability of a molecule between different structures to represent the change in the spatial distribution within a specific time period. The relative error between the predicted diffusion coefficient and the numerical experimental results is less than 20% in most cases.

This pb-PNM method leverages the characteristics of gas molecular movement in pores and throats, obtained from TPMC simulations, to scale up to larger core scales. As a result, this method can be used to analyze the effects of adsorption and heterogeneity on the diffusion of rarefied volatiles, and to simulate the unsteady diffusion process with significantly fewer computational resources compared to TPMC simulations.

References 1. Li, S., et al., Direct evidence of surface exposed water ice in the lunar polar regions. Proceedings of the National Academy of Sciences, 2018. 115(36): p. 8907-8912. 2. Jia, Y., et al., Research of Lunar Water-Ice and Exploration for China’s Future Lunar Water-Ice Exploration. Space: Science & Technology, 2023. 3: p. 0026. 3. Ding, C., et al., Moon-Based Ground Penetrating Radar Derivation of the Helium-3 Reservoir in the Regolith at the Chang'E-3 Landing Site. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2023. 16: p. 2764-2776. 4. Reshetnyk, V., et al., Transport Characteristics of the Near-Surface Layer of the Nucleus of Comet 67P/Churyumov–Gerasimenko. Solar System Research, 2021. 55(2): p. 106-123. 5. Rivera-Valentin, E.G., D.G. Blackburn, and R.K. Ulrich, Exploring the effects of overburden on the sublimation and transport of H2O on Iapetus. Icarus, 2012. 220(2): p. 808-820. 6. Heiken, G.H., D.T. Vaniman, and B.M. French, Lunar Sourcebook, a user's guide to the Moon. 1991. 7. Li, C., et al., Characteristics of the lunar samples returned by the Chang’E-5 mission. National Science Review, 2021. 9(2). 8. Steckelmacher, W., Knudsen flow 75 years on: the current state of the art for flow of rarefied gases in tubes and systems. Reports on Progress in Physics, 1986. 49(10): p. 1083-1107. 9. Li, Y., et al., The Mechanism for the Barrier of Lunar Regolith on the Migration of Water Molecules. Journal of Geophysical Research: Planets, 2023. 128(3). 10. Shi, Y., Y.T. Lee, and A.S. Kim, Knudsen Diffusion Through Cylindrical Tubes of Varying Radii: Theory and Monte Carlo Simulations. Transport in Porous Media, 2012. 93(3): p. 517-541. 11. Russ, S., et al., Lambert diffusion in porous media in the Knudsen regime: Equivalence of self-diffusion and transport diffusion. Physical Review E, 2005. 72(3): p. 030101. 12. Liu, J. and J. Wei, Knudsen diffusion in channels and networks. Chemical Engineering Science, 2014. 111: p. 1-14. 13. Zalc, J.M., S.C. Reyes, and E. Iglesia, The effects of diffusion mechanism and void structure on transport rates and tortuosity factors in complex porous structures. Chemical Engineering Science, 2004. 59(14): p. 2947-2960.
Country china
Acceptance of the Terms & Conditions Click here to agree

Primary authors

sunpeng zhou (College of Engineering, Peking University) Dr Ke Xu (Peking University)

Presentation materials

There are no materials yet.