Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Location Description”

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.

1,413 records · Page 79

IM3/HyperFACETS Thermodynamic Global Warming (TGW) Simulation Datasets

Publication For a thorough description of the methods, see the peer-reviewed paper: Jones, A.D., Rastogi, D., Vahmani, P. et al. Continental United States climate projections based on thermodynamic modification of historical weather. Sci Data 10, 664 (2023). https://doi.org/10.1038/s41597-023-02485-5 Overview The IM3 / HyperFACETS climate simulations provide 40-year historical (1980-2019) as well as four 80-year future simulations (2020-2099) over the U.S. The future simulations are split into near (2020-2059) and far future (2060-2099) segments. The future scenarios span a range of plausible changes in future climate (both Global Circulation Model (GCM) and Representative Concentration Pathways/Shared Socioeconomic Pathway (RCP/SSP) dimensions). The simulations provide climate variables with high spatiotemporal resolution (25 hourly variables and 207 3-hourly variables at 12 km2). The datasets are generated using dynamical downscaling with the WRF (Weather Research and Forecasting) model (version 4.2.1) and therefore preserve physical consistency across variables. WRF is a state-of-the-art, fully compressible, non-hydrostatic, mesoscale numerical weather prediction model. WRF is coupled with an urban canopy model (UCM), which resolves urban surfaces. The future scenarios were developed using a thermodynamic global warming approach where past events are replayed under a range of future warming conditions. These scenarios therefore provide a perspective on potential increases in extreme event intensity, geographic scope, and duration, with previously non-extreme conditions potentially crossing new thresholds to be considered extreme by today's standards. This approach is not intended to estimate future changes in extreme event frequency that might result from changes in large-scale atmospheric dynamics. This dataset has NOT been bias corrected. A bias corrected version of selected variables is under development and will be released here when available. Scenarios Files Data for each scenario is provided in weekly NetCDF files. 25 variables are available at hourly resolution, and 207 variables are available at three-hourly resolution. Spatial resolution is 12km and spans the conterminous United States (CONUS), including some areas of Canada and Mexico, resulting in a grid of 424 by 299 cells. The spatial projection is a Lambert Conformal Conic with the following proj-string: "+proj=lcc +lat_0=40.0000076293945 +lon_0=-97 +lat_1=30 +lat_2=45 +x_0=0 +y_0=0 +R=6370000 +units=m +no_defs". The available scenarios and simulation periods are listed below: historical | 1980 - 2019 rcp45cooler | 2020 - 2059 rcp45cooler | 2060 - 2099 rcp45hotter | 2020 - 2059 rcp45hotter | 2060 - 2099 rcp85cooler | 2020 - 2059 rcp85cooler | 2060 - 2099 rcp85hotter | 2020 - 2059 rcp85hotter | 2060 - 2099 * The first year (1979, 2019, and 2059) of data within each scenario represents a model warmup period and should not be used. These are located in the `spinup_files` directory. Historical year 2020 is considered an extra year of data beyond the simulation period and can be found in the `additional_files` directory. For information on specific variables and a more in-depth discussion of methodology, please refer to the data landing page at https://tgw-data.msdlive.org. Delta Warming Files The global and CONUS warming deltas for each scenario are provided in degrees Celsius annually and monthly. Restart Files Yearly restart files are provided for each scenario which can be used to restart the WRF model at a particular point in time. Spinup Files The first year of data within each simulation period represents a model warmup period and should not be used. The files are provided here for the sake of reproducibility. Additional Files Additional years of data are provided as an extension of the historic simulation.

Jones, Andrew D.↗

The Beam Dump eXperiment

