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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Effect of Unsaturated Flow Modes on Partitioning Dynamics of Gravity-Driven Flow at a Simple Fracture Intersection: Laboratory Study and Three-Dimensional Smoothed Particle Hydrodynamics Simulations

In this work, we study unsaturated gravity-driven free-surface flows on a synthetic surface intersected by a horizontal fracture and the importance of flow modes (droplet; rivulet) on the partitioning behavior at the fracture intersection. We present laboratory experiments, three-dimensional smoothed particle hydrodynamics (SPH) simulations using a heavily parallelized code, and an analytical solution. The flow-rate-dependent mode switching from droplets to rivulets is reproduced by the SPH model, and the transition ranges agree with the laboratory experiments. We show that flow modes heavily influence the bypass behavior. Flows favoring the formation of droplets exhibit a much stronger bypass capacity compared to rivulet flows, where nearly the whole fluid mass is initially stored within the horizontal fracture. This behavior is demonstrated for a multi-inlet laboratory setup where the inlet- specific flow rate is chosen so that either a droplet or rivulet flow persists. The effect of fluid buffering within the horizontal fracture is presented in terms of dimensionless fracture in- flow so that characteristic scaling regimes can be recovered. For both cases (rivulets and droplets), the flow within the horizontal fracture transitions into a Washburn regime until a critical threshold is reached and the bypass efficiency increases. For rivulet flows, the initial filling of the horizontal fracture is described by classical plug flow. Meanwhile, for droplet flows, a size-dependent partitioning behavior is observed, and the filling of the fracture takes longer. For the case of rivulet flow, we provide an analytical solution that demonstrates the existence of classical Washburn flow within the horizontal fracture.

unsaturated flow, fractured media, smoothed partic↗
DOE OSTI2024

Stochastic reconstruction of fracture network pattern using spatial point processes

Fracture spatial patterns can strongly affect fluid flow in the subsurface. Proximity and distribution of fractures control reservoir flow behavior over various length scales. In many studies, however, simplified geometrical patterns are generated for fractures which may lead to unrealistic subsurface models. Here we introduce a new method for characterization and modeling of fracture spatial patterns based on outcrop observations. We use Ripley's K-function to characterize the arrangement of fracture barycenters and intersection points over various length scales. In addition, we use semivariograms to quantify spatial correlation in fracture intensity maps. Using this information, we develop a stochastic algorithm that generates two-dimensional fracture network realizations with spatial properties similar to those of a real fracture network measured in the field. Numerical simulation models indicate that the generated fracture realizations exhibit similar flow behaviors as that of the original fracture network. Furthermore, such modeling tools expand and improve our capability in building representative fracture models and in quantification of uncertainty in naturally and hydraulically fractured reservoirs.

58 GEOSCIENCES↗
DOE OSTI2019

Characterizing the Influence of Fracture Density on Network Scale Transport

The topology of natural fracture networks is inherently linked to the structure of the fluid velocity field and transport therein. Here we study the impact of network density on flow and transport behaviors. We stochastically generate fracture networks of varying density and simulate flow and transport with a discrete fracture network model, which fully resolves network topology at the fracture scale. We study conservative solute trajectories with Lagrangian particle tracking and find that as fracture density decreases, solute channelization to large local fractures increases, thereby reducing plume spreading. Furthermore, in sparse networks mean particle travel distance increases and local network features, such as velocity zones where flow is counter to the primary pressure gradient, become increasingly important for transport. As the network density increases, network statistics homogenize and such local features have a reduced impact. We quantify local topological influence on transport behavior with an effective tortuosity parameter, which measures the ratio of total advective distance to linear distance at the fracture scale; large tortuosity values are correlated to slow-velocity regions. Additionally, these large tortuosity, slow-velocity regions delay downstream transport and enhance tailing on particle breakthrough curves. Finally, we predict transport with an upscaled, Bernoulli spatial Markov random walk model and parameterize local topological influences with a novel tortuosity parameter. Furthermore, Bernoulli model predictions improve when sampling from a tortuosity distribution, as opposed to a fixed value as has previously been done, suggesting that local network topological features must be carefully considered in upscaled modeling efforts of fracture network systems.

