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 991 records · Page 55

All sources NASA NTRS DOE OSTI

NASA NTRS1976

A two-dimensional finite difference solution for the transient thermal behavior of tubular solar collector

A numerical procedure was established using the finite-difference technique in the determination of the time-varying temperature distribution of a tubular solar collector under changing solar radiancy and ambient temperature. Three types of spatial discretization processes were considered and compared for their accuracy of computations and for selection of the shortest computer time and cost. The stability criteria of this technique was analyzed in detail to give the critical time increment to ensure stable computations. The results of the numerical analysis were in good agreement with the analytical solution previously reported. The numerical method proved to be a powerful tool in the investigation of the collector sensitivity to two different flow patterns and several flow control mechanisms.

Lansing, F. L.↗

NASA NTRS1980

A two-dimensional finite-difference solution for the transient thermal behavior of a tubular solar collector

A numerical procedure was established using the finite-difference technique in the determination of the time-varying temperature distribution of a tubular solar collector under changing solar radiancy and ambient temperature. Three types of spatial discretization processes were considered and compared for their accuracy of computations and for selection of the shortest computer time and cost. The stability criteria of this technique were analyzed in detail to give the critical time increment to ensure stable computations. The results of the numerical analysis were in good agreement with the analytical solution previously reported. The numerical method proved to be a powerful tool in the investigation of the collector sensitivity to two different flow patterns and several flow control mechanisms.

Lansing, F. L.↗

NASA NTRS1991

The asymptotic spectra of banded Toeplitz and quasi-Toeplitz matrices

Toeplitz matrices occur in many mathematical, as well as, scientific and engineering investigations. This paper considers the spectra of banded Toeplitz and quasi-Toeplitz matrices with emphasis on non-normal matrices of arbitrarily large order and relatively small bandwidth. These are the type of matrices that appear in the investigation of stability and convergence of difference approximations to partial differential equations. Quasi-Toeplitz matrices are the result of non-Dirichlet boundary conditions for the difference approximations. The eigenvalue problem for a banded Toeplitz or quasi-Toeplitz matrix of large order is, in general, analytically intractable and (for non-normal matrices) numerically unreliable. An asymptotic (matrix order approaches infinity) approach partitions the eigenvalue analysis of a quasi-Toeplitz matrix into two parts, namely the analysis for the boundary condition independent spectrum and the analysis for the boundary condition dependent spectrum. The boundary condition independent spectrum is the same as the pure Toeplitz matrix spectrum. Algorithms for computing both parts of the spectrum are presented. Examples are used to demonstrate the utility of the algorithms, to present some interesting spectra, and to point out some of the numerical difficulties encountered when conventional matrix eigenvalue routines are employed for non-normal matrices of large order. The analysis for the Toeplitz spectrum also leads to a diagonal similarity transformation that improves conventional numerical eigenvalue computations. Finally, the algorithm for the asymptotic spectrum is extended to the Toeplitz generalized eigenvalue problem which occurs, for example, in the stability of Pade type difference approximations to differential equations.

Beam, Richard M.↗

NASA NTRS1989

Numerical algorithms for computations of feedback laws arising in control of flexible systems

Several continuous models will be examined, which describe flexible structures with boundary or point control/observation. Issues related to the computation of feedback laws are examined (particularly stabilizing feedbacks) with sensors and actuators located either on the boundary or at specific point locations of the structure. One of the main difficulties is due to the great sensitivity of the system (hyperbolic systems with unbounded control actions), with respect to perturbations caused either by uncertainty of the model or by the errors introduced in implementing numerical algorithms. Thus, special care must be taken in the choice of the appropriate numerical schemes which eventually lead to implementable finite dimensional solutions. Finite dimensional algorithms are constructed on a basis of a priority analysis of the properties of the original, continuous (infinite diversional) systems with the following criteria in mind: (1) convergence and stability of the algorithms and (2) robustness (reasonable insensitivity with respect to the unknown parameters of the systems). Examples with mixed finite element methods and spectral methods are provided.

