Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 1,063 records · Page 59

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↗

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

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

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↗

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

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↗

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↗

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↗

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↗

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

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↗

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↗

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↗

Cis-Lunar Base Camp

Historically, when mounting expeditions into uncharted territories, explorers have established strategically positioned base camps to pre-position required equipment and consumables. These base camps are secure, safe positions from which expeditions can depart when conditions are favorable, at which technology and operations can be tested and validated, and facilitate timely access to more robust facilities in the event of an emergency. For human exploration missions into deep space, cis-lunar space is well suited to serve as such a base camp. The outer regions of cis-lunar space, such as the Earth-Moon Lagrange points, lie near the edge of Earth s gravity well, allowing equipment and consumables to be aggregated with easy access to deep space and to the lunar surface, as well as more distant destinations, such as near-Earth Asteroids (NEAs) and Mars and its moons. Several approaches to utilizing a cis-lunar base camp for sustainable human exploration, as well as some possible future applications are identified. The primary objective of the analysis presented in this paper is to identify options, show the macro trends, and provide information that can be used as a basis for more detailed mission development. Compared within are the high-level performance and cost of 15 preliminary cis-lunar exploration campaigns that establish the capability to conduct crewed missions of up to one year in duration, and then aggregate mass in cis-lunar space to facilitate an expedition from Cis-Lunar Base Camp. Launch vehicles, chemical propulsion stages, and electric propulsion stages are discussed and parametric sizing values are used to create architectures of in-space transportation elements that extend the existing in-space supply chain to cis-lunar space. The transportation options to cis-lunar space assessed vary in efficiency by almost 50%; from 0.16 to 0.68 kg of cargo in cis-lunar space for every kilogram of mass in Low Earth Orbit (LEO). For the 15 cases, 5-year campaign costs vary by only 15% from 0.36 to 0.51 on a normalized scale across all campaigns. Thus the development and first flight costs of assessed transportation options are similar. However, the cost of those options per flight beyond the initial operational capability varies by 70% from 0.3 to 1.0 on a normalized scale. The 10-year campaigns assessed begin to show the effect of this large range of cost beyond initial operational capability as they vary approximately 25% with values from 0.75 to 1.0 on the normalized campaign scale. Therefore, it is important to understand both the cost of implementation and first use as well as long term utilization. Finally, minimizing long term recurring costs is critical to the affordability of future human space exploration missions. Finally minimizing long term recurring costs is critical to the affordability of future human space exploration missions.

Merrill, Raymond G.↗

22” ADP Fan Rig Liners Design Report

