Abstract
CO2 storage and hydrocarbon recovery is modelled by mixed hyperbolic-elliptic PDEs, posed on unstructured, polyhedral grids with large aspect ratios. Such models are strongly nonlinear, have a truly multiscale behavior, and contain complex tabulated data or analytic submodels that makes linearization and solution nontrivial. Herein, I discuss challenges faced when developing new (multiscale) simulation methods for industry use. To overcome these challenges, we use the open-source MATLAB Reservoir Simulation Toolbox for rapid prototyping and bootstrapping of new methods, and to leverage the effort of others in a reproducible manner.