STOMP Subsurface Transport Over Multiple Phases, Version 4.0, User?s Guide

Abstract

This guide describes the general use, input file formatting, compilation and execution of the STOMP (Subsurface Transport Over Multiple Phases) simulator, a scientific tool for analyzing single and multiple phase subsurface flow and transport. A description of the simulator’s governing equations, constitutive functions and numerical solution algorithms are provided in a companion theory guide. In writing these guides for the STOMP simulator, the authors have assumed that the reader comprehends concepts and theories associated with multiple-phase hydrology, heat transfer, thermodynamics, radioactive chain decay, and relative permeability-saturation-capillary pressure constitutive relations. The authors further assume that the reader is familiar with the computing environment on which they plan to compile and execute the STOMP simulator. Source codes for the sequential versions of the simulator are available in pure FORTRAN 77 or mixed FORTRAN 77/90 forms. The pure FORTRAN 77 source code form requires a parameters file to define the memory requirements for the array elements. The mixed FORTRAN 77/90 form of the source code uses dynamic memory allocation to define memory requirements, based on a FORTRAN 90 preprocessor STEP, that reads the input files. The simulator utilizes a variable source code configuration, which allows the execution memory and speed to be tailored to the problem specifics, and essentially requires that the source code be assembled and compiled through a software maintenance utility. The memory requirements for executing the simulator are dependent on the complexity of physical system to be modeled and the size and dimensionality of the computational domain. Likewise execution speed depends on the problem complexity, size and dimensionality of the computational domain, and computer performance. Selected operational modes of the STOMP simulator are available for scalable execution on multiple processor (i.e., parallel) computers. These versions of the simulator are written in pure FORTRAN 90 with imbedded directives that are interpreted by a FORTRAN preprocessor. Without the preprocessor, the scalable version of the simulator can be executed sequentially on a single processor computer. The scalable versions of the STOMP modes carry the “-Sc” designator on the operational mode name. For example, STOMP-WCS-Sc is the scalable version of the STOMP-WCS (Water-CO2-Salt) mode. A separate mode containing an evaporation model as a boundary condition on the upper surface of the computation domain has also been included. This mode, STOMP-WAE-B (Water-Air-Energy-Barriers) can be viewed as an extension of the STOMP-WAE (Water-Air-Energy) mode. Details of this particular mode are outlined by Ward et al. (2005)(a). STOMP V4.0 includes the reactive transport module ECKEChem (Equilibrium-Conservation-Kinetic Equation Chemistry) for the STOMP-W (Water) and STOMP-WCS (Water-CO2-Salt) modes. For this particular module, the “-R” designator is included in the operational mode name (e.g., STOMP-W-R, STOMP-WCS-R-Sc). This mode is described in detail by White and McGrail (2005)(b). For all operational modes and processor implementations, the memory requirements for executing the simulator are dependent on the complexity of physical system to be modeled and the size and dimensionality of the computational domain. Likewise execution speed depends on the problem complexity, size and dimensionality of the computational domain, and computer performance. Additional information about the simulator can be found on the STOMP webpage: http://stomp.pnl.gov. The website includes an introductory short course with problems ranging from simple one-dimensional saturated flow to complex multiphase system computations.