Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Source Term Analysis”

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,369 records · Page 76

The Analysis and Design of Low Boom Configurations Using CFD and Numerical Optimization Techniques

The use of computational fluid dynamics (CFD) for the analysis of sonic booms generated by aircraft has been shown to increase the accuracy and reliability of predictions. CFD takes into account important three-dimensional and nonlinear effects that are generally neglected by modified linear theory (MLT) methods. Up to the present time, CFD methods have been primarily used for analysis or prediction. Some investigators have used CFD to impact the design of low boom configurations using trial and error methods. One investigator developed a hybrid design method using a combination of Modified Linear Theory (e.g. F-functions) and CFD to provide equivalent area due to lift driven by a numerical optimizer to redesign or modify an existing configuration to achieve a shaped sonic boom signature. A three-dimensional design methodology has not yet been developed that completely uses nonlinear methods or CFD. Constrained numerical optimization techniques have existed for some time. Many of these methods use gradients to search for the minimum of a specified objective function subject to a variety of design variable bounds, linear and nonlinear constraints. Gradient based design optimization methods require the determination of the objective function gradients with respect to each of the design variables. These optimization methods are efficient and work well if the gradients can be obtained analytically. If analytical gradients are not available, the objective gradients or derivatives with respect to the design variables must be obtained numerically. To obtain numerical gradients, say, for 10 design variables, might require anywhere from 10 to 20 objective function evaluations. Typically, 5-10 global iterations of the optimizer are required to minimize the objective function. In terms of using CFD as a design optimization tool, the numerical evaluation of gradients can require anywhere from 100 to 200 CFD computations per design for only 10 design variables. If one CFD computation requires an hour of computational time on a Cray computer, one can see that the use of constrained numerical optimization quickly becomes impractical.Hence, in order to practically couple a numerical design optimization technique with a CFD method, the CFD method must be extremely efficient with running times on the order of only minutes. The CFD Euler code developed under NASA sponsorship and referred to as MIM3D-SB for the most part fulfills these efficiency requirements. Analysis of wing- body configurations can be computed in a matter of a few minutes. The present study will concentrate on the feasibility of the use of this CFD code in conjunction with a numerical design optimization technique for the sonic boom reduction of candidate HSCT configurations. A preliminary supersonic aircraft design system has been established that utilizes the numerical design optimization code NPSOL developed at Stanford University coupled with the supersonic NUM3D-SB CFD code. Many questions still need to be answered in regard to using CFD and numerical optimizers as design tools. There are difficulties related to both the CFD codes and the numerical optimizers. Numerical optimizers can converge to a local minima rather than a global minima. This behavior is largely a function of the initial guess in the design space. The optimizer also is searching for a minimum of the function in terms of its derivative without any regard to the actual function value. Numerically (i.e. CFD) determined gradients can also generate spurious numerical local minima. In addition, for the sonic boom problem, grid fineness will also determine the accuracy of the final design solution. Design optimization methods work well on problems defined by continuous objective functions. The sonic boom signature design problem is not necessarily defined by a continuous objective function. The signature can have a variety of shapes; i.e. from N-wave to multiple shocks. The far-field or ground signature may not transition continuously from one shape to another and hence, may exhibit discontinuous behavior. This is also a source of difficulty in using design optimization methods.In the following sections, several low boom and one reference aircraft configuration will be analyzed to predict their sonic boom signature characteristics. Modifications to some of these designs will also be presented to demonstrate the feasibility of using CFD as a design tool and to demonstrate the feasibility of designing shaped sonic boom signatures. Design modifications to some configurations will be presented to demonstrate the feasibility of achieving shaped signatures with reduced levels and not necessarily to represent realistic or aerodynamically efficient design modifications. Fuselage volume or camber are used as design variables in order to have a minimal effect on the primary wing aerodynamics. The paper will also seek to demonstrate whether a hybrid or ramped signature is feasible to achieve. For the low-boom configurations, the CFD predicted signatures will be compared qualitatively to their MLT design signatures.

Siclari, Michael J.↗

Tracking Radioactive Isotopes in HVAC and Application for Hot Cell Analyses - 20155

