Search NASA⌕ Search

SEARCH · Search NASA

Results for “Entropy Stable”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS1976

Ice crystal growth in water vapor at high saturation

A simple technique is presented for estimating the energy of formation of monolayer icelike clusters at ice-vapor interfaces. Under the assumptions that the ice surfaces are smooth and sparsely covered with monomers, dimers, etc., in near equilibrium with the vapor, and that the bond energies and configurational entropy dominate the energy of formation, it is found that the basal surfaces prefer triangular embryos with an orientation which reverses from layer to layer, whereas the most stable clusters on the prism surfaces are rectangular in configuration. The preferred prism clusters are determined to have a significantly lower critical energy of formation than the basal clusters due to differences in both corner free energy and configurational entropy. This phenomenon provides a mechanism for strongly anisotropic crystal growth at high saturations.

Bartley, D. L.↗
NASA NTRS1988

Second order accurate finite difference approximations for the transonic small disturbance equation and the full potential equation

New shock-capturing finite difference approximations for solving two scalar conservation law nonlinear partial differential equations describing inviscid, isentropic, compressible flows of aerodynamics at transonic speeds are presented. A global linear stability theorem is applied to these schemes in order to derive a necessary and sufficient condition for the finite element method. A technique is proposed to render the described approximations total variation-stable by applying the flux limiters to the nonlinear terms of the difference equation dimension by dimension. An entropy theorem applying to the approximations is proved, and an implicit, forward Euler-type time discretization of the approximation is presented. Results of some numerical experiments using the approximations are reported.

Mostrel, M. M.↗
NASA NTRS2023

Chemical Thermodynamics and the Mathematical Integration of Reaction Kinetics

Key advances in the development of numerical methods for non-reacting compressible flows have been enabled by translating physical requirements into concrete numerical guidelines, such as the satisfaction of entropy inequalities for shock-capturing techniques [Lax, Contributions to Nonlinear Functional Analysis (1971) 603-634]. In the present work, we present nonlinear numerical analysis tools that draw from Chemical Thermodynamics , the branch of Nonequilibrium Thermodynamics that deals with chemical reactions. Through Gibbs formalism, chemical thermodynamics provides a well-known theoretical expression for the chemical equilibrium constant of a reaction in terms of reduced chemical potentials. A less-known, yet extremely valuable result, due to [Krambeck, Arch. Ration. Mech. Anal. , 38 (1970) 317], states that when this expression is implemented, mass-action kinetic models are consistent with the dynamical prescriptions of the 2nd law of thermodynamics. For fixed-temperature ordinary differential equations modeling constant-volume reacting gas mixtures, this leads to a decreasing Helmholtz free energy. If the temperature is allowed to vary in accordance with conservation of energy (1st law), this leads to the statement of increasing entropy. These nonlinear prescriptions can, and should be, used to further develop temporal integration techniques for reaction kinetics. We demonstrate that Krambeck's result holds even when the equilibrium constants are approximated from data. We prove this result by constructing the implicit free energy and the implicit entropy inherent to a given approximation. This is first done for a 5-species, 17-reaction model problem for air. With this structure established, elements of discrete entropy-stability theory [Tadmor, Acta Numer. , 12 (2003) 451] are leveraged to examine the consistency of time-integration schemes with these prescriptions. Using chemical potentials, one can compute the respective contributions of the kinetics model and of the temporal scheme to free energy/entropy variations. We introduce a nonlinear-stable version of the Discontinuous-Galerkin (DG) scheme in time which shows robustness improvements over the standard linearly-stable version. Most notably, the maximum timestep that can be resolved with the nonlinearly-stable variant tends to grow with polynomial order, in contrast to the linearly-stable variant. We generalize our constructions to arbitrary systems of reversible chemical reactions, ultimately showing that the compressible reacting Euler system admits the opposite of the implicitly constructed thermodynamic entropy as a mathematical entropy . This lays important theoretical foundations towards robust scheme development [Harten, J. Comput. Phys. 49 (1983) 151-164].

STMD↗
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 NTRS1985

Stable boundary conditions and difference schemes for Navier-Stokes equations

The Navier-Stokes equations can be viewed as an incompletely elliptic perturbation of the Euler equations. By using the entropy function for the Euler equations as a measure of energy for the Navier-Stokes equations, it was possible to obtain nonlinear energy estimates for the mixed initial boundary value problem. These estimates are used to derive boundary conditions which guarantee L2 boundedness even when the Reynolds number tends to infinity. Finally, a new difference scheme for modelling the Navier-Stokes equations in multidimensions for which it is possible to obtain discrete energy estimates exactly analogous to those we obtained for the differential equation was proposed.

