SEARCH · Engineering Papers

Results for “Moment equations”

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 721 records · Page 40

All sources NASA NTRS DOE OSTI

NASA NTRS1994

A note on the accuracy of spectral method applied to nonlinear conservation laws

Fourier spectral method can achieve exponential accuracy both on the approximation level and for solving partial differential equations if the solutions are analytic. For a linear partial differential equation with a discontinuous solution, Fourier spectral method produces poor point-wise accuracy without post-processing, but still maintains exponential accuracy for all moments against analytic functions. In this note we assess the accuracy of Fourier spectral method applied to nonlinear conservation laws through a numerical case study. We find that the moments with respect to analytic functions are no longer very accurate. However the numerical solution does contain accurate information which can be extracted by a post-processing based on Gegenbauer polynomials.

Shu, Chi-Wang↗

NASA NTRS1994

Sloshing dynamics modulated angular momentum fluctuations actuated by gravity gradient acceleration associated with spacecraft slew motion

The mathematical formulation of orbital spacecraft sloshing dynamics for a partially filled cryogenic superfluid liquid helium II in a dewar container actuated by the gravity gradient acceleration associated with slew motion is studied. The Advanced X-Ray Astrophysics Facility-Spectroscopy (AXAF-S) spacecraft is chosen as a practical example in this study. Explicit mathematical expressions which manage orbital galaxy gradient acceleration with slew motion which is acting on the spacecraft fluid systems are derived. The numerical computation of sloshing dynamics is based on the non-inertia frame spacecraft bound coordinates and the solution of time-dependent, three-dimensional formulations of partial differential equations subject to initial and boundary conditions. This study discloses the capillary effect of sloshing dynamics governed liquid-vapor interface fluctuations, angular momentum and moment fluctuations of fluid system, and also buble mass center fluctuations driven by the gravity gradient acceleration associated with slew motion which affects the stability of the orbital spacecraft fluid system in a microgravity environment.

Hung, R. J.↗

NASA NTRS2017

Two-Pendulum Model of Propellant Slosh in Europa Clipper PMD Tank

The objective of this fluids analysis is to model propellant slosh for the Europa Clipper mission using a two-pendulum model, such that controls engineers can predict slosh behavior during the mission. Propellant slosh causes shifts in center of mass and exerts forces and torques on the spacecraft which, if not adequately controlled, can lead to mission failure. The two-pendulum model provides a computationally simple model that can be used to predict slosh for the Europa Clipper tank geometry. The Europa Clipper tank is cylindrical with a domed top and bottom and includes a propellant management device (PMD). Due to the lack of experimental data in low gravity environments, computational fluid dynamics (CFD) simulation results were used as 'real' slosh behavior for two propellants at three fill fractions. Key pendulum parameters were derived that allow the pendulum model's center of mass, forces, and moments to closely match the CFD data. The parameter trends were examined as a function of tank fill fraction and compared with solutions to analytic equations that describe the frequency of slosh in tanks with simple geometries. The trends were monotonic as expected, and parameters resembled analytical predictions; any differences could be explained by the specific differences in the geometry of the tank. This paper summarizes the new method developed at Goddard Space Flight Center (GSFC) for deriving pendulum parameters for two-pendulum equivalent sloshing models. It presents the results of this method and discusses the validity of the results. This analysis is at a completed stage and will be applied in the immediate future to the evolving tank geometry as Europa Clipper moves past its preliminary design review (PDR) phase.

Propellant Slosh↗

NASA NTRSDate not supplied

Accuracy of Center of Pressure Determination via Motion Capture

