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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Characterizing the Impacts of Multi-Scale Heterogeneity on Solute Transport in Fracture Networks

We model flow and transport in fracture networks with varying degrees of fracture-to-fracture aperture heterogeneity and network intensity to show how changes in these properties can cause the emergence of anomalous flow and transport behavior. If fracture-to-fracture aperture heterogeneity is increased in sparse networks, velocity fluctuations inhibit high flow rates and solute transport can be delayed. Surprisingly, transport can be slowed even in cases where hydraulic aperture is monotonically increased. As the intensity of the networks is increased, more connected pathways allow for particles to bypass these effects. There exists transition behavior where with relatively few connected pathways in a network, first arrival times of particles are not heavily affected by fracture-to-fracture aperture heterogeneity, but the scaling behavior of the tails is strongly influenced. These results reinforce the importance of considering multi-scale effects in fractured systems and can inform flow and transport processes in both natural and engineered fractured systems.

58 GEOSCIENCES↗
DOE OSTI2022

A smeared crack modeling framework accommodating multi-directional fracture at finite strains

A generic smeared crack modeling framework predicated on the deformation gradient decomposition (DGD) approach is proposed for use in dynamic fracture problems at finite strains, accommodating failure along multiple mutually orthogonal fracture planes embedded within an independently defined bulk material model. Within this constitutive framework, the traction equilibrium conditions imposed at each failure surface are used to determine the associated crack displacements stored as internal state variables. In general, the enforcement of interfacial equilibrium entails the implicit solution of a non-linear system of equations within the constitutive update procedure. However, if inertial effects arising due to the relative motion of the fractured material are incorporated within the model, the traction equilibrium conditions are shown to give rise to corresponding dynamic equations of motion governing the time-evolution of the crack opening displacements. For dynamic problems, an explicit time-integration procedure is devised to efficiently update the material state, subject to a set of internal frictionless contact constraints to prevent material inversion. Finally, the efficacy of the proposed modeling framework is investigated through several benchmark dynamic fracture problems run within the explicit finite element code DYNA3D.

42 ENGINEERING↗
NASA NTRS1986

Composite interlaminar fracture toughness: Three-dimensional finite element modeling for mixed mode 1, 2 and 3 fracture

A computational method/procedure is described which can be used to simulate individual and mixed mode interlaminar fracture progression in fiber composite laminates. Different combinations of Modes 1, 2, and 3 fracture are simulated by varying the crack location through the specimen thickness and by selecting appropriate unsymmetric laminate configurations. The contribution of each fracture mode to strain energy release rate is determined by the local crack closure methods while the mixed mode is determined by global variables. The strain energy release rates are plotted versus extending crack length, where slow crack growth, stable crack growth, and rapid crack growth regions are easily identified. Graphical results are presented to illustrate the effectiveness and versatility of the computational simulation for: (1) evaluating mixed-mode interlaminar fracture, (2) for identifying respective dominant parameters, and (3) for selecting possible simple test methods.

Murthy, P. L. N.↗
NASA NTRS1988

Composite interlaminar fracture toughness - Three-dimensional finite-element modeling for mixed mode I, II, and fracture

A computational method/procedure is described which can be used to simulate individual and mixed mode interlaminar fracture progression in fiber composite laminates. Different combinations of Modes 1, 2, and 3 fracture are simulated by varying the crack location through the specimen thickness and by selecting appropriate unsymmetric laminate configurations. The contribution of each fracture mode to strain energy release rate is determined by the local crack closure methods while the mixed mode is determined by global variables. The strain energy release rates are plotted versus extending crack length, where slow crack growth, stable crack growth, and rapid crack growth regions are easily identified. Graphical results are presented to illustrate the effectiveness and versatility of the computational simulation for: (1) evaluating mixed-mode interlaminar fracture, (2) for identifying respective dominant parameters, and (3) for selecting possible simple test methods.

Murthy, Pappu L. N.↗
DOE OSTI2021

Predicting Fracture Network Development in Crystalline Rocks

