SEARCH · Search NASA

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

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 415 records · Page 23

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Thermodynamic Phase Transition in Magnetic Reconnection

By examining the entropy production in fully kinetic simulations of collisional plasmas, it is shown that the transition from collisional Sweet-Parker reconnection to collisionless Hall reconnection may be viewed as a thermodynamic phase transition. The phase transition occurs when the reconnection electric field satisfies E = $E_{D}\sqrt{m{e}/{m}_i}$ , where ${m{e}/{m}_i}$ is the electron-to-ion mass ratio and E D is the Dreicer electric field. This condition applies for all ${m{e}/{m}_i}$, including ${m{e}/{m}_i}$ = 1 , where the Hall regime vanishes and a direct phase transition from the collisional to the kinetic regime occurs. In the limit ${m{e}/{m}_i}$ → 0, this condition is equivalent to there being a critical electron temperature $Te ≈ m_iΩ^{2}_{i}δ^{2}$, where $Ω_i$ is the ion cyclotron frequency and δ is the current sheet half-thickness. The heat capacity of the current sheet changes discontinuously across the phase transition, and a critical power law is identified in an effective heat capacity. A model for the time-dependent evolution of an isolated current sheet in the collisional regime is derived.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

A high-order, localized-artificial-diffusivity method for Eulerian simulation of multi-material elastic-plastic deformation with strain hardening

A high-order method for Eulerian simulation of material undergoing large elastic–plastic deformation is developed. Thermodynamically consistent hyperelastic constitutive relations are assumed, facilitating the treatment of solids, liquids, and gases in a unified manner. Here, the method enables the simulation of multi-material interactions using a diffuse interface approach. Numerical capturing of material interfaces, shock waves, contact surfaces, and elastic-plastic strain discontinuities using high-order compact-difference schemes is assisted by Localized Artificial Diffusivity (LAD). In the new setting involving elastic–plastic deformation, the previously established terms for the artificial properties are verified to effectively regularize normal shocks. Additional LAD terms are introduced to the elastic and plastic kinematic equations to regularize shear shocks and other strain discontinuities, improving solution stability. Other important features of the method that improve robustness include the numerical treatment of compatibility terms in the kinematic equations, and the treatment of rotation. Particular emphasis is focused toward new advancements of the methods for plastic-deformation integration and the associated strain hardening of the material, including rate-dependent plasticity. The method is demonstrated on a variety of test problems, including 1-D impacts, a variant of the Shu-Osher problem, a Taylor impact, and a Richtmyer-Meshkov instability between two elastic–plastic solids with strain hardening.

42 ENGINEERING↗

DOE OSTI2021

Predicting large-scale pool fire dynamics using an unsteady flamelet- and large-eddy simulation-based model suite

A low-Mach, unstructured, large-eddy-simulation-based, unsteady flamelet approach with a generalized heat loss combustion methodology (including soot generation and consumption mechanisms) is deployed to support a large-scale, quiescent, 5-m JP-8 pool fire validation study. The quiescent pool fire validation study deploys solution sensitivity procedures, i.e., the effect of mesh and time step refinement on capturing key fire dynamics such as fingering and puffing, as mesh resolutions approach O(1) cm. A novel design-order, discrete-ordinate-method discretization methodology is established by use of an analytical thermal/participating media radiation solution on both low-order hexahedral and tetrahedral mesh topologies in addition to quadratic hexahedral elements. Here, the coupling between heat losses and the flamelet thermochemical state is achieved by augmenting the unsteady flamelet equation set with a heat loss source term. Soot and radiation source terms are determined using flamelet approaches for the full range of heat losses experienced in fire applications including radiative extinction. The proposed modeling and simulation paradigm are validated using pool surface radiative heat flux, maximum centerline temperature location, and puffing frequency data, all of which are predicted within 10% accuracy. Simulations demonstrate that under-resolved meshes predict an overly conservative radiative heat flux magnitude with improved comparisons as compared to a previously deployed hybrid Reynolds-averaged Navier–Stokes/eddy dissipation concept-based methodology.

