Search NASA⌕ Search

SEARCH · Search NASA

Results for “method”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
NASA NTRS1972

Interferometric spectropolarimetry - Alternate experimental methods.

Three alternate methods of obtaining spectra of the intensity and state of polarization of light are proposed. The methods make use of a two-beam amplitude division interferometer using the technique of Fourier spectroscopy. They do not require the presence of polarizers in the arms of the instrument. In one method a single analyzer is used in front of the detector with three successive orientations of its transmission axis azimuth. In another method, a (linear) polarizer assuming the same set of orientations is placed in the incident beam. A third method, a hybrid of the former two methods, makes use of both a polarizer and analyzer in the locations indicated. The first method permits the determination of all four Stokes parameters of polarization, whereas the other two methods cannot yield the ellipticity parameter. All methods require the recording of three interferograms. However, two interferograms can provide the intensity and degree of polarization in any of the methods described.

Fymat, A. L.↗
NASA NTRS1976

The Relation of Finite Element and Finite Difference Methods

Finite element and finite difference methods are examined in order to bring out their relationship. It is shown that both methods use two types of discrete representations of continuous functions. They differ in that finite difference methods emphasize the discretization of independent variable, while finite element methods emphasize the discretization of dependent variable (referred to as functional approximations). An important point is that finite element methods use global piecewise functional approximations, while finite difference methods normally use local functional approximations. A general conclusion is that finite element methods are best designed to handle complex boundaries, while finite difference methods are superior for complex equations. It is also shown that finite volume difference methods possess many of the advantages attributed to finite element methods.

Vinokur, M.↗
NASA NTRS1986

Methods for assessing wall interference in the 2- by 2-foot adaptive-wall wind tunnel

Discussed are two methods for assessing two-dimensional wall interference in the adaptive-wall test section of the NASA Ames 2 x 2-Foot Transonic Wind Tunnel: (1) a method for predicting free-air conditions near the walls of the test section (adaptive-wall methods); and (2) a method for estimating wall-induced velocities near the model (correction methods), both of which methods are based on measurements of either one or two components of flow velocity near the walls of the test section. Each method is demonstrated using simulated wind tunnel data and is compared with other methods of the same type. The two-component adaptive-wall and correction methods were found to be preferable to the corresponding one-component methods because: (1) they are more sensitive to, and give a more complete description of, wall interference; (2) they require measurements at fewer locations; (3) they can be used to establish free-stream conditions; and (4) they are independent of a description of the model and constants of integration.

Schairer, E. T.↗
NASA NTRS1987

Accuracies of southwell and force/stiffness methods in the prediction of buckling strength of hypersonic aircraft wing tubular panels

Accuracies of the Southwell method and the force/stiffness (F/S) method are examined when the methods were used in the prediction of buckling loads of hypersonic aircraft wing tubular panels, based on nondestructive buckling test data. Various factors affecting the accuracies of the two methods were discussed. Effects of load cutoff point in the nondestructive buckling tests on the accuracies of the two methods were discussed in great detail. For the tubular panels under pure compression, the F/S method was found to give more accurate buckling load predictions than the Southwell method, which excessively overpredicts the buckling load. It was found that the Southwell method required a higher load cutoff point, as compared with the F/S method. In using the F/S method for predicting the buckling load of tubular panels under pure compression, the load cutoff point of approximately 50 percent of the critical load could give reasonably accurate predictions.

Ko, William L.↗
NASA NTRS1988

An investigation of new methods for estimating parameter sensitivities

Parameter sensitivity is defined as the estimation of changes in the modeling functions and the design variables due to small changes in the fixed parameters of the formulation. There are currently several methods for estimating parameter sensitivities requiring either difficult to obtain second order information, or do not return reliable estimates for the derivatives. Additionally, all the methods assume that the set of active constraints does not change in a neighborhood of the estimation point. If the active set does in fact change, than any extrapolations based on these derivatives may be in error. The objective here is to investigate more efficient new methods for estimating parameter sensitivities when the active set changes. The new method is based on the recursive quadratic programming (RQP) method and in conjunction a differencing formula to produce estimates of the sensitivities. This is compared to existing methods and is shown to be very competitive in terms of the number of function evaluations required. In terms of accuracy, the method is shown to be equivalent to a modified version of the Kuhn-Tucker method, where the Hessian of the Lagrangian is estimated using the BFS method employed by the RPQ algorithm. Inital testing on a test set with known sensitivities demonstrates that the method can accurately calculate the parameter sensitivity. To handle changes in the active set, a deflection algorithm is proposed for those cases where the new set of active constraints remains linearly independent. For those cases where dependencies occur, a directional derivative is proposed. A few simple examples are included for the algorithm, but extensive testing has not yet been performed.

