Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “thermal fracturing”

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

Time-dependent THMC properties and microstructural evolution of damaged rocks in excavation damage zone

Modeling coupled thermo-hydro-mechanical-chemical (THMC) processes in host rocks near high-level nuclear waste (HLW) repositories at various time scales is an extremely challenging task. The current study integrates experimental, theoretical, and numerical methods in assessing the evolution of excavation damage zone (EDZ) over time and its implication on the long-term migration of hazardous species. Argillite and rock salt and are the focus of this study. The first part of the report presents a novel time-dependent directional microcrack damage theory for generic brittle rocks. It features detailed statistical description of the microcracks within a damaged solid, permitting a direct upscaling of microscale processes such as crack growth kinetics, crack closure/opening, sliding friction to explain the macroscopic creep, nonlinear elasticity, shear dilation, and loading-unloading hysteresis. This provides a basic platform for describing the anisotropic mechanical and transport properties of damaged rocks during excavation and subsequent THMC loadings. The model is validated and numerically implemented to Finite Element (FE) package ABAQUS through the user-defined material (UMAT) interface and have demonstrated great potential in resolving the time-dependent and anisotropic evolution of damage in the EDZ. The second part of the report presents a multi-scale experimental effort in characterizing the thermal, hydraulic, and mechanical properties of Mancos shale and Avery Island salt. For the Mancos shale, triaxial compression tests are performed at different confining pressures and temperatures to probe its thermomechanical properties relevant to HLW repositories. The obtained stress-strain data are interpreted using the proposed directional damage theory. Post-test specimens are subjected to gas permeability tests to reveal the correlation between permeability and the degree of microcracking. At microscale, temperature-controlled nanoindentation tests are performed and found a linear correlation between fracture toughness and elastic modulus from 25°C to 300°C. For Avery Island salt, we have designed and manufactured a novel relative-humidity controlled uniaxial creep device. Long-term creep tests at low stresses (< 5 MPa) are performed at different levels of relative humidity (RH). Besides confirming the much higher creep rates as one would expect through extrapolating the high-stress creep data, the results reveal that the steady-state creep rate of rock salt is strongly dependent on the ambient RH, an aspect that is often neglected in the literature. Both behaviors can be attributed to the pressure-solution creep mechanism which dominates at low stress and high RH levels. The third part of the report explores a set of numerical strategies in modeling the THMC behavior of porous geomaterials. A fully implicit, monolithic FE solution that can flexibly interface with different material models and coupling mechanisms for THMC problems is developed and verified through the ABAQUS user-defined element (UEL) interface. The scheme is used to study the THM response of a hypothetical HLW storage site with reference to an existing in-situ heater test. Strategies for integrating the UEL and the microcrack UMAT are suggested. Another numerical endeavor of this study is to implement a higher-order asymptotic homogenization method to account for the heterogeneous porous structures. The same method is then extended to perform microstructure-informed thermo-mechanical modeling of generalized continua. The above outcomes of this project provide a strong thrust towards enhancing the fundamental understanding and modeling capability of the long-term evolution of host rocks in EDZ, thus helping achieve the design goal of 1-million-year isolation of high-level nuclear wastes.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

DECOVALEX-2019 (Task G Final Report)

The DECOVALEX Project is an on-going international research collaboration, established in 1992, to advance the understanding and modeling of coupled Thermal (T), Hydrological (H), Mechanical (M) and Chemical (C) processes in geological in geological systems. This document is the final report of Task G which was proposed by the Swedish Radiation Safety Authority (SSM) and coordinated by geomecon GmbH, presenting the technical definitions of the problems studied, approaches applied, achievements made and outstanding issues for future research. Task G in DECOVALEX-2019 is primarily concerned with the evolution of transmissivity throughout the lifetime of a repository for radioactive waste, and in particular of spent nuclear fuel, in sparsely fractured and competent crystalline rock masses, and the representation of these processes via numerical simulation. The numerical simulations were to be evaluated for their suitability, and to be validated. Over the course of the work presented here, not only were transmissivity changes and the best approaches to simulate those change examined, but also strategies were to be developed on how to monitor transmissivity change. Guidelines on what a control program and a monitoring system for a repository should include and be able to measure, and how this can be done in practice, was a key objective of the research project.

