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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
NASA NTRS2004

Thermodynamic Vent System Test in a Low Earth Orbit Simulation

A thermodynamic vent system for a cryogenic nitrogen tank was tested in a vacuum chamber simulating oxygen storage in low earth orbit. The nitrogen tank was surrounded by a cryo-shroud at -40 F. The tank was insulated with two layers of multi-layer insulation. Heat transfer into cryogenic tanks causes phase change and increases tank pressure which must be controlled. A thermodynamic vent system was used to control pressure as the location of vapor is unknown in low gravity and direct venting would be wasteful. The thermodynamic vent system consists of a Joule-Thomson valve and heat exchanger installed on the inlet side of the tank mixer-pump. The combination is used to extract thermal energy from the tank fluid, reducing temperature and ullage pressure. The system was sized so that the tank mixer-pump operated a small fraction of the time to limit motor heating. Initially the mixer used sub-cooled liquid to cool the liquid-vapor interface inducing condensation and pressure reduction. Later, the thermodynamic vent system was used. Pressure cycles were performed until steady-state operation was demonstrated. Three test runs were conducted at tank fills of 97, 80, and 63 percent. Each test was begun with a boil-off test to determine heat transfer into the tank. The lower tank fills had time averaged vent rates very close to steady-state boil-off rates showing the thermodynamic vent system was nearly as efficient as direct venting in normal gravity.

VanOverbeke, Thomas J.↗
DOE OSTI2025

Flower‐Type Organized Trade‐Wind Cumulus: A Multi‐Day Lagrangian Large Eddy Simulation Intercomparison Study

Shallow cumulus cloud fields in subtropical marine trade wind environments, particularly over the tropical Atlantic Ocean, show distinct organizational patterns. Among these, Flower‐type clouds are characterized by expansive stratiform cloud patches surrounded by regions of scattered convection. The objectives of this study were (a) to construct a case study of a time period during the EUREC 4 A/ATOMIC field campaign when Flower‐type organization was observed, (b) to evaluate the fidelity of a multi‐model ensemble of large eddy simulations of that case, and (c) to analyze the interaction between cloud and precipitation processes and mesoscale organization in the simulations. The simulations follow a quasi‐Lagrangian trajectory, allowing mesoscale features to develop over time in a domain that follows the boundary‐layer airmass. The results show a broad agreement in simulated thermodynamic properties across different LES codes, with Flower‐type cloud patches appearing within hours of each other. The consensus among models is consistent with observations made during the EUREC 4 A/ATOMIC field campaign on the specific day of interest. The cloud structure reveals three distinct peaks in the joint probability densities of cloud base and cloud top height, with the dominant peak at any given time influenced by the stage of cloud organization. The simulated cloud system evolution reveals consistent occurrence of maxima in liquid water path and rain rate before Flower reaches its maximum length scale. Targeted sensitivity tests reveal a weak relationship between Cloud Droplet Number concentration and the extent/degree/type of organization.

EUREC4A↗
DOE OSTI2026

CLM5 Simulations of Soil Moisture and Gross Primary Productivity for CONUS at 0.125 degrees

This dataset provides 0.125-degree gridded simulations of soil moisture and gross primary productivity (GPP) for the Contiguous United States (CONUS), generated using the Community Land Model version 5 (CLM5) with the biogeochemistry module enabled. The data covers a historical baseline (1980-2015) and mid-century future projections (2020-2055). Future projections are organized into two sets of scenarios to distinguish the impacts of different drivers: (1) Atmospheric Only (ATM): These scenarios apply future atmospheric forcings while holding land use and land cover (LULC) at historical baseline levels. The atmospheric forcings represent moderately versus severely hotter/drier atmospheric conditions (dynamically downscaled perturbed thermodynamics simulations based on CMIP6 SSP245 and SSP585 warming signals), each with cooler versus hotter Earth System Model temperature sensitivity instantiations. These scenarios are identified in the folder names as rcp45_cooler_near, rcp45_hotter_near, rcp85_cooler_near, and rcp85_hotter_near. (2) Coupled Atmospheric and Land-Use (LAND+ATM): These scenarios apply future atmospheric forcing together with future LULC by pairing atmospheric pathways with lower versus higher population/economic growth scenarios representing Shared Socioeconomic Pathways 3 and 5 (SSP3 and SSP5). These scenarios are identified in the file names as ssp3_rcp45_cooler_near, ssp3_rcp45_hotter_near, ssp5_rcp85_cooler_near, and ssp5_rcp85_hotter_near. Please refer to the "README_first.md" file for detailed information on file structure, variables, units, and data formats.

