Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fracture modelling”

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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Peridynamic modeling of cementitious materials for nuclear waste management

Radioactive and hazardous waste generated from fuel processing plants, nuclear reactors, and hospitals, requires effective management strategies. Cementitious materials are widely applied for these needs, serving as structural materials, reactive barriers, or waste forms. In these applications cracking poses a significant risk to performance, driven by inconsistent shrinkage of the components and the varying strength and permeability of their interfaces. Traditional modeling approaches face challenges in representing the complex fracture behavior of cementitious materials due to the reliance on spatial derivatives and difficulties with mesh generation. Here, to overcome these limitations, we employ peridynamics, a novel continuum mechanics formulation that uses integrals to describe mechanical equilibrium, avoiding discontinuities associated with traditional methods. Through incorporation of a bi-linear softening model and quasistatics, an experimentally validated model for Portland cement concrete samples was created. Mechanical parameters, including compressive strength and elastic modulus were validated and variation due to aggregate packing was evaluated. Additionally, sensitivity analysis of the peridynamic parameters to the bulk material properties was established. The results lay the groundwork for evaluating the impact of unique conditions of cementitious waste forms that can be assessed to improve the reliability of waste management strategies.

Aggregates↗
DOE OSTI2021

The influence of preexisting host rock damage on fault network localization

The transition from stable to unstable fracture propagation occurs when fractures begin to interact and link. Thus, fracture network coalescence controls how rocks and engineered structures fail. Here, to constrain the factors that influence localization in shear zones under brittle conditions, we build discrete element method models with a rough fault embedded in a shear zone. We add varying numbers of diffuse, randomly-placed weaknesses to examine the influence of diffuse damage on fracture network localization. The number of weaknesses controls the localization behavior of the fault network and the final fault geometry. We quantify localization using the Gini coefficient of the fracture volume, which measures the nonuniformity in a population. Each model generally increases in localization toward failure. However, models with more diffuse damage experience delocalization phases that are superimposed on the overall trend of increasing localization. The observed link between delocalization and host rock damage may help explain the varying localization of low magnitude seismicity in southern California. Models with more diffuse damage produce more complex fault geometries comprised of several parallel strands of wing cracks. The propagation of these wing cracks reduces the shear stress acting on the model boundaries, indicating that this fracture development increases the mechanical efficiency of the system.

58 GEOSCIENCES↗
DOE OSTI2020

An energy-based coupling approach to nonlocal interface problems

Nonlocal models provide accurate representations of physical phenomena ranging from fracture mechanics to complex subsurface flows, where traditional partial differential equations fail to capture effects caused by long-range forces at the microscale and mesoscale. However, the application of nonlocal models to problems involving interfaces such as multimaterial simulations and fluid-structure interaction, is hampered by the lack of a rigorous nonlocal interface theory needed to support numerical developments. In this paper, we use an energy-based approach to develop a mathematically rigorous nonlocal interface theory which provides a physically consistent extension of the classical perfect interface PDE formulation. Numerical examples validate the proposed framework and demonstrate the scope of our theory.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Determination of Constraint-Independent Crack Tip Opening Angle for Stable Crack Growth in High-Strength Ductile Steels

The crack tip opening angle (CTOA) is one of fracture toughness parameters that has been used for decades in describing large stable crack growth in thin-walled aerospace structures under the low-constraint conditions. Recently, the pipeline industry has developed a growing interest in using the CTOA parameter to serve as the minimum required fracture toughness to arrest dynamic crack propagation in modern gas transmission pipelines made of high-strength ductile steel. To meet this industrial need, the CTOA test standard ASTM E3039 was therefore developed for measuring the constant critical CTOA. ASTM E3039 recommends a drop weight tearing test (DWTT) specimen with a shallow crack for standard CTOA testing, but its CTOA may depend on the low constraint condition of the DWTT specimen at the crack tip. Verifying the constraint independence of the DWTT-measured CTOA thus becomes indispensable for applying CTOA toughness to the running fracture control in the pipeline design. For this purpose, the present paper evaluates critical CTOA values in a set of fracture toughness tests on single-edge notched bend (SENB) specimens with shallow and deep cracks, based on four CTOA estimation models. Among these, the Ln(P)-LLD linear fit model is similar to that recommended by ASTM E3039 for CTOA calculation. Fracture test data for X80 pipeline steel and HY80 structural steel were considered in the CTOA evaluation. The results showed that the four CTOA models were able to determine a constraint independent CTOA value for stable crack growth in the SENB specimens. As a result, a single, reliable, constant CTOA value could be determined regardless of the specimen geometry or the crack-tip constraint conditions. Therefore, the CTOA measured using ASTM E3039 is constraint-independent and transferable to use in cases of actual cracks propagating in gas transmission pipelines.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Numerical and experimental analysis of mechanically induced failure in electric vehicle battery modules

