Search NASA⌕ Search

SEARCH · Search NASA

Results for “approximate dynamic programming”

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 505 records · Page 28

Towards a unified picture of the water self-ions at the air-water interface: a density functional theory perspective

The propensities of the water self ions, H3O+ and OH- , for the air-water interface has implications for interfacial acid-base chemistry. Despite numerous experimental and computational studies, no consensus has been reached on the question of whether or not H3O+ and/or OH- prefer to be at the water surface or in the bulk. Here we report a molecular dynamics simulation study of the bulk vs. interfacial behavior of H3O+ and OH- that employs forces derived from density functional theory with a generalized gradient approximation exchangecorrelation functional (specifically, BLYP) and empirical dispersion corrections. We computed the potential of mean force (PMF) for H3O+ as a function of the position of the ion in a 215-molecule water slab. The PMF is flat, suggesting that H3O+ has equal propensity for the air-water interface and the bulk. We compare the PMF for H3O+ to our previously computed PMF for OH- adsorption, which contains a shallow minimum at the interface, and we explore how differences in solvation of each ion at the interface vs. the bulk are connected with interfacial propensity. We find that the solvation shell of H3O+ is only slightly dependent on its position in the water slab, while OH- partially desolvates as it approaches the interface, and we examine how this difference in solvation behavior is manifested in the electronic structure and chemistry of the two ions. DJT was supported by National Science Foundation grant CHE-0909227. CJM was supported by the U.S. Department of Energy‘s (DOE) Office of Basic Energy Sciences, Division of Chemical Sciences, Geosciences and Biosciences. Pacific Northwest National Laboratory (PNNL) is operated for the Department of Energy by Battelle. The potential of mean force required resources of the Oak Ridge Leadership Computing Facility at the Oak Ridge National Laboratory, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DEAC05-00OR22725. The remaining simulations and analysis used resources of the National Energy Research Scientific Computing Center, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC02-05CH11231. at at Lawrence Berkeley National Laboratory. MDB is grateful for the support of the Linus Pauling Distinguished Postdoctoral Fellowship Program at PNNL.

Baer, Marcel D.↗

Conservation of Endo16 expression in sea urchins despite evolutionary divergence in both cis and trans-acting components of transcriptional regulation

Evolutionary changes in transcriptional regulation undoubtedly play an important role in creating morphological diversity. However, there is little information about the evolutionary dynamics of cis-regulatory sequences. This study examines the functional consequence of evolutionary changes in the Endo16 promoter of sea urchins. The Endo16 gene encodes a large extracellular protein that is expressed in the endoderm and may play a role in cell adhesion. Its promoter has been characterized in exceptional detail in the purple sea urchin, Strongylocentrotus purpuratus. We have characterized the structure and function of the Endo16 promoter from a second sea urchin species, Lytechinus variegatus. The Endo16 promoter sequences have evolved in a strongly mosaic manner since these species diverged approximately 35 million years ago: the most proximal region (module A) is conserved, but the remaining modules (B-G) are unalignable. Despite extensive divergence in promoter sequences, the pattern of Endo16 transcription is largely conserved during embryonic and larval development. Transient expression assays demonstrate that 2.2 kb of upstream sequence in either species is sufficient to drive GFP reporter expression that correctly mimics this pattern of Endo16 transcription. Reciprocal cross-species transient expression assays imply that changes have also evolved in the set of transcription factors that interact with the Endo16 promoter. Taken together, these results suggest that stabilizing selection on the transcriptional output may have operated to maintain a similar pattern of Endo16 expression in S. purpuratus and L. variegatus, despite dramatic divergence in promoter sequence and mechanisms of transcriptional regulation.

NASA Discipline Evolutionary Biology↗

Charge Recombination, Transport Dynamics, and Interfacial Effects in Organic Solar Cells

