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 469 records · Page 26

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Estimating ionization states and continuum lowering from ab initio path integral Monte Carlo simulations for warm dense hydrogen

Warm dense matter (WDM) is an active field of research, with applications ranging from astrophysics to inertial confinement fusion. Ionization degree and continuum lowering are important quantities to understand how materials behave under these conditions, but can be difficult to diagnose since experimental campaigns are limited and often require model-dependent analysis. This is especially true for hydrogen, which has a comparably low scattering cross section, making high-quality data particularly difficult to obtain. Consequently, building equation of state tables often relies on simulations in combination with untested approximations to extract properties from experiments. Here, we investigate an approach for extracting the ionization potential depression and ionization degree—quantities which are otherwise not directly accessible from the physical model—from first-principles path integral Monte Carlo (PIMC) simulations utilizing a chemical model. In contrast to experimental measurements, where noise and nonequilibrium effects add to the uncertainty of the inferred parameters, PIMC simulations provide a clean signal with well-defined thermodynamic conditions. Comparisons against commonly used models show a qualitative agreement, but we find deviations primarily for the high-density and high-temperature cases. We also demonstrate the decreasing sensitivity of the dynamic structure factor with respect to both ionization and continuum lowering for increasing scattering angles in x-ray Thomson scattering experiments. Our work has important implications for the design of future experiments, but also offers qualitative understanding of structure factors and the imaginary-time correlation function obtained from first-principles quantum Monte Carlo simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRSDate not supplied

Comparisons of CFD Simulations of Icing Wind Tunnel Clouds with Experiments Conducted at the NASA Propulsions Systems Laboratory

This paper evaluates simulation predictions against experimental test data of icing clouds that were produced during 2018 ice crystal icing physics tests conducted at the NASA Propulsion Systems Laboratory icing wind tunnel. Aero-thermal and cloud parameters are set and known upstream at the tunnel inlet and spray system, but change as the cloud and air thermodynamically interact as the flowing masses reach the tunnel test section. Utilizing the ANSYS Fluent Discrete Phase Model function, 3D computational fluid dynamics (CFD) simulations were performed, capturing the thermodynamic interactions between the test parameters, and providing predictions of the aero thermal and cloud conditions at the tunnel test section. Simulation predictions were compared with test data that were measured at the tunnel exit plane. Evaluations focused on the cloud concentration (total water content), humidity content, and air temperature. CFD simulation predictions showed areas of agreement and disagreement. Simulations showed that cloud concentration profiles at the test section are strongly related to the initial spray nozzle pattern used at the tunnel inlet. Experimental data suggest that greater dispersion of the cloud occurred as the simulated cloud predicted areas of high and low cloud concentration compared to test data profiles. Simulations, however, captured the magnitude and location of the change in humidity content and the change in air temperature due to the presence of the cloud reasonably well, when compared to test data. This result would suggest that while the ANSYS simulation did not fully predict the spreading of the cloud as measured during experiment, it did capture evaporation and the molecular movements of air and vapor relatively well.

CFD↗

NASA NTRSDate not supplied

Comparisons of CFD Simulations of Icing Wind Tunnel Clouds with Experiments Conducted at the NASA Propulsions Systems Laboratory

CFD↗

DOE OSTI2021

Vapor–liquid equilibrium of water with the MB-pol many-body potential

Among the many existing molecular models of water, the MB-pol many-body potential has emerged as a remarkably accurate model, capable of reproducing thermodynamic, structural, and dynamic properties across water's solid, liquid, and vapor phases. In this work, we assessed the performance of MB-pol with respect to an important set of properties related to vapor-liquid coexistence and interfacial behavior. Through direct coexistence classical molecular dynamics simulations at temperatures of 400 K < T < 600 K, we calculated properties such as equilibrium coexistence densities, vapor-liquid interfacial tension, vapor pressure, and enthalpy of vaporization and compared the MB-pol results to experimental data. We also compared rigid vs fully flexible variants of the MB-pol model and evaluated system size effects for the properties studied. We found that the MB-pol model predictions are in good agreement with experimental data, even for temperatures approaching the vapor-liquid critical point; this agreement was largely insensitive to system sizes or the rigid vs flexible treatment of the intramolecular degrees of freedom. Furthermore, these results attest to the chemical accuracy of MB-pol and its high degree of transferability, thus enabling MB-pol's application across a large swath of water's phase diagram.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Machine learning approaches for structural and thermodynamic properties of a Lennard-Jones fluid