In recent years, the surge in the number of isotope production facilities under design has increased the need to analyze isotope migration in facility Heating, Ventilation, and Air Conditioning (HVAC) during normal operations and accident scenarios. These facilities are used to produce isotopes for medical, security, and industrial applications and as such are subject to license and regulatory requirements 10CFR20, 10CFR30, 10CFR50, and 10CFR70. GOTHIC, a general-purpose thermal-hydraulics software package, includes the ability to model isotopic tracers, radioactive decay and isotope migration as well as HEPA and charcoal filters for isotope retention. A GOTHIC model was developed by Zachry Nuclear Engineering (ZNE) to examine the effect of negative room pressure, HEPA filtration, and HVAC fluctuations on radiation areas, hot cells and gloveboxes. Radioactive tracers were used to simulate the concentration of spills within contaminated areas, track the migration of isotopes of interest, and determine the isotopic retention and buildup on facility HEPA filters. A variety of isotopes with concern to dose (e.g., Kr-85, Sr-90, I-131, etc.), including their decay and progeny, are included in the analysis. Negative pressures are maintained in the regions of interest by a representative central HVAC system equipped with a volumetric fan that exhausts to the environment after a series of isolation valves and HEPA filters. GOTHIC is an industry trusted tool for providing engineering solutions for a variety of applications, including fission product tracking, aerosol and particulate transport and ventilation assessments. The software provides an integrated analysis environment that includes a graphical user interface (GUI) for constructing analysis models, a numerical solver that includes parallel processing capabilities and a post-processor for evaluating simulation results. It solves the conservation equations for mass, momentum and energy for multicomponent, multi-phase flow in lumped parameter and multi-dimensional geometries (1, 2, or full 3D), including the effects of turbulence, diffusion and buoyancy. It has been developed and maintained under a Quality Assurance program in compliance with the requirements of 10CFR50 Appendix B and applicable portions of ASME NQA-1 since 1995. GOTHIC has been used for assessing both forced and natural convection conditions for a wide range of applications, including: - Tracking concentration of hazardous gases and chemicals for habitability and safety assessments - Determining ventilation and filtration requirements and optimizing location and arrangement of these systems - Room heat-up, including diverse and Flexible coping strategies for Extended Loss of AC Power (FLEX/ELAP) - Equipment Qualification (EQ) A distinctive feature of GOTHIC is the ability to track many different fields/substances in a simulation, including user defined tracer elements, in the liquid, vapor and droplet fields as well as surfaces and filters. This capability allows GOTHIC to model fission product transport and release or the removal of particulates or harmful toxins from exhaust gases using a spray scrubber or other types of filtration systems. GOTHIC also includes models for engineered equipment, such as fans, filters, charcoal filters, dryers/demisters, dampers, etc. The aerosols and other filtered material are removed or accumulated in these components. The range of aerosol and radiological applications that GOTHIC has been used for includes: - Source Term: Primary Coolant (Equilibrium) Activity; Non-Water Coolant Source/Leakage. - Conditions for Iodine Re-evolution: Sump/Suppression Pool Conditions and pH; - RWST Conditions and pH. - Isotope Removal Mechanisms: Containment Sprayed and Unsprayed Region Mixing; Charcoal Filter Heating due to Iodine decay. - Radionuclide Transport and Decay: Post-LOCA Release in containment; Transport between connected Compartments and vent systems; Groundwater transport of radionuclides. - Non-Newtonian Fluid modeling for sludge, waste tanks, etc. ADAMS ML071581053 (titled 'Best Practice Guidelines for the use of CFD in Nuclear Reactor Safety Applications') poses guidelines for applying single phase CFD codes in nuclear reactor safety problems and GOTHIC is listed as a 'tool for 3D flows' and 'dispersal and deposition of radionuclides.' The Nuclear Quality Assurance (NQA) pedigree of GOTHIC is an important aspect for applications in the nuclear industry. The fundamental tracer models (convective transport, molecular and turbulent diffusion, removal mechanisms, etc.) have been verified using analytical solutions and validated against applicable separate effects tests. Also, GOTHIC has been benchmarked to many integrated effects tests, including Phebus FP (Fission Product). GOTHIC gives good agreement for the buildup and decay of fission products in Phebus Test 3. The model developed by ZNE demonstrates GOTHIC's applicability and acceptability for use in analyzing the migration and retention of radioactive isotopes and their progeny in normal operation and accident scenario analyses for isotope production facilities, hot cells, and gloveboxes. The tracer activities calculated by GOTHIC can then be used in downstream radiation transport and shielding codes like RADTRAD-NAI{sup C}, MCNP{sup R}, and MicroShield{sup R} to determine on-site and of-site doses. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Multitaper Magnitude‐Squared Coherence for Time Series With Missing Data: Understanding Oscillatory Processes Traced by Multiple Observables

To explore the hypothesis of a common source of variability in two time series, observers may estimate the magnitude-squared coherence (MSC), which is a frequency-domain view of the cross correlation. For time series that do not have uniform observing cadence, MSC can be estimated using Welch's overlapping segment averaging. However, multitaper has superior statistical properties to Welch's method in terms of the tradeoff between bias, variance, and bandwidth. The classical multitaper technique has recently been extended to accommodate time series with underlying uniform observing cadence from which some observations are missing. This situation is common for solar and geomagnetic data sets, which may have gaps due to breaks in satellite coverage, instrument downtime, or poor observing conditions. We demonstrate the scientific use of missing-data multitaper magnitude-squared coherence by detecting known solar mid-term oscillations in simultaneous, missing-data time series of solar Lyman α flux and geomagnetic Disturbance Storm Time index. Due to their superior statistical properties, we recommend that multitaper methods be used for all heliospheric time series with underlying uniform observing cadence.

Astro-statistics techniques (1886)↗

A comprehensive strategy for stable, high productivity cultivation of microalgae with controllable biomass composition

This project was performed primarily as a partnership between three Universities – The University of Toledo (UToledo), Montana State University (MSU) and The University of North Carolina at Chapel Hill (UNC). North Carolina State University also contributed to the project through a sub-subcontract. The project developed fundamental and applied knowhow for cultivation systems that are not restricted by (i) proximate availability of flue gas or other high concentration CO 2 sources and (ii) the energy and infrastructure burden to deliver CO 2 over long distances. To overcome these constraints, this project focused on cultivation of alkaliphilic algal cultures in high pH, high alkalinity media (pH>10, alkalinity >100mEq) and demonstrated high biomass productivity, improvement of biomass composition for fuel production, and culture stability. The high media alkalinity ensured availability of sufficient HCO 3 - in solution for high rates of carbon fixation while the high pH of the media facilitated rapid scavenging of atmospheric CO 2 . Overall, the direct air capture rates (DAC) in our system were sufficient to achieve high biomass productivities using atmospheric CO 2 alone. The project also developed and validated toolkits related to metabolic flux modeling and microbiome characterization.

