SEARCH · Engineering Papers

Results for “Location Of Events”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 937 records · Page 52

All sources NASA NTRS DOE OSTI

DOE OSTI2021

WRF-simulated low-level jets over Iowa: characterization and sensitivity studies

Abstract. Output from 6 months of high-resolution simulations with the Weather Research and Forecasting (WRF) model are analyzed to characterize local low-level jets (LLJs) over Iowa for winter and spring in the contemporary climate. Low-level jets affect rotor plane aerodynamic loading, turbine structural loading and turbine performance, and thus accurate characterization and identification are pertinent. Analyses using a detection algorithm wherein the wind speed above and below the jet maximum must be below 80 % of the jet wind speed within a vertical window of approximately 20–530 m a.g.l. (above ground level) indicate the presence of an LLJ in at least one of the 14 700 4 km×4 km grid cells over Iowa on 98 % of nights. Nocturnal LLJs are most frequently associated with stable stratification and low turbulent kinetic energy (TKE) and hence are more frequent during the winter months. The spatiotemporal mean LLJ maximum (jet core) wind speed is 9.55 m s−1, and the mean height is 182 m. Locations of high LLJ frequency and duration across the state are seasonally varying, with a mean duration of 3.5 h. The highest frequency occurs in the topographically complex northwest of the state in winter and in the flatter northeast of the state in spring. Sensitivity of LLJ characteristics to the (i) LLJ definition and (ii) vertical resolution at which the WRF output is sampled is examined. LLJ definitions commonly used in the literature are considered in the first sensitivity analysis. These sensitivity analyses indicate that LLJ characteristics are highly variable with definition. Use of different definitions identifies both different frequencies of LLJs and different LLJ events. Further, when the model output is down-sampled to lower vertical resolution, the mean jet core wind speed height decreases, but spatial distributions of regions of high frequency and duration are conserved. Implementation of a polynomial interpolation to extrapolate down-sampled output to full-resolution results in reduced sensitivity of LLJ characteristics to down-sampling.

17 WIND ENERGY↗

DOE OSTI2020

Theoretical analysis of energetic-ion-driven resistive interchange mode stabilization strategies using a Landau closure model

The aim of the present study is to perform a theoretical analysis of different strategies to stabilize energetic-ion-driven resistive interchange mode (EIC) in Large Helical Device (LHD) plasma. In this work, we use a reduced MHD for the thermal plasma coupled with a gyrofluid model for the energetic particles (EP) species. The hellically trapped EP component is introduced through a modification of the drift frequency to include their precessional drift. The stabilization trends of the 1/1 EIC observed experimentally with respect to the thermal plasma density and temperature are reproduced by the simulations, showing a reasonable agreement with the data. The LHD operation scenarios with stable 1/1 EIC are identified, leading to the stabilization of the 1/1 EIC if the thermal plasma density and temperature are above a given threshold. Additionally, the 1/1 EIC are also stabilized if the rotational transform is modified in a way that the 1/1 rational surface is located further away than 0.9 times the normalized radius, or the magnetic shear in the plasma periphery is enhanced. Also, LHD discharges with large magnetic fields show a higher EIC destabilization threshold with respect to the thermal plasma density. If the perpendicular NBI deposition region is moved further inward than 0.875 times the normalized radius the 1/1 EIC are also stabilized. Furthermore, increasing the perpendicular NBI voltage such that the EP energy is higher than 30 keV stabilizes the 1/1 EIC. Moreover, deuterium plasmas show a higher stability threshold for the 1/1 EIC than hydrogen plasmas. The experimental data shows a larger time interval between EIC events as the power of the tangential NBI is increased providing that the perpendicular NBI power is at least 13 MW. This implies a stabilizing effect of the tangential NBI.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Measurement of engine-oil fuel dilution using laser induced florescence spectroscopy

