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 1,027 records · Page 57

All sources NASA NTRS DOE OSTI

NASA NTRS1985

Symposium on Numerical and Physical Aspects of Aerodynamic Flows, 3rd, California State University, Long Beach, CA, January 21-24, 1985, Proceedings

Among the topics discussed are hydrodynamic stability and turbulent transition, separated flow simulations employing finite difference methods, a quasi-simultaneous finite difference approach for strongly interacting flow, the Newton solution of coupled Euler and boundary layer equations, multielement airfoil aerodynamics, vorticity dynamics, the significance of the thin layer Navier-Stokes approximation, a comparison of interactive and Navier-Stokes calculations of separating boundary layer flows, and a Navier-Stokes analysis for the mechanism of shock wave/boundary layer interactions. Also discussed are transonic interactions of unsteady vortical flows, boundary layers on oscillating airfoils, a laser-interferometer skin friction meter, flow over surface-mounted semibluff bodies, computations for spinning and spin-stabilized supersonic projectiles, transition calculations in three-dimensional flows, the effect of projectile boattailing in transonic flow, flow characteristics of a body of revolution at incidence, the influence of the boundary layer on the leading edge flow of a swept wing, and an interactive scheme for three-dimensional transonic flows.

Source record↗

NASA NTRS1985

Effects of stratospheric lapse rate on numerically simulated thunderstorm cloud top structure

The three-dimensional anelastic model of Schlesinger (1975, 1978, 1984) is utilized to study the effects of stratospheric stability on thunderstorm cloud top structure. Three experiments were conducted; in the first the stratosphere is assumed isothermal, in the second test there is a temperature lapse of 3 K/km, and in the third case there is a 3 K/km inversion. Time variations of vertical velocity, mature storm structure, and cloud top are examined. It is observed that the cloud summit height is lowered as stratospheric stability increases, a cloud-top thermal couplet occurs in all cases, the cold region has a U or V shape, and there is a height minimum in the warm regions.

Schlesinger, Robert E.↗

NASA NTRS1989

Simulation of collisional transport processes and the stability of planetary rings

The utility of the phase-space fluid method for the study of planetary ring dynamics is presently demonstrated through the numerical solution of a model kinetic equation for a flattened Keplerian disk. Attention is given to ringlets composed of single-sized particles, as well as to ringlets composed of two different-sized particles; in the latter case, the ringlets evolve in such a way that the lighter particles are confined by the heavier ones. The results obtained indicate that some natural process may sharpen the optical depth profile of edges even without an external forcing mechanism, and that intermediate optical depths are dynamically preferred in some cases.

Brophy, Thomas G.↗

NASA NTRS1990

WKBJ approximation of the stability of a frontal mean state

The baroclinic instability of a frontal mean state is investigated using the WKBJ approximation. The results are compared with numerical calculations performed on the same mean state. Excellent agreement (within 5 percent) is found for jets whose half-width is as a Rossby radius of deformation. For jets 20 percent broader, the agreement is almost perfect.

Lindzen, Richard S.↗

NASA NTRS2003

Structure and Stability of One-Dimensional Detonations in Ethylene-Air Mixtures

The propagation of one-dimensional detonations in ethylene-air mixtures is investigated numerically by solving the one-dimensional Euler equations with detailed finite-rate chemistry. The numerical method is based on a second-order spatially accurate total-variation-diminishing scheme and a point implicit, first-order-accurate, time marching algorithm. The ethylene-air combustion is modeled with a 20-species, 36-step reaction mechanism. A multi-level, dynamically adaptive grid is utilized, in order to resolve the structure of the detonation. Parametric studies over an equivalence ratio range of 0.5 less than phi less than 3 for different initial pressures and degrees of detonation overdrive demonstrate that the detonation is unstable for low degrees of overdrive, but the dynamics of wave propagation varies with fuel-air equivalence ratio. For equivalence ratios less than approximately 1.2 the detonation exhibits a short-period oscillatory mode, characterized by high-frequency, low-amplitude waves. Richer mixtures (phi greater than 1.2) exhibit a low-frequency mode that includes large fluctuations in the detonation wave speed; that is, a galloping propagation mode is established. At high degrees of overdrive, stable detonation wave propagation is obtained. A modified McVey-Toong short-period wave-interaction theory is in excellent agreement with the numerical simulations.

