Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reactive 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.

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

At least 649 records · Page 36

CO 2 -Responsive Fracturing Fluids for Enhanced Geothermal Systems (Final Report)

Our group has recently developed StimuFrac, a non-toxic stimuli-responsive fracturing fluid consisting of a CO 2 -reactive polymer which has shown at the lab-scale to consistently fracture rock cores at significantly lower net pressures in a range of representative geothermal pressure/temperature conditions. However, until now the mechanism/s responsible for more effective fracturing, of critical importance to optimize fracturing performance as well as strategize injection methodologies for field deployment, was not understood. In this document, we report (1) on the two main mechanisms responsible for fracturing rock at lower net pressures with StimuFrac; (2) the phase behavior of StimuFrac/CO 2 under geothermal wellbore conditions; and (3) based on high-temperature true triaxial stimulations, detailed evidence that StimuFrac/CO 2 is the best performing stimulation fluid under EGS T/P conditions as compared to water, waterless CO 2 , and CO 2 /water fracturing fluids. This is because i) it requires significantly lower volumes of CO 2 due to its reduced leak off into the formation by the crosslinked polymer solution; ii) large fractures can be generated reproducibly at both low and high CO 2 injection flow rates, and iii) the reversible (previously reported) viscosity increase of StimuFrac could be beneficial to transport proppants when they become available for EGS. These results were particularly evident for hot nearly dry rock as well as partially and nearly fully water-saturated granitic rock. GTO requires StimuFrac to be evaluated in fully water-saturated rock to determine whether the above-described performance applies under these conditions. Since (1) GTO considers there is no enough evidence for Sa different StimuFrac formulation to work below 90C (where full water-saturation in an open system is possible) and (2) no polyaxial loading frame larger than a few centimeters that can do hydraulic fracturing tests while maintaining the rock sample fully saturated with water at 200 C exists; PNNL concludes that the only way to determine if StimuFrac represents an advanced fracturing fluid alternative for EGS, is to perform a stimulation in an actual EGS reservoir.

15 GEOTHERMAL ENERGY↗

Heterolytic Scission of Hydrogen Within a Crystalline Frustrated Lewis Pair

We report the heterolysis of molecular hydrogen under ambient conditions by the crystalline frustrated Lewis pair (FLP) 1-{2-[bis(pentafluorophenyl)boryl]phenyl}-2,2,6,6-tetramethylpiperidine (KCAT). The gas–solid reaction provides an approach to prepare the solvent-free, polycrystalline ion pair KCATH2 through a single crystal to single crystal transformation. The crystal lattice of KCATH2 increases in size relative to the parent KCAT by approximately 2%. Microscopy was used to follow the transformation of the highly colored red/orange KCAT to the colorless KCATH2 over a period of 2 h at 300 K under a flow of H2 gas. There is no evidence of crystal decrepitation during hydrogen uptake. Inelastic neutron scattering employed over a temperature range from 4–200 K did not provide evidence for the formation of polarized H2 in a precursor complex within the crystal at low temperatures and high pressures. However, at 300 K, the INS spectrum of KCAT transformed to the INS spectrum of KCATH2. Calculations suggest that the driving force is more favorable in the solid state compared to the solution or gas phase, but the addition of H2 into the KCAT crystal is unfavorable. Ab Initio methods were used to calculate the INS spectra of KCAT, KCATH2, and a possible precursor complex of H2 in the pocket between the B and N of crystalline KCAT. Ex-situ NMR showed that the transformation from KCAT to KCATH2 is quantitative and our results suggest that the hydrogen heterolysis process occurs via H2 diffusion into the FLP crystal with a rate-limiting movement of H2 from inactive positions to reactive sites.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Optimization of artificial viscosity in production codes based on Gaussian Regression surrogate models

