Search NASA⌕ Search

SEARCH · Search NASA

Results for “Entropy stability”

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 73 records · Page 4

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.↗

Extrema principles of entrophy production and energy dissipation in fluid mechanics

A survey is presented of several extrema principles of energy dissipation as applied to problems in fluid mechanics. An exact equation is derived for the dissipation function of a homogeneous, isotropic, Newtonian fluid, with terms associated with irreversible compression or expansion, wave radiation, and the square of the vorticity. By using entropy extrema principles, simple flows such as the incompressible channel flow and the cylindrical vortex are identified as minimal dissipative distributions. The principal notions of stability of parallel shear flows appears to be associated with a maximum dissipation condition. These different conditions are consistent with Prigogine's classification of thermodynamic states into categories of equilibrium, linear nonequilibrium, and nonlinear nonequilibrium thermodynamics; vortices and acoustic waves appear as examples of dissipative structures. The measurements of a typical periodic shear flow, the rectangular wall jet, show that direct measurements of the dissipative terms are possible.

Horne, W. Clifton↗

Extrema principles of entropy production and energy dissipation in fluid mechanics

A survey is presented of several extrema principles of energy dissipation as applied to problems in fluid mechanics. An exact equation is derived for the dissipation function of a homogeneous, isotropic, Newtonian fluid, with terms associated with irreversible compression or expansion, wave radiation, and the square of the vorticity. By using entropy extrema principles, simple flows such as the incompressible channel flow and the cylindrical vortex are identified as minimal dissipative distributions. The principal notions of stability of parallel shear flows appear to be associated with a maximum dissipation condition. These different conditions are consistent with Prigogine's classification of thermodynamic states into categories of equilibrium, linear nonequilibrium, and nonlinear nonequilibrium thermodynamics; vortices and acoustic waves appear as examples of dissipative structures. The measurements of a typical periodic shear flow, the rectangular wall jet, show that direct measurements of the dissipative terms are possible.

Horne, W. Clifton↗

A Dynamical Model for The Jupiter Thermocline

A new conceptual model for the observed banded wind structure of the giant outer planets is proposed, assuming a smoothly distributed potential vorticity (PV) atop a layer of deeply seated, latitudinally variable stratification. With a vanishing horizontal entropy contrast presumably enforced at the bottom of the flow layer by the underlying convection, the thermal-wind balance of observed cloud-top motions implies a mapping of constant potential temperature surfaces, nearly vertical at upper tropospheric levels, down to a deeper, flatter, but variably sloped 'thermocline.' For a generally colder-poleward distribution of isentropes consistent with the strong equatorial jets of both Jupiter and Saturn, this kind of mapping would also imply a poleward decreasing static stability. The proposed temperature-stability distribution is just the reverse of the Earth's troposphere, where the strongest latitudinal potential temperature gradients are at the bottom instead of the top, with static stability generally increasing toward the pole. The warmer-stabler correlation for Jupiter would be consistent, however, with a nearly monotonic distribution of PV from low to high latitudes, over planetary scales for which the planetary vorticity dominates the relative vorticity of the jets. In this way efficient PV mixing for the isentropically bounded thermocline can account for the dynamical maintenance of the cyclonic flanks of the equatorial jet. Over latitudinal intervals comparable to the internal deformation scale, however, the gradient of the absolute vorticity is dominated by the flow curvature, with the local latitudinal gradient of the static stability proportional to the local potential temperature gradient and therefore the local variation in geostrophic velocity. The resulting correlation of velocities and vorticity gradients implies an alternation of the flow over the internal deformation radius set by the vertical entropy contrast.

Allison, Michael↗

Reversibility and stability of information processing systems

Classical and quantum models of dynamically reversible computers are considered. Instabilities in the evolution of the classical 'billiard ball computer' are analyzed and shown to result in a one-bit increase of entropy per step of computation. 'Quantum spin computers', on the other hand, are not only microscopically, but also operationally reversible. Readoff of the output of quantum computation is shown not to interfere with this reversibility. Dissipation, while avoidable in principle, can be used in practice along with redundancy to prevent errors.

Zurek, W. H.↗

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↗

On the Dynamics of TVD Schemes

The dynamics of a class of TVD schemes for model hyperbolic and parabolic equations is studied numerically using a highly parallel supercomputer (CM-5). The objective is to utilize the highly parallel property of the CM-5 to reveal the reliable time step and entropy parameter ranges, and the degree of compressible flux limiters to avoid slow convergence and the production of nonphysical numerical solutions. We choose to study the nonlinear stability property of TVD schemes numerically since it is otherwise not amenable analytically.

Yee, H. C.↗

Size and Structure of the Sequence Space of Repeat Proteins

