Search NASA⌕ Search

SEARCH · Search NASA

Results for “generator parameter estimation”

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,513 records · Page 84

Examination of L-Band Brightness Temperature Forecasts in the SMAP Level-4 Soil Moisture Analysis

The NASA Soil Moisture Active Passive (SMAP) mission [4] has been providing L-band (1.4 GHz) passive microwave brightness temperature (Tb) observations since April 2015. By assimilating the Tb observations into the NASA Catchment land surface model [5] using a spatially distributed ensemble Kalman filter, the NASA Global Modeling and Assimilation Office generates the SMAP Level-4 Soil Moisture (L4_SM) product, which provides global, 3-hourly, 9-km resolution estimates of surface (0-5 cm) and root-zone (0-100 cm) soil moisture with ~2.5-day latency for use in research and applications [6]. The L4_SM product also includes estimates of soil temperature, land surface fluxes, and assimilation diagnostics such as the model forecast and observed Tb values [7]. The output from the L4_SM system is routinely monitored by the L4_SM team. Such monitoring provides valuable information; instances of unusually large Tb observation-minus-forecast (O-F) residuals can indicate events during which soil moisture conditions are poorly described in the land modeling system [7,8]. For example, repeated occurrences of very large Tb O-F values in central Australia were traced back to deficiencies in the gauge-based precipitation product used in the land modeling system through L4_SM Version 5; this discovery prompted the use, beginning in Version 6, of satellite- and gauge-based precipitation observations outside of North America [9]. Recently, a statistical analysis of the Tb O-F residuals revealed systematic errors in the “tau-omega” L-band radiative transfer model that converts the land-model simulated soil moisture and temperature into the forecast Tb prior to the L4_SM analysis. The L4_SM analysis is built on the ensemble Kalman filter and assumes unbiased forecast errors. The presence of systematic Tb forecast errors could thus adversely impact the quality of the analyzed soil moisture. In this paper, we examine the Tb O-F residuals of the latest Version 7 L4_SM data (Science Version ID Vv7030 and Vv7032) [10,11,12]. In the L4_SM Version 7 algorithm, key parameters of the L-band radiative transfer model, including soil roughness, scattering albedo, and a (seasonally varying) climatology of vegetation opacity, are obtained from the SMAP Level-2 Radiometer retrieval product (Version 5) [1]. Additionally, we also use ground measurements of surface soil moisture and soil temperature from the SMAP core validation sites as in situ reference of soil conditions [2,3]. As part of the L4_SM system calibration, the seasonally-varying bias between the model forecast Tb and the observed values is removed prior to the assimilation of the SMAP Tb observations [7]. That is, the L4_SM system is designed to only correct errors in synoptic-scale and interannual variations from the long-term mean seasonal cycle while maintaining the model’s (potentially erroneous) climatology. The Tb O-F residuals examined here are thus computed using SMAP Tb observations after they are rescaled to the mean seasonal cycle of Tb from the modeling system without SMAP data assimilation. Consequently, the long-term average of the Tb O-F residuals has a global mean of only 0.13 K and locally small values, ranging from -1 to 3 K. Despite the small time-average values of the Tb O-F residuals, the model forecast Tb was nevertheless found to exhibit undesirable systematic errors. At some locations, the time-average Tb O-F values strongly depend on surface soil moisture (SM). At the Yanco SMAP core validation site, for example, the Tb O-F residuals typically range from 5 to 15 K under dry soil moisture conditions (SM < 0.15 m3 m-3) yet are predominantly negative under wet soil moisture conditions (SM > 0.25 m3 m-3), with values ranging from 0 to -40 K. This results in soil moisture analysis increments that persistently make the soil drier under dry SM conditions and persistently make the soil wetter under wet SM conditions, suggesting an error in the dynamic range of the simulated Tb, soil moisture or soil temperature. In this paper, we describe the higher-order systematic Tb forecast errors in more detail, examine their impact on the L4_SM product quality, and explore potential avenues to further improve the L4_SM algorithm.

Rolf Reichle↗

Acoustic sensors for fission gas characterization: R and D skills devoted to innovative instrumentation in MTR, non-destructive devices in hot lab facilities and specific transducers for measurements of LWR rods in nuclear plants

First of all, we will present the main principle of the method. A piezoelectric transducer, driven by a pulse generator, generates the acoustic waves in a cavity that may be the fuel rod or a chamber connected to an instrumented rod. The composition determination consists in measuring the time of flight of the acoustic signal emitted. The pressure can be estimated by a calibration process, above the measurement of the amplitude of the signal. Two projects will then be detailed. The first project consists in the development of advanced instrumentation for in-pile experiments in Material Testing Reactor. It constitutes a main goal for the improvement of the nuclear fuel behavior knowledge. This acoustic method was tested with success during a first experiment called REMORA 3, and the results were used to differentiate helium and fission gas release kinetics under transient operating conditions. This experiment was lead at OSIRIS reactor (CEA Saclay, France). As a first step of the development program, we performed in-pile tests on the most sensitive component, i.e., the piezoelectric transducer. For this purpose, the active part of this sensor has been qualified on gamma and neutron radiations and at high temperature. Various industrial piezo-ceramics were exposed to a high activity Cobalt source for few days. The cumulated dose was ranged from 50 kGy up to 2 MGy. Next, these devices were placed inside a Material Test Reactor to investigate their reliability towards neutron fluence. The final fluence after 150 days of irradiation was up to 1.6.10{sup 21}n/cm{sup 2} (for thermal neutron). Irreversible variations have been measured. Next, a specific sensor has been implemented on an instrumented fuel rod and tested in the frame of a REMORA 3 Irradiation test. It was the first experiment under high mixed, temperature neutron and gamma flux. A first irradiation phase took place in March 2010 in the OSIRIS reactor and in November 2010 for the second step of the irradiation. The instrumented fuel rod incorporating the ultrasonic gas composition sensor was finally irradiated during 2 weeks in nominal conditions. Neutronics calculation will be performed in order to calculate the thermal and fast neutron fluence and the gamma dose absorbed by acoustic sensor. A first evaluation gives a thermal fluence about 4,5.10{sup 19} n/cm{sup 2}, a fast neutrons fluence about 4,5.1018 n/cm{sup 2} and a total gamma dose up to 0,25 MGy The maximal temperature during the irradiation test was about 150 C. Although the ultrasonic sensor appears to be damaged, the optimization of the electrical attack parameters and the development of a new signal processing maintain the measurement feasibility up the end of the irradiation campaign. It was the first time that the composition of fission gas has been monitored all along an irradiation experiment in a MTR, giving access to the gas release kinetics. New researches involve thick film transducers produced by screen-printing process in order to propose piezoelectric structures for harsh temperature and irradiation measurements. The second project consists in the development of a non-destructive device that can be directly applied on a LWR fuel rod. The problem to be solved relates to the measurement of the fission gas pressure and composition in a fuel rod using a non-destructive method. Fuel rod internal pressure is one of the safety criteria applied in nuclear power analyses. This criterion must be verified in order to avoid any fuel-cladding gap reopening risk and therefore any local clad ballooning. Apart from the safety implications, this parameter is also a fuel behaviour indicator and reflects the overall fuel performance in operation, but also during shipping and long-term storage. Rod internal pressure is one criterion amongst others, like cladding corrosion, against which the acceptable fuel burn-up limit is set. A sensor has been achieved in 2007. A full-scale hot cell test of the internal gas pressure and composition measurement by an acoustic sensor was conducted successfully between 2008 and 2010 on 5 high burn-up MOX fuel rods and 2 very high burn-up UO{sub 2} fuel rods in LECA Facility at Cadarache Centre. An improvement of this sensor has been proposed, allowing us to divide by two the uncertainty on the pressure measurement. In the case of hot-cell measurements, viscous liquid can be used to couple the sensor with the rod. For gas content with a pressure exceeding 15 bars and a 10% Xe/Kr ratio, such coupling may reduce relative acoustic method accuracy by ±7% for pressure measurement result and ±0.25 % for the assessment of gas composition. These results make it possible to demonstrate the feasibility of the technique on LWR fuel rods. The transducer and the associated methodology are now operational for non-destructive measurements in hot lab facilities and allow characterising the fission gas without puncturing the fuel rods. Up to now, any other non-destructive method can be proposed. A next step will be the development of an industrial application in a fuel storage pool in order to perform a large number of measurements on a fuel assembly in nuclear plants.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Failure Analysis for Molten Salt Thermal Energy Storage Tanks for In-Service CSP Plants

