SEARCH · Engineering Papers

Results for “Level Set Method”

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 901 records · Page 50

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Systematically tracking the hourly progression of large wildfires using GOES satellite observations

Abstract. In the western United States, prolonged drought, a warming climate, and historical fuel buildup have contributed to larger and more intense wildfires as well as to longer fire seasons. As these costly wildfires become more common, new tools and methods are essential for improving our understanding of the evolution of fires and how extreme weather conditions, including heat waves, windstorms, droughts, and varying levels of active-fire suppression, influence fire spread. Here, we develop the Geostationary Operational Environmental Satellites (GOES)-Observed Fire Event Representation (GOFER) algorithm to derive the hourly fire progression of large wildfires and create a product of hourly fire perimeters, active-fire lines, and fire spread rates. Using GOES-East and GOES-West geostationary satellite detections of active fires, we test the GOFER algorithm on 28 large wildfires in California from 2019 to 2021. The GOFER algorithm includes parameter optimizations for defining the burned-to-unburned boundary and correcting for the parallax effect from elevated terrain. We evaluate GOFER perimeters using 12 h data from the Visible Infrared Imaging Radiometer Suite (VIIRS)-derived Fire Event Data Suite (FEDS) and final fire perimeters from the California's Fire and Resource Assessment Program (FRAP). Although the GOES imagery used to derive GOFER has a coarser resolution (2 km at the Equator), the final fire perimeters from GOFER correspond reasonably well to those obtained from FRAP, with a mean Intersection-over-Union (IoU) of 0.77, in comparison to 0.83 between FEDS and FRAP; the IoU indicates the area of overlap over the area of the union relative to the reference perimeters, in which 0 is no agreement and 1 is perfect agreement. GOFER fills a key temporal gap present in other fire tracking products that rely on low-Earth-orbit imagery, where perimeters are available at intervals of 12 h or longer or at ad hoc intervals from aircraft overflights. This is particularly relevant when a fire spreads rapidly, such as at maximum hourly spread rates of over 5 km h−1. Our GOFER algorithm for deriving the hourly fire progression using GOES can be applied to large wildfires across North and South America and reveals considerable variability in the rates of fire spread on diurnal timescales. The resulting GOFER product has a broad set of potential applications, including the development of predictive models for fire spread and the improvement of atmospheric transport models for surface smoke estimates. The resulting GOFER product has a broad set of potential applications, including the development of predictive models for fire spread and the improvement of atmospheric transport models for surface smoke estimates (https://doi.org/10.5281/zenodo.8327264, Liu et al., 2023).

Geology↗

DOE OSTI2024

Fox Trails

1. This software utilizes python pandas to pull data from P6 databases or XER files. The software transforms the datasets into multiple main tables by joining, filtering, iteratively flattening hierarchical structured data, and pivoting datasets to give simple flat output tables. The activity table includes all of the information related to an activity including activity codes, global, EPS, and project codes, UDFs, and WBS information as separate columns. This includes the code id, code value and sequence number for all levels in hierarchical codes. The resource table is similar to the activity table and includes all of the information related to resources on activities including UPFs and resource codes. The resource time phased table takes the resource information and time phases it for the budget, forecast, late, and actual dates/units/costs that closely matches P6's user interface's values as it implements the resource curve and calendars. The wbs table contains the WBS structure broken out by levels and includes UDFs, codes, and notebook topics. The final P6 data table is the relationships table which simply contains the relationships. 2. When a user updates the tool with data (via giving it P6 project names with database username/password information or XER files) the system creates the data in #1, then creates a networkx graph with the activity data imbedded in the node data and the relationships added as edges. Each edge also has it's float calculated (working time distance between the predecessor and successor) and attached to the edge. Activities are also tagged as a potential start of a path based on their constraints, constraint dates, remaining start date, and activity status. When a user enters an activity ID into the UI, it runs a shortest path calculation on the network graph between each node tagged as potential start to the entered activity id based on the float tagged on the edge. Each path returned by the algorithm contains all of the nodes on the path in order, as well as the total float of the edges that make the path. This data is then collected and returned to the user in the form of a gantt chart with groupings for each path that includes the total float for each group. 3. Similar to 2, if the user passes through a reference dataset each activity set in the path is checked to see if it had a path in the reference dataset, if that path was the primary path between the start and end activities, and what has changed regarding logic and durations. These changes are color coded and summarized before sent to the user to be displayed by the UI for simple discovery. 4. Utilizing the data from #1, the user can submit desired grouping code(s) and filters to the system. The system will then pull the activities, resources, and relationships and create a gantt chart based on the groupings sent and filtered based on the filters sent. 5. The system will produce a gantt chart in a similar method to #4, but allows interactivity with the data. As the user interacts with the gantt chart, the software captures the changes and stores it with the user making the change so that project controls and implement those changes in P6.

Fox, Ben↗

DOE OSTI2024

Active Aerodynamic Load Control for Wind Turbines

The goal of this project was to develop and demonstrate an advanced dielectric barrier discharge (DBD) plasma actuator technology. We set out to demonstrate the efficacy and impact of the new actuator technology as a key component of an active load control system for wind turbines. A DBD plasma actuator consists of a thin layer of dielectric material separating a pair of offset electrodes. One of those electrodes is embedded between the dielectric layer and a non-conductive substrate, while the other is exposed to air. When driven by an appropriate high voltage waveform, the device ionizes the air near the surface of the dielectric and adjacent to the exposed electrode. Collisions between the ions in the plasma and neutral air molecules result in a wall-jet – a region of induced air velocity that can be used to modify the flow around a lifting body. Use of the device near the trailing edge of a wind turbine blade, designed in such a way as to amplify the effects of the flow-modifying device, can result in large changes to the global forces experienced by the blade. Because of the fast response time of the device, it can allow the turbine to react in real time to changes in the wind associated with turbulence, wind shear, gusts, and the like, when paired with appropriate sensors and control algorithms. The goal of this project was to increase the capacity of the device to induce velocity on its surface at the levels required by large, utility-scale wind turbines. The first technical goal was to increase the induced velocity from the current industry-best of about 4 m/s to 10 m/s by modifying the electrical waveform used to drive the device and by introducing a semi-conductive surface coating to control electrostatic charge build-up on the surface. The second goal was to use the device to modify the lift on a representative airfoil in a wind tunnel, demonstrating a reduction in lift coefficient of 0.2 or better. The third goal was to produce an actuator-induced velocity of 20 m/s. Alongside our partners at the University of Texas at Dallas, we also applied modern, advanced design methods to optimize the impact of the technology on the design of wind turbines. We also had planned to install a segmented, active load control system on a test turbine to demonstrate the ability to reduce unsteady aerodynamic forces on the turbine associated with changes in the wind. The highest induced velocity achieved during this research was 11 m/s. However, practical design constraints limited the change in lift coefficient to about 0.12 for a representative airfoil in the wind tunnel at a Reynolds number of 400,000. The primary conclusion is that the plasma actuator control authority remains insufficient for practical purposes when extrapolated to Reynolds numbers over 1 million. The potential impact of the active lift control concept was evaluated through detailed simulations for three wind turbine sizes: a 3.4MW onshore turbine, a 10MW offshore turbine, and a 15MW offshore turbine. A feedback control system was designed for each turbine within two scenarios: one where the active lift control is used as a retrofit capability on the baseline design, and the other where the designers were permitted to “upscale” the turbines in order increase annual energy production. The levelized cost of energy was then evaluated for the range of turbine sizes and design configuration. It was found that the LCOE reduction associated with active lift control fell in the range of 0.7% to 7.2%, with the highest reduction associated with upscaling the 3.4MW turbine.

17 WIND ENERGY↗

DOE OSTI2021

VERA-Grizzly Ex-Core Calculations: Watts Bar Unit 1 Cycles 1-2

