Search NASASearch

SEARCH · Search NASA

Results for “Reservoir Simulation.”

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.

586 records · Page 19

An Overview of Experiments and Modeling of Polysiloxane-Coated Thermal Protection Systems for Missions to Mars, Titan, and Beyond.

Phenolic Impregnated Carbon Ablator (PICA) gained heritage during the historic Stardust mission, where it successfully returned samples from a comet’s tail and has since been instrumental in delivering payloads to the surface of Mars [1-3]. Most recently, PICA enabled the safe return of samples collected from asteroid Bennu as part of the OSIRIS-REx mission. This rich legacy underscores PICA’s critical role in allowing NASA’s most ambitious exploration missions. However, the friable nature of its phenolic phase presents challenges during handling and pre-launch activities. To mitigate this issue, PICA is coated with a polysiloxane resin system, which serves to suppress particulate dispersion and thereby safeguard spacecraft components. A comprehensive understanding of the polysiloxane resin’s behavior is imperative, as it profoundly shapes the material response of PICA during atmospheric entry by influencing its thermal and oxidative stability. This influence extends to thermocouple plugs embedded within thermal protection systems. These plugs have demonstrated their significance in missions such as Mars Science Laboratory (MSL) and Mars 2020, where the MEDLI and MEDLI2 instrumentation suites delivered in-valuable insights into the performance of thermal protection systems during entry into the Martian atmosphere [4]. Looking ahead, missions such as Dragonfly, set to descend into Titan’s dense atmosphere, aim to leverage advanced sensor technologies to further refine our understanding of thermal protection response [5]. Moreover, thermocouple plugs play an essential role in validating cutting-edge material response models, such as those pioneered under NASA’s Entry Systems Modeling Project (ESM), designed, in-part, to predict the operational integrity of thermal protection systems under the extreme stresses of atmospheric entry. To achieve these modeling goals, ground-based experiments are crucial to provide the foundational data necessary for developing and refining these predictive tools. To this end, an extensive test campaign was conducted at the Hypersonic Materials Environmental Test System (HyMETS) to investigate the high-temperature behavior of the polysiloxane resin in an air environment [6]. These experiments revealed critical phenomena, including the formation of a silicon oxycarbide layer that enhances oxidation resistance, moderates surface temperatures, and alters in-depth thermal response. Building on these findings, subsequent tests were designed to simulate atmospheric entry conditions in reactive gases, such as CO2 and N2, to mimic the environments of Mars and Titan, respectively, as well as non-reactive gases representing the atmospheres of the Ice Giants (Neptune and Uranus). A heating rate dependent decomposition mechanism has been identified for the polysiloxane resin under oxidizing conditions (Fig. 1). In the initial stage, the resin and the underlying thermal protection system undergo pyrolysis, rapidly generating a thin amorphous silicon oxycarbide interwoven with carbonaceous char and residual fibers from PICA. During the second stage, the nascent oxide layer establishes a robust, oxidation-resistant thermal barrier coating, which significantly impedes heat transfer to the underlying carbonaceous char, resulting in a stagnation of the surface temperature. A key factor contributing to this thermal resistance is the low recombination efficiency of atomic oxygen (γ), which further diminishes the heat load on the material’s interior layers [7]. Moreover, as the surface temperature stagnates, the silicon oxycarbide phase separates into distinct regions of silica and free graphite. Ultimately, when the heat flux reaches a critical threshold, a third stage is triggered, leading to the breakdown of the coating through carbothermal reduction, exposing the underlying char layer. This exposure leads to a dramatic surface temperature spike, driven by highly exothermic reactions between atomic oxygen and the char layer, further accelerating material degradation. A detailed mass and heat transfer model of PICA coated with polysiloxane resin was implemented in the Porous material Analysis Toolbox based on OpenFOAM, PATO [8]. The initial stage was considered negligible in this model because the resin decomposition occurs rapidly within a thin surface layer. Instead, the coating was directly considered as an oxygen-resistant thermal barrier coating. For the second stage, the thin amorphous silicon oxycarbide was treated as a pure silica surface to simplify the thermochemical behavior. The model ac-counts for surface equilibrium processes using representative elements of the coating-environment system. For the third stage, specific boundary conditions were developed to estimate the onset and progression of the coating removal. Two-dimensional material response simulations were conducted to compare uncoated and coated PICA using boundary conditions calibrated with HyMETS data. Fig. 2 illustrates that the simulations closely align with experimental data, successfully reproducing measured temperature profiles. This work will include the latest advancements in the coating model, including the calibration of recombination of atomic oxygen at the surface during the second phase. These simulated results will be further validated against additional CO2 data points from HyMETS, reinforcing the models’ predictive capabilities. These mechanisms and their effects on thermal protection systems, including thermochemical behavior and thermocouple probe performance in extreme environments, provide crucial insights for optimizing spacecraft designs that safeguard scientific payload and ensure mission success in future planetary exploration endeavors.