42 ENGINEERING↗

DOE OSTI2023

Ab initio Geochemistry of Hydrous Phases (Final Report)

Hydrous phases are among the most important Earth components. They are important for a broad suite of Earth processes, including the origin of life. From both thermodynamic and structural perspectives, however, they represent some of the most complex naturally occurring materials: their bonding often includes a combination of covalent solid, weak ionic, van der Waals, and hydrogen bonding, all within large unit cells. Most are solid solutions, and many are prone to variations in layer packing. Thermodynamic modeling of these materials is fundamental for understanding present and past natural processes, including those that shaped—and continue to shape—the structure and evolution of our planet. Yet, the thermodynamic properties of these materials at appropriate conditions are challenging to measure. To make significant progress and attain a deep understanding of these materials requires an atomistic theoretical approach.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Improvement of Atmospheric Objective Analysis Over Sloping Terrain and Its Impact on Shallow-Cumulus Clouds in Large-Eddy Simulations

Surface topography strongly impacts the regional atmospheric circulation dynamically and thermodynamically. Over the Great Plains in the United States, the gently tilted slope is an important factor that impacts clouds, convection and regional circulations. This study enhances an atmospheric constrained variational analysis by using a terrain-following sigma vertical coordinate and applies to the data collected at the Southern Great Plain (SGP) site of the Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) program. Sensitivity studies are performed to examine the impact of the terrain effects on the derived large-scale atmospheric forcing fields and the simulated shallow-cumulus clouds in large-eddy simulation (LES) driven by the forcing. We found that the terrain impacts on the large-scale forcing fields are mainly at lower levels and are strongly controlled by up/downslope winds. The response of the derived forcing to the slope of the terrain is monotonic, but the response of the simulated shallow-cumulus clouds is more complex and depends on several factors. Overall, the terrain impact is small over SGP due to the small slope angle. However, the flat-surface assumption may cause larger biases in the large-scale forcing fields at locations with steeper terrain. The new sigma-coordinate algorithm, with its consideration of surface slope, should be more suitable to derive large-scale objective analysis over regions with steep terrain for application to force single-column models, cloud resolving models and LES models.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Agreement of Analytical and Simulation-Based Estimates of the Required Land Depth in Climate Models

Previous analytical and simulation-based analyses suggest that deeper land surface models are needed to realistically simulate the terrestrial thermal state in climate models, with implications for land-atmosphere interactions. Analytical approaches mainly focused on the subsurface propagation of harmonics such as the annual temperature signal, and a direct comparison with climate-change model output has been elusive. Here, this study addresses the propagation of a harmonic pulse fitted to represent the timescale and amplitude of anthropogenic warming. Its comparison to land model simulations with stepwise increased bottom boundary depth leads to an agreement between the simulation-based and analytical frameworks for long-term climate trends. Any depth increase gradually decreases the relative error in the subsurface thermodynamics, and a minimum depth of 170 m is recommended to simulate the ground climate adequately. The approach provides an accurate estimate of the required land-model depth for climate-change simulations and assesses the relative bias in insufficiently deep land models.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRSDate not supplied

Validation of Ullage Collapse Due to Violent Lateral Slosh

Understanding the coupling between thermodynamics and propellant sloshing is important for cryogenic upper stages, especially during the lift-off phase and in missions including multiple engine restarts, as the fluid dynamic condition of the propellant is mainly influenced by various flight maneuvers. This paper describes the development and validation of a Volume of Fluid (VOF)-based CFD solver, Loci/STREAM-VOF, to the ullage collapse under violent later slosh conditions. Loci/STREAM-VOF is shown to be able to capture detailed slosh dynamics, including slosh frequency, non-linear slosh damping, and violent surface breakup due to forced acceleration. The comparison of surface shape and small drops formed during large amplitude slosh with experimental imagery shows excellent agreement throughout the sloshing period. Comparison of temperature histories in the ullage, at the gas-liquid interface, and inside liquid with the experiment conducted at JAXA also show very good agreement. The simulation points to a stratified temperature field inside the liquid near the interface, critical to the ullage pressure after the collapse. The predicted pressure collapse matches well with the experimental value due to sloshing. Parametric sensitivity studies suggest that high-fidelity Loci/STREM-VOF simulations can play an essential role in predicting the actual thermodynamics of ullage collapse when following the recommended simulation methodology.

