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 757 records · Page 42
All sourcesNASA NTRSDOE OSTI
NASA NTRS1999

Finite Differences and Collocation Methods for the Solution of the Two Dimensional Heat Equation

In this paper we combine finite difference approximations (for spatial derivatives) and collocation techniques (for the time component) to numerically solve the two dimensional heat equation. We employ respectively a second-order and a fourth-order schemes for the spatial derivatives and the discretization method gives rise to a linear system of equations. We show that the matrix of the system is non-singular. Numerical experiments carried out on serial computers, show the unconditional stability of the proposed method and the high accuracy achieved by the fourth-order scheme.

Kouatchou, Jules↗
NASA NTRS1999

Studies of Premixed Laminar and Turbulent Flames at Microgravity

Several topics relating to combustion limits in premixed flames at reduced gravity have been studied. These topics include: (1) flame balls; (2) numerical simulation of flame ball and planar flame structure and stability; (3) experimental simulation of buoyancy effects in premixed flames using aqueous autocatalytic reactions; and (4) premixed flame propagation in Hele-Shaw cells.

Abid, M.↗
NASA NTRS2002

Lattice Boltzmann Equation On a 2D Rectangular Grid

We construct a multi-relaxation lattice Boltzmann model on a two-dimensional rectangular grid. The model is partly inspired by a previous work of Koelman to construct a lattice BGK model on a two-dimensional rectangular grid. The linearized dispersion equation is analyzed to obtain the constraints on the isotropy of the transport coefficients and Galilean invariance for various wave propagations in the model. The linear stability of the model is also studied. The model is numerically tested for three cases: (a) a vortex moving with a constant velocity on a mesh periodic boundary conditions; (b) Poiseuille flow with an arbitras~y inclined angle with respect to the lattice orientation: and (c) a cylinder &symmetrically placed in a channel. The numerical results of these tests are compared with either analytic solutions or the results obtained by other methods. Satisfactory results are obtained for the numerical simulations.

Bouzidi, MHamed↗
NASA NTRS2006

Large-Scale Hybrid Motor Testing

Hybrid rocket motors can be successfully demonstrated at a small scale virtually anywhere. There have been many suitcase sized portable test stands assembled for demonstration of hybrids. They show the safety of hybrid rockets to the audiences. These small show motors and small laboratory scale motors can give comparative burn rate data for development of different fuel/oxidizer combinations, however questions that are always asked when hybrids are mentioned for large scale applications are - how do they scale and has it been shown in a large motor? To answer those questions, large scale motor testing is required to verify the hybrid motor at its true size. The necessity to conduct large-scale hybrid rocket motor tests to validate the burn rate from the small motors to application size has been documented in several place^'^^.^. Comparison of small scale hybrid data to that of larger scale data indicates that the fuel burn rate goes down with increasing port size, even with the same oxidizer flux. This trend holds for conventional hybrid motors with forward oxidizer injection and HTPB based fuels. While the reason this is occurring would make a great paper or study or thesis, it is not thoroughly understood at this time. Potential causes include the fact that since hybrid combustion is boundary layer driven, the larger port sizes reduce the interaction (radiation, mixing and heat transfer) from the core region of the port. This chapter focuses on some of the large, prototype sized testing of hybrid motors. The largest motors tested have been AMROC s 250K-lbf thrust motor at Edwards Air Force Base and the Hybrid Propulsion Demonstration Program s 250K-lbf thrust motor at Stennis Space Center. Numerous smaller tests were performed to support the burn rate, stability and scaling concepts that went into the development of those large motors.

Story, George↗
NASA NTRS2009

Parametric Optimization of Ares I Propellant Slosh Characteristics Using Frequency Response Criteria

