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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
NASA NTRS2003

Update on Area Production in Mixing of Supercritical Fluids

The focus of this research is on supercritical C7H16/N2 and O2/H2 mixing layers undergoing transitions to turbulence. The C7H16/N2 system serves as a simplified model of hydrocarbon/air systems in gas-turbine and diesel engines; the O2/H2 system is representative of liquid rocket engines. One goal of this research is to identify ways of controlling area production to increase disintegration of fluids and enhance combustion in such engines. As used in this research, "area production" signifies the fractional rate of change of surface area oriented perpendicular to the mass-fraction gradient of a mixing layer. In the study, a database of transitional states obtained from direct numerical simulations of the aforementioned mixing layers was analyzed to investigate global layer characteristics, phenomena in regions of high density-gradient magnitude (HDGM), irreversible entropy production and its relationship to the HDGM regions, and mechanisms leading to area production.

Okongo, Nora↗
NASA NTRS1984

Second- and third-order upwind difference schemes for hyperbolic conservation laws

Second- and third-order two time-level five-point explicit upwind-difference schemes are described for the numerical solution of hyperbolic systems of conservation laws and applied to the Euler equations of inviscid gas dynamics. Nonliner smoothing techniques are used to make the schemes total variation diminishing. In the method both hyperbolicity and conservation properties of the hyperbolic conservation laws are combined in a very natural way by introducing a normalized Jacobian matrix of the hyperbolic system. Entropy satisfying shock transition operators which are consistent with the upwind differencing are locally introduced when transonic shock transition is detected. Schemes thus constructed are suitable for shockcapturing calculations. The stability and the global order of accuracy of the proposed schemes are examined. Numerical experiments for the inviscid Burgers equation and the compressible Euler equations in one and two space dimensions involving various situations of aerodynamic interest are included and compared.

Yang, J. Y.↗
DOE OSTI2025

Higher Hall conductivity from a single wave function: Obstructions to symmetry-preserving gapped edge of (2+1)-dimensional topological order

A (2+1)D topologically ordered phase with U(1) symmetry may or may not have a symmetric gapped edge state, even if both thermal and electric Hall conductivity are vanishing. It has recently been discovered that there are “higher” versions of Hall conductivity valid for fermionic fractional quantum Hall (FQH) states that obstruct symmetry-preserving gapped edge states beyond thermal and electric Hall conductivity. In this paper, we show that one can extract higher Hall conductivity from a single wave function of an FQH state, by evaluating the expectation value of the “partial rotation” unitary, which is a combination of partial spatial rotation and a U(1) phase rotation. This result is verified numerically with the fermionic Laughlin state with 𝜈=1/3 and 1/5, as well as the non-Abelian Moore-Read state. Together with topological entanglement entropy, we prove that the expectation values of the partial rotation completely determine if a bosonic/fermionic Abelian topological order with U(1) symmetry has a symmetry-preserving gappable edge state or not. We also show that thermal and electric Hall conductivity of Abelian topological order can be extracted by partial rotations. Even in non-Abelian FQH states, partial rotation provides the Lieb-Schultz-Mattis type theorem constraining the low-energy spectrum of the bulk-boundary system. The generalization of higher Hall conductivity to the case with Lie group symmetry is also presented.

2-dimensional systems↗
NASA NTRS1989

Euler solvers for transonic applications

The 1980s may well be called the Euler era of applied aerodynamics. Computer codes based on discrete approximations of the Euler equations are now routinely used to obtain solutions of transonic flow problems in which the effects of entropy and vorticity production are significant. Such codes can even predict separation from a sharp edge, owing to the inclusion of artificial dissipation, intended to lend numerical stability to the calculation but at the same time enforcing the Kutta condition. One effect not correctly predictable by Euler codes is the separation from a smooth surface, and neither is viscous drag; for these some form of the Navier-Stokes equation is needed. It, therefore, comes as no surprise to observe that the Navier-Stokes has already begun before Euler solutions were fully exploited. Moreover, most numerical developments for the Euler equations are now constrained by the requirement that the techniques introduced, notably artificial dissipation, must not interfere with the new physics added when going from an Euler to a full Navier-Stokes approximation. In order to appreciate the contributions of Euler solvers to the understanding of transonic aerodynamics, it is useful to review the components of these computational tools. Space discretization, time- or pseudo-time marching and boundary procedures, the essential constituents are discussed. The subject of grid generation and grid adaptation to the solution are touched upon only where relevant. A list of unanswered questions and an outlook for the future are covered.