Dutt, P.↗
NASA NTRS1991

A new finite element formulation for computational fluid dynamics. X - The compressible Euler and Navier-Stokes equations

A space-time element method is presented for solving the compressible Euler and Navier-Stokes equations. The proposed formulation includes the variational equation, predictor multi-corrector algorithms and boundary conditions. The variational equation is based on the time-discontinuous Galerkin method, in which the physical entropy variables are employed. A least-squares operator and a discontinuity-capturing operator are added, resulting in a high-order accurate and unconditionally stable method. Implicit/explicit predictor multi-corrector algorithms, applicable to steady as well as unsteady problems, are presented; techniques are developed to enhance their efficiency. Implementation of boundary conditions is addressed; in particular, a technique is introduced to satisfy nonlinear essential boundary conditions, and a consistent method is presented to calculate boundary fluxes. Numerical results are presented to demonstrate the performance of the method.

Shakib, Farzin↗
NASA NTRS1989

Conservation equations and physical models for hypersonic air flows over the aeroassist flight experiment vehicle

The code development and application program for the Langley Aerothermodynamic Upwind Relaxation Algorithm (LAURA), with emphasis directed toward support of the Aeroassist Flight Experiment (AFE) in the near term and Aeroassisted Space Transfer Vehicle (ASTV) design in the long term is reviewed. LAURA is an upwind-biased, point-implicit relaxation algorithm for obtaining the numerical solution to the governing equations for 3-D, viscous, hypersonic flows in chemical and thermal nonequilibrium. The algorithm is derived using a finite volume formulation in which the inviscid components of flux across cell walls are described with Roe's averaging and Harten's entropy fix with second-order corrections based on Yee's Symmetric Total Variation Diminishing scheme. Because of the point-implicit relaxation strategy, the algorithm remains stable at large Courant numbers without the necessity of solving large, block tri-diagonal systems. A single relaxation step depends only on information from nearest neighbors. Predictions for pressure distributions, surface heating, and aerodynamic coefficients compare well with experimental data for Mach 10 flow over an AFE wind tunnel model. Predictions for the hypersonic flow of air in chemical and thermal nonequilibrium over the full scale AFE configuration obtained on a multi-domain grid are discussed.

Gnoffo, Peter A.↗
NASA NTRS1992

Discrete modes and continuous spectra in supersonic boundary layers

The disturbance field induced due to a harmonic point source consists of discrete eigenmodes and a continuous spectrum; these are studied by using generalized Fourier transform techniques. For a supersonic boundary layer, there exist seven branches of the continuous spectrum in the complex wavenumber space, four of which (two acoustic waves, one vorticity wave and one entropy wave) contribute to the flowfield downstream of the source. The discrete eigenmodes spring off from these branches at some critical Reynolds numbers. The results for Mach 2 and 4.5 boundary layers show that the receptivity coefficients for the stable discrete modes are much larger than that for the unstable mode. Therefore, the flow very near the source is dominated by the continuous spectrum and the stable discrete modes. However, the unstable mode takes over sufficiently far away from the source. It is shown that it is only necessary to consider the first few discrete modes to construct the solution. Calculations also show that, in a supersonic boundary layer, upstream influence from a localized disturbance is minimal.

Balakumar, P.↗
NASA NTRS1979

A reconsideration of nucleation phenomena in light of recent findings concerning the properties of small clusters, and a brief review of some other particle growth processes

The paper examines mechanisms of nucleation and growth by condensation and coagulation in the light of recent research on properties of small clusters. Homogeneous, hetero-molecular, and heterogeneous nucleation is analyzed, and expressions for the rate of formation of a stable condensed phase and evaluation of the free energy of formation of charged droplets are given. Application of high-pressure mass spectrometry which makes possible a direct determination of intensity spectra for cluster distributions, measurement of the thermodynamic properties of individual ion clusters and determination of cluster entropy and bond energy is discussed. Finally, coagulation of the condensed phase is considered, noting that concentration and mean particle size vary during coagulation, but the shape distribution is time independent, leading to the concept of a self-preserving aerosol size distribution.

Castleman, A. W., Jr.↗
NASA NTRS1986

Molecular microenvironments: Solvent interactions with nucleic acid bases and ions

