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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
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↗
DOE OSTI2025

Neuromorphic heat transport effects in a molecular junction

Understanding energy transport at the nanoscale is an open and fundamental challenge in the molecular sciences with direct implications for the design of new electronics, computing devices, and materials. While nanoscale energy transport under steady-state conditions has been studied extensively, there is much less known about energy transport under time-dependent driving forces, particularly in the far-from-equilibrium regime. In this work, we use nonequilibrium molecular dynamics simulations and stochastic thermodynamics to investigate energy transport in a well-studied nanoscale system—a molecular junction—subjected to a time-periodic temperature gradient. The primary observation is that molecular junctions can exhibit heat transport hysteresis, a phenomenon in which the heat flux through a system depends not only on the instantaneous value of a time-dependent temperature bias but also on the temporal history of that bias. The presented findings illustrate that molecular junctions can exhibit the specific memory effect—heat transport hysteresis—that is essential for the design of thermal neuromorphic computers. This work elucidates a potential pathway toward the realization of such devices.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

On the engulfment of antifreeze proteins by ice

Antifreeze proteins (AFPs) are remarkable biomolecules that suppress ice formation at trace concentrations. To inhibit ice growth, AFPs must not only bind to ice crystals, but also resist engulfment by ice. The highest supercooling, ΔT*, for which AFPs are able to resist engulfment is widely believed to scale as the inverse of the separation, L, between bound AFPs, whereas its dependence on the molecular characteristics of the AFP remains poorly understood. By using specialized molecular simulations and interfacial thermodynamics, here, we show that in contrast with conventional wisdom, ΔT* scales as L –2 and not as L –1 . We further show that ΔT* is proportional to AFP size and that diverse naturally occurring AFPs are optimal at resisting engulfment by ice. By facilitating the development of AFP structure–function relationships, we hope that our findings will pave the way for the rational design of AFPs.

Science & Technology - Other Topics↗
DOE OSTI2022

Single- and two-particle finite size effects in interacting lattice systems

Simulations of extended quantum systems are typically performed by extrapolating results of a sequence of finite-system-size simulations to the thermodynamic limit. In the quantum Monte Carlo community, twist-averaging was pioneered as an efficient strategy to eliminate one-body finite size effects. In the dynamical mean field community, cluster generalizations of the dynamical mean field theory were formulated to study systems with nonlocal correlations. In this work, we put the twist-averaging and the dynamical cluster approximation variant of the dynamical mean field theory onto equal footing, discuss commonalities and differences, and compare results from both techniques to the standard periodic boundary technique. At the example of Hubbard-type models with local, short-range and Yukawa-like longer range interactions we show that all methods converge to the same limit, but that the convergence speed differs in practice. We show that embedding theories are an effective tool for managing both one-body and two-body finite size effects, in particular if interactions are averaged over twist angles.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Dipole condensates in tilted Bose-Hubbard chains

Here we study the quantum phase diagram of a Bose-Hubbard chain whose dynamics conserves both a boson number and boson dipole moment, a situation which can arise in strongly tilted optical lattices. The conservation of the dipole moment has a dramatic effect on the phase diagram, which we analyze by combining a field theory analysis with DMRG simulations. In the thermodynamic limit, the phase diagram is dominated by various types of incompressible dipolar condensates. In finite-sized systems, however, it may be possible to stabilize a Bose-Einstein insulator: an exotic compressible phase which is insulating, despite the absence of a charge gap. We suggest several ways by which these exotic phases can be identified in near-term cold-atom experiments.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Studying Cloud and Radiative Impacts Through an Improved Physically Based Representation of Organic Aerosol in a Large-Scale Model (WRF-Chem) (Final Report)

The work performed under this award resulted in the development of a computationally-efficient kinetic model to simulate the chemistry, thermodynamics, and microphysics of organic aerosol. Versions of the kinetic model were used to study a suite of aerosol processes including, multigenerational gas-phase chemistry, phase-state influenced gas/particle partitioning, heterogeneous oxidation, oligomerization reactions, and experimental artifacts. The kinetic model has now been implemented into a regional chemistry-climate model (WRF-Chem) to further study the formation, evolution, properties, and spatiotemporal distribution of organic aerosol at urban, regional, and global scales.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Multimetallic Layered Composites (MMLCs) for Rapid, Economical Advanced Reactor Deployment (Final Report)

