Convergence analysis of fixed stress split iterative scheme for anisotropic poroelasticity with tensor Biot parameter

We perform a convergence analysis of the fixed stress split iterative scheme for the Biot system modeling coupled flow and deformation in anisotropic poroelastic media with tensor Biot parameter. The fixed stress split iterative scheme solves the flow subproblem with all components of the stress tensor frozen using a multipoint flux mixed finite element method, followed by the poromechanics subproblem using a conforming Galerkin method in every coupling iteration at each time step. The coupling iterations are repeated until convergence and Backward Euler is employed for time marching. The convergence analysis is based on studying the equations satisfied by the difference of iterates to show that the fixed stress split iterative scheme for anisotropic poroelasticity with Biot tensor is contractive. We also demonstrate that the scheme is numerically convergent using the classical Mandel’s problem solution for transverse isotropy.     

Robust iterative schemes for non-linear poromechanics

We consider a non-linear extension of Biot’s model for poromechanics, wherein both the fluid flow and mechanical deformation are allowed to be non-linear. Specifically, we study the case when the volumetric stress and the fluid density are non-linear functions satisfying certain assumptions. We perform an implicit discretization in time (backward Euler) and propose two iterative schemes for solving the non-linear problems appearing within each time step: a splitting algorithm extending the undrained split and fixed stress methods to non-linear problems, and a monolithic L-scheme. The convergence of both schemes are shown rigorously. Illustrative numerical examples are presented to confirm the applicability of the schemes and validate the theoretical results.     

A virtual element method for the two-phase flow of immiscible fluids in porous media

A primal C⁰-conforming virtual element discretization for the approximation of the bidimensional two-phase flow of immiscible fluids in porous media using general polygonal meshes is discussed. This work investigates the potentialities of the Virtual Element Method (VEM) in solving this specific problem of immiscible fluids in porous media involving a time-dependent coupled system of non-linear partial differential equations. The performance of the fully discrete scheme is thoroughly analysed testing it on general meshes considering both a regular problem and more realistic benchmark problems that are of interest for physical and engineering applications.

     

A two-scale poromechanical model for cohesive rocks

This paper presents a computational homogenization approach in the framework of poromechanics. A fully coupled hydromechanics problem is formulated at the macroscopic scale. The constitutive equations are replaced by results of numerical computations on a Representative Elementary Volume in order to obtain the overall stress of the mixture as well as its transmissivity properties. At the microscale, the material is assumed to be composed of an assembly of hyperelastic grains connected by cohesive interfaces. These interfaces are also channels of a network where the fluid can percolate. The fluid acts on the boundaries of the grains and influences the behavior of the cohesive interfaces. Conversely, the opening of the interfaces induces changes in the transmissivity properties of the corresponding channels. This yields a fully coupled hydromechanical problem at the microscopic scale. The finite element method is considered for the numerical solutions at both scales, the present approach extending the purely mechanical FE² scheme to the coupled hydromechanical framework. The local macroscopic behavior resulting from the homogenization scheme is illustrated on different numerical tests. The results clearly show the coupling between damage and fluid permeability in the overall response, as a consequence of the small-scale interaction between the action of the percolating fluid, the deformation of the solid skeleton, and the failure of the cohesive interfaces.     

A general method for simulating reactive dissolution in carbonate rocks with arbitrary geometry

The two-scale continuum model is widely used in simulating the reactive dissolution process and predicting the optimum injection rate for carbonate reservoir acidizing treatment. The numerical methods of this model are currently based on structured grids, which are not applicable for complicated geometries. In this study, a general numerical scheme for simulating a reactive flow problem on both structured and unstructured grids is presented based on the finite volume method (FVM). The convection and diffusion terms involved in the reactive flow model are discretized by using the upwind scheme and two-point flux approximation (TPFA), respectively. The location of the centroid node inside each control volume is moved by using an optimization algorithm to make the connections with the surrounding elements as orthogonal as possible, which systematically improves the accuracy of the TPFA scheme. Additionally, in order to avoid the computational complexity resulting from the discretization of the non-linear term, the mass balance equation is only discretized in the spatial domain to get a set of ordinary differential equations (ODEs). These ODEs are coupled with the reaction equations and then solved using the numerical algorithm on ODEs. The accuracy and efficiency of the proposed method are studied by comparing the results obtained from the proposed numerical method with previous experimental and numerical results. This comparison indicates that, compared with the previous methods, the proposed method predicts the wormhole structure more accurately. Finally, the presented method is used to check the effect of the domain geometry, and it is found that the geometry of the flow domain has no effect on the optimum injection velocity, but the radial domain requires a larger breakthrough volume than the linear domain when other parameters are fixed.     

