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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
NASA NTRS1988

An improved version of NCOREL: A computer program for 3-D nonlinear supersonic potential flow computations

A computer code called NCOREL (for Nonconical Relaxation) has been developed to solve for supersonic full potential flows over complex geometries. The method first solves for the conical at the apex and then marches downstream in a spherical coordinate system. Implicit relaxation techniques are used to numerically solve the full potential equation at each subsequent crossflow plane. Many improvements have been made to the original code including more reliable numerics for computing wing-body flows with multiple embedded shocks, inlet flow through simulation, wake model and entropy corrections. Line relaxation or approximate factorization schemes are optionally available. Improved internal grid generation using analytic conformal mappings, supported by a simple geometric Harris wave drag input that was originally developed for panel methods and internal geometry package are some of the new features.

Siclari, Michael J.↗
NASA NTRS1993

Energy considerations in computational aeroacoustics

A finite-volume multistage time-stepping Euler code is used to investigate the use of CFD algorithms for the direct calculation of acoustics. The 2D compressible inviscid flow about an accelerating or decelerating circular cylinder is used as a model problem. The time evolution of the energy transfer from the cylinder to the fluid, as the cylinder is moved from rest to some nonnegligible velocity, is clearly seen. By examining the temporal and spatial characteristics of the numerical solution, a distinction can be made between the propagating acoustic energy, the convecting energy associated with the entropy change in the fluid, and the energy contained in the local aerodynamic field. Systematic variation of the cylinder acceleration shows that the radiated acoustic energy depends strongly upon the rate of acceleration or deceleration. The computational grid has a large effect on the ratio of acoustic energy to nonphysical entropy associated energy, while the role of the explicit artificial viscosity seems to be of second order. The entropy term was nearly negligible in all cases the cylinder was started slowly.

Brentner, Kenneth S.↗
DOE OSTI2024

Random projection using random quantum circuits

The random sampling task performed by Google's Sycamore processor gave us a glimpse of the “quantum supremacy era.” This has definitely shed some light on the power of random quantum circuits in this abstract task of sampling outputs from the (pseudo)random circuits. In this paper, we explore a practical near-term use of local random quantum circuits in dimensional reduction of large low-rank data sets. We make use of the well-studied dimensionality reduction technique called the random projection method. This method has been extensively used in various applications such as image processing, logistic regression, entropy computation of low-rank matrices, etc. We prove that the matrix representations of local random quantum circuits with sufficiently shorter depths [ ∼ O ( n ) ] serve as good candidates for random projection. We demonstrate numerically that their projection abilities are not far off from the computationally expensive classical principal components analysis on MNIST and CIFAR-100 image datasets. We also benchmark the performance of quantum random projection against the commonly used classical random projection in the tasks of dimensionality reduction of image data sets and computing von Neumann entropies of large low-rank density matrices. And finally, using variational quantum singular value decomposition, we demonstrate a near-term implementation of extracting the singular vectors with dominant singular values after quantum random projecting a large low-rank matrix to lower dimensions. All such numerical experiments unequivocally demonstrate the ability of local random circuits to randomize a large Hilbert space at sufficiently shorter depths with robust retention of properties of large data sets in reduced dimensions. Published by the American Physical Society 2024

Kumaran, Keerthi (ORCID:0009000949125721)↗
DOE OSTI2023

Inferring Quantum Network Topology Using Local Measurements

Statistical correlations that can be generated across the nodes in a quantum network depend crucially on its topology. However, this topological information might not be known a priori, or it may need to be verified. In this paper, we propose an efficient protocol for distinguishing and inferring the topology of a quantum network. We leverage entropic quantities-namely, the von Neumann entropy and the measured mutual information-as well as measurement covariance to uniquely characterize the topology. We show that the entropic quantities are sufficient to distinguish two networks that prepare GHZ states. Moreover, if qubit measurements are available, both entropic quantities and covariance can be used to infer the network topology without state-preparation assumptions. We show that the protocol can be entirely robust to noise and can be implemented via quantum variational optimization. Numerical experiments on both classical simulators and quantum hardware show that covariance is generally more reliable for accurately and efficiently inferring the topology, whereas entropy-based methods are often better at identifying the absence of entanglement in the low-shot regime.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

The interpretation of temperature and salinity variables in numerical ocean model output and the calculation of heat fluxes and heat content