Thermal Energy Storage (TES) is a fundamental component in concentrating solar power (CSP) plants to increase the plant's dispatchability, capacity factor, while reducing the levelized cost of electricity. In central receivers CSP plants, nitrate molten salts have been used for several years for operation temperatures of up to 565 degrees C. Despite many efforts to advance nitrate salt to higher operation temperatures (even considering a replacement with molten chloride salts) to achieve higher energy conversion efficiencies, the 565 degrees C temperature is currently considered the state-of-the art. Although molten salt tanks have been broadly deployed in commercial CSP plants worldwide, several failures have been reported in these tanks after a few months or years of operation, causing significant economic loss and mistrust in CSP technologies. Most of these failures are associated with the infancy of the technology and multiple issues related to tank design, fabrication, commissioning, and aggressive operation. A technical standard dedicated to the design and fabrication of molten nitrate TES tanks does not exist today. Current in-service molten salt tanks have been generally designed based on the American Petroleum Institute's (API) 650 and ASME Section II standards. The API 650 code provides guidelines for dimensions and fabrication for oil storage tanks up to 260 degrees C. The ASME standard provides allowable stress values for various materials at a range of temperatures and conditions. Both standards seem to be limited for molten salt TES tanks where high temperatures, thermal cycling, and transient conditions are expected. In 2020, NREL released the Concentrating Solar Power Best Practices Study (NREL/TP-5500-75763) that summarized multiple issues in CSP plants, along with potential alternatives and recommendations to address those issues based on information collected from participants representing about 80% of operating CSP plants in the world. One of the recommendations from this study was the development of accurate and validated models to evaluate the plant's transient operation, capable of capturing the effect of short-term clouds and operator response, while being flexible in being adapted to various spatial and temporal resource data. The "Failure Analysis for Molten Salt Thermal Energy Tanks for In-Service CSP Plants" project was inspired on this recommendation and was focused on (1) the development and validation of a physics-based model for a representative, commercial-scale molten salt tank, (2) performing simulations to evaluate the behavior of the tank as a function of typical plant operation conditions, (3) understanding tank failures mechanisms, (4) determining the residual stress and distortion in the tank floor after welding fabrication and evaluating their impact in the stresses developed in the tank during operation, (5) assessing the impact of key operation parameters on the temperature and stress distribution, (6) conduct a preliminary evaluation of design features to reduce stress and improve tank's reliability, and (7) estimate the tank's service life based on the stress developed under diverse operation scenarios. From the analysis conducted in the project and presented in this report, it was found that maximum stresses surpassing the yield strength point of the stainless steel (SS) 347H are developed on the tank floor near the perimeter. These large stresses are strongly influenced by the initial residual stresses and distortion of the tank floor after welding fabrication. During operation, large stresses are developed in the tank floor at high operation temperatures with large salt inventory levels during transient operation. High stresses are also related to elevated temperature gradients in the tank floor that could be attributed to insufficient mixing within the salt inflow and the salt inventory. Based on the analysis, creep is the predominant failure mechanism. However, the large stress levels could favor the plastic deformation into buckles, and crack formation due to stress relaxation cracking during cycle operation. A lifetime below 3 years was estimated for the typical plant operation conditions and a specific initial residual stress and deformation distribution of the tank floor. The estimated life agrees with the service time to failure reported in several commercial molten salt tanks. Desing and operation guidelines can be extracted from the analysis presented in this report, which could be adopted by tank manufacturers and CSP operators to advance toward an ultimate solution for tank failures by reducing residual and operational stresses to achieve a tank service life of more than 30 years. Addressing failures in molten salt TES tanks is fundamental for the CSP industry's survivability, but it is also important for other industrial and power generation applications using this technology, including nuclear and concentrating solar thermal.

14 SOLAR ENERGY↗

The Acoustic Emission signal acquired by the microphones placed in the CABRI test device along the fourteen last R.I.A. experiments: an example of reproducible research in nuclear science

The CABRI facility is an experimental nuclear reactor of the French Atomic Energy Commission. It is located at the Cadarache Research Centre in southern France and it is designed to act as a support to the French nuclear infrastructure. The purpose of the new testing programme termed, 'CABRI International Programme' (CIP) is to study the behaviour of PWR-type fuel rods at high burnup, equipped with an 'advanced' cladding, under Reactivity Initiated Accident (RIA) conditions (such as the scenario of a control rod ejection). Within the framework of this programme, piloted and funded by the French Institute of Nuclear Radioprotection and Safety (IRSN), ten tests are to be conducted with a frequency of two tests per year. The LPRE laboratory of the CEA which is in charge of the Preparation, realisation and breakdown of the test results studies the possibility to set up a new test analysis based on the processing of signals coming from sensors placed within the test equipment. During the experimental phase, the behaviour of the fuel element generates acoustic waves which can be detected by two microphones placed upstream and downstream of the test device. Studies showed the interest to realize temporal and spectral analyses on these signals by showing the existence of signatures which can be correlated with physical phenomena as the rod failure or the test shutdown (i.e. scram). The work presented in this article results from the will to consolidate these studies. Since the main phenomenon to be tracked is the fuel rod failure, the aim would be to highlight specific events which would have been precursors of the rod failure in order to use in the future these signals for further interpretation. Such an antecedent works resumption leads to a better understanding of the experimental needs and constitutes a good initial state to prepare the new very fast digital data acquisition systems. Although all the raw data are accessible in the form of text files, analyses and graphics representations were not straightforward to reproduce from the ancient studies since that, on one hand, people who were in charge of the original work left the laboratory and on the other hand because it is not easy when the time passes, even with our own work, to be able to remember the steps of data manipulations and the exact setup: - During the ancient experiments the use of analog data acquisition systems required to digitize tapes to be able to realize computer treatments. That had had for consequence to lose the initial dating. This one must be correctly edited to do temporal comparisons. - Analyses require functions for calculations whose parameters has to be well-known to reach the same results. We thus wished to manage our workflow in the idea that it can be easily reproducible on all the experiments. The object of the work presented in this article was to put in practice this strong bind between the data, treatments and generation of the document in order not to hesitate to do the iteration principle in action. We do not have to be afraid by the data driven analyses. According to the philosophy of the literate programming, the text of the technical document is woven with the computer code that produces all the printed output as tables, graphs for the study eliminating hence the unrealistic cut and paste. This difficulty is not specific to the nuclear domain. For many years, researchers have been worked out solutions to this mundane issue. And, presently, new technologies and high-level programming languages offer us actual answers. We will firstly present the tools applied in our laboratory to implement this workflow, then we will describe the global perception carried out to continue the study of the Acoustic Emission signals recorded by the two microphones during the fourteen last CABRI R.I.A. test.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

