Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “thermodynamics of computing”

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 1,081 records · Page 60
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Interactive Exploration of High-Dimensional Phase Diagrams

High-dimensional thermodynamic phase stability databases are becoming increasingly common due to the convergence of three recent trends: (i) the widespread interest in so-called “high-entropy” alloys, (ii) the availability of high-throughput computational assessments of phase stability in broad composition spaces and (iii) the ongoing development of ever-increasingly broad, multicomponent, multiphase CALPHAD databases. Although automated computational tools can readily process such high-dimensional data, scientists are often unable to visualize the relevant phase relations, an ability that is crucial to gaining an intuitive understanding of the stability constraints governing materials design. The present work addresses this need by providing algorithms that enable the interactive exploration of phase equilibria in high-dimensional spaces. These algorithms concentrate the complex nonlinear nonsmooth optimization needed into a preprocessing step that generates a large number of high-dimensional yet elementary graphical primitives. Furthermore, these primitives can then be cross-sectioned to yield 3-dimensional views in a computationally efficient manner that enables an interactive exploration of high-dimensional spaces. All of these operations are highly parallelizable, thus facilitating scaling of this method to large data sets.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Recombination of 𝐵 𝑐 mesons in ultrarelativistic heavy-ion collisions

High-energy heavy-ion collisions have been suggested as a favorable environment for the production of 𝐵 𝑐 mesons, due to a much larger abundance of charm and bottom quarks compared to elementary reactions. Motivated by recent CMS data for 𝐵$^{+}_{𝑐}$ production in Pb-Pb(5.02 TeV) collisions at the Large Hadron Collider (LHC), we deploy a previously developed transport approach for charmonia and bottomonia to evaluate the kinetics of 𝐵 𝑐 mesons throughout the fireball formed in these reactions. The main inputs to our approach are two transport parameters: the 𝐵 𝑐 's reaction rate and equilibrium limit. Both quantities are determined by previous calculations via a combination of charm and bottom sectors. In-medium binding energies of 𝐵 𝑐 mesons are calculated from a thermodynamic 𝑇 matrix with a lattice-QCD constrained potential, and figure in their inelastic reaction rates. Temperature-dependent equilibrium limits include charm- and bottom-quark fugacities based on their initial production. We compute the centrality dependence of inclusive 𝐵 𝑐 production and transverse-momentum (𝑝 𝑇 ) spectra using two different recombination models: instantaneous coalescence and resonance recombination. The main uncertainty in the resulting nuclear modification factors, 𝑅 AA , is currently associated with the 𝐵 𝑐 cross section in elementary 𝑝⁢𝑝 collisions, caused by the uncertainty in the branching ratio for the 𝐵$^{−}_{𝑐}$ → 𝐽/𝜓⁢𝜇 − ⁢$\overline{𝜈}$ decay. Our results indicate a large enhancement of the 𝑅 AA at low 𝑝 𝑇 , with significant regeneration contributions up to 𝑝 𝑇 ≃ 20GeV. Comparisons to CMS data are carried out but firm conclusions will require a more accurate value of the branching ratio, or alternative channels to measure the 𝐵 𝑐 production in 𝑝⁢ collisions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1991

Computational modeling of the pressurization process in a NASP vehicle propellant tank experimental simulation

A multidimensional computational model of the pressurization process in a slush hydrogen propellant storage tank was developed and its accuracy evaluated by comparison to experimental data measured for a 5 ft diameter spherical tank. The fluid mechanic, thermodynamic, and heat transfer processes within the ullage are represented by a finite-volume model. The model was shown to be in reasonable agreement with the experiment data. A parameter study was undertaken to examine the dependence of the pressurization process on initial ullage temperature distribution and pressurant mass flow rate. It is shown that for a given heat flux rate at the ullage boundary, the pressurization process is nearly independent of initial temperature distribution. Significant differences were identified between the ullage temperature and velocity fields predicted for pressurization of slush and those predicted for pressurization of liquid hydrogen. A simplified model of the pressurization process was constructed in search of a dimensionless characterization of the pressurization process. It is shown that the relationship derived from this simplified model collapses all of the pressure history data generated during this study into a single curve.

Sasmal, G. P.↗
NASA NTRS1987

The validation of numerical computations

Information is given in viewgraph form on the validation of computational fluid dynamics codes. Topics covered include the types of validation required, aerodynamic heating to a slender 5 degree cone, heat transfer on cones with isentropic compression surfaces, the validation of turbulence data, modeling of thermodynamic and transport properties, real gas code validation, the Flight Dynamics Laboratory validation efforts, and reacting gas experimental data in low density flow.

