Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical integration”

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 901 records · Page 50

Modeling the Arecibo nighttime F2 layer. I - Overhead properties

The paper extends the servo model of Rishbeth et al. (1978) and applies it to fit the height and plasma density of the nighttime F2 layer, as measured from the Arecibo Observatory during solar maximum by Burnside (1984). The model equations are integrated numerically to fit the observed peak height and density. The model adequately reproduces the observed behavior of the F2 layer. The additional terms in the extended servo-height equation affect the peak height computation by no more than 10 percent on average. The applied current is generally eastward though it becomes westward during the postmidnight collapse or descent of the layer. Differences between model and optically measured meridional wind speeds appear to be related to the presence of large vertical shears in the wind. The observed peak density can be reproduced to within 20 to 40 percent accuracy. Variable plasmaspheric fluxes of the order of 10 exp 13/sq m s contribute to the maintenance and variability of the nighttime peak density.

Melendez-Alvira, D. J.↗

Hypersonic flutter of a curved shallow panel with aerodynamic heating

The general equations describing the nonlinear fluttering oscillations of shallow, curved, heated orthotropic panels have been derived. The formulation takes into account the location of the panel on the surface of a generic hypersonic vehicle, when calculating the aerodynamic loads. It is also shown that third order piston theory produces unsteady aerodynamic loading which is in close agreement with that based upon direct solution of the Euler equations. Results, for simply supported panels, are obtained using Galerkin's method combined with direct numerical integration in time to compute stable limit cycle amplitudes. These results illustrate the sensitivity of the aeroelastic behavior to the unsteady aerodynamic assumptions, temperature, orthotropicity and flow orientation.

Bein, T.↗

Radii effect on the translation spring constant of force transducer beams

Multi-component strain-gage force transducer design requires the designer to determine the spring constant of the numerous beams or flexures incorporated in the transducer. The classical beam deflection formulae that are used in calculating these spring constants typically assume that the beam has a uniform moment of inertia along the entire beam length. In practice all beams have a radius at the end where the beam interfaces with the shoulder of the transducer, and on short beams in particular this increases the beam spring constant considerably. A Basic computer program utilizing numerical integration is presented to determine this effect.

Scott, C. E.↗

ASDA - Advanced Suit Design Analyzer computer program

An ASDA model developed to evaluate the heat and mass transfer characteristics of advanced pressurized suit design concepts for low pressure or vacuum planetary applications is presented. The model is based on a generalized 3-layer suit that uses the Systems Integrated Numerical Differencing Analyzer '85 in conjunction with a 41-node FORTRAN routine. The latter simulates the transient heat transfer and respiratory processes of a human body in a suited environment. The user options for the suit encompass a liquid cooled garment, a removable jacket, a CO2/H2O permeable layer, and a phase change layer.

Bue, Grant C.↗

Stability of outer planetary orbits around binary stars - A comparison of Hill's and Laplace's stability criteria

A comparison is made between the stability criteria of Hill and that of Laplace to determine the stability of outer planetary orbits encircling binary stars. The restricted, analytically determined results of Hill's method by Szebehely and coworkers and the general, numerically integrated results of Laplace's method by Graziani and Black (1981) are compared for varying values of the mass parameter mu. For mu = 0 to 0.15, the closest orbit (lower limit of radius) an outer planet in a binary system can have and still remain stable is determined by Hill's stability criterion. For mu greater than 0.15, the critical radius is determined by Laplace's stability criterion. It appears that the Graziani-Black stability criterion describes the critical orbit within a few percent for all values of mu.

Kubala, A.↗

Pressure dither in venting analyses

Venting analyses simulate the time-dependent flow of air or another gas from one compartment to another, or to an ambient profile. The flow through the vent connecting the compartments is assumed to be one-dimensional, quasi-steady, and isentropic with an empirical discharge coefficient included to account for irreversible flow losses. When a venting analysis is performed numerically, an upper limit on the accuracy exists due to the effects of the numerical integration method. Because this manifests itself in a first-order method by a rapid oscillation about the true answer, this phenomenon is referred to as pressure dither. This paper explores the origin of the phenomenon and gives simple formulae for predicting when it will occur and the greatest achievable accuracy for a given time step.

Fay, John F.↗

A COMPACT model reduction methodology for articulated, multi-flexible bodies structures