To accurately model flows with shock waves using staggered-grid Lagrangian hydrodynamics, artificial viscosity has to be introduced to convert kinetic energy into internal energy, thereby increasing the entropy across shocks. Determining the appropriate strength of the artificial viscosity is an art and strongly depends on the particular problem and experience of the researcher. The objective of this study is to pose the problem of finding the appropriate strength of artificial viscosity as an optimization problem and solve this problem using machine learning (ML) tools, specifically using surrogate models based on Gaussian Process regression and Bayesian analysis. We describe the optimization method and discuss various practical details of its implementation. The shock-containing problems for which we apply this method all have been implemented in the LANL code FLAG. First, we apply ML to find optimal values to isolated shock problems of different strengths. Second, we apply ML to optimize viscosity for a 1D propagating detonation problem based on Zel’dovich-von Neumann-Doring (ZND) detonation theory using a reactive burn model. We compare results for default (currently used values in FLAG) and optimized values of artificial viscosity for these problems demonstrating the potential for significant improvement in the accuracy of computations.

42 ENGINEERING↗

Analysis of Microseismicity and Reactivated Fault Size to Assess the Potential for Felt Events by CO 2 Injection in the Illinois Basin

The results of monitoring of carbon dioxide (CO 2 ) injection at the Illinois Basin—Decatur Project (IBDP) and the companion Illinois Industrial Carbon Capture and Sequestration Sources (IL-ICCS) project—have shown that reservoir response to fluid pressure changes can vary significantly at different injection locations within the same reservoir. Predrill reservoir characterization is important to identify potentially seismogenic faults. However, interpretations of newly reprocessed 3D seismic reflection data illustrate the challenges related to their identification in a region dominated by faulting with small vertical offsets. Faults interpreted in the 3D seismic volume range from ~300 to 1200 m wide and are in the same size range as faults that could have been the source of historical events up to Mw 2.7 in central Illinois. The array of monitoring sensors that was installed for the IBDP continues to collect data, as injection operates in IL-ICCS, the second injection well. CO 2 injection rates for the IL-ICCS well are on average 1.7 times the rates injected in the IBDP well, but a significantly reduced rate of induced seismicity is observed. This article presents results of passive seismic monitoring for the duration of the project to date, integrating active and passive seismic data to develop a new interpretation of the subsurface structure at the Decatur site that explicitly identifies pathways for fluid flow into the basement leading to induced seismicity, and provides a geological explanation for the sharp reduction of induced seismicity during injection at higher rates into the second well. The use of seismic moment to estimate the length of seismogenic slip planes in the local subsurface suggests that faults large enough to produce felt seismicity are unlikely to be present at or near the Decatur site.

Geochemistry & Geophysics↗

EGS Collab Experiment 1: Time-series geochemistry data of the long-term circulation test

This submission presents the weekly geochemistry data of the long-term flow test performed within EGS Collab Experiment 1 from early 2019 to early 2020. The fluids from each producing borehole/interval (PI, PB, PDT and PST) along with the injectate were sampled roughly weekly from April 2019 to January 2020 for geochemistry analysis. The geochemical measurement was part of a long-term microbial profiling project (see details in the PNAS paper linked below). Additional background and methodologies are available in the PNAS paper linked bellow.

15 GEOTHERMAL ENERGY↗

Preferential adsorption of noble gases in zeolitic tuff with variable saturation: A modeling study of counter-intuitive diffusive-adsorptive behavior

