SEARCH · Engineering Papers

Results for “Conditional Normalizing Flow”

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 451 records · Page 25

All sources NASA NTRS DOE OSTI

DOE OSTI2015

Gas breakthrough and emission through unsaturated compacted clay in landfill final cover

Highlights: • Explore feasibility of unsaturated clay as a gas barrier in landfill cover. • Gas breakthrough pressure increases with clay thickness and degree of saturation. • Gas emission rate decreases with clay thickness and degree of saturation. • A 0.6 m-thick clay layer may be sufficient to meet gas emission rate limit. - Abstract: Determination of gas transport parameters in compacted clay plays a vital role for evaluating the effectiveness of soil barriers. The gas breakthrough pressure has been widely studied for saturated swelling clay buffer commonly used in high-level radioactive waste disposal facility where the generated gas pressure is very high (in the order of MPa). However, compacted clay in landfill cover is usually unsaturated and the generated landfill gas pressure is normally low (typically less than 10 kPa). Furthermore, effects of clay thickness and degree of saturation on gas breakthrough and emission rate in the context of unsaturated landfill cover has not been quantitatively investigated in previous studies. The feasibility of using unsaturated compacted clay as gas barrier in landfill covers is thus worthwhile to be explored over a wide range of landfill gas pressures under various degrees of saturation and clay thicknesses. In this study, to evaluate the effectiveness of unsaturated compacted clay to minimize gas emission, one-dimensional soil column tests were carried out on unsaturated compacted clay to determine gas breakthrough pressures at ultimate limit state (high pressure range) and gas emission rates at serviceability limit state (low pressure range). Various degrees of saturation and thicknesses of unsaturated clay sample were considered. Moreover, numerical simulations were carried out using a coupled gas–water flow finite element program (CODE-BRIGHT) to better understand the experimental results by extending the clay thickness and varying the degree of saturation to a broader range that is typical at different climate conditions. The results of experimental study and numerical simulation reveal that as the degree of saturation and thickness of clay increase, the gas breakthrough pressure increases but the gas emission rate decreases significantly. Under a gas pressure of 10 kPa (the upper bound limit of typical landfill gas pressure), a 0.6 m or thicker compacted clay is able to prevent gas breakthrough at degree of saturation of 60% or above (in humid regions). Furthermore, to meet the limit of gas emission rate set by the Australian guideline, a 0.6 m-thick clay layer may be sufficient even at low degree of saturation (i.e., 10% like in arid regions)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2021

AGR 5/6/7 Data Qualification Report for ATR Cycles 162B through 168A

This report provides the qualification status of experimental data for the Advanced Gas Reactor (AGR) 5/6/7 fuel irradiation. AGR-5/6/7 was conducted in the Advanced Test Reactor (ATR) at Idaho National Laboratory (INL) in support of development and qualification of tri-structural isotropic (TRISO) low-enriched fuel for use in high temperature gas-cooled reactors. The objectives of the AGR-5/6/7 experiments are to: (i) irradiate reference-design fuel particles to support fuel qualification, (ii) establish operating margins for the fuel beyond normal operating conditions, and (iii) provide irradiated-fuel performance data and irradiated-fuel samples for post-irradiation examination (PIE) and safety testing. The test train contains five separate capsules that were independently controlled and monitored. Each capsule contains multiple 12.51-mm-long compacts filled with low enriched uranium carbide/oxide (UCO) TRISO fuel particles. The primary objective of the AGR-5/6 test (Capsules 1, 2, 4, and 5) is to verify successful performance of the reference-design fuel under normal operating conditions. The AGR-7 test (Capsule 3) was designed to explore fuel performance at higher temperatures to demonstrate the capability of the fuel to withstand conditions beyond normal operating conditions in support of plant design and licensing. AGR 5/6/7 will also provide irradiated-fuel performance data on fission-gas release from failed particles during irradiation. The AGR-5/6/7 capsules were irradiated in the ATR northeast flux trap location. The experiment began on February 16, 2018 and ended on July 22, 2020, spanning nine ATR cycles over two and a half years. Thus, the AGR-5/6/7 fuel compacts were irradiated for a total of 360.9 effective full power days. The AGR 5/6/7 experiment was able to remain in the reactor core during all three Powered Axial Locator Mechanism (PALM) cycles (163A, 165A, and 167A) without overheating its fuel compacts. This report includes irradiation monitoring data from nine ATR Cycles: 162B, 163A, 164A, 164B, 165A, 166A, 166B, 167A, and 168A, as stored in the Nuclear Data Management and Analysis System (NDMAS). During irradiation, data records consisted of instantaneous measurements recorded every minute and provided by text files automatically every 2 hours. The AGR 5/6/7 data streams addressed in this report include thermocouple (TC) temperatures, sweep gas data (flow rates [capsule inlet, outlet, and downstream at detector], pressure, and moisture content), and Fission Product Monitoring System (FPMS) data (release rates and release to birth rate ratios [R/Bs]) for each of the five capsules. A total of 94,989,908 TC temperature and sweep gas data records were received and processed by NDMAS for AGR 5/6/7 irradiation. Of these records, 41,593,387 (or 43.7% of the total) met data collection and accuracy requirements and are labeled as Qualified. A total of 57,746,693 TC temperature readings were captured from 54 installed TCs. Among them, 10,034,676 TC temperature records (only 17.4%) were Qualified and 47,701,371 TC temperatures (or 82.6%) are Failed due to 48 TC failures (63.5%) and due to missing values (19.1%). To assess performance of the operational TCs, analysis of daily correlations between TCs found no evidence of virtual junction failure for any TCs. Analyses on control charts of TC temperature differences revealed trending in TC readings for TC2, 4, 5, and 13 in Capsule 3, but there is no conclusive indication of TC drift failure that caused those trends. Therefore, TC control charts are not used to disqualify TC data, but only for users’ consideration. For sweep gas flow rates, a total of 31,519,747 gas flow records (84.4%) are Qualified for use for AGR-5/6/7 experiment; 5,723,468 gas flow records (15.4%) are Failed due mostly to missing values; and 74,641 high sweep gas flow rates (0.2 %) are Trend. A large number of Failed missing TC temperature and gas flow values were caused by an error in the data output script that outputted a ‘NULL’ value when values were unchanged. This problem was fixed during the outage of Cycle 166B, which led to a substantially decreased number of missing values during the last three cycles. Nonetheless, a large amount of non-missing data remained because of the high data acquisition frequency (1-minute) and still provided sufficient data to effectively monitor the experiment as designed. For FPMS data, NDMAS received and processed fission product release and R/B data for nine ATR cycles, when ATR core reached full power during AGR 5/6/7 irradiation. These data consist of 110,388 release rate records and 110,388 R/B records for the twelve radionuclides (Kr 85m, Kr 87, Kr 88, Kr 89, Kr 90, Xe 131m, Xe 133, Xe 135, Xe 135m, Xe 137, Xe 138, and Xe 139) for each of the five capsules. Equivalent numbers of uncertainty records associated the release rates and R/B values were provided. To date, qualification status of the FPMS data stored in the NDMAS dat

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2019