The geometric properties of fractures influence whether they propagate, arrest, or coalesce with other fractures. Thus, quantifying the relationship between fracture network characteristics may help predict fracture network development, and perhaps precursors to catastrophic failure. To constrain the relationship and predictability of fracture characteristics, we deform eight one centimeter tall rock cores under triaxial compression while acquiring in situ X-ray tomograms. The tomograms reveal precise measurements of the fracture network characteristics above the spatial resolution of 6.5 µm. We develop machine learning models to predict the value of each characteristic using the other characteristics, and excluding the macroscopic stress or strain imposed on the rock. The models predict fracture development more accurately in the experiments performed on granite and monzonite, than the experiments on marble. Fracture network development may be more predictable in these igneous rocks because their microstructure is more mechanically homogeneous than the marble, producing more systematic fracture development that is not strongly impeded by grain contacts and cleavage planes. The varying performance of the models suggest that fracture volume, length, and aperture are the most predictable of the characteristics, while fracture orientation is the least predictable. Orientation does not correlate with length, as suggested by the idea that the orientation evolves with increasing differential stress and thus fracture length. This difference between the observed and expected relationship between orientation and length highlights the influence of mechanical heterogeneities and local stress perturbations on fracture growth as fractures propagate, link, and coalesce.

58 GEOSCIENCES↗
NASA NTRS2011

Surrogate Modeling of High-Fidelity Fracture Simulations for Real-Time Residual Strength Predictions

A surrogate model methodology is described for predicting, during flight, the residual strength of aircraft structures that sustain discrete-source damage. Starting with design of experiment, an artificial neural network is developed that takes as input discrete-source damage parameters and outputs a prediction of the structural residual strength. Target residual strength values used to train the artificial neural network are derived from 3D finite element-based fracture simulations. Two ductile fracture simulations are presented to show that crack growth and residual strength are determined more accurately in discrete-source damage cases by using an elastic-plastic fracture framework rather than a linear-elastic fracture mechanics-based method. Improving accuracy of the residual strength training data does, in turn, improve accuracy of the surrogate model. When combined, the surrogate model methodology and high fidelity fracture simulation framework provide useful tools for adaptive flight technology.

Spear, Ashley D.↗
DOE OSTI2020

Hydraulic characterization of a fault zone from fracture distribution

A quantitative assessment of how faults control the migration of geofluids is critical in many areas of geosciences. Here we integrated geological fieldwork, quantitative analysis of the fractures distribution and numerical modeling to build a geometrical representation of a fault zone and to characterize its hydraulic properties. Our target is a fault located in the Majella Mountain (Italy). We collected 21 scan lines across the fault profile in order to characterize its architecture. The numerical modeling of the fracture network of the damage zones and their hydraulic parameters was performed using both commercial (Move ® ) and open source software (dfnWorks and PFLOTRAN). Move ® was used to build a representative model of the fault zone using fracture spacing as a proxy, and to model the hydraulic parameters of the different fault domains. dfnWorks and PFLOTRAN were employed to infer the hydraulic parameters of the damage zones of the fault and then upscale these properties to an equivalent continuum domain, suitable for fluid flow simulations through the whole fault zone. Our findings show how even in a relatively small area it is possible to describe changes in terms of hydraulic properties of a fault zone and to build models capable to represent these variations.

58 GEOSCIENCES↗
NASA NTRS1992

Scale effects in the response and failure of fiber reinforced composite laminates loaded in tension and in flexure

The feasibility of using scale model testing for predicting the full-scale behavior of flat composite coupons loaded in tension and beam-columns loaded in flexure is examined. Classical laws of similitude are applied to fabricate and test replica model specimens to identify scaling effects in the load response, strength, and mode of failure. Experiments were performed on graphite-epoxy composite specimens having different laminate stacking sequences and a range of scaled sizes. From the experiments it was deduced that the elastic response of scaled composite specimens was independent of size. However, a significant scale effect in strength was observed. In addition, a transition in failure mode was observed among scaled specimens of certain laminate stacking sequences. A Weibull statistical model and a fracture mechanics based model were applied to predict the strength scale effect since standard failure criteria cannot account for the influence of absolute specimen size on strength.

Jackson, Karen E.↗
DOE OSTI2021

Hydraulic Modeling of Induced and Propagated Fractures: Analysis of Flow and Pressure Data From Hydromechanical Experiments in the COSC-1 Deep Borehole in Crystalline Rock Near Åre, Sweden

