Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Classical thermodynamics”

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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

High temperature melting of dense molecular hydrogen from machine-learning interatomic potentials trained on quantum Monte Carlo

We present results and discuss methods for computing the melting temperature of dense molecular hydrogen using a machine learned model trained on quantum Monte Carlo data. In this newly trained model, we emphasize the importance of accurate total energies in the training. We integrate a two phase method for estimating the melting temperature with estimates from the Clausius–Clapeyron relation to provide a more accurate melting curve from the model. We make detailed predictions of the melting temperature, solid and liquid volumes, latent heat, and internal energy from 50 to 180 GPa for both classical hydrogen and quantum hydrogen. At pressures of roughly 173 GPa and 1635 K, we observe molecular dissociation in the liquid phase. Here, we compare with previous simulations and experimental measurements.

08 HYDROGEN↗
DOE OSTI2020

Thermal interfaces in dynamic compression experiments

Temperature is a difficult thermodynamic variable to measure in dynamic compression experiments. Optical pyrometry is a general-purpose technique for measuring temperature from a radiant surface, but that surface is often the interface between distinct materials with temperatures that vary spatially along the loading direction. This leads to a fundamental problem: how does the measured interface temperature relate to this temperature profile along the compression axis? Numerical analysis of loading history and thermal diffusion at these interfaces shows that seemingly subtle changes in experiment geometry can lead to very different temperature profiles. We compare these results to laboratory temperature measurements of shock-compressed tin.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Metadynamics investigation of lanthanide solvation free energy landscapes and insights into separations energetics

Lanthanide ion solvation chemistry in nonaqueous phases is key to understanding and developing effective separation processes for these critical materials. Due to the complexity and inherent disorder of the solution phase, a comprehensive picture of the solvated metal ion is often difficult to generate solely from conventional spectroscopic approaches and electronic structure calculations, particularly in the extractant phase. In this work, we use classical molecular dynamics (MD) simulation with an advanced sampling technique, metadynamics, supplemented by experimental spectroscopy and speciation analysis, to measure lanthanide solvation free energy landscapes. We define coordination-based collective variables to probe the entire range of solvation configurations in the organic phase of lanthanum (La), europium (Eu), and lutetium (Lu) nitrate salts bound with a commonly used extractant, N,N′-dimethyl, N,N′-dioctylhexylethoxymalonamide (DMDOHEMA). The known lanthanide extraction trend of La ≈ Eu > Lu is readily explained by the measured free energy surfaces, which show consistent DMDOHEMA coordination from La to Eu, followed by loss of DMDOHEMA coordination from Eu to Lu. These simulations suggest how ligand crowding at the metal center can control selectivity, in this case resulting in the opposite extraction trend as observed with other conventional extractants, where the enthalpic contribution from increasing lanthanide charge density across the series dominates the extraction energetics. We also find that the presence of inner-sphere water, verified by time-resolved fluorescence, diversifies the accessible solvation structures. As a result, understanding solvation requires consideration of an entire thermodynamic ensemble, rather than the single dominant lowest-energy structure, as is often considered out of necessity in interpretation of spectroscopic data or in electronic structure-based ligand design approaches. In general, we demonstrate how metadynamics uniquely enables investigation of complex, multidimensional solvation energetic landscapes, and how it can explain selectivity trends where extraction is controlled by more complex mechanisms than simple charge density-based selectivity.

Wang, Xiaoyu↗
DOE OSTI2020

Solvent influence on non-adiabatic interfacial electron transfer at conductive oxide electrolyte interfaces