Experimental Investigation of Forced Convection and Natural Circulation Cooling of a VHTR Core under Normal Operation and Accident Scenarios

This NEUP project #15-8205 entitled “Experimental Investigation of Forced Convection and Natural Circulation Cooling of a VHTR Core under Normal Operation and Accident Scenarios”, has investigated important phenomena in the safety of Very High Temperature Reactors (VHTRs) using the experimental facilities at City College of New York (CCNY) and Kansas State University (KSU), in collaboration with the Idaho State University (ISU) and Idaho National Laboratory (INL). The main objectives of this project were to investigate the forced convection and bypass flow phenomena, natural circulation flow and heat transfer, and graphite oxidation due to air ingress which could occur in a Very High Temperature Reactor (VHTR) with a prismatic core. High pressure/high temperature experiments have been conducted to obtain data that could be used for validation of VHTR design and safety analysis codes. The focus of these experiments was on the generation of benchmark data for design and off-design heat transfer for forced and natural circulation in a VHTR core. This work has built on the previous NEUP project #11-3218 which investigated forced convection and flow laminarization phenomena in a prismatic core of a VHTR and natural circulation of a pure gas under Pressurized Conduction Cooldown (PCC) and Depressurized Conduction Cooldown (DCC) conditions. Under each Task, extensive literature reviews on natural circulation, bypass flow, air ingress and graphite oxidation phenomena were conducted. In Task 1, scaling analyses were conducted by Idaho State University followed by design and construction of experimental facilities needed to conduct both bypass flow and natural circulation experiments at CCNY. Natural circulation experiments were conducted in Task 2 to investigate the air ingress phenomena, and bypass flow phenomena in Task 3, both at CCNY. In Task 4, radiation heat transfer and graphite oxidation tests were performed at Kansas State University. The experimental data were then analyzed to obtain quantitative results on natural circulation of pure gases and gas mixtures under DCC conditions, bypass flow rates, radiation heat transfer and graphite oxidation. In addition, CFD and Direct Numerical Simulations were conducted by INL and CCNY to investigate the flow laminarization phenomenon in strongly heated gas flows. The scaling analyses have shown that the key nondimensional numbers and parameters all have acceptable distortion levels principally because the CCNY facility’s test sections have representative cooling channel dimensions and are using either helium or a scalable gas such as nitrogen at representative state conditions. The natural circulation experiments involving helium and nitrogen gas mixtures at CCNY have elucidated the effect of nitrogen (simulating air) ingress into a lower plenum on the onset of natural circulation, mass flow rates and heat transfer coefficients under different power input levels in the riser section. Also, the transport of nitrogen from the lower plenum to upper plenum due to natural circulation was quantified for the first time. Nusselt number correlations for the heat transfer coefficient under natural circulation as a function of Reynolds and Prandtl numbers was also developed. The experiments at KSU revealed that prior to the onset of natural circulationin a VHTR core. This work has built on the previous NEUP project #11-3218 which investigated forced convection and flow laminarization phenomena in a prismatic core of a VHTR and natural circulation of a pure gas under Pressurized Conduction Cooldown (PCC) and Depressurized Conduction Cooldown (DCC) conditions. Under each Task, extensive literature reviews on natural circulation, bypass flow, air ingress and graphite oxidation phenomena were conducted. In Task 1, scaling analyses were conducted by Idaho State University followed by design and construction of experimental facilities needed to conduct both bypass flow and natural circulation experiments at CCNY. Natural circulation experiments were conducted in Task 2 to investigate the air ingress phenomena, and bypass flow phenomena in Task 3, both at CCNY. In Task 4, radiation heat transfer and graphite oxidation tests were performed at Kansas State University. The experimental data were then analyzed to obtain quantitative results on natural circulation of pure gases and gas mixtures under DCC conditions, bypass flow rates, radiation heat transfer and graphite oxidation. In addition, CFD and Direct Numerical Simulations were conducted by INL and CCNY to investigate the flow laminarization phenomenon in strongly heated gas flows. The scaling analyses have shown that the key nondimensional numbers and parameters all have acceptable distortion levels principally because the CCNY facility’s test sections have representative cooling channel dimensions and are using either helium or a scalable gas such as nitrogen at representative state conditions.The natural circulation experiments involving helium and nitrogen gas mixtures at CCNY have elucidated the effect of nitrogen (simulating air) ingress into a lower plenum on the onset of natural circulation, mass flow rates and heat transfer coefficients under different power input levels in the riser section. Also, the transport of nitrogen from the lower plenum to upper plenum due to natural circulation was quantified for the first time. Nusselt number correlations for the heat transfer coefficient under natural circulation as a function of Reynolds and Prandtl numbers was also developed. The experiments at KSU revealed that prior to the onset of natural circulation(ONC), molecular diffusion plays a significant role in air-ingress, however, at higher temperatures convection currents may be influencing air ingress and in-turn ONC times. Bypass flow experiments faced difficulties in the measurement of gas velocity in the main channel due to the failure of a high temperature Hot Wire Anemometer probe, so an energy balance method was instead employed to estimate the flow distribution between the main and bypass flow channels. In the graphite oxidation experiments conducted at KSU, oxidation was found to increase the thermal diffusivity of nuclear grade graphite. This change is believed to be driven by oxidation decreasing graphite sample density rather than raising the thermal conductivity. Oxidation was also shown to greatly alter graphite emissivity. If large amounts of impurities have built up within the graphite over time and have significantly increased its effective ash content, then air ingress is likely to severely inhibit the radiative heat transfer abilities of the graphite due to the formation of low emissivity, ash layers on its surface. However, if the purity of the graphite at the time of the air ingress is still close to that of its virgin specifications, the oxidation induced roughening of the graphite’s outer surface would have the effect of slightly raising its emissivity. The flow laminarization phenomenon was investigated by INL by examining the results of Direct Numerical Simulations of previous experiments, and by CCNY conducting its own DNS simulations using DOE’s Nek5000 code

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2016

On the reversibility of the Meissner effect and the angular momentum puzzle