To characterize the coupled hydromechanical behavior of rock fractures, the step-rate injection method for fracture in-situ properties (SIMFIP) was conducted with a specialized downhole probe developed by Guglielmi et.al. In June 2019, a field campaign was carried out near Åre, Sweden, where the SIMFIP probe was applied in the Collisional Orogeny in the Scandinavian Caledonides-1 scientific borehole to understand the dynamics of injection-induced fracture initiation, fracture opening, and shearing due to water injection-withdrawal in a borehole interval isolated by two packers. Three intervals were investigated at ~500 m depth: (a) an unfractured section (intact rock), (b) a section with non-conductive fractures, and (c) a section with hydraulically conductive fractures. Pressure, injection flow rate, and borehole wall displacement were simultaneously measured during the tests. In the present study, the geometry of the induced fracture and deformation of existing fractures at different time stages of the tests are determined based on a hydrologic model by using the measured pressure and flow data during each time stage of the experiment. A numerical model for the fluid flow within the fracture and the packed-off borehole interval is implemented within COMSOL Multiphysics. By matching model simulations with observed data for all three sections, estimates of the induced and propagated fractures' radius and aperture at successive time stages have been obtained in each case. We could also determine the non-linear relationship between fracture aperture and pressure for values above fracture opening pressures. The model results provide insights for the understanding of pressure-induced fracture initiation and propagation in crystalline rock.

58 GEOSCIENCES↗
NASA NTRS2011

Surrogate Modeling of High-Fidelity Fracture Simulations for Real-Time Residual Strength Predictions

A surrogate model methodology is described for predicting in real time the residual strength of flight structures with discrete-source damage. Starting with design of experiment, an artificial neural network is developed that takes as input discrete-source damage parameters and outputs a prediction of the structural residual strength. Target residual strength values used to train the artificial neural network are derived from 3D finite element-based fracture simulations. A residual strength test of a metallic, integrally-stiffened panel is simulated to show that crack growth and residual strength are determined more accurately in discrete-source damage cases by using an elastic-plastic fracture framework rather than a linear-elastic fracture mechanics-based method. Improving accuracy of the residual strength training data would, in turn, improve accuracy of the surrogate model. When combined, the surrogate model methodology and high-fidelity fracture simulation framework provide useful tools for adaptive flight technology.

Spear, Ashley D.↗
DOE OSTI2020

A higher-order finite element reactive transport model for unstructured and fractured grids

Abstract This work presents a new reactive transport framework that combines a powerful geochemistry engine with advanced numerical methods for flow and transport in subsurface fractured porous media. Specifically, the PhreeqcRM interface (developed by the USGS) is used to take advantage of a large library of equilibrium and kinetic aqueous and fluid-rock reactions, which has been validated by numerous experiments and benchmark studies. Fluid flow is modeled by the Mixed Hybrid Finite Element (FE) method, which provides smooth velocity fields even in highly heterogenous formations with discrete fractures. A multilinear Discontinuous Galerkin FE method is used to solve the multicomponent transport problem. This method is locally mass conserving and its second order convergence significantly reduces numerical dispersion. In terms of thermodynamics, the aqueous phase is considered as a compressible fluid and its properties are derived from a Cubic Plus Association (CPA) equation of state. The new simulator is validated against several benchmark problems (involving, e.g., Fickian and Nernst-Planck diffusion, isotope fractionation, advection-dispersion transport, and rock-fluid reactions) before demonstrating the expanded capabilities offered by the underlying FE foundation, such as high computational efficiency, parallelizability, low numerical dispersion, unstructured 3D gridding, and discrete fraction modeling.

58 GEOSCIENCES↗
DOE OSTI2023

A multiscale phase field fracture approach based on the non-affine microsphere model for rubber-like materials