Yungster, S.↗

NASA NTRS1983

Numerical simulation of two-dimensional heat transfer in composite bodies with application to de-icing of aircraft components

Transient, numerical simulations of the de-icing of composite aircraft components by electrothermal heating were performed for a two dimensional rectangular geometry. The implicit Crank-Nicolson formulation was used to insure stability of the finite-difference heat conduction equations and the phase change in the ice layer was simulated using the Enthalpy method. The Gauss-Seidel point iterative method was used to solve the system of difference equations. Numerical solutions illustrating de-icer performance for various composite aircraft structures and environmental conditions are presented. Comparisons are made with previous studies. The simulation can also be used to solve a variety of other heat conduction problems involving composite bodies.

Chao, D. F. K.↗

NASA NTRS1990

A computer module used to calculate the horizontal control surface size of a conceptual aircraft design

The creation of a computer module used to calculate the size of the horizontal control surfaces of a conceptual aircraft design is discussed. The control surface size is determined by first calculating the size needed to rotate the aircraft during takeoff, and, second, by determining if the calculated size is large enough to maintain stability of the aircraft throughout any specified mission. The tail size needed to rotate during takeoff is calculated from a summation of forces about the main landing gear of the aircraft. The stability of the aircraft is determined from a summation of forces about the center of gravity during different phases of the aircraft's flight. Included in the horizontal control surface analysis are: downwash effects on an aft tail, upwash effects on a forward canard, and effects due to flight in close proximity to the ground. Comparisons of production aircraft with numerical models show good accuracy for control surface sizing. A modified canard design verified the accuracy of the module for canard configurations. Added to this stability and control module is a subroutine that determines one of the three design variables, for a stable vectored thrust aircraft. These include forward thrust nozzle position, aft thrust nozzle angle, and forward thrust split.

Sandlin, Doral R.↗

NASA NTRS2003

Identifying Head-Trunk and Lower Limb Contributions to Gaze Stabilization During Locomotion

The goal of the present study was to determine how the multiple, interdependent full-body sensorimotor subsystems respond to a change in gaze stabilization task constraints during locomotion. Nine subjects performed two 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 the level of their eyes. While subjects performed the tasks 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. We tested the hypothesis that with the increased demands placed on visual acuity during the number recognition task, subjects would modify full-body segmental kinematics in order to reduce perturbations to the head in order to successfully perform the task. We found that while reading numeral characters as - compared to the central point target: 1) compensatory head pitch movement was on average 22% greater despite the fact that the trunk pitch and trunk vertical translation movement control were not significantly changed; 2) coordination patterns between head and trunk as reflected by the peak cross correlation between the head pitch and trunk pitch motion as well as the peak cross correlation between the head pitch and vertical trunk translation motion were not significantly changed; 3) 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; 4) peak acceleration measured at the head was significantly reduced by an average of 13% in four of the six subjects. This was so even when the peak acceleration at the shank and the transmissibility of the shock wave at heel strike (measured by the peak acceleration ratio of the head/shank) remained unchanged. Taken together these results provide further evidence that the full body contributes to gaze stabilization during locomotion, and that its different functional elements can be modified online to contribute to gaze stabilization for different visual task constraints.

Mulavara, Ajitkumar P.↗

NASA NTRS1973

Survey of computational methods for three-dimensional supersonic inviscid flows with shocks

