Learn R Programming

marelac (version 2.1.11)

viscosity: Shear Viscosity of Water

Description

Calculates the shear viscosity of water, in centipoise (g/m/sec). Valid for 0 < t < 30 \(^\circ\)C, 0 < S < 36, 1 < P < 1000 bars.

Based on the code "CANDI" by B.P. Boudreau

Usage

viscosity(S = 35, t = 25, P = 1.013253)

Value

Shear visocisity in centipoise.

Arguments

S

salinity, -,

t

temperature, \(^\circ\)C,

P

True pressure, bar.

Author

Karline Soetaert <karline.soetaert@nioz.nl>

Details

The details given in the original code by B. Boudreau are repeated here:

Uses the equation given by Kukulka et al. (1987).

References

Based on the FORTRAN implementation of the diagenetic model "CANDI" of B.P. Boudreau:

Boudreau BP, 1996. A method-of-lines code for carbon and nutrient diagenesis in aquatic sediments. Computers & Geosciences 22 (5), 479-496.

Kulkula DJ, Gebhart B and Mollendorf JC, 1987. Thermodynamic and transport properties of pure and saline water. Adv. Heat transfer 18, 325-363.

See Also

coriolis, diffcoeff, ssd2rad, vertmean, gravity

Examples

Run this code
plot(0:30, viscosity(t = 0:30, S = 35, P = 1),
      xlab = "temperature", ylab = "g/m/s",
      main = "shear viscosity of water", type = "l")
lines(0:30, viscosity(t = 0:30, S = 0, P = 1), col = "red")
lines(0:30, viscosity(t = 0:30, S = 35, P = 100), col = "blue")
legend("topright", col = c("black","red","blue"), lty = 1,
        legend = c("S=35, P=1", "S=0, P=1", "S=35, P=100"))

Run the code above in your browser using DataLab