Abstract. The international Thermodynamic Equation of Seawater 2010 (TEOS-10) defined the enthalpy and entropy of seawater, thus enabling the global ocean heat content to be calculated as the volume integral of the product of in situ density, ρ, and potential enthalpy, h0 (with reference sea pressure of 0 dbar). In terms of Conservative Temperature, Θ, ocean heat content is the volume integral of ρcp0Θ, where cp0 is a constant “isobaric heat capacity”. However, many ocean models in the Coupled Model Intercomparison Project Phase 6 (CMIP6) as well as all models that contributed to earlier phases, such as CMIP5, CMIP3, CMIP2, and CMIP1, used EOS-80 (Equation of State – 1980) rather than the updated TEOS-10, so the question arises of how the salinity and temperature variables in these models should be physically interpreted, with a particular focus on comparison to TEOS-10-compliant observations. In this article we address how heat content, surface heat fluxes, and the meridional heat transport are best calculated using output from these models and how these quantities should be compared with those calculated from corresponding observations. We conclude that even though a model uses the EOS-80, which expects potential temperature as its input temperature, the most appropriate interpretation of the model's temperature variable is actually Conservative Temperature. This perhaps unexpected interpretation is needed to ensure that the air–sea heat flux that leaves and arrives in atmosphere and sea ice models is the same as that which arrives in and leaves the ocean model. We also show that the salinity variable carried by present TEOS-10-based models is Preformed Salinity, while the salinity variable of EOS-80-based models is also proportional to Preformed Salinity. These interpretations of the salinity and temperature variables in ocean models are an update on the comprehensive Griffies et al. (2016) paper that discusses the interpretation of many aspects of coupled Earth system models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

A multilayer multi-configurational approach to efficiently simulate large-scale circuit-based quantum computers on classical machines

Here, a multilayer multi-configurational theory framework is adapted to simulate circuit-based quantum computers. Quantum addition of superpositions of an exponential number of summands is performed in polynomial time with high accuracy. We demonstrate numerically accurate calculations including up to one million qubits for entangling benchmarks. Simulation cost can be assessed by entropy-based entanglement measures. For the considered systems, we show that the entanglement only grows weakly with the system size. The present simulations demonstrate how quantum algorithms in low-entropy regimes can be used efficiently on classically simulated quantum computers.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Conservative discontinuous Galerkin scheme of a gyro-averaged Dougherty collision operator

Here, a conservative discontinuous Galerkin scheme for a non-linear Dougherty collision operator in full-f long-wavelength gyrokinetics is presented. Analytically this model operator has the advective-diffusive form of Fokker-Planck operators, it has a non-decreasing entropy functional, and conserves particles, momentum and energy. Discretely these conservative properties are maintained exactly as well, independent of numerical resolution. In this work the phase space discretization is performed using a novel version of the discontinuous Galerkin scheme, carefully constructed using concepts of weak equality and recovery. Discrete time advancement is carried out with an explicit time-stepping algorithm, whose stability limits we explore. The formulation and implementation within the long-wavelength gyrokinetic solver of Gkeyll are validated with relaxation tests, collisional Landau-damping benchmarks and the study of 5D gyrokinetic turbulence on helical, open field lines.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1998

Self-Similar Compressible Free Vortices

Lie group methods are used to find both exact and numerical similarity solutions for compressible perturbations to all incompressible, two-dimensional, axisymmetric vortex reference flow. The reference flow vorticity satisfies an eigenvalue problem for which the solutions are a set of two-dimensional, self-similar, incompressible vortices. These solutions are augmented by deriving a conserved quantity for each eigenvalue, and identifying a Lie group which leaves the reference flow equations invariant. The partial differential equations governing the compressible perturbations to these reference flows are also invariant under the action of the same group. The similarity variables found with this group are used to determine the decay rates of the velocities and thermodynamic variables in the self-similar flows, and to reduce the governing partial differential equations to a set of ordinary differential equations. The ODE's are solved analytically and numerically for a Taylor vortex reference flow, and numerically for an Oseen vortex reference flow. The solutions are used to examine the dependencies of the temperature, density, entropy, dissipation and radial velocity on the Prandtl number. Also, experimental data on compressible free vortex flow are compared to the analytical results, the evolution of vortices from initial states which are not self-similar is discussed, and the energy transfer in a slightly-compressible vortex is considered.

vonEllenrieder, Karl↗
NASA NTRS2015

DNS of Flow in a Low-Pressure Turbine Cascade Using a Discontinuous-Galerkin Spectral-Element Method

A new computational capability under development for accurate and efficient high-fidelity direct numerical simulation (DNS) and large eddy simulation (LES) of turbomachinery is described. This capability is based on an entropy-stable Discontinuous-Galerkin spectral-element approach that extends to arbitrarily high orders of spatial and temporal accuracy and is implemented in a computationally efficient manner on a modern high performance computer architecture. A validation study using this method to perform DNS of flow in a low-pressure turbine airfoil cascade are presented. Preliminary results indicate that the method captures the main features of the flow. Discrepancies between the predicted results and the experiments are likely due to the effects of freestream turbulence not being included in the simulation and will be addressed in the final paper.