58 GEOSCIENCES↗
DOE OSTI2019

DECOVALEX2019 Task C: Results of Step 2 modeling

The work for Step 1 performed at Sandia National Laboratories and reported in Section 7 has been updated to incorporate new data and to conduct new simulations using a new larger base case domain. The new simulations also include statistical analysis for different fracture realizations. A sensitivity analysis was also conducted to the study of the effect of domain size. A much larger mesh was selected to minimize boundary effects. The DFN model was upscaled to the new base case domain and the much larger domain to generate relevant permeability and porosity fields for each case. The calculations updated for Step 2 are described in Section 12.1. New calculations have also been conducted to model the flooding of the CTD and the resulting pressure recovery. The modeling includes matching of pressure and chloride experimental data at the six observation locations in Well 12M133. The modeling was done for the 10 fracture realizations. The Step 2 recovery simulations are described in Section 12.2. The Step 2 work is summarized in Section 12.3.

42 ENGINEERING↗
DOE OSTI2018

Investigation on Cracking Behavior and Development of a Fracture Model of Ti-47Al-2Nb-2Cr Alloy During Hot Deformation

When alloys are strained to a certain limit under the plastic deformation, ductile fracture occurs. Whether the required deformation can be put into effect before fracture appears has become a serious concern. A combination of compression tests, ductile fracture criteria as well as finite element simulation (FEM) were employed in this study to establish the fracture initiation prediction model of Ti-47Al-2Nb-2Cr alloy. To obtain the critical fracture strain, a technique was proposed to calculate the data for deformed cylindrical samples through isothermal compression tests in the temperature range between 1223 and 1423 K with strain rates ranging from 0.001 to 1 s{sup −1}. Five commonly used ductile fracture criteria were compared by FE simulation of compressing tests. It was found that the Brozzo criterion was found to be suitable and accordingly a failure criterion, considering temperature and strain ratio was established. Hot working parameters (critical fracture value and critical strain) in deformation could be obtained through the finite element simulation, revealing that this newly established failure criterion could accurately predict cracking initiation in compressing deformation of Ti-47Al-2Nb-2Cr alloy ingot.

36 MATERIALS SCIENCE↗
DOE OSTI2015

dfnWorks: A discrete fracture network framework for modeling subsurface flow and transport

DFNWORKS is a parallelized computational suite to generate three-dimensional discrete fracture networks (DFN) and simulate flow and transport. Developed at Los Alamos National Laboratory over the past five years, it has been used to study flow and transport in fractured media at scales ranging from millimeters to kilometers. The networks are created and meshed using DFNGEN, which combines FRAM (the feature rejection algorithm for meshing) methodology to stochastically generate three-dimensional DFNs with the LaGriT meshing toolbox to create a high-quality computational mesh representation. The representation produces a conforming Delaunay triangulation suitable for high performance computing finite volume solvers in an intrinsically parallel fashion. Flow through the network is simulated in dfnFlow, which utilizes the massively parallel subsurface flow and reactive transport finite volume code PFLOTRAN. A Lagrangian approach to simulating transport through the DFN is adopted within DFNTRANS to determine pathlines and solute transport through the DFN. Example applications of this suite in the areas of nuclear waste repository science, hydraulic fracturing and CO 2 sequestration are also included.

42 ENGINEERING↗
DOE OSTI2014

A Nonlocal Peridynamic Plasticity Model for the Dynamic Flow and Fracture of Concrete.

