Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “constructive solid geometry”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

88 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Unsupervised Learning Based Interaction Force Model for Nonspherical Particles in Incompressible Flows

This project provides a neural network-based interaction force model for gas-solid flows from low to intermediate Reynolds numbers and concentration, which can be linked to MFiX-DEM. We have constructed a database of the interaction force between the irregular-shaped particles using a spherical harmonic method and the fluid phase based on the particle-resolved direct numerical simulation (PR-DNS) with immersed boundary-based gas kinetic scheme. Unsupervised learning method, i.e., variational auto-encoder (VAE) has been applied to extract the primitive shape factors determining the drag force, lifting forces, and torque. The interaction force model has been trained and validated with a simple but effective multi-layer feed-forward neural network: multi-layer perceptron (MLP), which will be concatenated after the encoder of the previously trained VAE for geometry feature extraction for single, irregular particles. We have trained transpose convolutional neural networks with the PR-DNS data to predict the velocity and pressure gradient of the single particle systems and utilized them to calculate drag force of multi-particle systems. This model can provide high computational efficiency because it does not require collecting multiparticle system data from PR-DNS.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2025

Design of a robot-automated flat plate/reflection geometry x-ray diffraction setup for accelerated materials discovery and structural screening

Here, we report the design, construction, and automation of a flat plate sample loading, alignment, and data acquisition system for X-ray diffraction measurements in reflection geometry implemented at the Stanford Synchrotron Radiation Lightsource. The system is built onto a single platform, enabling facile transferability, and is compartmentalized into sample storage, sample transfer, and sample position/alignment segments. The core feature of this system is a six-axis robotic arm that offers a large range of highly reproducible and programable movements. The degrees of freedom of the robot arm enable adaptability in which movements can be modified to fit various beamline environments and sample configurations. Samples are housed on 3D printed sample mounts, which are arranged onto a 6 × 2 array of sample cassettes capable of holding 7 samples. Using sample mounts designed for solid oxide electrolysis button cells (SOECs), the maximum tray capacity is 84 samples, which can be aligned and run in ~ 24 hours with long exposure scans. The sample array is additionally capable of accommodating a range of sample sizes and geometries due to the rapid 3D printed fabrication. The components of the setup will be described in detail and performance will be demonstrated with a set of representative SOEC and XRD standard samples. Opportunities for future developments and integration with the automated setup are summarized.

08 HYDROGEN↗
NASA NTRS2001

Study of Magnetic Damping Effect on Convection and Solidification Under G-Jitter Conditions

As shown in space flight experiments, g-jitter is a critical issue affecting solidification processing of materials in microgravity. This study aims to provide, through extensive numerical simulations and ground based experiments, an assessment of the use of magnetic fields in combination with microgravity to reduce the g-jitter induced convective flows in space processing systems. Analytical solutions and 2-D and 3-D numerical models for g-jitter driven flows in simple solidification systems with and without the presence of an applied magnetic field have been developed and extensive analyses were carried out. A physical model was also constructed and PIV measurements compared reasonably well with predictions from numerical models. Some key points may be summarized as follows: (1) the amplitude of the oscillating velocity decreases at a rate inversely proportional to the g-jitter frequency and with an increase in the applied magnetic field; (2) the induced flow oscillates at approximately the same frequency as the affecting g-jitter, but out of a phase angle; (3) the phase angle is a complicated function of geometry, applied magnetic field, temperature gradient and frequency; (4) g-jitter driven flows exhibit a complex fluid flow pattern evolving in time; (5) the damping effect is more effective for low frequency flows; and (6) the applied magnetic field helps to reduce the variation of solutal distribution along the solid-liquid interface. Work in progress includes developing numerical models for solidification phenomena with the presence of both g-jitter and magnetic fields and developing a ground-based physical model to verify numerical predictions.

Li, Ben Q.↗
NASA NTRS2004

Comparing morphologies of drainage basins on Mars and Earth using integral-geometry and neural maps

