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 829 records · Page 46

Variations in the modal characteristics of a telescopically deploying beam

The equations of motion for a two-segment deploying telescopic beam are derived through application of Lagrange's equation. The outer tube of the beam is fixed at one end and the inner tube slides freely relative to the fixed segment. The resulting nonlinear, non-autonomous set of equations is linearized and simplified to the standard Euler-Bernoulli partial differential equations for an elastic beam by freezing the deployment process at various stages of deployment, and examining the small amplitude and natural modes of vibration of the resulting configuration. Application of the natural boundary conditions and compatibility of motion relations for the two segments in their common region of overlap leads to a transcendental characteristic equation in the frequency parameter Beta(L). Numerical solution of the equation for the characteristic roots determines the modal frequencies, and the corresponding mode shapes are obtained from the general solution of the Euler-Bernoulli equation tailored to the natural boundary conditions. Sample results of modal frequencies and shapes are presented for various stages of deployment and discussed. It is shown that for all intermediate stages of deployment (between 0 and 100 percent) the spectral distribution is drastically altered by the appearance of regions of very closely spaced modal frequencies. The sources of this modal agglomeration are explored.

Amos, Anthony K.↗

Automatic Differentiation in MetaPhysicL and Its Applications in MOOSE

Efficient solution via Newton’s method of nonlinear systems of equations requires an accurate representation of the Jacobian, corresponding to the derivatives of the component residual equations with respect to the degrees of freedom. In practice these systems of equations often arise from spatial discretization of partial differential equations used to model physical phenomena. These equations may involve domain motion or material equations that are complex functions of the systems’ degrees of freedom. Computing the Jacobian by hand in these situations is arduous and prone to error. Finite difference approximations of the Jacobian or its action are prone to truncation error, especially in multiphysics settings. Symbolic differentiation packages may be used, but often result in an excessive number of terms in realistic model scenarios. An alternative to symbolic and numerical differentiation is automatic differentiation (AD), which propagates derivatives with every elementary operation of a computer program, corresponding to continual application of the chain rule. Automatic differentiation offers the guarantee of an exact Jacobian at a relatively small overhead cost. In this work, we outline the adoption of AD in the Multiphysics Object Oriented Simulation Environment (MOOSE) via the MetaPhysicL package. We describe the application of MOOSE’s AD capability to several sets of physics that were previously infeasible to model via hand-coded or Jacobian-free simulation techniques, including arbitrary Lagrangian-Eulerian and level-set simulations of laser melt pools, phase-field simulations with free energies provided through neural networks, and metallic nuclear fuel simulations that require inner Newton loop calculation of nonlinear material properties.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Quantum quenches of an SO(5) pseudospin reveal Higgs bosons

A controlled dynamical probe measurement of complex order parameter fluctuations of a system may reveal its massive collective excitations. Here, we design dynamical quench protocols to excite independently all ten midgap Higgs bosons in the isotropic Balian-Werthamer state of a spinful p -wave superfluid or superconductor. The analysis is based on microscopic equations of motion of an SO(5) pseudospin, an extension of the usual Bloch equation to a five-dimensional space. Key to these protocols is the realization of quenches that break the rotational symmetry of the kinetic energy and exploit the irreducible representation of the angular momentum J = 2. For perturbative quenches, we find ( nondecaying ) periodic oscillations in time of these Higgs modes. Finally, we present experimental protocols for superconductors (superfluids), with the intention of unveiling the nature of their order parameters.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Variational calculations for resonance oscillations of inhomogeneous plasmas

The electrostatic resonance properties of an inhomogeneous plasma column are treated by the Rayleigh-Ritz method. In contrast to Parker, Nickel & Gould (1964), who carried out an exact computation, the present treatment uses a description of the RF equation of motion and pressure term that allows one to express the system of equations in Euler-Lagrange form. The Rayleigh-Ritz procedure is then applied to the corresponding Lagrangian, to obtain approximate resonance frequences and eigenfunctions. An appropriate set of trial coordinate functions is defined, which leads to frequency and eigenfunction estimates in excellent agreement with the work of Parker et al. (1964).

Peng, Y.-K. M.↗

Flutter of articulated pipes at finite amplitude

The plane motion of an articulated pipe made of two segments is examined and the flow velocity at which flutter manifests itself is sought. The pressure in the reservoir feeding the pipe is kept constant. In contrast to previous works, the flow velocity is not taken as a prescribed parameter of the system but is left to follow the laws of motion. This approach requires a nonlinear formulation of the problem and the equations of motion are solved using Krylov-Bogoliubov's method. A graph of the amplitude of the limit cycles, as a function of the fluid-system mass ratio, is presented and conclusions are drawn as to the necessity of considering nonlinearities in the analysis.

