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 253 records · Page 14

Optimization and Evaluation of Stochastic Unified Convection Using Single-Column Model Simulations at Multiple Observation Sites

We extend the previously developed stochastic unified convection scheme (UNICON) for shallow convection to deep convection by parameterizing the impact of mesoscale organized flow on updraft properties. The extended stochastic UNICON parameterizes thermodynamic properties of updrafts at the near-surface as a multivariate Gaussian distribution, where the variances of the distribution are the summation of variances from non-organized turbulence and mesoscale organized flow. The distribution of updraft radius is parameterized as a power-law distribution with a scale break which is parameterized as a linear function of the strength of mesoscale organized flow. The proposed parameterization is validated using a series of large-eddy simulations of deep convection. The free parameters introduced in the formulation of stochastic UNICON are optimized using 10 cases of single-column model simulations over the ocean. Stochastic UNICON with the optimized parameters significantly reduces the biases of thermodynamic profiles and surface precipitation rates simulated in the original UNICON for tropical convection cases. The simulation of the variation in anomalies of temperature and moisture associated with the Madden-Julian oscillation is also improved. The overall improvements in simulated thermodynamic profiles are found to be due to the increased heating and drying tendencies by convective processes in stochastic UNICON. An additional simulation of an idealized deep convection case shows that stochastic UNICON produces enhanced cloud variabilities with dependency on updraft radius, indicating its ability to represent the coexistence of shallow and deep convection.

54 ENVIRONMENTAL SCIENCES↗

Thermomechanical response of thickly tamped targets and diamond anvil cells under pulsed hard x-ray irradiation

In the laboratory study of extreme conditions of temperature and density, the exposure of matter to high intensity radiation sources has been of central importance. Here, we interrogate the performance of multi-layered targets in experiments involving high intensity, hard x-ray irradiation, motivated by the advent of extremely high brightness hard x-ray sources, such as free electron lasers and 4th-generation synchrotron facilities. Intense hard x-ray beams can deliver significant energy in targets having thick x-ray transparent layers (tampers) around samples of interest for the study of novel states of matter and materials’ dynamics. Heated-state lifetimes in such targets can approach the microsecond level, regardless of radiation pulse duration, enabling the exploration of conditions of local thermal and thermodynamic equilibrium at extreme temperature in solid density matter. The thermal and mechanical responses of such thick layered targets following x-ray heating, including hydrodynamic relaxation and heat flow on picosecond to millisecond timescales, are modeled using radiation hydrocode simulation, finite element analysis, and thermodynamic calculations. Assessing the potential for target survival over one or more exposures and resistance to damage arising from heating and resulting mechanical stresses, this study doubles as an investigation into the performance of diamond anvil high pressure cells under high x-ray fluences. Long used in conjunction with synchrotron x-ray radiation and high power optical lasers, the strong confinement afforded by such cells suggests novel applications at emerging high intensity x-ray facilities and new routes to studying thermodynamic equilibrium states of warm, very dense matter.

Meza-Galvez, J. (ORCID:000000032699849X)↗

Simultaneous inference of equation of state parameters and unknown data errors with uncertainty quantification via hierarchical Bayesian posterior maximization

Equations of state (EOSs) are a key component in running hydrodynamic simulations as they relate the thermodynamic states for the material. The Davis reactants EOS is commonly used for modeling high explosives (HEs), and the EOS model parameters are calibrated using material specific data. The calibrations are often performed with uncertainty quantification via Bayesian inference to account for uncertainty in the data and generate ensembles of likely parameters. However, there are relatively few HE data sets to use for calibration and many are historical and lack error information. In this work, we simultaneously calibrate the Davis reactants EOS model parameters and unknown data error terms for the high explosive PBX 9501. To quantify the uncertainty in the models and the data, we use a Bayesian framework for the calibration and compute the hierarchical Bayesian posterior distribution with both a posteriori maximization approach and Markov Chain Monte Carlo. In general, we find that, given our assumptions, the two approaches result in similar calibrated parameters, posterior covariance matrices, and insights about the parameters but that the posterior maximization requires far less computational resources.

97 MATHEMATICS AND COMPUTING↗

Generating Skeletal Chemical Reaction Mechanisms for Post-Detonation Flows

