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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 55 records · Page 3

All sources NASA NTRS DOE OSTI

NASA NTRS1992

A quasi-static magnetospheric convection model in two-dimensions

A self-consistent 2D model of subsonic plasma sheet convection is presented. Specifically, time sequences of static equilibrium solutions for the 2D magnetospheric magnetic field are constructed consistent with adiabatic convection. This model self-consistently includes a dipole field and a reasonable accounting for the effects of inner magnetospheric shielding. Starting from a relaxed magnetospheric equilibrium, the earthward convection of plasma sheet flux tubes results in the stretching of inner plasma sheet field lines, the development of a local minimum in the equatorial magnetic field Be in the near-earth plasma sheet, and an increasing lobe magnetic field. This evolution in time occurs independent of the specific magnetopause or far-tail boundary conditions, provided the plasma sheet is not at marginal interchange stability. This behavior results solely from the convection of flux tubes of increasing entropy into the near-earth plasma sheet. These results are discussed in the general context of earthward convection in earth's plasma sheet and in the specific context of magnetospheric substorms.

Erickson, G. M.↗

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 NTRS2020

Toward Transition Modeling in a Hypersonic Boundary Layer at Flight Conditions

An accurate physics-based transition prediction method integrated with computational fluid dynamics (CFD) solvers is pursued for hypersonic boundary layer flows over slender hypersonic vehicles at flight conditions. The geometry and flow conditions are selected to match relevant trajectory locations from the ascent phase of the HIFiRE-1 flight experiment, namely, a 7-degree half-angle cone with 2.5 mm nose radius, freestream Mach numbers in the range of 3.8 – 5.5 and freestream unit Reynolds numbers in the range of 3.3 × 10(exp 6) – 21.4 × 10(exp 6) m(exp -1). Earlier research had shown that the onset of transition during the HIFiRE-1 flight experiment correlated with an amplification factor of N ≈ 13.5 for the planar Mack modes. However, to incorporate the N-factor correlations into a CFD code, we investigate surrogate models for disturbance amplification that avoid the direct computation of stability characteristics. A commonly used approach for low-speed flows is based on an a priori database of stability characteristics for locally similar profiles. However, the results presented in this paper demonstrate that the application of this approach to hypersonic boundary layers over blunt spherical nose-tip cones leads to large, unacceptable errors in the predictions of amplification factors, mainly due to its failure in accounting for the effects of the entropy layer on the boundary-layer profiles along the length of the model. We propose and demonstrate an alternate approach that employs the stability computations for a canonical set of blunt cone configurations to train a physics-informed convolutional neural network model that is shown to provide substantially improved transition predictions for hypersonic flow configurations with entropy-layer effects. Furthermore, the excellent performance of the neural network model is also confirmed for cone configurations with nose radius and half-angle values that do not correspond to those used to build the database. Finally, the convolutional neural network model is shown to outperform the linear stability calculations for underresolved basic states.

Pedro Paredes↗

NASA NTRS2024

Intrastep, Stage-Value Predictors for Diagonally-Implicit Runge-Kutta Methods

To better identify the necessary attributes of good stage-value predictors (SVPs), numerous SVPs are designed for an existing: ESDIRK4(3)7L[2]SA [26] and a new: ESDIRK4(3)8L[2]SA scheme.1 Both are stifflyaccurate, stage-order two, explicit, singly-diagonally implicit Runge–Kutta (ESDIRK) schemes. Tradeoffs are studied in the parameter spaces enforcing the constraints on accuracy, linear stability, nonlinear stability and coefficient size to determine which objectives correlate with effective predictors. The SVPs are tested on three challenging external aerodynamics problems [107 − 108 degrees of freedom (DoFs)], each with a different level of stiffness. The problems include two 3D airfoils simulations and one canonical turbulence simulation. All simulations use the compressible Navier-Stokes equations (CNSE). An entropy stable spectral collocation formulation is used for discretizing the spatial terms in the equations. Simulations are performed at a wide variety of temporal error tolerances. Problems that are sufficiently stiff (e.g., lax temporal error tolerances) benefit from SVPs designed with second-order accuracy and stability properties: A-stability, and L-stability, rather than high accuracy constraints. Simulations with modest stiffness (e.g., strict error tolerances) are better suited for SVPs designed using high accuracy constraints. Designing SVPs with enhanced stability properties is tedious but worthwhile. Simulation times are reduced with optimal SVPs by as much as 100% on some stages, with combined stepwise improvements of between 50 − 100% for both methods. A comparative study is performed with the two aforementioned methods as well as four other ESDIRKs. The newly designed ESDIRK4(3)8L[2]SA with γ ≈ 1/10, proves to be the most efficient of the six tested ESDIRK schemes simulating the CNSE.

