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.

631 records · Page 17

Introduction to Radiation Issues for International Space Station Extravehicular Activities: Chapter 1

The International Space Station (ISS) provides significant challenges for radiation protection of the crew due to a combination of circumstances including: the extended duration of missions for many crewmembers, the exceptionally dynamic nature of the radiation environment in ISS orbit, and the necessity for numerous planned extravehicular activities (EVA) for station construction and maintenance. Radiation protection requires accurate radiation dose measurements and precise risk modeling of the transmission of high fluxes of energetic electrons and protons through the relatively thin shielding provided by the space suits worn during EVA. Experiments and analyses have been performed due to the necessity to assure complete radiation safety for the EVA crew and thereby ensure mission success. The detailed characterization described of the material and topological properties of the ISS space suits can be used as a basis for design of space suits used in future exploration missions. In radiation protection practices, risk from exposure to ionizing radiation is determined analytically by the level of exposure, the detrimental quality of the radiation field, the inherent radiosensitivity of the tissues or organs irradiated, and the age and gender of the person at the time of exposure. During low Earth orbit (LEO) EVA, the relatively high fluxes of low-energy electrons and protons lead to large variations in exposure of the skin, lens of the eye, and tissues in other shallow anatomical locations. The technical papers in this publication describe a number of ground-based experiments that precisely measure the thickness of the NASA extravehicular mobility unit (EMU) and Russian Zvezda Orlan-M suits using medical computerized tomography (CT) X-ray analysis, and particle accelerator experiments that measure the minimum kinetic energy required by electrons and photons to penetrate major components of the suits. These studies provide information necessary for improving the understanding of the current ISS space suits and provide insights into improved approaches for the design of future suits. This chapter begins with a summary of the dynamic ionizing radiation environment in LEO space and introduces the concepts and quantities used to quantify exposure to space radiation in LEO. The space suits used for EVA and the experimental partial human phantom are described. Subsequent chapters report results from measured charged particle fields before and after incident protons and secondary particles are transported through the space suits and into organs and tissues.

M R Shavers

Controlled Parametric Forcing During Directional Solidification of a Bulk Organic Alloy Under Microgravity

The response of dendritic microstructures to step-like pulling velocity conditions is investigated using microgravity directional solidification experiments conducted on DECLIC-DSI combined with phase-field simulations. Under a constant pulling velocity of 1.5 µm/s, the evolution toward steady-state growth is characterized in terms of primary spacing, dendrite drift, and tip dynamics. For the first time, side-view observations enabled direct measurement of tip radius and sidebranching frequency. When step-like oscillations of the pulling velocity are imposed, the dendritic array exhibits a strongly period-dependent response: short periods lead to rapid tip adaptation, whereas longer periods induce a phase lag between tip position and morphology, resulting in progressive tip flattening and, above a critical period, interface destabilization and dendrite splitting. Quantitative phase-field simulations, including a realistic thermal field and stochastic noise, reproduce the experimental observations and provide insight into the governing mechanisms, highlighting the role of characteristic relaxation times, sequence-dependent effects, and the irreversible reorganization of the microstructure following splitting.

Microgravity

Mesh-based multiphysics coupling acceleration for fusion neutronics through clustering for fusion blanket applications

Accurate modeling of particle transport within fusion blankets is essential for predicting performance metrics such as heat deposition and the tritium breeding ratio (TBR). However, high-fidelity coupling of thermal fluids from computational fluid dynamics (CFD) to neutronics simulations often incurs significant computational costs due to the complexity of surface intersection calculations in Monte Carlo codes. This paper presents an accelerated multiphysics coupling method for neutronics that utilizes hierarchical agglomerative clustering to map complex material property distributions to a neutronics model. Implemented within the fusion reactor design and assessment (FREDA) framework, the method leverages existing Python packages to automate the creation of clustered geometries for OpenMC. The approach is demonstrated on a sector model of an ARC-class tokamak with an immersion molten salt blanket, and an simple geometry with varying isotopic concentrations. Results show that the clustering method significantly reduces computational burden without compromising fidelity, providing a foundation for agile iteration of neutronics simulations involving multiple coupled material properties.

Bae, Jin Whan [ORNL] (ORCID:0000000326548907)

Optimization of a 10–Step Reduced CH₄/O₂ Combustion Mechanism for RDRE Applications

