BZPEERReading, trust and collaboration for researchDiscover
Workspaces
Network
Opportunities
Account
Researcher profile
Beatrice Riviere contributes to research discovery and scholarly infrastructure.
Trust snapshot
Actions
Identity and collaboration
Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.
Log in to claimDirect collaboration
Claim this author entity first to unlock direct invitations.
Research graph
Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.
BZPEER is loading the nearby papers, people, topics and institutions for this page.
Published work
We formulate a numerical method for solving the two-phase flow poroelasticity equations. The scheme employs the interior penalty discontinuous Galerkin method and a sequential time-stepping method. The unknowns are the phase pressures and the displacement. Existence of the solution is proved. Three-dimensional numerical results show the accuracy and robustness of the proposed method.
The analyses of interior penalty discontinuous Galerkin methods of any order k for solving elliptic and parabolic problems with Dirac line sources are presented. For the steady state case, we prove convergence of the method by deriving a priori error estimates in the L2 norm and in weighted energy norms. In addition, we prove almost optimal local error estimates in the energy norm for any approximation order. Further, almost optimal local error estimates in the L2 norm are obtained for the case of piecewise linear approximations whereas suboptimal error bounds in the L2 norm are shown for any polynomial degree. For the time-dependent case, convergence of semi-discrete and of backward Euler fully discrete scheme is established by proving error estimates in L2 in time and in space. Numerical results for the elliptic problem are added to support the theoretical results.
This paper presents and analyzes a discontinuous Galerkin method for the incompressible three-phase flow problem in porous media. We use a first order time extrapolation which allows us to solve the equations implicitly and sequentially. We show that the discrete problem is well-posed, and obtain a priori error estimates. Our numerical results validate the theoretical results, i.e. the algorithm converges with first order.
A numerical method is formulated for the solution of the advective Cahn-Hilliard (CH) equation with constant and degenerate mobility in three-dimensional porous media with non-vanishing velocity on the exterior boundary. The CH equation describes phase separation of an immiscible binary mixture at constant temperature in the presence of a mass constraint and dissipation of free energy. Porous media/pore-scale problems specifically entail high-resolution images of rocks in which the solid matrix and pore spaces are fully resolved. The interior penalty discontinuous Galerkin method is used for the spatial discretization of the CH equation in mixed form, while a semi-implicit convex-concave splitting is utilized for temporal discretization. The spatial approximation order is arbitrary, while it reduces to a finite volume scheme for the choice of elementwise constants. The resulting nonlinear systems of equations are reduced using the Schur complement and solved via Newton's method. The numerical scheme is first validated using numerical convergence tests and then applied to a number of fundamental problems for validation and numerical experimentation purposes including the case of degenerate mobility. First-order physical applicability and robustness of the numerical method are shown in a breakthrough scenario on a voxel set obtained from a micro-CT scan of a real sandstone rock sample.
One-dimensional blood flow models take the general form of nonlinear hyperbolic systems but differ greatly in their formulation. One class of models considers the physically conserved quantities of mass and momentum, while another class describes mass and velocity. Further, the averaging process employed in the model derivation requires the specification of the axial velocity profile; this choice differentiates models within each class. Discrepancies among differing models have yet to be investigated. In this paper, we systematically compare several reduced models of blood flow for physiologically relevant vessel parameters, network topology, and boundary data. The models are discretized by a class of Runge-Kutta discontinuous Galerkin methods.
Mathematical modeling at the level of the full cardiovascular system requires the numerical approximation of solutions to a one-dimensional nonlinear hyperbolic system describing flow in a single vessel. This model is often simulated by computationally intensive methods like finite elements and discontinuous Galerkin, while some recent applications require more efficient approaches (e.g. for real-time clinical decision support, phenomena occurring over multiple cardiac cycles, iterative solutions to optimization/inverse problems, and uncertainty quantification). Further, the high speed of pressure waves in blood vessels greatly restricts the time step needed for stability in explicit schemes. We address both cost and stability by presenting an efficient and unconditionally stable method for approximating solutions to diagonal nonlinear hyperbolic systems. Theoretical analysis of the algorithm is given along with a comparison of our method to a discontinuous Galerkin implementation. Lastly, we demonstrate the utility of the proposed method by implementing it on small and large arterial networks of vessels whose elastic and geometrical parameters are physiologically relevant.