Rubber-like materials have a broad scope of applications due to their unique properties like high stretchability and increased toughness. Hence, computational models for simulating their fracture behavior are paramount for designing them against failures. In this study, the phase field fracture approach is integrated with a multiscale polymer model for predicting the fracture behavior in elastomers. At the microscale, damaged polymer chains are modeled to be made up of a number of elastic chain segments pinned together. Using the phase field approach, the damage in the chains is represented using a continuous variable. Both the bond stretch internal energy and the entropic free energy of the chain are assumed to drive the damage, and the advantages of this assumption are expounded. A framework for utilizing the non-affine microsphere model for damaged systems is proposed here by considering the minimization of a hypothetical undamaged free energy, ultimately connecting the chain stretch to the macroscale deformation gradient. At the macroscale, a thermodynamically consistent formulation is derived in which the total dissipation is assumed to be mainly due to the rupture of molecular bonds. Using a monolithic scheme, the proposed model is numerically implemented and the resulting three-dimensional simulation predictions are compared with existing experimental data. The capability of the model to qualitatively predict the propagation of complex crack paths and quantitatively estimate the overall fracture behavior is verified. Additionally, the effect of the length scale parameter on the predicted fracture behavior is studied for an inhomogeneous system.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Investigation of Coupled Processes in Fractures and the Bordering Matrix via a Micro‐Continuum Reactive Transport Model

Abstract In multi‐mineral fractured rocks, the altered porous layer on the fracture surface resulting from preferential dissolution of the fast‐reacting minerals can have profound impacts on subsequent chemical‐physical alteration of the fractures. This study adopts the micro‐continuum approach to provide further understanding of reactive transport processes in the altered layer (AL), and mass exchanges with the bordering matrix and fracture. The modeling framework couples the Darcy‐Brinkman‐Stokes (DBS) solver in COMSOL Multiphysics and the geochemical modeling capability of CrunchFlow. Three‐dimensional steady state simulations with systematically varied chemical‐physical parameters of the AL were performed to examine the impacts of individual factors and processes. Our simulation results confirm previous observations that dissolution of the fast‐reacting mineral (i.e., calcite) is largely controlled by diffusion across the AL. We also show that dissolution of the slow‐reacting mineral (i.e., dolomite), which controls AL development and fracture enlargement, increases with surface area and has a complex dependence on different local rate‐limiting processes. In particular, advection can result in evident spatial variations in the local dissolution rates of dolomite, although it does not affect the bulk chemistry significantly. The difference in the spatial patterns between simulations with and without advection in the AL is more noticeable in the locations with smaller apertures, with up to 20% difference in local reaction rates. Therefore, it is important to include a full depiction of advection, diffusion, and reactions for accurately capturing local dynamics that control long‐term fracture evolution.

58 GEOSCIENCES↗
NASA NTRS1993

Environment enhanced fatigue crack propagation in metals: Inputs to fracture mechanics life prediction models

This report is a critical review of both environment-enhanced fatigue crack propagation data and the predictive capabilities of crack growth rate models. This information provides the necessary foundation for incorporating environmental effects in NASA FLAGRO and will better enable predictions of aerospace component fatigue lives. The review presents extensive literature data on 'stress corrosion cracking and corrosion fatigue.' The linear elastic fracture mechanics approach, based on stress intensity range (Delta(K)) similitude with microscopic crack propagation threshold and growth rates, provides a basis for these data. Results are presented showing enhanced growth rates for gases (viz., H2 and H2O) and electrolytes (e.g. NaCl and H2O) in aerospace alloys including: C-Mn and heat treated alloy steels, aluminum alloys, nickel-based superalloys, and titanium alloys. Environment causes purely time-dependent accelerated fatigue crack growth above the monotonic load cracking threshold (KIEAC) and promotes cycle-time dependent cracking below (KIEAC). These phenomenon are discussed in terms of hydrogen embrittlement, dissolution, and film rupture crack tip damage mechanisms.

Gangloff, Richard P.↗
DOE OSTI2020

Mesoscale simulations of High Burnup Fuel Fragmentation: Applying phase-field modeling to understand high burnup fuel fragmentation during loss of coolant accident conditions