drought↗
DOE OSTI2026

CLM5 Simulations of Soil Moisture and Grain Carbon for CONUS at 0.125 degrees

This dataset provides 0.125° gridded simulations of soil moisture and crop grain carbon for the Contiguous United States (CONUS), generated using the Community Land Model version 5 (CLM5) with the biogeochemistry module enabled. The data covers a historical baseline (1980–2015) and mid-century future projections (2020–2055). Future projections are organized into two sets of scenarios to distinguish the impacts of different drivers: (1) Atmospheric Only (ATM-only): These scenarios apply future atmospheric forcings while holding land use and land cover change (LULCC) at historical baseline levels. The atmospheric forcings represent moderately versus severely hotter/drier atmospheric conditions (dynamically downscaled perturbed thermodynamics simulations based on CMIP6 SSP245 and SSP585 warming signals), each with cooler versus hotter Earth System Model temperature sensitivity instantiations. These scenarios are identified in the file names as atm45cooler, atm45hotter, atm85cooler, and atm85hotter. (2) Coupled Atmospheric and Land-Use (LAND+ATM): These scenarios apply future atmospheric forcing together with future LULCC by pairing atmospheric pathways with lower versus higher population/economic growth scenarios representing Shared Socioeconomic Pathways 3 and 5 (SSP3 and SSP5). These scenarios are identified in the file names as atm45cooler_ssp3, atm45hotter_ssp3, atm45cooler_ssp5, atm45hotter_ssp5, atm85cooler_ssp3, atm85hotter_ssp3, atm85cooler_ssp5, and atm85hotter_ssp5. Please refer to the README file for detailed information on file structure, variables, units, and data formats.

Yao, Lili [Pacific Northwest National Laboratory] ↗
NASA NTRS1992

Steady and transient performance calculation method for prediction, analysis, and identification

The detailed design and development of turbofans Involves the prediction and identification, by means of test analysis, of the performance of the engine and its components. The thermodynamic simulation and analysis codes integrate existing knowledge and Interpretations of the detailed operating procedure of the components of the engine being developed. The relevance of the predicted performance depends on the quality of the representation of the various physical phenomena affecting the characteristics of the components and, consequently, on the incorporation of experimental correlation in the modelling. In this context, the representation of compressor and turbine characteristics is particularly important. Firstly, we will analyze the ability of corrected parameters to represent MACH similitude at the component inlet under various conditions. The various measurements achieved on the engine during development are used to enhance the modelling. The methods of identifying the thermodynamic calculation code with the various measurements, considered here with their uncertainties, are then presented and described. The analysis of the tests performed on the powerplant, designed to Identify the real characteristics of the engine components, can be undertaken, considering one or more engine test points and incorporating knowledge acquired through experimentation or on the component test bench. There are many possible fields of use for these identification methods ranging from the rematching of components to the optimization of the control system.

Jean Pierre Duponchel↗
DOE OSTI2023

Microstructural refinement in ultrasonically modified A356 aluminum castings