A 10-step, 9-species reduced mechanism for methane/oxygen combustion is evaluated and optimized for use in rotating detonation rocket engine (RDRE) simulations at an operating pressure of 10 bar. The mechanism is assessed against GRI-Mech 3.0 for ignition delay times, adiabatic equilibrium temperature and species, Chapman-Jouguet (CJ) detonation parameters, and laminar flame speed. A sequential Nelder–Mead optimization procedure is applied to the rate constants of reactions R1, R4, and R5, targeting ignition delay agreement over the equivalence ratio range φ = 1.0–1.4 and temperature range T = 1400–2000 K. The optimized mechanism achieves ignition delay agreement within 28% of GRI-Mech 3.0 for the temperature range relevant to RDRE operation (τ ≥ 5 μs). Adiabatic equilibrium temperatures and CJ detonation parameters are reproduced within 0.15% of GRI-Mech 3.0 values, confirming thermodynamic consistency of the mechanism. Laminar flame speed predictions deviate significantly from GRI-Mech 3.0, which is attributed to the absence of HO₂ chemistry in the reduced mechanism; this limitation does not affect RDRE simulation fidelity since laminar flame propagation is of secondary importance in RDREs.

Computational Fluid Dynamics

Solar-Wind Bombardment of a Surface in Space

The solar wind is described and its effects on surfaces are reviewed. Experimental results are given for materials subjected to simulated solar-wind bombardment. Among the effects discussed are sputtering-type erosion, chemical reaction with paint vehicles and reduction of oxide pigments, and production of lunar-surface optical characteristics in a basalt powder.

Solar Wind

Assembly bias and local Primordial non-Gaussianity from DESI DR1 quasars

The analysis of the large-scale clustering of quasars (QSO) observed by the Dark Energy Spectroscopic Instrument (DESI) represents a promising avenue for constraining local Primordial non-Gaussianity (PNG), parameterized by f NL . The signal to be constrained is the scale-dependent bias induced in the 2-point clustering of the considered tracer sample. The resulting constraints on f NL , however, are fully degenerate with the local PNG bias parameter b ϕ , dependent on the assembly bias parameter p. Using IllustrisTNG hydrodynamical simulations, we select a QSO sample reflecting the selection criteria and properties of DESI QSOs, and provide a robust prior for p, and thus for b ϕ , building on the findings of Fondi et al. 2024. We find a distribution with mean p̅ ≃ 1.4 with weak redshift dependence, stable to selection noise and consistent with the expected recent merger history typical of quasar-hosting halos. By comparing with the CAMELS simulations we demonstrate that this prior is robust to astrophysical assumptions and cosmic variance. Finally, applying this prior to the DESI DR1 dataset, we derive updated constraints on local PNG, obtaining f NL = -3.3±9.2.

cosmological parameters from LSS

A Radiometric Consistent Spectral Fingerprinting Algorithm for Continuity Products of Hyperspectral Sounders

A radiometric consistent climate fingerprinting methodology has been developed to derive long-term temperature, water vapor, cloud, trace gases, and surface skin temperature anomaly time series from the hyper-spectral sounder measurements of multiple platforms. The spectral fingerprinting methodology requires the use of radiative kernels that are radiometrically consistent with observations. Radiative kernels are built using space-time averaged Jacobians that are physically retrieved from observations under all sky conditions. The physical retrieval algorithm uses the Principal Component based Radiative Transfer Model (PCRTM) for the forward simulation. The incorporation of multiple scattering simulation in PCRTM allows the direct radiative relationship between single field-of-view (FOV) radiance observations and corresponding thermal dynamic variables including cloud properties to be established. Therefore, radiance ?closure? can be achieved under all-sky conditions by the fingerprinting scheme. This methodology has been used to derive climate anomalies from the space-time averaged spectra of AIRS/AMSU and CrIS/ATMS. The use of a consistent fingerprinting scheme provides an effective mean of generating continuity product by merging observations from different platforms and therefore facilitating the long-term climate trend study.

Wan Wu

General Electric 20-ampere hour nickel-cadmium battery

The interaction, effect, and controllability of the performance parameters of the General Electric 20-ampere-hour, 24-cell nickel cadmium battery are investigated. The battery was cycled under simulated orbit conditions. The acquired data was analyzed and evaluated in terms of battery parameters and performance characteristics. Conclusions and tests results are presented along with recommendations for further study.

Kirsch, W. W.