Hadronic matter makes about 14% of the known universe. The remaining 86% is Dark Matter (DM). Since it does not interact with the ordinary matter via electromagnetic force, DM is not visible and, to date, it escaped detection. The search for Dark Matter (DM) is one of the hottest topic in modern physics. Despite the increasing number of astrophysical and cosmological observations proving the existence, so far no particle physics experiment has detected DM yet. Up to now, most of the experimental efforts have been focused on the so called WIMP (Weakly Interacting Massive Particle) paradigm, which predicts heavy DM (10 GeV-10 TeV mass range) interacting with Standard Model (SM) particles via the weak force mediators (W or Z bosons). More recently, due to the lack of clear evidence of WIMPs, other models of DM gained the interest of the physics community. These models consider Light DM particles (LDM), in MeV-GeV mass range. Among LDM theories, the Dark Photon theory predicts the existence of a Dark Sector interacting with SM particles via a new massive vector boson (Dark Photon, Heavy Photon or A0), mediator of a new force. This scenario, despite being theoretically well motivated, is remarkably experimentally unexplored. The Beam Dump eXperiment (BDX), is an approved experiment at Jefferson Lab (JLab), aiming to discover the DM predicted witihn the Dark Photon theory. The experiment uses the CEBAF (Continuous Electron Beam Accelerator Facility) 11 GeV electron beam, impinging on the JLab Hall-A beam-dump, to produce a beam of DM particles, detected by a ~ 1 m3 detector made of thallium doped cesium iodide (CsI(Tl)) crystals located ~ 20 m downstream. In order to achieve excellent background rejection, the experimental setup includes active vetos and passive shielding surrounding CsI(Tl) crystals. Given the weakness of the DM-SM interaction, the scattering of a DM particle in the BDX detector is a rare event. Therefore, the characterization of the expected background is a critical aspect of the experiment. Both cosmogenic and penetrating SM particles produced by the beam interaction in the dump contribute to the background of the experiment. While the cosmogenic contribution can be measured during the experiment when the beam is off, beam-related background can only be estimated via Monte Carlo (MC) simulations. Therefore, a careful assessment of possible systematics introduced by the MC is needed. The beam-correlated background characterization was made by measuring the muon flux produced by the interaction of the CEBAF beam with Hall-A dump in a dedicated experimental campaign in spring 2018 [1]. The comparison to the ?ux predicted by MC allowed us to validate the BDX simulation framework. The reach of BDX, i.e. the region that the experiment can probe in the LDM theory parameter space, depends critically on the background rejection capability and signal detection effciency. For this reason, the detector setup was fine-tuned through a dedicated study. The response to LDM and background events was evaluated for different setups and selection cuts. As a result of this procedure, the configuration resulting in the best sensitivity was selected [2]. The reach calculation was performed taking into account the different LDM production mechanisms, including the contribution of secondary particles produced in the beam-dump. In particular, the effect of the secondary positrons annihilation was found to be extremely significant. In the majority of the sensitivity studies, this contribution is neglected, but recent results [3] [4] demonstrated that this process significantly enhances the sensitivity of lepton beam-dump experiments. Currently, the BDX collaboration is focused on the deployment and operation of a small detector, called BDX-MINI, built to perform a preliminary physics measurement searching for LDM at JLab. This test will pave the way to the realization of the full BDX experiment. The measurement is currently ongoing but results are expected by the end of this year. During my PhD I was involved in all aspects of the BDX experiment: design, simulation, prototyping and data analysis. The main results of my work are reported in this thesis. This manuscript is organized as follows: the first Chapter provides an introduction to the theory of LDM, with particular attention to the Dark Photon paradigm; the second Chapter illustrates the BDX experimental setup, the LDM production and detection mechanisms and the expected backgrounds; Chapter 3 and 4 describe, respectively, the BDX-HODO measurement, with a detailed description of the simulations, and the BDX experimental setup and analysis cuts optimization. Chapter 5 reports about BDXMINI detector characterization, calibration and sensitivity estimate. Finally, Chapter 6 describes in detail the calculation of the secondary positron annihilation contribution to the sensitivity of BDX and other electron-beam thick-target experiment.

Marsicano, Luca↗

HarDWR - Harmonized Water Rights Records

For a detailed description of the database of which this record is only one part, please see the HarDWR meta-record. Here we present a new dataset of western U.S. water rights records. This dataset provides consistent unique identifiers for each spatial unit of water management across the domain, unique identifiers for each water right record, and a consistent categorization scheme that puts each water right record into one of 7 broad use categories. These data were instrumental in conducting a study of the multi-sector dynamics of intersectoral water allocation changes through water markets (Grogan et al., in review). Specifically, the data were formatted for use as input to a process-based hydrologic model, WBM, with a water rights module (Grogan et al., in review). While this specific study motivated the development of the database presented here, U.S. west water management is a rich area of study (e.g., Anderson and Woosly, 2005; Tidwell, 2014; Null and Prudencio, 2016; Carney et al, 2021) so releasing this database publicly with documentation and usage notes will enable other researchers to do further work on water management in the U.S. west. The raw downloaded data for each state is described in Lisk et al. (in review), as well as here. The dataset is a series of various files organized by state sub-directories. The first two characters of each file name is the abbreviation for the state the in which the file contains data for. After the abbreviation is the text which describes the contents of the file. Here is each file type described in detail: XXFullHarmonizedRights.csv: A file of the combined groundwater and surface water records for each state. Essentially, this file is the merging of XXGroundwaterHarmonizedRights.csv and XXSurfaceWaterHarmonizedRights.csv by state. The column headers for each of this type of file are: state - The name of the state the data comes from. FIPS - The two digit numeric state ID code. waterRightID - The unique identifying ID of the water right, the same identifier as its state uses. priorityDate - The priority date associated with the right. origWaterUse - The original stated water use from the state. waterUse - The water use category under the unified use categories established here. source - Whether the right is for surface water or groundwater. basinNum - The alpha-numeric identifier of the WMA the record belongs to. CFS - The maximum flow of the allocation in cubic feet per second (ft3s-1). Arizona is unique among the states, as its surface and groundwater resources are managed with two different sets of boundaries. So, for Arizona, the basinNum column is missing and instead there are two columns: surBasinNum - The alpha-numeric identifier of the surface water WMA the record belongs to. grdBasinNum - The alpha-numeric identifier of the groundwater WMA the record belongs to. XXStatePOD.shp: A shapefile which identifies the location of the Points of Diversion for the state's water rights. It should be noted that not all water right records in XXFullHarmonizedRights.csv have coordinates, and therefore may be missing from this file. XXStatePOU.shp: A shapefile which contains the area(s) in which each water right is claimed to be used. Currently, only Idaho and Washington provided valid data to included within this file. XXGroundwaterHarmonizedRights.csv: A file which contains only harmonized groundwater rights collected from each state. See XXFullHarmonizedRights.csv for more details on how the data is formatted. XXSurfaceWaterHarmonizedRights.csv: A file which contains only harmonized surface water rights collected from each state. See XXFullHarmonizedRights.csv for more details on how the data is formatted. Additionally, one file, stateWMALabels.csv, is not stored within a sub-directory. While we have referred to the spatial boundaries that each state uses to manage its water resources as WMAs, this term is not shared across all states. This file lists the proper name for each boundary set, by state.

Economics↗

