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 199 records · Page 11

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.↗

Tilted outer and inner structures in edge-on galaxies?

ABSTRACT Tilted and warped discs inside tilted dark matter haloes are predicted from numerical and semi-analytical studies. In this paper, we use deep imaging to demonstrate the likely existence of tilted outer structures in real galaxies. We consider two SB0 edge-on galaxies, NGC 4469 and NGC 4452, which exhibit apparent tilted outer discs with respect to the inner structure. In NGC 4469, this structure has a boxy shape, inclined by ΔPA ≈ 3° with respect to the inner disc, whereas NGC 4452 harbours a discy outer structure with ΔPA ≈ 6°. In spite of the different shapes, both structures have surface brightness profiles close to exponential and make a large contribution (∼30 per cent) to the total galaxy luminosity. In the case of NGC 4452, we propose that its tilted disc likely originates from a former fast tidal encounter (probably with IC 3381). For NGC 4469, a plausible explanation may also be galaxy harassment, which resulted in a tilted or even a tumbling dark matter halo. A less likely possibility is accretion of gas-rich satellites several Gyr ago. New deep observations may potentially reveal more such galaxies with tilted outer structures, especially in clusters. We also consider galaxies, mentioned in the literature, where a central component (a bar or a bulge) is tilted with respect to the stellar disc. According to our numerical simulations, one of the plausible explanations of such observed ‘tilts’ of the bulge/bar is a projection effect due to a not exactly edge-on orientation of the galaxy coupled with a skew angle of the triaxial bulge/bar.

Mosenkov, Aleksandr V.↗

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.↗

Hetero-Orbital Two-Component Fractional Quantum Hall States in Bilayer Graphene

A two-dimensional electron system exposed to a strong magnetic field produces a plethora of strongly interacting fractional quantum Hall (FQH) states, the complex topological orders of which are revealed through exotic emergent particles, such as composite fermions, and fractionally charged Abelian and non-Abelian anyons. Much insight has been gained by the study of multicomponent FQH states, where spin and pseudospin indices of the electron contribute additional correlation. Traditional multicomponent FQH states develop in situations where the components share the same orbital states and the resulting interactions are pseudospin independent; this homo-orbital nature is also crucial to their theoretical understanding. Here, we study “hetero-orbital” two-component FQH states, in which the orbital index is part of the pseudospin, rendering the multicomponent interactions strongly SU(2) anisotropic in the pseudospin space. Such states, obtained in bilayer graphene at the isospin transition between 𝑁 = 0 and 𝑁 = 1 electron Landau levels, are markedly different from previous homo-orbital two-component FQH states. In particular, we observe strikingly different behaviors for the parallel-vortex and reverse-vortex attachment composite fermion states, and an anomalously strong two-component 2/5 state over a wide range of magnetic field before it abruptly disappears at a high field. Our findings, combined with detailed theoretical calculations, reveal the surprising robustness of the hetero-orbital FQH effects, significantly enriching our understanding of FQH physics in this novel regime.

Bilayer graphene↗

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↗

Explicit physics-informed neural networks for nonlinear closure: The case of transport in tissues

In upscaling methods, closures for nonlinear problems present a well-known challenge. While a number of theoretical methods have been proposed for handling such closures, nonlinearities still remain a significant obstacle for many problems. In this work, we use a combination of formal upscaling and data-driven machine learning for explicitly closing a nonlinear transport and reaction process in multiscale tissues. The classical effectiveness factor model is used to formulate the macroscale reaction kinetics. We train a multilayer perceptron network using training data generated by direct numerical simulations over microscale examples. Once trained, the network is used in an algorithm for numerically solving the upscaled (coarse-grained) differential equation describing mass transport and reaction in two example tissues. The network is described as being explicit in the sense that the network is trained using macroscale concentrations and gradients of concentration as components of the feature space rather than incorporating them as part of a constraint in the optimization process. Network training and solutions to the macroscale transport equations were computed for two different tissues. The two tissue types (brain and liver) exhibit markedly different geometrical complexity and spatial scale (cell size and sample size). The upscaled solutions for the average concentration are compared with numerical solutions derived from the microscale concentration fields by a posteriori averaging. There are three outcomes of this work of particular note. 1) Our overall approach results in an upscaled nonlinear PDE. The PDE is closed using a neural network, and our approach results in the definition of the classical effectiveness factor for effecting closure. 2) We identify particular source terms for the closure problem that are important for representing the structure of the closure. These source terms involve macroscale concentrations and their gradients. We adopt these source terms to use as explicit features in the learning algorithm. We find the trained networks that include the macroscale source terms generate models that are able to predict the correction factor with increased fidelity over those that do not. 3) We find that the trained network exhibits good generalizability, and it is able to predict the effectiveness factor with high fidelity for realistically-structured tissues despite the significantly different scale and geometrical complexity of the two example tissue types. This latter result emphasizes our purposeful connection between conventional averaging methods with the use of machine learning for closure; this contrasts with some machine learning methods for upscaling where the exact form of the macroscale equation remains unknown.