A novel technique for developing propellant slosh damping requirements with respect to the stability characteristics of large flexible launch vehicles is presented. A numerical algorithm is devised which allows an automated software program to rapidly converge to pseudo-optimal solutions that minimize required propellant slosh damping for multiple tanks while maintaining constraints on the frequency response characteristics of a particular open-loop plant transfer function. An implementation of the algorithm using a high-order linear model of the Ares I plant dynamics considers all relevant dynamic interactions of flexible body modes, propellant slosh, and nozzle inertia effects. A high-resolution propellant damping requirements table is produced that can be used for baffle design. The method is demonstrated to provide exceptional speed and accuracy when compared with the alternative human-in-the-loop approach.

Orr, Jeb S.↗
NASA NTRS1990

Instability of a supersonic shock free elliptic jet

This paper presents a comparison of the measured and the computed spatial stability properties of an aspect ratio 2 supersonic shock free elliptic jet. The shock free nature of the elliptic jet provides an ideal test of validity of modeling the large scale coherent structures in the initial mixing region of noncircular supersonic jets with linear hydrodynamic stability theory. Both aerodynamic and acoustic data were measured. The data are used to compute the mean velocity profiles and to provide a description of the spatial composition of pressure waves in the elliptic jet. A hybrid numerical scheme is applied to solve the Rayleigh problem governing the inviscid linear spatial stability of the jet. The measured mean velocity profiles are used to provide a qualitative model for the cross sectional geometry and the smooth velocity profiles used in the stability analysis. Computational results are presented for several modes of instability at two jet cross sections. The acoustic measurements show that a varicose instability is the jet's perferred mode of motion. The stability analysis predicts that the Strouhal number varies linearly as a function of axial distance in the jet's initial mixing region, which is in good qualitative agreement with previous measurements.

Baty, Roy S.↗
NASA NTRS1991

Stability and bifurcation analysis of oscillators with piecewise-linear characteristics - A general approach

A general approach is developed for determining the periodic solutions and their stability of nonlinear oscillators with piecewise-smooth characteristics. A modified harmonic balance/Fourier transform procedure is devised for the analysis. The procedure avoids certain numerical differentiation employed previously in determining the periodic solutions, therefore enhancing the reliability and efficiency of the method. Stability of the solutions is determined via perturbations of their state variables. The method is applied to a forced oscillator interacting with a stop of finite stiffness. Flip and fold bifurcations are found to occur. This led to the identification of parameter ranges in which chaotic response occurred.

Noah, S. T.↗
NASA NTRS2003

Design Considerations of ISTAR Hydrocarbon Fueled Combustor Operating in Air Augmented Rocket, Ramjet and Scramjet Modes

The development and ground test of a rocket-based combined cycle (RBCC) propulsion system is being conducted as part of the NASA Marshall Space Flight Center (MSFC) Integrated System Test of an Airbreathing Rocket (ISTAR) program. The eventual flight vehicle (X-43B) is designed to support an air-launched self-powered Mach 0.7 to 7.0 demonstration of an RBCC engine through all of its airbreathing propulsion modes - air augmented rocket (AAR), ramjet (RJ), and scramjet (SJ). Through the use of analytical tools, numerical simulations, and experimental tests the ISTAR program is developing and validating a hydrocarbon-fueled RBCC combustor design methodology. This methodology will then be used to design an integrated RBCC propulsion system that produces robust ignition and combustion stability characteristics while maximizing combustion efficiency and minimizing drag losses. First order analytical and numerical methods used to design hydrocarbon-fueled combustors are discussed with emphasis on the methods and determination of requirements necessary to establish engine operability and performance characteristics.

Andreadis, Dean↗
NASA NTRS2002

Design Considerations of Istar Hydrocarbon Fueled Combustor Operating in Air Augmented Rocket, Ramjet and Scramjet Modes

