Search NASA⌕ Search

SEARCH · Search NASA

Results for “Numerical entropy”

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

Feasibility of measuring the speed of sound of the quark-gluon plasma from the multiplicity and mean 𝑝 𝑇 of ultracentral heavy-ion collisions

The mean transverse momentum ⟨𝑝 𝑇 ⟩ of hadrons has been observed experimentally and in numerical simulations to have a power-law dependence on the hadronic multiplicity 𝑁 in ultracentral relativistic heavy-ion collisions: ⟨𝑝 𝑇 ⟩∝𝑁 𝑏 UC . It has been put forward that this exponent 𝑏 UC is the speed of sound of quark-gluon plasma measured at a temperature determined from ⟨𝑝 𝑇 ⟩. We study step by step the connection between (i) the energy and entropy of hydrodynamic simulations and (ii) experimentally measurable observables. We show that an argument based on energy and entropy should yield an exponent equal to the pressure over energy density 𝑃/ɛ, rather than the speed of sound 𝑐$_s^2$; however, we also observe that ⟨𝑝 𝑇 ⟩ and 𝑁 are not sufficiently accurate proxies for the energy and entropy to make this possible in practice. From simulations, we find that the exponent 𝑏 UC is significantly different whether the “effective volume” is strictly constant or not, a condition that cannot be enforced experimentally. Additional tests using a modified equation of state find that the exponent 𝑏 UC exhibits a variable degree of correlations with the speed of sound and with 𝑃/ɛ, but is not an accurate measurement of either quantity in general.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1994

The r-process and neutrino-heated supernova ejecta

As a neutron star is formed by the collapse of the iron core of a massive star, its Kelvin-Helmholtz evolution is characterized by the release of gravitational binding energy as neutrinos. The interaction of these neutrinos with heated material above the neutron star generates a hot bubble in an atmosphere that is nearly in hydrostatic equilibrium and heated, after approximately 10 s, to an entropy of S/N(sub AS)k greater than or approximately = 400. The neutron-to-proton ratio for material moving outward through this bubble is set by the balance between neutrino and antineutrino capture on nucleons. Because the electron antineutrino spectrum at this time is hotter than the electron neutrino spectrum, the bubble is neutron-rich (0.38 less than or approximately = Y(sub e) less than or approximately = 0.47). Previous work using a schematic model has shown that these conditions are well suited to the production of heavy elements by the r-process. In this paper we have advanced the numerical modeling of a 20 solar mass 'delayed' supernova explosion to the point that we can follow the detailed evolution of material moving through the bubble at the late times appropiate to r-process nucleosynthesis. The supernova model predicts a final kinetic energy for the ejecta of 1.5 x 10(exp 51) ergs and leaves behind a remnant with a baryon mass of 1.50 solar mass (and a gravitational mass of 1.445 solar mass). We follow the thermodynamic and compositional evolution of 40 trajectories in rho(t), T(t), Y(sub e)(t) for a logarithmic grid of mass elements for the last approximately = 0.03 solar mass to be ejected by the proto-neutron star down to the last less than 10(exp -6) solar mass of material expelled at up to approximately = 18 s after core collapse. We find that an excellent fit to the solar r-process abundance distribution is obtained with no adjustable parameters in the nucleosynthesis calculations. Moreover, the abundances are produced in the quantities required to account for the present Galactic abundances. However, at earlier times, this one-dimensional model ejects too much material with entropies S/N(sub A)k approximately 50 and Y(sub e) approximately 0.46. This leads to an acceptable over production of N = 50 nuclei, particularly Sr-88, Y-89, and Zr-90, relative to their solar abundances. We speculate on various means to avoid the early overproduction and/or ejection of N = 50 isotonic nuclei while still producing and ejecting the correct amount of r-process material.

Woosley, S. E.↗
DOE OSTI2025

Performance of wave function and Green's function methods for non-equilibrium many-body dynamics

