Compressible single phase fluid model

Overview

This model represents a compressible single-phase fluid with constant compressibility and pressure-dependent viscosity. These assumptions are valid for slightly compressible fluids, such as water, and some types of oil with negligible amounts of dissolved gas.

Specifically, fluid density is computed as

$\rho(p) = \rho_0 e^{c_\rho(p - p_0)}$

where $c_\rho$ is compressibility, $p_0$ is reference pressure, $\rho_0$ is density at reference pressure. Similarly,

$\mu(p) = \mu_0 e^{c_\mu(p - p_0)}$

where $c_\mu$ is viscosibility (viscosity compressibility), $\mu_0$ is reference viscosity.

Either exponent may be approximated by linear (default) or quadratic terms of Taylor series expansion. Currently there is no temperature dependence in the model, although it may be added in future.

Parameters

The model is represented by <CompressibleSinglePhaseFluid> node in the input.

The following attributes are supported:

Name	Type	Default	Description
compressibility	real64	0	Fluid compressibility
defaultDensity	real64	required	Default value for density.
defaultViscosity	real64	required	Default value for viscosity.
densityModelType	geos_constitutive_ExponentApproximationType	linear	Type of density model. Valid options: * exponential * linear * quadratic
name	groupName	required	A name is required for any non-unique nodes
referenceDensity	real64	1000	Reference fluid density
referencePressure	real64	0	Reference pressure
referenceViscosity	real64	0.001	Reference fluid viscosity
viscosibility	real64	0	Fluid viscosity exponential coefficient
viscosityModelType	geos_constitutive_ExponentApproximationType	linear	Type of viscosity model. Valid options: * exponential * linear * quadratic

Example

<Constitutive>
  <CompressibleSinglePhaseFluid name="water"
                                referencePressure="2.125e6"
                                referenceDensity="1000"
                                compressibility="1e-19"
                                referenceViscosity="0.001"
                                viscosibility="0.0"/>
</Constitutive>