58 GEOSCIENCES↗
DOE OSTI2021

Creating Hydride Rim Structures in Zircaloy 4 Cladding Tubes to Mimic High Burnup

During their long service in commercial light water reactors, zirconium alloy cladding tubes are subject to oxidation and subsequent hydrogen pickup which can degrade the materials mechanical properties [1]. These degradations can lead to premature cladding ruptures during reactor power transients [2] as well as negatively impact the claddings post quench ductility following high temperature steam oxidation in thermal hydraulic transients [3]. Hydrogen has low solid solubility in zirconium and usually precipitates as brittle zirconium hydride platelets in the cladding metal matrix. Due to temperature driven diffusion (Soret effect) and the grain texture of stress relieved zirconium metal tubes, these hydrides precipitate with a circumferential orientation near the outer surface of the cladding tubes forming a hydride rim which can serve as a site for crack initiation during circumferential ‘hoop’ loading of the cladding [4rim formation][5]. In order to better understand the effect of hydride rims on the fracture behavior of stress relieved zirconium alloy cladding tubes a procedure to produce these rim-like structures is being developed at Idaho National Laboratory (INL).

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Compatibility of mesophase pitch and linear low-density polyethylene for low-cost carbon fiber

Mesophase pitch and linear low-density polyethylene (LLDPE) were blended by shear mixing and extruded producing precursor fibers which were then converted to carbon fibers (CFs). LLDPE additives agglomerate within the isotropic spheres rather than uniformly dispersing throughout the entire matrix. The blends possess different thermal properties and precursor tensile properties evident by two-phase tensile stress nature not present in precursor fibers produced from pure pitch. The resultant CFs were oxidized at 280 °C for 4 h and carbonized at 1000 °C for 15 min. This process created pitch CFs with a tensile strength of 1.23 GPa and a Young’s modulus of 99 GPa. All CFs produced from blends of pitch and LLDPE possessed a tensile strength of 0.63 GPa and less. These CFs showed a microstructure slightly different from the pure pitch CFs with fractured CF cross sections revealing voids present in all blended CFs due to burnt off LLDPE stemming from carbonization. Further study into compatibilizers and new conversion processes is required to successfully blend LLDPE with mesophase pitch for the purpose of producing a new low-cost CF.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Tribological behavior of spark plasma sintered Ti 3 SiC 2 MAX phase composites

Ti₃SiC₂ MAX phases are considered promising candidates for tribological applications; but, their low hardness (5-6 GPa) can lead to increased abrasive wear and higher wear rates. This study investigates the effect of incorporating hard SiC particles on microstructure, mechanical, and tribological behavior of Ti₃SiC₂-SiC-based MAX phase composites prepared using spark plasma sintering at 1400 °C, 40 MPa, 15-minute holding time. Phase and microstructural characterization of composites confirmed the formation of Ti₃SiC₂ MAX phase (90%) along with TiC as a minor phase (10%). For Ti₃SiC₂-SiC, homogeneous distribution of SiC grains within Ti₃SiC₂ matrix resulted in increased hardness from ~11.6 GPa to ~14.8 GPa; however, the flexural strength decreased from ~615 MPa to ~597 MPa due to coefficient of thermal expansion mismatch. Tribological behavior of Ti₃SiC₂-SiC MAX phase composites was assessed using unidirectional sliding ball-on-flat tests with a 52100 steel ball at different loads and a constant speed of 0.05 m/s. SiC reinforced composites exhibited a decreased friction coefficient from 0.39 to 0.29 and wear rates comparable to Ti₃SiC₂ (7-10×10⁻³ mm³/N·m). SEM and EDS analysis of the worn surfaces indicated material transfer from the steel ball counter body and its oxidation due to frictional heating, at lower loads (2N and 5N), whereas at higher load (10N), fracture of the transfer layer dominated, with the presence of microcracks, delamination, grain pull-outs, and wear debris. Improved mechanical properties and good adhesion of SiC with Ti₃SiC₂ matrix resulted in reduced microcracks and grain pull-outs, making Ti₃SiC₂-SiC/Steel tribo-pair more suitable for tribological applications.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Controlled Spalling of 4H Silicon Carbide with Investigated Spin Coherence for Quantum Engineering Integration