An alternative approach for quasi‐static large deformation analysis of saturated porous media using meshfree method

An alternative coupled large deformation formulation combined with a meshfree approach is proposed for flow–deformation analysis of saturated porous media. The formulation proposed is based on the Updated Lagrangian (UL) approach, except that the spatial derivatives are defined with respect to the configuration of the medium at the last time step rather than the configuration at the last iteration. In this way, the Cauchy stresses are calculated directly, rendering the second Piola–Kirchhoff stress tensor not necessary for the numerical solution of the equilibrium equations. Moreover, in contrast with the UL approach, the nodal shape function derivatives are calculated once in each time step and stored for use in subsequent iterations, which reduces the computational cost of the algorithm. Stress objectivity is satisfied using the Jaumann stress rate, and the spatial discretisation of the governing equations is achieved using the standard Galerkin method. The equations of equilibrium are satisfied directly, and the nonlinear parts of the system matrix are derived independent of the stresses of the medium resulting in a stable numerical algorithm. Temporal discretisation is effected based on a three‐point approximation technique that avoids spurious ripple effects and has second‐order accuracy. The radial point interpolation method is used to construct the shape functions. The application of the formulation and the significance of large deformation effects on the numerical results are demonstrated through several numerical examples. Copyright © 2015 John Wiley & Sons, Ltd.     

Thermo-hydro-mechanical modeling of fracturing porous media with two-phase fluid flow using X-FEM technique

In this paper, a fully coupled thermo-hydro-mechanical model is presented for two-phase fluid flow and heat transfer in fractured/fracturing porous media using the extended finite element method. In the fractured porous medium, the traction, heat, and mass transfer between the fracture space and the surrounding media are coupled. The wetting and nonwetting fluid phases are water and gas, which are assumed to be immiscible, and no phase-change is considered. The system of coupled equations consists of the linear momentum balance of solid phase, wetting and nonwetting fluid continuities, and thermal energy conservation. The main variables used to solve the system of equations are solid phase displacement, wetting fluid pressure, capillary pressure, and temperature. The fracture is assumed to impose the strong discontinuity in the displacement field and weak discontinuities in the fluid pressure, capillary pressure, and temperature fields. The mode I fracture propagation is employed using a cohesive fracture model. Finally, several numerical examples are solved to illustrate the capability of the proposed computational algorithm. It is shown that the effect of thermal expansion on the effective stress can influence the rate of fracture propagation and the injection pressure in hydraulic fracturing process. Moreover, the effect of thermal loading is investigated properly on fracture opening and fluids flow in unsaturated porous media, and the convective heat transfer within the fracture is captured successfully. It is shown how the proposed computational model is capable of modeling the fully coupled thermal fracture propagation in unsaturated porous media.     

Implementation of a hydromechanical elastoplastic constitutive model for fully coupled dynamic analysis of unsaturated soils and its validation using centrifuge test results

