Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Exact two component”

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 127 records · Page 7

Accurate interlaminar stress recovery from finite element analysis

The accuracy and robustness of a two-dimensional smoothing methodology is examined for the problem of recovering accurate interlaminar shear stress distributions in laminated composite and sandwich plates. The smoothing methodology is based on a variational formulation which combines discrete least-squares and penalty-constraint functionals in a single variational form. The smoothing analysis utilizes optimal strains computed at discrete locations in a finite element analysis. These discrete strain data are smoothed with a smoothing element discretization, producing superior accuracy strains and their first gradients. The approach enables the resulting smooth strain field to be practically C1-continuous throughout the domain of smoothing, exhibiting superconvergent properties of the smoothed quantity. The continuous strain gradients are also obtained directly from the solution. The recovered strain gradients are subsequently employed in the integration o equilibrium equations to obtain accurate interlaminar shear stresses. The problem is a simply-supported rectangular plate under a doubly sinusoidal load. The problem has an exact analytic solution which serves as a measure of goodness of the recovered interlaminar shear stresses. The method has the versatility of being applicable to the analysis of rather general and complex structures built of distinct components and materials, such as found in aircraft design. For these types of structures, the smoothing is achieved with 'patches', each patch covering the domain in which the smoothed quantity is physically continuous.

Tessler, Alexander↗

Radiative Impacts of Cloud Heterogeneity and Overlap in an Atmospheric General Circulation Model

The radiative impacts of introducing horizontal heterogeneity of layer cloud condensate, and vertical overlap of condensate and cloud fraction are examined with the aid of a new radiation package operating in the GEOS-5 Atmospheric General Circulation Model. The impacts are examined in terms of diagnostic top-of-the-atmosphere shortwave (SW) and longwave (LW) cloud radiative effect (CRE) calculations for a range of assumptions and parameter specifications about the overlap. The investigation is conducted for two distinct cloud schemes, the one that comes with the standard GEOS-5 distribution, and another which has been recently used experimentally for its enhanced GEOS-5 distribution, and another which has been recently used experimentally for its enhanced cloud microphysical capabilities; both are coupled to a cloud generator allowing arbitrary cloud overlap specification. We find that cloud overlap radiative impacts are significantly stronger for the operational cloud scheme for which a change of cloud fraction overlap from maximum-random to generalized results to global changes of SW and LW CRE of approximately 4 Watts per square meter, and zonal changes of up to approximately 10 Watts per square meter. This is because of fewer occurrences compared to the other scheme of large layer cloud fractions and of multi-layer situations with large numbers of atmospheric being simultaneously cloudy, conditions that make overlap details more important. The impact on CRE of the details of condensate distribution overlap is much weaker. Once generalized overlap is adopted, both cloud schemes are only modestly sensitive to the exact values of the overlap parameters. We also find that if one of the CRE components is overestimated and the other underestimated, both cannot be driven towards observed values by adjustments to cloud condensate heterogeneity and overlap alone.

Oreopoulos, L.↗

In Situ Resource Utilization (ISRU) Experiments for Mars Exploration

In situ resource utilization can best be described as living off the land. In our case the land is the planet Mars. ISRU is based on the idea that some fraction of the consumables, life support and propellant materials do not have to be flown from earth. Rather, they can be manufactured or extracted from resources already present on Mars. The primary resources on Mars are the atmosphere, polar caps and regolith. The atmosphere of Mars is mostly carbon dioxide as shown in the table below. The proportion of oxygen on the other hand is quite small. Still, there is quite a bit of oxygen in the Martian atmosphere, but it is unfortunately tied up with carbon. Thus, one of the goals of ISRU is the separation of breathable oxygen from the carbon dioxide. Several means of separation have been proposed. We have begun experiments on another approach for production of oxygen with carbon monoxide as a useful by product. Our work on a CO2 separator is described later in this report. Regolith melting is another means of obtaining materials. Two materials of interest are iron and silicon. Iron oxide is plentiful on Mars and is of obvious importance for structural components. Silicon is the foundation of solid state devices. Power generation on Mars may be accomplished using silicon solar cells. There is discussion of the feasibility of in situ production of solar cells. This would require a means of extracting silicon from the regolith. We have conducted several experiments concerning melting and glassification of the Mars soil simulant. Other summer faculty fellows have tried various means of processing the stimulant material. These include furnace melting, microwave melting and laser ablation. We have conducted several furnace melting experiments in both air and carbon dioxide environments. We have also carried out experiments to test spark melting in a carbon dioxide atmosphere. These experiments suggest the possibility of using arc melting in a reducing atmosphere. It is important to keep in mind that we are working with a soil stimulant. Any simulant, no matter how chemically similar it is to Martian regolith, may differ in mineralogy. The underlying assumption in this work is that once a glass is formed, any differences between simulant and regolith are unimportant. The exact means of forming the glass do, however, depend on the mineralogy of the regolith. A sample return mission is required to help answer these questions.