This report documents the generation of a skeletal chemical reaction mechanism for use with hemispherical pentaerythritol tetranitrate charges. Skeletal mechanisms can substantially reduce computation time while maintaining accuracy. The methodology within uses faster running sample simulations to build a representative thermodynamic state space. These thermodynamic states are used with a constant-volume reactor analysis and a reaction flow analysis to remove unimportant species and reactions from a full chemical reaction mechanism. For the given test case, this results in a 6x speedup in computation time for directly comparable simulations in 2D axisymmetric simulations. We see a 30x speedup in simulations in 3D Cartesian coordainates when compared to a prior full kinetics simulation. There is strong agreement between temperature and species mass fraction profiles between the full and skeletal chemical reaction mechanisms. These methodologies can be applied to any explosive, given the availability of sample simulations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Simulation of Delamination Under High Cycle Fatigue in Composite Materials Using Cohesive Models

A new thermodynamically consistent damage model is proposed for the simulation of high-cycle fatigue crack growth. The basis for the formulation is an interfacial degradation law that links Fracture Mechanics and Damage Mechanics to relate the evolution of the damage variable, d, with the crack growth rate da/dN. The damage state is a function of the loading conditions (R and (Delta)G) as well as the experimentally-determined crack growth rates for the material. The formulation ensures that the experimental results can be reproduced by the analysis without the need of additional adjustment parameters.

Camanho, Pedro P.↗

A Landau–Devonshire analysis of strain effects on ferroelectric Al 1-x Sc x N

We present a thermodynamic analysis of the recently discovered nitride ferroelectric materials using the classic Landau–Devonshire approach. Electrostrictive and dielectric stiffness coefficients of Al 1-x Sc x N with a wurtzite structure (6 mm) are determined using a free energy density function assuming a hexagonal parent phase (6/mmm), with the first-order phase transition based on the dielectric stiffness relationships. The results of this analysis show that the strain sensitivity of the energy barrier is one order of magnitude larger than that of the spontaneous polarization in these wurtzite ferroelectrics, yet both are less sensitive to strain compared to classic perovskite ferroelectrics. These analysis results reported here explain experimentally reported sensitivity of the coercive field to elastic strain/stress in Al 1-x Sc x N films and would enable further thermodynamic analysis via phase field simulation and related methods.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Kinematic and Thermodynamic Structures of Hurricane Bob (1991) Determined From A 1.3-KM Resolution Numerical Simulation

Hurricane Bob (1991) is simulated using the Penn State/NCAR mesoscale model MM5. The simulation is conducted for a 24-h period at 4-km resolution and for a 6-h period at 1.3-km resolution. The 4-km simulation is able to fairly realistically capture the intensity and structure of the storm. The 1.3-km simulation depicts very small-scale convective structures and produces a convective rain band outside of the eye wall that did not occur in the 4-km simulation. The 1.3-km results are used to characterize several kinematic structures in the storm, including low-level outflow in the eye wall, multiple outflows at upper levels, and the convective rain band structure outside of the eye wall. Thermodynamic characteristics of air parcels flowing into and rising within the eye wall will also be examined through trajectory calculations.

Braun, Scott A.↗

Accelerating Thermochemical Equilibrium Calculations for Nuclear Reactor Applications

Thermochemical properties play a key role in modeling and simulation of several key phenomena in nuclear reactors. There has been an increasing interest in incorporating CALPHAD-based formulations in multiphysics simulations including for Molten Salt Reactors where knowledge of phase evolution of the salt and the chemical potentials of various elements are of utmost importance in source term analyses and redox control. However, the size of such simulations is often limited by the high computational cost of full thermodynamic equilibrium calculations. This work discusses the current efforts aimed at accelerating thermochemical equilibrium calculations for multiphysics simulations performed using the open-source finite element / finite volume code Multiphysics Object Oriented Simulation Environment (MOOSE) [1]. While several methods have been proposed for accelerating phase equilibrium calculations [2], most focus on relatively small systems and often rely on a- priori knowledge of the state-space of the system. Nuclear materials, however, are often multi-component systems owing to the evolution of composition under irradiation and an approach based on a-priori mapping of phase diagram is often not enough. This work is aimed at demonstrating an on-the-fly surrogate modeling framework that uses active learning to reduce the number of full equilibrium calculations that must be performed. By combining with efficient coupling approaches, the surrogate framework helps in reducing the computational cost of thermodynamic equilibrium informed multiphysics simulations of nuclear materials. The performance is benchmarked against full coupling with the thermochemistry library Thermochimica [3]. This work uses a machine learning based approach for constructing surrogate models to predict the stable phases in a multicomponent system. The surrogates were constructed using neural networks and Gaussian process classification. In this work, we compare the relative performance of the two methods. We also demonstrate the use of caching previous calculations by interpolating the values from nearest neighbors. References [1] Lindsay, A.D., et al. "2.0 – MOOSE: Enabling massively parallel multiphysics simulation", SoftwareX, 20 (2022): 101202. [2] Roos, W.A. and Zietsman J.H. "Accelerating complex chemical equilibrium calculations – A Review", Calphad, 77 (2022): 102380. [3] Piro, M.H.A., et al. "The thermochemistry library Thermochimica", Computational Materials Science, 67 (2013): 266-272.

