cdmuir / tealeaves

Solve for leaf temperature using energy balance
Other
11 stars 1 forks source link

discontinuity in Tleaf with changing windspeed #12

Open CraigBarton9 opened 2 years ago

CraigBarton9 commented 2 years ago

I wanted to see how sensitive leaf temperature is to wind speed so kept everything as default and varied windspeed using script below.

require(tealeaves)

leaf_par <- make_leafpar() enviro_par <- make_enviropar(replace=list(wind=set_units(c(0.1,0.2,0.5,0.6,0.65, 0.7,1,2,3),m/s))) constants <- make_constants() pars <- c(leaf_par, enviro_par, constants)

tl<-tleaves(leaf_par, enviro_par, constants, quiet = FALSE, set_units = TRUE) with(tl,plot(T_leaf-T_air~wind))

there is a discontinuity in the leaf temperature as wind goes from 0.60 to 0.65 m/s. Not sure what is causing it.

cdmuir commented 2 years ago

Thanks for your inquiry. I am really busy with teaching right now, but I will try to help soon. My hunch is that it a switch from free to forced convection. The Nusselt number changes abruptly at a cutoff in Reynolds number. Details described in the Appendix of Muir 2019 (https://doi.org/10.1093/aobpla/plz054). I need to double-check this is correct.

CraigBarton9 commented 2 years ago

Yes I was thinking it would be the change from free to forced convection or from laminar to turbulent flow at the leaf surface. So Reynolds number. Somewhat stochastic and empirical but it does make a big step in estimates of leaf temperature. One almost needs to provide a confidence range for leaf temperature to show that at under certain conditions it is uncertain due to chaos creeping in.

From: Chris Muir @.> Sent: Friday, 28 January 2022 3:46 AM To: cdmuir/tealeaves @.> Cc: Craig Barton @.>; Author @.> Subject: Re: [cdmuir/tealeaves] discontinuity in Tleaf with changing windspeed (Issue #12)

Thanks for your inquiry. I am really busy with teaching right now, but I will try to help soon. My hunch is that it a switch from free to forced convection. The Nusselt number changes abruptly at a cutoff in Reynolds number. Details described in the Appendix of Muir 2019 (https://doi.org/10.1093/aobpla/plz054). I need to double-check this is correct.

- Reply to this email directly, view it on GitHubhttps://github.com/cdmuir/tealeaves/issues/12#issuecomment-1023427079, or unsubscribehttps://github.com/notifications/unsubscribe-auth/AEOTW4H6R6P2NCB2GJHMCNLUYFZCXANCNFSM5MSVEYPA. Triage notifications on the go with GitHub Mobile for iOShttps://apps.apple.com/app/apple-store/id1477376905?ct=notification-email&mt=8&pt=524675 or Androidhttps://play.google.com/store/apps/details?id=com.github.android&referrer=utm_campaign%3Dnotification-email%26utm_medium%3Demail%26utm_source%3Dgithub. You are receiving this because you authored the thread.Message ID: @.***>