Rousselet, J.↗

Evidence for large-scale velocity features on the sun

Giant velocity features on the solar surface are seen in line-of-sight velocity data from the Mount Wilson magnetograph. Velocity amplitudes around 40 m/s are seen in features with dimensions about 15 deg in latitude and 30-60 deg in longitude. These features are associated with solar activity. The largest such feature accompanied the largest complex of activity in this interval, lived for about 16 rotations, and predated the first magnetic fields in the area by about two rotations. A general pattern of upward motion at the equator is seen, and motions away from the earth at higher latitudes could represent a meridional flow toward the poles of the order of 20 m/s.

Howard, R.↗

Periodic solutions about the collinear Lagrangian solution in the general problem of three bodies

The article describes the solutions near Lagrange's circular collinear configuration in the planar problem of three bodies with three finite masses. The article begins with a detailed review of the properties of Lagrange's collinear solution. Lagrange's quintic equation is derived and several expressions are given for the angular velocity of the rotating frame. The equations of motion are then linearized near the circular collinear solution, and the characteristic equation is also derived in detail. The different types of roots and their corresponding solutions are discussed. The special case of two equal outer masses receives special attention, as well as the special case of two small outer masses. Finally, the fundamental family of periodic solutions is extended by numerical integration all the way up to and past a binary collision orbit. The stability and the bifurcations of this family are briefly enumerated.

Broucke, R.↗

Improved perfect-fluid energy-momentum tensor with spin in Einstein-Cartan space-time

The description of the spin given here is classical in that it is intrinsic but not quantized. The approach in this matter is similar to, for example, the work of Bailey and Israel (1973, 1975, 1979), where the fluid particles, which have intrinsic spin, may be galaxies or clusters of galaxies. The elementary particles of these objects and the 'ferromagnetic alignment' of their quantum spins are not resorted to in order to describe a fluid with spin. Physically this means that the equation of motion for the spin tensor is a modified Fermi-Walker transport equation (Misner et al., 1973), arising as a direct result of the inclusion of spin as an intrinsic variable in the thermodynamic description of the internal energy. The variables in this description are classical variables throughout and are not microscopic fields. An improved perfect-fluid energy-momentum tensor that includes spin and torsion is presented. Use is made of a Lagrangian variational principle based on the tetrad formalism of Halbwach (1960) and the method od constraints of Ray (1972).

Ray, J. R.↗

Normal mode study of the earth's rigid body motions

In this paper it is shown that the earth's rigid body (rb) motions can be represented by an analytical set of eigensolutions to the equation of motion for elastic-gravitational free oscillations. Thus each degree of freedom in the rb motion is associated with a rb normal mode. Cases of both nonrotating and rotating earth models are studied, and it is shown that the rb modes do incorporate neatly into the earth's system of normal modes of free oscillation. The excitation formula for the rb modes are also obtained, based on normal mode theory. Physical implications of the results are summarized and the fundamental differences between rb modes and seismic modes are emphasized. In particular, it is ascertained that the Chandler wobble, being one of the rb modes belonging to the rotating earth, can be studied using the established theory of normal modes.

Chao, B. F.↗

A vectorized algorithm for 3D dynamics of a tethered satellite

Equations of motion characterizing the three dimensional motion of a tethered satellite during the retrieval phase are studied. The mathematical model involves an arbitrary number of point masses connected by weightless cords. Motion occurs in a gravity gradient field. The formulation presented accounts for general functions describing support point motion, rate of tether retrieval, and arbitrary forces applied to the point masses. The matrix oriented program language MATLAB is used to produce an efficient vectorized formulation for computing natural frequencies and mode shapes for small oscillations about the static equilibrium configuration; and for integrating the nonlinear differential equations governing large amplitude motions. An example of time response pertaining to the skip rope effect is investigated.

Wilson, Howard B.↗

Spatial operator algebra framework for multibody system dynamics

The Spatial Operator Algebra framework for the dynamics of general multibody systems is described. The use of a spatial operator-based methodology permits the formulation of the dynamical equations of motion of multibody systems in a concise and systematic way. The dynamical equations of progressively more complex grid multibody systems are developed in an evolutionary manner beginning with a serial chain system, followed by a tree topology system and finally, systems with arbitrary closed loops. Operator factorizations and identities are used to develop novel recursive algorithms for the forward dynamics of systems with closed loops. Extensions required to deal with flexible elements are also discussed.