The need for renewable sources of energy is well known. Conversion of sunlight to electricity using solar cells is one of the most important opportunities for creating renewable energy sources. The research carried out under DE-FG02-08ER46535 focused on the science and technology of “Plastic” solar cells comprised of organic (i.e. carbon based) semiconductors. The Bulk Heterojunction concept involves a phase separated blend of two organic semiconductors each with dimensions in the nano-meter length scale --- one a material that functions as a donor for electrons and the other a material that functions as an acceptor for electrons. The nano-scale inter-penetrating network concept for “Plastic” solar cells was created at UC Santa Barbara. A simple measure of the impact of this concept can be obtained from a Google search which gives 244,000 “hits” for the Bulk Heterojunction solar cell. Research funded through this program focused on four major areas: 1. Interfacial effects in organic photovoltaics, 2. Charge transfer and photogeneration of mobile charge carriers in organic photovoltaics, 3. Transport and recombination of the photogenerated charge carriers in organic photovoltaics, 4. Synthesis of novel organic semiconducting polymers and semiconducting small molecules, including conjugated polyelectrolytes. Following the discovery of ultrafast charge transfer at UC Santa Barbara in 1992, the nano-organic (Bulk Heterojunction) concept was formulated. The need for a morphology comprising two interpenetrating bicontinuous networks was clear: one network to carry the photogenerated electrons (negative charge) to the cathode and one network to carry the photo-generated holes (positive charge) to the anode. This remarkable self-assembled network morphology has now been established using Transmission electron Microscopy (TEM) either in the Phase Contrast mode or via TEM-Tomography. The steps involved in delivering power from a solar cell to an external circuit are the following: • Photo-excitation of the donor (or the acceptor). • Charge transfer with holes in the donor domain and electrons in the acceptor domain. • Sweep-out to electrodes prior to recombination by the internal electric field. • Energy delivered to the external circuit. Each of these four steps was studied in detail using a wide variety of organic semiconductors with different molecular structures. This UC Santa Barbara group was the first to clarify the origin and the mechanism involved in the ultrafast charge transfer process. The ultrafast charge transfer (time scale approximately 100 times faster than the first step in the photo-synthesis of green plants) is the fundamental reason for the potential for high power conversion efficiency of sunlight to electricity from plastic solar cells. The UCSB group was the first to emphasize, clarify and demonstrate the need for sweep-out to electrodes prior to recombination by the internal electric field. The UCSB group was the first to synthesize small molecule organic semiconductors capable of high power conversion efficiencies. The results of this research were published in high impact peer-reviewed journals. Our published papers (40 in number) provide answers to fundamental questions that have been heavily discussed and debated in the field of Bulk Heterojunction Solar Cells; scientific questions that must be resolved before this technology can be ready for commercialization in large scale for production of renewable energy. Of the forty publications listed, nineteen were co-authored by two or more of the PIs, consistent with the multi-investigator approach described in the original proposal. The specific advantages of this “plastic” solar cell technology are the following: a. Manufacturing by low-cost printing technology using soluble organic semiconductors; this approach can be implemented in large scale by roll-to-roll printing on plastic substrates. b. Low energy cost in manufacturing; all steps carried out at room temperature (approx. a factor of ten less than the use of Silicon which requires high temperature processing). c. Low carbon footprint d. Lightweight, flexible and rugged Because of the resolution of many scientific issues, a significant fraction of which were addressed in the research results of DE-FG02-08ER46535, the power conversion efficiencies are improving at an ever increasing rate. During the funding period of DE-FG02-08ER46535, the power conversion efficiencies of plastic solar cells improved from just a few per cent to values greater than 11% with contributions from our group and from researchers all over the world.

14 SOLAR ENERGY↗

Planetary Atmosphere Dynamics and Radiative Transfer

This research program has dealt with two projects in the field of planetary atmosphere dynamics and radiative energy transfer, one theoretical and one experimental. The first project, in radiative energy transfer, incorporated the capability to isolate and quantify the contribution of individual atmospheric components to the Venus radiative balance and thermal structure to greatly improve the current understanding of the radiative processes occurring within the Venus atmosphere. This is possible by varying the mixing ratios of each gas species, and the location, number density and aerosol size distributions of the clouds. This project was a continuation of the work initiated under a 1992 University Consortium Agreement. Under the just completed grant, work has continued on the use of a convolution-based algorithm that provided the capability to calculate the k coefficients of a gas mixture at different temperatures, pressures and spectral intervals from the separate k-distributions of the individual gas species. The second primary goal of this research dealt with the Doppler wind retrieval for the Successful Galileo Jupiter probe mission in December, 1995. In anticipation of the arrival of Galileo at Jupiter, software development continued to read the radioscience and probe/orbiter trajectory data provided by the Galileo project and required for Jupiter zonal wind measurements. Sample experiment radioscience data records and probe/orbiter trajectory data files provided by the Galileo Radioscience and Navigation teams at the Jet Propulsion Laboratory, respectively, were used for the first phase of the software development. The software to read the necessary data records was completed in 1995. The procedure by which the wind retrieval takes place begins with initial consistency checks of the raw data, preliminary data reductions, wind recoveries, iterative reconstruction of the probe descent profile, and refined wind recoveries. At each stage of the wind recovery consistency is checked and maintained between the orbiter navigational data, the radioscience data, and the probe descent profile derived by the Atmospheric Instrument Team. Preliminary results show that the zonal winds at Jupiter increase with depth to approximately 150 m/s.

Atkinson, David H.↗

Moist Process Biases in Simulations of the Madden–Julian Oscillation Episodes Observed during the AMIE/DYNAMO Field Campaign

Two Madden–Julian oscillation (MJO) episodes observed during the 2011 Atmospheric Radiation Measurement Program MJO Investigation Experiment (AMIE)/DYNAMO field campaign are simulated using a regional model with various cumulus parameterizations, a regional cloud-permitting model, and a global variable-resolution model with a high-resolution region centered over the tropical Indian Ocean. Model biases in relationships relevant to existing instability theories of MJO are examined and their relative contributions to the overall model errors are quantified using a linear statistical model. The model simulations capture the observed approximately log-linear relationship between moisture saturation fraction and precipitation, but precipitation associated with the given saturation fraction is overestimated especially at low saturation fraction values. This bias is a major contributor to the excessive precipitation during the suppressed phase of MJO. After accounting for this bias using a linear statistical model, the spatial and temporal structures of the model-simulated MJO episodes are much improved, and what remains of the biases is strongly correlated with biases in saturation fraction. The excess precipitation bias during the suppressed phase of the MJO episodes is accompanied by excessive column-integrated radiative forcing and surface evaporation. A large portion of the bias in evaporation is related to biases in wind speed, which are correlated with those of precipitation. Furthermore, these findings suggest that the precipitation bias sustains itself at least partly by cloud radiative feedbacks and convection–surface wind interactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Dipolar Fermi gases in anisotropic traps

