Search NASA⌕ Search

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.

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

At least 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
NASA NTRS2004

High-Fidelity Three-Dimensional Simulation of the GE90

A full-engine simulation of the three-dimensional flow in the GE90 94B high-bypass ratio turbofan engine has been achieved. It would take less than 11 hr of wall clock time if starting from scratch through the exploitation of parallel processing. The simulation of the compressor components, the cooled high-pressure turbine, and the low-pressure turbine was performed using the APNASA turbomachinery flow code. The combustor flow and chemistry were simulated using the National Combustor Code (NCC). The engine simulation matches the engine thermodynamic cycle for a sea-level takeoff condition. The simulation is started at the inlet of the fan and progresses downstream. Comparisons with the cycle point are presented. A detailed look at the blockage in the turbomachinery is presented as one measure to assess and view the solution and the multistage interaction effects.

Turner, Mark G.↗
NASA NTRSDate not supplied

Implementation of Interactive Sea Ice Radiative Transfer into the GISS GCM and Its Impact on Sea Ice Simulation

To improve the representation of solar radiation in sea ice in the Goddard Institute for Space Studies (GISS) climate model (ModelE), we have implemented an advanced radiative transfer (RT) scheme to interactively calculate the surface albedo and transmittance of the ice layer and the solar absorption within. The ice physical properties (temperature, salinity, and density), which are simulated by the thermodynamic module in ModelE, are directly used as the input of the RT scheme to determine the inherent optical properties (IOPs) of brine pockets and air bubbles in the ice layers. This provides a more sophisticated, accurate, and complete treatment of sea ice RT for ModelE. The sea ice RT model is also validated with measured spectral albedo and transmittance for different ice types. Compared with ModelE simulations with current radiative parameterization for sea ice, the new RT implementation will slightly increase the ice surface albedo in the summer in the Arctic and greatly increase the radiative transmittance. Consequently, the ice surface melting is reduced and the base melting is increased but the difference in the simulated ice area is small. The effect of this new RT implementation on the Antarctic ice is negligible. The observed decrease in Arctic sea ice extent is well captured by the ModelE simulation.

Sea ice↗
DOE OSTI2026

Thermal exchange-correlation functionals: Capturing quantum electron behavior in warm, dense plasmas

We summarize and give perspective upon recent progress in developing non-empirical constraint-based thermal (i.e., free energy) exchange-correlation (XC) density functionals essential for accurate description of the quantum behavior of electrons in warm, dense plasmas. After delineating the critical role of ground-state functionals for zero-temperature, time-dependent DFT, we outline the underpinnings of local density approximation, generalized gradient approximation (GGA), and meta-GGA XC free-energy functionals. Two basic thermalization principles for upgrading ground-state XC functionals to successful thermal ones are emphasized. Then, we turn to a long-standing challenge, assessment of the accuracy of well-founded functionals. Unlike the ground state, there are a few exact results for large T and P. An exception is path integral Monte Carlo (PIMC) data for dense H/D and He plasmas. For those, we did ab initio molecular dynamics simulations under selected thermodynamic conditions employing five thermal XC functionals: two approximate thermal GGAs, fully thermal GGA, an approximate meta-GGA, and fully thermal meta-GGA. Comparisons with the PIMC data show that functionals thermalized by augmenting a non-thermal functional with a lower-level thermal contribution are inferior to functionals with thermal XC and spatial inhomogeneity effects taken into account at the same level of refinement. We believe this and similar evidence should be convincing to the high-energy density physics community of the necessity of use of proper thermal XC functionals in simulation studies of finite-temperature quantum effects in warm, dense plasmas.

Ab-initio molecular dynamics↗
NASA NTRS1980

Sea ice simulations based on fields generated by the GLAS GCM

The GLAS General Circulation Model (GCM) was applied to the four-month simulation of the thermodynamic part of the Parkinson-Washington sea ice model using atmospheric boundary conditions. The sea ice thickness and distribution were predicted for the Jan. 1-Apr. 30 period using the GCM-fields of solar and infrared radiation, specific humidity and air temperature at the surface, and snow accumulation; the sensible heat and evaporative surface fluxes were consistent with the ground temperatures produced by the ice model and the air temperatures determined by the atmospheric concept. It was concluded that the Parkinson-Washington sea ice model results in acceptable ice concentrations and thicknesses when used with GLAS GCM for the Jan.-Apr. period suggesting the feasibility of fully coupled ice-atmosphere simulations with these two approaches.