Vanleer, Bram↗
NASA NTRS1993

An Eulerian/Lagrangian coupling procedure for three-dimensional vortical flows

A coupled Eulerian/Lagrangian method is presented for the reduction of numerical diffusion observed in solutions of 3D vortical flows using standard Eulerian finite-volume time-marching procedures. A Lagrangian particle tracking method, added to the Eulerian time-marching procedure, provides a correction of the Eulerian solution. In turn, the Eulerian solution is used to integrate the Lagrangian state-vector along the particles trajectories. While the Eulerian solution ensures the conservation of mass and sets the pressure field, the particle markers describe accurately the convection properties and enhance the vorticity and entropy capturing capabilities of the Eulerian solver. The Eulerian/Lagrangian coupling strategies are discussed and the combined scheme is tested on a constant stagnation pressure flow in a 90 deg bend and on a swirling pipe flow. As the numerical diffusion is reduced when using the Lagrangian correction, a vorticity gradient augmentation is identified as a basic problem of this inviscid calculation.

Felici, Helene M.↗
DOE OSTI2021

Differentiating Hund from Mott physics in a three-band Hubbard-Hund model: Temperature dependence of spectral, transport, and thermodynamic properties

We study the interplay between Mott physics, driven by Coulomb repulsion U, and Hund physics, driven by Hund's coupling J, for a minimal model for Hund metals, the orbital-symmetric three-band Hubbard-Hund model (3HHM) for a lattice filling of 1/3. Hund-correlated metals are characterized by spin-orbital separation (SOS), a Hund's-rule-induced two-stage Kondo-type screening process, in which spin screening occurs at much lower energy scales than orbital screening. By contrast, in Mott-correlated metals, lying close to the phase boundary of a metal-insulator transition, the SOS window becomes negligibly small and the Hubbard bands are well separated. Using dynamical mean-field theory and the numerical renormalization group as real-frequency impurity solver, we identify numerous fingerprints distinguishing Hundness from Mottness in the temperature dependence of various physical quantities. Furthermore, These include ARPES-type spectra, the local self-energy, static local orbital and spin susceptibilities, resistivity, thermopower, and lattice and impurity entropies. Our detailed description of the behavior of these quantities within the context of a simple model Hamiltonian will be helpful for distinguishing Hundness from Mottness in experimental and theoretical studies of real materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Improved Models of R Coronae Borealis Stars

We present an improved numerical method to model subsolar He+CO-WD merger progenitors of R Corona Borealis stars that builds on our previous work. These improvements include a smooth entropy transition from the core to the envelope of the post-merger, inclusion of single-zone nucleosynthesis to mimic the effects of burning during the merger event, and post-processing the models with a larger nuclear network for analysis of s-process nucleosynthesis. We perform a parameter study to understand the effects of the entropy transition, peak temperature, and overshooting on our models. The models that best agree with observations of R Corona Borealis stars are processed with a much larger nuclear network to investigate s-process nucleosynthesis and the dredge-up of s-process products into the outer envelope in detail. We present a model with a significant enhancement in s-process elements, which also agrees with observed surface abundances and isotopic ratios of 16 O/ 16 O and C/O between 1 and 10. Finally, we find that the neutron exposure and initial neutron densities this model requires to obtain such an enhancement are much more consistent with i-process nucleosynthesis.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1951

The method of characteristics for the determination of supersonic flow over bodies of revolution at small angles of attack

The method of characteristics has been applied for the determination of the supersonic-flow properties around bodies of revolution at a small angle of attack. The system developed considers the effect of the variation of entropy due to the curved shock and determines a flow that exactly satisfies the boundary conditions in the limits of the simplifications assumed. Two practical methods for numerical calculations are given. (author)

Ferri, Antonio↗
NASA NTRS2001

Energy Stable Flux Formulas For The Discontinuous Galerkin Discretization Of First Order Nonlinear Conservation Laws

We consider the discontinuous Galerkin (DG) finite element discretization of first order systems of conservation laws derivable as moments of the kinetic Boltzmann equation. This includes well known conservation law systems such as the Euler For the class of first order nonlinear conservation laws equipped with an entropy extension, an energy analysis of the DG method for the Cauchy initial value problem is developed. Using this DG energy analysis, several new variants of existing numerical flux functions are derived and shown to be energy stable.

Barth, Timothy↗
NASA NTRS1994

Three-Dimensional Navier-Stokes Method with Two-Equation Turbulence Models for Efficient Numerical Simulation of Hypersonic Flows