Lasiecka, Irena↗

NASA NTRS1987

Nonparallel stability of boundary layers

The asymptotic formulations of the nonparallel linear stability of incompressible growing boundary layers are critically reviewed. These formulations can be divided into two approaches. The first approach combines a numerical method with either the method of multiple scales, or the method of averaging, of the Wentzel-Kramers-Brillouin (WKB) approximation; all these methods yield the same result. The second approach combined a multi-structure theory with the method of multiple scales. The first approach yields results that are in excellent agreement with all available experimental data, including the growth rates as well as the neutral stability curve. The derivation of the linear stability of the incompressible growing boundary layers is explained.

Nayfeh, Ali H.↗

NASA NTRS2003

Functional Coordination of a Full-Body Gaze Control Mechanisms Elicited During Locomotion

Control of locomotion requires precise interaction between several sensorimotor subsystems. Exposure to the microgravity environment of spaceflight leads to postflight adaptive alterations in these multiple subsystems leading to postural and gait disturbances. Countermeasures designed to mitigate these postflight gait alterations will need to be assessed with a new generation of functional tests that evaluate the interaction of various elements central to locomotor control. The goal of this study is to determine how the multiple, interdependent, full- body sensorimotor subsystems aiding gaze stabilization during locomotion are functionally coordinated. To explore this question two experiments were performed. In the first study (Study 1) we investigated how alteration in gaze tasking changes full-body locomotor control strategies. Subjects (n=9) performed two discreet gaze stabilization tasks while walking at 6.4 km/hr on a motorized treadmill: 1) focusing on a central point target; 2) reading numeral characters; both presented at 2m in front at eye level. The second study (Study 2) investigated the potential of adaptive remodeling of the full-body gaze control systems following exposure to visual-vestibular conflict. Subjects (n=14) walked (6.4 km/h) on the treadmill before and after they were exposed to 0.5X minifying lenses worn for 30 minutes during self-generated sinusoidal vertical head rotations performed while seated. In both studies we measured: temporal parameters of gait, full body sagittal plane segmental kinematics of the head, trunk, thigh, shank and foot, accelerations along the vertical axis at the head and the shank, and the vertical forces acting on the support surface. Results from Study 1 showed that while reading numeral characters as compared to the central point target: 1) compensatory head pitch movements were on average 22% greater 2) the peak acceleration measured at the head was significantly reduced by an average of 13% in four of the six subjects 3) the knee joint total movement was on average 11% greater during the period from the heel strike event to the peak knee flexion event in stance phase of the gait cycle. Results from Study 2 indicate that following exposure to visual-vestibular conflict changes in full-body strategies were observed consistent with the requirement to aid gaze stabilization during locomotion.

Bloomberg, Jacob J.↗

NASA NTRS1992

The effect of initial imperfections on shell stability: An updated review

The development of the Delft Interactive Shell Design Code (DISDECO) is described. The purpose is to make the accumulated theoretical, numerical, and practical knowledge of the last years readily accessible to users interested in the analysis of buckling sensitive structures. With this open ended, hierarchical, interactive computer code, the user can access workstation programs of increasing complexity. As a first update of DISDECO, the computational modules of the initial level were modified so as to handle combined loads consisting of axial compression, internal and external pressure, and clockwise or counter clockwise torque. The applied load is assumed to have a uniform spatial distribution and is divided into a fixed and a variable part. The magnitude of the variable part is allowed to vary in proportion to a load parameter Lambda. This leads to an eigenvalue problem for the critical lambda sub C.

Arbocz, J.↗

NASA NTRS2001

Driving Parameters for Distributed and Centralized Air Transportation Architectures