Parkinson, C. L.↗
DOE OSTI2020

Prediction of structure and cation ordering in an ordered normal-inverse double spinel

Abstract Spinels represent an important class of technologically relevant materials, used in diverse applications ranging from dielectrics, sensors and energy materials. While solid solutions combining two “single spinels” have been explored in a number of past studies, no ordered “double” spinels have been reported. Based on our first principles computations, here we predict the existence of such a double spinel compound MgAlGaO 4 , formed by an equimolar mixing of MgAl 2 O 4 normal and MgGa 2 O 4 inverse spinels. After studying the details of its atomic and electronic structure, we use a cluster expansion based effective Hamiltonian approach with Monte Carlo simulations to study the thermodynamic behavior and cation distribution as a function of temperature. Our simulations provide strong evidence for short-ranged cation order in the double spinel structure, even at significantly elevated temperatures. Finally, an attempt was made to synthesize the predicted double spinel compound. Energy Dispersive X-ray Spectrometry and X-ray diffraction Rietveld refinements were performed to characterize the single-phase chemical composition and local configurational environments, which showed a favorable agreement with the theoretical predictions. These findings suggest that a much larger number of compounds can potentially be realized within this chemical space, opening new avenues for the design of spinel-structured materials with tailored functionality.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Multiple Environmental Influences on the Lightning of Cold-Based Continental Convection. Part II: Sensitivity Tests for Its Charge Structure and Land–Ocean Contrast

In Part I, an electrification scheme was described and a simulation of an observed cold-based storm from the U.S. Great Plains was validated with electrical observations. Most charge in the storm was separated by rebounding collisions of secondary ice originating from prior graupel–snow collisions. In this Part II, sensitivity tests are performed with the control simulation (Part I) and influences from environmental factors (aerosols, temperature, moisture, and shear) on lightning are elucidated. Environmental factors [e.g., convective available potential energy (CAPE)] controlling updraft speed are salient. When ascent is reduced by 30% and 70%, flashes become 70% fewer and disappear, respectively; faster ascent promotes positive cloud-to-ground (+CGs) flashes. Since cloud base is too cold (1°C) for coalescence, cloud condensation nucleus aerosol concentrations do not influence the lightning appreciably. The electrical response to varying concentrations of active ice nuclei is limited by most ice particles being secondary and less sensitive—a natural “buffer.” Imposing a maritime sounding suggests that the land–sea contrast in lightning for such storms is due to the vertical structure of environmental temperature and humidity. Weak CAPE, and both entrainment and condensate weight from a low cloud base, suppress ascent and charging. Maritime thermodynamic conditions reduce simulated flash rates by two orders of magnitude. Reducing aerosol loadings from continental to maritime only slightly reinforces this suppression. Last, a conceptual model is provided for how any simulated storm is either normal because graupel/hail is mostly positively charged or else is inverted/anomalous because graupel/hail is mostly negatively charged, with environmental factors controlling the charging. Impacts from microphysical processes, including three processes of secondary ice production, on lightning are analyzed.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

As a Matter of State: The Role of Thermodynamics in Magnetohydrodynamic Turbulence

Turbulence simulations play a key role in advancing the general understanding of the physical properties of turbulence and in interpreting astrophysical observations of turbulent plasmas. For the sake of simplicity, however, turbulence simulations are often conducted in the isothermal limit. Given that the majority of astrophysical systems are not governed by isothermal dynamics, we aim to quantify the impact of thermodynamics on the physics of turbulence, through varying adiabatic index, γ, combined with a range of optically thin cooling functions. Here, we present a suite of ideal magnetohydrodynamics simulations of thermally balanced stationary turbulence in the subsonic, super-Alfvénic, high ${\beta }_{{\rm{p}}}$ (ratio of thermal to magnetic pressure) regime, where turbulent dissipation is balanced by two idealized cooling functions (approximating linear cooling and free–free emission) and examine the impact of the equation of state by considering cases that correspond to isothermal, monatomic, and diatomic gases. We find a strong anticorrelation between thermal and magnetic pressure independent of thermodynamics, whereas the strong anticorrelation between density and magnetic field found in the isothermal case weakens with increasing γ. Similarly, the linear relation between variations in density and thermal pressure with sonic Mach number becomes steeper with increasing γ. This suggests that there exists a degeneracy in these relations with respect to thermodynamics and Mach number in this regime, which is dominated by slow magnetosonic modes. These results have implications for attempts to infer (e.g.,) Mach numbers from (e.g.,) Faraday rotation measurements, without additional information regarding the thermodynamics of the plasma. However, our results suggest that this degeneracy can be broken by utilizing higher-order moments of observable distribution functions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Uncertainty Quantification in Atomistic Modeling of Metals and Its Effect on Mesoscale and Continuum Modeling: A Review

