SEARCH · Engineering Papers
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.
Airplane identification, equation error method, maximum likelihood method, parameter estimation in frequency domain, extended Kalman filter, aircraft equations of motion, aerodynamic model equations, criteria for the selection of a parsimonious model, and online aircraft identification are addressed.
Explore the source record for details and available documents.
Simulation of turbulence as it relates to the flight training environment is discussed in general terms. Simulators that put random white noise into the system and simulators that put random motion into the equations of motion are discussed, as are simulators that incorporate pitch and roll moments into vertical turbulence. Wind shear models and simplified models of microburst phenomena are covered.
Various Symbolic Manipulation Programs (SMP) were tested to check the functioning of their commands and suitability under various operating systems. Support systems for SMP were found to be relatively better than the one for MACSYMA. The graphics facilities for MACSYMA do not work as expected under the UNIX operating system. Not all commands for MACSYMA function as described in the manuals. Shape representation is a central issue in computer graphics and computer-aided design. Aside from appearance, there are other application dependent, desirable properties like continuity to certain order, symmetry, axis-independence, and variation-diminishing properties. Several shape representations are studied, which include the Osculatory Method, a Piecewise Cubic Polynomial Method using two different slope estimates, Piecewise Cubic Hermite Form, a method by Harry McLaughlin, and a Piecewise Bezier Method. They are applied to collected physical and chemical data. Relative merits and demerits of these methods are examined. Kinematics of a single link, non-dissipative robot arm is studied using MACSYMA. Lagranian is set-up and Lagrange's equations are derived. From there, Hamiltonian equations of motion are obtained. Equations suggest that bifurcation of solutions can occur, depending upon the value of a single parameter. Using the characteristic function W, the Hamilton-Jacobi equation is derived. It is shown that the H-J equation can be solved in closed form. Analytical solutions to the H-J equation are obtained.
A simple, approximate equation describing the velocity-density relationship (or velocity-altitude) has been derived from the flight of large ballistic coefficient projectiles launched at high speeds. The calculations obtained by using the approximate equation compared well with results for numerical integrations of the exact equations of motion. The flightpath equation was used to parametrically calculate maximum body decelerations and stagnation pressures for initial velocities from 2 to 6 km/s. Expressions were derived for the stagnation-point convective heating rates and total heat loads. The stagnation-point heating was parametrically calculated for a nonablating wall and an ablating carbon surface. Although the heating rates were very high, the pulse decayed quickly. The total nose-region heat shield weight was conservatively estimated to be only about 1 percent of the body mass.
This paper surveys computational issues for large-scale simulation of dynamics and controls of space structures, which involve structural elements that are capable of large combined rigid and flexible motions, accurate and efficient treatment of constraints, robust integration of both translational and large rotational motions in the equations of motion, modular interface with active control synthesis packages, and a capability from wave motions to slowly varying transient responses. A particular feature of the present survey is a partitioned solution (or divide-and-conquer) procedure that can handle the numerical solution of multidisciplinary simulation problems by relying on individual modular solution packages that treat each aspect of simulation requirements.
The Navier-Stokes equations of motion and the equation of continuity are transformed so as to apply to an orthogonal curvilinear coordinate system rotating with a uniform angular velocity about an arbitrary axis in space. A usual simplification of these equations as consistent with the accepted boundary-layer theory and an integration of these equations through the boundary layer result in boundary-layer momentum-integral equations for three-dimensional flows that are applicable to either rotating or nonrotating fluid boundaries. These equations are simplified and an approximate solution in closed integral form is obtained for a generalized boundary-layer momentum-loss thickness and flow deflection at the wall in the turbulent case. A numerical evaluation of this solution carried out for data obtained in a curving nonrotating duct shows a fair quantitative agreement with the measures values. The form in which the equations are presented is readily adaptable to cases of steady, three-dimensional, incompressible boundary-layer flow like that over curved ducts or yawed wings; and it also may be used to describe the boundary-layer flow over various rotating surfaces, thus applying to turbomachinery, propellers, and helicopter blades.
This is my final report for the grant Dusty Plasmas in Planetary Magnetospheres. The funding from this grant supported our research on dusty plasmas to study: a) dust plasma interactions in general plasma environments, and b) dusty plasma processes in planetary magnetospheres (Earth, Jupiter and Saturn). We have developed a general purpose transport code in order to follow the spatial and temporal evolution of dust density distributions in magnetized plasma environments. The code allows the central body to be represented by a multipole expansion of its gravitational and magnetic fields. The density and the temperature of the possibly many-component plasma environment can be pre-defined as a function of coordinates and, if necessary, the time as well. The code simultaneously integrates the equations of motion with the equations describing the charging processes. The charging currents are dependent not only on the instantaneous plasma parameters but on the velocity, as well as on the previous charging history of the dust grains.
A model for pulsed inductive plasma acceleration containing an energy equation to account for the various sources and sinks in such devices is presented. The model consists of a set of circuit equations coupled to an equation of motion and energy equation for the plasma. The latter two equations are obtained for the plasma current sheet by treating it as a one-element finite volume, integrating the equations over that volume, and then matching known terms or quantities already calculated in the model to the resulting current sheet-averaged terms in the equations. Calculations showing the time-evolution of the various sources and sinks in the system are presented to demonstrate the efficacy of the model, with two separate resistivity models employed to show an example of how the plasma transport properties can affect the calculation. While neither resistivity model is fully accurate, the demonstration shows that it is possible within this modeling framework to time-accurately update various plasma parameters.
A model for pulsed inductive plasma acceleration is presented that consists of a set of circuit equations coupled to both a one-dimensional equation of motion and an equation governing the partitioning of energy. The latter two equations are obtained for the plasma current sheet by treating it as a single element of finite volume and integrating the governing equations over that volume. The integrated terms are replaced where necessary by physically-equivalent quantities that are calculated through the solution of other parts of the governing equation set. The model improves upon previous one-dimensional performance models by permitting the time-evolution of the energy and state properties of the plasma, the latter allowing for the tailoring of the model to different gases that may be chosen as propellants. The time evolution of the various energy modes in the system and the associated plasma properties, calculated for argon propellant, are presented to demonstrate the efficacy of the model. The model produces a result where efficiency is maximized at a given value of the electrodynamic scaling term known as the dynamic impedance parameter. Qualitatively and quantitatively, the model compares favorably with performance measured for two separate inductive pulsed plasma thrusters, with disagreements attributable to simplifying assumptions employed in the generation of the model solution.
A seismic action which is characterized by the acceleration function of the seismic motion is analyzed. The equations of motion in the direction of the seismic action are determined for a completely rigid ground and for a ground with linear deformability characteristics (rotation deformability). The equation for the rotation motion is derived by considering elastic and viscoelastic properties of the foundation ground. The differential equations are written by treating the mechanical model of a many-storied structure as a fixed vertical beam with a uniformly distributed mass as well as with concentrated masses on each floor.
Classical equation of motion of extended monopole charged particle without preacceleration or runaway solutions
The Lagrange characteristics method and the constants of motion method are applied to various collisional kinetic equations. The relationship between the characteristics and the particle equations of motion is discussed. A special solution of Chandrasekhar equation is given.
A computational procedure suitable for the solution of equations of motion for multibody systems is presented. The present procedure adopts a differential partitioning of the translational motions and the rotational motions. The translational equations of motion are then treated by either a conventional explicit or an implicit direct integration method. A principle feature of this procedure is a nonlinearly implicit algorithm for updating rotations via the Euler four-parameter representation. This procedure is applied to the rolling of a sphere through a specific trajectory, which shows that it yields robust solutions.
A computational procedure suitable for the solution of equations of motion for multibody systems is presented. The present procedure adopts a differential partitioning of the translational motions and the rotational motions. The translational equations of motion are then treated by either a conventional explicit or implicit direct integration method. A principal feature of this procedure is a nonlinearly implicit algorithm for updating rotations via the Euler four-parameter representation. This procedure is applied to the rolling of a sphere through a specific trajectory, which shows that it yields robust solutions.
The presence of tip stores influences both aerodynamic and aeroelastic performances of wings. Such effects are more pronounced in the transonic regime. In this study, transonic aeroelasticity of wings with tip stores is studied for the first time by a theoretical method using the unsteady-small disturbance transonic aerodynamic equations coupled with modal structural equations of motion. The aerodynamic and structural equations of motion are simultaneously integrated by a time-accurate numerical scheme. To validate the tip store simulation, aeroelastic computations are made for a typical rectangular wing with a tip store and results are correlated with available wind tunnel data for the corresponding wing without a tip store at various flight conditions. Aeroelastic computations are also made for a typical fighter wing with a tip store. Present computations show that it is important to account for the aerodynamics of the tip store, particularly in the transonic regime where the tip store can make the wing aeroelastically unstable.
Transonic aeroelasticity of wings with active control surfaces is studied by using the unsteady-small disturbance transonic aerodynamic equations coupled with modal structural equations of motion. The aerodynamic and structural equations of motion are simultaneously integrated by a time-accurate numerical scheme. A procedure of synthesizing active controls with unsteady transonics is presented. Flutter suppression in the transonic regime using active controls is demonstrated for a rectangular wing. Characteristics of a selected control law in the transonic regime are studied. The results from this study are useful in the design of active control systems in the transonic regime.
A full set of equations of motion for the twin lift system, linearized about a hover trim condition were derived and presented. It is shown that this full set of equations of motion decouples into simpler sets of equations of motion if the aerodynamic coupling derivatives of the helicopters are neglected. One of these decoupled sets of equations of motion (referred to as the planar set) was studied at length. The other decoupled set (referred to as the nonplanar set) is examined here. It is shown that when the geometric configuration of the twin-lift is symmetric that a further decoupling is possible into antisymmetric and symmetric sets of equations. One set of these reduced equations of motion referred to as the antisymmetric set is directly equivalent to the longitudinal motion of a single helicopeter with a sling load. The second set, referred to as the symmetric set corresponds to the rotation of the entire system without load motion. As in the case of the planar symmetric motion, the location of the tether attachment point influences the stability of the nonplanar symmetric mode. The trend is opposite however in the nonplanar case in that an attachment point below the helicopter center of gravity gives a favorable effect on the stability. The effect is not as strong however as the unfavorable effect on the planar symmetric mode.