Active Oxidation

Challenges in Continuous In-Field Critical Current Testing of High-Temperature Superconducting Tapes: Thermal and Mechanical Perspectives

High-temperature superconductors (HTS) are essential for ultra-high-field applications requiring exceptional current-carrying capacity under extreme conditions. However, systematic characterization of critical current in long-length conductors remains challenging due to complex thermal, electromag netic, and mechanical interactions during continuous testing. This study reports the development of a continuous in-field magnetization testing system for position-dependent critical current measurement in HTS tapes at 20 K under 7.5 T fields applied normal to the tape plane, enabling identification of performance-limiting regions that could compromise magnet stability. Here, the system addresses two fundamental challenges inherent to cryogenic reel to-reel testing. First, thermal management requires continuous cooling of a moving conductor to 20 K, achieved through liquid nitrogen precooling combined with a 100 W@20 K Gifford McMahon cryocooler. Second, screening currents in high fields generate Lorentz forces that induce twisting, bowing, and potential delamination. To mitigate these risks, we propose mechanical reinforcement and active current density suppression strategies. Numerical simulations using the stream function formulation reveal four primary failure modes: frictional heating at guide interfaces, unstable equilibria causing deformation, transverse current-induced stresses at guide transitions, and unsupported forces in vertical spans. Our mitigation strategies include PTFE coated guides to minimize friction, spring-loaded stabilization mechanisms to maintain tape alignment, controlled pre-heating using the liquid nitrogen thermal jacket to suppress critical current at stress points, and optimized guide positioning to minimize force accumulation. The experimental system is nearing completion, with testing planned to commence within two months. Preliminary validation at 65 K under 0.5 T demonstrates strong correlation between simulation-predicted mechanical instabilities and observed critical current variations during conductor tran sitions through the measurement region. These findings establish a robust foundation for quality assurance protocols essential to next-generation superconducting magnet applications.