During a loss-of-coolant accident (LOCA) in a nuclear power plant, a rapid increase in temperature is experienced in the fuel. In commercial light water reactors (LWRs), the high burnup structure (HBS) forms in regions of the UO 2 fuel where local burnup is high, characterized by smaller grain sizes and large, overpressurized bubbles. During a LOCA transient, fuel in the HBS region is susceptible to fine fragmentation, where the fuel breaks up into micron-size fragments (sometimes referred to as pulverization). To better understand the mechanisms behind this phenomenon, mesoscale computational modeling using the Idaho National Laboratory code Marmot has been employed to simulate the process of HBS formation and response to a LOCA transient. The process of complete and partial HBS formation was demonstrated in prototypical LWR conditions. The simulated HBS microstructures from these simulations were passed to a phase-field model of fracture. To provide a more direct comparison with experiments, a set of experiments conducted by Studsvik on fine fragmentation/pulverization during LOCA conditions was simulated using the BISON fuel performance code. The temperature profiles from the BISON simulations were passed to a phase-field model to determine the gas bubble pressure as a function of time. These pressure histories were used in a phase-field model of fracture to determine under what conditions and in which types of bubbles fine fragmentation was likely to occur.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Recent Developments at the Raft River Geothermal Field

The Raft River geothermal field, located in Cassia County in southwestern Idaho, is the site of a Department of Energy Enhanced Geothermal System project. U.S. Geothermal, Inc. currently produces about 11 MWe from Precambrian metamorphic rocks. These lie beneath ~5,000 ft of Quaternary and Tertiary volcaniclastic and volcanic deposits. Maximum temperatures range from 271°F to 300°F. Well RRG-9 ST1, the well targeted for stimulation is located approximately 1 mile south of the main bore field. The open hole section of the well, from 5,551 to 5,900 ft MD, consists of Precambrian Elba Quartzite, the stimulation target, granite and minor diabase. Prior to setting the casing acoustic, gamma ray, and density logs were run. After completing the well, a step rate/step down test was conducted. The maximum injection rate achieved was 18 bpm at a wellhead pressure of 1,150 psig. A borehole televiewer run in the open hole section showed evidence of more than eighty fractures. The majority of these fractures trend from N20⁰W to N20°E and dip from 40° to 60°W. Permeable fractures were encountered in the Elba Quartzite at 5,640-5,660 ft MD. Analysis of the injection test indicates that the minimum in-situ principal stress in this zone is 3,050-3,200 psi, corresponding to a fracture gradient of 0.59-0.62 psi/ft. A discrete fracture network model was developed using measured and inferred fracture orientations, distributions and dimensions. A three-phase stimulation program is proposed for RRG-9 ST-1. During the first two stages, water at 140°F, and later 40°F, will be injected to pre-condition and thermally fracture the reservoir. Here, the third stage will consist of a high rate, large volume conventional hydraulic stimulation.

Bradford, Jacob↗
DOE OSTI2020

The Influence of an Interlayer on Dual Hydraulic Fractures Propagation

Multi-cluster hydraulic fracturing of long-range horizontal wells is an approach for enhancing the productivity of low-permeability shale reservoirs. In this study, RFPA-Petrol (rock failure process analysis on petroleum problems) is applied for modeling hydraulic fracture propagation in multilayered formations. RFPA-Petrol based on coupled hydraulic-mechanical-damage (HMD) modeling was first tested by modeling a laboratory scale experiment on a physical (cement) model with a single completion. The modeling demonstrated the capability of RFPA-Petrol for simulating hydraulic fracture propagation. Then, we used RFPA-Petrol to investigate how the difference in material properties between oil-bearing layers and interlayers and the fracturing fluid properties influence the propagation of dual fractures in multilayered laboratory-scale models. In this case, the models with geological discontinuities in the vertical direction are strongly heterogeneous and RFPA-Petrol simulations successfully modeled the fracture configurations.

58 GEOSCIENCES↗
NASA NTRS2003

Fracture Evaluation of In-Situ Sensors for High Temperature Applications

The feasibility of fabricating an in-situ crack sensor for real-time detection of surface cracks propagating in engine components was evaluated using a computational fracture mechanics model. The in-situ sensor system would be required to: (1) Be capable of sustaining normal function in a severe environment; (2) Transmit a signal if a detected crack in the component was above a predetermined length, but below a critical length that would lead to failure; (3) Act neutrally upon the overall performance of the engine system and not interfere with the engine maintenance operations. In this work, fracture mechanics methodologies are used to identify the requirements for an in-situ sensor system that could withstand the engine operating environment, foreign object damage, and minimally degrade engine performance. A computational fracture mechanics model was developed to evaluate the feasibility of fabricating an in-situ crack sensor for real-time damage propagation detection in engine components.

Forth, Scott C.↗