Search NASA⌕ Search

SEARCH · Search NASA

Results for “Outer approximation”

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 991 records · Page 55
All sourcesNASA NTRSDOE OSTI
NASA NTRS2004

Understanding Columbia's Reentry Problems

Soon after the Space Shuttle Columbia accident occurred last year, a group of CFD analysts from NASA centers and private industry was organized to help determine the cause of the accident. This group was under the direction of the Applied Aeroscience and CFD Branch of the Aeroscience and Flight Mechanics Division at the Johnson Space Center. For external flow simulations, noncommercial CFD codes that specialize in hypersonic or high Mach number flows were used. These tools were used to determine heating rates, pressures, and temperatures for a large number of vehicle damage scenarios. Lockheed Martin Space Operations was called upon to provide CFD support in the area of internal flows within the shuttle wing cavity, and or these simulations, FLUENT 6.1 was chosen. Two large-scale, simplified models were run to understand the flow patterns once a breach of the internal wing cavity was initiated. The results were primarily used to visualize flow patterns within the wing cavity. The first CFD model included the entire left wing without the wheel. well cavity. The purpose of the first model, which did not include the reinforced carbon-carbon (RCC) cavity along the wing leading edge, was to visualize the flow field within the wing cavity immediately after the leading edge spar breach, This model assumed that the flow coming into the wing cavity was nom1 to the spar. It included all o f the primary vents that allow for flow between the main cavities of the wing. A six-inch diameter hole was modeled in the spar at the approximate location where the spar breach was judged to have occurred, which was between RCC panels 8 and 9. The results of the modeling showed that at this location, the high temperature, high velocity gas stream entering the wing cavity impinged on the outboard wheel well cavity. Instrumentation in the Shuttle wheel well cavity registered abnormal temperatures during reentry, so the FLUENT results helped support the conclusion of the accident investigation team that the spar breach was in the RCC panel 8-9 area, and that the initial spar breach was likely entering the wing cavity normal to the spar. This model also showed that the flow entering the wing cavity tended to swirl within the cavity just outboard of the wheel well, and did not initially penetrate further into the rear cavities of the wing. The second CFD model was a 2-D simulation of the left wing cavity and the RCC cavity. It was used to visualize the flow through the RCC breach, through the wing spar breach, and into the wing cavity directly outboard of the wheel well. The purpose of this model was to verify whether or not it was possible for the flow to come into the wing cavity normal to the leading edge spar. The results from this 2-D model showed that the internal structure behind the RCC panels probably deflected the flow entering the RCC cavity so that it impinged normal to the spar. As in the 3-D model, this deflected flow stream was found to impinge on the wheel well outer wall. This model again supported the conclusions regarding the location of the spar breach and how the flow behaved inside the wing cavity.

Paul, Thomas H.↗
NASA NTRS2006

Understanding Columbia's Reentry Problems

Soon after the Space Shuttle Columbia accident occurred last year, a group of CFD analysts from NASA centers and private industry was organized to help determine the cause of the accident. This group was under the direction of the Applied Aeroscience and CFD Branch of the Aeroscience and Flight Mechanics Division at the Johnson Space Center. For external flow simulations, noncommercia2 CFD codes that specialize in hypersonic or high Mach number flows were used. These tools were used to determine heating rates, pressures, and temperatures for a large number of vehicle damage scenarios. Lockheed Martin Space Operations was called upon to provide CFD support in the area of internal flows within the shuttle wing cavity, and for these simulations, FLUENT 6.1 was chosen. Two large-scale, simplified models were m to understand the flow patterns once a breach of the internal wing cavity was initiated. The results were primarily used to visualize flow patterns within the wing cavity. The first CFD model included the entire lee wing without the wheel well cavity. The purpose of the first model, which did not include the reinforced carbon-carbon (RCC) cavity along the wing leading edge, was to visualize the flow field within the wing cavity immediately after the leading edge spar breach, This model assumed that the flow coming into the wing cavity was normal to the spar. It included all of the primary vents that allow for flow between the main cavities of the wing. A six-inch diameter hole was modeled in the spar at the approximate location where the spar breach was judged to have occurred, which was between RCC panels 8 and 9. The results of the modeling showed that at this location, the high temperature, high velocity gas stream entering the wing cavity impinged on the outboard wheel well cavity. Instrumentation in the Shuttle wheel well cavity registered abnormal temperatures during reentry, so the FLUENT results helped support the conclusion of the accident investigation team that the spar breach was in the RCC panel 8-9 area, and that the initial spar breach was likely entering the wing cavity normal to the spar. This model also showed that the flow entering the wing cavity tended to swirl within the cavity just outboard of the wheel well, and did not initially penetrate further into the rear cavities of the wing. The second CFD model was a 2-D simulation of the left wing cavity and the RCC cavity. It was used to visualize the flow through the RCC breach, through the wing spar breach, and into the wing cavity directly outboard of the wheel well. The purpose of this model was to verify whether or not it was possible for the flow to come into the wing cavity normal to the leading edge spar. The results from this 2-D model showed that the internal structure behind the RCC panels probably deflected the flow entering the RCC cavity so that it impinged normal to the spar. As in the 3-D model, this deflected flow stream was Found to impinge on the wheel well outer wall. This model again supported the conclusions regarding the location of the spar breach and how the flow behaved inside the wing cavity.

Paul, Thomas↗
NASA NTRS1996

Studies using IMP, Voyager and Pioneer cosmic ray data to determine the size of the heliosphere

