A Time-Continuation Solver for Hydraulic Fracture Propagation

2021 ◽  
Author(s):  
Guotong Ren ◽  
Rami M. Younis
Keyword(s):  
A Priori ◽  
Fracture Propagation ◽  
A Priori Estimate ◽  
Extended Finite Element ◽  
Time Step ◽  
Linear Elastic ◽  
Fully Implicit ◽  
Stress Shadow ◽  
Propagation Regime ◽  
Frame Work

Abstract We present an efficient time-continuation scheme for fluid-driven fracture propagation problems in the frame-work of the extended finite element method (XFEM). The fully coupled, fully implicit hydro-mechanical system is solved in conjunction with the linear elastic fracture propagation criterion by the Newton-Raphson method. Therefore, at the end of each time-step solve, the model ensures the energy release rate of weakest fracture tips within the equilibrium propagation regime. Besides, an initialization procedure for newly created fracture space as well as a priori estimate of stress intensity factor (SIF) growth rates are also developed to further improve the solver performance. We validate the model by the analytical solution and extend the problem to the multiple fracture propagation where stress shadow phenomenon occur.

scholarly journals A Fully Implicit Finite Volume Scheme for a Seawater Intrusion Problem in Coastal Aquifers

Water ◽  
2020 ◽  
Vol 12 (6) ◽  
pp. 1639
Author(s):  
Abdelkrim Aharmouch ◽  
Brahim Amaziane ◽  
Mustapha El Ossmani ◽  
Khadija Talali
Keyword(s):  
Finite Volume ◽  
Coastal Aquifers ◽  
Nonlinear Partial Differential Equations ◽  
Time Integration ◽  
Mass Conservation ◽  
Coupled System ◽  
Finite Volume Scheme ◽  
Time Step ◽  
Volume Scheme ◽  
Fully Implicit

We present a numerical framework for efficiently simulating seawater flow in coastal aquifers using a finite volume method. The mathematical model consists of coupled and nonlinear partial differential equations. Difficulties arise from the nonlinear structure of the system and the complexity of natural fields, which results in complex aquifer geometries and heterogeneity in the hydraulic parameters. When numerically solving such a model, due to the mentioned feature, attempts to explicitly perform the time integration result in an excessively restricted stability condition on time step. An implicit method, which calculates the flow dynamics at each time step, is needed to overcome the stability problem of the time integration and mass conservation. A fully implicit finite volume scheme is developed to discretize the coupled system that allows the use of much longer time steps than explicit schemes. We have developed and implemented this scheme in a new module in the context of the open source platform DuMu X . The accuracy and effectiveness of this new module are demonstrated through numerical investigation for simulating the displacement of the sharp interface between saltwater and freshwater in groundwater flow. Lastly, numerical results of a realistic test case are presented to prove the efficiency and the performance of the method.

scholarly journals Fatigue Crack Growth Analysis with Extended Finite Element for 3D Linear Elastic Material

Metals ◽  
2021 ◽  
Vol 11 (3) ◽  
pp. 397
Author(s):  
Yahya Ali Fageehi
Keyword(s):  
Finite Element ◽  
Fatigue Crack ◽  
Fatigue Life ◽  
Crack Propagation ◽  
Crack Growth ◽  
Mixed Mode ◽  
Propagation Path ◽  
Loading Conditions ◽  
Extended Finite Element ◽  
Linear Elastic

This paper presents computational modeling of a crack growth path under mixed-mode loadings in linear elastic materials and investigates the influence of a hole on both fatigue crack propagation and fatigue life when subjected to constant amplitude loading conditions. Though the crack propagation is inevitable, the simulation specified the crack propagation path such that the critical structure domain was not exceeded. ANSYS Mechanical APDL 19.2 was introduced with the aid of a new feature in ANSYS: Smart Crack growth technology. It predicts the propagation direction and subsequent fatigue life for structural components using the extended finite element method (XFEM). The Paris law model was used to evaluate the mixed-mode fatigue life for both a modified four-point bending beam and a cracked plate with three holes under the linear elastic fracture mechanics (LEFM) assumption. Precise estimates of the stress intensity factors (SIFs), the trajectory of crack growth, and the fatigue life by an incremental crack propagation analysis were recorded. The findings of this analysis are confirmed in published works in terms of crack propagation trajectories under mixed-mode loading conditions.

