Search NASASearch

SEARCH · Search NASA

Results for “simulation software”

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.

617 records · Page 16

Vidyut3d: A Gpu Accelerated Fluid Solver for Non-Equilibrium Plasmas on Adaptive Grids

We present the numerical methods, programming methodology, verification, and performance assessment of a non-equilibrium plasma fluid solver that can effectively utilize current and upcoming central processing and graphics processing unit (CPU+GPU) architectures, in this work. Our plasma fluid model solves the coupled conservation equations for species transport, electrostatic Poisson and electron temperature on adaptive Cartesian grids. Our solver is written using performance portable adaptive-grid/particle management library, AMReX, and is portable over widely available vendor specific GPU architectures. We present verification of our solver using method of manufactured solutions that indicate formal second order accuracy with central diffusion and fifth-order weighted-essentially-non-oscillatory (WENO) advection scheme. We also verify our solver with published literature on capacitive discharges and atmospheric pressure streamer propagation. We demonstrate the use of our solver on two 3D simulation cases: an atmospheric streamer propagation in Ar-H2 mixtures and a low pressure twin electrode radio frequency reactor. Our performance studies on three different CPU+GPU architectures indicate approximately 150-400X speed-up using AMD and NVIDIA GPUs per time step compared to a single CPU core for a 4 million cell simulation with 15 species.

Sitaraman, Hariswaran

Progress on Shape Memory Alloy Actuator Development for Active Clearance Control

Results of a numerical analysis evaluating the feasibility of high-temperature shape memory alloys (HTSMA) for active clearance control actuation in the high-pressure turbine section of a modern turbofan engine has been conducted. The prototype actuator concept considered here consists of parallel HTSMA wires attached to the shroud that is located on the exterior of the turbine case. A transient model of an HTSMA actuator was used to evaluate active clearance control at various operating points in a test bed aircraft engine simulation. For the engine under consideration, each actuator must be designed to counteract loads from 380 to 2000 lbf and displace at least 0.033 in. Design results show that an actuator comprised of 10 wires 2 in. in length is adequate for control at critical engine operating points and still exhibit acceptable failsafe operability and cycle life. A proportional-integral- derivative (PID) controller with integrator windup protection was implemented to control clearance amidst engine transients during a normal mission. Simulation results show that the control system exhibits minimal variability in clearance control performance across the operating envelope. The final actuator design is sufficiently small to fit within the limited space outside the high-pressure turbine case and is shown to consume only small amounts of bleed air to adequately regulate temperature.

Jonathan DeCastro

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]

Material Properties and Modeling of Room Temperature Vulcanizing Silicone

Room Temperature Vulcanizing silicone (RTV) is a high-temperature adhesive that has successfully been used as a gap-filler between Thermal Protection System (TPS) tiles for heatshields on numerous missions. It is also used to bond instrumentation plugs such as temperature and pressure sensors into the heatshields. While RTV has been traditionally assumed to be a non-porous and non-ablating material, numerous experiments have shown that RTV pyrolyzes and becomes highly porous as it is heated. Heating RTV has also shown swelling, or intumescence, which can pose unique problems that lead to roughness induced boundary-layer transition, surface oxide formation and contamination of heat shield sensors. Therefore, it is crucial to understand and model the intumescence phenomenon of RTV. As data for RTV material properties is limited, the first step in modeling RTV is to collect material properties such as pyrolysis mass-loss, microstructure change, virgin and char porosity, etc. which was performed in our initial study. Additionally, thermomechanical properties such as Young’s modulus and Poisson ratio are required for modeling the intumescence of RTV, which were taken from literature and the coefficient of thermal expansion was collected using in-situ heating and Micro Computed Tomography (µ-CT) in previous studies. Finally, numerous other properties such as pyrolysis gas properties, virgin and char thermal conductivity and specific heat were compiled from previous experiments and literature into a material database that can be used for simulations. In Porous Material Analysis Toolbox based on OpenFOAM (PATO) [4], structural mechanics coupled with material response was used for simulating the intumescence of RTV as it is heated. However, since the permeability of the material is very low, the pyrolysis gas creates an internal pressure build-up as the material is being heated, significantly contributing to the deformation of the material. To correctly characterize this phenomenon, additional physics models were implemented into PATO's stress analysis solver, and results were compared with RTV dilatometry test data as a preliminary verification case. Future work will include experiments of RTV at the Plasmatron X facility and the in-situ heating cell with µ-CT, and improvement of simulation tools to more accurately model RTV intumescence.