Rodriguez, G.↗

Dynamics of flare-driven chromospheric condensations

A general model of chromospheric condensation dynamics is developed under the assumption of sudden impulsive-phase flare energy release. The basic equation of motion for the downward-moving condensation is derived from the momentum equation and from general features seen in numerical simulations. This equation is solved for a number of different cases. The 'lifetime' of condensation and the time for the velocity to drop to half its peak value are computed. The peak downflow speed is examined as a function of the flare energy flux, and comparisons are made between downflow speeds predicted by both thick-target 'explosive evaporation' and evaporation driven conductively.

Fisher, George H.↗

Attitude control/momentum management of the Space Station Freedom for large angle torque-equilibrium-attitude configurations

An attitude-control and momentum-management (ACMM) system for the Space Station in a large-angle torque-equilibrium-attitude (TEA) configuration is developed analytically and demonstrated by means of numerical simulations. The equations of motion for a rigid-body Space Station model are outlined; linearized equations for an arbitrary TEA (resulting from misalignment of control and body axes) are derived; the general requirements for an ACMM are summarized; and a pole-placement linear-quadratic regulator solution based on scheduled gains is proposed. Results are presented in graphs for (1) simulations based on configuration MB3 (showing the importance of accounting for the cross-inertia terms in the TEA estimate) and (2) simulations of a stepwise change from configuration MB3 to the 'assembly complete' stage over 130 orbits (indicating that the present ACCM scheme maintains sufficient control over slowly varying Space Station dynamics).

Parlos, Alexander G.↗

A multidimensional finite element method for CFD

A finite element method is used to solve the equations of motion for 2- and 3-D fluid flow. The time-dependent equations are solved explicitly using quadrilateral (2-D) and hexahedral (3-D) elements, mass lumping, and reduced integration. A Petrov-Galerkin technique is applied to the advection terms. The method requires a minimum of computational storage, executes quickly, and is scalable for execution on computer systems ranging from PCs to supercomputers.

Pepper, Darrell W.↗

Navigation strategies for multiple autonomous mobile robots moving in formation

The problem of deriving navigation strategies for a fleet of autonomous mobile robots moving in formation is considered. Here, each robot is represented by a particle with a spherical effective spatial domain and a specified cone of visibility. The global motion of each robot in the world space is described by the equations of motion of the robot's center of mass. First, methods for formation generation are discussed. Then, simple navigation strategies for robots moving in formation are derived. A sufficient condition for the stability of a desired formation pattern for a fleet of robots each equipped with the navigation strategy based on nearest neighbor tracking is developed. The dynamic behavior of robot fleets consisting of three or more robots moving in formation in a plane is studied by means of computer simulation.

Wang, P. K. C.↗

Spatial Operator Algebra for multibody system dynamics

The Spatial Operator Algebra framework for the dynamics of general multibody systems is described. The use of a spatial operator-based methodology permits the formulation of the dynamical equations of motion of multibody systems in a concise and systematic way. The dynamical equations of progressively more complex grid multibody systems are developed in an evolutionary manner beginning with a serial chain system, followed by a tree topology system and finally, systems with arbitrary closed loops. Operator factorizations and identities are used to develop novel recursive algorithms for the forward dynamics of systems with closed loops. Extensions required to deal with flexible elements are also discussed.

Rodriguez, G.↗

A Sectored-One-Dimensional Model for Simulating Combustion Instabilities in Premix Combustors

A one-dimensional, CFD based combustor simulation has been developed that exhibits self-excited, thermoacoustic oscillations in premixed combustor geometries that typically have large, abrupt changes in cross sectional area. The combustor geometry is approximated by dividing it into a finite number of one-dimensional sectors. Within each sector, the equations of motion are integrated numerically, along with a species transport and a reaction equation. Across the sectors, mass and energy are conserved, and momentum loss is prescribed using appropriately compatible boundary conditions that account for the area change. The resulting simulation and associated boundary conditions essentially represent a one-dimensional, multi-block technique. Details of the simulation code are presented herein. Results are then shown comparing experimentally observed and simulated operation of a particular combustor rig that exhibited different instabilities at different operating points. It will be shown that the simulation closely matched the rig data in oscillation amplitudes, frequencies, and operating points at which the instabilities occurred. Finally, advantages and limitations of the simulation technique are discussed.

Paxson, Daniel E.↗