Prediction of unsaturated soil behavior during earthquake loading has received increasing attention in geotechnical engineering research and practice in recent years. Development of a fully coupled analysis procedure incorporating a coupled hydromechanical elastoplastic constitutive model for dynamic analysis of unsaturated soils has, however, been limited. This paper presents the implementation of a coupled hydromechanical elastoplastic constitutive model into a fully coupled dynamic analysis procedure and its validation using a centrifuge test. First, the fully coupled finite element equations governing the dynamic behavior of unsaturated soils with the solid skeleton displacement, pore water pressure, and pore air pressure as nodal unknowns are briefly presented. The closest point projection method is then utilized to implement the coupled hydromechanical elastoplastic constitutive model into the finite element equations. The constitutive model includes hysteresis in soil–water characteristic curves, cyclic elastoplasticity of the solid skeleton, and the coupling mechanisms between the SWCCs and the solid skeleton. Finally, the analysis procedure is validated using the results from a dynamic centrifuge test on an embankment constructed of compacted unsaturated silt subjected to base shaking. Reasonable comparisons between the predicted and measured accelerations, settlements, and deformed shapes are obtained.

     

Thermo-mechanical performance of layered transversely isotropic media around a cylindrical/tubular heat source

We develop a new numerical model based on a precise integration method to investigate the coupled thermo-mechanical performance of layered transversely isotropic media around a cylindrical/tubular heat source. To obtain the relational matrices of the extended precise integration method, we first convert the governing equations of the problem into ordinary differential matrix equations through the Laplace–Hankel transform. Then, the cylindrical heat source is divided into a series of plane heat sources, and the plane temperature load term is added to the state vector between layer elements. By combining the layer elements, we build a layered transversely isotropic numerical model containing a cylindrical heat source in the transformed domain. Finally, we solve the model in the transformed domain and obtain the solution of the problem in the real domain through the Laplace–Hankel transform inversion. The accuracy of this method is verified by comparing the solutions with the results of the analytical method and the finite element method. Then, we study the influence of the anisotropy of thermal parameters, the embedded depth, the length/radius ratio, the type of heat source and the stratification of the medium on the thermo-mechanical coupled performance.

     

The dilative intake of poroelastic inclusions an alternative to the Mandel–Cryer effect

This study examines the time-dependent poromechanics behavior of a fluid-saturated spherical inclusion embedded inside a fluid-saturated porous medium with different poroelastic properties. Both media comprise compressible constituents with distinctively defined poroelastic parameters. It is assumed that the inclusion is subjected to a fluid source at the center. The problem is formulated and solved using Biot theory of poromechanics. The contrasts in inclusion and the medium matrix stiffnesses and their respective hydraulic conductivities can be recognized as two competing factors, which affect the inclusion’s rate of expansion during fluid injection. Findings show a certain type of behavior that the inclusion exhibits at the onset of fluid injection when having greater stiffness than the medium matrix, where the inclusion experiences some decrease in the pore pressure. Compared to what announced as the stress redistribution due the Mandel–Cryer effect in earlier researches on dilation of free spheres, this study shows that the associated phenomenon would be likewise attributed to the coupled nature of pressures and deformations in the theory of poroelasticity. However, it is a consequence of the inclusion-matrix stiffness contrast where a dilating free sphere can be regarded as a special case of this new problem. The asymptotic expansions of pressure terms verify the existence of such an effect. The results of this study would put forward very good insight in some engineering applications.     

矩形移动荷载作用下路面-双层地基系统三维振动分析

考虑路面和地基之间的相互作用,建立了路面-双层地基的三维模型,将车辆荷载模拟成矩形移动荷载,利用Fourier变换方法对车辆荷载作用下路面-双层地基系统的三维振动问题进行了研究。假设地基为上部弹性土体和下部为饱和土体组成的双层混合结构,整个系统置于刚性基岩上。通过引入势函数,利用Lame分解理论和积分变换方法分别对弹性土层和饱和土层进行求解。在Fourier变换域内,联立路面和下卧双层地基系统的的运动方程,获得了车辆荷载作用下路面-双层地基系统三维振动的位移和孔压响应的积分形式解,并利用IFFT算法和自适应数值积分算法得到了数值计算结果。研究结果可为路基动力响应分析提供参考。     

Coupled formulation and algorithms for the simulation of non‐planar three‐dimensional hydraulic fractures using the generalized finite element method

