SEARCH · Engineering Papers
Results for “Effective stress”
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.
A study on texture stability and the biaxial creep behavior of as-hydrided CWSR Zircaloy-4 cladding at the effective stresses from 55 MPa to 65 MPa and temperatures from 300°C to 400°C
The creep rupture of high burnup used nuclear fuel (> 45 GWD/MTU) cladding is regarded as one of the failure mechanisms during long-term dry storage. A high amount of zirconium hydride in the cladding matrix would degrade the mechanical properties of the cladding especially leading to delayed hydride cracking. To better understand the influence of zirconium hydride on the biaxial thermal creep performance and the crystalline texture stability of Zircaloy-4 cladding, the pressurized tube technique is applied to test the durability of the as-hydrided material. Tests were performed on as-received Zircaloy-4 tubular specimens as well as-hydrided ones with targeted 300 wt parts per million (wppm) or 750 wppm hydrogen. The biaxial creep experiments were conducted at temperatures from 300 degrees C to 400 degrees C and at equivalent stresses at mid-wall from 55 MPa to 65 MPa. The hydridation process prior to creep tests induces the formation of FCC delta-hydride platelets along the circumferential direction of the tube. This alignment and phase structure of hydride show no significant change after biaxial creep tests. The creep strain-rate negatively depends on the hydrogen content. The synchrotron wide-angle X-ray diffraction (WAXD) technique and electron backscatter diffraction (EBSD) analysis were applied for the study of the crystallographic orientation relationship. zirconium-to-hydride grain orientation follows Shoji-Nishiyama crystallographic relationship. This relationship is stable before and after creep deformation. These results of creep performance and texture stability of Zircaloy-4 claddings can help support the design basis of interim and long-term dry storage facilities.
Initial Acoustoelastic Measurements in Olivine: Investigating the Effect of Stress on P- and S-Wave
Explore the source record for details and available documents.
Photo-oxidation of semicrystalline polymers: Effect of stress triaxiality on ductility
The effect of stress triaxiality on the strain-to-fracture of as-received and photo-oxidized polyamide-6 (PA-6) was investigated using mechanical testing, synchrotron X-ray tomography, and finite element analyses. Mechanical tests were conducted on cylindrical and round notched specimens, where different notch radii were used to vary the stress triaxiality. The specimens were aged by exposure to ultra-violet (UV) radiation at 60∘, causing photo-oxidation. As-received and so-aged specimens were loaded to failure (complete loss of load carrying capacity). For both unaged and aged specimens, a higher triaxiality led to a lower strain-to-fracture. To elucidate the micromechanical damage that mediates fracture in both conditions, specimens with an intermediate notch sharpness were loaded to the peak load, unloaded, and scanned ex situ using synchrotron X-ray tomography. Damage in the unaged bar was found to occur by cavitation and was concentrated at the center of the specimen, where the triaxiality is highest. In the UV-aged bar, a network of inter-connected chemical cracks were found on the notch surface, where the triaxiality is lowest. Finite element analyses were deployed to approximate the local triaxiality at damaged regions in the unaged and UV-aged specimens using a constitutive relation for semicrystalline polymers. From these analyses, the relationship between local triaxiality and strain-to-fracture was quantified for both unaged and photo-oxidized PA-6. Both unaged and photo-oxidized PA-6 showed similar decreases in ductility with triaxiality, hinting at common ductile fracture processes.
Effect of Stress on Irradiation Responses of Highly Oriented Pyrolytic Graphite
The effect of stress on irradiation responses of highly oriented pyrolytic graphite (HOPG) was studied by combing molecular dynamics (MD) simulation, proton irradiation, and Raman characterization. MD simulations of carbon knock-on at energies < 60 eV were used to obtain average threshold displacement energies (E¯d) as a function of strain ranging from 0 to 10%. Simulations at a higher irradiation energy of 2–5 keV were used to study the effect of strain on damage cascade evolution. With increasing tensile strain, E¯d was reduced from 35 eV at 0% strain to 31 eV at 10% strain. The strain-reduced E¯d led to a higher damage peak and more surviving defects (up to 1 ps). Furthermore, high strains induced local cleavage around the cavities, as one additional mechanism of damage enhancement. Experimentally, HOPG film was folded, and the folded region with the maximum tensile stress was irradiated by a 2 MeV proton beam. Raman characterization showed significantly enhanced D to G modes in comparison to the stress-free irradiation. Based on the strain dependence of E¯d and the Kinchin–Pease model, a formula for displacement estimation under different tensile strains is proposed. The stress effects need to be considered in graphite applications in a reactor’s harsh environment where both neutron damage and stress are present.
Insights into the soft brittle-to-ductile transition from discrete dislocation dynamics
The Brittle-to-ductile transition (BDT) in body centered cubic metals exhibits a soft transition wherein the fracture toughness gradually rises to before the onset of ductility. The resultant brittle-to-ductile transition temperature can be described with an Arrhenius relationship whose activation energy is related to plasticity in the material. To provide further insight into the nature of the BDTT, in this work we utilized a discrete dislo- cation dynamics model with a crack to simulate the BDT and how it depends on the thermally activated nature of plasticity. The interrelationship between the BDT activation energy and the dislocation mobility parameters were determined via the calculation of first order sensitivity coefficients. This analysis allows us to demonstrate that the activation energy for the BDT is directly related to the activation energy for plasticity through an effective stress that defines this relationship. This effective stress physically is the average stress on the dislocations that move out of the crack. Lastly, we are able to show that this effective stress is dictated by the low temperature fracture toughness or cleave energy of the material and the source position, the latter of which can be affected by processing. Collectively, these results provide new insight into what controls the thermal activation of the BDT and what are the important parameters to control it.
Integrating Experiments and Well Logs to Predict Caney Shale Static Mechanical Properties During Production with Supervised Machine Learning
Caney shale is one of the emerging oil reservoirs in Oklahoma. Understanding the impact of effective stress on its mechanical properties is critical for predicting hydraulic fracture geometry and overall hydrocarbon production. The objective of our study is to evaluate the impact of effective stress on the dynamic Young’s modulus using ultrasonic velocity measurements for Caney shale samples. A triaxial cell was utilized to measure ultrasonic (P-wave and S-wave) velocities for ten downhole Caney shale samples under various effective stresses to indirectly assess the impact of pore pressure change. The dynamic Young’s moduli estimated from these measurements were integrated with available conventional well logs (excluding sonic logs) and triaxial test results from Benge et al. (2021) to predict the static Young’s modulus using Random Forest (RF) and Extreme Gradient Boosting (XGBoost) models. The results showed that the estimated dynamic Young’s moduli from ultrasonic measurements were higher than the corresponding static Young’s modulus of cores from the same vertical well at similar depths. With increasing effective stress, the dynamic Young’s modulus increased for all samples. The estimated dynamic-to-static correction factor tended to be higher in zones of high neutron porosity (PHIN) and low density compared to other zones. Finally, SHapley Additive exPlanations (SHAP) for RF and XGBoost models identified depth, gamma ray (GR), and PHIN as key features for predicting the static Young’s modulus. This study enhances our understanding of the dynamic and static Young’s moduli for the Caney shale interval, as a function of effective stress and conventional well logs. The findings from this study can improve predictions of production throughout the well's lifespan by offering insights into the mechanical property degradation resulting from pore pressure depletion.
Anisotropic elastoplastic response of double-porosity media
Here we present a continuum framework to simulate fluid flow through anisotropic elastoplastic media with double porosity. Two effective stress measures σ' and δ" emerge from the thermodynamic formulation, which are energy-conjugate to the elastic and plastic components of strain, respectively. Both effective stress measures can be expressed as a combination of the total Cauchy stress σ and the average pore pressure $\bar{p}$ in the two pore scales. In the effective stress for elasticity, $\bar{p}$ is scaled with a rank-2 Biot tensor, whereas the effective stress for plasticity follows the Terzaghi form in which $\bar{p}$ is scaled by the Kronecker delta. The Biot tensor and storage coefficients are derived as functions of elasticity parameters and porosities. A mixed finite element formulation is introduced to discretize the domain and solve initial boundary-value problems. A stabilization scheme is employed on equal-order interpolation for both displacement and pressure fields throughout the entire range of drainage responses. Numerical simulations reproduce the hydromechanical response of Opalinus shale in one-dimensional consolidation tests throughout the range of primary and secondary consolidation under different external loads. Numerical simulations of the consolidation of a rectangular domain subjected to a strip load demonstrate the efficacy of the proposed stabilization scheme, as well as illustrate the impacts of stress history, mass transfer, and different pore systems on the hydromechanical response.
Development of Time Lapse VSP Integration Workflow: A Case Study at Farnsworth CO2-EOR Project
Abstract This study aims to develop a 4D Vertical Seismic Profile (VSP) integration workflow to improve the prediction of subsurface stress changes. The selected study site is a 5-spot pattern within the ongoing CO2-EOR operations at the Farnsworth Field Unit FWU in Ochiltree County, Texas. The specific pattern has undergone extensive geological and geomechanical characterization through the acquisition of 3D seismic data, geophysical well logs, and core. This workflow constrains a numerical hydromechanical model by applying a penalty function formed between "modeled" versus "observed" time-lapse compressional and shear seismic velocity changes. Analyses of geophysical logs and ultra-sonic measurements on core exhibit measurable sensitivities to changes in both fluid saturation and mean effective stress. These data are used to develop a site-specific rock physics model and stress-velocity relationship, which inform the numerical models used to generate the "modeled" portion of the penalty function. The "observed" portion of the penalty function is provided by a novel elastic full-waveform inversion of the available 3D baseline and three monitor surveys to produce high-quality estimates of time-lapse compressional and shear seismic velocity changes. The modeling workflow accounts sequentially for fluid substitution and stress impacts. Hydrodynamic and geomechanical properties of the 3D coupled numerical model are estimated through geostatistical integration of well log and core data with 3D seismic inversion products. Changes in seismic velocities due to fluid substitution are computed using the Biot-Gassmann workflow and site-specific rock physics. Stress impacts on time-lapse seismic velocity changes are modeled from the effective stress output of the hydromechanical model and are initially based on the velocity versus effective stress relationship extracted from core mechanical testing. Based on the principle of superposition of seismic wavefields, seismic velocity changes attributed to fluid substitution and that due to changes in mean effective stress are treated as linearly additive. The modeled results are upscaled using Backus averaging to reconcile scale discrepancies between the modeled and measured datasets to formulate the penalty function. This manuscript presents the forward modeling process and concludes that for the base case, the seismic velocity changes due to mean effective stress dominates over the seismic velocity changes attributed to fluid substitution because of the extensive range of the pressure perturbations. Successful minimization of this penalty function calibrates the coupled hydrodynamic geomechanical numerical model and affirms the suitability of acoustic time-lapse measurements such as 4D-VSP for geomechanical calibration.
Wellbore cement fracture permeability as a function of confining stress and pore pressure
Cemented annulus fractures are a major leakage path in a wellbore system, and their permeability plays an important role in the behavior of fluid flow through a leaky wellbore. The permeability of these fractures is affected by changing conditions including the external stresses acting on the fracture and the fluid pressure within the fracture. Laboratory gas flow experiments were conducted in a triaxial cell to evaluate the permeability of a wellbore cement fracture under a wide range of confining stress and pore pressure conditions. For the first time, an effective stress law that considers the simultaneous effect of confining stress and pore pressure was defined for the wellbore cement fracture permeability. Here the results showed that the effective stress coefficient (λ) for permeability increased linearly with the Terzaghi effective stress ( -p) with an average of λ = 1 in the range of applied pressures. The relationship between the effective stress and fracture permeability was examined using two physical-based models widely used for rock fractures. The results from the experimental work were incorporated into numerical simulations to estimate the impact of effective stress on the interpreted hydraulic aperture and leakage behavior through a fractured annular cement. Accounting for effective stress-dependent permeability through the wellbore length significantly increased the leakage rate at the wellhead compared with the assumption of a constant cemented annulus permeability.
Effects of Aluminum Plate Initial Residual Stress on Machined-Part Distortion
Dimensional tolerances for high-speed-machined aluminum products continue to tighten due to the demand for automated assembly of complex monolithic parts in aerospace and other industries. Understanding the contribution of inherent residual stress in wrought Al 7050-T7451 plate, common in aircraft manufacture, to distortion of high-aspect-ratio machined parts is critical but remains problematic due to the alloy's low residual stress magnitude over large geometries. Prior investigations into residual stress effects on machined part distortion suffer inadequate characterizations of the wrought material stress field, either because of low fidelity due to “slitting” methods, confounding effects in machined-layer removal methods, or small sample size when using neutron diffraction (ND). In this work, inherent residual stress is measured via ND at 860 locations in a 90.5 mm thick Al 7050-T7451 plate having dimensions 399 mm in the rolling direction and 335 mm in the transverse direction. Unlike prior studies, the ND residual stress is reconstructed using an iterative algorithm to ensure fully compatible, equilibrated 3D field prior to examining its effect on distortion. Further, the findings from simulations and experiments show that inherent residual stress alone could distort a high-aspect-ratio part beyond aerospace industry requirements, that slitting measurements may not sufficiently characterize residual stress for predicted distortion, and that parts machined from different plate thickness locations could exhibit reversed distortion patterns. Thus, research into distortion prediction that considers machining should carefully characterize and reconstruct inherent residual stress so that the coupled machining effects are accurately modeled.
The effects of stress, temperature and facet structure on growth of { 10 1 ¯ 2 } twins in Mg: a molecular dynamics and phase field study
The present study focuses on the roles of structure, stress and temperature on the mobility of facets bounding twin domains and on the overall impact on the growth kinetics of { 10 1 ¯ 2 } twins in Mg. We first use molecular dynamics simulations (MD) to study the relationship between interface structures, stress, temperature and mobility. As facets can exhibit distinct atomic scale structures, we also quantify the effect of stress relaxation mediated by misfit dislocations on the mobility of the basal-prismatic interface. The information gathered for ten different facets is used to calibrate an anisotropic phase field model predicting the effects of temperature and stress on the kinetics and shape of a growing { 10 1 ¯ 2 } 3D twin, by comparison with MD data. The comparative analysis suggests that, as compared to simulations of individual facets, the effective mobility of facets can be significantly different during growth of a 3D twin. This can be attributed to the complex three-dimensional internal stress state of the twin and to coupling in facet motion triggered by facet junctions. As such, by directly fitting both the twin growth kinetics and morphology as predicted by MD for different stresses and temperatures, we can identify a limited set of rate-limiting facets that control twin growth and morphology, and measure their effective mobility.
Stepwise Dynamic Calibration of a Hydromechanical Simulation Using Time-Lapse Vertical Seismic Profile
This study aims to develop a methodology for calibrating subsurface stress changes through time-lapse Vertical Seismic Profiling (VSP) integration. The selected study site is the 13-10A injector well within the ongoing CO2-EOR operation of the Farnsworth Field Unit. The Time-lapse VSP dataset carries the combined effects of fluid substitution and mean effective stress changes, thereby providing a dataset amenable for the calibration of production and injection-induced stress changes. The concept is similar to calibrating a reservoir simulation model in that the process honor real field data to set up an inverse problem. The solution optimizes the independent and impactful geomechanical parameters that replicate the observed time-lapse seismic velocity changes. This stress calibration is enabled by 4D geomechanical modeling and the VSP Integration workflow. This calibration benefits from extensive geological, geophysical and geomechanical characterization through 3D seismic data, geophysical well logs, and core assessed as part of the 1D MEM conducted on the 13-10A subject well. These data are used to develop a site-specific rock physics model. The Biot Gassmann workflow combines rock physics and reservoir simulation outputs to determine the fluid substitution contribution to seismic velocity change. Additionally, modeled seismic velocity attributed to mean effective stress are determined from the geomechanical simulation outputs, and the stress-velocity relationship developed from the ultrasonic seismic velocity measurements on the extracted Morrow B core. A penalty function is then formed between the modeled seismic velocities and the observed time-lapse VSP dataset. Four independent and impactful geomechanical parameters have been determined. These are the bulk modulus and shear modulus for zero porosity and the shear and compressional seismic velocity to mean effective stress derivatives. The dataset of numerous coupled hydromechanical- geomechanical simulation realizations is built by combining variations of the four stated geomechanical parameters. A machine learning-assisted workflow comprised of an artificial neural network and a particle swarm optimizer are used to converge on the optimal geomechanical parameters. The successful execution of this workflow has affirmed the suitability of acoustic time-lapse measurements for 4D-VSP geomechanical stress calibration pending measurable stress sensitivities within the anticipated effective stress changes and the availability of suitable and reliable datasets for petroelastic modeling.
DOE Award No.: DE-FE0023919 Phase 3 Scientific/Technical Report
This is the Phase 3 Report for the ‘Deepwater Methane Hydrate Characterization and Scientific Assessment or Genesis of Methane Hydrates in the Gulf of Mexico (GOM2)’ research project (DOE Award No. DE-FE0023919). The report summarizes activities from January 16, 2018 to September 30, 2019. The project is led by the University of Texas at Austin (UT). The objective is to gain insight into the nature, formation, occurrence and physical properties of methane hydrate-bearing sediments for the purpose of methane hydrate resource appraisal through the planning and execution of drilling, coring, logging, testing and analytical activities that assess the geologic occurrence, regional context, and characteristics of marine methane hydrate deposits in the Gulf of Mexico outer continental shelf (OCS). We determined that it would not be possible to pursue the project with the International Ocean Discovery Program (IODP) and then developed a revised science and operations plan to maximize the science using a commercial vessel within the budget that is available. We improved the scientific capability of the UT pressure core center by adding the ability to X-ray pressure cores, and the ability to cut and store multiple core samples for experimental analysis. We transferred pressure cores to peer institutions, including the National Energy Technology Lab (NETL) the United States Geological Survey (USGS Woods Hole). We made advances in understanding the composition and source of gasses locked in the methane hydrate from Green Canyon 955 (GC-955). We developed approaches to determine the in-situ salinity of hydrate-bearing samples, revealing that the in situ salinity of the GC 955 reservoir is just above that of seawater. We determined that the GC-955 hydrate reservoir is composed of sandy silt with a high concentration of hydrate interbedded with clayey silt with no hydrate present. We determined that some or all of the intervals bounding the hydrate reservoir are composed of material similar to the reservoir (sandy silt with hydrate and clayey silt with no hydrate), but with a lower net to gross. We determined that the reservoir effective permeability at a hydrate saturation of ~80-90 % is ~0.1 mD (or 1.0×10-16 m2) to ~0.5 mD (or 5.0×10-16 m2) and the intrinsic permeability is ~12 mD (or 1.2×10-14 m2). The GC-955 reservoir in-situ porosity of sandy silt is 0.38 to 0.40 and is largely independent of effective stress, and that the porosity of clayey silt is 0.33 at an in situ effective stress of 3.8 MPa to ~0.37 at zero effective stress. We developed a systematic, repeatable approach to studying hydrate reservoir properties by reconstituting individual lithofacies from dissociated pressure cores. We determined index properties of GC-955 reservoir, including liquid limit and plasticity, porosity, capillary behavior, and particle size distribution. We developed a more robust pressure coring technology. We finalized and published the UT-GOM2-1 Expedition Volume and we finalized a dedicated volume on the UT-GOM2-1 expedition that will be published in the American Association of Petroleum Geologists Bulletin (AAPG) in spring 2020.
A Novel Multiphysics Multiscale Multiporosity Shale Gas Transport Model for Geomechanics/Flow Coupling in Steady and Transient States
Summary A novel multiphysics multiscale multiporosity shale gas transport (M3ST) model was developed to investigate shale gas transport in both transient and steady states. The microscale model component contains a kerogen domain and an inorganic matrix domain, and each domain has its own geomechanical and gas transport properties. Permeabilities of various shale cores were measured in the laboratory using a pulse decay permeameter (PDP) with different pore pressure and confining stress combinations. The PDP-measured apparent permeability as a function of pore pressure under two effective stresses was fitted using the microscale M3ST model component based on nonlinear least squares fitting (NLSF), and the fitted model parameters were able to provide accurate model predictions for another effective stress. The parameters and petrophysical properties determined in the steady state were then used in the transient-state, continuum-scale M3ST model component, which performed history matching of the evolutions of the upstream and downstream gas pressures. In addition, a double-exponential empirical model was developed as a powerful alternative to the M3ST model to fit laboratory-measured apparent permeability under various effective stresses and pore pressures. The developed M3ST model and the research findings in this study provided critical insights into the role of the multiphysics mechanisms, including geomechanics, fluid dynamics and transport, and the Klinkenberg effect on shale gas transport across different spatial scales in both steady and transient states.
Fracturing in Wet Granular Media Illuminated by Photoporomechanics
Here, we study fluid-induced deformation and fracture of granular media and apply photoporomechanics to uncover the underpinning grain-scale mechanics. We fabricate spherical photoelastic particles of 2-mm diameter to form a monolayer granular pack in a circular Hele-Shaw cell that is initially filled with a viscous fluid. The key distinct feature of our system is that, with spherical particles, the granular pack has a connected pore space, thus allowing for pore-pressure diffusion and the study of effective stress in coupled poromechanical processes. We inject air into the fluid-filled photoelastic granular pack, varying the initial packing density and confining weight. With our recently developed experimental technique, photoporomechanics, we find two different modes of fluid invasion: fracturing in fluid-filled elastic media (with strong photoelastic response) and viscous fingering in frictional fluids (with weak or negligible photoelastic response). We directly visualize the evolving effective stress field and discover an effective stress shadow behind the propagating fracture tips, where the granular pack exhibits undrained behavior. We conceptualize the behavior of the system by means of a mechanistic model for a wedge of the granular pack bounded by two growing fractures. The model captures the pore-pressure build-up inside the stress shadow region and the grain compaction in the annular region outside. Our model reveals that a jamming transition determines the distinct rheological behavior of the wet granular pack, from a friction-dominated to an elasticity-dominated response.