Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Flow Conditions”

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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Interactions between molecular-scale processes and hyporheic exchange for understanding Fe-S-C cycling in riparian wetlands (Final Report)

Wetlands represent some of the most productive ecosystems on the planet and critically influence global environmental health. Specifically, wetlands promote water quality by transforming nutrients and organic compounds and sequestering metals and contaminants. Riparian wetland hyporheic zones, where toxic surface water and anoxic groundwater mix, exhibit dynamic conditions that drive steep redox gradients and promote hotspots of diverse and fluctuating microbial activity. Changes in climate, water quality, and water quantity can disturb hydrologic flow and biogeochemical processing in these environments. Understanding how sulfate loading, impacted by hydrologic flux and anthropogenic inputs, influences iron and carbon cycling in wetlands will be crucial for predicting water quality issues driven by iron mineral precipitation and sorption, such as the release of heavy metals and other toxic elements. We used a fully integrated multi-scale and multi-method approach to develop a mechanistic understanding of how hydrologic flow influences coupled iron and sulfur cycles in riparian wetlands. This entailed hydrological, geochemical, and microbial observations at two locations: an anthropogenic sulfate-impacted riparian wetland in northern Minnesota, and a low-sulfate Fe-rich riparian wetland in Tims Branch at the Savannah River Site (SRS). These sites were characterized by hydrologically dynamic conditions where the stream and wetland systems oscillated between gaining (upward flow) and losing (downward flow) conditions that recharged the system occasionally with oxidants that fueled a variety of biogeochemical reactions. Aqueous geochemical measurements of surface water, groundwater, and porewater samples were made alongside solid-phase geochemical analyses of sediment gravity cores. Bulk X-ray absorption spectroscopy at the Advanced Photon Source (APS; Argonne) interrogated the speciation and distribution of Fe and S mineral phases of the sediments. Interestingly, it was discovered that, despite strongly reducing conditions, Fe(III) compounds and a variety of intermediate valence S compounds were stable in the subsurface. Indeed, compounds like thiosulfate, S(0), and intermediate valence organosulfur compounds were more prevalent than FeS and pyrite. The composition of the sediments did change with changes in hydrologic flow, showing their reactivity in changing redox conditions. The abundance of these intermediate S compounds was likely formed as a result of anaerobic oxidation by aqueous and solid-phase Fe(III) compounds, fueling a cryptic S cycle that is driving the breakdown of organic matter in the hyporheic zone. Microbiome surveys showed a core community that seemed stable across the landscape, but changed with increasing depth into the sediment. The community composition did not seem to change dramatically with changes in season or hydrologic flow, except for some organisms that have the potential to contribute to S cycling. More work is needed to confirm their functional activity. These fine process-scale analyses were placed within a dynamic field context using physical flow parameters from surface water and groundwater level measurements. These data helped shed light on sulfur-driven biogeochemical processes in hydrologically dynamic riparian wetlands, addressing a gap in our understanding about the impacts of pollution and other anthropogenic changes on ecologically sensitive environments.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Directing the far-from-equilibrium assembly of nanoparticles in confined liquid crystals by hydrodynamic fields

The assembly of nematic colloids relies on long-range elastic interactions that can be manipulated through external stimuli. Confinement and the presence of a hydrodynamic field alter the defect structures and the energetic interactions between the particles. Here, the assembly landscape of nanoparticles embedded in a nematic liquid crystal confined in a nanochannel under a pressure-driven flow is determined. The dynamics of the liquid crystal tensor alignment field is determined through a Poisson-Bracket framework, namely the Stark–Lubensky equations, coupled with the zero-Reynolds momentum equations and the liquid crystal Landau-de Gennes free energy functional. A second order semi-implicit time integration and a three-dimensional Galerkin finite element method are used to resolve flow and nematic fields under several conditions. In general, the zero Reynolds flow displaces the defects around the particles in the upstream direction and renders the surface anchoring ineffective when the flow strength dominates over the nematic elasticity. More importantly, the potential of mean force for particle assembly is non-monotonic independent of surface anchoring. Our results show that the confinement length scale determines the repulsion/attraction transition between colloids, while the flow strength modifies the static defect structure surrounding the particles and determines the magnitude of the energetic barrier for successful assembly. In the attractive regime, the particles move at different rates through the nematic until one particle eventually catches up with the other. This process occurs against or along the direction of flow depending on the flow strength. Ultimately, these results provide a template for engineering and controlling the transport and assembly of nanoparticles under far-from equilibrium conditions in anisotropic media.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Reproducible benchmark for the SNAP 8 experimental reactor at operating conditions