BACKGROUND: This study was conducted to support the stability assessment for tasks in lunar gravity and exercises on a Vibration Isolation and Stabilization (VIS) system in microgravity based on the dynamic feasibility criterion of whether the calculated position of the center of pressure (COP) falls within the base of support (BOS) which outlines the subject’s feet. Motion capture data combined with biomechanical modeling and simulation allows the forces and moments between the human and the VIS platform to be computed and the position of the COP as well as the location and shape of the BOS to be determined. The goal of this study was to assess the accuracy of the COP trajectory calculated using motion capture-based data. METHODS AND RESULTS: To obtain the dynamic quantities from which COP is calculated, motion capture data is first collected in the 1g lab environment by recording the trajectories of passive retroreflective markers placed on a subject during exercise or performance of a given task. The OpenSim [1] inverse kinematics (IK) tool is used to fit a scaled subject model to recorded marker trajectories while minimizing marker error to obtain joint angles. Then, a custom OpenSim plugin [2] is used to determine the subject’s time-varying moment of inertia and its time derivative, center of mass (CM) position, velocity, and acceleration, as well as the angular momentum and its time derivative relative to the subject’s CM. Some of these quantities are not needed for modeling tasks performed on a stationary lunar surface but, due to the moving exercise platform, are needed to model VIS response to the subject’s motion. Hand positions, used in calculating a cable force if present, are recorded as well. These quantities are used to calculate the total force (F ⃗^((plate) )) and moment (M ⃗^((plate) )) exerted by the lunar surface or the VIS plate on the subject’s shoe soles. COP is then calculated from the following equations: r_x^((cop) )= M_z^((plate) )/F_y^((plate) ) and r_z^((cop) )= 〖-M〗_x^((plate) )/F_y^((plate) ), where the y axis is normal to the surface. COP accuracy for feasibility assessments is then determined by whether it falls within the BOS, which is also computed by the plugin. To study the accuracy of COP calculated from motion capture, we first investigated whether COP remained within the BOS, as it must, for exercises performed in the 1g lab environment. Standard exercises such as back squat and deadlift were analyzed, as well as more explosive exercises including hang clean and press. Cases in which the COP exited the BOS indicated that COP accuracy required further investigation. In this study, an exercise device with cables was used, so cable force modeling accuracy should also be considered. In a separate study, we collected motion capture and force plate data for twenty-seven motions not involving an exercise device. About a third were genuine countermeasures exercises (e.g., hang clean and press), some were relevant for lunar tasks (e.g., object pick up), and the rest were of a “unit test” nature (e.g., swaying back and forth or side to side). Motion capture-based COP positions were compared with force plate measured COP. We found that while force plate measured COP remained within the BOS, motion capture-based COP was observed to briefly exit the BOS on occasion. Techniques to mitigate IK artifacts and filtering of calculated data could be used to improve the agreement of calculated and measured results, resolving excursions from the BOS within this dataset. The mean error between calculated and measured COP was found to be less than 6 mm. Additionally, we derived and investigated equations for the COP in terms of the cable force, cable location, as well as the human CM position, acceleration, and angular momentum with respect to the CM, and analyzed them for sensitivity to errors in individual quantities. Several were found, but the most significant one was that when the vertical force on the feet approaches zero, indicating a near-detachment or ‘jump off’ condition, errors are amplified. This is consistent with the observation that in the absence of pressure, the concept of the center of pressure would become meaningless.

C A Bell↗

NASA NTRS2012

Low-degree Structure in Mercury's Planetary Magnetic Field

The structure of Mercury's internal magnetic field has been determined from analysis of orbital Magnetometer measurements by the MESSENGER spacecraft. We identified the magnetic equator on 531 low-altitude and 120 high-altitude equator crossings from the zero in the radial cylindrical magnetic field component, Beta (sub rho). The low-altitude crossings are offset 479 +/- 6 km northward, indicating an offset of the planetary dipole. The tilt of the magnetic pole relative to the planetary spin axis is less than 0.8 deg.. The high-altitude crossings yield a northward offset of the magnetic equator of 486 +/- 74 km. A field with only nonzero dipole and octupole coefficients also matches the low-altitude observations but cannot yield off-equatorial Beta (sub rho) = 0 at radial distances greater than 3520 km. We compared offset dipole and other descriptions of the field with vector field observations below 600 km for 13 longitudinally distributed, magnetically quiet orbits. An offset dipole with southward directed moment of 190 nT-R-cube (sub M) yields root-mean-square (RMS) residuals below 14 nT, whereas a field with only dipole and octupole terms tuned to match the polar field and the low-altitude magnetic equator crossings yields RMS residuals up to 68 nT. Attributing the residuals from the offset-dipole field to axial degree 3 and 4 contributions we estimate that the Gauss coefficient magnitudes for the additional terms are less than 4% and 7%, respectively, relative to the dipole. The axial alignment and prominent quadrupole are consistent with a non-convecting layer above a deep dynamo in Mercury's fluid outer core.