Low-Pressure Turbine Cascade↗
NASA NTRS1991

A two-dimensional Euler solution for an unbladed jet engine configuration

A two dimensional nonaxisymmetric Euler solution in a geometric representation of a jet engine configuration without blades is presented. The domain, including internal and external flow, is covered with a multiblock grid. To construct the grid, a domain decomposition technique is used to subdivide the domain and smooth grids are dimensioned and placed in each block. The grid contains 44 blocks which cover the external field, the inlet, bypass duct, core duct and nozzle of the nonaxisymmetric engine configuration. The geometry is symmetric about the meanline of the hub, but the grid is not since there is no symmetry condition applied to the grid between the two halves. With a symmetric grid at zero angle of attack, the measures of the solution would cancel exactly. With an asymmetric grid, the solution will not necessarily be symmetric and the lift coefficient will not necessarily be zero. Thus, grid asymmetry can be exploited to verify the resolution of the solution. The solution may be verified on the basis of five theoretical quantities: conservation of mass and energy, deviation of the lift coefficient from zero, deviation of the drag coefficient from zero, deviations from constant entropy, and deviations in the pressure distributions over the symmetric surfaces of the components. This technique is suitable for obtaining numerical solutions in complex geometries and provides a foundation for complete engine throughflow calculations.

Stewart, Mark E. M.↗
NASA NTRS1977

Sonic boom analysis for high-altitude flight at high Mach number

Numerical programs for the computation of the flow field from the airplane at the flight altitude to the ground are presented. They take into account the nonlinear effects of high Mach number, the entropy change across the shock, the entropy and enthalpy variations in the atmospheric layer, and the gravitational effect. Extension of the programs for the axisymmetric problems to handle nonaxisymmetric terms is described. The asymmetry can be caused by the geometry of the body and the lift, and also by the fact that the variations in the atmospheric layer are two-dimensional. Numerical results for ground level signatures of several configurations at various flight conditions are presented and compared with existing approximate theories to demonstrate the influences of these nonlinear effects.

Ferri, A.↗
NASA NTRS1985

Steady inviscid three-dimensional flows

The present analysis combines some of the theoretical concepts suggested by Hawthorne (1955) with a numerical integration procedure suggested by Martin (1978). The resulting algorithm is for inviscid subsonic flows. Thus, it is restricted to high Reynolds number flows. Chang and Adamczyk (1983) have provided a detailed derivation of the present algorithm along with a discussion of its stability bounds. The present paper represents a summary of this work. The integration of the continuity equation is considered along with an evaluation of the entropy, total temperature, and vorticity field. Attention is given to the shear-flow algorithm construction, and an application to a shear flow in a turning channel. A description of numerical results is also provided. The discussed algorithm represents a new procedure for solving inviscid subsonic three-dimensional rotational flows.

Adamczyk, J. J.↗
DOE OSTI2020

Quantifying Information without Entropy: Identifying Intermittent Disturbances in Dynamical Systems

A system’s response to disturbances in an internal or external driving signal can be characterized as performing an implicit computation, where the dynamics of the system are a manifestation of its new state holding some memory about those disturbances. Identifying small disturbances in the response signal requires detailed information about the dynamics of the inputs, which can be challenging. This paper presents a new method called the Information Impulse Function (IIF) for detecting and time-localizing small disturbances in system response data. The novelty of IIF is its ability to measure relative information content without using Boltzmann’s equation by modeling signal transmission as a series of dissipative steps. Since a detailed expression of the informational structure in the signal is achieved with IIF, it is ideal for detecting disturbances in the response signal, i.e., the system dynamics. Those findings are based on numerical studies of the topological structure of the dynamics of a nonlinear system due to perturbated driving signals. The IIF is compared to both the Permutation entropy and Shannon entropy to demonstrate its entropy-like relationship with system state and its degree of sensitivity to perturbations in a driving signal.

42 ENGINEERING↗
NASA NTRS1987

Vortex shedding in compressor blade wakes

The wakes of highly loaded axial compressor blades were often considered to be turbulent, unstructured flows. Recent work has suggested that the blade wakes are in fact dominated by a vortex street-like structure. The work on the wake structure at MIT is reviewed, the results of a viscous numerical simulation are presented, the blade wake vortices are compared to those shed from a cylinder, and the implications of the wake structure on compressor performance are discussed. In particular, a two-dimensional, time accurate, viscous calculation shows both a periodic wake structure and time variations in the passage shock strength. The numerical calculations are compared to laser anemometer and high frequency response probe data. The effect of the wake structure on the entropy production and apparent adiabatic efficiency of the compressor rotor is discussed.

