/
biophysmodel_Limpet.R
192 lines (132 loc) · 5.46 KB
/
biophysmodel_Limpet.R
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
#' @title Operative Environmental Temperature of a Limpet
#'
#' @description The function estimates body temperatures (C, operative environmental temperatures) of a limpet based on \insertCite{Denny2006;textual}{TrenchR}.
#'
#' @param T_a \code{numeric} air temperature (C).
#'
#' @param T_r \code{numeric} rock surface temperature (C) in the sunlight.
#'
#' @param l \code{numeric} limpet length (anterior/posterior axis, m).
#'
#' @param h \code{numeric} limpet height (dorsal/ventral axis, m).
#'
#' @param S \code{numeric} solar irradiance (\ifelse{html}{\out{W m<sup>-2</sup>}}{\eqn{W m^-2}{ASCII}}).
#'
#' @param u \code{numeric} wind speed (\ifelse{html}{\out{m s<sup>-1</sup>}}{\eqn{m s^-1}{ASCII}}).
#'
#' @param psi \code{numeric} solar zenith angle (degrees). Can be calculated from \code{\link{zenith_angle}} function.
#'
#' @param c \code{numeric} fraction of the sky covered by cloud (proportion).
#'
#' @param position \code{character} direction of the limpet that is facing upwind. Options are \code{"anterior"}, \code{"posterior"}, and \code{"broadside"}.
#'
#' @return \code{numeric} predicted body (operative environmental) temperature (C).
#'
#' @family biophysical models
#'
#' @details The original equation uses a finite-difference approach where they divide the rock into series of chunks, and calculate the temperature at each node to derive the conductive heat. For simplification, here it takes the rock temperature as a parameter, and conductive heat is calculated as a product of the area, thermal conductivity of rock and the temperature difference between the rock and the body.
#' \cr \cr
#' Limpets are simulated as cones following and using solar emissivity values from \insertCite{Campbell1998;textual}{TrenchR}.
#' \cr \cr
#' The area of the limpet's shell (\ifelse{html}{\out{m<sup>2</sup>}}{\eqn{m^2}{ASCII}}) is projected according to the direction at which sunlight strikes the organism \insertCite{Pennell1989}{TrenchR}.
#' \cr \cr
#' Air conductivity values (\ifelse{html}{\out{W m<sup>-1</sup> K<sup>-1</sup>}}{\eqn{W m^-1 K^-1}{ASCII}}) are calculated following \insertCite{Denny2006;textual}{TrenchR}.
#'
#' @export
#'
#' @references
#' \insertAllCited{}
#'
#' @examples
#' Tb_limpet(T_a = 25,
#' T_r = 30,
#' l = 0.0176,
#' h = 0.0122,
#' S = 1300,
#' u = 1,
#' psi = 30,
#' c = 1,
#' position = "anterior")
#'
Tb_limpet <- function (T_a,
T_r,
l,
h,
S,
u,
psi,
c,
position = "anterior") {
stopifnot(l > 0,
h > 0,
S > 0,
u >= 0,
psi >= 0,
psi <= 90,
c >= 0,
c <= 1,
length(position) == 1,
position %in% c("anterior", "posterior", "broadside"))
# Conversions
T_a <- celsius_to_kelvin(T_a)
T_r <- celsius_to_kelvin(T_r)
psi <- degrees_to_radians(psi)
r <- l / 2
# Calculations
# Short wave heat transfer
# Area of the limpet's shell (m^2) projected in the direction at which sunlight strikes the organism (Pennell and Deignan 1989)
if (tan(psi) < r / h) {
Ap <- pi * r^2 * cos(psi) + h * r * sin(psi) - pi * r^2 / 2 * cos(psi)
} else {
Ap <- pi * r^2 * cos(psi)
}
# short-wave absorptivity of the shell (the fraction of light energy that is absorbed)
alpha_sw <- 0.68
# Long-wave energy transfer
# View factor simulating limpets as a cone (Campbell and Norman 1998).
Vs <- cos(psi) * r / sqrt(r^2 + h^2)
# lateral area of a limpet shell (m^2)
Al <- pi * r * sqrt(h^2 + r^2)
# long-wave emissivity of the shell (proportion)
eps_ws <- 0.97
sigma <- stefan_boltzmann_constant()
# clear sky emissivity (Campbell and Norman 1998, 10.11)
eps_ac <- 9.2 * 10^-6 * T_a^2
# emissivity of air with clouds (same as above, 10.12)
eps_wa <- (1 - 0.84 * c) * eps_ac + 0.84 * c
# Convective heat transfer
# conductivity of air (W m^-1 K^-1)
Ka <- 0.00501 + 7.2 * 10^-5 * T_a
# kinematic viscosity of air (m^2 s^-1)
v <- -1.25 * 10^-5 + 9.2 * 10^-8 * T_a
# Reynolds number
Re <- u * l / v
if (position == "anterior") {
a <- 1.955
b <- 0.371
} else if (position == "posterior") {
a <- 1.881
b <- 0.376
} else if (position == "broadside"){
a <- 1.304
b <- 0.404
}
# Nusselt number
Nu <- a * Re^b
# Heat transfer coefficient (W m^-2 K^-1)
hc <- a * Ka * ((u / v)^b)*(l^(b - 1))
# area of the shell in convective contact with the air (m^2)
A_cv <- Al
# Conductive heat transfer
# area of conductive contact between the limpet's foot and the rock (m^2)
A_cd <- pi * r^2
# thermal conductivity of rock (W m^-1 K^-1)
Kr <- 3.06
q1 <- Ap * alpha_sw * S
q2 <- Vs * Al * eps_ws * sigma * T_a^4 * (eps_wa - 1)
q3 <- 4 * Vs * Al* eps_ws * sigma * T_a^3
q4 <- hc * A_cv
q5 <- A_cd * Kr
T_b <- (q1 + q2 + (q3 + q4) * T_a + q5 * T_r) / (q3 + q4 + q5)
kelvin_to_celsius(T_b)
}