The quest for quantum degenerate Fermi gases interacting through the anisotropic and long-range dipole-dipole interaction is an exciting and fast developing branch within the cold-atom research program. Recent experimental progress in trapping, cooling, and controlling polar molecules with large electric dipole moments has, therefore, motivated much theoretical effort. In a recent letter, we have briefly discussed the application of a variational time-dependent Hartree-Fock approach to study theoretically both the static and the dynamic properties of such a system in a cylinder-symmetric harmonic trap. We focused on the hydrodynamic regime, where collisions ensure the equilibrium locally. Here we present a detailed theory extended to encompass the general case of a harmonic trap geometry without any symmetry. After deriving the equations of motion for the gas, we explore their static solutions to investigate key properties like the aspect ratios in both real and momentum space as well as the stability diagram. We find that, despite the lack of symmetry of the trap, the momentum distribution remains cylinder symmetric. The equations of motion are then used to study the low-lying hydrodynamic excitations, where, apart from the quadrupole and monopole modes, the radial quadrupole mode is also investigated. Furthermore, we study the time-of-flight dynamics as it represents an important diagnostic tool for quantum gases. We find that the real-space aspect ratios are inverted during the expansion, while that in momentum space becomes asymptotically unity. In addition, anisotropic features of the dipole-dipole interaction are discussed in detail. These results could be particularly useful for future investigations of strongly dipolar heteronuclear polar molecules deep in the quantum degenerate regime.

74 ATOMIC AND MOLECULAR PHYSICS↗

NPP Simulators for Coupled Thermal and Electric Power Dispatch

The Light Water Reactor Sustainability (LWRS) program within the United States Department of Energy supports extending the operation of the U.S. commercial nuclear power plant (NPP) fleet. Within the LWRS program, the Flexible Plant Operation and Generation (FPOG) Pathway works to diversify the revenue streams of light water reactors (LWRs) by opening opportunities for the co-generation of non-electric products in addition to supplying electrical power to the grid. Recent events have added greater motivation to these efforts. For example, the recent Inflation Reduction Act (IRA) passed by the U.S. federal government offers substantial tax incentives for producing clean hydrogen, the technology readiness level of dispatchable and high-efficiency hydrogen production has dramatically increased in a short time, and societal response to world climate change is driving a transition away from fossil fuels. Producing hydrogen with maximum efficiency using nuclear power requires dispatching both electrical and thermal power from the nuclear plant to the hydrogen plant, so testing concepts of operations for combined electrical and thermal power dispatch (TPD) from an NNP to a hydrogen plant is of interest. This report documents achievement of the Light Water Reactor Sustainability (LWRS) program milestone “Install and demonstrate a vendor-developed simulator on the Human Systems Simulation Laboratory (HSS) for dispatch of LWR electrical power to a close-coupled electrolysis plant” with a due date of Dec. 22, 2022. Several factors provide motivation for this effort. Coupling the power generation deck of a nuclear power plant to a hydrogen production facility introduces new possibilities for operational transients that must be addressed. In particular, the performance of the integrated system during startup and shutdown of the hydrogen production facility, as well as offnormal conditions, need to be evaluated to ensure there are no adverse effects on the operation of the existing NPP. The concept of operations involving the NPP, the hydrogen plant, and the electric power grid must be tested using NPP simulators and operating procedures that have been modified for TPD operations. These tests must also include dynamic simulations of the coupled tertiary thermal and electric loads as well as coordinated activities with NPP operators, tertiary load operators and grid power coordinators. The report summarizes progress in developing and testing full-scope NPP simulators at the HSSL, including a generic BWR simulator from GSE Systems, Inc. and generic PWR simulator from Westinghouse. In the case of the TPD-GBWR Simulator from GSE Systems, Inc., a BWR is thermally coupled to a high temperature electrolysis (HTE) plant that produces hydrogen and oxygen from de-ionized water. The hydrogen plant is not explicitly simulated but only included as a transient heat sink. A thermal power dispatch (TPD) system transfers heat between the steam systems at the BWR and the hydrogen plant. Operational results from two versions of the modified simulator are presented. The first version uses synthetic oil as a heat transfer fluid in a closed delivery heat loop (DHL) that generates steam at the hydrogen plant. The second version uses steam as the heat transfer fluid in a delivery steam line (DSL) to provide steam to the hydrogen plant. For both versions, the estimated thermal power delivery distance is approximately one kilometer. The amount of thermal power dispatched in the simulators is 15% of the total reactor thermal power such that the simulators provide a tool to study the feasibility of coupling a BWR to industrial processes that benefit from a combination electrical and thermal power dispatch. Ongoing work within a CRADA is also developing a full-scope PWR simulator provided by Westinghouse for both thermal and electric power coupling. This simulator is based on a PWR plant with two three-loop Westinghouse reactors. Westinghouse PWRs are sufficiently similar that a simulator of a three-loop reactor is an appropriate representation for two-loop and four-loop PWR reactors. The three-loop simulator will initially be modified for close-coupling to a 100 MW HTE hydrogen production plant that will require approximately 25 MW of thermal power while operating at its maximum rated capacity. The simulator testing will include full coupling to dynamic simulations of a hydrogen production plant and a representative bulk electric grid. The simulator provided by Westinghouse is similar to the GPWR simulator that INL has already obtained from GSE Systems but has a few important added benefits. First, the Westinghouse simulator is based on digital controls and has additional screens that can be called up to show parameter trends to assist operators in decision-making. The Westinghouse simulator also has upgrades to the controls and hardware representations, such as valve actuators, that make it more realistic and flexible in terms of accurately sim

99 GENERAL AND MISCELLANEOUS↗

Finite Element Models for Electron Beam Freeform Fabrication Process