The purpose of this project was to use the cosmic ray data from the IMP, Voyager and Pioneer spacecraft in the heliosphere out to approximately 65 AU to estimate the size of the heliosphere. We used several techniques to develop a consistent picture of the size of the heliosphere. The first method used a response function approach which determined the intensity as a function of time by scaling the modulation effect as they move outward and eventually reach the boundary of the heliosphere. In this model the effects of transient cosmic ray disturbances is included. A second approach using the perturbation method in which drifts are considered as a perturbation to the standard diffusion-convection modulation models was not fully developed. In a third approach the location of the modulation boundary beyond the termination shock was estimated using observations of the intensity and radial gradients between Voyager 2 and Pioneer 10 along with new estimates of the interstellar intensity of more than 70 MeV galactic cosmic rays. Using this method we found that for 7 years, from 1983 to 1990, the modulation boundary remained constant at 83 +/- 5 AU. We infer from these studies that a modulation boundary can be estimated only by extrapolating the observed radial gradients when the solar magnetic field polarity is such that cosmic-ray particles are drifting in the heliosphere inward toward the Earth along the neutral sheet. The boundary distance is larger than the estimates of the location of the termination shock at 67 +/- 5 AU using the same method. Two other studies partially supported by this grant are attached. The first deals with the recovery period of the greater than 70 MeV cosmic rays in the outer heliosphere from 1992-1995. In the second paper we compare the rigidity dependence of the 11-year cosmic ray variation at the Earth in two cycles of opposite solar magnetic field polarity.

Lockwood, John A.↗
NASA NTRS2001

Estimating the Libration of Mercury by Remote Sensing of Gravity and Altimetry

One of the most intriguing scientific questions about Mercury is whether the planet possesses a liquid core. One way that this can be ascertained is through the measurement of the longitudinal libration of Mercury about its rotation axis. The MESSENGER (Mercury Surface, Space, Environment, Geochemistry, and Ranging) mission to Mercury will attempt this measurement by deriving the rotation of the planet from remotely sensed observations of the Mercury gravity field and the planet's shape. For a fully fluid core the libration is expected to be an approximately 350-m-amplitude oscillation (at the equator) about the mean rotation of the planet. Thus, the shape of the planet undergoes a librational oscillation in the same way that features on the surface participate in the motion. In this sense altimetric topography can be considered analogous to imaging, with the important exception that for altimetry, the long wavelengths corresponding to low degrees and orders in a spherical harmonic expansion, can be precisely determined from orbital measurements. The rotation of the planet thus introduces a time dependence of the topographic image in an inertial (i.e. fixed) reference frame and the libration is a variation in the planetary rotation rate, manifest as an oscillation of the topographic image. Similarly, gravity measures the distribution of mass within a planet and as such can be viewed as providing a vertically integrated 'image' of the internal density structure. For a librating planet the gravitational field will also reveal a variation in the rotation rate, manifest as an oscillation of the field with respect to an inertial reference frame. While the long-wavelength gravity field contains contributions from the radial distribution of mass with depth, the long wavelength terms are most sensitive to the mass distribution of the deep interior. Thus, both the topography and the gravity independently can be used to determine the (irregular) rotation of a planet, although the topographical method measures the rotation of the crust while the gravity measures the rotation of the distributed mass, and these are not necessarily identical. For example, for a planet that contains a fluid outer core that is effectively decoupled from the mantle the gravity field will not librate in the manner of a solid body. The libration obtained from the variable rotation of the gravity field will represent a combination of contributions due to a librating mantle and a differentially rotating fluid core. Differential rotation of the Earth's inner core had been theoretically predicted from three dimensional geodynamical models and subsequently reported from seismic observations. A difference between the libration of the lithosphere as determined from altimetry and the deep interior as determined from gravity, if it could be detected, might provide further insight into the nature of core-mantle coupling. Additional information is contained in the original extended abstract.

Smith, D. E.↗
NASA NTRS1994

Hydrogen and related materials at high density: Physics, chemistry and planetary implications

Recent studies of low-Z molecular materials including hydrogen to multimegabar pressures (less than 300 GPa) have uncovered a range of phenomena relevant to understanding the nature of the interiors of the outer planets and their satellites. Synchrotron x ray diffraction measurements (to 42 GPa) have been used to determine the crystal structure of the solid (hexagonal-close packed) and equation of state. Sound velocities in fluid and solid hydrogen (to 24 GPa) have been inverted to obtain elastic constants and aggregate bulk and shear moduli. In addition, an improved intermolecular potential has been determined which fits both static and shock-wave data. Use of the new potential for the molecular envelope of Jupiter suggests the need for major revisions of existing Jovian models or a reanalysis of reported free oscillations for the planet. Studies at higher pressures (greater than 100 GPa) reveal a sequence of pressure-induced symmetry-breaking transitions in molecular hydrogen, giving rise to three high-pressure phases (1, 2, and 3). Phase 1 is the rotationally disordered hcp phase which persists from low pressure to well above 100 GPa at high temperature (e.g., 300 K). Phase 2 is a low-temperature, high-pressure phase (transition at 100 GPa and 77 K in H2) with spectral features indicative of partial rotational ordering and crystallographic distortion. The transition to Phase 3 at 150 GPa is accompanied by a weakening of the molecular bond, gradual changes in orientational ordering, strong enhancement of the infrared intramolecular vibrational absorption, and strong intermolecular interactions similar to those of ambient-pressure network solids. Studies of the phase diagram reveal a triple point near 130 K and 160 GPa. Higher pressure measurements of vibrational spectra place a lower bound of approximately 250 GPa on the predicted transition pressure for dissociation of molecular hydrogen to form a monatomic metal.

Hemley, R. J.↗
DOE OSTI2024

AEOLUS: Advances in Experimental Design, Optimal Control, and Learning for Uncertain Complex Systems

