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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
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↗
DOE OSTI2021

Complement to thermodynamics of dyonic Taub-NUT-AdS spacetime

We examine the thermodynamics of Euclidean dyonic Taub-Nut/Bolt-AdS 4 black holes for a variety of horizon geometries to understand how gauge field regularity conditions influence the thermodynamic relations. We find several distinct features that distinguish the NUT-charged case from its dyonic Reissner-Nordstrom counterpart. For the Nut solution, the gauge field vanishes at the horizon and so regularity is ensured. For the Bolt solution we find that the norm of the gauge field is required to vanish at the horizon in order to satisfy both regularity and the first law of thermodynamics. This regularity condition yields a constraint on the electric and magnetic charges and so reduces cohomogeneity of the system; for spherical horizons, the regularity condition removing the Misner string singularity further reduces cohomogeneity. We observe that bolt solutions with increasing electric charge have positive heat capacity, but upon turning on the magnetic charge to make the solution dyonic, we find that the properties of the uncharged one are retained, having both positive and negative heat capacity. We also study the extremal Bolt solution, finding that Misner string disappears at the horizon in the zero temperature limit. We find that the extremal solution has finite-temperature-like behaviour, with the electric potential playing a role similar to temperature.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Optimization techniques in self-similar compressible flow

We investigate the one-dimensional (1D) inviscid compressible flow equations for an ideal gas through the lens of optimization techniques. It is the case that, to our knowledge, optimization analysis applied to the so-called “linear velocity” solutions of the Euler compressible flow equations has not been previously conducted. Through both gradient-based and variational techniques, new variants of well-studied flow scenarios, i.e., self-similar, 1D, linear velocity solution class to idealized inviscid compressible flow equations, are determined, as encoded in both the kinematic and thermodynamic properties of this self-similar solution class. With the kinematics of the said solutions being driven by a self-similar “scale radius” and the thermodynamics being driven separately through the appearance of an arbitrary function, a myriad of new solution classes is possible. Acting as a guide to more realistic physical circumstances as well as discovery, it is the hope that the presented cases serve as the framework for future investigations into the intersection of self-similarity and optimization techniques. Fields of study that may find this work to be of interest include aerodynamic design, flow control, inertial confinement fusion, physics-informed neural networks, and other related areas of interest.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Improved reference system for the corrected rigid spheres equation of state model

The Corrected Rigid Spheres (CRIS) EOS model, developed from fluid perturbation theory using a hard sphere reference system, has been successfully used to calculate the EOS of many materials, including gases and metals. The radial distribution function (RDF) plays a pivotal role in choosing the sphere diameter, through a variational principle, as well as the thermodynamic response. Despite its success, the CRIS model has some shortcomings in that it predicts too large a temperature for liquid-vapor critical points, can break down at large compression, and is computationally expensive. In this paper, we first demonstrate that an improved analytic representation of the hard sphere RDF does not alleviate these issues. Relaxing the strict adherence of the RDF to hard spheres allows an accurate fit to the isotherms and vapor dome of the Lennard-Jones fluid using an arbitrary reference system. Second order correction are eliminated, limiting the breakdown at large compression and significantly reducing the computation cost. The transferrability of the new model to real systems is demonstrated on argon, with an improved vapor dome compared to the original CRIS model.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Self-compensation in heavily Ge doped AlGaN: A comparison to Si doping

Self-compensation in Ge- and Si-doped Al 0.3 Ga 0.7 N has been investigated in terms of the formation of III vacancy and donor-vacancy complexes. Both Ge- and Si-doped AlGaN layers showed a compensation knee behavior with impurity compensation (low doping regime), compensation plateau (medium doping regime), and self-compensation (high doping regime). A maximum free carrier concentration of 4–5 × 10 19 cm –3 was obtained by Ge doping, whereas Si doping resulted in only half of that value, ~2 × 10 19 cm –3 . A DFT calculation with the grand canonical thermodynamics model was developed to support the hypothesis that the difference in self-compensation arises from the difference in the formation energies of the V III -n•donor complexes relative to their onsite configurations. The model suggested that the V III -2•donor and V III -3•donor complexes were responsible for self-compensation for both Ge- and Si-doped AlGaN. However, a lower free carrier concentration in Si-doped samples was due to a high V III -3•Si concentration, resulting from a lower energy of formation of V III -3•Si.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
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 OSTI2022