Predicting the functional properties of many molecular systems relies on understanding how atomistic interactions give rise to macroscale observables. However, current attempts to develop predictive models for the structural and thermodynamic properties of condensed-phase systems often rely on extensive parameter fitting to empirically selected functional forms whose effectiveness is limited to a narrow range of physical conditions. Here, we illustrate how these traditional fitting paradigms can be superseded using machine learning. Specifically, we use the results of molecular dynamics simulations to train machine learning protocols that are able to produce the radial distribution function, pressure, and internal energy of a Lennard-Jones fluid with increased accuracy in comparison to previous theoretical methods. The radial distribution function is determined using a variant of the segmented linear regression with the multivariate function decomposition approach developed by Craven et al. [J. Phys. Chem. Lett. 11, 4372 (2020)]. The pressure and internal energy are determined using expressions containing the learned radial distribution function and also a kernel ridge regression process that is trained directly on thermodynamic properties measured in simulation. The presented results suggest that the structural and thermodynamic properties of fluids may be determined more accurately through machine learning than through human-guided functional forms.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Development of Future Heatwaves for Different Hazard Thresholds

In 2018 and 2019, heatwaves set all-time temperature records around the world and caused adverse effects on human health, agriculture, natural ecosystems, and infrastructure. Often, severe impacts relate to the joint spatial and temporal extent of the heatwaves, but most research so far focuses either on spatial or temporal attributes of heatwaves. Furthermore, sensitivity of heatwaves characteristics to the choice of the heatwave thresholds in a warming climate are rarely discussed. Here, we analyze the largest spatiotemporal moderate heatwaves—that is, three-dimensional (space-time) clusters of hot days—in simulations of global climate models. We use three different hazard thresholds to define a hot day: fixed thresholds (time-invariant climatological thresholds), seasonally moving thresholds based on changes in the summer means, and fully moving thresholds (hot days defined relative to the future climatology). We find a substantial increase of spatiotemporally contiguous moderate heatwaves with global warming using fixed thresholds, whereas changes for the other two hazard thresholds are much less pronounced. In particular, no or very little changes in the overall magnitude, spatial extent, and duration are detected when heatwaves are defined relative to the future climatology using a temporally fully moving threshold. This suggests a dominant contribution of thermodynamic compared to dynamic effects in global climate model simulations. The similarity between seasonally moving and fully moving thresholds indicates that seasonal mean warming alone can explain large parts of the warming of extremes. The strong sensitivity of simulated future heatwaves to hazard thresholds should be considered in the projections of potential future heat-related impacts.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

The dynamics, mixing, and thermonuclear burn of compressed foams with varied gas fills