Electron beam freeform fabrication (EBF3) is a member of an emerging class of direct manufacturing processes known as solid freeform fabrication (SFF); another member of the class is the laser deposition process. Successful application of the EBF3 process requires precise control of a number of process parameters such as the EB power, speed, and metal feed rate in order to ensure thermal management; good fusion between the substrate and the first layer and between successive layers; minimize part distortion and residual stresses; and control the microstructure of the finished product. This is the only effort thus far that has addressed computer simulation of the EBF3 process. The models developed in this effort can assist in reducing the number of trials in the laboratory or on the shop floor while making high-quality parts. With some modifications, their use can be further extended to the simulation of laser, TIG (tungsten inert gas), and other deposition processes. A solid mechanics-based finite element code, ABAQUS, was chosen as the primary engine in developing these models whereas a computational fluid dynamics (CFD) code, Fluent, was used in a support role. Several innovative concepts were developed, some of which are highlighted below. These concepts were implemented in a number of new computer models either in the form of stand-alone programs or as user subroutines for ABAQUS and Fluent codes. A database of thermo-physical, mechanical, fluid, and metallurgical properties of stainless steel 304 was developed. Computing models for Gaussian and raster modes of the electron beam heat input were developed. Also, new schemes were devised to account for the heat sink effect during the deposition process. These innovations, and others, lead to improved models for thermal management and prediction of transient/residual stresses and distortions. Two approaches for the prediction of microstructure were pursued. The first was an empirical approach involving the computation of thermal gradient, solidification rate, and velocity (G,R,V) coupled with the use of a solidification map that should be known a priori. The second approach relies completely on computer simulation. For this purpose a criterion for the prediction of morphology was proposed, which was combined with three alternative models for the prediction of microstructure; one based on solidification kinetics, the second on phase diagram, and the third on differential scanning calorimetry data. The last was found to be the simplest and the most versatile; it can be used with multicomponent alloys and rapid solidification without any additional difficulty. For the purpose of (limited) experimental validation, finite element models developed in this effort were applied to three different shapes made of stainless steel 304 material, designed expressly for this effort with an increasing level of complexity. These finite element models require large computation time, especially when applied to deposits with multiple adjacent beads and layers. This problem can be overcome, to some extent, by the use of fast, multi-core computers. Also, due to their numerical nature coupled with the fact that solid mechanics- based models are being used to represent the material behavior in liquid and vapor phases as well, the models have some inherent approximations that become more pronounced when dealing with multi-bead and multi-layer deposits.

Chandra, Umesh↗

Magnetic Untwisting in Most Solar X-Ray Jets

From 54 X‐ray jets observed in the polar coronal holes by Hinode's X‐Ray Telescope (XRT) during coverage in movies from Solar Dynamic Observatory's Atmospheric Imaging Assembly (AIA) taken in its He II 304 Å band at a cadence of 12 s, we have established a basic characteristic of solar X‐ray jets: untwisting motion in the spire. In this presentation, we show the progression of few of these X‐ray jets in XRT images and track their untwisting in AIA He II images. From their structure displayed in their XRT movies, 19 jets were evidently standard jets made by interchange reconnection of the magnetic‐arcade base with ambient open field, 32 were evidently blowout jets made by blowout eruption of the base arcade, and 3 were of ambiguous form. As was anticipated from the >10,000 km span of the base arcade in most polar X‐ray jets and from the disparity of standard jets and blowout jets in their magnetic production, few of the standard X‐ray jets (3 of 19) but nearly all of the blowout X‐ray jets (29 of 32) carried enough cool (T is approximately 105 K) plasma to be seen in their He II movies. In the 32 X‐ray jets that showed a cool component, the He II movies show 10‐100 km/s untwisting motions about the axis of the spire in all 3 standard jets and in 26 of the 29 blowout jets. Evidently, the open magnetic field in nearly all blowout X‐ray jets and probably in most standard X‐ray jets carries transient twist. This twist apparently relaxes by propagating out along the open field as a torsional wave. High‐resolution spectrograms and Dopplergrams have shown that most Type‐II spicules have torsional motions of 10‐30 km/s. Our observation of similar torsional motion in X‐ray jets strengthens the case for Type‐II spicules being made in the same way as X‐ray jets, by blowout eruption of a twisted magnetic arcade in the spicule base and/or by interchange reconnection of the twisted base arcade with the ambient open field. This work was funded by NASA's Heliophysics Division through its LWS TRT Program and its Hinode Project, and by NSF through its Research Experience for Undergraduates Program.

Moore, Ronald↗

The Jefferson Lab High Power Broadband THz Facility