Mitigating thermal runaway and cell-to-cell propagation is essential for improving the safety of electric and hybrid vehicles. Enhancing digital twin capabilities to predict battery mechanical abuse is particularly critical for automotive and aerospace applications, where crashworthiness is a key concern. Understanding failure conditions and propagation in battery modules during mechanical abuse is complex due to interactions between structural deformation, heat transfer, electrochemical processes, exothermic reactions and mechanical fracture. While prior studies have focused on modeling cell-level behavior, extending these models to module or pack level is necessary for a system level understating of electric vehicle safety. This study develops coupled large deformation finite element models that simultaneously solve for electrochemistry, material failure, internal short circuit and thermal runaway propagation. The models account for mechanical and thermal interactions between lithium-ion cells and other battery components while the contact interfaces are evolving with time. Model-predicted voltage, temperature and force responses are compared with experimental data for validation. The results demonstrate that the approach captures key failure mechanisms, including thermal propagation through heat transfer, electrical propagation from short circuits in parallel-connected cells, and mechanical propagation via penetration and crack formation. These findings show that computational models are valuable tools for understanding battery module failure and providing insight that can reduce the need for extensive experimental testing.

25 ENERGY STORAGE↗
DOE OSTI2021

The homogenous alternative to biomineralization: Zn- and Mn-rich materials enable sharp organismal “tools” that reduce force requirements

We measured hardness, modulus of elasticity, and, for the first time, loss tangent, energy of fracture, abrasion resistance, and impact resistance of zinc- and manganese-enriched materials from fangs, stings and other “tools” of an ant, spider, scorpion and nereid worm. The mechanical properties of the Zn- and Mn-materials tended to cluster together between plain and biomineralized “tool” materials, with the hardness reaching, and most abrasion resistance values exceeding, those of calcified salmon teeth and crab claws. Atom probe tomography indicated that Zn was distributed homogeneously on a nanometer scale and likely bound as individual atoms to more than ¼ of the protein residues in ant mandibular teeth. This homogeneity appears to enable sharper, more precisely sculpted “tools” than materials with biomineral inclusions do, and also eliminates interfaces with the inclusions that could be susceptible to fracture. Based on contact mechanics and simplified models, we hypothesize that, relative to plain materials, the higher elastic modulus, hardness and abrasion resistance minimize temporary or permanent tool blunting, resulting in a roughly 2/3 reduction in the force, energy, and muscle mass required to initiate puncture of stiff materials, and even greater force reductions when the cumulative effects of abrasion are considered. We suggest that the sharpness-related force reductions lead to significant energy savings, and can also enable organisms, especially smaller ones, to puncture, cut, and grasp objects that would not be accessible with plain or biomineralized “tools”.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Characterization of electromigration-induced short-range stress development in Al(0.25 at. % Cu) conductor line

Scanning x-ray microbeam topography and fluorescence experiments were conducted in situ to study the electromigration behavior of a 0.5 μm thick, 10 μm wide, and 200 μm long Al(0.25 at. % Cu) conductor line with 1.5 μm-thick SiO2 passivation on a single crystal Si substrate. The strain sensitivity of x-ray topography measurement allowed detailed examination of the electromigration-induced stress distribution and evolution in the conductor line in response to the depletion of Cu solute early in the electromigration process. Upon electromigration at 0.4 MA/cm2 and 303 °C, a short-range stress gradient was quickly induced by Al migration in the Cu-depleted cathode region to counteract further Al flow. The stress gradient was fully developed during the 5.3 h incubation time, extending over the critical Blech length of about 66 μm from the cathode end. Plastic deformation then occurred at the downstream end of the Cu-depleted region. The preferential electromigration of Cu did not cause detectable stress change outside the Cu-depleted region, except for the significant stress development from the Al2Cu precipitation at the anode end which appeared to initiate the fracture in the passivation. Preliminary finite difference modeling was undertaken to simulate the experimental observations, from which important parameters dictating electromigration in Al(Cu) line were extracted: an apparent effective valence of −5.6 and −1.9 for Cu and Al in Al(Cu), respectively, and a critical Cu concentration of 0.16 at. % above which Al grain boundary diffusion is effectively blocked.

