The open Delft Advanced Research Terra Simulator (open-DARTS) framework is an open-source reservoir simulation software. The open-DARTS focused on energy transition applications, such as geothermal energy production and carbon sequestration. It enables the modeling of compositional thermal flow, coupled with a geomechanical solver based on the Finite Volume discretization and adjoints method for inverse modeling. The open-DARTS supports different grid types (structured, corner-point geometry, and unstructured), discrete fracture networks, contact mechanics, and various thermal-chemical interactions. The recently proposed generic nonlinear formulation supports the most general nonlinear PDEs designed for various energy transition applications. The open-DARTS has been implemented in C++ and Python to optimize hardware utilization while ensuring flexibility. The most computationally expensive part is written in C++ and compiled into libraries, which are subsequently exposed to Python using pybind11. This allows the extension and overriding of C++ functions by user-defined Python code. For example, using only a Python interface, one can adjust a timestep strategy, nonlinear solver, or properties output. Besides, the Python interface of open-DARTS provides straightforward coupling with other Python-based numerical modeling packages, including the meshing, file storage, caching, and visualization modules. The open-DARTS core uses the advantages of C++ language, such as efficient low-level memory management, object-oriented programming, compile-time polymorphism, and parallelization with OpenMP. One of the advantages of open-DARTS is the Operator-Based Linearization (OBL) technique, which can resolve challenges associated with complex physics and reduce the computation time, especially for ensemble-based simulations. We would also like to share our experience on the project, repository, and the development workflow configuration using gitlab.com, including the build system (cmake), handling merge requests, automated testing in CI/CD pipelines, documentation management (gitlab.io), wiki utilization, and release publishing. Additionally, Python’s integration into open-DARTS offers the advantage of straightforward installation via PyPI and simplifies defining requirements for users who prefer to avoid compiling code from source files.

In this work, we present a unified framework for the simulation of CO₂ sequestration problems at various time and space scales. The parametrization technique utilizes thermodynamic state-dependent operators expressing the governing equations for the thermal-compositional-reactive system to solve the nonlinear problem. This approach provides flexibility in the assembly of the Jacobian, which allows straightforward implementation of advanced thermodynamics. We validate our simulation framework through several simulation studies including complex physical phenomena relevant to CCUS. The proposed simulation framework is validated against a set of numerical and experimental benchmark tests, demonstrating the efficiency and accuracy of the modeling framework for CCUS-related subsurface applications. Important physical phenomena resulting from the complex thermodynamic interactions of CO₂ and impurities with reservoir fluids can be accurately captured now in detailed dynamic simulation. The investigated simulation scenarios include a reproduction of lab experiments at the core scale, investigation of macro-scale analog model and simulation of large-scale industrial application. The simulation time can also span from hours to years among various applications. Complex thermal-compositional-reactive phenomena can be addressed at each of these space and time scales. The unified thermodynamic description allows us to perform all these simulations for a reasonable CPU time due to advanced parametrization techniques and efficient GPU capabilities in our in-house reservoir simulator DARTS.