SEARCH · Search NASA

Results for “numerical stabilization”

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 325 records · Page 18

All sources NASA NTRS DOE OSTI

NASA NTRS1989

An Entropy-Based Approach to Nonlinear Stability

Many numerical methods used in computational fluid dynamics (CFD) incorporate an artificial dissipation term to suppress spurious oscillations and control nonlinear instabilities. The same effect can be accomplished by using upwind techniques, sometimes augmented with limiters to form Total Variation Diminishing (TVD) schemes. An analysis based on numerical satisfaction of the second law of thermodynamics allows many such methods to be compared and improved upon. A nonlinear stability proof is given for discrete scalar equations arising from a conservation law. Solutions to such equations are bounded in the L sub 2 norm if the second law of thermodynamics is satisfied in a global sense over a periodic domain. It is conjectured that an analogous statement is true for discrete equations arising from systems of conservation laws. Analysis and numerical experiments suggest that a more restrictive condition, a positive entropy production rate in each cell, is sufficient to exclude unphysical phenomena such as oscillations and expansion shocks. Construction of schemes which satisfy this condition is demonstrated for linear and nonlinear wave equations and for the one-dimensional Euler equations.

Merriam, Marshal L.↗

NASA NTRS1985

Hydrodynamic stability problem formulated by numerical solutions of the Navier-Stokes equations

The problem of hydrodynamic stability and the transition from laminar to turbulent flows are reformulated by seeking numerical solutions of the full, unapproximated Navier-Stokes equations. This method differs significantly from the well known Orr-Sommerfeld equation approach. The oncoming laminar flow is disturbed by forced, time dependent perturbations. The magnitudes of these perturbations are arbitrary. Then, the ensuing spatial and temporal development of the imposed perturbations on the basic flow is calculated by direct numerical solutions of the time dependent Navier-Stokes equations. One of the main advantages of this method is its ability to simulate nonlinear processes. As a specific application of this technique to SSME (Space Shuttle Main Engine) flow configurations, computer programs have been written for the two dimensional flow over a backward facing step. This numerical code will be tested for operational use as part of a continued research collaboration effort with the NASA/MSFC counterparts.

Hyun, J. M.↗

DOE OSTI2018

Direct numerical simulation of flame stabilization assisted by autoignition in a reheat gas turbine combustor

A three-dimensional direct numerical simulation (DNS) is performed for a turbulent hydrogen-air flame, represented with detailed chemistry, stabilized in a model gas-turbine combustor. The combustor geometry consists of a mixing duct followed by a sudden expansion and a combustion chamber, which represents a geometrically simplified version of Ansaldo Energia’s GT26/GT36 sequential combustor design. In this configuration, a very lean blend of hydrogen and vitiated air is prepared in the mixing duct and convected into the combustion chamber, where the residence time from the inlet of the mixing duct to the combustion chamber is designed to coincide with the ignition delay time of the mixture. The results show that when the flame is stabilized at its design position, combustion occurs due to both autoignition and flame propagation (deflagration) modes at different locations within the combustion chamber. A chemical explosive mode analysis (CEMA) reveals that most of the fuel is consumed due to autoignition in the bulk-flow along the centerline of the combustor, and lower amounts of fuel are consumed by flame propagation near the corners of the sudden expansion, where the unburnt temperature is reduced by the thermal wall boundary layers. An unstable operating condition is also identified, wherein periodic auto-ignition events occur within the mixing duct. These events appear upstream of the intended stabilization position, due to positive temperature fluctuations induced by pressure waves originating from within the combustion chamber. This DNS investigation represents the initial step of a comprehensive research effort aimed at gaining detailed physical insight into the rate-limiting processes that govern the sequential combustor behavior and avoid the insurgence of the off-design auto-ignition events.

20 FOSSIL-FUELED POWER PLANTS↗

DOE OSTI2018

Inclined porous medium convection at large Rayleigh number

