Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “characteristic 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 271 records · Page 15

Spacetimes with Killing tensors

The characteristics of the Killing equation and the Killing tensor are discussed. A conformal Killing tensor is of interest inasmuch as it gives rise to a quadratic first integral for null geodesic orbits. The Einstein-Maxwell equations are considered together with the Bianchi identity and the conformal Killing tensor. Two examples for the application of the considered relations are presented, giving attention to the charged Kerr solution and the charged C-metric.

Hughston, L. P.↗

Ray trajectories in a torus: An application of MACSYMA to a complex numerical computation

The study of ray trajectories of plasma waves in a torodial geometry using MACSYMA is an example of how symbolic, numerical, and graphical facilities can be used in concert to accomplish a complex computational goal. Computational features of this study which are of particular significance include: the derivation of code (i.e. writing functions to generate program fragments), the use of array functions to simplify the specification of a numerical iteration scheme, and the graphical presentation of the results. Mathematically, this study originates in the solution of a linear inhomogeneous partial differential equation in 3 dimensions by the method of characteristics. While it is possible to describe this equation compactly by using vector notation, and by specifying the spatial variation of the coefficients in terms of intermediate parameters, the transformation of the equation into a form amenable to solution is very tedious. A MACSYMA program is presented for obtaining description of the rf field structure excited by a waveguide located at the edge of a toroidal plasma confinement device.

Kulp, J. L.↗

Self-adjusting grid methods for one-dimensional hyperbolic conservation laws

The automatic adjustment of a grid which follows the dynamics of the numerical solution of hyperbolic conservation laws is given. The grid motion is determined by averaging the local characteristic velocities of the equations with respect to the amplitudes of the signals. The resulting algorithm is a simple extension of many currently popular Godunov-type methods. Computer codes using one of these methods can be easily modified to add the moving mesh as an option. Numerical examples are given that illustrate the improved accuracy of Godunov's and Roe's methods on a self-adjusting mesh. Previously announced in STAR as N83-15008

Harten, A.↗

Remote sensing of earth terrain

Progress on the investigation of the anisotropy of the terrain media, such as vegetation canopy and sea ice, and the study of the fluctuation-dissipation theorem in conjunction with the application of strong fluctuation theory for passive remote sensing of snowpacks is reported. The Feynman diagrammatic technique is used to derive the Dyson equation for the mean field and the Bethe-Salpeter equation for the correlation or the covariance of the field for electromagnetic wave propagation and scattering in an anisotropic random medium. With the random permittivity expressed in a general form, the bilocal and the nonlinear approximations are employed to solve the Dyson equation and the ladder approximation to the Bethe-Salpeter equation. The mean dyadic Green's function for a two layer anisotropic random medium with arbitrary three dimensional correlation function was investigated with the zeroth-order solutions to the Dyson equation under the four characteristic waves associated with the coherent vector fields propagating in an anisotropic random medium layer, which are the ordinary and extraordinary waves with upward and downward propagating vectors.

Kong, J. A.↗

Magnetized thermal conduction fronts

The evolution of planar thermal conduction fronts in the presence of a dynamically weak, but otherwise self-consistent, magnetic field is considered. The field is assumed to be connected and untangled. In the diffusion limit for the thermal conductivity, these fronts exhibit self-similar behavior, even in the presence of a field. The role of the field is restricted to channeling the heat flux along its lines of force, and it enters into the problem as a dimensionless angle variable. 'Combing' (or opening) of insulating field lines by the evaporative flow is explicitly demonstrated. Unless the field is nearly perpendicular to the front normal in the hot gas, insulating effects are not profound. Self-similarity breaks down if the front becomes saturated, and under certain conditions magnetized saturated conduction fronts cannot propagate: the solution characteristics of the wave equation form caustics. The physical resolution is the advent of two-fluid (nonlocal) heating. Such Coulomb-heated fronts are expected to be relatively rare in typical astrophysical systems. The large-scale effects of a magnetic field on cloud evaporation in the interstellar medium are briefly discussed, and it is suggested that these fields preclude the presence of time-independent evaporative solutions. Thermal interfaces may then continue to evolve until radiative cooling halts their development; large tracts of warm 10,000 K gas may result.

Balbus, S. A.↗

A linear quadratic tracker for Control Moment Gyro based attitude control of the Space Station

The paper discusses a design for an attitude control system for the Space Station which produces fast response, with minimal overshoot and cross-coupling with the use of Control Moment Gyros (CMG). The rigid body equations of motion are linearized and discretized and a Linear Quadratic Regulator (LQR) design and analysis study is performed. The resulting design is then modified such that integral and differential terms are added to the state equations to enhance response characteristics. Methods for reduction of computation time through channelization are discussed as well as the reduction of initial torque requirements.