The critical structures that comprise light-water reactor (LWR) nuclear power plants are subjected to operating environments that can challenge their integrity. Structures in close proximity to the reactor core, such as the reactor pressure vessel (RPV) and the biological shield wall, are subjected to high levels of radiation emanating from the core, as well as elevated temperatures. As the US fleet of operating LWRs ages, the effects of these operating environments on the integrity of these structures must be considered to ensure their continued safe operation. Extending the lifetime of commercial reactors and maintaining the aging reactor fleet require accurate prediction of the exposure of ex-core components to neutron and photon radiation. In particular, concrete degradation studies must be performed to evaluate the safety and long-term operation of reactors with lifetime extensions. The concrete reactor bioshield is important for providing radiological protection during operation and must last for the entire lifetime of the reactor. Recent interest in lifetime extensions furthers the need to accurately simulate concrete material degradation in the reactor bioshield. As a result of this need, the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program has funded this study to couple its tools, Virtual Environment for Reactor Applications (VERA) and Grizzly. VERA allows users to set up models to calculate time-dependent and fully coupled solutions (with thermal feedback) for ex-core quantities of interest such as vessel and coupon fluence and detector responses for multiple statepoints and cycles. Grizzly is a finite-element application based on the Multiphysics Object Oriented Simulation Environment (MOOSE) framework that is used to enable aging materials calculations. This report highlights the work performed to calculate the fluence in the vessel and concrete for Watts Bar Nuclear Plant Unit 1 (WBN1) Cycles 1 and 2. The fluences obtained from VERA were successfully transferred to Grizzly using a Python script. Four simulations were run with Grizzly: (1) the Mazars model with the initial Young’s modulus being the instantaneous modulus, (2) the Mazars model with the initial Young’s modulus being the delayed modulus, (3) the Mazars model with the initial Young’s modulus being the delayed modulus with the addition of the effects of micro-damage caused by irradiation, and (4) the Mazars model with the initial Young’s modulus being the instantaneous modulus, and with the addition of micro-damage and creep. Details regarding the methods used to obtain the fluence and the statistical errors associated with the VERA Monte Carlo Shift calculations are discussed in greater detail in this report. The results obtained from the four Grizzly models are also presented in this report.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2023

Light Water Sustainability Program: Optimizing Information Automation Using a New Method Based on System-Theoretic Process Analysis

This report describes the interim progress for research supporting the design and optimization of information automation systems for nuclear power plants. Much of the domestic nuclear fleet is currently focused on modernizing technologies and processes, including transitioning toward digitalization in the control room and elsewhere throughout the plant, along with a greater use of automation, artificial intelligence, robotics, and other emerging technologies. While there are significant opportunities to apply these technologies toward greater plant safety, efficiency, and overall cost-effectiveness, optimizing their design and avoiding potential safety and performance risks depends on ensuring that human-performance-related organizational and technical design issues are identified and addressed. This report describes modeling tools and techniques, based on sociotechnical system theory, to support these design goals and their application in the current research effort. The report is intended for senior nuclear energy stakeholders, including regulators, corporate management, and senior plant management. We have developed and employed a method to design an optimized information automation ecosystem (IAE) based on the systems-theoretic constructs underlying sociotechnical systems theory in general and the Systems-Theoretic Accident Modeling and Processes (STAMP) approach in particular. We argue that an IAE can be modeled as an interactive information control system whose behavior can be understood in terms of dynamic control and feedback relationships amongst the system’s technical and organizational components. Up to this point, we have employed a Causal Analysis based on STAMP (CAST) technique to examine a performance- and safety-related incident at an industry partner’s plant that involved the unintentional activation of an emergency diesel generator. This analysis provided insight into the behavior of the plant’s current information control structure within the context of a specific, significant event. Our ongoing analysis is focused on identifying near-term process improvements and longer-term design requirements for an optimized IAE system. The latter analyses will employ a second STAMP-derived technique, System-Theoretic Process Analysis (STPA). STPA is a useful modeling tool for generating and analyzing actual or potential information control structures. Finally, we have begun modeling plantwide organizational relationships and processes. Organizational system modeling will supplement our CAST and STPA findings and provide a basis for mapping out a plantwide information control architecture. CAST analysis findings indicate an important underlying contributor to the incident under investigation, and a significant risk to information automation system performance, was perceived schedule pressure, which exposed weaknesses in interdepartmental coordination between and within responsible plant organizations and challenged the resilience of established plant processes, until a human caused the initiating event. These findings are discussed in terms of their risk to overall system performance and their implications for information automation system resilience and brittleness. We present two preliminary information automation models. The proactive issue resolution model is a test case of an information automation concept with significant near-term potential for application and subsequent reduction in significant plant events. The IAE model is a more general representation of a broader, plantwide information automation system. From our results, we have generated a set of preliminary system-level requirements and safety constraints. These requirements will be further developed over the remainder of our project in collaboration with nuclear industry subject matter experts and specialists in the technical systems under consideration. Additionally, we will continue to pursue the system analyses initiated in the first part of our effort, with a particular emphasis on STPA as the main tool to identify weak or weakening control structures that affect the resilience of organizations and programs. Our intent is to broaden the scope of the analysis from an individual use case to a related set of use cases (e.g., maintenance tasks, compliance tasks) with similar human-system performance challenges. This will enable more generalized findings to refine the Proactive Issue Resolution and IAE models, as well as their system-level requirements and safety constraints. We will use organizational system modeling analyses to supplement STPA findings and model development. We conclude the report with a set of summary recommendations and an initial draft list of system-level requirements and safety constraints for optimized information automation systems.

99 GENERAL AND MISCELLANEOUS↗

DOE OSTI2023

Evaluating Impacts of the Inflation Reduction Act and Bipartisan Infrastructure Law on the U.S. Power System

The Inflation Reduction Act of 2022 (IRA) and the Infrastructure Investment and Jobs Act of 2021, commonly referred to as the 'Bipartisan Infrastructure Law (BIL),' collectively represent the largest commitment of the U.S. Federal Government to invest in the modernization and decarbonization of the U.S. energy system. The Congressional Budget Office (CBO) estimates that total support for the broad range of climate and clean energy programs, tax credits, and other incentives authorized through the two laws will exceed $430 billion from 2022 through 2031 (CRS 2022; CBO 2021, 2022). While the climate and clean energy provisions are numerous and have the potential to impact all aspects of the U.S. energy system from fuel and electricity production to final consumption in industry, transportation, and buildings, the provisions relevant to the electricity sector - in particular the suite of tax credits for clean generation, storage, and carbon dioxide ( CO 2 ) capture and storage - are expected to be some of the most consequential in terms of emissions reduction and clean energy deployment (Larsen et al. 2022; Jenkins, Mayfield, et al. 2022; Mahajan et al. 2022; Zhao et al. 2022). In this report, we detail the methods and results of a study estimating the potential impacts of key provisions of IRA and BIL on the contiguous U.S. power sector from present day through 2030. The analysis employs an advanced power system planning model, the Regional Energy Deployment System (ReEDS), to evaluate how major provisions from both laws impact investment in and operation of utility-scale generation, storage, and transmission, and, in turn, how those changes impact power system costs, emissions, and climate and health damages. While not exhaustive in capturing every provision, the analysis estimates the possible scale of power-sector impacts that could result from the modeled provisions in IRA and BIL. The study is structured around two scenarios to evaluate the potential impacts of both laws on the power sector: 1) No New Policy: A counter-factual scenario that reflects all Federal and state policies enacted as of September 2022, with exception to IRA and BIL, and assumes load growth consistent with the Energy Information Administration's Annual Energy Outlook 2022 (AEO22) Reference case (EIA 2022a); 2) IRA-BIL: A scenario reflecting all Federal and state policies enacted as of September 2022, including key IRA and BIL provisions, most notably the investment and production tax credits for zero-carbon emitting electricity generation and storage (ITC and PTC), the tax credit for CO 2 capture and storage (45Q), and the tax credit for existing nuclear plants (described further in Section 2.3). To account for the impacts of IRA and BIL on electrification, assumes increased load growth consistent with a scaled version of the Medium Electrification scenario from the Electrification Futures Study (Mai et al. 2018). These scenarios are simulated across seven sets of assumptions with varying projected future electricity market conditions, including technology costs and performance, natural gas prices, and the degree of availability, feasibility, and cost of development of renewable resources, electricity transmission, and CO 2 pipeline, injection, and storage infrastructure. In addition, we simulate two sensitivities on the 'policy' treatment in which we vary key assumptions pertaining to the realized value of the clean electricity ITC and PTC: 1) the cost of monetization of tax credits, and 2) the level of bonus crediting realized by project developers. We demonstrate that IRA and BIL have the collective potential to drive substantial growth in clean electricity by 2030, while reducing costs for consumers, mitigating climate change, and decreasing the human health impacts of power sector emissions. However, we also demonstrate that if expected cost improvements of clean technologies are not realized and/or constraints on deployment driven by factors such as supply-chain challenges, regulatory hurdles, and the social acceptability of energy infrastructure development limit the rate of clean energy and associated infrastructure deployment (such as transmission), then the share of clean generation achieved and the associated emissions benefits realized may be substantively reduced.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

DOE OSTI2026

From Modular ADMS to Plug-and-Play Ops: Distribution Grid Operations with Platform-Level Orchestration to Enable Ambitious App Hosting