Sustained advances in the mathematics of modeling and simulation have resulted in the capability today for routine simulation of a number of large scale complex DOE-relevant systems. As remarkable as this capability for solving the so-called forward problem is, it is typically only the first step-an inner loop within an outer loop that explores the simulation model's parameter space and decision space to characterize uncertainty in the model's predictions, learn unknown model parameters from data, design the most informative experiments, determine optimal control strategies, and create optimal designs. Broadly, what unifies all of these outer loop problems is that they are, in one form or another, optimization problems over parameter/control/design space that are constrained by complex uncertain models. To fully realize the power of scientific simulation as a basis for scientific discovery, technological innovation, and rational decision-making, it is imperative to move beyond simulation to tackle the outer loop of optimization for learning from data, experimental design, and control with complex uncertain models. When the models under consideration are large-scale and complex, and when the optimization variable and uncertain parameter spaces are high (or infinite) dimensional, this constitutes a grand challenge of the highest order, and is intractable with conventional methods. To overcome these challenges, the AEOLUS Center was established to develop a unified mathematical, computational, and statistical framework for (1) Learning predictive models from complex data via Bayesian inference and optimization, and (2) Optimizing experiments, processes, and designs using the resulting uncertain models. These problems are intractable with conventional methods, for several reasons: (1) The simulation problems that govern the inner loops of the optimization problems are expensive to execute (due to severe nonlinearity, heterogeneity, multiphysics/multiscale coupling); (2) The optimization variable and uncertain parameter spaces are high dimensional, often stemming from discretizations of infinite dimensional fields such as initial conditions, sources, or material properties. We argue that the key to overcoming these challenges is to develop new mathematical, computational, and statistical methods that exploit the structure of the Bayesian inference and optimization problems mediated by their underlying complex uncertain models. This structure includes the regularity, sparsity, geometry, low intrinsic dimensionality, and multifidelity nature of the maps from uncertain parameter/optimization variable spaces to the specific objectives targeted: Bayesian inference, optimal experimental design, and optimal control design. Black box methods developed as generic tools are incapable of exploiting this structure. To be successful, we must create, integrate, and cross-fertilize ideas across multiple areas of applied math--including approximation theory, Bayesian inference, data science, experimental design, information theory, machine learning, model reduction, optimal control theory, parallel algorithms, PDE-constrained optimization, randomized algorithms, stochastic optimization, and uncertainty quantification--all while exploiting the structure of the problems at hand. With this goal in mind, we have marshaled a team of leading authorities in these areas. While the methods we develop will be broadly applicable across a wide spectrum of DOE problems in which experiments inform models and the systems those models describe must be optimized under uncertainty, we have chosen a specific area, advanced manufacturing and materials, to drive our work. AMM is characterized by complex models across multiple scales, and is a rich source of challenging problems in inference, experimental design, and optimal control, requiring multifaceted and integrated advances in applied mathematics. As such, AMM serves as an excellent vehicle to motivate and demonstrate the advances in applied mathematics developed by our center.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1995

Low Intensity Low Temperature (LILT) Measurements and Coefficients on New Photovoltaic Structures

Past NASA missions to Mars, Jupiter and the outer planets were powered by radioisotope thermal generators (RTGs). Although these devices proved to be reliable, their high cost and highly toxic radioactive heat source has made them far less desirable for future planetary missions. This has resulted in a renewed search for alternate energy sources, some of them being photovoltaics (PV) and thermophotovoltaics (TPV). Both of these alternate energy sources convert light/thermal energy directly into electricity. In order to create a viable PV data base for planetary mission planners and cell designers, we have compiled low intensity low temperature (LILT) I-V data on single junction and multi-junction high efficiency solar cells. The cells tested here represent the latest photovoltaic technology. Using this LILT data to calculate Short Circuit Current (I(sub sc)), Open Circuit Voltage (V(sub os)), and Fill Factor (FF) as a function of temperature and intensity, an accurate prediction of cell performance under the AM0 spectrum can be determined. When combined with QUantum efficiency at Low Temperature (QULT) data, one can further enhance the data by adding spectral variations to the measurements. This paper presents an overview of LILT measurements and is only intended to be used as a guideline for material selection and performance predictions. As single junction and multi-junction cell technologies emerge, new test data must be collected. Cell materials included are Si, GaAs/Ge, GaInP/GaAs/GaAs, InP, InGaAs/InP, InP/InGaAs/InP, and GaInP. Temperatures range down to as low as -180 C and intensities range from 1 sun down to 0.02 suns. The coefficients presented in this paper represent experimental results and are intended to provide the user with approximate numbers.

Scheiman, David A.↗
NASA NTRS1993

Enhanced decoding for the Galileo S-band mission

A coding system under consideration for the Galileo S-band low-gain antenna mission is a concatenated system using a variable redundancy Reed-Solomon outer code and a (14,1/4) convolutional inner code. The 8-bit Reed-Solomon symbols are interleaved to depth 8, and the eight 255-symbol codewords in each interleaved block have redundancies 64, 20, 20, 20, 64, 20, 20, and 20, respectively (or equivalently, the codewords have 191, 235, 235, 235, 191, 235, 235, and 235 8-bit information symbols, respectively). This concatenated code is to be decoded by an enhanced decoder that utilizes a maximum likelihood (Viterbi) convolutional decoder; a Reed Solomon decoder capable of processing erasures; an algorithm for declaring erasures in undecoded codewords based on known erroneous symbols in neighboring decodable words; a second Viterbi decoding operation (redecoding) constrained to follow only paths consistent with the known symbols from previously decodable Reed-Solomon codewords; and a second Reed-Solomon decoding operation using the output from the Viterbi redecoder and additional erasure declarations to the extent possible. It is estimated that this code and decoder can achieve a decoded bit error rate of 1 x 10(exp 7) at a concatenated code signal-to-noise ratio of 0.76 dB. By comparison, a threshold of 1.17 dB is required for a baseline coding system consisting of the same (14,1/4) convolutional code, a (255,223) Reed-Solomon code with constant redundancy 32 also interleaved to depth 8, a one-pass Viterbi decoder, and a Reed Solomon decoder incapable of declaring or utilizing erasures. The relative gain of the enhanced system is thus 0.41 dB. It is predicted from analysis based on an assumption of infinite interleaving that the coding gain could be further improved by approximately 0.2 dB if four stages of Viterbi decoding and four levels of Reed-Solomon redundancy are permitted. Confirmation of this effect and specification of the optimum four-level redundancy profile for depth-8 interleaving is currently being done.

Dolinar, S.↗
NASA NTRS2011

Wheel-Based Ice Sensors for Road Vehicles

