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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
NASA NTRS1999

Driving Stresses in Mars Polar Ice Caps and Conditions for Ice Flow

Measurements of the topography of the North polar ice cap by the Mars Orbiter Laser Altimeter (MOLA) show that the ice cap is 2950 +/- 200 meters thick. The volume of the cap is about 1.2 x 10(exp 6) cu km covering an area of 1.04 x 10(exp 6) sq km, which is about 40 percent of the Greenland ice sheet in volume and 62 percent in area. The composition of the Northern cap was previously concluded to be predominately H2O, rather than CO2 ice, based on thermodynamic considerations of the insustainablity of CO2 during summer. Principal questions about the cap are: does the ice move and at what rate, is the cap currently growing or depleting in volume, and how and when was the cap formed? Recent research on terrestrial ice sheets indicates that rates of ice deformation at the low stress values characteristic of ice sheets are significantly higher than the rates given by the classic viscous-plastic flow laws commonly used.

Zwally, H. Jay↗
DOE OSTI2024

Floquet insulators and lattice fermions beyond naive time discretization

Periodically driven quantum systems known as Floquet insulators can host topologically protected bound states known as “ π modes” that exhibit response at half the frequency of the drive. Such states can also appear in undriven lattice field theories when time is discretized as a result of fermion doubling, raising the question of whether these two phenomena could be connected. Recently we demonstrated such a connection at the level of an explicit mapping between the spectra of a continuous-time Floquet model and a discrete-time undriven lattice fermion model. However, this mapping relied on a symmetry of the single-particle spectrum that is not present for generic drive parameters. Inspired by the example of the temporal Wilson term in lattice field theory, in this paper we extend this mapping to the full drive parameter space by allowing the parameters of the discrete-time model to be frequency-dependent. The spectra of the resulting lattice fermion models exactly match the quasienergy spectrum of the Floquet model in the thermodynamic limit. Our results demonstrate that spectral features characteristic of beyond-equilibrium physics in Floquet systems can be replicated in static systems with appropriate time discretization. Published by the American Physical Society 2024

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Thermodynamics of electrolyte solutions near charged surfaces: Constant surface charge vs constant surface potential

Electric double layers are ubiquitous in science and engineering and are of current interest, owing to their applications in the stabilization of colloidal suspensions and as supercapacitors. While the structure and properties of electric double layers in electrolyte solutions near a charged surface are well characterized, there are subtleties in calculating thermodynamic properties from the free energy of a system with charged surfaces. These subtleties arise from the difference in the free energy between systems with constant surface charge and constant surface potential. In this work, we present a systematic, pedagogical framework to properly account for the different specifications on charged bodies in electrolyte solutions. Our approach is fully variational—that is, all free energies, boundary conditions, relevant electrostatic equations, and thermodynamic quantities are systematically derived using variational principles of thermodynamics. We illustrate our approach by considering a simple electrolyte solution between two charged surfaces using the Poisson–Boltzmann theory. Our results highlight the importance of using the proper thermodynamic potential and provide a general framework for calculating thermodynamic properties of electrolyte solutions near charged surfaces. Specifically, we present the calculation of the pressure and the surface tension between two charged surfaces for different boundary conditions, including mixed boundary conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRSDate not supplied

Uncertainty Quantification of Classical Theories of Dendritic Growth Kinetics Applied to Nickel-Based Alloys

The solidification velocity in a model nickel-alloy single crystal during laser spot melting was recently characterized using synchrotron X-ray imaging. The measured solidification velocity was found to exceed the absolute stability threshold predicted by the Kurz-Giovanola-Trivedi (KGT) model. The discrepancies between the model and experiments motivate the further assessment of accurate material properties. This work quantifies the impact material property uncertainty has on model predictions of the absolute stability threshold velocities. Properties from the literature are reviewed and compared to those calculated using computational thermodynamics to provide uncertainty estimates on input properties to the KGT model. Global sensitivity analysis is used to quantify the influence of each uncertain input property on the predicted threshold velocity. This work supports the understanding of the nickel-alloy solidification during powder bed fusion additive manufacturing and identifies the solidification material properties that are the most important to assess from first-principles computations and experiments.