It is generally believed that the laws of thermodynamics govern superconductivity as an equilibrium state of matter, and hence that the normal-superconductor transition in a magnetic field is reversible under ideal conditions. Because eddy currents are generated during the transition as the magnetic flux changes, the transition has to proceed infinitely slowly to generate no entropy. Experiments showed that to a high degree of accuracy no entropy was generated in these transitions. However, in this paper we point out that for the length of times over which these experiments extended, a much higher degree of irreversibility due to decay of eddy currents should have been detected than was actually observed. We also point out that within the conventional theory of superconductivity no explanation exists for why no Joule heat is generated in the superconductor to normal transition when the supercurrent stops. In addition we point out that within the conventional theory of superconductivity no mechanism exists for the transfer of momentum between the supercurrent and the body as a whole, which is necessary to ensure that the transition in the presence of a magnetic field respects momentum conservation. We propose a solution to all these questions based on the alternative theory of hole superconductivity. The theory proposes that in the normal-superconductor transition there is a flow and backflow of charge in direction perpendicular to the phase boundary when the phase boundary moves. We show that this flow and backflow explains the absence of Joule heat generated by Faraday eddy currents, the absence of Joule heat generated in the process of the supercurrent stopping, and the reversible transfer of momentum between the supercurrent and the body, provided the current carriers in the normal state are holes. - Highlights: • The normal-superconductor phase transition is reversible. • Within the conventional theory, Foucault currents give rise to irreversibility. • To suppress Foucault currents, charge has to flow in direction perpendicular to the phase boundary. • The charge carriers have to be holes. • This solves also the angular momentum puzzle associated with the Meissner effect.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2024

Revisiting a Drag Partition Model For Canopy-Like Roughness Elements

Turbulent flows over a large surface area (S) covered by n obstacles experience an overall drag due to the presence of the ground and the protruding obstacles into the flow. The drag partition between the roughness obstacles and the ground is analyzed using an analytical model proposed by Raupach (Boundary-Layer Meteorol 60:375-395, 1992) and is hereafter referred to as R92. The R92 is based on the premise that the wake behind an isolated roughness element can be described by a shelter area A and a shelter volume V. The individual sizes of A and V without any interference from other obstacles can be determined from scaling analysis for the spread of wakes. To upscale from an individual roughness element to n/S elements where wakes may interact, R92 adopted a background stress re-normalizing instead of reducing A or V with each element addition. This work demonstrates that R92’s approach results in a linear background stress reduction in A and V only when the ratio of n/S is small, due to a low probability of wake interactions. This probabilistic nature suggests that up-scaling from individual to multiple roughness elements can be re-formulated using stochastic averaging methods proposed here. The two approaches are shown to recover R92 under plausible conditions. An alternative scaling for the shelter volume is also proposed here using thermodynamic arguments of work and dissipation though the final outcome remains similar to R92. Comparisons between R92 and available data spanning more than two decades after R92 on blocks and vegetation-like roughness elements confirm the practical utility of R92. The agreement between R92 and this updated databases of experiments and simulations confirm the potential use of R92 in large-scale models provided that the relevant parameters accommodate certain features of the roughness element type (cube versus vegetation-like) and, to a lesser extent, their configuration throughout S. Last, a comparison between R92 and models based on first-order closure principles with constant mixing length suggests that R92 can outperform such models when evaluated across a wide range of roughness densities.

canopy turbulence↗

DOE OSTI2020

Active- and transfer-learning applied to microscale-macroscale coupling to simulate viscoelastic flows

Active- and transfer-learning are applied to microscale dynamics of polymer flows for the multiscale discovery of effective constitutive approximations required in viscoelastic flow simulation. The result is macroscopic rheology directly connected to a microstructural model. Micro and macroscale simulations are adaptively coupled by means of Gaussian process regression (GPR) to run the expensive microscale computations only as necessary. This multiscale method is demonstrated with flows of a polymer solution as a model system. At the microscale level dissipative particle dynamics (DPD) is employed to model the fluid as a suspension of bead-spring micro-structures subjected to steady shear flow. The results yield the non-Newtonian viscosity and the first normal stress difference at strain rates as training data used in a GPR model. DPD parameters are calibrated with respect to experimental data for a real polymer solution. Compliance with these data requires adjustment of the DPD model's cutoff radius, which then becomes a function of the second invariant of the strain rate tensor. The FENE-P model is chosen for the macroscale description using the spectral element method (SEM) to simulate channel flow and flow past a circular cylinder. The DPD results at the lowest possible shear strain rate yield an estimate of the zero-shear rate viscosity, which allows the initiation of the macroscale flow by SEM as a Newtonian fluid. The resulting strain-rate field is surveyed to determine additional shear strain rate sampling points for the DPD system. This new information allows an initial fitting of parameters of the constitutive equation followed by new SEM simulations at the macroscale. Additionally, guided by active-learning GPR to select new sampling points, this process continues until convergence is achieved. The effectiveness of this new simulation paradigm for viscoelastic flows is tested with different macroscale operating conditions. The effective closure learned in the channel simulation is then transferred directly to the flow past a circular cylinder at low Reynolds number, where the results show that only two additional DPD simulations are required to achieve a satisfactory constitutive model. With an increase of the Reynolds number, the active-learning scheme automatically detects the inaccuracy of the learned constitutive model, and initiates additional DPD simulations for the extra data needed to once again close the microscale-macroscale coupled system. This new paradigm of active- and transfer-learning for multiscale modeling is readily applicable to other microscale-macroscale coupled simulations of complex fluids and other materials. Furthermore, the coupling between microscale and macroscale solvers can be seamlessly implemented with our open source multiscale universal interface (MUI) library.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2012

Advanced High-Temperature Reactor Dynamic System Model Development: April 2012 Status

The Advanced High-Temperature Reactor (AHTR) is a large-output fluoride-salt-cooled high-temperature reactor (FHR). An early-phase preconceptual design of a 1500 MW(e) power plant was developed in 2011 [Refs. 1 and 2]. An updated version of this plant is shown as Fig. 1. FHRs feature low-pressure liquid fluoride salt cooling, coated-particle fuel, a high-temperature power cycle, and fully passive decay heat rejection. The AHTR is designed to be a “walk away” reactor that requires no action to prevent large off-site releases following even severe reactor accidents. This report describes the development of dynamic system models used to further the AHTR design toward that goal. These models predict system response during warmup, startup, normal operation, and limited off-normal operating conditions. Severe accidents that include a loss-of-fluid inventory are not currently modeled. The scope of the models is limited to the plant power system, including the reactor, the primary and intermediate heat transport systems, the power conversion system, and safety-related or auxiliary heat removal systems. The primary coolant system, the intermediate heat transport system and the reactor building structure surrounding them are shown in Fig. 2. These systems are modeled in the most detail because the passive interaction of the primary system with the surrounding structure and heat removal systems, and ultimately the environment, protects the reactor fuel and the vessel from damage during severe reactor transients. The reactor silo also plays an important role during system warmup. The dynamic system modeling tools predict system performance and response. The goal is to accurately predict temperatures and pressures within the primary, intermediate, and power conversion systems and to study the impacts of design changes on those responses. The models are design tools and are not intended to be used in reactor qualification. The important details to capture in the primary system relate to flows within the reactor vessel during severe events and the resulting temperature profiles (temperature and duration) for major components. Critical components include the fuel, reactor vessel, primary piping, and the primary-to-intermediate heat exchangers (P-IHXs). The major AHTR power system loops are shown in Fig. 3. The intermediate heat transfer system is a group of three pumped salt loops that transports the energy produced in the primary system to the power conversion system. Two dynamic system models are used to analyze the AHTR. A Matlab/Simulink-based model initiated in 2011 has been updated to reflect the evolving design parameters related to the heat flows associated with the reactor vessel. The Matlab model utilizes simplified flow assumptions within the vessel and incorporates an empirical representation of the Direct Reactor Auxiliary Cooling System (DRACS). A Dymola/Modelica model incorporates a more sophisticated representation of primary coolant flow and a physics-based representation of the three-loop DRACS thermal hydraulics. This model is not currently operating in a fully integrated mode. The Matlab model serves as a prototype and provides verification for the Dymola model, and its use will be phased out as the Dymola model nears completion. The heat exchangers in the system are sized using spreadsheet-based, steady-state calculations. The detail features of the heat exchangers are programmed into the dynamic models, and the overall dimensions are used to generate realistic plant designs. For the modeling cases where the emphasis is on understanding responses within the intermediate and primary systems, the power conversion system may be modeled as a simple boundary condition at the intermediate-to-power conversion system heat exchangers.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2012

