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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
NASA NTRS1987

Stability/Instability Analysis of Rotating Machinery

Numerical index of stability calculated for nonlinear system. Technique for determining rotor stability or instability from analysis of measurements adapted for use with computer simulations of rotor motion. Involves calculation of log decrement or increment of vibration amplitude. Applicable to rotors mounted in loose bearings and to similar problems in which load-versus-deflection characteristics nonlinear. Developed for assessments of vibrational characteristics of turbopump rotors, technique also usable with such mechanisms as high-speed ball bearings.

Powers, Richard W.↗
NASA NTRS1987

The study of the physics of cometary nuclei

The numerical calculations of stability for many possible orbits of the double nucleus for P/Holmes showed that the likelihood of such a precollision history was quite high. A number of investigations were made of hypothetical orbits for particles about the asteroid Amphitrite to test for stability. The purpose was to establish more favorable fly-by orbits close to the asteroid for the Galileo missions en-route to Jupiter, reducing the collisional hazards. A statistical study was made of the orbits of long-period comets with small original semi-major axes recently perturbed from the great Opik-Oort Cloud. The results from the space missions to Halley's comet are partially reported in the two papers in the appendices.

Whipple, Fred L.↗
NASA NTRS1988

A delayed action oscillator for ENSO

A simple nonlinear model is proposed for the El Nino/Southern Oscillation (ENSO) phenomenon. Its key feature is the inclusion of oceanic wave transit effects through a negative, delayed feedback. A linear stability analysis and numerical results are presented, showing that the period of the oscillation is typically several times the delay. It is suggested that such an effect can account for the long time scale of ENSO.

Suarez, Max J.↗
NASA NTRS1989

Linear instability of supersonic plane wakes

In this paper we present a theoretical and numerical study of the growth of linear disturbances in the high-Reynolds-number and laminar compressible wake behind a flat plate which is aligned with a uniform stream. No ad hoc assumptions are made as to the nature of the undisturbed flow (in contrast to previous investigations) but instead the theory is developed rationally by use of proper wake-profiles which satisfy the steady equations of motion. The initial growth of near wake perturbation is governed by the compressible Rayleigh equation which is studied analytically for long- and short-waves. These solutions emphasize the asymptotic structures involved and provide a rational basis for a nonlinear development. The evolution of arbitrary wavelength perturbations is addressed numerically and spatial stability solutions are presented that account for the relative importance of the different physical mechanisms present, such as three-dimensionality, increasing Mach numbers enough (subsonic) Mach numbers, there exists a region of absolute instability very close to the trailing-edge with the majority of the wake being convectively unstable. At higher Mach numbers (but still not large-hypersonic) the absolute instability region seems to disappear and the maximum available growth-rates decrease considerably. Three-dimensional perturbations provide the highest spatial growth-rates.

Papageorgiou, D. T.↗
NASA NTRS1991

Dynamic dissipative compensator design for large space structures

Control system design is considered for attitude control and vibration suppression of flexible space structures. The problem addressed is that of controlling both the zero-frequency rigid-body modes and the elastic modes. Model-based compensators, which employ observers tuned to the plant parameters, are first investigated. Such compensators are shown to generally exhibit high sensitivity to the knowledge of the parameters, especially the elastic mode frequencies. To overcome this problem a class of dynamic dissipative compensators is next proposed, which robustly stabilize the plant in the presence of unmodeled dynamics and parametric uncertainties. An analytical proof of robust stability is given, and a method of implementing the controller as a strictly proper compensator is given. Methods of designing such controllers to obtain optimal performance and robust stability are presented. Numerical and experimental results of application of the methods are presented, which indicate that dynamic dissipative controllers can simultaneously provide excellent performance and robustness.

Joshi, S. M.↗
NASA NTRS1991

Linear and nonlinear dynamic analysis by boundary element method