Noble gas transport through geologic media has important applications in the prediction and characterization of measured gas signatures related to underground nuclear explosions (UNEs). Retarding processes such as adsorption can cause significant species fractionation of radionuclide gases, which has implications for measured and predicted signatures used to distinguish radioxenon originating from civilian nuclear facilities or from UNEs. Accounting for the effects of variable water saturation in geologic media on tracer transport is one of the most challenging aspects of modeling gas transport because there is no unifying relationship for the associated tortuosity changes between different rock types, and reactive transport processes such as adsorption that are affected by the presence of water likewise behave differently between gas species. In this study, we perform numerical diffusive-adsorptive transport simulations to estimate gas transport parameters associated with bench-scale laboratory diffusion cell experiments measuring breakthrough in zeolitic and non-zeolitic rocks for a gaseous mixture of xenon, krypton, and SF 6 at varying degrees of water saturation (S w ). Counter-intuitive transport behavior was observed in the zeolitic rock experiments whereby breakthrough concentrations were significantly higher when the core was partially saturated (S w = 17 %) than under dry (S w = 0 %) conditions. Breakthrough of xenon was especially retarded in the dry core – likely due to comparatively high affinity of xenon for zeolitic adsorption sites – and estimated effective diffusion coefficients for all gases were approximately an order of magnitude lower than what is predicted by porosity-tortuosity models. We propose the counter-intuitive behavior observed is because water infiltration of zeolite nanopores reduces both the adsorptive capacity of the rock and the tortuosity of connected flow paths. We developed a two-site competitive kinetic Langmuir adsorption reaction for the porous media transport simulator in order to constrain transport parameters within zeolitic tuff, where differential adsorption to zeolite and non-zeolite pores was observed. We determined that liquid saturation-dependent diffusive-adsorptive transport is affected by subtle and at times competing processes that are specific to different gases, which have a significant overall influence on effective transport parameters.

58 GEOSCIENCES↗

The Role of Retardation, Attachment and Detachment Processes during Microbial Coal-Bed Methane Production after Organic Amendment

Microbially enhanced coal-bed methane could allow for a more sustainable method of harvesting methane from un-mineable coaldbeds. The model presented here is based on a previously validated batch model; however, this model system is based on upflow reactor columns compared to previous experiments and now includes flow, transport and reactions of amendment as well as intermediate products. The model implements filtration and retardation effects, biofilm decay, and attachment and detachment processes of microbial cells due to shear stress. The model provides additional insights into processes that cannot be easily observed in experiments. This study improves the understanding of complex and strongly interacting processes involved in microbially enhanced coal-bed methane production and provides a powerful tool able to model the entire process of enhancing methane production and transport during microbial stimulation.

54 ENVIRONMENTAL SCIENCES↗

Comparison of Geochemical Reactivity of Marcellus and Caney Shale Based on Effluent Analysis

ABSTRACT: In this comparative study, we analyzed the changes in elemental concentrations of hydraulic fracturing fluids after interaction with the Marcellus and Caney Shale formations. The focus was on assessing the inherent risks and environmental implications associated with flowback waters, including their impact on soil, ground, and drinking water quality, and human health safety. The chemical compositions of effluents were determined through Inductively Coupled Plasma Optical Emission Spectroscopy (ICP-OES) and Mass spectrometry (ICP-MS). The Marcellus Shale, showed concentrations of Cd, averaging 0.380 ppm far exceeding safe water thresholds. Significant levels of As, Se, B, and Pb were detected in both shales, raising concerns about soil and water contamination. Analyzing the effluents from representative samples of sections of the Marcellus (S2 and S7) and Caney (R1 and R2) indicates different geochemical responses over 4 weeks of experiments. This comparison underscores the chemical changes and environmental considerations linked to hydraulic fracturing across shale formations, suggesting the value of tailored monitoring and regulatory measures for each type. 1. INTRODUCTION The Caney and Marcellus shales, differing in geology and geochemistry, represent distinct unconventional reservoirs. The Caney shale, is more ductile, with higher produced water volumes, (Smith et al., 2022) contrasts with the brittle Marcellus shale known for lower brine production but significant data availability. This study aims to elucidate the possible environmental health and safety issues that may arise from the hydraulic fracturing processes. The interaction between fracturing fluids and clays presents a significant challenge. The primary base of these fluids is water, which, when introduced to clay, can induce swelling and constrict flow pathways. This phenomenon is attributed to water molecules infiltrating the layers of clay, particularly in 2:1-type clays, leading to an increased distance between layers. To counteract this, clay stabilizers are employed (Awejori et al., 2021). Despite their effectiveness, these stabilizers are considered temporary solutions. Upon completion of the fracturing process, a concomitant amount of contaminated waters (flow back), with varied content is collected at the surface. From these, we can infer the geochemical reactions and the environmental challenges associated with these waters. Flowback waters can be reinjected or used for other purposes such as irrigation. This requires adequate screening and treatment for safe use.