A nonlocal, ordinary peridynamic constitutive model is formulated to numerically simulate the pressure-dependent flow and fracture of heterogeneous, quasi-brittle ma- terials, such as concrete. Classical mechanics and traditional computational modeling methods do not accurately model the distributed fracture observed within this family of materials. The peridynamic horizon, or range of influence, provides a characteristic length to the continuum and limits localization of fracture. Scaling laws are derived to relate the parameters of peridynamic constitutive model to the parameters of the classical Drucker-Prager plasticity model. Thermodynamic analysis of associated and non-associated plastic flow is performed. An implicit integration algorithm is formu- lated to calculate the accumulated plastic bond extension and force state. The gov- erning equations are linearized and the simulation of the quasi-static compression of a cylinder is compared to the classical theory. A dissipation-based peridynamic bond failure criteria is implemented to model fracture and the splitting of a concrete cylinder is numerically simulated. Finally, calculation of the impact and spallation of a con- crete structure is performed to assess the suitability of the material and failure models for simulating concrete during dynamic loadings. The peridynamic model is found to accurately simulate the inelastic deformation and fracture behavior of concrete during compression, splitting, and dynamically induced spall. The work expands the types of materials that can be modeled using peridynamics. A multi-scale methodology for simulating concrete to be used in conjunction with the plasticity model is presented. The work was funded by LDRD 158806.

Vogler, Tracy↗
NASA NTRS2012

Hypervelocity Impact Initiation of Explosive Transfer Lines

The Gemini, Apollo and Space Shuttle spacecraft utilized explosive transfer lines (ETL) in a number of applications. In each case the ETL was located behind substantial structure and the risk of impact initiation by micrometeoroids and orbital debris was negligible. A current NASA program is considering an ETL to synchronize the actuation of pyrobolts to release 12 capture latches in a contingency. The space constraints require placing the ETL 50 mm below the 1 mm thick 2024-T72 Whipple shield. The proximity of the ETL to the thin shield prompted analysts at NASA to perform a scoping analysis with a finite-difference hydrocode to calculate impact parameters that would initiate the ETL. The results suggest testing is required and a 12 shot test program with surplused Shuttle ETL is scheduled for February 2012 at the NASA White Sands Test Facility. Explosive initiation models are essential to the analysis and one exists in the CTH library for HNS I, but not the HNS II used in the Shuttle 2.5 gr/ft rigid shielded mild detonating cord (SMDC). HNS II is less sensitive than HNS I so it is anticipated that these results using the HNS I model are conservative. Until the hypervelocity impact test results are available, the only check on the analysis was comparison with the Shuttle qualification test result that a 22 long bullet would not initiate the SMDC. This result was reproduced by the hydrocode simulation. Simulations of the direct impact of a 7 km/s aluminum ball, impacting at 0 degree angle of incidence, onto the SMDC resulted in a 1.5 mm diameter ball initiating the SMDC and 1.0 mm ball failing to initiate it. Where one 1.0 mm ball could not initiate the SMDC, a cluster of six 1.0 mm diameter aluminum balls striking simultaneously could. Thus the impact parameters that will result in initiating SMDC located behind a Whipple shield will depend on how well the shield fragments the projectile and spreads the fragments. An end-to-end simulation of the impact of an aluminum ball onto a Whipple shield covering SMDC is problematic due to the hydrocode fracture models. Regardless, two simulations were performed resulting in a 5 mm ball initiating the SMDC and a 4 mm ball failing to initiate the SMDC.

Bjorkman, Michael D.↗
DOE OSTI2018

Hydraulic Fracture Modeling in Support of EGS Collab Treatment Designs

This paper presents ongoing hydraulic fracturing modeling work in support of the Department of Energy's EGS Collab project. CFRAC, a hydraulic fracturing research code with discrete fracture modeling capabilities was used in the modeling. The models are run under different geomechanical conditions to predict the hydraulic fracturing treatments to be pumped within EGS Collab's experiments in Sanford Underground Research Facility (SURF). Effect of the parameters such as a varying stress field and the presence of natural fractures are investigated. The simulations show that when injecting at 0.1 L/s for 180 s with uniform stress field conditions, a penny-shaped fracture with maximum aperture of 0.14 mm at radial extent of about 8 m is created. A disturbance in the stress field, caused by the presence of an offset mine drift, forced the fractures to grow asymmetrically. The presence of an intersecting natural fracture halted growth of the hydraulic fracture. These results agree with other modeling work done within the EGS Collab team.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2024

Fracture Network Influence on Rock Damage and Gas Transport following an Underground Explosion

