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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
NASA NTRS2012

The Outer Limits of Galaxy Clusters: Observations to the Virial Radius with Suzaku, XMM,and Chandra

The outskirts of galaxy clusters, near the virial radius, remain relatively unexplored territory and yet are vital to our understanding of cluster growth, structure, and mass. In this presentation, we show the first results from a program to constrain the sate of the outer intra-cluster medium (ICM) in a large sample of galaxy clusters, exploiting the strengths of three complementary X-ray observatories: Suzaku (low, stable background), XMM-Newton (high sensitivity),and Chandra (good spatial resolution). By carefully combining observations from the cluster core to beyond r200, we are able to identify and reduce systematic uncertainties that would impede our spatial and spectral analysis using a single telescope. Our sample comprises nine clusters at z is approximately 0.1-0.2 fully covered in azimuth to beyond r200, and our analysis indicates that the ICM is not in hydrostatic equilibrium in the cluster outskirts, where we see clear azimuthal variations in temperature and surface brightness. In one of the clusters, we are able to measure the diffuse X-ray emission well beyond r200, and we find that the entropy profile and the gas fraction are consistent with expectations from theory and numerical simulations. These results stand in contrast to recent studies which point to gas clumping in the outskirts; the extent to which differences of cluster environment or instrumental effects factor in this difference remains unclear. From a broader perspective, this project will produce a sizeable fiducial data set for detailed comparison with high-resolution numerical simulations.

Miller, Eric D.↗
NASA NTRS2022

Interaction of a Tunnel-like Acoustic Disturbance Field with a Blunt Cone Boundary Layer at Mach 8

The existing measurements of laminar-to-turbulent transition over circular cones in conventional (i.e., “noisy”) hypersonic wind tunnels have established that the transition location moves downstream when the nose radius is increased from zero. However, this initially downstream movement slows down and ultimately reverses beyond a critical value of the nose radius, and may be related to external forcing in the form of freestream disturbances and/or surface roughness. To understand the effects of freestream acoustic disturbances on transition reversal over a blunt body, hypersonic boundary-layer receptivity to broadband freestream acoustic disturbances from the nozzle wall of a digital conventional wind tunnel is investigated by both direct numerical simulations (DNS) and modal and nonmodal stability analysis. A Mach 8 flow over a 7 deg half-angle cone with a nose radius of 𝑅 𝑛 = 5.2mm and freestream Reynolds number of 12.2 × 10 6 m -1 is considered. The results show that the broadband tunnel noise in the free stream of a convectional hypersonic wind tunnel (i.e., outside of the nozzle-wall turbulent boundary layer) can be well represented by an acoustic model with an ansatz of slow acoustic waves. With successful calibration of the model parameters against the precursor tunnel DNS, such an acoustic ansatz can successfully reproduce both the frequency-wavenumber spectra and the temporal evolution of the broadband tunnel noise radiated from the nozzle wall. Additionally, the DNS of the Mach 8 blunt cone with tunnel-like acoustic input above the bow shock showed that the spectra of wall-pressure and heat-transfer fluctuations recovers the signature of the axisymmetric waves predicted by the nonmodal analysis. Furthermore, the azimuthal wavenumber and frequency spectrum of the temperature fluctuations as a function of the wall-normal distance show higher amplitudes for three-dimensional waves above the boundary-layer edge. The numerical schlieren contours show the inclined structures commonly observed in blunt cone experiments, demonstrating that they correspond to three-dimensional structures due to freestream disturbances in the presence of an entropy layer.

Boundary-layer transition↗
DOE OSTI2024

Magnetism in the two-dimensional dipolar XY model

Here, motivated by a recent experiment on a square-lattice Rydberg atom array realizing a long-range dipolar XY model [C. Chen et al., Nature (London) 616, 691 (2023)], we numerically study the model's equilibrium properties. We obtain the phase diagram, critical properties, entropies, variance of the magnetization, and site-resolved correlation functions. We consider both ferromagnetic and antiferromagnetic interactions and apply quantum Monte Carlo and pseudo-Majorana functional renormalization group techniques, generalizing the latter to a U⁡(1) symmetric setting. Our simulations perform extensive thermometry in dipolar Rydberg atom arrays and establish conditions for adiabaticity and thermodynamic equilibrium. On the ferromagnetic side of the experiment, we determine the entropy per particle S/N≈0.5, close to the one at the critical temperature, S c /N=0.585⁢(15). The simulations suggest the presence of an out-of-equilibrium plateau at large distances in the correlation function, thus motivating future studies on the nonequilibrium dynamics of the system.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Black-box statistical prediction of lossy compression ratios for scientific data