The possibility of reconstructing plausible sequences of events in prebiotic molecular evolution is limited by the lack of fossil remains. However, with hindsight, one goal of molecular evolution was obvious: the development of molecular systems that became constituents of living systems. By understanding the interactions among molecules that are likely to have been present in the prebiotic environment, and that could have served as components in protobiotic molecular systems, plausible evolutionary sequences can be suggested. When stable aggregations of molecules form, a net decrease in free energy is observed in the system. Such changes occur when solvent molecules interact among themselves, as well as when they interact with organic species. A significant decrease in free energy, in systems of solvent and organic molecules, is due to entropy changes in the solvent. Entropy-driven interactioins played a major role in the organization of prebiotic systems, and understanding the energetics of them is essential to understanding molecular evolution.

Macelroy, R. D.↗
NASA NTRSDate not supplied

Detecting Risk and Anomalies in Airplane Dynamics Through Entropic Analysis of Time Series Data

Despite recent efforts to move away from traditional threshold exceedance detection methods for aircraft state monitoring, modern aircraft still rely on safety thresholds to communicate to pilots the identification of an anomaly in the aircraft when a threshold is surpassed. Current anomaly detection methods mainly depend on uninterpretable machine learning models to learn complex patterns and relationships contained in the time series data of aircraft. Although these methods are capable of identifying known anomalies, their deficiency in interpretability presents a challenge when translating them to different aircraft. To overcome this deficiency, entropic analysis of aircraft dynamics seeks to characterize the complexity, or lack thereof, of the aircraft dynamics prior to the development of a risk scenario. This complexity characterization provides a more straightforward summary of state changes in the dynamics of flight variables. To build a foundation for entropic analysis, we analyzed the complexity of unstable approaches, an anomalous event present in many of today’s aviation accidents. The analysis revealed a statistically significant difference in the complexity distribution of flight variables under a stable approach versus an unstable approach. These differences in complexity were especially notable minutes before an approach was identified as unstable. Moreover, the multiscale entropic analysis revealed the presence of signal complexity at multiple time scales across multiple time windows before landing. By capturing state changes and corrections in the aircraft dynamics using entropy, advanced, yet still interpretable, sensor systems based on entropic frameworks from this study can be constructed in the future using classical machine learning approaches.

Risk detection↗
NASA NTRS1984

Space and Time Variations and Turbopause Dynamical Structure

Time spectra of the critical frequency f sub o E sub s and blanketing frequency f sub b E sub s for the records of a spaced chain of ionosphere vertical sounding stations were studied by means of a maximum entropy method, while coherence spectra were analyzed with the Blackman and Tukey method. It was found out that the frequency parameters spectrum of E sub s within the range of 1-10 cycles/h are discrete with one or two maxima. The first, a low frequency maximum with T=40 min, is of larger amplitude, stable and exists permanently; the second one with T = 10 min is less stable, occurs irregularly and more often within f sub o E sub s. Further, the first maximum corresponds to a higher level of coherence which decreases as the distance to the stations increases. It is concluded that the low-frequency maximum is induced by cellular eddies, usually interpreted as wind shears with horizontal dimensions of not more than 300 km. Irregularity, small amplitude, and low coherence in the range of the second maximum are indicative of the fact that turbulence must be its only source.

Ovezgeldiyev, O. G.↗
NASA NTRS2001

Entropy Splitting for High Order Numerical Simulation of Vortex Sound at Low Mach Numbers

A method of minimizing numerical errors, and improving nonlinear stability and accuracy associated with low Mach number computational aeroacoustics (CAA) is proposed. The method consists of two levels. From the governing equation level, we condition the Euler equations in two steps. The first step is to split the inviscid flux derivatives into a conservative and a non-conservative portion that satisfies a so called generalized energy estimate. This involves the symmetrization of the Euler equations via a transformation of variables that are functions of the physical entropy. Owing to the large disparity of acoustic and stagnation quantities in low Mach number aeroacoustics, the second step is to reformulate the split Euler equations in perturbation form with the new unknowns as the small changes of the conservative variables with respect to their large stagnation values. From the numerical scheme level, a stable sixth-order central interior scheme with a third-order boundary schemes that satisfies the discrete analogue of the integration-by-parts procedure used in the continuous energy estimate (summation-by-parts property) is employed.

Mueller, B.↗
NASA NTRS1991

Atomic force microscope studies of fullerene films - Highly stable C60 fcc (311) free surfaces

Atomic force microscopy and X-ray diffractometry were used to study 1500 A-thick films of pure C60 grown by sublimation in ultrahigh vacuum onto a CaF2 (111) substrte. Topographs of the films did not reveal the expected close-packed structures, but they showed instead large regions that correspond to a face-centered cubic (311) surface and distortions of this surface. The open (311) structure may have a relatively low free energy because the low packing density contributes to a high entropy of the exposed surface.