Diagonally-Implicit Runge-Kutta↗

NASA NTRS2002

Entropy Splitting for High Order Numerical Simulation of Compressible Turbulence

A stable high-order numerical scheme for direct numerical simulation (DNS) of shock-free compressible turbulence is presented. The method is applicable to general geometries. It contains no upwinding, artificial dissipation, or filtering. Instead the method relies on the stabilizing mechanisms of an appropriate conditioning of the governing equations and the use of compatible spatial difference operators for the interior points (interior scheme) as well as the boundary points (boundary scheme). An entropy-splitting approach splits the inviscid flux derivatives into conservative and nonconservative portions. The spatial difference operators satisfy a summation-by-parts condition, leading to a stable scheme (combined interior and boundary schemes) for the initial boundary value problem using a generalized energy estimate. A Laplacian formulation of the viscous and heat conduction terms on the right hand side of the Navier–Stokes equations is used to ensure that any tendency to odd–even decoupling associated with central schemes can be countered by the fluid viscosity. The resulting methods are able to minimize the spurious high-frequency oscillations associated with pure central schemes, especially for long time integration applications such as DNS. For validation purposes, the methods are tested in a DNS of compressible turbulent plane channel flow at low values of friction Mach number, where reference turbulence data bases exist. It is demonstrated that the methods are robust in terms of grid resolution, and in good agreement with published channel data. Accurate turbulence statistics can be obtained with moderate grid sizes. Stability limits on the range of the splitting parameter are determined from numerical tests.

N.D. Sandham↗

NASA NTRS2021

Hypersonic Boundary-Layer Transition on Blunted Cones at Angle of Attack

Experimental studies of cones at several high-speed facilities have demonstrated that, for small nosetip bluntness, transition onset over a circular cone moves upstream along the leeward side and downstream along the windward side, but this trend may be reversed at large bluntness values, where transition onset moves downstream along the leeward side and upstream along the windward side. A theoretical and numerical investigation is performed to characterize the effects of nose bluntness on disturbance amplification over the circular cone for several angles of attack, with the goal of understanding the potential physical mechanisms behind the experimental observations. The three-dimensional laminar basic states over a 1.5 m long, 7-degree half-angle cone with 9.525 mm nosetip radius are computed for selected angles of attack values and freestream conditions that are selected to match the Mach 10 experiments conducted within the Hypervelocity Wind Tunnel 9 at the Arnold Engineering Development Complex (AEDC). The solutions at a freestream unit Reynolds number of 17.1 million per meter are used to perform detailed instability analyses for angles of attack equal to 0, 1, 3, and 5 degrees. Results indicate that the linear amplification of stationary crossflow waves along inflection lines may begin to influence transition along the acreage of the cone for angles of attack equal to or larger than 5 degrees. The measured trend in transition front with respect to increasing angle of attack is found to be consistent with the predicted increase in the amplification factors for Mack mode disturbances along the streamline trajectories. The increase in Mack mode amplification along the windward ray for higher angles of attack is shown to be the result of a progressively earlier entropy-layer swallowing. Computations also indicate that the transition amplification factor along the windward ray is not constant and increases with the angle of attack and that the transition #-factors along the leeward ray are rather small. The nonmodal analysis for zero degrees angle of attack shows that entropy-layer disturbances with appreciably strong energy growth can coexist with Mack mode instabilities at the measured transition location.

Boundary Layer Transition↗

NASA NTRS1994

Stability of fireballs and gamma-ray bursts

Fireballs are an essential part of any cosmological gamma-ray burst. We derive a stability criterion for fireballs and show that fireballs are Rayleigh-Taylor unstable in any region in which the entropy decreases outward. The instability begins to operate when the fireball becomes matter dominated. Among the possible implications of the instability are the following: (1) Conversion of a fraction of the radiation energy to a convective energy expressed in the motion of bubbles relative to each other. (2) Penetration of fast bubbles through slower ones and creation of high-gamma regimes which are essential for efficient conversion of the energy to gamma-rays. (3) Formation of rapid time variation (of the scale of the bubbles) in the observed gamma-rays.