This project focused on the development of multi-metallic layered composites (MMLCs) for advanced fission reactor technologies. There are many instances where one alloy or material simply cannot meet all the demands thrown at it by a reactor system, or cannot allow it to perform as strongly as one would like. Instead of focusing all our effort on developing one perfect alloy, we seek to leverage the design principle of “separation of functionality,” used in many other arenas in design, to boost performance beyond single alloys alone. One illustrative example shows the power of this approach for molten salt-cooled reactors: A three meter tall, three meter diameter reactor vessel made of Incoloy 800 was quoted at $\$$500k in 2018. A Hastelloy N vessel was quoted at $\$$5M. An MMLC vessel, in which a layer of Hastelloy N would be weld-overlaid onto Incoloy 800, was quoted at $\$$700k, and it would achieve the same performance. The potential economic gains of leveraging this approach are therefore substantial. At a minimum, each MMLC would contain one core structural layer and one coolant-facing corrosion-resistant layer. Sometimes, MMLCs required buffer layers, as the structural and corrosion-resistant layers were metallurgically incompatible. In other words, they didn’t always play nice, thus separating layers compatible with both functioned as intermediaries to keep the composite together. However, in doing so we inevitably produce new interfaces, where new issues can arise. Therefore, this project focused on what happens at these interfaces from a combination of high temperatures, irradiation, corrosion, and time. After all, a reactor makes money when it is operating, and outages of any kind erode its economic viability. First, we set out to experimentally prove that MMLCs for at least two advanced reactor systems can be made, today, in US domestic facilities. In this respect we were successful – one MMLC (a Ni-201/Incoloy 800H composite) was successfully made and drawn into two-inch coolant piping. Others were attempted, though new issues relating to cracking in vanadium layers for one and radiation damage performance of the corrosion-resistant layer in another prevented us from moving further in those specific arenas – these are engineering problems which deserve continued focus after this project. Additional experimental work focused on long-term corrosion testing of the outermost layers of the salt-cooled and liquid lead-cooled MMLC concepts, which would then be fed into predictions of how long the MMLCs could last. Next, computational (thermodynamics and atomistic) simulation studies studied how much we expect the interfaces to “blend,” due to the mixing action of neutron irradiation. This eats into both the margin for the structural layer of each MMLC, as dilution from the corrosion-resistant layer into the structural layer would decrease the total load-bearing capacity of an MMLC of finite size. On the other hand, dilution of the corrosion-resistant layer into the structural layer further reduced the margin of corrodible material, reducing the lifetime of the MMLC or necessitating extra thickness to be imparted to the MMLC to meet its functional requirements. Work here focused on irradiation-induced segregation to predict new phases which may embrittle the MMLCs, as well as quantifying irradiation-induced mixing at each interface. The results showed that mixing is expected, but it is both steady and therefore predictable, and not lifetime-limiting for most MMLC concepts – it simply has to be accounted for in calculations of reactor performance when utilizing an MMLC. Then, full-core simulations using the experimentally-derived corrosion data, the computationally discovered irradiation-induced mixing data (partially validated by experiment), and existing, benchmarked core designs for large and small sized reactor concepts (one salt-cooled, one lead-cooled) were conducted to quantify any expansion of reactor operating envelopes achieved by utilizing these MMLCs. This new framework, called REX (Reactor Envelope Expansion), incorporates a combination of core neutronics, thermal hydraulics, and the material performance data derived from this project to see how using an MMLC expands advanced fission reactor operating envelopes. It was discovered that in some cases, MMLC utilization does indeed increase the maximum operating temperatures and cycle lengths of reactor concepts, while in other cases it does not. Finally, our tech-to-market (T2M) strategy was not necessarily to create specific embodiments of MMLCs for immediate sale (because getting into the nuclear market is incredibly slow and laden with regulation, this is a long-term goal), but rather immediate stimulation of US industry using the design approach of MMLCs derived from this project. In this respect we were successful, as one of the PhD students funded on this project co-founded Allium Engineering, Inc., which created a stainless steel / low-alloy steel MMLC to function as chloride corrosion-resistant rebar for embedding into concrete structures. Allium Engineering continues to be successful, having recently opened their first factory as of this writing.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Manufacturing of Al-Ce/Metal-Matrix-Composites (MMC) for Harsh Environments