The core function of the distribution grid is to provide electricity to consumers affordably, reliably, and securely. In pursuing these core objectives, distribution utilities are accountable to customers, regulators, and in some cases, shareholders. Other third parties such as aggregators and microgrids can also have a stake in the smooth operation of the grid. Each of these stakeholders has economic, business, and/or governance objectives that inform their expectations of the distribution grid. This multi-objective, multi-stakeholder environment creates tension that must be reconciled to successfully design and operate the distribution grid. Innovative companies are competing to bring high-tech solutions to electric utilities and their customers that address each of these objectives. Many developers of advanced distribution management systems (ADMS) and distributed energy resource management systems (DERMS) have adopted a modular architecture that allows grid operators to select functions and features according to their individual system needs. A modular platform also allows the solution provider to develop and integrate specific new product modules; however, the need to pursue multiple objectives with a fixed set of controllable devices makes integration expensive whether it is done at the product development stage or the deployment stage. This cost creates a significant barrier to adoption and can lengthen the product to market time of new solutions. To fundamentally address the complexity of system integration for distribution grid operations, the U.S. Department of Energy Office of Electricity has funded the GridAPPS-D project at PNNL, which streamlines integration by contributing to standards development, defining system architecture, applying advanced mathematics, and developing open-source software to demonstrate the concept of an open data-integration platform for distribution operations. The open data-integration platform concept enables system operators and solution providers to deploy ambitious, best-of-breed applications (or apps) without continually reengineering for integration. Ambitious apps developed by different solution providers will inevitably attempt to achieve different control objectives with the same set of controllable devices. If the open platform itself can resolve these conflicts in a way that achieves the best available outcomes for all apps, doesn’t restrict the ambitious design of apps, and ensures safe and secure operations, apps will be able to plug-and-play with the platform at the same time as other ambitious apps. In this paper, we describe a framework called App Deconfliction that empowers a platform to assign setpoints to controllable devices based on the values preferred by different apps (and even external stakeholder entities like customers or aggregators). The App Deconfliction framework is compatible with several methods for determining setpoint values. We present two methods based on game theory that provide a subtle built-in incentive structure for developers to adapt their apps to the fact that they will be operating in a moderated multi-app environment and to favor device setpoints that have the most effect on their objectives over those that have the least effect. Our simulation-based demonstrations have shown that game-theory-based deconfliction can lead to a 7% improvement in control space utilization compared to design-based methods.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2022

Findings on Subtask 3.1 - Bakken Rich Gas Enhanced Oil Recovery Project

Total in-place oil for the Bakken petroleum system (BPS) (which includes the Bakken and Three Forks Formations) has been estimated to be 600 billion barrels (bbl). However, BPS wells have decline rates as high as 85% over the first 3 years of their lives, and primary recovery factors typically range from 3% to 10% of original oil in place. Given the low initial recovery rates, even small incremental productivity improvements could dramatically increase technically recoverable oil in the BPS. One potential solution is enhanced oil recovery (EOR) using gas injection, such as carbon dioxide (CO2) or hydrocarbon (HC) gases. While commonly used in conventional reservoirs, CO2 EOR in unconventional tight oil reservoirs has been limited to pilot tests. EOR using rich gas (mixture of methane, ethane, and propane) has also been employed in numerous pilots in several unconventional plays and has recently been successfully applied in the Eagle Ford play. If successful, large-scale gas-based EOR in the BPS could dramatically increase oil productivity and recovery factors and extend the life of the play for decades. While CO2 may be a technically suitable working fluid for EOR in the BPS, supplies are limited and costs for using CO2 in EOR pilots are prohibitively high. Meanwhile, produced gas flaring has presented challenges for BPS operators in North Dakota. Analysis conducted by the North Dakota Pipeline Authority indicates that the current gas-gathering infrastructure in North Dakota is insufficient to accommodate all of the associated gas that is produced from the BPS. The geographically isolated location of North Dakota relative to large natural gas markets, combined with suppressed natural gas prices, has made it economically challenging for industry to invest capital in expanding gas-gathering infrastructure in the state. These circumstances led to a research program conducted by the Energy & Environmental Research Center (EERC) in partnership with Liberty Resources Management Company LLC (LR) to examine the potential to use rich gas injection for EOR and mitigate flaring. A rich gas EOR pilot test was designed and executed by LR at its Stomping Horse development area in Williams County, North Dakota. From July 2018 through May 2019, a total of 160 million standard cubic feet (MMscf) of rich produced gas was injected into the BPS using five different wells in a sequential injection strategy. LR’s Leon–Gohrick drill spacing unit (DSU) was used as the test site. Regulatory oversight was provided by the North Dakota Industrial Commission (NDIC). Technical support was provided by the EERC through a series of laboratory, modeling, and field-based activities, and additional post-pilot research activities incorporated learnings from the test, developed new laboratory data, improved fracture modeling methods, and developed machine learning and big data analytics. The results from the Stomping Horse rich gas EOR pilot activities indicate that developing an effective, economical EOR approach for the BPS will require more field tests. Another key lesson learned from the Stomping Horse tests is that detailed pre- and posttest data on reservoir conditions and fluids production are essential. Robust reservoir characterization provides information that is crucial to creating realistic geomodels and conducting valid dynamic simulations of potential EOR scenarios. A detailed understanding of the completions and production history of offset wells is also necessary for valid test result interpretations. This knowledge is essential to designing the operational parameters of injectivity tests and interpreting the results. A conformance control strategy is also essential to success. Laboratory-based examinations of rich gas interactions with reservoir fluids and rocks were conducted, with an emphasis on determining the ability to mobilize oil in the tight reservoir rocks and shales of the BPS. Injection fluid composition was shown to have a positive impact on reducing reservoir oil minimum miscibility pressure (MMP), reducing interfacial tension (IFT), and altering wettability. IFT and contact angle measurements demonstrated that wettability can be altered in the presence of rich gas, suggesting the potential to improve oil recovery. Iterative modeling of surface infrastructure and reservoir performance using data generated by the various project activities was conducted. A geologic model of the Stomping Horse area was built; history-matched oil, gas, and water production was used in simulations of various EOR scenarios. Early programmatic modeling results were used to support LR’s design and operation of the EOR pilot and to provide insight regarding optimization of future commercial-scale BPS EOR design and operations. Post-pilot modeling focused on alternative methods of understanding complex fracture networks and accelerating simulation time. These led to improved simulation run times and provide excellent history-matching results. Several of these iterative models were used as the bases for developing algorithms into machine learning and big data analytics. History matching in reservoir simulation is time-consuming and computer processing-intensive. Machine learning algorithms were created, and an automated history-matching tool was developed. A large set of synthetic reservoir simulations were created to generate well responses (oil, gas, and water production, well bottomhole pressure [BHP], and tracer or propane breakthrough) for a set of EOR operating parameters that included offset well status (open or closed), injectate (rich gas or propane), injection rate, and injection well BHP. A user interface was developed to provide real-time visualization. Machine learning-based models were developed to provide rapid forecasting of well performance given a set of user-defined EOR operating parameters. These predictive models allow the user to modify the offset well status, injection rate, and injection well BHP and rapidly forecast future production performance. The combination of real-time visualization tools with real-time forecasting tools provides a framework for real-time control—operational changes that the EOR site operator can enact (e.g., changing gas injection rates) to affect the observed performance and potentially improve the EOR outcome. There is great reason to be optimistic about the future of EOR in the Bakken. The results of the laboratory studies suggest significant potential for high rates of oil mobilization using produced field gas injection under the right conditions. The results of the lab studies, combined with rigorous statistical analysis of well production data and associated modeling efforts, confirm the notion that fluid mobility within the reservoir is controlled by fractures. As more knowledge is gained about the nature and distribution of fracture networks in the Bakken, the industry will be in a better position to predict and, ultimately, influence fluid mobility. New field tests are necessary to develop a more complete understanding of those conditions. Thoughtful and creatively engineered field tests within a well-characterized geologic setting will yield the fundamental knowledge needed to take Bakken oil production to the next level. This subtask was cofunded through the EERC–U.S. Department of Energy Joint Program on Research and Development for Fossil Energy-Related Resources Cooperative Agreement No. DE-FE0024233. Nonfederal funding was provided by the North Dakota Industrial Commission’s Oil and Gas Research Program and Computer Modelling Group.

04 OIL SHALES AND TAR SANDS↗

DOE OSTI2022

Westinghouse Accident Tolerant Fuel Phase 2B with Higher Enriched and Higher Burnup Add-On Project Final Technical Report Deliverable Volume 1

