Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “motion correction”

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 541 records · Page 30

TOWARD A METHOD FOR SCALING HUMAN BODY MODELS IN AN IMU-BASED WORKFLOW

BACKGROUND Scaled biomechanical models can more accurately inform crew health decisions when tailored to the wide range of astronaut sizes. One component to improve scaling of existing models to better represent each unique astronaut’s size is the individual length scaling of limbs. Traditionally, these lengths are determined by motion capture or manual measurement. A new method is herein proposed for length scaling which can be done by measuring linear and angular accelerations at a desired point during isolated motion around a point of rotation, then calculating the distance between the desired point and point of rotation. When an Inertial Measurement Unit (IMU) device is placed at the distal point of a limb, the isolated motion is about that limb’s proximal joint. For example, to measure forearm length, an IMU is placed at the wrist, the point of rotation is at the elbow, and the isolated motion is forearm flexion and extension. These calculated lengths are then used to scale models to each unique astronaut’s size, thereby improving the applicability of the model. This method of scaling limb segments can be used for any limb that has an easily defined proximal joint for the limb to rotate around including hands, arms, legs, feet. Utilizing IMUs for data collection also provides the synergistic ability to record data without a dedicated space in a room with many cameras, therefore reducing the data collection footprint, or record data where optical motion capture is not possible, such as inside a spacesuit. METHODS AND RESULTS To test this method, upper body data collection was performed with 5 Xsens DOT IMUs on a single subject. IMUs consist of an accelerometer, a gyroscope, and a magnetometer which collect linear acceleration, angular velocity, and magnetic fluctuations, respectively. Before any ground-based laboratory collection, the magnetic fluctuations are used to correct the heading of the IMU in space relative to the Earth’s magnetic field. The direct measurement of angular velocity is integrated to calculate angular acceleration. Then the linear acceleration ( a ) and angular acceleration (α) are solved using r = at/α to calculate the radius, which in this case is the distance between the IMU and the point of rotation (i.e., segment length). The IMU must be placed at the most distal point of the limb being measured (i.e., ankle if measuring lower leg length) and the test plan must consist of an isolated motion about that limb’s proximal joint (i.e., knee flexion and extension if measuring lower leg length). The distances (radii) calculated at every time interval were filtered (bandpass filter keeping 5th-90th percentile data) to eliminate outliers and spurious data that occur when the isolated motion was stopped or nearly stopped. The remaining distances were averaged, resulting in the calculated limb length. Scaling factors were then computed by dividing the calculated limb length by the unscaled model’s length. These scale factors are plugged into the Scale Tool in OpenSim [1,2] to apply the scaling to the OpenSim Full Body Rajagopal Model [3,4]. Manual measurements of limb lengths were taken before data collection started and used for comparing against the calculated lengths. The Anthropometric Survey of US Army Personnel (ANSUR II) [5] was also used as a third source of reference for limb length measurements. The following measurements were retrieved from the subject before data collection: 34.5 cm from L1 to C7 (thorax), 25.7 cm from C7 Joint Center (JC) to head vertex (neck and head), 36.3 cm from shoulder JC to elbow JC (humerus), 29.5 cm from elbow JC to wrist JC (forearm), and 16.2 cm from clavicle to acromion (clavicle). Of those five, forearm and humerus lengths were calculated using this proposed method to obtain preliminary results. The forearm length after filtering and averaging was calculated to be 37.6 cm. This is a 28% overestimation from the measured forearm length (29.5 cm). The humerus length after filtering and averaging was calculated to be 47.8 cm. This is a 31% difference from the measured subject length (36.3 cm). Sources of error include imperfect isolated motion (method currently expects that motion should be perfectly circular in a 2D plane, include no rotation of the IMU, and be relatively smooth; a more secure IMU attachment method will help), unrefined filter techniques (removed highest and lowest values with 20% high and low pass filters and no smoothing filters), arbitrary removal of stopped or near stopped data (kept data for only a short range before and after the angular velocity peaking), and a more representative method for removing gravitational acceleration is needed (current method is to zero all accelerations against a baseline taken just before the isolated motion which does not account for the gravitational acceleration changed due to IMU rotation during movement). Addressing these error sources will improve the accuracy of the limb length calculation. Next steps include creating a method for whole-body scaling estimation using individual limb scale factors. Continued pursuit of these techniques is expected to enable acquiring anthropometric information using only IMUs in real-time.