An integrated modeling/simulation and experimental approach will be used to design and synthesize Al-Ce/MMC for use in key technologies where improved energy efficiency and performance is harsh environments is required. Modeling and simulation work including thermodynamic modeling, DFT calculations, and continuum modeling of casting processes and mechanical behavior, will be combined with lab-scale synthesis to identify potential matrix and reinforcing phase compositions for large-scale synthesis. The experimental synthesis work will rely on lab-scale casting and gas atomization to identify the optimal matrix compositions for pilot/commercial-scale synthesis is Al-Ce/MMC. Selected Al-Ce/MMC will be synthesized by commercial scale melt processing, powder metallurgy and thermomechanical processing. The corrosion resistance, mechanical properties as a function of temperature, and thermal conductivity of the composites will be characterized to validate target properties are met or exceeded. This work will focus on two selected applications; however, the design framework will have broad applicability for composites for numerous harsh service condition applications. We will have synthesized Al-Ce/MMC that exhibit better performance in harsh environments than currently used materials.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Achieving High Strength and Creep Resistance in Inconel 740H Superalloy through Wire-Arc Additive Manufacturing and Thermodynamic-Guided Heat Treatment

Inconel 740H superalloy is commonly used in advanced ultra-supercritical power plants since it possesses excellent strength and creep resistance. This study investigates the microstructure and mechanical properties of Inconel 740H superalloy fabricated using wire-arc additive manufacturing. The as-printed microstructure consisted of columnar γ grains with the Laves phase and (Nb, Ti)C carbides as secondary phases. The anisotropy in grain structure increased from the bottom to the top regions, while the hardness was highest in the middle portion of the build. To guide the post-heat treatment design, thermodynamic and kinetic simulations were employed to predict the temperature and time. Complete recrystallization with the Laves phase dissolution occurred throughout the build after homogenization at 1200 °C for 2 h. The peak hardness was achieved after aging at 760 °C for 12 h with the M 23 C 6 carbides decorating the grain boundaries and γ’ precipitates in the grain interior. The yield strength (655 MPa) and ductility (29.5%) in the post-heat treated condition exceeded the design targets (620 MPa, 20%). Stress rupture tests at 750 °C showed that the high-temperature performance was at par with the wrought counterparts. The fracture mode after rupture was identified to be intergranular with the presence of grain boundary cavities along with grain boundary sliding.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Statistical Mechanics on Lattices (smol) v0.0.1

smol is a Python implementation of the Cluster Expansion Method and extensions of this methodology; which is used to fit applied lattice models from first principle calculations for subsequent use in Monte Carlo simulations for computing thermodynamic and statistical properties of atomic configuration. smol has a several notable advantages over other similar software packages. The first being its Python implementation which makes it easy to use even for users with little software and coding skills. Despite being implemented in Python, essential Monte Carlo routines are implemented in Cython such that performance is not compromised and is competitive even to similar software implemented in C/C++. Furthermore, the package has specific functionality for handling complex ionic materials (such as cathodes and electrolytes) that is not readily available in other packages. Lastly, the software has a flexible and modular design with the intention of making it fast and efficient to develop and extend the methodology.

Barroso-Luque, Luis↗
DOE OSTI2024

Lambda-PFLOTRAN 1.0: a workflow for incorporating organic matter chemistry informed by ultra high resolution mass spectrometry into biogeochemical modeling

Abstract. Organic matter (OM) composition plays a central role in microbial respiration of dissolved organic matter and subsequent biogeochemical reactions. Here, a direct connection of organic matter chemistry and thermodynamics to reactive transport simulators has been achieved through the newly developed Lambda-PFLOTRAN workflow tool that succinctly incorporates carbon chemistry data generated from Fourier transform ion cyclotron resonance mass spectrometry (FTICR-MS) into reaction networks to simulate organic matter degradation and the resulting biogeochemistry. Lambda-PFLOTRAN is a Python-based workflow, executed through a Jupyter notebook interface, that digests raw FTICR-MS data, develops a representative reaction network based on substrate-explicit thermodynamic modeling (also termed lambda modeling due to its key thermodynamic parameter λ used therein), and completes a biogeochemical simulation with the open source, reactive flow and transport code PFLOTRAN. The workflow consists of the following five steps: configuration, thermodynamic (lambda) analysis, sensitivity analysis, parameter estimation, and simulation output and visualization. Two test cases are provided to demonstrate the functionality of the Lambda-PFLOTRAN workflow. The first test case uses laboratory incubation data of temporal oxygen depletion to fit lambda parameters (i.e., maximum utilization rate and microbial carrying capacity). A slightly more complex second test case fits multiple lambda formulation and soil organic matter release parameters to temporal greenhouse gas generation measured during a soil incubation. Overall, the Lambda-PFLOTRAN workflow facilitates upscaling by using molecular-scale characterization to inform biogeochemical processes occurring at larger scales.

58 GEOSCIENCES↗