09 BIOMASS FUELS↗

Evolution of the Antarctic Ice Sheet from 2000–2300 and beyond: model sensitivity and uncertainty analysis using MPAS-Albany Land Ice

We present a description of the Antarctic Ice Sheet model configuration submitted to the ISMIP6-Antarctica-2300 experiment using the MPAS-Albany Land Ice model, along with three new sets of simulations: (1) a set of extended simulations to 2500 for three forced experiments and to 2775 for the control experiment; (2) a sensitivity analysis of our model configuration to parameters controlling basal sliding and sub-shelf melt, and to model structural choices including the choice of the energy and stress balances; and (3) a 72-member ensemble run on graphics processing units (GPUs) and analysis of variance to determine the primary sources of uncertainty in our ice-sheet model projections. Our extended simulations predict rapid retreat beginning after 2300 for SSP1-2.6 forcing and after 2500 for present-day (control) forcing, primarily in the Amundsen Sea Embayment. We find that varying the sub-shelf melt parameter between the 5th to 95th percentile values for a mean-Antarctic calibration target results in an up to ∼ ± 40 % change in sea-level contribution relative to our baseline simulations that used the median value. Using a linear basal sliding law reduces sea-level contribution by 51 %–73 % relative to our baseline nonlinear sliding law with an exponent of 1/5. When using basal sliding law exponents of 1/3 and 1/10, the overall difference from our baseline simulations at 2300 is on the order of 10 %. The Amundsen Sea Embayment region displays a strongly non-linear dependence of mass loss on the sliding law exponent, with no discernible relationship between the sliding law exponent and the mass loss by 2300, while the sectors feeding the Ross and Filchner-Ronne ice shelves exhibit more mass loss with a more-plastic sliding law. Our model fidelity sensitivity experiments reveal a 9 %–31 % increase in sea-level contribution when using a depth-integrated stress balance approximation relative to our three-dimensional solver, while using a fixed-in-time temperature field increases sea-level contribution by 14 %–88 % relative to two thermomechanically coupled configurations. Our 72-member ensemble and analysis of variance show that the uncertainty in long-term projections is dominated by the choice of Earth system model forcing and the presence or absence of hydrofracture forcing, rather than uncertainty in sliding and sub-shelf melt parameters.

58 GEOSCIENCES↗

Unveiling the Journey of A Highly Inclined CME: Insights from the March 13, 2012, Event With 110° Longitudinal Separation

Context. A fast (∼2000 km s -1 ) and wide (>110°) coronal mass ejection (CME) erupted from the Sun on March 13, 2012. Its interplanetary counterpart was detected in situ two days later by STEREO-A and near-Earth spacecraft, such as ACE, Wind, and Cluster. We suggest that at 1 au the CME extended at least 110° in longitude, with Earth crossing its east flank and STEREO-A crossing its- west flank. Despite their separation, measurements from both positions showed very similar in situ CME signatures. The solar source region where the CME erupted was surrounded by three coronal holes (CHs). Their locations with respect to the CME launch site were east (negative polarity), southwest (positive polarity) and west (positive polarity). The solar magnetic field polarity of the area covered by each CH matches that observed at 1 au in situ. Suprathermal electrons at each location showed mixed signatures with only some intervals presenting clear counter streaming flows as the CME transits both locations. The strahl population coming from the shortest magnetic connection of the structure to the Sun showed more intensity. Aims. The aim of this work is to understand the propagation and evolution of the CME and its interaction with the surrounding CHs, to explain the similarities and differences between the observations at each spacecraft, and report what one of the most longitudinal expanded CME structures measured in situ would be. Methods. Known properties of the large-scale structures from a variety of catalogues and previous studies were used to have a better overview of this particular event. In addition, multipoint observations were used to reconstruct the 3D geometry of the CME and determine the context of the solar and heliospheric conditions before the CME eruption and during its propagation. The graduated cylindrical shell model (GCS) was used to reproduce the orientation, size and speed of the structure with a simple geometry. Also, the Drag-Based Model (DBM) was utilised to understand the conditions of the interplanetary medium better in terms of the drag undergone by the structure while propagating in different directions. Finally, a comparative analysis of the different regions of the structure through the different observatories was carried out in order to directly compare the in situ plasma and magnetic field properties at each location. Results. The study presents important findings regarding the in situ measured CME on March 15, 2012, detected at a longitudinal separation of 110° in the ecliptic plane despite its initial inclination being around 45° when erupted (March 13). This suggests that the CME may have deformed and/or rotated, allowing it to be observed near its legs with spacecraft at a separation angle greater than 100°. The CME structure interacted with high-speed streams generated by the surrounding CHs. The piled-up plasma in the sheath region exhibited an unexpected correlation in magnetic field strength despite the large separation in longitude. In situ observations reveal that at both locations there was a flank encounter – where the spacecraft crossed the first part of the CME – then encountered ambient solar wind, and finally passed near the legs of the structure. Conclusions. A scenario covering all evidence is proposed for both locations with a general view of the whole structure and solar wind conditions. Also, the study shows the necessity of having multipoint observations of large-scale structures in the heliosphere.

methods: data analysis↗

