Accurate and efficient simulation of multiphase flow in heterogeneous porous media motivates the development of space-time multiscale strategies for the coupled nonlinear flow (pressure) and saturation transport equations. The flow equation entails heterogeneous high-resolution (fine-scale) coefficients and is global (elliptic or parabolic). The time-dependent saturation profile, on the other hand, may exhibit sharp local gradients or discontinuities (fronts) where the solution accuracy is highly sensitive to the time-step size. Therefore, accurate flow solvers need to address the multiscale spatial scales, while advanced transport solvers need to also tackle multiple time scales. This paper presents the first multirate multiscale method for space-time conservative multiscale simulation of sequentially coupled flow and transport equations. The method computes the pressure equation at the coarse spatial scale with a multiscale finite volume technique, while the transport equation is solved by taking variable time-step sizes at different locations of the domain. At each coarse time step, the developed local time-stepping technique employs an adaptive recursive time step refinement to capture the fronts accurately. The applicability (accuracy and efficiency) of the method is investigated for a wide range of two-phase flow simulations in heterogeneous porous media. For the studied cases, the proposed method is found to provide 3 to 4 times faster simulations. Therefore, it provides a promising strategy to minimize the tradeoff between accuracy and efficiency for field-scale applications.

@en

We present ADM-LTS, an adaptive multilevel space-time-stepping scheme for transport in heterogeneous porous media. At each time step, firstly, the flow (pressure) solution is obtained. Then, the transport equation is solved using the ADM-LTS method, which consists of two stages. In the first stage, an initial solution is obtained by imposing the coarsest space-time grid. This initial solution is then improved, in the second stage, by imposing a space-time adaptive grid on the cells where the solution does not satisfy the desired quality. The quality control is based on error estimators with user-defined threshold values. The time-integration procedure, in which the coarsest-scale solution provides local flux boundary conditions for sub-domains with local time refinement, is strictly mass conservative. In addition, the method employs space-time fine grid cells only at the moving saturation fronts. In order to ensure local mass conservation at all levels, finite-volume restriction operators and unity prolongation operators are developed. Several numerical experiments have been performed to analyze the efficiency and accuracy of the proposed ADM-LTS method for both homogeneous and heterogeneous permeability fields on two and three dimensional domains. The results show that the method provides accurate solutions, at the same time it maintains the computational efficiency. The ADM-LTS implementation is publicly available at https://gitlab.com/darsim2simulator.@en

We propose a mathematical model and a discretization strategy for the simulation of pressurized fractures in porous media accounting for the poroelastic effects due to the interaction of pressure and flow with rock deformations. The aim of the work is to develop a numerical scheme suitable to model the interplay among several fractures subject to fluid injection in different geometric configurations, in view of the application of this technique to hydraulic fracturing. The eXtended Finite Element Method, here employed for both the mechanical and fluid-dynamic problems, is particularly useful to analyze different configurations without remeshing. In particular, we adopt an ad hoc enrichment for the displacement at the fracture tip and a hybrid dimensional approach for the fluid. After the presentation of the model and discretization details we discuss some test cases to assess the impact of fracture spacing on aperture during injection.

@en

The present work proposes a novel method for the simulation of crack propagation in brittle elastic materials that combines two of the most popular approaches in literature. A large scale displacement solution is obtained with the well known extended finite elements method (XFEM), while propagation is governed by the solution of a local phase field problem at the tip scale. The method, which we will refer to as Xfield, is here introduced and tested in 2D under mixed modes I and II loads. The main features and the capability of the Xfield to efficiently simulate crack propagation are shown in some numerical benchmarks.

@en

The process by which rocks are formed from the burial of a fresh sediment involves the coupled effects of mechanical compaction and geochemical reactions. Both of them affect the porosity and permeability of the rock and, in particular, geochemical reactions can significantly alter them, since dissolution and precipitation processes may cause a structural transformation of the solid matrix. Often, the differential problems that arise from the modeling of these chemical reactions may present a discontinuous right hand side, where the discontinuity depends on the solution itself. In this work we have developed a numerical model to simulate this complex multi-physics problem by treating the discontinuous right hand side with a specially tailored event-driven numerical scheme. We show the performance of this strategy in terms of positivity and mass conservation, also in comparison with a more traditional approach that relies on a regularization of the discontinuous terms.@en

The present work deals with the numerical simulation of porous media subject to the coupled effects of mechanical compaction and reactive flows that can significantly alter the porosity due to dissolution, precipitation or transformation of the solid matrix. These chemical processes can be effectively modelled as ODEs with discontinuous right hand side, where the discontinuity depends on time and on the solution itself. Filippov theory can be applied to prove existence and to determine the solution behaviour at the discontinuities. From the numerical point of view, tailored numerical schemes are needed to guarantee positivity, mass conservation and accuracy. In particular, we rely on an event-driven approach such that, if the trajectory crosses a discontinuity, the transition point is localized exactly and integration is restarted accordingly.@en

We propose a mathematical model and a numerical scheme to describe compaction processes in a sedimentary rock layer undergoing both mechanical and geochemical processes. We simulate the sedimentation process by providing a sedimentation rate, and we account for chemical reactions using simplified kinetics describing either the conversion of a solid matrix into a fluid, as in the case of kerogen degradation into oil, or the precipitation of a mineral solute on the solid matrix of the rock. We use a Lagrangian description that enables to recast the equations in a fixed frame of reference. We present an iterative splitting scheme that allows solving the set of governing equations efficiently in a sequential manner. We assess the performances of this strategy in terms of convergence and mass conservation. Some numerical experiments show the capability of the scheme to treat two test cases, one concerning the precipitation of a mineral, the other the dissolution of kerogen.@en