Physics↗
DOE OSTI2022

Stress Intensity Thresholds for Development of Reliable Brittle Materials

Brittle material failure in high consequence systems can appear random and unpredictable at subcritical stresses. Gaps in our understanding of how structural flaws and environmental factors (humidity, temperature) impact fracture propagation need to be addressed to circumvent this issue. A combined experimental and computational approach composed of molecular dynamics (MD) simulations, numerical modeling, and atomic force microscopy (AFM) has been undertaken to identify mechanisms of slow crack growth in silicate glasses. AFM characterization of crack growth as slow as 10 -13 m/s was observed, with some stepwise crack growth. MD simulations have identified the critical role of inelastic relaxation in crack propagation, including evolution of the structure during relaxation. A numerical model for the existence of a stress intensity threshold, a stress intensity below which a fracture will not propagate, was developed. This transferrable model for predicting slow crack growth is being incorporated into mission-based programs.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Geochemical Modeling of Celestite (SrSO 4 ) Precipitation and Reactive Transport in Shales

Celestite (SrSO 4 ) precipitation is a prevalent example of secondary sulfate mineral scaling issues in hydraulic fracturing systems, particularly in basins where large concentrations of naturally occurring strontium are present. Herein, we present a validated and flexible geochemical model capable of predicting celestite formation under such unconventional environments. Simulations were built using CrunchFlow and guided by experimental data derived from batch reactors. These data allowed the constraint of key kinetic and thermodynamic parameters for celestite precipitation under relevant synthetic hydraulic fracturing fluid conditions. Effects of ionic strength, saturation index, and the presence of additives were considered in the combined experimental and modeling construction. This geochemical model was then expanded into a more complex system where interactions between hydraulic fracturing fluids and shale rocks were allowed to occur subject to diffusive transport. We find that the carbonate content of a given shale and the presence of persulfate breaker in the system strongly impact the location and extent of celestite formation. The results of this study provide a novel multicomponent reactive transport model that may be used to guide future experimental design in the pursuit of celestite and other sulfate mineral scale mitigation under extreme conditions typical of hydraulic fracturing in shale formations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Scale bridging damage model for quasi-brittle metals informed with crack evolution statistics

Computationally efficient methods for bridging length scales, from highly resolved micro/meso-scale models that can explicitly model crack growth, to macro-scale continuum models that are more suitable for modeling large parts, have been of interest to researchers for decades. In this work, an improved brittle damage model is presented for the simulation of dynamic fracture in continuum scale quasi-brittle metal components. Crack evolution statistics, including the number, length, and orientation of individual cracks, are extracted from high-fidelity, finite discrete element method (FDEM) simulations and used to generate effective material moduli that reflect the material’s damaged state over time. This strategy allows for the retention of small-scale physical behaviors such as crack growth and coalescence in continuum scale hydrodynamic simulations. However, the high-fidelity simulations required to generate the crack statistics are computationally expensive. Thus, steps were taken to produce a flexible constitutive model to reduce the number of costly high-fidelity simulations needed to produce accurate results. A new stress based degradation criterion is introduced for the degradation of individual material zones. This allows for the development of a heterogeneous damage distribution within the bulk material. Then a flow stress model is added to the hydrodynamic simulation to account for plasticity in quasi-brittle materials. As a result, the effective moduli model can be applied to a larger range of materials. The effective moduli constitutive model is used to simulate beryllium flyer plate experiments. The results from the continuum scale simulations using statistics from a single high-fidelity simulation are found to be in excellent agreement with numerical and experimental velocity interferometer data. The same set of crack statistics are used to extrapolate the results of a higher rate flyer plate case using the effective moduli model. In conclusion, the extension of this model to higher rate cases shows promise for further reducing the number of costly high-fidelity simulations needed to generate crack statistics.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Multilevel Graph Partitioning for Three-Dimensional Discrete Fracture Network Flow Simulations