High-Rayleigh-number ( $Ra$ ) convection in an inclined two-dimensional porous layer is investigated using direct numerical simulations (DNS) and stability and variational upper-bound analyses. When the inclination angle $\unicode[STIX]{x1D719}$ of the layer satisfies $0^{\circ }<\unicode[STIX]{x1D719}\lesssim 25^{\circ }$ , DNS confirm that the flow exhibits a three-region wall-normal asymptotic structure in accord with the strictly horizontal ( $\unicode[STIX]{x1D719}=0^{\circ }$ ) case, except that as $\unicode[STIX]{x1D719}$ is increased the time-mean spacing between neighbouring interior plumes also increases substantially. Both DNS and upper-bound analysis indicate that the heat transport enhancement factor (i.e. the Nusselt number) $Nu\sim CRa$ with a $\unicode[STIX]{x1D719}$ -dependent prefactor $C$ . When $\unicode[STIX]{x1D719}>\unicode[STIX]{x1D719}_{t}$ , however, where $30^{\circ }<\unicode[STIX]{x1D719}_{t}<32^{\circ }$ independently of $Ra$ , the columnar flow structure is completely broken down: the flow transitions to a large-scale travelling-wave convective roll state, and the heat transport is significantly reduced. To better understand the physics of inclined porous medium convection at large $Ra$ and modest inclination angles, a spatial Floquet analysis is performed, yielding predictions of the linear stability of numerically computed, fully nonlinear steady convective states. The results show that there exist two types of instability when $\unicode[STIX]{x1D719}\neq 0^{\circ }$ : a bulk-mode instability and a wall-mode instability, consistent with previous findings for $\unicode[STIX]{x1D719}=0^{\circ }$ (Wen et al. , J. Fluid Mech. , vol. 772, 2015, pp. 197–224). The background flow induced by the inclination of the layer intensifies the bulk-mode instability during its subsequent nonlinear evolution, thereby favouring increased spacing between the interior plumes relative to that observed in convection in a horizontal porous layer.

Mechanics↗

DOE OSTI2014

Stability analysis of the Euler discretization for SIR epidemic model

In this paper we consider a discrete SIR epidemic model obtained by the Euler method. For that discrete model, existence of disease free equilibrium and endemic equilibrium is established. Sufficient conditions on the local asymptotical stability of both disease free equilibrium and endemic equilibrium are also derived. It is found that the local asymptotical stability of the existing equilibrium is achieved only for a small time step size h. If h is further increased and passes the critical value, then both equilibriums will lose their stability. Our numerical simulations show that a complex dynamical behavior such as bifurcation or chaos phenomenon will appear for relatively large h. Both analytical and numerical results show that the discrete SIR model has a richer dynamical behavior than its continuous counterpart.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1990

Numerical prediction of turbulent flame stability in premixed/prevaporized (HSCT) combustors

A numerical analysis of combustion instabilities that induce flashback in a lean, premixed, prevaporized dump combustor is performed. KIVA-II, a finite volume CFD code for the modeling of transient, multidimensional, chemically reactive flows, serves as the principal analytical tool. The experiment of Proctor and T'ien is used as a reference for developing the computational model. An experimentally derived combustion instability mechanism is presented on the basis of the observations of Proctor and T'ien and other investigators of instabilities in low speed (M less than 0.1) dump combustors. The analysis comprises two independent procedures that begin from a calculated stable flame: The first is a linear increase of the equivalence ratio and the second is the linear decrease of the inflow velocity. The objective is to observe changes in the aerothermochemical features of the flow field prior to flashback. It was found that only the linear increase of the equivalence ratio elicits a calculated flashback result. Though this result did not exhibit large scale coherent vortices in the turbulent shear layer coincident with a flame flickering mode as was observed experimentally, there were interesting acoustic effects which were resolved quite well in the calculation. A discussion of the k-e turbulence model used by KIVA-II is prompted by the absence of combustion instabilities in the model as the inflow velocity is linearly decreased. Finally, recommendations are made for further numerical analysis that may improve correlation with experimentally observed combustion instabilities.