Jefferson Lab operates a broadband terahertz (THz) facility, whose light is generated by 300 fs FWHM 135 pC bunches of electrons at an energy of 125 MeV as they traverse a magnetic compression chicane prior to entering the optical cavity of the Jefferson Lab Free Electron Laser (FEL)[1-3]. The repetition rate of the electron bunches, and hence that of the light pulses, can be as high as 75 MHz continuous. Approximately 1 microJoule of energy is emitted for each pulse, which implies an average power of up to 75W, a peak power of several MW and electric fields of MV/m. Niche applications of this unique source are pump-probe dynamics, both in and out of equilibrium, and large-scale imaging. The beamline itself is designed to deliver an f/8 beam 2mm vertical × 3mm horizontal in size and 60% vertically polarized from a subtended angle of ~200 milliradians horizontal by ~135 milliradians vertical in the accelerator. The first mirror is copper-braid cooled. The philosophy of the optical design is to use a relay optics configuration to transport the beam via reflection off metal mirrors through a series of focal points, with the final one being at a diamond window. The optical scheme was developed using the Synchrotron Radiation Workshop (SRW)[4]. This code performs a full calculation of the electric field from a relativistic electron. It does not handle multiparticle coherent enhancement[5], but it is significantly different from all the other synchrotron radiation calculations because it retains a term called the Coulomb term in which the electric field falls off as 1/R2. In addition to calculating the electric field and the intensity, the code also allows one to propagate the field through the optical focusing elements and apertures. The beam exiting the final diamond window in the THz user lab passes into a series of vacuum chambers with considerable flexibility for attenuation, collimation, and/or beam splitting for imaging, spectroscopy and 2 photon pump-probe measurements. In this talk we will discuss the details of the philosophy, optical design and operating experience of this facility and also present ideas for discussion. We deeply appreciate the excellent support of our colleagues without which these facilities would not have been built and operated. This work was supported by the U.S. Dept. of Energy under contract DE-AC05-06-OR23177, and the Commonwealth of Virginia. The FEL facility is supported by the Office of Naval Research. 1. S. V. Benson et al, ?The VUV/IR/THz Free Electron Laser Program at Jefferson Lab?, Nucl. Instr. and Meth. A649 9 (2011). 2. G.L. Carr, M.C. Martin, W.R. McKinney, K. Jordan, G.R. Neil and G.P. Williams ?High Power Terahertz Radiation from Relativistic Electrons?, Nature 420 153-156 (2002). 3. J.M. Klopf et al. ?The Jefferson Lab high power THz user facility?, Nucl. Instrum. Methods A582, 114-116 (2007). 4. O. Chubar, P. Elleaume, "Accurate And Efficient Computation Of Synchrotron Radiation In The Near Field Region", proc. of the EPAC98 Conference, 22-26 June 1998, p.1177-1179. 5. C. J. Hirschmugl, M. Sagurton and G. P. Williams, ?Multiparticle coherence calculations for synchrotron radiation emission? Physical Review A44, 1316 (1991).

Klopf, John↗

Advancements in Validation of TSLs through Inelastic Neutron Scattering and Transmission Measurements [Abstract]

Historically, the free gas approximation has been used to treat the thermal scattering of neutrons with energies below a few electron-volts (eV) in unevaluated materials. However, this method inadequately reproduces neutron scattering at these energies. Until recently, only a limited number of materials had available thermal scattering law (TSL) files/libraries in the ENDF nuclear data libraries in this energy range. With advancements in atomistic modeling techniques, such as molecular dynamics, ab-initio molecular dynamics, and density functional theory, TSL libraries have become available for many more materials. This is particularly relevant due to the rising interest in several advanced reactor systems that require novel moderator and reflector materials. While quasi-integral and integral benchmarks have been designed to validate historically important moderator materials (such as light water and polyethylene), there is currently a lack of standard validation methods for TSLs, especially when multiple conflicting TSLs exist. To address this issue, the Oak Ridge National Lab Nuclear Data group has been working on utilizing inelastic neutron scattering (INS) measurements combined with transmission (i.e., total cross section) measurements to evaluate and validate TSLs for different materials. We plan to demonstrate how this method has worked on materials such as polyethylene, lucite, and polystyrene. In addition, we will compare the newly created libraries to ENDF libraries for these materials and explain why integral benchmarks should not be used for validation when multiple TSLs exist.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Demonstration of Capability to Simulate Particle Irregular Shape and Poly-Disperse Mixtures Within Lunar Lander Plume-Surface Interaction

Plume-Surface Interaction (PSI) between lander engine plumes and regolith soil creates hazards in obscuration and contamination by particle clouds, high-energy ejecta streams, and landing area cratering damage. The MSFC Fluid Dynamics Branch is developing simulation tools to offer a predictive PSI capability to NASA customers such as the Human Lander System (HLS) and Commercial Lunar Payload Services (CLPS). The Gas-Granular Flow Solver (GGFS) is the main application tool for coupled gas-particle two-phase flow simulations to predict the range of PSI effects from onset of surface erosion to deep crater formation. GGFS features an Eulerian-Eulerian modeling approach, treating both gas and granular material as interacting continuum phases. Modeling the lunar regolith granular material fluidic characteristics poses special challenges due to complex particle shapes and mixture composition. The lunar regolith is poorly sorted with broad particle size distributions and large fines content. It has significant cohesion, due to interlocking jagged particle shapes. Eulerian granular material flow modeling requires closure formulations for the granular material constitutive models (stress, friction, collisional and kinetic energy dissipation, drag, etc.). While closure models for spherical particles are available from particle kinetic theory, closure models for realistic non-spherical particles must be extracted from unit physics Discrete Element Model (DEM) particle interaction simulations and provided in the form of tabular datasets. The effects of particle irregular shape (non-spherical shape factors, angular particle surface roughness, and interlocking features) are simulated by approximating the particle features in the form of grouped elemental spheres to form composite particles in the DEM simulations. The effects of the wide range of regolith mixture particle sizes and the strong effects of the presence of the small particle sizes results in high cohesion and low porosity of the regolith mixture. The range of particle sizes is simulated by binning the particle sizes into an appropriate finite number of particle-size species and solving the problem as a species mixture. Combining these two modeling approaches enables simulations to capture both, the contributions of the irregular particle shape and the particle size distribution. The integration and maturation of the DEM-based constitutive model database generation process and poly-disperse mixture binning approach into the GGFS simulation framework are proceeding under funding by the NASA Game Changing Development program. The status of current capabilities will be presented in comparisons of crater characteristics resulting for spherical and irregular shape particles, and for mono-, bi-, and tri-disperse mixture simulations of Apollo LM plume-surface interaction. The computational results confirm the significance of including the particle shape and mixture effects. Going forward plans for the full implementation of the general poly-disperse regolith modeling capability and maturation towards NASA project application readiness under the GCD program will be presented.

