Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “EQUATIONS OF MOTION”

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,081 records · Page 60
All sourcesNASA NTRSDOE OSTI
NASA NTRS1988

Control strategies

Techniques for improving the performance of robotic-manipulator path-control systems comprising independent SISO feedback controllers for each joint are discussed and illustrated with block diagrams, reviewing the results of recent analytical investigations. Topics examined include the servo design for a single link, the equations of motion for manipulators, SISO servo design for multiple links, inverse methods, pole placement with compensation of the gravity terms, linear state-feedback control based on the perturbation equations, and adaptive control methods. Consideration is given to variable-structure systems, suboptimal controllers, and the optimal-design-strategy approach.

Wang, J. C.↗
NASA NTRS1989

Multi-flexible-body dynamics capturing motion-induced stiffness

A multi-flexible-body dynamics formulation incorporating a recently developed theory for capturing motion induced stiffness for a arbitrary structure undergoing large rotation and translation accompanied by small vibrations is presented. In essence, the method consists of correcting prematurely linearized dynamical equations for an arbitrary flexible body with generalized active forces due to geometric stiffness corresponding to a system of twelve inertia forces and nine inertia couples distributed over the body. Equations of motion are derived by means of Kane's method. A useful feature of the formulation is its treatment of prescribed motions and interaction forces. Results of simulations of motions of three flexible spacecraft, involving stiffening during spinup motion, dynamic buckling, and a repositioning maneuver, demonstrate the validity and generality of the theory.

Banerjee, Arun K.↗
NASA NTRS1991

Simulation of cryogenic liquid flows with vapor bubbles

Liquid flows in rocket engine components (such as bearings, seals, and pumps) often involve the formation of vapor bubbles due to local superheating of the fluid (either boiling or cavitation). Under the present effort, an analysis has been developed for liquid flows with vapor bubbles, based on a combined Eulerian-Lagrangian technique, in which the continuous (liquid) phase is treated by solving a system of Eulerian conservation equations, while the discrete (vapor bubble) phase is dealt with by integrating Lagrangian equations of motion in computational coordinates. Vapor bubbles of changing size can be accommodated easily by this analysis, and models for the simulation of bubble formation, growth, and motion have been included. The effect of bubble motion and other bubble processes on the continuous (liquid) phase has been accounted for by appropriate bubble mass, momentum, and energy interchange source terms in the Eulerian conservation equations. To demonstrate the viability of the resulting procedure, the cavitating flow of liquid oxygen through a simplified model of a labyrinth seal has been successfully calculated.

De Jong, Frederik J.↗
NASA NTRS1937

Motion of the two-control airplane in rectilinear flight after initial disturbances with introduction of controls following an exponential law

An airplane in steady rectilinear flight was assumed to experience an initial disturbance in rolling or yawing velocity. The equations of motion were solved to see if it was possible to hasten recovery of a stable airplane or to secure recovery of an unstable airplane by the application of a single lateral control following an exponential law. The sample computations indicate that, for initial disturbances complex in character, it would be difficult to secure correlation with any type of exponential control. The possibility is visualized that the two-control operation may seriously impair the ability to hasten recovery or counteract instability.

Klemin, Alexander↗
NASA NTRS1995

An Analytic Approximation to Very High Specific Impulse and Specific Power Interplanetary Space Mission Analysis

A simple, analytic approximation is derived to calculate trip time and performance for propulsion systems of very high specific impulse (50,000 to 200,000 seconds) and very high specific power (10 to 1000 kW/kg) for human interplanetary space missions. The approach assumed field-free space, constant thrust/constant specific power, and near straight line (radial) trajectories between the planets. Closed form, one dimensional equations of motion for two-burn rendezvous and four-burn round trip missions are derived as a function of specific impulse, specific power, and propellant mass ratio. The equations are coupled to an optimizing parameter that maximizes performance and minimizes trip time. Data generated for hypothetical one-way and round trip human missions to Jupiter were found to be within 1% and 6% accuracy of integrated solutions respectively, verifying that for these systems, credible analysis does not require computationally intensive numerical techniques.

Williams, Craig Hamilton↗
NASA NTRS2004

Formation Control for the MAXIM Mission

Over the next twenty years, a wave of change is occurring in the space-based scientific remote sensing community. While the fundamental limits in the spatial and angular resolution achievable in spacecraft have been reached, based on today s technology, an expansive new technology base has appeared over the past decade in the area of Distributed Space Systems (DSS). A key subset of the DSS technology area is that which covers precision formation flying of space vehicles. Through precision formation flying, the baselines, previously defined by the largest monolithic structure which could fit in the largest launch vehicle fairing, are now virtually unlimited. Several missions including the Micro-Arcsecond X-ray Imaging Mission (MAXIM), and the Stellar Imager will drive the formation flying challenges to achieve unprecedented baselines for high resolution, extended-scene, interferometry in the ultraviolet and X-ray regimes. This paper focuses on establishing the feasibility for the formation control of the MAXIM mission. MAXIM formation flying requirements are on the order of microns, while Stellar Imager mission requirements are on the order of nanometers. This paper specifically addresses: (1) high-level science requirements for these missions and how they evolve into engineering requirements; and (2) the development of linearized equations of relative motion for a formation operating in an n-body gravitational field. Linearized equations of motion provide the ground work for linear formation control designs.

