9.9 Appendix I

Calculation of the Right-Hand Limit of the Climate Space
In their original paper, Porter and Gates (1969) state with regard to Equation (9.3), “… An estimate was made of ¯Qa as generally related to air temperature for value of absorptivity from 0.2 to 1.0.” Gates (1977) shows how to calculate S and s as a function of latitude, time of year, and time of day. A representative value of 40° was chosen for latitude and then values of S, s were calculated under clear sky conditions in the late morning and afternoon. For r, a representative value of 0.15 was probably used. Curves of Ta and Tg as a function of time of day such as shown in Figure 9.2 were then taken from weather bureau statistics. Ra and Rg could then be calculated using the Stefan-Boltzmann law. All the numbers necessary to estimate ¯Qa are then available. The final step is to choose pairs of Ta and Qa that are to be used. In Figure 9.3, 9:00-10:00 are hours of the day when this is true. This procedure is repeated at several latitudes and times of the year, from which the right-hand boundary can be derived.

Campbell (1977 pp.89-92) presents simplified equations to calculate the left- and right-hand boundaries of the climate space. He includes a correction factor to average the longwave radiation from the ground and the sky for the left-hand boundary. The direct beam and diffuse shortwave radiation fluxes are simply given for the right-hand boundary. The reflected shortwave component seems to be included in the diffuse term which at 25% of the direct flux is higher than Gates (1978) gives. The absorbed shortwave radiation is then added to the left-hand boundary values. Therefore as Ta increases the shortwave component is constant using Campbell’s equations but using Gates’ method the shortwave flux increases.

9.9.1 Appendix II

Derivation of the Constant ks for the Convection Coefficient
Mitchell (1976) reported that the best overall relationship between the Reynolds and Nusselt numbers is given by

Nu=0.34Re0.6 Recalling that Reynolds number Re is the ratio of interim forces to viscous forces in the fluid, we write

Re=VLν

where
V = the fluid velocity (m s-1)
L = the characteristic length (m), and
ν = the kinematic viscosity (m2 s-1)

The Nusselt number is a way to scale the rate of heat transfer as a function of wind velocity, size of the organism and fluid thermal diffusivity.
It can also be expressed as

Nu=hcLk

where
hc = heat transfer coefficient (W m-2 °C-1),
L = characteristic length (m), and
k = thermal conductivity W m-1 °C-1).

Mitchell (1976) defined the characteristic length as

L=(Mbρ)1/3

where Mb is the mass kg
ρ is the mass density kg m-3.

Using these expressions we can solve for the heat transfer coefficient as a function of weight and wind velociy. Rearranging Equation (9.21) and substituting the Equation (9.19) for Nu, we have hc=kL(0.34Re0.6) We can use Equations (9.20) and (9.22) to incorporate wind velocity and weight, respectively.

hc=0.34kL(VLν)0.6=0.34kν0.6V0.6L0.4=0.34×kν0.6V0.6(Mbρ)0.133=0.34×k×ρ0.133ν0.6V0.6M0.133b

Letting ks=0.34×k×ρ0.133ν0.6 if ρ=1×103 kg m-3 which is the density of water k=2.57×102Wm1K1 ν=1.51×105m2s1 at 20 °C then ks=17.24.

To see the error of assumng k and ν at 20 °C we can compare the ratio of kν0.6.

Air temperature (°C) kν0.6
-10 20.79
20 20.07
50 19.51

The difference over the 20 °C value is 1.2820.79=6.4 percent.