Wheel-based sensors for detection of ice on roads and approximate measurement of the thickness of the ice are under development. These sensors could be used to alert drivers to hazardous local icing conditions in real time. In addition, local ice-thickness measurements by these sensors could serve as guidance for the minimum amount of sand and salt required to be dispensed locally onto road surfaces to ensure safety, thereby helping road crews to utilize their total supplies of sand and salt more efficiently. Like some aircraft wing-surface ice sensors described in a number of previous NASA Tech Briefs articles, the wheelbased ice sensors are based, variously, on measurements of changes in capacitance and/or in radio-frequency impedance as affected by ice on surfaces. In the case of ice on road surfaces, the measurable changes in capacitance and/or impedance are attributable to differences among the electric permittivities of air, ice, water, concrete, and soil. In addition, a related phenomenon that can be useful for distinguishing between ice and water is a specific transition in the permittivity of ice at a temperature- dependent frequency. This feature also provides a continuous calibration of the sensor to allow for changing road conditions. Several configurations of wheel-based ice sensors are under consideration. For example, in a simple two-electrode capacitor configuration, one of the electrodes would be a circumferential electrode within a tire, and the ground would be used as the second electrode. Optionally, the steel belts that are already standard parts of many tires could be used as the circumferential electrodes. In another example (see figure), multiple electrodes would be embedded in rubber between the steel belt and the outer tire surface. These electrodes would be excited in alternating polarities at one or more suitable audio or radio frequencies to provide nearly continuous monitoring of the road surface under the tire. In still another example, one or more microwave stripline(s) or coplanar waveguide(s) would be embedded in a tire near its outer surface; in comparison with lower-frequency capacitive devices, a device of this type could be more sensitive.

Arndt, G. Dickey↗
NASA NTRS1992

Search for evidence of a clock related to the solar 154 day complex of periodicities

Evidence that has recently been compiled (Bai and Sturrock 1991) indicates that the enigmatic 154-day periodicity in solar activity may be viewed as part of a complex of periodicities that are approximate multiples of 25.8 days, suggesting that the Sun contains a 'clock' with frequency in the range 440 to 463 nano Hz. The clock may comprise either an oscillator or a rotator, each of which may be either real or virtual. We have reconsidered a previous spectrum analysis of the Zurich sunspot-number sequence by Knight, Schatten, and Sturrock (1979) which revealed a sharp, persistent and significant periodicity with a period of 12.072 days, corresponding to a frequency of about 958.8 nano Hz. This periodicity may be regarded as the (second) upper sideband of the second harmonic (2nu(sub R) + 2nu(sub E)) of a fundamental frequency of 447.7 nano Hz that is clearly within the search band. In this expression, nu(sub R) is the sidereal frequency of the hypothetical rotator and nu(sub E) is the frequency (31.69 nano Hz) of the Earth in its orbital motion around the Sun. In analyzing sunspot area data derived from the Greenwich data set, and on noting that any frequency is defined only to within the Nyquist frequency, we find clear evidence not only for the upper sideband of the second harmonic, but also for the second harmonic (2nu(sub R)) and the lower sideband of the second harmonic (2nu(sub R) - 2nu(sub E)). There is no strong peak at the fundamental frequency in the Greenwich data, but there is in the Zurich sunspot data. The effect of a linear oscillator is, to the lowest order in the amplitude, the same as the combined effect of two rotators of opposite polarities. A rotator that has arbitrary orientation with respect to the ecliptic may influence the outer layers of the Sun and thereby modulate the occurrence of solar activity such as sunspots. By analyzing a simple model, we find that such a rotator would influence surface activity in such a way that the spectrum of a 'signal' (such as the record of sunspots), as seen from the Earth, would contain components with frequencies that are certain integral combinations of nu(sub R) and nu(sub E). The amplitudes of the various components depend sensitively on theta, the angle between the axis of the rotator and the axis of the Earth's orbital motion. This simple model therefore offers a kinematical (but not dynamical) interpretation of the sunspot spectrum. The present analysis, while offering support of our conjectures that the Sun contains a clock that regulates the 154-day complex of periodicities, cannot distinguish between an osillator or a rotator (that might be a traveling wave), nor between a real rotator or a virtual rotator (that might be an apparent traveling wave due to the aliasing effect of an oscillator in a rotating system). Further analysis of sunspot and other data sets will be required to confirm the existence of such clock and (if it is real) to determine its physical nature.

Sturrock, P. A.↗
NASA NTRS2015

The Compositional Evolution of C/2012 S1 (ISON) from Ground-Based High-Resolution Infrared Spectroscopy as Part of a Worldwide Observing Campaign

Volatile production rates, relative abundances, rotational temperatures, and spatial distributions in the coma were measured in C/2012 S1 (ISON) using long-slit high-dispersion (lambda/delta lambda approximately 2.5 times 10 (sup 4)) infrared spectroscopy as part of a worldwide observing campaign. Spectra were obtained on Universal Time 2013 October 26 and 28 with NIRSPEC (Near Infrared Spectrometer) at the W.M. Keck Observatory, and Universal Time 2013 November 19 and 20 with CSHELL (Cryogenic Echelle Spectrograph) at the NASA IRTF (Infrared Telescope Facility). H2O was detected on all dates, with production rates increasing markedly from (8.7 plus or minus 1.5) times 10 (sup 27) molecules per second on October 26 (Heliocentric Distance = 1.12 Astronomical Units) to (3.7 plus or minus 0.4) times 10 (sup 29) molecules per second on November 20 (Heliocentric Distance = 0.43 Astronomical Units). Short-term variability of H2O production is also seen as observations on November 19 show an increase in H2O production rate of nearly a factor of two over a period of about 6 hours. C2H6, CH3OH and CH4 abundances in ISON (International Scientific Optical Network) are slightly depleted relative to H2O when compared to mean values for comets measured at infrared wavelengths. On the November dates, C2H2, HCN and OCS abundances relative to H2O appear to be within the range of mean values, whereas H2CO and NH3 were significantly enhanced. There is evidence that the abundances with respect to H2O increased for some species but not others between October 28 (Heliocentric Distance = 1.07 Astronomical Units) and November 19 (Heliocentric Distance = 0.46 Astronomical Units). The high mixing ratios of H2CO to CH3OH and C2H2 to C2H6 on November 19, and changes in the mixing ratios of some species with respect to H2O between October 28 to November 19, indicates compositional changes that may be the result of a transition from sampling radiation-processed outer layers in this dynamically new comet to sampling more pristine natal material as the outer processed layer was increasingly eroded and the thermal wave propagated into the nucleus as the comet approached perihelion for the first time. On November 19 and 20, the spatial distribution for dust appears asymmetric and enhanced in the antisolar direction, whereas spatial distributions for volatiles (excepting CN) appear symmetric with their peaks slightly offset in the sunward direction compared to the dust. Spatial distributions for H2O, HCN, C2H6, C2H2, and H2CO on November 19 show no definitive evidence for significant contributions from extended sources; however, broader spatial distributions for NH3 and OCS may be consistent with extended sources for these species. Abundances of HCN and C2H2 on November 19 and 20 are insufficient to account for reported abundances of CN and C2 in ISON near this time. Differences in HCN and CN spatial distributions are also consistent with HCN as only a minor source of CN in ISON on November 19 as the spatial distribution of CN in the coma suggests a dominant distributed source that is correlated with dust and not volatile release. The spatial distributions for NH3 and NH2 are similar, suggesting that NH3 is the primary source of NH2 with no evidence of a significant dust source of NH2; however, the higher production rates derived for NH3 compared to NH2 on November 19 and 20 remain unexplained. This suggests a more complete analysis that treats NH2 as a distributed source and accounts for its emission mechanism is needed for future work.