Luquette, Richard J.↗
DOE OSTI2024

A Particle-in-Cell Method for Plasmas with a Generalized Momentum Formulation, Part II: Enforcing the Lorenz Gauge Condition

In a previous paper Christlieb et al. (A particle-in-cell method for plasmas with a generalized momentum formulation, part I: Model formulation, 2024), we developed a new particle-in-cell (PIC) method for the relativistic Vlasov–Maxwell system in which the electromagnetic fields and the equations of motion for the particles were cast in terms of scalar and vector potentials through a Hamiltonian formulation. This new method evolved the potentials under the Lorenz gauge using integral equation methods. New methods to construct spatial derivatives of the potentials that converge at the same rates as the fields were also presented. The new particle method was compared against standard explicit discretizations, including the well-known FDTD-PIC method, for a range of applications involving sheaths and particle beams. Here, this paper extends this new class of methods by focusing on the enforcement the Lorenz gauge condition in both exact and approximate forms using co-located meshes. A time-consistency property of the proposed field solver for the vector potential form of Maxwell’s equations is established, which is shown to preserve the equivalence between the semi-discrete Lorenz gauge condition and the analogous semi-discrete continuity equation. Using this property, we present three methods to enforce a semi-discrete gauge condition. The first method introduces an update for the continuity equation that is consistent with the discretization of the Lorenz gauge condition. Both the finite difference and spectral implementations satisfy this discrete gauge condition to machine precision. The second approach we propose enforces a semi-discrete continuity equation using the boundary integral solution to the field equations. The potential benefit of this approach is that it eliminates spatial derivatives that appear on the particle data, namely the current density, which is often calculated by linear combinations of low-order spline basis functions. This method is ideally suited to boundary integral equation methods that invert multi-dimensional operators without dimensional splitting techniques and will be the subject of future work. The third approach introduces a gauge correcting method that makes direct use of the gauge condition to modify the scalar potential and uses local maps for both the charge and current densities. This results in a gauge error, as the maps do not enforce the continuity equation. The vector potential coming from the current density is taken to be exact, and using the Lorenz gauge, we compute a correction to the scalar potential that makes the two potentials satisfy the gauge condition. This method also enforces the gauge condition to machine precision. We demonstrate two of the proposed methods in the context of periodic domains. Problems defined on bounded domains, including those with complex geometric features remain an ongoing effort. However, this work shows that it is possible to design computationally efficient methods that can effectively enforce the Lorenz gauge condition in a non-staggered PIC formulation.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1978

Interaction of a strong blast wave with a free surface

When a point source explosion is initiated at the ocean surface, the shock propagated into the water is reflected at the surface as a centered expansion wave. The solution in the neighborhood of the interaction point is obtained by writing the equations of motion in the appropriate similarity variables and then changing the independent variables to polar coordinates based at the interaction point. From the zero-order solution of the resulting equations the slopes of boundaries at the interaction point are obtained. A first-order perturbation of this solution provides more accurate representation of the flow variables and the curvature of the shock surface near the interaction point.

Falade, A.↗
NASA NTRS1980

Longitudinal stability and control in wind shear with energy height rate feedback

The longitudinal linearized equations of motion in wind shear were derived for the NASA Terminal Configured Vehicle, a modified Boeing 737 airplane. In addition to the apparent acceleration terms resulting from wind shear, the equations included altitude dependent stability derivatives. A linear analysis of these equations indicates a first order divergence type of instability due to wind shear in which head wind decreased with altitude. Furthermore, this instability cannot be stabilized by attitude control alone. However, attitude control used in combination with an addition feedback loop which consisted of the energy height rate feedback to the throttle proved to be effective in suppressing instability due to wind shear. A brief piloted, real time, nonlinear simulation indicated the desirability of using a display based on the rate of change of energy height rate and of commanded thrust.

Gera, J.↗
NASA NTRS1981

Two-dimensional integrating matrices on rectangular grids

The use of integrating matrices in solving differential equations associated with rotating beam configurations is examined. In vibration problems, by expressing the equations of motion of the beam in matrix notation, utilizing the integrating matrix as an operator, and applying the boundary conditions, the spatial dependence is removed from the governing partial differential equations and the resulting ordinary differential equations can be cast into standard eigenvalue form. Integrating matrices are derived based on two dimensional rectangular grids with arbitrary grid spacings allowed in one direction. The derivation of higher dimensional integrating matrices is the initial step in the generalization of the integrating matrix methodology to vibration and stability problems involving plates and shells.

Lakin, W. D.↗
NASA NTRS1982

Fluid dynamical description of relativistic nuclear collisions