An advanced implementation of the direct boundary element method (BEM) applicable to free-vibration, periodic (steady-state) vibration and linear and nonlinear transient dynamic problems involving two and three-dimensional isotropic solids of arbitrary shape is presented. Interior, exterior, and half-space problems can all be solved by the present formulation. For the free-vibration analysis, a new real variable BEM formulation is presented which solves the free-vibration problem in the form of algebraic equations (formed from the static kernels) and needs only surface discretization. In the area of time-domain transient analysis, the BEM is well suited because it gives an implicit formulation. Although the integral formulations are elegant, because of the complexity of the formulation it has never been implemented in exact form. In the present work, linear and nonlinear time domain transient analysis for three-dimensional solids has been implemented in a general and complete manner. The formulation and implementation of the nonlinear, transient, dynamic analysis presented here is the first ever in the field of boundary element analysis. Almost all the existing formulation of BEM in dynamics use the constant variation of the variables in space and time which is very unrealistic for engineering problems and, in some cases, it leads to unacceptably inaccurate results. In the present work, linear and quadratic isoparametric boundary elements are used for discretization of geometry and functional variations in space. In addition, higher order variations in time are used. These methods of analysis are applicable to piecewise-homogeneous materials, such that not only problems of the layered media and the soil-structure interaction can be analyzed but also a large problem can be solved by the usual sub-structuring technique. The analyses have been incorporated in a versatile, general-purpose computer program. Some numerical problems are solved and, through comparisons with available analytical and numerical results, the stability and high accuracy of these dynamic analysis techniques are established.

Ahmad, Shahid↗
NASA NTRS1995

Convergence acceleration of the Proteus computer code with multigrid methods

This report presents the results of a study to implement convergence acceleration techniques based on the multigrid concept in the two-dimensional and three-dimensional versions of the Proteus computer code. The first section presents a review of the relevant literature on the implementation of the multigrid methods in computer codes for compressible flow analysis. The next two sections present detailed stability analysis of numerical schemes for solving the Euler and Navier-Stokes equations, based on conventional von Neumann analysis and the bi-grid analysis, respectively. The next section presents details of the computational method used in the Proteus computer code. Finally, the multigrid implementation and applications to several two-dimensional and three-dimensional test problems are presented. The results of the present study show that the multigrid method always leads to a reduction in the number of iterations (or time steps) required for convergence. However, there is an overhead associated with the use of multigrid acceleration. The overhead is higher in 2-D problems than in 3-D problems, thus overall multigrid savings in CPU time are in general better in the latter. Savings of about 40-50 percent are typical in 3-D problems, but they are about 20-30 percent in large 2-D problems. The present multigrid method is applicable to steady-state problems and is therefore ineffective in problems with inherently unstable solutions.

Demuren, A. O.↗
NASA NTRS1995

Reliability enhancement of Navier-Stokes codes through convergence acceleration

Methods for enhancing the reliability of Navier-Stokes computer codes through improving convergence characteristics are presented. The improving of these characteristics decreases the likelihood of code unreliability and user interventions in a design environment. The problem referred to as a 'stiffness' in the governing equations for propulsion-related flowfields is investigated, particularly in regard to common sources of equation stiffness that lead to convergence degradation of CFD algorithms. Von Neumann stability theory is employed as a tool to study the convergence difficulties involved. Based on the stability results, improved algorithms are devised to ensure efficient convergence in different situations. A number of test cases are considered to confirm a correlation between stability theory and numerical convergence. The examples of turbulent and reacting flow are presented, and a generalized form of the preconditioning matrix is derived to handle these problems, i.e., the problems involving additional differential equations for describing the transport of turbulent kinetic energy, dissipation rate and chemical species. Algorithms for unsteady computations are considered. The extension of the preconditioning techniques and algorithms derived for Navier-Stokes computations to three-dimensional flow problems is discussed. New methods to accelerate the convergence of iterative schemes for the numerical integration of systems of partial differential equtions are developed, with a special emphasis on the acceleration of convergence on highly clustered grids.