We compare morphologies of drainage basins on Mars and Earth in order to confine the formation process of Martian valley networks. Basins on both planets are computationally extracted from digital topography. Integral-geometry methods are used to represent each basin by a circularity function that encapsulates its internal structure. The shape of such a function is an indicator of the style of fluvial erosion. We use the self-organizing map technique to construct a similarity graph for all basins. The graph reveals systematic differences between morphologies of basins on the two planets. This dichotomy indicates that terrestrial and Martian surfaces were eroded differently. We argue that morphologies of Martian basins are incompatible with runoff from sustained, homogeneous rainfall. Fluvial environments compatible with observed morphologies are discussed. We also construct a similarity graph based on the comparison of basins hypsometric curves to demonstrate that hypsometry is incapable of discriminating between terrestrial and Martian basins. INDEX TERMS: 1824 Hydrology: Geomorphology (1625); 1886 Hydrology: Weathering (1625); 5415 Planetology: Solid Surface Planets: Erosion and weathering; 6225 Planetology: Solar System Objects Mars. Citation: Stepinski, T. F., and S. Coradetti (2004), Comparing morphologies of drainage basins on Mars and Earth using integral-ge

Stepinski, T. F.↗
DOE OSTI2021

Closures and Simulation for Thermal Radiation Transport in Stochastic Media with Nonlinear Temperature Dependence

Because of the practical challenge of rendering very complex realistic spatial structures for numerical work, it is common practice to resort to characterizing such media as stochastic mixtures of materials, ideally parametrized with low order statistics such as the mean, variance, and correlation functions of the now random material properties. This enables realizations of the medium to be repeatedly generated and radiation transport computations to, in principle, be performed for a large ensemble of these realizations to obtain a statistically well-characterized radiation field. Statistical post-processing yields desired quantities such as conditional and unconditional mean radiation flux and probability distributions of transmitted radiation. However, such computations prove expensive for all but the simplest stochastic geometries and are most suited for benchmarking approximate models. The most common approximations lead to homogenized media so that transport computations are required only on a single medium realization but by construct provide only limited statistical information on the radiation field. Almost all approximate approaches to this problem attempt to develop equations for low order moments of the radiation intensity (mean, second moment, correlation function) but inevitably encounter a closure problem: the equation for any statistical moment will contain terms depending on unknown higher-order moments. Thus, the challenge shifts to one of developing closure relations that relate the unknown moments to the lower order moments. Under very special conditions, an exact closure can be derived but in general closures are heuristically stated constitutive relations. Also, closure approaches depend on whether the mixing statistics are spatially and/or temporally continuous as in fluctuating turbulent fields, or discontinuous as in randomly mixed solid chunks of material. Thus, unconditional averaging is generally applied in the former case but conditional averaging is more appropriate when the mixing is discontinuous. In this work, the emphasis is on binary statistical mixtures of immiscible fluids as as such the quantities of interest are averages (flux, temperature) conditioned on the material type.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2011

Local Leak Detection and Health Monitoring of Pressurized Tanks

An optical gas-detection sensor safely monitors pressurized systems (such as cryogenic tanks) and distribution systems for leaks. This sensor system is a fiber-coupled, solid optical body interferometer that allows for the miniaturized sensing element of the device to be placed in the smallest of recesses, and measures a wide range of gas species and densities (leaks). The deflection of the fringe pattern is detected and recorded to yield the time-varying gas density in the gap. This technology can be used by manufacturers or storage facilities with toxic, hazardous, or explosive gases. The approach is to monitor the change in the index of refraction associated with low-level gas leaks into a vacuum environment. The completion of this work will provide NASA with an enabling capability to detect gas system leaks in space, and to verify that pressurized systems are in a safe (i.e. non-leaking) condition during manned docking and transit operations. By recording the output of the sensor, a time-history of the leak can be constructed to indicate its severity. Project risk is mitigated by having several interferometric geometries and detection techniques available, each potentially leveraging hardware and lessons learned to enhance detectability.

Polzin, Kurt↗
DOE OSTI2023

Theoretical estimates of flammability bounds for thin condensed fuel diffusion flames in microgravity using detailed models of chemistry and radiation