Lossy compressors are increasingly adopted in scientific research, tackling volumes of data from experiments or parallel numerical simulations and facilitating data storage and movement. In contrast with the notion of entropy in lossless compression, no theoretical or data-based quantification of lossy compressibility exists for scientific data. Users rely on trial and error to assess lossy compression performance. As a strong data-driven effort toward quantifying lossy compressibility of scientific datasets, we provide a statistical framework to predict compression ratios of lossy compressors. Our method is a two-step framework where (i) compressor-agnostic predictors are computed and (ii) statistical prediction models relying on these predictors are trained on observed compression ratios. Proposed predictors exploit spatial correlations and notions of entropy and lossyness via the quantized entropy. We study 8+ compressors on 6 scientific datasets and achieve a median percentage prediction error less than 12%, which is substantially smaller than that of other methods while achieving at least a 8.8× speedup for searching for a specific compression ratio and 7.8× speedup for determining the best compressor out of a collection.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1995

Investigation of computational and spectral analysis methods for aeroacoustic wave propagation

Most computational fluid dynamics (CFD) schemes are not adequately accurate for solving aeroacoustics problems, which have wave amplitudes several orders of magnitude smaller yet with frequencies larger than the flow field variations generating the sound. Hence, a computational aeroacoustics (CAA) algorithm should have minimal dispersion and dissipation features. A dispersion relation preserving (DRP) scheme is, therefore, applied to solve the linearized Euler equations in order to simulate the propagation of three types of waves, namely: acoustic, vorticity, and entropy waves. The scheme is derived using an optimization procedure to ensure that the numerical derivatives preserve the wave number and angular frequency of the partial differential equations being discretized. Consequently, simulated waves propagate with the correct wave speeds and exhibit their appropriate properties. A set of radiation and outflow boundary conditions, compatible with the DRP scheme and derived from the asymptotic solutions of the governing equations, are also implemented. Numerical simulations are performed to test the effectiveness of the DRP scheme and its boundary conditions. The computed solutions are shown to agree favorably with the exact solutions. The major restriction appears to be that the dispersion relations can be preserved only for waves with wave lengths longer than four or five spacings. The boundary conditions are found to be transparent to the outgoing disturbances. However, when the disturbance source is placed closer to a boundary, small acoustic reflections start appearing. CAA generates enormous amounts of temporal data which needs to be reduced to understand the physical problem being simulated. Spectral analysis is one approach that helps us in extracting information which often can not be easily interpreted in the time domain. Thus, three different methods for the spectral analysis of numerically generated aeroacoustic data are studied. First, the capabilities of two traditional methods for spectral analysis, namely, the Blackman-Tukey method and periodogram method, are compared in estimating the spectra of a simple-periodic process. The periodogram is then applied to analyze transitory-deterministic processes. Finally, these two methods are compared with a more recent method, referred as the Weighted-Overlapped-Segment-Averaging (WOSA) method, in estimating the spectra of a chaotic (random-like) process. From the demonstrative case for the spectral analyses of data generated by simple-periodic process, the periodogram method is found to give a better estimate of the steep-sloped spectra than the Blackman-Tukey method. Also, for this problem, the Hanning window is found to perform better with the periodogram method than with the Blackman-Tukey method. Finally, for the spectral analysis of data generated by the chaotic process, the periodogram method does not perform well, whereas, the WOSA and Blackman-Tukey methods give equivalently good results.

Vanel, Florence O.↗
NASA NTRS2001

High Order Difference Method for Low Mach Number Aeroacoustics

