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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,135 records · Page 63

All sources NASA NTRS DOE OSTI

NASA NTRS1973

Advanced control concepts

The selection of a trim solution that provides the space shuttle with the highest level of performance and dynamic control in the presense of wind disturbances and bias torques due to misalignment of rocket engines is described. It was determined that engine gimballing is insufficient to provide control to trim the vehicle for headwind and sidewind disturbances, and that it is necessary to use aerodynamic surfaces in conjunction with engine gimballing to achieve trim. The algebraic equations for computing the trim solution were derived from the differential equations describing the motion of the vehicle by substituting the desired trim conditions. The general problem of showing how the trim equations are derived from the equations of motion and the mathematical forms of the performance criterion is discussed in detail, along with the general equations for studying the dynamic response of the trim solution.

Hutton, M. F.↗

DOE OSTI2021

Gravitational wave propagation beyond general relativity: waveform distortions and echoes

We study the cosmological propagation of gravitational waves (GWs) beyond general relativity (GR) across homogeneous and isotropic backgrounds. We consider scenarios in which GWs interact with an additional tensor field and use a parametrized phenomenological approach that generically describes their coupled equations of motion. We analyze four distinct classes of derivative and non-derivative interactions: mass, friction, velocity, and chiral. We apply the WKB formalism to account for the cosmological evolution and obtain analytical solutions to these equations. We corroborate these results by analyzing numerically the propagation of a toy GW signal. We then proceed to use the analytical results to study the modified propagation of realistic GWs from merging compact binaries, assuming that the GW signal emitted is the same as in GR. We generically find that tensor interactions lead to copies of the originally emitted GW signal, each one with its own possibly modified dispersion relation. These copies can travel coherently and interfere with each other leading to a scrambled GW signal, or propagate decoherently and lead to echoes arriving at different times at the observer that could be misidentified as independent GW events. Depending on the type of tensor interaction, the detected GW signal may exhibit amplitude and phase distortions with respect to a GW waveform in GR, as well as birefringence effects. Furthermore, we discuss observational probes of these tensor interactions with both individual GW events, as well as population studies for both ground- and space-based detectors.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1972

Numerical calculations of flow fields

The solutions to the equations of motion for inviscid fluid flow around a pointed elliptic cone at incidence are presented. The numerical method used, MacCormack's second order preferential predictor-corrector finite difference approximation, is applied to the fluid flow equations derived in conservation-law form. The entropy boundary condition, hitherto unused for elliptic cone problems, is investigated and compared to reflection boundary condition solutions. The stagnation streamline movement of the inclined elliptic cone is noted and surface pressure coefficients are plotted. Also presented are solutions for an elliptic cone and a circular cone at zero incidence and a circular cone at a small angle of attack. Comparisons are made between these present solutions and previously published theory.

Anderson, D. M.↗

NASA NTRS1974

Estimation of elastic aircraft parameters using the maximum likelihood method

The application of the maximum likelihood method to estimate the aerodynamic parameters of elastic flight vehicles in a symmetric flight condition is discussed. In this application, particular attention is directed toward the center of mass, elastic deformation, and sensor equations of motion. It is shown that the two major computational problems to be overcome are the inversion of large-sized matrices and the time-wise integration of a large number of linear, ordinary, differential equations.

Schwanz, R. C.↗

NASA NTRS1974

Numerical solution of the steady-state Navier-Stokes equations for hypersonic flow about blunt axisymmetric bodies

The steady-state Navier-Stokes equations are solved for hypersonic flow about blunt axisymmetric bodies. The equations of motion are solved by successive approximations using an implicit finite-difference scheme. The results are compared with viscous shock-layer theory, experimental data, and time-dependent solutions of the Navier-Stokes equations. It is demonstrated that viscous shock-layer theory is sufficiently accurate for the range of flight conditions normally encountered by entry vehicles.

Anderson, E. C.↗

NASA NTRS1977

An analytical method for ride quality of flexible airplanes