An Integrated Supercritical System for Efficient Produced Water Treatment and Power Generation

The main goal of this project was to evaluate the feasibility of an integrated supercritical (SC) water desalination and power generation system for treatment of produced waters from CO 2 sequestration, oil production, and coal-bed methane recovery. This system may offer a transformative approach to generating power from coal or natural gas (NG) and purifying saline water or produced water with high total dissolved solids (TDS) in a zero-liquid discharge plant. The closed-loop steam cycle of conventional power plants was replaced with an open-loop steam cycle that heated the pretreated produced water above the SC point of water (374 °C and 221 bar) to precipitate and separate the salts dissolved in the water. After separation of the precipitated salt, the SC steam required further treatment to reduce the salt content from ~100 mg/L (ppm) to <1 ppm to make it suitable for use in turbines for power generation. A SC membrane distillation unit was used to polish the SC steam, which could be sent through a series of turbines for power generation, and the steam exiting the low-pressure turbine could be condensed by cooling water to produce pure water. Project tasks included a process simulation and a techno-economic evaluation of the integrated system; design and assembly of SC salt precipitation and membrane distillation systems; development and characterization of advanced carbon membranes for SC membrane distillation; and pretreatment, desalination, and purification of different produced water samples with salt concentrations of 30,000–200,000 ppm to a high-purity (<1 ppm) level. Process simulations were performed for 550 MW SC coal-fired or NG baseline cases and various scenarios for the proposed integrated SC system using produced water as the water source and coal or NG as the fuel source. Simulation results indicated that the net plant efficiency loss of the proposed SC cogeneration process would be in the range of <1% to ~4%, in which NG and lower TDS produced water cases would have higher efficiencies. The proposed process for coal-fired cases could provide all water needs for the power generation process and generate additional pure water for beneficial use. A techno-economic analysis was performed to estimate the cost of produced water treatment (dollars /kgal of treated water) by conventional treatment processes compared with the proposed cogeneration process. The cost estimation for the proposed cogeneration system included the cost of pretreatment, the cost of power loss to integrate the water desalination process into the power generation system, and the cost of plant modification or retrofitting. The cost of produced water treatment was calculated by assuming three different scenarios in which converting a baseline power plant to the proposed cogeneration plant would increase the cost of the boiler and feed water system by 10%, 50%, and 100%. Results indicated that the cost of produced water treatment, even with a 100% cost increase, was in the range of $8.6–11.6 dollars per kgal, which is about half of the treatment cost by conventional methods. The cost of produced water treatment by the proposed cogeneration system could be further reduced or a considerable income could be generated if the values of saltwater and freshwater products are considered. Overall, the techno-economic analysis results indicated that the proposed integrated SC system could potentially treat produced water at a cost significantly lower than the cost of produced water treatment by conventional technologies. To evaluate the feasibility of the proposed treatment system at a laboratory scale, a SC salt precipitation and membrane distillation system was designed and fabricated to treat a simulated brine or pretreated produced water by precipitating salts and further polishing the SC steam containing ~100 ppm salts to a high-purity steam. Through baseline tests with NaCl solutions, with concentrations ranging from 30,000 to 200,000 ppm, it was demonstrated that concentrated feed solutions could be treated by precipitating salt at SC conditions to a TDS level of 130 ppm. Produced water samples were collected from a potential CO 2 sequestration site, a coal-bed methane site, and an oilfield in Illinois. Representative samples were characterized for various water quality parameters, including pH; turbidity; conductivity; and concentrations of selected anions, cations, alkalinity, dissolved organic carbon, total suspended solids, TDS, and total petroleum hydrocarbons. The produced water samples were first treated by conventional pretreatment methods to remove residual oil, dissolved organic matter, and scale-forming species. Pretreatment methods included coagulation, sedimentation, filtration, activated carbon adsorption, and ion-exchange adsorption. Pretreated produced waters were treated at SC conditions in which salts were precipitated. Analysis of treated samples revealed that the concentration of Na (the main cation) in produced waters was reduced dramatically from 14,000–47,400 ppm in the feed to 21–80 ppm in treated waters. Concentrations of other main cations (Ca, Mg, Sr, K) also decreased significantly from initial values of ~1,000–22,000 ppm to ~1 ppm or less. However, a significant increase in the concentrations of Fe, Ni, Cu, and Co was observed, mainly caused by leaching of these elements from the system (vessels, filters, tubing), which is manufactured from alloys containing Fe, Ni, Cu, and Co. To polish the SC steam containing a salt content of ~100 ppm to a high-purity level of <1 ppm, a SC membrane distillation treatment using carbon membranes was employed. Different types of self-standing or supported carbon membranes were developed from graphite, graphene, carbon nanotubes, carbon nanofibers, and pyrolytic carbon. Screening tests revealed that only a limited number of developed carbon membranes were suitable for the purification of steam at SC conditions. Membranes that were either fragile and lacked mechanical strength or had large pores and showed limited salt rejection (<20%) were found to be unsuitable for the purification of SC steam. Promising carbon membranes that were successfully tested for purifying SC steam included two graphite membranes (i.e., AF-5, a porous graphite material, and FGS, a flexible graphite sheet) and a carbon nanotube membrane grown on a Hastelloy substrate and then further modified using a proprietary method (i.e., membrane C). The AF-5 membrane exhibited a stable water flux of ~50 kg/m 2 ∙h and ~40% salt rejection. The FGS membrane had a high salt rejection of 95%–98% but a lower water flux of ~15 kg/m 2 ∙h. Membrane C showed an extremely high flux of 200-800 kg/m 2 ∙h, which is one order of magnitude higher than the water flux of conventional membranes. Membrane C also exhibited 87%–93% salt rejection. Overall, a two-stage or multistage polishing treatment by these membranes could lower the salt content of the SC steam from initial values of ~100–200 ppm to the steam quality of <1 ppm or to parts per billion levels required for the power cycle.

20 FOSSIL-FUELED POWER PLANTS↗

Investigating Low-Altitude Constellations of Ad-Hoc Lunar PNT System for Distributed Spacecraft Autonomy

