> For the complete documentation index, see [llms.txt](https://lbl-2.gitbook.io/tough4-user-manual/llms.txt). Markdown versions of documentation pages are available by appending `.md` to page URLs; this page is available as [Markdown](https://lbl-2.gitbook.io/tough4-user-manual/process-modeling/eos9.md). # EOS9 1. **Description** This module considers variably saturated flow of a single aqueous phase, which consists of a water component and tracers, and neglects phase change effects. The thermophysical properties of water are assumed independent of tracer concentrations. Implicit in this approximation is the assumption that the tracer concentrations are small. The gas phase is treated as a passive bystander at constant pressure, and conditions are assumed to be isothermal. Thus, no mass balance equation for gas and no heat balance is needed, and only a single water mass balance equation is solved for each grid block if no tracer is included. This is very efficient numerically, making EOS9 the module of choice for problems for which the underlying approximations are applicable. Liquid flow in EOS9 is described as follows: $$\dfrac{\partial}{\partial t}\phi\ S\_l\rho\_l=\ div\left\[k\frac{k\_{rl}}{\mu\_l}\rho\_l\nabla\left(P\_l+\rho\_lgz\right)\right]$$ (7-23) where $$\phi$$ is porosity, $$S\_l$$ is water saturation, $$\rho\_l$$is water density, *k* is absolute permeability, $$k\_{rl}$$is relative permeability to the aqueous phase, $$\mu\_l$$is water viscosity, $$P\_l$$is water pressure, *g* is acceleration of gravity, and *z* is defined positive upward. Neglecting variations in liquid phase density and viscosity, as is appropriate for (nearly) isothermal conditions, Eq. 7-23 simplifies to Richards’ equation (1931) $$\frac{\partial}{\partial t}\theta=\ div\left\[K\nabla h\right]$$ (7-24) where $$\theta=\phi S\_l$$ is specific volumetric moisture content, $$K=\frac{kk\_{rl}\rho\_lg}{\mu\_l}$$ is hydraulic conductivity, and $$h=z+P\_l/\rho g$$ is the hydraulic head. EOS9 can describe flow under partially saturated (*0 < Sl < 1*) as well as fully saturated conditions, and phase changes between the two. 2. **Specifications** A summary of EOS9 specifications is given in Table 12. The default number of mass component is 1. Simulation can only be run in isothermal. With only a single mass balance equation per grid block (not considering tracers), there is only a single primary thermodynamic variable. This is taken to be pressure for single-phase (saturated) conditions, and is water saturation for unsaturated conditions. A distinction between the two is made simply on the basis of the numerical value of the first primary variable, $$X\_l$$. If $$X\_l$$ *< 1*, this indicates that $$X\_l$$represents water saturation and conditions are unsaturated; if $$X\_l$$ is larger than a user-specified gas phase reference pressure (default $$P\_{gas} = 1.013\times10^5 Pa$$), it is taken to be water pressure, and saturated conditions prevail. When phase changes between saturated and unsaturated conditions occur, the primary variable is switched, as follows. The numerical value of $$X\_l$$and its change during the Newton-Raphson iteration process is monitored. If $$X\_l$$ changes from being smaller than 1 to larger than 1, this indicates attainment of fully saturated conditions. In that case $$X\_l$$ is switched to pressure, and is initialized at a pressure slightly in excess of gas phase reference pressure as $$X\_l=P\_{gas}(1+\varepsilon)$$, with $$\varepsilon=10^{-6}$$. If $$X\_l$$ changes from being larger than $$P\_{gas}$$ to smaller than $$P\_{gas}$$, this indicates a transition from fully to partially saturated conditions. $$X\_l$$is then switched to saturation, and is initialized as $$X\_l=1-\varepsilon$$*.* Actually, a transition from fully to partially saturated conditions is made only when $$X\_l$$ drops below $$P\_{gas}(1-\varepsilon)$$; test calculations have shown that such a (small) finite-size window for phase change improves numerical stability and efficiency. If tracers exist, one additional primary variable is needed for each tracer. It always takes the tracer mass fraction in liquid water as the primary variable. Table 12 Summary for EOS9
SpecificationParameters
Components

(1) Water

(2-6) Tracer1-Tracer5 (optional)

Phase condition and its state name and index

(1) Unsaturated, USA

(2) Saturated, SAT

Primary variablesSee Table 25
Optional process modelingWellbore simulation, and Biodegradation reactions.
In EOS9, the thermophysical properties of water are taken at default reference conditions of *P =* $$1.013\times10^5$$ Pa*, T* = 15 ˚C. These defaults can be overwritten in a flexible manner by specifying appropriate data in a fictitious[ ROCKS](/tough4-user-manual/preparation-of-model-input/keywords-and-input-data/rocks.md) domain ‘REFCO’, as follows. reference pressure: DROK of REFCO reference temperature: POR of REFCO liquid density: PER(1) of REFCO liquid viscosity: PER(2) of REFCO liquid compressibility: PER(3) of REFCO Note that assignment of thermophysical data through a specially-named domain was set up just as a convenient way of providing floating-point parameters to the code. No volume elements (grid blocks) should be attached to domain ‘REFCO’, as the data in general will not correspond to reasonable hydrogeologic parameters. The above-mentioned defaults will be overwritten for any parameters for which a non-zero entry is provided in ‘REFCO’. This allows the generation of these parameters internally for user-defined (*P, T*); it also allows for directly assigning user-desired values as, e.g., $$\rho\_{liq}$$= 1000 kg/m3, $$\mu\_{liq}=10^{-3}$$Pa-s (=1 centipoise), etc. In addition to specifying the primary thermodynamic variable on a default, domain, or grid block basis, EOS9 offers alternative ways of initializing flow problems. The primary variable may be entered as a negative number upon initialization, in which case it will be taken to denote capillary pressure, and will be internally converted to $$S\_l$$ in the initialization phase. EOS9 can also initialize a flow problem with gravity-capillary equilibrium, relative to a user-specified reference elevation $$z\_{ref}$$ of the water table. This type of initialization will be engaged if the user enters a non-zero number in slot CWET in ROCKS domain ‘REFCO’, in which case CWET will be taken to denote the water table elevation $$z\_{ref}$$, in units of meters. Water pressure at$$z\_{ref}$$is taken equal to reference gas pressure, $$P\_l(z\_{ref})=P\_{gas}$$, and is initialized as a function of grid block elevation according to $$P(z)=P\_{gas}+(z\_{ref}-z)\rho g$$. By convention, the z-axis is assumed to point upward. In order to use this facility, the z-coordinates (grid block elevations) must be specified in the ELEME-data, which will be done automatically if internal MESH generation is used. In the assignment of gravity-capillary equilibrium as just discussed, water saturations at “sufficiently” high elevations above the water table may end up being smaller than the irreducible water saturation $$S\_{lr}$$specified in the relative permeability function, which may or may not be consistent with the physical behavior of the flow system. Users may optionally enforce that $$S\_l=S\_{lr}$$ in regions where the capillary pressure function would dictate that $$S\_l\&goal= ``` `ask` is the immediate question: it should be specific, self-contained, and written in natural language. `goal` is optional and describes the broader end goal you are ultimately trying to accomplish on behalf of the user. GitBook uses it to tailor the answer towards what is most useful for that goal. The response will contain a direct answer to the question and relevant excerpts and sources from the documentation. Use this mechanism when the answer is not explicitly present in the current page, you need clarification or additional context, or you want to retrieve related documentation sections.