An optical diagnostic, based on laser-induced fluorescence (LIF), has been developed for on-engine measurements of real-time fuel dilution of engine oil, or fuel in oil (FiO). Fuel dilution of oil is believed to be one of the underlying mechanisms that promotes low speed pre-ignition (LSPI) during high-load engine operations. While standard (e.g., ASTM D3524-90) methods are not capable of transient measurements, the LIF technique resolves transient dilution on the minutes time scale, making it possible to study the relation between FiO and LSPI. The measurement is based on adding a dye (ca. 500 ppm) to the fuel and monitoring for its presence in the oil via LIF. The instrument contains a 532-nm laser diode, and a spectrometer to resolve the LIF spectra. A fiber-optic assembly links the instrument to the engine measurement location; the probe is in a Y-configuration with the six fibers around one on a common probe end that connects to the engine; the Y-assembly legs are connected to the instrument; the two legs are the six fibers for transporting the laser light, and the single fiber for the fluorescence. We have expanded on our previous work (Parks et al, SAE paper 2007-01-4108, 2007), by demonstrating the diagnostic technique with to measure transient oil dilution on a 2.0-liter GM ECOTEC LNF single cylinder SI engine. The proposed LIF diagnostic will be critical in studying oil-dilution correlation with LSPI events.

Neupane, Sneha↗

DOE OSTI2022

Use of Machine Learning on PMU Data for Transmission System Fault Analysis

Synchrophasor technology has been used for monitoring, control, and protection of bulk power system for over 10 years. Deployment of phasor measurement units (PMUs) in the USA power system has surpassed 3000 units installed in the transmission substations as stand-alone intelligent electronic devices (IEDs) or as a software add-on to other devices such as digital protective relays (DPRs) or digital fault recorders (DFRs). By now, thousands of terabytes of PMU data may have been captured and stored by various transmission system operators (TSOs) and independent system operators (ISOs). This creates an opportunity to deploy advanced machine learning (ML) techniques to detect and classify faults recorded by PMUs automatically to be used by the system operators for rapid, critical decision-making when manual analysis of the past or unfolding events is not feasible. In this paper we offer a brief background on how the automated fault analysis may be done using DPR and/or DFR data, and compare some of the legacy approaches to the new ML approaches in the context of the system-wide PMU recordings. We then offer insights from developing practical ML solutions that have been applied on field recordings captured by close to 450 PMUs from all three US interconnections (Western, Eastern and ERCOT) over two years (2016-2017). We identify and illustrate ML challenges we addressed: inaccurate data, data with scarce and temporally imprecise fault labels, data recorded by PMUs sparsely located at substations resulting in the fault records taken afar from the ends of the faulted lines, data containing only positive sequence values, and data taken at different voltage levels. We then illustrate the ML model results for fault analysis under different application scenarios. The novelty of this study is not only in the design, implementation, and performance analysis of the ML algorithms, but also in the use of advanced fault modelling and simulation approaches to improve the training results when developing supervised ML models for fault detection and classification. Extensive simulations of faults were conducted on a 14-bus power system to create a training dataset with over 1400 accurately labelled faults. This dataset was applied to enhance the accuracy of fault detection and classification of machine learning-based models trained with small number of labelled faults in large datasets recorded in the grid interconnections ranging from 5,000 to 70,000 buses.

Synchrophasors, Machine Learning, Fault Analysis, ↗

DOE OSTI2024

Deep desert aquifers as an archive for Mid- to Late Pleistocene hydroclimate: An example from the southeastern Mediterranean