Simulations of rock damage and gas transport following underground explosions that omit preexisting fracture networks in the subsurface cannot fully characterize the influence of geo-structural variability on gas transport. Previous studies do not consider the impact that fracture network structure and variability have on gas seepage. In this study, we develop a sequentially coupled, axi-symmetric model to look at the damage pattern and resulting gas breakthrough curves following an underground explosion given different fracture network realizations. We simulate 0.327 and 0.164 kT chemical explosives with burial depths of 100 m for 90 stochastically generated fracture networks. Gases quickly reach the surface in 30% of the higher yield simulations and 5% of the lower yield simulations. The fast breakthrough can be attributed to the formation of connected pathways between fractures to the surface. The formation of a connected damage pathway to the surface is not clearly correlated with the fracture intensity (P32) in our simulations. Breakthrough curves with slower transport are highly variable depending on the fracture network sample. The variability in the breakthrough behavior indicates that ignoring the influence of fracture networks on rock damage, which strongly influences the hydraulic properties following an underground explosion, will likely lead to a large underestimation of the uncertainty in the gas transport to the surface. This work highlights the need for incorporation of fracture networks into models for accurately predicting gas seepage following underground explosions.

58 GEOSCIENCES↗
DOE OSTI2016

Examining the microtexture evolution in a hole-edge punched into 780 MPa grade hot-rolled steel

The deformation behavior in the hole-edge of 780 MPa grade hot-rolled steel during the punching process was investigated via microstructure characterization and computational simulation. Microstructure characterization was conducted to observe the edges of punched holes through the thickness direction, and electron back-scattered diffraction (EBSD) was used to analyze the heterogeneity of the deformation. Finite element analysis (FEA) that could account for a ductile fracture criterion was conducted to simulate the deformation and fracture behaviors of 780 MPa grade hot-rolled steel during the punching process. Calculation of rotation rate fields at the edges of the punched holes during the punching process revealed that metastable orientations in Euler space were confined to specific orientation groups. Rotation-rate fields effectively explained the stability of the initial texture components in the hole-edge region during the punching process. A visco-plastic self-consistent (VPSC) polycrystal model was used to calculate the microtexture evolution in the hole-edge region during the punching process. FEA revealed that the heterogeneous effective strain was closely related to the heterogeneity of the Kernel average misorientation (KAM) distribution in the hole-edge region. A simulation of the deformation microtexture evolution in the hole-edge region using a VPSC model was in good agreement with the experimental results. - Highlights: •We analyzed the microstructure in a hole-edge punched in HR 780HB steel. •Rotation rate fields revealed the stability of the initial texture components. •Heterogeneous effective stain was closely related to the KAM distribution. •VPSC model successfully simulated the deformation microtexture evolution.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Efficient Monte Carlo With Graph-Based Subsurface Flow and Transport Models

Simulating flow and transport in fractured porous media frequently involves solving numerical discretizations of partial differential equations with a large number of degrees of freedom using discrete fracture network (DFN) models. Uncertainty in the properties of the fracture network that controls flow and transport requires a large number of DFN simulations to statistically describe quantities of interest. However, the computational cost of solving more than a few realizations of a large DFN can be intractable. As a means of circumventing this problem, we utilize both a high-fidelity DFN model and a graphbased model of flow and transport in combination with a multifidelity Monte Carlo (MC) method to reduce the number of high-fidelity simulations that are needed to obtain an accurate estimate of the quantity of interest. We demonstrate the approach by estimating quantiles of the breakthrough time for a conservative tracer in an ensemble of fractured porous media. Our results demonstrate that a multifidelity MC estimate, whose computational cost is equal to the cost of 10 DFN simulations, can be as accurate as a standard MC estimate that utilizes 1,000 DFN simulations. Thus the combination of our graphbased model with multifidelity MC estimates effectively reduces the computational cost of the problem by a factor of approximately 100.

58 GEOSCIENCES↗
DOE OSTI2023

Impact of artificial topological changes on flow and transport through fractured media due to mesh resolution