Russo, N. Dello↗
NASA NTRS1997

Applying Space Technology to Enhance Control of an Artificial Arm

At the present time, myoelectric prostheses perform only one function of the hand: open and close with the thumb, index and middle finger coming together to grasp various shaped objects. To better understand the limitations of the current single-function prostheses and the needs of the individuals who use them, The Institute for Rehabilitation and Research (TIRR), sponsored by the National Institutes of Health (August 1992 - November 1994), surveyed approximately 2500 individuals with upper limb loss. When asked to identify specific features of their current electric prosthesis that needed improvement, the survey respondents overwhelmingly identified the lack of wrist and finger movement as well as poor control capability. Simply building a mechanism with individual finger and wrist motion is not enough. Individuals with upper limb loss tend to reject prostheses that require continuous visual monitoring and concentration to control. Robotics researchers at NASA's Johnson Space Center (JSC) and Rice University have made substantial progress in myoelectric teleoperation. A myoelectric teleoperation system translates signals generated by an able-bodied robot operator's muscles during hand motions into commands that drive a robot's hand through identical motions. Farry's early work in myoelectric teleoperation used variations over time in the myoelectric spectrum as inputs to neural networks to discriminate grasp types and thumb motions. The resulting schemes yielded up to 93% correct classification on thumb motions. More recently, Fernandez achieved 100% correct non-realtime classification of thumb abduction, extension, and flexion on the same myoelectric data. Fernandez used genetic programming to develop functions that discriminate between thumb motions using myoelectric signal parameters. Genetic programming (GP) is an evolutionary programming method where the computer can modify the discriminating functions' form to improve its performance, not just adjust numerical coefficients or weights. Although the function development may require much computational time and many training cases, the resulting discrimination functions can run in realtime on modest computers. These results suggest that myoelectric signals might be a feasible teleoperation medium, allowing an operator to use his or her own hand and arm as a master to intuitively control an anthropomorphic robot in a remote location such as outer space.

Atkins, Diane↗
DOE OSTI2020

Fast and Accurate Core Analysis by the Full-Immersion Pressure-Pulse Decay: Part 2—Practice and Demonstration

Summary This study continues the work of presenting a novel approach for making petrophysical assessments of tight core samples. This method, the full-immersion pressure-pulse decay (Hannon 2019), involves applying a rapid increase in pressure in a chamber surrounding the entire outer surface area of a cylindrical sample, shutting the system in, and monitoring the pressure decay in the chamber as it reaches a new equilibrium. A precursor article covered the numerical simulator designed to model flow through the sample, demonstrating its performance and accuracy in addition to providing a first-order comparison between the speed and shape of the pressure-decay responses of the full-immersion method with those of other similar transient methods. This study covers the parameter-estimation procedure and experimental verification through a proof-of-concept laboratory investigation. The investigations provided here demonstrate that under appropriate, achievable experimental conditions, the pressure data can be analyzed in such a way that returns an estimate of the porosity and apparent permeabilities both parallel and perpendicular to bedding from a single test performed on a single cylindrical sample. After determining these experimental conditions (the uniqueness window), this report outlines a data-inversion strategy to estimate the petrophysical properties (porosity, horizontal permeability, and vertical permeability) from each test. This strategy is put to the test through comparisons with measurements performed by a commercial core laboratory. A common set of samples recovered from an outcrop of a tight-gas sandstone formation were investigated using the full-immersion method, and their results are compared with those from conventional steady-state measurement procedures performed by the commercial laboratory. Comparisons between petrophysical characterizations of these samples, which had permeabilities between 25 nd and 2.3 μd, demonstrated close agreement in most cases. However, whereas steady-state measurements performed at the professional laboratory required 4 to 5 hours of testing time per measurement of a single permeability, similar assessments using the full-immersion technique, requiring approximately 5 to 10 minutes to complete, returned estimates of the horizontal and vertical permeability simultaneously. Additional analyses are provided to determine principal reasons of discrepancies in instances where agreement was not as strong. Based on lessons learned from these experiences, the report closes with suggestions on areas of improvement in the experimental approach. Once complete, these developments should propel this technology to fill a critical need to determine petrophysical properties (porosity and permeability) of tight rocks in a time-efficient manner and in a way that does not compromise their accuracy.

Energy & Fuels↗
DOE OSTI2025

DIF3D-VARIANT 12.0: Updates and New Features