Merkle, Charles L.↗
NASA NTRS1998

Dynamics of Liquids in Edges and Corners (DYLCO): IML-2 Experiment for the BDPU

Knowledge of the behavior of fluids possessing free surfaces is important to many fluid systems, particularly in space, where the normally subtle effects of surface wettability play a more dramatic and often surprising role. DYLCO for the IML-2 mission was proposed as a simple experiment to probe the particular behavior of capillary surfaces in containers of irregular cross section. Temperature control was utilized to vary the fluid-solid contact angle, a questionable thermodynamic parameter of the system, small changes in which can dramatically influence the configuration, stability, and flow of a capillary surface. Container shapes, test fluid, and temperature ranges were selected for observing both local changes in interface curvature as well as a global change in fluid orientation due to a critical wetting phenomenon. The experiment hardware performed beyond what was expected and fluid interfaces could be readily digitized post flight to show the dependence of the interface curvature on temperature. For each of the containers tested surfaces were observed which did not satisfy the classic equations for the prediction of interface shape with constant contact angle boundary condition. This is explained by the presence of contact angle hysteresis arising from expansion and contraction of the liquid during the heating and cooling steps of the test procedure. More importantly, surfaces exceeding the critical surface curvature required for critical wetting were measured, yet no wetting was observed. These findings are indeed curious and pose key questions concerning the role of hysteresis for this critical wetting phenomena. The stability of such surfaces was determined numerically and it is shown that stability is enhance (reduced) when a surface is in its 'advancing' ('receding') state, The analysis shows complete instability as the critical wetting condition is reached. The case of ideal dynamic wetting is addressed analytically in detail with results of significant flow characteristics presented in closed form. The solutions indicate a square root of T dependence of the capillary 'rise' rate which is corroborated by drop tower tests. The analysis clearly shows that infinite time is necessary for surfaces to reorient at the critical wetting transition.

Langbein, D.↗
NASA NTRS1999

Unified CFD Methods Via Flowfield-Dependent Variation Theory

This paper addresses the flowfield-dependent variation (FDV) methods in which complex physical phenomena are taken into account in the final form of partial differential equations to be solved so that finite difference methods (FDM) or finite element methods (FEM) themselves will not dictate the physics, but rather are no more than simply the options how to discretize between adjacent nodal points or within an element. The variation parameters introduced in the formulation are calculated from the current flowfield based on changes of Mach numbers, Reynolds numbers, Peclet numbers, and Damkohler numbers between adjacent nodal points, which play many significant roles such as adjusting the governing equations (hyperbolic, parabolic, and/or e!liptic), resolving various physical phenomena, and controlling the accuracy and stability of the numerical solution. The theory is verified by a number of example problems addressing the physical implications of the variation parameters which resemble the flowfield itself, shock capturing mechanism, transitions and interactions between inviscid/viscous, compressibility/incompressibility, and laminar/turbulent flows.

Chung, T. J.↗
NASA NTRS2004

GPC-Based Stable Reconfigurable Control

This paper presents development of multi-input multi-output (MIMO) Generalized Pre-dictive Control (GPC) law and its application to reconfigurable control design in the event of actuator saturation. A Controlled Auto-Regressive Integrating Moving Average (CARIMA) model is used to describe the plant dynamics. The control law is derived using input-output description of the system and is also related to the state-space form of the model. The stability of the GPC control law without reconfiguration is first established using Riccati-based approach and state-space formulation. A novel reconfiguration strategy is developed for the systems which have actuator redundancy and are faced with actuator saturation type failure. An elegant reconfigurable control design is presented with stability proof. Several numerical examples are presented to demonstrate the application of various results.

Soloway, Don↗
NASA NTRS2013

Comparative Study on High-Order Positivity-preserving WENO Schemes

