[Ibpsausa] Groundwater temperature estimating
Joe Huang
yjhuang at whiteboxtechnologies.com
Wed Jan 6 12:11:44 PST 2016
Here's an older but more compact version of the Kusuda algorithm written in Fortran 77,
just 20 lines of code but you have to give it the following information:
DIFF = soil diffusivity
DBTMIN = minimum average monthly air temperature
DBTMAX = maximum average monthly air temperature
DBTYAV -= average annual air temperature
The output is TG, the average monthly ground temperature.
I stripped this from the DOE-2.1E program, and gave it to the EnergyPlus Team back in 2002
(?) for use in their weather processing program. Although the subroutine is structured to
produce average monthly ground temperatures, I suggest that you code it to do daily or at
least weekly calculations (but using a "running" month, i.e., calculate DBTMIN and DBTMAX
from -15 to +15 days around the day) to avoid the monthly staircase of temperatures using
the algorithm at the monthly timestep.
Joe
Joe Huang
White Box Technologies, Inc.
346 Rheem Blvd., Suite 205A
Moraga CA 94556
yjhuang at whiteboxtechnologies.com
http://weather.whiteboxtechnologies.com for simulation-ready weather data
(o) (925)388-0265
(c) (510)928-2683
"building energy simulations at your fingertips"
On 1/6/2016 10:57 AM, Neal Kruis wrote:
> Andrew,
>
> EnergyPlus has source for the Kusuda algorithm here
> <https://github.com/NREL/EnergyPlus/blob/develop/src/EnergyPlus/GroundTemperatureModeling/KusudaAchenbachGroundTemperatureModel.cc#L188>.
> As Joe mentioned, this is a pretty crude estimate. For a more accurate estimate you
> might want to use a finite difference calculation that takes into account the surface
> boundary conditions (solar, wind, IR radiation to the sky, etc.).
>
> Neal Kruis
> 720 583-5587 <tel:720%20583-5587>
>
> Vice President
> Big Ladder Software
> 1624 Market Street, Suite 304
> Denver, CO 80202
> bigladdersoftware.com <http://bigladdersoftware.com/>
>
>
> On Mon, Jan 4, 2016 at 8:50 AM, Andrew Corney <andrew.corney at sefaira.com
> <mailto:andrew.corney at sefaira.com>> wrote:
>
> Hi Nick and everyone who replied - thank you all for your help. It's great to see so
> much support in the community :-)
>
> So it turns out what I was actually after was a way to generate ground temperatures
> for a random set of weather data. Someone kindly send me the fortran script that E+
> uses for its EPW files which is exactly what I needed.
>
> Have a great new year all look forward to catching up in 2016.
>
> Andrew
>
> Andrew Corney • PE, M.CIBSE, M.ASHRAE
>
> Product Director
>
> Sefaira - Sustainability • Performance • Design
>
> M +44 (0)7887 059 722 <tel:%2B44%20%280%297887%20059%20722>
>
> T +44 (0)2037 147 619 <tel:%2B44%20%280%292037%20147%20619>
>
> Skype andrew.corney.sefaira
>
> Readthis article <http://continuingeducation.bnpmedia.com/article.php?L=480&C=1446>
> in Architecture Record and complete a quiz to earl AIA and GBCI credit
>
>
> This email was sent from Sefaira UK Ltd.
>
> For company contact and registration information, please
> visit:http://www.sefaira.com/info/contact
>
>
>
> On Mon, Jan 4, 2016 at 2:47 PM, Nicholas Caton <ncaton at catonenergy.com
> <mailto:ncaton at catonenergy.com>> wrote:
>
> Hi Andrew!
>
> Oklahoma State University in Stillwater, OK produces quite a bit of research
> concerning the developing science of ground source wellfield system design and
> verification. Also, individuals I personally consider leaders in our circles of
> energy simulation when it comes to anything ground source spent time with that
> department.
>
> If there exists an established best practice for extrapolating ground
> temperatures (short of direct measurement), there's a good chance someone at
> that institution knows and possibly even derived the answer ;-).
>
> Of course if you're lucky someone will reply with a direct answer - here's hoping!
>
> ~Nick
>
> On Dec 23, 2015 2:38 PM, "Andrew Corney" <andrew.corney at sefaira.com
> <mailto:andrew.corney at sefaira.com>> wrote:
>
> Hi All,
>
> Does anyone know of a peer-reviewed paper / methodology they could share
> with me that describes a process for converting annual weather data into
> monthly groundwater temperatures?
>
> Obviously a lot of EPW files already have ground water temperatures but many
> don't and a lot of other weather data doesn't have that information so it'd
> be great to know if there was a well-recognised process for making those
> estimates.
>
> Thanks in advance for any help anyone can provide.
>
> best regards,
>
> Andrew
>
> Andrew Corney • PE, M.CIBSE, M.ASHRAE
>
> Product Director
>
> Sefaira - Sustainability • Performance • Design
>
> M +44 (0)7887 059 722 <tel:%2B44%20%280%297887%20059%20722>
>
> T +44 (0)2037 147 619 <tel:%2B44%20%280%292037%20147%20619>
>
> Skype andrew.corney.sefaira
>
> Readthis article
> <http://continuingeducation.bnpmedia.com/article.php?L=480&C=1446> in
> Architecture Record and complete a quiz to earl AIA and GBCI credit
>
>
> This email was sent from Sefaira UK Ltd.
>
> For company contact and registration information, please
> visit:http://www.sefaira.com/info/contact
>
>
>
> _______________________________________________
> Ibpsausa mailing list
> Ibpsausa at lists.onebuilding.org <mailto:Ibpsausa at lists.onebuilding.org>
> http://lists.onebuilding.org/listinfo.cgi/ibpsausa-onebuilding.org
> To unsubscribe from this mailing list send a blank message to
> IBPSAUSA-UNSUBSCRIBE at ONEBUILDING.ORG
> <mailto:IBPSAUSA-UNSUBSCRIBE at ONEBUILDING.ORG>
>
>
>
> _______________________________________________
> Ibpsausa mailing list
> Ibpsausa at lists.onebuilding.org <mailto:Ibpsausa at lists.onebuilding.org>
> http://lists.onebuilding.org/listinfo.cgi/ibpsausa-onebuilding.org
> To unsubscribe from this mailing list send a blank message to
> IBPSAUSA-UNSUBSCRIBE at ONEBUILDING.ORG <mailto:IBPSAUSA-UNSUBSCRIBE at ONEBUILDING.ORG>
>
>
>
>
> _______________________________________________
> Ibpsausa mailing list
> Ibpsausa at lists.onebuilding.org
> http://lists.onebuilding.org/listinfo.cgi/ibpsausa-onebuilding.org
> To unsubscribe from this mailing list send a blank message to IBPSAUSA-UNSUBSCRIBE at ONEBUILDING.ORG
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.onebuilding.org/pipermail/ibpsausa-onebuilding.org/attachments/20160106/f07994cc/attachment.htm>
-------------- next part --------------
Lines 19790-2000 of wth.f
SUBROUTINE GTEMP(DIF,TMIN,TMAX,TAV,TG) GTEMP 2
DIMENSION AMON(12),TG(12) GTEMP 3
DATA AMON / 15.,46.,74.,95.,135.,166.,196.,227.,258.,288., GTEMP 4
1 319.,349. / GTEMP 5
DATA P,PI,PO / 8760.,3.14159265,0.6 / GTEMP 6
C GTEMP 7
BETA = SQRT(PI/(DIF*P))*10. GTEMP 8
X = EXP(-BETA) GTEMP 9
X2 = X*X GTEMP 10
C = COS(BETA) GTEMP 11
S = SIN(BETA) GTEMP 12
Y = X2 - 2.*X*C + 1. GTEMP 13
Y = Y / (2.*BETA*BETA) GTEMP 14
GM = SQRT(Y) GTEMP 15
Z = (1.-X*(C+S)) / (1.-X*(C-S)) GTEMP 16
PHI = ATAN(Z) GTEMP 17
BO = (TMAX-TMIN)*0.5 GTEMP 18
DO 40 I=1,12 GTEMP 19
THETA = AMON(I)*24. GTEMP 20
40 TG(I) = TAV - BO*COS(2.*(PI/P)*THETA-PO-PHI)*GM + 460. GTEMP 21
RETURN GTEMP 22
END GTEMP 23
Line 2990
DIFF = .025 -089 31
Line 3430
CALL GTEMP(DIFF,DBTMIN,DBTMAX,DBTYAV,TG) PACKER 458
More information about the Ibpsausa
mailing list