Valorization of consolidated bioprocessing residues for bioplastics

This study demonstrates an organic solvent-free processing strategy to valorize consolidated bioprocessing (CBP) residues, from switchgrass and poplar biomass, into functional poly(butylene succinate) (PBS)-based biocomposites using high-shear homogenization (HSH). HSH transformed the switchgrass and poplar CBP residues (CBP-R) into fine, uniformly distributed particles and microfibers. The composites of PBS with homogenized switchgrass residues (H-CBP-R-SG) or homogenized poplar residues (H-CBP-R-P) at a 70/30 weight ratio exhibited improved processability and mechanical integrity, with the Young's modulus for the PBS/H-CBP-R-SG and PBS/H-CBP-R-P nearly doubling to 0.66 ± 0.07 GPa and 0.65 ± 0.04 GPa, respectively, compared to neat PBS (0.36 ± 0.02 GPa). Dynamic Mechanical Analysis (DMA) reveals a significant suppression of the tan δ peak magnitude, indicating that HSH-mediated physical activation facilitates stress transfer in composites typical of covalent chemical grafting systems. While the transition to a stiffness-dominated profile reduces ductility, the resulting composites exhibit the dimensional stability and resistance to thermal warping required for high-fidelity FDM 3D printing and injection molding. Beyond material performance, comprehensive techno-economic analysis (TEA) and life cycle assessment (LCA) confirmed that diverting CBP residues into composite products can improve the economic viability of the biorefinery without substantially increasing biorefinery global warming potential (GWP). At a 30 wt% blend ratio, incorporating residuals into PBS yielded a minimum selling price for the composite of $\$4.07$ per kg compared to the conventional bioplastic price of $\$5.00$ per kg. This approach aligns with circular bioeconomy principles by converting waste streams into value-added products. Furthermore, this innovative strategy addresses key challenges in bioplastic development, including cost, compatibility, and performance, while simultaneously advancing waste minimization strategies for sustainable manufacturing systems.

09 BIOMASS FUELS

Socio-Techno-Economic Feasibility of Deep Geothermal with EGS for Residential District Heating in the Northeastern United States

In this paper, we present socio-techno-economic simulation results for enhanced geothermal system (EGS) district heating systems in the northeastern United States. We applied the geospatial Distributed Geothermal Market Demand (dGeo) tool coupled with GEOPHIRES to evaluate the feasibility of geothermal deep direct-use, utilizing EGS reservoirs for providing residential heating with district heating systems. The results were combined with JEDI results to assess social impacts. Simulation results assuming doublet EGS reservoirs with target production temperature of 80 degrees C and newly constructed district heating systems indicate a wide range in system levelized cost of heat (LCOH) values ranging from ~$15/MMBtu to over $1,000/MMBtu, with most the attractive regions urban areas with high total thermal demand and high thermal demand density. We estimate about 60 GWth in potential installed capacity with LCOH values under $50/MMBtu and 100 GWth in potential capacity with LCOH values under $100/MMBtu, respectively.

15 GEOTHERMAL ENERGY

Ejecta Generation and Redistribution on 433 Eros: Modeling Ejecta Launch Conditions

The NEAR-Shoemaker mission to asteroid 433 Eros presents an unprecedented opportunity to gain fundamental new knowledge about the processes governing regolith formation and redistribution on small bodies. NEAR-Shoemaker’s high-resolution imaging of the surface of Eros makes the asteroid a valuable and heretofore unparalleled laboratory for the detailed study of impact ejecta reaccretion and regolith redistribution on low-gravity (of order 10 -3 g) objects. Regolith is produced on asteroids by impact cratering, and the existence of regolith on the smallest solar system bodies supports the view that some of the ejecta from impact events on such objects may be retained. Impact craters and retained ejecta on low-gravity objects like Eros represent valuable natural laboratories for evaluating various models of impact cratering processes, since they may present crater structures or ejecta features that either do not form or are hidden on higher-gravity bodies like the Moon. Further, quantifying the extent to which impact processes generate and redistribute regoliths on small body surfaces (excavation depths, retained fraction, turnover timescales, etc.) is pivotal to the issue of how to relate meteoritical samples to their asteroidal parent bodies when surficial processes ( i.e., “space weathering”) may disguise or cover up underlying material and confound the ability of remote sensing techniques to provide reliable mineralogical assays of the parent objects. The rich variety of data on Eros’ regolith properties and distribution returned by NEAR-Shoemaker now require detailed analysis in order to take full advantage of the clues these observations offer for elucidating details of the impact cratering process on small bodies. Complicating simple interpretations of crater and ejecta morphology are dynamical effects on ejecta emplacement resulting from Eros’ irregular shape, rapid (5.27 hr) rotation, and low gravity. Figure 1 shows the very different ejecta deposit morphology that can result if the effects of rotation alone are neglected. Considering the additional complicating factors of Eros’ irregular shape and complex gravitational field, simple calculations of the extent and thickness of ejecta blankets and the spatial distribution of ejecta blocks from basic crater scaling laws or numerical hydrocodes alone do not suffice. In order to fully interpret the suite of NEAR-Shoemaker observations of regolith features across the surface of Eros and to evaluate various impact models for specific craters on the asteroid, detailed dynamical modeling of the deposition of crater ejecta from those craters is required . Here, I describe some modifications and improvements to the dynamical model being used for these studies.