Marone, Matt↗

A Study of the Free Oscillations of the Earth

Published observations on the toroidal oscillations of the earth are critically reviewed. A supplementary analysis of the record obtained by the Lamont strain seismometer is presented. Eleven toroidal modes are identified, and it is concluded that the periods are known to within 1 per cent. A perturbation scheme involving the ratio of the angular velocity of the earth to the resonant frequency is used in calculating the effects due to the rotation of the earth on the resonant frequency. The free oscillations are viewed as a superposition of traveling waves. In a nonrotating system two traveling waves combine to produce a stationary standing wave. In a rotating system, the rotation distinguishes between waves that travel in the direction of rotation and those that travel in the opposite direction. Rotation removes a degeneracy and results in a splitting of a spectral peak of order l into 2 times l plus 1 peaks. The fractional displacement in frequency for the lowest-order toroidal oscillations is 1/206 and of the same order as the Q of the peak, so that splitting will probably not be observed in the toroidal oscillations. Viewed locally, rotation causes a particle to precess about a direction parallel to the axis of rotation. This precession will cause a variation of amplitude with time if the motion is recorded by an instrument with an anisotropic response function. Care is therefore needed in studying the time decay of a given spectral peak. Rotation also couples the normal coordinates so that a motion that is initially purely horizontal will develop a vertical component. It is expected that vertical seismometers should record particle motion with the toroidal frequencies. The perturbations of the toroidal oscillations due to core-mantle interaction are treated in detail. An exact expression is obtained for the rate of energy dissipated by a finitely conducting plate oscillating across a magnetic field. The energy dissipated at the core-mantle boundary due to viscous and hydromagnetic coupling is shown to be insignificant as compared with the energy dissipated within the mantle. The toroidal magnetic field leaking into the lower mantle combines with the dipole field, resulting in a stress on the mantle, tending to stiffen the lower boundary. The stress is of sufficient magnitude to produce a displacement toward higher frequency in the lower-order toroidal oscillations. Observations on the (sub 0) T (sub 2) oscillations lead to an estimate of the toroidal magnetic field in the lower mantle. A calculation of elastic energy in the low-order oscillations suggests a value of 10 (sup 18) ergs per cycles per hour for the energy density at low frequencies in the Chilean earthquake. Each mode of oscillation has a characteristic radial distribution of elastic energy associated with it. This distribution determines which parts of the earth contribute most heavily in determining a particular resonant frequency. The distribution of energy for the lower 17 modes for a homogeneous and a Gutenberg model earth is calculated. The resonant frequencies for models of the earth based on the Gutenberg and Lehmann distribution of elastic properties are presented. It is shown that the Gutenberg model earth fits the observations more closely than the Lehmann model and that a slight alteration of the Gutenberg model gives a significantly better fit to the observations. The alteration involves a lower shear-wave velocity in the lower mantle while the Gutenberg velocity distribution is maintained in the upper mantle. Various studies of the earth's oscillations coupled with surface-wave investigations substantiate Gutenberg's hypothesis of a layer of low velocity in the upper mantle. The physical conditions required for the formation of a region of low velocity are examined in detail. The results confirm Birch's earlier statement that a temperature gradient in excess of 6 degrees to 7 degrees per kilometer is needed to produce a decrease in velocity. The low-velocity layer does not require that the temperature approach or exceed the melting temperature. If tile upper mantle is homogeneous, the region of lower velocity should commence at the base of the crust and extend to 150 kilometers under the oceans and about 100 kilometers under continental regions. The distribution of thermal conductivity and radioactivity consistent with the low-velocity layer is also considered.

MacDonald, Gordon J. F.↗

Discrete-Layer Piezoelectric Plate and Shell Models for Active Tip-Clearance Control