We detail scientific and engineering advances which enable the controlled spalling and layer transfer of single crystal 4H silicon carbide (4H-SiC) from bulk substrates. 4HSiC’s properties, including high thermal conductivity and a wide bandgap, make it an ideal semiconductor for power electronics. Moreover, 4H-SiC is an excellent host of solid-state atomic defect qubits for quantum computing and quantum networking. Because 4H-SiC substrates are expensive (due to long growth times and limited yield), techniques for removal and transfer of bulk-quality films are desirable for substrate reuse and integration of the separated films. In this work, we utilize updated approaches for stressor layer thickness control and spalling crack initiation to demonstrate controlled spalling of 4H-SiC, the highest fracture toughness crystal spalled to date. We achieve coherent spin control of neutral divacancy (VV 0 ) qubit ensembles and measure a quasi-bulk spin T 2 of 79.7 μs in the spalled films.

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

Methane Plume Emissions Associated With Puget Sound Faults in the Cascadia Forearc

Methane gas plumes have been discovered to issue from the seafloor in the Puget Sound estuary. These gas emission sites are co-located over traces of three major fault zones that fracture the entire forearc crust of the Cascadia Subduction Zone. Multibeam and single-beam sonar data from cruises conducted in 2011, 2018, 2019, 2020 and 2021 identified the acoustic signature of over 330 bubble plumes. Dissolved gas from the plumes combines to elevate seawater methane concentrations of the entire Puget Sound estuary. Fluid samples from adjacent terrestrial hot springs and deep-water wells surrounding the estuary contain a helium-3 isotope signature, indicating a deep fluid source located near the underlying Cascadia Subduction Zone. However, Puget Sound seawater emission sites lack chemical isotope signatures and elevated thermal anomalies expected from association with a deep plate-interface reservoir. The existence of vigorous marine methane plumes located in areas of thin sediment cover and above the traces of deeply-penetrating forearc fault zones – but still presenting no thermal or chemical anomalies that are found in many other similar forearc environments and remain an unresolved paradox.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Mutually Reinforced Polymer–Graphene Bilayer Membranes for Energy‐Efficient Acoustic Transduction

Abstract Graphene holds promise for thin, ultralightweight, and high‐performance nanoelectromechanical transducers. However, graphene‐only devices are limited in size due to fatigue and fracture of suspended graphene membranes. Here, a lightweight, flexible, transparent, and conductive bilayer composite of polyetherimide and single‐layer graphene is prepared and suspended on the centimeter scale with an unprecedentedly high aspect ratio of 10 5 . The coupling of the two components leads to mutual reinforcement and creates an ultrastrong membrane that supports 30 000 times its own weight. Upon electromechanical actuation, the membrane pushes a massive amount of air and generates high‐quality acoustic sound. The energy efficiency is ≈ 10–100 times better than state‐of‐the‐art electrodynamic speakers. The bilayer membrane's combined properties of electrical conductivity, mechanical strength, optical transparency, thermal stability, and chemical resistance will promote applications in electronics, mechanics, and optics.

Khan, Assad U.↗
DOE OSTI2022

Penta-SiCN: A Highly Auxetic Monolayer