Many efforts have been made to illuminate the nature of past hydroclimates in semi-arid and arid regions, where current and future shifts in water availability have enormous consequences on human subsistence. Deep desert aquifers, where groundwater is stored for prolonged periods, might serve as a direct record of major paleo-recharge events. To date, groundwater-based paleoclimate reconstructions have mainly focused on a relatively narrow timescale (up to ∼40 kyr), limited by the relatively short half-life of the widely used radiocarbon (5.73 kyr). Here we demonstrate the usage of deep regional aquifers in the arid southeastern Mediterranean as a hydroclimate archive for earlier Mid-to-Late Pleistocene epochs. State-of-the-art dating tools, primarily the 81 Kr radioisotope (t 1/2 = 229 kyr), were combined with other atmosphere-derived tracers to illuminate the impact of four distinguishable wetter episodes over the past 400 kyr, with differences in climatic conditions and paleo-recharge locations. Variations in stable water isotope composition suggest moisture transport from more proximal (Mediterranean) and distal (Atlantic) sources to different parts of the region at distinct times. Large variability in the computed noble gas-based recharge temperature (NGT), ranging ~15–30 °C, cannot be explained by climate variations solely, and points to different recharge pathways, including geothermal heating in the deep unsaturated zone and recharge from high-elevation (colder) regions. The obtained groundwater record complements and enhances the interpretation of other terrestrial archives in the arid region, including a contribution of valuable information regarding the moisture source origin as reflected in the deuterium-excess values, which is unattainable from the common practice analysis of calcitic cave deposits. We conclude that similar applications in other deep (hundred-m-order) regional groundwater systems (e.g., the Sahara desert aquifers) can significantly advance our understanding of long-term (up to 1 Myr) paleo-hydroclimate in arid regions, including places where no terrestrial remnants, such as cave, lake, and spring sediments, are available.

Atom Trap Trace Analysis↗

DOE OSTI2024

Offshore low-level jet observations and model representation using lidar buoy data off the California coast

Abstract. Low-level jets (LLJs) occur under a variety of atmospheric conditions and influence the available wind resource for wind energy projects. In 2020, lidar-mounted buoys owned by the US Department of Energy (DOE) were deployed off the California coast in two wind energy lease areas administered by the Bureau of Ocean Energy Management: Humboldt and Morro Bay. The wind profile observations from the lidars and collocated near-surface meteorological stations (4–240 m) provide valuable year-long analyses of offshore LLJ characteristics at heights relevant to wind turbines. At Humboldt, LLJs were associated with flow reversals and north-northeasterly winds, directions that are more aligned with terrain influences than the predominant northerly flow. At Morro Bay, coastal LLJs were observed primarily during northerly flow as opposed to the predominant north-northwesterly flow. LLJs were observed more frequently in colder seasons within the lowest 250 m a.s.l. (above sea level), in contrast with the summertime occurrence of the higher-altitude California coastal jet influenced by the North Pacific High, which typically occurs at heights of 300–400 m. The lidar buoy observations also validate LLJ representation in atmospheric models that estimate potential energy yield of offshore wind farms. The European Centre for Medium-Range Weather Forecasts Reanalysis version 5 (ERA5) was unsuccessful at identifying all observed LLJs at both buoy locations within the lowest 200 m. An extension of the National Renewable Energy Laboratory (NREL) 20-year wind resource dataset for the Outer Continental Shelf off the coast of California (CA20-Ext) yielded marginally greater captures of observed LLJs using the Mellor–Yamada–Nakanishi–Niino (MYNN) planetary boundary layer (PBL) scheme than the 2023 National Offshore Wind dataset (NOW-23), which uses the Yonsei University (YSU) scheme. However, CA20-Ext also produced the most LLJ false alarms, which are instances when a model identified an LLJ but no LLJ was observed. CA20-Ext and NOW-23 exhibited a tendency to overestimate the duration of LLJ events and underestimate LLJ core heights.

17 WIND ENERGY↗

DOE OSTI2024

Spatial variability in Arctic–boreal fire regimes influenced by environmental and human factors

