Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linear equations”

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 1,279 records · Page 71

Prediction of blade stresses due to gust loading

An analysis is developed for investigating the response of a rotor-fuselage system in a three-dimensional gust field wherein the gust velocity conponents can have arbitrary variation in space and time. Each rotor blade undergoes flap bending, lag bending and torsional deflections. The blades are divided into beam elements and each element consists of fifteen nodal degrees of freedom. Quasi-steady strip theory is used to obtain the aerodynamic loads. Unsteady aerodynamic effects are introduced through dynamic inflow modeling. Dynamic stall and reverse flow effects are also included. The fuselage is allowed five degrees of freedom: vertical, longitudinal, lateral, pitch and roll motions. The gust response equations are linearized about the vehicle trim state and the blade steady-state deflected position, and then solved by time integration. The blade bending moments, which determine blade stresses, are evaluated using the force summation technique. Systematic studies are made to identify the importance of several parameters including dynamic stall, forward speed, lag stiffness, gust profile, gust penetration rate and gust velocity direction.

Bir, Gunjit↗

Stability analysis of double-peaked proton distribution functions in the solar wind

Twelve proton velocity distribution functions as measured for different solar wind speeds by the plasma analyzer at the Helios spacecraft are modeled analytically by a superposition of two bi-Maxwellians. The stability of these distributions is investigated by solving the linear dispersion equation for parallel propagation of both left- and right-hand polarized modes. Due to the temperature anisotropy of the main proton component the presence of the resonant left-hand ion cyclotron instability is found to be typical for high-speed streams and is independent of an ion beam. The results indicate that the ratio of the beam density to total particle density and the beam drift velocity are generally close to the threshold values necessary to drive a highly unstable nonresonant right-hand mode. This instability is found for one proton distribution and is a consequence of the beam particles. It is concluded that both the ion cyclotron instability and the ion beam instability act as a regulating mechanism on the shape of solar wind double humped proton distribution functions.

Leubner, Manfred P.↗

Down-to-earth supercomputing at NASA Ames

The requirements for computers for aerospace vehicle design problems, in particular fluid dynamics, are examined. The levels of approximation to Navier-Stokes equations are: (1) linearized inviscid form, (2) nonlinear inviscid equation, (3) Reynolds-averaged Navier-Stokes equation, (4) large eddy simulation with models for subgrid scale turbulence, and (5) full Navier-Stokes equations to calculate laminar/turbulent transition and turbulence dissipation. The Cray 2 with 256 megaword memory and 4.1 nsec clock and the Cray X-MP/48 with a 128 megaword memory and a 9.5 nsec clock of the Ames facility are described, and an example revealing their applicability is presented. The use of zonal methods to simulate aerodynamic flow fields numerically is discussed. Consideration is given to the development of a program for simulation of aerodynamic flow fields, the compilers, symbolic debuggers, the scheduling of multiple CPUs, and the analysis and storage of the data.

Stevens, Kenneth G., Jr.↗

Stationary disturbances in three-dimensional boundary layers over concave surfaces

A two-dimensional boundary layer on a concave surface is known to be susceptible to centrifugal instability which manifests itself in the form of stationary streamwise counter-rotating vortices commonly known as Goertler vortices. In this paper, the problem of the stability of a three-dimensional boundary layer on a concave surface is considered. Linear stability equations, including streamline and surface curvature effects, are solved for an infinitely swept wing. The results indicate that the Goetler vortex structure begins to transform into co-rotating vortices when the crossflow Reynolds number is increased. This transition is almost complete when the crossflow Reynolds number is in excess of 45. It is shown that the centrifugal effects destabilize the crossflow vortices. Some comparisons with available experimental results are made.

Collier, F. S., Jr.↗

Monte Carlo modeling of ionospheric oxygen acceleration by cyclotron resonance with broad-band electromagnetic turbulence

It is demonstrated that cyclotron resonance with observed electric field fluctuations is responsible for production of the oxygen-ion conics that are observed by the Dynamics Explorer 1 satellite in the central plasma-sheet region of the earth's magnetosphere. The ion-velocity distribution is described by a quasi-linear diffusion equation which is solved by the Monte Carlo technique. The acceleration produced by the observed wave spectrum agrees well with the ion observations, in both form and magnitude. This is believed to represent the first successful comparison of an observed conic with any theoretical model.

Retterer, John M.↗

Simulation of the electron acoustic instability for a finite-size electron beam system