Winowich, Nicholas S.↗

NASA NTRS1977

Efficient numerical treatment of periodic systems with application to stability problems

Two efficient numerical methods for dealing with the stability of linear periodic systems are presented. Both methods combine the use of multivariable Floquet-Liapunov theory with an efficient numerical scheme for computing the transition matrix at the end of one period. The numerical properties of these methods are illustrated by applying them to the simple parametric excitation problem of a fixed end column. The practical value of these methods is shown by applying them to some helicopter rotor blade aeroelastic and structural dynamics problems. It is concluded that these methods are numerically efficient, general and practical for dealing with the stability of large periodic systems.

Friedmann, P.↗

NASA NTRS1979

Explicit numerical formulas of improved stability and accuracy for the solution of parabolic equations

A class of explicit numerical formulas which involve next nearest neighbor as well as nearest neighbor points are explored in this paper. These formulas are formal approximations to the linear parabolic partial-differential equation of first order in time and second order in distance. It was found that some of these formulas can employ time steps as much as four times that for the conventional explicit technique without becoming unstable. Others showed improved accuracy for a given time step and spatial grid spacing. One formula achieved a steady-state solution of specified accuracy for an example problem in less than 4 percent of the total computational time required by the conventional explicit technique.

Olstad, W. B.↗

DOE OSTI2016

Stability of exact solutions of the nonlinear Schrödinger equation in an external potential having supersymmetry and parity-time symmetry

We discuss the stability properties of the solutions of the general nonlinear Schrödinger equation (NLSE) in 1+1 dimensions in an external potential derivable from a parity-time (${ \mathcal P }{ \mathcal T }$) symmetric superpotential $W(x)$ that we considered earlier, Kevrekidis et al (2015 Phys. Rev. E 92 042901). In particular we consider the nonlinear partial differential equation $\{{\rm{i}}\,{\partial }_{t}+{\partial }_{x}^{2}-{V}^{-}(x)+| \psi (x,t){| }^{2\kappa }\}\,\psi (x,t)=0,$ for arbitrary nonlinearity parameter κ. We study the bound state solutions when ${V}^{-}(x)\,=(1/4-{b}^{2}){\text{sech}}^{2}(x)$, which can be derived from two different superpotentials $W(x)$, one of which is complex and ${ \mathcal P }{ \mathcal T }$ symmetric. Using Derrick's theorem, as well as a time dependent variational approximation, we derive exact analytic results for the domain of stability of the trapped solution as a function of the depth b 2 of the external potential. We compare the regime of stability found from these analytic approaches with a numerical linear stability analysis using a variant of the Vakhitov–Kolokolov (V–K) stability criterion. The numerical results of applying the V–K condition give the same answer for the domain of stability as the analytic result obtained from applying Derrick's theorem. Our main result is that for $\kappa \gt 2$ a new regime of stability for the exact solutions appears as long as $b\gt {b}_{{\rm{crit}}}$, where ${b}_{{\rm{crit}}}$ is a function of the nonlinearity parameter κ. In the absence of the potential the related solitary wave solutions of the NLSE are unstable for $\kappa \gt 2$.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1985

A simple model for the interaction between vertical eddy heat fluxes and static stability

A numerical model for studying the interaction of vertical eddy heat fluxes and vertical temperature structure in midlatitude regions is described. The temperature profile for the model was derived from calculations of the equilibrium among heating rates in simplified representations of large-scale vertical eddy heat flux, moist convection and radiation. An eddy flux profile is calculated based on the quasi-geostrophic, liner baroclinic instability of a single wave. Model equilibrium states for summer and winter conditions are compared with observations, and the results are discussed in detail.

Gutowski, W. J., Jr.↗

NASA NTRS1995

Unsteady numerical simulations of the stability and dynamics of flames

