Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fracture simulation”

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 865 records · Page 48

Time-temperature-stress capabilities of composite materials for advanced supersonic technology application, phase 1

Implementation of metal and resin matrix composites into supersonic vehicle usage is contingent upon accelerating the demonstration of service capacity and design technology. Because of the added material complexity and lack of extensive service data, laboratory replication of the flight service will provide the most rapid method of documenting the airworthiness of advanced composite systems. A program in progress to determine the time temperature stress capabilities of several high temperature composite materials includes thermal aging, environmental aging, fatigue, creep, fracture, and tensile tests as well as real time flight simulation exposure. The program has two parts. The first includes all the material property determinations and aging and simulation exposures up through 10,000 hours. The second continues these tests up to 50,000 cumulative hours. Results are presented of the 10,000 hour phase, which has now been completed.

Kerr, J. R.↗

An elastic-plastic finite element analysis of crack initiation, stable crack growth, and instability

Studies have been conducted to develop efficient techniques to simulate crack extension and to examine various local and global fracture criteria. Of the considered criteria, the crack-tip-opening angle (CTOA) or displacement (CTOD) at a specified distance from the crack tip was shown to be most suited for modeling stable crack growth and instability during the fracture process. The results obtained in a number of studies show the necessity for studying different crack configurations when assessing the validity of any fracture criteria. One of the objectives of the present investigation is related to a critical evaluation of the CTOD growth criterion using an elastic-plastic finite element analysis under monotonic loading to failure. The analysis was found to predict three stages of crack growth behavior under monotonic loading to failure. Calculated CTOD values agreed well with experimental values for crack growth initiation.

Newman, J. C., Jr.↗

Progressive delamination in polymer matrix composite laminates: A new approach

A new approach independent of stress intensity factors and fracture toughness parameters has been developed and is described for the computational simulation of progressive delamination in polymer matrix composite laminates. The damage stages are quantified based on physics via composite mechanics while the degradation of the laminate behavior is quantified via the finite element method. The approach accounts for all types of composite behavior, laminate configuration, load conditions, and delamination processes starting from damage initiation, to unstable propagation, and to laminate fracture. Results of laminate fracture in composite beams, panels, plates, and shells are presented to demonstrate the effectiveness and versatility of this new approach.

Chamis, C. C.↗

Residual Strength Predictions with Crack Buckling

Fracture tests were conducted on middle crack tension, M(T), and compact tension, C(T), specimens of varying widths, constructed from 0.063 inch thick sheets of 2024-T3 aluminum alloy. Guide plates were used to restrict out-of-plane displacements in about half of the tests. Analyses using the three-dimensional, elastic-plastic finite element code WARP3D simulated the tests with and without guide plates using a critical CTOA fracture criterion. The experimental results indicate that crack buckling reduced the failure loads by up to 40%. Using a critical CTOA value of 5.5 deg., the WARP3D analyses predicted the failure loads for the tests with guide plates within +/- 10% of the experimentally measured values. For the M(T) tests without guide plates, the WARP3D analyses predicted the failure loads for the 12 and 24 inch tests within 10%, while over predicting the failure loads for the 40 inch wide tests by about 20%.

Dawicke, D. S.↗

hifem_transport

SAND2022-14657 O The hifem-transport software is used to simulate transport processes in large geologic environments characterized by complex fractures and man-made infrastructure at an affordable computational cost. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

Beskardes, Gungor↗

Finite Element Analysis (FEA) for Water-Foam Fracturing of Granite Rock

