# viscosity: Shear Viscosity of Water In marelac: Tools for Aquatic Sciences

 viscosity R Documentation

## 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)
``````

### Arguments

 `S ` salinity, -, `t ` temperature, `^\circ`C, `P ` True pressure, bar.

### Details

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

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

### Value

Shear visocisity in centipoise.

### Author(s)

Karline Soetaert <karline.soetaert@nioz.nl>

### 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.

`coriolis`, `diffcoeff`, `ssd2rad`, `vertmean`, `gravity`

### Examples

``````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"))
``````

marelac documentation built on Sept. 25, 2023, 5:06 p.m.