In this study, we examine a low-altitude Lunar Position, Navigation, and Timing (LPNT) constellations and the localization performance of Centralized Extended Kalman Filter (CEKF) and Decentralized Extended Kalman Filter (DEKF) algorithms. The primary investigation involves a 100-node swarm operating at a 100 km altitude, in contrast to previous studies that examined a 21-node asset in a frozen-orbit at 5,500 km. The autonomous operation of large-scale swarm is based on two-way Inter-Satellite Link (ISL) measurements, which involve pseudoranges and relative velocities among swarm nodes. We perform a numerical assessment of the two filtering approaches, utilizing ‘fully sampled’ measurements from all available assets as well as ‘two ISL’ measurements where each spacecraft is restricted to only two antennas. This research includes an analysis of CEKF under 2-ISL constraints and evaluates the performance of DEKF in a 100-node swarm, which has not been explored in previous studies. In addition, we examine the impact of increasing the sampling frequency for DEKF, showing that the update cycle can be shortened from a 10-minute interval. A novel approach for ‘2-ISL limited’ DEKF will also be introduced, using a matching formulation that exhaustively enumerates all potential matches. This study provides valuable insights into large-scale distributed swarm operations, considering various filter configurations, sampling frequencies, matching strategies, and scalability of CEKF and DEKF for low-altitude LPNT applications. The Lunar PNT technology plays a key role in providing reliable and robust navigation services on the Moon's surface and the South pole, where the primary Lunar missions are planned. To support upcoming Lunar missions, including small satellites from NASA's Commercial Lunar Payload Services program, the Lunar PNT system must be adaptable to smaller platforms like CubeSats. Driven by the growing involvement of public and private exploration partnerships, the traditional low Earth orbit missions are shifting to beyond geosynchronous orbit [1]. These upcoming missions aim to foster a sustainable and innovative exploration program, in collaboration with commercial and international partners, to facilitate human expansion throughout the solar system and return new knowledge and opportunities to Earth [2]. As part of this trend, there are increasing efforts to utilize science missions in Lunar orbit to develop a non-dedicated and ad-hoc PNT network system. Two traditional approaches, the Deep Space Network (DSN) and the weak signal Global Positioning System (GPS), are established deep-space navigation technologies for missions beyond the geosynchronous orbit. Beginning in 1958, the DSN was developed to communicate with the Explorer 1 spacecraft based on the use of radiometric tracking in spacecraft navigation [3]. The DSN is capable of providing nearly unfettered coverage to spacecraft beyond low-Earth orbit (LEO), however, increased space mission volume has created concerns about future expectations of DSN usage for spacecraft navigation [4]. For cislunar mission applications, the position accuracy using DSN achieves 100 m (3σ) with at least three geometrically diverse ground stations when using radiometric tracking alone [5]. The DSN's dependence on Earth-based ground stations restricts its operational capabilities to periods of Earth visibility. This limitation, coupled with its poor localization performance, renders the DSN unsuitable for future lunar missions that demand continuous tracking and precise positioning. To satisfy the increasing requirements of DSN in Lunar applications, spacecrafts are also required to improve their onboard antenna power and efficiency of the transmission. However, there is an important aggregate cost trade between adding capabilities to every spacecraft and adding to a capacity on the ground that serves multiple spacecraft [6]. A weak GPS system can provide PNT service while the user spacecraft is bound to the Moon, leveraging a single, steerable high gain antenna with the relatively narrow beam which includes all the sources in its field of view [7]. However, the higher the altitude the receiver is above the GPS constellations, the poorer and the weaker are the relative geometry and the received signal powers, respectively, leading to a significant navigation accuracy reduction [8]. The transmitted power becomes weaker with increasing distance from the Earth as well as signals tracked from one of the side lobes of the GPS antenna pattern. As a results, the number of visible satellites and relative geometric condition of the GPS satellites at very high altitude drops dramatically and reduces the navigation solution accuracy. Therefore, the weak GPS system is also not an ideal way to provide PNT service to upcoming Lunar missions when considering its limited geometric condition and the recued navigation accuracy. Another navigation approach on the Moon is being developed, similar to the Global Navigation Satellite System (GNSS) on Earth, aiming to offer navigation service with continuous 24/7 coverage across the entire Lunar surface. For example, lunar communications relay and navigation systems (LCRNS) by NASA and Lunar navigation satellite systems (LNSS) by JAXA are designed to serve as dedicated Position, Navigation, and Timing (PNT) systems for the Moon. However, designing a dedicated LNSS and PNT service involves additional challenges, which are unique to the lunar environment, including limited payload capacity for the CubeSat platform, i.e., the size, weight, and power (SWaP) of the onboard clock, limited lunar ground monitoring stations, and limited financial investment as compared to the legacy Earth-GPS [9]. NASA’s focus on utilizing CubeSat platforms on the Moon leads to an alternative Lunar navigation platform that leverages the existing Lunar science and exploration assets. The small satellites used in Lunar missions can be used to create a low-cost, autonomous, ad-hoc, and on-demand mission-centric Lunar PNT swarm capable of providing PNT services to these low-cost lunar missions [10]. As upcoming Lunar missions will often operate at low-altitude about 30 km to 100 km for scientific observations and mapping purposes, the low-altitude orbital constellations could be employed to create an ad-hoc Lunar PNT system. However, several issues must be addressed, such as the instability of these orbits, which often require maintenance or are only suitable for short-duration missions, operating for fewer than 90 days. Additionally, at an altitude of 100 km, the satellites have a limited period during which they are above the horizon and capable of providing PNT service to users. The implementation of a non-dedicated, ad-hoc Lunar navigation constellation facilitates on-demand PNT services. A preliminary study of ad-hoc Lunar PNT system was conducted using 21 spacecraft in 5,5000 km altitude frozen orbits to test its feasibility and a basic performance of orbital asset localization among ad-hoc Lunar constellations in small satellites format [10]. These swarm assets are designed for autonomous localization with minimal Earth interaction, reducing dependency on bandwidth and ground resources. The design in [10] demonstrated the feasibility of a decentralized PNT approach, specifically employing a DEKF approach for state estimation, which helps minimize onboard operating costs. The DEKF method distributes computation across individual satellites, which lightens the computational load while maintaining accuracy in orbit ephemeris and clock offsets, similar to centralized systems [11]. In a follow-on study [12], each spacecraft was limited to 2 communications antennae, forcing the selection of measurements and scheduling spacecraft activities to perform the measurements. A matching algorithm is implemented to select the best measurements and schedule position estimation updates. The decentralized localization performance is also investigated with increasing levels of network degradation for swarm assets considering the impact of intermittent and permanent communication failure, to demonstrate the robustness and fidelity of the decentralized Lunar PNT service [13]. This study confirmed that the ad-hoc PNT constellations in frozen orbit are highly robust and resilient to communication failures. However, unlike frozen orbit swarm assets, the low-altitude satellites have a limited ground view at an altitude of 100 km, where the ad-hoc Lunar constellation consists of 98 low-altitude satellites, evenly distributed across seven circular polar orbital planes, alongside two satellites in a frozen orbit at an altitude of 5,500 km (Figure 1). Therefore, the number of satellites visible to ground users is significantly limited in low-altitude orbit constellations. As each visibility of a spacecraft remains intact for only a few ticks before it moves out of the field of view, the ground user encounters challenges in maintaining continuous navigation service, resulting in sparse availability and provision of Lunar PNT system. Consequently, service availability is primarily restricted to the Lunar South Pole region (Figure 2). Given these limitations and concerns, the localization performance of low-altitude swarm assets will be assessed in this study. We focus on the investigation of the localization performance of low-altitude swarm assets and ground users near the Lunar South Pole. The overall flow of the Lunar PNT simulation incorporates the DEKF approach of asset localization and the weighted least-squares approach in user localization (Figure 3). The autonomous Lunar PNT simulation is primarily implemented in MATLAB, where the DEKF based on the matching scheduler is implemented with Google’s OR-tools as a model builder and Gurobi optimization tool as a backend solver. The General Mission Analysis Tool (GMAT) is utilized to generate ephemeris data for swarm assets, and accounts for satellite orbital details, mass, and perturbations like solar radiation pressure and drag coefficients. Each ephemeris dataset is produced in the Moon International Celestial Reference Frame (ICRF) inertial coordinate system. For state estimation, the distributed swarm assets rely on two-way Inter-Satellite Link (ISL) measurements, which involve tracking pseudoranges and relative velocities between visible satellites and anchor nodes during each observation. Numerical evaluations of the decentralized localization process are conducted to demonstrate the feasibility of the low-altitude PNT system in providing reliable navigation services. The main approach involves using DEKF and CEKF to localize 100 satellites in low-altitude constellations, where the CEKF is implemented to serve as a baseline for comparing the performance of distributed algorithms. In both cases, we evaluate ‘fully sampled’ measurements from all available assets, and ‘two ISL’ measurements when spacecraft are constrained to have only two antennas. We test four estimation techniques: CEKF fully sampled, CEKF two ISL, DEKF fully sampled, and DEKF two ISL filters. As the DEKF update cycle is comprised of network setup, communication, and computations, a global broadcast network and 2-way ISL network setup will take from 4 to 6 minutes as maximum [12]. In this simulation, the DEKF update cycle is set to 10 minutes, including a 4-minute latency for obtaining and computing the actual measurement updates. We experiment an increased update cycle to demonstrate the feasibility and evaluate the impact on localization performance using various tuning values for measurement noise covariances (Figures 4 and 5). By comparing centralized and decentralized approaches using a matching algorithm, we analyze the influence of cross-correlation factors in the covariance matrix, assuming 100% reliability of all assets and measurements. The increased frequency and the adjustments of tuning parameters reveal distinct error patterns between the two scenarios. The localization accuracy of the swarm assets and ground users is assessed by taking the median error across 100 assets and one ground user (84.9°S, 137.5°E) over 7-day simulation period (Table 1). Since the user localization accuracy is significantly affected by the performance of the swarm assets, it is crucial to maintain high localization accuracy within the swarm. This study will continue to explore decentralized filtering for autonomous LPNT operations, with further investigation of an 'iterative' matching approach which enumerates every valid matching pair, planned for the following month.