In addition to the foam data that were obtained from literature and that were collected from the current study, simulation data was also generated from finite element analysis (FEA) conducted in this study using COMSOL Multiphysics software. The FEA models were built to simulate the experiments conducted at Oak Ridge National Laboratory (ORNL) on cement and granite samples. In these FEA models, temperature was kept at ambient while the pressure profile resembled the loading conditions during the ORNL experiments, where pressure was either monotonically increased or applied cyclically. The cement material was used as a model material and was used to study Von Mises stress and tensile stress distribution for different bore hole length geometry using a parametric sweep with water as fracturing fluid using solid-fluid interaction module. For the granite material, FEA models were developed for stress analysis of cylindrical samples with water or foam fluids. The solid mechanics module in COMSOL was implemented to solve for Von Mises stress and tensile stress. The fluid-structure interaction module was implemented to solve for water-foam interaction on granite cylinder with addition of fluid-loading on structure, i.e., large deformation in solid mechanics with no impact on fluid deformation. Foam was considered as a pseudo single-phase compressible fluid for which material properties were calculated from water and gas (nitrogen) phases. The density of foam is calculated as a function of the densities of water and nitrogen, while viscosity is a function of temperature. Four types of FEA analyses were modelled: 1. Monotonic injection with water 2. Monotonic injection with foam 3. Cyclic injection with water 4. Cyclic injection with foam All the COMSOL files are converted to a zip file which is save in .mph.

15 GEOTHERMAL ENERGY↗

Numerical Simulation Applications in the Design of EGS Collab Experiment 1

The United States Department of Energy, Geothermal Technologies Office (GTO) is funding a collaborative investigation of enhanced geothermal systems (EGS) processes at the meso-scale. This study, referred to as the EGS Collab project, is a unique opportunity for scientists and engineers to investigate the creation of fracture networks and circulation of fluids across those networks under in-situ stress conditions. The EGS Collab project is envisioned to comprise three experiments and the site for the first experiment is on the 4850 Level in phyllite of the Precambrian Poorman formation, at the Sanford Underground Research Facility, located at the former Homestake Gold Mine, in Lead, South Dakota. Principal objectives of the project are to develop a number of intermediate-scale field sites and to conduct well-controlled in situ experiments focused on rock fracture behavior and permeability enhancement. Data generated during these experiments will be compared against predictions of a suite of computer codes specifically designed to solve problems involving coupled thermal, hydrological, geomechanical, and geochemical processes. Comparisons between experimental and numerical simulation results will provide code developers with direction for improvements and verification of process models, build confidence in the suite of available numerical tools, and ultimately identify critical future development needs for the geothermal modeling community. Moreover, conducting thorough comparisons of models, modelling approaches, measurement approaches and measured data, via the EGS Collab project, will serve to identify techniques that are most likely to succeed at the Frontier Observatory for Research in Geothermal Energy (FORGE), the GTO’s flagship EGS research effort. As noted, outcomes from the EGS Collab project experiments will serve as benchmarks for computer code verification, but numerical simulation additionally plays an essential role in designing these meso-scale experiments. This paper describes specific numerical simulations supporting the design of Experiment 1, a field test involving hydraulic stimulation of two fractures from notched sections of the injection borehole and fluid circulation between sub-horizontal injection and production boreholes in each fracture individually and collectively, including the circulation of chilled water. Whereas the mine drift allows for accurate and close placement of monitoring instrumentation to the developed fractures, active ventilation in the drift cooled the rock mass within the experimental volume. Numerical simulations were executed to predict seismic events and magnitudes during stimulation, initial fracture orientations for smooth horizontal wellbores, pressure requirements for fracture initiation from notched wellbores, fracture propagation during stimulation between the injection and production boreholes, tracer travel times between the injection and production boreholes, and produced fluid temperatures with chilled water injections, and pressure limits on fluid circulation to avoid fracture growth.

Numerical Simulation, EGS Collab, EGS, Sanford Und↗

Analysis of the Effects of Surface Pitting and Wear on the Vibrations of a Gear Transmission System

A comprehensive procedure to simulate and analyze the vibrations in a gear transmission system with surface pitting, 'wear' and partial tooth fracture of the gear teeth is presented. An analytical model was developed where the effects of surface pitting and wear of the gear tooth were simulated by phase and magnitude changes in the gear mesh stiffness. Changes in the gear mesh stiffness were incorporated into each gear-shaft model during the global dynamic simulation of the system. The overall dynamics of the system were evaluated by solving for the transient dynamics of each shaft system simultaneously with the vibration of the gearbox structure. In order to reduce the number of degrees-of-freedom in the system, a modal synthesis procedure was used in the global transient dynamic analysis of the overall transmission system. An FFT procedure was used to transform the averaged time signal into the frequency domain for signature analysis. In addition, the Wigner-Ville distribution was also introduced to examine the gear vibration in the joint time frequency domain for vibration pattern recognition. Experimental results obtained from a gear fatigue test rig at NASA Lewis Research Center were used to evaluate the analytical model.