We report the negative Poisson’s (NPR) ratio in a two-dimensional (2D) material is a counterintuitive mechanical property that facilitates the development of nanoscale devices with sophisticated functionality. Inspired by the peculiar buckled lower-symmetric, trilayered geometry of pentagonal monolayers, we theoretically predict penta-SiCN, a ternary auxetic metallic monolayer with highly tunable NPR. The penta-SiCN is structurally, thermally, dynamically, and mechanically stable, and sustainable at and beyond room temperature with experimental feasibility. It possesses nontrivial geometrical and mechanical isotropy and relatively moderate thickness. Remarkably, the shorter and quasi sp3-hybridized C–N bond and the rigidity against the strain allow the monolayer to possess a high value of NPR (-0.136), even higher than that of black phosphorene, extendable up to -0.639 by 4% of biaxial stretching. On the other hand, the 2D Young’s modulus of 129.88 N/m decreases to 41.34 N/m at equivalent stretching, indicating relative softening and flexibility. Interestingly, a buckled-to-planar phase transition is identified at 10% biaxial strain before it suffers the fracture at 16%. Additionally, the strong optical anisotropy, absorbance (up to 6.51 x 105 cm -1 ), and presence of plasmon frequency demonstrate its potential application in optomechanical and plasmonics.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Guidelines for predicting stress in cemented doublets undergoing temperature change

This work explores quick predictive methods for calculating potentially risky stresses and deflections in cemented doublets experiencing temperature change that agree well with finite element analysis. There are three failure modes of interest: cohesive failure of the adhesive, delamination (surface bond failure or debonding), and glass fracture. Adhesive theory, confirmed by finite element analysis, predicts stress singularities that complicate interpretation of the stress calculations. The presence of a stress singularity indicates the breakdown of linear elastic assumptions, but damage initiation and stress singularities are related. The authors find that geometry details near a bond edge can exacerbate or minimize damage initiation and stress concentrations. Because the interpretation of the stress results is complicated, the authors investigated predicted stresses in doublets that have been successfully tested between -40 and 85 °C. This study found that the product (ΔT ∙ Δα) should be less than 189 ppm, where ΔT is the temperature excursion and Δα is the difference in glass coefficient of thermal expansion. If the product (ΔT ∙ Δα) is equal to or greater than 189 ppm, further analysis and testing is warranted. But the authors also show that the fabrication process can significantly influence stress failure, particularly with large diameter doublets.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Improving Ni-Based SOFC Anode Resilience and Durability through Secondary Phase Formation

The project used a multi-disciplinary approach to test test if 2° materials, added to traditional cermet components, enhanced SOFC mechanical strength, anode resilience and durability when subjected to electrochemical and environmental redox cycling. Findings from fracture testing, electrochemical studies and operando Raman spectroscopy measurements demonstrated that aluminum titanate (Al 2 TiO 5 or ALT) mechanically mixed with NiO-YSZ cermets confers superior mechanical strength and slower degradation rates under polarization. No deleterious effects from ALT in conversion efficiency were observed. Specific accomplishments included the following: (1) Refined methods used to fabricate ALT enhanced anodes into bi-layer anode supports to achieve high power densities; (2) Compared the effects of adding ALT mechanically to Ni-YSZ powders prior to anode fabrication with adding ALT through infiltration and co-infiltration of YSZ scaffolds; (3) Tested the resilience of these novel materials to electrochemical and environmental redox cycling and thermal stresses commonly encountered in functioning SOFCs; and (4) Characterized the ability of ALT to improve the carbon tolerance of Ni-YSZ anodes operating with dry methane.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2022

Impact Response of Control Atmosphere Plasma Spray Deposited Materials

Thermal spray processing of metals and respective blends is becoming increasingly attractive due to the unique properties such as increased yield strength, low ductility, and differences in tensile and compressive strengths that result from microstructural features due to the spray process compared to other additive manufacturing methods. Here we report the results of plate impact experiments applied to Controlled Atmosphere Plasma Spray deposits of tantalum (Ta), niobium (Nb), and a tantalum-niobium blend (TaNb). These methods allowed for definition of the Hugoniot for each material type and the assessment of the Hugoniot Elastic Limit (HEL). Spallation experiments were conducted, and soft recovery of each material type allowed for scanning electron microscopy to characterize the fracture mechanism during tensile loading.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Postweld Heat Treatment on Toughness of Electric- Resistance Welded X70 Steel