The design of next-generation alloys through the integrated computational materials engineering (ICME) approach relies on multiscale computer simulations to provide thermodynamic properties when experiments are difficult to conduct. Atomistic methods such as density functional theory (DFT) and molecular dynamics (MD) have been successful in predicting properties of never before studied compounds or phases. However, uncertainty quantification (UQ) of DFT and MD results is rarely reported due to computational and UQ methodology challenges. Over the past decade, studies that mitigate this gap have emerged. These advances are reviewed in the context of thermodynamic modeling and information exchange with mesoscale methods such as the phase-field method (PFM) and calculation of phase diagrams (CALPHAD). The importance of UQ is illustrated using properties of metals, with aluminum as an example, and highlighting deterministic, frequentist, and Bayesian methodologies. Finally, challenges facing routine uncertainty quantification and an outlook on addressing them are also presented.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Hybrid machine learning/physics-based approach for predicting oxide glass-forming ability

Predicting the liquid compositions that will vitrify at experimentally accessible quench rates remains one of the grand challenges in the field of condensed matter physics. This glass-forming ability can be quantified as the critical quench rate needed to suppress crystallization. Knowledge of this critical quench rate also informs which glass composition could be used for new applications. There have been several physical and empirical models presented in the literature to predict the critical quench rate/glass forming ability. These models range from those theoretically derived to those quantified only through experimental characterization. In this work, we instead propose a new method to calculate the critical quench rate using the recently developed toy landscape model combined with machine learning. The toy landscape model accesses the underlying physics that control the vitrification behavior by directly simulating the liquid thermodynamics and kinetics. In conclusion, the results are discussed in terms of industrial impact, physical insights, and how the glass science community can develop improved predictions of glass-forming ability.

Crystallization↗
DOE OSTI2021

Thermodynamic modeling of hydrogen fueling process from high-pressure storage tank to vehicle tank

This study develops a hydrogen fueling station (HFS) thermodynamic model that simulates the actual fueling process in which hydrogen is supplied from a high-pressure (HP) storage tank into a fuel cell electric vehicle (FCEV) tank. To make the model as accurate as possible, we use the same components and specifications as in actual HFSs, such as a pressure control valve, a pre-cooling system, and an FCEV tank. After the components and their specifications are set, pressure and temperature profiles are set as the HP tank supply conditions. Based on the pressure and temperature profiles, the model solves for the temperature, pressure, and mass flow rate of hydrogen at each downstream position, including the inside of the vehicle tank. The values predicted by the model are compared with experimental data, and we show that the developed model makes it possible to accurately simulate those values at any position during the fueling process.

08 HYDROGEN↗
DOE OSTI2025

Directed energy deposition of radially functionally graded nickel based superalloys IN625 and IN738: Microstructure and mechanical properties

This study investigates the microstructure and mechanical properties of radially functionally graded (RFG) materials fabricated from nickel-based superalloys IN625 and IN738 using directed energy deposition (DED). RFG cylindrical deposits with radial gradients of IN625-IN738 and IN738-IN625 were successfully fabricated and confirmed to be crack-free through tailored composition grading, scan strategy and process parameters that minimized thermal gradients and solidification velocities. Microstructural analysis revealed distinct variations across the gradient, characterized by differences in dendrite arm spacing, microsegregation, and the formation of distinct γ' precipitates and carbide/boride phases (MC, M 6 C, M 2 C, and MB 2 ) in as-deposited (AD) conditions, which evolved into M 23 C 6 and M 2 B upon subsequent heat treatment (HT). EBSD showed predominant columnar grains with a [100] texture, while microhardness profiles demonstrated consistent variation across the gradients. Mechanical testing highlighted the synergy of the graded structure, combining IN738's strength with IN625's ductility, resulting in enhanced properties such as a yield strength of 1131 MPa, ultimate tensile strength of 1382 MPa, and elongation of 12 % for HT IN625-IN738 grading. Numerical simulations employing CALPHAD thermodynamic calculations and finite element thermal modeling are used to establish phase evolution and microstructural feature correlation with thermal profiles, aiding in the optimization of process-structure-property relationship.

