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.
Equations of motion of shell are solved using Houbolt's numerical procedure with nonlinear terms being moved to right-hand side of equilibrium equations and treated as generalized loads. Program was written in FORTRAN IV for IBM 360 or CDC 6000 series computers.
The derivation of minimum-dimension sets of discrete-coordinate and hybrid-coordinate equations of motion of a system consisting of an arbitrary number of hinge-connected rigid bodies assembled in tree topology is presented. These equations are useful for the simulation of dynamical systems that can be idealized as tree-like arrangements of substructures, with each substructure consisting of either a rigid body or a collection of elastically interconnected rigid bodies restricted to small relative rotations at each connection. Thus, some of the substructures represent elastic bodies subjected to small strains or local deformations, but possibly large gross deformations, in the hybrid formulation, distributed coordinates referred to herein as large-deformation modal coordinates, are used for the deformations of these substructures. The equations are in a form suitable for incorporation into one or more computer programs to be used as multipurpose tools in the simulation of spacecraft and other complex electromechanical systems.
The gravitational collapse of a rotating cloud or vortex is analyzed by expanding the dependent variables in the equations of motion in two-dimensional Taylor series in the space variables. It is shown that the gravitation and rotation terms in the equations are of first order in the space variables, the pressure gradient terms are of second order, and the turbulent viscosity term is of third order. The presence of a turbulent viscosity insures that the initial rotation is solid-body-like near the origin. The effect of pressure on the collapse process is found to depend on the shape of the initial density disturbance at the origin. Dimensionless collapse times, as well as the evolution of density and velocity, are calculated by solving numerically the system of nonlinear ordinary differential equations resulting from the series expansions. The axial inflow plays an important role and allows collapse to occur even when the rotation is large. An approximate solution of the governing partial differential equations is also given; the equations are used to study the spacial distributions of the density and velocity.
The gravitational collapse of a rotating cloud or vortex is analyzed by expanding the dependent variables in the equations of motion in two-dimensional Taylor series in the space variables. It is shown that the gravitation and rotation terms in the equations are of first order in the space variables, the pressure gradient terms are of second order, and the turbulent viscosity term is of third order. The presence of a turbulent viscosity insures that the initial rotation is solid-body-like near the origin. The effect of pressure on the collapse process is found to depend on the shape of the intial density disturbance at the origin. Dimensionless collapse times, as well as the evolution of density and velocity, are calculated by solving numerically the system of nonlinear ordinary differential equations resulting from the series expansions. The axial inflow plays an important role and allows collapse to occur even when the rotation is large. An approximate solution of the governing partial differential equations is also given, in order to study the spacial distributions of the density and velocity.
The gravitational collapse of a rotating cloud or vortex is analyzed by expanding the dependent variables in the equations of motion in two-dimensional Taylor series in the space variables. It is shown that the gravitational and rotational terms in the equations are of first order in the space variables, the pressure-gradient terms are of second order, and the turbulent-viscosity term is of third order. The presence of turbulent viscosity ensures that the initial rotation is solid-body-like near the origin. The effect of pressure on the collapse process is found to depend on the shape of the initial density disturbance at the origin. Dimensionless collapse times, as well as the evolution of density and velocity, are calculated by solving numerically the system of nonlinear ordinary differential equations resulting from the series expansions. The axial flow is always inward and allows collapse to occur (axially) even when the rotation is large. An approximate solution of the governing partial differential equations is also given in order to study the spatial distributions of the density and velocity.
Lubricants, usually Newtonian fluids, are assumed to experience laminar flow. The basic equations used to describe the flow are the Navier-Stokes equation of motion. The study of hydrodynamic lubrication is, from a mathematical standpoint, the application of a reduced form of these Navier-Stokes equations in association with the continuity equation. The Reynolds equation can also be derived from first principles, provided of course that the same basic assumptions are adopted in each case. Both methods are used in deriving the Reynolds equation, and the assumptions inherent in reducing the Navier-Stokes equations are specified. Because the Reynolds equation contains viscosity and density terms and these properties depend on temperature and pressure, it is often necessary to couple the Reynolds with energy equation. The lubricant properties and the energy equation are presented. Film thickness, a parameter of the Reynolds equation, is a function of the elastic behavior of the bearing surface. The governing elasticity equation is therefore presented.
A design procedure is presented for the installation of simple pendulums on the blades of a helicopter rotor to suppress the root reactions. The procedure consists of a frequency response analysis for a hingeless rotor blade excited by a harmonic variation of spanwise airload distributions during forward flight, as well as a concentrated load at the tip. The structural modeling of the blade provides for elastic degrees of freedom in flap and lead-lag bending plus torsion. Simple flap and lead-lag pendulums are considered individually. Using a rational order scheme, the general nonlinear equations of motion are linearized. A quasi-steady aerodynamic representation is used in the formation of the airloads. The solution of the system equations derives from their representation as a transfer matrix. The results include the effect of pendulum tuning on the minimization of the hub reactions.
Methods of using rotor vacuum whirl data to improve the ability to model helicopter rotors were developed. The work consisted of the formulation of the equations of motion of elastic blades on a hub using a Galerkin method; the development of a general computer program for simulation of these equations; the study and implementation of a procedure for determining physical parameters based on measured data; and the application of a method for computing the normal modes and natural frequencies based on test data.
The corotating magnetospheric convection model is analyzed with the objective of deriving the fundamental wave modes present in the Io torus. The time-dependent equations of motion are linearized for the case of the absence of plasma mass and momentum sources. By separating the magnetosphere-ionosphere coupling equation into two equations and by adopting a specific analytic form for the equilibrium plasma density, a cubic equation is derived for the dispersion relation defining the angular frequencies of the fundamental wave modes. A detailed analysis of the dispersion relation yields three fundamental wave modes. One of the modes is identified as the previously studied inertial interchange (unstable) mode, while the other two are newly derived decaying modes. The results obtained are discussed in the context of plasma transport in the Io torus.
The dynamical equations-of-motion (EOM) for cruising airships require nonconventional terms to account for buoyancy and apparent-mass-effects, but systematic validation of these equations against flight data is not available. Using a candidate set of EOM, three comparisons are made with carefully-measured describing functions derived from frequency-sweep flight tests on the Skyship-500 airship. The first compares the pre-flight predictions to the data; the second compares the 'best-fit' equations to data at each of two airspeeds and the third compared the ability to extrapolate from one condition to another via airship-specific scaling laws. Two transient responses are also compared. The generally good results demonstrate that fairly simple, perturbation equation models are adequate for many types of flight-control analysis and flying quality evaluations of cruising airships.
This paper describes a multibody dynamics algorithm formulated for parallel implementation on multiprocessor computing platforms using the divide-and-conquer approach. The system of interest is a general topology of rigid and elastic articulated bodies with or without loops. The algorithm divides the multibody system into a number of smaller sets of bodies in chain or tree structures, called "branches" at convenient joints called "connection points", and uses an Order-N (O (N)) approach to formulate the dynamics of each branch in terms of the unknown spatial connection forces. The equations of motion for the branches, leaving the connection forces as unknowns, are implemented in separate processors in parallel for computational efficiency, and the equations for all the unknown connection forces are synthesized and solved in one or several processors. The performances of two implementations of this divide-and-conquer algorithm in multiple processors are compared with an existing method implemented on a single processor.
Thermal fluctuations are a fundamental feature of dissipative systems that are essential for understanding physics near the expected critical point of QCD and in small systems. When such fluctuations are modeled naively in relativistic systems, strange features can appear such as negative self-correlation functions. We construct an effective theory for nonlinear stochastic relativistic hydrodynamics that ensure a well-posed mathematical formulation. Using Crooks fluctuation theorem, we derive a symmetry of the effective action that incorporates fluctuations through a suitable free energy functional. For divergence type theories, the action can then be fully specified using a single vector generating current. The equations of motion obtained using this procedure are guaranteed to be flux conservative and symmetric hyperbolic when the dynamics is causal. This ensures that these equations are well-posed (for suitable initial data) and are in a form that can easily be simulated, including with Metropolis techniques.
The density matrix quantum Monte Carlo (DMQMC) set of methods stochastically samples the exact N-body density matrix for interacting electrons at finite temperature. We introduce a simple modification to the interaction picture DMQMC (IP-DMQMC) method that overcomes the limitation of only sampling one inverse temperature point at a time, instead allowing for the sampling of a temperature range within a single calculation, thereby reducing the computational cost. At the target inverse temperature, instead of ending the simulation, we incorporate a change of picture away from the interaction picture. The resulting equations of motion have piecewise functions and use the interaction picture in the first phase of a simulation, followed by the application of the Bloch equation once the target inverse temperature is reached. We find that the performance of this method is similar to or better than the DMQMC and IP-DMQMC algorithms in a variety of molecular test systems.
The concept of the Lagrangian displacement of a balloon is introduced. It is shown that the general balloon response problem is extremely complicated because the wind-forcing functions in the balloon equations of motion are functions of the wind velocity vector and its Eulerian first derivatives evaluated at the location of the balloon. The linear perturbation equations for a spherical balloon are derived by perturbing the components of velocity of the balloon about a terminal velocity state which is in equilibrium with a space-time invariant mean horizontal flow. The atmospheric flow is also perturbed such that the resulting equations can be used to analyze the responses of spherical balloons to three-dimensional time-dependent flows. The wind field is represented in terms of a four-fold Fourier integral that involves three orthogonal wave numbers and a frequency, while the balloon components of velocity are represented as Fourier integrals involving a frequency which, in turn, is a function of the wind field wave numbers and frequency and the unperturbed flow components of velocity.
The results of 18 months of investigations are reported. A movable mass control system to convert the tumbling motion of a spacecraft into simple spin was studied along with the optimization techniques for generating displacement profiles for a tumbling asymmetrical body. Equations of motion are discussed for two asymmetrical vehicles with flexible beams and one spacecraft with flexible solar arrays. The characteristics which allow reasonable safety and reliability in bailout are also discussed.
Coupled flap-lag-torsion equations of motion of an isolated horizontal axis wind turbine (HAWT) blade have been formulated. The analysis neglects blade-tower coupling. The final nonlinear equations have periodic coefficients. A new and convenient method of generating an appropriate time-dependent equilibrium position, required for the stability analysis, has been implemented and found to be computationally efficient. Steady-state response and stability boundaries for an existing (typical) HAWT blade are presented. Such stability boundaries have never been published in the literature. The results show that the isolated blade under study is basically stable. The tower shadow (wake) has a considerable effect on the out-of-plane response but leaves blade stability unchanged. Nonlinear terms can significantly affect linearized stability boundaries; however, they have a negligible effect on response, thus implying that a time-dependent equilibrium position (or steady-state response), based completely on the linear system, is appropriate for the type of HAWT blades under study.
Selected space structure configurations were reviewed in order to define dynamic analysis problems associated with variable geometry. The dynamics of a beam being constructed from a flexible base and the relocation of the completed beam by rotating the remote manipulator system about the shoulder joint were selected. Equations of motion were formulated in physical coordinates for both of these problems, and FORTRAN programs were developed to generate solutions by numerically integrating the equations. These solutions served as a standard of comparison to gauge the accuracy of approximate solution techniques that were developed and studied. Good control was achieved in both problems. Unstable control system coupling with the system flexibility did not occur. An approximate method was developed for each problem to enable the analyst to investigate variable geometry effects during a short time span using standard fixed geometry programs such as NASTRAN. The average angle and average length techniques are discussed.
The dynamic equations of motion for a spinning airplane are simplified to facilitate the analysis of spin entry and recovery maneuvers. Solutions of a nonlinear, first order equation for spin rate enable calculation of spin entry and recovery times and control required for recovery. The analysis is applied to a light airplane for which both aerodynamic data, rotory balance wind tunnel tests, and spin flight test data have been obtained. A comparison of predicted and actual transient spin responses is made which supports the validity of the analytical approach taken but which also illustrates the difficulty in obtaining accurate aerodynamic data for spinning airplanes. Certain normalized parameters of the reduced state spin equation suggest criteria for assessing the transient spin characteristics of light airplanes.