Search NASASearch

SEARCH · Search NASA

Results for “numerical simulations”

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.

282 records · Page 12

Three-Dimensional, Multi-Phase Unsteady Flow of Water Undergoing Multiple Flow Regimes (Free Surface, Droplet, and Fluid Film)

Upgrades to the Mobile Launcher Rainbird nozzles were implemented for the Artemis II launch to increase water flow rates and improve the mobile launcher deck water coverage while minimizing the water throw distance to avoid impacting the Space Launch System (SLS) rocket nozzles. This analysis describes the use of a three-dimensional, multi-phase, free surface transient computational fluid dynamics (CFD) model to simulate the water flowing out through the Rainbird nozzle and ensure compliance with current Artemis requirements. Flow rates, pressures, velocities, nozzle throw distances and coverage areas were analyzed to ensure that system requirements are met. The completed results were used to make necessary modifications to the Ignition Overpressure Protection and Sound Suppression (IOP/SS) system Rainbird design by reshaping the nozzle used in the system.

CFD

The Origin of Palimpsests

Palimpsests are the most enigmatic impact features observed on icy worlds of the Solar System. Palimpsests are large impact features which are bright, old, mostly flat, and roughly circular. They range from >150 to ~350 km in diameter and sometimes exhibit recognizable associated radial secondary ejecta crater strings beyond their rims. Curiously, among all the Solar System's icy worlds, palimpsests are exclusively seen on the icy Galilean moons. Palimpsests are obviously different than similarly sized impact features seen on other icy satellites. For example, Evander, a 350 km diameter peak-ring impact basin on Saturn’s moon Dione, has a very low depth/diameter ratio of 0.008, indicating significant viscous relaxation. Nevertheless, it still retains short wavelength topography including a prominent central peak, peak ring, and a continuous rim, none of which are seen at palimpsests. Palimpsests display an approximately circular, and dominating, expanse of undulating plains, which is level with surrounding terrain and shows no recognizable rim. At the center of palimpsests are smooth central plains that are enclosed by a low scarp. Central plains material superposes the surrounding undulating plains and is littered with small blocks reminiscent of “small chaos” on Europa. Palimpsests have yet to be explored using modern impact simulation tools. We have begun to investigate two competing hypotheses: (A) Palimpsests form on targets which have a preexisting near-subsurface layer that is liquid or solid but very weak, and is excavated as proximal ejecta; or (B) They form in a solid ice target underlain by a steep thermal gradient that strongly influences the resulting impact feature, including the production of impact-induced melting at depth that is then a major component of proximal ejecta. Our results will have profound implications for the early thermal history of the icy Galilean satellites. We will test our two hypotheses for the formation of palimpsests on Ganymede (where they are both best observed and preserved) by means of 2D iSALE hydrodynamic simulations of impacts into ice targets that will be structured according to the scenarios for each hypothesis, with select runs in 3D. We have already conducted a pilot study using 2D iSALE that has yielded promising results, which we will report out at the meeting.

palimpsests

Development of the METMAN Thermal Desktop Human Model for System Integration

Human thermal analysis and model development inform NASA’s space suit development, vehicle/habitat design, and survivability studies. Historically, human thermal models such as the 41-node metabolic man (METMAN) and the Wissler model have been standalone software tools developed with FORTRAN, a programming language known for its high performance in computationally intensive applications. Though efficient, these standalone programs pose challenges to coupled human-system analysis with detailed life support and thermal control subsystem models in other commercial software. This paper describes the conversion of the METMAN human thermal model from a standalone FORTRAN program to a model in Thermal Desktop (Ansys), a commonly used CAD-based simulation software for engineering that specializes in heat transfer, thermal radiation, and fluid flow analysis. This format was chosen to best facilitate model sharing and compatibility, enabling the direct integration of METMAN human thermal analysis with subsystem models across NASA programs and commercial partners.

Elizabeth Thurston

Dynamic Analysis of Reynolds Number Effects on Trailing Edge Transonic Vortex Shedding and Its Impact on Turbine Blade Aerodynamic Performance