Beltracchi, Todd J.↗
NASA NTRS1991

Development of advanced modal methods for calculating transient thermal and structural response

Higher-order modal methods for predicting thermal and structural response are evaluated. More accurate methods or ones which can significantly reduce the size of complex, transient thermal and structural problems are desirable for analysis and are required for synthesis of real structures subjected to thermal and mechanical loading. A unified method is presented for deriving successively higher-order modal solutions related to previously-developed, lower-order methods such as the mode displacement and mode-acceleration methods. A new method, called the force-derivative method, is used to obtain higher-order modal solutions for both uncoupled (proportionally-damped) structural problems as well as thermal problems and coupled (non-proportionally damped) structural problems. The new method is called the force-derivative method because, analogous to the mode-acceleration method, it produces a term that depends on the forcing function and additional terms that depend on the time derivatives of the forcing function.

Camarda, Charles J.↗
NASA NTRS1997

Development of Improved Surface Integral Methods for Jet Aeroacoustic Predictions

The accurate prediction of aerodynamically generated noise has become an important goal over the past decade. Aeroacoustics must now be an integral part of the aircraft design process. The direct calculation of aerodynamically generated noise with CFD-like algorithms is plausible. However, large computer time and memory requirements often make these predictions impractical. It is therefore necessary to separate the aeroacoustics problem into two parts, one in which aerodynamic sound sources are determined, and another in which the propagating sound is calculated. This idea is applied in acoustic analogy methods. However, in the acoustic analogy, the determination of far-field sound requires the solution of a volume integral. This volume integration again leads to impractical computer requirements. An alternative to the volume integrations can be found in the Kirchhoff method. In this method, Green's theorem for the linear wave equation is used to determine sound propagation based on quantities on a surface surrounding the source region. The change from volume to surface integrals represents a tremendous savings in the computer resources required for an accurate prediction. This work is concerned with the development of enhancements of the Kirchhoff method for use in a wide variety of aeroacoustics problems. This enhanced method, the modified Kirchhoff method, is shown to be a Green's function solution of Lighthill's equation. It is also shown rigorously to be identical to the methods of Ffowcs Williams and Hawkings. This allows for development of versatile computer codes which can easily alternate between the different Kirchhoff and Ffowcs Williams-Hawkings formulations, using the most appropriate method for the problem at hand. The modified Kirchhoff method is developed primarily for use in jet aeroacoustics predictions. Applications of the method are shown for two dimensional and three dimensional jet flows. Additionally, the enhancements are generalized so that they may be used in any aeroacoustics problem.

Pilon, Anthony R.↗
NASA NTRS2000

Evaluation of Non-Ozone-Depleting-Chemical Cleaning Methods for Space Mechanisms Using a Vacuum Spiral Orbit Rolling Contact Tribometer

