Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “peak estimation”

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 1,297 records · Page 72

OGLE-2016-BLG-0156: Microlensing Event with Pronounced Microlens-parallax Effects Yielding a Precise Lens Mass Measurement

We analyze the gravitational binary-lensing event OGLE-2016-BLG-0156, for which the lensing light curve displays pronounced deviations induced by microlens-parallax effects. The light curve exhibits three distinctive widely separated peaks and we find that the multiple-peak feature provides a very tight constraint on the microlens-parallax effect, enabling us to precisely measure the microlens parallax . All the peaks are densely and continuously covered from high-cadence survey observations using globally located telescopes and the analysis of the peaks leads to the precise measurement of the angular Einstein radius . From the combination of the measured and , we determine the physical parameters of the lens. It is found that the lens is a binary composed of two M dwarfs with masses M(1) = 0.18 ± 0.01 M(⊙) and M(2) = 0.16 ± 0.01 M(⊙) located at a distance D(L) = 1.35 ± 0.09 kpc. According to the estimated lens mass and distance, the flux from the lens comprises an important fraction, ~25%, of the blended flux. The bright nature of the lens combined with the high relative lens-source motion, μ = 6.94 ± 0.50 mas/yr, suggests that the lens can be directly observed from future high-resolution follow-up observations.

Youn Kil Jung↗

Exploring the Predictability of the High-Energy Tail of MEE Precipitation Based on Solar Wind Properties

Medium Energy Electron (MEE) precipitation (≳30 keV) ionizes the mesosphere and initiates chemical reactions, which ultimately can reduce mesospheric and stratospheric ozone. Currently, there are considerable differences in how existing parameterizations represent flux response, timing, and duration of MEE precipitation, especially considering its high-energy tail (≳300 keV). This study compares the nature of ≳300 to ≳30 keV electron fluxes to better understand differences within MEE precipitation. The MEE fluxes are estimated from measurements by the Medium Energy Proton and Electron Detector (MEPED) onboard the Polar Orbiting Environmental Satellite (POES) from 2004 to 2014. The fluxes are explored in the context of solar wind drivers: corotating high-speed solar wind streams (HSSs) and coronal mass ejections (CMEs) alongside their associated solar wind properties. Three key aspects of ≳300 keV electron fluxes are investigated: maximum response, peak timing, and duration. The results reveal a structure-dependent correlation (0.89) between the peak fluxes of ≳30 and ≳300 keV electrons. The epsilon coupling function correlates well (0.84) with the ≳300 keV peak flux, independent of solar wind structure. The ≳300 keV flux peaks 0–3 days after the ≳30 keV flux peaks. The highest probability (∼42%) occurs for a 1-day delay, while predictive capabilities increase when accounting for solar wind speed. The ≳300 keV flux response has the highest probability of lasting 4 days for both CMEs and HSSs. The results form a base for a stochastic MEE parameterization that goes beyond the average picture, enabling realistic flux variability on both daily and decadal scales.

J A Salice↗

Trunk muscle forces and spinal loads while walking in persons with lower limb amputation: Influences of chronic low back pain

