Search NASA⌕ Search

SEARCH · Search NASA

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Pseudospectral simulation of compressible turbulence using logarithmic variables

The direct numerical simulation of dissipative, highly compressible turbulent flow is performed using a pseudospectral Fourier technique. The governing equations are cast in a form where the important physical variables are the fluid velocity and the natural logarithms of the fluid density and temperature. Bulk viscosity is utilized to model polyatomic gases more accurately and to ensure numerical stability in the presence of strong shocks. Numerical examples include three-dimensional supersonic homogeneous turbulence and two-dimensional shock-turbulence interactions.

Shebalin, John V.↗
NASA NTRS1985

Characteristics pertaining to a stiffness cross-coupled Jeffcott model

Rotordynamic studies of complex systems utilizing multiple degree-of-freedom analysis have been performed to understand response, loads, and stability. In order to understand the fundamental nature of rotordynamic response, the Jeffcott rotor model has received wide attention. The purpose of this paper is to provide a generic rotordynamic analysis of a stiffness cross-coupled Jeffcott rotor model to illustrate characteristics of a second order stiffness-coupled linear system. The particular characteristics investigated were forced response, force vector diagrams, response orbits, and stability. Numerical results were achieved through a fourth order Runge-Kutta method for solving differential equations and the Routh Hurwitz stability criterion. The numerical results were verified to an exact mathematical solution for the steady state response.

Spanyer, K. L.↗
NASA NTRS1984

Stability Test for Transient-Temperature Calculations

Graphical test helps assure numerical stability of calculations of transient temperature or diffusion in composite medium. Rectangular grid forms basis of two-dimensional finite-difference model for heat conduction or other diffusion like phenomena. Model enables calculation of transient heat transfer among up to four different materials that meet at grid point.

Campbell, W.↗
NASA NTRS1982

Implicit numerical integration for periodic solutions of autonomous nonlinear systems

A change of variables that stabilizes numerical computations for periodic solutions of autonomous systems is derived. Computation of the period is decoupled from the rest of the problem for conservative systems of any order and for any second-order system. Numerical results are included for a second-order conservative system under a suddenly applied constant load. Near the critical load for the system, a small increment in load amplitude results in a large increase in amplitude of the response.

Thurston, G. A.↗
NASA NTRS1984

Stability and Convergence of Underintegrated Finite Element Approximations

The effects of underintegration on the numerical stability and convergence characteristics of certain classes of finite element approximations were analyzed. Particular attention is given to hourglassing instabilities that arise from underintegrating the stiffness matrix entries and checkerboard instabilities that arise from underintegrating constrain terms such as those arising from incompressibility conditions. A fundamental result reported here is the proof that the fully integrated stiffness is restored in some cases through a post-processing operation.

Oden, J. T.↗
NASA NTRS1978

Filtering and error analysis via the UDU super T covariance factorization

Kalman filter algorithms based on the UDU super T covariance factorization are discussed, with special attention given to algorithm implementation efficiency. A U-D-factored covariance error-analysis algorithm is formulated, and its efficiency and numerical stability are demonstrated in a representative orbit determination problem. The numerical results are compared with those obtained using covariance error-analysis formulas, and the comparison highlights the numerical superiority of the present algorithm. A byproduct of the U-D analysis is a highly efficient algorithm mechanization of the arbitrary gain covariance update formula.

Thornton, C. L.↗
NASA NTRS1986

Analysis of coaxial spray combustion flames and related numerical issues

An approach to the simulation of strongly coupled multiphase flows in combustion hardware is sketched and its unique requirements highlighted. An example of a successful application to a coaxial injector flame is presented. Furthermore, several numerical issues that tend to interact with the physics of the problem are discussed with special regard to their potential impact on the choices of numerical parameters by the analyst. These include the issues of stability, numerical diffusivity, stiffness, and boundary conditions. The theme of this paper focuses on the intriguing relationships among the grid, the solution algorithm, and the actual physical mechanisms themselves.

Liang, P. Y.↗
NASA NTRS1990

Numerical Methods For Chemically Reacting Flows