Peter A Liever↗

Demonstration of Capability to Simulate Particle Irregular Shape and Poly-Disperse Mixtures Within Lunar Lander Plume-Surface Interaction

Plume-Surface Interaction (PSI) between lander engine plumes and regolith soil creates hazards in obscuration and contamination by particle clouds, high-energy ejecta streams, and landing area cratering damage. The MSFC Fluid Dynamics Branch is developing simulation tools to offer a predictive PSI capability to NASA customers such as the Human Lander System (HLS) and Commercial Lunar Payload Services (CLPS). The Gas-Granular Flow Solver (GGFS) is the main application tool for coupled gas-particle two-phase flow simulations to predict the range of PSI effects from onset of surface erosion to deep crater formation. GGFS features an Eulerian-Eulerian modeling approach, treating both gas and granular material as interacting continuum phases. Modeling the lunar regolith granular material fluidic characteristics poses special challenges due to complex particle shapes and mixture composition. The lunar regolith is poorly sorted with broad particle size distributions and large fines content. It has significant cohesion, due to interlocking jagged particle shapes. Eulerian granular material flow modeling requires closure formulations for the granular material constitutive models (stress, friction, collisional and kinetic energy dissipation, drag, etc.). While closure models for spherical particles are available from particle kinetic theory, closure models for realistic non-spherical particles must be extracted from unit physics Discrete Element Model (DEM) particle interaction simulations and provided in the form of tabular datasets. The effects of particle irregular shape (non-spherical shape factors, angular particle surface roughness, and interlocking features) are simulated by approximating the particle features in the form of grouped elemental spheres to form composite particles in the DEM simulations. The effects of the wide range of regolith mixture particle sizes and the strong effects of the presence of the small particle sizes results in high cohesion and low porosity of the regolith mixture. The range of particle sizes is simulated by binning the particle sizes into an appropriate finite number of particle-size species and solving the problem as a species mixture. Combining these two modeling approaches enables simulations to capture both, the contributions of the irregular particle shape and the particle size distribution. The integration and maturation of the DEM-based constitutive model database generation process and poly-disperse mixture binning approach into the GGFS simulation framework are proceeding under funding by the NASA Game Changing Development program. The status of current capabilities will be presented in comparisons of crater characteristics resulting for spherical and irregular shape particles, and for mono-, bi-, and tri-disperse mixture simulations of Apollo LM plume-surface interaction. The computational results confirm the significance of including the particle shape and mixture effects. Going forward plans for the full implementation of the general poly-disperse regolith modeling capability and maturation towards NASA project application readiness under the GCD program will be presented.

Peter A. Liever↗

Orbiter Return-To-Flight Entry Aeroheating

The Columbia accident on February 1, 2003 began an unprecedented level of effort within the hypersonic aerothermodynamic community to support the Space Shuttle Program. During the approximately six month time frame of the primary Columbia Accident Investigation Board activity, many technical disciplines were involved in a concerted effort to reconstruct the last moments of the Columbia and her crew, and understand the critical events that led to that loss. Significant contributions to the CAIB activity were made by the hypersonic aerothermodynamic community(REF CAIB) in understanding the re-entry environments that led to the propagation of an ascent foam induced wing leading edge damage to a subsequent breech of the wing spar of Columbia, and the subsequent breakup of the vehicle. A core of the NASA hypersonic aerothermodynamics team that was involved in the CAIB investigation has been combined with the United Space Alliance and Boeing Orbiter engineering team in order to position the Space Shuttle Program with a process to perform in-flight Thermal Protection System damage assessments. This damage assessment process is now part of the baselined plan for Shuttle support, and is a direct out-growth of the Columbia accident and NASAs response. Multiple re-entry aeroheating tools are involved in this damage assessment process, many of which have been developed during the Return To Flight activity. In addition, because these aeroheating tools are part of an overall damage assessment process that also involves the thermal and stress analyses community, in addition to a much broader mission support team, an integrated process for performing the damage assessment activities has been developed by the Space Shuttle Program and the Orbiter engineering community. Several subsets of activity in the Orbiter aeroheating communities support to the Return To Flight effort have been described in previous publications (CFD?, Cavity Heating? Any BLT? Grid Generation?). This work will provide a description of the integrated process utilized to perform Orbiter tile damage assessment, and in particular will seek to provide a description of the integrated aeroheating tools utilized to perform these assessments. Individual aeroheating tools will be described which provide the nominal re-entry heating environment characterization for the Orbiter, the heating environments for tile damage, heating effects due to exposed Thermal Protection System substrates, the application of Computational Fluid Dynamics for the description of tile cavity heating, and boundary layer transition prediction. This paper is meant to provide an overall view of the integrated aeroheating assessment process for tile damage assessment as one of a sequence of papers on the development of the boundary layer transition prediction capability in support of Space Shuttle Return To Flight efforts.

Campbell, Charles H.↗

The International Space Station: A Unique Platform for Remote Sensing of Natural Disasters