Recently, U-shape flammability maps have been constructed showing minimal oxygen vs. flame strain for opposed flame spread in micro-gravity by Olson and Ferkul. The U-shape defines the limiting flammability bounds from radiative extinction and flame blow-off. Here, the minimum of the U corresponds to the minimum possible oxidizer concentration where burning can occur, and is an important quantity of interest for fire safety. While high strain extinction bounds have been well analyzed, low strain radiative extinction has not. To estimate low strain extinction, in this study an analytical theory is developed based on thin flame theory coupled with a heat and mass transfer model for solid fuels. A reaction progress variable based on the Damköhler number is adapted in the theory to account for incomplete combustion at high strain rates and enable the capturing of the full flammability map. The analytical model is compared to a one dimensional numerical model w/ detailed chemical kinetics and coupled radiation heat transfer in planar and spherical geometries. The flammability maps are then qualitatively compared to experimental extinguishment data compiled by Olson and Ferkul for cylindrical rods of PMMA showing similar trends. The results show the newly developed analytics capture the radiative extinction bound compared to the numerical model and qualitatively agrees with microgravity data.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2024

Quarterly Research Performance Progress Report (Q8)

As part of Task 1, we have started by testing our modeling capabilities by reproducing isothermal DFIT simulations presented in the literature. Once satisfied with the results we have started by targeting the modeling of the DFITs at conducted at well 58-32. We have a identified a specific test (cycle 4 in zone 2) as the most interesting to be model with GEOS hydraulic fracturing module. Thus, we have first produced results with an isothermal model and adjusted model parameters to get a satisfying match with field pressure data. The, we have added thermal effects and compared the modeling results with and without thermal effects to estimate how thermal effects may influence test interpretation. Models seem to suggest that, for small volumes of fluid, thermal effects are moderate. In Task 2, we have adapted GEOS phase-field formulation to be able to simulate near-wellbore hydraulic fracture nucleation and propagation. We have devised a novel formulation that, compared to other existing ones, incorporates rock strengths. We have submitted a journal publication about our work. We are currently employing this phase-field formulation to model the experiments taking place at U Pitt and help us understand the effect of various parameters. In Task 3, we have built a model of the region surround well 16A and have started modeling stage 3 stimulation because of its simpler planar geometry. After calibrating simulation parameters using known analytical solutions, we have simulated the stage 3 stimulation using our isothermal hydraulic fracturing module, varying the permeability field, the stress conditions including different physics to get a better understanding of the numerical challenges and of the effects of varying these parameters on the simulation results. In Task 4 laboratory experimentation, a set of specialized drilling and injection tools has been customized and constructed to accommodate an inclined well with an orientation of up to 30 degrees relative to material anisotropy or principal stress axes. These inclined samples have also undergone thermal stress and hydraulic fracturing at a temperature of 190 degrees Celsius. Furthermore, both vertical and deviated sampling testing setups enable an extended analysis of post-peak pressure behaviors, facilitating post-test pressure analyses such as the G-function, step rate, and fracture reopening measurements. Thus, the key components of in-situ stress estimation can be extracted and validated through our experiment, providing a solid foundation for validating existing in-situ stress estimation theories or proposing new ones. Simultaneously, we are integrating computer vision techniques with traditional experimental fracture observation methods such as multi-overcore/slicing and water-penetration fracture observation. This combination will prove beneficial in populating the hydraulic fracture patterns database, generated under challenging EGS conditions. This approach aims to deepen our understanding of the complexities in EGS reservoirs and pave the way for future data-driven investigations. Additionally, PITT has also equipped the ELE International compression machine, which is now prepared for conducting indirect tensile and fracture toughness tests. These tests will aid in characterizing how rock fabrics influence the resulting fracture patterns. Additionally, we have completed the required personnel training and gained access to Scanning Electron Microscopy (SEM) and Energy Dispersive Spectroscopy (EDS) for conducting more detailed characterization and analysis of rock fabrics, as well as the examination of thermal and hydraulically induced fracture patterns. Thus, the PITT team has effectively demonstrated the capabilities of our experimental apparatuses in exploring the thermal effects, well deviation angles, material anisotropy, and operational choices (such as circulation rate, injection fluid viscosity, and injection rate) and their impact on pressure responses and fracture trajectories under the Utah FORGE conditions.