Experimental study of heat transfer in a 7-element bundle cooled with supercritical Freon-12

Experimental data on Supercritical-Water (SCW) cooled bundles are very limited. Major problems with performing such experiments are technical difficulties in testing and experimental costs at high pressures, temperatures and heat fluxes. Also, there are only a few SCW experimental setups currently in the world capable of providing data. Supercritical Water-cooled nuclear Reactors (SCWRs), as one of the six concepts of Generation IV reactors, cannot be designed without such data. Therefore, a preliminary approach uses modeling fluids such as carbon dioxide and refrigerants instead of water is practical. In particularly, experiments in supercritical refrigerant-cooled bundles can be used. One of the SC modeling fluids typically used is Freon-12 (R-12) with the critical pressure of 4.136 MPa and the critical temperature of 111.97 deg. C. These conditions correspond to the critical pressure of 22.064 MPa and critical temperature of 373.95 deg. C in water. A set of experimental data obtained at the Inst. of Physics and Power Engineering (IPPE, Obninsk, Russia) in a vertically-oriented bundle cooled with supercritical R-12 was analyzed. This dataset consisted of 20 runs. The test section was 7-element bundle installed in a hexagonal flow channel with 3 grid spacers. Data was collected at pressures of approximately 4.65 MPa for several different combinations of wall and bulk-fluid temperatures that were below, at, or above the pseudo-critical temperature. The values of mass flux were ranged from 400 to 1320 kg/m{sup 2}s and inlet temperatures ranged from 72 to 120 deg. C. The test section consisted of fuel-element simulators that were 9.5 mm in OD with the total heated length of about 1 m. Bulk-fluid and wall temperature profiles were recorded using a combination of 8 different thermocouples. Analysis of the data has confirmed that there are three distinct heat-transfer regimes for forced convention in supercritical fluids: 1) Normal heat transfer; 2) Deteriorated heat transfer characterized with higher than expected wall temperatures; and 3) Enhanced heat transfer characterized with lower than expected wall temperatures. It was also confirmed that the effects of spacers are evident which was previously observed in sub-critical experimental data. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2014

PMU Data Integrity Evaluation through Analytics on a Virtual Test-Bed

Power systems are rapidly becoming populated by phasor measurement units (PMUs) in ever increasing numbers. PMUs are critical components of today s energy management systems, designed to enable near real-time wide area monitoring and control of the electric power system. They are able to measure highly accurate bus voltage phasors as well as branch current phasors incident to the buses at which PMUs are equipped. Synchrophasor data is used for applications varying from state estimation, islanding control, identifying outages, voltage stability detection and correction, disturbance recording, and others. However, PMU-measured readings may suffer from errors due to meter biases or drifts, incorrect configurations, or even cyber-attacks. Furthermore, the testing of early PMUs showed a large disparity between the reported values from PMUs provided by different manufacturers, particularly when frequency was off-nominal, during dynamic events, and when harmonic/inter-harmonic content was present. Detection and identification of PMU gross measurement errors are thus crucial in maintaining highly accurate phasor readings throughout the system. In this paper, we present our work in conducting analytics to determine the trustworthiness and worth of the PMU readings collected across an electric network system. By implementing the IEEE 118 bus test case on a virtual test bed (VTB) , we are able to emulate PMU readings (bus voltage and branch current phasors in addition to bus frequencies) under normal and abnormal conditions using (virtual) PMU sensors deployed across major substations in the network. We emulate a variety of failures such as bus, line, transformer, generator, and/or load failures. Data analytics on the voltage phase angles and frequencies collected from the PMUs show that specious (or compromised) PMU device(s) can be identified through abnormal behaviour by comparing the trend of its frequency and phase angle reading with the ensemble of all other PMU readings in the network. If the reading trend of a particular PMU deviates from the weighted average of the reading trends of other PMUs at nearby substations, then it is likely that the PMU is malfunctioning. We assign a weight to each PMU denoting how electric-topology-wise close it is from where the PMU under consideration is located. The closer a PMU is, the higher the weight it has. To compute the closeness between two nodes in the power network, we employ a form of the resistance distance metric. It computes the electrical distance by taking into consideration the underlying topology as well as the physical laws that govern the electrical connections or flows between the network components. The detection accuracy of erroneous PMUs should be improved by employing this metric. We present results to validate the proposed approach. We also discuss the effectiveness of using an end-to-end VTB approach that allows us to investigate different types of failures and their responses as seen by the ensemble of PMUs. The collected data on certain types of events may be amenable to certain types of analysis (e.g., alerting for sudden changes can be done on a small window of data) and hence determine the data analytics architectures is required to evaluate the streaming PMU data.

and data analytics architecture.↗

DOE OSTI2023

How does uncertainty of soil organic carbon stock affect the calculation of carbon budgets and soil carbon credits for croplands in the U.S. Midwest?