Computational thermodynamics↗
DOE OSTI2025

Roadmap on thermodynamics and thermal metamaterials

Thermal metamaterials represent a transformative paradigm in modern physics, synergizing thermodynamic principles with metamaterial engineering to master heat flow at will. As next-generation technologies demand multi-scale thermal control, this field urgently requires systematic frameworks to unify its multidisciplinary advances. Curated through a global collaboration involving over 50 specialists across 25 subdisciplines, this review primarily summarizes two decades of advancements, ranging from theoretical breakthroughs to functional implementations. The review reveals groundbreaking innovations in heat manipulation through the exploration of both classical and non-classical transport regimes, topological thermal control mechanisms, and quantum-informed phonon engineering strategies. By bridging physical insights like non-Hermitian thermal dynamics and valleytronic phonon transport with cutting-edge applications, we demonstrate paradigm-shifting capabilities: environment-adaptive thermal cloaks, AI-optimized metamaterials, and nonlinear thermal circuits enabling heat-based computation. Experimental milestones include 3D thermal null media with reconfigurable invisibility and thermal designs breaking classical conductivity limits. Here, this collaborative effort establishes an indispensable roadmap for physicists, highlighting pathways to quantum thermal management, entropy-controlled energy systems, and topological devices. As thermal metamaterials transition from laboratory marvels to technological cornerstones, this work provides the foundational lexicon and design principles for the coming era of intelligent thermal matter.

heat conduction control↗
DOE OSTI2020

Dissipation of electron-beam-driven plasma wakes

Metre-scale plasma wakefield accelerators have imparted energy gain approaching 10 gigaelectronvolts to single nano-Coulomb electron bunches. To reach useful average currents, however, the enormous energy density that the driver deposits into the wake must be removed efficiently between shots. Yet mechanisms by which wakes dissipate their energy into surrounding plasma remain poorly understood. Here, we report picosecond-time-resolved, grazing-angle optical shadowgraphic measurements and large-scale particle-in-cell simulations of ion channels emerging from broken wakes that electron bunches from the SLAC linac generate in tenuous lithium plasma. Measurements show the channel boundary expands radially at 1 million metres-per-second for over a nanosecond. Simulations show that ions and electrons that the original wake propels outward, carrying 90 percent of its energy, drive this expansion by impact-ionizing surrounding neutral lithium. The results provide a basis for understanding global thermodynamics of multi-GeV plasma accelerators, which underlie their viability for applications demanding high average beam current.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

FSISPH: An SPH formulation for impacts between dissimilar materials

Here, we present an SPH formulation with several new features designed to better model the fully-compressible interaction of dissimilar materials. We developed the new method to simulate the atmospheric entry and break-up of small celestial bodies in planetary atmospheres. The formulation uses a unity-based, density-energy discretization of the hydrodynamic conservation laws with linear-corrected kernel gradients. To account for variations in compressibility, we use an HLLC approximate Riemann solver to adjust the velocity gradient at material interfaces. To handle large transverse velocity discontinuities, we introduce a simple slip interface model that limits the artificial viscosity at material interfaces. Diffusion is optionally applied through the velocity gradient and this allows the density and specific thermal energy to evolve in a manner more consistent with the first law of thermodynamics in comparison to other more direct diffusion schemes. We also introduce a material-local second-order artificial conduction scheme used to smooth the specific thermal energy field. Material damage fits neatly under this framework by treating the damage front as a material interface. The method has been implemented as a solver, FSISPH, within the code, Spheral++, and is publicly available on github. We test our new solver on a number of classic shock, mixing, and multi-material problem. The components we outline can significantly improve accuracy of SPH for problems with sharp contact discontinuities.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1992

Calorimetric determination of the thermoneutral potential of Li/BCX and Li/SOCl2 cells

