Flexible and rigorous numerical modelling of multiphysics processes in fractured porous media using PorePy
The paper argues the importance of two central features of simulation tools for research into mathematical and numerical models. The requirements are rigour, ensuring that the governing equations are accurately solved, and flexibility, allowing experimentation with different physical models and solution techniques.
Main content
Corresponding author: Ivar Stefansson (ivar.stefansson@uib.no)
Complementing the general discussion on code design, we present how the features manifest in the PorePy toolbox in the context of coupled processes in fractured porous media. We start out from high-level design principles including structured testing and a high level of modularity. Then, we provide more detail in a suite of examples illustrated by code snippets. The examples include samples from the PorePy test suite as well as simulations highlighting modelling flexibility and demonstrating how to achieve complicated, tailored models with modest implementational effort.
Flexible and rigorous numerical modelling of multiphysics processes in fractured porous media using PorePy
https://doi.org/10.1016/j.rinam.2023.100428
Ivar Stefansson, Jhabriel Varela, Eirik Keilegavlen, Inga Berre