A new computational efficient 3-D compressible Reynolds-averaged implicit Navier-Stokes method with advanced two equation turbulence models for high speed flows is presented. All convective terms are modeled using an entropy satisfying higher-order Total Variation Diminishing (TVD) scheme based on implicit upwind flux-difference split approximations and arithmetic averaging procedure of primitive variables. This method combines the best features of data management and computational efficiency of space marching procedures with the generality and stability of time dependent Navier-Stokes procedures to solve flows with mixed supersonic and subsonic zones, including streamwise separated flows. Its robust stability derives from a combination of conservative implicit upwind flux-difference splitting with Roe's property U to provide accurate shock capturing capability that non-conservative schemes do not guarantee, alternating symmetric Gauss-Seidel 'method of planes' relaxation procedure coupled with a three-dimensional two-factor diagonal-dominant approximate factorization scheme, TVD flux limiters of higher-order flux differences satisfying realizability, and well-posed characteristic-based implicit boundary-point a'pproximations consistent with the local characteristics domain of dependence. The efficiency of the method is highly increased with Newton Raphson acceleration which allows convergence in essentially one forward sweep for supersonic flows. The method is verified by comparing with experiment and other Navier-Stokes methods. Here, results of adiabatic and cooled flat plate flows, compression corner flow, and 3-D hypersonic shock-wave/turbulent boundary layer interaction flows are presented. The robust 3-D method achieves a better computational efficiency of at least one order of magnitude over the CNS Navier-Stokes code. It provides cost-effective aerodynamic predictions in agreement with experiment, and the capability of predicting complex flow structures in complex geometries with good accuracy.

Bardina, J. E.↗
DOE OSTI2021

Analytic continuation of noisy data using Adams Bashforth residual neural network

We propose a data-driven learning framework for the analytic continuation problem in numerical quantum many-body physics. Designing an accurate and efficient framework for the analytic continuation of imaginary time using computational data is a grand challenge that has hindered meaningful links with experimental data. The standard Maximum Entropy (MaxEnt)-based method is limited by the quality of the computational data and the availability of prior information. Also, the MaxEnt is not able to solve the inversion problem under high level of noise in the data. Here we introduce a novel learning model for the analytic continuation problem using a Adams-Bashforth residual neural network (AB-ResNet). Additionally, the advantage of this deep learning network is that it is model independent and, therefore, does not require prior information concerning the quantity of interest given by the spectral function. More importantly, the ResNet-based model achieves higher accuracy than MaxEnt for data with higher level of noise. Finally, numerical examples show that the developed AB-ResNet is able to recover the spectral function with accuracy comparable to MaxEnt where the noise level is relatively small.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Graph neural networks predict energetic and mechanical properties for models of solid solution metal alloy phases

Here, we developed a PyTorch-based architecture called HydraGNN that implements graph convolutional neural networks (GCNNs) to predict the formation energy and the bulk modulus for models of solid solution alloys for various atomic crystal structures and relaxed volumes. We trained the GCNN surrogate model on a dataset for nickel–niobium (NiNb) generated by the embedded atom model (EAM) empirical interatomic potential for demonstration purposes. The dataset was generated by calculating the formation energy and the bulk modulus as a prototypical elastic property for optimized geometries starting from initial body-centered cubic (BCC), face-centered cubic (FCC), and hexagonal compact packed (HCP) crystal structures, with configurations spanning the possible compositional range for each of the three types of initial crystal structures. Numerical results show that the GCNN model effectively predicts both the formation energy and the bulk modulus as function of the optimized crystal structure, relaxed volume, and configurational entropy of the model structures for solid solution alloys.

36 MATERIALS SCIENCE↗
NASA NTRS1972

Numerical calculations of flow fields

The solutions to the equations of motion for inviscid fluid flow around a pointed elliptic cone at incidence are presented. The numerical method used, MacCormack's second order preferential predictor-corrector finite difference approximation, is applied to the fluid flow equations derived in conservation-law form. The entropy boundary condition, hitherto unused for elliptic cone problems, is investigated and compared to reflection boundary condition solutions. The stagnation streamline movement of the inclined elliptic cone is noted and surface pressure coefficients are plotted. Also presented are solutions for an elliptic cone and a circular cone at zero incidence and a circular cone at a small angle of attack. Comparisons are made between these present solutions and previously published theory.

Anderson, D. M.↗
DOE OSTI2023