Hyperbolic systems of partial differential equations governing supersonic inviscid flows are discussed and analyzed. Finite-difference analogues for integrating these systems in the interior of fluid domains are described from two points of view: a differential form approach and an integral form approach. The algorithms presented are analyzed for stability and accuracy. The concept of time splitting is discussed and supplied to these methods to achieve increased numerical efficiency. Techniques for treating conditions at the boundaries of the fluid domain and shock-wave discontinuities at surfaces within the domain are described.

Maccormack, R. W.↗

NASA NTRS1983

Viscous, resistive MHD stability computed by spectral techniques

Expansions in Chebyshev polynomials are used to study the linear stability of one dimensional magnetohydrodynamic (MHD) quasi-equilibria, in the presence of finite resistivity and viscosity. The method is modeled on the one used by Orszag in accurate computation of solutions of the Orr-Sommerfeld equation. Two Reynolds like numbers involving Alfven speeds, length scales, kinematic viscosity, and magnetic diffusivity govern the stability boundaries, which are determined by the geometric mean of the two Reynolds like numbers. Marginal stability curves, growth rates versus Reynolds like numbers, and growth rates versus parallel wave numbers are exhibited. A numerical result which appears general is that instability was found to be associated with inflection points in the current profile, though no general analytical proof has emerged. It is possible that nonlinear subcritical three dimensional instabilities may exist, similar to those in Poiseuille and Couette flow.

Dahlburg, R. B.↗

NASA NTRS1983

Heavy-lift airship dynamics

The basic aerodynamic and dynamic properties of an example heavy-lift airship (HLA) configuration are analyzed using a nonlinear, multibody, 6-degrees-of-freedom digital simulation. The slung-payload model is described, and a preliminary analysis of the coupled vehicle-payload dynamics is presented. Trim calculations show the importance of control mixing selection and suggest performance deficiencies in crosswind stationkeeping for the unloaded example HLA. Numerically linearized dynamics of the unloaded vehicle exhibit a divergent yaw mode and an oscillatory pitch mode whose stability characteristic is sensitive to flight speed. An analysis of the vehicle-payload dynamics shows significant coupling of the payload dynamics with those of the basic HLA. It is shown that significant improvement in the vehicle's dynamic behavior can be achieved with the incorporation of a simple flight controller having proportional, rate, and integral-error feedbacks.

Tischler, M. B.↗

NASA NTRS1983

Aeroelastic stability of forward swept composite winged aircraft

This paper reviews the author's past and present aeroelastic stability and performance studies related to forward swept, composite wing aircraft. The influence of laminate elastic bend/twist coupling upon wing divergence, lateral control, and lift effectiveness will be illustrated by means of closed-form solutions, numerical analysis and simple wind-tunnel experiments. In addition, results of analyses of a freely flying flexible FSW aircraft are discussed to indicate the possible effects of the flexible forward swept wing on aircraft dynamic stability. These studies show, both theoretically and experimentally, that, if the aircraft is not carefully designed, a phenomenon referred to as body freedom flutter may appear.

Weisshaar, T. A.↗

NASA NTRS1984

Optimum design of rotor blades for vibration reduction in forward flight

Modern structural optimization techniques are applied to vibration reduction of helicopter rotor blades in forward flight. The objective function minimized consists of the oscillatory vertical hub shears or the hub rolling moments at one particular advance ratio. The behavior constraints are the frequency placements of the blade and the requirement that aeroelastic stability margins, in hover, remain unaffected by the optimization process. The aeroelastic stability and response analysis is based on a fully coupled flap-lag-torsional analysis of the blade. Numerical results are presented for some typical soft-in-plane hingeless rotor configurations indicating a 15-40 percent reduction in vibration levels, as well as a blade which is 20 percent lighter than the initial design. These results imply that structural optimization techniques can yield substantial practical benefits in the design process of rotor systems.

Friedmann, P. P.↗

NASA NTRS1985

Stability of limit cycles in frictionally damped and aerodynamically unstable rotor stages