The Westinghouse Electric Company LLC (Westinghouse) accident tolerant fuel (ATF) program funded by Westinghouse and the US Department of Energy (DOE) utilized chromium (Cr) coated zirconium alloy cladding with doped UO 2 (ADOPT(TM)) and uranium silicide (U 3 Si 2 ) high density/high thermal conductivity fuel for its lead test rod (LTR) program with irradiation beginning in 2019. ADOPT is a Cr 2 O 3 +Al 2 O 3 doped UO 2 pellet with increased oxidation resistance and density, and increased resistance to fission gas release. Due to the issues with U 3 Si 2 fuel reaction in pressurized water reactor (PWR) environments, uranium nitride (U 15 N) has been substituted for U 3 Si 2 as a long-term fuel option. Cr coated cladding with ADOPT fuel is the near-term Westinghouse EnCore fuel product. The long-term EnCore ® fuel product is SiGA ® SiC/SiC composite cladding with high density/high thermal conductivity UN fuel. In 2020, the higher burnup and higher enriched and (HBHE) program was integrated into the ATF program. The objective of this expanded program is to extend ATF burnups to at least 75 MWd/kgU to economically facilitate 24-month cycles in PWRs. The ATF program now includes ADOPT fuel with enrichments >5% 235 U. Over the past several years, Westinghouse has tested the Cr coated zirconium (Zr) and silicon carbide (SiC) claddings in Westinghouse Churchill autoclaves and the Massachusetts Institute of Technology (MIT) reactor. Additionally, the Cr coated cladding has been tested in the Advanced Test Reactor (ATR) at Idaho National Laboratory (INL) and in LTR programs at the Doel 4 and Bryon 2 commercial reactors. High temperature tests at the state-of-the-art Westinghouse facilities in Churchill, PA and at Karlsruhe Institute of Technology (KIT) have been carried out to determine the time and temperature limits for the Cr coated zirconium claddings. These tests indicate that Cr coated Zr cladding can take temperatures up to about 1500°C for short periods of time without becoming totally oxidized. The main issue is the formation of the Cr-Zr eutectic at 1333°C resulting in the migration of this eutectic inward with the formation of ZrCr 2 which rapidly oxidizes to ZrO 2 on the outside of the tube. The Cr coated Zr cladding also delays the bursting of the tube, reduces the burst area, and reduces hydriding and embrittlement of the Zr which is a major benefit for reducing fuel fragmentation, redistribution, and dispersal (FFRD), the major licensing issue faced by HBHE. The manufacturing parameters for the SiC have been found to have a significant effect on the corrosion rate of the SiC in light water reactor (LWR) conditions and significant improvements have been made in the yield, quality, and oxidation resistance of the SiC cladding by identifying and tightening manufacturing process parameters. Current autoclave results for SiC composite claddings indicate that a corrosion rate of fewer than 2 micrometers per year can be achieved which meets corrosion requirements under normal operating conditions. High temperature steam oxidation tests at the state-of-the-art Westinghouse facilities in Churchill, PA and at Karlsruhe Institute of Technology (KIT) have been carried and determined that SiC cladding can withstand temperatures up to about 1800°C to 1900°C before excessive corrosion reactions begin, and that the SiC will not balloon and burst. ADOPT fuel pellets were incorporated into the Westinghouse ATF program because of their increased density and oxidation resistance. A topical report supporting the use of ADOPT pellets has been submitted to the NRC and will be approved shortly. ADOPT pellets have been a product in Europe for over 15 years. The additives work to make larger grain sizes which appear to reduce fission gas release during transients as well as increase the density of the pellet. Fuel rod and assembly design in preparation for the lead test rod (LTR) and lead test assembly (LTA) programs is underway as well as licensing efforts with the Nuclear Regulatory Commission (NRC). Finally, accident analyses coupled with economic evaluations for both operating savings as well as fuel savings have been initiated. Modular Accident Analysis Program 5 (MAAP5) calculations indicate that both Cr coated cladding and SiC composite cladding can increase the time to fuel rod loss of geometry for up to two hours longer than current Zr based cladding in a station blackout scenario. This additional time is due to the much higher oxidation resistance of the SiC and Cr coated cladding. These two hours can be used to implement additional response options instituted through the Diverse and Flexible Mitigation Capability (FLEX) program by reactor operators. Both ATF cladding options also reduce hydrogen production which dramatically reduces primary system and containment pressure and the risk of fission product release beyond containment in the unlikely event of an accident. The lower pressure in the system allows more time to feed cooling water to the core, resulting in the avoidance of fuel melting. The coping time is extended indefinitely if the modest water flow provided by FLEX continues. Experimental work on methods to rapidly fabricate SiC composite structures with high density and reduce the fabrication price of SiC fibers while maintaining a high level of performance is needed. Methods for mitigating or stopping the Cr-Zr eutectic are also needed to further improve the oxidation resistance of Cr coated Zr. Minimal (<1%) swelling of U 3 Si 2 and subsequent fission gas release has been demonstrated up to 20 MWd/kgU. Irradiation experiments with U 3 Si 2 fuel in ATR to determine these properties at 40 MWd/kgU were completed but post irradiation examinations were not done since U 3 Si 2 has been replaced by UN. UN has three issues that are being addressed. The first is increasing the oxidation resistance so that there is not excessive reaction up to ~1500°C. While UN increases the pellet density and additives to the UN postponed the temperature at which rapid oxidation occurred by ~100°C to 200°C, this is not enough, and pellet coating options are now being pursued. The second is developing a method to manufacture that does not require the multi-step process of UF 6 to UO 2 to UC to UN. Efforts are underway looking at UF 6 to UN 2 to UN and UF 6 to UF 4 to UN 2 to UN reactions. Finally, an economically acceptable process for potentially enriching the 15 N content of the nitrogen used to make UN to >95% 15N was identified though not experimentally validated. In addition to the work supported by the DOE, research and testing activities are being carried on in a world-wide effort funded by many countries such as Sweden, United Kingdom, Belgium, Netherlands, Spain, Germany, Japan, and France. This work is being facilitated through the Westinghouse led Collaboration for Advanced Research on Accident Tolerant Fuel (CARAT) program. Annual meetings were organized in 2018 and 2019 by Westinghouse as a venue for presentation of this work and to provide for the cross-fertilization of ideas among the many outstanding researchers in the ATF area. No meetings were held in 2020 and 2021 due to the COVID-19 related restrictions on traveling and gatherings. Since SiC, coated cladding, and high-density fuel options are not currently used in the nuclear industry; support from the government and industry members is needed to further the significant effort of setting new standards. The same is true for the Nuclear Regulatory Commission (NRC) which must review and approve the commercial use of these new fuels since all current regulations are oriented toward Zr/UO 2 fuel. Several meetings have been held with the NRC to generate a fast-track approach to licensing ATF using a combination of atomic modeling, in-rod sensors, and in-reactor testing. This approach was memorialized in an “Accelerated Fuel Qualification White Paper” prepared by the Accelerated Fuel Qualification Working Group lead by General Atomics.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2022

Westinghouse Accident Tolerant Fuel Phase 2B with Higher Enriched and Higher Burnup Add-On Project Final Technical Report Deliverable Volume 2