Modeling atomically mixed graded density impactors

Graded density impactors (GDIs) are multi-material composite impactors used in gas gun experiments to tailor the drive conditions imparted to a sample test material. Previous graded density impactors generally rely on thin, but discrete, layers of different materials. The thinner and the greater number of layers will result in smoother compression. Taken to the limit of very thin layers would be pure material 1 at one surface, such as the front surface of an impactor, smoothly transitioning at the atomic scale to pure material 2 on the back surface. Such an impactor can initially shock, then smoothly compresses a material during a dynamic experiment. This type of experiment can serve to explore a larger region of thermodynamic space than a single or even multi-shock experiments. An overview of how graded density impactors are made is reviewed and sample results are given. Additionally, a strategy for modeling these kinds of impactors is presented. The length scales of constituent mixing are given from the experimental build through electrochemical-deposition. Equation of state models for pure constituents and their subsequent mixtures are presented. It is demonstrated that the time scales for pressure and temperature equilibration, for atomically mixed GDIs, are short enough to be a justifiable closure for the resulting multiphase flow. Furthermore, we present simulation results of dynamic shock followed by a ramp compression, utilizing a silver/gold graded density impactor, onto a tantalum sample.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Rigorous incorporation of pH effects into ab initio electrochemical models

This report summarizes the work carried out with support of the United States Department of Energy under Award DE-SC0023410. The theme of this project was to update electronically grand canonical atomistic simulations to properly incorporate pH effects, on a rigorous thermodynamic basis, and to compare this to experimental findings. This project was primarily conducted by PI Andrew Peterson and his group at Brown University, with an unfunded collaboration with the Quantum Simulations Group at Lawrence Livermore National Laboratory with Dr. Joel Varley, who hosted the primary graduate student (Alexander Lackey) funded by this project and contributed to all major intellectual pursuits, utilizing the expertise and capabilities of the national laboratory to advance the project goals. Postdoctoral associate Dr. Juye Kim was also a major contributor to this project, and implementation of the final section included cooperation with several electronic-structure open-source software developers, including Colin Baker (Brown), Sandeep Sharma (CalTech), Georg Kastlunger (DTU), and Jens Jørgen Mortensen (DTU).

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

Beyond Thermodynamics: Assessing the Dynamical Softness of Hydrated Ions from First Principles

Ion (de)hydration is a key rate-determining step in interfacial processes from corrosion to electrochemical energy storage. However, predicting the kinetics of ion (de)hydration remains challenging, prompting the use of static proxies such as hydration energy and valence. While useful for assessing thermodynamic preferences, such descriptors cannot fully capture the dynamical softness of the hydration shell that dictates kinetics. Accordingly, we use first-principles molecular dynamics to analyze hydration shell softness for a diverse set of metal cations. Additionally, three dynamic metrics are introduced to intuitively describe the bond rigidity, shape deformability, and exchange fluidity of the solvation shell. Together, these metrics capture the relevant physics in the static descriptors, while offering a far more complete and efficient representation for the overall propensity for (de)hydration. Application to the hydrated ion set demonstrates a weak connection between dynamical softness and hydration energy, confirming that dynamical descriptors of hydration are key for correctly describing ion transfer processes.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Equilibrium spin polarization arising from chirality

Chirality-induced spin selectivity (CISS) describes how chiral molecules and materials generate spin polarization even at thermal equilibrium. This observation has challenged established principles of microscopic reversibility and Onsager reciprocity. We resolve this paradox by formulating a pseudo-Hermitian quantum framework that separates thermodynamic equilibrium from time-reversal symmetry. Within this approach, structural chirality and electron correlations, irrespective of their microscopic origin, are sufficient to produce CISS observables. Chirality enters through a non-local metric η that couples spin and spatial motion, leading to real spectra, unitary evolution, and thermodynamic consistency. The framework predicts a chirality-induced spin magnetic ordering characterized by a spin-displacement order, which reconciles equilibrium spin polarization with detailed balance and explains the persistence of CISS in materials composed of light elements. We derive generalized Onsager-Casimir relations that respect the observed CISS symmetry, i.e., parity $(\mathscr{P})$-odd and time-reversal $(\mathscr{T})$-odd, but exhibiting $\mathscr{PT}$-even symmetry. This approach establishes a coherent foundation for equilibrium CISS and provides a route to link chemical chirality with measurable spin-to-charge conversion effects.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Fast estimation of ion-pairing for screening electrolytes: A cluster can approximate a bulk liquid