The kinetics for interfacial electron transfer (ET) from a transparent conductive oxide (tin-doped indium oxide, ITO, Sn:In 2 O 3 ) to molecular acceptors 4-[N,N-di(p-tolyl)amino]benzylphosphonic acid, TPA, and [Ru II (bpy) 2 (4,4' -(PO 3 H 2 ) 2 -bpy)] 2+ , RuP, positioned at variable distances within and beyond the electric double layer (EDL), were quantified in benzonitrile and methanol by nanosecond absorption spectroscopy as a function of the thermodynamic driving force, -ΔG ° . Relevant ET parameters such as the rate constant, k et , reorganization energy, λ, and electronic coupling, Hab, were extracted from the kinetic data. Overall, k et increased as the distance between the molecular acceptor and the conductor decreased. For redox active molecules within the Helmholtz planes of the EDL, k et was nearly independent of -ΔG° , consistent with a negligibly small λ value. Rips–Jortner analysis revealed a non-adiabatic electron transfer mechanism consistent with H ab < 1 cm -1 . The data indicate that the barrier for electron transfer is greatly diminished at the conductor–electrolyte interface.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1994

The chromospheric structure of the cool giant star g Herculis

Non-Local Thermodynamic Equilibrium (LTE) calculations of semiempirical chromospheric models are presented for 30 g Her (M6 III). This star is one of the coolest (T(sub eff) = 3250 K) SRb (semiregular) variable stars and has a mass perhaps as great as 4 solar mass. Chromospheric features we have observed in its spectrum include Mg II h and k; C II) UV0.01, which is sensitive to electron density; Mg I lambda 2852; Ca II H, K, and IRT; Ca I lambda 4227 and lambda 6573; Al II) UV 1; and H alpha. We pay special attention to fitting the C II intersystem lines and the Mg II resonance lines but use all the other features as constraints to some extent. The equations of radiative transfer and statistical equilibrium are solved self-consistently for H I, H(-), H2, He I, C I, C II, Na I, Mg I, Mg II, Al I, Al II, Ca I, and Ca II with the equivalent two-level technique. To simplify these calculations, a one-dimensional hydrostatic, plane-parallel atmosphere is assumed. We investigate 10 separate 'classical' chromospheric models, differing most importantly in total mass column density above the temperature minimum. Synthetic spectra from these models fit some but not all of the observations. These comparisons are discussed in detail. However, we find that no single-component classical model in hydrostatic equilibrium is able to reproduce both the Mg II line profiles and the relative strengths of the CII) lines. In all these models, chromospheric emission features are formed relatively close to the star (approximately less than 0.05 R(sub *). The circumstellar environment has a thick, cool component overlying the Mg II emission region, which is relatively static and very turbulent. Finally, we find that thermalization in the Mg II h and k lines in the coolest giant stars is controlled by continuum absorption from Ca I 4p 4p3 P0 bound-free opacity and not collisional de-excitation as is the case for warmer K giants.

Luttermoser, Donald G.↗
DOE OSTI2023

Light-trapping structures fabricated in situ for ultrathin III-V solar cells

Here, we describe a fully in situ method of fabricating light-scattering structures on III-V materials that generates a rough morphology via vapor phase etching and redeposition. Fully in situ methods support higher industrial throughput by utilizing the growth reactor to generate the light-trapping structures after device growth without removal from the reactor. We use HCl and PH 3 to etch and redeposit scattering morphologies on Ga 0.5 In 0.5 P in a dynamic hydride vapor phase epitaxy (D-HVPE) reactor. We show that the addition of PH 3 leads to redeposition during the vapor phase HCl etching of Ga 0.5 In 0.5 P and that HCl flow rate and time exposed to HCl-PH 3 each independently cause a linear increase in the redeposited feature size, indicating that redeposition proceeds by island growth in a III-Cl-limited, hydride-enhanced HVPE regime. Auger electron spectroscopy and scanning transmission electron microscopy with energy dispersive spectroscopy (STEM-EDS) reveal redeposition to be highly Ga-rich GaInP, i.e., Ga(In)P. The Ga-rich nature of the redeposition results from the higher thermodynamic driving force for Ga incorporation than for In during HVPE growth and the difference in the volatility of the III-Cl etch products. The resulting morphologies have high broadband scattering, as determined by normal specular reflectance and integrating sphere measurements, indicating effectiveness as light-scattering structures. In a 270-nm-thick GaAs photovoltaic device with a textured back surface, we achieve a 4.9% increase in short circuit current density (J SC ) without any loss in open-circuit voltage (V OC ) relative to a planar control using only a 60 s in situ texturing treatment.

14 SOLAR ENERGY↗
DOE OSTI2023

Real-time latent heat emission during dynamic-compression freezing of water

Abstract Dynamic compression studies have been used to study the nucleation kinetics of water to ice VII for decades. Diagnostics such as photon Doppler velocimetry, transmission loss, and imaging have been used to measure pressure/density, and phase fraction, while temperature has remained the difficult thermodynamic property to quantify. In this work, we measured pressure/density and implemented a diagnostic to measure the temperature. In doing so the temperature shows quasi-isentropically compressed liquid water forms ice at pressures below the previously defined metastable limit, and the liquid phase is not hypercoooled as previously thought above that limit. Instead, the latent heat raises the temperature to the liquid-ice-VII melt line, where it remains with increasing pressure. We propose a hypothesis to corroborate these results with previous work on dynamic compression freezing. These results provide constraints for nucleation models, and suggest this technique be used to investigate phase transitions in other materials.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Biasing topological charge injection in topological matter

Here we explore the interplay between topologies in the momentum and real spaces to formulate a thermodynamic description of nonequilibrium injection of topological charges under external bias. We show that the edge modes engendered by the momentum-space topology can play a functional role of connecting the external reservoirs to the bulk transport of topological charges in the real space. We illustrate our general results with two examples: the spin-torque injection of skyrmions in an electrically biased integer quantum Hall system, and the vortex injection in a topological $p + ip$ superconductor coupled to heat reservoirs. Based on the universal fractional entropy of the Majorana zero modes bound to the vortices, their controllable injection proposed in this work could provide a route for creating and manipulating Majorana fermions.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Adsorbate Free Energies from DFT-Derived Translational Energy Landscapes

Adsorption free energies are fundamental to surface chemistry and catalysis. Standard models combine some assumed analytical form of the translational potential energy surface, often parameterized against density functional theory (DFT) calculations, with an analytical expression for the resultant translational densities of states (DOS), free energy, and entropy. Here we compare the performance of such models against numerical evaluations of the DOS and thermodynamic functions derived from solutions to the translational Schrodinger equation. We compare results for a translational potential energy surface (PES) derived from nudged eleastic band calculations with those obtained from adsorbate rastering across a series of monatomic (O, S, C, N, H) and polyatomic (NH x ) adsorbates on (100) Pt and Au facets. We find that analytical models as commonly parameterized have mixed performance for describing the translational PES and that the consequences for computed free energies are modest but potentially significant in microkinetic models. Here, numerical solutions are possible for modest to no additional computational cost over analytical models and thus should be considered when reliable free energy estimates are needed or translational potential energy surfaces are available.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Transfer Learning of High-Fidelity Opacity Spectra in Autoencoders and Surrogate Models

Simulations of high energy density physics are expensive, largely in part for the need to produce nonlocal thermodynamic equilibrium opacities. High-fidelity spectra may reveal new physics in the simulations not seen with low-fidelity spectra, but the cost of these simulations also scales with the level of fidelity of the opacities being used. Neural networks are capable of reproducing these spectra, but neural networks need data to train them, which limits the level of fidelity of the training data. Here this article demonstrates that it is possible to reproduce high-fidelity spectra with median errors in the realm of 3%–4% using as few as 50 samples of high-fidelity Krypton data by performing transfer learning on a neural network trained on many times more low-fidelity data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Large Nernst Effect and Thermodynamics Properties in Weyl Antiferromagnet

The noncollinear antiferromagnet Mn3Sn exhibits a large anomalous Hall effect (AHE) driven by a large Berry curvature despite the vanishingly small magnetization. Theoretically, the large Berry curvature in antiferromagnets is predicted to be induced by the time-reversal breaking Weyl state in momentum space. Here, we comprehensively study the thermal conductivity, specific heat, thermoelectric power S(T), and thermomagnetic power S ji (T) in single crystals of Mn 3 Sn. Mn 3 Sn exhibited a large spontaneous anomalous Nernst effect (ANE) at room temperature; in particular, the signal magnitude of Mn 3 Sn exceeded 0.6 µV/K, which is comparable to that of ordinary ferromagnets. The large ANE results from the magnetic Weyl states near the Fermi level E F and likely shares its origin with the AHE. Moreover, the possibility of Nernst-type thermoelectric modules was evaluated by analyzing the basic thermal and electric physical properties.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Polyconvex neural network models of thermoelasticity

Machine-learning function representations such as neural networks have proven to be excellent constructs for constitutive modeling due to their flexibility to represent highly nonlinear data and their ability to incorporate constitutive constraints, which also allows them to generalize well to unseen data. Here, in this work, we extend a polyconvex hyperelastic neural network framework to (isotropic) thermo-hyperelasticity by specifying the thermodynamic and material theoretic requirements for an expansion of the Helmholtz free energy expressed in terms of deformation invariants and temperature. Different formulations which a priori ensure polyconvexity with respect to deformation and concavity with respect to temperature are proposed and discussed. The physics-augmented neural networks are furthermore calibrated with a recently proposed sparsification algorithm that not only aims to fit the training data but also penalizes the number of active parameters, which prevents overfitting in the low data regime and promotes generalization. The performance of the proposed framework is demonstrated on synthetic data, which illustrate the expected thermomechanical phenomena, and existing temperature-dependent uniaxial tension and tension-torsion experimental datasets.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Hydro-Code Implementation and Testing of a Kinetic Phase Transition Framework

In this report we describe the Kinetic Phase Transition (KPT) framework that has been worked out over the last 10 years (from around 2014) and the implementation of it into three different codes, the one-dimensional hydro- LASLO and the three-dimensional magneto-hydro- ALEGRA, Sandia codes, via subroutines in the LAMBDA Equations of State and constitutive models package, and Flag, an arbitrary Lagrangian-Eulerian multiphysics code developed within the Lagrangian Applications project (LAP) at LANL. We discuss the introduction of phase mass (and/or volume) fractions that are needed in a code for it to be ‘phase aware’, that is, not only the thermodynamic state is known in each point but also the mixture of the materials’ phases in that point. Further we point to the need of a full Equations of State for each phase in a material to achieve phase awareness and we review the equilibrium phase model, where a phase mixture is at its lowest Gibbs free energy state, to make this point clear. Contrasting the kinetic phase transition to this equilibrium model seamlessly introduce us to the KPT framework that is subsequently thoroughly discussed. While the determination of the total state and the states and mass fractions of phases in each point is a problem that can borrow many of its numerical details from Eulerian codes and mixture of materials (not phases), the update of mass fractions with time in a KPT framework needs a new set of considerations. General for any update model is that we need to prevent mass fractions from becoming unphysical (negative or their sum to be larger than one). We have solved this problem by implementing a subdivision of the hydro time step that prevents the phase from being fully present to not present at all in one subdivided time step by limiting the size of the subdivided time step. This scheme also corrects numerical problems from abrupt changes in parameter values, the so called Gibbs phenomena, that gives rise to slushing between phases in the KPT framework. Interspersed throughout the report are discussions on different thermodynamics considerations. EOS validity windows, limitations on the EOS phase space, are needed for the KPT framework and are discussed separately and exemplified. The KPT framework described in this report has been verified by code comparison, but validation is still an active area of research. There is room for improvement in the update model, both in the model for determination of rates and in how to prevent the mass fractions from becoming unphysical. In addition, the parameters in the KPT update model and the placement of the phase boundary in the EOS phase space, and interactions with other constitutive models, are closely related and interfering with each other. One possible way forward is to simultaneously develop KPT parameters, EOS, and constitutive models for each material.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Machine learning accelerated turbulence modeling of transient flashing jets

We report modeling the sudden depressurization of superheated liquids through nozzles is a challenge because the pressure drop causes rapid flash boiling of the liquid. The resulting jet usually demonstrates a wide range of structures, including ligaments and droplets, due to both mechanical and thermodynamic effects. As the simulation comprises increasingly numerous phenomena, the computational cost begins to increase. One way to moderate the additional cost is to use machine learning surrogacy for specific elements of the calculation. This study presents a machine learning-assisted computational fluid dynamics approach for simulating the atomization of flashing liquids accounting for distinct stages, from primary atomization to secondary breakup to small droplets using the Σ - Y model coupled with the homogeneous relaxation model. Notably, the models for thermodynamic non-equilibrium (HRM) and Σ - Y are coupled, for the first time, with a deep neural network that simulates the turbulence quantities, which are then used in the prediction of superheated liquid jet atomization. The data-driven component of this method is used for turbulence modeling, avoiding the solution of the two-equation turbulence model typically used for Reynolds-averaged Navier-Stokes simulations for these problems. Both the accuracy and speed of the hybrid approach are evaluated, demonstrating adequate accuracy and at least 25% faster computational fluid dynamics simulations than the traditional approach. This acceleration suggests that perhaps additional components of the calculation could be replaced for even further benefit. Published under an exclusive license by AIP Publishing.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Magnetic dynamics with Weyl fermions

Transport of charge and valley (axial) degrees of freedom coupled to order-parameter dynamics in magnetic Weyl semimetals is studied in the framework of nonequilibrium thermodynamics. In addition to the established valley-related transport anomalies that are rooted in band-structure topology, we construct dissipative couplings between the three dynamic constituents of the system driven out of equilibrium by electromagnetic perturbations. We show how the valley degree of freedom mediates an effective coupling between the charge and magnetic sectors of the system, through a combination of the chiral anomaly, on the electric side, and the Onsager-paired valley torque and pumping, on the magnetic side. This work complements previous studies of magnetic Weyl semimetals by a more systematic analysis of collective dissipation. We discuss several concrete examples of the valley-mediated current-driven magnetic instabilities and charge pumping, and extend the theory to the antiferromagnetic case.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Nonequilibrium dynamics of spontaneous symmetry breaking into a hidden state of charge-density wave

Abstract Nonequilibrium phase transitions play a pivotal role in broad physical contexts, from condensed matter to cosmology. Tracking the formation of nonequilibrium phases in condensed matter requires a resolution of the long-range cooperativity on ultra-short timescales. Here, we study the spontaneous transformation of a charge-density wave in CeTe 3 from a stripe order into a bi-directional state inaccessible thermodynamically but is induced by intense laser pulses. With ≈100 fs resolution coherent electron diffraction, we capture the entire course of this transformation and show self-organization that defines a nonthermal critical point, unveiling the nonequilibrium energy landscape. We discuss the generation of instabilities by a swift interaction quench that changes the system symmetry preference, and the phase ordering dynamics orchestrated over a nonadiabatic timescale to allow new order parameter fluctuations to gain long-range correlations. Remarkably, the subsequent thermalization locks the remnants of the transient order into longer-lived topological defects for more than 2 ns.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Coordination and thermophysical properties of select trivalent lanthanides in LiCl–KCl

The coordination chemistry of various fission and decay products, such as actinides and lanthanides, are crucial to the commercial deployment of molten salt reactors as they can affect the thermophysical properties. In this report we examined the structure, coordination environment, and physical properties such as the density and the vibrational density of states for three lanthanide species, namely Ce, Eu, and Sm in the LiCl–KCl eutectic system using a combination of quantum mechanics simulations and spectroscopic experiments. Quantum mechanics molecular dynamics (QM-MD) modelling was employed to determine the physical properties of each system resulting in accurate local coordination of each species. Then, the vibrational density of states (DOS) was determined using a two-phase thermodynamic modelling which was then compared to the experimentally obtained Raman spectra of the species in molten LiCl–KCl having the eutectic composition. We find that Ce 3+ , Eu 3+ and Sm 3+ all adopt octahedral local coordination environments in the eutectic salt composition in good agreement with experimental results. Ce3+ is found to fluctuate between an octahedral six-coordinated and a seven-coordinated structure due to the increased local proximity of Cl in the eutectic salt, resulting in a lower fluidicity/diffusivity than the other trivalent lanthanides studied. The thermophysical properties of the eutectic composition with trivalent lanthanides were not significantly different from the pure eutectic salt composition, but several changes were noted.

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

A molecular dynamics study of laser-excited gold

The structural evolution of laser-excited systems of gold has previously been measured through ultrafast MeV electron diffraction. However, there has been a long-standing inability of atomistic simulations to provide a consistent picture of the melting process, leading to large discrepancies between the predicted threshold energy density for complete melting, as well as the transition between heterogeneous and homogeneous melting. We make use of two-temperature classical molecular dynamics simulations utilizing three highly successful interatomic potentials and reproduce electron diffraction data presented by Mo et al. [Science 360, 1451–1455 (2018)]. We recreate the experimental electron diffraction data, employing both a constant and temperature-dependent electron–ion equilibration rate. In all cases, we are able to match time-resolved electron diffraction data, and find consistency between atomistic simulations and experiments, only by allowing laser energy to be transported away from the interaction region. This additional energy-loss pathway, which scales strongly with laser fluence, we attribute to hot electrons leaving the target on a timescale commensurate with melting.

36 MATERIALS SCIENCE↗