On the basis of both a conventional relativistic nuclear fluid dynamic model and a two fluid generalization that takes into account the interpenetration of the target and projectile upon contact, collisions between heavy nuclei moving at relativistic speeds are calculated. This is done by solving the relevant equations of motion numerically in three spatial dimensions by use of particle in cell finite difference computing techniques. The effect of incorporating a density isomer, or quasistable state, in the nuclear equation of state at three times normal nuclear density, and the effect of doubling the nuclear compressibility coefficient are studied. For the reaction 20Ne + 238U at a laboratory bombarding energy per nucleon of 393 MeV, the calculated distributions in energy and angle of outgoing charged particles are compared with recent experimental data both integrated over all impact parameters and for nearly central collisions.

Nix, J. R.↗
NASA NTRS1984

Description of quantum noise by a Langevin equation

General features of the quantum noise problem expressed as the equations of motion for a particle coupled to a set of oscillators are investigated analytically. Account is taken of the properties of the companion oscillators by formulating quantum statistical correlation Langevin equations (QSLE). The frequency of the oscillators is then retained as a natural cut-off for the quantum noise. The QSLE is further extended to encompass the particle trajectory and is bounded by initial and final states of the oscillator. The states are expressed as the probability of existence at the moment of particle collision that takes the oscillator into a final state. Two noise sources then exist: a statistical uncertainty of the initial state and the quantum dynamical uncertainty associated with a transition from the initial to final state. Feynman's path-integral formulation is used to characterize the functional of the particle trajectory, which slows the particle. It is shown that the energy loss may be attributed to friction, which satisfies energy conservation laws.

Metiu, H.↗
NASA NTRS2009

Mean Element Propagations Using Numerical Averaging

The long-term evolution characteristics (and stability) of an orbit are best characterized using a mean element propagation of the perturbed two body variational equations of motion. The averaging process eliminates short period terms leaving only secular and long period effects. In this study, a non-traditional approach is taken that averages the variational equations using adaptive numerical techniques and then numerically integrating the resulting EOMs. Doing this avoids the Fourier series expansions and truncations required by the traditional analytic methods. The resultant numerical techniques can be easily adapted to propagations at most solar system bodies.

orbit propagations↗
NASA NTRS2016

Simulation of an Electromechanical Spin Motor System of a Control Moment Gyroscope

A two-phase brushless DC motor (BDCM) with pulse-width modulated (PWM) voltage drive is simulated to control the flywheel speed of a control moment gyroscope (CMG). An overview of a double-gimballed control moment gyroscope (DGCMG) assembly is presented along with the CMG torque effects on the spacecraft. The operating principles of a two-phase brushless DC motor are presented and the system's electro-mechanical equations of motion are developed for the root-mean-square (RMS) currents and wheel speed. It is shown that the system is an extremely "stiff" set of first-order equations for which an implicit Euler integrator is required for a stable solution. An adaptive proportional voltage controller is presented which adjusts the PWM voltages depending on several control modes for speed, current, and torque. The simulation results illustrate the interaction between the electrical system and the load dynamics and how these influence the overall performance of the system. As will be shown, the CMG spin motor model can directly provide electrical power use and thermal power output to spacecraft subsystems for effective (average) calculations of CMG power consumption.

Inampudi, Ravi↗
NASA NTRS1989

Determination of joint drives for stable end-effector motion in flexible robotic systems

The prescribed tasks in high speed robotic systems are severely deteriorated because of their manipulator dynamic deflections. On the other hand conventional dynamic modeling techniques fail to reveal appropriate control forces in flexible systems. The conventional dynamic equations of motion for systems subject to kinematical constraints are modified by a new concept of control force representation. The directions of the control forces are selected such that they correspond to the joint degrees of freedom. Then the joint control forces and torques that yield unperturbed prescribed motions are solved simultaneously with the system motion. A flexible manipulator is presented to illustrate the methods proposed.

Ider, Sitki Kemal↗
NASA NTRS1971

Feasibility study of scanning celestial Attitude System (SCADS) for Earth Resources Technology Satellite (ERTS)

The feasibility of using the Scanning Celestial Attitude Determination System (SCADS) during Earth Resources Technology Satellite (ERTS) missions to compute an accurate spacecraft attitude by use of stellar measurements is considered. The spacecraft is local-vertical-stabilized. A heuristic discussion of the SCADS concept is first given. Two concepts are introduced: a passive system which contains no moving parts, and an active system in which the reticle is caused to rotate about the sensor's axis. A quite complete development of the equations of attitude motions is then given. These equations are used to generate the true attitude which in turn is used to compute the transit times of detectable stars and to determine the errors associated with the SCADS attitude. A more complete discussion of the analytical foundation of SCADS concept and its use for the geometries particular to this study, as well as salient design parameters for the passive and active systems are included.

Source record↗
NASA NTRS1975

Numerical implications of stabilization by the use of integrals

Liapunov or energy restraint methods for dynamic stabilization in two body motion perturbation problems are considered. Results of computerized orbital stabilization estimates show that the application of energy restraint prevents the occurrence of consistent timing errors in the stepwise integration of equations of motion for a nearly circular orbit.

Beaudet, P. R.↗