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 811 records · Page 45
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

Counterrotating prop-fan simulations which feature a relative-motion multiblock grid decomposition enabling arbitrary time-steps

Improvements are presented of a computer algorithm developed for the time-accurate flow analysis of rotating machines. The flow model is a finite volume method utilizing a high-resolution approximate Riemann solver for interface flux definitions. The numerical scheme is a block LU implicit iterative-refinement method which possesses apparent unconditional stability. Multiblock composite gridding is used to orderly partition the field into a specified arrangement of blocks exhibiting varying degrees of similarity. Block-block relative motion is achieved using local grid distortion to reduce grid skewness and accommodate arbitrary time step selection. A general high-order numerical scheme is applied to satisfy the geometric conservation law. An even-blade-count counterrotating unducted fan configuration is chosen for a computational study comparing solutions resulting from altering parameters such as time step size and iteration count. The solutions are compared with measured data.

Janus, J. Mark↗
NASA NTRS2000

Improving Tiltrotor Whirl-Mode Stability with Rotor Design Variations

Further increases in tiltrotor speeds are limited by coupled wing/rotor whirl-mode aeroelastic instability. Increased power, thrust, and rotor efficiency are not enough: the whirl-mode stability boundary must also be improved. With current technology, very stiff, thick wings of limited aspect ratio are essential to meet the stability requirements, which severely limits cruise efficiency and maximum speed. Larger and more efficient tiltrotors will need longer and lighter wings, for which whirl-mode flutter is a serious design issue. Numerous approaches to improving the whirl-mode airspeed boundary have been investigated, including tailored stiffness wings, active stability augmentation, variable geometry rotors, highly swept tips, and at one extreme, folding rotors. The research reported herein began with the much simpler approach of adjusting the chordwise positions of the rotor blade aerodynamic center and center of gravity, effected by offsetting the airfoil quarter chord or structural mass with respect to the elastic axis. The research was recently extended to include variations in blade sweep, control system stiffness, and pitch-flap coupling (delta(sub 3)). As an introduction to the subject, and to establish a baseline against which to measure stability improvements, this report will first summarize results. The paper will then discuss more advanced studies of swept blades and control-system modifications.

Acree, C. W., Jr.↗
NASA NTRS1994

Anisotropic Alfven-ballooning modes in Earth's magnetosphere

We have carried out a theoretical analysis of the stability and parallel structure of coupled shear Alfven and slow magnetosonic waves in Earth's inner magnetopause (i.e., at equatorial distances between about five and ten Earth radii) including effects of finite anisotropic Grad-Shafranov equation yields an approximate self-consistent magnetohydrodynamic (MHD) equilibrium. This MHD equilibrium is used in the numerical solution of a set of eigenmode equations which describe the field line eigenfrequency, linear stability, and parallel eigenmode structure. We call these modes anisotropic Alfven-ballooning modes. The main results are: (1) The field line eigenfrequency can be significantly lowered by finite pressure effects. (2) The parallel mode structure of the transverse wave components is fairly insensitive to changes in the plasma pressure, but the compressional magnetic component can become highly peaked near the magnetic equator as a result of increased pressure, especially when P(sub perpendicular to) is greater than P(sub parallel) (here P(sub perpendicular to) and P(sub parallel) are the perpendicular and parallel plasma pressure). (3) For the isotropic (P(sub parallel) = P(sub perpendicular to) = P) case ballooning instability can occur when the ratio of the plasma presure to the magnetic pressure, beta = P/(B squared/8 pi), exceeds a critical value beta(sup B)(sub 0) is approximately equal to 3.5 at the equator. (4) Compared to the isotropic case the critical beta value is lowered by anisotropy, either due to decreased field line bending stabilization when P(sub parallel) is greater than P(sub perpendicular to) or due to increased ballooning-mirror destabilization when P(sub perpendicular to) is greater than P(sub parallel). (5) We use a beta-delta stability diagram to display the regions of instability with respect to the equatorial values of the parameters bar beta and delta, where bar beta = (1/3)(beta(sub parallel) + 2 beta(sub perpendicular to)) is an average beta value and delta = 1 - P(sub parallel)/P(sub perpendicular to) is a measure of the plasma anisotropy. The diagram is divided into regions corresponding to the firehose, mirror and ballooning instabilities. It appears that observed values of the plasma pressure are below the critical value for the isotropic ballooning instability but it may be possible to approach a ballooning-mirror instability when P(sub perpendicular to)/P(sub parallel) is greater than or approximately 2.

Chan, Anthony A.↗
NASA NTRS2012

Receptivity of Hypersonic Boundary Layers to Distributed Roughness and Acoustic Disturbances