Dje, L. B.↗

Wound Field and Hybrid Synchronous Machines for EV Traction with Brushless Capacitive Rotor Field Excitation (Final Report)

This project focused on the development of wound field synchronous machines (WFSMs) and hybrid excitation synchronous machines (HESMs) with brushless capacitive power transfer for the field excitation. The target application for the machines developed is the main traction motor in electric vehicles. The magnetization in these types of machines is provided by a field winding on the rotor which is excited with DC current. The magnetization level in the machine can be varied by changing the magnitude of the field current. The variable magnetization or field is one of the key features of WFSMs. WFSMs are commonly used as generators however they have several attractive features for automotive traction applications. 1) No use of permanent magnet: Rare earth permanent magnets are primarily mined and processed in China. They have been subject to large price and supply variations and their export may be restricted during times of geopolitical tension. 2) Easy field weakening: Wound field synchronous machines have complete control of their field excitation. With proper design, this type of machine can electromagnetically have an infinite constant power speed range. 3) High power factor: With proper choice of the field excitation, WFSMs may be operated with high or even unity power factor. This potentially allows for the inverter connected to the stator winding to be downsized. In comparison induction machines and interior permanent magnet synchronous machines must supply reactive power to the stator increasing the kVA rating and cost of the inverter. 4) Reduced iron losses at high speed: By reducing the field excitation, the iron loss in the stator can be reduced. This is in comparison to interior permanent magnet synchronous machines which must use stator current to buck or reduce the flux produced by the permanent magnets. Generally, WFSMs have their highest efficiency at high speed. 5) Torque output at high temperatures: The magnetization provided by the field winding only depends on the field current and not on the field winding temperature. This is in contrast to permanent magnet machines where the permanent magnet flux decreases as the magnet temperature increases. Historically a number of approaches have been developed to provide DC current to the rotating field winding in WFSM’s including brushes and slip rings, low frequency brushless exciters, and high frequency rotary transformers and rectifiers. In this project, a different approach was used: brushless capacitive power transfer. Brushless capacitive power transfer uses two sets of rotating capacitors or electrodes in which an AC electric field is established by a high frequency inverter. A displacement current can flow through the airgap in the rotating capacitors which is rectified on the rotor using a diode bridge. The potential advantage of capacitive power transfer is that there is no need for heavy iron to guide magnetic flux. The electric flux lines terminate on the charges on the rotary capacitor surfaces. This should also limit the electric field outside the rotary capacitor airgaps. The main challenge with capacitive power transfer is that the capacitance and surface area of the rotating capacitors is small. Because capacitive power transfer systems are essentially a dual of a magnetic system, an Ampere per Hertz relationship is characteristic versus a Volts per Hertz relationship. A very high frequency power inverter must be used to provide sufficient excitation to the field winding. The concept of using capacitive power transfer to excite a high performance WFSM was initially developed in a previous U.S.A. Dept. of Energy project, DE-EE0006829. This project focused on increasing the power density of the WFSMs and reducing the cost and manufacturing complexity of the capacitive power transfer system. This project has demonstrated that WFSMs with brushless capacitive field power transfer can provide a high-power density and low-cost automotive powertrain technology.

33 ADVANCED PROPULSION SYSTEMS↗

UO 2 + 5 vol% ZrB 2 nano composite nuclear fuels with full boron retention and enhanced oxidation resistance