Inertial confinement fusion (ICF) implosions involve highly coupled physics and complex hydrodynamics that are challenging to model computationally. Due to the sensitivity of such implosions to small features, detailed simulations require accurate accounting of the geometry and dimensionality of the initial conditions, including capsule defects and engineering features such as fill tubes used to insert gas into the capsule, yet this is computationally prohibitive. It is therefore difficult to evaluate whether discrepancies between the simulation and experiment arise from inadequate fidelity to the capsule geometry and drive conditions, uncertainties in physical data used by simulations, or inadequate physics. We present results from detailed high-resolution three-dimensional simulations of ICF implosions performed as part of the MARBLE campaign on the National Ignition Facility [Albright et al., Phys. Plasmas 29, 022702 (2022)]. These experiments are foam-filled separated-reactant experiments, where deuterons reside in the foam and tritons reside in the capsule gas fill and deuterium–tritium (DT) fusion reactions only occur in the presence of mixing between these materials. Material mixing in these experiments is primarily seeded by shock interaction with the complex geometry of the foam and gas fill, which induces the Richtmyer–Meshkov instability. We compare results for experiments with two different gas fills (ArT and HT), which lead to significant differences in the hydrodynamic and thermodynamic developments of the materials in the implosion. Our simulation results show generally good agreement with experiments and demonstrate a substantial impact of hydrodynamic flows on measured ion temperatures. The results suggest that viscosity, which was not included in our simulations, is the most important unmodeled physics and qualitatively explains the few discrepancies between the simulation and experiment. The results also suggest that the hydrodynamic treatment of shocks is inadequate to predict the heating and yield produced during shock flash, when the shock converges at the center of the implosion. Alternatively, underestimation of the level of radiative preheat from the shock front could explain many of the differences between the experiment and simulation. Nevertheless, simulations are able to reproduce many experimental observables within the level of experimental reproducibility, including most yields, time-resolved X-ray self-emission images, and an increase in burn-weighted ion temperature and neutron down-scattered ratio in the line of sight that includes a jet seeded by the glue spot that joins capsule hemispheres.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Hydrogen Localization and Cluster Formation in α-Zr from First-Principles Investigations

We report alpha-zirconium (α-Zr) is a hydrogen atom host material used in nuclear energy and semiconducting applications due to Zr's excellent properties such as corrosion resistance, low thermal neutrons cross-section, and thermodynamically favorable interstitial sequestering of H atoms. Herein we use density functional theory simulations to systematically analyze the initial stages of H insertion in the Zr matrix, and determine optimal hydrogen localization, distribution, and clustering tendencies from a thermodynamics and kinetics perspective. Our results suggest that H occupation of tetrahedral sites is strongly favored, consistent with previous experimental and theoretical studies. Further, a cooperative effect in H localization is shown where H atoms favor clustering. The H pair is found to be the most stable by occupying tetrahedral sites with pairs at second nearest-neighbor positions (TTn2). Among the four H atom clusters, the most favorable configuration has two TTn2 pairs forming "box " or "zig-zag" shaped-cluster with the pairs separated by 3.03 or 3.25 Å, respectively. Nudged elastic band simulations suggest that the formation of TTn2 is facile and is limited by the small diffusion barrier of the monomer of less than 0.1 eV. Further, the formation routes for the 4H atom clusters are also presented. The current study can serve as a foundation to systematically model H atom clustering in metals and alloys.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Use of Remote Sensing and In-Situ Observations to Develop and Evaluate Improved Representations of Convection and Clouds for the ACME Model

The overachieving goal of the whole CMDV-MCS project is to improve understanding of warm season continental convection and to develop treatments of convection and microphysics capable of representing mesoscale convective systems (MCSs) features in large-scale models. Our tasks for this project contributing to the overachieving goal include: (1) Improve the ice nucleation formulation for MG2 and P3 cloud microphysics schemes; (2) Improve the treatment of subgrid dynamics and thermodynamics driving the ice nucleation in E3SM; and (3) Test the performance of improved ice microphysics in E3SM with observation data. In this project, we have (1) Improved the ice nucleation parameterization for MG2 and P3 in E3SM by implementing two advanced empirical parameterizations with connection to aerosols. The two deterministic heterogeneous ice nucleation parameterizations (i.e., DeMott et al., 2015; Niemand et al., 2012) were merged with the MG2 and P3 cloud microphysics schemes in E3SM. Long-term simulations were conducted to examine the impacts of these new parameterizations on simulated cloud properties; (2) Improved the treatment of subgrid dynamics and thermodynamics driving the ice nucleation in E3SM. We evaluated the double Gaussian PDF of vertical velocity simulated by the Cloud Layers Unified By Binormals (CLUBB) and the sub-column vertical velocity sampled from the Subgrid Importance Latin Hypercube Sampler (SILHS) in E3SM. We introduced the vertical velocity variance induced by topographic gravity waves for ice nucleation and droplet activation; and (3) Tested the performance of improved ice microphysics in E3SM with observation data. We tested the new treatments of ice nucleation in the single column model (SCM) mode for the stratiform mixed-phase clouds observed during 9-10 October 2004 in the DOE ARM Mixed-Phase Arctic Cloud Experiment (M-PACE) and for the convective clouds observed on 20 May 2011 in the Midlatitude Continental Convective Clouds Experiment (MC3E). Modeled ice nucleating particles (INPs) concentrations were compared against observations collected around the globe.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Towards the Accurate Thermodynamic Characterization of Enzyme Reaction Mechanisms