This work presents fully reproducible multiphysics benchmark models of the Systems for Nuclear Auxiliary Power (SNAP) 8 Experimental Reactor at operating conditions with coolant flow. Wet experiment (with coolant, at power) validation benchmarks are presented using both deterministic (Serpent-Griffin) and Monte-Carlo (OpenMC-Cardinal) multiphysics frameworks coupled with thermal-hydraulic solvers in MOOSE. Reactivity coefficient measurements including fuel temperature, isothermal temperature, and power coefficients show good agreement with experiments, with discrepancies within experimental uncertainty. Reactivity worth experiments for coolant, samarium, and xenon poisoning are reproduced with differences under 200 pcm. Comparison between Serpent-Griffin and OpenMC-Cardinal frameworks reveal multiphysics coupling introduces positive reactivity effects (100-200 pcm) compared to uniform temperature and density fields at nominal operating conditions. Comparison between Serpent-Griffin and reference Serpent solution shows that power distributions maintain consistent radial and axial peaking behavior. All models, assumptions, thermophysical and thermomechanical properties, and material definitions are thoroughly documented with cited references; model inputs and model generating scripts are stored in the snapReactors GitHub repository.

SNAP↗
DOE OSTI2025

Computational Fluid Dynamics Analysis of the Molten Salt Tritium Transport Experiment Test Section

Tritium, a radionuclide produced through neutron capture by lithium and other elements (beryllium and fluoride) in molten salts, presents unique challenges to radionuclide release. This is true for both fusion energy breeder blankets and molten salt fission reactors. The fundamental understanding of tritium transport is crucial to the safe design and operation of these reactors. The Molten Salt Tritium Transport Experiment (MSTTE), currently under construction at Idaho National Laboratory, aims to investigate tritium transport phenomena using a forced-convection fluoride salt loop. This loop is designed to study various transport mechanisms, such as permeation through metals and gas-liquid interactions, and is intended to support future research on tritium extraction units. A critical aspect of the MSTTE loop design is ensuring a fully developed velocity profile before the fluid reaches the permeation test section where measurements are made. This study employs computational fluid dynamics to model the salt flow behavior within the MSTTE permeation test section. A realizable k-ε turbulent model with enhanced wall treatment is used to simulate the single-phase, vertical upward flow of molten salt FLiNaK under isothermal conditions. The simulation results indicated flow distortion and underdeveloped profiles at all planned flow rates within the test section due to the 85-deg sharp bend. To address this issue, a reduced diameter with a reducer and expander and a flow conditioner are investigated to achieve fully developed flow. The analysis showed that the flow conditioner successfully corrected the flow profile, achieving fully developed behavior at a flow rate of 50 liters per minute (LPM). In conclusion, this research enhances our understanding of flow dynamics in molten salt systems and contributes to optimizing tritium transport control technologies.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Thermal-hydraulic and Fuel Performance Scoping Studies of a Flowing Water Capsule in TREAT

The restart of the Transient Reactor Test (TREAT) facility has provided a much needed capability for integral transient testing of nuclear fuel. This testing is necessary to qualify new fuel concepts such as those developed under the Accident Tolerant Fuel (ATF) program, increasing the burnup limit of light water reactor (LWR) fuels, or a variety of other programs under the Advanced Fuels Campaign (AFC). The ATF campaign has been the main driver behind the development and implementation of a variety of capsules for TREAT experiments. The Separate-Effect Test Holder (SETH) was a dry capsule that enabled testing of ATF concepts under reactivity-initiated accident (RIA) heating conditions. Following SETH, the Static Environment Rodlet Transient Test Apparatus (SERTTA) was developed to enable RIA testing in a static water environment. In efforts to support the need for future Loss-of-Coolant Accident (LOCA) tests, the Transient Water Irradiation System for TREAT (TWIST) capsule has been developed that allows for water to drain from around the fuel rod and lower the pressure to simulate LOCA conditions. All these capsules that have been developed for LWR fuel testing all lack the capability for forced convective cooling which in some applications may limit their ability to test under prototypic conditions. A design and modeling effort has been started to modify the TWIST capsule by adding a flow tube and impeller that can create flowing coolant conditions for the fuel rod. RELAP5-3D and BISON models are being used to study the differences between RIA, LOCA, Anticipated Operational Occurrences (AOO), LWR power cycling, and other transient scenarios under flowing conditions capable in the flowing capsule and the current stagnant water capsules (SERTTA and TWIST). The scoping study will provide guidance on the needed capabilities for the flowing capsule and the limitations of the currently developed capsule for LWR testing in TREAT.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