Buck, Melvin L.↗
DOE OSTI2025

Framework for phase transitions between the Maxwell and Gibbs constructions at finite temperature

The characteristics of the hadron-to-quark first-order phase transition differ depending on whether charge neutrality is locally or globally fulfilled. In 𝛽-equilibrated matter, these two possibilities correspond to the Maxwell and Gibbs constructions. Recently, we presented a new framework in which a continuously varying parameter allows one to describe a first-order phase transition in intermediate scenarios to the two extremes of fully local and fully global charge neutrality. In this work, we extend the previous framework to finite temperatures and out-of-𝛽 equilibrium conditions, making it available for simulations of core-collapse supernovae and binary neutron star mergers. We investigate its impact on key thermodynamic quantities across a range of baryon densities, temperatures, and electron fractions. We find that when matter is not in 𝛽 equilibrium, the pressure in the mixed phase is not constant even for the case of fully local charge neutrality. Moreover, we compute the thermal index using three different approaches, demonstrating that the finite-temperature extension of an equation of state using a constant thermal index can be ill defined when applied to the mixed phase.

QCD phase transitions↗
DOE OSTI2021

Watch and learn—a generalized approach for transferrable learning in deep neural networks via physical principles

Transfer learning refers to the use of knowledge gained while solving a machine learning task and applying it to the solution of a closely related problem. Such an approach has enabled scientific breakthroughs in computer vision and natural language processing where the weights learned in state-of-the-art models can be used to initialize models for other tasks which dramatically improve their performance and save computational time. Here we demonstrate an unsupervised learning approach augmented with basic physical principles that achieves fully transferrable learning for problems in statistical physics across different physical regimes. By coupling a sequence model based on a recurrent neural network to an extensive deep neural network, we are able to learn the equilibrium probability distributions and inter-particle interaction models of classical statistical mechanical systems. Our approach, distribution-consistent learning, DCL, is a general strategy that works for a variety of canonical statistical mechanical models (Ising and Potts) as well as disordered interaction potentials. Using data collected from a single set of observation conditions, DCL successfully extrapolates across all temperatures, thermodynamic phases, and can be applied to different length-scales. This constitutes a fully transferrable physics-based learning in a generalizable approach.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Atomic structure considerations for the low-temperature opacity of Xenon

Here, we have begun a preliminary investigation into the opacity of Xe at low temperatures (<50 eV). The emissivity and opacity of Xe is a crucial factor in determining the utility of Xe in EUV lithography, with numerous industrial applications. To this end, we have been exploring the accuracy of some approximations used in opacity models for the relevant ion stages of Xe (6 times ionized through 20 times ionized). Due to the complexity of Xe atomic structure, one needs to use full configuration-interaction to properly describe the strong mixing in the various n=4 sub-shells that give rise to the Δn=0 and Δn=1 transitions that dominate the opacity spectrum at low temperatures. Since calculations that include full configuration-interaction for large numbers of configurations quickly become computationally prohibitive, we have explored hybrid calculations, in which full configuration-interaction is retained for the most important transitions, while intermediate coupling is employed for all other transitions. After verifying atomic structure properties, local-thermodynamic-equilibrium (LTE) opacities are generated using the ATOMIC code at selected temperatures and densities and compared to various experimental results.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

SAP-X2C: Optimally-Simple Two-Component Relativistic Hamiltonian with Size-Intensive Picture Change

We present a simple relativistic exact 2-component (X2C) Hamiltonian that models two-electron picture-change effects using Lehtola’s superposition of atomic potentials (SAP) [S. Lehtola, J. Chem. Theory Comput. 15, 1593−1604 (2019)]. The SAP-X2C approach retains the low cost and technical simplicity of the popular 1-electron X2C (1eX2C) predecessor but is significantly more accurate and has a well-defined thermodynamic limit, making it applicable to extended systems (such as large molecules and periodic crystals). The assessment of the SAP-X2C-based Hartree−Fock total and spinor energies, spin−orbit splittings, equilibrium bond distances, and harmonic vibrational frequencies suggests that SAP-X2C is similar to the more complex atomic meanfield (AMF) X2C counterparts in its ability to approximate the 4-component Dirac−Hartree−Fock reference.

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

“Magic” Orientation Angles to Suppress Spin-Driven Hall Currents in Anisotropic 2D Materials with an Ideal Skyrmion Gas

