Engineering Papers⌕ Search

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.

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

At least 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
NASA NTRS1989

Numerical modeling of a vortex stabilized arcjet thruster

A numerical method to solve the equations governing a vortex stabilized arcjet thruster is being developed. The model will allow the effects of swirling flow and geometry on arcjet thruster performance to be determined. The propellant flow equations, i.e., the axisymmetric, thin layer, Navier-Stokes equations, are solved using a Gauss-Siedel line-relaxation procedure. An implicit FTCS method is used to solve the electromagnetic field equations. A grid generation scheme was developed for an arbitrary arcjet geometry. The model will allow the radial and axial components of velocity and current distributions to be determined from a region upstream of the cathode, through the constrictor, to the exit plane of the nozzle.

Keith, Theo G., Jr.↗
DOE OSTI2026

Hourglass control in staggered-grid hydrodynamics using virtual element stabilization techniques

Numerical simulations using the staggered-grid hydrodynamics (SGH) discretization suffer from hourglass instabilities. In this work, we develop a stabilization method to suppress the hourglass instabilities using techniques from the virtual element method (VEM). The stiffness matrix of the VEM consists of two terms: the consistency matrix which is rank deficient and the stability matrix. Here, we first show that in two dimensions and on general polygons, the stiffness matrix of the SGH is identical to the consistency matrix of the linear VEM for both the diffusion equation and the linear elasticity equation. These analyses explain the origin of the hourglass instabilities of the SGH discretization method, and establish a theoretical foundation for our proposed stabilization method by augmenting the stiffness matrix of the SGH discretization using the VEM stability matrix. Then, we present numerical examples using Lagrangian SGH simulations. The numerical experiments demonstrate that the proposed VEM stabilization method is effective at eliminating hourglass modes in the SGH discretization.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Comparing numerical accuracy and stability for different horizontal discretizations in MPAS-Ocean

This manuscript investigates the effectiveness of two possible horizontal discretizations for the global ocean model MPAS-Ocean, both applied to Spherical Centroidal Voronoi Tessellations (SCVTs). The first discretization is TRiSK, a C-grid, finite-volume method, that possesses many desirable mimetic properties, but has a low order accuracy. The second discretization was introduced for the first time by Peixoto (2016), and consists of modifications to the TRiSK scheme designed to achieve at least first-order accuracy in the L ∞ norm, with the loss of some of the mimetic properties. Tests on shallow-water and primitive-equation models show that the scheme due to Peixoto is indeed more accurate, but presents stability issues with respect to TRiSK. Here, TRiSK is indeed found to be often more stable in time and more robust with respect to errors in the geometric properties of the grid.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1975

Numerical implications of stabilization by the use of integrals

Liapunov or energy restraint methods for dynamic stabilization in two body motion perturbation problems are considered. Results of computerized orbital stabilization estimates show that the application of energy restraint prevents the occurrence of consistent timing errors in the stepwise integration of equations of motion for a nearly circular orbit.

Beaudet, P. R.↗
NASA NTRS2018

Far-field Lorenz-Mie Scattering in an Absorbing Host Medium. II: Improved Stability of the Numerical Algorithm

A recently developed FORTRAN program computing far-field optical observables for spherical particles in an absorbing medium has exhibited numerical instability arising when the product of the particle vacuum size parameter and the imaginary part of the refractive index of the host becomes sufficiently large. We offer a simple analytical explanation of this instability and propose a compact numerical algorithm for the stable computation of Lorenz–Mie coefficients based on an upward recursion formula for spherical Hankel functions of a complex argument. Extensive tests confirm an excellent accuracy of this algorithm approaching machine precision. The improved public-domain FORTRAN program is available at https://www.giss.nasa.gov/staff/mmishchenko/Lorenz-Mie.html.

Lorenz-Mie theory↗
NASA NTRS1987

Ponderomotive stabilization of flute modes in mirrors Feedback control and numerical results