CFD-based design optimization of ducted hydrokinetic turbines

Abstract Hydrokinetic turbines extract kinetic energy from moving water to generate renewable electricity, thus contributing to sustainable energy production and reducing reliance on fossil fuels. It has been hypothesized that a duct can accelerate and condition the fluid flow passing the turbine blades, improving the overall energy extraction efficiency. However, no substantial evidence has been provided so far for hydrokinetic turbines. To investigate this problem, we perform a CFD-based optimization study with a blade-resolved Reynolds-averaged Navier–Stokes (RANS) solver to explore the design of a ducted hydrokinetic turbine that maximizes the efficiency of energy extraction. A gradient-based optimization approach is utilized to effectively deal with the high-dimensional design space of the blade and duct geometry, with gradients being calculated through the adjoint method. The final design is re-evaluated through higher-fidelity unsteady RANS (URANS) simulations. Our optimized ducted turbine achieves an efficiency of about 54% over a range of operating conditions, higher than the typical 46% efficiency of unducted turbines.

13 HYDRO ENERGY↗
DOE OSTI2022

Data-driven wind turbine wake modeling via probabilistic machine learning

Wind farm design primarily depends on the variability of the wind turbine wake flows to the atmospheric wind conditions and the interaction between wakes. Physics-based models that capture the wake flow field with high-fidelity are computationally very expensive to perform layout optimization of wind farms, and, thus, data-driven reduced-order models can represent an efficient alternative for simulating wind farms. In this work, we use real-world light detection and ranging (LiDAR) measurements of wind-turbine wakes to construct predictive surrogate models using machine learning. Specifically, we first demonstrate the use of deep autoencoders to find a low-dimensional latent space that gives a computationally tractable approximation of the wake LiDAR measurements. Then, we learn the mapping between the parameter space and the (latent space) wake flow fields using a deep neural network. Additionally, we also demonstrate the use of a probabilistic machine learning technique, namely, Gaussian process modeling, to learn the parameter-space-latent-space mapping in addition to the epistemic and aleatoric uncertainty in the data. Finally, to cope with training large datasets, we demonstrate the use of variational Gaussian process models that provide a tractable alternative to the conventional Gaussian process models for large datasets. Furthermore, we introduce the use of active learning to adaptively build and improve a conventional Gaussian process model predictive capability. Overall, we find that our approach provides accurate approximations of the wind-turbine wake flow field that can be queried at an orders-of-magnitude cheaper cost than those generated with high-fidelity physics-based simulations.

Deep neural networks↗
DOE OSTI2026

Hybrid combustion modeling approach for turbulent jet ignition in natural-gas pre-chamber spark-ignition engines at high EGR

Here, this study presented a hybrid modeling approach for simulating turbulent jet ignition and combustion processes in a natural-gas pre-chamber spark-ignition engine operating under exhaust gas recirculation (EGR) diluted conditions. In-depth analyses of experimental data and simulation results from previous work [Chinnathambi et al., ICEF2021-67836; Kim et al., Fuel 409: 137815, 2026] revealed two key findings: (i) the magnitude of pressure difference between the pre-chamber and main chamber ($∆P_{PC-MC}$) was positively correlated with the combustion duration from the moment of $∆P_{PC-MC}=0$ to the point of 5% mass fraction burned, with larger $∆P_{PC-MC}$ associated with longer duration; and (ii) the turbulent combustion regime in the main chamber transitioned from the broken reaction zone to the corrugated flamelet regime, with the Karlovitz number exceeding 100 immediately after turbulent hot jets were ejected from nozzles, coinciding with observed local extinction events. To accurately simulate the entire combustion process, a hybrid approach was developed under Reynolds-Averaged Navier Stokes framework, combining the G-equation model for pre-chamber combustion with the multi-zone well-stirred reactor approach and a turbulence-chemistry interaction (TCI) submodel for main chamber combustion. The TCI submodel accounted for the attenuation of reaction rates due to turbulent strain and modeled local extinction by suppressing reaction rates under certain flow and flame conditions. When applied to three EGR rate conditions toward the dilution limit, the hybrid modeling approach accurately reproduced experimental data in terms of cylinder pressure, apparent heat release rate, and the observed positive correlation, including the delayed onset of main chamber combustion—a feature not captured by existing combustion models.