The objectives of this work were to develop computational tools for the analysis of active-sensory composite structures with added or embedded piezoelectric layers. The targeted application for this class of smart composite laminates and the analytical development is the accomplishment of active tip-clearance control in turbomachinery components. Two distinct theories and analytical models were developed and explored under this contract: (1) a discrete-layer plate theory and corresponding computational models, and (2) a three dimensional general discrete-layer element generated in curvilinear coordinates for modeling laminated composite piezoelectric shells. Both models were developed from the complete electromechanical constitutive relations of piezoelectric materials, and incorporate both displacements and potentials as state variables. This report describes the development and results of these models. The discrete-layer theories imply that the displacement field and electrostatic potential through-the-thickness of the laminate are described over an individual layer rather than as a smeared function over the thickness of the entire plate or shell thickness. This is especially crucial for composites with embedded piezoelectric layers, as the actuating and sensing elements within these layers are poorly represented by effective or smeared properties. Linear Lagrange interpolation polynomials were used to describe the through-thickness laminate behavior. Both analytic and finite element approximations were used in the plane or surface of the structure. In this context, theoretical developments are presented for the discrete-layer plate theory, the discrete-layer shell theory, and the formulation of an exact solution for simply-supported piezoelectric plates. Finally, evaluations and results from a number of separate examples are presented for the static and dynamic analysis of the plate geometry. Comparisons between the different approaches are provided when possible, and initial conclusions regarding the accuracy and limitations of these models are given.

Heyliger, P. R.↗

Magnetohydrodynamic Riemann problem and the structure of the magnetic reconnection layer

We present a complete solution for a set of magnetohydrodynamic (MHD) Riemann problems in which the upstream and downstream states have the same total pressure, and in which the normal component of the magnetic field is very small. These solutions are pertinent to subfast flows in the earth's magnetic tail and near the magnetopause. In a coplanar situation a family of solutions exists that depend on two parameters as well as on dissapation mechanisms. In the parallel case the transverse magnetic field either does not change direction or changes the direction twice by involving two intermediate shocks. In the antiparallel case an intermediate shock is always required, except when the solution consists of two switch-off shocks. In a noncoplanar case the solution is not self-similar as a function of x/t, but continues to evolve. At early times the evolution is similar to the coplanar case. In general two time-dependent intermediate shocks are required to rotate the magnetic fields. The velocity shear has a strong effect on the Riemann solution. In some cases no Riemann solution can exist because of the cavitation caused by the slow refraction waves. The calculated magnetopause structure resembles the observed structure for northward interplanetary magnetic field (IMF). However, for southward IMF, the MHD result shows the existence of a depletion layer, which is not supported by observations. We also show that on the magnetosheath side, the Walen relation, which is exact for a rotational discontinuity, can also be well satisfied by a slow shock, an intermediate shock, or the head of a slow rarefaction wave.

Wu, C. C.↗

Nonlinear structure formation in flat cosmological models

This paper describes the formation of nonlinear structure in flat (zero curvature) Friedmann cosmological models. We consider models with two components: the usual nonrelativistic component that evolves under gravity and eventually forms the large-scale structure of the universe, and a uniform dark matter component that does not clump under gravity, and whose energy density varies with the scale factor a(t) like a(t)(sup -n), where n is a free parameter. Each model is characterized by two parameters: the exponent n and the present density parameter Omega(sub 0) of the nonrelativistic component. The linear perturbation equations are derived and solved for these models, for the three different cases n = 3, n is greater than 3, and n is less than 3. The case n = 3 is relevant to model with massive neutrinos. The presence of the uniform component strongly reduces the growth of the perturbation compared with the Einstein-de Sitter model. We show that the Meszaros effect (suppression of growth at high redshift) holds not only for n = 4, radiation-dominated models, but for all models with n is greater than 3. This essentially rules out any such model. For the case n is less than 3, we numerically integrate the perturbation equations from the big bang to the present, for 620 different models with various values of Omega(sub 0) and n. Using these solutions, we show that the function f(Omega(sub 0), n) = (a/delta(sub +))d(delta)(sub +)/da, which enters in the relationship between the present density contrast delta(sub 0) and peculiar velocity field u(sub 0) is essentially independent of n. We derive approximate solutions for the second-order perturbation equations. These second-order solutions are tested against the exact solutions and the Zel'dovich approximation for spherically symmetric perturbations in the marginally nonlinear regime (the absolute value of delta is less than or approximately 1). The second-order and Zel'dovich solutions have comparable accuracy, significantly higher than the accuracy of the linear solutions. We then investigate the dependence of the delta(sub 0) - u(sub 0) relationship upon the value of n in the nonlinear regime using the second-order solutions for marginally nonlinear, general perturbations, and the exact solutions for strongly nonlinear, spherically symmetric perturbations. In both cases, we find that the delta(sub 0) - u(sub 0) relationship remains independent of n. We speculate that this result extends to strongly nonlinear, general perturbations as well. This eliminates any hope to determine the presence of the uniform component or the value of n using dynamical methods. Finally, we compute the nonlinear evolution of the skewness of the distribution of values of delta, assuming Gaussian initial conditions. We find that the skewness is not only independent of n, but also of Omega(sub 0). Thus the skewness cannot be used to discriminate among various models with Gaussian initial conditions. However, it can be used for testing the Gaussianity of the initial conditions themselves. We conclude that the uniform component leaves no observable signature in the present large-scale structure of the universe. To determine its presence and nature, we must investigate the relationship between the past and present universe, using redshift-dependent tests.