The coding space of protein sequences is shaped by evolutionary constraints set by requirements of function and stability. We show that the coding space of a given protein family— the total number of sequences in that family—can be estimated using models of maximum entropy trained on multiple sequence alignments of naturally occurring amino acid sequences. We analyzed and calculated the size of three abundant repeat proteins families, whose members are large proteins made of many repetitions of conserved portions of *30 amino acids. While amino acid conservation at each position of the alignment explains most of the reduction of diversity relative to completely random sequences, we found that correlations between amino acid usage at different positions significantly impact that diversity. We quantified the impact of different types of correlations, functional and evolutionary, on sequence diversity. Analysis of the detailed structure of the coding space of the families revealed a rugged landscape, with many local energy minima of varying sizes with a hierarchical structure, reminiscent of frustrated energy landscapes of spin glass in physics. This clustered structure indicates a multiplicity of subtypes within each family and suggests new strategies for protein design.

Jacopo Marchi↗

A Brief Review of Chaotropic vs. Kosmotropic Salts on the Degree of Saturation of Membrane Fatty Acids: Implications for Preservation of Environmental Biomarkers

Water activity is well known to have a significant impact on habitability, due to the presence of salts. However, salt composition is also important for habitability, as it can impact microbial activity and biomolecular interactions since certain salts are known to disrupt hydrogen bonding. Chaotropic agents disrupt hydrogen bonding in aqueous solutions leading to increased entropy (chaos) whereas kosmotropic agents are compounds that promote hydrogen bonding (order). The Hofmeister series (Fig. 1) shows the ordering of ions in terms of their ability to stabilize and induce precipitation or aggregation of proteins.

M B Wilhelm↗

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.↗

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.↗

Nonlinear Nonmodal Analysis of Hypersonic Flow over Blunt Cones

The linear amplification of modal disturbances that lead to boundary-layer transition in two-dimensional/axisymmetric hypersonic configurations is strongly reduced by the presence of a blunt nosetip, and the mechanisms underlying the observed onset of transition over the cone frustum are currently unknown. Linear nonmodal analysis has shown that both planar and oblique traveling disturbances that peak within the entropy layer experience appreciable energy amplification for moderate to large nosetip bluntness. The present study extends the previous linear analysis by including the nonlinear effects. Specifically, the perturbation form of the 2D, harmonic Navier-Stokes equations (HNSE) are solved with a fully implicit formulation and the Newton-Raphson method. The increased number of degrees of freedom for the nonlinear system presents difficulties for solution strategies based on direct solution of the linearized system. Such difficulties are overcome by using the GMRES iterative method with a preconditioner corresponding to a simplified Jacobian without the cross derivative terms. The HNSE solver is verified by comparing with nonlinear parabolized stability equation (NPSE) results for the nonlinear evolution of planar waves in an incompressible Blasius boundary layer and in a Mach 6 flow over a blunt cone. Finally, nonlinear nonmodal results are presented for planar traveling disturbances over the blunt cone. The nonmodal analysis demonstrates that entropy-layer disturbances generated close to the nose tip can seed the amplification of higher frequency Mack’s second-mode instabilities further downstream.

boundary layer transition↗

Nonlinear Nonmodal Analysis of Hypersonic Flow over Blunt Cones

The linear amplification of modal disturbances that lead to boundary-layer transition in two-dimensional/axisymmetric hypersonic configurations is strongly reduced by the presence of a blunt nosetip, and the mechanisms underlying the observed onset of transition over the cone frustum are currently unknown. Linear nonmodal analysis has shown that both planar and oblique traveling disturbances that peak within the entropy layer experience appreciable energy amplification for moderate to large nosetip bluntness. The present study extends the previous linear analysis by including the nonlinear effects. Specifically, the perturbation form of the 2D, harmonic Navier-Stokes equations (HNSE) are solved with a fully implicit formulation and the Newton-Raphson method. The increased number of degrees of freedom for the nonlinear system presents difficulties for solution strategies based on direct solution of the linearized system. Such difficulties are overcome by using the GMRES iterative method with a preconditioner corresponding to a simplified Jacobian without the cross derivative terms. The HNSE solver is verified by comparing with nonlinear parabolized stability equation (NPSE) results for the nonlinear evolution of planar waves in an incompressible Blasius boundary layer and in a Mach 6 flow over a blunt cone. Finally, nonlinear nonmodal results are presented for planar traveling disturbances over the blunt cone. The nonmodal analysis demonstrates that entropy-layer disturbances generated close to the nose tip can seed the amplification of higher frequency Mack’s second-mode instabilities further downstream.

boundary layer transition↗

Convection without eddy viscosity: An attempt to model the interiors of giant planets

In the theory of hydrostatic quasi-geostrophic flow in the Earth's atmosphere the principal results do not depend on the eddy viscosity. This contrasts with published theories of convection in deep rotating fluid spheres, where the wavelength of the fastest growing disturbance varies as E sup 1/3, where E, the Ekman number, is proportional to the eddy viscosity. A new theory of quasi-columnar motions in stably stratified fluid spheres attempts to capture the luck of the meteorologists. The theory allows one to investigate the stability of barotropic and baroclinic zonal flows that extend into the planetary interior. It is hypothesized that the internal heat Jupiter and Saturn comes out not radially but on sloping surfaces defined by the internal entropy distribution. To test the hypothesis one searches for basic states in which the wavelength of the fastest-growing disturbance remains finite as E tends to zero, and is which the heat flux vector is radially outward and poleward.