E. K. Marecki↗

Comparison of analysis and flight test data for a drone aircraft with active flutter suppression

A drone aircraft equipped with an active flutter suppression system is considered with emphasis on the comparison of modal dampings and frequencies as a function of Mach number. Results are presented for both symmetric and antisymmetric motion with flutter suppression off. Only symmetric results are given for flutter suppression on. Frequency response functions of the vehicle are presented from both flight test data and analysis. The analysis correlation is improved by using an empirical aerodynamic correction factor which is proportional to the ratio of experimental to analytical steady-state lift curve slope. The mathematical models are included and existing analytical techniques are described as well as an alternative analytical technique for obtaining closed-loop results.

Newsom, J. R.↗

Comparison of analysis and flight test data for a drone aircraft with active flutter suppression

This paper presents a comparison of analysis and flight test data for a drone aircraft equipped with an active flutter suppression system. Emphasis is placed on the comparison of modal dampings and frequencies as a function of Mach number. Results are presented for both symmetric and antisymmetric motion with flutter suppression off. Only symmetric results are presented for flutter suppression on. Frequency response functions of the vehicle are presented from both flight test data and analysis. The analysis correlation is improved by using an empirical aerodynamic correction factor which is proportional to the ratio of experimental to analytical steady-state lift curve slope. In addition to presenting the mathematical models and a brief description of existing analytical techniques, an alternative analytical technique for obtaining closed-loop results is presented.

Newsom, J. R.↗

Mass distribution and performance of free flight models

This note deals with the mass distribution and performance of free flight models. An airplane model which is to be used in free flight tests must be balanced dynamically as well as statically, e.g., it must not only have a given weight and the proper center of gravity but also a given ellipsoid of inertia. Equations which relate the motions of an airplane and its model are given. Neglecting scale effect, these equations may be used to predict the performance of an airplane, under the action of gravity alone, from data obtained in making dropping tests of a correctly balanced model.

Scherberg, Max↗

Physics-based, neural network force fields for reactive molecular dynamics: Investigation of carbene formation from [EMIM + ][OAc - ]

Reactive molecular dynamics simulations enable detailed understanding of solvent effects on chemical reaction mechanisms and reaction rates. While classical molecular dynamics using reactive force fields allows significantly longer simulation time scales and larger system sizes compared with ab initio molecular dynamics, constructing reactive force fields is a difficult and complex task. In this work, we describe a general approach following the Empirical Valence Bond (EVB) framework for constructing ab initio reactive force fields for condensed phase simulations by combining physics-based methods with neural networks (PB/NN). The physics-based terms ensure correct asymptotic behavior of electrostatic, polarization, and dispersion interactions, and are compatible with existing solvent force fields. Neural networks are utilized for versatile description of short-range orbital interactions within the transition state region, and accurate rendering of vibrational motion of the reacting complex. Herein, we demonstrate our methodology for a simple deprotonation reaction of the 1-ethyl-3-methylimidazolium (EMIM+) cation with acetate to form 1-ethyl- 3-methylimidazol-2-ylidene and acetic acid. Our PB/NN force field exhibits ~ 1 kJ/mol MAE accuracy within the transition state region for the gas-phase complex. To characterize solvent modulation of the reaction profile, we compute potentials of mean force (PMFs) for the gas-phase reaction as well as the reaction within a four ion cluster, and benchmark against ab initio molecular dynamics simulations. We find that the surrounding ionic environment significantly destabilizes formation of the carbene product, and we show that this effect is accurately captured by the reactive force field. By construction, the PB/NN potential may be directly employed for simulations of other solvents/chemical environments without additional parameterization.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Orbit improvement from satellite imaging data obtainable from the Mariner/Jupiter/Saturn 1977 missions

Equations of motion are established for a dynamic system in which a spacecraft flies close to and interacts with an outer planet and one or more of its satellites. For the computation of the state and mass partials needed in the orbit corrections, a set of variational equations is derived. These assume a notably compact form through the introduction of a matrix operator D. The above system of differential equations is integrated numerically on a computer.

Aksnes, K.↗

Current collection in a magnetoplasma

