SEARCH · Engineering Papers

Results for “Numerical 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 451 records · Page 25

All sources NASA NTRS DOE OSTI

NASA NTRS1989

On some numerical schemes for transonic flow problems

New second-order-accurate finite-difference approximations for a class of nonlinear PDEs of mixed type, which includes the two-dimensional low-frequency transonic small-disturbance equation (TSD) and full-potential equation (FP), are presented. For the TSD equation, the scheme is implemented via a time-splitting algorithm; the inclusion of flux limiters keeps the total variation nonincreasing and eliminates spurious oscillations near shocks. Global-linear-stability, TVD, and entropy-stability results are proved. Numerical results for the flow over a thin airfoil are presented. Current techniques used to solve the TSD equation may easily be extended to second-order accuracy by this method. For the FP equation, the new scheme requires no subsonic/supersonic switching and no numerical flux biasing. Global linear stability for all values of the Mach number is proved.

Mostrel, Marco Mosche↗

NASA NTRS1974

Static, stability, and dynamic analysis of shells of revolution by numerical integration - A comparison

Recent innovations in digital computer technology have enabled engineers to analyze shell structures of complex configurations without unduly restrictive approximations. An attempt is made to compare the various programs now generally available from the point of view of the advantages of the relative technique utilized, as well as the programmed state of the art. Many of the comparisons are based on the sample problems solved by the STARS-2 system of programs. These examples indicate both the structural detail which can be analyzed by, and the analytical capabilities available in, the numerical shell-of-revolution programs. All advantages and differences are demonstrated by use of solutions for realistic shell problems in the areas of statics, stability, vibrations, and dynamic response of shells subjected to time-dependent loadings.

Svalbonas, V.↗

DOE OSTI2025

Physics-based stabilized finite element approximations of the Poisson–Nernst–Planck equations

We present and analyze two stabilized finite element methods for solving numerically the Poisson–Nernst–Planck equations. The stabilization we consider is carried out by using a shock detector and a discrete graph Laplacian operator for the ion equations, whereas the discrete equation for the electric potential need not be stabilized. Discrete solutions stemmed from the first algorithm preserve both maximum and minimum discrete principles. For the second algorithm, its discrete solutions are conceived so that they hold discrete principles and obey an entropy law provided that an acuteness condition is imposed for meshes. Remarkably the latter is found to be unconditionally stable. We validate our methodology through transient numerical experiments that show convergence toward steady-state solutions.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1999

Influence of Buoyant Convection on the Stability of Enclosed Laminar Flames

Enclosed diffusion flames are commonly found in practical combustion systems, such as the power-plant combustor, gas turbine combustor, and jet engine after-burner. In these systems, fuel is injected into a duct with a co-flowing or cross-flowing air stream. In combustors, this flame is anchored at the burner (i.e., fuel jet inlet) unless adverse conditions cause the flame to lift off or blow out. Investigations of burner stability study the lift off, reattachment, and blow out of the flame. There have been numerous studies of flame stability. Relatively few studies have investigated the stability of flames with an oxidizer co-flow, compared with the number of studies on (nearly) free jet diffusion flames. The air flow around the fuel jet can significantly alter the lift off, reattachment and blow out of the jet diffusion flame. In normal gravity, however, the effects of the air flow on flame stability are often complicated by the presence of buoyant convection. A comparison of normal-gravity and microgravity flames can provide clear indication of the influence of forced and buoyant flows on the flame stability. The overall goal of the Enclosed Laminar Flames (ELF) research, described at the following URL site: http://zeta.lerc.nasa.gov/expr/elf.htm, is to improve our understanding of the effects of buoyant convection on the structure and stability of co-flow diffusion flames.

Brooker, John E.↗

NASA NTRS1998

Stability of Capillary Surfaces in Rectangular Containers: The Right Square Cylinder

The linearized governing equations for an ideal fluid are presented for numerical analysis for the stability of free capillary surfaces in rectangular containers against unfavorable disturbances (accelerations,i.e. Rayleigh-Taylor instability). The equations are solved for the case of the right square cylinder. The results are expressed graphically in term of a critical Bond number as a function of system contact angle. A critical wetting phenomena in the corners is shown to significantly alter the region of stability for such containers in contrast to simpler geometries such as the right circular cylinder or the infinite rectangular slot. Such computational results provide additional constraints for the design of fluids systems for space-based applications.

Weislogel, M. M.↗