Ingersoll, A. P.↗

The convergence rate of approximate solutions for nonlinear scalar conservation laws

The convergence rate is discussed of approximate solutions for the nonlinear scalar conservation law. The linear convergence theory is extended into a weak regime. The extension is based on the usual two ingredients of stability and consistency. On the one hand, the counterexamples show that one must strengthen the linearized L(sup 2)-stability requirement. It is assumed that the approximate solutions are Lip(sup +)-stable in the sense that they satisfy a one-sided Lipschitz condition, in agreement with Oleinik's E-condition for the entropy solution. On the other hand, the lack of smoothness requires to weaken the consistency requirement, which is measured in the Lip'-(semi)norm. It is proved for Lip(sup +)-stable approximate solutions, that their Lip'convergence rate to the entropy solution is of the same order as their Lip'-consistency. The Lip'-convergence rate is then converted into stronger L(sup p) convergence rate estimates.

Nessyahu, Haim↗

Towards Accurate Predictions of Martensitic Transition Temperatures for Shape Memory Alloys from Ab Initio Simulations

Experimentally, NiTi undergoes a single martensitic phase transition around 341 K from the lowtemperature (T) monoclinic B19’ phase (P21/m) to the high-temperature cubic B2 phase (Pm3m). Theoretically, an orthorhombic B33 (Cmcm) has also been proposed as the T=0 ground state structure, although this phase has never been observed. Accurate predictions of martensitic transition temperatures (MTT) have remained elusive in part due to several well-known theoretical complexities of these systems including low temperature instabilities of the B2 phase. Recently, we proposed a rigorous thermodynamic integration approach based on ab initio simulations to resolve many of these difficulties [1,2]. However, an unsatisfying overprediction of the MTT relative to experiment (by ~100 K) means a fully quantitative theory is still lacking. In this work, we report several new developments to our method that bring first principles theory and experiment much closer into agreement. Our calculations indicate that phonon free energies at low temperature stabilizes B19’ over B33, rationalizing B19’ as the ground state down to T=0. We also find that accurate computations of the electronic free energy, i.e. the change in energy and the appearance of electronic configurational entropy due to finite temperature, is crucial to obtain accurate MTT. Incorporating these corrections results in an MTT prediction of 365 K for binary NiTi, which is in very close agreement with experiment. Our theoretical approach is expected to be a broadly applicable and predictive theory for MTT of SMAs.

Zhigang Wu↗

Theoretical prediction of phase relationships in planetary mantles

Thermodynamic and phase equilibrium data are used to generate an internally consistent set of enthalpies and entropies for important components of the CaO-MgO-Al2O3-SiO2 (CMAS) system in silica-undersaturated compositions. The addition of Na and Fe(2+) to the CMAS system produces shifts in the plagioclase, spinel and garnet stability fields. While the Morgan and Anders (1979) model Martian composition has stability fields of plagioclase and garnet lherzolite, and a small spinel lherzolite field at temperatures below 900 C, the Martian mantle composition of McGetchin and Smyth (1978) would not contain orthopyroxene, and a low pressure assemblage of plagioclase-spinel wehrlite would be replaced by garnet-spinel wehrlite at higher pressure. In both cases, the Fe-Mg ratio would be substantially greater than that found in primitive terrestrial basalts.

Wood, B. J.↗

Sloshing of the Magnetized Cool Gas in the Cores of Galaxy Clusters

X-ray observations of many clusters of galaxies reveal the presence of edges in surface brightness and temperature, known as "cold fronts". In relaxed clusters with cool cores, these edges have been interpreted as evidence for the "sloshing" of the core gas in the cluster's gravitational potential. The smoothness of these edges has been interpreted as evidence for the stabilizing effect of magnetic fields "draped" around the front surfaces. To check this hypothesis, we perform high-resolution magnetohydrodynamics simulations of magnetized gas sloshing in galaxy clusters initiated by encounters with subclusters. We go beyond previous works on the simulation of cold fronts in a magnetized intracluster medium by simulating their formation in realistic, idealized mergers with high resolution ((Delta)x approx. 2 kpc). Our simulations sample a parameter space of plausible initial magnetic field strengths and field configurations. In the simulations, we observe strong velocity shears associated with the cold fronts amplifying the magnetic field along the cold front surfaces, increasing the magnetic field strength in these layers by up to an order of magnitude, and boosting the magnetic pressure up to near-equipartition with thermal pressure in some cases. In these layers, the magnetic field becomes strong enough to stabilize the cold fronts against Kelvin-Helmholtz instabilities, resulting in sharp, smooth fronts as those seen in observations of real clusters. These magnetic fields also result in strong suppression of mixing of high and low-entropy gas in the cluster, seen in our simulations of mergers in the absence of a magnetic field. As a result, the heating of the core due to sloshing is very modest and is unable to stave off a cooling catastrophe.

ZuHone, J. A.↗