Waxman, Eli↗

NASA NTRS2002

Interchange Method in Compressible Magnetized Couette Flow: Magnetorotational and Magnetoconvective Instabilities

We obtain the general forms of the axisymmetric stability criteria in a magnetized compressible Couette flow using an energy variational principle, the so-called interchange or Chandrasekhar s met hod, which we applied successfully in the incompressible case. This formulation accounts for the simultaneous presence of gravity, rotation, a toroidal magnetic field, a weak axial magnetic field, entropy gradients, and density gradients in the initial equilibrium state. The power of the method lies in its simplicity which allows us to derive extremely compact and physically clear expressions for the relevant stability criteria despite the inclusion of so many physical effects. In the implementation of the method, all the applicable conservation laws are explicitly taken into account during the variations of a quantity with dimensions of energy which we call the free energy function. As in the incompressible case, the presence of an axial field invalidates the conservation laws of angular momentum and azimuthal magnetic flux and introduces instead isorotation and axial current conservation along field lines. Our results are therefore markedly different depending on whether an axial magnetic field is present, and generalize in two simple expressions all previously known, partial stability criteria for the appearance of magnetorotational instability. Furthermore, the coupling between magnetic tension and buoyancy and its influence to the dynamics of nonhomoentropic magnetized flows becomes quite clear from our results. In the limits of plane-parallel atmospheres and homoentropic flows, our formulation easily recovers the stability criteria for suppression of convective and Parker instabilities, as well as some related special cases studied over 40 years ago by Newcomb and Tserkovnikov via laborious variational techniques.

Christodoulou, Dimitris M.↗

NASA NTRS1993

Instabilities of two-dimensional inviscid compressible vortices

An investigation of the linear stability and subsequent nonlinear evolution and acoustic radiation of a planar inviscid compressible vortex is presented. The effects of the entropy gradient are investigated, and for the particular entropy profile chosen, the positive average entropy in the vortex core is destabilizing, while the opposite is true for the negative average entropy gradient. Finite-difference methods are used to study the linear initial value problem. These methods are found to be capable of accurately computing the frequencies and weak growth rates of the normal modes. When the initial condition consists of random perturbations, the long-time behavior is found to correspond to the most unstable normal mode in all cases. The numerical solution of the Euler equations is used to study the nonlinear development of an elliptic-mode perturbation.

Chan, W. M.↗

NASA NTRS2018

Nosetip Bluntness Effects on Transition at Hypersonic Speeds: Experimental and Numerical Analysis

The existing database of transition measurements in hypersonic ground facilities has established that the onset of boundary layer transition over a circular cone at zero angle of attack shifts downstream as the nosetip bluntness is increased with respect to a sharp cone. However, this trend is reversed at sufficiently large values of the nosetip Reynolds number, so that the transition onset location eventually moves upstream with a further increase in nosetip bluntness. This transition reversal phenomenon, which cannot be explained on the basis of linear stability theory, was the focus of a collaborative investigation under the NATO STO group AVT-240 on Hypersonic Boundary-Layer Transition Prediction. The current paper provides an overview of that effort, which included wind tunnel measurements in three different facilities and theoretical analysis related to modal and nonmodal amplification of boundary layer disturbances. Because neither first and second-mode waves nor entropy-layer instabilities are found to be substantially amplified to initiate transition at large bluntness values, transient (i.e., nonmodal) disturbance growth has been investigated as the potential basis for a physics based model for the transition reversal phenomenon. Results of the transient growth analysis indicate that stationary disturbances that are initiated within the nosetip or in the vicinity of the juncture between the nosetip and the frustum can undergo relatively significant nonmodal amplification and that the maximum energy gain increases nonlinearly with the nose radius of the cone. This finding does not provide a definitive link between transient growth and the onset of transition, but it is qualitatively consistent with the experimental observations that frustum transition during the reversal regime was highly sensitive to wall roughness, and furthermore, was dominated by disturbances that originated near the nosetip. Furthermore, the present analysis shows significant nonmodal growth of traveling disturbances that peak within the entropy layer and could also play a role in the transition reversal phenomenon.

Freestream Mach Number↗

NASA NTRS1994

Rapidly rotating polytropes in general relativity