Results are presented of the calorimetric determination of the effective thermoneutral potential, Eetp, of Li/BCS and Li/SOCl2 cells in the temperature range 0-60 C through a continuous recording of the cell voltage, heat flow, and current. The average effective thermoneutral potential at 25 C was 4.0 and 3.84 V for BCX and Li/SOCl2 cells, respectively. Based on the classical approach, the reversible cell potential, Er, and temperature dependence of reversible cell potential, dEr/dT, for BCX cell were 3.74 V and -0.952 mV/K, respectively, and for Li/SOCl2, Er = 3.67 V and dEr/dT = 0.567 mV/K. The thermal polarization (Eetp-E1), where E1 is the load voltage, for both cells, showed that they are the most thermally efficient near 40 C. An overall reaction proposed for the BCX chemistry is supported by the calculated thermodynamic parameters.

Kalu, E. E.↗
DOE OSTI2026

Molecular dynamics simulation of hydrodynamic transport coefficients in plasmas

Molecular dynamics simulations are used to calculate transport coefficients in a two-component plasma interacting through a repulsive Coulomb potential. The thermal conductivity, electrical conductivity, electrothermal coefficient, thermoelectric coefficient, and shear viscosity are computed using the Green–Kubo formalism over a broad range of Coulomb coupling strength, 0.01 ≤ Γ ≤ 140. Emphasis is placed on testing standard results of the Chapman–Enskog solution in the weakly coupled regime (Γ ≪ 1) using these first-principles simulations. As expected, the results show good agreement for Γ ≲ 0.1. However, this agreement is only possible if careful attention is paid to the definitions of linear constitutive relations in each of the theoretical models, a point that is often overlooked. For example, the standard Green–Kubo expression for thermal conductivity is a linear combination of thermal conductivity, electrothermal, and thermoelectric coefficients computed in the Chapman–Enskog formalism. Meaningful results for electrical conductivity are obtained over the full range of coupling strengths explored, but it is shown that potential and virial components of the other transport coefficients diverge in the strongly coupled regime (Γ ≫ 1). In this regime, only the kinetic components of the transport coefficients are meaningful for a classical plasma.

Electrical conductivity↗
DOE OSTI2020

Thermodynamics of rotating quantum matter in the virial expansion

We characterize the high-temperature thermodynamics of rotating bosons and fermions in two-dimensional (2D) and three-dimensional (3D) isotropic harmonic trapping potentials. In this work, we begin by calculating analytically the conventional virial coefficients $b_n$ for all $\textit{n}$ in the noninteracting case, as functions of the trapping and rotational frequencies. We also report on the virial coefficients for the angular momentum and associated moment of inertia. Using the $b_n$ coefficients, we analyze the deconfined limit (in which the angular frequency matches the trapping frequency) and derive explicitly the limiting form of the partition function, showing from the thermodynamic standpoint how both the 2D and 3D cases become effectively homogeneous 2D systems. To tackle the virial coefficients in the presence of weak interactions, we implement a coarse temporal lattice approximation and obtain virial coefficients up to third order.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Scaling functions for eigenstate entanglement crossovers in harmonic lattices

For quantum matter, eigenstate entanglement entropies obey an area law or log-area law at low energies and small subsystem sizes and cross over to volume laws for high energies and large subsystems. This transition is captured by crossover functions, which assume a universal scaling form in quantum critical regimes. We demonstrate this for the harmonic lattice model, which describes quantized lattice vibrations and is a regularization for free scalar field theories, modeling, e.g., spin-0 bosonic particles. In one dimension, the ground-state entanglement obeys a log-area law. For dimensions d ≥ 2, it displays area laws, even at criticality. The distribution of excited-state entanglement entropies is found to be sharply peaked around subsystem entropies of corresponding thermodynamic ensembles in accordance with the eigenstate thermalization hypothesis. Numerically, we determine crossover scaling functions for the quantum critical regime of the model and do a large-deviation analysis. We show how infrared singularities of the system can be handled and how to access the thermodynamic limit using a perturbative trick for the covariance matrix. Eigenstates for quasi-free bosonic systems are not Gaussian. Here, we resolve this problem by considering appropriate squeezed states instead. For these, entanglement entropies can be evaluated efficiently.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1983