A high order finite difference method with improved accuracy and stability properties for computational aeroacoustics (CAA) at low Mach numbers is proposed. The Euler equations are split into a conservative and a symmetric non- conservative portion to allow the derivation of a generalized energy estimate. Since the symmetrization is based on entropy variables, that splitting of the flux derivatives is referred to as entropy splitting. Its discretization by high order central differences was found to need less numerical dissipation than conventional conservative schemes. Owing to the large disparity of acoustic and stagnation quantities in low Mach number aeroacoustics, the split Euler equations are formulated in perturbation form. The unknowns are the small changes of the conservative variables with respect to their large stagnation values. All nonlinearities and the conservation form of the conservative portion of the split flux derivatives can be retained, while cancellation errors are avoided with its discretization opposed to the conventional conservative form. The finite difference method is third-order accurate at the boundary and the conventional central sixth-order accurate stencil in the interior. The difference operator satisfies the summation by parts property analogous to the integration by parts in the continuous energy estimate. Thus, strict stability of the difference method follows automatically. Spurious high frequency oscillations are suppressed by a characteristic-based filter similar to but without limiter. The time derivative is approximated by a 4-stage low-storage second-order explicit Runge-Kutta method. The method has been applied to simulate vortex sound at low Mach numbers. We consider the Kirchhoff vortex, which is an elliptical patch of constant vorticity rotating with constant angular frequency in irrotational flow. The acoustic pressure generated by the Kirchhoff vortex is governed by the 2D Helmholtz equation, which can be solved analytically using separation of variables.

Mueller, B.↗
NASA NTRS2019

Efficient Entropy Stable Gauss Collocation Methods

The construction of high-order entropy stable collocation schemes on quadrilateral and hexahedral elements has relied on the use of Gauss-Legendre-Lobatto collocation points and their equivalence with summation-by-parts (SBP) finite difference operators. In this work, we show how to efficiently generalize the construction of semidiscrete, entropy stable schemes on tensor product elements to Gauss points and generalized SBP operators. Numerical experiments suggest 8 that the use of Gauss points significantly improves accuracy on curved meshes.

Jesse Chan↗
DOE OSTI2025

The superconformal index and black hole instabilities

The superconformal index of $\mathcal{N}$ = 4 supersymmetric Yang-Mills theory with gauge group U(N) has provided powerful insights into the entropy of supersymmetric black holes in AdS 5 × S 5 , including some sub-leading logarithmic and non-perturbative corrections. Recently, the phase space of supersymmetric solutions has been argued to contain configurations other than the asymptotically AdS 5 black hole. Such configurations include the so-called grey galaxies where the black hole at the center is surrounded by a gas of gravitons. By numerically evaluating the superconformal index of $\mathcal{N}$ = 4 supersymmetric Yang-Mills at small values of N, we detect systematic deviations from the entropy of black holes with two distinct angular momenta. We find that the giant graviton expansion of the index is a numerically efficient way of evaluating the index that complements the direct character evaluation and allows for explicit access to N ≤ 15 with up to two giant gravitons in the expansion. We find it remarkable that a supersymmetric quantity in field theory, usually thought of as a rigid counting observable, indeed contains information about different phases in the space of supersymmetric solutions on the gravity side.

AdS-CFT correspondence↗
DOE OSTI2020

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

Benchmark data are presented for the zeroth- through third-order many-body perturbation corrections to the electronic Helmholtz energy, internal energy, and entropy in the canonical ensemble in a wide range of temperature. They are determined as numerical λ-derivatives of the respective quantities computed by thermal full configuration interaction with a perturbation-scaled Hamiltonian, $\hat{H} = \hat{H}_0 + λ\hat{V}$. Sum-over-states analytical formulas for up to the third-order corrections to these properties are also derived as analytical λ-derivatives. These formulas, which are verified by exact numerical agreement with the benchmark data, are given in terms of the Hirschfelder–Certain degenerate perturbation energies and should be valid for both degenerate and nondegenerate reference states at any temperature down to zero. Further, the results in the canonical ensemble are compared with the same in the grand canonical ensemble.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Biased degenerate ground-state sampling of small Ising models with converged quantum approximate optimization algorithm

The quantum alternating operator ansatz, a generalization of the quantum approximate optimization algorithm (QAOA), is a quantum algorithm used for approximately solving combinatorial optimization problems. QAOA typically uses the transverse field mixer as the driving Hamiltonian. One of the interesting properties of the transverse field driving Hamiltonian is that it results in nonuniform sampling of degenerate ground states of optimization problems. In this study, we numerically examine the fair sampling properties of the transverse field mixer QAOA, and Grover mixer QAOA (GM-QAOA), which provides theoretical guarantees of fair sampling of degenerate optimal solutions, up to a large enough p such that the mean expectation value converges to an optimal approximation ratio of 1. This comparison is performed with high-quality heuristically computed, but not necessarily optimal, QAOA angles, which give strictly monotonically improving solution quality as p increases. These angles are computed using the Julia based numerical simulation software JuliQAOA. Fair sampling of degenerate ground states is quantified using the Shannon entropy of the ground-state amplitudes distribution. The fair sampling properties are reported on several quantum signature Hamiltonians from previous quantum annealing fair sampling studies. Small random fully connected spin glasses are shown, which exhibit exponential suppression of some degenerate ground states with transverse field mixer QAOA. The transverse field mixer QAOA simulations show that some problem instances clearly saturate the Shannon entropy of 0 with a maximally biased distribution that occurs when the learning converges to an approximation ratio of 1 while other problem instances never deviate from a maximum Shannon entropy (uniform distribution) at any p step. Published by the American Physical Society 2025

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Enhanced Estimation of Quantum Properties with Common Randomized Measurements