Abstract Wildfire activity in Arctic and boreal regions is rapidly increasing, with severe consequences for climate and human health. Regional long-term variations in fire frequency and intensity characterize fire regimes. The spatial variability in Arctic–boreal fire regimes and their environmental and anthropogenic drivers, however, remain poorly understood. Here we present a fire tracking system to map the sub-daily evolution of all circumpolar Arctic–boreal fires between 2012 and 2023 using 375 m Visible Infrared Imaging Radiometer Suite active fire detections and the resulting dataset of the ignition time, location, size, duration, spread and intensity of individual fires. We use this dataset to classify the Arctic–boreal biomes into seven distinct ‘pyroregions’ with unique climatic and geographic environments. We find that these pyroregions exhibit varying responses to environmental drivers, with boreal North America, eastern Siberia and northern tundra regions showing the highest sensitivity to climate and lightning density. In addition, anthropogenic factors play an important role in influencing fire number and size, interacting with other factors. Understanding the spatial variability of fire regimes and its interconnected drivers in the Arctic–boreal domain is important for improving future predictions of fire activity and identifying areas at risk for extreme events.

Geology↗

DOE OSTI2023

Diffusion–convection model of runaway electrons due to large magnetohydrodynamic perturbations in post-thermal quench plasmas

Systematic test particle tracing simulations for runaway electrons (REs) are performed for six post-thermal quench equilibria from DIII-D and ITER, where large scale, kink-like n = 1 (n is the toroidal mode number) magnetohydrodynamic (MHD) instabilities are found. The modeled particle guiding center orbits allow extraction of the effective diffusion–convection coefficients of REs in the presence of large three-dimensional (3D) perturbations up to 10% of the equilibrium toroidal field. With a fixed spatial distribution of the field perturbation, the RE transport coefficients along the plasma radial coordinate track reasonably well with the surface-averaged perturbation level. A substantial variation in the value of the transport coefficients—by three orders of magnitude in most cases, however, occurs with varying launching location of REs along the plasma radius. Large 3D perturbations almost always lead to comparable diffusion and convection processes, meaning that diffusion alone is insufficient to describe the particle motion. At lower (but still high) level of perturbation, the RE convection is found to be dominant over diffusion. A similar observation is made when the perturbation is too strong. In the presence of large perturbation, the dependence of the RE transport on the particle energy is sensitive to the spatial distribution of the perturbation. Based on numerically obtained RE transport coefficients, an analytic fitting model is proposed to quantify the particle diffusion and convection processes due to large MHD events in post-thermal quench plasmas. The model is shown to reasonably well reproduce the direct test particle tracing results for the RE loss fraction and can, thus, be useful for incorporating into other kinetic RE codes in order to simulate the RE beam evolution in the presence of large 3D perturbations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Monitoring natural gas storage using Synthetic Aperture Radar: are the residuals informative?

Estimates of line-of-sight (LOS) displacements from Interferometric Synthetic Aperture Radar (InSAR) observations serve as the basis of the long-term monitoring of an operating natural gas storage site at Honor Rancho in California. Here, an inversion algorithm is used to estimate the portion of the signal that is attributable to deformation within the gas storage reservoir, located at a depth of around 3 km. Removing this contribution produces residuals that are used to characterize the background variation is surface deformation at the gas storage facility and to determine a threshold that can signify unusually large residuals. An application to almost 7 yr of InSAR data, from 2011 until 2018, indicates that there are intervals of heightened residuals as well as brief episodes of anomalously large misfits. An examination of the spatial distributions of the individual residual LOS displacements indicates larger displacements in an alluvial valley just south of the reservoir, with rapid spatial variations in sign, indicating a rather shallow origin. Furthermore, the two anomalous events also involve rapid spatial variations in the LOS displacement residuals directly above the storage facility. The results demonstrate that the technique of extracting residuals after removing the reservoir signal is a useful approach, even in the case of this deep reservoir, and is a promising method for long-term monitoring.

03 NATURAL GAS↗

DOE OSTI2022

Scaling Sustainable Biochar Research & Commercialization for Agriculture & Conservation: A Summary from a Stakeholder Convening