Yeji Kim↗

Microstructural analysis of undoped and moderately Sc-doped TiO2 anatase nanoparticles using Scherrer equation and Debye function analysis

Highlights: • UV light sensitive undoped and moderate Sc-doped anatase nanocrystallites were prepared without post-synthesis annealing. • Elongated nanocrystallites were found by Scherrer equation and Debye Function Analysis applied to XRD data and also HRTEM. • 4 at.% Sc doping reduced the amount of Ti vacancies by half. • Doping showed ~0.3% growth of lattice parameters. - Abstract: UV light sensitive undoped and moderate (4 at.%) Sc-doped anatase nanocrystallites were prepared by an efficient and environmentally benign method based on homogeneous hydrolysis of TiOSO{sub 4} and Sc{sub 2}(C{sub 2}O{sub 4}){sub 3.5}H{sub 2}O aqueous solutions using urea as a precipitation agent, without post-synthesis calcination or annealing. X-ray diffraction study of the obtained powders revealed single phase anatase in both samples. Two different methods were applied to estimate the anatase crystallite size and shape from the X-ray diffraction data. The first one is based on the more conventional Scherrer equation (as implanted in GSAS-II program) while the second one is based on Debye function analysis (as implanted in DUBUSSY suite version 2.2), which is considered to be more adequate for nanocrystallites. Although each program has its own shape models, both methods indicated similar results of elongated nanocrystallite geometry for Sc-doped TiO{sub 2}: (i) GSAS-II confirmed ellipsoid shape with equatorial size of 5.1(0) nm and axial size of 6.8(0), (ii) DEBUSSY suite 2.2 showed cylinder shape with diameter of equivalent circle in the ab-plane, D{sub ab} = 4.9(1.0) nm and crystallite length along the c-axis, L{sub c} = 5.4(3.2) nm. For the undoped TiO{sub 2} sample, equatorial size of 5.5(0) nm and axial size of 7.3(0) nm for the ellipsoid shape was determined according to GSAS-II, while D{sub ab} = 5.5(1.3) nm and L{sub c} = 6.0(2.2) nm for the cylinder shape was determined by DEBUSSY suite 2.2. These outcomes were found to corroborate with HRTEM analysis. Doping with 4 at.% Sc caused ~0.3% growth of lattice parameters from a = 3.7979 [4] Å, c = 9.4995 [9] Å for the undoped sample to a = 3.8110 [7] Å, c = 9.5274 [16] Å for the moderately Sc doped sample, while reduction in the Ti vacancies by half was confirmed by both methods. The introducing of Sc as electronically active secondary species into the crystal lattice of TiO{sub 2} can greatly alter its optical absorption. The incorporation of Sc into the crystal lattice of TiO{sub 2} favors the substitution of Ti, as both have approximately the same ionic size, and generates surface oxygen vacancies which can act as trapping center for the photogenerated electrons and reduce the recombination of electron–hole pairs. We can expect Sc doped TiO{sub 2} nanomaterials to have applications as effective photocatalysts for selective photocatalytic oxidation of various environmental organic pollutants.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Prognostics for Systems Health Management - Model and Hybrid Based Approaches. Where are We Heading?

To facilitate and solve the prediction problem, awareness of the current state and health of the system is key, since it is necessary to perform condition-based system health predictions. To accurately predict the future state of any system, it is required to possess knowledge of its current health state and future operational conditional. In case of next generation electric aircrafts, computing remaining flying time is safety-critical, since an aircraft that runs out of power (battery charge) while in the air will eventually lose control leading to catastrophe. In order to tackle and solve the prediction problem, it is essential to have awareness of the current health state of the system, especially since it is necessary to perform condition-based predictions. To be able to predict the future state of the system, it is also required to possess knowledge of the current and future operational conditions and flight profiles for accurate estimation of end-of-discharge (EOD) for the batteries. Similar framework can be implemented to other complex systems and subsystems. Our research approach is to develop a system level health monitoring safety indicator which runs estimation and prediction algorithms to estimate remaining useful life predictions at system, subsystem swell as component levels. Given models of the current and future system behavior, a general approach of model-based prognostics is discussed as a solution to the prediction problem and further for decision making. Data driven prognostics approaches have been equally used with good results in the past, where respective approaches have their own challenges to tackle. This limits their applicability to complex real-world domains: (a) high complexity or incompleteness of physics-based models and (b) limited representativeness of the training dataset for data-driven models. With the advent of internet of things for data collection and increased use of ML algorithms, hybrid approaches are the next avenue to reduce the challenges and achieve better results. An hybrid framework for fusing information from physics-based performance models along with deep learning algorithms for prognostics of complex safety critical systems is presented. In this framework, we use physics-based performance models to infer unobservable model parameters related to the system's components health solving a calibration problem.

Prognostics↗

Cost Risk Analysis Based on Perception of the Engineering Process