Choy, F. K.↗

An integrated numerical model for coupled poro-hydro-mechanics and fracture propagation using embedded meshes

Integrated models for fluid-driven fracture propagation and general multiphase flow in porous media are valuable to the study and engineering of several systems, including hydraulic fracturing, underground disposal of waste, and geohazard mitigation across such applications. This work extends the coupled model multiphase flow and poromechanical model of Ren et al. (2018) to admit fracture propagation (FP). The coupled XFEM-EDFM scheme utilizes a separate fracture mesh that is embedded on a static background mesh. The onset and dynamics of fracture propagation are governed by the equivalent stress intensity factor (SIF) criterion. A domain-integral method (J integral) is applied to compute this information. Additionally, an adaptive time-marching scheme is proposed to rapidly restrict and grow temporal resolution to match the underlying time-scales. The proposed model is verified with analytical solutions, and shows the capability to accurately and adaptively co-simulate fluid transport and deformation as well as the propagation of multiple fractures.

42 ENGINEERING↗

Numerical Simulation of Delamination Growth in Composite Materials

The use of decohesion elements for the simulation of delamination in composite materials is reviewed. The test methods available to measure the interfacial fracture toughness used in the formulation of decohesion elements are described initially. After a brief presentation of the virtual crack closure technique, the technique most widely used to simulate delamination growth, the formulation of interfacial decohesion elements is described. Problems related with decohesion element constitutive equations, mixed-mode crack growth, element numerical integration and solution procedures are discussed. Based on these investigations, it is concluded that the use of interfacial decohesion elements is a promising technique that avoids the need for a pre-existing crack and pre-defined crack paths, and that these elements can be used to simulate both delamination onset and growth.

Camanho, P. P.↗

VCCT with Progressive Nodal Release for Simulating Mixed-Mode Delamination: Formulation, Algorithmic Improvements and Implications

To simulate the propagation of cracks using the Virtual Crack Closure Technique (VCCT), intermediate crack positions between existing node pairs can be modeled by progressively releasing the nodes using fracture mechanics principles. This approach has been implemented in commercial software and can be applied to simulate delamination growth. Recent research has provided insight into the synchronized mixed mode nodal release, extending and formalizing the methodology for delamination growth. This paper documents the development of the Jacobians associated with delamination propagation under quasi-static mixed-mode conditions, leading to improvements in convergence of the implicit nonlinear solvers used.

Mabson, G. E.↗

VCCT with Progressive Nodal Release for Simulating Mixed-Mode Delamination: Formulation, Algorithmic Improvements and Implications

To simulate the propagation of cracks using the Virtual Crack Closure Technique (VCCT), intermediate crack positions between existing node pairs can be modeled by progressively releasing the nodes using fracture mechanics principles. This approach has been implemented in commercial software and can be applied to simulate delamination growth. Recent research has provided insight into the synchronized mixed mode nodal release, extending and formalizing the methodology for delamination growth. This paper documents the development of the Jacobians associated with delamination propagation under quasi-static mixed-mode conditions, leading to improvements in convergence of the implicit nonlinear solvers used.

G E Mabson↗

Modeling shockwaves and impact phenomena with Eulerian peridynamics