Model and Standard Operating Procedures Supporting Signal Variation Flow Graph Analysis

This document will describe the principles of the Signal Variation Flow Technique, and the uncertainty models generated using it. We focus on the capture of the variation between the ideal signal and the measured signal. The ideal signal is defined to represent the signal output of a system whose full state behavior is known, with no variation in environment or during operation, and whose photons trajectories are not modified by the object. A CT uncertainty analysis is the result of two main steps. First, the radiography regime extends from the source to the collected image, which is 2D in the case of standard CT. This maps all upstream uncertainties into the variation observed on the radiograph and captures all variation in the physical domain. Second, the reconstruction regime extends from the captured images to the reconstructed 3D image. This regime is purely in the mathematical domain and corrects reconstruction algorithm artifacts/anomalies. The present work focuses on building the model through the radiography regime. The reconstruction regime is expected to be largely a study of algorithmic sensitivity, requiring the definition of a range of standardized tests through which the algorithms would be run. Radiographic variation would then be mapped through reconstruction sensitivities to predict the signal variation in the final image. An incomplete list for the reconstructed image variation output basis includes voxel density, edge blur and length variation, in analogous fashion to the basis functions presented in this work. The signal variation occurs in several forms at the radiograph. These forms are gathered into a complete basis set of functions describing all variation on the radiograph. The scale of each basis function is calculated independently via a specific SVFG. This includes 0D pixel noise (0DI), 0D energy noise (0DE), 1D blur (1DB), 1D length (1DL) or 2D position (2DP). These models are orthogonal in that they each explore a space in the signal variation domain that cannot be reached by the other basis functions. The variation basis functions, and associated SVFG models are split by output dimensionality, a term loosely used for categorization, and explained further below.

42 ENGINEERING↗

Carbon monoxide chemistry of α-V70I Mo-nitrogenase: Evidence from EPR- and IR-monitored photolysis – or, what a difference a methyl makes

A critical step in the global nitrogen cycle is the conversion of dinitrogen into biologically accessible ammonia. In Nature this is accomplished by the nitrogenase (N 2 ase) family of enzymes. Carbon monoxide (CO) has long been known as an inhibitor of dinitrogen reduction by N 2 ase, but it can also be a substrate of the enzyme, when it is catalytically reduced to hydrocarbons. Understanding the CO interactions with N 2 ases are thus relevant to both dinitrogen fixation and Fischer-Tropsch-like chemistry. Here, in this work, the interaction of CO with the α-V70I variant of Azotobacter vinelandii MoFe N 2 ase was investigated using electron paramagnetic resonance (EPR) and infrared (IR) monitored photolysis of bound CO under cryogenic conditions. This was supplemented by further analysis of stopped-flow Fourier transform IR (SF-FT-IR) data under turnover conditions. The α-V70I variant adds a single methyl group close to the FeMo-cofactor active site, and the results show that this inhibits and slows, but does not substantially chemically change, the binding of CO to the FeMocofactor. The EPR spectra of both the hi-CO and lo-CO states closely resemble those from the wild-type enzyme. Similarly, the SF-FT-IR spectra of CO inhibited α-V70I and wild-type enzyme are strikingly similar, showing only small shifts in band energies which allow better interpretation of the published wild-type spectra. The extra carbon does, however, impact and inhibit the photochemical release and migration of CO at cryogenic temperatures, resulting in novel CO-bound species. These include a product species, termed Lo-1*, which may involve CO photochemically migrating on the FeMo-cofactor.

Carbon monoxide↗

Integrated Mid-Continent Stacked Carbon Storage Hub Project Phase II (Final Summary Report)

