It is occasionally useful to be able to calculate thermodynamic properties of the fluids that FEHM simulates. For convenience, the functions that calculate these properties, for twophase water and CO2, are supplied here. Density (dens()), enthalpy (enth()) and viscosity (visc()), and their corresponding derivatives with respect to pressure and temperature, can be calculated for given pressure and temperature conditions. Furthermore, the position and derivatives of the saturation line for water can be calculated for a given temperature or pressure (sat() and tsat()).
Note, these functions are vectorised; that is, the user can supply vectors of temperature and pressure to output a vector of the corresponding variable.
Return liquid water, vapor water and CO2 density, or derivatives with respect to temperature or pressure, for specified temperature and pressure.
Parameters: 


Returns:  Three element tuple containing (liquid, vapor, CO2) density or derivatives if requested. 
Return liquid water, vapor water and CO2 enthalpy, or derivatives with respect to temperature or pressure, for specified temperature and pressure.
Parameters: 


Returns:  Three element tuple containing (liquid, vapor, CO2) enthalpy or derivatives if requested. 
Return liquid water, vapor water and CO2 viscosity, or derivatives with respect to temperature or pressure, for specified temperature and pressure.
Parameters: 


Returns:  Three element tuple containing (liquid, vapor, CO2) viscosity or derivatives if requested. 
Return saturation pressure and first derivative for given temperature.
Parameters:  T (fl64) – Temperature (degC) 

Returns:  Two element tuple containing (saturation pressure, derivative). 
Return saturation temperature and first derivative for given pressure.
Parameters:  P (fl64) – Pressure (degC) 

Returns:  Two element tuple containing (saturation temperature, derivative). 
Calculate thermodynamic properties of a column of fluid.
Parameters: 


Returns:  Three element tuple containing (liquid, vapor, CO2) properties. Each contains a three column array corresponding to pressure, temperature, density, enthalpy and viscosity of the fluid. 
z=np.linspace(0,4e3,200)
Pco2=fluid_column(z,0.025,20,8)[2][:,0]