36 MATERIALS SCIENCE↗

Low global warming potential (GWP) refrigerant supermarket refrigeration system modeling and its application

As an environmentally friendly low global warming potential (GWP) refrigerant, Carbon dioxide (CO 2 ) has continuously gained popularity and research attention as alternative refrigerant for supermarket refrigeration system. In this paper, to fulfill the increasing need of accurate Low-GWP supermarket refrigeration models for development of supervisory level control and optimization strategies, a high fidelity model is developed for CO 2 transcritical supermarket refrigeration system which includes compressor rack of low temperature (LT) compressors and medium temperature (MT) compressors, air-cooled gas cooler, evaporator, expansion valves and other auxiliary equipment. A resistance-capacity model structure is proposed to simulate the display cases. Semi-thermodynamic models are proposed to estimate reciprocating compressors volumetric efficiency and power consumption. Furthermore the zone modeling approach is used for evaporator simulation, and air-cooled gas cooler is modeled with distributed modeling method. The expansion valve simulation is based on orifice flow model. To calibrate these models, both manufacture performance data and experimental data are used. The experiments are conducted with a full instrumental CO 2 supermarket refrigeration system installed in Oak Ridge National Lab Flexible Research Platform (FRP). The simulation model can predict the system performance, including power consumption, cooling capacity, mass flow rate, temperature, and pressure, with high accuracy (within ±4%) compared to experimental data. In addition, this developed model has been used to create the system optimum high side pressure for high side expansion valve control, and to generate wide operating range simulation data for developing the batter-equivalent commercial refrigeration system model which can be used in grid interactive control development.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

High-temperature lean Cu alloys with Cr-to-Nb atomic ratio of 2

Two Cu-Cr-Nb alloys, denoted as alloy 1 (comprising Cu-0.89 at% Cr-0.42 at% Nb) and alloy 2 (comprising Cu-1.84 at% Cr-0.99 at% Nb), were produced through a series of manufacturing processes including vacuum induction melting, melt spinning, consolidation, brazing, and baking, with both alloys aimed at achieving a nominal Cr-to-Nb atomic ratio of 2. Microstructural characterization using transmission electron microscopy and X-ray diffraction identified the cubic C15 Laves-phase Cr 2 Nb as the dominant precipitate in both alloys, cross-validated by thermodynamic calculations and atomistic simulation-based density functional theory (DFT). Besides cubic C15 Cr 2 Nb, hexagonal C14-phase Cr 2 Nb and α-BiF3 cubic structured Cr 3 Nb were also observed in the alloys, including a coherent interface formed between the Cr 3 Nb precipitate and the Cu matrix. The hardness of the alloys increases, and the electrical conductivity decreases with increasing alloying addition content; two practical equations described the trends. Further DFT simulations revealed that the electrical conductivity (conductance) of the Cu/Cr 2 Nb interface is an order of magnitude higher than the intrinsic Cu high-angle grain boundaries.

36 MATERIALS SCIENCE↗

Entropy Pair Functional Theory: Direct Entropy Evaluation Spanning Phase Transitions