We construct an extensive set of equilibrium sequences of rotating polytropes in general relativity. We determine a number of important physical parameters of such stars, including maximum mass and maximum spin rate. The stability of the configurations against quasi-radial perturbations is diagnosed. Two classes of evolutionary sequences of fixed rest mass and entropy are explored: normal sequences which behave very much like Newtonian evolutionary sequences, and supramassive sequences which exist solely because of relativistic effects. Dissipation leading to loss of angular momentum causes a star to evolve in a quasi-stationary fashion along an evolutionary sequence. Supramassive sequences evolve towards eventual catastrophic collapse to a black hole. Prior to collapse, the star must spin up as it loses angular momentum, an effect which may provide an observational precursor to gravitational collapse to a black hole.

Cook, Gregory B.↗

NASA NTRS2012

Crystallographic Stability of Metastable Phase Formed by Containerless Processing in REFeO3 (RE: Rare-Earth Element)

Undercooling a melt often facilitates a metastable phase to nucleate preferentially. Although the classical nucleation theory shows that the most critical factor for forming a metastable phase is the interface free energy, the crystallographic stability is also indispensable for the phase to be frozen at ambient temperature. In compound materials such as oxides, authors have suggested that the decisive factors for forming a critical nucleus are not only the free energy difference but also the difference of the entropy of fusion between stable and metastable phases. In the present study, using REFeO3 (RE: rare-earth element) as a model material, we investigate the formation of a metastable phase from undercooled melts with respect to the competitive nucleation and crystallographical stabilities of both phases.

Kuribayashi, Kazuhiko↗

NASA NTRS1982

On the instability of shock waves attached to wedges and cones

A numerical study is made of the stability of weak and strong shocks attached to wedges and cones. The stability is shown to depend on the boundary conditions imposed, in good agreement with the experimental observations and the principle of minimum entropy. The behavior of the solutions is explained by quasi-static considerations.

Salas, M. D.↗

NASA NTRS2022

Linear Disturbance Amplification Over Blunted Flat Plates in High-Speed Flows

Modal and nonmodal instability characteristics of cylindrically blunted flat plates with varying leading edge radii are described for Mach 4 and Mach 6 freestream conditions. The selection of leading edge radii and freestream parameters is informed by experimental conditions. The investigation of this 2D problem provides a slow entropy layer swallowing which allows for an isolated development of perturbations seeded upstream within different wall-normal regions of the flow. At both Mach numbers, a decrease in modal instability amplification was seen as the leading edge radius was increased. Nonmodal analysis reveals amplifying perturbations in the boundary layer as well as in the entropy layer. The medium bluntness regime exhibits the strongest amplification of nonmodal disturbances that is nonmonotonic in character. Small amplitude boundary forcing at the plate surface or volumetric forcing at various wall-normal heights was used to account for receptivity effects. While wall forcing effectively induced modal instabilities, only an actuation above the boundary layer captured disturbances that amplify within the entropy layer. The optimal nonmodal theory’s entropy-layer disturbance evolution exhibited outstanding agreement with controlled forcing, including receptivity effects. The evolution of entropy-layer disturbances from the optimal nonmodal theory showed excellent agreement with the results of receptivity to controlled forcing. Therefore, the nonmodal optimal growth analysis may provide a useful as well as efficient technique to identify the complete disturbance spectrum in blunt hypersonic configurations, where both modal and nonmodal disturbances can amplify in the boundary-layer and entropy-layer regions.

Boundary layer transition↗

NASA NTRS2019

Laminar-Turbulent Transition Upstream of the Entropy-Layer Swallowing Location in Hypersonic Boundary Layers