Martel, Hugo↗

Implementation of Speed Variation in the Structural Dynamic Assessment of Turbomachinery Flow-Path Components

During the design of turbomachinery flow path components, the assessment of possible structural resonant conditions is critical. Higher frequency modes of these structures are frequently found to be subject to resonance, and in these cases, design criteria require a forced response analysis of the structure with the assumption that the excitation speed exactly equals the resonant frequency. The design becomes problematic if the response analysis shows a violation of the HCF criteria. One possible solution is to perform "finite-life" analysis, where Miner's rule is used to calculate the actual life in seconds in comparison to the required life. In this situation, it is beneficial to incorporate the fact that, for a variety of turbomachinery control reasons, the speed of the rotor does not actually dwell at a single value but instead dithers about a nominal mean speed and during the time that the excitation frequency is not equal to the resonant frequency, the damage accumulated by the structure is diminished significantly. Building on previous investigations into this process, we show that a steady-state assumption of the response is extremely accurate for this typical case, resulting in the ability to quickly account for speed variation in the finite-life analysis of a component which has previously had its peak dynamic stress at resonance calculated. A technique using Monte Carlo simulation is also presented which can be used when specific speed time histories are not available. The implementation of these techniques can prove critical for successful turbopump design, as the improvement in life when speed variation is considered is shown to be greater than a factor of two.

Brown, Andrew M.↗

Implementation of Speed Variation in the Structural Dynamic Assessment of Turbomachinery Flow-Path Components

During the design of turbomachinery flow path components, the assessment of possible structural resonant conditions is critical. Higher frequency modes of these structures are frequently found to be subject to resonance, and in these cases, design criteria require a forced response analysis of the structure with the assumption that the excitation speed exactly equals the resonant frequency. The design becomes problematic if the response analysis shows a violation of the HCF criteria. One possible solution is to perform "finite-life" analysis, where Miner's rule is used to calculate the actual life in seconds in comparison to the required life. In this situation, it is beneficial to incorporate the fact that, for a variety of turbomachinery control reasons, the speed of the rotor does not actually dwell at a single value but instead dithers about a nominal mean speed and during the time that the excitation frequency is not equal to the resonant frequency, the damage accumulated by the structure is diminished significantly. Building on previous investigations into this process, we show that a steady-state assumption of the response is extremely accurate for this typical case, resulting in the ability to quickly account for speed variation in the finite-life analysis of a component which has previously had its peak dynamic stress at resonance calculated. A technique using Monte Carlo simulation is also presented which can be used when specific speed time histories are not available. The implementation of these techniques can prove critical for successful turbopump design, as the improvement in life when speed variation is considered is shown to be greater than a factor of two

Brown, Andrew M.↗

Blind Study Validating Parametric Costing Tools PRICE TruePlanning and SEER-H for NASA Science Missions

Two of the primary parametric costing tools used to estimate the development and production cost of future spacecraft hardware are PRICE TruePlanning - Space Missions Catalog by PRICE Systems, and SEER-H by Galorath. These are standard tools used by NASA and industry to estimate the cost of new aerospace hardware. However, no independent verification of the accuracy of these tools is publicly available. Both PRICE Systems and Galorath have completed internal validation studies of their parametric cost estimating tools; however, they only provided the results of the studies and did not detail the exact methods used to perform the validation. In the present study, cost estimators used PRICE TruePlanning and SEER-H to estimate the cost of twelve different past NASA science missions. The estimators were prevented from knowing the actual cost of the missions in an effort to minimize cognitive biases. In the present study, SEER had an average error of 23%, median error of -0.3%, with a standard deviation of 43%. PRICE had an average error of 52%, median error of 50%, and standard deviation of 45%. There were several factors independent of PRICE and SEER which may have affected the accuracy of the results in the present study including: uncertainty in the technical data used for the estimates, the methods used to estimate uncertainty in spacecraft component mass and numbers of prototypes, and the experience of the estimators.

Paul D Friz↗

Blind Study Validating Parametric Costing Tools PRICE True Planning and SEER-H for NASA Science Missions