15 GEOTHERMAL ENERGY↗
DOE OSTI2022

Crystallography of hyperbolic lattices

Hyperbolic lattices are a revolutionary platform for tabletop simulations of holography and quantum physics in curved space and facilitate efficient quantum error correcting codes. Their underlying geometry is non-Euclidean, and the absence of Bloch's theorem precludes the straightforward application of the often indispensable energy band theory to study model Hamiltonians on hyperbolic lattices. Motivated by recent insights into hyperbolic band theory, we initiate a crystallography of hyperbolic lattices. Here we show that many hyperbolic lattices feature a hidden crystal structure characterized by unit cells, hyperbolic Bravais lattices, and associated symmetry groups. Using the mathematical framework of higher-genus Riemann surfaces and Fuchsian groups, we derive a list of example hyperbolic {p,q} lattices and their hyperbolic Bravais lattices, including five infinite families and several graphs relevant for experiments in circuit quantum electrodynamics and topolectrical circuits. This dramatically simplifies the computation of energy spectra of tight-binding Hamiltonians on hyperbolic lattices, from exact diagonalization on the graph to solving a finite set of equations in terms of irreducible representations. The significance of this achievement needs to be compared to the all-important role played by conventional Euclidean crystallography in the study of solids. We exemplify the high potential of this approach by constructing and diagonalizing finite-dimensional Bloch wave Hamiltonians.

2-dimensional systems↗
NASA NTRS1995

Marshall Space Flight Center surface modeling and grid generation applications

The Solid Rocket Motors (SRM) used by NASA to propel the Space Shuttle employ gimballing nozzles as a means for vehicular guidance during launch and ascent. Gimballing a nozzle renders the pressure field of the exhaust gases nonaxisymmetric. This has two effects: (1) it exerts a torque and side load on the nozzle; and (2) the exhaust gases flow circumferentially in the aft-dome region, thermally loading the flexible boot, case-to-nozzle joint, and casing insulation. The use of CFD models to simulate such flows is imperative in order to assess SRM design. The grids for these problems were constructed by obtaining information from drawings and tabulated coordinates. The 2D axisymmetric grids were designed and generated using the EZ-Surf and GEN2D surface and grid generation codes. These 2D grids were solved using codes such as FDNS, GASP, and MINT. These axisymmetric grids were rotated around the center-line to form 3D nongimballed grids. These were then gimballed around the pivot point and the gaps or overlaps resurfaced to obtain the final domains, which contained approximately 366,000 grid points. The 2D solutions were then rotated and manipulated as appropriate for geometry and used as initial guesses in the final solution. The analyses were used in answering questions about flight criteria.

Williams, Robert W.↗
DOE OSTI2026

AnisONet: A deep neural operator-based anisotropic permeability upscaler from pore to Darcy scale

Directional permeability variations, which govern directional fluid flow in porous media with anisotropy, are important to accurately predict flow behavior, reactive transport, and fluid–solid interactions for various processes such as enhanced geothermal systems, energy storage devices, and biological systems. However, the intricate architecture of porous media makes it difficult to predict directional permeabilities. In this work, we present a novel machine learning (ML) framework, AnisONet, built upon an integration of a convolutional neural network, Swin transformer, and the deep operator network architecture, designed to predict anisotropic permeability and upscale predictions to larger spatial domains. First, AnisONet was evaluated with three classes of two-dimensional (2D) porous media, including synthetic circular and elliptical grains and natural sandstone grains from micro-computed tomography images. A lattice Boltzmann model (LBM) was used to calculate directional permeabilities at every 10° angle, producing 19 data points per image of porous media. AnisONet is then trained to predict permeability as a function of rotation angle. AnisONet showed strong predictive capability of directional permeability. Second, we tested our model for five upscaling cases with a large image size in the finite-element method (FEM) for 2D Darcy flow with various permeability tensor construction methods. Overall, upscaled permeability tensors in FEM simulations produce a reasonably good match with LBM results, highlighting the importance of selecting appropriate tensor formation strategies for accurate permeability upscaling. AnisONet, as a directional permeability estimator, could be further developed for more complex geometries, with the potential to develop a foundational ML model for various applications in porous media.