PATO

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

Numerical Modeling of the Transient Chilldown Process of a Cryogenic Propellant Transfer Line

Before cryogenic fuel depots can be fully realized, efficient methods with which to chill down the spacecraft transfer line and receiver tank are required. This paper presents numerical modeling of the chilldown of a liquid hydrogen tank-to-tank propellant transfer line using the Generalized Fluid System Simulation Program (GFSSP). To compare with data from recently concluded turbulent LH2 chill down experiments, seven different cases were run across a range of inlet liquid temperatures and mass flow rates. Both trickle and pulse chill down methods were simulated. The GFSSP model qualitatively matches external skin mounted temperature readings, but large differences are shown between measured and predicted internal stream temperatures. Discrepancies are attributed to the simplified model correlation used to compute two-phase flow boiling heat transfer. Flow visualization from testing shows that the initial bottoming out of skin mounted sensors corresponds to annular flow, but that considerable time is required for the stream sensor to achieve steady state as the system moves through annular, churn, and bubbly flow. The GFSSP model does adequately well in tracking trends in the data but further work is needed to refine the two-phase flow modeling to better match observed test data.

Temperature Sensors

Development and Experimental Validation of a Path-Dependent Spin Forming Finite Element Model

Spin forming is an advanced manufacturing process widely used in the aerospace and defense sectors to produce lightweight, high-strength cylindrical components with tight dimensional tolerances. This study explores the applicability of the path-dependent Mechanical Threshold Stress (MTS) constitutive model by simulating the evolution of geometry, machining forces, and plastic deformation during the spin forming of a 10-mm thick 6061-O aluminum cylinder. While numerical modeling of spin forming has advanced substantially over the past decade, systematic verification and experimental validation of material models remain limited, particularly in predicting through-thickness process evolution. The MTS model, incorporating a Voce hardening rule, is employed for its ability to represent cyclic loading, rapidly varying temperature fields, and strain rates characteristic of spin forming. Numerical convergence analysis indicates discretization uncertainties between 0.3% and 9.2% for key quantities of interest. Experimental validation demonstrates that the MTS model, when implemented with a verified mesh, accurately reproduces both elastic and plastic behavior of 6061-O aluminum, predicting peak roller loads within 11–18% of measurements, geometric tolerances within 3%, and plastic strain distributions within 10% of experimental values. Collectively, these results establish a validated computational framework for predictive spin-forming simulations with quantified confidence, providing a foundation for extension to other alloys, geometries, and forming conditions.

Spin forming

Synthetic battery cycling techniques

The group of techniques that as a class are referred to as synthetic battery cycling are described with reference to spacecraft battery systems. Synthetic battery cycling makes use of the capability of computer graphics to illustrate some of the basic characteristics of operation of individual electrodes within an operating electrochemical cell. It can also simulate the operation of an entire string of cells that are used as the energy storage subsystem of a power system.

Leibecki, H.

Modeling Ni-Cd performance. Planned alterations to the Goddard battery model

The Goddard Space Flight Center (GSFC) currently has a preliminary computer model to simulate a Nickel Cadmium (Ni-Cd) performance. The basic methodology of the model was described in the paper entitled Fundamental Algorithms of the Goddard Battery Model. At present, the model is undergoing alterations to increase its efficiency, accuracy, and generality. A review of the present battery model is given, and the planned charges of the model are described.

Jagielski, J. M.

Hypervelocity Impact Testing of Nickel Hydrogen Battery Cells

Nickel-Hydrogen (Ni/H2) battery cells have been used on several satellites and are planned for use on the International Space Station. In January 1992, the NASA Lewis Research Center (LeRC) conducted hypervelocity impact testing on Ni/H2 cells to characterize their failure modes. The cell's outer construction was a 24 mil-thick Inconel 718 pressure vessel. A sheet of 1.27 cm thick honeycomb was placed in front of the battery cells during testing to simulate the on-orbit box enclosure. Testing was conducted at the NASA White Sands Test Facility (WSTF). The hypervelocity gun used was a 7.6 mm (0.30 caliber) two-stage light gas gun. Test were performed at speeds of 3, 6, and 7 km/sec using aluminum 2017 spherical particles of either 4.8 or 6.4 mm diameter as the projectile. The battery cells were electrically charged to about 75 percent of capacity, then back-filled with hydrogen gas to 900 psi simulating the full charge condition. High speed film at 10,000 frames/sec was taken of the impacts. Impacts in the dome area (top) and the electrode area (middle) of the battery cells were investigated. Five tests on battery cells were performed. The results revealed that in all of the test conditions investigated, the battery cells simply vented their hydrogen gas and some electrolyte, but did not burst or generate any large debris fragments.