Theoretical descriptions of the non-equilibrium dynamics of quantum many-body systems essentially employ either (i) explicit treatments, relying on the truncation of the expansion of the many-body wave function, (ii) compressed representations of the many-body wave function, or (iii) evolution of an effective (downfolded) representation through Green's functions. In this work, we select representative cases of each of the methods and address how these complementary approaches capture the dynamics driven by intense field perturbations to non-equilibrium states. Under strong driving, the systems are characterized by strong entanglement of the single-particle density matrix and natural populations approaching those of a strongly interacting equilibrium system. We generate a representative set of results that are numerically exact and form a basis for a critical comparison of the distinct families of methods. We demonstrate that the compressed formulation based on similarity-transformed Hamiltonians (coupled-cluster approach) is practically exact in weak fields and, hence, weakly or moderately correlated systems. Coupled cluster, however, struggles for strong driving fields, under which the system exhibits strongly correlated behavior, as measured by the von Neumann entropy of the single-particle density matrix. The dynamics predicted by Green's functions in the (widely popular) G W approximation are less accurate, but improve significantly upon the mean-field results in the strongly driven regime. Published by the American Physical Society 2025

Reeves, Cian C. (ORCID:0009000642581845)↗
NASA NTRS1989

Advanced rotor computations with a corrected potential method

An unsteady Full-Potential Rotor code (FPR) has been enhanced with modifications directed at improving its drag prediction capability. The potential code has been rewritten with modifications to increase the code accuracy. Also, the shock generated entropy has been included to provide solutions comparable to the Euler equations. Two different weakly interacted boundary layer models have also been coupled to FPR in order to estimate skin-friction drag. One is a two-dimensional integral method and the other is a three-dimensional finite-difference scheme. The new flow solver is able to find accurate inviscid drags without recourse to numerical error tares. This permits the resolution of drag distributions resulting from rotor geometric variations. Good comparisons have been obtained between computed and measured torque for a rectangular and a highly swept model rotor.

Bridgeman, John O.↗
NASA NTRS1976

Development of a computer code for calculating the steady super/hypersonic inviscid flow around real configurations. Volume 1: Computational technique

A numerical procedure has been developed to compute the inviscid super/hypersonic flow field about complex vehicle geometries accurately and efficiently. A second order accurate finite difference scheme is used to integrate the three dimensional Euler equations in regions of continuous flow, while all shock waves are computed as discontinuities via the Rankine Hugoniot jump conditions. Conformal mappings are used to develop a computational grid. The effects of blunt nose entropy layers are computed in detail. Real gas effects for equilibrium air are included using curve fits of Mollier charts. Typical calculated results for shuttle orbiter, hypersonic transport, and supersonic aircraft configurations are included to demonstrate the usefulness of this tool.

Marconi, F.↗
DOE OSTI2023

Deformation and failure of the CrCoNi medium-entropy alloy subjected to extreme shock loading

The extraordinary work hardening ability and fracture toughness of the face-centered cubic (fcc) high-entropy alloys render them ideal candidates for many structural applications. Here, the deformation and failure mechanisms of an equiatomic CrCoNi medium-entropyalloy (MEA) were investigated by powerful laser-driven shock experiments. Multiscale characterization demonstrates that profuse planar defects including stacking faults, nanotwins, and hexagonal nanolamella were generated during shock compression, forming a three-dimensional network. During shock release, the MEA fractured by strong tensile deformation and numerous voids was observed in the vicinity of the fracture plane. High defect populations, nanorecrystallization, and amorphization were found adjacent to these areas of localized deformation. Molecular dynamics simulations corroborate the experimental results and suggest that deformation-induced defects formed before void nucleation govern the geometry of void growth and delay their coalescence. Our results indicate that the CrCoNi-based alloys are impact resistant, damage tolerant, and potentially suitable in applications under extreme conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Towards absolutely stable ergodicity breaking in two and three dimensions

We propose physically reasonable systems capable of avoiding ergodicity at infinite time in the thermodynamic limit, even with generic perturbations and when coupled to a heat bath. In two dimensions, the rainbow loop soup has (stretched) exponentially numerous absolutely stable nonergodic states with diverging energy but vanishing energy density. In three dimensions the rainbow membrane soup has (stretched) exponentially numerous nonergodic states with diverging energy barriers, leading to infinite-time robust ergodicity breaking that even survives coupling to a nonzero temperature heat bath. We describe our results in the language of exact emergent symmetries and demonstrate how the systems avoid common instabilities. Furthermore, our construction naturally connects to quantum dimer models, topologically ordered systems, the group word construction, and Hamiltonians whose low-energy eigenstates exhibit anomalous entanglement entropy.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1976