Persons with lower limb amputation (LLA) are at high risk for developing chronic low back pain (LBP), often with biomechanical factors considered as likely contributors. Here, trunk and pelvis kinematics, muscle forces, and resultant spinal loads were characterized in persons with LLA, with and without chronic LBP. Thirty-five persons with unilateral LLA – 19 with chronic LBP (“LLA-cLBP”), 16 without LBP (“LLA-nLBP”) – and 15 (uninjured) persons without LBP (“CTR-nLBP”) walked overground (1.3 m/s) while thorax and pelvis kinematics were tracked (and ranges of motion [ROM] computed), and used as inputs for a non-linear finite element model of the spine to estimate global and local muscle forces, and resultant spinal loads. In the frontal and transverse planes, thorax ROM were up to 66.6% smaller in LLA-nLBP versus LLA-cLBP (P < 0.001) and CTR-nLBP (P < 0.001). In the sagittal plane, pelvis ROM was 50.4% smaller in LLA-nLBP versus LLA-cLBP (P = 0.014). LLA-cLBP exhibited 45.5% and 34.2% greater peak local and global muscle forces, respectively, versus CTR-nLBP (P < 0.011). Up to 48.1% greater spinal loads were observed in LLA-cLBP versus CTR-nLBP (P < 0.013); peak compression and local muscle forces were respectively 20.2% and 41.0% larger in LLA-nLBP versus CTR-nLBP (P < 0.005). Despite differences in trunk and pelvis kinematics between LLA-cLBP and LLA-nLBP, trunk muscle forces and spinal loads were similar (P > 0.101) between these groups. Similar loading parameters regardless of LBP presence, while highly dependent on trunk muscle activation strategies, may mitigate further accumulation of mechanical fatigue. Finally, it remains important to understand the temporality of loading with respect to LBP onset following LLA.

60 APPLIED LIFE SCIENCES↗

Time-Resolved X-ray Operando Observations of Lithiation Gradients across the Cathode Matrix and Individual Oxide Particles during Fast Cycling of a Li-Ion Cell

Lithiated transition metal oxides serve as active materials in the positive electrode (cathode) of lithium-ion cells. During electrochemical cycling, lithium ions intercalate and deintercalate into these oxide particles. This behavior causes two types of lithiation gradients to emerge: (i) a bulk gradient across the depth of the cathode matrix (averaged over individual oxide particles) and (ii) a microscopic gradient across the particles themselves, which also depends on their location in the electrode. Here we show how both gradients can be studied using operando X-ray diffraction during 4C charge and 4C discharge. The oxide (de)lithiation is estimated from the unit cell parameters by indexing the X-ray diffraction spectra. By fitting the lithiation profiles with orthogonal polynomials, the bulk gradients across the electrode thickness are quantified. These gradients develop as the current flows through the cell and dissipate during open-circuit and potentiostatic-hold periods. Further details of lithiation dynamics can be obtained through shape analysis of the Bragg peaks. In particular, from electrochemical model simulations, we show that the width and skewness of the (003) peak track (de)lithiation fronts moving across the individual oxide particles.

(003) peak↗

Design of infrasound-detection system via adaptive LMSTDE algorithm

A proposed solution to an aviation safety problem is based on passive detection of turbulent weather phenomena through their infrasonic emission. This thesis describes a system design that is adequate for detection and bearing evaluation of infrasounds. An array of four sensors, with the appropriate hardware, is used for the detection part. Bearing evaluation is based on estimates of time delays between sensor outputs. The generalized cross correlation (GCC), as the conventional time-delay estimation (TDE) method, is first reviewed. An adaptive TDE approach, using the least mean square (LMS) algorithm, is then discussed. A comparison between the two techniques is made and the advantages of the adaptive approach are listed. The behavior of the GCC, as a Roth processor, is examined for the anticipated signals. It is shown that the Roth processor has the desired effect of sharpening the peak of the correlation function. It is also shown that the LMSTDE technique is an equivalent implementation of the Roth processor in the time domain. A LMSTDE lead-lag model, with a variable stability coefficient and a convergence criterion, is designed.

Khalaf, C. S.↗

Comparison of Southern Hemisphere radiosonde and LIMS temperatures at 100 mb

Radiosonde (RS)and satellite-derived (Nimbus-7 LIMS) 100-mb temperatures over New Zealand at 12 GMT are compared for the 1978-79 summer. The colocated LIMS temperature information consists of synoptically mapped values (for 12 GMT), as well as the primary nighttime orbital retrievals valid at about 1030 GMT. The RS time series of temperature is dominated by temporal fluctuations associated mainly with the eastward passage of waves which have characteristic periods of 4-5 and 11-12 days and peak-to-peak amplitudes of 10-15 K. The LIMS temperatures and the corresponding temperature time series are also found to exhibit quite close agreement (in terms of temporal phase for the latter) with the RS data. However, the LIMS-mapped temperature fluctuations suffer from a noticeable attenuation in amplitude (approaching 50 percent for higher-frequency fluctuations), which will affect the accuracy of LIMS-derived estimates of dynamical quantities such as wind velocity and relative vorticity in the lower stratosphere.