Kaidy, J. T.↗

Evaluation of Navier-Stokes and Euler solutions for leading-edge separation vortices

Extensive study on the numerical simulation of the vortical flow over a double delta wing is carried out using the thin layer Navier-Stokes and Euler equations. Two important flow characteristics, vortex interaction and vortex breakdown, are successfully simulated. Grid resolution is one of the most important factors associated with the vortex problem. Computations were performed on a series of grids with various levels of refinement, coarse, medium, and fine. Computations using either the coarse or medium grids fail to capture the proper physical phenomena. The computed result using a fine grid shows flow unsteadiness once the vortex breakdown takes place. The C sub L - alpha characteristics are well predicted up to the breakdown angle of attack for all the grid distributions. The Euler solutions show fairly good agreement with the experiment on the C sub L - alpha characteristics. However, other aspects of the solution at each angle of attack, such as the locus of the leading edge separation vortex, are not consistent with the experiment. Even for the fine grid Navier-Stokes computations, further grid resolution is required to obtain good quantitative agreement with the experiment.

Fujii, K.↗

A one-dimensional shock capturing finite element method and multi-dimensional generalizations

Multi-dimensional generalizations of a one-dimensional finite element shock capturing scheme are proposed. A scalar model problem is used to emphasize that 'preferred directions' are important in multi-dimensional applications. Schemes are developed for the two-dimensional Euler equations. One, based upon characteristics, employs the Mach lines and streamlines as preferred directions.

Hughes, T. J. R.↗

Relationship of Thematic Mapper simulator data to leaf area index of temperate coniferous forests

Regional relationships between remote sensing data and the leaf area index (LAI) of coniferous forests were analyzed using data acquired by an Airborne Thematic Mapper. Eighteen coniferous forest stands with a range of projected leaf area index of 0.6-16.1 were sampled from an environmental gradient in moisture and temperature across west-central Oregon. Spectral radiance measurements to account for atmospheric effects were acquired above the canopies from a radiometer mounted on a helicopter. A strong positive relationship was observed between LAI of closed canopy forest stands and the ratio of near-infrared and red spectral bands. A linear regression based on LAI explained 83 percent of the variation in the ratio of the atmospherically corrected bands. A log-linear equation fit the asymptotic characteristic of the relationship better, explaining 91 percent of the variance. The positive relationship is explained by a strong asymptotic inverse relationship between LAI and red radiation and a relatively flat response between LAI and near-infrared radiation.

Peterson, David L.↗

Evaluation of Navier-Stokes and Euler solutions for leading-edge separation vortices

Extensive study on the numerical simulation of the vortical flow over a double delta wing is carried out using the thin layer Navier-Stokes and Euler equations. Two important flow characteristics, vortex interaction and vortex breakdown, are successfully simulated. Grid resolution is one of the most important factors associated with the vortex problem. Computations were performed on a series of grids with various levels of refinement, coarse, medium, and fine. Computations using either the coarse or medium grids fail to capture the proper physical phenomena. The computed result using a fine grid shows flow unsteadiness once the vortex breakdown takes place. The C sub L - alpha characteristics are well predicted up to the breakdown angle of attack for all the grid distributions. The Euler solutions show fairly good agreement with the experiment on the C sub L - alpha characteristics. However, other aspects of the solution at each angle of attack, such as the locus of the leading edge separation vortex, are not consistent with the experiment. Even for the fine grid Navier-Stokes computations, further grid resolution is required to obtain good quantitative agreement with the experiment.

Fujii, K.↗

Performance limitations in parallel processor simulations

A jet-engine model is partitioned and simulated on a parallel processor system consisting of five 8086/8087 floating-point computers. The simulation uses Heun's integration method. A near-optimal parallel simulation (in the sense of minimum execution time) achieves speedup of only 2.13 and efficiency of 42.6 percent, in effect wasting 57.4 percent of the available processing power. A detailed analysis identifies and graphically demonstrates why the system fails to achieve ideal performance (viz., speedup of 5 and efficiency of 100 percent). Inherent characteristics of the problem equations and solution algorithm account for the loss of nearly half of the available processing power. Overheads associated with interprocessor communication and processor synchronization account for only a small fraction of the lost processing power. The effects of these and other factors which limit parallel processor performance are illustrated through real-time timing-analyzer tracers describing the run/idle status of the parallel processors during the simulation.

O'Grady, E. Pearse↗

Global surface-water-induced seasonal variations in the earth's rotation and gravitational field