Two of the primary parametric costing tools used to estimate the development and production cost of future spacecraft hardware are PRICE True Planning - Space Missions Catalog by PRICE Systems, and SEER-H by Galorath. These are standard tools used by NASA and industry to estimate the cost of new aerospace hardware. However, no independent verification of the accuracy of these tools is publicly available. Both PRICE Systems and Galorath have completed internal validation studies of their parametric cost estimating tools; however, they only provided the results of the studies and did not detail the exact methods used to perform the validation. In the present study, cost estimators used PRICE True Planning and SEER-H to estimate the cost of twelve different past NASA science missions. The estimators were prevented from knowing the actual cost of the missions in an effort to minimize cognitive biases. In the present study, SEER had an average error of 23%, median error of -0.3%, with a standard deviation of 43%. PRICE had an average error of 52%, median error of 50%, and standard deviation of 45%. There were several factors independent of PRICE and SEER which may have affected the accuracy of the results in the present study including: uncertainty in the technical data used for the estimates, the methods used to estimate uncertainty in spacecraft component mass and numbers of prototypes, and the experience of the estimators.

Friz, Paul D.↗

Parallel 3D Multi-Stage Simulation of a Turbofan Engine

A 3D multistage simulation of each component of a modern GE Turbofan engine has been made. An axisymmetric view of this engine is presented in the document. This includes a fan, booster rig, high pressure compressor rig, high pressure turbine rig and a low pressure turbine rig. In the near future, all components will be run in a single calculation for a solution of 49 blade rows. The simulation exploits the use of parallel computations by using two levels of parallelism. Each blade row is run in parallel and each blade row grid is decomposed into several domains and run in parallel. 20 processors are used for the 4 blade row analysis. The average passage approach developed by John Adamczyk at NASA Lewis Research Center has been further developed and parallelized. This is APNASA Version A. It is a Navier-Stokes solver using a 4-stage explicit Runge-Kutta time marching scheme with variable time steps and residual smoothing for convergence acceleration. It has an implicit K-E turbulence model which uses an ADI solver to factor the matrix. Between 50 and 100 explicit time steps are solved before a blade row body force is calculated and exchanged with the other blade rows. This outer iteration has been coined a "flip." Efforts have been made to make the solver linearly scaleable with the number of blade rows. Enough flips are run (between 50 and 200) so the solution in the entire machine is not changing. The K-E equations are generally solved every other explicit time step. One of the key requirements in the development of the parallel code was to make the parallel solution exactly (bit for bit) match the serial solution. This has helped isolate many small parallel bugs and guarantee the parallelization was done correctly. The domain decomposition is done only in the axial direction since the number of points axially is much larger than the other two directions. This code uses MPI for message passing. The parallel speed up of the solver portion (no 1/0 or body force calculation) for a grid which has 227 points axially.

Turner, Mark G.↗

Comparing Theoretically Scaled Biomechanical Models