In most cost estimating applications at the NASA Langley Research Center (LaRC), it is desirable to present predicted cost as a range of possible costs rather than a single predicted cost. A cost risk analysis generates a range of cost for a project and assigns a probability level to each cost value in the range. Constructing a cost risk curve requires a good estimate of the expected cost of a project. It must also include a good estimate of expected variance of the cost. Many cost risk analyses are based upon an expert's knowledge of the cost of similar projects in the past. In a common scenario, a manager or engineer, asked to estimate the cost of a project in his area of expertise, will gather historical cost data from a similar completed project. The cost of the completed project is adjusted using the perceived technical and economic differences between the two projects. This allows errors from at least three sources. The historical cost data may be in error by some unknown amount. The managers' evaluation of the new project and its similarity to the old project may be in error. The factors used to adjust the cost of the old project may not correctly reflect the differences. Some risk analyses are based on untested hypotheses about the form of the statistical distribution that underlies the distribution of possible cost. The usual problem is not just to come up with an estimate of the cost of a project, but to predict the range of values into which the cost may fall and with what level of confidence the prediction is made. Risk analysis techniques that assume the shape of the underlying cost distribution and derive the risk curve from a single estimate plus and minus some amount usually fail to take into account the actual magnitude of the uncertainty in cost due to technical factors in the project itself. This paper addresses a cost risk method that is based on parametric estimates of the technical factors involved in the project being costed. The engineering process parameters are elicited from the engineer/expert on the project and are based on that expert's technical knowledge. These are converted by a parametric cost model into a cost estimate. The method discussed makes no assumptions about the distribution underlying the distribution of possible costs, and is not tied to the analysis of previous projects, except through the expert calibrations performed by the parametric cost analyst.

Dean, Edwin B.↗

Accuracy Assessment of Global Internal-Tide Models Using Satellite Altimetry Models Using Satellite Altimetry

