System Modeling of a Lunar Molten Regolith Electrolysis Plant

Introduction: In-Situ Resource Utilization (ISRU) is the process of extracting local resources to produce commodities for propulsion, life support systems, and off-planet construction rather than transporting consumables from Earth. Molten Regolith Electrolysis (MRE) is a novel ISRU method of extracting oxygen gas and metal alloy from lunar regolith. The MRE process involves placing lunar regolith between two electrodes, through which current is passed, to melt the regolith and reduce the metal oxide constituents by direct electrolysis (e.g. FeO, SiO2, MgO, Al2O3) into oxygen gas and metal alloys. The oxygen is liquefied and used as propellant for landers, while the metals (e.g. Ferro-alloys) are further processed and used in structural building materials and parts manufacturing.

A system model was developed that accounted for the major processes of an MRE plant (from excavation of raw materials to storage of products) to assess the feasibility of a lunar MRE plant. The System Engineering and Integration (SE&I) ISRU Modeling and Analysis (SIMA) team utilized its previously documented system sizing model, the Mission Analysis and Integration Tool (MAIT) [1] as framework of the system model. MAIT uses MATLAB/Simulink to integrate subsystem models into a complete system model of the MRE plant. Total mass, volume, and power requirements were computed for numerous iterations of a MRE plant.

System Model:
Figure 1: MRE Plant Block Diagram
The regolith excavation model determines the mass and power needed to excavate sufficient regolith. The preheating auger initiates the regolith heating process before regolith enters the MRE re-actor to reduce the energy required to turn the solid into a molten liquid. The MRE reactor is modeled in COMSOL Multiphysics and based on the research by Dominguez, Sibille, and Schreiner [2, 3, 4]. This preliminary reactor model provides an accurate calculation of thermal equilibrium during electrochemical operation of the reactor system to assess the optimal mass and power required to process the inlet flow of regolith. The model also computes the outlet flowrates of oxygen and molten products. For this analysis, the primary components of the metal alloy considered were iron and silicon. The oxygen is then purified using an Yttrium Stabilized Zirconia (YSZ) electrode, followed by liquefaction using a 90K cryocooler to be stored as liquid oxygen in insulated cylindrical tanks. In future iterations of the system model, the molten metal tapped from the MRE reactor will undergo additional processing or refinement. However, downstream handling of metals is currently a technology gap that is missing a high TRL subsystem model. Therefore, for this analysis, the accumulated metal alloy stream terminates after leaving the MRE reactor.
Study Goals: This analysis investigates multiple input variables to the system to determine the sensitivity of a (near) complete plant at full-scale. This preliminary investigation ran parametric sweeps on the MRE reactor geometry, electrical current supply, layers of multi-layer insulation (MLI) on the reactor, size of the electrodes in the oxygen purification model, and regolith composition (based on landing site location). Three production targets of oxygen (1,000, 10,000, and 50,000 kg/yr) were investigated for this analysis.

The parametric sweeps conducted in this analysis provide valuable insight into the expected impact of the various model inputs on plant size. This information can be used to identify the most critical components of the plant and guide future decisions on allocating funding for research and development, providing subsystem developers with appropriate interfaces with downstream and upstream processes, and assessing the overall feasibility of MRE when compared to other ISRU plants.

References:
[1] Carlson, A. et al. (2024) ICES, ICES-2024-53. [2] Dominguez, D.A., and Sibille, L. (2011) AIAA, AIAA-2011-700. [3] Schreiner, S.S. (2015) MIT, Dissertation. [4] Schreiner, S.S. et al. (2016) ASR, 57(7), pp.1585-1603.