This report considers the problem of intersecting aircraft flows under decentralized conflict avoidance rules. Using an Eulerian standpoint (aircraft flow through a fixed control volume), new air traffic control models and scenarios are defined that enable the study of long-term airspace stability problems. Considering a class of two intersecting aircraft flows, it is shown that airspace stability, defined both in terms of safety and performance, is preserved under decentralized conflict resolution algorithms. Performance bounds are derived for the aircraft flow problem under different maneuver models. Besides analytical approaches, numerical examples are presented to test the theoretical results, as well as to generate some insight about the structure of the traffic flow after resolution. Considering more than two intersecting aircraft flows, simulations indicate that flow stability may not be guaranteed under simple conflict avoidance rules. Finally, a comparison is made with centralized strategies to conflict resolution.

Feron, Eric↗

NASA NTRS1974

The three dimensional motion and stability of a rotating space station-cable-counterweight configuration

The three dimensional equations of motion for a cable connected space station-counterweight system are developed using a Lagrangian formulation. The system model employed allows for cable and end body damping and restoring effects. The equations are then linearized about the equilibrium motion and nondimensionalized. To first degree, the out-of-plane equations uncouple from the in-plane equations. From the general in-plane characteristic equation, necessary conditions for stability are obtained. The Routh-Hurwitz necessary and sufficient conditions for stability are derived for the general out-of-plane characteristic equation. Special cases of the in-plane and out-of-plane equations are then examined for stability. Time constants for the least damped mode are obtained for a range of system parameters by numerical examination of the roots of the in-plane and out-of-plane characteristic polynomials.

Bainum, P. M.↗

NASA NTRS1976

On the computation and updating of the modified Cholesky decomposition of a covariance matrix

Methods for obtaining and updating the modified Cholesky decomposition (MCD) for the particular case of a covariance matrix when one is given only the original data are described. These methods are the standard method of forming the covariance matrix K then solving for the MCD, L and D (where K=LDLT); a method based on Householder reflections; and lastly, a method employing the composite-t algorithm. For many cases in the analysis of remotely sensed data, the composite-t method is the superior method despite the fact that it is the slowest one, since (1) the relative amount of time computing MCD's is often quite small, (2) the stability properties of it are the best of the three, and (3) it affords an efficient and numerically stable procedure for updating the MCD. The properties of these methods are discussed and FORTRAN programs implementing these algorithms are listed.

Vanrooy, D. L.↗

NASA NTRS1975

Effect of halo component on bar-formation in disk galaxies

Numerical experiments are performed to determine the effect of a fixed halo component of mass on the stability of purely stellar disks. The rotation curve of the fixed halo component corresponds to the Schmidt model of the galaxy. It is found that when the stellar disk contains less than 50% of the total mass, the large-scale bar-making instability is effectively suppressed. For disks containing 50% or more of the total mass, a bar structure quickly forms.

Hohl, F.↗

NASA NTRS1978

Dynamics of spinning spacecraft during the deployment of flexible boom

The effect of flexibility on the dynamics of a spinning symmetrical spacecraft during the extension (or retraction) of boom-type appendages is treated by modeling the appendages as compound spherical pendula with varying lengths. Both boom stiffness and structural damping are included. For constant boom length the resulting linearized equations of motion contain periodic coefficients involving the spin frequency. A bounded transformation converts this system into a kinematically equivalent one (with only constant coefficients), whose stability is analyzed using the Kelvin-Tait-Chetaev theorem. The dynamics during extension is predicted by numerical simulation of the nonlinear equations.

James, P. K.↗

NASA NTRS1981

Stable boundary approximations for a class of implicit schemes for the one-dimensional inviscid equations of gas dynamics

The applicability to practical calculations of recent theoretical developments in the stability analysis of difference approximations for initial-boundary-value problems of the hyperbolic type. For the numerical experiments, select the one-dimensional inviscid gas-dynamic equations in conservation-law form is selected. A class of implicit schemes based on linear multistep methods for ordinary differential equations is chosen and the use of space or space-time extrapolations as implicit or explicit boundary schemes is emphasized. Some numerical examples with various inflow-outflow conditions highlight the commonly discussed issues: explicit versus implicit boundary schemes, unconditionally stable schemes, and underspecification or overspecification of boundary conditions.