D D Durda

A Probabilistic Approach to Load Modeling for Central HVAC Systems in Large Commercial Buildings for Retrofit Decisions Under Uncertainty

Retrofitting central HVAC systems in large commercial buildings with advanced technologies like heat recovery chillers (HRCs) offers a significant opportunity to enhance energy efficiency. However, analyzing these retrofits is challenging with traditional whole-building simulation tools, which require intensive calibration and struggle to model innovative system configurations and controls. To overcome these limitations, this study proposes a load profilebased retrofit analysis framework that provides better decisions under uncertainty. The main focus of this paper is the development of a probabilistic load profile model that can be used in the framework by using exploratory data analysis (EDA) of measured building data to properly quantify its inherent variability. A non-parametric Gaussian Process (GP) model was employed to capture the time- and weather-dependent characteristics of the heating load while explicitly modeling its uncertainty. The model's effectiveness is demonstrated through strong predictive performance on unseen data and physically interpretable insights into load behavior. This data-driven, probabilistic load profile serves as a robust and flexible input for subsequent system simulations, enabling a more confident and statistically sound analysis of retrofit potential.

Ham, S W

From Exploration Flight Test-1 to Artemis II--A NASA Langley's Orion Aerosciences Overview

The Orion Aerosciences program at NASA Langley has played a central role in advancing the aerodynamic and aeroheating prediction capabilities required for the Orion crew vehicle’s return from deep space. This presentation provides a technical overview of aerosciences contributions spanning Exploration Flight Test-1 (EFT-1), Artemis I, and the ongoing post-flight analysis of Artemis II. EFT-1 provided the first high-energy entry dataset for Orion, enabling critical validation of aerodynamic force and moment predictions, static and dynamic stability characteristics, and aeroheating environments at relevant flight Mach and Reynolds numbers. Flight-derived pressure data were used to refine the Flush Air Data System (FADS) methodology for atmospheric density reconstruction and to improve Best Estimated Trajectory (BET) solutions. The EFT-1 data also offered key insights into heat shield performance, including char layer recession, in-depth thermal response, and material retention behavior under flight conditions, informing updates to both thermal response models and uncertainty quantification practices. Building on EFT-1, Artemis I extended the database to true lunar-return conditions. Observations of heat shield performance, including localized char loss, bondline response, and recession variability, provided an unprecedented opportunity to reassess Thermal Protection System (TPS) and aeroheating modeling assumptions. Aerodynamic reconstruction efforts incorporated improved FADS calibration, enhanced atmospheric modeling, and refined force and moment databases to reduce trajectory and load uncertainties. Aeroheating comparisons between pre-flight predictions and flight data enabled targeted model updates, particularly in transitional flow environments and wake heating regions. For Artemis II, these lessons were systematically incorporated into the pre-flight prediction process. Updates included refined aerodynamic databases anchored to flight-validated corrections, improved density estimation and BET methodologies using enhanced database interpolation algorithm and FADS modeling, and revised aeroheating design environments informed by Artemis I material response observations. By the time of the workshop, Artemis II post-flight analysis will be underway, and preliminary findings will be presented where available, including early comparisons of aerodynamic reconstruction, atmospheric density estimation, and thermal protection system performance relative to updated predictions. Collectively, this body of work is a testament to the dedicated and multidisciplinary team whose sustained efforts have contributed to the program’s success and to the progressive maturation of Orion aerosciences modeling through numerical modeling, ground and flight data assimilation. The integrated advancement of aerodynamics, trajectory reconstruction, FADS-based density estimation, and aeroheating analysis has reduced predictive uncertainty and strengthened confidence for future crewed lunar and deep-space missions.