Cropland carbon budget depicts the amount of carbon flowing in and out of agroecosystems and the changes in carbon stocks of soil and living biomass during the same period. Soil carbon credit is the additional change in soil carbon stock under certain farming practices compared with the business-as-usual practices. Accurately calculating cropland carbon budget and soil carbon credit is critical to assessing climate change mitigation potential in agroecosystems. The calculation of cropland carbon budget and soil carbon credit is sensitive to local soil and climatic conditions, especially initial soil organic carbon (SOC) stock, which is determined by both SOC concentration (SOC%) and bulk density (Bulk_Density). SOC stock data are either from soil sampling or gridded public survey data. In agroecosystem models, SOC stock data are a key model input for quantifying cropland carbon budget and soil carbon credit. However, various types and degrees of uncertainties exist in SOC stock datasets, which propagate to the quantification of SOC stock change. In particular, a large discrepancy is found in two widely used SOC stock datasets — Rapid Carbon Assessment dataset (RaCA) and Gridded Soil Survey Geographic Database (gSSURGO) — in the U.S. Midwest, with a relative difference (quantified using Normalized Root Mean Square Error, NRMSE) of 48.0% for 0–30 cm SOC stock between the two datasets. It remains largely unclear how uncertainty in SOC stocks affects the calculation of cropland carbon budget and soil carbon credit. To address this question, we used a well-validated process-based agroecosystem model, ecosys, to assess the impacts of SOC stock uncertainty on carbon budget and soil carbon credit calculation in the U.S. Midwestern corn-soybean rotation systems. Our results reveal the following findings: (1) A sizable discrepancy exists in simulated cropland carbon budget between using gSSURGO and using RaCA for their SOC% and Bulk_Density as model inputs, with a Pearson correlation coefficient (r) of only 0.4 for simulated change of SOC stock (ΔSOC) using these two different soil datasets. (2) Simulated cropland carbon budget components were more sensitive to initial SOC% than to Bulk_Density. For example, the upper and lower quartiles of multi-year averaged ΔSOC were –29.8 and 4.8 gC/m 2 /year for the selected counties respectively, with an uncertainty of 13.7 and 0.7 gC/m 2 /year induced by uncertainties in initial SOC% and Bulk_Density, respectively. (3) Both simulated ΔSOC and its uncertainty were negatively correlated with initial SOC%, whereas ΔSOC was negatively correlated with air temperature, and ΔSOC uncertainty was positively correlated with air temperature. (4) The uncertainty of calculated soil carbon credits was much smaller compared with the uncertainty of calculated absolute carbon budgets assuming the same SOC stock uncertainty level in the inputs. Specifically, in our assessment comparing planting cover crops vs no cover crop, the uncertainty of calculated soil carbon credits induced by initial SOC% uncertainty was less than 4% (relative to the quantified value of the soil carbon credits) for 90% of the cases. Our analysis highlights that high accuracy measurement of SOC% as inputs is needed for the calculation of cropland carbon budgets; however, soil carbon credit quantification is much less sensitive to the initial SOC% inputs, and the current publicly available soil datasets (e.g., gSSURGO) are largely suitable for the calculation of soil

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2012

Review of the WECC EDT phase 2 EIM benefits analysis and results report.

A region-wide Energy Imbalance Market (EIM) was recently proposed by the Western Electricity Coordinating Council (WECC). In order for the Western Area Power Administration (Western) to make more informed decisions regarding its involvement in the EIM, Western asked Argonne National Laboratory (Argonne) to review the EIM benefits study (the October 2011 revision) performed by Energy and Environmental Economics, Inc. (E3). Key components of the E3 analysis made use of results from a study conducted by the National Renewable Energy Laboratory (NREL); therefore, we also reviewed the NREL work. This report examines E3 and NREL methods and models used in the EIM study. Estimating EIM benefits is very challenging because of the complex nature of the Western Interconnection (WI), the variability and uncertainty of renewable energy resources, and the complex decisions and potentially strategic bidding of market participants. Furthermore, methodologies used for some of the more challenging aspects of the EIM have not yet matured. This review is complimentary of several components of the EIM study. Analysts and modelers clearly took great care when conducting detailed simulations of the WI using well-established industry tools under stringent time and budget constraints. However, it is our opinion that the following aspects of the study and the interpretation of model results could be improved upon in future analyses. The hurdle rate methodology used to estimate current market inefficiencies does not directly model the underlying causes of sub-optimal dispatch and power flows. It assumes that differences between historical flows and modeled flows can be attributed solely to market inefficiencies. However, flow differences between model results and historical data can be attributed to numerous simplifying assumptions used in the model and in the input data. We suggest that alternative approaches be explored in order to better estimate the benefits of introducing market structures like the EIM. In addition to more efficient energy transactions in the WI, the EIM would reduce the amount of flexibility reserves needed to accommodate forecast errors associated with variable production from wind and solar energy resources. The modeling approach takes full advantage of variable resource diversity over the entire market footprint, but the projected reduction in flexibility reserves may be overly optimistic. While some reduction would undoubtedly occur, the EIM is only an energy market and would therefore not realize the same reduction in reserves as an ancillary services market. In our opinion the methodology does not adequately capture the impact of transmission constraints on the deployment of flexibility reserves. Estimates of flexibility reserves assume that forecast errors follow a normal distribution. Improved estimates could be obtained by using other probability distributions to estimate up and down reserves to capture the underlying uncertainty of these resources under specific operating conditions. Also, the use of a persistence forecast method for solar is questionable, because solar insolation follows a deterministic pattern dictated by the sun's path through the sky. We suggest a more rigorous method for forecasting solar insolation using the sun's relatively predictable daily pattern at specific locations. The EIM study considered only one scenario for hydropower resources. While this scenario is within the normal range over the WI footprint, it represents a severe drought condition in the Colorado River Basin from which Western schedules power. Given hydropower's prominent role in the WI, we recommend simulating a range of hydropower conditions since the relationship between water availability and WI dispatch costs is nonlinear. Also, the representation of specific operational constraints faced by hydropower operators in the WI needs improvements. The model used in the study cannot fully capture all of the EIM impacts and complexities of power system operations. In particular, a primary benefit of the EIM is a shorter dispatch interval; namely, 5 minutes. However, the model simulates the dispatch hourly. Therefore it cannot adequately measure the benefits of a more frequent dispatch. A tool with a finer time resolution would significantly improve simulation accuracy. When the study was conducted, the rules for the EIM were not clearly defined and it was appropriate to estimate societal benefits of the EIM assuming a perfect market without a detailed specification of the market design. However, incorporating a more complete description of market rules will allow for better estimates of EIM benefits. Furthermore, performing analyses using specific market rules can identify potential design flaws that may be difficult and expensive to correct after the market is established. Estimated cost savings from a more efficient dispatch are less than one percent of the total cost of electricity production.

14 SOLAR ENERGY↗

DOE OSTI2025

Operation of helium sub-atmospheric multistage cryogenic centrifugal compressor trains: Part 2 – Transient modeling and pump-down path selection

Low-pressure conditions required for operation of helium cryogenic systems below the normal boiling point (i.e. 4.2 K) are established through a transient process, commonly referred to as ‘pump-down’. This process is defined as the transition from pressures above atmospheric conditions to the saturation pressure which corresponds to a specified operational temperature. The FRIB 2 K system consists of five cryogenic centrifugal compressors which are operated in series. Historically, the pump-down process path has been established through empirical methods and system operator experience. Investigation into the pump-down process at FRIB aimed to develop a pump-down methodology which relies on theoretical model predictions rather than empirically developed process paths. Ensuring stable operation during the pump-down process involved application of a centrifugal compressor performance prediction model, which is described in Part 1 of this paper. Compressor performance maps can be directly used to evaluate the stability of a selected pump-down path and anticipate the overall reliability of the selected path. In conjunction with the compressor performance maps, a system pressure model was developed to estimate the transient pressure response during the pump-down process. Lastly, an explicit equation was developed to establish a mass flow rate profile for the pump-down process. Implementation of the presented methodology (including the developed models) allows for the system operator to determine a continuous pump-down path which maintains compressor stability while conforming to overall system capabilities. Altogether, the methodology presented has resulted in simplification of transient pump-down operations and increased the reliability, stability and efficiency of the pump-down process.