The authors present a survey of a very incomplete subject, current collection in a magnetoplasma. The best-developed and simplest theories for current collection are steady-state collisionless theories, and these must be understood before departures from them can be analyzed usefully. Thus, the authors begin with a review of them. The authors include some recent numerical results which indicate that steady-state collisionless Laplace-limit currents remain substantially below the Parker-Murphy (1967) canonical upper bound out to very large electrode potentials, and approach it as a limit only very slowly if at all. Attempts to correct this theory for space-charge effects lead to potential disturbances which extend to infinite distance along the electrode's magnetic shadow, unless collisional effects are also taken into account. However, even a small amount of relative plasma drift motion, such as that involved in a typical rocket experiment, can change this conclusion fundamentally. It is widely believed that time-averaged current collection may be increased by effects of plasma turbulence, and the authors review the available evidence for and against this contention. Steady-state collisionless particle dynamics predicts the existence of a toroidal region of trapped orbits which surrounds the electrode. Light emissions from this region have been photographed, indicating that collisional ionization may also occur there, and this, and/or scattering by collisions or possibly turbulent fluctuations in this region, may also increase current collection by the electrode. The authors also discuss effects on particle motions near the electrode, associated with breakdown of magnetic insulation in the region of large electric fields near it.

Laframboise, Jim G.↗

On motion of fluid in boundary layer near line of intersection of two planes

In the paper "The Mutual Interference of Boundary Layers," the authors investigated the problem of the interference of two planes intersecting at right angles on the boundary layers formed by the motion of fluid along the line of intersection of these planes. In the present paper, the results of the preceding one are generalized to the case of planes intersecting at any angle. The motion of a fluid in an angle less than 180 degrees is discussed and the enlargement of the boundary layers near the line of intersection of the planes, the limits of the interference effects of the boundary layers, and the corrections on the drag are determined. All computations are conducted by the Karman-Pohlhausen method for laminar and turbulent boundary layers. The results are reduced to tabulated form.

BOUNDARY LAYER - PLATES↗

High-Performance Parallel Analysis of Coupled Problems for Aircraft Propulsion

This research program dealt with the application of high-performance computing methods to the numerical simulation of complete jet engines. The program was initiated in January 1993 by applying two-dimensional parallel aeroelastic codes to the interior gas flow problem of a bypass jet engine. The fluid mesh generation, domain decomposition and solution capabilities were successfully tested. Attention was then focused on methodology for the partitioned analysis of the interaction of the gas flow with a flexible structure and with the fluid mesh motion driven by these structural displacements. The latter is treated by a ALE technique that models the fluid mesh motion as that of a fictitious mechanical network laid along the edges of near-field fluid elements. New partitioned analysis procedures to treat this coupled three-component problem were developed during 1994 and 1995. These procedures involved delayed corrections and subcycling, and have been successfully tested on several massively parallel computers, including the iPSC-860, Paragon XP/S and the IBM SP2. For the global steady-state axisymmetric analysis of a complete engine we have decided to use the NASA-sponsored ENG10 program, which uses a regular FV-multiblock-grid discretization in conjunction with circumferential averaging to include effects of blade forces, loss, combustor heat addition, blockage, bleeds and convective mixing. A load-balancing preprocessor tor parallel versions of ENG10 was developed. During 1995 and 1996 we developed the capability tor the first full 3D aeroelastic simulation of a multirow engine stage. This capability was tested on the IBM SP2 parallel supercomputer at NASA Ames. Benchmark results were presented at the 1196 Computational Aeroscience meeting.

Felippa, C. A.↗

Visual Acuity Using Head-fixed Displays During Passive Self and Surround Motion

The ability to read head-fixed displays on various motion platforms requires the suppression of vestibulo-ocular reflexes. This study examined dynamic visual acuity while viewing a head-fixed display during different self and surround rotation conditions. Twelve healthy subjects were asked to report the orientation of Landolt C optotypes presented on a micro-display fixed to a rotating chair at 50 cm distance. Acuity thresholds were determined by the lowest size at which the subjects correctly identified 3 of 5 optotype orientations at peak velocity. Visual acuity was compared across four different conditions, each tested at 0.05 and 0.4 Hz (peak amplitude of 57 deg/s). The four conditions included: subject rotated in semi-darkness (i.e., limited to background illumination of the display), subject stationary while visual scene rotated, subject rotated around a stationary visual background, and both subject and visual scene rotated together. Visual acuity performance was greatest when the subject rotated around a stationary visual background; i.e., when both vestibular and visual inputs provided concordant information about the motion. Visual acuity performance was most reduced when the subject and visual scene rotated together; i.e., when the visual scene provided discordant information about the motion. Ranges of 4-5 logMAR step sizes across the conditions indicated the acuity task was sufficient to discriminate visual performance levels. The background visual scene can influence the ability to read head-fixed displays during passive motion disturbances. Dynamic visual acuity using head-fixed displays can provide an operationally relevant screening tool for visual performance during exposure to novel acceleration environments.