Miles, T.↗

Magnitude and Seasonality of Wetland Methane Emissions from the Hudson Bay Lowlands (Canada)

The Hudson Bay Lowlands (HBL) is the second largest boreal wetland ecosystem in the world and an important natural source of global atmospheric methane. We quantify the HBL methane emissions by using the GEOS-Chem chemical transport model to simulate aircraft measurements over the HBL from the ARCTAS and pre-HIPPO campaigns in May-July 2008, together with continuous 2004-2008 surface observations at Fraserdale (southern edge of HBL) and Alert (Arctic background). The difference in methane concentrations between Fraserdale and Alert is shown to be a good indicator of HBL emissions, and implies a sharp seasonal onset of emissions in late May (consistent with the aircraft data), a peak in July-August, and a seasonal shut-off in September. The model, in which seasonal variation of emission is mainly driven by surface temperature, reproduces well the observations in summer but its seasonal shoulders are too broad. We suggest that this reflects the suppression of emissions by snow cover and greatly improve the model simulation by accounting for this effect. Our resulting best estimate for HBL methane emissions is 2.3 Tg/a, several-fold higher than previous estimates (Roulet et al., 1994; Worthy et al., 2000).

Pickett-Heaps, C. A.↗

Role of MODIS Vegetation Phenology Products in the ForWarn System for Monitoring of Forest Disturbances in the Conterminous United States

This presentation discusses MODIS vegetation phenology products used in the ForWarn Early Warning System (EWS) tool for near real time regional forest disturbance detection and surveillance at regional to national scales. The ForWarn EWS is being developed by the USDA Forest Service NASA, ORNL, and USGS to aid federal and state forest health management activities. ForWarn employs multiple historical land surface phenology products that are derived from MODIS MOD13 Normalized Difference Vegetation Index (NDVI) data. The latter is temporally processed into phenology products with the Time Series Product Tool (TSPT) and the Phenological Parameter Estimation Tool (PPET) software produced at NASA Stennis Space Center. TSPT is used to effectively noise reduce, fuse, and void interpolate MODIS NDVI data. PPET employs TSPT-processed NDVI time series data as an input, outputting multiple vegetation phenology products at a 232 meter resolution for 2000 to 2011, including NDVI magnitude and day of year products for seven key points along the growing season (peak of growing season and the minima, 20%, and 80% of the peak NDVI for both the left and right side of growing season), cumulative NDVI integral products for the most active part of the growing season and sequentially across the growing season at 8 day intervals, and maximum value NDVI products composited at 24 day intervals in which each product date has 8 days of overlap between the previous and following product dates. MODIS NDVI phenology products are also used to compute nationwide NRT forest change products refreshed every 8 days. These include percent change in forest NDVI products that compare the current NDVI from USGS eMODIS products to historical MODIS MOD13 NDVI. For each date, three forest change products are produced using three different maximum value NDVI baselines (from the previous year, three previous years, and all previous years). All change products are output with a rainbow color table in which forests with the most severe NDVI decreases are assigned hot colors (yellow to red) and forests with prominent NDVI increases are assigned cold colors (blue tones). All mentioned products have been integrated as data layers into ForWarn s geospatial data viewer known as the U.S. Forest Change Assessment Viewer (FCAV). The latter is used to view and assess the context of the mentioned forest change products with respect to ancillary data layers, such as land cover, elevation, hydrologic features, climatic data, storm data, aerial disturbance surveys, fire data, and land ownership. The FCAV also includes a temporal NDVI profiler for viewing phenological change in multi-year NDVI associated with known or suspected regionally apparent forest disturbances (e.g., from fire and insects). ForWarn forest change products have been used to detect, track, and assess several biotic and abiotic regional forest disturbance events across the country, including ephemeral and longer lasting damage from storms, drought, and insects. Such change products are most effective for viewing severe disturbances affecting multiple MODIS pixels. MODIS vegetation phenology products contribute vital current information on forest conditions to the ForWarn system and this role is expected to grow as these products are refined and derivative products are added.