Additive manufacturing↗
DOE OSTI2022

Depletion-Driven Assembly of Polymer-Coated Nanocrystals

Depletion-driven assembly has been widely studied for micron-sized colloids, but questions remain at the nanoscale where the governing physics are impacted by stabilizing surface ligands or wrapping polymers, whose dimensions are similar to those of the colloidal core and the depletant. Here, we probe how wrapping colloidal tin-doped indium oxide nanocrystals with polymers affects their depletion-induced interactions and assembly in solutions of polyethylene glycol. Copolymers of polyacrylic acid grafted with polyethylene oxide provide nanocrystal wrappings with different effective polymer graft densities and molecular weights. (Ultra) small angle X-ray scattering, coarse-grained molecular dynamics simulation, and molecular thermodynamic theory were combined to analyze how depletant size and polymer wrapping characteristics affect depletion interactions, structure, and phase behavior. The results show how depletant molecular weight, as well as surface density and molecular weight of polymer grafts, set thresholds for assembly. These signatures are unique to depletion-driven assembly of nanoscale colloids and mirror phase behaviors of grafted nanoparticle–polymer composites. Optical and rheological responses of depletion-driven assemblies of nanocrystals with different polymer shell architectures were probed to learn how their structural differences impact properties. Finally, we discuss how these wrinkles for depletion-driven assembly at the nanoscale may provide fresh opportunities for designing responsive depletion interactions and dynamically reconfigurable materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Observed and Bin Model Simulated Evolution of Drop Size Distributions in High-Based Cumulus Congestus Over the United Arab Emirates

This study examines microphysical processes in developing high-based cumulus congestus over the United Arab Emirates using aircraft observations and a large-eddy-simulation model with bin microphysics. A notable feature of this case is the lack of mm-sized drops despite having a liquid cloud layer >3.5 km deep, contrasting with copious large drops observed in maritime tropical cumulus congestus having a similar vertical extent. Modeled drop size distributions are similar to observations at various temperatures between 9.5 and –12°C, including the lack of mm-sized drops. Cloud dilution leads to low-to-moderate liquid water contents (~0.5–1.5 g m –3 ) in most of the cloud core, several times smaller than adiabatic values. Dilution is enhanced in the inflowing branch of the toroidal circulations associated with individual cloud thermals, which are favored regions for secondary droplet activation (activation above cloud base). Secondary activation in general contributes substantially to the droplet population. Turning it off leads to a sharp decrease in droplet concentration and increase in mean size aloft, but does little to increase rain drop production. Warm rain generation (or lack thereof) in this case is therefore determined more by the sub-cloud aerosol and the cloud base droplet size distribution (DSD) than DSD evolution aloft from secondary droplet activation. Decreasing the aerosol concentration by a factor of 10 greatly increases production of large drops via collision-coalescence. Thus, despite its high base (low temperatures) and substantial dilution, the simulated cloud is thermodynamically and dynamically capable of rapidly producing copious mm-sized drops from collision-coalescence under pristine aerosol conditions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Effects of Urbanization in China on the East Asian Summer Monsoon as Revealed by Two Global Climate Models