Biochar is a carbon-rich substance with the potential to increase soil carbon sequestration, reduce greenhouse gas (GHG) emissions and improve soil health. When used as part of soil health management systems, biochar can provide multiple synergistic benefits for the agricultural sector and society, while also supporting the production of renewable, climate-neutral biofuels. Establishing a pyrolysis biochar bioenergy industry (PBBI) is necessary to harness the environmental benefits of biochar. To galvanize the necessary coordination to support a PBBI, the Foundation for Food & Agriculture Research (FFAR), the National Center for Appropriate Technology (NCAT), and American Farmland Trust (AFT) co-hosted a virtual convening event on biochar research and commercialization in March 2022. Participants represented a diverse set of stakeholders, including biochar producers, agricultural producers, and members of nonprofit, industry and government agencies. This whitepaper offers a summary of the convening, as well as insights from subsequent stakeholder engagement. Sustainably sourced, fit-for-purpose biochar can be a powerful tool in the soil health management systems and climate-smart toolbox, with decades of research supporting the soil, environment, production, and climate benefits. However, because biochar represents a range of feedstocks, processes and products, there is a critical need to: • Coordinate efforts to characterize and test the application of biochar types across locations with diverse soil and management conditions • Address research and decision-support gaps • Support the growing industry These objectives can be addressed or advanced through a coordinated strategy. To support establishment of this strategy, the convening participants identified key gaps in implementation and research. The experiences of foresters, ranchers, and farmers already using biochar, some of whom participated in the convening, helped inform priorities. Key gaps include: • Practical, applied on-farm research • Supporting management for agricultural productivity • Maximum benefits for conservation • Decision support • Life Cycle Assessment (LCA) • Quantification of climate mitigation potential Existing research gaps include: • Potential synergies among biochar, soil organic matter, and enhanced rock weathering • Development of new storage reservoirs for biochar in industrial products Convening participants further provided next steps for rapidly developing a PBBI to produce sustainable, fit-for-purpose biochars that increase soil health and mitigate climate change. Participants also stressed the need for actionable, coordinated, large-scale research relevant to commercial production in the next five years. To have a chance to reduce greenhouse gas emissions, we need an array of tools that would aid in mitigation; a coordinated strategy to enable biochar to be an effective tool in the climate-smart toolbox must be created; We cannot wait 50 years.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Reversed magnetic shear scenario development in NSTX-U using TRANSP

Abstract Understanding and control of electron thermal transport is a critical point of research in magnetic fusion experiments. Previous experiments have shown that operation with reversed magnetic shear (RMS) can suppress electron thermal transport, resulting in the generation of internal transport barriers (ITBs), with the location of the ITB correlated with the location of minimum magnetic shear. The recent upgrades to NSTX—increased magnetic field up to 1 T, increased plasma current up to 2 MA, 2nd neutral beam—present an increased operating space in which to explore electron thermal transport in RMS plasmas. Utilizing TRANSP, we have developed operating scenarios by which to generate RMS in NSTX-U. The results suggest that RMS in NSTX-U can be generated through fast current ramp and early beam injection into a large plasma volume. This is very similar to the procedure that was followed in both TFTR and NSTX to generate RMS. Sustainment of RMS, disregarding non-( q min = 1) MHD events, requires maintaining a large plasma volume, and increasing the core T e , either via increased plasma current and/or adding heating power. Using this procedure, RMS was sustained for ∼1 s, with q min > 1 for that period.

Galante, M. E. (ORCID:0009000098149425)↗

DOE OSTI2020

Cybersecurity for Distance Relay Protection