Abstract We performed a set of numerical simulations to characterize the interplay of fracture network topology, upscaling, and mesh refinement on flow and transport properties in fractured porous media. We generated a set of generic three-dimensional discrete fracture networks at various densities, where the radii of the fractures were sampled from a truncated power-law distribution, and whose parameters were loosely based on field site characterizations. We also considered five network densities, which were defined using a dimensionless version of density based on percolation theory. Once the networks were generated, we upscaled them into a single continuum model using the upscaled discrete fracture matrix model presented by Sweeney et al. (2019). We considered steady, isothermal pressure-driven flow through each domain and then simulated conservative, decaying, and adsorbing tracers using a pulse injection into the domain. For each simulation, we calculated the effective permeability and solute breakthrough curves as quantities of interest to compare between network realizations. We found that selecting a mesh resolution such that the global topology of the upscaled mesh matches the fracture network is essential. If the upscaled mesh has a connected pathway of fracture (higher permeability) cells but the fracture network does not, then the estimates for effective permeability and solute breakthrough will be incorrect. False connections cannot be eliminated entirely, but they can be managed by choosing appropriate mesh resolution and refinement for a given network. Adopting octree meshing to obtain sufficient levels of refinement leads to fewer computational cells (up to a 90% reduction in overall cell count) when compared to using a uniform resolution grid and can result in a more accurate continuum representation of the true fracture network.

58 GEOSCIENCES↗
DOE OSTI2019

Stress-Induced Anomalous Transport in Natural Fracture Networks

In this work we investigate the effects of geological stress on fluid flow and tracer transport in natural fracture networks. We show the emergence of non-Fickian (anomalous) transport from the interplay among fracture network geometry, aperture heterogeneity, and geological stress. In this study, we extract the fracture network geometry from the geological map of an actual rock outcrop, and we simulate the geomechanical behavior of fractured rock using a hybrid finite-discrete element method. We analyze the impact of stress on the aperture distribution, fluid flow field, and tracer transport properties. Both stress magnitude and orientation have strong effects on the fracture aperture field, which in turn affects fluid flow and tracer transport through the system. We observe that stress anisotropy may cause significant shear dilation along long, curved fractures that are preferentially oriented to the stress loading. This, in turn, induces preferential flow paths and anomalous early arrival of tracers. An increase in stress magnitude enhances aperture heterogeneity by introducing more small apertures, which exacerbates late-time tailing. This effect is stronger when there is higher heterogeneity in the initial aperture field. To honor the flow field with strong preferential flow paths, we extend the Bernoulli Continuous Time Random Walk model to incorporate dual velocity correlation length scales. The proposed upscaled transport model captures anomalous transport through stressed fracture networks and agrees quantitatively with the high-fidelity numerical simulations.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2002

Probabilistic Fatigue: Computational Simulation

Fatigue is a primary consideration in the design of aerospace structures for long term durability and reliability. There are several types of fatigue that must be considered in the design. These include low cycle, high cycle, combined for different cyclic loading conditions - for example, mechanical, thermal, erosion, etc. The traditional approach to evaluate fatigue has been to conduct many tests in the various service-environment conditions that the component will be subjected to in a specific design. This approach is reasonable and robust for that specific design. However, it is time consuming, costly and needs to be repeated for designs in different operating conditions in general. Recent research has demonstrated that fatigue of structural components/structures can be evaluated by computational simulation based on a novel paradigm. Main features in this novel paradigm are progressive telescoping scale mechanics, progressive scale substructuring and progressive structural fracture, encompassed with probabilistic simulation. These generic features of this approach are to probabilistically telescope scale local material point damage all the way up to the structural component and to probabilistically scale decompose structural loads and boundary conditions all the way down to material point. Additional features include a multifactor interaction model that probabilistically describes material properties evolution, any changes due to various cyclic load and other mutually interacting effects. The objective of the proposed paper is to describe this novel paradigm of computational simulation and present typical fatigue results for structural components. Additionally, advantages, versatility and inclusiveness of computational simulation versus testing are discussed. Guidelines for complementing simulated results with strategic testing are outlined. Typical results are shown for computational simulation of fatigue in metallic composite structures to demonstrate the versatility of this novel paradigm in predicting a priori fatigue life.

Chamis, Christos C.↗