NASA NTRS1985

A numerical study of global axisymmetric states and their stability to quasi-geostrophic disturbances

A steady, axisymmetric model of the general circulation is developed as a basis for climate stability studies. The model includes the effects of heating, rotation, and internal friction, but neglects topography. It is assumed that the axisymmetric flow may be modeled by making the Boussinesq and deep convection approximations. The hydrostatic assumption is not made, thus permitting the advective terms to be included in the vertical equation of motion. The initial set of five primitive equations is reduced to three equations in terms of the zonal velocity, meridional streamfunction, and the potential temperature perturbation. The application of the Boussinesq, deep convection, and quasi-geostrophic assumptions limits the ranges of the heating and rotation rates. For values not too far from typical atmospheric values, the model produces a stability boundary separating Hadley from Rossby flow. The boundary is characterized by a particular value of vertical wind shear, which suggests that baroclinic instability is the primary mechanism for the loss of stability. The initial growth rates are largest for longitudinal waves 4-7, also in agreement with studies of baroclinic instability.

Henderson, H. W.↗

NASA NTRS1987

Experimental and numerical analyses of laminar boundary-layer flow stability over an aircraft fuselage forebody

Fuelled by a need to reduce viscous drag of airframes, significant advances have been made in the last decade to design lifting surface geometries with considerable amounts of laminar flow. In contrast to the present understanding of practical limits for natural laminar flow over lifting surfaces, limited experimental results are available examining applicability of natural laminar flow over axisymmetric and nonaxisymmetric fuselage shapes at relevantly high length Reynolds numbers. The drag benefits attainable by realizing laminar flow over nonlifting aircraft components such as fuselages and nacelles are shown. A flight experiment to investigate transition location and transition mode over the forward fuselage of a light twin engine propeller driven airplane is examined.

Vijgen, Paul M. H. W.↗

DOE OSTI2021

On the Stochastic Stability of Deep Markov Models

Deep Markov models (DMM) are generative models which are scalable and expressive generalization of Markov models for representation, learning, and inference problems. However, the fundamental stochastic stability guarantees of such models have not been thoroughly investigated. In this paper, we present a novel stability analysis method and provide sufficient conditions of DMM's stochastic stability. The proposed stability analysis is based on the contraction of probabilistic maps modeled by deep neural networks. We make connections between the spectral properties of neural network's weights and different types of used activation function on the stability and overall dynamic behavior of DMMs with Gaussian distributions. Based on the theory, we propose a few practical methods for designing constrained DMMs with guaranteed stability. We empirically substantiate our theoretical results via intuitive numerical experiments using the proposed stability constraints.

Drgona, Jan↗

NASA NTRS2002

Local Discontinuous Galerkin Methods for Partial Differential Equations with Higher Order Derivatives

In this paper we review the existing and develop new continuous Galerkin methods for solving time dependent partial differential equations with higher order derivatives in one and multiple space dimensions. We review local discontinuous Galerkin methods for convection diffusion equations involving second derivatives and for KdV type equations involving third derivatives. We then develop new local discontinuous Galerkin methods for the time dependent bi-harmonic type equations involving fourth derivatives, and partial differential equations involving fifth derivatives. For these new methods we present correct interface numerical fluxes and prove L(exp 2) stability for general nonlinear problems. Preliminary numerical examples are shown to illustrate these methods. Finally, we present new results on a post-processing technique, originally designed for methods with good negative-order error estimates, on the local discontinuous Galerkin methods applied to equations with higher derivatives. Numerical experiments show that this technique works as well for the new higher derivative cases, in effectively doubling the rate of convergence with negligible additional computational cost, for linear as well as some nonlinear problems, with a local uniform mesh.

Yan, Jue↗

DOE OSTI2023

Comparison of MHD stability properties between QH-mode and ELMy H-mode plasmas by considering plasma rotation and ion diamagnetic drift effects