BACKGROUND An investigation into incorporating space suit aspects required a 50th-percentile male model which was not part of the existing dataset of models. Previously, theoretical models for 5th-percentile female and 95th-percentile male were created utilizing a scaled test subject as close as possible to the target height and weight. Scaling factors were created by taking the height ratio and applying it uniformly to the model body segments, then fine tuning to ensure the theoretical model’s height is as expected. This study was initiated to examine the existing method of isometrically scaling in OpenSim [1,2] and to create alternative methods which do not rely on an existing subject being close in height and weight to the theoretical model of interest. Generating a theoretical biomechanical model provides additional abilities without the reliance on available OpenSim models or real subjects. Gained abilities include creating different percentile models and attaining representative anthropometry for any target height/weight, such as targeting specific crew populations or gaps within the current dataset. METHODS AND RESULTS This study includes 3 methods for generating the theoretical model; utilizing the modified unscaled OpenSim Full Body Rajagopal Model (FBRM) [3,4], a dataset of scaled OpenSim models, and existing scaling factors from Dumas et al. [5]. All methods use the Anthropometric Survey of US Army Personnel (ANSUR II) [6] as an input of necessary anthropometric measurements. The following measurements are retrieved from the ANSUR II collection: mass, stature, cervicale height, acromial height, axilla height, waist height, trochanterion height, lateral femoral epicondyle height, lateral malleolus height, acromion-radiale length, radiale-stylion length, palm length, ball of foot length, bicristal breadth, bimalleolar breadth. Some of the measurements are direct segment lengths and others are utilized to derive segment lengths. Height, weight, and age ranges are the required inputs to parse the collection for a mean value of the measurements. A simple iterative process may be necessary for the output of mean mass and stature consistent with the theoretical model of interest. In some cases, one may specify exact values instead of a range but that is dependent upon whether those exact values pertain to a single subject from the ANSUR II collection. Equations from Dumas et al. were used to calculate scaling factors for the two methods that utilize OpenSim scaled and unscaled models. In the third method, Dumas’ scaling factors were used directly instead of generating our own. After scaling factors are achieved, they are applied along with the mass and segment lengths retrieved from ANSUR II in order to calculate the segment mass, Center of Mass (CoM), mass Moment of Inertia (MoI), and the joint location in the parent frame. The theoretical models from the different methods were compared to the previous existing method of isometrically scaling in OpenSim, as well as the standards found in NASA STD-3000 [7] and the NASA Human Integration Design Handbook [8]. The whole-body center of mass is the main comparison performed between the models and the standards. The body segment mass properties were also compared when possible. In some cases, the standards do not provide complete data, or the center of mass location is not provided with respect to the appropriate joint coordinate system. For a 50th-percentile male subject, results from isometric scaling in OpenSim and the method utilizing the unscaled FBRM were compared to the standards in NASA STD-3000. According to the standard, the whole-body vertical CoM location was taken from the head vertex to the CoM. The standards provide 80.2 cm for the CoM vertical component whereas isometric scaling in OpenSim and utilizing the unscaled FBRM provide 82.2 and 82.3 cm, respectively. Through definition, the whole-body CoM in the medial/lateral direction agreed across the models. Investigation is on-going to clearly identify the reference point pertaining to dorsal/ventral whole-body CoM position in standards and relate that to the theoretical models. REFERENCES [1] Delp, S.L., et al., “OpenSim: Open-source Software to Create and Analyze Dynamic Simulations of Movement”, IEEE Transactions on Biomedical Engineering, (2007). [2] Seth, A., Hicks J.L., Uchida, T.K., Habib, A., Dembia, C.L., Dunne, J.J., Ong, C.F., DeMers, M.S., Rajagopal, A., Millard, M., Hamner, S.R., Arnold, E.M., Yong, J.R., Lakshmikanth, S.K., Sherman,n M.A., Delp, S.L. OpenSim: Simulating musculoskeletal dynamics and neuromuscular control to study human and animal movement. Plos Computational Biology, 14(7). (2018) [3] R. K. Huffman, W. K. Thompson, C. A. Gallo, L. J. Quiocho, “Improvement of Scaling and Inverse Kinematic Results with Additional Upper Body Joints Added to Opensim Rajagopal Model”, NASA Human Research Program Investigator’s Workshop, (2019). [4] Huffman R.K., Thompson W., Gallo C., “Modified OpenSim Rajagopal Full Body Model”, New Technology Report, MSC-26872-1. (August 2020). [5] Dumas R., Cheze L. and Verriest J., 2007. “Adjustments to McConville et al. and Young et al. body segment inertial parameters”. Journal of Biomechanics 40, pp. 543-553 [6] Gordon C.C., 2012 “Anthropometric Survey of U.S. Army Personnel: Methods and Summary Statistics”, US Army Natick Soldier RD&E Center. [7] NASA STD-3000/REV-B, 1995. “The Man-System Integration Standards” [8] NASA/SP-2010-3407/REV-1, 2014. “The Human Integration Design Handbook (HIDH)”

F N Matari↗

Use of Otoacousticemission Phase Change to Evaluate Countermeasures for Spaceflight-Associated Neuro-Ocular Syndrome