We present a technique for enhancing the estimation of quantum state properties by incorporating approximate prior knowledge about the quantum state. This consists in performing randomized measurements on a quantum processor and comparing the results with those obtained from a classical computer that stores an approximation of the quantum state. We provide unbiased estimators for expectation values of multicopy observables and present performance guarantees in terms of variance bounds that depend on the prior knowledge accuracy. We demonstrate the effectiveness of our approach through experimental and numerical examples detecting mixed-state entanglement, and estimating polynomial approximations of the von Neumann entropy and state fidelities. Published by the American Physical Society 2024

Physics↗
NASA NTRS1986

The discrete one-sided Lipschitz condition for convex scalar conservation laws

Physical solutions to convex scalar conservation laws satisfy a one-sided Lipschitz condition (OSLC) that enforces both the entropy condition and their variation boundedness. Consistency with this condition is therefore desirable for a numerical scheme and was proved for both the Godunov and the Lax-Friedrichs scheme--also, in a weakened version, for the Roe scheme, all of them being only first order accurate. A new, fully second order scheme is introduced here, which is consistent with the OSLC. The modified equation is considered and shows interesting features. Another second order scheme is then considered and numerical results are discussed.

Brenier, Yann↗
NASA NTRS2020

Entropy Stable Nonconforming Discretizations with the Summation-By-Parts Property for Curvilinear Coordinates

The entropy conservative/stable algorithm of Friedrichet al.(2018) for hyperbolic conservation laws on nonconforming p-refined/coarsened Cartesian grids, is extended to curvilinear grids for the compressible Euler equations. The primary focus is on constructing appropriate coupling procedures across the curvilinear nonconforming interfaces. A simple and flexible approach is proposed that uses interpolation operators from one element to the other. On the element faces,the analytic metrics are used to construct coupling terms, while metric terms in the volume are approximated to satisfy a discretization of the geometric conservation laws. The resulting scheme is entropy conservative/stable, elementwise conservative, and freestream preserving. The accuracy and stability properties of the resulting numerical algorithm are shown to be comparable to those ofthe original conforming scheme (∼p+ 1 convergence) in the context of the isentropic Euler vortex and the inviscid Taylor–Green vortex problems on manufactured high order grids.

David Del Rey Fernandez↗
DOE OSTI2021

First-principles calculation of lattice distortions in four single phase high entropy alloys with experimental validation

Exceptional properties of high entropy alloys (HEAs) are attributed to the disordered random solid solution of multiple alloying elements. Despite numerous studies, the fundamental understanding at the electronic and atomic levels is still missing. We report a comparative study of four fcc HEAs of NiFeCoCr and NiFeCoCrX (X = Mn, Cu, or Pd) based on ab initio calculations using large supercells with 500 atoms in equal composition. After fully optimizing their structures using the VASP package, their electronic structure, interatomic bonding, partial charge distribution, and mechanical properties are calculated and compared, revealing the intricate interdependence among them. A novel parameter based on the quantum mechanical metric for internal cohesion, the total bond order density (TBOD), is used to interpret the calculated properties. The highest TBOD is found in Cantor alloy NiFeCoCrMn but lower in NiFeCoCrPd. The atomic radii vary depending on their local chemical environment. The resulting lattice distortions is validated experimentally in NiFeCoCrMn and NiFeCoCrPd. Moreover, modeling of Cu and Pd clustering in the supercell shows they have lower total energy in agreement with the observation of Cu-enhanced nano-participates in NiFeCoCrCu and Pd-induced concentration wave in NiFeCoCrPd HEAs.

36 MATERIALS SCIENCE↗
NASA NTRS1979

