Search NASA⌕ Search

SEARCH · Search NASA

Results for “finite element 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.

739 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

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↗
DOE OSTI2025

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↗
NASA NTRS2006

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↗
DOE OSTI2026

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]↗
NASA NTRSDate not supplied

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↗
NASA NTRS1983

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.↗
DOE OSTI2026

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 NTRS1996

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.↗
NASA NTRS1997

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↗
NASA NTRS2015

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↗
NASA NTRS1982

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.↗
NASA NTRS1986

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.↗
NASA NTRS1996

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.↗
DOE OSTI2025

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↗
DOE OSTI2025

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↗
NASA NTRSDate not supplied

Experimental Characterization of Additively Manufactured Nickel-Titanium Shape Memory Alloy Heat Pipes

Shape memory alloys (SMA) have been identified for use in spacecraft components as replacement for conventional deployment mechanisms. They may be used in thermal management components such as radiators to create self-deploying radiators. One SMA, NiTi, has also been developed for additive manufacturing processes. Heat pipes are a common way to create highly effective and lightweight spaceflight radiators, and heat pipes can also be made from NiTi and related alloys. The wick is the critical element of a functioning heat pipe, and recent progress over the past years has led to the development of additively manufactured heat pipe wicks in various materials. The combination of these efforts is the focus of this project: creating an additively manufactured, shape memory alloy self-deploying heat pipe radiator. This paper will focus on the experimental characterization of these additively manufactured NiTi heat pipes. The heat pipe coupons were additively manufactured by direct metal laser sintering (DMLS), with an integral liquid cooled condenser. Heat is input to the heat pipe via a thin film heater. Thermocouples were spot welded to the heat pipes to measure temperature at several axial locations. The heat pipes were tested with two working fluids: water and ethanol. Ethanol is not an ideal working fluid for heat pipes but is useful in characterizing them because it wets well to a wide variety of surfaces. Water is in general a superior working fluid for heat pipes, but its contact angle and therefore wicking performance strongly depends on the surface chemistry of the surface it is in contact with. A particular measurement of interest in this test is the evaporator to condenser thermal conductance, which will be compared in the full paper to recently published correlations for additively manufactured heat pipes. Experimental results for two straight geometry and one bellows geometry heat pipe will be presented. The bellows geometry is of interest for condenser of the self-deploying radiator design.

Additive manufacturing↗
DOE OSTI2026

The ABCs of phase retrieval: Connecting the acronyms of scanning transmission electron microscopy

High-resolution scanning transmission electron microscopy (S/TEM) is an indispensable tool for characterizing the structure and properties of materials down to the atomic scale. Conventional S/TEM imaging, however, is limited by the phase problem, whereby the phase of the electron exit wave is lost upon detection. Recent advances in diffractive imaging and 4D-STEM have enabled a range of phase-retrieval techniques that computationally reconstruct the missing information encoded in the phase of the transmission function. These approaches offer improved dose efficiency and enhanced sensitivity to weakly scattering signals, extending quantitative imaging to beam-sensitive materials composed of light elements. In this work, we introduce the phase problem in electron microscopy and survey the diverse landscape of phase-retrieval techniques used in the field. Despite their many acronyms and algorithmic variations, these techniques share a common physical and mathematical foundation. We present a unified framework that connects these seemingly distinct methods, from parallax imaging and tilt-corrected bright-field (tcBF-STEM), to aberration-corrected bright-field (acBF-STEM), optimum bright-field (OBF-STEM) and single-sideband (SSB) ptychography, as well as first-moment integrated center of mass techniques (iCOM) and iterative ptychographic algorithms. Based on these insights, we discuss the opportunities and practical limitations of applying these methods across different materials systems, detector designs, and microscope configurations.Graphical abstractRepresentative electron microscopy configurations used for phase retrieval and diffractive imaging in S/TEM: (a) Zernike phase-contrast transmission electron microscopy (TEM), (b) small-convergence-angle four-dimensional scanning transmission electron microscopy (4D-STEM) for nanobeam-based phase reconstruction methods, and (c) large-convergence-angle 4D-STEM for ptychographic and related diffractive imaging techniques reviewed in this work.

36 MATERIALS SCIENCE↗