magnetometer↗

NASA NTRS2017

Evaluation of Full Reynolds Stress Turbulence Models in FUN3D

Full seven-equation Reynolds stress turbulence models are a relatively new and promising tool for todays aerospace technology challenges. This paper uses two stress-omega full Reynolds stress models to evaluate challenging flows including shock-wave boundary layer interactions, separation and mixing layers. The Wilcox and the SSGLRR full second-moment Reynolds stress models are evaluated for four problems: a transonic two-dimensional diffuser, a supersonic axisymmetric compression corner, a compressible planar shear layer, and a subsonic axisymmetric jet. Simulation results are compared with experimental data and results using the more commonly used Spalart-Allmaras (SA) one-equation and the Menter Shear Stress Transport (SST) two-equation models.

computational fluid dynamics↗

DOE OSTI2023

On the Cauchy problem for Boltzmann equation modeling a polyatomic gas

In the present article, we consider the Boltzmann equation that models a polyatomic gas by introducing one additional continuous variable, referred to as microscopic internal energy. We establish existence and uniqueness theory in the space homogeneous setting for the full non-linear case, under an extended Grad-type assumption on transition probability rates, which comprises hard potentials for both the relative speed and internal energy with the rate in the interval 0,2, multiplied by an integrable angular part and integrable partition functions. The Cauchy problem is resolved by means of an abstract ordinary differential equation (ODE) theory in Banach spaces for the initial data with finite and strictly positive gas mass and energy, finite momentum, and additionally finite K* polynomial moment, with K* depending on the rate of the transition probability and the structure of a polyatomic molecule or its internal degrees of freedom. Moreover, we prove that polynomially and exponentially weighted Banach space norms associated with the solution are both generated and propagated uniformly in time.

Physics↗

NASA NTRS1973

V/STOL tilt rotor aircraft study mathematical model for a real time simulation of a tilt rotor aircraft (Boeing Vertol Model 222), volume 8

This report documents the development of a real time mathematical model of a tilt rotor aircraft. This mathematical model is to be used in conjunction with the NASA Flight Simulator for Advanced Aircraft (FSAA) at Ames Research Center for evaluation of aircraft performance and handling qualities. In addition to developing the mathematical model, a parallel programming effort was conducted utilizing Boeing-Vertol's Hybrid Simulation Laboratory for the purpose of developing and evaluating model simplification. The mathematical model is an eleven degree of freedom total force model. This model includes the basic six degree of freedom rigid body .outer loop equations written about the instantaneous center of gravity with the inertial and aerodynamic terms included. The rotor is treated as a point source of forces and moments with appropriate response time lags and actuator dynamics. The wing has one vertical bending and one wing torsion degree of freedom. These structural degrees of freedom are treated on a "quasistatic" basis; i.e., the natural frequencies of vibration of the structure are much higher than the· frequencies of the rigid body motion, and the coupling is in the aerodynamic terms. Each nacelle has an independent pitch degree of freedom about the wing pivot. The aerodynamics of the wing, tail, rotors, landing gear and fuselage are included. Wing and tail mutual interference effects and turbine engine performance and dynamic responses' are represented.

H Rosenstein↗

NASA NTRS1986

Kalman filtering, smoothing and recursive robot arm forward and inverse dynamics

The inverse and forward dynamics problems for multi-link serial manipulators are solved by using recursive techniques from linear filtering and smoothing theory. The pivotal step is to cast the system dynamics and kinematics as a two-point boundary-value problem. Solution of this problem leads to filtering and smoothing techniques identical to the equations of Kalman filtering and Bryson-Frazier fixed time-interval smoothing. The solutions prescribe an inward filtering recursion to compute a sequence of constraint moments and forces followed by an outward recursion to determine a corresponding sequence of angular and linear accelerations. In addition to providing techniques to compute joint accelerations from applied joint moments (and vice versa), the report provides an approach to evaluate recursively the composite multi-link system inertia matrix and its inverse. The report lays the foundation for the potential use of filtering and smoothing techniques in robot inverse and forward dynamics and in robot control design.

Rodriguez, G.↗

NASA NTRS1995

The estimation of material and patch parameters in a PDE-based circular plate model