Assembly of the International Space Station (ISS) was completed in 2012, and the station is now fully operational as a platform for remote sensing instruments tasked with collecting scientific data about the Earth system. Remote sensing systems are mounted inside the ISS, primarily in the U.S. Destiny Module's Window Observational Research Facility (WORF), or are located on the outside of the ISS on any of several attachment points. While NASA and other space agencies have had remote sensing systems orbiting Earth and collecting publicly available data since the early 1970s, these sensors are carried onboard free-flying, unmanned satellites. These satellites are traditionally placed into Sun-synchronous polar orbits that allow imaging of the entire surface of the Earth to be repeated with approximately the same Sun illumination (typically local solar noon) over specific areas, with set revisit times that allow uniform data to be taken over long time periods and enable straightforward analysis of change over time. In contrast, the ISS has an inclined, Sun-asynchronous orbit (the solar illumination for data collections over any location changes as the orbit precesses) that carries it over locations on the Earth between approximately 52degnorth and 52deg south latitudes (figure 1). The ISS is also unique among NASA orbital platforms in that it has a human crew. The presence of a crew provides options not available to robotic sensors and platforms, such as the ability to collect unscheduled data of an unfolding event using handheld digital cameras as part of the Crew Earth Observations (CEO) facility and on-the-fly assessment of environmental conditions, such as cloud cover, to determine whether conditions are favorable for data collection. The crew can also swap out internal sensor systems installed in the WORF as needed. The ISS orbit covers more than 90 percent of the inhabited surface of the Earth, allowing the ISS to pass over the same ground locations at different times of the day and night. This is important for two reasons: 1) certain surface processes (i.e., development of coastal fog banks) occur at times other than local solar noon, making it difficult to collect relevant data from traditional satellite platforms, and 2) it provides opportunities for the ISS to collect data for short-duration events, such as natural disasters, that polar-orbiting satellites may miss due to their orbital dynamics - in essence, the ISS can be "in the right place at the right time" to collect data. An immediate application of ISS remote sensing data collection is that the data can be used to provide information for humanitarian aid after a natural disaster. This activity contributes directly to the station's Benefits to Humanity mission. The International Charter, Space and Major Disasters (also known as the International Disaster Charter, or IDC) is an agreement between agencies of several countries to provide - on a best-effort basis - remotely sensed data related to natural disasters to requesting countries in support of disaster response. In the United States, the lead agency for interaction with the IDC is the United States Geological Survey (USGS); when an IDC request, or activation, is received, the USGS notifies the science teams for NASA instruments with targeting information for data collection. In the case of the ISS, Earth scientists in the JSC ARES Directorate, in association with the ISS Program Science Office, coordinate targeting and data collection with the USGS. If data is collected, it is passed back to the USGS for posting on its Hazards Data Distribution System and made available for download. The ISS was added to the USGS's list of NASA remote sensing assets that could respond to IDC activations in May 2012. Initially, the NASA ISS sensor systems available to respond to IDC activations included the ISS Agricultural Camera (ISSAC), an internal multispectral visible-near infrared wavelength system mounted in the WORF; CEO, a project that collects imagery through the ISS windows using off-the-shelf handheld digital visible-wavelength cameras; and the Hyperspectral Imager for the Coastal Oceans (HICO), a visible to near-infrared system mounted externally on the Japanese Experiment Module - Exposed Facility. Since May 2012, there have been 37 IDC activations; ISS sensor systems have collected data for 10 of these events.

Stefanov, William L.↗

Titan's stratospheric temperature asymmetry: a radiative origin?

During the 1981 Voyager encounter, Titan's stratosphere exhibited a large thermal asymmetry, with high northern latitudes being colder than comparable southern latitudes. Given the short radiative time constant, this asymmetry would not be expected at the season of the Voyager observations (spring equinox), if the infrared and solar opacity sources were distributed symmetrically. We have investigated the radiative budget of Titan's stratosphere, using two selections of Voyager IRIS spectra recorded at symmetric northern and southern latitudes. In the region 0.1-1 mbar, temperatures are 7 K colder at 50 degrees N than at 53 degrees S and the difference reaches approximately 13 K at 5 mbar. On the other hand, the northern region is strongly enriched in nitriles and hydrocarbons, and the haze optical depth derived from the continuum emission between 8 and 15 micrometers is twice as large as in the south. Cooling rate profiles have been computed at the two locations, using the gas and haze abundances derived from the IRIS measurements. We find that, despite lower temperatures, the cooling rate profiles in the pressure range 0.15-5 mbar are 20 to 40% larger in the north than in the south, because of the enhanced concentrations of infrared radiators. Because the northern hemisphere appears darker than the southern one in the Voyager images, enhanced solar heating is also expected to take place at 50 degrees N. Solar heating rate profiles have been calculated, with two different assumptions on the origin of the hemispheric asymmetry. In the most likely case where it results from a variation in the absorbance of the haze material, the heating rates are found to be 12-15% larger at the northern location than at the southern one, a smaller increase than that in the cooling rates. If the lower albedo in the north results from an increase in the particle number density, a 55 to 75% difference is found for the pressure range 0.15-5 mbar, thus larger than that calculated for the cooling rates. Considering the uncertainties in the haze model, dynamical heat transport may significantly contribute to the meridional temperature gradients observed in the stratosphere. On the other hand, the latitudinal variation in gas and haze composition may be sufficient to explain the entire temperature asymmetry observed, without invoking a lag in the thermal response of the atmosphere due to dynamical inertia.

NASA Program Exobiology↗

Drop Size Distribution Measurements Supporting the NASA Global Precipitation Measurement Mission: Infrastructure and Preliminary Results