Ponderomotive stabilization of rigid plasma flute modes is numerically investigated by use of a variational principle, for a simple geometry, without eikonal approximation. While the near field of the studied antenna can be stabilizing, the far field has a small contribution only, because of large cancellation by quasi mode-coupling terms. The field energy for stabilization is evaluated and is a nonnegligible fraction of the plasma thermal energy. A new antenna design is proposed, and feedback stabilization is investigated. Their use drastically reduces power requirements.

Similon, P. L.↗
DOE OSTI2026

Hybrid Symbolic-Numerical Modeling and Parametric Stability Analysis of DC–AC Power Systems

Hybrid DC-AC power systems integrating diverse inverter-based resources (IBRs) and multi-terminal high-voltage direct current (MTDC) networks represent a promising architecture for future power grids, while introducing challenges for modeling, stability analysis, and control design. This paper develops a hybrid symbolic-numerical modeling framework and tool to characterize the parametric small-signal stability of DC-AC coupled power systems. The proposed approach constructs parametric state-space models to enable efficient representation of system dynamics under varying control parameters and network configurations, with target parameters retained as symbolic variables and the remainder treated numerically. The stability analysis framework covers eigenvalue, sensitivity, and stability boundary and region characterization. Enhanced linear matrix inequality (LMI) techniques are proposed to directly certify small-signal stability over regions of parameter space while also reducing the conservativeness and computational burden. The resulting tools and frameworks enable rapid parametric model construction across diverse grid conditions, thereby facilitating stability-informed control and operation in complex DC–AC power systems.

DC–AC power systems↗
NASA NTRS1986

Numerical-perturbation technique for stability of flat-plate boundary layers with suction

A numerical-perturbation scheme is proposed for determining the stability of flows over plates with suction through a finite number of porous suction strips. The basic flow is calculated as the sum of the Blasius flow and closed-form linearized triple-deck solutions of the flow due to the strips. A perturbation technique is used to determine the increment a(ij) in the complex wavenumber at a given location x(j) due to the presence of a strip centered at x(i). The end result is a set of influence coefficients that can be used to determine the growth rates and amplification factors for any suction levels without repeating the calculations. The numerical-perturbation results are verified by comparison with interacting boundary layers for the case of six strips and the experimental data of Reynolds and Saric for single- and multiple-strip configurations. The influence coefficient form of the solution suggests a scheme for optimizing the strip configuration. The results show that one should concentrate the suction near branch I of the neutral stability curve, a conclusion verified by the experiments.

Reed, H. L.↗
DOE OSTI2020

F4‐TCNQ as an Additive to Impart Stretchable Semiconductors with High Mobility and Stability

Abstract Numerous strategies are developed to impart stretchability to polymer semiconductors. Although these methods improve the ductility, mobility, and stability of such stretchable semiconductors, they nonetheless still need further improvement. Here, it is shown that 2,3,5,6‐tetrafluoro‐7,7,8,8‐tetracyanoquinodimethane (F4‐TCNQ) is an effective molecular additive to tune the properties of a diketopyrrolopyrrole‐based (DPP‐based) semiconductor. Specifically, the addition of F4‐TCNQ is observed to improve the ductility of the semiconductor by altering the polymer’s microstructures and dynamic motions. As a p‐type dopant additive, F4‐TCNQ can also effectively enhance the mobility and stability of the semiconductor through changing the host polymer’s packing structures and charge trap passivation. Upon fabricating fully stretchable transistors with F4‐TCNQ‐DPP blended semiconductor films, it is observed that the resulting stretchable transistors possess one of the highest initial mobility of 1.03 cm 2 V −1 s −1 . The fabricated transistors also exhibit higher stability (both bias and environmental) and mobility retention under repeated strain, compared to those without F4‐TCNQ additive. These findings offer a new direction of research on stretchable semiconductors to facilitate future practical applications.

Mun, Jaewan↗
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.↗
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.↗