Modeling the Thermal Extraction of Water Ice from Regolith

Modeling has been developed to support the development of volatile extraction technologies for the Moon, Mars, asteroids, or other bodies. This type of modeling capability is important to avoid the high cost of multiple test campaigns in simulated lunar conditions as the hardware design is iterated. The modeling uses the Crank-Nicholson algorithm applied in a two dimensional (2D) axisymmetric (extendable to 3D) finite difference formalism. It uses soil thermal parameters developed from Apollo soil measurements with adaptations for asteroid regolith. Simulations show that it successfully replicates thermal measurements on the surfaces of asteroids and the Moon and helps interpret those measurements to provide insight into the subsurface properties of those bodies. The 2D simulations have provided insight into the cooling of a lunar drill bit and provide a method to determine the original subsurface temperature despite the presence of the warm bit.