The boron isotope ( 10 B) can be used as a neutron absorber in UO 2 to control the reactivity of nuclear fuel pellets, however, the boron source can react with oxygen source in UO 2 to form B 2 O 3 that vaporize readily at temperatures above 1200 °C. Unfortunately, the sintering of UO 2 fuel requires hours holding at high temperature (>1700 °C), resulting in an inevitable B loss during sintering and unpredictable B concentration in final product. It is challenging to incorporate boron through a conventional sintering method. In this work, we demonstrated that spark plasma sintering (SPS), a field assisted sintering technology, can effectively densify UO 2 + 5 vol% ZrB 2 composite fuel pellets by rapid consolidation at 1600 °C for a short duration of 5 min under an applied pressure of 40 MPa. Further, thermogravimetric analysis (TGA) measurements confirm that ZrB 2 is fully retained inside the composite fuel pellets. Inside the composite fuel pellets, nano sized ZrB 2 particles are uniformly distributed along the grain boundaries of the UO 2 matrix. The ZrB 2 particle transforms to a glassy B 2 O 3 phase covering the sample surface and grain boundaries of UO 2 matrix after a simple post-sintering annealing at 1000 °C in flowing Argon gas for 4 h. The formed glassy B 2 O 3 slows down the diffusion of oxygen ions and postpones the onset temperature for oxidation of UO 2 from 400 °C to 550 °C. This study demonstrates the capability of SPS, an advanced fuel manufacturing technique, to achieve a full retention of ZrB 2 in UO 2 oxide fuel and increase oxidation resistance through a simple post-sintering annealing. The reported work holds great engineering potential for development of advanced oxide fuel for nuclear application.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Oxidative coupling of methane (OCM) conversion into C 2 products through a CO 2 /O 2 co-transport membrane reactor

Oxidative coupling of methane (OCM), which transforms CH 4 into C 2 products (C 2 H 6 and C 2 H 4 ) with molecular O 2 as the oxidant, is one of the most studied direct methane conversions (DMCs). However, a major technical hurdle to the OCM process is to achieve high CH 4 conversion at high C 2 (C 2 H 6 and C 2 H 4 ) selectivity. One rudimentary cause for this “tradeoff” behavior is the high chemical reactivity of the products (C 2 H 6 or C 2 H 4 ), which can be re-oxidized by O 2 . To overcome this thermodynamic challenge, minimizing the oxidizing power of the oxidant and lowering the local oxygen partial pressure are keys. Here, In this work, we demonstrate a new membrane reactor that capture CO 2 /O 2 from a flue gas and uses it for OCM conversion. The results show that the co-captured CO 2 /O 2 mixture converts CH 4 into C 2 H 6 in the presence of a 2%Mn–5%Na 2 WO 4 /SiO 2 catalyst, followed by thermal cracking of C 2 H 6 into C 2 H 4 and H 2 . The presence of CO 2 decreases the local partial pressure of O 2 , thus reducing the propensity of C 2 -products re-oxidation and leading to a higher C 2 selectivity. We show that a small button-type membrane reactor can achieve 12% C 2 yield with ~57% C 2 -selectivity using a diluted CH 4 -Ar mixture as the feedstock. We expect higher C 2 yield with tubular plug-flow membrane reactors in the future. We also highlight the unique advantage of the membrane reactor in intensifying CO 2 capture from both flue gas and OCM purification process into one single step.

42 ENGINEERING↗

Geomechanical in situ testing of fault reactivation in argillite repositories