A new and easily used method of aircraft ride quality analysis is developed, and its use is illustrated by a numerical example using DC-8 airplane equations of motion. The method readily allows any type of stability augmentation system to be included for good handling qualities. A particularly powerful method of specifying and achieving all roots of a desired closed-loop characteristic equation by use of full state variable feedback is described and applied to the DC-8 example.

Swaim, R. L.↗

NASA NTRS1976

Numerical solution of the steady-state Navier-Stokes equations for hypersonic flow about blunt axisymmetric bodies

Anderson, E. C.↗

NASA NTRS1994

Computationally efficient multibody simulations

Computationally efficient approaches to the solution of the dynamics of multibody systems are presented in this work. The computational efficiency is derived from both the algorithmic and implementational standpoint. Order(n) approaches provide a new formulation of the equations of motion eliminating the assembly and numerical inversion of a system mass matrix as required by conventional algorithms. Computational efficiency is also gained in the implementation phase by the symbolic processing and parallel implementation of these equations. Comparison of this algorithm with existing multibody simulation programs illustrates the increased computational efficiency.

Ramakrishnan, Jayant↗

NASA NTRS2003

A Flight Dynamics Model for a Small Glider in Ambient Winds

In this paper we describe the equations of motion developed for a point-mass zero-thrust (gliding) aircraft model operating in an environment of spatially varying atmospheric winds. The wind effects are included as an integral part of the flight dynamics equations, and the model is controlled through the three aerodynamic control angles. Formulas for the aerodynamic coefficients for this model are constructed to include the effects of several different aspects contributing to the aerodynamic performance of the vehicle. Characteristic parameter values of the model are compared with those found in a different set of small glider simulations. We execute a set of example problems which solve the glider dynamics equations to find aircraft trajectory given specified control inputs. The ambient wind conditions and glider characteristics are varied to compare the simulation results under these different circumstances.

Beeler, Scott C.↗

NASA NTRS2003

A Flight Dynamics Model for a Small Glider in Ambient Winds

In this paper we describe the equations of motion developed for a point-mass zero-thrust (gliding) aircraft model operating in an environment of spatially varying atmospheric winds. The wind effects are included as an integral part of the flight dynamics equations, and the model is controlled through the three aerodynamic control angles. Formulas for the aerodynamic coefficients for this model are constructed to include the effects of several different aspects contributing to the aerodynamic performance of the vehicle. Characteristic parameter values of the model are compared with those found in a different set of small glider simulations. We execute a set of example problems which solve the glider dynamics equations to find the aircraft trajectory given specified control inputs. The ambient wind conditions and glider characteristics are varied to compare the simulation results under these different circumstances.

Beeler, Scott C.↗

NASA NTRS2015

Multi-Objective Hybrid Optimal Control for Multiple-Flyby Interplanetary Mission Design using Chemical Propulsion

Preliminary design of high‐thrust interplanetary missions is a highly complex process. The mission designer must choose discrete parameters such as the number of flybys and the bodies at which those flybys are performed. For some missions, such as surveys of small bodies, the mission designer also contributes to target selection. In addition, real‐valued decision variables, such as launch epoch, flight times, maneuver and flyby epochs, and flyby altitudes must be chosen. There are often many thousands of possible trajectories to be evaluated. The customer who commissions a trajectory design is not usually interested in a point solution, but rather the exploration of the trade space of trajectories between several different objective functions. This can be a very expensive process in terms of the number of human analyst hours required. An automated approach is therefore very desirable. This work presents such an approach by posing the impulsive mission design problem as a multi‐objective hybrid optimal control problem. The method is demonstrated on several real‐world problems. Two assumptions are frequently made to simplify the modeling of an interplanetary high‐thrust trajectory during the preliminary design phase. The first assumption is that because the available thrust is high, any maneuvers performed by the spacecraft can be modeled as discrete changes in velocity. This assumption removes the need to integrate the equations of motion governing the motion of a spacecraft under thrust and allows the change in velocity to be modeled as an impulse and the expenditure of propellant to be modeled using the time‐independent solution to Tsiolkovsky's rocket equation [1]. The second assumption is that the spacecraft moves primarily under the influence of the central body, i.e. the sun, and all other perturbing forces may be neglected in preliminary design. The path of the spacecraft may then be modeled as a series of conic sections. When a spacecraft performs a close approach to a planet, the central body switches from the sun to that planet and the trajectory is modeled as a hyperbola with respect to the planet. This is known as the method of patched conics. The impulsive and patched‐conic assumptions significantly simplify the preliminary design problem.

