Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userA preconditioned mixed-FE scheme with stabilized Lagrange multiplier for frictional contact mechanics of crossing fractures in porous media
64
0
0.0
(
0
)
Added by
Hadi Hajibeygi
Publication date
2021
fields
Informatics Engineering
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
Simulation of contact mechanics in fractured media is of paramount important in the scope of computational mechanics. In this work, a preconditioned mixed-finite element scheme with Lagrange multipliers is proposed in the framework of constrained variational principle, which has the capability to handle frictional contact mechanics of the multi-crossing fractures. The slippage, opening and contact traction on fractures are calculated by the resulted saddle-point algebraic system. A novel treatment is devised to guarantee physical solutions at the intersected position of crossing fractures. A preconditioning technique is introduced to re-scale the resulting saddle-point algebraic system, to preserve the robustness of the system. An iteration strategy, namely monolithic-updated contact algorithm, is then designed to update the two primary unknowns (displacement and Lagrange multiplier) in one algebraic block. A series of numerical tests is conducted to study the contact mechanics of single- and multi-crossing fractures. Benchmark study is presented to verify the presented numerical method. Two tests with crossing fractures are studied, in which the slippage and opening can be calculated. The effects of crossing fractures on the deformation field can be observed in the calculated results, in which the variation of slippage/opening is analyzed by different loading conditions.
rate research
Read More
Accurate numerical simulations of interaction between fluid and solid play an important role in applications. The task is challenging in practical scenarios as the media are usually highly heterogeneous with very large contrast. To overcome this computational challenge, various multiscale methods are developed. In this paper, we consider a class of linear poroelasticity problems in high contrast heterogeneous porous media, and develop a mixed generalized multiscale finite element method (GMsFEM) to obtain a fast computational method. Our aim is to develop a multiscale method that is robust with respect to the heterogeneities and contrast of the media, and gives a mass conservative fluid velocity field. We will construct decoupled multiscale basis functions for the elastic displacement as well as fluid velocity. Our multiscale basis functions are local. The construction is based on some suitable choices of local snapshot spaces and local spectral decomposition, with the goal of extracting dominant modes of the solutions. For the pressure, we will use piecewise constant approximation. We will present several numerical examples to illustrate the performance of our method. Our results indicate that the proposed method is able to give accurate numerical solutions with a small degree of freedoms.
We introduce a phase-field method for continuous modeling of cracks with frictional contacts. Compared with standard discrete methods for frictional contacts, the phase-field method has two attractive features: (1) it can represent arbitrary crack geometry without an explicit function or basis enrichment, and (2) it does not require an algorithm for imposing contact constraints. The first feature, which is common in phase-field models of fracture, is attained by regularizing a sharp interface geometry using a surface density functional. The second feature, which is a unique advantage for contact problems, is achieved by a new approach that calculates the stress tensor in the regularized interface region depending on the contact condition of the interface. Particularly, under a slip condition, this approach updates stress components in the slip direction using a standard contact constitutive law, while making other stress components compatible with stress in the bulk region to ensure non-penetrating deformation in other directions. We verify the proposed phase-field method using stationary interface problems simulated by discrete methods in the literature. Subsequently, by allowing the phase field to evolve according to brittle fracture theory, we demonstrate the proposed methods capability for modeling crack growth with frictional contact.
The purposes of this work are to study the $L^{2}$-stability of a Navier-Stokes type model for non-stationary flow in porous media proposed by Hsu and Cheng in 1989 and to develop a Lagrange-Galerkin scheme with the Adams-Bashforth method to solve that model numerically.The stability estimate is obtained thanks to the presence of a nonlinear drag force term in the model which corresponds to the Forchheimer term. We derive the Lagrange-Galerkin scheme by extending the idea of the method of characteristics to overcome the difficulty which comes from the non-homogeneous porosity. Numerical experiments are conducted to investigate the experimental order of convergence of the scheme. For both simple and complex designs of porosities, our numerical simulations exhibit natural flow profiles which well describe the flow in non-homogeneous porous media.
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 present a barrier method for treating frictional contact on interfaces embedded in finite elements. The barrier treatment has several attractive features, including: (i) it does not introduce any additional degrees of freedom or iterative steps, (ii) it is free of inter-penetration, (iii) it avoids an ill-conditioned matrix system, and (iv) it allows one to control the solution accuracy directly. We derive the contact pressure from a smooth barrier energy function that is designed to satisfy the non-penetration constraint. Likewise, we make use of a smoothed friction law in which the stick-slip transition is described by a continuous function of the slip displacement. We discretize the formulation using the extended finite element method to embed interfaces inside elements, and devise an averaged surface integration scheme that effectively provides stable solutions without traction oscillations. Subsequently, we develop a way to tailor the parameters of the barrier method to embedded interfaces, such that the method can be used without parameter tuning. We verify and investigate the proposed method through numerical examples with various levels of complexity. The numerical results demonstrate that the proposed method is remarkably robust for challenging frictional contact problems, while requiring low cost comparable to that of the penalty method.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