The Phase II Integrated Midcontinent Stacked Carbon Storage Hub (IMSCS-HUB) is part of the Carbon Storage Assurance Facility Enterprise (CarbonSAFE) established by the United States Department of Energy (DOE) National Energy Technology Laboratory (NETL). CarbonSAFE is phased to support the development of commercial-scale (50 million metric tonnes [Mt] over a 30-year period) carbon capture, utilization, and storage (CCUS) in the United States. The IMSCS-HUB study area comprises carbon dioxide (CO 2 ) sources in Iowa, Kansas, and Nebraska (the source corridor), and CO 2 sinks in Kansas and Nebraska (the storage corridor), representing the first large-scale project for the Midcontinent region. The stacked storage corridor is characterized by alternating sequences of deep saline formations, oil-bearing reservoirs, shale, and evaporite units that are conducive to vertically stacked CO 2 injection for geologic storage and enhanced oil recovery (EOR). Three sites within the IMSCS-HUB stacked storage corridor were evaluated in Phase II for commercial CCUS feasibility: one in southwest-central Nebraska, Sleepy Hollow Field (SHF), a second in southwestern Nebraska near Madrid (Madrid), and a third in southwestern Kansas, the Patterson Site (composed of the Patterson, Heinitz, Hartland, and Oslo fields). In Phase II, the team assessed the feasibility of storage complexes at the potential storage sites in Nebraska and Kansas to support a commercial-scale storage hub that integrates proven CO 2 capture technology and transport from nearby ethanol sources. Building on lessons learned from the DOE-NETL Regional Carbon Sequestration Partnerships (RCSPs), the Project Team has identified a clear strategy to meet DOE’s 2025 objective of commercial carbon capture and storage (CCS) implementation by developing a CO 2 market and infrastructure that relies on multiple ethanol-based CO 2 sources in the short term and the incorporation of multiple coal-fired power plant CO 2 sources when commercial capture is economically viable. The team also leveraged the updated 45Q tax credit to develop capture and transport infrastructure. Commercial-scale CCUS is feasible at two candidate storage sites studied, the Madrid, Nebraska Site and the Patterson Site in Kearny County, Kansas. The Sleepy Hollow Field in Nebraska was found to be an attractive candidate for stacked storage with CO 2 -EOR (Battelle 2020e). Outreach efforts facilitated engagement from industry, government, and research sectors (Battelle and GPI, 2020) and an outreach plan for future phases of the project was developed to address issues that are of concern in the IMSCS-HUB project area (Battelle, 2020f). All components of a CCUS project were determined to be feasible in the IMSCS-HUB region and Risk Mitigation Plan was developed and includes strategies to mitigate risks associated with each project component (Battelle, 2020j). A roadmap was developed to obtain the required UIC permits for an integrated CCUS project (Battelle, 2020k). The regional storage resource characterization demonstrated significant opportunity for commercial-scale projects in the IMSCS-HUB storage corridor with 577.4 Mt of stacked CO 2 storage capacity and the potential to produce 181.9 MMbbls of oil via EOR across 17 individual storage areas (Battelle and ARI, 2020). The pipeline assessment study found viable pipeline routes that connected 45Q-eligible ethanol plants, coal fired power plants, and other sources in the IMSCS-HUB corridor. The comprehensive results of subsurface characterization, modeling efforts, outreach assessment, and regulatory analysis from were integrated to develop a Detailed Commercial Development Plan for the IMSCS-HUB (Battelle, 2020n). Commercialization efforts will involve obtaining Class VI UIC permits, establishing and finalizing the pipeline route, and evaluating capture projects at participating CO 2 sources. Phases I and II of the IMSCS-HUB CarbonSAFE provide a strong foundation for safely, efficiently, and cost-effectively characterizing and permitting commercial-scale project sites in the region. The plan for implementation of commercial-scale CCUS projects in the IMSCS-HUB is aligned with the objectives of CarbonSAFE Phase III: Site Characterization and CO 2 Capture Assessment.

20 FOSSIL-FUELED POWER PLANTS↗

Applications of MODIS Fluorescent Line Height Measurements to Monitor Water Quality Trends and Algal Bloom Activity

Recent advances in satellite and airborne remote sensing, such as improvements in sensor and algorithm calibrations, processing techniques and atmospheric correction procedures have provided for increased coverage of remote-sensing, ocean-color products for coastal regions. In particular, for the Moderate Resolution Imaging Spectrometer (MODIS) sensor calibration updates, improved aerosol retrievals and new aerosol models has led to improved atmospheric correction algorithms for turbid waters and have improved the retrieval of ocean color in coastal waters. This has opened the way for studying ocean phenomena and processes at finer spatial scales, such as the interactions at the land-sea interface, trends in coastal water quality and algal blooms. Human population growth and changes in coastal management practices have brought about significant changes in the concentrations of organic and inorganic, particulate and dissolved substances entering the coastal ocean. There is increasing concern that these inputs have led to declines in water quality and have increase local concentrations of phytoplankton, which cause harmful algal blooms. In two case studies we present MODIS observations of fluorescence line height (FLH) to 1) assess trends in water quality for Tampa Bay, Florida and 2) illustrate seasonal and annual variability of algal bloom activity in Monterey Bay, California as well as document estuarine/riverine plume induced red tide events. In a comprehensive analysis of long term (2003-2011) in situ monitoring data and satellite imagery from Tampa Bay we assess the validity of the MODIS FLH product against chlorophyll-a and a suite of water quality parameters taken in a variety of conditions throughout a large optically complex estuarine system. A systematic analysis of sampling sites throughout the bay is undertaken to understand how the relationship between FLH and in situ chlorophyll-a responds to varying conditions and to develop a near decadal trend in water quality changes. In situ monitoring locations that correlated well with satellite imagery were in depths greater than seven meters and were located over five kilometers from shore. Water quality parameter of total nitrogen, phosphorous, turbidity and biological oxygen demand had high correlations with these sites, as well. Satellite FLH estimates show improving water quality from 2003-2007 with a slight decline up through 2011. Dinoflagellate blooms in Monterey Bay, California (USA) have recently increased in frequency and intensity. Nine years of MODIS FLH observations are used to describe the annual and seasonal variability of bloom activity within the Bay. Three classes of MODIS algorithms were correlated against in situ chlorophyll measurements. The FLH algorithm provided the most robust estimate of bloom activity. Elevated concentrations of phytoplankton were evident during the months of August-November, a period during which increased occurrences of dinoflagellate blooms have been observed in situ. Seasonal patterns of FLH show the on- and offshore movement of areas of high phytoplankton biomass between oceanographic seasons. Higher concentrations of phytoplankton are also evident in the vicinity of the land-based nutrient sources and outflows, and the cyclonic bay-wide circulation can transport these nutrients to the northern Bay bloom incubation region. Both of these case studies illustrate the utility MODIS FLH observations in supporting management decisions in coastal and estuarine waters.

Fischer, Andrew↗

Stochastic Optimization and Uncertainty Quantification of Natrium-based Nuclear-Renewable Energy Systems for Flexible Power Applications in Deregulated Markets