Because CFC 113, an ozone depleting chemical (ODC), can no longer be produced, alternative bearing cleaning methods must be studied. The objective of this work was to study the effect of the new cleaning methods on lubricant lifetime using a vacuum bearing simulator (spiral orbit rolling contact tribometer). Four alternative cleaning methods were studied: ultra-violet (UV) ozone, aqueous levigated alumina slurry (ALAS), super critical fluid (SCF) CO2 and aqueous Brulin 815GD. Baseline tests were done using CFC 113. Test conditions were the following: a vacuum of at least 1.3 x 10(exp -6) Pa, 440C steel components, a rotational speed of 10 RPM, a lubricant charge of between 60-75 micrograms, a perfluoropolyalkylether lubricant (Z-25), and a load of 200N (44.6 lbs., a mean Hertzian stress of 1.5 GPa). Normalized lubricant lifetime was determined by dividing the total number of ball orbits by the amount of lubricant. The failure condition was a friction coefficient of 0.38. Post-test XPS analysis was also performed, showing slight variations in post-cleaning surface chemistry. Statistical analysis of the resultant data was conducted and it was determined that the data sets were most directly comparable when subjected to a natural log transformation. The natural log life (NL-Life) data for each cleaning method were reasonably normally (statistically) distributed and yielded standard deviations that were not significantly different among the five cleaning methods investigated. This made comparison of their NL-Life means very straightforward using a Bonferroni multiple comparison of means procedure. This procedure showed that the ALAS, UV-ozone and CFC 113 methods were not statistically significantly different from one another with respect to mean NL-Life. It also found that the SCF CO2 method yielded a significantly higher mean NL-Life than the mean NL-Lives of the ALAS, UV-ozone and CFC 113 methods. It also determined that the aqueous Brulin 815GD method yielded a mean NL-Life that was statistically significantly higher than the mean NL-Lives of each of the other four methods. Baseline tests using CFC 113 cleaned parts yielded a mean NL-Life 3.62 orbits/micro-g. ALAS and UV-ozone yielded similar mean NL-Life (3.31 orbits/mg and 3.33 orbits/micro-g, respectively). SCF CO2, gave a mean NL-Life of 4.08 orbits/mg and aqueous Brulin 8l5GD data yielded the longest mean NL-Life (4.66 orbits/micro-g).

Jansen, Mark J.↗
NASA NTRS2001

A Unified Development of Basis Reduction Methods for Rotor Blade Analysis

The axial foreshortening effect plays a key role in rotor blade dynamics, but approximating it accurately in reduced basis models has long posed a difficult problem for analysts. Recently, though, several methods have been shown to be effective in obtaining accurate,reduced basis models for rotor blades. These methods are the axial elongation method,the mixed finite element method, and the nonlinear normal mode method. The main objective of this paper is to demonstrate the close relationships among these methods, which are seemingly disparate at first glance. First, the difficulties inherent in obtaining reduced basis models of rotor blades are illustrated by examining the modal reduction accuracy of several blade analysis formulations. It is shown that classical, displacement-based finite elements are ill-suited for rotor blade analysis because they can't accurately represent the axial strain in modal space, and that this problem may be solved by employing the axial force as a variable in the analysis. It is shown that the mixed finite element method is a convenient means for accomplishing this, and the derivation of a mixed finite element for rotor blade analysis is outlined. A shortcoming of the mixed finite element method is that is that it increases the number of variables in the analysis. It is demonstrated that this problem may be rectified by solving for the axial displacements in terms of the axial forces and the bending displacements. Effectively, this procedure constitutes a generalization of the widely used axial elongation method to blades of arbitrary topology. The procedure is developed first for a single element, and then extended to an arbitrary assemblage of elements of arbitrary type. Finally, it is shown that the generalized axial elongation method is essentially an approximate solution for an invariant manifold that can be used as the basis for a nonlinear normal mode.

Ruzicka, Gene C.↗
NASA NTRS2001

Development of an Innovative Algorithm for Aerodynamics-Structure Interaction Using Lattice Boltzmann Method

The lattice Boltzmann equation (LBE) is a kinetic formulation which offers an alternative computational method capable of solving fluid dynamics for various systems. Major advantages of the method are owing to the fact that the solution for the particle distribution functions is explicit, easy to implement, and the algorithm is natural to parallelize. In this final report, we summarize the works accomplished in the past three years. Since most works have been published, the technical details can be found in the literature. Brief summary will be provided in this report. In this project, a second-order accurate treatment of boundary condition in the LBE method is developed for a curved boundary and tested successfully in various 2-D and 3-D configurations. To evaluate the aerodynamic force on a body in the context of LBE method, several force evaluation schemes have been investigated. A simple momentum exchange method is shown to give reliable and accurate values for the force on a body in both 2-D and 3-D cases. Various 3-D LBE models have been assessed in terms of efficiency, accuracy, and robustness. In general, accurate 3-D results can be obtained using LBE methods. The 3-D 19-bit model is found to be the best one among the 15-bit, 19-bit, and 27-bit LBE models. To achieve desired grid resolution and to accommodate the far field boundary conditions in aerodynamics computations, a multi-block LBE method is developed by dividing the flow field into various blocks each having constant lattice spacing. Substantial contribution to the LBE method is also made through the development of a new, generalized lattice Boltzmann equation constructed in the moment space in order to improve the computational stability, detailed theoretical analysis on the stability, dispersion, and dissipation characteristics of the LBE method, and computational studies of high Reynolds number flows with singular gradients. Finally, a finite difference-based lattice Boltzmann method is developed for inviscid compressible flows.