Issues related to numerical stability, accuracy, and resolution discussed. Technical memorandum presents issues in numerical solution of hyperbolic conservation laws containing "stiff" (relatively large and rapidly changing) source terms. Such equations often used to represent chemically reacting flows. Usually solved by finite-difference numerical methods. Source terms generally necessitate use of small time and/or space steps to obtain sufficient resolution, especially at discontinuities, where incorrect mathematical modeling results in unphysical solutions.

Leveque, R. J.↗
NASA NTRS1992

Runge-Kutta methods combined with compact difference schemes for the unsteady Euler equations

Recent development using compact difference schemes to solve the Navier-Stokes equations show spectral-like accuracy. A study was made of the numerical characteristics of various combinations of the Runge-Kutta (RK) methods and compact difference schemes to calculate the unsteady Euler equations. The accuracy of finite difference schemes is assessed based on the evaluations of dissipative error. The objectives are reducing the numerical damping and, at the same time, preserving numerical stability. While this approach has tremendous success solving steady flows, numerical characteristics of unsteady calculations remain largely unclear. For unsteady flows, in addition to the dissipative errors, phase velocity and harmonic content of the numerical results are of concern. As a result of the discretization procedure, the simulated unsteady flow motions actually propagate in a dispersive numerical medium. Consequently, the dispersion characteristics of the numerical schemes which relate the phase velocity and wave number may greatly impact the numerical accuracy. The aim is to assess the numerical accuracy of the simulated results. To this end, the Fourier analysis is to provide the dispersive correlations of various numerical schemes. First, a detailed investigation of the existing RK methods is carried out. A generalized form of an N-step RK method is derived. With this generalized form, the criteria are derived for the three and four-step RK methods to be third and fourth-order time accurate for the non-linear equations, e.g., flow equations. These criteria are then applied to commonly used RK methods such as Jameson's 3-step and 4-step schemes and Wray's algorithm to identify the accuracy of the methods. For the spatial discretization, compact difference schemes are presented. The schemes are formulated in the operator-type to render themselves suitable for the Fourier analyses. The performance of the numerical methods is shown by numerical examples. These examples are detailed. described. The third case is a two-dimensional simulation of a Lamb vortex in an uniform flow. This calculation provides a realistic assessment of various finite difference schemes in terms of the conservation of the vortex strength and the harmonic content after travelling a substantial distance. The numerical implementation of Giles' non-refelctive equations coupled with the characteristic equations as the boundary condition is discussed in detail. Finally, the single vortex calculation is extended to simulate vortex pairing. For the distance between two vortices less than a threshold value, numerical results show crisp resolution of the vortex merging.

Yu, Sheng-Tao↗
NASA NTRS1993

Stability of mixing layers

The research program for the first year of this project (see the original research proposal) consists of developing an explicit marching scheme for solving the parabolized stability equations (PSE). Performing mathematical analysis of the computational algorithm including numerical stability analysis and the determination of the proper boundary conditions needed at the boundary of the computation domain are implicit in the task. Before one can solve the parabolized stability equations for high-speed mixing layers, the mean flow must first be found. In the past, instability analysis of high-speed mixing layer has mostly been performed on mean flow profiles calculated by the boundary layer equations. In carrying out this project, it is believed that the boundary layer equations might not give an accurate enough nonparallel, nonlinear mean flow needed for parabolized stability analysis. A more accurate mean flow can, however, be found by solving the parabolized Navier-Stokes equations. The advantage of the parabolized Navier-Stokes equations is that its accuracy is consistent with the PSE method. Furthermore, the method of solution is similar. Hence, the major part of the effort of the work of this year has been devoted to the development of an explicit numerical marching scheme for the solution of the Parabolized Navier-Stokes equation as applied to the high-seed mixing layer problem.

Tam, Christopher↗
NASA NTRS2015

Higher-Order Methods for Compressible Turbulent Flows Using Entropy Variables

Turbulent flows have a large range of spatial and temporal scales which need to be resolved in order to obtain accurate predictions. Higher-order methods can provide greater efficiency for simulations requiring high spatial and temporal resolution, allowing for solutions with fewer degrees of freedom and lower computational cost than traditional second-order computational fluid dynamics (CFD) methods.1 Higher-order methods have been widely used for turbulent flows. However, the reduced numerical stabilization present in higher-order schemes implies that special care needs to be taken in the development of numerical methods to suppress nonlinear instabilities.2–6 In this work we present the development of a higher-order space-time discontinuous Galerkin method with a focus on the aspects of our numerical scheme required for ensuring nonlinear stability for turbulent simulations at high Reynolds numbers.