Rapid integration of variable renewable energy sources (VRES) has made modeling and stochastic optimization of hybrid energy systems crucial for studying their long-term performance and viability. However, most studies have focused on just historical data, which may be unreliable for capturing short-term fluctuations, rare events, and long-term patterns of energy demand, price, and the variability of renewable energy sources. For this study, optimal synthetic time series models were developed using Wasserstein distance. The models were validated by comparing the key statistical measures against those of the historical data. They were then used to optimize the integrated Natrium-style advanced energy systems and their long-term (30 years) economics. The stochastic model performs bi-level optimization to find the optimal sizes for the balance of plant and thermal energy storage, while also optimizing energy dispatch to achieve the maximum net present value. In studies of two deregulated markets (California ISO and the Electric Reliability Council of Texas), the integrated Natrium-style system performed better in CAISO than in ERCOT, given higher and more consistent electricity prices during peak-demand periods. The potentially enlarged cost associated with the variable operation and maintenance of the TES system also plays a significant role in driving the system sizing, thus its impacts on the system are investigated in detail through comparison against a baseline case. The study also finds that the bi-level optimization results based on stochastic gradient descent closely match the grid search results. The uncertainty quantification of the stochastic signals provides further NPV-related insights and probability distributions for the case studies. The normal standard error of the mean of NPV for the case with and without TES VOM for CAISO were found to be 7.73M (plus-minus sign) 1.09M USD and 104.99M (plus-minus sign) 1.25M USD, respectively based on a 95% confidence. Given the relatively small NPV variance based on 150 samples, the analysis affords the most robust possible prediction of the techno-economic performance of the integrated Natrium-style energy systems.

25 ENERGY STORAGE↗

Task 12 Sustainability - Methodological Guidelines on Net Energy Analysis of Photovoltaic Electricity (2nd Edition)

Net Energy Analysis (NEA) is a structured, comprehensive method of quantifying the extent to which a given energy source is able to provide a net energy gain (i.e., an energy surplus) to the end user, after accounting for all the energy losses occurring along the chain of processes that are required to exploit it (i.e., for its extraction, processing and transformation into a usable energy carrier, and delivery to the end user), as well as for all the additional energy 'investments' that are required in order to carry out the same chain of processes. However, this general framework leaves the individual practitioner with a range of choices that can affect the results and thus, the conclusions of a NEA study. The current IEA PVPS guidelines were developed to provide guidance on assuring consistency, balance, and quality to enhance the credibility and reliability of the results from photovoltaic (PV) NEAs. The guidelines represent a consensus among the authors - PV NEA experts in North America and Europe - for assumptions made on PV performance, process inputs and outputs, methods of analysis, and reporting of the results. Guidance is given on photovoltaic-specific parameters used as inputs in NEA and on choices and assumptions in inventory data analysis and on implementation of modelling approaches. A consistent approach towards system modelling, the functional unit, the system boundaries and allocation aspects enhance the credibility of PV electricity NEA studies and enables balanced NEA-based comparisons. Specifically, "apples-to-oranges" comparisons of different energy carriers (e.g., fuels vs. electricity) are not methodologically sound and are to be avoided in all cases; also, any comparison across renewable and non-renewable electricity generation technologies must clearly point out the intrinsically short-term nature of the NEA viewpoint, which does not capture the long-term sustainability implications of renewable vs. non-renewable primary energy harvesting and use: non-renewable primary energy resources are depleted and finally exhausted (irrespective of the size of the EROI), while renewable primary energy resources are not. This document provides an in-depth discussion of a common metric of NEA, namely the energy return on investment (EROI), and how this is to be interpreted vis-a-vis the deceptively similar-sounding metrics in the field of Life Cycle Assessment (LCA): cumulative energy demand (CED) and non-renewable cumulative energy demand (nr-CED) per unit output. Specifically, a number of key differences are highlighted between these metrics as applied to electricity production systems, which are listed in Table S-1.

14 SOLAR ENERGY↗

Alpha Decay Chains as Thermal Power Sources: Analysis and Applications for RTGs

Radioactive sources can provide power in remote and environmentally harsh locations such as the arctic or space. The generators powered by such sources are rugged and can withstand extreme temperatures, lack of sunlight, and require no human intervention for multiple years. Radioisotopes are used in thermoelectric generators to provide power at remote sites and deep in space. Isotopes like Pu-238, Cm-244, and Am-241 are used in these generators by NASA for power in space probes and spacecrafts. These power sources deliver a steady supply of energy over extended periods of time. Alpha particles created during decay do not travel far in a material. Their kinetic energy is transferred to heat that we can then convert into energy. Unlike beta and gamma decay, the slower-moving alpha particles stop in the material, making their energy available for use. Energy from these natural decay processes provides a reliable source of power. Spontaneous fission is rare and unreliable, and unlike induced fission processes, alpha decay occurs naturally and does not require external management or ignition. The ideal properties of an isotope for use as a power source depend upon the intended use. For use in an Arctic research base over a period of several years, but less than a decade, an isotope that provides high power output over a shorter lifespan may be the most suitable option. Whereas, for deep space missions where a consistent power source for decades or perhaps more than 100 years is needed that would require a very different isotope. One with a much longer half-life that would provide consistent power throughout that time and survive in that state in for these extended periods of time. These examples represent two extreme sides in terms of time frames. By analyzing the power produced by different radioactive decay processes over time, we can evaluate the suitability of various isotope decay chains for specific uses. Some unstable isotopes undergo a series of radioactive decays, transforming into different isotopes at each step and resulting in a stable isotope. The lists of isotopes in these decay processes are known as decay chains. Some of these chains, illustrated in the figures below, are currently being investigated for use in radioisotope thermoelectric generators (RTGs) designed for a range of operational durations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Development Testing of the Gateway Integrated Bipropellant Refueling Subsystem

