Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userIterative coupling of flow, geomechanics and adaptive phase-field fracture including level-set crack width approaches
117
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
In this work, we present numerical studies of fixed-stress iterative coupling for solving flow and geomechanics with propagating fractures in a porous medium. Specifically, fracture propagations are described by employing a phase-field approach. The extension to fixed-stress splitting to propagating phase-field fractures and systematic investigation of its properties are important enhancements to existing studies. Moreover, we provide an accurate computation of the fracture opening using level-set approaches and a subsequent finite element interpolation of the width. The latter enters as fracture permeability into the pressure diffraction problem which is crucial for fluid filled fractures. Our developments are substantiated with several numerical tests that include comparisons of computational cost for iterative coupling and nonlinear and linear iterations as well as convergence studies in space and time.
rate research
Read More
This paper concerns the analysis and implementation of a novel iterative staggered scheme for quasi-static brittle fracture propagation models, where the fracture evolution is tracked by a phase field variable. The model we consider is a two-field variational inequality system, with the phase field function and the elastic displacements of the solid material as independent variables. Using a penalization strategy, this variational inequality system is transformed into a variational equality system, which is the formulation we take as the starting point for our algorithmic developments. The proposed scheme involves a partitioning of this model into two subproblems; phase field and mechanics, with added stabilization terms to both subproblems for improved efficiency and robustness. We analyze the convergence of the proposed scheme using a fixed point argument, and find that under a natural condition, the elastic mechanical energy remains bounded, and, if the diffusive zone around crack surfaces is sufficiently thick, monotonic convergence is achieved. Finally, the proposed scheme is validated numerically with several bench-mark problems.
We compare the accuracy, convergence rate and computational cost of eigenerosion (EE) and phase-field (PF) methods. For purposes of comparison, we specifically consider the standard test case of a center-crack panel loaded in biaxial tension and assess the convergence of the energy error as the length scale parameter and mesh size tend to zero simultaneously. The panel is discretized by means of a regular mesh consisting of standard bilinear or Q1 elements. The exact stresses from the known analytical linear elastic solution are applied to the boundary. All element integrals over the interior and the boundary of the domain are evaluated exactly using the symbolic computation program Mathematica. When the EE inelastic energy is enhanced by means of Richardson extrapolation, EE is found to converge at twice the rate of PF and to exhibit much better accuracy. In addition, EE affords a one-order-of-magnitude computational speed-up over PF.
Mineral precipitation and dissolution processes in a porous medium can alter the structure of the medium at the scale of pores. Such changes make numerical simulations a challenging task as the geometry of the pores changes in time in an apriori unknown manner. To deal with such aspects, we here adopt a two-scale phase-field model, and propose a robust scheme for the numerical approximation of the solution. The scheme takes into account both the scale separation in the model, as well as the non-linear character of the model. After proving the convergence of the scheme, an adaptive two-scale strategy is incorporated, which improves the efficiency of the simulations. Numerical tests are presented, showing the efficiency and accuracy of the scheme in the presence of anisotropies and heterogeneities.
We consider a phase-field fracture propagation model, which consists of two (nonlinear) coupled partial differential equations. The first equation describes the displacement evolution, and the second is a smoothed indicator variable, describing the crack position. We propose an iterative scheme, the so-called $L$-scheme, with a dynamic update of the stabilization parameters during the iterations. Our algorithmic improvements are substantiated with two numerical tests. The dynamic adjustments of the stabilization parameters lead to a significant reduction of iteration numbers in comparison to constant stabilization values.
A new fast algebraic method for obtaining an $mathcal{H}^2$-approximation of a matrix from its entries is presented. The main idea behind the method is based on the nested representation and the maximum-volume principle to select submatrices in low-rank matrices. A special iterative approach for the computation of so-called representing sets is established. The main advantage of the method is that it uses only the hierarchical partitioning of the matrix and does not require special proxy surfaces to be selected in advance. The numerical experiments for the electrostatic problem and for the boundary integral operator confirm the effectiveness and robustness of the approach. The complexity is linear in the matrix size and polynomial in the ranks. The algorithm is implemented as an open-source Python package that is available online.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