Boundary-layer receptivity and stability of Mach 6 flow over smooth and rough 7 half-angle sharp-tipped cones are numerically investigated. The receptivity of the boundary layer to slow acoustic disturbances, fast acoustic disturbances, and vortical disturbances are considered. The effects of two-dimensional isolated and distributed roughness on the receptivity and stability are also simulated. The results show that the instability waves are generated in the leading edge region and that the boundary layer is much more receptive to slow acoustic waves than to the fast waves. Vortical disturbances also generate unstable second modes, however the receptivity coefficients are smaller than that of the slow acoustic wave. An isolated two-dimensional roughness element of height h/delta =1/4 did not produce any difference in the receptivity or in the stability of the boundary layer. Distributed roughness elements produced a small decrease in the receptivity coefficient and also stabilized the boundary layer by small amounts.

Balakumar, Ponnampalam↗
NASA NTRS1996

The effects of complex chemistry on triple flames

The structure, ignition, and stabilization mechanisms for a methanol (CH3OH)-air triple flame are studied using Direct Numerical Simulations (DNS). The methanol (CH3OH)-air triple flame is found to burn with an asymmetric shape due to the different chemical and transport processes characterizing the mixture. The excess fuel, methanol (CH3OH), on the rich premixed flame branch is replaced by more stable fuels CO and H2, which burn at the diffusion flame. On the lean premixed flame side, a higher concentration of O2 leaks through to the diffusion flame. The general structure of the triple point features the contribution of both differential diffusion of radicals and heat. A mixture fraction-temperature phase plane description of the triple flame structure is proposed to highlight some interesting features in partially premixed combustion. The effects of differential diffusion at the triple point add to the contribution of hydrodynamic effects in the stabilization of the triple flame. Differential diffusion effects are measured using two methods: a direct computation using diffusion velocities and an indirect computation based on the difference between the normalized mixture fractions of C and H. The mixture fraction approach does not clearly identify the effects of differential diffusion, in particular at the curved triple point, because of ambiguities in the contribution of carbon and hydrogen atoms' carrying species.

Echekki, T.↗
NASA NTRS2001

Rotor Design Options for Improving V-22 Whirl-Mode Stability

With current tilt-rotor technology, very stiff, thick wings of limited aspect ratio are essential to meet aeroelastic stability requirements, which severely limits cruise efficiency and maximum speed. Larger and more efficient tiltrotors; will need longer and lighter wings, for which whirl-mode flutter is a serious design issue. Numerous approaches to improving the whirl-mode airspeed boundary have been investigated, typically focusing on wing design, active stability augmentation, and variable geometry rotors. The research to be reported applies the much simpler approaches of sweeping the outboard blade sections and extending balance weights at the tip. Improving proprotor whirl-mode stability margins is an ongoing research activity at NASA Ames Research Center. Previous efforts focussed on the XV-15 rotor. This is the first application of these design principles to a full V-22 rotor model in CAMRAD II (a previous publication presented results for a preliminary analytical model). The analysis to be discussed in the proposed paper includes a dual load-path blade model, a compressible-flow aerodynamic model, a representation of swashplate/rotor modal coupling, and other new features. The paper will include discussions of the mechanisms by which blade sweep and tip-mass offsets improve whirl-flutter margins. The V-22 was modeled with CAMRAD II; the rotor model is shown. It is also shows an example swept blade derived from the V-22 rotor system. For this rotor, the primary significance of sweep is the improved whirl-flutter boundary, not the reduced Mach-number effects. An offset tip mass is also shown.

Acree, C. W., Jr.↗
NASA NTRS2000

Effects of Rotor Design Variations on Tiltrotor Whirl-Mode Stability

Further increases in tiltrotor speeds are limited by coupled wing/rotor whirl-mode aeroelastic instability. Increased power, thrust, and efficiency are not enough: the whirl-mode stability boundary must also be improved. With current technology, very stiff, thick wings of limited aspect ratio are essential to meet the stability requirements, which severely limits cruise efficiency and maximum speed. Larger and more efficient tiltrotors will need longer and lighter wings, for which whirl-mode flutter is a serious design issue. Numerous approaches to improving the whirl-mode airspeed boundary have been investigated, typically focusing on wing design, active stability augmentation, and variable geometry rotors. The research reported herein applies the much simpler approach of sweeping the outboard blade sections. The effects-of control-system stiffness were also studied. Reference covers the first phase of the current research and discusses the evolution of the concept.

Acree, C. W., Jr.↗
NASA NTRS1977

The use of a movable telescoping end mass system for the time-optimal control of spinning spacecraft