Global Precipitation Measurement Mission (GPM) retrieval algorithm validation requires datasets that characterize the 4-D structure, variability, and correlation properties of hydrometeor particle size distributions (PSD) and accumulations over satellite fields of view (5 -- 50 km). Key to this process is the combined use of disdrometer and polarimetric radar platforms. Here the disdrometer measurements serve as a reference for up-scaling dual-polarimetric radar observations of the PSD to the much larger volumetric sampling domain of the radar. The PSD observations thus derived provide a much larger data set for assessing DSD variability, and satellite-based precipitation retrieval algorithm assumptions, in all three spatial dimensions for a range of storm types and seasons. As one component of this effort, the GPM Ground Validation program recently acquired five 3rd generation 2D Video disdrometers as part of its Disdrometer and Radar Observations of Precipitation Facility (DROP), currently hosted in northern Alabama by the NASA Marshall Space Flight Center and the University of Alabama in Huntsville. These next-generation 2DVDs were operated and evaluated in different phases of data collection under the scanning domain of the UAH ARMOR C-band dual-polarimetric radar. During this period approximately 7500 minutes of PSD data were collected and processed to create gamma size distribution parameters using a truncated method of moments approach. After creating the gamma parameter datasets the DSDs were then used as input to T-matrix code for computation of polarimetric radar moments at C-band. The combined dataset was then analyzed with two basic objectives in mind: 1) the investigation of seasonal variability in the rain PSD parameters as observed by the 2DVDs; 2) the use of combined polarimetric moments and observed gamma distribution parameters in a functional form to retrieve PSD parameters in 4-D using the ARMOR radar for precipitation occurring in different seasons and for different rain system types. Preliminary results suggest that seasonal variations in the DSD parameters do occur, but are most pronounced when comparing tropical PSDs to either winter or summer convective precipitation. For example the previously documented shift to relatively smaller drop diameters in higher number concentrations for equivalent rain rate bins was observed in tropical storm rainbands occurring over Huntsville. On a more inter seasonal basis empirical fits between parameters such as D0 and ZDR do not appear to exhibit robust seasonal biases- i.e., one fit seems to work for all seasons within acceptable standard error (O[10%]) for estimates of D0. In polarimetric retrievals of the vertical variability in PSD (rain layer) for a tropical rainband we find that the Do varies with height when partitioned by specified precipitation categories (e.g., convective or stratiform, heavy and light stratiform etc.) but this variation is of order 10-20% and is smaller than the difference in D0 observed between the basic delineation of convective and stratiform precipitation types. Currently we are expanding our analysis of the vertical structure of the PSD to include several seasonally and/or dynamically-different storm system types (e.g., winter convection and stratiform events; summer mid-latitude convective etc.) sampled by ARMOR. The study will present the results of our combined analyses.

Petersen, Walter A.↗

The Potential Role of Concentrating Solar Power within the Context of DOE's 2030 Solar Cost Targets

For solar electricity generating technologies to be cost competitive at a large scale with conventionally generated electricity, cost reductions are needed for both concentrating solar power (CSP) and solar photovoltaic (PV) systems. In 2011, the U.S. Department of Energy (DOE) established solar cost targets that corresponded to reducing CSP and PV prices by approximately 75% in order to achieve a levelized cost of electricity (LCOE) of $0.06 per kilowatt-hour (kWh) for both utility-scale PV and high-capacity factor CSP-TES systems in 2020. Utility-scale PV achieved its 2020 cost target in 2017, and recent estimates for the LCOE of CSP-TES with a molten-salt power tower system are approximately $0.10/kWh for projects that are expected to come online in 2020, which represents a substantial reduction since 2010 - when the LCOE for CSP-TES was around $0.21/kWh. To continue the momentum for cost reductions in solar technologies, DOE recently established cost targets for 2030 that would make solar one of the lowest-cost sources of new electricity in the United States. For CSP-based systems, the new targets correspond to an LCOE in 2030 of $0.05/kWh for a dispatchable, high-capacity factor CSP-TES plant configuration. This aggressive target would have been unimaginable a decade ago. However, building on the previously described reduction in CSP-TES costs over the past decade, recent announcements suggest the next phase of projects will continue this downward trend through lower installation costs, attractive financing, longer-duration PPAs, and the ability to capitalize on the value that the flexibility of storage brings CSP. This presentation summarizes a recent report that evaluates the potential impacts of simultaneously achieving the 2030 cost targets for PV and CSP-TES, and it includes a detailed evaluation of the role that CSP-TES could play in realizing those impacts. The scenarios in this analysis are designed to isolate and assess the potential impacts of achieving DOE's 2030 cost targets for CSP-TES and PV, and they do not reflect the potential benefits or system impacts associated with success in other DOE research programs. The low-cost solar scenarios include DOE's 2030 solar cost targets--which are represented via a roughly 50% reduction in LCOE by 2030 (from current levels) with additional cost reductions thereafter representing technology learning and/or improvements that could result from innovation - coupled with a variety of market, technology, and demand assumptions. With these assumptions, the evolution of the contiguous U.S. electricity system is evaluated with NREL's Renewable Energy Deployment System (ReEDS) model, which was specifically designed to represent the temporal and locational value of renewable generation technologies in the U.S. power system. ReEDS relies on system-wide least-cost optimization to estimate the type and location of future generation and transmission capacity. In addition, it accounts for the locational and temporal variations in variable renewable technologies, including the need for new transmission, curtailment, dynamic capacity value, and the need to hold operating reserves to account for the uncertainty and variability of these technologies. This presentation summarizes the key findings that arise from a detailed evaluation of the impacts of achieving DOE's 2030 cost targets for CSP-TES and PV systems, noting again the inherent challenges associated with modeling future scenarios of the large, complex system electricity system in the contiguous United States within the context of various market features.

14 SOLAR ENERGY↗