Snyder, Eric J.↗
NASA NTRS1990

An upwind-biased, point-implicit relaxation algorithm for viscous, compressible perfect-gas flows

An upwind-biased, point-implicit relaxation algorithm for obtaining the numerical solution to the governing equations for three-dimensional, viscous, compressible, perfect-gas flows is described. The algorithm is derived using a finite-volume formulation in which the inviscid components of flux across cell walls are described with Roe's averaging and Harten's entropy fix with second-order corrections based on Yee's Symmetric Total Variation Diminishing scheme. Viscous terms are discretized using central differences. The relaxation strategy is well suited for computers employing either vector or parallel architectures. It is also well suited to the numerical solution of the governing equations on unstructured grids. Because of the point-implicit relaxation strategy, the algorithm remains stable at large Courant numbers without the necessity of solving large, block tri-diagonal systems. Convergence rates and grid refinement studies are conducted for Mach 5 flow through an inlet with a 10 deg compression ramp and Mach 14 flow over a 15 deg ramp. Predictions for pressure distributions, surface heating, and aerodynamics coefficients compare well with experiment data for Mach 10 flow over a blunt body.

Gnoffo, Peter A.↗
NASA NTRS2004

A Novel Liquid-Liquid Transition in Undercooled Ti-Zr-Ni Liquids

If crystallization can be avoided, liquids enter a metastable (undercooled) state below their equilibrium liquidus temperatures, T(sub l), finally 'freezing' into a glass below a characteristic temperature called the glass transition temperature, T(sub g). In rare cases, the undercooled liquid may undergo a liquid-liquid phase transition (liquid polymorphism) before entering the glassy state. This has been suggested from experimental studies of H2O and Si. Such phase transitions have been predicted in some stable liquids, ie. above T(sub l) at atmospheric pressure, for SiO2 and BeF2, but these have not been verified experimentally. They have been observed in liquids of P, Si and C, but only under high pressure. In this letter we present the first experimental evidence for a phase transition in a low viscosity metallic liquid that is driven by an approach to a constant entropy configuration state and correlated with a growing icosahedral order in the liquid. A maximum in the specific heat at constant pressure, similar to what is normally observed near T(sub g), is reported for undercooled liquids of quasicrystal-forming Ti-Zr-Ni alloys. A two-state excitation model that includes cooperativity by incorporating a temperature-dependent excitation energy, fits the specific heat data well, signaling a phase transition. An inflection in the liquid density with decreasing temperature instead of a discontinuity indicates that this is not a typical first order phase transition; it could be a weakly first order or higher order transition. While showing many similarities to a glass transition, this liquid-liquid phase transition occurs in a mobile liquid, making it novel.

Lee, G. W.↗
NASA NTRS1989

Upwind-biased, point-implicit relaxation strategies for hypersonic flowfield simulations on supercomputers

An upwind-biased, point-implicit relaxation algorithm for obtaining the numerical solution to the governing equations for three-dimensional, viscous, hypersonic flows in chemical and thermal nonequilibrium is described. The algorithm is derived using a finite-volume formulation in which the inviscid components of flux across cell walls are described with Roe's averaging and Harten's entropy fix with second-order corrections based on Yee's Symmetric Total Variation Diminishing scheme. The relaxation strategy is well suited for computers employing either vector or parallel architectures, and the relation between computer architecture and algorithm is emphasized. It is also well suited to the numerical solution of the governing equations on unstructured grids. Because of the point-implicit relaxation strategy, the algorithm remains stable at large Courant numbers without the necessity of solving large. block tri-diagonal systems. A single relaxation step depends only on information from nearest neighbors. Predictions for pressure distributions, surface heating, and aerodynamic coefficients compare well with experimental data for Mach 10 flow over a blunt body. Predictions for the hypersonic flow of air in chemical and thermal nonequilibrium (velocity = 8917 m/s, altitude = 78 km.) over the Aeroassist Flight Experiment (AFE) configuration obtained on a multi-domain grid are discussed.

Gnoffo, Peter A.↗
NASA NTRS1992

Maximum entropy approach to fuzzy control

For the same expert knowledge, if one uses different &- and V-operations in a fuzzy control methodology, one ends up with different control strategies. Each choice of these operations restricts the set of possible control strategies. Since a wrong choice can lead to a low quality control, it is reasonable to try to loose as few possibilities as possible. This idea is formalized and it is shown that it leads to the choice of min(a + b,1) for V and min(a,b) for &. This choice was tried on NASA Shuttle simulator; it leads to a maximally stable control.

Ramer, Arthur↗