This paper presents a coupled hydro‐mechanical formulation for the simulation of non‐planar three‐dimensional hydraulic fractures. Deformation in the rock is modeled using linear elasticity, and the lubrication theory is adopted for the fluid flow in the fracture. The governing equations of the fluid flow and elasticity and the subsequent discretization are fully coupled. A Generalized/eXtended Finite Element Method (G/XFEM) is adopted for the discretization of the coupled system of equations. A Newton–Raphson method is used to solve the resulting system of nonlinear equations. A discretization strategy for the fluid flow problem on non‐planar three‐dimensional surfaces and a computationally efficient strategy for handling time integration combined with mesh adaptivity are also presented. Several three‐dimensional numerical verification examples are solved. The examples illustrate the generality and accuracy of the proposed coupled formulation and discretization strategies. Copyright © 2015 John Wiley & Sons, Ltd.     

Algorithm of Coupled Normal Stress and Fluid Flow in Fractured Rock Mass by the Composite Element Method

This paper presents a composite element algorithm of coupled normal stress and fluid flow process for fractured rock mass, developed from the composite element method (CEM). The coupled relation between the fracture flow and normal stress makes use of the “filled model”, which examines the asperities in the fracture as a layer of granular medium having high porosity and being clipped by the two parallel plates. The existence of fractures is not considered in the mesh generation, but it will be considered explicitly in the mapped composite element. The coupled normal stress and fluid flow process has been simulated by applying a cross iterative algorithm between the two fields. The proposed algorithm considers not only the flow through the fractures, but also the flow exchange between fractures and the surrounding rock blocks. In addition, it can be used for both the filled and non-filled fractures. The verification of the proposed algorithm has been conducted through the illustration of three examples by comparison with the conventional finite element method (FEM), from which the advantages and reliability of the proposed algorithm have been shown clearly.     

Integration of a continuum damage model for shale with the cutting plane algorithm

The stability of integration is essential to numerical simulations especially when solving nonlinear problems. In this work, a continuum damage mechanics model proposed by the first author is implemented with an integration method named cutting plane algorithm (CPA) to improve the robustness of the simulation. This integration method is one type of return mapping algorithm that bypasses the need for computing the gradients. We compare the current integration method with the previous direct method, and the result shows that the cutting plane algorithm exhibits excellent performance under large loading rate conditions. To enhance accuracy of the new method, a control procedure is utilized in the implementation of the algorithm based on error analysis. Thereafter, the theory of poromechanics is utilized with the damage model to account for the effects of fluid diffusion. Laboratory tests simulated with finite element method illustrate distinct behaviors of shale with different loading rates and indicate the development of microcrack propagation under triaxial compression. Copyright © 2016 John Wiley & Sons, Ltd.     

基于热-流-固耦合模型的石油钻井施工过程数值分析

岩土介质多场耦合问题需考虑诸多因素,温度、渗流及应力之间的耦合关系复杂,试验条件不易控制,且难以实现,因此,辅以数值模拟手段具有重要的意义。基于混合物理论,推导出岩土介质温度、渗流和应力耦合的数学模型及其控制方程,提出该数学模型的求解方法,以MATLAB语言为平台,将Abaqus程序作为一个模块嵌入迭代算法程序中,编制了多场耦合分析程序,并给出了2个典型算例验证该方法的有效性和实用性。然后,将建立的多场耦合模型和计算程序应用于石油钻井施工过程的模拟,重点分析井壁围岩内温度场、渗流场和应力场的变化规律,以及钻井液温度的变化对井壁稳定性的影响。研究成果对我国地下石油、核废料储存等工程设计和施工具有一定的指导意义。     

Analysis of Pore Pressure and Stress Distribution around a Wellbore Drilled in Chemically Active Elastoplastic Formations