Development of a computer code for calculating the steady super/hypersonic inviscid flow around real configurations. Volume 2: Code description

A numerical procedure was developed to compute the inviscid super/hypersonic flow field about complex vehicle geometries accurately and efficiently. A second-order accurate finite difference scheme is used to integrate the three-dimensional Euler equations in regions of continuous flow, while all shock waves are computed as discontinuities via the Rankine-Hugoniot jump conditions. Conformal mappings are used to develop a computational grid. The effects of blunt nose entropy layers are computed in detail. Real gas effects for equilibrium air are included using curve fits of Mollier charts. Typical calculated results for shuttle orbiter, hypersonic transport, and supersonic aircraft configurations are included to demonstrate the usefulness of this tool.

Marconi, F.↗
NASA NTRS1975

Computational procedure for three-dimensional boundary layers on aircraft and aerospace vehicles

An implicit numerical procedure is developed to solve the full three-dimensional laminar, boundary layer equations over realistic aircraft and spacecraft configurations. This procedure utilizes splined functions in the direction normal to the surface and finite differences in the cross flow direction. Calculation of off-pitch plane attachment and detachment lines can be calculated. The solution procedure permits the incorporation of a wide variety of turbulent models and boundary conditions including; entropy layer procedures, specified surface temperatures or heat flux distributions, surface catalysis, mass addition, and the coupling to both three-dimensional inviscid codes and surface ablation codes. In addition, reasonable computing times are obtained.

Kendall, R. M.↗
NASA NTRS1989

MHD intermediate shock discontinuities. I - Rankine-Hugoniot conditions

Recent numerical investigations have focused attention once more on the role of intermediate shocks in MHD. Four types of intermediate shock are identified using a graphical representation of the MHD Rankine-Hugoniot conditions. This same representation can be used to exhibit the close relationship of intermediate shocks to switch-on shocks and rotational discontinuities. The conditions under which intermediate discontinuities can be found are elucidated. The variations in velocity, pressure, entropy and magnetic-field jumps with upstream parameters in intermediate shocks are exhibited graphically. The evolutionary arguments traditionally advanced against intermediate shocks may fail because the equations of classical MHD are not strictly hyperbolic.

Kennel, C. F.↗
DOE OSTI2021

Elaboration of Metallic Materials by SPS: Processing, Microstructures, Properties, and Shaping

After a few decades of increasing interest, spark plasma sintering (SPS) has now become a mature powder metallurgy technique, which allows assessing its performances toward fabricating enhanced materials. Here, the case of metals and alloys will be presented. The main advantage of SPS lies in its rapid heating capability enabled by the application of high intensity electric currents to a metallic powder. This presents numerous advantages balanced by some limitations that will be addressed in this review. The first section will be devoted to sintering issues, with an emphasis on the effect of the electric current on the densification mechanisms. Then, typical as-SPS microstructures and properties will be presented. In some cases, they will be compared with that of materials processed by conventional techniques. As such, examples of nanostructured materials, intermetallics, metallic glasses, and high entropy alloys, will be presented. Finally, the implementation of SPS as a technique to manufacture complex, near-net shape industrial parts will be discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Accelerating the Discovery of New, Single Phase High Entropy Ceramics via Active Learning