This project is a DOE follow-up effort on the CREDC workshop held on September 13, 2018 in Cambridge, MA to discuss cybersecurity of distance relays, which considered the benefits, vulnerabilities and risk mitigations for the use of communication systems in power system protection. The objectives of this project are to define the taxonomy of relay protection and associated communications; define use cases describing approaches to reduce the cyber-attack surface on those protective relays; and evaluate the loss of operational functional capability from changes to communication coverage. Mitigating controls will also be evaluated to understand if there are other approaches to reduce attack surfaces while maintaining communications or partial communications. Distance relays are used to protect transmission lines of approximately 10 to 300 miles in length, by detecting short circuits (i.e., faults) on the lines and then tripping circuit breakers in the substation. Such protection systems are a subset of the power system and they incorporate sensing, logic and communication functions. Protection system exposure to cyberattack could be drastically limited by disconnecting relays from all vulnerable communication systems, but this may adversely impact overall power system performance in the absence of cyberattack. This project began with a use case analysis of protection systems with communications, as summarized in this report. It continued with modeling, testing and evaluation in a miniature power system (MPS), located in the Western Area Power Administration (WAPA) Electric Power Training Center (EPTC). The project also incorporated feedback from two industry meetings held in February and September 2019. The suggested next steps account for and complement the work already underway with DOE/CESER funding: 1. Study the performance of LCD and PC vs. PUTT, which is less reliant on communication system performance and GPS timing references. The PUTT scheme could prove to be more resilient to cyberattack or communications-related disruption. It could also be more tolerant of message re-routing with SDN/SDR communication systems. On the other hand, it will be more vulnerable to false tripping during dynamic events or to loss of the voltage signal. The optimum choice of scheme may depend on the specific power system and risk assessment. This study could provide a new template for evaluation based on business functions. 2. Research and develop new methods to detect and monitor distributed physical attacks, possibly using drones, video sensors, thermal sensors, machine learning and other advanced techniques. This will help mitigate the impact of cyberattack on the protection system, and will also help mitigate the impact of wild fires. 3. Implement a scalable PKI for use in electric utility protection systems. This will encourage widespread adoption of secure authentication methods that are already available, but not widely used at present. This will help secure engineering access to the relays. 4. Investigate the use of SDN in combination with SDR to achieve better cybersecurity and electromagnetic security of the network, incorporating path variability. This would help secure both engineering access and peer-to-peer GOOSE messaging. 5. Perform additional testing, with operator evaluation of “red button” scenarios, PUTT vs. LCD, relay mis-operations, and other cyberattacks in the EPTC. This is an important advantage of testing in the EPTC rather than by computer simulation or even hardware-in-the-loop simulation; the EPTC is already dedicated to managing the situational awareness, operator response times and other human impacts. One of the project objectives was to settle on a common nomenclature for this problem space. We have concluded that the OSI layer model, supplemented by ANSI device numbers and other IEEE standards, is already well-accepted by the industry. The IEEE PSRC knowledge base provides a great deal of public information

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2025

The Short-Baseline Near Detector at Fermilab: Input to the European Strategy for Particle Physics 2026 Update

SBND is a 112 ton liquid argon time projection chamber (LArTPC) neutrino detector located 110 meters from the Booster Neutrino Beam (BNB) target at Fermilab. Its main goals include searches for eV-scale sterile neutrinos as part of the Short-Baseline Neutrino (SBN) program, other searches for physics beyond the Standard Model, and precision studies of neutrino-argon interactions. In addition, SBND is providing a platform for LArTPC neutrino detector technology development and is an excellent training ground for the international group of scientists and engineers working towards the upcoming flagship Deep Underground Neutrino Experiment (DUNE). SBND began operation in July 2024, and started collecting stable neutrino beam data in December 2024 with an unprecedented rate of ~7,000 neutrino events per day. During its currently approved operation plans (2024-2027), SBND is expected to accumulate nearly 10 million neutrino interactions. The near detector dataset will be instrumental in testing the sterile neutrino hypothesis with unprecedented sensitivity in SBN and in probing signals of beyond the Standard Model physics. It will also be used to significantly advance our understanding of the physics of neutrino-argon interactions ahead of DUNE. After the planned accelerator restart at Fermilab (2029+), opportunities are being explored to operate SBND in antineutrino mode in order to address the scarcity of antineutrino-argon scattering data, or in a dedicated beam-dump mode to significantly enhance sensitivity to searches for new physics. SBND is an international effort, with approximately 40% of institutions from Europe, contributing to detector construction, commissioning, software development, and data analysis. Continued European involvement and leadership are essential during SBND's operations and analysis phase for both the success of SBND, SBN and its role leading up to DUNE.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2024