The effects of seasonal changes in continental surface-water storage on the low-degree gravitational-field coefficients (J), the annual wobble excitation (Psi), and the seasonal length-of-day (LOD) variations are investigated by means of numerical simulations based on compiled meteorological data (Willmott et al., 1985) and satellite snow-load estimates (Chao et al., 1987). The formulation of the model equations and the overall characteristics of the data sets are discussed in detail, and the computation results are presented in tables and graphs. The effect on Psi is found to be relatively small due to longitudinal cancellation, but those on LOD and J are considered significant.

Chao, B. F.↗

Spectral simulation of unsteady compressible flow past a circular cylinder

An unsteady compressible viscous wake flow past a circular cylinder was successfully simulated using spectral methods. A new approach in using the Chebyshev collocation method for periodic problems is introduced. It was further proved that the eigenvalues associated with the differentiation matrix are purely imaginary, reflecting the periodicity of the problem. It was been shown that the solution of a model problem has exponential growth in time if improper boundary conditions are used. A characteristic boundary condition, which is based on the characteristics of the Euler equations of gas dynamics, was derived for the spectral code. The primary vortex shedding frequency computed agrees well with the results in the literature for Mach = 0.4, Re = 80. No secondary frequency is observed in the power spectrum analysis of the pressure data.

Don, Wai-Sun↗

Inverse dynamics of a 3 degree of freedom spatial flexible manipulator

A technique is presented for solving the inverse dynamics and kinematics of 3 degree of freedom spatial flexible manipulator. The proposed method finds the joint torques necessary to produce a specified end effector motion. Since the inverse dynamic problem in elastic manipulators is closely coupled to the inverse kinematic problem, the solution of the first also renders the displacements and rotations at any point of the manipulator, including the joints. Furthermore the formulation is complete in the sense that it includes all the nonlinear terms due to the large rotation of the links. The Timoshenko beam theory is used to model the elastic characteristics, and the resulting equations of motion are discretized using the finite element method. An iterative solution scheme is proposed that relies on local linearization of the problem. The solution of each linearization is carried out in the frequency domain. The performance and capabilities of this technique are tested through simulation analysis. Results show the potential use of this method for the smooth motion control of space telerobots.

Bayo, Eduardo↗

Semi-empirical model for prediction of unsteady forces on an airfoil with application to flutter

A semi-empirical model is described for predicting unsteady aerodynamic forces on arbitrary airfoils under mildly stalled and unstalled conditions. Aerodynamic forces are modeled using second order ordinary differential equations for lift and moment with airfoil motion as the input. This model is simultaneously integrated with structural dynamics equations to determine flutter characteristics for a two degrees-of-freedom system. Results for a number of cases are presented to demonstrate the suitability of this model to predict flutter. Comparison is made to the flutter characteristics determined by a Navier-Stokes solver and also the classical incompressible potential flow theory.

Mahajan, Aparajit J.↗

The 26-day oscillation observed in the satellite sea surface temperature measurements in the equatorial western Indian Ocean

A 26-d oscillation in SST data is observed in the western Indian Ocean, from 52 to 60 deg E and in the vicinity of the equator. The SST data used in this study are obtained from the NOAA 9 satellite and are for the years 1987 and 1988. This fluctuation of SST at a period near 26-d is found to be antisymmetric about the equator and is trapped within the equatorial waveguide latitude; furthermore, the variance decreases at a faster rate toward the equator than poleward. These characteristics are consistent with the latitudinal structure for the mixed Rossby-gravity (or Yanai) waves as predicted from linear wave theory. The temporal variation of this 26-d oscillation is most energetic during the summer season (July to September), with maximum values of 0.4 C and 0.8 C found during August of 1987 and 1988, respectively. This observation agrees with temporal variation of Yanai waves inferred from drifting buoy observations and numerical studies of the Indian Ocean. Thus it is concluded that the Yanai wave is responsible for the 26-d fluctuation observed in the SST data in this region.

Tsai, Pedro T. H.↗

Semi-empirical model for prediction of unsteady forces on an airfoil with application to flutter

A semi-empirical model is described for predicting unsteady aerodynamic forces on arbitrary airfoils under mildly stalled and unstalled conditions. Aerodynamic forces are modeled using second order ordinary differential equations for lift and moment with airfoil motion as the input. This model is simultaneously integrated with structural dynamics equations to determine flutter characteristics for a two degrees-of-freedom system. Results for a number of cases are presented to demonstrate the suitability of this model to predict flutter. Comparison is made to the flutter characteristics determined by a Navier-Stokes solver and also the classical incompressible potential flow theory.

Mahajan, A. J.↗