42 ENGINEERING↗
NASA NTRS2003

Flow Visualization of Low Prandtl Number Fluids using Electrochemical Measurements

It is well established that residual flows exist in contained liquid metal processes. In 1-g processing, buoyancy forces often drive these flows and their magnitudes can be substantial. It is also known that residual flows can exist during microgravity processing, and although greatly reduced in magnitude, they can influence the properties of the processed materials. Unfortunately, there are very few techniques to visualize flows in opaque, high temperature liquid metals, and those available are not easily adapted to flight investigation. In this study, a novel technique is developed that uses liquid tin as the model fluid and solid-state electrochemical cells constructed from Yttria-Stabilized Zirconia (YSZ) to establish and measure dissolved oxygen boundary conditions. The melt serves as a common electrode for each of the electrochemical cells in this design, while independent reference electrodes are maintained at the outside surfaces of the electrolyte. By constructing isolated electrochemical cells at various locations along the container walls, oxygen is introduced or extracted by imposing a known electrical potential or passing a given current between the melt and the reference electrode. This programmed titration then establishes a known oxygen concentration boundary condition at the selected electrolyte-melt interface. Using the other cells, the concentration of oxygen at the electrolyte-melt interface is also monitored by measuring the open-circuit potentials developed between the melt and reference electrodes. Thus the electrochemical cells serve to both establish boundary conditions for the passive tracer and sense its path. Rayleigh-Benard convection was used to validate the electrochemical approach to flow visualization. Thus, a numerical characterization of the second critical Rayleigh numbers in liquid tin was conducted for a variety of Cartesian aspect ratios. The extremely low Prandtl number of tin represents the lowest value studied numerically. Additionally, flow field oscillations are visualized and the effect of tilt on convecting systems is quantified. Experimental studies of the effect of convection in liquid tin are presented. Three geometries are studied: (1) double electrochemical cell with vertical concentration gradients; (2) double cell with horizontal concentration gradients; and (3) multiple cells with vertical temperature gradients. The first critical Rayleigh number transition is detected with geometry (1) and it is concluded that current measurements are not as affected by convection as EMF measurements. The system is compared with numerical simulations in geometry (2), and oscillating convection is detected with geometry (3).

Crunkleton, D.↗
NASA NTRS2004

Vibration Characteristics Determined for Stainless Steel Sandwich Panels With a Metal Foam Core for Lightweight Fan Blade Design

The goal of this project at the NASA Glenn Research Center is to provide fan materials that are safer, weigh less, and cost less than the currently used titanium alloy or polymer matrix composite fans. The proposed material system is a sandwich fan construction made up of thin solid face sheets and a lightweight metal foam core. The stiffness of the sandwich structure is increased by separating the two face sheets by the foam layer. The resulting structure has a high stiffness and lighter weight in comparison to the solid facesheet material alone. The face sheets carry the applied in-plane and bending loads (ref. 1). The metal foam core must resist the transverse shear and transverse normal loads, as well as keep the facings supported and working as a single unit. Metal foams have ranges of mechanical properties, such as light weight, impact resistance, and vibration suppression (ref. 2), which makes them more suitable for use in lightweight fan structures. Metal foams have been available for decades (refs. 3 and 4), but the difficulties in the original processes and high costs have prevented their widespread use. However, advances in production techniques and cost reduction have created a new interest in this class of materials (ref. 5). The material chosen for the face sheet and the metal foam for this study was the aerospace-grade stainless steel 17-4PH. This steel was chosen because of its attractive mechanical properties and the ease with which it can be made through the powder metallurgy process (ref. 6). The advantages of a metal foam core, in comparison to a typical honeycomb core, are material isotropy and the ease of forming complex geometries, such as fan blades. A section of a 17-4PH sandwich structure is shown in the following photograph. Part of process of designing any blade is to determine the natural frequencies of the particular blade shape. A designer needs to predict the resonance frequencies of a new blade design to properly identify a useful operating range. Operating a blade at or near the resonance frequencies leads to high-cycle fatigue, which ultimately limits the blade's durability and life. So the aim of this study is to determine the variation of the resonance frequencies for an idealized sandwich blade as a function of its face-sheet thickness, core thickness, and foam density. The finite element method is used to determine the natural frequencies for an idealized rectangular sandwich blade. The proven Lanczos method (ref. 7) is used in the study to extract the natural frequency.