Mei, Ren-Wei↗
NASA NTRS2002

The Fractional Step Method Applied to Simulations of Natural Convective Flows

This paper describes research done to apply the Fractional Step Method to finite-element simulations of natural convective flows in pure liquids, permeable media, and in a directionally solidified metal alloy casting. The Fractional Step Method has been applied commonly to high Reynold's number flow simulations, but is less common for low Reynold's number flows, such as natural convection in liquids and in permeable media. The Fractional Step Method offers increased speed and reduced memory requirements by allowing non-coupled solution of the pressure and the velocity components. The Fractional Step Method has particular benefits for predicting flows in a directionally solidified alloy, since other methods presently employed are not very efficient. Previously, the most suitable method for predicting flows in a directionally solidified binary alloy was the penalty method. The penalty method requires direct matrix solvers, due to the penalty term. The Fractional Step Method allows iterative solution of the finite element stiffness matrices, thereby allowing more efficient solution of the matrices. The Fractional Step Method also lends itself to parallel processing, since the velocity component stiffness matrices can be built and solved independently of each other. The finite-element simulations of a directionally solidified casting are used to predict macrosegregation in directionally solidified castings. In particular, the finite-element simulations predict the existence of 'channels' within the processing mushy zone and subsequently 'freckles' within the fully processed solid, which are known to result from macrosegregation, or what is often referred to as thermo-solutal convection. These freckles cause material property non-uniformities in directionally solidified castings; therefore many of these castings are scrapped. The phenomenon of natural convection in an alloy under-going directional solidification, or thermo-solutal convection, will be explained. The development of the momentum and continuity equations for natural convection in a fluid, a permeable medium, and in a binary alloy undergoing directional solidification will be presented. Finally, results for natural convection in a pure liquid, natural convection in a medium with a constant permeability, and for directional solidification will be presented.

Westra, Douglas G.↗
NASA NTRS2003

Finite Element Method for Capturing Ultra-relativistic Shocks

While finite element methods are used extensively by researchers solving computational fluid dynamics in fields other than astrophysics, their use in astrophysical fluid simulations has been predominantly overlooked. Current simulations using other methods such as finite difference and finite volume (based on finite difference) have shown remarkable results, but these methods are limited by their fundamental properties in aspects that are important for simulations with complex geometries and widely varying spatial and temporal scale differences. We have explored the use of finite element methods for astrophysical fluids in order to establish the validity of using such methods in astrophysical environments. We present our numerical technique applied to solving ultra-relativistic (Lorentz Factor Gamma >> 1) shocks which are prevalent in astrophysical studies including relativistic jets and gamma-ray burst studies. We show our finite element formulation applied to simulations where the Lorentz factor ranges up to 2236 and demonstrate its stability in solving ultra-relativistic flows. Our numerical method is based on the Flowfield Dependent Variation (FDV) Method, unique in that numerical diffusion is derived from physical parameters rather than traditional artificial viscosity methods. Numerical instabilities account for most of the difficulties when capturing shocks in this regime. Our method results in stable solutions and accurate results as compared with other methods.

Richardson, G. A.↗
NASA NTRS2004

Discontinuous Spectral Difference Method for Conservation Laws on Unstructured Grids