A liner design study was conducted as part of a cooperative effort between five (5) government/industry teams that together seek to demonstrate the technology for designing and manufacturing acoustic liners that are twenty five percent (25%) more efficient than 1992 technology liners. The study emphasized teaming collaboration. The improved liners were designed by Pratt & Whitney and Boeing Airplane Company and built by Rohr Inc, and will be tested in the NASA/P& W 22-inch ADP fan rig at NASA Lewis Research Center's 9' x 15' wind tunnel. Design guidelines and decisions were made collectively during monthly design review telecons and at the formal final design review. The tools that were used to design the new improved liners were not new, but the process that was developed as part of this study that led into the evaluation and selection of the final and best liner designs is new. Until now, the liner design process as practiced by different industry teams varied greatly. Some procedures were based on empirical liner attenuation databases which could not adequately account for engine-to-engine hardwall fan noise spectral differences, while others were entirely theoretical. And regardless of which method one chose, the major difficulty was still the lack of understanding and knowledge of the actual hardwall fan source noise modal structure and farfield SPL spectra. The NASA-led effort to attempt actual measurements of the fan source noise modal structure by means of the rotating microphone array is expected to contribute significantly to the understanding of the nature of the fan tone noise modes, but the application to broadband noise is still a long way off. The new design process included consideration of the measured fan tone modes, but for the majority of the spectra a separate, systematic process was used. It is a common practice in the engine/nacelle industry that acoustic liners for new products are designed before actual measured hardwall engine far-field noise spectra are available. Target noise spectra are normally derived from existing engine noise databases with some adjustments to absolute SPLs and frequencies. This practice has been acceptable as long as the new engine was a derivative of the base engine. However, in the case of the ADP, the transition from a current engine base is too - great, and the adjustments could not adequately account for the quantum changes in SPL spectral differences. Examples were the 1992 single degree of freedom (SDOF) inlet and aft liners (designed by P&W) that would be used as the baseline liners against which the new improved liners' efficien¬cies would -be measured. As will be shown, these baseline liners have been found to be deeper than the desire optimum depths. The new liner design process begins with the selection of the target hard wall fan noise spectra. These spectra were obtained by scaling up (5.91 scale factor) the measured hardwall fan spectral data from the 22” ADP fan rig. Next, the modal energy contents for each 1/3-octave band center frequency of these hardwall fan noise spectra are estimated. (Within each 1/3-octave band center frequency, the model energy distribution approximation is for both tone and broadband noise). For the inlet noise, P&W uses a derivative of the NASA Lewis (Ed Rice's) method of classifying; and grouping propagating modes by their cutoff ratios. The next step is to assign energy level to each group of modes having the same cutoff ratio. In Ed Rice's model, the modes were grouped according to ten (10) cutoff ratio intervals with center values located at 1.026, 1.085, 1.155, 1.24, 1.35, 1.49, 1.69, 2.0, and 4.47 (with equal number of modes in each interval). The modes were assumed to have equal energy. P&W’s model expands the cutoff ratio-mode grouping into two hundred (200) smaller cutoff ratio intervals with center values located f1".'m 1.003 to 11.5 in 199 increasing incremental intervals. Next, P&W's model uses a "2-parameter" normal distribution as a template to assign energy levels to these 200 pre-determined cutoff ratio values (for each 1/3-octave band center frequency). In the past, P&W had conducted an extensive study to determine what "2-parameter" values are appropriate for fullscale inlet liners, and had developed a set of twenty-four (24) "2-pararneter" values (i.e. one for each 1/3-octave frequency band) that when used in Ed Rice's Inlet Attenuation Prediction Method produced predicted liner attenuations that closely matched measured liner attenuations from several P&W's engines. In the absence of actual measured tone and broadband modal data from the 22" ADP rig, the process will use P&W's proprietary set of "2-parameter" values for this liner design study. For the aft noise, Boeing uses a modal energy approximation that the "transport energy" of each propagating mode is equal. This approximation is almost the same as the equal energy per mode approximation, except for modes that are near cutoff. Boeing's model forces these modes to have lower energy levels. Both P&W and Boeing agreed that the "transport energy" approximation should work well for the ADP aft fan noise which appears to be dominated by broadband noise. The design process then proceeds to calculate the optimum liner impedances for each frequency in both the inlet and the aft. These optimum impedances represent the target impedances that the designed liners should have. Liners with impedances matching the optimum impedances at all frequencies are "ideal" liners. These ideal liners are theoretically the best liners. Unfortunately, it has been showed that it is impossible to design and build such ideal liners. The next best liners are ones that have impedances matching the optimum impedances over some frequency range (not all frequencies as for the ideal). This is accomplished by the use of "frequency weightings". Several of P&W's and Boeing's existing computer decks were used for optimizing and matching the designed liner impedances to the target optimum values (with the various frequency weightings specified). The optimization produces liner candidates with predicted liner impedances and descriptions of their physical liner characteristics. These candidate liner impedances are then used to predict their spectral attenuation characteristics which are then used together with the target hardwall fan noise spectra to determine the resulting treated noise spectra and PNLT values. Further optimization around the selected candidate designs yield final designs that are best in PNLT attenuations. Use of the optimization decks allowed a large number of liner candidates to be screened in a relatively short period of time. This design process is systematic and is efficient. The new inlet SDOF liner design obtained from the improved process was predicted to be 34% more effective (per unit area) than the 1992 baseline inlet liner. This inlet liner design is a 112 rayl wovenwiremesh facesheet over a 0.312-inch deep honeycomb core. The new aft SDOF liner design was predicted to be 52% more efficient than the aft baseline liner. The new aft SDOF liner is "segmented" with a shallower liner on the core cowl (inner duct wall) and a deeper liner on the fan cowl (outer duct wall). The shallower liner is a 70.6 rayl woven-wiremesh facesheet over a 0.141-inch deep honeycomb core. The deeper liner is a 68 rayl woven-wiremesh facesheet over a 0.309-inch deep honeycomb core. All liner dimensions are for the 22-inch model-scale ADP fan rig liners. The selected advanced liners are "segmented" double-layer (DDOF) liners for the inlet and aft locations. Also, for the inlet, a bulk liner with ceramic foam for wider broadband noise absorption was also selected. Triple-layer liner designs were not considered since the model scaled liners ( 1/5. 91) were dimensionally too small to be built correctly, and irrin earlier concept study, Boeing found a triple-layer to have only very small benefits over a double-layer. The inlet DDOF liner was predicted to be 83% more effective than the baseline. The inlet DDOF design is a 78 rayl facesheet over a 0.080 top cavity depth, a 68 rayl septum and a 0.227-inch bottom cavity depth. The inlet bulk liner is a 60 rayl facesheet over a 0.33-inch deep honeycomb filled with high temperature (HTP) ceramic foam with a density of 4.8 lb/cu.ft and a flow resistivity of 167 rayl/cm. The inlet bulk liner was predicted to be 83% more effective than the baseline liner. The aft DDOF liners are segmented with a shallower liner on the core cowl and a deeper liner on the fan cowl. The shallow DDOF liner is a 49.8 rayl facesheet over a 0.093-inch top cavity depth, a septum of 88.2 rayls over a 0.181-inch bottom cavity depth. The deep DDOF liner is a 12.9 rayl facesheet over a 0.140-inch top cavity depth, a septum of 53.1 ray ls over a 0.258-inch botom cavity depth. The segmented aft DDOF liners were predicted to be 86% more efficient than the baseline liner.

turbofan acoustic treatment↗

Multifrequency studies of bright radio supernova remnants. 2: W49B

We report radio observations of the supernova remnant W49B using the Very Large Array (VLA) of the NRAO at 90, 20, and 6 cm. This work continues the study of the properties of young, bright supernova remnants (SNRs) begun with 3C 391 (Moffett & Reynolds 1994). Here we present high-resolution images of total intensity, polarization, and spectral index of W49B. In large-scale morphology it is basically a box-shaped remnant expanding into an apparently inhomogeneous medium, but we also find small-scale structures, arcs and filaments, which lie interior to the outer shell in projection. The actual spatial position of the filaments is unknown, though we suspect that they are in the remnant interior rather than on the front or back face of the blast wave. In addition, their distribution in two dimensions suggests the projection of a helical structure. The shell edge, sharply rising in brightness for at least 3/4 of its circumference, still remains unresolved at our highest resolution (4 sec). This for our assumed distance of 8 kpc, the width of the region in which the emission appears is less than 0.16 pc, indicating short mean free paths for shock-accelerated electrons and high levels of MHD turbulence presumably causing the scattering. We find no polarized flux at 90 or 20 cm, with 3 sigma upper limits in polarized intensity of 36 mJy at 90 cm (45 sec resolution) and 7.2 mJy at 20 cm (5 sec resolution), or 22 micro-Jy/sq arcsec and 370 micro-jy/sq arcsec, respectively, for any emission extended on those scales. Polarized flux is present at 6 cm, but at a very low mean polarized fraction (total polarized flux divided by total flux) of 0.44% +/- 0.06%, which, as for 3C 391, is much lower than typical for bright SNRs at this frequency. The morphology in polarized intensity is poorly correlated with that in total intensity. We see excursions in the polarized fraction up to at least 10% in a few locations, but even where polarization is seen, typical levels are a few percent. Tangled or disordered magnetic fields in the emitting region of the radio shell may be responsible for depolarizing the radio synchrotron radiation, but some form of internal Faraday depolarization may also occur. We estimate the foreground Faraday rotation measure to be about -450 rad/sq m, similar to that found for 3C 391, which is also the Galactic plane and just about as distant. Spectral index images created from the total intensity images show that the spectral index across W49B is constant to within about Delta alpha approximately 0.1 in bright regions. This result weakly supports a common origin of the radio-emitting electrons, as in the blast wave, rather than in inhomogeneous turbulent regions of differing properties due to the stochastic (second-order Fermi) acceleration process. Variations at the level of Delta alpha approximately 0.1 are seen, but their significance is doubtful. New observations at 90 cm, using experimental three-dimensional imaging technology, may improve on this limit.

Moffett, David A.↗

Z-Pinch Pulsed Plasma Propulsion Technology Development

Fusion-based propulsion can enable fast interplanetary transportation. Magneto-inertial fusion (MIF) is an approach which has been shown to potentially lead to a low cost, small reactor for fusion break even. The Z-Pinch/dense plasma focus method is an MIF concept in which a column of gas is compressed to thermonuclear conditions by an axial current (I approximates 100 MA). Recent advancements in experiments and the theoretical understanding of this concept suggest favorable scaling of fusion power output yield as I(sup 4). This document presents a conceptual design of a Z-Pinch fusion propulsion system and a vehicle for human exploration. The purpose of this study is to apply Z-Pinch fusion principles to the design of a propulsion system for an interplanetary spacecraft. This study took four steps in service of that objective; these steps are identified below. 1. Z-Pinch Modeling and Analysis: There is a wealth of literature characterizing Z-Pinch physics and existing Z-Pinch physics models. In order to be useful in engineering analysis, simplified Z-Pinch fusion thermodynamic models are required to give propulsion engineers the quantity of plasma, plasma temperature, rate of expansion, etc. The study team developed these models in this study. 2. Propulsion Modeling and Analysis: While the Z-Pinch models characterize the fusion process itself, propulsion models calculate the parameters that characterize the propulsion system (thrust, specific impulse, etc.) The study team developed a Z-Pinch propulsion model and used it to determine the best values for pulse rate, amount of propellant per pulse, and mixture ratio of the D-T and liner materials as well as the resulting thrust and specific impulse of the system. 3. Mission Analysis: Several potential missions were studied. Trajectory analysis using data from the propulsion model was used to determine the duration of the propulsion burns, the amount of propellant expended to complete each mission considered. 4. Vehicle Design: To understand the applicability of Z-Pinch propulsion to interplanetary travel, it is necessary to design a concept vehicle that uses it -- the propulsion system significantly impacts the design of the electrical, thermal control, avionics and structural subsystems of a vehicle. The study team developed a conceptual design of an interplanetary vehicle that transports crew and cargo to Mars and back and can be reused for other missions. Several aspects of this vehicle are based on a previous crewed fusion vehicle study -- the Human Outer Planet Exploration (HOPE) Magnetized Target Fusion (MTF) vehicle. Portions of the vehicle design were used outright and others were modified from the MTF design in order to maintain comparability.

Polsgrove, Tara↗

Organic Nano-Grains in Comet 103P/Hartley 2: The Organic Glue of Porous Aggregate Grains?

The GNIRS instrument on the Gemini 8-­‐m telescope observed comet 103P/Hartley on 2010-­‐ Dec-­‐04UT, a month after the EPOXI Mission encounter, and detected the 3.3 and 3.4 um bands in emission. The 3.3/3.4 ratio and the broad band widths are consistent with experiments of heated (approximately 600 K) aliphatic carbon (-CH3, -CH2) thin films. For the 3.4 micron band to be in emission, the aliphatic bonds must be attached to a carrier possessing the strongly UV-­‐absorbing C=C aromatic rings, and these rings have to be less than 50-­‐100 carbon atoms (4-6 Angstrom) for attached -CH bonds to also generate a 3.3 micron-band in emission. Slightly larger (≥10Å) Very Small Grains (VSGs) can absorb single UV photons comparable to or exceeding their heat capacity, thermally fluctuate and release IR photon(s). The 3.3 micron and 3.4 micron bands observed by GNIRS suggest that organic macromolecules/ nano-­‐grains with both aliphatic and aromatic bonds are fluorescing/thermally fluctuating in the coma. Aliphatic and aromatic materials have been seen in Stardust samples and the primitive carbonaceous chondrite 'Tagish Lake'. The larger the ratio of the -CH2/-CH3 components of the aliphatic 3.4 micron band, the more 'primitive' the organic material. In a Stardust organic globule, some aliphatic bonds were transformed into aromatic bonds during the low dosage of Transmission Electron Microscope imaging. Conversely, lab experiments show irradiation of ices containing small PAHs generates aliphatic organics. Photo-­‐processing of ices also likely forms the ubiquitous aliphatic coatings that appear on the surfaces of all silicate subgrains constituting nine cometary interplanetary dust particles. The aliphatic coatings, dominated by -CH2, likely were important in sticking the aggregates together, and existed prior to incorporation of dust aggregates into comet nuclei. These comet aliphatics may be some of the sought-­‐after precursors to the more robust and complex organics studied as Insoluble Organic Matter in carbonaceous chondrites. Aliphatic coatings on submicron grains, however, will not be observable in absorption because they are fairly transparent, nor do the aliphatic carbonaceous coatings produce the 3.4 micron emission band because the particles they are attached to are too large (too many vibration modes). We must probe the nano-­‐sized organic carriers that undergo substantive thermal fluctuations in cometary comae and emit at 3.3 3.4 micron. Observations of the 3.3 and 3.4 micron emission features contribute to characterizing the evolution of organics prior to their incorporation into cometary nuclei as well as their rapid evolution in cometary comae, which in turn contributes to deepening our understanding of the evolution of organics on the surfaces of asteroids and outer icy bodies in our solar system. Studying organics in comets contributes to understanding the formation and evolution pathways of ISM organics through to the formation of the robust insoluble organic matter in meteorites. A'Hearn, M.F., et al. 2011, Science, 332, 1396; Bockelee-­‐Morvan, D. et al. 1995, Icarus, 116, 18; De Gregorio, B.T., et al. 2010, GCA, 74, 4454; Dello Russo, N., et al. 2011, ApJ, 734, L8; Dischler et al. 1983, Solid State Communications, 48, 105; Flynn, G., et al. 2010a, LPSC, 41, #1079; Flynn, G., et al. 2010b, COSPAR, 38, F31-­‐0012-­‐10; Flynn, G., Wirick, S. 2011, LPSC, 42, #1856; Fomenkova, et al. 1994, GCA 58, 4503; Matrajt, G., et al. 2013, ApJ, 765, 145; Schutte, et al. 1993, ApJ, 415, 397; Wooden, D.H. et al. 2011, EPSC-­‐DPS, 1557; Wooden, D.H. et al. 2013, submitted.

Stardust samples↗