We employed QM/MM molecular dynamics (MD) simulations to characterize the rate-limiting step of the glycosylation reaction of pancreatic α-amylase with combined DFT/molecular dynamics methods (PBE/def2-SVP : AMBER). Upon careful choice of four starting active site conformations based on thorough reactivity criteria, Gibbs energy profiles were calculated with umbrella sampling simulations within a statistical convergence of 1 -2 kcal · mol -1 . Nevertheless, Gibbs activation barriers and reaction energies still varied from 11.0 to 16.8 kcal · mol -1 and -6.3 to +3.8 kcal · mol -1 depending on the starting conformations, showing that despite significant state-of-the-art QM/MM MD sampling (0.5 ns/profile) the result still depends on the starting structure. In conclusion, the results supported the one step dissociative mechanism of Asp197 glycosylation preceded by an acid-base reaction by the Glu233, which are qualitatively similar to those from multi-PES QM/MM studies, and thus support the use of the latter to determine enzyme reaction mechanisms.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Assessment of Precipitating Marine Stratocumulus Clouds in the E3SMv1 Atmosphere Model: A Case Study from the ARM MAGIC Field Campaign

This paper presents a process-oriented evaluation of precipitating stratocumulus and its transition to cumulus in version 1 of the Energy Exascale Earth System Model (E3SMv1) using comprehensive case-study observations from a field campaign of the Atmospheric Radiation Measurement program (ARM). The E3SMv1 single-column model (SCM) of the marine boundary layer and its low clouds and precipitation are compared to observations including subcloud drizzle retrievals from a combination of Doppler radar and lidar backscatter measurements. Additionally, the SCM is also compared to a large-eddy simulation (LES) of the same case. The combination of advanced remote sensing observations and LES is a powerful framework to evaluate the physical parameterizations of large-scale models. Given the observed large-scale environment, the E3SMv1 SCM realistically represents the evolution of clouds and boundary layer structure during the stratocumulus-to-cumulus transition. The model well simulates the liquid water path and its diurnal cycle in the stratocumulus period as well as the two-layer vertical thermodynamic structure and lower cloud fraction in the transition period. E3SMv1’s success in simulating the cloud in the stratocumulus period permitted examination of its precipitation processes. Here problems were identified with E3SMv1 producing an unrealistically small subcloud precipitation fraction, an unrealistic double peak in the vertical profiles of precipitation mass, and drizzle that evaporates too close to the surface. Further model diagnostics determined that these unrealistic characteristics resulted from an overly long microphysics time step and an unrealistic parameterization of the precipitation fraction. These results imply that careful consideration of these issues is needed in order to better simulate precipitation processes in marine stratocumulus.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1985

Simulation of mesoscale convective response

Results are reported from a transfer of thermodynamic and dynamic data downscale from a two-dimensional mesoscale model to a two-dimensional cloud model and a subsequent examination of the differential convective response of the cloud model to a mesoscale structure. The mesoscale model included a high resolution PBL formulation, with convective forces expressed in a profile of exchange coefficients over the height of the PBL and the magnitude of the surface heat flux. Account was also taken of forcing by long- and short-wave radiation, surface forcing by the surface energy budget, the terrain shape, and possible wave reflection at the top boundary. The slab-symmetric cloud model possessed subgrid-scale features, five types of moisture terms, and a subroutine for accretion processes. The mesoscale environment was observed to have a significant impact on convective response, i.e., the vertical velocity and the cloud water. Various factors which were not included in the study and which must be considered are discussed.