Spaceflight-associated neuro-ocular syndrome (SANS) is a human spaceflight risk recognized by NASA. Elevated intracranial pressure (ICP) has been implicated as a root cause of many SANS signs and symptoms, yet there is no reliable noninvasive means of monitoring ICP. We have developed a noninvasive method of monitoring ICP change that exploits ear canal acoustic and otoacoustic emission (OAE) measurements. Changed ICP alters pressure in the inner ear, leading to changes in the tension and position of middle ear (ME) components; tension of these components determines the phase of the stimulus in the ear canal and the OAE response sound transmission back through the ME. The OAE method has been validated in several studies, including our own experiments as part of the NASA Fluid Shifts study. Systematic OAE phase changes demonstrating increased ME tension (an ICP indicator) are observed as posture is changed from seated to supine to head-down tilt (HDT). This effect can be substantially mitigated by lower body negative pressure (LBNP). The OAE technique has also been used on International Space Station (ISS) crewmembers, providing evidence that ICP in microgravity is similar to that seen on the ground in the supine position. The OAE method is also a rapid and noninvasive means of assessing the effectiveness of SANS countermeasures. Here we report results from two studies which used OAEs. In the most recent study (Venous Congestion Countermeasures - VCCM), three promising countermeasures [LBNP, an impedance threshold device (ITD), and veno-occlusive thigh cuffs (VTC)] were applied individually and in combination. In our previous ITD-only study, ITD was investigated for its ability to reduce ICP and cephalic venous congestion in supine and various HDT postures. Internal jugular vein (IJV) ultrasound showed a clear decongestive effect at all postures, however OAE data showed that ITD only caused a phase decrease (tension decrease) in HDT postures. In supine, ITD appeared to INCREASE tension. This paradox leads us to hypothesize that the OAE method is not accurately representing ICP changes with countermeasures (CM), which can alter ME tension through other means, such as ME pressure (MEP) changes. More generally, the exact mechanism for observed OAE response and stimulus phase shifts are not clearly understood, specifically with regard to the effects of MEP. The VCCM study examined the effects of externally-applied MEP on OAE recordings to document the relationship between these parameters. Analysis of these data provide new insights to these OAE mechanisms, in addition to results on CM effectiveness.

Kemp, D.↗

Theoretical studies of the radar properties of the icy Galilean moons of Jupiter

The icy Galilean satellites of Jupiter - Europa, Ganymede, and Callisto - have unusual radar scattering properties compared with those of the terrestrial planets or Earth's Moon. There are three main features of the data that distinguish these targets: (1) the radar cross-section normalized by the geometrical cross-section is an order of magnitude larger than that of any terrestrial planet; (2) the reflected power is almost evenly distributed between two orthogonal polarizations with more power being returned in the same circular polarization as was transmitted whereas virtually all of the power returned from the terrestrial planets is contained in the opposite circular polarization to the one that was transmitted; and (3) the echo power spectra have a broad shape indicating a nearly uniformly radar-bright surface in contrast to the spectra from the terrestrial planets that contain a strong quasi-specular component from the vicinity of the sub-radar point and very little reflected power from the rest of the surface. The normalized radar cross-sections decrease as the areal water ice coverage decreases from Europa to Ganymede to Callisto. Recently, radar echoes from the polar caps of Mars and Mercury, and from Saturn's satellite Titan imply similarly strong cross-sections and have classically unexpected polarization properties and it is also thought that this is due to the presence of ice on the surface. A model called the radar glory model is analyzed and it is shown that the main features of the radar echoes calculated from this model agree well with the observations from all three icy Galilean satellites. This model involves long radar paths in the ice below the surface and special structures in which the refractive index decreases abruptly at a hemispherical boundary. It is not known whether such structures exist or how they could be created, but possible scenarios can be imagined such as the formation of an impact crater followed by deposition of a frost layer followed by a resurfacing event in which a layer of solid ice is placed above the layer of frost. Regardless of the exact geophysical processes required to create such structures, the superior ability of this model to account for all of the important observations with very few adjustable parameters and with no ad hoc assumptions is a compelling argument in support of at least the electromagnetic model. The key features of the electromagnetic model are multiple subsurface scattering events, total internal reflection, and a low degree of randomness imposed on a deterministic geometry that strongly favors backscattering.

Eshleman, Von R.↗

Blind Validation Study of PRICE TruePlanning and SEER-H

Two of the primary parametric costing tools used to estimate the development and production cost of future spacecraft hardware are PRICE TruePlanning by PRICE Systems, and System Estimation and Evaluation of Resources-Hardware (SEER-H) by Galorath. These are standard tools used by NASA and industry to estimate the cost of new aerospace hardware. This study is an independent verification of the accuracy of these tools and was originally published in 2018. The purpose of this presentation is to circulate the findings of that study among the NASA cost community. Both PRICE Systems and Galorath have completed internal validation studies of their parametric cost estimating tools; however, they only provided the results of the studies and did not detail the exact methods used to perform the validation. In the present study, cost estimators used PRICE TruePlanning and SEER-H to estimate the cost of twelve different past NASA science missions. The estimators were prevented from knowing the actual cost of the missions in an effort to minimize cognitive biases. In the present study, SEER-H had an average error of 23%, median error of -0.3%, with a standard deviation of 43%. PRICE had an average error of 52%, median error of 50%, and standard deviation of 45%. Nine of the twelve mission's actual costs fell within the 80% confidence interval of SEER's probabilistic estimates; however, none of the missions fell within PRICE’s confidence intervals. There were several factors independent of PRICE and SEER-H which may have affected the accuracy of the results in the present study including: uncertainty in the technical data used for the estimates, the methods used to estimate uncertainty in spacecraft component mass and numbers of prototypes, and the experience of the estimators.