Frate, David T.

Thermally induced mimicry of quantum cluster excitations and implications for the magnetic transition in FePSe 3

In two dimensional magnets, the interplay of thermal fluctuations and spin anisotropy control the existence of long-range magnetic order. In the van der Waals antiferromagnets FePX 3 , orbital degeneracy in the 𝑡 2⁢𝑔 levels of the Fe 2+ ions in octahedral coordination yields strong uniaxial anisotropy, which stabilizes magnetic order up to 𝑇 ≈ 100 K. Recent inelastic neutron scattering measurements around the magnetic ordering transition have shown the existence of a broad spectrum of magnetic fluctuations with nontrivial momentum dependence, which has been interpreted as evidence for localized entangled cluster excitations. In this paper, we offer an alternative interpretation using classical nonlinear spin dynamics simulations. We present stochastic Landau Lifshitz dynamics simulations that reproduce the neutron scattering measurements of Chen et al. [npj Quantum Mater. 9, 40 (2024)] on FePSe 3 . These calculations faithfully explain the dynamical structure factor's momentum and energy dependence and point to a classical origin for the excitations observed in neutron spectroscopy and that the order-disorder transition can be understood in terms of thermal fluctuations overcoming the anisotropy energy.

Landau-Lifschitz-Gilbert equation

Transforming Windows from Energy Liabilities to Zero-Energy Assets: Next-Generation Solutions for Buildings

Windows have traditionally contributed to a building's HVAC load, but they can also become a source of net energy gain or even operate as zero-energy components. For heating applications, highly insulating windows can harness more solar heat than the energy lost through them, transforming windows from energy liabilities to assets. Dynamic glazings provide further benefits by regulating solar heat gain, reducing cooling loads in summer and heating demands in winter. This simulation study focuses on developing the next generation of zero-energy windows (ZEW) for residential new construction. Through annual energy simulations across climate zones 1-8, ZEW performance benchmarks were established based on current code-level buildings, and we've identified the regions where meeting ZEW standards are most achievable. This work evaluates both static and dynamic window technologies, assessing their effects on annual energy use and cost. Key findings demonstrate that ZEW performance is achievable across diverse climate zones, with specific regional requirements. Most climate zones from 3-8 can achieve ZEW with specific configurations, while some warm climates (1-2) appear challenging for ZEW implementation. Climate zones 4-6 consistently allow for zero energy window implementation, offering multiple pathways through either static or dynamic window technologies. Colder climate zones (7-8) ZEW products allow for higher SHGC values while requiring low U-values.

Yu, Lili

Decoupling thermal and irradiation effects on grain boundary segregation

Radiation-induced segregation (RIS) is most often measured by peak solute concentration at a boundary. However, this may give an incomplete picture of segregation quantity and phenomena. Radiation-induced and thermal segregation at grain boundaries was investigated in Fe-9.6 at.% Cr after 9 MeV Fe 3+ ion irradiation at 400 °C. The experimental results were compared to kinetic Monte Carlo (KMC) simulations. The study revealed that Cr enrichment (peak segregation) at the grain boundaries was comparable in both the irradiated and non-irradiated conditions, although irradiation resulted in broader segregation profiles in both experiments and simulations, indicating an overall increase in grain boundary segregation due to irradiation. This broadening is attributed to back diffusion into the grain interior. While it is an established phenomenon, this study offers a quantitative evaluation using experimental data and KMC modeling. Further, these results emphasize the importance of analyzing the entire segregation profile and decoupling the thermal and irradiation contributions to solute segregation.

Grain boundary segregation

Predicting Atomistic Transitions with Transformers

Accurate knowledge of the atomistic transition pathways in materials and material surfaces is crucial for many material science problems. However, conventional simulation techniques used to find these transitions are extremely computationally intensive. Even with large-scale, accelerated material simulations, the computational cost constrains the applicable domain in practice. Machine learning models, with the potential to learn the complex emergent behaviors governing atomistic transitions as a fast surrogate model, have great promise to predict transitions with a vastly reduced computational cost. Here, we demonstrate how transformers can be trained to predict atomistic transitions in nano-clusters. We show how we evaluate physical validity of the predictions and how a multitude of additional, different microstates can be generated by slightly varying the data provided to the model.

36 MATERIALS SCIENCE