Acoustic Phased Array Quantification of Quiet Technology Demonstrator 3 Advanced Inlet Liner Noise Component

Acoustic phased array flyover noise measurements were acquired as part of the Boeing 737 MAX-7 NASA Advanced Inlet Liner segment of the Quiet Technology Demonstrator 3 (QTD3) flight test program. This paper reports on the processes used for separating and quantifying the engine inlet, exhaust and airframe noise source components and provides sample phased array-based comparisons of the component noise source levels associated with the inlet liner treatment configurations. Full scale flyover noise testing of NASA advanced inlet liners was conducted as part of the Quiet Technology Demonstrator 3 flight test program in July and August of 2018. Details on the inlet designs and testing are provided in the companion paper of Reference 1. The present paper provides supplemental details relating to the acoustic phased array portion of the analyses provided in Ref. 1. In brief, the test article was a Boeing 737MAX-7 aircraft with a modified right hand (starboard side) engine inlet, which consisted of either a production inlet liner, a NASA designed inlet liner or a simulated hard wall configuration (accomplished by applying speed tape over the inlet acoustic treatment areas). In all three configurations, the engine forward fan case acoustic panel was replaced with a unperforated (hardwall) panel. No other modifications to any other acoustic treatment areas were made. The left hand (port side) engine was a production engine and was flown at idle thrust for all measurements in order to isolate the effects of the inlet liners to the right hand engine. As described in Ref. 1, the NASA inlet treatment consists of laterally cut slots (cut perpendicular to the flow direction) which are designed to reduce excrescence drag while maintaining or exceeding the liner acoustic noise reduction capabilities. The NASA inlet liner consists of a Multi-Degree of Freedom (MDOF) design with two breathable septum layers inserted into each honeycomb cell [1]. The aircraft noise measurements were acquired for both takeoff (flaps 1 setting, gear up) and approach (flaps 30 gear up and gear down) configurations. The inlet and flight test configurations are summarized in Table 1. Table 1: Inlet Treatment and Flight Configurations Inlet Forward Fan Case Aircraft Production Hardwall Flaps 1, gear up; flaps 30 gear up; flaps 30 gear down NASA Hardwall Flaps 1, gear up; flaps 30 gear up; flaps 30 gear down Hardwall Hardwall Flaps 1, gear up; flaps 30 gear up; flaps 30 gear down III.Test Description and Hardware The flight testing was conducted at the Grant County airport in Moses Lake, WA, between 27 July and 6 August 2018. The noise measurement instrumentation included 8 flush dish microphones arranged in a noise certification configuration as well as an 840 microphone phased array. The flush dish microphones were used to quantify the levels and differences in levels between the various inlet treatments. The phased array was used to separate and quantify the narrowband (tonal) and broadband noise component levels from the engine inlet/exhaust and from the airframe. Phased array extraction of the broadband component was critical to this study because it allowed for the separation of the inlet component from the total airplane level noise even when it was significantly below the total level. Figure 1 provides an overview of the phased array microphone layout as well as a detailed image of an individual phased array microphone mounted in a plate holder (the microphone sensor is the dot in the center of the plate). The ground plane ensemble array microphones (referred to as “ensemble array” in this paper) were mounted in plates with “flower petal” edges designed to minimize edge scattering effects. Fig. 1 Flyover test microphone layout. The phased array configuration was the result of a progressive development of concepts originally implemented in Ref. 2 and refined over the following years, consisting namely of multiple multi-arm logarithmic spiral subarrays designed to cover overlapping frequency ranges and optimized for various aircraft emission angles. For the present case, the signals from all 840 microphones were acquired on a single system. The 840 microphones were parsed into 11 primary subarray sets spanning from smallest to largest aperture size and labeled accordingly as “a, b, …, k”, where “a” corresponds to the smallest fielded subarray and “k” corresponds to the largest aperture subarray. The apertures ranged from approximately 10 ft to 427 ft in size (in the flight direction) with the subarrays consisting of between 215 and 312 microphones. Figure 2 shows three such subarrays, k, h and a. As done in Ref. 2, microphones were shared between subarrays in order to reduce total channel count. Fig. 2 Sample subarray sizes (20° from overhead – refer to Figure 3a discussion). In addition to the above, each of the 11 primary subarray sets consisted of four subarrays optimized to provide near equivalent array spatial resolution in both the flight and lateral directions within 30 degrees of overhead (i.e., airplane directly above the center of the array), namely, at angles  of 0, ±10, ±20 and ±30 degrees relative to overhead where angle  is defined as shown in Figure 3a. This allowed for optimized aircraft noise measurements from 60 to 120 degree emission angle.6 An example of this pletharray design is shown in Figure 3b for the k subarray. When the aircraft is at overhead, the microphones indicated by the blue markers are used for beamforming. When the aircraft is at angles ±10 degrees from overhead, both the blue and red colored microphones are used, and so on for the ±20 and ±30 degree aircraft locations. See Ref. 3 for extensive details on pletharray design for aeroacoustic phased array testing. 6 In the discussions that follow, emission angle values are used. These are the angles at the time sound is emitted relative to the engine axis and are calculated based on flight path angle, body aircraft body angle with respect to the relative wind direction, and engine axis angle relative to aircraft body angle.

Brusniak, Leon↗

TRACER-Coastal Urban Boundary-Layer Interactions with Convection (TRACER-CUBIC) Field Campaign Report

