On the upstream mobility scheme for two-phase flow in porous media

When neglecting capillarity, two-phase incompressible flow in porous media is modelled as a scalar nonlinear hyperbolic conservation law. A change in the rock type results in a change of the flux function. Discretizing in one-dimensional with a finite volume method, we investigate two numerical fluxes, an extension of the Godunov flux and the upstream mobility flux, the latter being widely used in hydrogeology and petroleum engineering. Then, in the case of a changing rock type, one can give examples when the upstream mobility flux does not give the right answer.Comment: A preprint to be published in Computational Geoscience

A multiscale flux basis for mortar mixed discretizations of reduced Darcy-Forchheimer fracture models

In this paper, a multiscale flux basis algorithm is developed to efficiently solve a flow problem in fractured porous media. Here, we take into account a mixed-dimensional setting of the discrete fracture matrix model, where the fracture network is represented as lower-dimensional object. We assume the linear Darcy model in the rock matrix and the non-linear Forchheimer model in the fractures. In our formulation, we are able to reformulate the matrix-fracture problem to only the fracture network problem and, therefore, significantly reduce the computational cost. The resulting problem is then a non-linear interface problem that can be solved using a fixed-point or Newton-Krylov methods, which in each iteration require several solves of Robin problems in the surrounding rock matrices. To achieve this, the flux exchange (a linear Robin-to-Neumann co-dimensional mapping) between the porous medium and the fracture network is done offline by pre-computing a multiscale flux basis that consists of the flux response from each degree of freedom on the fracture network. This delivers a conserve for the basis that handles the solutions in the rock matrices for each degree of freedom in the fractures pressure space. Then, any Robin sub-domain problems are replaced by linear combinations of the multiscale flux basis during the interface iteration. The proposed approach is, thus, agnostic to the physical model in the fracture network. Numerical experiments demonstrate the computational gains of pre-computing the flux exchange between the porous medium and the fracture network against standard non-linear domain decomposition approaches

Eruptive Event Generator Based on the Gibson-Low Magnetic Configuration

Coronal Mass Ejections (CMEs), a kind of energetic solar eruptions, are an integral subject of space weather research. Numerical magnetohydrodynamic (MHD) modeling, which requires powerful computational resources, is one of the primary means of studying the phenomenon. With increasing accessibility of such resources, grows the demand for user-friendly tools that would facilitate the process of simulating CMEs for scientific and operational purposes. The Eruptive Event Generator based on Gibson-Low flux rope (EEGGL), a new publicly available computational model presented in this paper, is an effort to meet this demand. EEGGL allows one to compute the parameters of a model flux rope driving a CME via an intuitive graphical user interface (GUI). We provide a brief overview of the physical principles behind EEGGL and its functionality. Ways towards future improvements of the tool are outlined