The reflected entanglement spectrum for free fermions

We consider the reflected entropy and the associated entanglement spectrum for free fermions reduced to two intervals in 1 + 1 dimensions. Working directly in the continuum theory the reflected entropy can be extracted from the spectrum of a singular integral equation whose kernel is determined by the known free fermion modular evolved correlation function. We find the spectrum numerically and analytically in certain limits. For intervals that almost touch the reflected entanglement spectrum approaches the spectrum of the thermal density matrix. This suggests that the reflected entanglement spectrum is well suited to the task of extracting physical data of the theory directly from the ground state wave function.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1994

Inside the supernova: A powerful convective engine

We present an extensive study of the inception of supernova explosions by following the evolution of the cores of two massive stars (15 and 25 Solar mass) in multidimension. Our calculations begin at the onset of core collapse and stop several hundred milliseconds after the bounce, at which time successful explosions of the appropriate magnitude have been obtained. Similar to the classical delayed explosion mechanism of Wilson, the explosion is powered by the heating of the envelope due to neutrinos emitted by the protoneutron star as it radiates the gravitational energy liberated by the collapse. However, as was shown by Herant, Benz, & Colgate, this heating generates strong convection outside the neutrinosphere, which we demonstrate to be critical to the explosion. By breaking a purely stratified hydrostatic equilibrium, convection moves the nascent supernova away from a delicate radiative equilibrium between neutrino emission and absorption, Thus, unlike what has been observed in one-dimensional calculations, explosions are rendered quite insensitive to the details of the physical input parameters such as neutrino cross sections or nuclear equation of state parameters. As a confirmation, our comparative one-dimensional calculations with identical microphysics, but in which convection cannot occur, lead to dramatic failures. Guided by our numerical results, we have developed a paradigm for the supernova explosion mechanism. We view a supernova as an open cycle thermodynamic engine in which a reservoir of low-entropy matter (the envelope) is thermally coupled and physically connected to a hot bath (the protoneutron star) by a neutrino flux, and by hydrodynamic instabilities. This paradigm does not invoke new or modified physics over previous treatments, but relies on compellingly straightforward thermodynamic arguments. It provides a robust and self-regulated explosion mechanism to power supernovae that is effective under a wide range of physical parameters.

Herant, Marc↗
NASA NTRS1988

Numerical simulation of underexpanded plumes using upwind algorithms

The solution of the flow characteristics of highly underexpanded plumes is presently sought by means of a high-resolution total variation-diminishing upwind scheme featuring finite-volume metrics, Roe averaging with entropy correction, and monotonic upstream schemes for conservation laws. The existing approximate factorization schemes are used both as preconditioners and in order to apply the conjugate residual algorithm. The plumes treated are axisymmetric and laminar, and for an ideal gas that is allowed to be thermally incomplete. Simulations of the SSME and SRB plumes have been successfully conducted.

Obayashi, Shigeru↗
DOE OSTI2022

Numerical schemes for a multi-species BGK model with velocity-dependent collision frequency

Here, we consider a kinetic description of multi-species gas mixture modeled with Bhatnagar-Gross-Krook (BGK) collision operators, in which the collision frequency varies not only in time and space but also with the microscopic velocity. In this model, the Maxwellians typically used in standard BGK operators are replaced by a generalization of such target functions, which are defined by a variational procedure. In this paper we present a numerical method for simulating this model, which uses an Implicit-Explicit (IMEX) scheme to minimize a certain potential function, mimicking the Lagrange functional that appears in the theoretical derivation. We show that theoretical properties such as conservation of mass, total momentum and total energy as well as positivity of the distribution functions are preserved by the numerical method, and illustrate its usefulness and effectiveness with numerical examples.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Universal Bound to the Amplitude of the Vortex Nernst Signal in Superconductors

A liquid of superconducting vortices generates a transverse thermoelectric response. This Nernst signal has a tail deep in the normal state due to superconducting fluctuations. Here, we present a study of the Nernst effect in two-dimensional heterostructures of Nb-doped strontium titanate (STO) and in amorphous MoGe. Additionally, the Nernst signal generated by ephemeral Cooper pairs above the critical temperature has the magnitude expected by theory in STO. On the other hand, the peak amplitude of the vortex Nernst signal below Tc is comparable in both and in numerous other superconductors despite the large distribution of the critical temperature and the critical magnetic fields. In four superconductors belonging to different families, the maximum Nernst signal corresponds to an entropy per vortex per layer of ≈kBln 2 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