The DIF3D code has been a workhorse of fast reactor analysis work at Argonne National Laboratory for over 40 years. In 1995, a transport option called VARIANT was added to DIF3D to improve the flux solutions for fast reactor problems which we term DIF3D-VARIANT today. DIF3D-VARIANT performs nodal neutron transport calculations using P N or SP N theory in Cartesian and hexagonal two- and three-dimensional geometries. The limited computing capabilities of the time restricted DIF3D-VARIANT to use at most a 6 th order spatial approximation combined with a P3 flux approximation and P1 scattering kernel for a 33 group structure on most studied reactor problems. Computer capabilities have increased steadily since 1995 and today much larger space-angle-energy approximations are possible. This manuscript serves as an update to the theory section of the original DIF3D-VARIANT manual and details more than twenty years of changes made to DIF3D to make version 12 which was released on November 1 st , 2024. The primary focus of the initial work was to extend the space-angle approximations available in DIF3D-VARIANT such that the error due to transport approximations could be better understood. This work was started and completed in 2002 and marked the official version 10. Unfortunately, those higher order approximations could not be used at that time due to the memory constraints of the BPOINTER part of DIF3D (limited to 2 GB). In version 11, completed in 2012, BPOINTER was circumvented in DIF3D-VARIANT for the largest arrays by introducing a Fortran 90 module called LMA (Large Memory Array). This seamlessly replaces all of the functionality of the BPOINTER concept, but it allows 64 bit addressing for every array such that they can be larger than 2 GB. It is now common for DIF3D-VARIANT jobs to consume 50 GB of memory on modern workstations when using high order space-angle approximations and a large number of groups. Many improvements were made to version 11 from 2012 to 2022 when work to create version 12 started. For version 12, several parts of DIF3D were updated to improve performance and thread parallelism was introduced to further reduce the runtime. Numerous minor bugs were discovered in DIF3D-VARIANT as part of the process of creating the perturbation and sensitivity code PERSENT. All of these algorithmic problems were identified in the transition from version 10 to version 11 which prevented DIF3D-VARIANT from running efficiently and reliably. Firstly, the coarse mesh rebalance scheme would routinely diverge and a study detailed in this report demonstrates how it was also typically not effective. This is not a failure of the coarse mesh rebalance methodology, but a failure of its implementation in DIF3D-VARIANT for hexagonal geometries. The fission source extrapolation algorithm was also found to be unreliable on larger group structure problems, leading to divergence in some cases and a negligible improvement in performance overall. Finally, the “Omega” acceleration applied to the partial current solver routine of DIF3D-VARIANT was found to cause DIF3D-VARIANT to converge to the wrong answer. To resolve these issues, both the coarse mesh rebalance and fission source extrapolation were permanently disabled in version 11. The Tchebychev acceleration was put in as a temporary reliable alternative but it is generally inferior to coarse mesh rebalance or coarse mesh finite difference. For the Omega acceleration, the factor was restricted to guarantee that it would not cause follow-on errors in PERSENT. Due to limited funding to support maintenance and development of DIF3D in the last 10 years, no effort was spent since to resolve the outer iteration acceleration. Except for the threading work, all of the changes discussed in this manuscript refer to changes made between version 10 and version 11. Performance comparisons are done to demonstrate the improvements from version 9 to version 12. As will be demonstrated, the updated versi

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS2016

Precision Oxygen Isotope Measurements of Two C-Rich Hydrated Interplanetary Dust Particles