To better understand the complicated web of processes governing convective cloud life cycle and aerosol-convection interactions, the U.S. Department of Energy (DOE)’s Atmospheric Radiation Measurement (ARM) user facility supported deployment of a variety of advanced atmospheric measurement systems to the greater Houston, Texas, area from 1 October 2021 to 30 September 2022 as part of the Tracking Aerosol Convection Interactions Experiment (TRACER). Houston was selected as a study area because isolated convection and a variety of aerosol conditions are common in this region. This one-year ARM Mobile Facility (AMF) deployment featured a four-month intensive operational period (IOP) during summer 2022 (1 June–30 September). The ARM instrumentation was deployed at three sites along an east-west transect from La Porte, Texas to an ancillary site in a less-polluted rural region southwest of downtown Houston (Figure 1). At the La Porte Site, which is located near the Houston ship channel in an area that experiences significant pollution, the first ARM Mobile Facility (AMF1) was deployed. During the IOP, the ARM tethered balloon system (TBS) operated at the ancillary site. The second-generation C-Band Scanning ARM Precipitation Radar (CSAPR) operated near Pearland, Texas, roughly halfway between the Laporte and ancillary sites. As part of the TRACER- Coastal Urban Boundary-Layer Interactions with Convection (CUBIC) project, three boundary-layer profiling systems) were deployed along a north-south transect spanning from the University of Houston Coastal Center to the Aldine site north of downtown Houston (also blue dot in Figure 1) during the TRACER IOP. These systems included the National Oceanic and Atmospheric Administration (NOAA) National Severe Storms Laboratory CLAMPS2 (C2), which was deployed at the UHCC, the University of Wisconsin SPARC, which was deployed at the ARM CSAPR site near Pearland (orange diamond in middle of map in Figure 1), and the University of Oklahoma CLAMPS1 (C1), which was deployed at Aldine. These three systems have been successfully operated in various field campaigns, providing data sets that collectively offer new insights into atmospheric-boundary-layer (ABL) processes, sea-breeze (SB) circulations, and convection initiation (CI). For the TRACER IOP window, these systems ran continuously between 1 June and 26 September, 2022. Due to commitments to NOAA projects, the Doppler lidar at the UHCC site was not available until 24 June 2022. The CLAMPS and SPARC profiling systems are self-contained platforms that have benefited from several years of development and deployment. Instruments and data processing were maintained remotely, which made the 4-month deployment for the TRACER-CUBIC IOP period possible. The same basic instrument configuration comprises each system: a scanning Doppler wind lidar for flow characterization and passive profiler(s) for characterizing planetary-boundary-layer (PBL) thermodynamic properties. Each platform includes a Halo Streamline Doppler wind lidar, an Atmospheric Emitted Radiance Interferometer (AERI), and a surface meteorology station. CLAMPS1 and CLAMPS2 each also include a microwave radiometer (MWR, Figure 1d). The TRACER-CUBIC hypotheses included (i) Interactions of SB and urban circulations and how they affect the PBL structure in the Houston environment, causing spatially (horizontally and vertically) and temporally highly variable flow patterns, (ii) heat, moisture, and aerosol transport and mixing depend on these flow dynamics, and (iii) an improved understanding of the flow patterns and PBL structure are critical for investigating the processes leading to CI. To test these hypotheses, the project aimed at (i) characterizing SB circulations and their impacts on the diurnal evolution of the structure of the ABL, (ii) studying the evolution of Houston’s complex urban boundary layer, and (iii) identifying effects of urban-induced circulations on pre-convective environments. TRACER-CUBIC observations generally provide good coverage during the summer IOP. Initial screening of the data indicates a good number of cases with bay-breeze (BB) and/or SB signatures, local CI, and interesting boundary-layer features such as strong nocturnal low-level jets (LLJs, Table 1). The numbers listed in rows 3-5 in this table will be further updated as part of ongoing in-depth analyses and systematic identification of local circulations and CI events. More detailed information about the data availability and quality for each instrument is provided in the “readme” files that were submitted to the ARM Data Center along with each archived data sets. These “readme” files also provide instrument descriptions, information about the data collection and processing procedures, data formats, and any additional information relevant for further data analysis.

54 ENVIRONMENTAL SCIENCES↗

HarDWR - Harmonized Water Rights Records