Ghosn, Louis J.↗
NASA NTRS1993

Flow visualization study in high aspect ratio cooling channels for rocket engines

The structural integrity of high pressure liquid propellant rocket engine thrust chambers is typically maintained through regenerative cooling. The coolant flows through passages formed either by constructing the chamber liner from tubes or by milling channels in a solid liner. Recently, Carlile and Quentmeyer showed life extending advantages (by lowering hot gas wall temperatures) of milling channels with larger height to width aspect ratios (AR is greater than 4) than the traditional, approximately square cross section, passages. Further, the total coolant pressure drop in the thrust chamber could also be reduced, resulting in lower turbomachinery power requirements. High aspect ratio cooling channels could offer many benefits to designers developing new high performance engines, such as the European Vulcain engine (which uses an aspect ratio up to 9). With platelet manufacturing technology, channel aspect ratios up to 15 could be formed offering potentially greater benefits. Some issues still exist with the high aspect ratio coolant channels. In a coolant passage of circular or square cross section, strong secondary vortices develop as the fluid passes through the curved throat region. These vortices mix the fluid and bring lower temperature coolant to the hot wall. Typically, the circulation enhances the heat transfer at the hot gas wall by about 40 percent over a straight channel. The effect that increasing channel aspect ratio has on the curvature heat transfer enhancement has not been sufficiently studied. If the increase in aspect ratio degrades the secondary flow, the fluid mixing will be reduced. Analysis has shown that reduced coolant mixing will result in significantly higher wall temperatures, due to thermal stratification in the coolant, thus decreasing the benefits of the high aspect ratio geometry. A better understanding of the fundamental flow phenomena in high aspect ratio channels with curvature is needed to fully evaluate the benefits of this geometry. The fluid dynamic and conjugate heat transfer problem of high aspect ratio rocket engine coolant channels are being investigated numerically, but these efforts have been hampered by a lack of validating data. Wall temperature data is available for the conjugate problem for channels without curvature and aspect ratio = 5.0, and unheated fluid dynamic data are available for square and circular cross section channels with curvature at Reynold's numbers up to 40,000. But the effects of aspect ratio on secondary flow development have not been experimentally studied. To provide some insight into the effects of channel aspect ratio on secondary flow and to qualitatively provide anchoring for the numerical codes, a flow visualization experiment was initiated at the NASA Lewis Research Center.

Meyer, Michael L.↗
DOE OSTI2025

High-resolution 3D-printed plastic scintillators with tertiary dye

Additive manufacturing offers efficient production of plastic scintillators with nontrivial geometries using vat polymerization, allowing fabrication of geometries which would be difficult or even impossible to produce using conventional subtractive manufacturing. This work presents a novel photocurable scintillator formula that includes coumarin 450 as a tertiary dye to enable high-resolution 3D printing via the manipulation of the 405 nm cure light. Bulk photocured and 3D printed (with and without tertiary dye) samples were compared through observational assessment and spectral response. All samples showed pulse shape discrimination between neutron and gamma events. Inclusion of the tertiary dye has minimal impact on emission spectrum and light output, but significant impact on print resolution as shown by comparison of printed high-complexity geometries and feature resolution test objects. With the use of a cure-limiting dye, unsupported features — such as freestanding pillars — were resolvable down to 0.7 mm. Even finer resolution at or below 0.1 mm was achieved in fully supported, integrated structures printed with off-the-shelf 405 nm desktop 3D printer. Scintillators demonstrated a light output up to 50% of EJ-200 with a PSD figure of merit up to 1.35 at 0.9–1.1 MeVee.

36 MATERIALS SCIENCE↗
NASA NTRSDate not supplied

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.

ISRU↗