In this report we describe the research performed at the Naval Research Laboratory in support of the NASA Microgravity Science and Applications Program over the past three years (from Feb. 1992) with emphasis on the work performed since the last microgravity combustion workshop. The primary objective of our research is to develop an understanding of the differences in the structure, stability, dynamics and extinction of flames in earth gravity and in microgravity environments. Numerical simulations, in which the various physical and chemical processes can be independently controlled, can significantly advance our understanding of these differences. Therefore, our approach is to use detailed time-dependent, multi-dimensional, multispecies numerical models to perform carefully designed computational experiments. The basic issues we have addressed, a general description of the numerical approach, and a summary of the results are described in this report. More detailed discussions are available in the papers published which are referenced herein. Some of the basic issues we have addressed recently are (1) the relative importance of wall losses and gravity on the extinguishment of downward-propagating flames; (2) the role of hydrodynamic instabilities in the formation of cellular flames; (3) effects of gravity on burner-stabilized flames, and (4) effects of radiative losses and chemical-kinetics on flames near flammability limits. We have also expanded our efforts to include hydrocarbon flames in addition to hydrogen flames and to perform simulations in support of other on-going efforts in the microgravity combustion sciences program. Modeling hydrocarbon flames typically involves a larger number of species and a much larger number of reactions when compared to hydrogen. In addition, more complex radiation models may also be needed. In order to efficiently compute such complex flames recent developments in parallel computing have been utilized to develop a state-of-the-art parallel flame code. This is discussed below in some detail after a brief discussion of the numerical models.

Kailasanath, K.↗

DOE OSTI2015

Numerical study on the stability of weakly collisional plasma in E×B fields

Plasma stability in weakly collisional plasmas in the presence of E×B fields is studied with numerical simulations. Different types of ion-neutral collisions are considered in a fully magnetized regime. We study the influence of ion-neutral collisions and the role of collision types on the stability of plasma. It is found that the stability of plasma depends on the type of ion-neutral collisions, with the plasma being unstable for charge exchange collisions, and stable for the elastic scattering. The analysis focuses on the temporal evolution of the velocity phase space, RMS values of the potential fluctuations, and coherent structures in potential densities. For the unstable case, we observe growth and propagation of electrostatic waves. Simulations are performed with a three-dimensional electrostatic particle in cell code.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS2009

Towards an "All Speed" Unstructured Upwind Scheme

In the authors previous studies [1], a time-accurate, upwind finite volume method (ETAU scheme) for computing compressible flows on unstructured grids was proposed. The scheme is second order accurate in space and time and yields high resolution in the presence of discontinuities. The scheme features a multidimensional limiter and multidimensional numerical dissipation. These help to stabilize the numerical process and to overcome the annoying pathological behaviors of upwind schemes. In the present paper, it will be further shown that such multidimensional treatments also lead to a nearly all-speed or Mach number insensitive upwind scheme. For flows at very high Mach number, e.g., 10, local numerical instabilities or the pathological behaviors are suppressed, while for flows at very low Mach number, e.g., 0.02, computation can be directly carried out without invoking preconditioning. For flows in different Mach number regimes, i.e., low, medium, and high Mach numbers, one only needs to adjust one or two parameters in the scheme. Several examples with low and high Mach numbers are demonstrated in this paper. Thus, the ETAU scheme is applicable to a broad spectrum of flow regimes ranging from high supersonic to low subsonic, appropriate for both CFD (computational fluid dynamics) and CAA (computational aeroacoustics).

Loh, Ching Y.↗

NASA NTRS1984

A numerical study of the stability of radiative shocks

Attention is given to the oscillatory instability of optically thin radiative shocks in time-dependent numerical calculations of accretion flows onto degenerate dwarfs. The present nonlinear calculations yield good quantitative agreement with the linear results obtained for oscillation frequencies, damping rates, and critical alpha-values. The fundamental mode and the first overtone in the shock radius and luminosity variations can be clearly identified, and evidence is sometimes seen for the second overtone. Time-dependent calculations are also performed which include additional physics relevant to degenerate dwarf accretion, such as electron thermal conduction, unequal electron and ion temperatures, Compton cooling, and relativistic corrections to the bremsstrahlung cooling law. All oscillatory modes are found to be damped, and hence stable, in the case of a 1-solar mass white dwarf accreting in spherical symmetry.