The time-optimal control of a spin-stabilized spacecraft with a movable telescoping appendage (boom) is considered analytically and numerically. The motion of a control mass at the end of the boom is determined such that the terminal time will be minimized for two-axis control of a symmetric spacecraft. The equations of rotational motion are linearized about the desired state of spin about the symmetry axis. The equations for the transverse angular velocity components have the form of a coupled two dimensional harmonic oscillator with boom motion as a control force. The control function which brings the system to the desired state is known to be a series of positive and negative pulses. If the initial state is such that the system can be driven to rest in a single switch, the responses, switching and final times, and required boom motion may be determined analytically. Some typical numerical results based on these solutions are discussed.

Bainum, P. M.↗
NASA NTRS1977

Finite element analysis of low speed viscous and inviscid aerodynamic flows

A weak interaction solution algorithm was established for aerodynamic flow about an isolated airfoil. Finite element numerical methodology was applied to solution of each of differential equations governing potential flow, and viscous and turbulent boundary layer and wake flow downstream of the sharp trailing edge. The algorithm accounts for computed viscous displacement effects on the potential flow. Closure for turbulence was accomplished using both first and second order models. The COMOC finite element fluid mechanics computer program was modified to solve the identified equation systems for two dimensional flows. A numerical program was completed to determine factors affecting solution accuracy, convergence and stability for the combined potential, boundary layer, and parabolic Navier-Stokes equation systems. Good accuracy and convergence are demonstrated. Each solution is obtained within the identical finite element framework of COMOC.

Baker, A. J.↗
NASA NTRS1977

Computer simulation of plasma and N-body problems

Any spiral structure in computer-generated galaxies is generally short lived and the final state is a rotating bar. The bar thus obtained rotates more slowly than the stars. It has been argued that core/halo components have a stabilizing effect on galaxies and result in longer lived spiral structure. However, numerical experiments with large fixed stellar components representing the core/halo component show that multiarmed spiral structure develops and persists for many rotations but only in an evolving manner. That is, the spiral structure is either wound up into a tight pattern or it is wound up and then reappears again. A recent study of the effect of fixed core/halo components does show that the bar instability is indeed inhibited by a sufficiently large fixed component. The present study determines the effect of a self-consistent (rather than fixed) core/halo component in order to show whether there are any instabilities (such as two-stream) or other important interactions present that may be suppressed with a fixed core. Also studied were the effects of finite thickness of the disk and of three-dimensional essentially spherical core/halo components.

Harries, W. L.↗
NASA NTRS1983

Periodic motions (close to stationary) of an axisymmetric satellite with magnetic damping

Close to stationary periodic motions of an axisymmetric satellite in a circular orbit are considered. The satellite was equipped with a spherical magnetic damper. The investigation was conducted on the assumption that a strong magnet was installed on the damper float. Stationary rotations of the satellite around the axis of symmetry are selected as the generating solutions. The solutions are constructed in the form of power series of the smaller parameter, and they are extended numerically to the region of random values of the damping coefficient. The stability of the resulting solutions was investigated.

Ovchinnikov, M. Y.↗
NASA NTRS1983

Inelastic stress analyses at finite deformation through complementary energy approaches

A new hybrid-stress finite element algorithm, suitable for analyses of large, quasistatic, inelastic deformations, is presented. The algorithm is based upon a generalization of de Veubeke's (1972) complementary energy principle. The principal variables in the formulation are the nominal stress rate and spin, and the resulting finite element equations are discrete versions of the equations of compatibility and angular momentum balance. The algorithm produces true rates, time derivatives, as opposed to 'increments'. There results a boundary value problem (for stress rate and velocity) and an initial value problem (for total stress and deformation). A discussion of the numerical treatment of the boundary value problem is followed by a detailed examination of the numerical treatment of the initial value problem, covering the topics of efficiency, stability, and objectivity. The paper is closed with a set of examples, finite homogeneous deformation problems, which serve to bring out important aspects of the algorithm.

Atluri, S. N.↗
NASA NTRS1983

Analyses of large quasistatic deformations of inelastic bodies by a new hybrid-stress finite element algorithm

A new hybrid-stress finite element algorithm, suitable for analyses of large, quasistatic, inelastic deformations, is presented. The algorithm is base upon a generalization of de Veubeke's complementary energy principle. The principal variables in the formulation are the nominal stress rate and spin, and thg resulting finite element equations are discrete versions of the equations of compatibility and angular momentum balance. The algorithm produces true rates, time derivatives, as opposed to 'increments'. There results a complete separation of the boundary value problem (for stress rate and velocity) and the initial value problem (for total stress and deformation); hence, their numerical treatments are essentially independent. After a fairly comprehensive discussion of the numerical treatment of the boundary value problem, we launch into a detailed examination of the numerical treatment of the initial value problem, covering the topics of efficiency, stability and objectivity. The paper is closed with a set of examples, finite homogeneous deformation problems, which serve to bring out important aspects of the algorithm.