DECOVALEX-2023: Task F1 Final Report

DECOVALEX-2023 Task F is a comparison of models and methods for post-closure performance assessment (PA) of a deep geologic repository for radioactive waste. The general aims of Task F are to build confidence in the models, methods, and software used for PA and to stimulate additional research and development in PA methodologies. The task objectives are to motivate development of PA modelling skills and capabilities, to examine the influence of model choices on calculated repository performance, and to compare the uncertainties introduced by model choices to other sources of uncertainty. Task F involves no actual experiment or site. It is a PA modelling exercise that requires the conceptual development of hypothetical repository designs and geologic settings. Because three of the teams were interested in salt and the rest of the teams were interested in crystalline rock, Task F was split into two branches: Task F1 for crystalline rock and Task F2 for salt. This report is for Task F1, crystalline rock. Teams from seven countries (Canada, Czech Republic, Germany, Korea, Sweden, Taiwan, and United States) participated in Task F1. The teams worked together to define the features, events, and processes of the reference case repository and established a set of performance measures. In addition, they defined a set of benchmark problems designed to test and compare modelling capabilities for fracture flow and transport at different scales. The repository design and benchmark problems are documented in a Task Specification that evolved over time as the group honed the specifications. The benchmark problems verified that each team can aptly model flow and transport in fractured media in 1-, 2-, and 3-dimensions. Two general approaches were used for the 3-dimensional benchmarks: discrete fracture network (DFN) and equivalent continuous porous medium (ECPM). DFN modelling involves explicit meshing of each fracture while ECPM modelling aims to capture the effective porosity and directional permeability of each cell in a space-filling mesh as affected by intersecting fractures. In some models, a combination of the two is used, i.e., DFN for large known fractures and ECPM for the rest of the domain. Transport is solved by using either the advection-dispersion equation or particle tracking. Although some variation is observed among model breakthrough curves in the benchmark problems, there is strong agreement in breakthrough behaviour up to at least the 75 th percentile for all benchmarks. At the 90 th percentile, breakthrough results show larger differences, suggesting several models retain substantially higher fractions of tracer in regions of slower moving water. In addition to the flow and transport benchmarks, several teams completed the source term benchmark, verifying capabilities for modelling radionuclide decay and ingrowth, waste package breach, instant release fractions, fuel matrix degradation rates, and radionuclide solubility limitations. The reference case is conceptualized as a generic spent fuel repository at a depth of 450 m in fractured crystalline rock. The repository has 50 parallel backfilled drifts, each with 50 deposition holes 6 m apart. Each deposition hole contains a 4-PWR waste package and bentonite buffer. The rock domain is 5 km in length, 2 km in width, and 1 km in depth. It has 6 deterministic fractured deformation zones and a multitude of stochastic fractures. Teams generally used the ECPM approach for the entire rock or a hybrid approach in which the deterministic fracture zones are modelled with a DFN and the rest of the rock is modelled by ECPM. Of the reference case problems specified, only the results of the initial reference case problem are compared in this report. The initial problem focuses on transport from the deposition holes to the surface, i.e., it neglects waste package performance. Tracers are released at all waste package locations at time zero and tracked for their releases to the near field and ground surface. The water fluxes calculated at the ground surface entry and exit regions of the domain are similar for all models except for two that have considerably lower fluxes. For tracer transport, large differences are observed among models in the magnitude of tracer transported. Much of the difference appears to be due to how the repository is implemented and hence the different degrees of repository simplification. Models that exclude the drifts, buffer, and backfill from the domain tend to show greater release of tracers and radionuclides from the repository. The initial study presented here indicates that major differences in modelling important processes within the repository (e.g., diffusion through buffer and backfill) can produce broadly different release and transport results, especially when those processes are excluded. Even for the models that included all specified features, events, and processes, the results show significant differences and demonstrate the importance of examining multiple modelling approaches in performance assessment. The differences in results observed in this study are expected to motivate teams to either increase complexity in future versions of the reference case models or to improve methods to account for the effects of simplified features and processes. Either way, future improvements in these models are expected to produce results that more closely agree.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2023