In this report the effect of postweld heat treatment (PWHT) on Charpy V-notch impact toughness of a highfrequency electric-resistance welded (HF-ERW) grade X70 pipeline steel is investigated. PWHT thermal cycles are simulated using the Gleeble on the as-welded specimens. Microstructure, along with crystallographic texture and microhardness, is characterized in specimens that are impact tested at –5°C (23°F), –30°C (–22°F), and –45°C (–49°F). The impact toughness values show a wide scatter band and decrease with the increasing peak temperature of the PWHT. For higher peak temperatures, the microstructure of the heat-affected zone (HAZ) gradually changed from equiaxed ferrite to bainitic ferrite. Furthermore, the prior austenite grain size (PAGS) increases with the increasing peak temperature. The density of high-angle grain boundaries decreases, and the fraction of cleavage planes {100} parallel to the impact fracture plane increases for higher peak temperatures of the PWHT

36 MATERIALS SCIENCE↗
DOE OSTI2024

Stress Relaxation Cracking in 347H Stainless Steel Arc Welds: Susceptibility Evaluation of Heat-Affected Zone

Stress relaxation cracking (SRC) is considered one of the major failure mechanisms for 347H stainless steel welds at elevated service temperatures or during post weld heat treatment (PWHT), especially within the heat-affected zone (HAZ). This work focuses on the characterization of SRC susceptibility within 347H physically simulated arc welded HAZ at elevated temperatures. A four-step SRC thermomechanical test in combination with finite element modeling (FEM) of the welding and testing processes is developed to establish a susceptibility map for HAZ. The test first runs a thermal cycle with three different peak temperatures (1335, 1275, and 1150 °C) to duplicate representative HAZ subzone microstructures, followed by time-to-failure examination under a variety of pre-stress (260–600 MPa) and pre-strain conditions (0.03–0.19) as a function of reheat temperatures between 750 and 1050 °C. With the aid of FEM, SRC susceptibility maps are generated to identify the threshold stress, plastic strain, and creep strain as a function of test temperature. It was found out that HAZ subzone with a lower peak temperature (1150 °C) appears to be slightly less susceptible to SRC than the other two subzones that experienced higher peak temperatures. Generally, time-to-fracture reduces with increasing initially applied stress and strain for all test temperatures. The pre-stress thresholds decrease from about 500 to 330 MPa as the testing temperature increases from 800 to 1050 °C, while the corresponding initial plastic strain thresholds reduces from 0.15 to 0.06. The SRC susceptibility was also evaluated through the Larson–Miller Parameter (LMP) analysis as a function of plastic strain, initial stress and starting stress upon reaching the testing temperature, respectively. The 1050 °C test with a high pre-applied strain (0.1) exhibits an extremely short time to failure (t = 3 s) that lies outside the general trend in LMP analysis. Additionally, it was identified that a plastic strain above 0.07 is identified to significantly reduce the bulk creep strain tolerance to fracture and therefore increases SRC susceptibility. Hardness measurement and fractography analysis indicated that the strain aging of niobium carbonitrides and other potential phases in conjunction with intergranular precipitates contributes to an increase in microhardness and increased intergranular cracking susceptibility.

14 SOLAR ENERGY↗
DOE OSTI2020

A Dual Approach in Direct Ink Writing of Thermally Cured Shape Memory Rubber Toughened Epoxy

Bisphenol A-based epoxies are much used in a wide range of composite and coating applications due to their excellent thermomechanical properties. However, their 3D printability remains a challenge with most reported materials suffering from high brittleness and low toughness. In this work, we have described especially modified epoxy resins that enable 3D printing with both fast and slow curing rates. These materials exhibit greatly enhanced toughness, tunable thermomechanical properties, and excellent shape memory behavior. Two different printing systems, including a two-part static mixing printhead and a single extrusion printhead, were developed for fast- and slow-curing epoxies, respectively. The rheology of inks in both systems has been modified into printable thixotropic fluids with the aid of silica nanoparticles and other additives. Epoxide-functionalized telechelic polybutadiene was added into the resins, which are then introduced inside the epoxy network after cross-linking. The addition of polybutadiene rubber significantly improves the toughness (over 135%), fracture strain (over 200%), and shape memory behavior. By adding different amounts of the rubber telechelic, thermomechanical properties, including modulus, elongation, and Tg of epoxy, can be well controlled in a wide range to satisfy different applications.

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