The goal of this study is to compare the results obtained by non-positivity-preserving methods with the recently developed positivity-preserving schemes for representative test cases. In particular the more di cult 3D Noh and Sedov problems are considered. These test cases are chosen because of the negative pressure/density most often exhibited by standard high-order shock-capturing schemes. The simulation of a hypersonic nonequilibrium viscous shock tube that is related to the NASA Electric Arc Shock Tube (EAST) is also included. EAST is a high-temperature and high Mach number viscous nonequilibrium ow consisting of 13 species. In addition, as most common shock-capturing schemes have been developed for problems without source terms, when applied to problems with nonlinear and/or sti source terms these methods can result in spurious solutions, even when solving a conservative system of equations with a conservative scheme. This kind of behavior can be observed even for a scalar case (LeVeque & Yee 1990) as well as for the case consisting of two species and one reaction (Wang et al. 2012). For further information concerning this issue see (LeVeque & Yee 1990; Griffiths et al. 1992; Lafon & Yee 1996; Yee et al. 2012). This EAST example indicated that standard high-order shock-capturing methods exhibit instability of density/pressure in addition to grid-dependent discontinuity locations with insufficient grid points. The evaluation of these test cases is based on the stability of the numerical schemes together with the accuracy of the obtained solutions.

Kotov, Dmitry V.↗
NASA NTRS2014

Elastic Differential Cross Sections

The eikonal, partial wave (PW) Lippmann-Schwinger, and three-dimensional Lippmann-Schwinger (LS3D) methods are compared for nuclear reactions that are relevant for space radiation applications. Numerical convergence of the eikonal method is readily achieved when exact formulas of the optical potential are used for light nuclei (A less than or equal to 16) and the momentum-space optical potential is used for heavier nuclei. The PW solution method is known to be numerically unstable for systems that require a large number of partial waves, and, as a result, the LS3D method is employed. The effect of relativistic kinematics is studied with the PW and LS3D methods and is compared to eikonal results. It is recommended that the LS3D method be used for high energy nucleon- nucleus reactions and nucleus-nucleus reactions at all energies because of its rapid numerical convergence and stability.

Werneth, Charles M.↗
NASA NTRS2015

Nuclear Cross Sections for Space Radiation Applications

The eikonal, partial wave (PW) Lippmann-Schwinger, and three-dimensional Lippmann-Schwinger (LS3D) methods are compared for nuclear reactions that are relevant for space radiation applications. Numerical convergence of the eikonal method is readily achieved when exact formulas of the optical potential are used for light nuclei (A = 16) and the momentum-space optical potential is used for heavier nuclei. The PW solution method is known to be numerically unstable for systems that require a large number of partial waves, and, as a result, the LS3D method is employed. The effect of relativistic kinematics is studied with the PW and LS3D methods and is compared to eikonal results. It is recommended that the LS3D method be used for high energy nucleon-nucleus reactions and nucleus-nucleus reactions at all energies because of its rapid numerical convergence and stability for both non-relativistic and relativistic kinematics.

Werneth, C. M.↗
NASA NTRS2019

Pulsed Plasma Acceleration Modeling in Detonation and Deflagration Modes

Gas-fed electromagnetic pulsed plasma accelerators discharge electrical energy into a gas to ionize and electromagnetically accelerate propellant in the domain. Efforts to model pulsed accelerators have either assumed that the discharge is short and completely transient, accelerating the gas by entraining it in a moving current sheet, or that the discharge is relatively long, establishing a stable quasi-steady current distribution pattern that accelerates a plasma flowing through it. There have been pulsed plasma accelerator tests that appear to fit somewhere between these two bounds, exhibiting some properties that are associated with the purely transient devices while also showing others that are associated with quasi-steady-state plasma acceleration. A model is presented based upon the premise that all pulsed plasma accelerators first form an accelerating current sheet (detonation mode accelerator). Depending upon the pulse length and the gas conditions in the dis- charge channel, the plasma sheet may reach the end of the accelerator before the discharge has completed a full half-cycle, wherein the proposed model transitions to a quasi-steady description of the acceleration process (deflagration mode accelerator). A review of the entire model is presented, highlighting improvements and upgrades implemented to aid in the stability of the numerical scheme used to model the gas flow in the channel and an improved treatment of current sheet mass shedding, which adds gas to the wake of the sheet. The assumptions employed to model the transition from detonation to deflagration mode are presented and used to generate solutions to the governing equations. The modeling of the deflagration-mode under the present assumptions results in very low deflagration impulse bits relative to those obtained at the end of the detonation mode, implying that the present assumptions and modeling of the deflagration mode may not represent a fruitful approach.