Simulations of the electron acoustic instability for a finite-size electron beam system are performed with a particle-in-cell code to investigate the heating phenomena associated with the instability and the width of the heating region. The results show that the beam radiates electrostatic electron acoustic waves whose decay time outside the beam agrees with the spatial decay length derived from the linear dispersion equation. The ambient cold electrons in a diffusion region surrounding the beam are heated to a higher temperature by absorbing the radiated electron acoustic waves, with the heating occurring mainly in the parallel direction. In the heat diffusion region, the temperature of the cold electrons decreases with distance from the beam with a temperature gradient length smaller than the decay length of the wave energy. The results are discussed with respect to the DE 1 plasma and wave observations in the polar cusp region.

Lin, C. S.↗

The dynamics and control of the Space Station polar platform

The Space Station polar platform will carry a variety of earth observation instruments for NASA and the National Oceanic and Atmospheric Administration. In this paper, the asymmetrical platform is modeled as three connected rigid bodies. A generalized angular momentum equation is employed to derive the rotational equations of motion. These equations are linearized and used for preliminary sizing of control devices using a classical control approach. Two control systems are considered to stabilize the platform and satisfy the pointing requirements. The first system is composed of a single variable-speed, double-gimbaled momentum wheel and the second consists of three-reaction wheels. The performance of each system is assessed using a linear optimal control approach.

Wahbah, M. M.↗

Equations of motion for maneuvering flexible spacecraft

This paper is concerned with the derivation of the equations of motion for maneuvering flexible spacecraft both in orbit and in an earth-based laboratory. The structure is assumed to undergo large rigid-body maneuvers and small elastic deformations. A perturbation approach is presented in which the quantities defining the rigid-body maneuver are regarded as the unperturbed motion and the elastic motions and deviations from the rigid-body motions are regarded as the perturbed motion. The perturbation equations are linear, non-self-adjoint, and with time-dependent coefficients. A maneuver force distribution exciting the least amount of elastic deformation of the spacecraft is developed. Numerical results highlight the vibration caused by rotational maneuvers.

Meirovitch, L.↗

The dynamics and control of large-flexible space structures, part 10

A mathematical model is developed to predict the dynamics of the proposed orbiting Spacecraft Control Laboratory Experiment (SCOLE) during the station keeping phase. The equations of motion are derived using a Newton-Euler formulation. The model includes the effects of gravity, flexibility, and orbital dynamics. The control is assumed to be provided to the system through the Shuttle's three torquers, and through six actuators located by pairs at two points on the mast and at the mass center of the reflector. The modal shape functions are derived using the fourth order beam equation. The generic mode equations are derived to account for the effects of the control forces on the modal shape and frequencies. The equations are linearized about a nominal equilibrium position. The linear regulator theory is used to derive control laws for both the linear model of the rigidized SCOLE as well as that of the actual SCOLE including the first four flexible modes. The control strategy previously derived for the linear model of the rigidized SCOLE is applied to the nonlinear model of the same configuration of the system and preliminary single axis slewing maneuvers conducted. The results obtained confirm the applicability of the intuitive and appealing two-stage control strategy which would slew the SCOLE system, as if rigid to its desired position and then concentrate on damping out the residual flexible motions.

Bainum, Peter M.↗

Creep behavior of tungsten/niobium and tungsten/niobium-1 percent zirconium composites

The creep behavior and microstructural stability of tungsten fiber reinforced niobium and niobium 1 percent zirconium was determined at 1400 and 1500 K in order to assess the potential of this material for use in advanced space power systems. The creep behavior of the composite materials could be described by a power law creep equation. A linear relationship was found to exist between the minimum creep rate of the composite and the inverse of the composite creep rupture life. The composite materials had an order of magnitude increase in stress to achieve 1 percent creep strain and in rupture strength at test temperatures of 1400 and 1500 K compared to unreinforced material. The composite materials were also stronger than the unreinforced materials by an order of magnitude when density was taken into consideration. Results obtained on the creep behavior and microstructural stability of the composites show significant potential improvement in high temperature properties and mass reduction for space power system components.

Petrasek, Donald W.↗

Unstructured multigrid methods

The use of the multigrid procedure with a sequence of unnested coarser grids is discussed. Validity of the procedure is assessed by considering the solution of a single linear elliptic equation. It is demonstrated how a scheme with the optimum order of operations can be constructed. Application to the solution of the Euler equations is considered.

Loehner, R.↗

Parallel computer methods for eigenvalue extraction

A new numerical algorithm for the solution of large-order eigenproblems typically encountered in linear elastic finite element systems is presented. The architecture of parallel processing is used in the algorithm to achieve increased speed and efficiency of calculations. The algorithm is based on the frontal technique for the solution of linear simultaneous equations and the modified subspace eigenanalysis method for the solution of the eigenproblem. The advantages of this new algorithm in parallel computer architecture are discussed.