For a detailed description of the database of which this record is only one part, please see the HarDWR meta-record. Here we present a new dataset of western U.S. water rights records. This dataset provides consistent unique identifiers for each spatial unit of water management across the domain, unique identifiers for each water right record, and a consistent categorization scheme that puts each water right record into one of 7 broad use categories. These data were instrumental in conducting a study of the multi-sector dynamics of intersectoral water allocation changes through water markets (Grogan et al., in review). Specifically, the data were formatted for use as input to a process-based hydrologic model, WBM, with a water rights module (Grogan et al., in review). While this specific study motivated the development of the database presented here, U.S. west water management is a rich area of study (e.g., Anderson and Woosly, 2005; Tidwell, 2014; Null and Prudencio, 2016; Carney et al, 2021) so releasing this database publicly with documentation and usage notes will enable other researchers to do further work on water management in the U.S. west. The raw downloaded data for each state is described in Lisk et al. (in review), as well as here. The dataset is a collection of objects stored within an RDate file, stateWaterRightsHarmonized.RData. Each object either describes how one should interact with the other objects, or contains the harmonized water rights data. Here is each object described in detail: states: A character vector containing the state names for those states in which data was collected for. More importantly, the index of the state name is also the index in which that state's data can be found in the various following list objects. For example, if California is the third index in this object, the data for California will also be in the third index for each accompanying list. rightsByState_ground: A list of data frames with the cleaned ground water rights collected from each state. The column headers for each data frame are: waterRightID - The unique identifying ID of the water right, the same identifier as its state uses. basinNum - The alpha-numeric identifier of the WMA the record belongs to. priorityDate - The priority date associated with the right. allocatedFlow - The maximum flow of the allocation in cubic feet per second (ft3s-1). source - Whether the right is for surface water or groundwater. origWaterUse - The original stated water use from the state. waterUse - The water use category under the unified use categories established here. rightsByState_surface: A list of the cleaned surface water rights collected from each state. See rightsByState_ground for more details on how the data is formatted. fullRightsRecs: A list of the combined groundwater and surface water records for each state. Essentially, these lists are the merging of rightsByState_ground and rightsByState_surface by state. See rightsByState_ground for more details on how the data is formatted. projProj: The spatial projection used for map creation in the beginning of the project. Specifically, the World Geodetic System (WGS84) as a coordinate reference system (CRS) string in PROJ.4 format. wmaStateLabel: The name and/or abbreviation for what each state legally calls their WMAs. h2oUseByState: A list of spatial polygon data frames which contain the area(s) in which each water right is claimed to be used. It should be noted that not all water right records have a listed area(s) of use in this object. Currently, only Idaho and Washington provided valid data to be included in this object. h2oDivByState: A list of spatial points data frames which identifies the location of the Point of Diversion for the state's water rights. It should be noted that not all water right records have a listed Point of Diversion in this object. spatialWMAByState: A list of spatial polygon data frames which contain the spatial WMA boundaries for each state. The only data contained within the table are identifiers for each polygon. It is worth reiterating that Arizona is the only state in which the surface and groundwater WMA boundaries are not the same. wmaIDByState: A list which contains the unique ID values of the WMAs for each state. plottingDim: A character vector used to inform mapping functions for internal map making. Each state is classified as either "tall" or "wide", to maximize space on a typical 8x11 page.

Economics↗

The Pursuit of K: Reflections on the Current State-of-the-Art in Stress Intensity Factor Solutions for Practical Aerospace Applications

The elastic stress intensity factor (SIF, commonly denoted as K) is the foundation of practical fracture mechanics (FM) analysis for aircraft structures. This single parameter describes the first-order effects of stress magnitude and distribution as well as the geometry of both structure/component and crack. Hence, the calculation of K is often the most significant step in fatigue analysis based on FM. This presentation will provide several reflections on the current state-of-the-art in SIF solution methods used for practical aerospace applications, including a brief historical perspective, descriptions of some recent and ongoing advances, and comments on some remaining challenges. Newman and Raju made significant early contributions to practical structural analysis by developing closed-form SIF equations for surface and corner cracks in simplified geometries, often based on empirical fits of finite element (FE) solutions. Those solutions (and others like them) were sometimes revised as new analyses were conducted or limitations discovered. The foundational solutions have exhibited striking longevity, despite the relatively "coarse" FE models employed many decades ago. However, in recent years, the accumulation of different generations of solutions for the same nominal geometry has led to some confusion (which solution is correct?), and steady increases in computational capabilities have facilitated the discovery of inaccuracies in some (not all!) of the legacy solutions. Some examples of problems and solutions are presented and discussed, including the challenge of maintaining consistency with legacy design applications. As computational power has increased, the prospect of calculating large numbers of SIF solutions for specific complex geometries with advanced numerical methods has grown more attractive. Fawaz and Andersson, for example, have been generating literally millions of new SIF solutions for different combinations of multiple cracks under simplified loading schemes using p-version FE methods. These data are invaluable, but questions remain about their practical use, because the tabular databases of key results needed to support practical life analysis can occupy gigabytes of storage for only a few classes of geometries. The prospect of using such advanced numerical methods to calculate in real time only those K solutions actually needed to support a specific crack growth analysis is also tempting, but the stark reality is that the computational cost is still so high that the approach is not practical except for specific, critical application problems. Some thoughts are offered about alternative paradigms. Compounding approaches are some of the earliest building blocks of SIF development for more complex geometries. These approaches are especially attractive because of their very low computational cost and their conceptual robustness; they are, in some ways, an intriguing contrast and complement to the brute-force numerical methods. In recent years, researchers at NRC-Canada have published remarkable results showing how compounding approaches can be used to generate accurate solutions for very difficult problems. Examples are provided of some successes--and some limitations--using this approach. These closed-form, tabulated numerical, and compounding approaches have typically been used for simple remote loading with simple load paths to the crack. However, many significant cracks occur in complex stress gradient fields. This is a job for weight function (WF) methods, where the arbitrary stress distribution on the crack plane in the corresponding uncracked body (typically determined using FE methods) is used to determine K. Several significant recent advances in WF methods and solutions are highlighted here. Fueled by advanced 3D numerical methods, many new solutions have been generated for classic geometries such as surface and corner cracks with wide ranges of geometrical validity. A new WF formulation has also be developed for part-through cracks considering the arbitrary stress gradients in all directions in the crack plane (so-called bivariant solutions). Basic WF methods have recently been combined with analytical expressions for crack plane stresses to develop a large family of accurate SIF solutions for corner, surface, and through cracks at internal or external notches with very wide ranges of shapes, sizes, acuities, and offsets. Finally, WF solutions are much faster than FE or boundary element solutions, but can still be much slower than simple closed-form solutions, especially for bivariant solutions that can require 2D numerical integration. Novel pre-integration and dynamic tabular methods have been developed that substantially increase the speed of these advanced WF solutions. The practical utility of advanced SIF methods, including both WF and direct numerical methods, is greatly enhanced if the FM life analysis can be directly and efficiently linked with digital models of the actual structure or component (e.g., FE models for stress analysis). Two recent advances of this type will be described. One approach directly interfaces the FM life analysis with the FE model of the uncracked component (including stress results). Through a powerful graphical user interface, simplified FM life models can be constructed (and visualized) directly on the component model, with the computer collecting the geometry and stress gradient information needed for the life calculation. An even more powerful paradigm uses expert logic to automatically build an optimum simple fracture model at any and every desired location in the component model, perform the life calculation, and even generate fatigue crack growth life contour maps, all with minimal user intervention. This paradigm has also been extended to the automatic calculation of fracture risk, considering uncertainty or variability in key input parameters such as initial crack size or location. Another new integrated approach links the engineering life analysis, the component model, and a 3D numerical fracture analysis built with the same component model to generate a table of SIF values at a specific location that can then be employed efficiently to perform the life calculation. Some attention must be given to verification and validation (V&V) issues and challenges: how good are these SIF solutions, how good is good enough, and does anyone believe the life answer? It is important to think critically about the different sources of error or uncertainty and to perform V&V in a hierarchal, building-block manner. Some accuracy issues for SIF solutions, for example, may actually involve independent material behavior issues, such as constraint loss effects for crack fronts near component surfaces, and can be a source of confusion. Recommendations are proposed for improved V&V approaches. This presentation will briefly but critically survey the range of issues and advances mentioned above, with a particular view towards assembling an integrated approach that combines different methods to create practical tools for real-world design and analysis problems. Examples will be selectively drawn from the recent literature, from recent enhancements in the NASGRO and DARWIN computer codes, and from previously unpublished research