Spruce, Joseph P.↗

Isolating Added Mass Load Components of CPAS Main Clusters

The current simulation for the Capsule Parachute Assembly System (CPAS) lacks fidelity in representing added mass for the 116 ft Do ringsail Main parachute. The availability of 3-D models of inflating Main canopies allowed for better estimation the enclosed air volume as a function of time. This was combined with trajectory state information to estimate the components making up measured axial loads. A proof-of-concept for an alternate simulation algorithm was developed based on enclosed volume as the primary independent variable rather than drag area growth. Databases of volume growth and parachute drag area vs. volume were developed for several flight tests. Other state information was read directly from test data, rather than numerically propagated. The resulting simulated peak loads were close in timing and magnitude to the measured loads data. However, results are very sensitive to data curve fitting and may not be suitable for Monte Carlo simulations. It was assumed that apparent mass was either negligible or a small fraction of enclosed mass, with little difference in results.

Ray, Eric S.↗

Comparison of scattering ratio profiles retrieved from ALADIN/Aeolus and CALIOP/CALIPSO observations and preliminary estimates of cloud fraction profiles

The space-borne active sounders have been contributing invaluable vertically resolved information of atmospheric optical properties since the launch of Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observation (CALIPSO) in 2006. To build long-term records from space-borne lidars useful for climate studies, one has to understand the differences between successive space lidars operating at different wavelengths, flying on different orbits, and using different viewing geometries, receiving paths, and detectors. In this article, we compare the results of Atmospheric Laser Doppler INstrument (ALADIN) and Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) lidars for the period from 28 June to 31 December 2019. First, we build a dataset of ALADIN–CALIOP collocated profiles (Δdist<1°; Δtime<6 h). Then we convert ALADIN's 355 nm particulate backscatter and extinction profiles into the scattering ratio vertical profiles SR(z) at 532 nm using molecular density profiles from Goddard Earth Observing System Data Assimilation System, version 5 (GEOS-5 DAS). And finally, we build the CALIOP and ALADIN globally gridded cloud fraction profiles CF(z) by applying the same cloud detection threshold to the SR(z) profiles of both lidars at the same spatial resolution. Before comparing the SR(z) and CF(z) profiles retrieved from the two analyzed lidar missions, we performed a numerical experiment to estimate the best achievable cloud detection agreement CDA norm (z) considering the differences between the instruments. We define CDA norm (z) in each latitude–altitude bin as the occurrence frequency of cloud layers detected by both lidars, divided by a cloud fraction value for the same latitude–altitude bin. We simulated the SR(z) and CF(z) profiles that would be observed by these two lidars if they were flying over the same atmosphere predicted by a global model. By analyzing these simulations, we show that the theoretical limit for CDA$_{norm}^{theor}$(z) for a combination of ALADIN and CALIOP instruments is equal to 0.81±0.07 at all altitudes. In other words, 19 % of the clouds cannot be detected simultaneously by two instruments due to said differences. The analyses of the actual observed CALIOP–ALADIN collocated dataset containing ~78 000 pairs of nighttime SR(z) profiles revealed the following points: (a) the values of SR(z) agree well up to ~3 km height. (b) The CF(z) profiles show agreement below ~3 km, where ~80 % of the clouds detected by CALIOP are detected by ALADIN as expected from the numerical experiment. (c) Above this height, the CDA$_{norm}^{obs}$(z) reduces to ~50 %. (d) On average, better sensitivity to lower clouds skews ALADIN's cloud peak height in pairs of ALADIN–CALIOP profiles by ~0.5±0.6 km downwards, but this effect does not alter the heights of polar stratospheric clouds and high tropical clouds thanks to their strong backscatter signals. (e) The temporal evolution of the observed CDA$_{norm}^{obs}$(z) does not reveal any statistically significant change during the considered period. This indicates that the instrument-related issues in ALADIN L0/L1 have been mitigated, at least down to the uncertainties of the following CDA$_{norm}^{obs}$(z) values: 68±12 %, 55±14 %, 34±14 %, 39±13 %, and 42±14 % estimated at 0.75, 2.25, 6.75, 8.75, and 10.25 km, respectively.