Mcnider, R. T.↗

NASA NTRS2020

Experimental Validation of Nozzle Flow Simulations for Rotating Detonation Rocket Engines

Rotating detonation engines (RDEs) promise increased thermodynamic performance that may significantly enhance the capabilities of current rocket platforms. Little work has been conducted thus far to characterize the effect of nozzle design on high chamber pressure RDEs. Previous computational work was completed to understand the nozzle performance of the Purdue methane/oxygen rocket RDE. A new experimental study on a similar kerosene/oxygen RDE was conducted to validate the results of the computational study. Both a nozzleless geometry and several aerospike designs were hot-fire tested. New pressure instrumentation on the different nozzle surfaces were included to better understand the flow physics unique to RDE chamber exit conditions. Both the previous computational study and the new experimental results confirmed that for a nozzleless geometry, the RDE cycle enhances suction on the base region, an important result to determining RDE engine performance separate from nozzle effects. While the aerospike experiments showed agreement with computational results, delay of flow separation due to the RDE cycle could not be confirmed. Two aerospike geometries with different nozzle pressure ratios were experimentally evaluated. This study showed both the strength and the necessity for 3D transient computations to better understand the RDE flow field with nozzle geometries.

Alexis J Harroun↗

DOE OSTI2022

Thermodynamic modeling of adsorption at the liquid-solid interface

Adsorptive separation techniques are significantly energy efficient in comparison to conventional thermal separation techniques such as distillation. Despite extensive research and development activities undertaken for mixed-gas adsorption, the use of adsorption techniques for the separation of multicomponent liquid mixtures is still limited. A major barrier is the lack of accurate adsorption thermodynamic models, which form the scientific foundation of process simulation of such systems, making the translation to industrial scale challenging. In this work, we have rigorously computed the surface excess of adsorption for six binary liquid mixtures on silica gel at 303 K using the frameworks of the generalized Langmuir isotherm (gL) and the adsorbed solution theory (AST). The six binary liquid mixtures were formed by the pair-wise combinations of four components: benzene, 1,2-dichloroethane, cyclohexane, and n-heptane. We have based our calculations by considering simultaneous equilibria of three phases: saturated vapor phase, bulk liquid phase, and adsorbed phase. Further, the composition of the corresponding saturated vapor phase was estimated by the Nonrandom Two-Liquid activity coefficient model and experimental vapor-liquid equilibria data. The activity coefficients of the adsorbed phase, the central issue of multicomponent adsorption thermodynamics, were calculated using the adsorption Nonrandom Two-Liquid activity coefficient model. Devoid of simplifying assumptions, gL and AST provide rigorous thermodynamic frameworks for adsorption equilibria of multicomponent liquid mixtures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Comparative Energetics of Various Membrane Distillation Configurations and Guidelines for Design and Operation

This paper presents a comparative performance study of single-stage desalination processes with major configurations of membrane distillation (MD) modules. MD modules covered in this study are (a) direct contact MD (DCMD), (b) vacuum MD (VMD), (c) sweeping gas MD (SGMD), and (d) air gap MD (AGMD). MD-based desalination processes are simulated with rigorous theoretical MD models supported by molecular thermodynamic property models for the accurate calculation of performance metrics. The performance metrics considered in MD systems are permeate flux and energy efficiency, i.e., gained output ratio (GOR). A general criterion is established to determine the critical length of these four MDs (at fixed width) for the feasible operation of desalination in a wide range of feed salinities. The length of DCMD and VMD is restricted by the feed salinity and permeate flux, respectively, while relatively large AGMD and SGMD are allowed. The sensitivity of GOR flux with respect to permeate conditions is investigated for different MD configurations. AGMD outperforms other configurations in terms of energy efficiency, while VMD reveals the highest permeate production. With larger MD modules, utilization of thermal energy supplied by the hot feed for evaporation is in the order of VMD > AGMD > SGMD > DCMD. Simulation results highlight that energy efficiency of the overall desalination process relies on the efficient recovery of spent for evaporation, suggesting potential improvement in energy efficiency for VMD-based desalination.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Using a Coarse-Grained Modeling Framework to Identify Oligomeric Motifs with Tunable Secondary Structure