We present a topology-based method for mesh-partitioning in three-dimensional discrete fracture network (DFN) simulations that takes advantage of the intrinsic multi-level nature of a DFN. DFN models are used to simulate flow and transport through low-permeability fractured media in the subsurface by explicitly representing fractures as discrete entities. The governing equations for flow and transport are numerically integrated on computational meshes generated on the interconnected fracture networks. Modern high-fidelity DFN simulations require high-performance computing on multiple processors where performance and scalability depends partially on obtaining a high-quality partition of the mesh to balance work-loads and minimize communication across all processors. The discrete structure of a DFN naturally lends itself to various graph representations, which can be thought of as coarse-scale representations of the computational mesh. Using this concept, we develop two applications of the multilevel graph partitioning algorithm to partition the mesh of a DFN. In the first, we project a partition of the graph based on the DFN topology onto the mesh of the DFN and in the second, this DFN-based projection is used as the initial condition for further partitioning refinement of the mesh. We compare the performance of these methods with standard multi-level graph partitioning using graph-based metrics (cut, imbalance, partitioning time), computational-based metrics (FLOPS, iterations, solver time), and total run time. The DFN-based and the mesh-based partitioning methods are comparable in terms of the graph-based metrics, but the time required to obtain the partition is several orders of magnitude faster using the DFN-based partitions. The computation-based metrics show comparable performance between both methods so, in combination, the DFN-based partitions are several orders of magnitude faster than the mesh-based partition. Furthermore, the method which uses the DFN-partition solution as the initial condition of the mesh partition provided cut and imbalance values that were close to the mesh-based partition but in a fraction of the time. In turn, this hybrid method outperformed both of the other methods in terms of the total run time.

58 GEOSCIENCES↗
DOE OSTI2020

NNSA-IAEC Science Area V Environmental ISR: Waste Management and Subsurface Science (Final Report FY18-FY20)

In August 2017, the U.S. National Nuclear Security Administration (NNSA) and the Israel Atomic Energy Commission (IAEC) launched a joint waste management program to evaluate the feasibility of siting a radioactive waste repository at intermediate depths in Israel. The bilateral collaboration involves three National Laboratories (LLNL, LANL, and Sandia), and the Nuclear Research Center, Negev, Geological Survey of Israel, and Ben Gurion University in Israel. This research is designed to assist with the evaluation of locations for an intermediate nuclear waste repository in Israel. This final report is a high-level summary of the results from the past three years by the team that worked on the WM-2: Radionuclide colloid-facilitated transport in fractured carbonate rock project for Science Area V, Environmental ISR: Subsurface Science and Waste Management. We have completed all of our tasks and met all of our milestones. MW-2 had three integrated tasks. Task 1 was laboratory experiments investigating colloid transport of radionuclides in natural fractured carbonate rocks from the Avdat formation, Israel. Task 2 was field experiments injecting radionuclide analogues pre-sorbed to clay colloids into fractured carbonate rock, and task 3 was numerical modeling of laboratory and field experimental results to assess the importance of colloids in the migration of relevant radionuclides. This task was incorporated into tasks 1 and 2. A 4 th task was initiated in FY20, to investigate the role of organics in facilitating radionuclide transport under the same conditions explored in task 1 due to the high organic content of the rocks being investigated for the immediate borehole locations. Preliminary results will be summarized here and the work will continue in FY21-23. The overall objective of this research is to evaluate the role of colloids (naturally occurring < 1 micron particles) in facilitating the transport of long-lived radionuclides from a nuclear waste repository situated in fractured carbonate rocks. Currently little data exists on radionuclide transport in carbonate rocks, but this is the main rock type available for siting a repository in Israel. Laboratory experiments were carried out in both Israel and the U.S., field experiments have taken place in Israel, and reactive transport modeling involved LLNL, LANL and Israel partner institutions. LLNL hosted a graduate student, Emily Tran, from Ben Gurion University in FY18 and FY19 and much of the work presented here was part of her PhD research.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2024

CONSTRAINT-INDEPENDENT CONSTANT CTOA DETERMINATION FOR DUCTILE STABLE CRACK GROWTH

Crack tip opening angle (CTOA) has been used as a reliable fracture toughness parameter for decades to characterize stable ductile crack growth for thin-walled aerospace structures in the low-constraint conditions. Recently, the CTOA parameter was also applied to the pipeline industry, and a CTOA test standard ASTM E3039 was thus developed for testing a critical constant CTOA. Research showed that the constant CTOA can reasonably describe fracture toughness required to arrest a dynamic crack propagation for a modern gas pipeline. However, the CTOA fracture criterion requires constraint-independent CTOA toughness against stable ductile crack growth. ASTM E3039 recommends a drop weight tearing test (DWTT) specimen for CTOA testing. Since a shallow crack is used, DWTT measured CTOA may depend on constraint level at the crack tip. To understand if it is the case, this paper evaluates the critical CTOA for a set of fracture toughness tests on single edge notched bend (SENB) specimens with shallow and deep cracks based on four CTOA estimation models. In which, the Ln(P)-LLD linear fit model is similar to that used by ASTM E3039 in the CTOA calculation. Fracture test data for X80 pipeline steel and HY80 structural steel are considered in the CTOA evaluation. The results show that the four CTOA models can determine a crack size-independent constant CTOA over stable ductile crack growth for the SENB specimens. As a result, CTOA determined by ASTM E3039 is constraint&#x2;independent and transferable to use for an actual crack propagating in a gas pipeline.