To simulate the dynamical motion of articulated, multiflexible body structures, one can use multibody simulation packages such as DISCOS. To this end, one need appropriate reduced-order models for all the flexible components involved. The component modes projection and assembly model reduction methodology is one way to construct these reduced-order component models. However, the generated reduced-order system models typically contain high frequency extraneous modes. The presence of high frequency extraneous modes in a reduced-order system model forces us to use a small stepsize in the numerical integration of that model, resulting in an increase in the elapsed time of the integration process. To eliminate these high-frequency modes, we first identify and discard the component modes that 'contributed' to the generations of these modes. A sensitivity analysis is then used to determine how the remaining component modes should be perturbed in order to still closely capture the selected system modes. The effectiveness of the proposed methodology has been successfully verified using a finite-element model of the Galileo spacecraft.

Lee, Allan Y.↗

Correcting for diffraction in the far-infrared reflectance measurement of rough surfaces

The manner in which diffraction can increase the reflectance measured from very rough surfaces is shown by observations of diffracted light within the instrument profile of a far-IR reflectometer system. A correction to the calibration signal based on numerical integration of the diffracted part of the instrument profile is described. Diffraction correction factors as large as 2.94 have been found with small optics at long wavelength (630 microns). The effect of diffraction on diffuse reflectance measurements of a very rough perfect reflector is shown at wavelengths from 56 to 200 microns.

Smith, Sheldon M.↗

Dynamic analysis and control of lightweight manipulators with flexible parallel link mechanisms

The flexible parallel link mechanism is designed for increased rigidity to sustain the buckling when it carries a heavy payload. Compared to a one link flexible manipulator, a two link flexible manipulator, especially the flexible parallel mechanism, has more complicated characteristics in dynamics and control. The objective of this research is the theoretical analysis and the experimental verification of dynamics and control of a two link flexible manipulator with a flexible parallel link mechanism. Nonlinear equations of motion of the lightweight manipulator are derived by the Lagrangian method in symbolic form to better understand the structure of the dynamic model. A manipulator with a flexible parallel link mechanism is a constrained dynamic system whose equations are sensitive to numerical integration error. This constrained system is solved using singular value decomposition of the constraint Jacobian matrix. The discrepancies between the analytical model and the experiment are explained using a simplified and a detailed finite element model. The step response of the analytical model and the TREETOPS model match each other well. The nonlinear dynamics is studied using a sinusoidal excitation. The actuator dynamic effect on a flexible robot was investigated. The effects are explained by the root loci and the Bode plot theoretically and experimentally. For the base performance for the advanced control scheme, a simple decoupled feedback scheme is applied.

Lee, Jeh Won↗

Landing-gear impact

Report deals with the impact forces in landing gears. Both the landing impact and the taxiing impact have been considered, but drag forces have so far been excluded. The differential equations are developed and their numerical integration is shown, considering the nonlinear properties of the oleo shock strut. A way is shown for determining the dimensions of the metering pin from a given load-time diagram. A review of German literature on landing-gear impact is also presented.

Flugge, W↗

Determination of elastic stresses in gas-turbine disks

A method is presented for the calculation of elastic stresses in symmetrical disks typical of those of a high-temperature gas turbine. The method is essentially a finite-difference solution of the equilibrium and compatibility equations for elastic stresses in a symmetrical disk. Account can be taken of point-to-point variations in disk thickness, in temperature, in elastic modulus, in coefficient of thermal expansion, in material density, and in Poisson's ratio. No numerical integration or trial-and-error procedures are involved and the computations can be performed in rapid and routine fashion by nontechnical computers with little engineering supervision. Checks on problems for which exact mathematical solutions are known indicate that the method yields results of high accuracy. Illustrative examples are presented to show the manner of treating solid disks, disks with central holes, and disks constructed either of a single material or two or more welded materials. The effect of shrink fitting is taken into account by a very simple device.

Manson, S S↗

On flow of electrically conducting fluids over a flat plate in the presence of a transverse magnetic field

The use of a magnetic field to control the motion of electrically conducting fluids is studied. The incompressible boundary-layer solutions are found for flow over a flat plate when the magnetic field is fixed relative to the plate or to the fluid. The equations are integrated numerically for the effect of the transverse magnetic field on the velocity and temperature profiles, and hence, the skin friction and rate of heat transfer. It is concluded that the skin friction and the heat-transfer rate are reduced when the transverse magnetic field is fixed relative to the plate and increased when fixed relative to the fluid. The total drag is increased in all of the areas.