Coarse-grained modeling can be used to explore general theories that are independent of specific chemical detail. In this paper, we present cg_openmm, a Python-based simulation framework for modeling coarse-grained hetero-oligomers and screening them for structural and thermodynamic characteristics of cooperative secondary structures. cg_openmm facilitates the building of coarse-grained topology and random starting configurations, setup of GPU-accelerated replica exchange molecular dynamics simulations with the OpenMM software package, and features a suite of postprocessing thermodynamic and structural analysis tools. In particular, native contact analysis, heat capacity calculations, and free energy of folding calculations are used to identify and characterize cooperative folding transitions and stable secondary structures. In this work, we demonstrate the capabilities of cg_openmm on a simple 1–1 Lennard-Jones coarse-grained model, in which each residue contains 1 backbone and 1 side-chain bead. By scanning both nonbonded and bonded force-field parameter spaces at the coarse-grained level, we identify and characterize sets of parameters which result in the formation of stable helices through cooperative folding transitions. Furthermore, we show that the geometries and stabilities of these helices can be tuned by manipulating the force-field parameters.

36 MATERIALS SCIENCE↗

NASA NTRS1987

Nonequilibrium effects for hypersonic transitional flows

Presented are the results of numerical simulations of hypersonic flow about blunt cones and hemispherical nose configurations for reentry velocities of 7.5 and 10 km/s. Cone half angles 0, 5, and 10 deg are considered at zero angle of incidence; however, the focus is for the 5 deg cone. The body size and altitude ranges considered (70 to 110 km) are such that the flow is in the transitional regime. Translational, thermodynamic, and chemical nonequilibrium effects are considered in the numerical simulation by utilizing the direct simulation Monte Carlo (DSMC) method of Bird. The DSMC results are compared with those obtained with viscous shock-layer and Navier-Stokes methods. Comparisons between the DSMC and continuum calculations show the altitude range where differences in flowfield structure and surface quantities become significant. The current calculations show that the binary scaling similitude provides a means of correlating the blunt body surface quantities in the hypersonic, transitional regime. Furthermore, for the higher velocity entry conditions, the results highlight some of the concerns in the application of multitemperature continuum formulations, particularly the use of some proposed functional relations for the chemical rate constants under thermodynamic nonequilibrium conditions.

Moss, James N.↗

NASA NTRSDate not supplied

Computational Fluid Dynamics Simulation of Methane Slosh and Drain Experiments

NASA possesses a wealth of historical cryogenic experiments that provide valuable insights into the design, troubleshooting, and understanding involved in the complex fluid and thermodynamics of managing cryogenic propellants. One approach to leveraging these historical datasets is by simulating these experiments to validate the accuracy of simulation environments and constituent models. This study focused on simulating a selection of the K site test series for both static and sloshing pressurized liquid methane draining experiments conducted at NASA in the 1970’s, utilizing computational fluid dynamics. The simulations were performed in the ANSYS FLUENT environment using the Volume of Fluid (VOF) numerical approach. A k-omega turbulence model was used with interfacial turbulence damping, and accurately predicted the amount of pressurant needed to maintain the required tank pressure throughout the static drain. The static simulation predicted the temperature stratification in the ullage observed at the end of the drain. During the methane expulsion with sloshing test, many features were successfully captured using the k-omega turbulence model with interfacial turbulence damping included. The rate of phase change and liquid temperature was overpredicted compared to the experimental measurements. The overprediction may be attributed to uncertainties in the vessel geometry, methane pressurant temperature and composition, and methodological differences in how the sloshing frequency was adjusted during the expulsion.

Cryogenic Propellants↗