The Westinghouse Electric Company LLC (Westinghouse) accident tolerant fuel (ATF) program funded by Westinghouse and the US Department of Energy (DOE) utilized chromium (Cr) coated zirconium alloy cladding with doped UO 2 (ADOPT(TM)) and uranium silicide (U 3 Si 2 ) high density/high thermal conductivity fuel for its lead test rod (LTR) program with irradiation beginning in 2019. ADOPT is a Cr 2 O 3 +Al 2 O 3 doped UO 2 pellet with increased oxidation resistance and density, and increased resistance to fission gas release. Due to the issues with U 3 Si 2 fuel reaction in pressurized water reactor (PWR) environments, uranium nitride (U 15 N) has been substituted for U 3 Si 2 as a long-term fuel option. Cr coated cladding with ADOPT fuel is the near-term Westinghouse EnCore ® fuel product. The long-term EnCore fuel product is SiGA ® SiC/SiC composite cladding with high density/high thermal conductivity UN fuel. In 2020, the higher burnup and higher enriched and (HBHE) program was integrated into the ATF program. The objective of this expanded program is to extend ATF burnups to at least 75 MWd/kgU to economically facilitate 24-month cycles in PWRs. The ATF program now includes ADOPT fuel with enrichments >5% 235 U. Over the past several years, Westinghouse has tested the Cr coated zirconium (Zr) and silicon carbide (SiC) claddings in Westinghouse Churchill autoclaves and the Massachusetts Institute of Technology (MIT) reactor. Additionally, the Cr coated cladding has been tested in the Advanced Test Reactor (ATR) at Idaho National Laboratory (INL) and in LTR programs at the Doel 4 and Bryon 2 commercial reactors. High temperature tests at the state-of-the-art Westinghouse facilities in Churchill, PA and at Karlsruhe Institute of Technology (KIT) have been carried out to determine the time and temperature limits for the Cr coated zirconium claddings. These tests indicate that Cr coated Zr cladding can take temperatures up to about 1500°C for short periods of time without becoming totally oxidized. The main issue is the formation of the Cr-Zr eutectic at 1333°C resulting in the migration of this eutectic inward with the formation of ZrCr 2 which rapidly oxidizes to ZrO 2 on the outside of the tube. The Cr coated Zr cladding also delays the bursting of the tube, reduces the burst area, and reduces hydriding and embrittlement of the Zr which is a major benefit for reducing fuel fragmentation, redistribution, and dispersal (FFRD), the major licensing issue faced by HBHE. The manufacturing parameters for the SiC have been found to have a significant effect on the corrosion rate of the SiC in light water reactor (LWR) conditions and significant improvements have been made in the yield, quality, and oxidation resistance of the SiC cladding by identifying and tightening manufacturing process parameters. Current autoclave results for SiC composite claddings indicate that a corrosion rate of fewer than 2 micrometers per year can be achieved which meets corrosion requirements under normal operating conditions. High temperature steam oxidation tests at the state-of-the-art Westinghouse facilities in Churchill, PA and at Karlsruhe Institute of Technology (KIT) have been carried and determined that SiC cladding can withstand temperatures up to about 1800°C to 1900°C before excessive corrosion reactions begin, and that the SiC will not balloon and burst. ADOPT fuel pellets were incorporated into the Westinghouse ATF program because of their increased density and oxidation resistance. A topical report supporting the use of ADOPT pellets has been submitted to the NRC and will be approved shortly. ADOPT pellets have been a product in Europe for over 15 years. The additives work to make larger grain sizes which appear to reduce fission gas release during transients as well as increase the density of the pellet. Fuel rod and assembly design in preparation for the lead test rod (LTR) and lead test assembly (LTA) programs is underway as well as licensing efforts with the Nuclear Regulatory Commission (NRC). Finally, accident analyses coupled with economic evaluations for both operating savings as well as fuel savings have been initiated. Modular Accident Analysis Program 5 (MAAP5) calculations indicate that both Cr coated cladding and SiC composite cladding can increase the time to fuel rod loss of geometry for up to two hours longer than current Zr based cladding in a station blackout scenario. This additional time is due to the much higher oxidation resistance of the SiC and Cr coated cladding. These two hours can be used to implement additional response options instituted through the Diverse and Flexible Mitigation Capability (FLEX) program by reactor operators. Both ATF cladding options also reduce hydrogen production which dramatically reduces primary system and containment pressure and the risk of fission product release beyond containment in the unlikely event of an accident. The lower pressure in the system allows more time to feed cooling water to the core, resulting in the avoidance of fuel melting. The coping time is extended indefinitely if the modest water flow provided by FLEX continues. Experimental work on methods to rapidly fabricate SiC composite structures with high density and reduce the fabrication price of SiC fibers while maintaining a high level of performance is needed. Methods for mitigating or stopping the Cr-Zr eutectic are also needed to further improve the oxidation resistance of Cr coated Zr. Minimal (<1%) swelling of U 3 Si 2 and subsequent fission gas release has been demonstrated up to 20 MWd/kgU. Irradiation experiments with U 3 Si 2 fuel in ATR to determine these properties at 40 MWd/kgU were completed but post irradiation examinations were not done since U 3 Si 2 has been replaced by UN. UN has three issues that are being addressed. The first is increasing the oxidation resistance so that there is not excessive reaction up to ~1500°C. While UN increases the pellet density and additives to the UN postponed the temperature at which rapid oxidation occurred by ~100°C to 200°C, this is not enough, and pellet coating options are now being pursued. The second is developing a method to manufacture that does not require the multi-step process of UF 6 to UO 2 to UC to UN. Efforts are underway looking at UF 6 to UN 2 to UN and UF 6 to UF 4 to UN 2 to UN reactions. Finally, an economically acceptable process for potentially enriching the 15 N content of the nitrogen used to make UN to >95% 15 N was identified though not experimentally validated. In addition to the work supported by the DOE, research and testing activities are being carried on in a world-wide effort funded by many countries such as Sweden, United Kingdom, Belgium, Netherlands, Spain, Germany, Japan, and France. This work is being facilitated through the Westinghouse led Collaboration for Advanced Research on Accident Tolerant Fuel (CARAT) program. Annual meetings were organized in 2018 and 2019 by Westinghouse as a venue for presentation of this work and to provide for the cross-fertilization of ideas among the many outstanding researchers in the ATF area. No meetings were held in 2020 and 2021 due to the COVID-19 related restrictions on traveling and gatherings. Since SiC, coated cladding, and high-density fuel options are not currently used in the nuclear industry; support from the government and industry members is needed to further the significant effort of setting new standards. The same is true for the Nuclear Regulatory Commission (NRC) which must review and approve the commercial use of these new fuels since all current regulations are oriented toward Zr/UO 2 fuel. Several meetings have been held with the NRC to generate a fast-track approach to licensing ATF using a combination of atomic modeling, in-rod sensors, and in-reactor testing. This approach was memorialized in an “Accelerated Fuel Qualification White Paper” prepared by the Accelerated Fuel Qualification Working Group lead by General Atomics.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2023

Optimizing Information Automation Using a New Method Based on System-Theoretic Process Analysis: Tool Development and Method Evaluation

This report is an update to a prior report that describes progress and findings for a program of research supporting the design and optimization of information automation systems for nuclear power plants. Much of the domestic nuclear fleet is currently focused on modernizing technologies and processes, including transitioning toward digitalization in the control room and throughout the plant, along with a greater interest in the use of automation, artificial intelligence, robotics, and other emerging technologies. While there are significant opportunities to apply these technologies toward greater plant safety, efficiency, and overall cost-effectiveness, optimizing their design and avoiding potential safety and performance risks depends on ensuring that human performance-related organizational and technical design issues are identified and addressed early in the design process. This report describes modeling tools and techniques, based on sociotechnical systems theory, to support these design goals and their application in the current research effort. The report is primarily intended for senior nuclear energy stakeholders, including regulators, corporate management, and senior plant management. We have developed and employed a method to design an optimized information automation ecosystem (IAE) based on the systems-theoretic constructs underlying sociotechnical systems theory in general and the Systems-Theoretic Accident Modeling and Processes (STAMP) approach in particular. We argue that an IAE can be modeled as an interactive information control system whose behavior can be understood in terms of dynamic control, feedback, and communication relationships amongst the system’s technical and organizational components. We have employed two STAMP-based tools in this effort. The first is Causal Analysis based on STAMP (CAST), an accident and incident analysis technique that was used to examine a performance- and safety-related incident at an industry partner’s plant involving the unintentional activation of an emergency diesel generator. This analysis provided insight into the behavior of the plant’s current information control structure within the context of a specific, significant event. The second tool is Systems Theoretic Process Analysis (STPA) which is a proactive risk analysis tool used to examine existing and potential, planned sociotechnical systems. STPA was used to identify risk factors in the current design of a generic nuclear power plant (NPP) preventive maintenance system. Our analyses focused on identifying near-term system improvements and longer-term design requirements for an optimized IAE system. CAST analysis findings indicate an important underlying contributor to the incident under investigation, and a significant risk to information automation system performance, was perceived time and schedule pressure, which exposed weaknesses in interdepartmental coordination between and within responsible plant organizations and challenged the resilience of established plant processes, until a human caused the eventual event. These findings are discussed in terms of their risk to overall system performance and their implications for information automation system resilience and brittleness. STPA findings exposed several areas of concern in the design of current preventive maintenance systems. We also present two preliminary information automation models. The proactive issue resolution (PIR) model is a test case of an information automation concept with significant near-term potential for application and subsequent reduction in significant plant events. The IAE model is a more general representation of a broader, plantwide information automation system and represents an end-state vision for our work. From our results, we have generated an initial set of preliminary system-level requirements and safety constraints for these models. We have also focused on early development of easy to learn, easy to use “transportable” tools for sociotechnical systems analysis. We intend these to be used by NPP personnel as a means of gaining reliable and relatively quick insight into (1) sociotechnical systems factors impacting incidents and accidents, (2) potential sociotechnical risk factors in existing or planned system designs, and (3) potential weaknesses in a system’s safety and/or information control structure. We conclude the report with a set of summary recommendations, a discussion of planned and potential follow-on research and development, and a draft list of system-level requirements and safety constraints for optimized information automation systems.

99 GENERAL AND MISCELLANEOUS↗

DOE OSTI2020