computational fluid dynamics simulation↗
DOE OSTI2023

OC6 project Phase III: validation of the aerodynamic loading on a wind turbine rotor undergoing large motion caused by a floating support structure

Abstract. This paper provides a summary of the work done within Phase III of the Offshore Code Comparison Collaboration, Continued, with Correlation and unCertainty (OC6) project, under the International Energy Agency Wind Technology Collaboration Programme Task 30. This phase focused on validating the aerodynamic loading on a wind turbine rotor undergoing large motion caused by a floating support structure. Numerical models of the Technical University of Denmark 10 MW reference wind turbine were validated using measurement data from a 1:75 scale test performed during the UNsteady Aerodynamics for FLOating Wind (UNAFLOW) project and a follow-on experimental campaign, both performed at the Politecnico di Milano wind tunnel. Validation of the models was performed by comparing the loads for steady (fixed platform) and unsteady (harmonic motion of the platform) wind conditions. For the unsteady wind conditions, the platform was forced to oscillate in the surge and pitch directions under several frequencies and amplitudes. These oscillations result in a wind variation that impacts the rotor loads (e.g., thrust and torque). For the conditions studied in these tests, the system aerodynamic response was almost steady. Only a small hysteresis in airfoil performance undergoing angle of attack variations in attached flow was observed. During the experiments, the rotor speed and blade pitch angle were held constant. However, in real wind turbine operating conditions, the surge and pitch variations would result in rotor speed variations and/or blade pitch actuations, depending on the wind turbine controller region that the system is operating. Additional simulations with these control parameters were conducted to verify the fidelity of different models. Participant results showed, in general, a good agreement with the experimental measurements and the need to account for dynamic inflow when there are changes in the flow conditions due to the rotor speed variations or blade pitch actuations in response to surge and pitch motion. Numerical models not accounting for dynamic inflow effects predicted rotor loads that were 9 % lower in amplitude during rotor speed variations and 18 % higher in amplitude during blade pitch actuations.

17 WIND ENERGY↗
DOE OSTI2022

Stochastic boundary condition effects on supersonic leading edge blowing

The present manuscript studies the robustness of leading edge blowing for active control of the upstream bow shock and overall aerothermal characteristics of a supersonic airfoil subject to variability in the injection conditions. The injection boundary conditions were considered stochastic and their randomness was propagated to the flow topology and aerothermal quantities of interest by means of non-intrusive Polynomial Chaos Expansions. The evaluations needed to build these expansions were carried out via steady two-dimensional Reynolds-Averaged Navier Stokes simulations. Two different injection port sizes were considered, with the smallest promoting a Coanda effect that gives rise to non-symmetric flow topologies. The uncertain boundary conditions lead to significant variations in the flow topology, particularly at the injection exit and the bow shock position, with variations of up to 75% of the nominal value. Here this study identifies the key operating conditions to ensure the optimal application of leading-edge injection for active control of supersonic airfoils.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Visco-inertial gas flow through wellbore cement fractures

Understanding the nature of fluid flow through fractured wellbore cement is fundamental for evaluating the leakage potential and risk assessments of leaky wellbores. In this study, the conditions that require considering visco-inertial flow for describing the gas flow through wellbore cement fractures were investigated. Nitrogen gas flow tests were conducted on fractured cement samples under varying pressure conditions and flow rates, covering both viscous and visco-inertial flow regimes. The data substantially deviated from Darcy's law at higher flowrates and were well-fit to Forchheimer's equation for visco-inertial flow. The inertial coefficient and critical Reynolds number were expressed as a function of the hydraulic aperture. Finally, the empirical function obtained from the experiments was used as an input to numerical simulations which showed the significant role of visco-inertial flow in wellhead pressure build-up and leakage rates, and demonstrated the importance of visco-inertial flow when modeling gas flow through wellbore cement fractures.