SEER↗

Progress of the NASA/USGS Lunar Regolith Simulant Project

Beginning in 2004 personnel at MSFC began serious efforts to develop a new generation of lunar simulants. The first two products were a replication of the previous JSC-1 simulant under a contract to Orbitec and a major workshop in 2005 on future simulant development. Beginning in 2006 the project refocused its efforts and approached simulant development in a new and more comprehensive manner, examining new approaches in simulant development and ways to more accurately compare simulants to actual lunar materials. This led to a multi-year effort with five major tasks running in parallel. The five tasks are Requirements, Lunar Analysis, Process Development, Feed Stocks, and Standards. Major progress has been made in all five areas. A substantial draft of a formal requirements document now exists and has been largely stable since 2007. It does evolve as specific details of the standards and Lunar Analysis efforts proceed. Lunar Analysis has turned out to be vastly more difficult than anticipated. After great effort to mine existing published and gray literature, the team has realized the necessity of making new measurements of the Apollo samples, an effort that is currently in progress. Process development is substantially ahead of expectations in 2006. It is now practical to synthesize glasses of appropriate composition and purity. It is also possible to make agglutinate particles in significant quantities. A series of minerals commonly found on the Moon has been synthesized. Separation of mineral constituents from starting rock material is also proceeding. Customized grinding and mixing processes have been developed and tested are now being documented. Identification and development of appropriate feedstocks has been both easier and more difficult than anticipated. The Stillwater Mining Company, operating in the Stillwater layered mafic intrusive complex of Montana, has been an amazing resource for the project, but finding adequate sources for some of the components remains a difficult problem. For example the ratio of clino- to ortho-pyroxenes in the Stillwater is not an exact match for lunar materials. One of the sources being examined as an alternative pyroxene source is the Bushveld Complex in South Africa. Standards have been a major success for the project. The Figure of Merit algorithms have been created, tested, and are being considered for an ISO standard. Agreement has been reached in the community about how to make many of the critical measurements. There remains much work to do: (1) driving down the cost of simulants remains a major obstacle; (2) documentation and cost data analysis have not kept up with progress; (3) educating users in the complexity of the lunar regolith and the use of simulants remains a major task. In summary the project has made enormous progress and is successfully placing simulant development and use on a rigorous, scientifically defensible, engineering basis.

Rickman, Doug↗

Autonomous Medical Care for Exploration

The goal of Autonomous Medical Care (AMC) is to ensure a healthy, well-performing crew which is a primary need for exploration. The end result of this effort will be the requirements and design for medical systems for the CEV, lunar operations, and Martian operations as well as a ground-based crew health optimization plan. Without such systems, we increase the risk of medical events occurring during a mission and we risk being unable to deal with contingencies of illness and injury, potentially threatening mission success. AMC has two major components: 1) pre-flight crew health optimization and 2) in-flight medical care. The goal of pre-flight crew health optimization is to reduce the risk of illness occurring during a mission by primary prevention and prophylactic measures. In-flight autonomous medical care is the capability to provide medical care during a mission with little or no real-time support from Earth. Crew medical officers or other crew members provide routine medical care as well as medical care to ill or injured crew members using resources available in their location. Ground support becomes telemedical consultation on-board systems/people collect relevant data for ground support to review. The AMC system provides capabilities to incorporate new procedures and training and advice as required. The on-board resources in an autonomous system should be as intelligent and integrated as is feasible, but autonomous does not mean that no human will be involved. The medical field is changing rapidly, and so a challenge is to determine which items to pursue now, which to leverage other efforts (e.g. military), and which to wait for commercial forces to mature. Given that what is used for the CEV or the Moon will likely be updated before going to Mars, a critical piece of the system design will be an architecture that provides for easy incorporation of new technologies into the system. Another challenge is to determine the level of care to provide for each mission type. The level of care refers to the amount and type of care one will render based on perceived need and ability. This is in contrast to the standard of care which is the benchmark by which that care is provided. There are certainly some devices and procedures that have unique microgravity or partial gravity requirements such that terrestrial methods will not work. For example, performing CPR on Mars cannot be done in exactly the same way as on Earth because the reduced gravity causes too large a reduction in the forces available for effective compression of the chest. Likewise, fluid behavior in microgravity may require a specialized water filtration and mixing system for the creation of intravenous fluids. This paper will outline the drivers for the design of the medical care systems, prioritization and planning techniques, key system components, and long term goals.

Johnson-Throop, Kathy A.↗