Two A356 aluminum alloys (Al-Si-Mg), one with 0.09 wt.% Fe and one with 0.91 wt.% Fe, were cast in a graphite mold with the simultaneous application of local ultrasonic intensification to refine the as-cast microstructure. Ultrasonication during casting transformed the morphology of primary Al grains from dendritic (~140-290 microns in size) to globular (~33-36 microns in size). The alloy with high Fe exhibited globular grains at distances up to 45 mm away from the ultrasound probe, while the alloy with low Fe exhibited globular grains at distances only up to 6 mm away from the ultrasound probe. Near the location of the ultrasound probe (< 2 mm away), a second non-dendritic microstructural morphology was observed with fine aluminum grains (~9-25 microns in size). This unique fine-grained morphology has not been previously reported, contains a greater concentration of Si relative to the globular microstructure, and may be a large, fully eutectic region. Ultrasonication during casting also transformed the morphology of the ß-Al 5 FeSi phase particles (which are deleterious to the strength and ductility of the alloy) in the high Fe alloy from needle-like to rectangular, which could enable the greater use of secondary Al alloys. Thermodynamic simulations conducted to calculate the solidification paths of the two alloys studied predict that the ß-Al 5 FeSi phase begins to form earlier in the alloy with high Fe. Finally, data suggest that the ß-Al 5 FeSi phase (which is more abundant in alloys with high Fe content) may enhance ultrasonically-induced grain refinement.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Strontium Effects on the Formation of Iron-Intermetallic Phases in Secondary Al–9Si–0.6Fe Alloys

The presence and morphology of Fe-containing intermetallic phases affect the mechanical properties of aluminum alloys, especially in secondary Al–Si-based cast alloys. Although strontium (Sr) addition of 50 to 500 ppm is known to refine the needle-type eutectic silicon structure, the influence of Sr on the formation of Fe-intermetallic phases remains unclear. The present work investigates the combined additions of Sr and Mn to Al–9Si–0.6Fe–0.35Mg (All compositions are in wt pct except otherwise stated.) alloys on the formation of Fe-intermetallic phases at different solidification rates from ~ 1.5 to ~ 60 °C/s. Long and branched-type AlFeSi phase with size ranging from 50 to 120 µm are more common when solidified at the rate of 1.5 °C/s regardless of Sr and Mn additions. However, at the fast solidification rate of 60 °C/s, a 60 ppm Sr addition significantly reduced the average length of needle-shaped AlFeSi phase to less than 3 to 5 µm. Thermodynamic simulations have been performed using CALculation of PHAse Diagrams (CALPHAD) models to predict the formation of various phases and their possible interactions during solidification. The results indicated that the combination of a high solidification rate and about 60 ppm of Sr is beneficial to refining the δ-Al 3 FeSi 2 phase in Al–Si–Mg alloys containing 0.6 pctFe. As a result, this unexpected finding of Fe-intermetallic refinement by low Sr addition (~60 ppm) provides an important guide in designing secondary alloys for sustainable casting applications.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Functionalized benzylamines from commercial kraft lignin

Benzylamines are key intermediates in pharmaceuticals, agrochemicals, and polymers, but their conventional production relies on benzyl chloride - a petroleum-derived compound with high toxicity and energy demands. Lignin, accounting for up to 30% of plant biomass, is the largest renewable source of aromatic carbon on Earth. However, its highly complex and recalcitrant structure poses a major barrier to efficient conversion into high-value chemicals. Here, in this study, we developed a catalytic approach to convert commercial kraft lignin into phenolic benzylamines through selective depolymerization and subsequent functionalization. We systematically evaluated the effects of three alcohol solvents, formic acid (FA), and a ruthenium-on‑carbon (Ru/C) catalyst on monophenol yield and selectivity. Up to 6.5 wt% monophenol yield was achieved using methanol (MeOH), FA, and Ru/C at 300 °C for 2 h. Quantum thermodynamic simulations based on the COSMO-RS model confirmed the superior solvation and reactivity of the MeOH + FA system, rationalizing observed product yield. The purified monophenolic products, primarily guaiacol and alkyl guaiacols, were then converted into functionalized benzylamines with >90% yield via a multicomponent Mannich reaction under mild conditions. Techno economic analysis (TEA) and life cycle assessment (LCA) underscore the importance of improving lignin depolymerization yields and expanding biorefinery scale. Solvent-only configurations outperform other options in both cost and emissions, with the methanol-only case performing the best ($\$$105 /kg and 26 kg CO 2e /kg) at a large-scale facility. This study establishes a scalable, bio-based pathway for producing benzylamines from commercial kraft lignin, advancing lignin valorization and offering a sustainable alternative to produce petrochemical-based benzylamines.