Reed, K. W.↗
NASA NTRS1984

Optimization and mechanisms of mistuning in cascades

In the present inverse design procedure for the optimum mistuning of a high bypass ratio shroudless fan that is modeled as a cascade of blades (each with a single torsional degree-of-freedom), linearized supersonic aerodynamic theory is used to compute the unsteady aerodynamic forces in the 'influence coefficient' form at a typical blade section. The mistuning pattern is then numerically optimized in order to achieve a specified increase in the aeroelastic stability margin with a minimum amount of mistuning. If the blades are self-damped, an optimized mistuning pattern can be found that achieves a given stability margin for a much lower level of mistuning than required for the alternate mistuning pattern, which requires only two blade frequencies and is relatively insensitive to implementation errors.

Crawley, E. F.↗
NASA NTRS1987

Stability of streamwise vortices

A brief overview of some theoretical and computational studies of the stability of streamwise vortices is given. The local induction model and classical hydrodynamic vortex stability theories are discussed in some detail. The importance of the three-dimensionality of the mean velocity profile to the results of stability calculations is discussed briefly. The mean velocity profile is provided by employing the similarity solution of Donaldson and Sullivan. The global method of Bridges and Morris was chosen for the spatial stability calculations for the nonlinear eigenvalue problem. In order to test the numerical method, a second order accurate central difference scheme was used to obtain the coefficient matrices. It was shown that a second order finite difference method lacks the required accuracy for global eigenvalue calculations. Finally the problem was formulated using spectral methods and a truncated Chebyshev series.

Khorrami, M. K.↗
NASA NTRS1987

The dynamics and control of large flexible space structures X, part 1

The effect of delay in the control system input on the stability of a continuously acting controller which is designed without considering the delay is studied. The stability analysis of a second order plant is studied analytically and verified numerically. For this example it is found that the system becomes unstable for a delay which is equivalent to only 16 percent of its natural period of motion. It is also observed that even a small amount of natural damping in the system can increase the amount of delay that can be tolerated before the onset of instability. The delay problem is formulated in the discrete time domain and an analysis procedure suggested. The maximum principle from optimal control theory is applied to minimize the time required for the slewing of a general rigid spacecraft. The slewing motion need not be restricted to a single axis maneuver. The minimum slewing time is calculated based on a quasi-linearization algorithm for the resulting two point boundary value problem. Numerical examples based on the rigidized in-orbit model of the SCOLE also include the more general reflector line-of-sight slewing maneuvers.

Bainum, Peter M.↗
NASA NTRS1988

Design for performance enhancement in feedback control systems with multiple saturating nonlinearities

A systematic control design methodology is introduced for multi-input/multi-output systems with multiple saturations. The methodology can be applied to stable and unstable open loop plants with magnitude and/or rate control saturations and to systems in which state limitations are desired. This new methodology is a substantial improvement over previous heuristic single-input/single-output approaches. The idea is to introduce a supervisor loop so that when the references and/or disturbances are sufficiently small, the control system operates linearly as designed. For signals large enough to cause saturations, the control law is modified in such a way to ensure stability and to preserve, to the extent possible, the behavior of the linear control design. Key benefits of this methodology are: the modified compensator never produces saturating control signals, integrators and/or slow dynamics in the compensator never windup, the directional properties of the controls are maintained, and the closed loop system has certain guaranteed stability properties. The advantages of the new design methodology are illustrated by numerous simulations, including the multivariable longitudinal control of modified models of the F-8 (stable) and F-16 (unstable) aircraft.

Kapasouris, Petros↗
NASA NTRS1989

Modeling of high speed chemically reacting flow-fields

The SPARK3D and SPARK3D-PNS computer programs were developed to model 3-D supersonic, chemically reacting flow-fields. The SPARK3D code is a full Navier-Stokes solver, and is suitable for use in scramjet combustors and other regions where recirculation may be present. The SPARK3D-PNS is a parabolized Navier-Stokes solver and provides an efficient means of calculating steady-state combustor far-fields and nozzles. Each code has a generalized chemistry package, making modeling of any chemically reacting flow possible. Research activities by the Langley group range from addressing fundamental theoretical issues to simulating problems of practical importance. Algorithmic development includes work on higher order and upwind spatial difference schemes. Direct numerical simulations employ these algorithms to address the fundamental issues of flow stability and transition, and the chemical reaction of supersonic mixing layers and jets. It is believed that this work will lend greater insight into phenomenological model development for simulating supersonic chemically reacting flows in practical combustors. Currently, the SPARK3D and SPARK3D-PNS codes are used to study problems of engineering interest, including various injector designs and 3-D combustor-nozzle configurations. Examples, which demonstrate the capabilities of each code are presented.

Drummond, J. P.↗