The Hall effect depending on conduction electron spin projection becomes very different for anisotropic 2D crystals. In this case the spin-dependent electron current strongly determined by the orientation angle, θ, of the sample with respect to an applied electric eld. The spin-up and -down components of the direct and Hall charge currents oscillate with the angle 2&theta;. The direct and Hall components of the current have the structure where there are the angle-independent part, oscillation amplitude, and phase shift. All three quantitates strongly depend on the electron spin projection, electron mass ratio, and skyrmion size. We nd that there are "magic" orientation angles where the spin-up, spin-down, and total Hall currents vanish. There is a great interest in computations based on 2D materials with skyrmions. Such properties can be useful for computer logic operations based on skyrmions.

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

An Atomistic Study of Reactivity in Solid-State Electrolyte Interphase Formation for Li/Li7P3S11

Lithium metal batteries offer superior volumetric and gravimetric specific capacities compared to those based on traditional graphite anodes. Although advancements in solid-state electrolytes address safety concerns, challenges remain, particularly regarding interphase formation in lithium metal anodes. This work presents a computational framework based on high-throughput first-principles density functional theory and machine-learning interatomic potentials (MLIPs) including automated iterative, active learning to enable robust computational exploration of interphase formation between lithium metal anodes and an inorganic solid-state electrolyte. As a demonstration, we apply the framework to a Li/Li7P3S11 interface and find that it accurately identifies the experimentally observed, thermodynamically stable interphase products as well as their overall spatial arrangement within a heterogeneous, amorphous layered structure, with Li2S domains of nanocrystallinity. Our simulations show two stages, a fast and slow diffusion reaction regime, that corroborate the relative phase formation rate of Li x P, Li2S, and Li3P. Using the Onsager transport theory, we capture time-dependent ionic diffusion within the reacting interface, including cross-correlation effects. We found that cross-correlation effects between Li-P and P-S ionic motion significantly influence P-ion diffusion, making it highly sensitive to the local environment and potentially leading to "kinetic trapping" of Li-P phases. The passivation of the interface is shown as the ionic fluxes all approach zero, effectively halting interphase growth.

Diffusion↗
DOE OSTI2025

Theoretical and experimental quantification of Suzuki segregation enthalpy and strengthening mechanisms in a binary alloy

Solute segregation to planar defects in metallic alloys has been shown to drastically alter mechanical properties. While various works using first-principles and thermodynamic calculations have studied the fundamental driving forces for solute segregation via the Suzuki criterion, planar defect energy, or a comparison of energies of the HCP-like phase and FCC matrix, a quantitative experimental and computational comparison of equilibrium composition and segregation enthalpies has not yet been reported. In this work, we predict the equilibrium composition and segregation enthalpy to intrinsic stacking faults in a Ni-60Co (at.%) alloy and compare the results to two independent experimental methods. We observed that Co segregates to the innermost two planes of the intrinsic stacking fault, and we found that the experimental segregation enrichment, measured from transmission electron microscopy energy dispersive X-ray spectroscopy, of the faults is 6.8 at.% Co, which is 2.2 at.% less than the predicted value at the same temperature. We also find that the segregation enthalpy measured from the composition profile is −21.1 ± 6.4 meV/atom and separately from differential scanning calorimetry segregation enthalpy is −33.2 meV/atom, whereas the predicted enthalpy is −31 ± 1 meV/atom. Based on these results, we determine that segregation occurs very rapidly, within 8 min at temperatures as low as 36% of the homologous solidus temperature. Furthermore, this analysis provides an overview of the possible dislocation mechanisms responsible for strengthening effects due to solute segregation, and concludes that changes in room temperature hardness from local phase transformation is likely tied to post-segregation room temperature equilibrium partial separation distance.

Ab initio calculation↗
DOE OSTI2024

Wave steepening and shock formation in ultracold neutral plasmas

We present observations of wave steepening and signatures of shock formation during expansion of ultracold neutral plasmas formed with an initial density distribution that is centrally peaked and decays exponentially with distance. The plasma acceleration and velocity decrease at large distance from the plasma center, leading to central ions overtaking ions in the outer regions and the development of a steepening front that is narrow compared to the size of the plasma. The density and velocity change dramatically across the front, and significant heating of the ions is observed in the region of steepest gradients. For a reasonable estimate of electron temperature, the relative velocity of ions on either side of the front modestly exceeds the local sound speed (Mach number M≳1). This indicates that by sculpting steep density gradients, it is possible to create the conditions for shock formation, or very close to it, opening a new avenue of research for ultracold neutral plasmas.

Computational fluid dynamics↗
NASA NTRS1973

New thermodynamic functions for the C3 molecule.