97 MATHEMATICS AND COMPUTING↗

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↗

Understanding the effect of density functional choice and van der Waals treatment on predicting the binding configuration, loading, and stability of amine-grafted metal organic frameworks

Metal organic frameworks (MOFs) are crystalline, three-dimensional structures with high surface areas and tunable porosities. Made from metal nodes connected by organic linkers, the exact properties of a given MOF are determined by node and linker choice. MOFs hold promise for numerous applications, including gas capture and storage. M2(4,4′-dioxidobiphenyl-3,3′-dicarboxylate)—henceforth simply M2(dobpdc), with M = Mg, Mn, Fe, Co, Ni, Cu, or Zn—is regarded as one of the most promising structures for CO2 capture applications. Further modification of the MOF with diamines or tetramines can significantly boost gas species selectivity, a necessity for the ultra-dilute CO2 concentrations in the direct-air capture of CO2. There are countless potential diamines and tetramines, paving the way for a vast number of potential sorbents to be probed for CO2 adsorption properties. The number of amines and their configuration in the MOF pore are key drivers of CO2 adsorption capacity and kinetics, and so a validation of computational prediction of these quantities is required to suitably use computational methods in the discovery and screening of amine-functionalized sorbents. In this work, we study the predictive accuracy of density functional theory and related calculations on amine loading and configuration for one diamine and two tetramines. In particular, we explore the Perdew–Burke–Ernzerhof (PBE) functional and its formulation for solids (PBEsol) with and without the Grimme-D2 and Grimme-D3 pairwise corrections (PBE+D2/3 and PBEsol+D2/3), two revised PBE functionals with the Grimme-D2 and Grimme-D3 pairwise corrections (RPBE+D2/3 and revPBE+D2/3), and the nonlocal van der Waals correlation (vdW-DF2) functional. We also investigate a universal graph deep learning interatomic potential’s (M3GNet) predictive accuracy for loading and configuration. These results allow us to identify a useful screening procedure for configuration prediction that has a coarse component for quick evaluation and a higher accuracy component for detailed analysis. Our general observation is that the neural network-based potential can be used as a high-level and rapid screening tool, whereas PBEsol+D3 gives a completely qualitatively predictive picture across all systems studied, and can thus be used for high accuracy motif predictions. We close by briefly exploring the predictions of relative thermal stability for the different functionals and dispersion corrections.

Chemistry↗

Transporting Non-Compliant TRU Wastes in the OPTIMUS{sup TM} Packages - 20493

The OPTIMUS{sup TM} product line of packages was developed with the purpose of providing a versatile and modular packaging option for shipping problematic wastes and fissile material contents. The product line includes a high-activity design in the OPTIMUS-H and a low activity design in the OPTIMUS-L. The two variants utilize a singular containment vessel design, that is enclosed in different protective components that offer different levels of shielding along with thermal and impact protection. The thick shielding provided by the OPTIMUS-H packaging allows for higher activity contents than the OPTIMUS-L to be transported. But the lightweight design of the OPTIMUS-L packaging allows for more packages to be transported in a single shipment. The containment vessel design shared between the two packages provides leak-tight containment of all radioactive contents along with the capabilities to easily leak test and the option to inert and backfill the contents, as necessary, prior to each shipment of the package. Together the OPTIMUS package designs offer significant flexibility and are capable of handling a wide variety of waste materials and other radioactive contents. One of the primary contents covered in the initial design effort for these packages is TRU waste contents beyond the standard waste materials that are compliant with the WIPP Waste Acceptance Criteria. More specifically, the initial content of interest is TRU Waste drums containing sealed containers with potentially flammable gases. While there are multiple other packaging options available for transporting standard WIPP compliant TRU wastes, the goal of the OPTIMUS packages is to offer a superior option for transporting both standard TRU wastes as well as the more problematic non-compliant wastes, among other contents. The non-compliant TRU waste contents explicitly included for the OPTIMUS packages are standard aerosol cans and DOT 3E lecture bottles. However, the methods applied for the package containment can be easily adjusted to cover other sealed container types with potentially flammable gases. The primary challenges with including these items are both the obvious issue of potential for flammable gases (e.g. aerosol propellants) in the package, but also the uncertainty in the exact state of the contents. Because these desired contents are waste materials, the characterization of the materials present in the waste may include some uncertainties. For example, though it may be known that there are one or more aerosol cans in a TRU waste drum, it may not be known if this can is full, spent, or anywhere in between. Also there may be equal uncertainty in the contents of the aerosol can, specifically in the potentially flammable propellant remaining in the can. To include these types of non-compliant items as acceptable contents for the package, the methods utilized in the safety analyses of the OPTIMUS packages must consider all of the uncertainties in the characterization of the waste. The methods utilized cover a range of concerns with transporting contents of this nature including pressure buildup, gas generation, and gaseous combustion for demonstrating the containment of the package. The ability to transport these problematic contents in the OPTIMUS packages gives TRU waste generating sites the option to relocate these drums to an offsite location where the non-compliant items can be properly managed through removal or destruction. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

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↗