The Lunar Gateway is a deep space orbiting outpost being developed by the National Aeronautics and Space Administration (NASA) in partnership with the European Space Agency (ESA) and other domestic and international partners. Because Gateway is a vital component of NASA’s Artemis program supporting long-term human exploration of the moon, it is designed for on-orbit refueling to allow for a longer life performance. The reaction control system (RCS) bipropellant refueling system onboard the station will have the capability to transfer propellants, monomethyl hydrazine (MMH) and mixed oxides of nitrogen-3 (MON-3), under controlled conditions and will span across three modules: the European System Providing Refueling Infrastructure and Telecommunication Refueling Module (ESPRIT-RM or ERM), the Habitation and Logistics Outpost (HALO), and the Propulsion and Power Element (PPE). These propellant transfers are complex operations with known hazards, and one of the risk areas is exceeding the flight system’s maximum design pressure (MDP) in priming and refueling pause/stop operations. In priming, liquid propellant is transferred from a pressurized source tank to evacuated transfer lines, which could result in excessive transient surge pressures and potentially damage hardware. In refueling pause/stop operations, fast-acting isolation valves (IV) are closed, which could lead to damaging water hammer. To mitigate these risks and develop the system, a collaborative NASA/ESA/Thales Alenia Space – United Kingdom (TAS-UK) test program was completed at TAS-UK on a simplified refueling breadboard representing ERM, HALO, and PPE bipropellant transfer systems. The objectives of the integrated breadboard testing were to (1) gather performance data to characterize and demonstrate critical refueling operations, (2) help inform flight designs, and (3) to validate numerical models that can be extended to predicting flight system performance. Integrated breadboard test data has shown the architected system performance is closing initial design assumptions. Further testing with propellant on a higher fidelity fluid simulator and analysis are planned.

Adela D Han↗

Development Testing of the Gateway Integrated Bipropellant Refueling Subsystem

The Lunar Gateway is a deep space orbiting outpost being developed by the National Aeronautics and Space Administration (NASA) in partnership with the European Space Agency (ESA) and other domestic and international partners. Because Gateway is a vital component of NASA’s Artemis program supporting long-term human exploration of the moon, it is designed for on-orbit refueling to allow for a longer life performance. The reaction control system (RCS) bipropellant refueling system onboard the station will have the capability to transfer propellants, monomethyl hydrazine (MMH) and mixed oxides of nitrogen-3 (MON-3), under controlled conditions and will span across three modules: the European System Providing Refueling Infrastructure and Telecommunication Refueling Module (ESPRIT-RM or ERM), the Habitation and Logistics Outpost (HALO), and the Propulsion and Power Element (PPE). These propellant transfers are complex operations with known hazards, and one of the risk areas is exceeding the flight system’s maximum design pressure (MDP) in priming and refueling pause/stop operations. In priming, liquid propellant is transferred from a pressurized source tank to evacuated transfer lines, which could result in excessive transient surge pressures and potentially damage hardware. In refueling pause/stop operations, fast-acting isolation valves (IV) are closed, which could lead to damaging water hammer. To mitigate these risks and develop the system, a collaborative NASA/ESA/Thales Alenia Space – United Kingdom (TAS-UK) test program was completed at TAS-UK on a simplified refueling breadboard representing ERM, HALO, and PPE bipropellant transfer systems. The objectives of the integrated breadboard testing were to (1) gather performance data to characterize and demonstrate critical refueling operations, (2) help inform flight designs, and (3) to validate numerical models that can be extended to predicting flight system performance. Integrated breadboard test data has shown the architected system performance is closing initial design assumptions. Further testing with propellant on a higher fidelity fluid simulator and analysis are planned.

Adela D Han↗

Transitioning from Stage-III to Stage-IV: Cosmology from galaxy×CMB lensing and shear×CMB lensing (Final Report for Grant SC0021429 "Cosomology with Galaxy and CMB Surveys: Mitigating Systematic Effects Through Combining Datasets")

We examine the cosmological constraining power from two cross-correlation probes between galaxy and cosmic microwave background (CMB) surveys: the cross-correlation of lens galaxy density with CMB lensing convergence $\langle δ_g$$^κ\text{CMB}\rangle$, and source galaxy weak lensing shear with CMB lensing convergence $\langle^{γκ}$$_{\text{CMB}}\rangle$. These two cross-correlation probes provide an independent cross-check of other large-scale structure constraints and are insensitive to galaxy-only or CMB-only systematic effects. In addition, when combined with other large-scale structure probes, the cross-correlations can break degeneracies in cosmological and nuisance parameters, improving both the precision and robustness of the analysis. In this work, we study how the constraining power of $\langle δ_g$$^κ\text{CMB}\rangle + \langle^{γκ}$$_{\text{CMB}}\rangle$ changes from Stage-III (ongoing) to Stage-IV (future) surveys. Given the flexibility in selecting the lens galaxy sample, we also explore systematically the impact on cosmological constraints when we vary the redshift range and magnitude limit of the lens galaxies using mock galaxy catalogs. We find that in our setup, the contribution to cosmological constraints from $\langle δ_g$$^κ\text{CMB}\rangle$ and $\langle^{γκ}$$_{\text{CMB}}\rangle$ are comparable in the Stage-III data sets; but in Stage-IV surveys, the noise in$\langle δ_g$$^κ\text{CMB}\rangle$ becomes subdominant to cosmic variance, preventing $\langle δ_g$$^κ\text{CMB}\rangle$ to further improve the constraints. This implies that to maximize the cosmological constraints from future $\langle δ_g$$^κ\text{CMB}\rangle$ + $\langle^{γκ}$$_{\text{CMB}}\rangle$ analyses, we should focus more on the requirements on$^{γκ}$$_{\text{CMB}}\rangle$ instead of$\langle δ_g$$^κ\text{CMB}\rangle$. Furthermore, the selection of the lens sample should be optimized in terms of our ability to characterize its redshift or galaxy bias instead of its number density.