Optimization↗

NASA NTRS2019

Modal Analysis of the Orion Capsule Two Parachute System

As discussed in Ref [1], it is apparent from flight tests that the system made up of two main parachutes and a capsule can undergo several distinct dynamical behaviors. The most significant and problematic of these is the pendulum mode in which the system develops a pronounced swinging motion with an amplitude of up to 24 deg. Large excursions away from vertical by the capsule could cause it to strike the ground at a large horizontal or vertical speed and jeopardize the safety of the astronauts during a crewed mission. In reference [1], Ali et al. summarized a series of efforts taken by the Capsule Parachute Assembly System (CPAS) Program to understand and mitigate the pendulum issue. The period of oscillation and location of the system's pivot point are determined from post-flight analysis. Other noticeable but benign modes include: 1) flyout (scissors) mode, where the parachutes move back and forth symmetrically with respect to the vertical axis similar to the motion of a pair of scissors; 2) maypole mode, where the two parachutes circle around the vertical axis at a nearly constant radius and period; and 3) breathing mode, in which deformation of the non-rigid canopies affects the axial acceleration of the system in an oscillatory manner. Because these modes are relatively harm- less, little effort has been devoted to analyzing them in comparison with the pendulum motion. Motions of the actual system made up of two parachutes and a capsule are extremely complicated due to nonlinearities and flexibility effects. Often it is difficult to obtain insight into the fundamental dynamics of the system by examining results from a multi-body simulation based on nonlinear equations of motion (EOMs). As a part of this study, the dynamics of each mode observed during flight is derived from first principles on an individual basis by making numerous simplifications along the way. The intent is to gain a better understanding into the behavior of the complex multi-body system by studying the reduced set of differential equations associated with each mode. This approach is analogous to the traditional modal analysis technique used to study airplane flight dynamics, in which the full nonlinear behavior of the airframe is decomposed into the phugoid and short period modes for the longitudinal dynamics and the spiral, roll-subsidence, and dutch-roll modes for the lateral dynamics. It is important to note that the study does not address the mechanisms that cause the system to transition from one mode to another, nor does it discuss motions during which two or more modes occur simultaneously.

Pei, J.↗

NASA NTRS1982

Astrometric observations and the motion of comet Halley

Beginning with studies related to the 1835 appearance of comet Halley, it became apparent that comets may undergo substantial nongravitational perturbations due to the rocket-like effect of an outgassing cometary nucleus. Marsden (1968, 1969) modeled these nongravitational effects with a semiempirical term in the cometary equations of motion. The most successful orbital solution was consistent with an outgassing water ice nucleus in direct rotation. For the conduction and coordination of accurate astrometric observations of comet Halley, an Astrometry Network is being established within the International Halley Watch program. Attention is given to the support of the Halley flight projects provided by the Astrometry Network, and the support of ground based observers by this network.

Yeomans, D. K.↗

DOE OSTI2022

Simulation of crumpled sheets via alternating quasistatic and dynamic representations