Most previous development of the peridynamic theory has assumed a Lagrangian formulation, in which the material model refers to an undeformed reference configuration. Here, an Eulerian form of material modeling is developed, in which bond forces depend only on the positions of material points in the deformed configuration. The formulation is consistent with the thermodynamic form of the peridynamic model and is derivable from a suitable expression for the free energy of a material. We show that the resulting formulation of peridynamic material models can be used to simulate strong shock waves and fluid response in which very large deformations make the Lagrangian form unsuitable. The Eulerian capability is demonstrated in numerical simulations of ejecta from a wavy free surface on a metal subjected to strong shock wave loading. The Eulerian and Lagrangian contributions to bond force can be combined in a single material model, allowing strength and fracture under tensile or shear loading to be modeled consistently with high compressive stresses. Furthermore, we demonstrate this capability in numerical simulation of bird strike against an aircraft, in which both tensile fracture and high pressure response are important.

36 MATERIALS SCIENCE↗

Coupled Flow and Deformation Modeling of Carbon Dioxide Migration in the Presence of a Caprock Fracture during Injection

Understanding the transport of carbon dioxide (CO{sub 2}) during long-term CO{sub 2} injection into a typical geologic reservoir, such as a saline aquifer, could be complicated because of changes in geochemical, hydrogeological, and hydromechanical behavior. While the caprock layer overlying the target aquifer is intended to provide a tight, impermeable seal in securing injected CO{sub 2}, the presence of geologic uncertainties, such as a caprock fracture or fault, may provide a channel for CO{sub 2} leakage. There could also be a possibility of the activation of a new or existing dormant fault or fracture, which could act as a leakage pathway. Such a leakage event during CO{sub 2} injection may lead to a different pressure and ground response over a period of time. In the present study, multiphase fluid flow simulations in porous media coupled with geomechanics were used to investigate the overburden geologic response and plume behavior during CO{sub 2} injection in the presence of a hypothetical permeable fractured zone in a caprock, existing or activated. Both single-phase and multiphase fluid flow simulations were performed. The CO{sub 2} migration through an existing fractured zone leads to changes in the fluid pressure in the overburden geologic layers and could have a significant impact on ground deformation behavior. Results of the study show that pressure signatures and displacement patterns are significantly different in the presence of a fractured zone in the caprock layer. The variation in pressure and displacement signatures because of the presence of a fractured zone in the caprock at different locations may be useful in identifying the presence of a fault/fractured zone in the caprock. The pressure signatures can also serve as a mechanism to identify the activation of leakage pathways through the caprock during CO{sub 2} injection. Pressure response and ground deformation behavior from sequestration modeling could be useful in the development of smart technologies to monitor safe CO{sub 2} storage and understand CO{sub 2} transport, with limited field instrumentation.

54 ENVIRONMENTAL SCIENCES↗

Numerical estimations of storage efficiency for the prospective CO 2 storage resource of shales

Here, hydrocarbon-bearing shale formations might be an attractive geologic reservoir for permanent carbon dioxide (CO 2 ) storage. Shale formations applicable for storage require previous hydrocarbon production to deplete a reservoir, hydraulic fracturing to provide a highly-permeable stimulated zone, sufficient depths (generally > 800 m or 2600 feet) to maintain CO 2 in a supercritical state, and an overlying seal to prevent CO 2 migration to underground sources of drinking water and into the atmosphere. CO 2 is stored in shale as a free phase within fractures and matrix pores, and as a sorbed phase on organic matter and clay minerals. Recently in our previous work, we presented a screening-level assessment methodology for CO 2 storage in shales using a volumetric approach. The approach deals with reduction of CO 2 storage through estimations of efficiency factors based on petrophysical properties (i.e., bulk volume, porosity, sorption, etc.) and their limitations on fluid transport. Here, we conducted numerical simulations using the FRACGEN/NFFLOW simulator to study the CO 2 injection into a depleted hydro-fractured shale reservoir to estimate storage efficiencies using a range of reservoir parameters and injection scenarios. Specifically, the ranges for two efficiency factors, E Φ and E S , measure the effectiveness of free and sorbed CO 2 storage. These efficiency factors were estimated to have P 10 to P 90 probability ranges of 0.15 to 0.36 for and 0.11 to 0.24 for E S , reported after 60 years of CO 2 injection.

58 GEOSCIENCES↗

Composition-dependent fracture energy in metallic glasses