Magnetohydrodynamic (MHD) stability at tokamak edge pedestal in a quiescent H-mode (QH-mode) and type-I ELMy H-mode plasmas in DIII-D experiment was analyzed by considering plasma rotation and ion diamagnetic drift effects. QH-mode plasma is marginally stable to kink/peeling mode (K/PM), but ELMy H-mode one is almost unstable to peeling-ballooning mode (PBM). It was identified that there are three physics features responsible for the difference in the MHD stability properties between QH-mode plasma and ELMy H-mode one. These are the distance of pedestal foot from the last closed flux surface (LCFS), the amount of the ion diamagnetic drift frequency at pedestal, and impact of coupled rotation and ion diamagnetic drift effects. These features were confirmed through the numerical experiments that the stability properties of the QH-mode plasma can be changed to that of the ELMy H-mode one by shifting the plasma profiles inward in the radial direction and halving the ion diamagnetic drift frequency. The reasons of the change in the stability properties are thought as that K/PM is stabilized due to the inward shift of the bootstrap current profile, and PBM is destabilized due to the reduction of the coupled rotation and ion diamagnetic drift stabilizing effect. Importance of these features was validated through numerical experiments with experimental data of other QH-mode plasmas in DIII-D. All the results show that MHD stability properties of QH-mode plasma can be obtained in case that pedestal foot is close to LCFS, ion diamagnetic drift frequency is large due to high ion temperature, and strong rotation shear exists near pedestal.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1987

A finite volume Euler calculation of the aerodynamics of transonic airfoil-vortex interaction

Unsteady inviscid transonic airfoil-vortex interaction is numerically analyzed by solving the two-dimensional unsteady Euler equations in integral form using a finite volume scheme. The solution procedure is based on an explicit Runge-Kutta time-stepping scheme wherein the spatial terms are central-differenced and a combination of second- and fourth-differences in the flow variables is used to form the numerical dissipation terms to stabilize the scheme. A velocity decomposition technique is applied to alleviate the problem of vortex diffusion by the numerical dissipation terms and to treat the interaction of a Rankine vortex with an airfoil accurately. Results obtained are compared with available numerical data.

Damodaran, Murali↗

NASA NTRS2005

The Apollo Capsule Optimization for Improved Stability and Computational/Experimental Data Comparisons

Numerical optimization was employed on the Apollo Command Module to modify its external shape. The Apollo Command Module (CM) that was used on all NASA human space flights during the Apollo Space Program is stable and trimmed in an apex forward (alpha of approximately 40 to 80 degrees) position. This poses a safety risk if the CM separates from the launch tower during abort. Optimization was employed on the Apollo CM to remedy the undesirable stability characteristics of the configuration. Geometric shape changes were limited to axisymmetric modifications that altered the radius of the apex (R(sub A)), base radius (R(sub O)), corner radius (R(sub C)), and the cone half angle (theta), while the maximum diameter of the CM was held constant. The results of multipoint optimization on the CM indicated that the cross-range performance can be improved while maintaining robust apex-aft stability with a single trim point. Navier-Stokes computations were performed on the baseline and optimized configurations and confirmed the Euler-based optimization results. Euler Analysis of ten alternative CM vehicles with different values of the above four parameters are compared with the published experimental results of numerous wind tunnel tests during the late 1960's. These comparisons cover a wide Mach number range and a full 180-degree pitch range and show that the Euler methods are capable of fairly accurate force and moment computations and can separate the vehicle characteristics of these ten alternative configurations.

Cliff, Susan E.↗

NASA NTRS1993

Effects of gas on the global stability of galactic disks - Radial flows

We study numerically the effect of gas on the global stability of a two-component self-gravitating galactic disk embedded in a live halo. The stars are evolved by using a 3D collisionless N-body code, and the gas is represented by an ensemble of finite size inelastic particles. The gravitational interaction of stars and gas is calculated using a TREE method. We find that the evolution of the gaseous distribution in the globally unstable disks can be described by two different regimes. When the gas mass fraction is less than about 10 percent, the gas is channeled toward the galactic center by a growing stellar bar. For higher gas fractions, the gas becomes highly inhomogeneous, and the bar instability in the disk is heavily damped. The gas falls toward the inner kpc due to dynamical friction. Domains of both regimes depend on the efficiency of dissipation in the gas. We also discuss the relevance of the Jeans instability and give an empirical criterion for the global bar instability in a two-component self-gravitating disk.

Shlosman, Isaac↗

DOE OSTI2023

On the relationship between the multi-region relaxed variational principle and resistive inner-layer theory