Epstein, A. H.↗
NASA NTRS1992

Essentially nonoscillatory postprocessing filtering methods

High order accurate centered flux approximations used in the computation of numerical solutions to nonlinear partial differential equations produce large oscillations in regions of sharp transitions. Here, we present a new class of filtering methods denoted by Essentially Nonoscillatory Least Squares (ENOLS), which constructs an upgraded filtered solution that is close to the physically correct weak solution of the original evolution equation. Our method relies on the evaluation of a least squares polynomial approximation to oscillatory data using a set of points which is determined via the ENO network. Numerical results are given in one and two space dimensions for both scalar and systems of hyperbolic conservation laws. Computational running time, efficiency, and robustness of method are illustrated in various examples such as Riemann initial data for both Burgers' and Euler's equations of gas dynamics. In all standard cases, the filtered solution appears to converge numerically to the correct solution of the original problem. Some interesting results based on nonstandard central difference schemes, which exactly preserve entropy, and have been recently shown generally not to be weakly convergent to a solution of the conservation law, are also obtained using our filters.

Lafon, F.↗
NASA NTRS2005

Shock Melting of Permafrost on Mars: Water Ice Multiphase Equation of State for Numerical Modeling and Its Testing

The presence of water/ice/brine in upper layers of Martian crust affects many processes of impact cratering. Modeling of these effects promises better understanding of Martian cratering records. We present here the new ANEOS-based multiphase equation of state for water/ice constructed for usage in hydrocodes and first numerical experiments on permafrost shock melting. Preliminary results show that due to multiple shock compression of ice inclusions in rocks the entropy jump in shocked ice is smaller than in pure ice for the same shock pressure. Hence previous estimates of ice melting during impact cratering on Mars should be re-evaluated. Additional information is included in the original extended abstract.

Ivanov, B. A.↗
NASA NTRSDate not supplied

NASA HECC Geometry and Performance Review Part 3: A Numerical and Experimental Investigation of Tip Clearance Effects on the Vaneless Diffuser Configuration

Tip clearance effects in centrifugal compressors have been extensively investigated to understand the losses associated with the flow in the impeller tip clearance gap between the rotating blades and the stationary shroud. In Part 3 of this multipart investigation, experimental data and validated numerical simulations from Parts 1 and 2 were used to analyze the effects of the size of the tip clearance gap on the High Efficiency Centrifugal Compressor performance and aerodynamics. Four tip clearance gaps ranging from 0.012-in to 0.030-in (2.0% to 4.9% of the impeller exit blade height) were considered at both design and off-design operating conditions. The total pressure ratio and efficiency of the stage are found to decay linearly with increasing tip gaps. The sensitivity of the impeller performance to the tip gap was found to vary with the rotational speed of the compressor. Spanwise surveys of flow angle, total pressure, and total temperature collected at the impeller exit at design speed are used to validate numerical simulations at each experimental tip gap condition. Numerical simulations show that increased turbulence kinetic energy near the shroud at larger tip gaps leads to spanwise mixing of high entropy fluid near the impeller trailing edge which decreases the useful work input by the impeller. The data presented have been made available to the public in the HECC Data Archive located at https://storage.googleapis.com/hecc-data/NASA-HECC-Data-Archive.zip.

centrifugal compressor↗
NASA NTRSDate not supplied

NASA HECC Geometry and Performance Review Part 3: A Numerical and Experimental Investigation of Tip Clearance Effects on the Vaneless Diffuser Configuration

Tip clearance effects in centrifugal compressors have been extensively investigated to understand the losses associated with the flow in the impeller tip clearance gap between the rotating blades and the stationary shroud. In Part 3 of this multipart investigation, experimental data and validated numerical simulations from Parts 1 and 2 were used to analyze the effects of the size of the tip clearance gap on the High Efficiency Centrifugal Compressor performance and aerodynamics. Four tip clearance gaps ranging from 0.012-in to 0.030-in (2.0% to 4.9% of the impeller exit blade height) were considered at both design and off-design operating conditions. The total pressure ratio and efficiency of the stage are found to decay linearly with increasing tip gaps. The sensitivity of the impeller performance to the tip gap was found to vary with the rotational speed of the compressor. Spanwise surveys of flow angle, total pressure, and total temperature collected at the impeller exit at design speed are used to validate numerical simulations at each experimental tip gap condition. Numerical simulations show that increased turbulence kinetic energy near the shroud at larger tip gaps leads to spanwise mixing of high entropy fluid near the impeller trailing edge which decreases the useful work input by the impeller. The data presented have been made available to the public in the HECC Data Archive located at https://storage.googleapis.com/hecc-data/NASA-HECC-Data-Archive.zip.

centrifugal compressor↗