A new, high-order, conservative, and efficient discontinuous spectral finite difference (SD) method for conservation laws on unstructured grids is developed. The concept of discontinuous and high-order local representations to achieve conservation and high accuracy is utilized in a manner similar to the Discontinuous Galerkin (DG) and the Spectral Volume (SV) methods, but while these methods are based on the integrated forms of the equations, the new method is based on the differential form to attain a simpler formulation and higher efficiency. Conventional unstructured finite-difference and finite-volume methods require data reconstruction based on the least-squares formulation using neighboring point or cell data. Since each unknown employs a different stencil, one must repeat the least-squares inversion for every point or cell at each time step, or to store the inversion coefficients. In a high-order, three-dimensional computation, the former would involve impractically large CPU time, while for the latter the memory requirement becomes prohibitive. In addition, the finite-difference method does not satisfy the integral conservation in general. By contrast, the DG and SV methods employ a local, universal reconstruction of a given order of accuracy in each cell in terms of internally defined conservative unknowns. Since the solution is discontinuous across cell boundaries, a Riemann solver is necessary to evaluate boundary flux terms and maintain conservation. In the DG method, a Galerkin finite-element method is employed to update the nodal unknowns within each cell. This requires the inversion of a mass matrix, and the use of quadratures of twice the order of accuracy of the reconstruction to evaluate the surface integrals and additional volume integrals for nonlinear flux functions. In the SV method, the integral conservation law is used to update volume averages over subcells defined by a geometrically similar partition of each grid cell. As the order of accuracy increases, the partitioning for 3D requires the introduction of a large number of parameters, whose optimization to achieve convergence becomes increasingly more difficult. Also, the number of interior facets required to subdivide non-planar faces, and the additional increase in the number of quadrature points for each facet, increases the computational cost greatly.

Liu, Yen↗
NASA NTRS2010

The ICARE Method

The ICARE method is a flexible, widely applicable method for systems engineers to solve problems and resolve issues in a complete and comprehensive manner. The method can be tailored by diverse users for direct application to their function (e.g. system integrators, design engineers, technical discipline leads, analysts, etc.). The clever acronym, ICARE, instills the attitude of accountability, safety, technical rigor and engagement in the problem resolution: Identify, Communicate, Assess, Report, Execute (ICARE). This method was developed through observation of Space Shuttle Propulsion Systems Engineering and Integration (PSE&I) office personnel approach in an attempt to succinctly describe the actions of an effective systems engineer. Additionally it evolved from an effort to make a broadly-defined checklist for a PSE&I worker to perform their responsibilities in an iterative and recursive manner. The National Aeronautics and Space Administration (NASA) Systems Engineering Handbook states, engineering of NASA systems requires a systematic and disciplined set of processes that are applied recursively and iteratively for the design, development, operation, maintenance, and closeout of systems throughout the life cycle of the programs and projects. ICARE is a method that can be applied within the boundaries and requirements of NASA s systems engineering set of processes to provide an elevated sense of duty and responsibility to crew and vehicle safety. The importance of a disciplined set of processes and a safety-conscious mindset increases with the complexity of the system. Moreover, the larger the system and the larger the workforce, the more important it is to encourage the usage of the ICARE method as widely as possible. According to the NASA Systems Engineering Handbook, elements of a system can include people, hardware, software, facilities, policies and documents; all things required to produce system-level results, qualities, properties, characteristics, functions, behavior and performance. The ICARE method can be used to improve all elements of a system and, consequently, the system-level functional, physical and operational performance. Even though ICARE was specifically designed for a systems engineer, any person whose job is to examine another person, product, or process can use the ICARE method to improve effectiveness, implementation, usefulness, value, capability, efficiency, integration, design, and/or marketability. This paper provides the details of the ICARE method, emphasizing the method s application to systems engineering. In addition, a sample of other, non-systems engineering applications are briefly discussed to demonstrate how ICARE can be tailored to a variety of diverse jobs (from project management to parenting).

Henke, Luke↗
NASA NTRS2004

Improved Method of Locating Defects in Wiring Insulation