Homogeneous crystal nucleation in binary metallic melts

A method for calculating the homogeneous crystal nucleation frequency in binary metallic melts is developed. The free energy of crystallization is derived from regular solution models for the liquid and solid and is used, together with model-based estimates of the interfacial tension, to calculate the nucleation frequency from the classical theory. The method can account for the composition dependence of the maximum undercooling observed in a number of experiments on small droplet dispersions. It can also be used to calculate the driving force for crystal growth and to obtain more precise estimates of the homogeneous crystal nucleation frequency in glass-forming alloys. This method, although approximate, is simple to apply, and requires only knowledge of the phase diagram and a few readily available thermodynamic quantities as input data.

Thompson, C. V.↗
DOE OSTI2023

On the black hole/string transition

We discuss aspects of the possible transition between small black holes and highly excited fundamental strings. We focus on the connection between black holes and the self gravitating string solution of Horowitz and Polchinski. This solution is interesting because it has non-zero entropy at the classical level and it is natural to suspect that it might be continuously connected to the black hole. Surprisingly, we find a different behavior for heterotic and type II cases. For the type II case we find an obstruction to the idea that the two are connected as classical solutions of string theory, while no such obstruction exists for the heterotic case. We further provide a linear sigma model analysis that suggests a continuous connection for the heterotic case. We also describe a solution generating transformation that produces a charged version of the self gravitating string. This provides a fuzzball-like construction of near extremal configurations carrying fundamental string momentum and winding charges. We provide formulas which are exact in α' relating the thermodynamic properties of the charged and the uncharged solutions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

First-principle investigation on catalytic hydrogenation of benzaldehyde over Pt-group metals

We report that understanding the hydrogenation of organic compounds in the aqueous phase has always been fundamentally important for improving carbon neutral pathways to fuels and value-added chemicals. In this study, we investigated both thermodynamic and kinetic profiles of benzaldehyde hydrogenation over the Pd(111) and Pt(111) metal surfaces using density functional theory (DFT) and ab initio molecular dynamic (AIMD) simulations. The adsorption of H 2 shows the mixed preference of H adsorption sites on the Pt(111), while the fcc adsorption site is dominant for H on the Pd(111). When benzaldehyde is added to the systems, we observe a strong reduction of benzaldehyde on charged Pd (111) surface compared with that on neutral surface. In contrast, charged state of the Pt(111) surface does not change their interaction. Subsequent hydrogenation reaction of benzaldehyde over Pd(111), proceeding via Langmuir-Hinshelwood mechanism, is affected by two major factors: the presence of H 2 O solvent and surface charge. The presence of H 2 O solvent greatly reduces the activation energy of C-H and O-H bond formation during the hydrogenation process. Furthermore, the hydrogenation step via C-H bond formation is preferred thermodynamically and kinetically over O-H bond formation during thermocatalytic hydrogenation, while the opposite trend holds true during electrocatalytic hydrogenation.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2022

(U) Marker Lagrangian Simulations using Geometry from CAD Assembly Files [Slides]

The Eulerian code PAGOSA now contains an option for the forward time integration of any user selected material in the Lagrangian frame of reference. This option is accomplished using a Multifield version of Brackbill’s FLIP plus Sulsky’s MPM. FLIP+MPM is a Marker Lagrangian scheme that is not to be confused with well known Grid Lagrangian schemes. This Marker Lagrangian method allows materials to undergo completely arbitrary deformation, up to and including fracture and separation. Using this option makes PAGOSA into a mixed frame simulation tool such that fluids can be integrated in the Eulerian frame and solids integrated in the Lagrangian frame. Boundaries between the fields are said to be ’immersed’ in the computational domain, and their nature is formulated using Multifield Theory. Each Marker in a field of FLIP+MPM material tracks the thermodynamic state of a small (but finite) piece of mass. Initiation of the Markers using the geometry and material information from CAD Assembly files is the subject of the talk. The ease and speed with which the dynamics of the assembly can be studied is demonstrated.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Topology stabilized fluctuations in a magnetic nodal semimetal