Abstract The East Asian summer monsoon (EASM), vulnerable to anthropogenic activities, exerts profound impacts on precipitation in East Asia. In recent decades, China has been experiencing rapid urbanization. It has not yet reached a consensus on its influences on EASM. Here, we introduce satellite‐derived urban impervious areas of China during 1985–2017 into two global climate models with contrasting biases of climatological EASM simulation. Both models consistently show that urbanization weakens EASM by increasing surface friction drag. The increased surface roughness reduces climatological southerly winds in the concentrated urbanization area. As a result, moisture convergence occurs south of the concentrated urbanization area with ascending motion, while moisture divergence with descending motion appears to the north, resulting in the precipitation anomalies of “southern flood and northern drought.” The dominant role of urban dynamic forcing in regulating EASM is confirmed by reducing urban roughness to the level of croplands in sensitivity experiments. As urban dynamic forcing is largely mitigated, the urbanization impact mainly driven by urban thermal forcing conversely strengthens EASM due to enhanced land‐sea thermal contrast via urbanization‐induced warming. Spatial patterns and magnitudes of the urbanization‐induced EASM weakening vary between the two models due to disparities in their simulations of urban thermodynamic forcing and atmospheric adjustments. These results underscore the robust weakening effect of urbanization on EASM by urban friction drag and document uncertainties in the related spatial pattern and magnitude resulting from model discrepancies.

Meteorology & Atmospheric Sciences↗
DOE OSTI2020

A magma ocean origin to divergent redox evolutions of rocky planetary bodies and early atmospheres

Magma oceans were once ubiquitous in the early solar system, setting up the initial conditions for different evolutionary paths of planetary bodies. In particular, the redox conditions of magma oceans may have profound influence on the redox state of subsequently formed mantles and the overlying atmospheres. The relevant redox buffering reactions, however, remain poorly constrained. Using first-principles simulations combined with thermodynamic modeling, we show that magma oceans of Earth, Mars, and the Moon are likely characterized with a vertical gradient in oxygen fugacity with deeper magma oceans invoking more oxidizing surface conditions. This redox zonation may be the major cause for the Earth’s upper mantle being more oxidized than Mars’ and the Moon’s. These contrasting redox profiles also suggest that Earth’s early atmosphere was dominated by CO 2 and H 2 O, in contrast to those enriched in H 2 O and H 2 for Mars, and H 2 and CO for the Moon.

58 GEOSCIENCES↗
DOE OSTI2025

Using magnetic dynamics to measure the spin gap in a candidate Kitaev material

Spin-orbit entangled materials have attracted widespread interest due to the novel magnetic phenomena arising from the interplay between spin-orbit coupling and electronic correlations. However, the intricate nature of spin interactions within Kiteav materials complicates the precise measurement of low-energy spin excitations. Using Na 2 Co 2 TeO 6 as an example, we study these low-energy spin excitations using the time-resolved resonant elastic x-ray scattering (tr-REXS). Our observations unveil remarkably slow spin dynamics at the magnetic peak, whose recovery timescale is several nanoseconds. This timescale aligns with the extrapolated spin gap of ~1 μeV, obtained by density matrix renormalization group (DMRG) simulations in the thermodynamic limit. The consistency demonstrates the efficacy of tr-REXS in discerning low-energy spin gaps inaccessible to conventional spectroscopic techniques.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Elucidating the role of reactive nitrogen intermediates in hetero-cyclization during hydrothermal liquefaction of food waste

Hydrothermal liquefaction (HTL) has emerged as a promising strategy for converting abundant, water-rich organic streams into an energy-dense, sustainable, biocrude. Despite major strides in improving oil yields and process intensification, a key underlying challenge persists in relating the composition of the feedstock to the ultimate fate of nitrogen, which is especially important for biocrude quality. To elucidate how nitrogen partitions into the respective gas, aqueous, oil and solid hydrothermal products, food waste was treated under HTL conditions (15 wt% solids, 575 K reaction temperature, and 1 hour reaction time) and elemental, chromatographic and mass spectroscopy analyses were performed on the products to establish dominant product classes and to close mass, carbon, and nitrogen balances. Liquid products (aqueous and biocrude) were found to contain nearly 60% of the nitrogen, with a majority of the biocrude-phase nitrogen in the form of amides. Ab initio thermodynamic and kinetic simulations on model reactants were used to evaluate potential reaction pathways involving reactive nitrogen intermediates. A subsequent kinetic model was evaluated at reaction temperatures ranging 300–1000 K and for a range of feedstock compositions to identify primary reaction pathways. The Maillard reaction is revealed to be the dominant pathway for converting reactive nitrogen intermediates into observed products, especially those in which nitrogen is present as a heterocyclic aromatic. Subsequent product analysis by Fourier Transform Ion Cyclotron Resonance Mass Spectrometry (FT-ICR MS) confirms strong agreement among the experimental data, kinetic model and proposed Maillard pathways.

09 BIOMASS FUELS↗