Chen, Siwei [Princeton Plasma Physics Laboratory (

Computational Fluid Dynamics at the Lewis Research Center: An Overview

Lewis is a multidisciplinary Center with strong research and development programs in aeronautical and space propulsion, power, space communications, space experiments and materials. Computational fluid dynamics (CFD) is playing an important and growing role in most of these areas. Described here is how CFD is integrated into these programs and highlights elements of the CFD activities. Examples are presented of codes developed to predict flow fields in advanced propulsion systems and several of the code validation experiments are described. The CFD effort at Lewis ranges from basic research on new and improved algorithms through code development to the application of these codes to specific engineering problems. Because of the substantial improvement in CFD's predictive capability, its use at Lewis is on a steep growth path, spreading rapidly into new areas which had not traditionally taken advantage of the techniques of numerical simulation. Multidisciplinary codes and the future direction of CFD at Lewis are discussed.

Robert M Stubbs

Lunar and Planetary Science XXXV: Icy Worlds: Moving and Grooving

Reports from the conference session entitled Icy Worlds: Moving and Grooving, include:Mass Anomalies on Ganymede; Europan Chaos and Lenticulae: A Synthesis of Size, Spacing, and Areal Density Analyses; Thermal and Topographic Tests of Europa Chaos Formation Models; Flexure of Europa s Lithosphere Due to Ridge-Loading; Ridges on Europa: Origin by Incremental Ice-Wedging ; Convergent Boundaries on Europa: a Numerical Approach to Euler Pole Analysis and Its' Implications for Plate Reconstruction; Numerical Simulations of Subsolidus Convection in the Ice Shell of Europa: Implications for the Thermal Evolution and Present State; Effects of Plasticity on Convection in an Ice Shell: Implications for Europa; Non-Newtonian Convection and Compositional Buoyancy: Advances in Modeling Convection and Dome Formation on Europa; Convective Instability in Ice I: Application to Callisto and Ganymede; Crater Size Distributions on Callisto: A Galileo SSI Summary; Neutron Diffraction Studies of Planetary Ices; and H2O2 Synthesis Induced by Irradiation of H2O with Energetic H+ and Ar+ Ions at Various Temperatures.

Source record

Magnetically Arrested Circumbinary Accretion Flows

Abstract Binary systems with comparable masses and a surrounding accretion disk can accrete gas through spiral accretion streams penetrating the central cavity formed by tidal interactions. Using three-dimensional Newtonian magnetohydrodynamics simulations, we investigate the possibility of a magnetically arrested accretion flow through the cavity. Rather than solely continuously feeding the binary through spiral accretion streams, the accretion is regulated by the strong magnetic field inside the cavity. Transport of mass and angular momentum onto the binary then proceeds largely periodically in magnetic flux eruption episodes. The ejected flux tubes carry angular momentum outward and away from the binary, inject hot plasma into the disk, and can launch flares. This likely intermittent scenario could have potential implications for the emission signatures of supermassive black hole binaries and shed light onto the role magnetic fields play in the binary’s orbital evolution.

Astronomy & Astrophysics

Numerical Investigation of Fluid Flow and Space Charge in Liquid Argon Time Projection Chamber (LArTPC) Detectors

Overview This project focused on developing a high-fidelity numerical framework to simulate the multiphysics environment within Liquid Argon Time Projection Chamber (LArTPC) detectors. The primary objective was to characterize the complex interplay between ion transport, background fluid dynamics, and electric field distortions—a critical factor for the calibration and sensitivity of next-generation High Energy Physics experiments, such as DUNE. Technical Achievements The research successfully yielded a hybrid numerical space-charge solver utilizing a Cell-Centered Finite Volume Method (FVM) for ion transport coupled with a Finite Element Method (FEM) for electric potential. Key accomplishments include: • Verification & Validation: The 3-D solver was rigorously verified against 1-D analytical solutions, demonstrating high numerical accuracy in predicting space-charge-induced field deviations. • Field Distortion Analysis: 3D simulations revealed that space charge effects introduce significant non-uniformities in the electric field. Critically, the research identified that background LAr flow velocities, when comparable to ion drift velocities, markedly exacerbate these distortions. • Technology Transfer: The resulting source code and comprehensive user manuals were successfully transferred to collaborators at Fermilab, providing a portable computational tool for the broader scientific community. Challenges and Future Directions While the space-charge solver achieved all performance metrics, the integrated fluid dynamics modeling encountered convergence challenges stemming from the extreme 200-fold disparity in length scales between the detector's 37 mm inlet pipes and the 8-meter global domain. To address this, the project has identified a clear technical pivot toward Hierarchical Geometric Adaptive Mesh Refinement (HG-AMR). By implementing an h-type refinement strategy with hanging nodes, future iterations of this solver will be capable of resolving localized high-gradient inlet flows without the prohibitive computational costs of regular grids. This advancement, combined with data-driven uncertainty quantification based on MicroBooNE-style calibration, will enable the precise modeling of detector responses in large-scale cryogenic environments where direct measurement remains difficult. Impact The computational tools developed under this award provide a foundation for enhancing the energy resolution and spatial reconstruction of noble liquid detectors. By bridging the gap between theoretical fluid dynamics and experimental field calibration, this work supports the DOE’s mission to advance the frontiers of neutrino physics and dark matter detection.

42 ENGINEERING

Discharge Transient Coupling in Large Space Power Systems

Experiments have shown that plasma environments can induce discharges in solar arrays. These plasmas simulate the environments found in low earth orbits where current plans call for operation of very large power systems. The discharges could be large enough to couple into the power system and possibly disrupt operations. Here, the general concepts of the discharge mechanism and the techniques of coupling are discussed. Data from both ground and flight experiments are reviewed to obtain an expected basis for the interactions. These concepts were applied to the Space Station solar array and distribution system as an example of the large space power system. The effect of discharges was found to be a function of the discharge site. For most sites in the array discharges would not seriously impact performance. One location at the negative end of the array was identified as a position where discharges could couple to charge stored in system capacitors. This latter case could impact performance.

N John Stevens

Opportunities to Intercalibrate Radiometric Sensors from International Space Station

Highly accurate measurements of Earth's thermal infrared and reflected solar radiation are required for detecting and predicting long-term climate change. Consideration is given to the concept of using the International Space Station to test instruments and techniques that would eventually be used on a dedicated mission, such as the Climate Absolute Radiance and Refractivity Observatory (CLARREO). In particular, a quantitative investigation is performed to determine whether it is possible to use measurements obtained with a highly accurate (0.3%,with 95%confidence) reflected solar radiation spectrometer to calibrate similar, less accurate instruments in other low Earth orbits. Estimates of numbers of samples useful for intercalibration are made with the aid of yearlong simulations of orbital motion. Results of this study support the conclusion that the International Space Station orbit is ideally suited for the purpose of intercalibration between spaceborne sensors.

Intercalibrate

Topological prethermal strong zero modes on superconducting processors

Abstract Symmetry-protected topological phases 1–4 cannot be described by any local order parameter and are beyond the conventional symmetry-breaking model 5 . They are characterized by topological boundary modes that remain stable under symmetry respecting perturbations 1–4,6–8 . In clean, gapped systems without disorder, the stability of these edge modes is restricted to the zero-temperature manifold; at finite temperatures, interactions with mobile thermal excitations lead to their decay 9–11 . Here we report the observation of a distinct type of topological edge mode 12–14 , which is protected by emergent symmetries and persists across the entire spectrum, in an array of 100 programmable superconducting qubits. Through digital quantum simulation of a one-dimensional disorder-free stabilizer Hamiltonian, we observe robust long-lived topological edge modes over up to 30 cycles for a wide range of initial states. We show that the interaction between these edge modes and bulk excitations can be suppressed by dimerizing the stabilizer strength, leading to an emergent U(1) × U(1) symmetry in the prethermal regime of the system. Furthermore, we exploit these topological edge modes as logical qubits and prepare a logical Bell state, which exhibits persistent coherence, despite the system being disorder-free and at finite temperature. Our results establish a viable digital simulation approach 15–18 to experimentally study topological matter at finite temperature and demonstrate a potential route to construct long-lived, robust boundary qubits in disorder-free systems.

Science & Technology - Other Topics

The Development of Kinetic and Radiation Hydrodynamics Modeling of Thermonuclear Burn Propagation in Isochoric p - 11 B Through the Support of the INFUSE Program

The report summarizes DOE INFUSE-supported work between HB11 Energy and the University of Rochester’s TriForce Institute to improve computational modeling of advanced fusion fuels, especially proton–boron-11 (p- 11 B). The project extended the TriForce particle-in-cell/Monte Carlo collision code to include physics needed for dense, high-temperature p- 11 B burn studies, including p- 11 B fusion reactions, three-alpha-particle reaction products, relativistic Coulomb collisions, large-angle nuclear scattering, bremsstrahlung radiation, inverse bremsstrahlung absorption, and photon transport. The upgraded models were verified against focused physics tests and against known deuterium–tritium burn behavior. The study then used one-dimensional spherical simulations to estimate the conditions required for thermonuclear burn propagation in isochoric p- 11 B fuel. The calculations found that burn propagation is possible in the model, but only under very extreme hot-spot conditions, such as about 7000 g/cm 3 at 500 keV or 9000 g/cm 3 at 300 keV for a 20-micron hot spot. These conditions are much more demanding than current demonstrated inertial confinement fusion hot spots. The report concludes that the INFUSE collaboration successfully advanced kinetic and radiation modeling capabilities for p- 11 B fusion and provided useful estimates of ignition requirements. However, the simulated fuel gains remain below what would be needed for practical inertial fusion energy, and further work is needed to reconcile differences among kinetic, radiation-hydrodynamic, and analytic models and to identify more achievable target designs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY

Design Criteria for Integrated Flight/Propulsion Control Systems for STOVL Fighter Aircraft

As part of NASA's program to develop technology for short takeoff and vertical landing (STOVL) fighter aircraft, control system designs have been developed for a conceptual STOVL aircraft. This aircraft is representative of the class of mixed-flow remote-lift concepts that was identified as the preferred design approach by the US/UK STOVL Joint Assessment and Ranking Team. The control system designs have been evaluated throughout the powered-lift flight envelope on Ames Research Center's Vertical Motion Simulator. Items assessed in the control system evaluation were: maximum control power used in transition and vertical flight, control system dynamic response associated with thrust transfer for attitude control, thrust margin in the presence of ground effect and hot gas ingestion, and dynamic thrust response for the engine core. Effects of wind, turbulence, and ship airwake disturbances are incorporated in the evaluation. Results provide the basis for a reassessment of existing flying qualities design criteria applied to STOVL aircraft.

James A Franklin

Analysis of Surface and Bulk Behavior in Ni-Pd Alloys

The most salient features of the surface structure and bulk behavior of Ni-Pd alloys have been studied using the BFS method for alloys. Large-scale atomistic simulations were performed to investigate surface segregation profiles as a function of temperature, crystal face, and composition. Pd enrichment of the first layer was observed in (111) and (100) surfaces, and enrichment of the top two layers occurred for (110) surfaces. In all cases, the segregation profile shows alternate planes enriched and depleted in Pd. In addition, the phase structure of bulk Ni-Pd alloys as a function of temperature and composition was studied. A weak ordering tendency was observed at low temperatures, which helps explain the compositional oscillations in the segregation profiles. Finally, based on atom-by-atom static energy calculations, a comprehensive explanation for the observed surface and bulk features will be presented in terms of competing chemical and strain energy effects.

Ordering

Two-mode bosonic quoctit for high energy physics

In this work, we study a two-mode bosonic encoding of a quoctit inside a non-Abelian group-structured constellation of coherent states. This work is motivated by the importance of non-Abelian symmetry in particle physics and the desire to have transversal non-Abelian logical gates. We use the previously developed 2 T constellation of states used to encode a so-called 2 T qutrit. The fidelity of the 2 T quoctit is benchmarked against other bosonic qudits for different noise models and find it compare favorably when power constraints are considered. This paves the way for the construction of higher-dimensional qudits (e.g., a quicosotetrit with 2 T group structure) in bosonic systems with practical applications in quantum simulations of particle physics.

Kürkçüoglu, Doga Murat [Fermilab] (ORCID:000000031

Thermoplastic Matrix Composite Design for Cryotanks Using Multiscale Modeling and Bayesian Optimization

Designing lightweight, robust cryogenic storage tanks is critical for future launch vehicles, in-space propellant storage, and hydrogen powered aircraft. This work presents a multiscale modeling and Bayesian optimization framework for the design of thermoplastic matrix composite cryotanks. Molecular dynamics simulations are first used to determine temperature-dependent constituent properties for candidate thermoplastic matrices, which are homogenized to the lamina scale using NASA’s Multiscale Analysis Tool (NASMAT). These lamina properties, in combination with laminate family generation rules, are evaluated in HyperX structural optimization software to identify stacking sequences that meet all cryogenic load requirements. A Bayesian optimization framework is applied, with HyperX in the loop (via the HyperX API) to efficiently search across material and laminate design variables, yielding an optimized cryotank configuration with significant reductions in design cycle time compared to exhaustive search approaches.

thermoplastics

Development of a high fidelity CFD model for solvent evaporation and transport in porous structure during battery electrode drying

An efficient battery manufacturing process is the key to the mass production of Electric Vehicles (EV), in which drying is one of the most energy-intensive steps significantly influencing the battery cell performance. An accurate 3D CFD model for drying is essential for predicting the drying mechanism and optimizing its parameters. By optimizing the drying process, it is possible to reduce energy consumption and cost during battery manufacturing, minimize binder loading and maximize active material loading to achieve superior electrochemical performances and facilitate wider and faster public adoption of EV. This project aims to optimize the drying process during electrode manufacturing by leveraging high-fidelity, porous electrode simulations for solvent evaporation. By optimizing this process, we seek to reduce energy consumption during battery manufacturing, while minimizing binder loading and maximizing active material loading, with the overall goal of enhancing electrical vehicle performance.

Horner, Jeffrey Scott [Sandia National Laboratorie

Digital Twin Based Condition Monitoring of LCC-LCC Inductive Power Transfer Systems

Inductive power transfer (IPT) systems provide a flexible, hands-free charging opportunity to electric vehicles (EV). The resonant network components and the transmitter and receiver coils are often subjected to high voltages or currents. Component aging in the compensation network and coils of resonant IPT systems is detrimental to the reliability and power transfer efficiency of the IPT system. Monitoring the component health of such multi-element complex systems requires robust optimization algorithms. This paper discusses condition monitoring of a resonant IPT system for an EV charger using a digital twin model. A hybrid estimation algorithm based on genetic algorithms and adaptive particle swarm optimization is developed to estimate the parameters of the digital twin model. Simulation results are used to verify the monitoring capabilities of the developed algorithm under various operating conditions of the IPT system.

Weldehawaryat, Lidya Mussie [graduate research ass

A nickel-cadmium battery reconditioning circuit

The circuit presented is simple and small enough to be included in a typical battery charge/power control assembly, yet provides the advantage of a complete ground-type battery reconditioning discharge. Test results on the circuit when used to recondition two 24 cell, 20 A-h nickel-cadmium batteries are given. These results show that a battery reconditioned with this circuit returns to greater than 90 percent of its original capacity (greater than nameplate capacity) and follows a typical new battery degradation curve even after over 20,000 simulated orbital cycles for a 4 year period. Applications of the circuit are considered along with recommendations relative to its use. Its application in low voltage (22 to 36 Vdc) power systems and in high voltage (100 to 150 Vdc) power systems is discussed. The implications are that the high voltage systems have a greater need for battery reconditioning than their low voltage counterparts, and that using these circuit techniques, the expected life of a battery in low Earth orbit can be up to 5 years.

Lanier, R.

Temperature-dependent spin dynamics in Cr 2 O 3

Cr 2 ⁢O 3 emerges as a prominent candidate material for spintronics and magnetoelectronic applications. However, a comprehensive understanding of its temperature-dependent spin dynamics remains elusive, impeding the engineering of novel spintronics based on this material. Here, we delve into this through a combination of inelastic neutron scattering experiments and atomistic simulations. Our results unveiled the emergence of paramagnons above and below 𝑇 𝑁 . We demonstrated a significant softening of linear magnons upon heating in the antiferromagnetic state. Further analysis revealed that this softening primarily originated from four-magnon interactions, while thermal expansion played a minor role.

antiferromagnetism