Constraint conditions for the partial saturations and the pressure fractions of CO2 and brine are defined.
To characterize the stress state in the solid matrix, the effective stress principle is applied.
The coupled problem is solved using the in-house scientific code OpenGeoSys (an open source finite element code) and verified with benchmarks.