Zhu, Xian-Kui↗
DOE OSTI2025

Prediction and Analysis of Utah FORGE Injection Activities using a Coupled Thermo-hydro-mechanical and Earthquake (THM+E) Modeling Workflow

A coupled thermo-hydro-mechanical (THM) numerical workflow that is capable of modeling seismic slip is critical for the successful development of enhanced geothermal systems (EGS). By integrating key physical processes, this workflow enables accurate simulation of temperature and pressure diffusions, stress changes, and induced seismicity. As a result, it serves as a vital tool for predicting induced seismicity and optimizing reservoir stimulation strategies. The Utah FORGE (Frontier Observatory for Research in Geothermal Energy) project, located near Milford, Utah, is a U.S. Department of Energy initiative aimed at advancing EGS technology. In April 2024, eight new stimulation stages (Stages 3R-10) were conducted in well 16A (injection well) subsequent to the first series of stimulation (Stages 1-3) performed in April, 2022. To monitor the induced seismicity, geophones were deployed in wells 58-32, 56-32, and 78B-32, while fiber optic cables were also installed in wells 16B, 78-32, and 78B-32 to collect microseismic data and detect frac hits Preliminary analyses of microseismic catalogs and fiber optic data suggest that the stimulated fractures in Stages 3R–6 closely align with that generated during Stage 3, indicating that the new stimulations were likely reactivating the previously stimulated fracture. To better understand the underlying process, a comprehensive modeling approach that can accurately capture thermal, hydrological, mechanical, and seismic responses is essential. In this work, we propose and utilize a coupled thermo-hydro-mechanical and earthquake (THM+E) simulation workflow to numerically investigate the stimulation activities on well 16A. The specific objective is to confirm whether the new stimulation stages (Stages 3R–6) reactivated fractures previously stimulated during Stage 3. For this purpose, we perform THM+E simulations individually for Stages 3, 3R, 4, and 5, incorporating the discrete fracture networks (DFNs) created by the plane-fitting technique based on the microseismic catalogs. The simulation workflow consists of two separate models: a THM model and an earthquake model, coupled in a one-way manner. Detailed descriptions of the workflow are provided in Section 3. Simulation results are presented in terms of injection pressure, permeability evolution, and predicted seismic catalogs, which are then compared with field data for further analyses. This report is structured as follows. In Section 2, we present detailed analyses of the field data and propose the hypothesis that the new stimulation stages (Stages 3R–6) were probably reactivating the previously stimulated fractures in Stage 3. In Section 3, we introduce the coupled THM+E workflow and the problem setup to validate our hypothesis, followed by the simulation results for each stage in Section 4. Meanwhile, discussions are included to analyze the model predictions and their comparison with field data. Lastly, we conclude the report and outline future plans in Section 5.

15 GEOTHERMAL ENERGY↗
DOE OSTI2020

Tailoring ductile-phase toughened tungsten hierarchical microstructures for plasma-facing materials

Biological materials, such as bones, nacre, etc., have been found to show attractive and unique combinations of stiffness, strength and fracture toughness. Research has been conducted in various engineering areas to mimic the naturally hierarchical microstructures or nanostructures of these materials to produce materials with optimum stiffness, high strength and fracture toughness for their intended structural applications. In the same objective, this work applies a multiscale microstructural approach recently developed to investigate the deformation and fracture behavior of ductile phase toughened tungsten (W) materials such as tungsten-nickel/iron (W-Ni/Fe) composites that possess lamellar-like and hierarchical “brick-and-mortar” (BAM) microstructures. First, the approach is used to simulate tensile loading of W-Ni/Fe specimens cut out from hot-rolled W-Ni/Fe plates which possess lamellar-like microstructures. The finite element model of the gage section of the specimen consists of a W-Ni/Fe dual-phase microstructural domain in which the constitutive behaviors of W and Ni/Fe phases are described by an elastic-plastic damage model. The predicted material stress-strain response and crack pattern development as a function of loading are then compared to the corresponding experimental results to determine the constitutive model parameters for W and Ni/Fe. Subsequently, the model parameters are used to analyze W-Ni/Fe BAM microstructures that are artificially created to investigate the effects of microstructural features and morphologies on the composite stress-strain response, damage, and fracture patterns. The modeling shows that regular BAM microstructures that experience bridging mechanism combined with crack penetration across W-phase regions exhibit significantly higher strengths than the rather random lamellar-like microstructures. More significantly, however adjusting the brick’s length-to-height ratio, the BAM microstructure can be designed to allow a more distributed damage zone which leads to increased strength, ductility and fracture energy.