H Q Yang↗

DOE OSTI2024

A platform for non-local thermodynamic equilibrium atomic physics research using the buried layer target approach at a kJ level laser facility

Here we present a design and first use of a kJ level laser facility for research of non-local thermodynamic equilibrium atomic physics using the buried layer target method. The target design included a metal layer buried inside a plastic tamper with thicknesses tailored to the expected laser intensities. The target was illuminated from each side by two laser beams with intensities of 0.5–5 × 10 14 W/cm 2 . The advanced diagnostic suite included static and time-resolved imagers and spectrometers with various spectral resolutions. A 3D printed dual elliptically curved spectrometer is presented, and its results are compared to a traditional crystal spectrometer. Experimental results and radiation hydrodynamic simulations demonstrate that the target achieved the desired thermodynamic conditions of n e ≈ 10 21 –10 22 cm -3 and T e ≈ 1–2 keV.

3D printing↗

DOE OSTI2025

Thermodynamic Modeling of Complex Solid Solutions in the Lu-H-N System via Graph Neural Network Accelerated Monte Carlo Simulations

Metal hydrides are important across diverse applications, such as hydrogen storage, batteries, gas sensors, nuclear reactions, and high-temperature superconductivity. Previous computational studies of metal hydrides under extreme pressures, e.g., 𝑂⁡(10 2 ) ⁢GPa, usually treat them as stoichiometric compounds without considering interstitial lattice disorder. As pressures become more moderate in the 𝑂⁡(10 0 ) ⁢GPa and below range, hydrogen disorder at interstitial lattice sites becomes prominent, e.g., in the N-doped Lu hydride that was recently claimed superconducting near 1 GPa. Further adding compositional complexity from alloying and/or multielement interstitial occupation makes elucidating pressure- and temperature-dependent observables intractable by first-principles calculations alone. We therefore propose a lattice graph neural-network surrogate modeling approach to predict configuration- and pressure-dependent equation-of-state properties. Their efficiency permits Monte Carlo simulations to calculate Gibbs energies and pressure-dependent phase diagrams, thereby revealing insights into the synthesis conditions required for achieving desired phase equilibria. We demonstrate this concept for the compositionally complex cubic Lu(H,N,Va) 3 system where three constituents (hydrogen, nitrogen and vacancy) have disordered multielement interstitial occupancies and insights into pressure-dependent phase equilibria are critically needed, e.g., N-doping levels can significantly lower dehydrogenation temperatures and provide a new strategy to optimize hydrogen-storage alloys. This work can improve the thermodynamic understanding of the Lu-H-N system and help rational synthesis of N-doped Lu hydrides, but more generally demonstrates an efficient approach to model pressure-dependent thermodynamics of multicomponent solid solutions.

Monte Carlo methods↗

NASA NTRS1993

Technologies for Visualization in Computational Aerosciences

State-of-the-art research in computational aerosciences produces' complex, time-dependent datasets. Simulations can also be multidisciplinary in nature, coupling two or more physical disciplines such as fluid dynamics, structural dynamics, thermodynamics, and acoustics. Many diverse technologies are necessary for visualizing computational aerosciences simulations. This paper describes these technologies and how they contribute to building effective tools for use by domain scientists. These technologies include data management, distributed environments, advanced user interfaces, rapid prototyping environments, parallel computation, and methods to visualize the scalar and vector fields associated with computational aerosciences datasets.

Miceli, Kristina D.↗

NASA NTRS2004

Multi-d CFD Modeling of a Free-piston Stirling Convertor at NASA Glenn