79 ASTRONOMY AND ASTROPHYSICS↗

Interannual Variability in Soil Trace Gas (CO2, N2O, NO) Fluxes and Analysis of Controllers

Interannual variability in flux rates of biogenic trace gases must be quantified in order to understand the differences between short-term trends and actual long-term change in biosphere-atmosphere interactions. We simulated interannual patterns (1983-1988) of global trace gas fluxes from soils using the NASA Ames model version of CASA (Carnegie-Ames-Stanford Approach) in a transient simulation mode. This ecosystem model has been recalibrated for simulations driven by satellite vegetation index data from the NOAA Advanced Very High Resolution Radiometer (AVHRR) over the mid-1980s. The predicted interannual pattern of soil heterotropic CO2 emissions indicates that relatively large increases in global carbon flux from soils occurred about three years following the strong El Nino Southern Oscillation (ENSO) event of 1983. Results for the years 1986 and 1987 showed an annual increment of +1 Pg (1015 g) C-CO2 emitted from soils, which tended to dampen the estimated global increase in net ecosystem production with about a two year lag period relative to plant carbon fixation. Zonal discrimination of model results implies that 80-90 percent of the yearly positive increments in soil CO2 emission during 1986-87 were attributable to soil organic matter decomposition in the low-latitudes (between 30 N and 30 S). Soils of the northern middle-latitude zone (between 30 N and 60 N) accounted for the residual of these annual increments. Total annual emissions of nitrogen trace gases (N2O and NO) from soils were estimated to vary from 2-4 percent over the time period modeled, a level of variability which is consistent with predicted interannual fluctuations in global soil CO2 fluxes. Interannual variability of precipitation in tropical and subtropical zones (30 N to 20 S appeared to drive the dynamic inverse relationship between higher annual emissions of NO versus emissions of N2O. Global mean emission rates from natural (heterotrophic) soil sources over the period modeled (1983-1988) were estimated at 57.1 Pg C-CO2yr-1, 9.8Tg (1012 g) N-NO yr-1, and 9.7 Tg N-N2O yr-1. Chemical fertilizer contributions to global soil N gas fluxes were estimated at between 1.3 to 7.3 Tg N-NO yr-1, and 1.2 to 4.0 Tg N-N2O yr-1.

Potter, C.↗

Radiative Effect of Springtime Biomass-Burning Aerosols over Northern Indochina During 7-SEAS Baseline 2013 Campaign

The direct aerosol radiative effects of biomass-burning (BB) aerosols over northern Indochina were estimated by using aerosol properties (physical, chemical, and optical) along with the vertical profile measurements from ground-based measurements with integration of an optical and a radiative transfer model during the Seven South East Asian Studies Biomass-Burning Aerosols Stratocumulus Environment: Lifecycles Interactions Experiment (7-SEASBASELInE) conducted in spring 2013. Cluster analysis of backward trajectories showed the air masses arriving at mountainous background site (Doi Ang Khang; 19.93degN, 99.05degE, 1536 m above mean sea level) in northern Indochina, mainly from near-source inland BB activities and being confined in the planetary boundary layer. The PM(sub10) and black carbon (BC)mass were 87 +/- 28 and 7 +/- 2 micrograms m(exp -3), respectively. The aerosol optical depth (AOD (sub 500) was found to be 0.26--1.13 (0.71 +/- 0.24). Finer (fine mode fraction is approximately or equal to 0.95, angstrom-exponent at 440-870 nm is approximately or equal to 1.77) and significantly absorbing aerosols(single scattering albedo is approximately or equal to 0.89, asymmetry-parameter is approximately or equal to 0.67, and absorption AOD 0.1 at 440 nm) dominated over this region. BB aerosols (water soluble and BC) were the main contributor to the aerosol radiative forcing (ARF), while others (water insoluble, sea salt and mineral dust) were negligible mainly due to their low extinction efficiency. BC contributed only 6 to the surface aerosol mass but its contribution to AOD was 12 (2 times higher). The overall mean ARF was 8.0 and -31.4 W m(exp -2) at top-of-atmosphere (TOA) and at the surface (SFC), respectively. Likely, ARF due to BC was +10.7 and -18.1 W m(exp -2) at TOA and SFC, respectively. BC imposed the heating rate of +1.4 K d(exp -1) within the atmosphere and highlighting its pivotal role in modifying the radiation budget. We propose that to upgrade our knowledge on BB aerosol radiative effects in BB source region, the long-term and extensive field measurements are needed.

aerosol optical properties↗