36 MATERIALS SCIENCE↗
DOE OSTI2026

One Earth Energy Seismic Interpretation

The objectives of the Illinois Storage Corridor (ISC) project are to accelerate commercial deployment of carbon capture utilization and storage at two individual sites and receive approvals for Underground Injection Control (UIC) Class VI permits for construction at each site (ISC Project Narrative, 2020). As part of this project, and as part of the subsurface geologic characterization, 2D seismic data was acquired at both sites. This report summarizes the findings from the 2D and 3D seismic interpretation at the One Earth Energy site near Gibson City, Illinois. The seismic data confirms the stratigraphic continuity of the Mt. Simon Arkose Zone storage interval and the Eau Claire confining unit across the project area. The seismic data also indicates that there are faults that transect the Mt. Simon Arkose Zone Sandstone storage reservoir within the modeled CO 2 plume (for more detailed information, see Faults and Fractures section of One Earth Energy Class VI Permit applications). However, the seismic data also shows that there are no faults within the modeled CO 2 plume that transect the confining unit Eau Claire Formation. The faults that transect the Mt. Simon Arkose Zone Sandstone storage reservoir all tip out in the Lower Mt. Simon Formation and do not reach the overlying Eau Claire confining unit. A small 3D survey acquired around the One Earth Energy #1 characterization well confirms these findings.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2020

Predicting the potential for mineral scale precipitation in unconventional reservoirs due to fluid-rock and fluid mixing geochemical reactions

Mineral precipitation within hydraulically fractured shale may affect fluid flow pathways and impact longterm hydrocarbon production. The ability to predict geochemical reactions that lead to problematic mineral precipitation will lead to active reservoir management strategies to improve overall production. Using the Marcellus Shale as a case study, a combination of laboratory experiments and reaction path modeling was applied to determine which reactions are likely to occur upon introduction of hydraulic fracturing fluid into the shale reservoir. Experimental results indicate that contact between fracturing fluid and shale will result in dissolution of primary minerals (quartz, feldspars, kaolinite, chlorite, pyrite) and secondary mineral precipitation over time periods of less than one week. Precipitation of barite, Fe-oxides, feldspars, amorphous silica and clay is likely to occur within the reservoir during shut in and early flowback due to mixing between fracturing fluid and reservoir brine as based on modeling saturation indices using experimental fluid data. Reaction path modeling of the experimental scenarios corroborates the specific dissolution and precipitation reactions observed experimentally. Comparison of the results to injected and produced waters from a Marcellus Shale well pad in Greene County, PA, USA, shows that the mineral reactions occur during the hydraulic fracturing, shut in, and early flowback periods. The results presented in this paper demonstrate the value in applying experimental approaches to identify mineral precipitation/dissolution reactions that may significantly impact reservoir performance. The good agreement between geochemical models and experimental results provides confidence that numerical models can be applied to screen the potential fluid-mineral and fluid-mixing reactions in unconventional reservoirs that result in undesired mineral scale precipitation.

58 GEOSCIENCES↗
DOE OSTI2021

Fatigue-Like Behavior of Silver Metallization Gridlines and Proposed Damage Mechanics Model

In this work, we present a novel method to evaluate the fracture and fatigue-like behavior of crystalline silicon PV cell metallization. Initial results demonstrate that silver metal gridlines are capable of bridging a ~10–20 µm gap prior to becoming electrically open and do demonstrate fatigue-like behavior under cyclic straining conditions. This later phenomenon was well fit with a traditional power law equation for high cycle fatigue and enabled the development of a damage mechanics model for gridline failure.

41 EE - Solar Energy Technologies Office (EE-4S)↗