scholarly journals The Effect of Perforation Spacing on the Variation of Stress Shadow

Energies ◽  
2021 ◽  
Vol 14 (13) ◽  
pp. 4040
Author(s):  
Weige Han ◽  
Zhendong Cui ◽  
Zhengguo Zhu
Keyword(s):  
Compressive Stress ◽  
Reservoir Rock ◽  
Hydraulic Fractures ◽  
Extended Finite Element ◽  
Gas Reservoir ◽  
Fracture Length ◽  
Shadow Effect ◽  
Stress Shadow ◽  
Initiation Pressure ◽  
Small Perforation

When the shale gas reservoir is fractured, stress shadows can cause reorientation of hydraulic fractures and affect the complexity. To reveal the variation of stress shadow with perforation spacing, the numerical model between different perforation spacing was simulated by the extended finite element method (XFEM). The variation of stress shadows was analyzed from the stress of two perforation centers, the fracture path, and the ratio of fracture length to spacing. The simulations showed that the reservoir rock at the two perforation centers is always in a state of compressive stress, and the smaller the perforation spacing, the higher the maximum compressive stress. Moreover, the compressive stress value can directly reflect the size of the stress shadow effect, which changes with the fracture propagation. When the fracture length extends to 2.5 times the perforation spacing, the stress shadow effect is the strongest. In addition, small perforation spacing leads to backward-spreading of hydraulic fractures, and the smaller the perforation spacing, the greater the deflection degree of hydraulic fractures. Additionally, the deflection angle of the fracture decreases with the expansion of the fracture. Furthermore, the perforation spacing has an important influence on the initiation pressure, and the smaller the perforation spacing, the greater the initiation pressure. At the same time, there is also a perforation spacing which minimizes the initiation pressure. However, when the perforation spacing increases to a certain value (the result of this work is about 14 m), the initiation pressure will not change. This study will be useful in guiding the design of programs in simultaneous fracturing.

scholarly journals Statistical Uncertainty of DNS in Geometries without Homogeneous Directions

2021 ◽  
Vol 11 (4) ◽  
pp. 1399
Author(s):  
Jure Oder ◽  
Cédric Flageul ◽  
Iztok Tiselj
Keyword(s):  
Channel Flow ◽  
A Priori ◽  
Time Integration ◽  
Navier Stokes ◽  
Statistical Uncertainty ◽  
Time Step ◽  
Order Of Magnitude ◽  
Averaging Time ◽  
The Individual ◽  
Time Required

In this paper, we present uncertainties of statistical quantities of direct numerical simulations (DNS) with small numerical errors. The uncertainties are analysed for channel flow and a flow separation case in a confined backward facing step (BFS) geometry. The infinite channel flow case has two homogeneous directions and this is usually exploited to speed-up the convergence of the results. As we show, such a procedure reduces statistical uncertainties of the results by up to an order of magnitude. This effect is strongest in the near wall regions. In the case of flow over a confined BFS, there are no such directions and thus very long integration times are required. The individual statistical quantities converge with the square root of time integration so, in order to improve the uncertainty by a factor of two, the simulation has to be prolonged by a factor of four. We provide an estimator that can be used to evaluate a priori the DNS relative statistical uncertainties from results obtained with a Reynolds Averaged Navier Stokes simulation. In the DNS, the estimator can be used to predict the averaging time and with it the simulation time required to achieve a certain relative statistical uncertainty of results. For accurate evaluation of averages and their uncertainties, it is not required to use every time step of the DNS. We observe that statistical uncertainty of the results is uninfluenced by reducing the number of samples to the point where the period between two consecutive samples measured in Courant–Friedrichss–Levy (CFL) condition units is below one. Nevertheless, crossing this limit, the estimates of uncertainties start to exhibit significant growth.