Rossow, Vernon J↗

Optimal design of control systems for prevention of aircraft flight departure

Control system design for general nonlinear flight dynamic models is considered through numerical simulation. The design is accomplished through a numerical optimizer coupled with analysis of flight dynamic equations. In the analysis, the general nonlinear flight dynamic equations with nonlinear aerodynamics are numerically integrated; and the dynamic characteristics needed in the optimization process are then identified from the dynamic response. To assure a reasonable solution, the initial input values of the design variables are estimated through a sensitivity analysis. To demonstrate the method, the pitch departure prevention for an F-16 configuration is demonstrated.

Lan, C. Edward↗

Spaceborne autonomous multiprocessor systems

The goal of this task is to provide technology for the specification and integration of advanced processors into the Space Station Freedom data management system environment through computer performance measurement tools, simulators, and an extended testbed facility. The approach focuses on five categories: (1) user requirements--determine the suitability of existing computer technologies and systems for real-time requirements of NASA missions; (2) system performance analysis--characterize the effects of languages, architectures, and commercially available hardware on real-time benchmarks; (3) system architecture--expand NASA's capability to solve problems with integrated numeric and symbolic requirements using advanced multiprocessor architectures; (4) parallel Ada technology--extend Ada software technology to utilize parallel architectures more efficiently; and (5) testbed--extend in-house testbed to support system performance and system analysis studies.

Fernquist, Alan↗

Precision GPS ephemerides and baselines

The emphasis of this grant was focused on precision ephemerides for the Global Positioning System (GPS) satellites for geodynamics applications. During the period of this grant, major activities were in the areas of thermal force modeling, numerical integration accuracy improvement for eclipsing satellites, analysis of GIG '91 campaign data, and the Southwest Pacific campaign data analysis.

Source record↗

The long-term dynamical behavior of short-period comets

The orbits of the known short-period comets under the influence of the Sun and all the planets except Mercury and Pluto are numerically integrated. The calculation was undertaken in order to determine the dynamical lifetimes for these objects as well as explaining the current orbital element distribution. It is found that a comet can move between Jupiter-family and Halley-family comets several times in its dynamical lifetime. The median lifetime of the known short-period comets from the time they are first injected into a short-period comet orbit to ultimate ejection is approximately 50,000 years. The very flat inclination distribution of Jupiter-family comets is observed to become more distended as it ages. The only possible explanation for the observed flat distribution is that the comets become extinct before their inclination distribution can change significantly. It is shown that the anomalous concentration of the argument of perihelion of Jupiter-family comets near 0 and 180 deg is a direct result of their aphelion distance being close to 5.2AU and the comet being recently perturbed onto a Jupiter-family orbit. Also the concentration of their aphelion near Jupiter's orbit is a result of the conservation of the Tisserand invariant during the capture process.

Levison, Harold F.↗

Orbital simulations of satellite escape/capture and the origin of satellites such as Triton

We investigate satellite escape/capture in the context of the restricted, circular three body problem as applied to the Sun, Neptune, and Triton. We have computed a large number of coplanar prograde and retrograde orbital simulations over a range of initial distances and velocities. The satellite starts at superior conjunction within approximately 2 Hill radii of Neptune and has a velocity orthogonal to the Sun-planet line. Orbits with these initial conditions can be reflected with respect to time, so an escape is simply the reverse of a capture. We numerically integrate the equations of motion to compute the satellite's position until it escapes, collides with Neptune, or after 100 planetary years fails to escape, when computations cease. The initial distance x and velocity v in the restricted problem uniquely define the Jacobi constant C, a conserved energy-like quantity. Plots of the simulation outcomes in the prograde and retrograde C, x phase spaces reveal distinct zones in which temporary satellites approach the planet closely enough that permanent capture can be effected by gas drag with a protoplanetary nebula or by collision with a pre-existing satellite. Single and double close-flybys constitute the most common possible capture orbits. Long term multiple flyby orbits occur near the stability limits between bound and unbound orbits, and are more common among retrograde captures.

Benner, Lance A. M.↗