The estimation of material and patch parameters for a system involving a circular plate, to which piezoceramic patches are bonded, is considered. A partial differential equation (PDE) model for the thin circular plate is used with the passive and active contributions form the patches included in the internal and external bending moments. This model contains piecewise constant parameters describing the density, flexural rigidity, Poisson ratio, and Kelvin-Voigt damping for the system as well as patch constants and a coefficient for viscous air damping. Examples demonstrating the estimation of these parameters with experimental acceleration data and a variety of inputs to the experimental plate are presented. By using a physically-derived PDE model to describe the system, parameter sets consistent across experiments are obtained, even when phenomena such as damping due to electric circuits affect the system dynamics.

Banks, H. T.↗

NASA NTRS2006

Free Body Dynamics of a Spinning Cylinder With Planar Restraint (a.k.a. Barrel of Fun)

The dynamic motion of a cylinder on a floor or hard surface is both entertaining and instructive. With maintenance torques, motion can be sustained and controlled as illustrated in a video clip that can be viewed in the PDF file of this document. The analysis of such a cylinder with and without end caps is burned on rotation about its center of mass and restrained by a plane normal to the axis passing through its center of mass at an angle alpha. For small values of alpha, the governing equations are simplified, and for symmetric bodies, stability requires rotation greater than [2 square root of(JWL*)]/J(sub X), where J is the transverse mass moment of inertia, W is the weight of the cylinder, L* is the cylinder length from the base to the center of mass, and JX is the mass moment of inertia about the longitudinal axis OX of the barrel. Comparisons to data are made and some applications are discussed.

Moraru, Laurentiu↗

NASA NTRSDate not supplied

Transonic Correction Method for Flight Dynamic Stability Analysis of Mach 0.745 Transonic Truss-Braced Wing

This paper presents a transonic correction method for obtaining dynamic stability derivatives for flight dynamic stability analysis. The method provides a transonic correction to the Theodorsen's theory of unsteady aerodynamics using FUN3D CFD solver of unsteady Reynolds-averaged Navier-Stokes equations (RANS) for a series of wing sections of the Mach 0.745 Transonic Truss-Braced Wing in pitch and plunge oscillations. Unsteady lift and pitching moment coefficients are obtained and used to develop the correction terms in the Theodorsen's theory to account for transonic aerodynamics. The unsteady lift and pitching moment derivatives with respect to the unsteady angle of attack are obtained as frequency response functions of the reduced frequency. These frequency response functions are used to compute the dynamic stability derivatives of lift and pitching moment due to the angle of attack and pitch rate and the dynamic stability derivatives for the rolling moment and yawing moment with respect to the roll rate and yaw rate. A transonic correction is applied to steady-state stability derivatives computed by VSPAERO solver using transonic small disturbance code TSFOIL coupled to an integral boundary method. A dynamic stability analysis is conducted for longitudinal and lateral-directional motions. Without transonic corrections and dynamic stability derivatives, the analysis indicates an unstable phugoid mode. The transonic correction applied to the steady-state stability derivatives computed by VSPAERO shows a stable phugoid mode. This is due to the increase of the drag stability derivatives as a result of the additional wave drag contribution in transonic flow. The effect of the transonic dynamic stability derivatives is observed to be a significant contributor to the increase in the damping values of all the flight dynamic modes of the Mach 0.745 Transonic-Truss Braced Wing.

Transonic↗

NASA NTRS1977

A mathematical force and moment model of a UH-1H helicopter for flight dynamics simulations

A model of a UH-1H helicopter was developed to support flight simulations and for developmental work on an avionics system known as V/STOLAND system. Equations and numerical values of constants used to represent the helicopter are presented. Responses to stop inputs of the cyclic and collective controls are shown and compared with flight test data for a UH-1H. The model coefficients were adjusted in an attempt to get a consistant match with the flight time histories at hover and 60 knots. Response matching was obtained at 60 knots, but the matching at hover was not as successful. Pilot evaluations of the model, both fixed and moving base, were made.

Talbot, P. D.↗

NASA NTRS1994

Superfluid helium sloshing dynamics induced oscillations and fluctuations of angular momentum, force and moment actuated on spacecraft driven by gravity gradient or jitter acceleration associated with slew motion