High-entropy ceramics have garnered interest due to their remarkable hardness, compressive strength, thermal stability, and fracture toughness; yet the discovery of new high-entropy ceramics (out of a tremendous number of possible elemental permutations) still largely requires costly, inefficient, trial-and-error experimental and computational approaches. The entropy forming ability (EFA) factor was recently proposed as a computational descriptor that positively correlates with the likelihood that a 5-metal high-entropy carbide (HECs) will form the desired single phase, homogeneous solid solution; however, discovery of new compositions is computationally expensive. If you consider 8 candidate metals, the HEC EFA approach uses 49 optimizations for each of the 56 unique 5-metal carbides, requiring a total of 2744 costly density functional theory calculations. Here, we describe an orders-of-magnitude more efficient active learning (AL) approach for identifying novel HECs. To begin, we compared numerous methods for generating composition-based feature vectors (e.g., magpie and mat2vec), deployed an ensemble of machine learning (ML) models to generate an average and distribution of predictions, and then utilized the distribution as an uncertainty. Here we then deployed an AL approach to extract new training data points where the ensemble of ML models predicted a high EFA value or was uncertain of the prediction. Our approach has the combined benefit of decreasing the amount of training data required to reach acceptable prediction qualities and biases the predictions toward identifying HECs with the desired high EFA values, which are tentatively correlated with the formation of single phase HECs. Using this approach, we increased the number of 5-metal carbides screened from 56 to 15,504, revealing 4 compositions with record-high EFA values that were previously unreported in the literature. Our AL framework is also generalizable and could be modified to rationally predict optimized candidate materials/combinations with a wide range of desired properties (e.g., mechanical stability, thermal conductivity).

36 MATERIALS SCIENCE↗
NASA NTRS1993

Navier-Stokes relaxation to sinh-Poisson states at finite Reynolds numbers

A mathematical framework is proposed in which it seems possible to justify the computationally-observed relaxation of a two-dimensional Navier-Stokes fluid to a 'most probable', or maximum entropy, state. The relaxation occurs at large but finite Reynolds numbers, and involves substantial decay of higher-order ideal invariants such as enstrophy. A two-fluid formulation, involving interpenetrating positive and negative vorticity fluxes (continuous and square integrable) is developed, and is shown to be intimately related to the passive scalar decay problem. Increasing interpenetration of the two fluids corresponds to the decay of vorticity flux due to viscosity. It is demonstrated numerically that, in two dimensions, passive scalars decay rapidly, relative to mean-square vorticity (enstrophy). This observation provides a basis for assigning initial data to the two-fluid field variables.

Montgomery, David↗
NASA NTRS2018

Nosetip Bluntness Effects on Transition at Hypersonic Speeds: Experimental and Numerical Analysis Under NATO STO AVT-240

The existing database of transition measurements in hypersonic ground facilities has established that the onset of boundary layer transition over a circular cone at zero angle of attack shifts downstream as the nosetip bluntness is increased with respect to a sharp cone. However, this trend is reversed at sufficiently large values of the nosetip Reynolds number, so that the transition onset location eventually moves upstream with a further increase in nosetip bluntness. This transition reversal phenomenon, which cannot be explained on the basis of linear stability theory, was the focus of a collaborative investigation under the NATO STO group AVT-240 on Hypersonic Boundary-Layer Transition Prediction. The current paper provides an overview of that effort, which included wind tunnel measurements in three different facilities and theoretical analysis related to modal and nonmodal amplification of boundary layer disturbances. Because neither first and second-mode waves nor entropy-layer instabilities are found to be substantially amplified to initiate transition at large bluntness values, transient (i.e., nonmodal) disturbance growth has been investigated as the potential basis for a physics-based model for the transition reversal phenomenon. Results of the transient growth analysis indicate that disturbances that are initiated within the nosetip or in the vicinity of the juncture between the nosetip and the frustum can undergo relatively significant nonmodal amplification and that the maximum energy gain increases nonlinearly with the nose radius of the cone. This finding does not provide a definitive link between transient growth and the onset of transition, but it is qualitatively consistent with the experimental observations that frustum transition during the reversal regime was highly sensitive to wall roughness, and furthermore, was dominated by disturbances that originated near the nosetip.

Numerical analysis↗
NASA NTRS1990

Development of a linearized unsteady aerodynamic analysis for cascade gust response predictions

A method for predicting the unsteady aerodynamic response of a cascade of airfoils to entropic, vortical, and acoustic gust excitations is being developed. Here, the unsteady flow is regarded as a small perturbation of a nonuniform isentropic and irrotational steady background flow. A splitting technique is used to decompose the linearized unsteady velocity into rotational and irrotational parts leading to equations for the complex amplitudes of the linearized unsteady entropy, rotational velocity, and velocity potential that are coupled only sequentially. The entropic and rotational velocity fluctuations are described by transport equations for which closed-form solutions in terms of the mean-flow drift and stream functions can be determined. The potential fluctuation is described by an inhomogeneous convected wave equation in which the source term depends on the rotational velocity field, and is determined using finite-difference procedures. The analytical and numerical techniques used to determine the linearized unsteady flow are outlined. Results are presented to indicate the status of the solution procedure and to demonstrate the impact of blade geometry and mean blade loading on the aerodynamic response of cascades to vortical gust excitations. The analysis described herein leads to very efficient predictions of cascade unsteady aerodynamic response phenomena making it useful for turbomachinery aeroelastic and aeroacoustic design applications.