Catalytic reduction↗
DOE OSTI2023

Effect of Cr and Ni concentrations on resilience of cast Nb- alloyed heat resistant austenitic steels at extreme high temperatures

Austenitic Cr–Ni alloyed heat resistant steels with Nb additions are used for intensively thermo-mechanically loaded cast components working in extreme high temperature oxidizing environment. Their performance during static oxidation and transient thermo-mechanical loading was investigated to recommend an optimal cost-effective Cr/Ni composition of Nb-alloyed austenitic class steels. Here, the static oxidation and transient thermo-mechanical behavior of three austenitic steels with different Cr/Ni alloying levels were investigated and compared for variety of working conditions. Static oxidation was performed between 900 °C and 1000 °C in air for 400 h. The critical temperature which increases spallation during static oxidation was determined for each of the steel alloying levels. In addition, thermal cycling of a constrained specimen was done with varying upper cycling temperatures between 850 °C and 1000 °C. SEM and TEM analyses were supported by thermodynamic simulation of the phases precipitated in the metal matrix and the structure of formed oxide layers. These studies were used to determine the mechanisms of degradation of thermo-mechanically loaded Cr/Ni austenitic steels at extreme high temperatures. A recommendation for a cost-effective Cr/Ni alloying level for different working conditions was determined.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Sigma phase evolution and nucleation mechanisms revealed by atom probe tomography in a 347H stainless steel

The size evolution and phase composition of sigma phases were investigated for a 347H stainless steel subject to isothermal aging at 750 °C for up to 10,000 h. Scanning electron microscopy reveals that sigma phases are present before 336 h aging and they continue to grow for all aging times studied. Atom probe tomography shows little variation in the sigma phase composition with aging time up to 10,000 h, and only the C content decreases with aging time. The comparison of C concentration suggests that the sigma phase does not nucleate within the C-depleted matrix or precipitate free grain boundaries, leaving dissolving metastable M 23 C 6 as the only possible nucleation site. Based on the experimental results and thermodynamic simulations, the desirable Nb concentration is discussed to reduce the kinetics of sigma phase growth, and possible design directions are suggested to improve 347H.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Sustainable bioleaching of lithium-ion batteries for critical metal recovery: Process optimization through design of experiments and thermodynamic modeling

Recycling spent lithium-ion batteries (LIBs) could alleviate supply risks for critical metals and be less harmful to the environment compared to new production of metals from mining. Developing a cost-effective LIB bioleaching process could be a promising alternative to traditional energy-intensive recycling technologies. Here, this study aimed to optimize bioleaching conditions for maximum economic competitiveness through design of experiments using iterative response surface methodology (RSM), assisted by thermodynamic modeling. The optimal condition was identified as 2.5% pulp density in 75 mM gluconic acid biolixiviant at 55°C for 30 h which could recover 57%–84% of nickel, 71%–86% of cobalt, and 100% of lithium and manganese, yielding a 17%–26% net profit margin. The recommended pulp density and acid concentrations, together with the observed metal solubilization, were supported by thermodynamic modeling predictions. Our study demonstrated that combining RSM with thermodynamic simulations could be a powerful tool for optimizing bioleaching conditions.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Enhanced Convective Microphysics Scheme and Its Impacts on Mean Climate in E3SM