The interplay between metallic glasses (MGs) mechanical properties, fracture energy (G), and glass-forming ability (GFA) and their dependence on alloy composition remain poorly understood. Here, we perform molecular dynamics simulations to investigate the intrinsic composition-dependence of G in Cu x Zr 100-x MGs (x = 20, 30, 40, 50, 64). The results indicate that the value of G increases with Cu content. In addition, it is revealed that MGs with higher G values display higher Poisson’s ratio (ν) and GFA, suggesting a close correlation between fracture toughness, mechanical properties, and GFA. This correlation between G, ν and GFA can be understood based on the fragility (m) of supercooled liquids, which is directly related to the structural heterogeneity in MGs. Larger m values are related to dynamic slowdown and supercooled liquid stabilization, which enhance GFA and the formation of pronounced structural heterogeneity, comprised of loosely packed regions that favor β-relaxation and the activation of shear transformation zones. Those concurrently promote the expansion of the plastic zone at the crack tip, enhancing the observed value of G. Furthermore, these simulation results shed light on the intrinsic relationship between fracture toughness, mechanical properties, and alloy composition in MGs.

36 MATERIALS SCIENCE↗

Deformation Mechanisms of Vertically Stacked WS 2 /MoS 2 Heterostructures: The Role of Interfaces

Here, the mechanical and optical properties generated due to the stacking of different atomically thin materials have made it possible to tune and engineer these materials for next-generation electronics. The understanding of the interlayer interactions in such stacked structures is of fundamental interest for structure and property correlation. Here, a combined approach of in situ Raman spectroscopy and mechanical straining along with molecular dynamics (MD) simulations has been used to probe one such interface, namely, the WS 2 /MoS 2 heterostructure. Vertical heterostructures on poly(methyl methacrylate), when flexed, showed signs of decoupling at 1.2% strain. Theoretical calculations showed strain-induced stacking changes at 1.75% strain. The sliding characteristics of layers were also investigated using scanning probe microscopy based nanoscratch testing, and the results are further supported by MD simulations. The present study could be used to design future optoelectronic devices based on WS 2 /MoS 2 heterostructures.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Findings on Subtask 3.1 - Bakken Rich Gas Enhanced Oil Recovery Project