Akl, Fred↗

A mapped finite difference study of noise propagation in nonuniform ducts with mean flow

The primary objective of this work is to study noise propagation in acoustically lined variable area ducts with mean fluid flow. The method of study is numerical in nature and involves a body-fitted grid mapping procedure in conjunction with a factored-implicit finite difference technique. The mean fluid flow model used is two-dimensional, inviscid, irrotational, incompressible, and nonheat conducting. Fully-coupled solutions of the linearized gasdynamic equations are obtained for both positive and negative Mach numbers as well as for hard and soft wall conditions. The factored-implicit finite difference technique used did give rise to short wavelength perturbations, but these were dampened by the introduction of higher order artificial dissipation terms into the scheme. Results compared favorably with available numerical and experimental data.

Raad, Peter E.↗

Vacillations induced by interference of stationary and traveling planetary waves

The interference pattern produced when a traveling planetary wave propagates over a stationary forced wave is explored, examining the interference signature in a variety of diagnostics. The wave field is first restricted to a diatomic spectrum consisting of two components: a single stationary wave and a single monochromatic traveling wave. A simple barotropic normal mode propagating over a simple stationary plane wave is considered, and closed form solutions are obtained. The wave fields are then restricted spatially, providing more realistic structures without sacrificing the advantages of an analytical solution. Both stationary and traveling wave fields are calculated numerically with the linearized Primitive Equations in a realistic basic state. The mean flow reaction to the fluctuating eddy forcing which results from interference is derived. Synoptic geopotential behavior corresponding to the combined wave and mean flow fields is presented, and the synoptic signature in potential vorticity on isentropic surfaces is examined.

Salby, Murry L.↗

Electrostatic charge on a dust size distribution in a plasma

The capacitance of a grain immersed in a steady state plasma containing a size distribution of dust particles is studied. The grain charge is determined by assuming the equilibrium potential has been obtained by a simple balance of electron and ion collection currents. It is shown that the validity of the analytical treatment given here for the linearized Poisson equation is confined to a certain region of space. Within this region and starting at very small plasma Debye length lambda(D), the capacitance at first exhibits a monotonic increase with increasing lambda(D). The capacitance eventually reaches a maximum, followed by a monotonic decrease. The charge density of the dust in the plasma is found to be only a function of the lambda(D); there is no significant dependence on the interparticle spacing.

Houpis, Harry L. F.↗

Determination of electron density and temperature in non-LTE plasmas from spectral lines of impurity ions

A novel approach for calculating the populations of the excited Li-like ions C IV, N V, O VI, and Ne VIII is presented. The populations of the 2(2P), 3(2S), 3(2P), and 3(2D) electronic levels in these ions in optically thin plasmas with a broad range of electron density, N(e), and temperature, T(e), are determined from the collisional-radiative model by solving the system of rate equations for the production of excited ions; the equations are linear with respect to the excited ion populations, and the N(e) and T(e) are taken as independent variables. These populations are used to determine the ratios of line intensities for dipole allowed transitions between various energy levels. This approach can be applied to impurities other than the lithiumlike ions and is especially useful for diagnostics of systems where nonintrusive spectroscopic techniques must be used.

Kunc, Joseph A.↗

Cosmic-ray transport in accelerating flows

The quasi-linear transport equation of energetic charged particles is derived, including scattering by Alfven waves propagating parallel and antiparallel to a uniform magnetic field, losses by synchrotron radiation, and acceleration of the cold background fluid supporting the waves. As in the comoving frame equations of radiative transfer, scattering, and loss terms are evaluated in the frame in which the fluid is locally at rest. The diffusion approximation is applied to the resulting equation, yielding a transport equation for the isotropic part of the distribution function and a first-order approximation to the anisotropy. The conditions under which this approximation may be applied are discussed. In the nonrelativistic regime, the standard equation of the diffusion approximation is recovered. As contributions previously investigated, new effects appear from the spatial and temporal variations of the fluid velocity as well as the synchrotron radiation terms. These are discussed in detail.

Kirk, John G.↗

Direct finite element solution on an optical laboratory matrix-vector processor

The first optical laboratory system results employing a direct LU decomposition solution of a system of linear algebraic equations are presented for a finite element problem solution. This also represents the first laboratory demonstration of the use of sign-magnitude negative number representation as well as new bit partitioning techniques to increase the accuracy of an optical encoded processor beyond the number of bit channels available.

Casasent, David↗