When graphite is used as ablation material in heat shields, very often a prominent carbon species vaporized into the stream is C3. The thermodynamic properties of the C3 molecule are, therefore, important in calculating transport phenomena in the ablation flow field. The nature of the C3 thermodynamic functions has been in doubt because of the uncertain contribution of the bending mode vibrations to the total internal energy of the molecule. An approach for overcoming these difficulties is considered. The results of the computations are presented in the form of graphs and approximating functions.

Pearson, W. E.↗
NASA NTRS1985

Laplace-transform technique for deriving thermodynamic equations from the classical microcanonical ensemble

A direct and convenient method is presented for deriving expressions which equate any thermodynamic state function to averages of specific dynamical functions and their fluctuations over the classical microcanonical distribution. Specific expressions are obtained for a variety of thermodynamic quantities. The effect of various entropy definitions on the results are assessed, and the latter are compared to previous work in the literature. The derived formulas are applied to the analysis of molecular-dynamics computer simulations.

Pearson, E. M.↗
NASA NTRS1994

Microgravity sciences application visiting scientist program

Contract NAS8-38785, Microgravity Experimental and Theoretical Research, is a project involving a large number of individual research programs related to: determination of the structure of human serum albumin and other biomedically important proteins; analysis of thermodynamic properties of various proteins and models of protein nucleation; development of experimental techniques for the growth of protein crystals in space; study of the physics of electrical double layers in the mechanics of liquid interfaces; computational analysis of vapor crystal growth processes in microgravity; analysis of the influence of magnetic fields in damping residual flows in directional solidification processes; crystal growth and characterization of II-VI semiconductor alloys; and production of thin films for nonlinear optics. It is not intended that the programs will be necessarily limited to this set at any one time. The visiting scientists accomplishing these programs shall serve on-site at MSFC to take advantage of existing laboratory facilities and the daily opportunities for technical communications with various senior scientists.

Source record↗
NASA NTRS1994

NLTE model stellar atmospheres with line blanketing near the series limits

In this paper we study the influence of line-merging regions at the intermediate long-wavelength side of a continuum threshold on the computed model atmosphere structure and predicted spectrum. In order to model these regions sufficiently accurately, we have developed two concepts. First, we have extended the occupation probability formalism of Hummer and Mihalas to non-local thermodynamic equilibrium (non-LTE) plasmas. Second, in order to treat the very complicated opacity in the line merging region, we have generalized the concept of opacity distribution functions to treat non-LTE situations. All Rydberg states are consistently included within this framework, so that no arbitrary cutoff of high (LTE) levels is made. We have calculated several pure hydrogen models atmospheres for two effective temperatures, T(sub eff) = 20000 and 35000 K, and discussed the differences between models calculated with various treatments of the line merging. In particular, we have shown that the error in the predicted profiles of Balmer lines resulting from the neglect of line merging is typically of the order of 3-4%, while the errors in the far-UV portion of the Balmer continuum reaches 15-35%. The errors generally decrease with increasing effective temperature. At the same time, the internal accuracy of the models is shown to be about or below 0.5% for all predicted spectral features. We conclude that for interpreting current high-accuracy spectrophotometric observations models including the line merging are necessary, and that the formalism developed in this paper is capable of providing a sufficiently accurate and robust modeling technique.

Hubeny, I.↗
DOE OSTI2021

Computational discovery of stable and metastable ternary oxynitrides

Materials design from first principles enables exploration of uncharted chemical spaces. Extensive computational searches have been performed for mixed-cation ternary compounds, but mixed-anion systems are gaining increased interest as well. Central to computational discovery is the crystal structure prediction, where the trade-off between reliance on prototype structures and size limitations of unconstrained sampling has to be navigated. We approach this challenge by letting two complementary structure sampling approaches compete. We use the kinetically limited minimization approach for high-throughput unconstrained crystal structure prediction in smaller cells up to 21 atoms. On the other hand, ternary—and, more generally, multinary—systems often assume structures formed by atomic ordering on a lattice derived from a binary parent structure. Thus, we additionally sample atomic configurations on prototype lattices with cells up to 56 atoms. Using this approach, we searched 65 different charge-balanced oxide–nitride stoichiometries, including six known systems as the control sample. The convex hull analysis is performed both for the thermodynamic limit and for the case of synthesis with activated nitrogen sources. We identified 34 phases that are either on the convex hull or within a viable energy window for potentially metastable phases. We further performed structure sampling for “missing” binary nitrides whose energies are needed for the convex hull analysis. Among these, we discovered metastable Ce3N4 as a nitride analog of the tetravalent cerium oxide, which becomes stable under slightly activated nitrogen condition µN > +0.07 eV. Given the outsize role of CeO2 in research and application, Ce3N4 is a potentially important discovery.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