Total in-place oil for the Bakken petroleum system (BPS) (which includes the Bakken and Three Forks Formations) has been estimated to be 600 billion barrels (bbl). However, BPS wells have decline rates as high as 85% over the first 3 years of their lives, and primary recovery factors typically range from 3% to 10% of original oil in place. Given the low initial recovery rates, even small incremental productivity improvements could dramatically increase technically recoverable oil in the BPS. One potential solution is enhanced oil recovery (EOR) using gas injection, such as carbon dioxide (CO2) or hydrocarbon (HC) gases. While commonly used in conventional reservoirs, CO2 EOR in unconventional tight oil reservoirs has been limited to pilot tests. EOR using rich gas (mixture of methane, ethane, and propane) has also been employed in numerous pilots in several unconventional plays and has recently been successfully applied in the Eagle Ford play. If successful, large-scale gas-based EOR in the BPS could dramatically increase oil productivity and recovery factors and extend the life of the play for decades. While CO2 may be a technically suitable working fluid for EOR in the BPS, supplies are limited and costs for using CO2 in EOR pilots are prohibitively high. Meanwhile, produced gas flaring has presented challenges for BPS operators in North Dakota. Analysis conducted by the North Dakota Pipeline Authority indicates that the current gas-gathering infrastructure in North Dakota is insufficient to accommodate all of the associated gas that is produced from the BPS. The geographically isolated location of North Dakota relative to large natural gas markets, combined with suppressed natural gas prices, has made it economically challenging for industry to invest capital in expanding gas-gathering infrastructure in the state. These circumstances led to a research program conducted by the Energy & Environmental Research Center (EERC) in partnership with Liberty Resources Management Company LLC (LR) to examine the potential to use rich gas injection for EOR and mitigate flaring. A rich gas EOR pilot test was designed and executed by LR at its Stomping Horse development area in Williams County, North Dakota. From July 2018 through May 2019, a total of 160 million standard cubic feet (MMscf) of rich produced gas was injected into the BPS using five different wells in a sequential injection strategy. LR’s Leon–Gohrick drill spacing unit (DSU) was used as the test site. Regulatory oversight was provided by the North Dakota Industrial Commission (NDIC). Technical support was provided by the EERC through a series of laboratory, modeling, and field-based activities, and additional post-pilot research activities incorporated learnings from the test, developed new laboratory data, improved fracture modeling methods, and developed machine learning and big data analytics. The results from the Stomping Horse rich gas EOR pilot activities indicate that developing an effective, economical EOR approach for the BPS will require more field tests. Another key lesson learned from the Stomping Horse tests is that detailed pre- and posttest data on reservoir conditions and fluids production are essential. Robust reservoir characterization provides information that is crucial to creating realistic geomodels and conducting valid dynamic simulations of potential EOR scenarios. A detailed understanding of the completions and production history of offset wells is also necessary for valid test result interpretations. This knowledge is essential to designing the operational parameters of injectivity tests and interpreting the results. A conformance control strategy is also essential to success. Laboratory-based examinations of rich gas interactions with reservoir fluids and rocks were conducted, with an emphasis on determining the ability to mobilize oil in the tight reservoir rocks and shales of the BPS. Injection fluid composition was shown to have a positive impact on reducing reservoir oil minimum miscibility pressure (MMP), reducing interfacial tension (IFT), and altering wettability. IFT and contact angle measurements demonstrated that wettability can be altered in the presence of rich gas, suggesting the potential to improve oil recovery. Iterative modeling of surface infrastructure and reservoir performance using data generated by the various project activities was conducted. A geologic model of the Stomping Horse area was built; history-matched oil, gas, and water production was used in simulations of various EOR scenarios. Early programmatic modeling results were used to support LR’s design and operation of the EOR pilot and to provide insight regarding optimization of future commercial-scale BPS EOR design and operations. Post-pilot modeling focused on alternative methods of understanding complex fracture networks and accelerating simulation time. These led to improved simulation run times and provide excellent history-matching results. Several of these iterative models were used as the bases for developing algorithms into machine learning and big data analytics. History matching in reservoir simulation is time-consuming and computer processing-intensive. Machine learning algorithms were created, and an automated history-matching tool was developed. A large set of synthetic reservoir simulations were created to generate well responses (oil, gas, and water production, well bottomhole pressure [BHP], and tracer or propane breakthrough) for a set of EOR operating parameters that included offset well status (open or closed), injectate (rich gas or propane), injection rate, and injection well BHP. A user interface was developed to provide real-time visualization. Machine learning-based models were developed to provide rapid forecasting of well performance given a set of user-defined EOR operating parameters. These predictive models allow the user to modify the offset well status, injection rate, and injection well BHP and rapidly forecast future production performance. The combination of real-time visualization tools with real-time forecasting tools provides a framework for real-time control—operational changes that the EOR site operator can enact (e.g., changing gas injection rates) to affect the observed performance and potentially improve the EOR outcome. There is great reason to be optimistic about the future of EOR in the Bakken. The results of the laboratory studies suggest significant potential for high rates of oil mobilization using produced field gas injection under the right conditions. The results of the lab studies, combined with rigorous statistical analysis of well production data and associated modeling efforts, confirm the notion that fluid mobility within the reservoir is controlled by fractures. As more knowledge is gained about the nature and distribution of fracture networks in the Bakken, the industry will be in a better position to predict and, ultimately, influence fluid mobility. New field tests are necessary to develop a more complete understanding of those conditions. Thoughtful and creatively engineered field tests within a well-characterized geologic setting will yield the fundamental knowledge needed to take Bakken oil production to the next level. This subtask was cofunded through the EERC–U.S. Department of Energy Joint Program on Research and Development for Fossil Energy-Related Resources Cooperative Agreement No. DE-FE0024233. Nonfederal funding was provided by the North Dakota Industrial Commission’s Oil and Gas Research Program and Computer Modelling Group.

04 OIL SHALES AND TAR SANDS↗