Altimeter measurements are corrected for several geophysical parameters in order to access ocean signals of interest, like mesoscale or sub-mesoscale variability. The ocean tide is one of the most critical corrections due to the amplitude of the tidal elevations and to the aliasing phenomena of high-frequency signals into the lower-frequency band, but the internal-tide signatures at the ocean surface are not yet corrected globally. Internal tides can have a signature of several centimeters at the surface with wavelengths of about 50–250 km for the first mode and even smaller scales for higher-order modes. The goals of the upcoming Surface Water Ocean Topography (SWOT) mission and other high-resolution ocean measurements make the correction of these small-scale signals a challenge, as the correction of all tidal variability becomes mandatory to access accurate measurements of other oceanic signals. In this context, several scientific teams are working on the development of new internal-tide models, taking advantage of the very long altimeter time series now available, which represent an unprecedented and valuable global ocean database. The internal-tide models presented here focus on the coherent internal-tide signal and they are of three types: empirical models based upon analysis of existing altimeter missions, an assimilative model and a three-dimensional hydrodynamic model. A detailed comparison and validation of these internal-tide models is proposed using existing satellite altimeter databases. The analysis focuses on the four main tidal constituents: M2, K1, O1 and S2. The validation process is based on a statistical analysis of multi-mission altimetry including Jason-2 and Cryosphere Satellite-2 data. The results show a significant altimeter variance reduction when using internal-tide corrections in all ocean regions where internal tides are generating or propagating. A complementary spectral analysis also gives some estimation of the performance of each model as a function of wavelength and some insight into the residual non-stationary part of internal tides in the different regions of interest. This work led to the implementation of a new internal-tide correction (ZARON'one) in the next geophysical data records version-F (GDR-F) standards.

Loren Carrere↗

Simulation of Physics-Based 0-10Hz Strong Motion Using High Performance Computing Supporting Refinements to Regional Ground Motion Models for the Central Eastern US

In collaboration with the U.S. Nuclear Regulatory Commission (NRC) the LLNL has developed a computationally efficient simulation platform designed to perform physics-based ground motion simulations for crustal earthquakes in the Stable Continental Regions of Central and Eastern US (CEUS), using high-performance computing. The main objective of the earthquake simulations was to use synthetic ground motion to provide constrains to refinements of existing ergodic Ground Motion Models (GMMs), for large magnitude earthquakes and near-fault distances, for which these models are less reliable. Physics-based broadband (0-10Hz) ground motion simulations were used to estimate the near-fault ground motion amplitudes and within event and between-event variabilities associated with fault rupture characteristics. In our simulations we used a 3D regional velocity model that was based on Saikia’s 1D velocity model (1994). In simulations performed during the first stage of this project the Saikia’s velocity model demonstrated better performance in modelling high frequency regional wave propagation for the CEUS region recorded during the Mw5.0 November 7, 2016, Cushing Oklahoma (Taylor et al., 2017), and Mw5.8 September 3, 2016, Pawnee Oklahoma earthquakes. The proposed regional 3D model includes random perturbations to the 1D background model using the stochastic scheme of Pitarka and Mellors (2021). In addition, validation analysis of the rupture generator and regional wave propagation models, using comparisons with different GMMs for Mw6.5 and Mw7.0 scenario earthquakes in the CEUS region resulted in a very good match between the simulated and empirical ground motion models. For the purposes of seismic hazard assessment at the existing and planned nuclear power plants, NRC is interested in studies aimed at improving the current ground motion models (GMM) for both Stable Continental Regions (SCR) in the Central and Eastern US and Active Crustal Regions (ACR) in the Western US. Due to lack of recorded data, these improvements require synthetic data for short fault distances and large magnitude earthquakes for which the existing recorded data is not enough to uniquely constrain the GMMs. The need for simulations and strong motion data is especially critical for the CEUS region where we do not have recorded data from potentially large damaging earthquakes with moment magnitudes 6.0 and higher. In this the project, we focused on 10Hz simulations of Mw7.0 scenario earthquakes with strike slip and thrust faulting mechanisms. We used more than 50 Mw7.0 earthquake rupture scenarios to investigate the ground motion uncertainty due to unknown earthquake rupture parameters, in particular, the slip distribution, rupture velocity, and faulting mechanism, and their implication on ground motion amplification due to forward rupture directivity effects.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

A comparison of cosmological hydrodynamic codes

We present a detailed comparison of the simulation results of various hydrodynamic codes. Starting with identical initial conditions based on the cold dark matter scenario for the growth of structure, with parameters h = 0.5 Omega = Omega(sub b) = 1, and sigma(sub 8) = 1, we integrate from redshift z = 20 to z = O to determine the physical state within a representative volume of size L(exp 3) where L = 64 h(exp -1) Mpc. Five indenpendent codes are compared: three of them Eulerian mesh-based and two variants of the smooth particle hydrodynamics 'SPH' Lagrangian approach. The Eulerian codes were run at N(exp 3) = (32(exp 3), 64(exp 3), 128(exp 3), and 256(exp 3)) cells, the SPH codes at N(exp 3) = 32(exp 3) and 64(exp 3) particles. Results were then rebinned to a 16(exp 3) grid with the exception that the rebinned data should converge, by all techniques, to a common and correct result as N approaches infinity. We find that global averages of various physical quantities do, as expected, tend to converge in the rebinned model, but that uncertainites in even primitive quantities such as (T), (rho(exp 2))(exp 1/2) persists at the 3%-17% level achieve comparable and satisfactory accuracy for comparable computer time in their treatment of the high-density, high-temeprature regions as measured in the rebinned data; the variance among the five codes (at highest resolution) for the mean temperature (as weighted by rho(exp 2) is only 4.5%. Examined at high resolution we suspect that the density resolution is better in the SPH codes and the thermal accuracy in low-density regions better in the Eulerian codes. In the low-density, low-temperature regions the SPH codes have poor accuracy due to statiscal effects, and the Jameson code gives the temperatures which are too high, due to overuse of artificial viscosity in these high Mach number regions. Overall the comparison allows us to better estimate errors; it points to ways of improving this current generation ofhydrodynamic codes and of suiting their use to problems which exploit their best individual features.

Kang, Hyesung↗

Electronics Upgrades to the Green is Clean Phoswich Detector Systems and Programmatic Implementation at LANL - Phase II Completion - 18403

Los Alamos National Laboratory (LANL) radiological facilities produce low-density room trash that, in many cases, is not contaminated with radioactivity. It has been estimated that 50 to 90% of low-density room trash is free of radioactive contamination and eligible for inclusion in LANL's Green is Clean (GIC) program. The GIC program is a verification program for non-regulated waste from radiation controlled areas that has been actively segregated as nonradioactive through the use of the waste generator's acceptable knowledge. The GIC program has employed two large area Phoswich gamma ray detector systems (HERCULES and ZEUS) to verify that GIC candidate low-density waste contains no radioactivity above background. The large area (127 cm{sup 2}) Phoswich detectors consist of a thin NaI front crystal (3 mm thick) that is optically coupled to a thick CsI back crystal (50 mm thick). The scintillation properties of the two crystals are different enough to allow the system's electronic components to identify the origin of any pulse. The original systems utilize Nuclear Instrumentation Module (NIM) electronics to process the detector signals. The electronics currently in use by the systems are out of date and some components are no longer manufactured. In 2015, phase I upgrades to the GIC detector systems were completed. Phase I objectives included identifying more modern commercially available electronics that could replace the NIM electronics. A pulse shape discriminator (PSD) module, the Digital Gamma Finder (DGF) Pixie-4 made by XIA LLC, that fits in a compact National Instruments{sup TM} PXI chassis was identified to replace the majority of the original NIM bin electronics - a single Pixie-4 module replaces four NIM amplifiers and four NIM pulse shape analyzers (PSAs). Another phase I objective was to show that that new PSD could separate out the NaI and CsI detector signals in a manner equivalent to or better than the NIM electronics. Proof of concept was achieved in 2015 when a prototype upgrade system on ZEUS was built and tested and gating parameter settings were determined that successfully separated the NaI and CsI signals. Phase II required establishing a LabVIEW software interface with Pixie-4 electronics in lieu of the manufacturer provided Pixie-4 software called the Pixie Viewer. In addition, phase II included 1) final optimization of the gating parameters on the prototype Pixie-4 ZEUS system, 2) finalizing Requirements, Design and Test Plan documents, 3) LabVIEW software interface quality assurance (QA) testing, 4) setting up the HERCULES system with the Pixie-4 upgrade electronics, 5) check source and background calibration of both systems, 6) absolute efficiency determination on both systems and 7) procedure revisions and implementation of both systems into routine operations. In 2017 phase II was completed including the preparation of detailed Requirements, Design and Test Plan documents. The LabVIEW interface software that communicates and controls the Pixie-4 has been completed, merged with the preexisting LabVIEW operations software and fully tested. Gating parameter optimization, configuration of HERCULES with the Pixie-4 electronics, check source and background calibrations and absolute efficiency determinations have all been performed - quantitative results are presented in the full paper. For performance validation, comparison measurements of actual GIC bags on both the HERCULES system using the old NIM Bin electronics and the ZEUS system using the upgraded Pixie-4 electronics have been done and evaluated. Finally, procedure revisions and implementation of both systems into routine operations is currently in process and planned to be complete by the end of calendar year 2017. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Modern Era Retrospective Restrospective-Analysis for Research and Applications (MERRA) Data and Services at the GES DISC

The Modern Era Retrospective-analysis for Research and Applications (MERRA) dataset is a NASA satellite era, 30 year (1979 - present), reanalysis using the Goddard Earth Observing System Data Assimilation System, Version 5 (GEOS-5). The project, run out of NASA's Global Modeling and Assimilation Office at Goddard Space Flight Center, provides the science and application communities with a state-of-the-art global analysis with emphasis on improved estimates of the hydrological cycle over a broad range of weather and climate time scales. MERRA products are generated as a long-term synthesis that places the NASA EOS suite of observations in a climate context. The MERRA analysis is performed at a horizontal resolution of 2/3 longitude x 1/2 latitude (540x361 global gridpoints) with observational analyses every 6 hours. The MERRA output data will include 3 dimensional state fields for every 6 hourly analysis cycle on 42 pressure levels (or 72 terrain following model coordinate levels) from the surface through the stratosphere. Several data products are specifically designed to support chemistry and stratosphere transport modeling. The 2 dimensional surface and atmospheric diagnostics (numbering 259) are being stored on the native grid at 1 hourly intervals. These include radiation and vertical integrals of the atmosphere for water and energy budget studies and also surface diagnostics where the diurnal cycle is important. The one hourly surface and near surface data product will also facilitate research on the integrated analysis of Earth system observations in the land, ocean and cryosphere. The MERRA products are archived and distributed by the Goddard Earth Sciences Data and Information Services Center (GES DISC) through its Modeling DISC Web (MDISC) portal. Multiple data access methods and services are available for MERRA data through MDISC: (1) Mirador offers a quick, comprehensive search of MERRA and all GES DISC archived data holdings, allowing searches on keywords, location names or latitude/longitude box, and date/time, with responses within a few seconds. (2) Giovanni is a GES DISC developed Web application that provides data visualization and analysis online. Giovanni features popular visualizations such as latitude-longitude maps, animations, cross sections, profiles, time series, etc. and some basic statistical analysis functions such as scatter plots and correlation coefficient maps. Users are able to download results in several different formats, including Google Earth. (3) On-the-fly parameter subsetting of data within a spatial/temporal window is provided through a simple select and click Web page. (4) MERRA data are also available via OPeNDAP, GrADS Data Server (GDS) and can be converted to netCDF on the fly.

Berrick, Stephen W.↗

Are ultra-faint galaxies at z = 6–8 responsible for cosmic reionization? Combined constraints from the Hubble Frontier Fields clusters and parallels

We use deep Hubble Space Telescope imaging of the Frontier Fields to accurately measure the galaxy rest-frame ultraviolet luminosity function (UV LF) in the redshift range z ∼ 6–8. We combine observations in three lensing clusters, A2744, MACS 0416, and MACS 0717, and their associated parallel fields to select high-redshift dropout candidates. We use the latest lensing models to estimate the flux magnification and the effective survey volume in combination with completeness simulations performed in the source plane. We report the detection of 227 galaxy candidates at z = 6–7 and 25 candidates at z ∼ 8. While the total survey area is about 4 arcmin{sup 2} in each parallel field, it drops to about 0.6–1 arcmin{sup 2} in the cluster core fields because of the strong lensing. We compute the UV LF at z ∼ 7 using the combined galaxy sample and perform Monte Carlo simulations to determine the best-fit Schechter parameters. We are able to reliably constrain the LF down to an absolute magnitude of M{sub UV} = −15.25, which corresponds to 0.005 L{sup ⋆}. More importantly, we find that the faint-end slope remains steep down to this magnitude limit with α=−2.04{sub −0.17}{sup +0.13}. We find a characteristic magnitude of M{sup ⋆}=−20.89{sub −0.72}{sup +0.60} and log(ϕ{sup ⋆}) = −3.54{sub −0.45}{sup +0.48}. Our results confirm the most recent results in deep blank fields but extend the LF measurements more than two magnitudes deeper. The UV LF at z ∼ 8 is not very well constrained below M{sub UV} = −18 owing to the small number statistics and incompleteness uncertainties. To assess the contribution of galaxies to cosmic reionization, we derive the UV luminosity density at z ∼ 7 by integrating the UV LF down to an observational limit of M{sub UV} = −15. We show that our determination of log(ρ{sub UV}) = 26.2 ± 0.13 (erg s{sup −1} Hz{sup −1} Mpc{sup −3}) can be sufficient to maintain reionization with an escape fraction of ionizing radiation of f{sub esc} = 10%–15%. Future Hubble Frontier Fields observations will certainly improve the constraints on the UV LF at the epoch of reionization, paving the way to more ambitious programs using cosmic telescopes with the next generation of large aperture telescopes such as the James Webb Space Telescope and the European Extremely Large Telescope.

79 ASTRONOMY AND ASTROPHYSICS↗

Diagnostics, Control and Performance Parameters for the BELLA High Repetition Rate Petawatt Class Laser

A laser system producing controllable and stable pulses with high power and ultrashort duration at high repetition rate is a key component of a high energy laser-plasma accelerator (LPA). Precise characterization and control of laser properties are essential to understanding laser-plasma interactions required to build a 10-GeV class LPA. This study discusses the diagnostics, control and performance parameters of a 1 Hz, 1 petawatt (PW) class laser at the Berkeley Lab Laser Accelerator (BELLA) facility. The BELLA PW laser provided up to 46 J on target with a 1% level energy fluctuation and 1.3-μrad pointing stability. The spatial profile was measured and optimized by using a camera, wavefront sensor, and deformable mirror (ILAO system). The focus waist was measured to be r 0 = 53 μm and the fraction of energy within the circular area defined by the first minimum of the diffraction pattern (r = 67 μm) was 0.75. The temporal profile was controlled via the angle of incidence on a stretcher and a compressor, as well as an acousto-optic programmable dispersive. The temporal pulse shape was measured to be about 33 fs in full width at half maximum (WIZZLER and GRENOUILLE diagnostics). In order to accurately evaluate peak intensity, the energy-normalized peak fluence, and energy-normalized peak power were analyzed for the measured spatial and temporal mode profiles, and were found to be 15 kJ/(cm 2 J) with 6% fluctuation (standard deviation) and 25 TW/J with 5% fluctuation for 46-J on-target energy, respectively. This yielded a peak power of 1.2 PW and a peak intensity of 17×10 18 W/cm 2 with 8% fluctuation. A method to model the pulse shape for arbitrary compressor grating distance with high accuracy was developed. The pulse contrast above the amplified spontaneous emission pedestal was measured by SEQUOIA and found to be better than 10 9 . The first order spatiotemporal couplings (STCs) were measured with GRENOUILLE, and a simulation of the pulse's evolution at the vicinity of the target was presented. A maximum pulse front tilt angle of less than 7 mrad was achieved. The reduction of the peak power caused by the first order STCs was estimated to be less than 1%. Finally, the capabilities described in this paper are essential for generation of high quality electron beams.

47 OTHER INSTRUMENTATION↗

First multi-institutional systematic comparison of the neutron ambient dose equivalent produced by proton therapy systems

Objective. Isochronous cyclotrons, synchrocyclotrons, and synchrotrons are used to accelerate protons for proton therapy. An accurate measurement of neutron doses generated by these accelerators and associated delivery systems and its clinical relevance requires systematic protocols and proper neutron dosimetry for a meaningful assessment. We present the first comprehensive comparison of neutron ambient dose equivalent (H*(10)) produced by clinically operational proton therapy systems. Approach. Treatment plans with 10 cm modulation-depth and ranges of 10 cm (R10M10) and 25 cm (R25M10) were created to cover a 10 × 10 × 10 cm 3 water target. The pencil beam scanning proton therapy machines studied were: two gantry-mounted synchrocyclotrons (Hyperscan, Mevion, half-gantry), two isochronous cyclotrons (ProBeam, Varian, full-gantry), one isochronous cyclotron (Proteus, IBA, full-gantry), and two synchrotrons (PROBEAT, Hitachi, full- and half-gantry). Proton beams were delivered to 30 × 30 × 40 cm 3 plastic water phantoms. WENDI-II and LUPIN-BF3-NP neutron rem-meters were positioned at three angles (0°, 45°, 90°) relative to the beam direction to measure the neutron H*(10) at distances between 50–300 cm from the isocenter. Main results. H*(10) showed dependence on beam energy, machine type, and measurement location. The highest reading was for the gantry-mounted synchrocyclotron, whereas other systems produced approximately comparable neutron doses. In all cases, the H*(10) reduced with distance from the isocenter. The H*(10) drop at 2 m distance compared to that at 0.5 m was a factor of ∼5 for the gantry-mounted synchrocyclotron whereas in other systems the decrease was a factor of 10. The WENDI-II device suffered from dead-time-associated under-estimation of the dose by a factor of ∼2–3 under the synchrocyclotron beam due to its high dose-per-pulse. However, WENDI-II and LUPIN-BF3-NP results were within reasonable agreement in isochronous cyclotron and synchrotron beams, indicating that both devices are suitable for those systems. Significance. Neutron H*(10) is dependent on various parameters including beam energy, measurement location, as well as machine design. Caution must be exercised in choosing the appropriate neutron-dose-measurement device to be used for low-duty-factor, particularly in high-instantaneous-rate proton delivery systems. By delivering the same volumetric proton dose across different machines, this work provides a benchmark for inter-system comparisons and serves as a foundation for future studies.

LUPIN↗

GAMA/DEVILS: constraining the cosmic star formation history from improved measurements of the 0.3–2.2 μm extragalactic background light

We present a revised measurement of the optical extragalactic background light (EBL), based on the contribution of resolved galaxies to the integrated galaxy light (IGL). The cosmic optical background radiation (COB), encodes the light generated by star formation, and provides a wealth of information about the cosmic star formation history (CSFH). We combine wide and deep galaxy number counts from the Galaxy And Mass Assembly survey (GAMA) and Deep Extragalactic VIsible Legacy Survey (DEVILS), along with the Hubble Space Telescope (HST) archive and other deep survey data sets, in nine multiwavelength filters to measure the COB in the range from 0.35 μm to 2.2 μm. We derive the luminosity density in each band independently and show good agreement with recent and complementary estimates of the optical-EBL from very high-energy (VHE) experiments. Our error analysis suggests that the IGL and γ-ray measurements are now fully consistent to within $\sim 10{{\ \rm per\ cent}}$, suggesting little need for any additional source of diffuse light beyond the known galaxy population. We use our revised IGL measurements to constrain the CSFH, and place amplitude constraints on a number of recent estimates. As a consistency check, we can now demonstrate convincingly, that the CSFH, stellar mass growth, and the optical-EBL provide a fully consistent picture of galaxy evolution. Here, we conclude that the peak of star formation rate lies in the range 0.066–0.076 M ⊙ yr -1 Mpc -3 at a lookback time of 9.1 to 10.9 Gyr.

79 ASTRONOMY AND ASTROPHYSICS↗