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

Stability, transition and turbulence

A glimpse is provided of the research program in stability, transition and turbulence based on numerical simulations. This program includes both the so-called abrupt and the restrained transition processes. Attention is confined to the prototype problems of channel flow and the parallel boundary layer in the former category and the Taylor-Couette flow in the latter category. It covers both incompressible flows and supersonic flows. Some representative results are presented.

Hussaini, M. Y.↗
NASA NTRS1987

Stability, transition and turbulence

A glimpse is provided of the research program in stability, transition, and turbulence based on numerical simulations. This program includes both the so-called abrupt and the restrained transition processes. Attention is confined to the prototype problems of channel flow and the parallel boundary layer in the former category and the Taylor-Couette flow in the latter category. It covers both incompressible flows and supersonic flows. Some representative results are presented.

Hussaini, M. Y.↗
NASA NTRS1989

A simple time-accurate turbomachinery algorithm with numerical solutions of an uneven blade count configuration

The present computer algorithm for the time-accurate flow analysis of rotating turbomachines is based on the finite-volume method and employs a high-resolution approximate Riemann solver for interface flux definitions and an implicit numerical scheme that possesses apparent unconditional stability. Block-block interfaces, including dynamic ones, are treated in such a way as to mimic interior block communication. The turbomachine configurations treated by way of illustration are 8-8-bladed and 11-9-bladed versions of a contrarotating unducted fan engine.

Janus, J. Mark↗
NASA NTRS1989

A study of the temporal stability of multiple cell vortices

The effect of initial mean velocity field on the stability characteristics of longitudinal vortices is documented in detail. The temporal stability of isolated multiple cell vortices is considered. The types of vortices studied include single cell as well as two and three cell vortices. It is shown that cell multiplicity in the vortex core has drastic effects on the stability characteristics. On the basis of numerical calculations, it is concluded that the growth rates of instabilities in multiple cell vortices are substantially larger (two to threefold increases are observed) than those of a single cell vortex. It is also determined that there is a substantial increase in the effective range of axial and azimuthal wavenumbers where instabilities are present. But most importantly, there is the appearance of a variety of viscous modes of instability. In the case of vortices, these latter instabilities which highlight the importance of viscous forces have never been reported before. These effects are discussed in detail for the case of a two cell vortex.

Khorrami, Mehdi R.↗
NASA NTRS1989

A controller design for multi-body large angle maneuvers

Active large angle slewing maneuvers of a multi-body flexible dynamic system are investigated. An appropriate state variable transformation and a feedback linearization technique are employed to transform the dynamics of the nonlinear system to a new state that is more amenable to control design procedures. Closed-loop feedback algorithms are implemented to perform slewing maneuvers, while simultaneously suppressing flexural vibrations of the system. Stability of this class of nonlinear systems is also investigated, whereby a sufficient condition for asymptotic stability of the system is established. Numerical examples are presented to demonstrate the proposed active control algorithms.

Ghaemmaghami, Peiman↗
NASA NTRS1991

Frequency effects on the stability of a journal bearing for periodic loading

The stability of a journal bearing is numerically predicted when a unidirectional periodic external load is applied. The analysis is performed using a cavitation algorithm, which mimics the Jakobsson-Floberg and Olsson (JFO) theory by accounting for the mass balance through the complete bearing. Hence, the history of the film is taken into consideration. The loading pattern is taken to be sinusoidal and the frequency of the load cycle is varied. The results are compared with the predictions using Reynolds boundary conditions for both film rupture and reformation. With such comparisons, the need for accurately predicting the cavitation regions for complex loading patterns is clearly demonstrated. For a particular frequency of loading, the effects of mass, amplitude of load variation and frequency of journal speed are also investigated. The journal trajectories, transient variations in fluid film forces, net surface velocity and minimum film thickness, and pressure profiles are also presented.

Vijayaraghavan, D.↗
NASA NTRS1992

Frequency effects on the stability of a journal bearing for periodic loading

The stability of a journal bearing is numerically predicted when an unidirectional periodic external load is applied. The analysis is performed using a cavitation algorithm, which mimics the JFO theory by accounting for the mass balance through the complete bearing. Hence, the history of the film is taken into consideration. The loading pattern is taken to be sinusoidal and the frequency of the load cycle is varied. The results are compared with the predictions using Reynolds boundary conditions for both film rupture and reformation. With such comparisons, the need for accurately predicting the cavitation regions for complex loading patterns is clearly demonstrated. For a particular frequency of loading, the effects of mass, amplitude of load vibration and frequency of journal speed are also investigated.

Vijayaraghavan, D.↗
NASA NTRS1992

Debris about asteroids: Where and how much?

