Numerical Investigation on the Fixed-Stress Splitting Scheme for Biot’s Equations: Optimality of the Tuning Parameter

We study the numerical solution of the quasi-static linear Biot equations solved iteratively by the fixed-stress splitting scheme. In each iteration the mechanical and flow problems are decoupled, where the flow problem is solved by keeping an artificial mean stress fixed. This introduces a numerical tuning parameter which can be optimized. We investigate numerically the optimality of the parameter and compare our results with physically and mathematically motivated values from the literature, which commonly only depend on mechanical material parameters. We demonstrate, that the optimal value of the tuning parameter is also affected by the boundary conditions and material parameters associated to the fluid flow problem suggesting the need for the integration of those in further mathematical analyses optimizing the tuning parameter.acceptedVersio

Multigrid particle-in-cell simulations of plasma microturbulence

A new scheme to accurately retain kinetic electron effects in particle-in-cell (PIC) simulations for the case of electrostatic drift waves is presented. The splitting scheme, which is based on exact separation between adiabatic and on adiabatic electron responses, is shown to yield more accurate linear growth rates than the standard df scheme. The linear and nonlinear elliptic problems that arise in the splitting scheme are solved using a multi-grid solver. The multi-grid particle-in-cell approach offers an attractive path, both from the physics and numerical points of view, to simulate kinetic electron dynamics in global toroidal plasmas

The role of rock joint frictional strength in the containment of fracture propagation

The fracturing phenomenon within the reservoir environment is a complex process that is controlled by several factors and may occur either naturally or by artificial drivers. Even when deliberately induced, the fracturing behaviour is greatly influenced by the subsurface architecture and existing features. The presence of discontinuities such as joints, artificial and naturally occurring faults and interfaces between rock layers and microfractures plays an important role in the fracturing process and has been known to significantly alter the course of fracture growth. In this paper, an important property (joint friction) that governs the shear behaviour of discontinuities is considered. The applied numerical procedure entails the implementation of the discrete element method to enable a more dynamic monitoring of the fracturing process, where the joint frictional property is considered in isolation. Whereas fracture propagation is constrained by joints of low frictional resistance, in non-frictional joints, the unrestricted sliding of the joint plane increases the tendency for reinitiation and proliferation of fractures at other locations. The ability of a frictional joint to suppress fracture growth decreases as the frictional resistance increases; however, this phenomenon exacerbates the influence of other factors including in situ stresses and overburden conditions. The effect of the joint frictional property is not limited to the strength of rock formations; it also impacts on fracturing processes, which could be particularly evident in jointed rock masses or formations with prominent faults and/or discontinuities