Abstract. Pressurization of natural faults as a result of repository-induced effects can lead to their reactivation and permeability generation in case such features are present near disposal tunnels. Potential driving forces for such pressurization are the temperature increase caused by heat-producing high-level radioactive waste and the generation of hydrogen and other gases due to corrosion of engineered materials. We are primarily concerned about pressurization and fault activation in host rocks and faults that have very low natural permeability for pore pressure increases to dissipate, such as the argillite rocks currently investigated in Switzerland, France and other countries. This presentation discusses a series of in situ experiments of fault activation by fluid injection conducted in the argillite rock (Opalinus clay) at the Mont Terri underground research laboratory in Switzerland. A multi-model monitoring test bed was installed at Mont Terri that includes distributed fiber optics for strain, temperature, and acoustics; local fault pore pressure and three-dimensional displacement sensors; active seismic imaging; and passive seismic monitoring. The fault experiments (and their subsequent analysis via hydromechanical modeling) provide a new fundamental understanding of the coupling between pore pressure, fault deformation and permeability generation as a function of time and allow exploring how seismic and aseismic events may impact the integrity of faulted argillite host rock. Our experimental observations demonstrate that significant flow (and transport) can occur along the initially impermeable argillite fault when rupture is activated. However, the fault permeability decreases to almost its pre-activation value when fluid injection ceases and fluid pressure drops. Dilatant slip on the fault plane alone does not explain the observed increase in fault permeability; pressure-induced fault opening also plays a role, favored by the softness of the shale along with the fact that the structure of the fault zone prevents fluids from diffusing into the adjacent damage zone. Rupture initiation and permeability generation is initially aseismic, associated only with an increase in noise level and emerging tremors. Micro-earthquakes are initiated later in the experiments and typically occur away from the fluid-pressurized area. Hydromechanical models show that stress transferred from the initial aseismic deformation can build up to stress criticality and later induce seismic rupture. After presenting the experimental results, we will close the presentation with an outlook to future experimental campaigns using the fault test bed at Mont Terri. We are currently planning a controlled thermal stimulation of the fault, where instead of fluid injection as a trigger mechanism we will heat up the nearby rock volume and measure potential effects on fault stability. Lessons learned from our past and future experiments will help inform the safety assessment of geologic disposal in argillite host rock.

Birkholzer, Jens T.↗

Reactive Carbide‐Based Synthesis and Microstructure of NASICON Sodium Metal All Solid‐State Electrolyte

Reactive carbide precursor-based synthesis of NASICON-type NZSP (Na 1+x Zr 2 Si x P 3-x O 12 ) solid-state electrolyte (SSE) is demonstrated, in contrast to the established oxide-based approach. Exothermic decomposition of ZrC and SiC in air homogenizes microstructure, yielding 98% compact density after conventional sintering at 1200 °C. Quantitative stereology demonstrates that significant microstructural differences are present. Compacts of carbide-derived Carb-NZSP are 98% dense with a secondary zirconium oxide (ZrO 2 ) volume fraction of 0.2% ± 0.3%, versus 93% dense and 3% ± 1% for oxide-derived baseline. For Carb-NZSP, the secondary glassy phosphate phase is agglomerated, while for baseline, it is dispersed and percolated. Electrochemical testing combined with post-mortem analysis demonstrates how microstructural control of secondary phases is critical for dendrite suppression: Carb-NZSP critical current density (CCD) is 3.1 ± 0.8 mA cm −2 at 0.1 mAh cm −2 , versus 1.0 ± 0.7 mA cm −2 at 0.1 mAh cm −2 . Cryogenic focused ion beam (cryo-FIB) analysis demonstrates that in both materials, the porous 2D sheet-like sodium metal dendrites propagate around and subsume NZSP grains, likely following a path enriched with glassy phase and with porosity. Dendrites also flow around isolated zirconia particles. Phase field simulation reveals deflection of dendrites by mechanically tough zirconia, while brittle glassy phase accelerates dendrite growth, especially when finely distributed.

36 MATERIALS SCIENCE↗

Review of Power System Support Functions for Inverter-Based Distributed Energy Resources - Standards, Control Algorithms, and Trends

Penetration of renewable energy in power systems has been increasing in the past decades in response to increased global electricity demand and concerns for the environment. Distributed energy resources (DERs) based on renewables have experienced rapid growth thanks to the incentive programs and broad-based participation. With the growing prevalence of DERs, the risk of grid instability and vulnerability increases due to the intermittent nature of renewable energy. At the same time, the voltage and frequency deviation problems emerge more often when the reverse power flow occurs under supply-demand imbalance in distributed power systems. Standards and grid codes have been issued for DER inverters to interconnect with the distribution grid. The updated standard and grid codes expect DERs to provide a variety of power system support functions in order to incorporate higher DER penetration and to maximize DER value to the grid. This paper provides an overview of the power system support functions from renewable DER inverters, which are categorized as: voltage regulation by active/reactive power control, frequency regulation by active power control, voltage ride-through, and frequency ride-through. The benefits and drawbacks of each algorithm are presented and compared with its predecessor, manifesting the logic in the evolution of the algorithms.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Mechanically Accelerated Depolymerization of Entangled Linear Polymer Melts