Introduction: Chondritic-smooth IDPs (Interplanetary Dust Particles) are low porosity objects whose mineralogy is dominated by aqueous alteration products such as Mg-rich phyllosilicates (smectite and serpentine group) and Mg-Fe carbonate minerals. Their hydrated mineralogy combined with low atmospheric entry velocities have been used to infer an origin largely from asteroidal sources. Spectroscopic studies show that the types and abundance of organic matter in CS IDPs is similar to that in CP IDPs. Although CS IDPs show broad similarities to primitive carbonaceous chondrites, only a few particles have been directly linked to specific meteorite groups such as CM and CI chondrites based on the presence of diagnostic minerals. Many CS IDPs however, have carbon contents that greatly exceed that of known meteorite groups suggesting that they either may derive from comets or represent samples of more primitive parent bodies than do meteorites. It is now recognized that many large, dark primitive asteroids in the outer main belt, as well as some trans-Neptunian objects, show spectroscopic evidence for aqueous alteration products on their surfaces. Some CS IDPs exhibit large bulk D enrichments similar to those observed in the cometary CP IDPs. While hydrated minerals in comets have not been unambiguously identified to date, the presence of the smectite group mineral nontronite has been inferred from infrared spectra obtained from the ejecta from comet 9P/Tempel 1 during the Deep Impact mission. Recent observations of low temperature sulfide minerals in Stardust mission samples suggest that limited aqueous activity occurred on comet Wild-2. All of these observations, taken together, suggest that the high-carbon hydrated IDPs are abundant and important samples of primitive solar system objects not represented in meteorite collections. Oxygen isotopic compositions of chondrites reflect mixing between a 16O-rich reservoir and a 17O,18O-rich reservoir produced via mass-independent fractionation. The composition of the 16O-rich reservoir is well constrained but material representing the 17O,18O-rich end-member is rare. Self-shielding models predict that cometary water, presumed to represent this reservoir, should be enriched in 17O and 18O by greater than 200 per mille. The high-carbon hydrated IDPs may be among the best materials available to search for preserved "cometary" H2O signatures. In order to better understand the origin and evolution of these particles, we have obtained 10 hydrated interplanetary dust particles for coordinated mineralogical, isotopic and organic analyses. We have previously reported the results of mineralogical and O isotopic measurements of two hydrated IDPs; here we present results of O isotopic measurements of three additional IDPs. Samples and Methods: Three interplanetary dust particles (L2079C35, L2083D46 and L2083E46) were embedded in S and partially ultramicrotomed into approximately 70 nanometer sections for analysis via transmission electron microscopy (TEM). The remainders of the unsliced particles were removed from S and pressed into high purity Au foil that was cleaned with HF acid and annealed at 800 degrees Centigrade. The pressed IDPs were analyzed via electron microprobe analysis (EPMA) for quantitative bulk chemical analysis. After EPMA analysis, the IDPs were subjected to precision O isotope analysis with the UCLA Cameca IMS-1270 ion probe. A 20 kiloelectronvolt, 0.5 nanoangstrom Cs+ primary beam of approximately 15 micrometers diameter was used for each measurement. Small particles of San Carlos olivine and Burma spinel were pressed into the Au foil for use as standards to correct for instrumental mass fractionation. The detection system was configured for multicollection, with 16O measured on a Faraday cup, and 17O and 18O measured on electron multipliers (EMs). Individual analyses consisted of 15 cycles of 10 seconds per cycle. Additionally, two microtome thin sections were measured for H isotopic compositions with the JSC NanoSIMS 50L ion probe. An 8 picoangstrom, 16 kiloelectronvolt Cs plus primary beam was used. Measurements consisted of H, C, 12C, 16O, and 18O collected with EMs in multicollection. Terrestrial biotite and kerogen were used for isotopic standards. A significant challenge in O isotope measurement of hydrated minerals is the interference from 16OH at mass 17O. We ensured that the 17O and 16OH peaks were fully resolved by using a mass resolution of greater than 7000 and by careful analyses of San Carlos olivine, Burma spinel and chlorite hydrated mineral standards. The hydride was further suppressed with a cold finger attached to an LN2 dewar to trap volatiles in the sample chamber. All sputtered ions were counted (i.e. presputtering was not used); after applying background, yield and deadtime corrections, we performed a change-point analysis on our data via R in order to determine when the sample reached sputtering equilibrium; data points collected prior to the change point were excluded. Change-point analysis was also used to determine whether the IDP had completely sputtered. Results: Mineralogy. IDPs C35 and E46 exhibited hydrated mineralogies, Fe-Ni sulfide grains, nanoglobules and occasional enstatite grains distributed throughout a fine-grained Mg-Fe saponite matrix. C35 also contained breunnerite (Mg,Fe)CO3; solar flare tracks were observed in enstatite, indicating minimal atmospheric entry heating. The mineralogy of D46 is dominated by a large FeS grain with a minor component of adhering silicate material. D46 was strongly heated during atmorpheric entry as evidenced by a well-developed magnetite rim. EPMA analyses show that both C35 and E46 have high carbon contents of 20 weight percentage (approximately 6X CI). D46 contains approximately 6 wt.weight percentage C. A significant amount of the carbon is present as carbon nanoglobules. Results: Hydrogen isotopes: Although the bulk delta D values of both sections of L2079C35 were within error of SMOW (minus 33 plus or minus 19 per mille, 1 plus or minus 14 per mille 1 sigma), several delta D-rich hotspots were also identified, reaching 2000 per mille. As shown in Fig. 1, these hotspots are clearly associated with discrete carbonaceous inclusions that are akin to nanoglobules found in many meteorites and other IDPs. Oxygen Isotopes. Results of the oxygen isotope measurements are shown in Figure 1. The oxygen isotope composition for L2079C35 was delta 18O equals plus11.6 plus or minus 1.9per mille, delta 17O equals plus 7.9 plus or minus 1.9per mille (2 standard errors). The oxygen isotope composition for L2083D46 was delta 18O equals minus 8.1 plus or minus 1.9 per mille, delta 17O equals minus 6.4 plus or minus 3.1 per mille (2 standard errors). The oxygen isotope composition for L2083E46 was delta 18O equals plus 12.0 plus or minus 1.9 per mille, delta 17O equals plus 9.2 plus or minus 2.0 per mille (2 standard errors). Discussion: Despite mineralogical similarities to highly aqueously altered carbonaceous chondrites, the hydrated IDPs we analyzed have oxygen isotopic compositions that are distinct from matrix materials in the CI, CM, and CR chondrites. The IDPs plot along the Young-Russell line, with delta 17O values for C35 and E46 suggestive of interaction with a 16O-poor reservoir. However, we have thus far not observed evidence of extreme 16O-poor reservoirs expected from self-shielding models and observed in Acfer 094 simplectite. The high carbon contents of the IDPs also set them apart from known meteoritic samples. The lack of atmospheric entry heating effects are consistent with low encounter velocities and suggest either an asteroidal source, or a low inclination, low eccentricity cometary origin. Conclusions: The unusual oxygen isotopic compositions, high carbon contents, and the abundance of Drich nanoglobules, together, suggest that the high-carbon, hydrated IDPs are derived from a primitive source that is not yet represented in meteorite collections.

Snead, C. J.↗
NASA NTRS2004

Interaction of a Boundary Layer with a Turbulent Wake

The objective of this grant was to study the transition mechanisms on a flat-plate boundary layer interacting with the wake of a bluff body. This is a simplified configuration presented and designed to exemplify the phenomena that occur in multi-element airfoils, in which the wake of an upstream element impinges on a downstream one. Some experimental data is available for this configuration at various Reynolds numbers. The first task carried out was the implementation and validation of the immersed-boundary method. This was achieved by performing calculations of the flow over a cylinder at low and moderate Reynolds numbers. The low-Reynolds number results are discussed, which is enclosed as Appendix A. The high-Reynolds number results are presented in a paper in preparation for the Journal of Fluid Mechanics. We performed calculations of the wake-boundary-layer interaction at two Reynolds numbers, Re approximately equal to 385 and 1155. The first case is discussed and a comparison of the two calculations is reported. The simulations indicate that at the lower Reynolds number the boundary layer is buffeted by the unsteady Karman vortex street shed by the cylinder. This is shown: long streaky structures appear in the boundary layer in correspondence of the three-dimensionalities in the rollers. The fluctuations, however, cannot be self-sustained due to the low Reynolds-number, and the flow does not reach a turbulent state within the computational domain. In contrast, in the higher Reynolds-number case, boundary-layer fluctuations persist after the wake has decayed (due, in part, to the higher values of the local Reynolds number Re achieved in this case); some evidence could be observed that a self-sustaining turbulence generation cycle was beginning to be established. A third simulation was subsequently carried out at a higher Reynolds number, Re=3900. This calculation gave results similar to those of the Re=l155 case. Turbulence was established at fairly low Reynolds number, as a consequence of the high level of the free-stream perturbation. An instantaneous flow visualization for that case is shown. A detailed examination of flow statistics in the transitional and turbulent regions, including the evolution of the turbulent kinetic energy (TKE) budget and frequency spectra showed the formation and evolution of turbulent spots characteristic of the bypass transition mechanism. It was also observed that the turbulent eddies achieved an equilibrium, fully developed turbulent states first, as evidenced by the early agreement achieved by the terms in the TKE budget with those observed in turbulent flows. Once a turbulent Reynolds stress profile had been established, the velocity profile began to resemble a turbulent one, first in the inner region and later in the outer region of the wall layer. An extensive comparison of the three cases, including budgets, mean velocity and Reynolds stress profiles and flow visualization, is included. The results obtained are also presented.