We prove that, within the class of pair potential Hamiltonians, the excess entropy is a universal, temperature-independent functional of the density and pair correlation function. This result extends Henderson’s theorem, which states that the free energy is a temperature dependent functional of the density and pair correlation. The stationarity and concavity of the excess entropy functional are discussed and related to the Gibbs–Bugoliubov inequality and to the free energy. We apply the Kirkwood approximation, which is commonly used for fluids, to both fluids and solids. Approximate excess entropy functionals are developed and compared to results from thermodynamic integration. The pair functional approach gives the absolute entropy and free energy based on simulation output at a single temperature without thermodynamic integration. We argue that a functional of the type, which is strictly applicable to pair potentials, is also suitable for first principles calculation of free energies from Born–Oppenheimer molecular dynamics performed at a single temperature. In conclusion, this advancement has the potential to reduce the evaluation the free energy to a simple modification to any procedure that evaluates the energy and the pair correlation function.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Predictive coarse-grained (CG) modeling of morphologies in polymer nanocomposites with specific and directional intermolecular interactions (Final Report)

The overarching goal of the proposed work was to develop predictive models for investigating structure and dynamics in soft materials with chemistries that have specific and directional molecular interactions. The motivation behind studying materials with specific and directional interactions lies in the many desirable features these interactions provide when designing novel soft materials. Soft materials with specific and direction interactions (such as hydrogen bonds or H-bonds) can have a) thermally reversible phase behavior with different functions with varying temperature, b) precisely tuned nanostructure with desirable geometries that afford unique physical properties (e.g., color response, mechanical properties) and c) well-mixed/blended morphologies that are useful for variety of applications in energy field (e.g., materials for batteries require use of blended polymers where one domain gives superior mechanical properties and one domain promotes electrical conduction). Notably, biology makes extensive use of specific and directional interactions, in many cases based on H-bonds, to construct materials with precisely defined architectures and properties. Engineered soft materials with precisely tuned nanostructures and improved processiblity through thermoresponsive phase behavior are useful in numerous applications that are relevant to the Department of Energy (DOE) including high efficiency electronic devices, light-weight high-strength composite materials for batteries and fuel cells, and polymer membranes for separations, etc. While past computational studies have been tremendously useful in understanding molecular phenomena and guiding synthesis of new macromolecular soft materials for a wide variety of applications, the inability to capture small scale specific and directional interactions alongside macromolecular length and time scales represented a key limitation of most studies to date. Our work in this project addressed this grand challenge in computational materials chemistry, i.e., ability to model the anisotropic, directional, and specific interactions that govern the behavior of many macromolecular soft matter systems of interest, thus, has greatly expanded the predictive potential of simulations. Specifically the key outcomes were: successful development of new coarse-grained (CG) polymer models to study generic and specific polymer chemistries in which hydrogen-bonding interactions are dominant. These CG models were then used in molecular simulations to study structure and thermodynamics in polymer nanocomposites and blends; some studies were conducted in collaboration with experimentalists. We also published a perspective and a viewpoint which included some of the work we completed in this DOE project; we believe these perspective and viewpoint articles guide other researchers in the soft materials community on how to extend the computational approaches and models we have developed for the purposes of their studies.

36 MATERIALS SCIENCE↗

Quantification of the Land Surface and Brown Ocean Influence on Tropical Cyclone Intensification Over Land

An investigation of TC Kelvin in February 2018 over northeast Australia was conducted to understand the mechanisms of the Brown Ocean effect (BOE), and to develop a comprehensive analysis framework for land falling tropical cyclones (TCs) in the process. NASA’s Land Information System (LIS) coupled to the NASA Unified WRF (NU-WRF) system was employed as the numerical model framework for 12 land/soil moisture perturbation experiments. Impacts of soil moisture and surface enthalpy flux conditions on TC Kelvin were investigated by closely evaluating simulated track and intensity, midlevel atmospheric thermodynamic properties, vertical wind shear, total precipitable water (TPW), and surface moisture flux. The results suggest that there were recognized differentiations among the sensitivity simulations as a result of land surface (e.g., soil moisture and texture) conditions. However, the intensification of TC Kelvin over land was more strongly related to atmospheric moisture advection and the diurnal cycle of solar radiation (i.e., radiative cooling) than to overall soil moisture conditions or surface fluxes. The analysis framework employed here for TC Kelvin can serve as a foundation to specifically quantify the factors governing the BOE. It also demonstrates that the BOE is not a binary influence (i.e. all or nothing), but instead operates in a continuum from largely to minimally influential such that it could be utilized to help improve prediction of inland effects for all land falling TCs.

Jinwoong Yoo↗