An improved method of locating small breaches in insulation on electrical wires combines aspects of the prior dielectric withstand voltage (DWV) and time-domain reflectometry (TDR) methods. The method was invented to satisfy a need for reliably and quickly locating insulation defects in spacecraft, aircraft, ships, and other complex systems that contain large amounts of wiring, much of it enclosed in structures that make it difficult to inspect. In the DWV method, one applies a predetermined potential (usually 1.5 kV DC) to the wiring and notes whether the voltage causes any arcing between the wiring and ground. The DWV method does not provide an indication of the location of the defect (unless, in an exceptional case, the arc happens to be visible). In addition, if there is no electrically conductive component at ground potential within about 0.010 in. (approximately equal to 0.254 mm) of the wire at the location of an insulation defect, then the DWV method does not provide an indication of the defect. Moreover, one does not have the option to raise the potential in an effort to increase the detectability of such a defect because doing so can harm previously undamaged insulation. In the TDR method as practiced heretofore, one applies a pulse of electricity having an amplitude of less than 25 V to a wire and measures the round-trip travel time for the reflection of the pulse from a defect. The distance along the wire from the point of application of the pulse to the defect is then calculated as the product of half the round-trip travel time and the characteristic speed of a propagation of an electromagnetic signal in the wire. While the TDR method as practiced heretofore can be used to locate a short or open circuit, it does not ordinarily enable one to locate a small breach in insulation because the pulse voltage is too low to cause arcing and thus too low to induce an impedance discontinuity large enough to generate a measurable reflection. The present improved method overcomes the weaknesses of both the prior DWV and the prior TDR method.

Greulich, Owen R.↗
NASA NTRS2016

Diagonally Implicit Runge-Kutta Methods for Ordinary Differential Equations. A Review

A review of diagonally implicit Runge-Kutta (DIRK) methods applied to rst-order ordinary di erential equations (ODEs) is undertaken. The goal of this review is to summarize the characteristics, assess the potential, and then design several nearly optimal, general purpose, DIRK-type methods. Over 20 important aspects of DIRKtype methods are reviewed. A design study is then conducted on DIRK-type methods having from two to seven implicit stages. From this, 15 schemes are selected for general purpose application. Testing of the 15 chosen methods is done on three singular perturbation problems. Based on the review of method characteristics, these methods focus on having a stage order of two, sti accuracy, L-stability, high quality embedded and dense-output methods, small magnitudes of the algebraic stability matrix eigenvalues, small values of aii, and small or vanishing values of the internal stability function for large eigenvalues of the Jacobian. Among the 15 new methods, ESDIRK4(3)6L[2]SA is recommended as a good default method for solving sti problems at moderate error tolerances.

Kennedy, Christopher A.↗
NASA NTRSDate not supplied

A METHOD TO REDUCE BIOBURDEN IN ASTROMATERIALS CURATION FACILITIES WITHOUT INTRODUCING UNWANTED CONTAMINATION