The generalized mathematical formulation of sloshing dynamics for partially filled liquid of cryogenic superfluid helium II in dewar containers driven by the gravity gradient and jitter accelerations associated with slew motion for the purpose to perform scientific observation during the normal spacecraft operation are investigated. An example is given with the Advanced X-Ray Astrophysics Facility-Spectroscopy (AXAF-S) for slew motion which is responsible for the sloshing dynamics. The jitter accelerations include slew motion, spinning motion, atmospheric drag on the spacecraft, spacecraft attitude motions arising from machinery vibrations, thruster firing, pointing control of spacecraft, crew motion, etc. Explicit mathematical expressions to cover these forces acting on the spacecraft fluid systems are derived. The numerical computation of sloshing dynamics is based on the non-inertia frame spacecraft bound coordinate, and solve time-dependent, three-dimensional formulations of partial differential equations subject to initial and boundary conditions. The explicit mathematical expressions of boundary conditions to cover capillary force effect on the liquid-vapor interface in microgravity environments are also derived. The formulations of fluid moment and angular moment fluctuations in fluid profiles induced by the sloshing dynamics, together with fluid stress and moment fluctuations exerted on the spacecraft dewar containers have also been derived. Examples are also given for cases applicable to the AXAF-S spacecraft sloshing dynamics associated with slew motion.

Hung, R. J.↗

NASA NTRS2005

Stochastic Representation of Chaos using Terminal Attractors

A nonlinear version of the Liouville equation based upon terminal attractors is proposed for describing post-instability motions of dynamical systems with exponential divergence of trajectories such as those leading to chaos and turbulence. As a result, the post-instability motions are represented by expectations, variances, and higher moments of the state variables as functions of time. The proposed approach can be applied to conservative chaos, and in particular, to n-bodies problem, as well as to dissipative systems, and in particular, to chaotic attractors and turbulence.

stochastic processes↗

DOE OSTI2024

Leading order track functions in a hot and dense QGP

We study the modifications to the fragmentation pattern of partons into charged particles in the presence of a hot and dense quark gluon plasma. To this end, we analyze the perturbative renormalization group equations of the track functions, which describe the energy fraction carried by charged hadrons. Focusing on pure Yang-Mills theory, we compute the lowest-order moments of the medium-modified track functions, which are found to be sensitive to the reduced phase space for emissions in the medium and to energy loss. We use the extracted moments to calculate the energy energy correlator (EEC) on tracks in the collinear limit. The EEC on medium-evolved tracks does not differ qualitatively from the EEC on vacuum tracks despite being sensitive to the color decoherence transition and suppressing the distribution due to quenching, as seen in other jet observables. Published by the American Physical Society 2024

Barata, João (ORCID:0000000342864555)↗

NASA NTRS1944

A theoretical investigation of the rolling oscillations of an airplane with ailerons free

Cohen, Doris↗

NASA NTRS1944

A Theoretical Investigation of the Rolling Oscillations of an Airplane with Ailerons Free

An analysis is made of the stability of an airplane with ailerons free, with particular attention to the motions when the ailerons have a tendency to float against the wind. The present analysis supersedes the aileron investigation contained in NACA Report No. 709. The equations of motion are first written to include yawing and sideslipping, and it is demonstrated that the principal effects of freeing the ailerons can be determined without regard to these motions. If the ailerons tend to float against the wind and have a high degree of aerodynamic balance, rolling oscillations, in addition to the normal lateral oscillations, are likely to occur. On the basis of the equations including only the rolling motion and the aileron deflection, formulas are derived for the stability and damping of the rolling oscillations in terms of the hinge moment derivatives and other characteristics of the ailerons and airplane. Charts are also presented showing the oscillatory regions and stability boundaries for a fictitious airplane of conventional proportions. The effects of friction in the control system are investigated and discussed. If the ailerons tend to trail with the wind, the condition for stable variation of stick force with aileron deflection is found to determine the amount of aerodynamic balance that may be used. If the ailerons tend to float against the wind, the period and damping of the rolling oscillations are found to be satisfactory (in a mass-balanced system) so long as the restoring moment is not completely balanced out. Unbalanced mass behind the hinge, however, has an unfavorable effect on the damping of the oscillations and so shifts the boundary that close aerodynamic balance may not be attainable.

Cohen, Doris↗