Wood, Scott J.↗

The straw tracking detector for the Fermilab Muon g-2 Experiment

The Muon g-2 Experiment at Fermilab uses a gaseous straw tracking detector to make detailed measurements of the stored muon beam profile, which are essential for the experiment to achieve its uncertainty goals. Positrons from muon decays spiral inward and pass through the tracking detector before striking an electromagnetic calorimeter. The tracking detector is therefore located inside the vacuum chamber in a region where the magnetic field is large and non-uniform. As such, the tracking detector must have a low leak rate to maintain a high-quality vacuum, must be non-magnetic so as not to perturb the magnetic field and, to minimize energy loss, must have a low radiation length. The performance of the tracking detector has met or surpassed the design requirements, with adequate electronic noise levels, an average straw hit resolution of (110 ± 20) μm, a detection efficiency of 97% or higher, and no performance degradation or signs of aging. The tracking detector's measurements result in an otherwise unachievable understanding of the muon's beam motion, particularly at early times in the experiment's measurement period when there are a significantly greater number of muons decaying. This is vital to the statistical power of the experiment, as well as facilitating the precise extraction of several systematic corrections and uncertainties. This paper describes the design, construction, testing, commissioning, and performance of the tracking detector.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Orbit improvement from satellite imaging data obtainable from outer planet missions.

Equations of motion are established for a dynamical system in which a spacecraft flies close to and interacts with an outer planet and one or more of its satellites. For the computation of the state and mass partials needed in a simultaneous orbit correction of n interacting bodies, a notably compact set of variational equations is derived. The above system of differential equations is integrated numerically on a computer. Spacecraft-satellite direction measurements accurate to plus or minus 10 sec were simulated along three representative trajectories (Mariner/Jupiter/Saturn 1977 missions) approaching Io, Titan, and Iapetus to within 41000, 13000, and 7000 km, respectively. The paper concludes with a brief discussion of the need for future work on the orbits of the satellites of the outer planets.

Aksnes, K.↗

Unsteady Newton-Busemann flow theory. Part 2: Bodies of revolution

Newtonian flow theory for unsteady flow past oscillating bodies of revolution at very high Mach numbers is completed by adding a centrifugal force correction to the impact pressures. Exact formulas for the unsteady pressure and the stability derivatives are obtained in closed form and are applicable to bodies of revolution that have arbitrary shapes, arbitrary thicknesses, and either sharp or blunt noses. The centrifugal force correction arising from the curved trajectories followed by the fluid particles in unsteady flow cannot be neglected even for the case of a circular cone. With this correction, the present theory is in excellent agreement with experimental results for sharp cones and for cones with small nose bluntness; gives poor agreement with the results of experiments in air for bodies with moderate or large nose bluntness. The pitching motions of slender power-law bodies of revulution are shown to be always dynamically stable according to Newton-Busemann theory.

Hui, W. H.↗

The compressible mixing layer

The capability of turbulence modeling correctly to handle natural unsteadiness appearing in compressible turbulent flows is investigated. Physical aspects linked to the unsteadiness problem and the role of various flow parameters are analyzed. It is found that unsteady turbulent flows can be simulated by dividing these motions into an 'organized' part for which equations of motion are solved and a remaining 'incoherent' part represented by a turbulence model. Two-equation turbulence models and second-order turbulence models can yield reasonable results. For specific compressible unsteady turbulent flow, graphic presentations of different quantities may reveal complementary physical features. Strong compression zones are observed in rapid flow parts but shocklets do not yet occur.

Vandromme, Dany↗

Short-term adaptation of the VOR: non-retinal-slip error signals and saccade substitution