Crowd cluster data in the USA for analysis of human response to COVID-19 events and policies

We provide data on daily social contact intensity of clusters of people at different types of Points of Interest (POI) by zip code in Florida and California. This data is obtained by aggregating fine-scaled details of interactions of people at the spatial resolution of 10 m, which is then normalized as a social contact index. We also provide the distribution of cluster sizes and average time spent in a cluster by POI type. This data will help researchers perform fine-scaled, privacy-preserving analysis of human interaction patterns to understand the drivers of the COVID-19 epidemic spread and mitigation. Current mobility datasets either provide coarse-level metrics of social distancing, such as radius of gyration at the county or province level, or traffic at a finer scale, neither of which is a direct measure of contacts between people. We use anonymized, de-identified, and privacy-enhanced location-based services (LBS) data from opted-in cell phone apps, suitably reweighted to correct for geographic heterogeneities, and identify clusters of people at non-sensitive public areas to estimate fine-scaled contacts.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2020

A Multi-Function AAA Algorithm Applied to Frequency Dependent Line Modeling

Modeling of power distribution system components that are valid for a wide range of frequencies are crucial for highly accurate modeling of electromagnetic transient (EMT) events. This has recently become of interest due to the improvements needed for the resilient operation of distribution systems. Vector fitting (VF) is a very popular and commonly used algorithm for wide band representations of power system components in EMT simulations. In this research, we present a new multi-input rational approximation algorithm (MIAAA) and illustrate its advantages with respect to VF using examples of approximations of admittance matrices discussed in the literature. We show that MIAAA not only outperforms VF in terms of achieving better accuracy using lesser number of poles, but also has no numerical issues achieving convergence. In contrast to VF, MIAAA is not sensitive to the location of input sample points and it does not require good estimates for the location of the desired approximation poles. The novelty of this research work is the use of recent mathematical results to solve existing challenges in distribution system modeling and to develop rational approximations for power system models that intend to be optimal in terms of accuracy and performance.

41 EE - Solar Energy Technologies Office (EE-4S)↗

DOE OSTI2020

Developing a sustainability science approach for water systems

We convened a workshop to enable scientists who study water systems from both social science and physical science perspectives to develop a shared language. This shared language is necessary to bridge a divide between these disciplines’ different conceptual frameworks. As a result of this workshop, we argue that we should view socio-hydrological systems as structurally co-constituted of social, engineered, and natural elements and study the “characteristic management challenges” that emerge from this structure and reoccur across time, space, and socioeconomic contexts. This approach is in contrast to theories that view these systems as separately conceptualized natural and social domains connected by bi-directional feedbacks, as is prevalent in much of the water systems research arising from the physical sciences. A focus on emergent characteristic management challenges encourages us to go beyond searching for evidence of feedbacks and instead ask questions such as: What types of innovations have successfully been used to address these challenges? What structural components of the system affect its resilience to hydrological events and through what mechanisms? Are there differences between successful and unsuccessful strategies to solve one of the characteristic management challenges? If so, how are these differences affected by institutional structure and ecological and economic contexts? To answer these questions, social processes must now take center stage in the study and practice of water management. We also argue that water systems are an important class of coupled systems with relevance for sustainability science because they are particularly amenable to the kinds of systematic comparisons that allow knowledge to accumulate. Indeed, the characteristic management challenges we identify are few in number and recur over most of human history and in most geographical locations. This recurrence should allow us to accumulate knowledge to answer the above questions by studying the long historical record of institutional innovations to manage water systems.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

A MIMO AAA Agorithm for Frequency Dependent Line Modeling

AAA algorithm↗