A high efficiency Stirling Radioisotope Generator (SRG) is being developed for possible use in long duration space science missions. NASA s advanced technology goals for next generation Stirling convertors include increasing the Carnot efficiency and percent of Carnot efficiency. To help achieve these goals, a multidimensional Computational Fluid Dynamics (CFD) code is being developed to numerically model unsteady fluid flow and heat transfer phenomena of the oscillating working gas inside Stirling convertors. Simulations of the Stirling convertors for the SRG will help characterize the thermodynamic losses resulting from fluid flow and heat transfer between the working gas and solid walls. The current CFD simulation represents approximated 2-dimensional convertor geometry. The simulation solves the Navier Stokes equations for an ideal helium gas oscillating at low speeds. The current simulation results are discussed.

Wilson, Scott D.↗

DOE OSTI2023

Transfer Learning of High-Fidelity Opacity Spectra in Autoencoders and Surrogate Models

Simulations of high energy density physics are expensive, largely in part for the need to produce nonlocal thermodynamic equilibrium opacities. High-fidelity spectra may reveal new physics in the simulations not seen with low-fidelity spectra, but the cost of these simulations also scales with the level of fidelity of the opacities being used. Neural networks are capable of reproducing these spectra, but neural networks need data to train them, which limits the level of fidelity of the training data. Here this article demonstrates that it is possible to reproduce high-fidelity spectra with median errors in the realm of 3%–4% using as few as 50 samples of high-fidelity Krypton data by performing transfer learning on a neural network trained on many times more low-fidelity data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS2007

CFD Modeling of Helium Pressurant Effects on Cryogenic Tank Pressure Rise Rates in Normal Gravity

A recently developed computational fluid dynamics modeling capability for cryogenic tanks is used to simulate both self-pressurization from external heating and also depressurization from thermodynamic vent operation. Axisymmetric models using a modified version of the commercially available FLOW-3D software are used to simulate actual physical tests. The models assume an incompressible liquid phase with density that is a function of temperature only. A fully compressible formulation is used for the ullage gas mixture that contains both condensable vapor and a noncondensable gas component. The tests, conducted at the NASA Marshall Space Flight Center, include both liquid hydrogen and nitrogen in tanks with ullage gas mixtures of each liquid's vapor and helium. Pressure and temperature predictions from the model are compared to sensor measurements from the tests and a good agreement is achieved. This further establishes the accuracy of the developed FLOW-3D based modeling approach for cryogenic systems.

Grayson, Gary↗

DOE OSTI2020

Climatology of Severe Local Storm Environments and Synoptic-Scale Features over North America in ERA5 Reanalysis and CAM6 Simulation

Severe local storm (SLS) activity is known to occur within specific thermodynamic and kinematic environments. These environments are commonly associated with key synoptic-scale features—including southerly Great Plains low-level jets, drylines, elevated mixed layers, and extratropical cyclones—that link the large-scale climate to SLS environments. This work analyzes spatiotemporal distributions of both extreme values of SLS environmental parameters and synoptic-scale features in the ERA5 reanalysis and in the Community Atmosphere Model, version 6 (CAM6), historical simulation during 1980–2014 over North America. Compared to radiosondes, ERA5 successfully reproduces SLS environments, with strong spatiotemporal correlations and low biases, especially over the Great Plains. Both ERA5 and CAM6 reproduce the climatology of SLS environments over the central United States as well as its strong seasonal and diurnal cycles. ERA5 and CAM6 also reproduce the climatological occurrence of the synoptic-scale features, with the distribution pattern similar to that of SLS environments. Compared to ERA5, CAM6 exhibits a high bias in convective available potential energy over the eastern United States primarily due to a high bias in surface moisture and, to a lesser extent, storm-relative helicity due to enhanced low-level winds. Additionally, composite analysis indicates consistent synoptic anomaly patterns favorable for significant SLS environments over much of the eastern half of the United States in both ERA5 and CAM6, though the pattern differs for the southeastern United States. Overall, our results indicate that both ERA5 and CAM6 are capable of reproducing SLS environments as well as the synoptic-scale features and transient events that generate them.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1992