Intermediate-Term Thermal Aging Effect Evaluation for Grade 92 and 316L at The LWR Relevant Temperature

Life extension of the existing nuclear reactors imposes accumulated damages, such as higher fluences and longer period of corrosion, to structural materials, which would result in significant challenges to the traditional reactor materials such as type 304 and 316 stainless steels. Advanced alloys with superior radiation resistance will increase safety margins, design flexibility, and economics for not only the life extension of the existing fleet but also new builds with advanced reactor designs. The Electric Power Research Institute (EPRI) teamed up with Department of Energy (DOE) on the Advanced Radiation Resistant Materials (ARRM) program, aiming to develop and test degradation resistant alloys from current commercial alloy specifications by 2021 to a new advanced alloy with superior degradation resistance in light water reactor (LWR)-relevant environments by 2024. Based on a comprehensive microstructure and property screening performed in Phase-1 of the ARRM program, a total of five alloys, together with 316L and X-750 as references, were down-selected for Phase-2 neutron irradiation studies. Because thermal aging could exert a synergistic effect on neutron irradiation due to the low neutron damage rate on the order of 10 –7 displacements per atom per second (dpa/s), Grade 92 (one of the five down-selected alloys) and two heats of 316L were selected in this task to study the effect of aging at 350°C for 12.5–12.7 kh on microstructure and mechanical properties. The aging time is approximately corresponding to 5 dpa neutron irradiation, which can be served as a reference for the 5-dpa-irradiated samples of the alloys to help understand the independent neutron damage influence on microstructure and mechanical properties. Optical microscopy and scanning electron microscopy were used for microstructural characterization. Hardness, tensile, Charpy impact toughness, and fracture toughness in the ductile regime of the aged samples were examined. The aging of Grade 92 led to the formation of many Laves phase in sizes of ~100–200 nm and resulted in some reduction in hardness and yield/ultimate tensile strength with some increases in uniform and total plastic elongations, which may have helped the enhancement of Charpy impact toughness, e.g., ~4 J increase in upper-shelf energy and 20.5°C reduction in ductile-brittle transition temperature compared with the unaged condition. The fracture toughness of the aged Grade 92 showed decent toughness of ~303 MPa√m (K Jq ) with ~87 tearing modulus at 22°C, which reduced to ~241 MPa√m (~20% reduction) with ~75 tearing modulus at 300°C. The two heats of 316L are differentiated by their amounts of δ-ferrite, i.e., ~1 vol% in heat T1103 and ~4 vol% in heat N5B8. The two heats of 316L were subjected to 15% cold work (CW) as this condition is often used in nuclear reactors to help trapping radiation-induced defects and thus postpone the steady state swelling stage. The aging did not result in any noticeable microstructure changes, except for possible segregations at grain boundaries that need to be further investigated. The aging resulted in slight reduction in hardness of 316L-T1103 and slight increase in hardness of 316L-N5B8, which is consistent with their tensile testing results. The 15%CW might have introduced some inhomogeneity, resulting in large standard deviations in hardness, which were reduced after the aging. The aged 316L followed the same trend as the unaged condition for the strength and elongation results from the tensile testing. Minimum elongations appeared within ~300–600°C, above which the elongations seem to increase. Unlike the negligible or minor changes in microstructure, hardness and tensile results of the two heats of 316L, their fracture toughness showed noticeable difference. The aged 316L-T1103 showed good toughness of ~350 MPa√m with ~80 tearing modulus at 22°C, which reduced to ~271 MPa√m (~23% reduction) with ~83 tearing modulus at 300°C. In contrast, the aged 316L-N5B8 showed lower toughness of ~244 MPa√m with ~53 tearing modulus at 22°C, which reduced to ~180 MPa√m (~26% reduction) with ~31 tearing modulus at 300°C. The preliminary results indicate that the presence of high volume of δ-ferrite would noticeable impair fracture toughness although it may not influence hardness and tensile properties.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Are we modeling the properties of unconventional shales correctly?