Abstract To improve the representation of microphysical processes in convective clouds and their interaction with aerosol and stratiform clouds, a two‐moment convective microphysics parameterization (CMP) scheme developed by Song and Zhang (2011, https://doi.org/10.1029/2010jd014833 ) is upgraded and implemented in E3SM. The new developments include: (a) implementing a parameterization for graupel to enhance the representation of ice‐phase microphysical processes; (b) representing the impact of spatial inhomogeneity of cloud droplets in cumulus ensembles on autoconversion and accretion processes to improve the representation of warm‐rain microphysical processes; (c) implementing a comprehensive Bergeron process parameterization to better represent mixed‐phase microphysical processes; and (d) representing the interactions between ice‐phase microphysics and cloud thermodynamics. Simulations show that the cloud microphysical properties simulated by the CMP are generally in good agreement with observations. It reasonably simulates the changes in droplets effective radius related to precipitation formation in convective clouds, as identified from satellite observations. It also successfully simulates the contrast in these processes between maritime and continental clouds, demonstrating its capability to simulate the impact of aerosols on convection. Analyses of the impact of CMP on climate mean state simulation demonstrate that the CMP slightly improves the simulations of precipitation, cloud macrophysical properties, longwave cloud radiative forcing, zonal wind, and temperature. However, a degradation in shortwave cloud radiative forcing occurs.

GCM↗
DOE OSTI2021

In situ synchrotron diffraction and modeling of non-equilibrium solidification of a MnFeCoNiCu alloy

The solidification mechanism and segregation behavior of laser-melted Mn 35 Fe 5 Co 20 Ni 20 Cu 20 was firstly investigated via in situ synchrotron x-ray diffraction at millisecond temporal resolution. The transient composition evolution of the random solid solution during sequential solidification of dendritic and interdendritic regions complicates the analysis of synchrotron diffraction data via any single conventional tool, such as Rietveld refinement. Therefore, a novel approach combining a hard-sphere approximation model, thermodynamic simulation, thermal expansion measurement and microstructural characterization was developed to assist in a fundamental understanding of the evolution of local composition, lattice parameter, and dendrite volume fraction corresponding to the diffraction data. This methodology yields self-consistent results across different methods. Via this approach, four distinct stages were identified, including: (I) FCC dendrite solidification, (II) solidification of FCC interdendritic region, (III) solid-state interdiffusion and (IV) final cooling with marginal diffusion. It was found out that in Stage I, Cu and Mn were rejected into liquid as Mn 35 Fe 5 Co 20 Ni 20 Cu 20 solidified dendritically. During Stage II, the lattice parameter disparity between dendrite and interdendritic region escalated as Cu and Mn continued segregating into the interdendritic region. After complete solidification, during Stage III, the lattice parameter disparity gradually decreases, demonstrating a degree of composition homogenization. The volume fraction of dendrites slightly grew from 58.3 to 65.5%, based on the evolving composition profile across a dendrite/interdendritic interface in diffusion calculations. Postmortem metallography further confirmed that dendrites have a volume fraction of 64.7 ± 5.3% in the final microstructure.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Multiphysics modeling of accelerators through code integration

This work aims to improve the ability of particle accelerator researchers to develop high-performance accelerator cavity designs by creating an overall multiphysics framework that integrates and couples existing application codes. This framework will allow accelerator researchers to build multiphysics models that will optimize cavity design, improve understanding of whole-device performance, and reduce the development and fabrication costs of accelerator research. We utilize the open-source VizSchema data standard as an intermediate data structure interface layer to standardize interfaces between individual application codes. VizScema is extensively documented online, and plugins for VizSchema are available for popular visualization packages, including VisIt and ParaView. Currently, the work focuses on coupling the EM field solver COMSOL and the electron gun code MICHELLE to allow COMSOL field-solve results to be seamlessly used by MICHELLE for particle-solve. Later work will extend this integration to include other fields, particles, and thermodynamics simulation codes.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Machine-Learned Force Field Modeling of Metal Organic Frameworks for CO2 Direct Air Capture

To cope with legacy greenhouse gas emissions and to achieve net-zero emissions by 2050, the U.S. Department of Energy (DOE) is funding efforts to develop direct air capture (DAC), a method for removing CO2 directly from air. Metal organic frameworks (MOFs) have been well studied as DAC sorbent materials due to their tunable structural and compositional properties. Thermodynamic simulations using force fields are often used to provide predictions of a material’s performance in many separations. However, these force fields often make assumptions about bonds and the physics of the adsorption process. A new class of force fields called machine-learned force fields (MLFFs) use machine learning to form quantitative relationships between a material’s chemical structure and the forces and energies predicted by more accurate quantum mechanical calculations, such as dispersion-corrected density functional theory (DFT). In this work, models were developed to achieve DFT-level accuracy for the forces and energies associated with MOF flexibility and CO2 adsorption using MLFFs. These methods were parametrized based on thousands of DFT calculations of CO2 in flexible MOFs and used to predict MOF structural properties as well as CO2 adsorption in several MOFs.

Findley, John↗
DOE OSTI2025

IM3 Data Center Driven Grid Stress Dataset for the U.S. Western Interconnection

This dataset provides projected grid stress and reliability results (including all model inputs and outputs from an open-source grid operations modeling framework - GO), for the Integrated Multisector, Multiscale Modeling (IM3) project, under varying levels of data center demand growth between 2025 and 2035 in the U.S. Western Interconnection. The scenarios and sensitivity experiments are combinations of different data center demand growth rates and energy, weather, population and economic pathways. Data center demand growth projections were sourced from the Electric Power Research Institute (EPRI). The data center demand growth projection names are: Low (3.71% annual data center demand growth) Moderate (5% annual data center demand growth) High (10% annual data center demand growth) Higher (15% annual data center demand growth) Energy, weather, population and economic pathways are informed by two Shared Socioeconomic Pathways (SSP3 and SSP5) and two Representative Concentration Pathways (RCP4.5 and RCP8.5) following the hotter general circulation model (GCM) forcing group from a set of perturbed thermodynamics simulations. The resulting pathway names are: rcp45hotter_ssp3 rcp45hotter_ssp5 rcp85hotter_ssp3 rcp85hotter_ssp5 The main scenarios and sensitivity experiments are detailed below. Reference scenario: The projected grid stress and reliability results for the U.S. Western Interconnection from a previous study. This scenario does not consider data center demand growth explicitly. Data center scenario: Building on the reference scenario, this scenario considers various data center growth rates and how they impact the U.S. Western Interconnection. Data center loads are modeled as flat 8760-hr profiles. This scenario does not consider new generation and transmission capacities specifically designed to meet the new data center demands. The related folder is named "flat". Delayed generator retirements sensitivity experiment: Building on the data center scenario, this experiment explores the impact of different levels of natural gas and nuclear generator retirement delays. The resulting scenario names are: (1) postponing 100% nuclear retirements; (2) postponing 100% nuclear and 25% natural gas retirements; (3) postponing only 50% natural gas retirements; (4) postponing 100% nuclear and 50% natural gas retirements; (5) postponing 100% nuclear and 75% natural gas retirements; and (6) postponing 100% nuclear and 100% natural gas retirements. The related folder names are: no_gen_retire_0_gas, no_gen_retire_25_gas, no_gen_retire_50_gas, no_gen_retire_50_gas_only, no_gen_retire_75_gas, and no_gen_retire_100_gas. Demand response through curtailment sensitivity experiment: Building on the data center scenario, this experiment explores the impact of different participation and compensation levels of data center demand response. The resulting scenario names are: (1) 5% demand available for curtailment with 750 $/MWh compensation; (2) 5% demand available for curtailment with 500 $/MWh compensation; (3) 5% demand available for curtailment with 250 $/MWh compensation; (4) 15% demand available for curtailment with 750 $/MWh compensation; (5) 15% demand available for curtailment with 500 $/MWh compensation; and (6) 15% demand available for curtailment with 250 $/MWh compensation. The related folder names are: dr_cost_250_drup_0_drdown_5, dr_cost_250_drup_0_drdown_15, dr_cost_500_drup_0_drdown_5, dr_cost_500_drup_0_drdown_15, dr_cost_750_drup_0_drdown_5, and dr_cost_750_drup_0_drdown_15. Combination of delayed generator retirements and demand response through curtailment sensitivity experiment: The impact of combining postponing 100% nuclear and 25% natural gas retirements with 5% demand available for curtailment with 750 $/MWh compensation is simulated. The related folder is named "dr_cost_750_drup_0_drdown_5_nuc_100_gas_25". Please refer to the README file for a detailed description of the dataset including individual files and references.

Artificial Intelligence↗