We studied short-term (30 min) adaptation of the vestibulo-ocular reflex (VOR) in five normal humans using a "position error" stimulus without retinal image motion. Both before and after adaptation a velocity gain (peak slow-phase eye velocity/peak head velocity) and a position gain (total eye movement during chair rotation/amplitude of chair motion) were measured in darkness using search coils. The vestibular stimulus was a brief ( approximately 700 ms), 15 degrees chair rotation in darkness (peak velocity 43 degrees /s). To elicit adaptation, a straight-ahead fixation target disappeared during chair movement and when the chair stopped the target reappeared at a new location in front of the subject for gain-decrease (x0) adaptation, or 10 degrees opposite to chair motion for gain-increase (x1.67) adaptation. This position-error stimulus was effective at inducing VOR adaptation, though for gain-increase adaptation the primary strategy was to substitute augmenting saccades during rotation while for gain-decrease adaptation both corrective saccades and a decrease in slow-phase velocity occurred. Finally, the presence of the position-error signal alone, at the end of head rotation, without any attempt to fix upon it, was not sufficient to induce adaptation. Adaptation did occur, however, if the subject did make a saccade to the target after head rotation, or even if the subject paid attention to the new location of the target without actually looking at it.

Non-NASA Center↗

MISR Level 2 TOA/Cloud Stereo parameters (MIL2TCST_V1)

The MISR Top-of-Atmosphere (TOA)/Cloud Stereo geophysical parameters include stereoscopically-derived cloud motion vectors (winds), cloud-top heights, and an accompanying cloud mask. The Stereo product geophysical parameters include a stereoscopically-derived cloud mask and cloud height on a 1.1 km grid. It also includes cloud motion vectors on a 70.4 km grid. The three types of stereo heights are: the BestWind heights are only calculated for those regions where the associated wind vectors passed the quality tests. Therefore, they have sparse coverage but since the wind correction is included, these contain our 'best guess' as to what the true heights are. The WithoutWind heights are calculated assuming a constant wind vector of zero. They have almost complete coverage and therefore form a nice 'pretty picture' of the relative cloud heights over small areas. The RawWind heights are a diagnostic product as they are calculated using all available wind vectors (even the bad ones). It is therefore recommended that one only use the Best and Without wind products. It is important to remember that the stereo matchers pick up the layer of maximum contrast, which is not necessarily the same as the highest cloud so all the stereo heights are keyed to this level of maximum contrast. Therefore, higher and thinner cirrus layers may not be detected by any of the height fields. [Location=GLOBAL] [Temporal_Coverage: Start_Date=2000-02-24; Stop_Date=] [Spatial_Coverage: Southernmost_Latitude=-90; Northernmost_Latitude=90; Westernmost_Longitude=-180; Easternmost_Longitude=180] [Data_Resolution: Latitude_Resolution=1.1 km; Longitude_Resolution=1.1 km; Horizontal_Resolution_Range=1 km - < 10 km or approximately .01 degree - < .09 degree; Temporal_Resolution=about 15 orbits/day; Temporal_Resolution_Range=Daily - < Weekly].

EOSDIS↗

MISR Level 2 TOA/Cloud Stereo parameters (MIL2TCST_V2)

The MISR Top-of-Atmosphere (TOA)/Cloud Stereo geophysical parameters include stereoscopically-derived cloud motion vectors (winds), cloud-top heights, and an accompanying cloud mask. The Stereo product geophysical parameters include a stereoscopically-derived cloud mask and cloud height on a 1.1 km grid. It also includes cloud motion vectors on a 70.4 km grid. The three types of stereo heights are: the BestWind heights are only calculated for those regions where the associated wind vectors passed the quality tests. Therefore, they have sparse coverage but since the wind correction is included, these contain our 'best guess' as to what the true heights are. The WithoutWind heights are calculated assuming a constant wind vector of zero. They have almost complete coverage and therefore form a nice 'pretty picture' of the relative cloud heights over small areas. The RawWind heights are a diagnostic product as they are calculated using all available wind vectors (even the bad ones). It is therefore recommended that one only use the Best and Without wind products. It is important to remember that the stereo matchers pick up the layer of maximum contrast, which is not necessarily the same as the highest cloud so all the stereo heights are keyed to this level of maximum contrast. Therefore, higher and thinner cirrus layers may not be detected by any of the height fields. [Temporal_Coverage: Start_Date=2000-02-24; Stop_Date=] [Spatial_Coverage: Southernmost_Latitude=-90; Northernmost_Latitude=90; Westernmost_Longitude=-180; Easternmost_Longitude=180] [Data_Resolution: Latitude_Resolution=1.1 km; Longitude_Resolution=1.1 km; Temporal_Resolution=about 15 orbits/day].

CLOUD MASK↗