A Quantum-Based Approach to Predict Primary Radiation Damage in Polymeric Networks

Initial atomistic-level radiation damage in chemically reactive materials is thought to induce reaction cascades that can result in undesirable degradation of macroscale properties. Ensembles of quantum-based molecular dynamics (QMD) simulations can accurately predict these cascades, but extracting chemical insights from the many underlying trajectories is a labor-intensive process that can require substantial a priori intuition. We develop here a general and automated graph-based approach to extract all chemically distinct structures sampled in QMD simulations and apply our approach to predict primary radiation damage of polydimethylsiloxane (PDMS), the main constituent of silicones. A postprocessing protocol is developed to identify underlying polymer backbone structures as connected components in QMD trajectories. Furthermore, these backbones form a repository of radiation-damaged structures. A scheme for extracting and updating a library of isomorphically distinct structures is proposed to identify the spanning set and aid chemical interpretation of the repository. The analyses are applied to ensembles of cascade QMD simulations in which the four element types in PDMS are selectively excited in primary knock-on atom events. Our approach reveals a much higher degree of combinatorial complexity in this system than was inferred through radiolysis experiments. Probabilities are extracted for radiation-induced network changes including formation of branch points, carbon linkages, cycles, bond scissions, and carbon uptake into the Si–O siloxane backbone network. The general analysis framework presented here is readily extendable to modeling chemical degradation of other polymers and molecular materials and provides a basis for future quantum-informed multiscale modeling of radiation damage.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Kimberlina 1.2 CCUS Geophysical Models and Synthetic Data Sets

This synthetic multi-scale and multi-physics data set was produced in collaboration with teams at the Lawrence Berkeley National Laboratory, National Energy Technology Laboratory, Los Alamos National Laboratory, and Colorado School of Mines through the Science-informed Machine Learning for Accelerating Real-Time Decisions in Subsurface Applications (SMART) Initiative. Data are associated with the following publication: Alumbaugh, D., Gasperikova, E., Crandall, D., Commer, M., Feng, S., Harbert, W., Li, Y., Lin, Y., and Samarasinghe, S., “The Kimberlina Synthetic Geophysical Model and Data Set for CO2 Monitoring Investigations”, The Geoscience Data Journal, 2023, DOI: 10.1002/gdj3.191. The dataset uses the Kimberlina 1.2 CO2 reservoir flow model simulations based on a hypothetical CO2 storage site in California (Birkholzer et al., 2011; Wainwright et al., 2013). Geophysical properties models (P- and S-wave seismic velocities, saturated density, and electrical resistivity) were produced with an approach similar to that of Yang et al. (2019) and Gasperikova et al. (2022) for 100 Kimberlina 1.2 reservoir models. Links to individual resources are provided below: [CO2 Saturation Models](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-co2-saturation-models); Resistivity Models – [part 1](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-resistivity-models-part-1), [part 2](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-resistivity-models-part-2), and [part 3](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-resistivity-models-part-3); [Vp Velocity Models](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-vp-velocity-models); [Vs Velocity Models](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-vs-velocity-models); [Density Models](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-density-models). The 3D distributions of geophysical properties for the 33 time stamps of the SIM001 model were used to generate synthetic seismic, gravity, and electromagnetic (EM) responses for 33 times between zero and 200 years. Synthetic surface seismic data were generated using 2D and 3D finite-difference codes that simulate the acoustic wave equation (Moczo et al., 2007). 2D data were simulated for six point-pressure sources along a 2D line with 10 m receiver spacing and a time spacing of 0.0005 s. 3D simulations were completed for 25 surface pressure sources using a source separation of 1 km in both the x and y directions and a time spacing of 0.001 s. Links to individual resources are provided below: [2D velocity models](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-2d-velocity-models) and [2D surface seismic data](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-2d-surface-seismic-data). [3D velocity models](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-velocity-models), and 3D seismic data [year0](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year0), [year1](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year1), [year2](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year2), [year5](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year5), [year10](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year10), [year15](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year15), [year20](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year20), [year25](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year25), [year30](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year30), [year35](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year35), [year40](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year40), [year45](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year45), [year49](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year49), [year50](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year50), [year51](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year51), [year52](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year52), [year55](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year55), [year60](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year60), [year65](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year65), [year70](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year70), [year75](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year75), [year80](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year80), [year85](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year85), [year90](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year90), [year95](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year95), [year100](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year100), [year110](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year110), [year120](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year120), [year130](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year130), [year140](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year140), [year150](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year150), [year175](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year175), [year200](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year200). The Python scripts to read these models and data are provided [here](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-python-scripts). EM simulations used a borehole-to-surface survey configuration, with the source located near the reservoir level and receivers on the surface using the code developed by Commer and Newman (2008). Pseudo-2D data for the source at [2500 m](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-pseudo-2d-csem-data-tz2500m) and [3025 m](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-pseudo-2d-csem-data-tz3025m), used a 2D inline receiver configuration to simulate a response over 3D resistivity models. The [3D data](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-csem-data) contain electric fields generated by borehole sources at monitoring well locations and measured over a surface receiver grid. Vector gravity data, both on the surface and in boreholes, were simulated using a modeling code developed by Rim and Li (2015). The simulation scenarios were parallel to those used for the EM: [pseudo-2D data](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-gravity-data) were calculated along the same lines and within the same boreholes, and [3D data](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-gravity-data) were simulated over 3D models on the surface and in three monitoring wells. A series of [synthetic well logs](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-well-logs) of CO2 saturation, acoustic velocity, density, and induction resistivity in the injection well and three monitoring wells are also provided at 0, 1, 2, 5, 10, 15, and 20 years after the initiation of injection. These were constructed by combining the low-frequency trend of the geophysical models with the high-frequency variations of actual well logs collected in the Kimberlina 1 well that was drilled at the proposed site. Measurements of permeability and pore connectivity were made on cores of Vedder Sandstone, which forms the primary reservoir unit: [CT micro scans](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-ct-micro-scans-of-vedder-formation) and [Industrial CT Images](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-industrial-ct-images-vedder-formation). These measurements provide the range of scales in the otherwise synthetic data set to be as close to a real-world situation as possible. References: Birkholzer, J.T., Zhou, Q., Cortis, A. and Finsterle, S., 2011. A sensitivity study on regional pressure buildup from large-scale CO2 storage projects. Energy Procedia, 4, 4371-4378. Commer, M., and Newman, G.A., 2008. New advances in three-dimensional controlled-source electromagnetic inversion, Geophysical Journal International, 172, 513-535. Gasperikova, E., Appriou, D., Bonneville, A., Feng, Z., Huang, L., Gao, K., Yang, X., Daley, T., 2022, Sensitivity of geophysical techniques for monitoring secondary CO2 storage plumes, Int. J. Greenh. Gas Control, Volume 114, 103585, ISSN 1750-5836, https://doi.org/10.1016/j.ijggc.2022.103585. Moczo, P., J.O. Robertsson and L. Eisner, 2007, The finite-difference time-domain method for modeling of seismic wave propagation: Advances in geophysics, 48, 421-516. Rim, H., and Y. Li, 2015, Advantages of borehole vector gravity in density imaging, Geophysics, 80, G1-G13. Wainwright, H. M.; Finsterle, S.; Zhou, Q.; Birkholzer, J. T., 2013. Modeling the Performance of Large-Scale CO2 Storage Systems: A Comparison of Different Sensitivity Analysis Methods. International Journal of Greenhouse Gas Control, 17, 189205. https://doi.org/10.1016/j.ijggc.2013.05.007, DOI: 10.18141/1603331. Yang, X., Buscheck, T.A., Mansoor, K., Wang, Z., Gao, K., Huang, L., Appriou, D., and Carroll, S.A., 2019. Assessment of geophysical monitoring methods for detection of brine and CO2 leakage in drinking water aquifers, International Journal of Greenhouse Gas Control, 90, 102803, https://doi.org/10.1016/j.ijggc.2019.102803.

CCUS↗

DOE OSTI2020

Austin Sustainable and Holistic Integration of Energy Storage and Solar PV [Austin SHINES]. Final Report, Version 2