Compressor train control↗

DOE OSTI2019

Enabling Staged Pressurized Oxy-Combustion: Improving Flexibility and Performance at Reduced Cost

Introduction: The staged pressurized oxy-combustion (SPOC) process, developed by Washington University in St. Louis (WUSTL), is an oxy-combustion carbon capture technology that employs a steam power cycle to generate electricity. As with atmospheric-pressure oxy-combustion processes, carbon dioxide (CO 2 ) capture rates of 90% or higher can be achieved. The staged-combustion approach of the SPOC process, when operated at elevated gas pressure, allows for near elimination of flue gas recycle (FGR) and leads to significant improvements in efficiency and reduced costs. The primary goal of this project was to investigate the potential for SPOC for flexible operation beyond the capabilities of conventional coal-fired power plants, particularly those employing carbon capture and storage (CCS). Oxy-combustion differs from conventional coal combustion in that the combustion of coal is carried out using oxygen as opposed to air. The resultant flue gas is a mixture of primarily CO 2 and water, greatly simplifying CO 2 capture. As SPOC does the oxy-combustion process under pressure, it is higher efficiency and lower cost than atmospheric oxy-combustion. The project was led by the Electric Power Research Institute, Inc., with the assistance of WUSTL, Doosan Babcock (DBL – boiler plant original equipment manufacturer), and Air Liquide (AL – air separation experts). The specific objectives of the project were to: (1) Evaluate the SPOC concept and develop a risk-based approach to the heating surface layout ensuring that performance (gas and steam side), manufacturing, transportation, and plant erection considerations were fully accounted for in the system design. (2) Improve the technology to ensure its performance and cost potential are substantially better than today’s baseline pulverized coal power plant with post-combustion capture (PCC) or atmospheric oxy-combustion, and show progress toward performance commensurate with projected commercial operation, including 90% or more CO 2 capture. (3) Address critical technology gaps and improving overall system performance for the technology. (4) Perform combustion tests at scale under commensurate pressure to the commercial operating system validate combustor and advance the SPOC combustion modeling tools that will facilitate full-scale design. Combustion Testing: Testing was initially carried out at atmospheric pressure with methane to ensure that all systems were correctly operating. Following this, the testing then proceeded to moderate pressure operation at 4 bara, where different stoichiometric conditions were used across an extended heat input range, showing that a stable flame could be maintained. Testing then proceeded at 10 and 15 bara where it was shown that the ignition system needed to be replaced at these operating conditions. When testing recommenced, a similar stable flame shape was achieved, without any methane support, and the full 100 kWth coal oxy-combustion was demonstrated. Heat flux measurements showed 400 kW/m 2 (0.126 MMBtu/hr-ft 2 ) from the flame, within the 450 kW/m 2 (0.142 MMBtu/hr-ft 2 ) preferred limit defined by DBL for the full-scale design case. Carbon monoxide and carbon-in-ash measurements showed that complete combustion was possible with ultra-low excess oxygen at 1 vol % in the product flue gas. This allowed the full-scale models to be calculated based on this level of excess oxygen, improving the performance of the system as lower feed oxygen is produced in the air separation unit (ASU), saving auxiliary power. Full-Scale Design: The full-scale design is sized to deliver 550 MWe net electrical power output to allow direct comparison with the National Energy Technology Laboratory (NETL) baseline cases using the same Montana Powder River Basin (PRB) fuel. The cases used for direct comparison are: Case S12A (supercritical [SC] coal without CO 2 capture), Case S12B (SC coal with PCC at 90% CO 2 capture, and Case S12F (atmospheric oxy-combustion SC coal with 90% CO 2 capture). The SPOC stages were configured to be identical arrangements, allowing for a more straightforward design, construction, and control strategy. Each stage consists of two PVs, the first being the combustor PV module and the second being the convective PV module. The combustor module is the tallest vessel, with local coal feeding equipment and the burner FGR and oxygen plenums at the top and ash management at the bottom. The combustor module is effectively an open cross section with membrane tubes around the circumference – this unobstructed slender profile ensures that the relatively slower burnout mechanism needed for the low-speed turbulent mixing burner can proceed unimpeded and that molten ash will not attach to any surfaces at an angle to the predominant flue gas flow. The steam turbine used in the SPOC design is essentially the same unit utilized in the atmospheric oxy-combustion case except the degree of heat recovery possible allows the low-pressure (LP) feedwater heater train to be eliminated. This is a result of being able to capture the latent heat of evaporation of the moisture produced in the combustion (normally lost to the stack). The SPOC system outperforms the NETL atmospheric oxy-combustion baseline case by 3.3% points and the PCC case by 7.5% points. This performance increase is mainly a result of the improved heat recovery possible from the product flue gas and from the reduced auxiliary power consumption from the ASU and the CO 2 compression and purification unit over the atmospheric equivalent. The combined improvement in boiler efficiency and reduction in the steam turbine heat rate resulted in a net improvement in the overall plant efficiency by 0.61%, to 35.6% on a HHV basis. This relatively small improvement in plant efficiency between sub-bituminous and bituminous coals shows the value of useful moisture latent heat recovery offered by the SPOC system when using low rank fuels. Flexibility and Turndown: The SPOC system has inherent high turndown capabilities because of a demonstrated high range of combustion stability and the ability to bypass stages, thereby reducing steam generation rates proportionally without altering firing rates in operating stages. ASU units do not generally offer significant turndown opportunities as they are limited by the performance of the compression plant. The baseline SPOC design has a 2-train ASU configuration, needed to be able to meet the 10,500 tonnes per day duty, that can deliver efficient turndown points between 85–100% load with both trains operating, and at 45–50% load with a single train operating. Outside these operating points, the flexibility and efficiency are limited and may be uncompetitive at loads below 45% as the ASU systems are unable to turndown beyond 85% for the base case and 70% for the flexible case. The ASU cold box for the flexible case was designed to operate at loads between 40-100%, below this operating range the purity of the produced oxygen cannot be guaranteed, resulting in additional load for the CPU equipment. The additional power consumption of these cases yields an uneconomic operating condition with poor efficiency when operating at low loads such as 25% and 12%. To address this flexibility constraint, AL evaluated how the ASU configuration could be altered to deliver efficient oxygen supply over the load range. This involved integrating the compressor duties across both ASU trains into a shared duty arrangement with multiple compressors, strategically sized to deliver the duty required across all load ranges between 40–100% by enabling different combinations of the main air compressors and the booster air compressors. The oxygen purity cannot be guaranteed below 40% load for each train – the power consumption at 25% and 12% loads were estimated based on the multi-compressor configuration at maximum turndown as these load cases were not considered in the flexible ASU design. With a different cold box configuration, i.e. number of trains, the entire load range could be delivered within oxygen purity specifications. Although there is a small efficiency reduction at full load, there are significant improvements in overall plant performance during low-load operation. The exact ASU configuration needed will therefore be dependent on the expected operating profile of the unit. Economic Analysis: Cost estimates for the baseline and the flexible cases were developed using NETL baseline case data for common BOP items and an Association for the Advancement of Cost Engineering International Class 5 (conceptual/screening study) assessment of key components unique to the SPOC system. The capital, operating, and maintenance costs were assessed along with the first-year power cost, levelized cost of electricity, and CO 2 captured and avoided cost for the SPOC cases were compared against the relevant NETL baseline cases in January 2019 dollars. Both SPOC cases achieve a lower cost that the alternative NETL baseline cases (with the flexible SPOC case being slightly higher than the baseline SPOC case due to the compounded impact of higher capital costs and lower efficiency at full load). The SPOC cases deliver a lower cost of CO 2 avoided than the NETL baseline cases due to greater efficiency in generating low-carbon power. However, the cost of CO 2 captured is higher because the higher efficiency of the SPOC process generates less CO 2 , hence less needs to be captured, resulting in a high capital cost being shared over a reduced CO 2 quantity. The flexible ASU adds 1.6% overall to the plant cost in comparison to the baseline case, but the efficiency improvements at loads below 50% are significant. Subsequently, depending on the plant load-profile expectation, the flexible ASU would be beneficial if the plant spends a significant portion of its operating life below 50% load. This kind of operating profile is likely to be required for all fossil plants when more renewable electricity generators are installed on the local grid, particularly solar power that can be predicted ahead of time, allowing for appropriate pricing signals to be incorporated into the diurnal cycle. R&D Recommendations: SPOC, while a promising technology, is a relatively recent concept and, as such, operability issues of combustor design and steam-side integration for such systems at a scale relevant to commercial deployment have not been evaluated. WUSTL has conducted extensive small pilot-scale (100 kWth) research to understand and advance pressurized oxy-combustion processes, including investigation of combustion and flame characteristics, radiative heat flux, burner operability, turn down, char burnout, ash characteristics, water-wash column operation, etc. for pressurized oxy-combustion systems. Nonetheless, at this stage in the development of the SPOC process, what is needed is a large-scale pilot plant that can serve to study pressurized oxy-combustion systems and components at a scale commensurate with the maturity of the technology. A scale of 10 MWth, which includes steam-side integration and two stages would yield essential information with respect to heat transfer characteristics both in the radiative and convective sections of the pressure vessels, and the ability to operate the fuel staging process. In addition, while modeling results indicate that combustion and flame characteristics improve with scale, direct studies of the combustor at this scale will ensure that the models can be relied upon for scale up to commercial scale. Furthermore, a detailed analysis must be performed to understand the scaling aspects of key components and systems.

20 FOSSIL-FUELED POWER PLANTS↗

DOE OSTI2011

Reliable algorithms for power system analysis in the presence of data uncertainties

A robust and reliable power flow analysis represents an essential requirement for many power systems applications as far as network optimization, voltage control, state estimation, and service restoration are concerned. The most common power flow approach, referred to here as a deterministic power flow (PLF), requires precise or 'crisp' values chosen by the analyst for each input variable. The solution provides precise network voltages and flows through each line. The specified values rest upon assumptions about the operating condition derived from historical measurements or predictions about future conditions and thus, cannot be considered accurate. Even in the case where the inputs are based on measurements, inaccuracies arise from time-skew problems, three-phase unbalance, static modeling approximations of dynamic components (e.g., transformer tap changers), variations in line parameters, and so on. The advent of deregulation and competitive power markets will only exacerbate this problem as well-known generation patterns change, loading becomes less predictable and the transmission paths grow more diverse. Conventional methodologies proposed in literature address tolerance analysis of power flow solution by means of detailed probabilistic methods, accounting for the variability and stochastic nature of the input data, and sampling based approaches. In particular uncertainty propagation using sampling based methods, such as the Monte Carlo, requires several model runs that sample various combinations of input values. Since the number of model runs can sometimes be very large, the required computer resources can sometimes be prohibitively expensive resulting in substantial computational demands. As far as probabilistic methods are concerned, they represent a useful tool, especially for planning studies, but, as evidenced by the many discussions reported in literature, they could reveal some shortcomings principally arising from: (1) the non-normal distribution and the statistical dependence of the input data; and (2) the difficulty arising in accurately identifying probability distributions for some input data, such as the power generated by wind or photovoltaic generators. All these could result in time consuming computations with several limitations in practical applications especially in power flow analysis of complex power networks. In order to try and overcome some of these limitations, obtaining thereby comprehensive power flow solution tolerance analysis at adequate computational costs, self validated computation could play a crucial role. Armed with such a vision, this chapter will analyze two advanced techniques for power flow analysis in the presence of data uncertainty namely the boundary power flow and the affine arithmetic power flow.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2021

Experimental and Numerical Evaluation of the Thermal-Fluid Characteristics of the Helium-Cooled Modular Divertor with Multiple Jets Using a Reversed Heat Flux Approach

A single-finger unit of the Helium-Cooled Modular Divertor with Multiple Jets (HEMJ) with a plasma-facing surface (PFS) area of about 2 cm 2 has been studied in a helium (He) loop at He mass flow rates m˙ H ≤ 8 g/s and nearly prototypical conditions. Based on previous studies of the single finger of the HEMJ, our Georgia Institute of Technology group is planning to experimentally study larger divertors. Given that the HEMJ test section was heated with an induction heater and that it is impractical to scale this up to divertors with larger PFS areas, a reversed heat flux approach is being considered to measure heat transfer coefficients (HTCs). In this approach, the direction of the heat flux is reversed with water cooling and high-temperature He heating of the outer shell attached to the PFS. This work presents an initial experimental and numerical evaluation of this approach for a single HEMJ finger. Experiments with brass and copper-chromium-zirconium outer shells were conducted at dimensionless He mass flow rates or Reynolds numbers Re = 1 × 10 4 to 4.7 × 10 4 , an inlet pressure of 10 MPa, temperatures as great as 673 K, and maximum heat flux of 8.4 MW/m 2 . The experiments verify that the He-side HTCs are independent of the direction of the heat flux. Here, the results agree well with previous Nusselt number correlation and pressure loss coefficients for the HEMJ obtained using the normal heating approach.

42 ENGINEERING↗

DOE OSTI2024

Causes of and Solutions to Wind Speed Bias in NREL's 2020 Offshore Wind Resource Assessment for the California Pacific Outer Continental Shelf

This report provides the results of a detailed analysis into the causes of high wind speed bias in the 20-year wind resource data set for offshore California the National Renewable Energy Laboratory (NREL) released in 2020, herein called CA20. The data set was developed using the state-of-the-art Weather Research and Forecasting (WRF) model. Notably, no floating lidars were available at the time in offshore California to validate offshore hub-height wind speeds. In late 2020, the Pacific Northwest National Laboratory (PNNL) deployed two floating lidars in the California outer continental shelf (OCS), near the Bureau of Ocean Energy Management (BOEM) call areas of Humboldt and Morro Bay. Using these observations through 2021, NREL found considerable bias in modeled hub-height winds at both locations: up to +2 m/s at Humboldt over a 6-month period, and up to +1 m/s at Morro Bay over a one-year period. Upon the discovery of this bias, the Department of Energy (DOE) and BOEM funded NREL and PNNL to investigate the causes of, impacts of, and solutions to the bias in the CA20 data set. This report summarizes the findings of this research. We first investigated whether different WRF model setups could lead to reduced bias. We found that the choice of planetary boundary layer (PBL) scheme - which controls the vertical turbulent mixing of momentum, heat, and moisture in the lowermost part of the atmosphere - greatly affected hub-height wind speeds in the region. Specifically, switching from the Mellor-Yamada-Nakanishi-Niino (MYNN) scheme used in CA20 (and widely used across a range of operational and research weather models) to the less common Yonsei University (YSU) scheme nearly eliminated the bias at both the Humboldt and Morro Bay lidar locations. The large discrepancy between the MYNN- and YSU-modeled hub-height winds pointed towards the role of atmospheric stability. In general, PBL schemes agree well in conditions of high turbulence and mixing, normally referred to as "unstable" conditions. By contrast, PBL schemes start to diverge in "stable" conditions, where turbulence is low and thermal stratification (i.e., higher temperature air sitting on top of colder air) greatly suppresses vertical mixing. Under such conditions, winds aloft can decouple from surface effects and greatly accelerate, causing high wind speeds at hub-height and frequent low-level jets (LLJs). We determined that these stable conditions are in fact dominant in offshore California. The region is characterized by moderate-to-extreme stable stratification with a LLJ on average around 200 meters above sea-level. To our knowledge, no wind energy area globally has as strongly stable stratification as offshore California. Under these extreme conditions, we determined that the MYNN scheme models higher stability than YSU, resulting in less vertical turbulent mixing than YSU, allowing for the acceleration of hub-height winds, more intense LLJs, and higher-amplitude inertial oscillations. Using surface observations, we found that MYNN overestimates near-surface stability, whereas YSU tends to model stability better. We then considered several short-term case studies to assess additional meteorological drivers of the bias at Humboldt. We found that during synoptic scale northerly flows driven by the North Pacific High and inland thermal low, a coastal warm bias in the MYNN case studies contributes to the modeled wind speed bias by altering the boundary layer thermodynamics via a thermal wind mechanism. Given the strong performance of the YSU-based runs in offshore California, NREL has produced and published an updated version of the CA20 data set with YSU as the PBL scheme. This updated data set is now part of NREL's 2023 National Offshore Wind (NOW-23) data set, which covers all the U.S. offshore waters. The development and final validation of the NOW-23 data set in offshore California is documented in this report.

17 WIND ENERGY↗

DOE OSTI2019

Transition from ion-coupled to electron-only reconnection: Basic physics and implications for plasma turbulence

Here using 2.5 dimensional kinetic particle-in-cell simulations, we simulate reconnection conditions appropriate for the magnetosheath and solar wind, i.e., plasma beta (ratio of gas pressure to magnetic pressure) greater than 1 and low magnetic shear (strong guide field). Changing the simulation domain size, we find that the ion response varies greatly. For reconnecting regions with scales comparable to the ion inertial length, the ions do not respond to the reconnection dynamics leading to “electron-only” reconnection with very large quasisteady reconnection rates. Note that in these simulations, the ion Larmor radius is comparable to the ion inertial length. Moreover, the transition to a more traditional “ion-coupled” reconnection is gradual as the reconnection domain size increases, with the ions becoming frozen-in in the exhaust when the magnetic island width in the normal direction reaches many ion inertial lengths. During this transition, the quasisteady reconnection rate decreases until the ions are fully coupled, ultimately reaching an asymptotic value. The scaling of the ion outflow velocity with the exhaust width during this electron-only to ion-coupled transition is found to be consistent with a theoretical model of a newly reconnected field line. In order to have a fully frozen-in ion exhaust with ion flows comparable to the reconnection Alfvén speed, an exhaust width of at least several ion inertial lengths is needed. In turbulent systems with reconnection occurring between magnetic bubbles associated with fluctuations, using geometric arguments, we estimate that fully ion-coupled reconnection requires magnetic bubble length scales of at least several tens of ion inertial lengths.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Utah FORGE: Fault Shear Reactivation Experimental Data for Fluid Injection-Rate Controls on Seismic Moment

Included are experimental data recorded from shear experiments that specifically explore the link between fluid-injection rate and seismic moment resulting from shear reactivation of laboratory faults. Raw mechanical data from three experiments are included alongside corresponding MATLAB scripts that import and plot the data, as well as use it to calculate shear and normal stress. Experiments are performed on 2.5-3 inch long granitoid cores from the Utah FORGE EGS demonstration site, containing a single inclined fracture with small-scale roughness added to the fracture surface. The raw data included here were recorded from an aluminum triaxial pressure vessel (TEMCO) configured with three independent servo-controlled pumps, with distilled water used as the working fluid. The pumps control confining pressure, upstream pore pressure, and axial pressure, with each pump connected to a LabView interface to record applied pressures, cumulative injected water volumes, and pump flow rates. The downstream outlet from the fracture is closed to allow pressurization, which is measured by an external pressure transducer. A linear variable differential transformer (LVDT) attached to the axial piston measures axial displacement, from which we calculate shear displacement along the fracture. Additionally, P-wave transducers are used to record acoustic signatures, where acoustic emission events and maximum amplitudes are compared against seismic moment and shear slip velocity. Fluid injection rates range between 0.05 mL/min, 0.25 mL/min, and 0.75 mL/min for each experiment. Along-fault pressure distributions are progressively less uniform as injection rates increase, representing a switch from steady-state to transient conditions. Triggered shear displacement is used as a proxy for seismic moment and is indexed against cumulative injection volume and rate. Each experiment is performed under constant shear stress conditions, and the sample is fully saturated with DI water. Axial and confining stresses are applied to 3 MPa through pressure-stepping in 500 kPa increments. The pore pressure is held constant at 200 kPa prior to initiating the experiment, and initial axial displacement is recorded. The axial stress is then increased to initiate shear mobilization during the loading phase (run-in) until a peak steady state is achieved. The initial shear stress is reduced to approximately 80% of the peak shear stress by decreasing the axial stress, then held constant for the duration of each experiment.

15 GEOTHERMAL ENERGY↗