Diosady, Laslo T.↗
NASA NTRS1973

The use of the modified Cholesky decomposition in divergence and classification calculations

The use of the Cholesky decomposition technique is analyzed as applied to the feature selection and classification algorithms used in the analysis of remote sensing data (e.g. as in LARSYS). This technique is approximately 30% faster in classification and a factor of 2-3 faster in divergence, as compared with LARSYS. Also numerical stability and accuracy are slightly improved. Other methods necessary to deal with numerical stablity problems are briefly discussed.

Vanroony, D. L.↗
NASA NTRS1976

Numerical computation of transient coaxial entry tube flows

A numerical program was developed to compute transient laminar flows in two dimensions including multicomponent mixing and chemical reaction. The program can compute both incompressible flows and compressible flows at all speeds, and it is applied to describe transient and steady state solutions for low subsonic, coaxial entry, tue flows. Single component, nonreacting flows comprise most of the solutions, but one steady state solution is presented for trace concentration constituents engaging in a second order reaction. Numerical stability was obtained by adding at each calculation point a correction for numerical diffusion errors caused by truncation of the Taylor series used to finite difference the conservation equations. Transient computations were made for fluids initially at rest, then subjected to step velocity inputs that were uniform across each region of the entry plane and were held constant throughout the computation period. For center tube to annulus velocity ratios of 0.5 and 2.0, the bulk fluid in the tube initially moved in plug flow, but strong radial flows developed near the injection plane which moved the fluid into the high shear region between the jets and away from the tube wall.

Wieber, P. R.↗
NASA NTRS1972

On the numerical solution of two-dimensional, laminar compressible flows with imbedded shock waves.

The complete, time-dependent Navier-Stokes equations are expressed in conservation form and solved by employing an explicit finite difference numerical technique which incorporates artificial viscosity terms of the form first suggested by Rusanov for numerical stability in the vicinity of shock waves. Surface boundary conditions are developed in a consistent and unique manner through the use of a physically oriented extrapolation procedure. From numerical experimentation an extended range for the explicit stability parameter is established. Also employed is an additional convergence parameter which relates incremental spatial steps. Convergence of the transient solution to a steady state flow was obtained after 400 to 500 time steps.

Goodrich, W. D.↗
NASA NTRS1981

Stability analysis of numerical boundary conditions and implicit difference approximations for hyperbolic equations

Implicit, noniterative, finite difference schemes were recently developed by several authors for multidimensional systems of nonlinear hyperbolic partial differential equations. When applied to linear model equations with periodic boundary conditions those schemes are unconditionally stable (A-stable). As applied in practice the algorithms often face a severe time step restriction. A major source of the difficulty is the treatment of the numerical boundary conditions. One conjecture was that unconditional stability requires implicit numerical boundary conditions. An apparent counter example was the space time extrapolation considered by Gustafsson, Kreiss, and Sunstrom. Spatial (implicit) and space time (explicit) extrapolation using normal mode analysis for a finite and infinite number of spatial mesh intervals are examined. The results indicate that for unconditional stability with a finite number of spatial mesh intervals, the numerical boundary conditions must be implicit.

Beam, R. M.↗
NASA NTRS1974

Study report on guidelines and test procedures for investigating stability of nonlinear cardiovascular control system models

A general study of the stability of nonlinear as compared to linear control systems is presented. The analysis is general and, therefore, applies to other types of nonlinear biological control systems as well as the cardiovascular control system models. Both inherent and numerical stability are discussed for corresponding analytical and graphic methods and numerical methods.

Fitzjerrell, D. G.↗
NASA NTRS1972

Generation of the invariant coefficients of the characteristic polynomial for an nxn matrix

In theories of numerical stability, roots to a characteristic polynomial are sought, which, in the case of the predictor with iterative correction method of numerical integration, are eigenvalues of a matrix whose elements depend on the coefficients used in the integration process. The characteristic polynomial is displayed explicitly in terms of the elements of the characteristic matrix.

Beaudet, P. R.↗