Mechanical forces can enhance the chemical depolymerization of synthetic polymers when shear flow accelerates chain scission. To quantify the extent of mechanically-accelerated scission, the effect of simple shear flow (duration and strength) with low Weissenberg and Deborah numbers was investigated by considering the impact of applied work in both simple shear and shear dominated mixed flows. Hydrogenated polyisoprene was chosen as a model linear, entangled system. The conditions (strain amplitude, frequency, and shearing time) necessary to increase chain scission were assessed in the rubbery melt. Shear flow accelerated chain scission at higher temperatures, suggesting an activated process. Isothermal scission versus work curves were superposed by applying shift factors a T,S , whose Arrhenius-like temperature dependence gave an apparent activation energy for chain scission of ~ 110 kJ/mol, which is likely a combination of the activation energy of viscosity and bond energy. This work provides a base for quantifying the impact of shear on depolymerization of polymer melts and highlight the connection between viscous dissipation and scission chemistry.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Relating the depth of the water table to the depth of weathering

Weathering of bedrock creates and occludes permeability, affecting subsurface water flow. Often, weathering intensifies above the water table. On the other hand, weathering can also commence below the water table. To explore relationships between weathering and the water table, a simplified weathering model for an eroding hillslope was formulated that takes into account both saturated and unsaturated subsurface water flow (but does not fully account for changes in dissolved gas chemistry). The phreatic line was calculated using solutions to mathematical treatments for both zones. In the model, the infiltration rate at the hill surface sets both the original and the eventual steady-state position of the water table with respect to the weathering reaction front. Depending upon parameters, the weathering front can locate either above or below the water table at steady state. Erosion also affects the water table position by changing porosity and permeability even when other hydrological conditions (e.g. hydraulic conductivity of parent material, infiltration rate at the surface) do not change. The total porosity in a hill (water storage capacity) was found to increase with infiltration rate (all else held constant). Furthermore, this effect was diminished by increasing the erosion rate. We also show examples of how the infiltration rate affects the position of the water table and how infiltration rate affects weathering advance.

58 GEOSCIENCES↗

Impact of nuclear data covariance libraries on uncertainty quantification of sodium cooled fast reactor simulation

Uncertainty quantification in core modelling to assess reactor performance is a regular practice to identify design safety limits. Best-estimate plus uncertainty (BEPU) methods are commonly used to propagate uncertainties to core responses from input parameters, characterized by variance-covariance libraries. This study evaluates the impact of different covariance libraries on the uncertainty associated with core safety parameters and their propagation through reactor transient modelling for a metallic core (ABR-1000). The COMAC-1.0, COMMARA-2.0, and COMMARA-2.1 covariance libraries are used to propagate nuclear data uncertainties to neutron feedback coefficients using the Generalized Perturbation Theory, then through transient safety analyses to establish the confidence interval for safety performance of the reactor using stochastic sampling. The considered transients include Unprotected Transient Over Power and Unprotected Loss Of Flow. Most influential nuclide and reaction pairs in different variance-covariance libraries are identified and their impact on the uncertainty of safety parameters are evaluated. Significant differences were observed between uncertainties from COMAC compared to COMMARA libraries. Uncertainties from COMMARA libraries were consistently larger than COMAC for all cases. For structure, coolant density feedback coefficients and control rod worth, total nuclear data uncertainties were approximately 50% larger for COMMARA than COMAC. Differences were also observed in identification of top contributors of uncertainty in terms of reaction-cross sections. For instance, largest contribution to total uncertainty of radial expansion coefficient from COMMARA-2.0 and 2.1 originates from Na{sup 23} P1 Elastic cross section (0.81%) but from Pu{sup 239} Fission (0.53%) with COM C-1.0. Such differences are evaluated in this research for all quantities of interest. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