Piomelli, Ugo↗
NASA NTRS2017

HST WFC3 Observations of Uranus' 2014 Storm Clouds and Comparison with VLT/SINFONI and IRTF/Spex Observations

In November 2014 Uranus was observed with the Wide Field Camera 3 (WFC3) instrument of the Hubble Space Telescope as part of the Hubble 2020: Outer Planet Atmospheres Legacy program, OPAL. OPAL annually maps Jupiter, Uranus and Neptune (and will also map Saturn from 2018) in several visible near- infrared wavelength filters. The Uranus 2014 OPAL observations were made on the 89th November at a time when a huge cloud complex, first observed by de Pater et al. (2015) and subsequently tracked by professional and amateur astronomers (Sayanagi et al., 2016), was present at 30-40deg N. We imaged the entire visible atmosphere, including the storm system, in seven filters spanning 467924 nm, capturing variations in the coloration of Uranus clouds and also vertical distribution due to wavelength dependent changes in Rayleigh scattering and methane absorption optical depth. Here we analyse these new HST observations with the NEMESIS radiative-transfer and retrieval code in multiple-scattering mode to determine the vertical cloud structure in and around the storm cloud system. The same storm system was also observed in the H-band (1.4-1.8 micrometers) with the SINFONI Integral Field Unit Spectrometer on the Very Large Telescope (VLT) on 31st October and 11th November, reported by Irwin et al. (2016, 10.1016j.icarus.2015.09.010). To constrain better the cloud particle sizes and scattering properties over a wide wavelength range we also conducted a limb-darkening analysis of the background cloud structure in the 30-40deg N latitude band by simultaneously fitting: a) these HSTOPAL observations at a range of zenith angles; b) the VLTSINFONI observations at a range of zenith angles; and c) IRTFSpeX observations of this latitude band made in 2009 at a single zenith angle of 23deg, spanning the wavelength range 0.8-1.8 micrometers (Irwin et al., 2015, 10.1016j.icarus.2014.12.020). We find that the HST observations, and the combined HSTVLTIRTF observations at all locations are well modelled with a three-component cloud comprised of: 1) a vertically thin, but optically thick deep tropospheric cloud at a pressure of approximately 2 bars; 2) a methane-ice cloud based at the methane-condensation level of approximately 1.23 bar, with variable vertical extent; and 3) a vertically extended tropospheric haze, also based at the methane-condensation level of 1.23 bar. We find that modelling both haze and tropospheric cloud with particles having an effective radius of approximately 0.1 micron provides a good fit the observations, although for the tropospheric cloud, particles with an effective radius as large as 1.0 micron provide a similarly good fit. We find that the particles in both the tropospheric cloud and haze are more scattering at short wave- lengths, giving them a blue color, but are more absorbing at longer wavelengths, especially for the tropospheric haze. We find that the spectra of the storm clouds are well modelled by localized thickening and vertical extension of the methane-ice cloud. For the particles in the storm clouds, which we assume to be composed of methane ice particles, we find that their mean radii must lie somewhere in the range 0. 1 1. 0 m. We find that the high clouds have low integrated opacity, and that streamers reminiscent of convective thunderstorm anvils are confined to levels deeper than 1 bar. These results argue against vigorous moist convective origins for the cloud features.

atmospheres↗
NASA NTRS1993

Flexure and the role of inplane force around coronae on Venus

Large coronae on Venus, such as Artemis and Latona, are rimmed by conspicuous trenches and associated outer rises. Sandwell and Schubert have observed that these systems resemble terrestrial subduction zones in planform and have succeeded in fitting an elastic plate bending equation to the inferred flexural topography. However, the first zero crossing bending moments required are -2.5 x 10(exp 17) N for Artemis and -5.0 x 10(exp 16) N for Latona. Since these moments are similar in magnitude to those of subducting slabs on Earth, a rollback subduction mechanism was proposed to explain the flexure around the largest coronae, although a differential thermal subsidence model is sufficient to account for the topography around some coronae. The purpose is to investigate the effect of inplane force as a possible alternative to large applied moments in producing flexure at Artemis and Latona. The close correlation of gravity to topography on Venus implies the absence of a low viscosity zone and the strong coupling of the lithosphere to mantle convection. If coronae are the surface manifestations of mantle plumes, they may be the sites of active convective stress coupling. As the upwelling reaches the lithosphere, it spreads radially outward, inducing shear tractions on the base of the plate. In addition, the hot, expanding corona may load the surrounding plate horizonally. Both the basal shear stresses and radial loading can be treated as an equivalent compressive inplane force in the mechanical lithosphere, which contributes to the bending of the outlying plate. Using a model that relates inplane force to the measured gravity anomalies, a rough value of the inplane force at Artemis was calculated. Recent Pioneer Venus spherical harmonic gravity models indicate a geoid anomaly of about 75 m over Artemis, which corresponds to an estimated inplane force on the order of -1x10(exp 13) N/m. The gravity model is unable to resolve Latona, but an inplane force of similar dimensions is assumed. The maximum possible inplane force based on the expected rheology can be constrained by using the approximate 5 K/km thermal gradient inferred from the best fit 30 km elastic plate at Artemis and Latona. For a dry olivine flow law in the upper mantle, the compressional load limit of the 60 km thick mechanical lithosphere is -4 x 10(exp 13) N/m. This value is equivalent to a load of -8 x 10(exp 13) N/m on a 30 km thick elastic plate.

Brown, C. David↗