We summarize several recent findings on the size and shape of the region within which material can stably orbit an asteroid. If the asteroid (with assumed density 2.38 g/cu cm) circles the Sun at 2.55 AU, co-planar prograde material will remain trapped whenever started on unperturbed circular orbits at less than about 220 R(sub A) (asteroid radii); co-planar retrograde particles are stable out twice as far. Our 3-D stability surface, which encloses several hundred numerically calculated orbits that start with various inclinations, is shaped like a sphere with its top and bottom sliced off; its dimensions scale like the Hill radius =(mu/3)(exp 1/3)R, where mu is the asteroid-to-solar mass ratio and R is the asteroid's orbital radius. If the asteroid moves along an elliptical orbit, a fairly reliable indicator of the dimensions of the hazard zone is the size of its Hill sphere at the orbit's pericenter. Grains with radii less than a few mm will be lost through the action of radiation forces which can induce escape or cause collisions with the asteroid on times scales of a few years; interplanetary micrometeoroids produce collisional break-up of these particles in approximately 10(exp 4) yrs. The effects of Jupiter and of asteroids that pass close to the target asteroid allow particles to diffuse from the system, again shrinking the hazard zone. None of the considered sources-primordial formation, debris spalled off the asteroid during micrometeoroid impact, captured interplanetary particles, feeder satellites, etc., seem capable of densely populating distant orbits from the asteroid. No certain detections of debris clouds or of binary asteroids have been made. Thus, it seems highly unlikely that a spacecraft fly-by targeted at 100 R(sub A) from the asteroid over its orbital pole would encounter any material.

Burns, Joseph A.↗
NASA NTRS1992

Finite element computation of incompressible flows

Among the solution techniques presented for FEM computation of incompressible flows are stabilized formulations designed to prevent the numerical oscillations that may be generated by (1) the presence of dominant advection terms in the governing equations, or (2) the failure to use an acceptable combination of interpolation functions to represent velocity and pressure fields. It is shown that the space-time formulation, in conjunction with the Galerkin-least squares methods, can be used to compute steady incompressible flows involving moving boundaries and interfaces; it is recommended that iteration methods be used to reduce the associated computational costs.

Tezduyar, T. E.↗
NASA NTRS1994

Semigroup approximation and robust stabilization of distributed parameter systems

Theoretical results that enable rigorous statements of convergence and exponential stability of Galerkin approximations of LQR controls for infinite dimensional, or distributed parameter, systems have proliferated over the past ten years. In addition, extensive progress has been made over the same time period in the derivation of robust control design strategies for finite dimensional systems. However, the study of the convergence of robust finite dimensional controllers to robust controllers for infinite dimensional systems remains an active area of research. We consider a class of soft-constrained differential games evolving in a Hilbert space. Under certain conditions, a saddle point control can be given in feedback form in terms of a solution to a Riccati equation. By considering a related LQR problem, we can show a convergence result for finite dimensional approximations of this differential game. This yields a computational algorithm for the feedback gain that can be derived from similar strategies employed in infinite dimensional LQR control design problems. The approach described in this paper also inherits the additional properties of stability robustness common to game theoretic methods in finite dimensional analysis. These theoretical convergence and stability results are verified in several numerical experiments.

Kurdila, A. J.↗
NASA NTRS1994

Effect of nonzero surface admittance on receptivity and stability of compressible boundary layer

The effect of small-amplitude short-scale variations in surface admittance on the acoustic receptivity and stability of two-dimensional compressible boundary layers is examined. In the linearized limit, the two problems are shown to be related both physically and mathematically. This connection between the two problems is used, in conjunction with some previously reported receptivity results, to infer the modification of stability properties due to surface permeability. Numerical calculations are carried out for a self-similar flat-plate boundary layer at subsonic and low supersonic speeds. Variations in mean suction velocity at the perforated admittance surface can also induce receptivity to an acoustic wave. For a subsonic boundary layer, the dependence of admittance-induced receptivity on the acoustic-wave orientation is significantly different from that of the receptivity produced via mean suction variation. The admittance-induced receptivity is generally independent of the angle of acoustic incidence, except in a relatively narrow range of upstream-traveling waves for which the receptivity becomes weaker. However, this range of angles is precisely that for which the suction-induced receptivity tends to be large. At supersonic Mach numbers, the admittance-induced receptivity to slow acoustic models is relatively weaker than that in the case of the fast acoustic modes. We also find that purely real values for the surface admittance tend to have a destabilizing effect on the evolution of an instability wave over a slightly permeable surface. The limits on the validity of the linearized approximation are also assessed in one specific case.

Choudhari, Meelan↗
NASA NTRS1994

Linear and nonlinear dynamic analysis of redundant load path bearingless rotor systems