Polzin, Kurt A.↗
NASA NTRS2019

Pulsed Plasma Acceleration Modeling in Detonation and Deflagration Modes

Gas-fed electromagnetic pulsed plasma accelerators discharge electrical energy into a gas to ionize and electromagnetically accelerate propellant in the domain. Efforts to model pulsed accelerators have either assumed that the discharge is short and completely tran- sient, accelerating the gas by entraining it in a moving current sheet, or that the discharge is relatively long, establishing a stable quasi-steady current distribution pattern that accel- erates a plasma flowing through it. There have been pulsed plasma accelerator tests that appear to fit somewhere between these two bounds, exhibiting some properties that are associated with the purely transient devices while also showing others that are associated with quasi-steady-state plasma acceleration. A model is presented based upon the premise that all pulsed plasma accelerators first form an accelerating current sheet (detonation mode accelerator). Depending upon the pulse length and the gas conditions in the dis- charge channel, the plasma sheet may reach the end of the accelerator before the discharge has completed a full half-cycle, wherein the proposed model transitions to a quasi-steady description of the acceleration process (deflagration mode accelerator). A review of the entire model is presented, highlighting improvements and upgrades implemented to aid in the stability of the numerical scheme used to model the gas flow in the channel and an im- proved treatment of current sheet mass shedding, which adds gas to the wake of the sheet. The assumptions employed to model the transition from detonation to deflagration mode are presented and used to generate solutions to the governing equations. The modeling of the deflagration-mode under the present assumptions results in very low deflagration im- pulse bits relative to those obtained at the end of the detonation mode, implying that the present assumptions and modeling of the deflagration mode may not represent a fruitful approach.

Polzin, Kurt A.↗
NASA NTRS1996

Multigrid Acceleration of Time-Accurate Navier-Stokes Calculations

A numerical scheme to solve the unsteady Navier-Stokes equations is described. The scheme is fully implicit in time and is unconditionally stable (at least for first- and second-order discretizations of the physical time derivatives). With unconditional stability, the choice of the time step is based on the physical phenomena to be resolved rather than limited by numerical stability. This is especially important for high Reynolds number viscous flows, where the spatial variation of grid cell size can be as much as six orders of magnitude. A multigrid-multiblock, steady-state, three-dimensional Navier-Stokes solver, TLNS3D, was modified to iteratively invert the equations at each physical time step. The implementation of this procedure in TLNS3D is discussed. The implications of applying several popular turbulence models to unsteady flow are also considered. Numerical results are presented to show the application of the scheme to various two-dimensional turbulent flows. The results of a three-dimensional laminar flow calculation are also given.

Melson, N. Duane↗
NASA NTRS1994

On the Dynamics of TVD Schemes

The dynamics of a class of TVD schemes for model hyperbolic and parabolic equations is studied numerically using a highly parallel supercomputer (CM-5). The objective is to utilize the highly parallel property of the CM-5 to reveal the reliable time step and entropy parameter ranges, and the degree of compressible flux limiters to avoid slow convergence and the production of nonphysical numerical solutions. We choose to study the nonlinear stability property of TVD schemes numerically since it is otherwise not amenable analytically.

Yee, H. C.↗