In this work, we will discuss observations from images acquired from a time-resolved, high-speed self-aligned focusing Schlieren campaign that was performed at the CW-22 facility at NASA Glenn Research Center to understand the dynamic behavior of thick trailing-edge high-pressure turbine blades simulating a ceramic matrix composite (CMC) construction at high inlet turbulence conditions. For the CMC-9 blade, which has a trailing edge thickness of 9% of the axial chord, we identified a regime where an excessive total pressure loss (loss anomaly) is observed only for a narrow range of Reynolds numbers at a fixed exit Mach number of 0.74. The loss anomaly is qualitatively explained by our images, which were taken at 10 distinct blade Reynolds numbers spanning a factor of 6. The images show a significant increase in energy related to the oscillations due to transonic vortex shedding at the Reynolds numbers related to the high loss conditions. This increased energy leads to the formation of strong acoustic waves that turn into shock waves at the highest loss conditions. From our observations stemming from Spectral POD analysis of the high-speed images, we see the acoustic waves produced by the trailing edge vortex shedding exist in all conditions tested; but the shedding frequency has a very slight trend upwards as the Reynolds number is increased, varying about 6% in the range tested. Considering this variation of shedding frequency as a function of Reynolds number, which is well-established for other bluff-body flows, we stipulate there may be a potential feedback mechanism involving an acoustic information transfer path across neighboring blades that may explain why only a narrow range of Reynolds numbers displays strong, shock-forming vortex shedding. We consider a few feedback paths and examine the timing based on the mean flow field from a high-resolution LES simulation. It appears that all feedback mechanisms are viable, presenting an integer number of delay cycles with respect to disturbances generated at the trailing edge. Most noteworthy, however, is the acoustic/shock-boundary layer interaction mechanism at the blade geometric throat. Based on our analysis, this mechanism likely plays an important role in realistic, curved turbine blade passages.

Aerodynamics

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

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

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

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

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.

Evaluation of Cycle Life and Characterization of YTP 45 Ah Li-Ion Battery for EMU

Li-ion batteries, with longer cycle life and higher energy density features, are now more and more attractive and applied in multiple fields. The YTP 45 Ah Li-ion battery has been evaluated here and may be employed in EMU in the future. Evaluations were on: (1) Cycle life tests - 500 cycles total (completed 40 cycles in simulated shuttle use mode and 460 cycles in an accelerated use mode, and recorded differential voltage of individual cell in battery); (2) Characterization test - discharge capacity measurement in environment temperature of -10, 25, 50 C before and after 500 cycles; and (3) Thermal testing - charge and discharge at 50 C and -10 C before and after 500 cycles. The battery showed less than a 9% drop of initial discharge capacity and energy within 500 cycles with 475 cycles 59% DOD plus 25 cycles 100% DOD. The EOD voltage ranged from 16.0 to 18.0 V, which fits the requirement for operating the EMU.

Deng, Yi