In this work, we present a method for simulating the large-scale deformation and crumpling of thin, elastoplastic sheets. Motivated by the physical behavior of thin sheets during crumpling, two different formulations of the governing equations of motion are used: (1) a quasistatic formulation that effectively describes smooth deformations, and (2) a fully dynamic formulation that captures large changes in the sheet's velocity. The former is a differential-algebraic system of equations integrated implicitly in time, while the latter is a set of ordinary differential equations (ODEs) integrated explicitly. We adopt a hybrid integration scheme to adaptively alternate between the quasistatic and dynamic representations as appropriate. Further, we demonstrate the capacity of this method to effectively simulate a variety of crumpling phenomena. Finally, we show that statistical properties, notably the accumulation of creases under repeated loading, as well as the area distribution of facets, are consistent with experimental observations.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

Real-time coupled-cluster approach for the cumulant Green's function

Green’s function methods within many-body perturbation theory provide a general framework for treating electronic correlations in excited states. Here we investigate the cumulant form of the one-electron Green’s function based on the coupled-cluster equation of motion approach in an extension of our previous study. The approach yields a non-perturbative expression for the cumulant in terms of the solution to a set of coupled first order, non-linear differential equations. The method thereby adds non-linear corrections to traditional cumulant methods linear in the self energy. The approach is applied to the core-hole Green’s function and illustrated for a number of small molecular systems. For these systems we find that the non-linear contributions lead to significant improvements both for quasiparticle properties such as core-level binding energies, as well as the satellites corresponding to inelastic losses observed in photoemission spectra.

Vila, Fernando D.↗

DOE OSTI2022

Geometric integration of classical spin dynamics via a mean-field Schrödinger equation

The Landau-Lifshitz equation describes the time evolution of magnetic dipoles and can be derived by taking the classical limit of a quantum mechanical spin Hamiltonian. To take this limit, one constrains the many-body quantum state to a tensor product of coherent states, thereby neglecting entanglement between sites. Expectation values of the quantum spin operators produce the usual classical spin dipoles. One may also consider expectation values of polynomials of the spin operators, leading to quadrupole and higher-order spin moments, which satisfy a dynamical equation of motion that generalizes the Landau-Lifshitz dynamics [Zhang and Batista, Phys. Rev. B 104, 104409 (2021)]. Here we reformulate the dynamics of these N 2 –1 generalized spin components as a mean-field Schrödinger equation on the N-dimensional coherent state. Furthermore, this viewpoint suggests efficient integration methods that respect the local symplectic structure of the classical spin dynamics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NASA NTRS1977

Transient wave propagation in composite plates due to impact

Stress-wave propagation in a composite plate due to an impact has been examined using anisotropic elasticity theory. The plate is modelled as a number of identical anisotropic layers whose elastic moduli are given. Mindlin's (1959) approximate theory of plates is then applied to each layer to obtain a set of difference-differential equations of motion by using the interlaminar stresses and displacements as explicit variables. Dispersion relationships for harmonic waves are found when traction-free boundary conditions are applied to both surfaces of the plate. The difference-differential equations are reduced to difference equations via integral transforms. With given impact boundary conditions these equations are solved for an arbitrary number of layers in the plate, and the transient propagation of a stress wave is calculated by means of a fast Fourier transform algorithm.

Sung Kim, B.↗

NASA NTRS1986

Sharply curved turn around duct flow predictions using spectral partitioning of the turbulent kinetic energy and a pressure modified wall law

Computational predictions of turbulent flow in sharply curved 180 degree turn around ducts are presented. The CNS2D computer code is used to solve the equations of motion for two-dimensional incompressible flows transformed to a nonorthogonal body-fitted coordinate system. This procedure incorporates the pressure velocity correction algorithm SIMPLE-C to iteratively solve a discretized form of the transformed equations. A multiple scale turbulence model based on simplified spectral partitioning is employed to obtain closure. Flow field predictions utilizing the multiple scale model are compared to features predicted by the traditional single scale k-epsilon model. Tuning parameter sensitivities of the multiple scale model applied to turn around duct flows are also determined. In addition, a wall function approach based on a wall law suitable for incompressible turbulent boundary layers under strong adverse pressure gradients is tested. Turn around duct flow characteristics utilizing this modified wall law are presented and compared to results based on a standard wall treatment.

Santi, L. Michael↗