CraigMcClung, R.↗

Enhancement of Optical Efficiency of CSP Mirrors for Reducing O&M Cost via Near-Continuous Operation of Self-Cleaning Electrodynamic Screens (EDS). Final Report

Over the past decade, techno-economical advancements in solar energy systems, particularly the improvement in the conversion efficiency of PV modules made with mono-crystalline silicon solar cells from 12% to 20%, as well as the cost reduction in manufacturing by a factor of about 10%, have made it possible to achieve a levelized cost of electricity (LCOE) in PV plants that is comparable to, or less than, the cost of deriving electricity from fossil fuels. Operating PV plants in mid-latitude sun-belt regions, where the solar irradiance level is highest, provides a high annual energy-yield (kWh/kWp) due to two factors: (1) the availability of predictable high solar irradiance throughout the year with the fewest interruptions in solar flux from clouds and rain, and (2) the increased conversion efficiency of crystalline solar cells, as recombination loss has decreased with increased intensity of the sunlight that illuminates the silicon solar cells. Semi-arid and desert regions, however, are plagued by high atmospheric dust concentrations and frequent sand storms. The deposition of a layer of dust on the optical surfaces of solar collectors such as PV modules and concentrating mirrors reduces the transmission efficiency of sunlight that actually reaches the solar cells or receivers, resulting in high energy-yield soiling loss. There are two major cost components to operating a solar plant: (1) installation costs, and (2) operation-and-maintenance (O&M) costs. There is no fuel cost; hence operating a solar plant in a semi-arid or desert region provides high returns on investment if soiling losses are mitigated via efficient cleaning methods and optimized cleaning frequency. If solar collectors are not cleaned, the accumulation of dust layers on solar collectors may cause the operation of such plants in arid regions to become economically unviable. Washing solar collectors with water and detergent, as is most commonly done now, is an efficient method for cleaning. The conventional approach in utility-scale solar plants is to use a large truck with a water tank and pump system for spraying deionized water on the surface of the solar collectors. Robotic cleaning with brushes, used for many solar plants, requires lesser water for cleaning. The water consumed using semi-automated cleaning of PV modules in utility-scale solar plants is approximately 2 liters/m 2 per cleaning cycle. The total optical surface area of the solar collectors in 1 TW-scale solar installation will be more than 3 × 109 m 2 ; hence an enormous amount of water be needed for cleaning. There simply is not enough fresh water in the sun-belt areas of the world for predicted cleaning needs. In solar power plants, the estimated cost of cleaning solar collectors includes expenses related to the equipment used, labor, cost of transportation of water, the energy required for cleaning, plus ancillary costs whereas the cost of water is not considered. The water used is obtained from sources located close to plant sites, unmindful of the environmental and societal impact as the power plants are oftentimes located in regions that face severe drought. This practice is very similar to the cost calculations in deriving the levelized cost of electricity (LCOE) in conventional power plants based on burning fossil fuel such as coal and gas, while disregarding the cost of climate change and health effects. Unless a water-free or low-water cleaning method is established, the expansion of solar plants may lose public support in areas suffering long intervals of drought. The goal of this research project has been the development and application of the Electrodynamic Screen (EDS) as a means for a water-free, scalable cleaning process applicable to solar-power installations, including rooftop applications. We describe here the development of an EDS film-based cleaning process as an emerging method for use on PV modules, parabolic troughs, and heliostats. This report aims to show the feasibility of integrating or retrofitting EDS films onto the optical surfaces of solar collectors (both PV and CSP) while maintaining high transmission or reflection efficiency. The cleaning action provided by the EDS film is an active method to remove dust deposits by electrodynamic force. Current lab-scale prototype EDS films, retrofitted onto solar panels and mirrors, have shown to be capable of maintaining optical transmission or specular-reflection efficiencies higher than 90% of initial values under clean conditions. The optical surfaces of solar collectors laminated with EDS films can remove more than 90% of deposited dust when the EDS is activated for less than two minutes. As an electrodynamic dust removal process, the EDS film-based method is designed primarily for the removal of dust in solar installations located in semi-arid and desert areas, where the atmosphere is often dry and dusty and rainfall is infrequent. While EDS film application minimizes water consumption and facilitates cleaning as frequently as needed, it has limitations in removing contaminants such as soot, organic pollutants deposited as fine films on the surface, and bird droppings. The dust removal efficiency of the EDS is maximum at relative humidity RH is 40 to 50% and decreases at when RH > 65%. Many solar plant sites undergo diurnal and seasonal cycles of high ambient, early afternoon temperatures and an RH that reaches the dew point early in the morning. These variations in atmospheric conditions do not limit the operation of the EDS film. This report presents a brief review of the progress and the potential of EDS film technology for mitigating the impact of dust on solar collectors via water-free cleaning, as well as current technical challenges regarding efficiency and durability. Our experimental data on the performance of EDS films show that: (1) the dust-removal efficiency (DRE) can reach levels higher than 90%, (2) the specular reflectivity (SR) of EDS film-laminated second-surface mirrors reach levels in excess of 90%, (3) the specular reflectivity restoration (SRR) can exceed 90%, (4) the output-power restoration (OPR) of PV modules can exceed 95%, and (5) the optical transmission efficiency (TE) of the EDS films can be greater than 90%. Working with Sandia National Laboratories (NM), Corning Research and Development Corporation (NY), Eastman Kodak (NY), Tomark-Worthen Industries (NH), and EDS Chile SPA (Chile), we have produced EDS film-laminated PV modules and demonstrated their self-cleaning functions without requiring water. We have demonstrated that the operational range of EDS films will cover the expected ambient temperature of solar fields at RH cycling varying from 20 to 95% as long as the EDS films are activated in the RH range 20 to 50%. (Typical solar-field climates in deserts and semi-arid lands often reach near dew point in coastal areas.) Our experiments on the application of hydrophobic-fluorinated nanoparticle coatings on EDS film surfaces show that the EDS operational range can be extended to higher RH levels that approach the dew point. At Eastman Kodak, as one of our industrial partners, we were able to establish a process for manufacturing EDS films using flexographic printing of the electrodes onto transparent polymer films. This process utilizes an existing manufacturing line at Eastman Kodak that allows fabrication of medium-scale EDS films (26 cm × 30 cm). The manufacturing process has the capacity to produce EDS films at high production speeds. The medium-scale EDS films that have been printed at Kodak were evaluated in the lab at Boston University. The EDS films produced at Kodak were laminated at Tomark-Worthen using an industrial scale vacuum laminator to produce EDS film stacks, which can be affixed onto the optical surfaces of PV modules or concentrating mirrors. The EDS film stack consists of the EDS films that have Willow® Glass (WG) which has a thickness of 100 μm as the front surface. The WG sheets obtained from Corning Research and Development Corporation are customized in size and shape to cover the active area of the EDS films. The back surface of the EDS film stack is integrated onto the optical surface of the solar panel or mirror using optically clear adhesive (OCA) films or silicone adhesives. The OCA films (thickness 25 μm) are produced by 3M. The EDS film stacks have the architecture: WG/OCA/EDS Film/OCA/ over PV module or solar mirror. The power-supply units needed for activating the electrodes of the EDS film were designed and produced at Boston University. These power supply units provide three-phase, 1.2 kV voltage pulses at a very low current (micro-ampere) level and at a low frequency (≈ 5 Hz). The voltage pulses are applied to the electrodes in a sequence such that the train of pulses resembles a unidirectional traveling wave of electrical field on the surface of the EDS film. The dust particles on the surface become charged electrostatically and are levitated by the Coulomb force. The lateral sweeping action of the traveling electric field created by the three-phase voltages pulses then sweeps the dust off the surface. The energy consumed by the EDS electrodes is less than 0.2 Wh/m 2 /cleaning cycle, enabling energy-efficient restoration of output power (OPR) of PV module or specular reflectivity restoration (SRR) for solar mirrors. The EDS system consists of (1) an EDS film stack laminated onto the solar collectors, (2) connection of the EDS film to its power supply unit and (3) Interconnection of the power supply to PV modules or solar mirrors. Design, construction and assessment of field-testing units that have EDS stack laminated PV modules for evaluating the performance of the EDS films in solar fields is being carried out at BU. Our progress under this project, aimed at the advancement of EDS film technology, has reached DOE Technology Readiness Level (TRL) 6. Based on the extensive laboratory evaluations and limited field trials, as well as contacts with potential users, we believe that the technology has reached its commercial stage. We are conducting a cost analysis using the National Renewable Energy Laboratory (NREL) System Advisory Model (SAM) and are preparing for field trials of EDS films in different solar fields in the US, Chile, India and in the Middle East. A brief description of EDS film performance, construction and testing of the field-test unit, autonomous operation of the field-test unit for evaluating EDS performance in increasing energy yield, associated revenue savings, and water conservation is presented.