Technological advancements have led to impressive growth in hydrocarbon (HC) production from unconventional shale reservoirs. However, major limiting factors in shale gas development are i) low recovery of gas in place (<20%), (ii) a rapid decline in well productivity, iii) release of contaminants, despite the use of advanced hydraulic fracturing fluids and multiple fracturing stages. This is mainly due to a lack of understanding of the nanoscale physicochemical properties of shales, especially that of kerogen, the macromolecule that is not only the source of the majority of the HC’s in shales but also holds most of the HC’s in adsorbed state. Over the last few years, a large number of studies have been published on simulating the physicochemical properties of shale utilizing molecular models of kerogen. However, the molecular models of kerogen input into these simulations are based on the kerogen “type” derived from a very limited number of shale samples. In this paper, we examine the variations that can exist in kerogen structure within a particular kerogen “type” across different shale basins and single shale basin at similar thermal maturity levels. We propose that using kerogen “type” based structural models for molecular simulations could lead to inaccurate estimation of HC reserves, HC recovery, HC production, frackability, and quantity and quality of produced water. Our study highlights the need for developing a better classification of kerogen based on its “molecular structure” instead of “type” for a more accurate prediction of physicochemical properties of shales.

04 OIL SHALES AND TAR SANDS↗
DOE OSTI2023

A statistical mechanics framework for polymer chain scission, based on the concepts of distorted bond potential and asymptotic matching

To design increasingly tough, resilient, and fatigue-resistant elastomers and hydrogels, the relationship between controllable network parameters at the molecular level (bond type, non-uniform chain length, entanglement density, etc.) to macroscopic quantities that govern damage and failure must be established. Many of the most successful constitutive models for elastomers have been rooted in statistical mechanical treatments of polymer chains. Typically, such constitutive models have used variants of the freely jointed chain model with rigid links. However, since the free energy state of a polymer chain is dominated by enthalpic bond distortion effects as the chain approaches its rupture point, bond extensibility ought to be accounted for if the model is intended to capture chain rupture. To that end, a new bond potential is supplemented to the freely jointed chain model (as derived in the u FJC framework of Buche and Silberstein (2021) and Buche et al. (2022)), which we have extended to yield a tractable, closed-form model of single chain behavior that should be amenable to continuum-level constitutive model development. Inspired by the asymptotically matched u FJC model response in both the low/intermediate chain force and high chain force regimes, a simple, quasi-polynomial bond potential energy function is derived. This bond potential exhibits harmonic behavior near the equilibrium state and anharmonic behavior for large bond stretches tending to a characteristic energy plateau (akin to the Lennard-Jones and Morse bond potentials). Using this bond potential, approximate yet highly-accurate analytical functions for bond stretch and chain force dependent upon chain stretch are established. Then, using this polymer chain model, a stochastic thermal fluctuation-driven chain rupture framework is developed. This framework is based upon a force-modified tilted bond potential that accounts for distortional bond potential energy, allowing for the derivation and subsequent calculation of the dissipated chain scission energy. Here, the cases of rate-dependent and rate-independent scission are accounted for throughout the rupture framework. The impact of Kuhn segment number on chain rupture behavior is also investigated. The model is fit to single chain mechanical response data collected from atomic force microscopy tensile tests for validation and to glean deeper insight into the molecular physics taking place. Due to their analytical nature, this polymer chain model and the associated rupture framework can in the future be implemented in finite element models accounting for fracture and fatigue in polydisperse elastomer networks.

36 MATERIALS SCIENCE↗