The Austin SHINES project and solution is a software management platform, for an electric grid with a high penetration of dispersed photovoltaic (PV) solar generation sites, which maintains the traditional power quality and reliability associated with grid service. This project developed and deployed the platform as a Distributed Energy Resource Management System (DERMS), engaging multiple advanced controls, to evaluate operation and optimization of a fleet of diverse DER assets, installed at several locations among Austin Energy’s customers and distribution system. The project also produced a methodology to create a replicable DERMS template, adaptable to other regions and market structures. Last, Austin SHINES aimed to demonstrate the solution’s methodology would enable the DER grid ecosystem to serve load at a technical cost (System Levelized Cost of Electricity, or System LCOE) of less than the U.S. Department of Energy SHINES program metric of $0.14/kWh, in a defined boundary, while enabling a high penetration of distributed PV. Research was categorized in 6 reports (Final Deliverables = FD) listed below, with titles and descriptions indicating which area of understanding was investigated: FD-1: System Levelized Cost of Electricity (System LCOE) Methodology The creation and use of the System LCOE to Serve Load metric that encompasses the holistic, system-level costs and benefits of all resources, and enables them to be evaluated based on their ability to support an efficient and low-cost integrated grid ecosystem. FD-2: Software Platform Product Description The creation of new DER control methodologies deployable within a utility-grade software platform that enable DER's to maximize their benefit within a grid, that is capable of serving load enabling a high penetration of distributed PV generation. FD-3: Optimal Design Methodology Optimal design methodologies for individual DER installations that enable utilities to determine the optimal combinations and sizing for individual DER sites. FD-4: Austin SHINES Ownership and Operation Models for DER System Performance A comparison of multiple DER aggregation and ownership methodologies including direct utility control, third-party aggregator, and autonomous. FD-5: Economic Modeling & Optimization A comparison of multiple DER technology mixes and configurations within the distribution system, providing insight into an optimal blend of technologies that best enable the distribution system to serve load at the lowest cost at high penetrations of solar. FD-6: Fielded Assets Deployed DER assets within the Austin Energy SHINES circuits. Austin SHINES provided an opening for state-of-the-art technology products to be deployed, providing a rich opportunity for improving how each of the products perform as stand-alone products, and in concert with other complementary products. The Austin SHINES project comprised of two key metrics for System LCOE: SystemLCOE_SHINES<$0.14/kWh Modeled ΔSystemLCOE_SHINES/ΔSystemLCOE_Base≥20% at same solar penetration The System LCOE calculation uses the costs of the utility-owned infrastructure as it exists today, the cost of the DERs that exist in the system today, and the cost of the purchase of energy from ERCOT wholesale markets over the course of the calendar year. All costs are on an annualized basis. The capital and operating costs are derived from the rate case, which produces a yearly cost. The net cost of energy and services imported to the system is integrated over the test year, as is the load served and solar penetration. The first metric was easily achieved by every scenario considered. The goal was set when the Department of Energy’s SHINES Funding Opportunity Announcement was written in 2015 and was a more difficult target at the time. Due mostly to rapidly declining costs for DERs and the significant decrease in the Electric Reliability Council of Texas (ERCOT) energy market prices, which results in lower net cost of energy purchases, the System LCOE is well below this target for all scenarios considered. A fleet of DERs can assume different mixtures, each of which serves the load at a different LCOE. The optimal mixture of DERs serves load at the smallest System LCOE. The second metric (hereinafter %delta metric) asks that the holistic DERMS controls reduce the incremental cost above the baseline of going to a high solar penetration future by at least 20% as compared to the case of a DER deployment with no sophisticated controls (autonomous). Many comparison sets were created throughout this project. Physical technology was installed for informing utility engineering and testing several types of operational control schemes, through the DERMS. The types of operational control which were compared for valuation of the System LCOE Metric were: Holistic control = using the full suite of the DERMS platform to decide and optimize how/why the systems operate depending on weather, market, and reliability signal input. Autonomous control = a local mode at the asset site, wherein a schedule operates the asset, with visibility into performance only No control = the baseline for comparing value against the other two types of control The types of ownership control included: Direct Utility control = the utility dispatches a signal to each asset Third-Party Aggregator = a third party aggregates a fleet of assets and the utility dispatches one signal for all Autonomous = a local mode is set for operation at the asset site, wherein a schedule operates the asset, with visibility into performance only The types of control methodologies deployable within a utility-grade software platform included: Utility Peak Load Reduction = Lower transmission cost obligation Day-Ahead Energy Arbitrage = Realize economic value through price differential Real-Time Price Dispatch = Realize economic value from real-time price spikes Voltage support = Reduce losses and increase solar generation Distribution Congestion Management = Increase local grid reliability Demand Charge Reduction = Lower customer bills and realize system benefit The fielded assets deployed for the project were: Utility Scale Kingsbery Energy Storage System: 1.5 MW / 3 MWh Li-Ion battery storage Mueller Energy Storage System: 1.75 MW / 3.2 MWh Li-Ion battery storage, 7 Energy Storage Units (250 kW each) La Loma Community Solar: 2.6 MW Commercial Scale Aggregated storage installations at 3 sites, with existing solar (300+ kW): One 18 kW / 36 kWh Li-Ion battery storage Two 72 kW / 144 kWh Li-Ion battery storage Residential Scale Aggregated storage installations: -Six stationary battery storage systems (10 kWh each) at homes with existing solar -One Electric Vehicle installed as Vehicle-to-Grid (V2G) Utility-Controlled Solar via Smart Inverters at 12 homes Autonomously-Controlled Smart Inverters at 6 homes Over the course of the project, Austin SHINES undertook installing more than 3 MW of distributed battery energy storage, smart PV inverters, a DER control platform, and other enabling technologies utilizing customer and utility locations and aggregation models. All of these resources were to be integrated and optimized at the utility level. DER assets and control methodologies were designed to achieve a credible pathway to a System LCOE for energy delivered to load of $0.14//kWh or less by 2020, while maximizing distributed solar generation and maintaining acceptable standards of power quality. The project also established a template for other regions to follow, to maximize the adoption of distributed solar PV in support of an economic and efficient grid. In total, the Austin SHINES project added value to the DER subject area in each layer of integration. From utility, to commercial to residential scales, the sheer hierarchy of communication and coordination was a significant accomplishment in addition to learnings from what these communications revealed was unique to each. Economically, the most effective method demonstrated was the criticality of planning phases. Contingencies and multiple projection scenarios helped guide the project to deploy optimal design as close as feasible, in real world conditions. The project and reports will serve public benefit by outlining specific areas of DER strategy and installation where many stakeholders and needs can be addressed with improved efficiency. Overall, communities and utilities should use the results to guide the increasing options available for powering the grid with DER, renewables, and carbon considerate energy.

14 SOLAR ENERGY↗

DOE OSTI2020

What can the GRACE seasonal cycle tell us about lake-aquifer interactions?

Lake-aquifer interactions have been the subject of investigation and debate for decades. Traditional investigation methods include measurement of water flux across the groundwater-surface water interface, application of heat and environmental tracer methods, conducting numerical simulations of the water flow, and mass balance-based approaches. We first review, evaluate and describe the different approaches that have been applied for examining lake-aquifer interactions and propose an additional complementary approach. While the leakage from lakes and its slow movement through porous media (laminar flow) is well characterized by existing methods, modeling the rapid and turbulent flow through preferred pathways (faults/karst) remains a challenging task. To better understand the nature, and the full scale, of the lake aquifer interactions, and using Lake Nasser and the underlying Dakhla subbasin of the Nubian Sandstone Aquifer System (NSAS) in the Western Desert of Egypt (area: 0.66 x 10 6 km 2 ) as a test site, we examine, model, and correlate temporally and spatially, the variations in Gravity Recovery and Climate Experiment terrestrial water storage (GRACE TWS ) solutions, with precipitation, lake levels, area, and water volume. We review current GRACE applications in hydrology and present our novel approach. Findings include: firstly, large seasonal fluctuations (peak: Nov./Dec.; trough: July/Aug.) in Lake Nasser's surface water levels (average annual fluctuations: from 3.7 m to 7.5 m), area (average area: 3622 km 2 but up to 4530 km 2 ), and volume (annual average: 13.4 km 3 , and up to 34 km 3 ) are observed throughout years 2006 to 2015. These fluctuations are accompanied by an increase in GRACE TWS (average: 50 ± 13 mm/yr, up to 77 ± 18 mm/yr) over Lake Nasser and by a progression of a front of increasing GRACE TWS values (> 50 ± 13 mm) away from the lake reaching distances of up to 700 km some 3 to 5 months following peak lake level periods. The areas witnessing the seasonal increase in GRACE TWS display a progressive increase in phase and decrease in amplitude with distance from the lake. Secondly, the negligible precipitation over the Dakhla subbasin cannot account for the observed seasonal GRACE TWS patterns and neither can the leakage signal from Lake Nasser. Thirdly, overall similarities in the distribution of GRACE TWS seasonal spatial patterns are observed. During high lake level periods (e.g., Dec. 2007; Feb. 2015) additional, and more distant, areas from Lake Nasser saw increases in GRACE TWS and vice versa during low lake level periods (e.g., July 2006; Aug. 2010). These observations are consistent with Lake Nasser being the main source of modern recharge for the Dakhla subbasin and suggest a new conceptual model for the subbasin: a slow groundwater flow through a porous matrix and a rapid groundwater flow along a network of faults, fractures, and karst topography across the matrix. We suggest that the proposed conceptual model for the interaction between Lake Nasser and the Dakhla subbasin could be applicable to aquifers of similar geologic, climatic, and hydrologic settings worldwide and that approaches similar to those advocated here could be used to investigate the validity of this suggestion.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Th 2 O – , Th 2 Au – , and Th 2 AuO 1,2 – Anions: Photoelectron Spectroscopic and Computational Characterization of Energetics and Bonding