03 NATURAL GAS↗
DOE OSTI2021

Quantification of energy and cost reduction from decreasing dissolved oxygen levels in full‐scale water resource recovery facilities

Abstract Aeration systems often lack the efficiency to maintain a desired residual dissolved oxygen (DO) concentration in the tank in part because little consideration is given to the dynamic daily and seasonal loading conditions. Although advanced aeration controllers exist, the majority of plants have DO set points typically based on common practice and literature values rather than site‐specific conditions, which can result in DO set points higher than those necessary to meet treatment objectives. DO set point reduction strategies have primarily been proposed through either static or dynamic simulations. In this study, the substantial improvements associated with DO set point reduction are demonstrated at full scale. A yearlong characterization of full‐scale aeration dynamics captured the effect of diurnal and seasonal fluctuations on oxygen transfer and energy demand and so facilitated the estimation of the potential savings of DO reduction strategies. Full‐scale validation provided direct evidence of DO reduction strategies inducing an overall enhancement of oxygen transfer efficiency along the different bioreactors, while confirming that energy savings as high as 20% were feasible. This study quantifies the influence of oxygen transfer efficiency on operating choices and site‐specific conditions (control strategy, loading conditions, and influent flow variability). Practitioner points We quantified the energy reduction and cost savings associated with a DO reduction in an aeration tank. For each 0.2 mg/L of DO decreased, the average power demand reduction per unit water treated exceeded 17%. Field measurements of dynamic alpha values eliminate the uncertainty in estimating aeration energy and cost savings from DO variations.

Pasini, Federico↗
DOE OSTI2021

Analysis and prediction of intersection traffic violations using automated enforcement system data

We report that the automated enforcement system (AES) is an effective way of supplementing traditional traffic enforcement, and the traffic violation data from AES can also be effectively used for safety research. In this study, traffic violation data were used to analyze the influencing factors associated with traffic violations and to predict the probability of violations at intersections. The potential factors influencing violations include 24 independent factors related to time, space, traffic and weather. Results from a logistic model showed that the midday period, weekends, residential districts, collector roads, congested traffic conditions, high traffic flow, lower wind speed and low temperature would increase the probability of traffic violations. The probability of violations was predicted by the random forest algorithm, which was proven to be the best traffic violation prediction model among logistic regression, Gaussian naive Bayes, and support vector machine. Moreover, the proximity weighted synthetic oversampling technique (ProWSyn) method was applied to reduce the impact of the imbalance ratio (IR) and improve the model’s prediction performance. The receiver operating characteristics (ROC) curves and Precision-Recall (PR) curves illustrated that the random forest algorithm using oversampling data had the best classifier prediction performance than undersampling data. The area under curve (AUC) and out-of-bag (OOB) error with IR = 1 reached 0.914 and 0.0787, which showed the better performance of the random forest algorithm using ProWSyn in dealing with imbalanced traffic violation data.

42 ENGINEERING↗
DOE OSTI2020

Molecular dynamics and continuum analyses of the electrokinetic zeta potential in nanostructured slit channels

This work presents a theoretical and numerical study of electrokinetic flow and the zeta potential for the case of slit channels with nanoscale surface features of dimensions comparable to the Debye length, by employing molecular dynamics simulations and continuum-level analyses. Additionally, a simple analytical model for considering the average effect of such surface features producing nanoscale roughness is proposed by employing matched asymptotic solutions for the charge density and fluid flow field, and matching conditions that satisfy electroneutrality and the Onsager reciprocal relation between the electroosmotic flow rate and streaming current. The proposed analytical model quantitatively accounts for results from molecular dynamics simulations that consider the presence of ion solvation shells and surface hydration layers. Our analysis indicates that a simultaneous knowledge of the electroosmotic and pressure-driven flow rate or streaming current can be instrumental to unambiguously determine the zeta potential in the presence of nanoscale surface features.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

TRITON: A Multi-GPU open source 2D hydrodynamic flood model