The development and ground test of a rocket-based combined cycle (RBCC) propulsion system is being conducted as part of the NASA Marshall Space Flight Center (MSFC) Integrated System Test of an Airbreathing Rocket (ISTAR) program. The eventual flight vehicle (X-43B) is designed to support an air-launched self-powered Mach 0.7 to 7.0 demonstration of an RBCC engine through all of its airbreathing propulsion modes - air augmented rocket (AAR), ramjet (RJ), and scramjet (SJ). Through the use of analytical tools, numerical simulations, and experimental tests the ISTAR program is developing and validating a hydrocarbon-fueled RBCC combustor design methodology. This methodology will then be used to design an integrated RBCC propulsion system thai: produces robust ignition and combustion stability characteristics while maximizing combustion efficiency and minimizing drag losses. First order analytical and numerical methods used to design hydrocarbon-fueled combustors are discussed with emphasis on the methods and determination of requirements necessary to establish engine operability and performance characteristics.

Andreadis, Dean↗
NASA NTRS1980

Strongly turbulent stabilization of electron beam-plasma interactions

The stabilization of electron beam interactions due to strongly turbulent nonlinearities is studied analytically and numerically for a wide range of plasma parameters. A fluid mode coupling code is described in which the effects of electron and ion Landau damping and linear growth due to the energetic electron beam are included in a phenomenological manner. Stabilization of the instability is found to occur when the amplitudes of the unstable modes exceed the threshold of the oscillating two-stream instability. The coordinate space structure of the turbulent spectrum which results clearly shows that soliton-like structures are formed by this process. Phenomenological models of both the initial stabilization and the asymptotic states are developed. Scaling laws between the beam-plasma growth rate and the fluctuations in the fields and plasma density are found in both cases, and shown to be in good agreement with the results of the simulation.

Freund, H. P.↗
NASA NTRSDate not supplied

Fluid System Stability Analysis Techniques

Fluid systems, or networks, consist of multiple components that work together to achieve some desired thermofluid state. For any generic application, this fluid state can be a combination of the fluid pressure, flow rate, enthalpy, or species concentration. Fluid system components, such as pumps and valves, are often governed by nonlinear differential equations, resulting in complex component-to-component interactions. System-level fluid network stability occurs when the flow through the system can maintain a steady-state solution in the presence of small perturbations, which depends on these component interactions. System instability, however, can go undetected until issues arise during integrated system testing. This presentation explores a method for system designers to think of the fluid network as an assembly of components, each with their own thermofluid surfaces of partial stability, called nullclines. The intersections of all nullclines yields system-level solutions, called equilibrium points. When designers define operating points, they are tuning system parameters so that these equilibrium points move to the desired location in the thermofluid state plane. However, linearization theory shows us that the dynamic behavior around these points can be unstable. The local stability of these equilibrium points can be assessed analytically with eigen-analysis, or numerically by propagating state-plane samples to construct a phase portrait. Investigating a phase portrait can help designers gain a qualitative understanding of a system’s dynamic performance. This understanding can then help inform requirement definitions, component selection, and operational procedures. This presentation includes an example of the phase portrait technique on a system featuring a centrifugal pump and a back-pressure regulator (BPR). Numerical modeling of this system suggests that equilibrium points on the left-hand side of the pump curve are dynamically unstable.

Fluid System↗
NASA NTRS1986

Variable-Conductance Heat Pipes