scholarly journals Positive solutions of higher order quasilinear elliptic equations

2002 ◽  
Vol 7 (8) ◽  
pp. 423-452
Author(s):  
Marcelo Montenegro
Keyword(s):  
A Priori Estimates ◽  
A Priori ◽  
A Priori Estimate ◽  
Higher Order ◽  
Dirichlet Boundary Conditions ◽  
Quasilinear Elliptic Equations ◽  
Radial Solutions ◽  
The Maximum Principle ◽  
Priori Estimates ◽  
Quasilinear Elliptic

The higher order quasilinear elliptic equation−Δ(Δp(Δu))=f(x,u)subject to Dirichlet boundary conditions may have unique and regular positive solution. If the domain is a ball, we obtain a priori estimate to the radial solutions via blowup. Extensions to systems and general domains are also presented. The basic ingredients are the maximum principle, Moser iterative scheme, an eigenvalue problem, a priori estimates by rescalings, sub/supersolutions, and Krasnosel'skiĭ fixed point theorem.

scholarly journals Extracting real-crack properties from non-linear elastic behaviour of rocks: abundance of cracks with dominating normal compliance and rocks with negative Poisson ratios

2017 ◽  
Vol 24 (3) ◽  
pp. 543-551 ◽  
Author(s):  
Vladimir Y. Zaitsev ◽  
Andrey V. Radostin ◽  
Elena Pasternak ◽  
Arcady Dyskin
Keyword(s):  
Experimental Data ◽  
Poisson Ratio ◽  
A Priori ◽  
Elastic Behaviour ◽  
Crack Model ◽  
Normal Compliance ◽  
S Wave ◽  
Linear Elastic ◽  
Non Linear ◽  
Linear Elastic Behaviour

Abstract. Results of examination of experimental data on non-linear elasticity of rocks using experimentally determined pressure dependences of P- and S-wave velocities from various literature sources are presented. Overall, over 90 rock samples are considered. Interpretation of the data is performed using an effective-medium description in which cracks are considered as compliant defects with explicitly introduced shear and normal compliances without specifying a particular crack model with an a priori given ratio of the compliances. Comparison with the experimental data indicated abundance (∼ 80 %) of cracks with the normal-to-shear compliance ratios that significantly exceed the values typical of conventionally used crack models (such as penny-shaped cuts or thin ellipsoidal cracks). Correspondingly, rocks with such cracks demonstrate a strongly decreased Poisson ratio including a significant (∼ 45 %) portion of rocks exhibiting negative Poisson ratios at lower pressures, for which the concentration of not yet closed cracks is maximal. The obtained results indicate the necessity for further development of crack models to account for the revealed numerous examples of cracks with strong domination of normal compliance. Discovering such a significant number of naturally auxetic rocks is in contrast to the conventional viewpoint that occurrence of a negative Poisson ratio is an exotic fact that is mostly discussed for artificial structures.

scholarly journals A PRIORI ESTIMATE FOR AN EQUATION WITH FRACTIONAL DERIVATIVES WITH DIFFERENT ORIGINS

2019 ◽  
pp. 41-47
Author(s):  
Л.М. Энеева
Keyword(s):  
Fractional Derivatives ◽  
Model Equation ◽  
A Priori ◽  
A Priori Estimate ◽  
Equation Of Motion ◽  
Point Boundary ◽  
Priori Estimate ◽  
Fractal Media ◽  
Different Origins ◽  
Variable Potential

В работе исследуется обыкновенное дифференциальное уравнение дробного порядка, содержащее композицию дробных производных с различными началами, с переменным потенциалом. Рассматриваемое уравнение выступает модельным уравнением движения во фрактальной среде. Для исследуемого уравнения доказана априорная оценка решения смешанной двухточечной краевой задачи. We consider an ordinary differential equation of fractional order with the composition of leftand right-sided fractional derivatives, and with variable potential. The considered equation is a model equation of motion in fractal media. We prove an a priori estimate for solutions of a mixed two-point boundary value problem for the equation under study.

Sign in / Sign up

Export Citation Format

Share Document