This paper deals with the stability of limit cycles (Steady-State Oscillations) associated with the multi-degree-of-freedom model of a frictionally damped and aerodynamically unstable rotor stage. By using the first order averaging technique, a generalized criterion has been established to sort out those unstable limit cycles which govern the maximum transient amplitude beyond which the rotor stage becomes unstable.The stability of the remaining steady-state solutions is analyzed by linearizing the averaged system of differential equations. Numerical results are discussed for three-, four- and five-bladed disks.

Sinha, A.↗

NASA NTRS1986

Motion of an electron bunch through a plasma

In space research, pulse-modulated beams, or bunches, of energetic electrons have been used for tracing magnetic field lines and for remote sounding of magnetospheric electric fields. Doubts have been expressed regarding the possibility that an electron bunch would be able to travel through the plasma in space. However, experiments in the ionosphere have shown that, at least on some occasions, bunches of electrons with energies of a few keV have traveled over long distances with no apparent energy loss. The present study is concerned with the theoretical explanation for this phenomenon. The method employed in this study makes use of a direct numerical integration of the Vlasov equation. The results obtained provide the basis for a provisional conclusion concerning the stability of a moving electron bunch in the considered case. It appears that the system is mildly unstable.

Shoucri, M. M.↗

NASA NTRS1989

Numerical Study of Three-dimensional Spatial Instability of a Supersonic Flat Plate Boundary Layer

The behavior of spatially growing three-dimensional waves in a supersonic boundary layer was studied numerically by solving the complete Navier-Stokes equations. Satisfactory comparison with linear parallel and non-parallel stability theories, and experiment are obtained when a small amplitude inflow disturbance is used. The three-dimensional unsteady Navier-Stokes equations are solved by a finite difference method which is fourth-order and second-order accurate in the convection and viscous terms respectively, and second-order accurate in time. Spanwise periodicity is assumed. The inflow disturbance is composed of eigenfunctions from linear stability theory. By increasing the amplitude of the inflow disturbance, nonlinear effects in the form of a relaxation type oscillation of the time signal of rho(u) are observed.

Maestrello, Lucio↗

NASA NTRS1990

H2, fixed architecture, control design for large scale systems

The H2, fixed architecture, control problem is a classic linear quadratic Gaussian (LQG) problem whose solution is constrained to be a linear time invariant compensator with a decentralized processing structure. The compensator can be made of p independent subcontrollers, each of which has a fixed order and connects selected sensors to selected actuators. The H2, fixed architecture, control problem allows the design of simplified feedback systems needed to control large scale systems. Its solution becomes more complicated, however, as more constraints are introduced. This work derives the necessary conditions for optimality for the problem and studies their properties. It is found that the filter and control problems couple when the architecture constraints are introduced, and that the different subcontrollers must be coordinated in order to achieve global system performance. The problem requires the simultaneous solution of highly coupled matrix equations. The use of homotopy is investigated as a numerical tool, and its convergence properties studied. It is found that the general constrained problem may have multiple stabilizing solutions, and that these solutions may be local minima or saddle points for the quadratic cost. The nature of the solution is not invariant when the parameters of the system are changed. Bifurcations occur, and a solution may continuously transform into a nonstabilizing compensator. Using a modified homotopy procedure, fixed architecture compensators are derived for models of large flexible structures to help understand the properties of the constrained solutions and compare them to the corresponding unconstrained ones.

Mercadal, Mathieu↗

NASA NTRS1991

On the effects of viscosity on the stability of a trailing-line vortex

The linear stability of the Batchelor (1964) vortex is investigated. Particular emphasis is placed on modes found recently in a numerical study by Khorrami (1991). These modes have a number of features very distinct from those found previously for this vortex, including exhibiting small growth rates at large Reynolds numbers and susceptibility to destabilization by viscosity. These modes are described using asymptotic techniques, producing results which compare favorably with fully numerical results at large Reynolds numbers.

Duck, Peter W.↗