Numerical and experimental studies have demonstrated that modal growth of planar Mack modes is responsible for laminar-turbulent transition on sharp cones at hypersonic speeds. However, the physical mechanisms that lead to transition onset upstream of the entropy-layer swallowing location over sufficiently blunt geometries are not well understood as yet. Modal amplification is too weak or nonexistent to initiate transition at moderate-to-large bluntness values. Nonmodal analysis shows that, with increasing nose bluntness, both planar and oblique traveling disturbances that peak within the entropy layer experience appreciable energy amplification. However, because of the relatively weak signature of the nonmodal traveling disturbances within the boundary-layer region, the route to transition onset subsequent to the nonmodal growth remains unclear. Thus, nonlinear parabolized stability equations (NPSE) and direct numerical simulations (DNS) have been used to investigate the potential transition mechanisms over a 7-degree blunt cone that was tested in the AFRL Mach-6 high-Reynoldsnumber facility. Computations are performed to separately follow the nonlinear development of two classes of inflow disturbances, namely, a pair of oblique traveling waves with equal but opposite angles with respect to the mean flow direction and a planar traveling wave. Results in both cases show an excellent agreement between the NPSE and DNS predictions, establishing that the NPSE is an accurate and efficient technique for predicting the nonlinear development for these particular nonmodal traveling disturbances. Computations reveal that the oblique mode interactions lead to the generation of stationary streaks inside the boundary layer that, in turn, facilitate the growth of a subharmonic sinuous disturbance. For relatively modest amplitudes of the inflow disturbance, the oblique-mode breakdown can lead to transition at the measured location of transition onset during the experiment. On the other hand, the nonlinear development of a planar traveling wave leads to the formation of inclined structures just above the boundary-layer edge and these structures are strongly reminiscent of the transitional events observed during blunt cone experiments by using schlieren flow visualizations.

Entropy↗

NASA NTRS2021

Transition Analysis for Isolated Trips on BOLT-II Wind-Tunnel and Flight Configuration

The BOLT-II configuration designed for an upcoming flight experiment exemplifies the complexities of hypersonic boundary layer transition in the presence of multiple and potentially interacting instability mechanisms. The present paper is related to the NASA-led roughness experiment on the secondary side of the BOLT-II configuration. First, computations are performed to aid mean-flow based correlations for the transition measurements acquired during an experimental campaign in the Ludwieg Tube facility at the United States Air Force Academy. To that end, both laminar and turbulent flow computations are performed for selected conditions of the experiment with a 1:3 scale model of the slightly shorter BOLT configuration. The variation in well-known correlation parameters for transition due to isolated roughness elements with respect to element position, height, and flow Reynolds number is delineated. Furthermore, the wake instability theory is used to examine the amplification of unstable disturbances behind the centerline trip at the design conditions of the BOLT-II flight experiment and for its dynamically scaled version on the 1:3 scale model of the wind tunnel experiment. To our knowledge, this study represents the first analysis of the wake instabilities in the presence of a prominent entropy layer outside the boundary layer. At the nominal flight design condition of Re ∞ ) = 5.44 and Re ∞ = 2.5 x 10 6 /m, a diamond planform (“pizza-box”) trip with k/δ = 0.70 and planform-halfwidth-to-height ratio of b/k = 3.0 at X/L = 0.254 induces a strong distortion of the boundary-layer flow in the vicinity of the symmetry plane, such that the wake structure includes a centerline streak produced by the flow retardation behind the trip along with multiple additional, prominent vorticity structures representing the necklace vortex system. Stability results reveal that each of these streaks can support a strong amplification of wake instabilities within a frequency band of 95 kHz–255 kHz, with a peak N-factor of greater than 15 beyond X/L ≈ 0.42. The dynamically scaled wind tunnel configuration with Re ∞ = 20.25 x 10 6 /m is also found to support N > 15 at similar axial stations.

BOLT↗

NASA NTRS1974

The electromagnetic interchange mode in a partially ionized collisional plasma

A collisional electromagnetic dispersion relation is derived from two-fluid theory for the interchange mode coupled to the Alfven, acoustic, drift and entropy modes in a partially ionized plasma. The fundamental electromagnetic nature of the interchange model is noted; coupling to the intermediate Alfven mode is strongly stabilizing for finite k sub z. Both ion viscous and ion-neutral stabilization are included, and it was found that collisions destroy the ion finite Larmor radius cutoff at short perpendicular wavelengths.

Hudson, M. K.↗

NASA NTRS1975

The electromagnetic interchange mode in a partly-ionized collisional plasma

A collisional electromagnetic dispersion relation is derived from two-fluid theory for the interchange mode coupled to the Alfven, acoustic, drift, and entropy modes in a partially ionized plasma. The fundamental electromagnetic nature of the interchange mode is noted: coupling to the intermediate Alfven mode is strongly stabilizing for finite perturbations of the magnetic field. Both ion-viscous and ion-neutral stabilization are included; and it is found that collisions destroy the FLR (finite Larmor radius) cutoff at short perpendicular wavelengths.

Hudson, M. K.↗