The observation and characterization of the anions: Th 2 O – , Th 2 Au – , and Th 2 AuO 1,2 – is reported. These species were studied through a synergetic combination of anion photoelectron spectroscopy and ab initio correlated molecular orbital theory calculations at the CCSD(T) level with large correlation-consistent basis sets. To better understand the energetics and bonding in these anions and their corresponding neutrals, a range of smaller diatomic to tetratomic species were studied computationally. Correlated molecular orbital theory calculations at the CCSD(T) level showed that in most of these cases, there are close-lying anions and neutral clusters with different geometries and spin states and are consistent with the experimentally observed spectra. Thus, comparison of experimentally determined and computationally predicted vertical detachment energies and electron affinities for different optimized geometries and spin states shows excellent agreement to within 0.1 eV. The structures for both the neutrals and anions have a significant ionic component to the bonding because of the large electron affinity of the Au atom and modest ionization potentials for Th 2 , Th 2 O, and Th 2 O 2 . The analysis of the bonding for the Th–Th bonds from the molecular orbitals is consistent with this ionic model. The results show that there is a wide variation in the bond distance from 2.7 to 3.5 Å for the Th–Th bonds all of which are less than twice the atomic radius of Th of 3.6 Å. The bond distances encompass bond orders from 4 to 0. Furthermore, there can be different bond orders for the same bond distance depending on the nature of the ionic bonding suggesting that one may not be able to correlate the bond order with the bond distance in these types of clusters. In addition, the presence of an Au atom may provide a unique probe of the bonding in such clusters because of its ability to accept an electron from clusters with modest ionization potentials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Using Separation-Enhanced Isotope Ratio Mass Spectrometry to Enable Increased Renewable Carbon Content in Transportation Fuels (CRADA 525)

Stable isotope ratio measurements of carbon atoms using isotope ratio mass spectrometry (IRMS) can be an effective tool for quantifying biogenic carbon in co-processed fuels, with results approaching the precision and accuracy of accelerator mass spectrometry (AMS). The lower cost of an IRMS may enable deployment to refineries, improving access and analysis turnaround times (≤2 hours), and, by extension, provide data that can allow process optimization to maximize renewable carbon in desired refinery products. This project explored the integration of chemical separation with IRMS analyses to enable highly detailed tracking of biogenic carbon into fuel product streams separated by boiling point range, chemical class, or specific compound. Forty-nine fuels and fuel components of fossil and biogenic origin, spanning gasoline and diesel boiling point ranges, were received from three refiners and were analyzed for their δ 13 C values via IRMS. Results spanned a 13 C range from ca. 10‰ to 44‰ and reflect materials derived from sustainable sources (e.g., C4 or C3 plants, animal-based pathways, syngas) or from fossil-derived fuels. Common ranges are approximately 18‰ to 9‰ and approximately 30‰ to 20‰ for C4 and C3 plants, respectively, and approximately 34‰ to 24‰ and approximately 70‰ to 33‰ for petroleum-derived fuels and methane, respectively. Fuel-like standards were developed and tested using direct-injection elemental analyzer (EA) IRMS for liquid fuels. This method was compared with the published methods, yielding statistically similar results. Four blend curve sets were produced ranging from 0% to 100% of a fuel containing biogenic carbon, focusing on 0% to 10% biogenic carbon. Linear fits were the most applicable for two of the four blend curve sets; however, two sets were found to exhibit slightly quadratic behavior, which was more pronounced in low biogenic blend samples, necessitating second-order fits. The origin of the slight quadratic behavior remains unclear; however, the discussion points to possible interpretations. CanmetENERGY thoroughly characterized a majority of the samples using one- and two-dimensional gas chromatography (GC and GC×GC, respectively) and other analyses. Selected samples were subjected to solid phase extraction (SPE) for saturate, olefin, aromatic, and polar (SOAP) analysis, and the resulting solvent-diluted fractions containing saturates and aromatics were returned to Pacific Northwest National Laboratory (PNNL), where the solvent was removed via evaporation or physical separation using GC techniques. Characterization and separations provided an understanding of saturate and aromatic content, as well as boiling point ranges for each sample and sample fraction. Samples resulting from SPE were examined using EA-IRMS and gas chromatography combustion IRMS (GC-C-IRMS) analyses. Both approaches suggest that the range in values between end-members can be increased by selecting the paraffinic or aromatic fraction of the end-member or by selecting among individual compounds resulting from GC separation of the paraffinic fractions. Considerable work remains to put these approaches into practice and statistically validate the benefit for using a fraction or individual compound over bulk analysis of a sample. However, initial results suggest that separations provide advantages for samples having blend ratios of less than 10% biogenic blendstocks. 13 C results showed statistically similar biofuel blend results to those obtained at PNNL, although additional work is needed to obtain better reproducibility. Select samples were sent to Los Alamos National Laboratory (LANL) for IRMS measurements and Beta Analytics for AMS measurements. This work suggests that IRMS and AMS yield closely comparable results and in some circumstances, IRMS could serve as a surrogate for AMS. While additional work is needed to better resolve statistical advantages for separations and better show the comparable nature of IRMS and AMS in both the biogenic carbon analysis of bulk chemical classes, initial results from this study suggest that these should be pursued in order to proliferate this approach for quantifying biogenic carbon in transportation fuels to the refinery level, thereby potentially enabling process optimization in co-processing scenarios.

09 BIOMASS FUELS↗

DOE OSTI2022

Nafion Passivation of c-Si Surface and Edge for Electron Paramagnetic Resonance

Effective surface passivation of crystalline silicon (c-Si) surface by reducing the carrier recombination rate has led to modern c-Si solar cells with efficiencies > 25% in both laboratory and industrial settings. Typical mainstream surface passivation techniques include high-temperature silicon oxide (SiOx), amorphous silicon (a-Si:H), hydrogen-rich silicon nitride (SiNx), and aluminum oxide (Al2O3) [1]. They have demonstrated excellent surface recombination velocity of < 1 cm/s owing to both chemical passivation (via the hydrogen saturation of Si dangling bonds at the c-Si surface), and field-effect passivation (via the band bending from the fixed charge of the dielectric layers). Recently, several groups have studied solution-based organic materials for c-Si passivation, including bis(trifluoromethane)sulfonimide (TFSI), polystyrenesulfonate, and Nafion [2, 3] via spin or dip coating. All films were processed at ambient room temperature, and a high lifetime of 12 ms, as well as a low saturation current density (J0) of 16 fA/cm2, have been demonstrated using Nafion passivation [2]. However, despite the air instability, Nafion has several advantages when introduced to PV applications: (a) Wafer quality can be measured at high throughput after several process stages. (b) It is compatible with PL mapping, compared to HF liquid passivation which cannot be performed in room ambient. (c) Nafion process is fast and poses fewer constraints on process complexity, and cleanness, compared to Al2O3 and a-Si:H passivation. (d) Nafion is the ideal room temperature passivation, which can be applied onto a small fragment of a degraded module to investigate microscopic mechanisms when studying degradation mechanisms. Edge passivation is needed for advanced characterization tests such as Electron Paramagnetic Resonance (EPR), Electrically Detected Magnetic Resonance (EDMR), Deep Level Transient Spectroscopy (DLTS), local cell current-voltage (J-V), and Suns-Voc. In such cases, the laser-damaged minicell edges will obscure the true degradation mechanisms, and Nafion can effectively passivate them without causing further degradation through elevated-temperature processes. In this contribution, we show the passivation results of Nafion on bare nCz and compare it with a Al2O3 witness. We highlight the importance of edge pre-treatment before Nafion to achieve good passivation. We show that Nafion can reach decent passivation without undergoing high-temperature processes. We also explore the temperature dependency of Nafion through photoluminescence (PL) study and demonstrate the application of Nafion under cryogenic temperature (~6 K) through EPR. Our results reveal that Nafion can reduce surface and edge Si dangling bonds at low temperatures. Thus, it can be used as an effective room temperature passivation technique for advanced characterization methods.

c-Si↗