The propensity for ion-pairing can often dictate the thermodynamic and kinetic properties of electrolyte solutions. Fast and accurate estimates of ion-pairing can thus be extremely valuable for supplementing design and screening efforts for novel electrolytes. We introduce an efficient cluster model to estimate the local ion-pair potential-of-mean-force between ionic solutes in electrolytes. The model incorporates an ion-pair and a few layers of explicit solvent in a gas-phase cluster and leverages an enhanced sampling approach to achieve high efficiency and accuracy. We employ harmonic restraints to prevent solvent escape from the cluster and restrict sampling of large inter-ion distances. We develop a cluster ion-pair sampling tool that implements our cluster model and demonstrate its potential utility for screening simple and poly-electrolyte systems.

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

Sequence-selective dynamic covalent assembly of information-bearing oligomers

Relatively robust dynamic covalent interactions have been employed extensively to mediate molecular self-assembly reactions; however, these assembly processes often do not converge to a thermodynamic equilibrium, instead yielding mixtures of kinetically-trapped species. Here, we report a dynamic covalent self-assembly process that mitigates kinetic trapping such that multiple unique oligomers bearing covalently coreactive pendant groups are able to undergo simultaneous, sequence-selective hybridization with their complementary strands to afford biomimetic, in-registry molecular ladders with covalent rungs. Analogous to the thermal cycling commonly employed for nucleic acid melting and annealing, this is achieved by raising and lowering the concentration of a multi-role reagent to effect quantitative dissociation and subsequently catalyze covalent bond rearrangement, affording selective assembly of the oligomeric sequences. The hybridization specificity afforded by this process further enabled information encoded in oligomers to be retrieved through selective hybridization with complementary, mass-labeled sequences.

36 MATERIALS SCIENCE↗
DOE OSTI2024

A Free-Standing Boron-Doped Diamond Grid Electrode for Fundamental Spectroelectrochemistry

Spectroelectrochemistry (SEC) is a powerful technique that enables a variety of redox properties to be studied, including formal potential (E o ), thermodynamic values (ΔG, ΔH, ΔS), diffusion coefficient (D), electron transfer stoichiometry (n), and others. SEC requires an electrode which light can pass through while maintaining sufficient electrical conductivity. This has been traditionally composed of metal or metal oxide films atop transparent substrates like glass, quartz, or metallic mesh. Robust electrode materials like boron-doped diamond (BDD) could help expand the environments in which SEC can be performed, but most designs are limited to thin films (~100–200 nm) on transparent substrates less resilient than free-standing BDD. Here, this work presents a free-standing BDD grid electrode (G-BDD) for fundamental SEC measurements, using the well-characterized Fe(CN) 6 3–/4– redox couple as proof-of-concept. With a combination of cyclic voltammetry (CV), thin-layer SEC, and chronoabsorptometry, several of the redox properties mentioned above were calculated and compared. For E o' , n, and D, similar results were obtained when comparing the CV [E o' = +0.279 (±0.002) V vs Ag/AgCl; n = 0.97; D = 4.1 × 10 –6 cm 2 ·s –1 ] and SEC [E o' = +0.278 (±0.001) V vs Ag/AgCl; n = 0.91; D = 5.2 × 10 –6 cm 2 ·s –1 ] techniques. Both values align with what has been previously reported. To calculate D from the SEC data, modification of the classical equation used in chronoabsorptometry was required to accommodate the G-BDD electrode geometry. Overall, this work expands on the applicability of SEC techniques and BDD as a versatile electrode material.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