The goal of this research is to develop the transfer matrix method to treat nonlinear autonomous boundary value problems with multiple branches. The application is the complete nonlinear aeroelastic analysis of multiple-branched rotor blades. Once the development is complete, it can be incorporated into the existing transfer matrix analyses. There are several difficulties to be overcome in reaching this objective. The conventional transfer matrix method is limited in that it is applicable only to linear branch chain-like structures, but consideration of multiple branch modeling is important for bearingless rotors. Also, hingeless and bearingless rotor blade dynamic characteristics (particularly their aeroelasticity problems) are inherently nonlinear. The nonlinear equations of motion and the multiple-branched boundary value problem are treated together using a direct transfer matrix method. First, the formulation is applied to a nonlinear single-branch blade to validate the nonlinear portion of the formulation. The nonlinear system of equations is iteratively solved using a form of Newton-Raphson iteration scheme developed for differential equations of continuous systems. The formulation is then applied to determine the nonlinear steady state trim and aeroelastic stability of a rotor blade in hover with two branches at the root. A comprehensive computer program is developed and is used to obtain numerical results for the (1) free vibration, (2) nonlinearly deformed steady state, (3) free vibration about the nonlinearly deformed steady state, and (4) aeroelastic stability tasks. The numerical results obtained by the present method agree with results from other methods.

Murthy, V. R.↗
NASA NTRS1995

Hydrodynamics of binary coalescence. 2: Polytropes with gamma = 5/3

We present a new numerical study of the equilibrium and stability properties of close binary systems. We use the smoothed-particle hydrodynamics (SPH) technique both to construct accurate equilibrium configurations in three dimensions and to follow their hydrodynamic evolution. We adopt a simple polytropic equation of state p = K(sub rho)(exp gamma) with gamma = 5/3 and K = constant within each star, applicable to low-mass degenerate dwarfs as well as low-mass main-sequence stars. For degenerate configurations, we set the two polytropic constants equal, K = K prime, independent of the mass ratio. For main-sequence stars, we adjust K and K prime so as to obtain a simple mass-radius relation of the form R/R prime = M/M prime, where R prime and M prime are the radius and mass of the secondary. Along a sequence of binary equilibrium configurations for two identical stars, we demonstrate the existence of both secular and dynamical instabilities, confirming directly the results of recent analytic work. We use the SPH method to calculate the nonlinear development of the dynamical instability and to determine the final fate of the system. We find that the two stars merge together into a single, rapidly rotating object in just a few orbital periods. Equilibrium sequences are also constructed for systems containing two nonidentical stars. These sequences terminate at a Roche limit, which we can determine very accurately using SPH. For two low-mass main-sequence stars with mass ratio q approximately less than 0.4 we find that the (synchronized) Roche limit configuration is secularly unstable. We discuss the implications of our results for the evolution of double white-dwarf systems and W Ursae Majoris binaries.

Rasio, Frederic A.↗
NASA NTRS1999

Implementation and Validation of the Chien k-epsilon Turbulence Model in the Wind Navier-Stokes Code

The two equation k-epsilon turbulence model of Chien has been implemented in the WIND Navier-Stokes flow solver. Details of the numerical solution algorithm, initialization procedure, and stability enhancements are described. Results obtained with this version of the model are compared with those from the Chien k-epsilon model in the NPARC Navier-Stokes code and from the WIND SST model for three validation cases: the incompressible flow over a smooth flat plate, the incompressible flow over a backward facing step, and the shock-induced flow separation inside a transonic diffuser. The k-epsilon model results indicate that the WIND model functions very similarly to that in NPARC, though the WIND code appears to he slightly more accurate in the treatment of the near-wall region. Comparisons of the k-epsilon model results with those from the SST model were less definitive, as each model exhibited strengths and weaknesses for each particular case.

Yoder, Dennis A.↗
NASA NTRS2000

Model Refinement Using Eigensystem Assignment

IA novel approach for the refinement of finite-element-based analytical models of flexible structures is presented. The proposed approach models the possible refinements in the mass, damping, and stiffness matrices of the finite element model in the form of a constant gain feedback with acceleration, velocity, and displacement measurements, respectively. Once the free elements of the structural matrices have been defined, the problem of model refinement reduces to obtaining position, velocity, and acceleration gain matrices with appropriate sparsity that reassign a desired subset of the eigenvalues of the model, along with partial mode shapes, from their baseline values to those obtained from system identification test data. A sequential procedure is used to assign one conjugate pair of eigenvalues at each step using symmetric output feedback gain matrices, and the eigenvectors are partially assigned, while ensuring that the eigenvalues assigned in the previous steps are not disturbed. The procedure can also impose that gain matrices be dissipative to guarantee the stability of the refined model. A numerical example, involving finite element model refinement for a structural testbed at NASA Langley Research Center (Controls-Structures-Interaction Evolutionary model) is presented to demonstrate the feasibility of the proposed approach.