54 ENVIRONMENTAL SCIENCES↗

Efficient excimer-based white OLEDs with reduced efficiency roll-off

A tetradentate Pt(II) complex, Pt2O2-p2m, which exhibits both efficient monomer and excimer emission, was synthesized. Through rational device design to broaden the recombination zone, an OLED device employing Pt2O2-p2m demonstrated an impressively high peak external quantum efficiency (EQE) of 21.5% and a warm white emission with a color rendering index of 71 and Commission Internationale de L'Eclairage coordinates of (0.409, 0.458). Moreover, a reduced efficiency roll-off was also observed to retain high EQE values of 20% EQE at 1000 cd/m 2 and 17% EQE at 5000 cd/m 2 , and an LT70 (time to 70% initial luminance) was estimated to be over 150 h at 1000 cd/m2 under constant operation at accelerated driving conditions.

36 MATERIALS SCIENCE↗

Wind-driven emissions of coarse-mode particles in an urban environment

Abstract. Quantifying surface–atmosphere exchange rates of particles is important for understanding the role of suspended particulate matter in radiative transfer, clouds, precipitation, and climate change. Emissions of coarse-mode particles with a diameter greater than 0.5 µm provide giant cloud condensation nuclei and ice nuclei. These emissions are critical for understanding the evolution of cloud microphysical properties yet remain poorly understood. Here we introduce a new method that uses lidar retrievals of the elastic backscatter and Doppler velocity to obtain surface number emissions of particles with a diameter greater than 0.53 µm. The technique is applied to study particle number fluxes over a 2-month period from 1 June to 10 August 2022 during the TRACER campaign at an urban site near Houston, TX, USA. We found that all the observed fluxes were positive (upwards), indicating particle emission from the surface. The fluxes followed a diurnal pattern and peaked near noon local time. Flux intensity varied through the 2 months with multi-day periods of strong fluxes and multi-day periods of weak fluxes. Emission particle number fluxes peaked near ∼ 100 cm−2 s−1. The daily averaged emission fluxes correlated with friction velocity and were anticorrelated with surface relative humidity. The emission flux can be parameterized as F= 3000 u*4, where u* is the friction velocity in m s−1 and the emission flux F is in cm−2 s−1. The u* dependence is consistent with emission from wind-driven erosion. Estimated values for the mass flux are in the lower range of literature values from non-urban sites. These results demonstrate that urban environments may play an important role in supplying coarse-mode particles to the boundary layer. We anticipate that quantification of these emissions will help constrain aerosol–cloud interaction models that use prognostic aerosol schemes.

54 ENVIRONMENTAL SCIENCES↗

The energetics of molecular clouds. I - Methods of analysis and application to the S255 molecular cloud

Observations of the S255 molecular cloud in the radio lines of CO, (C-13)O, and H2O as well as in the IR region from 12 to 20 microns are presented. These results show that an extended molecular cloud is associated with the optically visible H II regions S255 - S257, both the kinetic temperature and molecular density are enhanced in an area centered on a compact near-IR source, this source is also coincident with an OH and an H2O maser (within the errors), and the IR source exhibits the silicon absorption feature. First-order techniques which allow physical properties of molecular clouds to be estimated and the energetics of the gas and dust in the clouds to be analyzed are developed and applied to S255. It is concluded that: (1) the gas cooling rate is much less than the dust cooling rate; (2) the primary energy flow is in the IR through dust emission over a range of temperatures; (3) the near-IR source is apparently responsible for the local peak in gas kinetic temperature; and (4) the exciting stars of nearby H II regions appear to be the primary heat sources for the overall molecular cloud.