Effects of Adverse Weather on Aerodynamics [Les Effets des Conditions Météorologiques Adverses sur I'Aérodynamique]

The 19 technical papers developed for the AGARD Fluid Dynamics Panel (FDP) Specialists' Meeting on "Effects of Adverse Weather on Aerodynamics" held from 28th April—i May 1991 in Toulouse, France are documented in this Conference Proceedings. In addition, introductory material from the Program Chairman and the results of the Round Table Discussion held after the meeting are also included. The FDP organized this meeting to provide a timely review of the progress being made in advancing the state-of-the-art of predicting, simulating, and measuring the effects of icing, anti-icing fluids, and various forms of precipitation on the aerodynamic characteristics of flight vehicles. Topics included results from both theoretical and experimental programs and material related to procedures and regulations for certification and operation. International participation for the meeting included authors from eight nations and representatives from most of the 16 NATO nations.

Precipitation - meteorology

The dynamical origins of the dark comets and a proposed evolutionary track

So-called ‘dark comets’ are small, morphologically inactive near-Earth objects (NEOs) that exhibit nongravitational accelerations inconsistent with radiative effects. These objects exhibit short rotational periods (minutes to hours), where measured. We find that the strengths required to prevent catastrophic disintegration are consistent with those measured in cometary nuclei and expected in rubble pile objects. We hypothesize that these dark comets are the end result of a rotational fragmentation cascade, which is consistent with their measured physical properties. We calculate the predicted size-frequency distribution for objects evolving under this model. Using dynamical simulations, we further demonstrate that the majority of these bodies originated from the 𝜈6 resonance, implying the existence of volatiles in the current inner main belt. Moreover, one of the dark comets, (523599) 2003 RM, likely originated from the outer main belt, although a JFC origin is also plausible. These results provide strong evidence that volatiles from a reservoir in the inner main belt are present in the near-Earth environment.

Aster G Taylor

Updated global and regional trends of stratospheric ozone profiles

We present updated evaluation of stratospheric ozone profile trends in the 60° S–60°N latitude range using long-term ground-based and satellite climate data records, as well as simulations by chemistry-climate models. The trends are evaluated using the LOTUS (Long-term Ozone Trends and Uncertainties in the Stratosphere) regression model. Analyses of satellite data confirm the statistically significant positive ozone trends in the period 2000–2024 in the upper stratosphere of ~1–3% per decade, with larger trends at mid-latitudes compared to the tropics. The trends are slightly positive or close to zero in the middle stratosphere, and mostly negative, -1 to -2% per decade, in the lower stratosphere, but they are not statistically significant. The morphology and magnitude of ozone trends are similar to previous analyses (2000–2020 trends). Ozone trends in 2000–2024 predicted by chemistry-climate model simulations are in good agreement with combined satellite trends. In the upper stratosphere, models predict a slightly stronger ozone recovery than observations. In the lower stratosphere, both models and satellite observations report negative trends in the tropics, while modelled ozone trends are slightly positive at mid-latitudes. Ozone profile trends over several stations estimated from ground-based records capture the same overall vertical pattern of ozone trends as merged gridded satellite datasets. Analyses of regional ozone profile trends in 2003–2024 using merged satellite datasets confirmed the previous observations of a longitudinal structure in ozone trends in the NH mid-latitude stratosphere, with positive trends over Scandinavia and negative trends over Siberia. However, the magnitude of this dipole-like structure is reduced compared to previous analyses.

trends

Powering the Lunar Surface: Managing Dust, Extreme Environments, and Power Needs

Power availability remains one of the primary constraints for lunar surface science. This talk reviews power requirements from previously flown instruments to help prepare future payloads for upcoming CLPS opportunities and highlights the testing and environmental simulation capabilities at NASA JSC that enable reliable lunar payload development. It also outlines the power needs, environmental challenges, and emerging technologies required to support sustained human and robotic operations on the lunar surface as part of NASA’s Moon to Mars strategy. Key challenges include variable solar illumination at polar and equatorial regions, extreme thermal environments, and dust driven degradation that limit current surface power systems. The science data needed for resource identification and landing site planning will allow for the successful preparation of crewed Artemis activities and long-term presence. Building on recent missions, current test infrastructure, and emerging power technology pathways, this presentation equips industry, academia, and government teams with the information needed to design robust lunar payloads, reduce development risk, and fully leverage the increasing cadence of CLPS missions. These developments will form a critical technical foundation for long duration lunar presence and future Mars exploration.

Anastasia Ford

Hubble Space Telescope nickel-hydrogen battery testing: An update

The Marshall Space Flight Center (MSFC) began testing the HST Ni-H2 Six Battery Test and the 'Flight Spare Battery' Tests approximately one year before the launch of the HST. These tests are operated and reported on by the MSFC, but are managed and funded by Goddard Space Flight Center in direct support of the HST program. The HST Ni-H2 batteries are built from Eagle Picher RNH-90-3 cells. The HST EPS (electrical power system) is a direct energy transfer power system. The HST Ni-H2 Six Battery Test is a breadboard of the HST EPS. The batteries in the test are composed of test module cells and packaged into three battery modules identical to the flight modules. This test is the HST EPS testbed. The 'Flight Spare Battery' Test is a simulation of one of the six battery channels on the HST. The cells in the test are from the flight spare lot of cells, which are the same lot of cells that three of the six HST flight batteries are made from. This test is the battery life test for the HST program.

Whitt, Thomas H.

MEAD In-situ Sampling and Testing at Haughton Crater

The Mars Exploration through Analog-site Drilling (MEAD) project is designed to demonstrate the feasibility and scientific value of robotic drilling missions to Mars in the search for evidence of past or extant life. MEAD conducts high-fidelity field simulations integrating drilling, autonomous operations, and in-situ life-detection and mineralogical instruments. The project tests the Signs of Life Detector (SOLID), proposed for the 2019 Discovery Mars Icebreaker mission, and ARIA (Astronaut Raman Instrument for ISRU and Astrobiology), which provides mineral identification and detects trace organics and volatiles, including indicators of organic thermal maturity. Because current and planned Mars drilling missions operate largely “blind,” MEAD also evaluates the use of drill-induced vibrations recorded by deployed geophones to enable local subsurface mapping for improved targeting of drilling and sampling. In August 2025, MEAD completed its first-year field deployment to Haughton Crater in the Canadian High Arctic, an established Mars impact-crater terrestrial analog site, bringing a Honeybee Robotics TRIDENT drill with applied automation, alongside ARIA, SOLID, and FleetSpace geophones.

Haughton Crater