A new open source multi-GPU 2D flood model called TRITON is presented in this work. The model solves the 2D shallow water equations with source terms using a time-explicit first order upwind scheme based on an Augmented Roe's solver that incorporates a careful estimation of bed strengths and a local implicit formulation of friction terms. Here, the scheme is demonstrated to be first order accurate, robust and able to solve for flows under various conditions. TRITON is implemented such that the model effectively utilizes heterogeneous architectures, from single to multiple CPUs and GPUs. Different test cases are shown to illustrate the capabilities and performance of the model, showing promising runtimes for large spatial and temporal scales when leveraging the computer power of GPUs. Under this hardware configuration, communication and input/output subroutines may impact the scalability. The code is developed under an open source license and can be freely downloaded in https://code.ornl.gov/hydro/triton.

2D flood model↗
DOE OSTI2022

Manganese(II) complexes of 1,1'-bis[(pyridin-2-yl)methyl)]-2,2'-bipiperidine (PYBP): Synthesis, structure, catalytic properties in alkene epoxidation with hydrogen peroxide, and related mechanistic studies

In this study, several manganese(II) complexes with the stereoisomers of ligand PYBP (1,1'-bis[(pyridin-2-yl)methyl]-2,2'-bipiperidine) and different anions were prepared and characterized by X-ray diffractometry. Complex [Mn II (rac-PYBP)] 2+ (1) was found to be an efficient catalyst of alkene epoxidation by hydrogen peroxide in the presence of acetic acid in acetonitrile at room temperature. Cyclooctene was converted to its epoxide with up to 91 % yield, 99.6 % selectivity, and the turnover number of 180 within 5 min. Fast epoxidations of cyclohexene, 1-decene, styrene, and cis-stilbene were also achieved. Isomeric complex [Mn II (meso-PYBP)] 2+ (2) was catalytically inactive under the same experimental conditions. Stopped-flow spectrophotometry and freeze-quenched EPR spectra show that complex 2 is not oxidized by H 2 O 2 in the presence of acetic acid (AcOH) but instead undergoes partial ligand protonation and liberation of the Mn 2+ cations due to the relatively poor chelating ability of ligand meso-PYBP. The rac-PYBP isomer acts as a better ligand and retains the coordinated Mn center when complex 1 is treated with the H 2 O 2 /AcOH mixture in acetonitrile solution yielding a mixture of intensely colored intermediates likely involving Mn III , Mn IV , and Mn V complexes. Magnetic susceptibility measurements, UV–vis and EPR spectra suggest that dinuclear complexes [Mn III 2 (μ-O)(μ-OAc)(rac-PYBP) 2 ] 3+ and [Mn III Mn IV (μ-O) 2 (rac-PYBP) 2 ] 3+ gradually accumulate in the reaction mixture as inactivated states of the catalyst. Complex 1 also causes fast decomposition of hydrogen peroxide into O 2 gas and H 2 O, which competes with the epoxidation of alkenes and requires gradual addition of H 2 O 2 for its efficient use. The catalytic activity of complex 1 is strongly influenced by its counterions and decreases in order ClO 4 - ≈ SbF 6 - > NO 3 – > Cl - indicating that labile ligands in the coordination sphere of Mn are required for the activation of H 2 O 2 .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Shear strength and permeability of the cement-casing interface

Here, the shear strength and hydraulic permeability of the interface between well cement and casing was investigated using a triaxial direct shear apparatus. For the first time, these experiments provide measurements under controlled stress conditions with fluid flow measurements along the interface. The low cohesion (1.1 ± 1.1 MPa) and the high friction angle (43.4 ± 2.0°) indicates that the shear strength of the interface is provided by friction. This implies that the state of stress of the cement is critical to well integrity. The hydraulic aperture of the undamaged cement-steel samples was 6.8 ± 1.0 microns. Shear damage to the interface caused a decrease (-20 %) in hydraulic aperture for samples aged up to 1 month, and an increase (+300 %) for samples cured for two years. We performed numerical simulations to estimate the leakage potential from a carbon storage operation. This model predicts negligible leakage amounts (47 tonnes) in a shear-damaged well for the modeled injection of ~1.26 million tonnes of CO 2 . Thus, our measurements indicate that the cement-casing interface is not a significant leakage pathway in its intact or damaged state, and that shear-driven failure scenarios for this interface are not a significant risk to CO 2 storage security.

54 ENVIRONMENTAL SCIENCES↗