Imamura, J. N.↗

NASA NTRS1981

Numerical experiments on the stability of preplanetary disks

Gravitational stability of gaseous protostellar disks is relevant to theories of planetary formation. Stable gas disks favor formation of planetesimals by the accumulation of solid material; unstable disks allow the possibility of direct condensation of gaseous protoplanets. This paper presents the results of numerical experiments designed to test the stability of thin disks against large-scale, self-gravitational disruption. It is found that a disk as massive as 1 solar mass, surrounding a 1 solar mass protostar, can be stable against long-wavelength gravitational disruption if its temperature is about 300 K or greater. Stability of a cooler disk requires that it be less massive, but even at 100 K a stable disk can have an appreciable fraction (about 1/3) of a solar mass.

Cassen, P. M.↗

DOE OSTI2016

Numerical study on the Welander oscillatory natural circulation problem using high-order numerical methods

In this study, the classical Welander’s oscillatory natural circulation problem is investigated using high-order numerical methods. As originally studied by Welander, the fluid motion in a differentially heated fluid loop can exhibit stable, weakly instable, and strongly instable modes. A theoretical stability map has also been originally derived from the stability analysis. Numerical results obtained in this paper show very good agreement with Welander’s theoretical derivations. For stable cases, numerical results from both the high-order and low-order numerical methods agree well with the non-dimensional flow rate analytically derived. The high-order numerical methods give much less numerical errors compared to the low-order methods. For stability analysis, the high-order numerical methods could perfectly predict the stability map, while the low-order numerical methods failed to do so. For all theoretically unstable cases, the low-order methods predicted them to be stable. The result obtained in this paper is a strong evidence to show the benefits of using high-order numerical methods over the low-order ones, when they are applied to simulate natural circulation phenomenon that has already gain increasing interests in many future nuclear reactor designs.

42 ENGINEERING↗

DOE OSTI2020

Stability boundary and optimal operating parameter identification in milling using Bayesian learning

Here, this paper describes a novel Bayesian learning approach for stability boundary and optimal parameter identification in milling without the knowledge of the underlying tool dynamics or material cutting force coefficients. The paper is divided into two parts. First, a Bayesian learning method for stability lobe identification using test results is described. Each axial depth and spindle speed combination is characterized by a probability of stability which is updated using Bayes’ rule when a test result (stable or unstable) is made available. A novel likelihood function is defined which incorporates knowledge of the stability behavior. Numerical results show convergence to the analytical stability lobe diagram. Second, an adaptive experimental strategy to identify optimal operating parameters that maximize material removal rate is described. Numerical evaluation shows convergence to the optimal operating point with error less than 15 % within ten tests on average. The approach is validated using experimental results. Results show that the proposed method is an efficient and robust learning method to identify the stability lobe diagram and optimal operating parameters with a limited number of tests/data points.

Bayesian machine learning↗

NASA NTRS1978

Numerical experiments on the stability of disklike galaxies

It is noted that the stability properties of available galactic models are not consistent with presumptions based on observation and that axisymmetric disk systems with velocity dispersions like those of the Galaxy display major changes in form on a dynamical time scale. A report is given on a series of numerical experiments carried out as part of a systematic search for purely self-consistent disk galaxy models that might undergo little change over the time of several galactic rotations. The stability problem is reviewed, the method and calculations used are described, and initial conditions are outlined. The problem of estimating growth rates is discussed, and growth rates are evaluated for various disturbances. Experimental results are presented concerning growths of axisymmetric disturbances in 'cold' systems as well as disturbances with m = 1, 2, and 3. It is found that m = 2 disturbances are the most virulent and should dominate when all m values are allowed, but that the most rapidly growing disturbances are inhibited when several m values are allowed.

Miller, R. H.↗