Introduction: NASA curates its Astromaterials collections in cleanrooms that are carefully monitored for particulate, inorganic and trace metal contamination. Current sample collections are not particularly susceptible to organic contamination or biological alteration. However, new collections like those from the OSIRIS-REx and Hayabusa2 missions will have organic contamination requirements and are susceptible to biodegradation. It will be necessary sterilize or at least disinfect curation labs, as well as tools and equipment in a manner that does not introduce additional contamination and does not affect the samples 1. Current curation cleaning procedures utilize isopropyl alcohol which offers some bioburden reduction, but is not effective against spore-forming bacteria or fungal spores 2. We present a modified disinfection method that uses ultrapure hydrogen peroxide to reduce bioburden inside curation labs and glove boxes without introducing contamination or damaging curation equipment. We tested this method in the meteorite processing lab as well as on a glovebox being cleaned for use in processing ANGSA (Apollo Next Generation Sample Analysis) samples and present the results of those tests. We discuss the limitations of this method and describe potential situations in which it will not be applicable. The CDC guidelines for disinfection andsterilization in healthcare facilities discusses over 15different methods for reducing bioburden in hospitalsettings 3. The most common method, steamsterilization, is well suited to sterilizing curationprocessing tools but cannot easily be used to sterilizecleanroom surfaces or large equipment likegloveboxes. Chemical sterilization with bleach(NaOCl) is also a common strategy in healthcare andpharmaceutical settings that presents materialcompatibility issues as well as serious inorganiccontamination concerns for curation facilities.Introducing a new source of Na and Cl into curationlabs is not acceptable. Other chemical methods likeethylene oxide, formaldehyde, iodophors andquaternary ammonium compounds could introduceorganic and inorganic contamination. We chose tofocus on hydrogen peroxide because it is generallycompatible with commonly used curation materialslike stainless steel, aluminum and Teflon and becauseit decomposes to oxygen and water. The CDCguidelines for hydrogen peroxide specify using a 7.5wt% solution at 25 ̊C with a contact time of 30 minutesfor high level disinfection and 6 hours for sterilization.High level disinfection is defined as a technique thatwill kill all microorganisms except large numbers ofbacterial spores 3. Methods: We prepared a solution of 7.5 wt%hydrogen peroxide from a stock solution of ultrapure30 wt% peroxide (JT Baker) and curation gradeultrapure water. This ultrapure water is already used incuration cleaning procedures and thus is not consideredand additional source of contamination. We conducteda materials compatibility test by exposing unanodizedand anodized 6061 T6 Al alloy to the peroxide solutionfor up to six hours and periodically inspecting thesurfaces for visible defects. We used this peroxide todisinfect the floor of the meteorite processing lab andthe interior of a curation glovebox by exposing thesesurfaces to the peroxide solution for 30 min. Thesurfaces were swabbed with a dry macrofoam swabbefore (Puritan Brand 2518051PFRNDFD) and afterperoxide treatment to collect microbes present on thesurfaces. Microbes were extracted by sonication fromthe swab into 15 ml of PBS (phosphate buffered saline)and inoculated onto the following media: TSA (trypticsoy agar) BA (blood agar), R2A (Reasoners 2 agar),Potato Dextrose Agar, Saboraud Dextrose Agar andSaboraud Dextrose Agar with 0.1 mg/ mlchloramphenicol. Four TSA plates and two BA plateswere inoculated with 0.1 ml of PBS each andincubated at 35 and 37 for 48 hours. Two R2A°C°Cplates (0.1 ml of PBS each) were incubated at 25 .°CThe remaining plates were inoculated with 0.2ml ofPBS and incubated at 30 ̊C for seven days. Afterincubation bacterial and fungal isolates were countedand transferred to new plates for identification usingthe VITEK24 automated system or by sequencing aportion of the barcode gene (16S rRNA for bacteria,small subunit gene for fungi) on an ABI 3500 Sangersequencer. Negative controls consisted of swabs thatwere opened in the sampling environment andanalyzed alongside the experimental samples.Results: A 6 hour exposure to hydrogen peroxideresulted in visible pitting on un-anodized 6061 Al, butnot on anodized surfaces. No visible pitting occurredafter a 30 minute exposure. Therefore, we decided tolimit our experimental tests to 30 min. exposures. 17bacterial CFU (colony forming units) representing 4distinct organisms were isolated from the meteorite processing lab floor prior to hydrogen peroxidetreatment. We were unable culture any organisms afterperoxide treatment. In the glovebox we were able toculture three bacterial CFU representing three distinctspecies, including a spore forming bacterium prior todisinfection with peroxide. After the peroxidetreatment we were unable to culture any organisms.Routine monitoring of the meteorite processing lab andthe glovebox did not indicate any increase in unwantedinorganic contamination after these peroxidetreatments. Discussion: A 30 minute treatment with 7.5 wt%peroxide appears to be an effective method forreducing bioburden on typical cleanroom surfaces. Themethod does not introduce unwanted organic orinorganic contamination and is compatible withcommonly used curation materials like stainless steel,Teflon and anodized aluminum alloys. Special careshould be taken with un-anodized aluminum.Prolonged exposure to hydrogen peroxide can causepitting on this material. We recommend using thismethod to disinfect curation labs and equipment whenbiological alteration is a concern. This method iseffective at room temperature and cannot be used todisinfect labs and equipment where the ambienttemperature is < 0 ̊C. Astromaterials samples shouldbe removed from the area where disinfection is tooccur. Hydrogen peroxide is a powerful oxidizingagent and will react with any organic carbon present inthe sample. References: [1.] Mccubbin, F. M. et al.Sp. Sci Rev(2019) doi:10.1007/s11214-019-0615-9. [2.] Mogul, R.et al.Astrobiology 18, ast.2017.1814 (2018). [3.]Rutala, W. A. & Weber, D. J. Guideline for Disinfection and Sterilization in Healthcare Facilities, 2008. [4.] Pincus, D. H. in Encyclopedia of Rapid Microbiological Methods (2005).

A. B. Regberg↗
NASA NTRSDate not supplied

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↗