Verdon, Joseph M.↗
DOE OSTI2025

A simple introduction to the SiMPL method for density-based topology optimization

We introduce a novel method for solving density-based topology optimization problems: Sigmoidal Mirror descent with a Projected Latent variable (SiMPL). The SiMPL method (pronounced as “the simple method”) optimizes a design using only first-order derivative information of the objective function. The bound constraints on the density field are enforced with the help of the (negative) Fermi–Dirac entropy, which is also used to define a non-symmetric distance function called a Bregman divergence on the set of admissible designs. This Bregman divergence leads to a simple update rule that is further simplified with the help of a so-called latent variable. Because the SiMPL method involves discretizing the latent variable, it produces a sequence of pointwise-feasible iterates, even when high-order finite elements are used in the discretization. Numerical experiments demonstrate that the method outperforms other popular first-order optimization algorithms. In conclusion, to outline the general applicability of the technique, we include examples with (self-load) compliance minimization and compliant mechanism optimization problems.

Calculus of Variations and Optimization↗
DOE OSTI2022

Deep-learning-based workflow for boundary and small target segmentation in digital rock images using UNet++ and IK-EBM

We report three-dimensional (3D) X-ray micro-computed tomography (μCT) has been widely used in petroleum engineering because it can provide detailed pore structural information for a reservoir rock, which can be imported into a pore-scale numerical model to simulate the transport and distribution of multiple fluids in the pore space. The partial volume blurring (PVB) problem is a major challenge in segmenting raw μCT images of rock samples, which impacts boundaries and small targets near the resolution limit. We developed a deep-learning (DL)-based workflow for accurate and fast partial volume segmentation. The DL model's performance depends primarily on the training data quality and model architecture. This study employed the entropy-based-masking indicator kriging (IK-EBM) to segment 3D Berea sandstone images as training datasets. The comparison between IK-EBM and manual segmentation using a 3D synthetic sphere pack, which had a known ground truth, showed that IK-EBM had higher accuracy on partial volume segmentation. We then trained and tested the UNet++ model, a state-of-the-art supervised encoder-decoder model, for binary (i.e., void and solid) and four-class segmentation. We compared the UNet++ with the commonly used U-Net and wide U-Net models and showed that the UNet++ had the best performance in terms of pixel-wise and physics-based evaluation metrics. Specifically, boundary-scaled accuracy demonstrated that the UNet++ architecture outperformed the regular U-Net architecture in the segmentation of pixels near boundaries and small targets, which were subjected to the PVB effect. Feature map visualization illustrated that the UNet++ bridged the semantic gaps between the feature maps extracted at different depths of the network, thereby enabling faster convergence and more accurate extraction of fine-scale features. The developed workflow significantly enhances the performance of supervised encoder-decoder models in partial volume segmentation, which has extensive applications in fundamental studies of subsurface energy, water, and environmental systems.

02 PETROLEUM↗
DOE OSTI2020

Finite-temperature many-body perturbation theory in the grand canonical ensemble

A finite-temperature many-body perturbation theory is presented, which expands in power series the electronic grand potential, chemical potential, internal energy, and entropy on an equal footing. Sum-over-states and sum-over-orbitals analytical formulas for the second-order perturbation corrections to these thermodynamic properties are obtained in a time-independent, nondiagrammatic, algebraic derivation, relying on the sum rules of the Hirschfelder–Certain degenerate perturbation energies in a degenerate subspace as well as nine algebraic identities for the zeroth-order thermal averages of one- through four-indexed quantities and products thereof. Lastly, they reproduce numerically exactly the benchmark data obtained as the numerical derivatives of the thermal-full-configuration-interaction results for a wide range of temperatures.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