In response to need to accurately and efficiently predict performance of variable-conductance heat pipes (VCHP's) incorporated in spacecraft thermalcontrol systems, computer code VCHPDA developed to interact with thermal analyzer programs such as SINDA (Systems Improved Numerical Differencing Analyzer). Calculates length of gas-blocked region and vapor temperature in active portion. Advantages of VCHPDA over prior programs improved accuracy, unconditional stability, and increased efficiency of solution resulting from novel approach and use of state-of-the-art numerical techniques for solving VCHP mathematical model. Code valuable tool in design and evaluation of advanced thermal-control systems using variable-conductance heat pipes. Written in FORTRAN IV for use on CDC 600 computers.

Antoniuk, D.↗
NASA NTRS1993

Numerical study of cathode emission constraints on cylindrical self-field MPD thruster performance

A stability equation relating thruster discharge current, argon propellant mass flow rate, and electrode geometry has been solved for a variety of cylindrical self-field MPD thruster configurations and discharge currents. Realistic values for cathode emission current densities were used to provide additional constraints on the thruster geometries. A two-dimensional MPD code was employed to provide better estimates of the maximum achievable specific impulse, thrust, and flow efficiency for cases of interest. The model results indicate that long life, cylindrical self-field MPD thrusters operated with argon propellant may not be able to provide specific impulse values in excess of 2100s. Alternate electrode geometries, applied magnetic fields, and/or low molecular weight propellants may be necessary to achieve higher values of specific impulse.

Lapointe, Michael R.↗
NASA NTRSDate not supplied

Fluid System Stability Analysis Techniques

Fluid systems, or networks, consist of multiple components that work together to achieve some desired thermofluid state. For any generic application, this fluid state can be a combination of the fluid pressure, flow rate, enthalpy, or species concentration. Fluid system components, such as pumps and valves, are often governed by nonlinear differential equations, resulting in complex component-to-component interactions. System-level fluid network stability occurs when the flow through the system can maintain a steady-state solution in the presence of small perturbations, which depends on these component interactions. System instability, however, can go undetected until issues arise during integrated system testing. This presentation explores a method for system designers to think of the fluid network as an assembly of components, each with their own thermofluid surfaces of partial stability, called nullclines. The intersections of all nullclines yields system-level solutions, called equilibrium points. When designers define operating points, they are tuning system parameters so that these equilibrium points move to the desired location in the thermofluid state plane. However, linearization theory shows us that the dynamic behavior around these points can be unstable. The local stability of these equilibrium points can be assessed analytically with eigen-analysis, or numerically by propagating state-plane samples to construct a phase portrait. Investigating a phase portrait can help designers gain a qualitative understanding of a system’s dynamic performance. This understanding can then help inform requirement definitions, component selection, and operational procedures. This presentation includes an example of the phase portrait technique on a system featuring a centrifugal pump and a back-pressure regulator (BPR). Numerical modeling of this system suggests that equilibrium points on the left-hand side of the pump curve are dynamically unstable.

Fluid System↗
NASA NTRS1990

Numerical evaluation of the performance of active noise control systems

This paper presents a generalized numerical technique for evaluating the optimal performance of active noise controllers. In this technique, the indirect BEM numerical procedures are used to derive the active noise controllers for optimal control of enclosed harmonic sound fields where the strength of the noise sources or the description of the enclosure boundary may not be known. The performance prediction for a single-input single-output system is presented, together with the analysis of the stability and observability of an active noise-control system employing detectors. The numerical procedures presented can be used for the design of both the physical configuration and the electronic components of the optimal active noise controller.

Mollo, C. G.↗
NASA NTRS1994

Dynamics of Numerics and CFD

This lecture attempts to illustrate the basic ideas of how the recent advances in nonlinear dynamical systems theory (dynamics) can provide new insights into the understanding of numerical algorithms used in solving nonlinear differential equations (DEs). Examples will be given of the use of dynamics to explain unusual phenomena that occur in numerics. The inadequacy of the use of linearized analysis for the understanding of long time behavior of nonlinear problems will be illustrated, and the role of dynamics in studying the nonlinear stability, accuracy, convergence property and efficiency of using time- dependent approaches to obtaining steady-state numerical solutions in computational fluid dynamics (CFD) will briefly be explained.

Yee, H. C.↗
NASA NTRS1985

Numerical solution of the vertical structure equation in the normal mode method

In the present model of multilayered stability stratification, aimed at obtaining the analytic eigensolutions of the vertical structure equation, each layer is characterized by its own static stability value. By requiring continuity of pressure and vertical velocity across each interface level, and by imposing suitable upper and lower boundary conditions, matching eigensolutions are obtained in terms of the Bessel functions. Attention is given to an explicit example of a double-layered stratified atmosphere which demonstrates the mathematical manipulations involved; the resultant vertical structure functions are used to check the accuracy of the numerical solutions by the finite difference and finite element methods.

Sasaki, Y. K.↗