Feature selection via entropy minimization: An example using LANDSAT satellite data

The author has identified the following significant results. The minimum entropy model may provide several useful advantages over traditional techniques for processing LANDSAT data. Total computer time to conduct a complete pattern recognition process is reduced. Subjective (transformed image), as well as statistically derived information is made available to the analyst/user much earlier in the analysis process. A rapid feedback loop in which numerous training set combinations can be tested for difference and representativeness is available. Additional tests of LANDSAT data processing using the minimum entropy model are clearly justified.

Zandonella, A.↗
NASA NTRS2017

The Root Cause of the Overheating Problem

Previously we identified the receding flow, where two fluid streams recede from each other, as an open numerical problem, because all well-known numerical fluxes give an anomalous temperature rise, thus called the overheating problem. This phenomenon, although presented in several textbooks, and many previous publications, has scarcely been satisfactorily addressed and the root cause of the overheating problem not well understood. We found that this temperature rise was solely connected to entropy rise and proposed to use the method of characteristics to eradicate the problem. However, the root cause of the entropy production was still unclear. In the present study, we identify the cause of this problem: the entropy rise is rooted in the pressure flux in a finite volume formulation and is implanted at the first time step. It is found theoretically inevitable for all existing numerical flux schemes used in the finite volume setting, as confirmed by numerical tests. This difficulty cannot be eliminated by manipulating time step, grid size, spatial accuracy, etc, although the rate of overheating depends on the flux scheme used. Finally, we incorporate the entropy transport equation, in place of the energy equation, to ensure preservation of entropy, thus correcting this temperature anomaly. Its applicability is demonstrated for some relevant 1D and 2D problems. Thus, the present study validates that the entropy generated ab initio is the genesis of the overheating problem.

Liou, Meng-Sing↗
NASA NTRS1993

Aspects of a high-resolution scheme for the Navier-Stokes equations

In this paper we emphasize the importance of the form of the numerical dissipation model in computing accurate viscous flow solutions. A high-resolution scheme for viscous flows based on three-point central differencing and a matrix dissipation is considered. The various components of this scheme, including 'entropy fix', limiter function, and boundary-point dissipation are discussed. By analyzing boundary-point dissipation stencils, we confirm that with the matrix dissipation model the normal numerical dissipation terms in the streamwise momentum equation are independent of the Reynolds number. Such independence is not achieved with a scalar dissipation form. The accuracy of the central-difference scheme, with and without matrix dissipation, and the flux-difference split scheme of Roe, which is classified as a high-resolution scheme, is compared. For this comparison, three high Reynolds number laminar flows are considered. Solutions of the Navier-Stokes equations are obtained for low-speed flow over a flat plate, transonic flow over an airfoil with transition near the leading edge, and hypersonic flow over a compression ramp. The emphasis of the comparison is primarily on the details of the viscous flows. The necessity of the high-resolution property is revealed.

Swanson, R. C.↗
NASA NTRSDate not supplied

Hypersonic Boundary-Layer Transition over a Blunt Circular Cone in a Mach 8 Digital Wind Tunnel

To understand the effects of freestream acoustic disturbances on transition reversal over a blunt body, three-dimensional direct numerical simulations (DNS) were conducted to investigate the interaction of a tunnel-like acoustic disturbance field with a Mach 8 laminar boundary layer over a 7 deg half-angle blunt cone with a nose radius of 𝑅 𝑛 = 5.2 mm. The flow conditions and the axisymmetric cone geometry matched those measured in the Sandia Hypersonic Wind Tunnel at Mach 8 (Sandia HWT-8), and the incident freestream acoustic disturbances in the DNS were generated by acoustic radiation from the nozzle-wall turbulent boundary layer in the same tunnel. Consistent with the predictions of the nonmodal instability analysis, the study shows the emergence of a second spectral peak for moderately oblique waves between 150 kHz and 200 kHz in the aft part of the cone. The spectral energy of temperature and streamwise velocity fluctuations is initially concentrated within the entropy layer. These fluctiations are eventually “swallowed” by the boundary layer as the entropy-layer thickness 𝛿 𝑆 becomes similar to the boundary-layer thickness 𝛿h𝑡 at downstream locations. The numerical schlieren contours show the same inclined structures commonly observed in blunt cone experiments, and these structures are found to be three dimensional and advect with a constant speed approximately equal to the mean flow speed near the edge of the boundary layer.

Hypersonic Flow↗