Orion

Superdiffusion resilience in Heisenberg chains with two-dimensional interactions on a quantum processor

Superdiffusive spin transport in the one-dimensional (1D) Heisenberg model is a key theoretical discovery in nonequilibrium quantum many-body physics. Although extensively studied in 1D systems, the breakdown and sustenance of superdiffusion in two-dimensional (2D) lattices with integrability-breaking terms, as found in real materials, remains an open question. To address this, we develop a toy model that extends the 1D Heisenberg model with a representative set of 2D interaction types and tunable strengths. Our model exhibits varying degrees of superdiffusion breakdown depending on the interaction type, spanning ballistic to diffusive regimes. We establish and justify a hierarchy of 2D interactions based on their resilience against superdiffusion breakdown: Heisenberg >𝑋⁢𝑋 > Ising. This precise control over the superdiffusive behavior also enables rigorous benchmarking of quantum hardware, and our simulations on IBM's Heron devices confirm the hardware's ability to accurately capture these many-body nonequilibrium phenomena. Overall, our results are relevant not only to simulating superdiffusion in real materials, such as the 1D Heisenberg compound KCuF3, which contains modest nonintegrable 2D terms, but also to extending superdiffusive behavior to larger 2D qubit lattices and other 2D materials.

Alagarsamy Manikandan, Keerthi Kumaran [ORNL]

Battery model for electrical power system energy balance

A model to simulate nickel-cadmium battery performance and response in a spacecraft electrical power system energy balance calculation was developed. The voltage of the battery is given as a function of temperature, operating depth-of-charge (DOD), and battery state-of-charge. Also accounted for is charge inefficiency. A battery is modeled by analysis of the results of a multiparameter battery cycling test at various temperatures and DOD's.

Hafen, D. P.

Temperature as a control knob on spin-orbit coupling

Spin-orbit coupling (SOC) governs many physical phenomena. Through ab initio molecular dynamics simulations, Lu and Sun demonstrated that structural disorder at elevated temperatures substantially reduces the effective SOC contribution that stabilizes band inversion, driving a topological-to-normal-insulator transition in Bi2Se3. Their work identifies temperature as a meaningful control parameter for SOC-mediated topology and other properties.

Liang, Liangbo [ORNL] (ORCID:0000000311990049)

NASA battery testbed: A designed experiment for the optimization of LEO battery operational parameters

Simulation of spacecraft battery operation is implemented. Robust design experiment to obtain optimum battery operational parameters is performed. It is found that short term tests using robust design of experiments can provide guidelines for optimum battery operation. It is decided to use robust design approach to provide guidelines for battery operation on current spacecraft in orbit as batteries age (GRO, UARS, EUVE, TOPEX).

Deligiannis, F.

Studies of Flame Structure in Microgravity

The present research endeavor is concerned with gaining fundamental understanding of the configuration, structure, and dynamics of laminar premixed and diffusion flames under conditions of negligible effects of gravity. Of particular interest is the potential to establish and hence study the properties of spherically- and cylindrically-symmetric flames and their response to external forces not related to gravity. For example, in an earlier experimental study of the burner-stabilized cylindrical premixed flames, the possibility of flame stabilization through flow divergence was established, while the resulting one-dimensional, adiabatic, stretchless flame also allowed an accurate means of determining the laminar flame speeds of combustible mixtures. We have recently extended our studies of the flame structure in microgravity along the following directions: (1) Analysis of the dynamics of spherical premixed flames; (2) Analysis of the spreading of cylindrical diffusion flames; (3) Experimental observation of an interesting dual luminous zone structure of a steady-state, microbuoyancy, spherical diffusion flame of air burning in a hydrogen/methane mixture environment, and its subsequent quantification through computational simulation with detailed chemistry and transport; (4) Experimental quantification of the unsteady growth of a spherical diffusion flame; and (5) Computational simulation of stretched, diffusionally-imbalanced premixed flames near and beyond the conventional limits of flammability, and the substantiation of the concept of extended limits of flammability. Motivation and results of these investigations are individually discussed.

C K Law