Yee, H. C.↗

NASA NTRS1982

On current aspects of finite element computational fluid mechanics for turbulent flows

A set of nonlinear partial differential equations suitable for the description of a class of turbulent three-dimensional flow fields in select geometries is identified. On the basis of the concept of enforcing a penalty constraint to ensure accurate accounting of ordering effects, a finite element numerical solution algorithm is established for the equation set and the theoretical aspects of accuracy, convergence and stability are identified and quantized. Hypermatrix constructions are used to formulate the reduction of the computational aspects of the theory to practice. The robustness of the algorithm, and the computer program embodiment, have been verified for pertinent flow configurations.

Baker, A. J.↗

NASA NTRS1983

Quasi-axisymmetric circulation and superrotation in planetary atmospheres

Atmospheric superrotation is ubiquitous in the solar system, and numerous theoretical models have been proposed to describe the effect. However, no model is fully accepted. Problems are primarily related to an understanding of the angular momentum and the energy budgets. In connection with the present investigation, attention is given to hybrid models, where three-dimensional models may describe eddy processes and two-dimensional symmetric models may describe the global scale circulation in the framework of a viscous medium. It is believed that the concept of an axisymmetric circulation is still useful. It is the main purpose of the current investigation to explore some of the properties of the circulation in planetary atmospheres. It is assumed that the eddy transports are downgradient or diffusive in nature. Energy and momentum are perceived to cascade from the global scale circulation driven by the sun to the smaller scale eddies which in turn determine the diffusive properties. A numerical model of the Venusian circulation is presented, and results are discussed for different conditions of static stability.

Mayr, H. G.↗

NASA NTRS1987

Chemistry-split techniques for viscous reactive blunt body flow computations

The weak-coupling structure between the fluid and species equations has been exploited and resulted in three, closely related, time-iterative implicit techniques. While the primitive variables are solved in two separated groups and each by an Alternating Direction Implicit (ADI) factorization scheme, the rate-species Jacobian can be treated in either full or diagonal matrix form, or simply ignored. The latter two versions render the split technique to solving for species as scalar rather than vector variables. The solution is completed at the end of each iteration after determining temperature and pressure from the flow density, energy and species concentrations. Numerical experimentation has shown that the split scalar technique, using partial rate Jacobian, yields the best overall stability and consistency. Satisfactory viscous solutions were obtained for an ellipsoidal body of axis ratio 3:1 at Mach 35 and an angle of attack of 20 degrees.

Li, C. P.↗

NASA NTRS1987

Mutual interactions of redox couples via electron exchange in silicate melts - Models for geochemical melt systems

The mutual interactions via electron exchange of redox couples in glass-forming melts were investigated both theoretically and experimentally. A thermodynamic approach for considering the mutual interactions leads to conclusion that the degree of mutual interaction in the melt should be proportional in part to the difference in relative reduction potentials of the interacting redox couples. Experimental studies verify this conclusion for numerous redox couples in several composition/temperature/oxygen fugacity regimes. Geochemical systems simultaneously possess many potentially multivalent elements; the stabilized redox states in the resulting magmas can be explained in part by mutual interactions and by redox buffering through the central Fe(III)- Fe(II) couples in the melts. The significance of these results for basaltic magmas of the earth, moon, and meteorites is addressed.

Schreiber, Henry D.↗

NASA NTRS1987

Resistive-heating instability of a compressible, force-free plasma

The thermal stability of a resistive, force-free plasma in a two-dimensional, slab geometry is investigated numerically. Linearly unstable modes that arise as an energetic consequence of Ohmic heating are examined in relation to the dynamic magnetic-tearing instability. Plasma compression is found to give rise to a set of unstable thermal modes not seen in previous incompressible studies. Such modes can occur in parameter regimes where no purely growing tearing mode exists.

Sparks, L.↗