We show that the variational energy principle of the multi-region relaxed magnetohydrodynamic (MRxMHD) model can be used to predict finite-pressure linear tearing instabilities. In this model, the plasma volume is sliced into sub-volumes separated by 'ideal interfaces', and in each volume the magnetic field relaxes to a Taylor state, where the pressure gradient $\nabla p = 0$. The MRxMHD model is implemented in the Stepped-Pressure Equilibrium Code (SPEC) so that the equilibrium solution in each region is computed while preserving the force balance across the interfaces. As SPEC computes the Hessian matrix (a discretized stability matrix), the stability of an MRxMHD equilibrium can also be computed with SPEC. In this article, using SPEC, we investigate the effect of local pressure gradients and the $\nabla p = 0$ in the vicinity of the resonant surface of a tearing mode. For low-beta plasma, we have been able to illustrate a relationship between the resistive singular-layer theory (Coppi et al 1966 Nucl. Fusion6 101; Glasser et al 1975 Phys. Fluids18 875–88) and the MRxMHD model. Within the singular layer, the volume-averaged magnetic helicity and the flux-averaged toroidal flux are shown to be the invariants for the linear tearing modes in SPEC simulations. Our technique to compute MRxMHD stability is first tested numerically in a cylindrical tokamak and its application in toroidal geometry is demonstrated. We demonstrate an agreement between the stability boundary obtained with SPEC simulation and the resistive inner-layer theories.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

A novel approach for tetrahedral-element-based finite element simulations of anisotropic hyperelastic intervertebral disc behavior

Intervertebral discs are microstructurally complex spinal tissues that add greatly to the flexibility and mechanical strength of the human spine. Attempting to provide an adjustable basis for capturing a wide range of mechanical characteristics and to better address known challenges of numerical modeling of the disc, we present a robust finite-element-based model formulation for spinal segments in a hyperelastic framework using tetrahedral elements. We evaluate the model stability and accuracy using numerical simulations, with particular attention to the degenerated intervertebral discs and their likely skewed and narrowed geometry. To this end, 1) annulus fibrosus is modeled as a fiber-reinforced Mooney-Rivlin type solid for numerical analysis. 2) An adaptive state-variable dependent explicit time step is proposed and utilized here as a computationally efficient alternative to theoretical estimates. 3) Tetrahedral-element-based FE models for spinal segments under various loading conditions are evaluated for their use in robust numerical simulations. For flexion, extension, lateral bending, and axial rotation load cases, numerical simulations reveal that a suitable framework based on tetrahedral elements can provide greater stability and flexibility concerning geometrical meshing over commonly employed hexahedral-element-based ones for representation and study of spinal segments in various stages of degeneration.

59 BASIC BIOLOGICAL SCIENCES↗

NASA NTRS1984

Stability of a laminar premixed supersonic free shear layer with chemical reactions

The stability of a two-dimensional compressible supersonic flow in the wake of a flat plate is discussed. The fluid is a multi-species mixture which is undergoing finite rate chemical reactions. The spatial stability of an infinitesimal disturbance in the fluid is considered. Numerical solutions of the eigenvalue stability equations for both reactive and nonreactive supersonic flows are presented and discussed. The chemical reactions have significant influence on the stability behavior. For instance, a neutral eigenvalue is observed near the freestream Mach number of 2.375 for the nonreactive case, but disappears when the reaction is turned on. For reactive flows, the eigenvalues are not very dependent on the free stream Mach number.

Menon, S.↗

NASA NTRS1988

Cage stability analysis for SSME HPOTP bearings

A numerical model of cage motion (CAGEDYN) was used to analyze the stability of bearing cages in the Space Shuttle main engine (SSME) high pressure oxygen turbopump (HPOTP). The stability of existing bearing geometries, as well as perturbations of these geometries, was analyzed for various operating conditions. Results of the analyses show that some combinations of operating parameters, exacerbated by the sparse lubrication that exist in the HPOTP bearings, can cause unstable cage oscillations. Frequencies of cage oscillations were predicted by the CAGEDYN numerical model by Fourier analysis of predicted cage motions. Under conditions that cause unstable cage motion, high frequency oscillations were predicted that could cause premature cage failures.

Merriman, T. L.↗

NASA NTRS1983

Finite element formulation for transient heat treat problems

The macrothermomechanical behavior of materials subjected to rapid thermal or mechanical loading such as occurs in most heat treatments is described. The equations are developed for Lagrangian, Eulerian, and intermediary kinematic descriptions and are independent of the constitutive laws and the equation of state; they can be solved numerically for a specified material and boundary conditions. The coupled transport effects between dissipation and energy are included. The conventional linearized stability approach indicates the numerical procedure to be stable, with certain restriction on the time step size.

Mullen, R. L.↗