14 SOLAR ENERGY↗

Autonomous Ocean World Exploration: Advancement of a Virtual Testbed

The search for life (extinct or extant) and potentially habitable bodies in our solar system and beyond is one of the 12 priority science questions outlined in the National Acadamies’ 2022 decadal survey [5]. Extraterrestrial destinations containing liquid water present an opportunity to search for life as we know it, and in recent years an increasing number of such locations have been discovered within our solar system. Several Jovian moons—Europa, Ganymede, and Callisto [10]—and the Saturnian moons Enceladus [8] and Titan [9] are known or suspected to harbor massive subsurface oceans. Of these "ocean worlds", Europa is the focus of at least one planned NASA orbiter mission, Europa Clipper [4], and an early lander mission concept, the Europa Lander [2, 3]. Whereas most robotic missions to the Moon and Mars (e.g. orbiters, rovers, landers) to date have had ground controllers on Earth tightly involved in mission operations, missions to more distant worlds will require a high degree of onboard autonomy due to long communication lags and blackouts, harsh environments (radiation, cold), and more limited battery and hardware life. The past decade has seen great advances in both AI technologies and computing scalability and performance that offer promising solutions for spacecraft autonomy and motivate the software system and research programs described in this paper. The Ocean Worlds Autonomy Testbed for Exploration, Research, and Simulation (OceanWATERS) [1], which has been in development at the NASA Ames Research Center since 2018, is a virtual environment for testing lander autonomy solutions. It is built on the Robot Operating System (ROS), runs on consumer-grade Linux workstations, and was released as open source in 2020. OceanWATERS provides a physical and visual simulation of a prototypical lander in a Europa-like environment (Figure 1). The lander was modeled after requirements and specifications made in JPL’s Europa Lander Study of 2016 [3]. Simulated lander systems include stereo cameras and spotlights mounted on an antenna mast that pans and tilts, a 6 degrees of freedom (DoF) robotic arm with a force-torque sensor and two interchangeable end effectors, and a battery pack power system. The environment consists of multiple terrain models including a highly detailed model sourced from the FROST dataset [11], simulation of surrounding planetary bodies based on an ephemeris model, and lighting from the sun with associated surface illumination, reflectance, and shadows. Operations supported by OceanWATERS include panoramic and directed imaging of the environment and lander workspace, Cartesian and joint-level arm commanding, grinding of the terrain surface (e.g. digging a trench), and scooping of ground material (Figure 2) which can be discarded or collected as science samples in a receptacle that can be emptied (science operations themselves are not simulated). These operations are realized as ROS Actions and are complimented by a wide selection of telemetry that is continually produced by each lander subsystem. The power system model is driven by the open-source Generic Software Architecture for Prognostics (GSAP) [11] that predicts the battery’s remaining useful life and other characteristics. As a testbed for high-level autonomy, OceanWATERS provides an execution framework based on PLEXIL [12], an open-source plan specification language and execution engine developed largely at Ames. NASA's initial development of OceanWATERS, as well the Ocean Worlds Lander Autonomy Testbed (OWLAT) [6], a complimentary physical testbed developed at JPL, was the first step in a plan for realizing candidate onboard autonomy solutions for such planetary landers. In 2020 NASA solicited applications for its Autonomous Robotics Research for Ocean Worlds (ARROW) program, and in 2021 the similar Concepts for Ocean worlds Life Detection Technology (COLDTech) program. Collectively six research teams, based in universities and companies across the United States, were awarded grants to develop and demonstrate autonomy solutions on OceanWATERS and OWLAT. These 1–2-year projects have now finished or are nearing completion, and a wide variety of autonomy challenges in ocean world surface missions were addressed. Prototyped and demonstrated solutions have included autonomous discovery, response and adaptation to system faults and unexpected environmental events, world model synthesis through perception, plan synthesis using learned models, methods to optimize sample target selection and prioritize science data transmission, extension of PLEXIL for stochastic decision-making, and an integration of a model of JPL’s mission-ready COLDArm [7]. Technologies used in these projects include many forms of machine learning, causal reasoning, automated planning, Markov decision processes, formal methods, and other advanced techniques. A more detailed summary of the ARROW and COLDTech projects is given herein. OceanWATERS has had significant enhancements since its open-source release in 2020. Many of its new features were driven or shaped by feedback from the ARROW and COLDTech teams and requirements of their projects. In support of enabling autonomous adaptation to spacecraft faults (a specific capability solicited by both programs), a fault injection and detection framework was developed that supports a wide and growing range of fault types such as locked joints, image loss, and battery failures. The power system model was completed and integrated into the simulator, starting as a single-cell battery model and later upgraded to a multi-cell model with associated faults such as cell disconnection. Arm/terrain interaction was improved by adding a force-torque sensor and associated faults, and an analytic dig force model based on the Balovnev bucket force equations. Environment fidelity was increased by modeling terrain deformation resulting from digging and scooping; visual improvements were made in textures, lighting, and shadows. To facilitate interoperation with OWLAT, a unified command and telemetry interface between the testbeds was developed at the ROS level, along with a PLEXIL interface. The number of lander operations was greatly expanded (e.g. with Cartesian-based arm and antenna movement), and a framework was designed for users to build their own lander actions. A GUI for PLEXIL plan selection was created (Figure 3), and an expansive set of plans were added, such as those that illustrate patterns for fault handling. This paper provides a self-contained high-level description of OceanWATERS, focusing on more detailed coverage of the aforementioned enhancements. It provides a high-level summary of the projects undertaken by participants in the ARROW and COLDTech programs and how these efforts have helped shape OceanWATERS. Finally, potential future work and directions for the testbed are listed, as likely informed by the recent planetary science decadal survey [5].

K Michael Dalal↗