viscosity {marelac} | R Documentation |
Calculates the shear viscosity of water, in centipoise (g/m/sec). Valid for 0 < t < 30 degrees C, 0 < S < 36, 1 < P < 1000 bars.
Based on the code "CANDI" by B.P. Boudreau
viscosity(S = 35, t = 25, P = 1.013253)
S |
salinity, -, |
t |
temperature, degrees C, |
P |
True pressure, bar. |
The details given in the original code by B. Boudreau are repeated here:
Uses the equation given by Kukulka et al. (1987).
Shear visocisity in centipoise.
Karline Soetaert <karline.soetaert@nioz.nl>
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.
coriolis
,
diffcoeff
, ssd2rad
, vertmean
,
gravity
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"))