Abstract The interplay between magnetism and electronic band topology enriches topological phases and has promising applications. However, the role of topology in magnetic fluctuations has been elusive. Here, we report evidence for topology stabilized magnetism above the magnetic transition temperature in magnetic Weyl semimetal candidate CeAlGe. Electrical transport, thermal transport, resonant elastic X-ray scattering, and dilatometry consistently indicate the presence of locally correlated magnetism within a narrow temperature window well above the thermodynamic magnetic transition temperature. The wavevector of this short-range order is consistent with the nesting condition of topological Weyl nodes, suggesting that it arises from the interaction between magnetic fluctuations and the emergent Weyl fermions. Effective field theory shows that this topology stabilized order is wavevector dependent and can be stabilized when the interband Weyl fermion scattering is dominant. Our work highlights the role of electronic band topology in stabilizing magnetic order even in the classically disordered regime.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Unraveling the solvation phenomena of Ln 3+ cations in room-temperature ionic liquids: A computational study

Rare-earth elements (REEs), classified as critical materials, are difficult to separate due to their similar chemical properties and slight differences in ionic radius. Room-temperature ionic liquids (RTILs) have gained significant attention for REE separation because of their unique physicochemical properties and environmental advantages. In this study, we have investigated the solvation mechanisms of two RE cations (Ln 3+ ), Nd 3+ (light REE) and Yb 3+ (heavy REE), in two 1-butyl-3-methylimidazolium ([BMIM] + )-based RTILs using bis(trifluoromethylsulfonyl)imide ([NTf 2 ] − ) and acetate ([OAc] − ) as the respective anions, employing classical molecular dynamics (MD) simulations and density functional theory (DFT) techniques. From MD simulations, it is evident that Ln 3+ cations are primarily solvated by RTIL anions in the first solvation shell, while [BMIM] + cations form the second solvation shell with their numbers depending on the size and composition of the first solvation shell. The neutralizing [NO 3 − ] counterions are mainly solvated by [BMIM] + cations in their first solvation shell. Relative free energy change of solvation calculations using thermodynamic integration (TI) method indicate that Yb 3+ is more strongly solvated than Nd 3+ in both RTILs, a trend further supported by DFT calculations. While both methods predict consistent qualitative behavior, differences in models and energy evaluations lead to variations in absolute values. Overall, this study provides a comprehensive understanding of the solvation behavior of Ln 3+ cations in RTILs, demonstrating a stronger solvation preference for the heavy Ln 3+ cation (Yb 3+ ). In conclusion, these findings have implications for the design of RTIL-based separation processes for REEs.

Ash, Tamalika [Ames Lab., and Iowa State Univ., Am↗
NASA NTRS1979

Magnetic Cepheids

The role of a magnetic field in the pulsational behavior of classical Cepheids has been studied by computing linearized models of pulsating stellar envelopes pervaded by a well-tangled magnetic field. It is found that the 'pulsational' masses of Cepheids that are implied by the theoretical values of the quantities Q(0), P(1)/P(0), and (possibly) P(2)/P(0) can be brought into line with the large 'evolutionary' masses, if the pressure due to the postulated magnetic field is assumed to be comparable to the thermodynamic pressure everywhere in the pulsating layers. The field strengths that are required - several hundred gauss at the surface and several ten thousand gauss at the base - are quite reasonable from the point of view of the observed magnetic fields in Cepheids. The influence of the magnetic field on the location of the blue edge of the instability strip in the H-R diagram is estimated to be very slight.

Stothers, R.↗