Maghami, Peiman G.↗
NASA NTRS2001

Application of Aeroelastic Solvers Based on Navier Stokes Equations

The propulsion element of the NASA Advanced Subsonic Technology (AST) initiative is directed towards increasing the overall efficiency of current aircraft engines. This effort requires an increase in the efficiency of various components, such as fans, compressors, turbines etc. Improvement in engine efficiency can be accomplished through the use of lighter materials, larger diameter fans and/or higher-pressure ratio compressors. However, each of these has the potential to result in aeroelastic problems such as flutter or forced response. To address the aeroelastic problems, the Structural Dynamics Branch of NASA Glenn has been involved in the development of numerical capabilities for analyzing the aeroelastic stability characteristics and forced response of wide chord fans, multi-stage compressors and turbines. In order to design an engine to safely perform a set of desired tasks, accurate information of the stresses on the blade during the entire cycle of blade motion is required. This requirement in turn demands that accurate knowledge of steady and unsteady blade loading is available. To obtain the steady and unsteady aerodynamic forces for the complex flows around the engine components, for the flow regimes encountered by the rotor, an advanced compressible Navier-Stokes solver is required. A finite volume based Navier-Stokes solver has been developed at Mississippi State University (MSU) for solving the flow field around multistage rotors. The focus of the current research effort, under NASA Cooperative Agreement NCC3- 596 was on developing an aeroelastic analysis code (entitled TURBO-AE) based on the Navier-Stokes solver developed by MSU. The TURBO-AE code has been developed for flutter analysis of turbomachine components and delivered to NASA and its industry partners. The code has been verified. validated and is being applied by NASA Glenn and by aircraft engine manufacturers to analyze the aeroelastic stability characteristics of modem fans, compressors and turbines.

Keith, Theo G., Jr.↗
NASA NTRS2001

Discontinuous Galerkin Methods for NonLinear Differential Systems

This talk considers simplified finite element discretization techniques for first-order systems of conservation laws equipped with a convex (entropy) extension. Using newly developed techniques in entropy symmetrization theory, simplified forms of the discontinuous Galerkin (DG) finite element method have been developed and analyzed. The use of symmetrization variables yields numerical schemes which inherit global entropy stability properties of the PDE (partial differential equation) system. Central to the development of the simplified DG methods is the Eigenvalue Scaling Theorem which characterizes right symmetrizers of an arbitrary first-order hyperbolic system in terms of scaled eigenvectors of the corresponding flux Jacobian matrices. A constructive proof is provided for the Eigenvalue Scaling Theorem with detailed consideration given to the Euler equations of gas dynamics and extended conservation law systems derivable as moments of the Boltzmann equation. Using results from kinetic Boltzmann moment closure theory, we then derive and prove energy stability for several approximate DG fluxes which have practical and theoretical merit.

Barth, Timothy↗
NASA NTRS2003

Gravity Effects Observed In Partially Premixed Flames

Partially premixed flames (PPFs) contain a rich premixed fuel air mixture in a pocket or stream, and, for complete combustion to occur, they require the transport of oxidizer from an appropriately oxidizer-rich (or fuel-lean) mixture that is present in another pocket or stream. Partial oxidation reactions occur in fuel-rich portions of the mixture and any remaining unburned fuel and/or intermediate species are consumed in the oxidizer-rich portions. Partial premixing, therefore, represents that condition when the equivalence ratio (phi) in one portion of the flowfield is greater than unity, and in another section its value is less than unity. In general, for combustion to occur efficiently, the global equivalence ratio is in the range fuel-lean to stoichiometric. These flames can be established by design by placing a fuel-rich mixture in contact with a fuel-lean mixture, but they also occur otherwise in many practical systems, which include nonpremixed lifted flames, turbulent nonpremixed combustion, spray flames, and unwanted fires. Other practical applications of PPFs are reported elsewhere. Although extensive experimental studies have been conducted on premixed and nonpremixed flames under microgravity, there is a absence of previous experimental work on burner stabilized PPFs in this regard. Previous numerical studies by our group employing a detailed numerical model showed gravity effects to be significant on the PPF structure. We report on the results of microgravity experiments conducted on two-dimensional (established on a Wolfhard-Parker slot burner) and axisymmetric flames (on a coannular burner) that were investigated in a self-contained multipurpose rig. Thermocouple and radiometer data were also used to characterize the thermal transport in the flame.

Puri, Ishwar K.↗