A PC-based collapse factor analysis for high pressure cryogenic systems

A collapse factor simulation computer program, suitable for use on personal computers, for supercritical cryogen system pressurization is examined. The available literature is previewed, individual component models are developed, and the overall simulation program logic is examined. Extensive use is made of existing thermodynamic and transport property evaluation packages. Heat transfer is estimated using natural convection correlations, and transient wall heat conduction is considered. Examples of the simulation are demonstrated for two facilities. Partial validation is achieved by comparing the simulation output with data from an existing test facility. Preliminary results indicate that the collapse factors will seldom exceed 1.15 during supercritical operation of high-pressure cryogenic systems. Recommendations for improvements in the component models and the simulation are made.

Hodge, B. K.↗

NASA NTRS1983

Numerical simulation of the Weddell Sea pack ice

In order to determine the degree to which a coupled dynamic thermodynamic model can reproduce the seasonal cycle of the Antarctic Sea ice in the Weddell Sea, a series of two-year simulations of the Weddell Sea ice pack were carried out. The simulations employed Hibler's (1979) model and used a one-day time step on an 18 x 15 grid with a resolution of 222 km. Daily atmospheric data from 1979 were used to drive the simulations. The simulations yielded a seasonal cycle of ice with maximum and minimum extents close to that observed. Except for portions of the western Weddell, the advance of the ice is found to be primarily thermodynamic in nature, while the rapid decay depends critically on the presence of both leads and lateral ice advection. In early summer the ice motion causes a residual tongue of ice to extend eastward from the Antarctic Peninsula in agreement with observations. Mean ice thicknesses are consistent with observations, as are the mean drift rates of about 5 km/day.

Hibler, W. D., III↗

DOE OSTI2023

Atomistic simulations and machine learning of solute grain boundary segregation in Mg alloys at finite temperatures

Understanding solute segregation thermodynamics is the first step in investigating grain boundary (GB) properties, such as strong yttrium (Y) effects on grain growth and texture evolution in micro-scale polycrystalline magnesium (Mg) alloys. To estimate the average GB segregation behavior in low-solute-concentration Mg alloys (e.g., 2 at.% Y), a state-of-the-art spectral approach is applied based on a per-site segregation energy spectrum for Y solute atoms at zero K obtained from molecular statistics (MS) simulations of ~10 4 GB sites in Mg symmetric tilt GBs (STGBs). Although selected MS simulation results are consistent with verification by density functional theory (DFT) calculations, estimates of average segregation tendency based on the zero-K energy spectrum deviate from experimental observations. To resolve this problem, thermodynamic integration (TI) methods based on molecular dynamics (MD) simulations are used to determine the per-site segregation free energies of Y at representative GB sites, which show contributions beyond harmonic approximations can be important for certain GB sites at high temperatures. A surrogate model of per-site segregation free energy is constructed from a small set of TI data points using stacking cross-validation regressors and physics-informed descriptors. This model is applied to predict the Y segregation free energy spectra for thousands of GB sites in Mg STGBs with uncertainty quantification. Finally, the average segregation tendency predicted by the spectral approach based on the free energy spectra agrees well (within the uncertainty range) with experimental observations for micro-scale polycrystalline Mg alloys at typical thermomechanical processing temperatures (500 ~ 800 K), where thermodynamic equilibrium states are likely to be achieved due to fast diffusion.

Atomistic simulations↗

DOE OSTI2025

Ab initio equation of state for vanadium

Accurate materials’ equations of state (EOS) are essential for understanding materials properties as well as for use in multiphysics simulations. In particular, hydrodynamics simulations are based on three fundamental conservation laws (mass, momentum and energy) that form an under-determined system of equations. The equation of state serves as an additional closure relation between thermodynamic variables for a given material that enables numerical hydrodynamics simulation. In this report, we focus on the development of an ab initio EOS for the body centered cubic (BCC) phase of Vanadium (V) for eventual integration into a multiphase EOS in the OpenSesame EOS database.

36 MATERIALS SCIENCE↗