Drilling in low-permeable reactive shale formations with water-based drilling mud presents significant challenges, particularly in high-pressure and high-temperature environments. In previous studies, several models were proposed to describe the thermodynamic behaviour of shale. Most shale formations under high pressure are expected to undergo plastic deformation. An innovative algorithm including work hardening is proposed in the framework of thermo-chemo-poroelasticity to investigate the effect of plasticity on stresses around the wellbore. For this purpose a finite-element model of coupled thermo-chemo-poro-elastoplasticity is developed. The governing equations are based on the concept of thermodynamics of irreversible processes in discontinuous systems. In order to solve the plastic problem, a single-step backward Euler algorithm containing a yield surface-correction scheme is used to integrate the plastic stress–strain relation. An initial stress method is employed to solve the non-linearity of the plastic equation. In addition, super convergent patch recovery is used to accurately evaluate the time-dependent stress tensor from nodal displacement. The results of this study reveal that thermal and chemical osmosis can significantly affect the fluid flow in low-permeable shale formations. When the salinity of drilling mud is higher than that of pore fluid, fluid is pulled out of the formation by chemical osmotic back flow. Similar results are observed when the temperature of drilling mud is lower than that of the formation fluid. It is found that linear elastic approaches to wellbore stability analysis appear to overestimate the tangential stress around the wellbore and produce more conservative stresses compared to the results of field observation. Therefore, the drilling mud properties obtained from the elastoplastic wellbore stability in shales provide a safer mud weight window and reduce drilling cost.     

u‐p semi‐Lagrangian reproducing kernel formulation for landslide modeling

This paper presents a u‐p (displacement‐pressure) semi‐Lagrangian reproducing kernel (RK) formulation to effectively analyze landslide processes. The semi‐Lagrangian RK approximation is constructed based on Lagrangian discretization points with fixed kernel supports in the current configuration. As a result, it tracks state variables at discretization points while allowing extreme deformation and material separation that is beyond the capability of Lagrangian formulations. The u‐p formulation following Biot theory is incorporated into the formulation to describe poromechanics of saturated geomaterials. In addition, a stabilized nodal integration method to ensure stability of the domain integration and kernel contact algorithms to model contact between bodies are introduced in the u‐p semi‐Lagrangian RK formulation. The proposed method is verified with several numerical examples and validated with an experimental result and the field data of an actual landslide.     

Unsteady flow over offshore wind turbine airfoils and aerodynamic loads with computational fluid dynamic simulations

The first notable megawatt class wind turbine, which was the pioneer of improvement in the blade performance in large wind turbines, appeared in Vermont. Nowadays, modern wind turbines are using blades with multi-airfoils at different sections. In this study, in order to indicate the best airfoil profile for the optimum performance in different sections of a blade, five popular airfoils, including S8xx, FFA and AH series, were studied. On the large-scale profile, shear stress transport K–ω model was applied for the simulation of horizontal axis wind turbines for different wind speeds. The aerodynamic simulation was accomplished using computational fluid dynamic method, which in turn is based on the finite volume method, and semi-implicit method for pressure-linked equations algorithm is used for pressure–velocity coupling. The governing equations applied in this simulation are the unsteady Reynolds-averaged Navier–Stokes equations. The aerodynamic coefficients of lift and drag were calculated at different angle of attacks and different wind speeds. The results were validated by EPPLER code, XFOIL and experimental data of the US National Renewable Energy Laboratory. The results showed that S818 profile is the best profile in terms of gaining the highest lift coefficient with the lowest angle of attack at the root of the blades. The findings also indicated that the selected model can predict the exact geometry with a high precision.     

Accelerating the convergence of coupled geomechanical‐reservoir simulations

The pressure variations during the production of petroleum reservoir induce stress changes in and around the reservoir. Such changes of the stress state can induce marked deformation of geological structures for stress sensitive reservoirs as chalk or unconsolidated sand reservoirs. The compaction of those reservoirs during depletion affects the pressure field and so the reservoir productivity. Therefore, the evaluation of the geomechanical effects requires to solve in a coupling way the geomechanical problem and the reservoir multiphase fluid flow problem. In this paper, we formulate the coupled geomechanical‐reservoir problem as a non‐linear fixed point problem and improve the resolution of the coupling problem by comparing in terms of robustness and convergence different algorithms. We study two accelerated algorithms which are much more robust and faster than the conventional staggered algorithm and we conclude that they should be used for the iterative resolution of coupled reservoir‐geomechanical problem. Copyright © 2006 John Wiley & Sons, Ltd.