Evans, N. J., II↗

Solar potential inventory and modeling

Image processing procedures for calculating the energy that roof-mounted solar collectors can potentially supply in a metropolitan area are presented. Satellite multispectral imagery from which land cover types can be determined digitally was sampled in order to estimate the percentage of land area occupied by flat or south-facing roof tops in residential and commercial/industrial areas. Procedures were applied to the various power subdistricts of the western San Fernando valley of California, and it was found that on the average 120% of the existing power demand could be met if only half the useable rooftop area were utilized, amounting to 385 MW of peak power and indicating the applicability of solar cells to power generation in urban areas.

Angelici, G. L.↗

A comparison of ship and Coastal Zone Color Scanner mapped distribution of phytoplankton in the southeastern Bering Sea

Twenty-one Coastal Zone Color Scanner (CZCS) images of the southeastern Bering Sea are examined in order to map the near-surface distribution of phytoplankton during 1979 and 1980. The information is compared with the mesoscale (100-1000 km) distribution of phytoplankton inferred from pooled ship sampling obtained during the Processes and Resources of the Bering Shelf (PROBES) intensive field study during the late 1970s and early 1980s. The imagery indicates that open-water phytoplankton blooms occur first in late April in coastal waters, peak in early May over the middle shelf, and decay rapidly afterwards, reaching concentration minima in June in both regions. These patterns show that the earlier ship observations are valid for most of the eastern Bering shelf. A very tight correlation is found between the PROBES surface chlorophyll a concentrations and mean mixed-layer chlorophyll concentrations. The significant discrepancies between CZCS and ship-based chlorophyll estimates may be due to aliasing in time by the CZCS. It is concluded that neither satellite nor ship alone can do an adequate job of characterizing the physics or biological dynamics of the ocean.

Mcclain, C. R.↗

Convective lofting links Indian Ocean air pollution to paradoxical South Atlantic ozone maxima

We describe a broad resolution of the Atlantic Parado concerning the seasonal and geographic distribution, of tropical tropospheric ozone. We highlight periods of significant maximum tropospheric O3 for Jan.- April, 1999, exploiting satellite estimates and SHADOZ (Southern Hemisphere Additional Ozonesondes). Trajectory analyses connecting sondes and Total Tropospheric Ozone (TTO) maps suggest a complex influence from the Indian Ocean: beginning with mixed combustion sources, then low level transport, cumulonimbus venting, possible stratospheric input, and finally high-level transport to the west, with possible mixing over Africa. For the Jan.-March highest column-O3 periods in the Atlantic, distinct sounding peaks trace to specific NO sources, especially lightning, while in the same episodes, recurring every 20-50 days, more diffuse buildups of Indian-to-Atlantic pollution make important contributions.

Chatfield, R. B.↗

Siderophile Element Depletion in the Angrite Parent Body (APB) Mantle: Due to Core Formation?

The origin of angrites has evaded scientists due in part to unusual mineralogy, oxidized character, and small numbers of samples. Increased interest in the origin of angrites has stemmed from the recovery of approximately 10 new angrites in the past decade. These new samples have allowed meteoriticists to recognize that angrites are compositionally diverse, old, and record very early differentiation. Also, a magma ocean has been proposed to have been involved in APB early differentiation, but this remains untested for siderophile elements which are commonly cited as one of the main lines of evidence for magma oceans on the early Earth, Moon, Mars and eucrite parent body (e.g., [6]). And recent suggestions that angrites may or may not be from Mercury have also peaked interest in these achondrites. Given all of this background, a detailed understanding of the early differentiation process is desired. Previous efforts at examining siderophile element (SE) concentrations with respect to core formation processes in the APB have not resulted in any definite conclusions regarding segregation of a metallic core. The goal of this study is to summarize what is known about SE concentrations in the suite, estimate depletions of SE compared to chondrites, and apply metal/silicate experimental partition coefficients to assess whether the APB had a core.

Righter, K.↗

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↗