Search NASA⌕ Search

SEARCH · Search NASA

Results for “Methods”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
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↗
NASA NTRSDate not supplied

Comparison of Two Load Prediction Methods for Strain-Gage Balances

Data from a five-component semi-span balance is used to perform a systematic comparison of the load prediction accuracy of two load prediction methods. Both methods independently obtain the load prediction equations from multivariate least squares fits of balance calibration data. The first method is called the Non-Iterative Method. This approach directly uses regression models of the individual load components of a balance for the load prediction. The second method is called the Iterative Method. This alternate approach uses a load iteration equation for the load prediction that is constructed from the regression coefficients of the gage outputs of the balance. Basic characteristics of the two methods are reviewed. Afterwards, both methods are applied to calibration, check load, and wind tunnel test data of a five-component semi-span balance. Selected analysis results are compared. These comparisons confirm that the accuracy of the two methods is the same for all practical purposes.

strain-gage balance↗
NASA NTRS1972

A comparative study of stability methods for flexible satellites.

This paper compares three approaches to the stability of hybrid dynamical systems, all three methods being based on the Liapunov direct method. The first method uses testing density functions, whereas the second involves defining certain integral coordinates. Both the method using testing density functions and the method of integral coordinates lead to closed-form stability criteria in terms of the system parameters. Criteria obtained using the method of integral coordinates are in general less restrictive than those derived by the method using testing density functions. On the other hand, the latter method is easier to apply and requires less work than the former. The third method is the standard modal analysis. The modal analysis generally yields more involved criteria, depending on the number of modes used to represent the elastic displacements. As an application, the attitude stability of an earth-pointing satellite with multi-elastic domains is investigated.

Meirovitch, L.↗
NASA NTRS1973

Analysis and optimization of cyclic methods in orbit computation

The mathematical analysis and computation of the K=3, order 4; K=4, order 6; and K=5, order 7 cyclic methods and the K=5, order 6 Cowell method and some results of optimizing the 3 backpoint cyclic multistep methods for solving ordinary differential equations are presented. Cyclic methods have the advantage over traditional methods of having higher order for a given number of backpoints while at the same time having more free parameters. After considering several error sources the primary source for the cyclic methods has been isolated. The free parameters for three backpoint methods were used to minimize the effects of some of these error sources. They now yield more accuracy with the same computing time as Cowell's method on selected problems. This work is being extended to the five backpoint methods. The analysis and optimization are more difficult here since the matrices are larger and the dimension of the optimizing space is larger. Indications are that the primary error source can be reduced. This will still leave several parameters free to minimize other sources.

Pierce, S.↗
NASA NTRS1974

Comparison of several methods for predicting separation in a compressible turbulent boundary layer

Several methods for predicting the separation point for a compressible turbulent boundary layer were applied to the flow over a bump on a wind-tunnel wall. Measured pressure distributions were used as input. Two integral boundary-layer methods, three finite-difference boundary-layer methods, and three simple methods were applied at five free-stream Mach numbers ranging from 0.354 to 0.7325. Each of the boundary-layer methods failed to explicitly predict separation. However, by relaxing the theoretical separation criteria, several boundary-layer methods were made to yield reasonable separation predictions, but none of the methods accurately predicted the important boundary-layer parameters at separation. Only one of the simple methods consistently predicted separation with reasonable accuracy in a manner consistent with the theory. The other methods either indicated several possible separation locations or only sometimes predicted separation.

Gerhart, P. M.↗
NASA NTRS1978

Application of higher-order numerical methods to the boundary-layer equations

A fourth-order method is presented for calculating numerical solutions to parabolic, partial differential equations in two variables or ordinary differential equations. The method is the natural extension of the second-order Keller Box Scheme to fourth order and is demonstrated with application to the incompressible, laminar and turbulent boundary-layer equations for both attached and separated flows. The efficiency of the present method is compared with other higher-order methods; namely, the Keller Box Scheme with Richardson extrapolation, the method of deferred corrections, the three-point spline methods, and a modified finite-element method. For equivalent accuracy, numerical results show the present method to be more efficient than the other higher-order methods for both laminar and turbulent flows.

Wornom, S. F.↗
NASA NTRS1978

Critical study of higher order numerical methods for solving the boundary-layer equations

A fourth order box method is presented for calculating numerical solutions to parabolic, partial differential equations in two variables or ordinary differential equations. The method, which is the natural extension of the second order box scheme to fourth order, was demonstrated with application to the incompressible, laminar and turbulent, boundary layer equations. The efficiency of the present method is compared with two point and three point higher order methods, namely, the Keller box scheme with Richardson extrapolation, the method of deferred corrections, a three point spline method, and a modified finite element method. For equivalent accuracy, numerical results show the present method to be more efficient than higher order methods for both laminar and turbulent flows.

Wornom, S. F.↗
NASA NTRS1981

Methods for the calculation of axial wave numbers in lined ducts with mean flow

A survey is made of the methods available for the calculation of axial wave numbers in lined ducts. Rectangular and circular ducts with both uniform and non-uniform flow are considered as are ducts with peripherally varying liners. A historical perspective is provided by a discussion of the classical methods for computing attenuation when no mean flow is present. When flow is present these techniques become either impractical or impossible. A number of direct eigenvalue determination schemes which have been used when flow is present are discussed. Methods described are extensions of the classical no-flow technique, perturbation methods based on the no-flow technique, direct integration methods for solution of the eigenvalue equation, an integration-iteration method based on the governing differential equation for acoustic transmission, Galerkin methods, finite difference methods, and finite element methods.

Eversman, W.↗