Sneak Peak Video of the 
New Solar Hydrogen Home DVD
Coming SOON!

Download Over 100Meg of
FREE Hydrogen Video
Ride in the Famous H2 Geo
Click Here

0 then ssgain=ssgain-tgain else tgain=0 360 ts=ts+(is+ssgain)/cs'find new shelf temp (f) 370 if ts>120 then ts=120'limit upper shelf temp (f) 380 dhw=8000/24'water heating load (btu/h) 390 tt=tt+(it+tgain-dhw)/ct'find new tank temp (f) 400 if tt>170 then tt=170'limit upper tank temp (f) 410 if httmax then ttmax=tt'find min and max tank temps 440 x=(h-hs)*640/(he-hs) 450 pset(x,349-v*(ta))'plot ambient temp 460 'pset(x,349-v*(trs))'plot room temp setpoint 470 pset(x,349-v*(tr))'plot room temp 480 'pset(x,349-v*(ts))'plot shelf temp 490 pset(x,349-v*(tt))'plot tank temp 500 if hour=.5 then line (x,349)-(x,345)'mark midnights 510 h=h+1'increment hours 520 if h75f? -----x----------x---------------x--------- duct fan | | | shelfair<75? --------|--------- cloudy day pump it would be nice to change the shelf air setpoint to 55 in setback times, but that complicates the controls... nick
re: deployable doubt dispellers
18 mar 2006
the virtual 20' cube in rochester ny is now more of a virtual house.
the sizing calc below (based on an average december day) has twinwall
polycarbonate sunspace glazing and water shelves under the ceiling to
store heat from warm sunspace air on an average day and an air heater
inside the sunspace with some fin-tube pipe near the top to heat water
for showers and a large cloudy-day heat storage tank...

20 aceil=800:rceil=40'ceiling area (ft^2) and r-value
30 gceil=aceil/rceil'ceiling thermal conductance (btu/h-f)
40 awn=0:awe=24:aww=24:aws=48'n,e,w,s window areas (ft^2)
50 aw=awn+awe+aww+aws'total window area (ft^2)
60 rw=1/.25:rwall=34'window (1/u) and wall r-values
70 gwind=aw/rw'window thermal conductance (btu/h-f)
80 asg=240:gsg=.58:tsg=.8'sspace glazing area (ft^2) and u-value and solar xmsn
90 gsun=asg/(rwall+1/gsg)'nighttime sunspace conductance (btu/h-f)
100 wgain=.5*(160*awn+260*awe+250*aww+540*aws)'window solar gain (btu/day)
110 egain=226/27'indoor electrical use (kwh/day)
120 pgain=12*300'heat gain from 1 12h/day occupant(s) (btu/day)
130 igain=(wgain+3412*egain+pgain)'internal heat gain (btu/day)
140 perim1=120:perim2=80'1st and 2nd floor perimeters (ft)
150 awall=(perim1+perim2)*8-aw-asg'wall area (ft^2)
160 gwall=awall/rwall'wall thermal conductance (btu/h-f)
170 area1=800:area2=400'1st and 2nd floor areas (ft^2)
180 vol=(area1+area2)*8'volume of heated space (ft^3)
190 ach=.2:cfm=vol*ach/60'natural air leakage (cfm)
200 g=gceil+gwind+gsun+gwall+cfm'house thermal conductance (btu/h-f)
210 print"1000'conductances (btu/h-f):"
220 print"1010'ceiling","window","wall","air leaks","total"
230 print"1020'";gceil,gwind,gwall,cfm,g
240 tavg=29.1:tpek=35.8:tday=(tavg+tpek)/2'average dec outdoor temps (f)
250 ssgain=asg*tsg*540'sunspace gain (btu/day)
260 sigain=4*8*(.8*540-6*(70-tday)*gsg)'solar siding gain (btu/day)
270 gnosun=g-gsun'remainder of house (non-sunspace) conductance (btu/h-f)
280 snloss=18*(70-tavg)*gsun'sunspace night loss (btu/night)
290 nsnloss=24*(65-tavg)*gnosun'non-sunspace house loss (btu/day)
300 twi=47.6:twg=100:twh=110'incoming, gwhx and hot water temps (f)
310 dhw=50000!'water heating (btu/day from aug 2004 gas bill)
320 chw=dhw/(twh-twi)'hot water quantity (lb/day)
330 shwreq=chw*(twh-twg)'solar water heating requirement (btu/day)
340 print"1050'solar water heating:",shwreq,"(btu/day)"
350 netgain=ssgain+sigain+igain-shwreq-snloss-nsnloss'net gain (btu/day)
360 ts=tday+netgain/(6*asg*gsg)'daytime sunspace temp (f)
370 imass=2000'inherent structural thermal mass (btu/f)
380 nightheat=(70-60)*imass+.75*(igain)'heat supplied by mass and gains (btu)
390 shelfheat=snloss+.75*nsnloss-nightheat'heat supplied by shelves (btu)
400 n=4:ash=n*4*8'number of 4'x8' shelves and one-sided shelf area (ft^2)
410 t=1'initial shelf water depth (feet)
420 dloss=(60-(tavg-(tpek-tavg)))*g'dawn heat loss (btu/h)
430 tmin=60+dloss/(3*ash)'min water temp (f)
440 rc=62.33*t/3'time constant (hours)
450 tmax=ts+(tmin-ts)*exp(-6/rc)'dusk temp (f)
460 th=shelfheat/(ash*62.33*(tmax-tmin))'new water depth (feet)
470 if abs(th-t)>.01 then t=th: goto 440'iterate
480 print"1030'shelves","depth (in)","tmax (f)","tmin (f)","weight (lb)"
490 print"1040'";n,12*t,tmax,tmin,62.33*32*t
500 aah=8*12'air heater glazing area (ft^2)
510 ahgain=aah*.72*540-shwreq'air heater gain (btu/day)
520 tah=ts+ahgain/(6*aah)'air heater temp (f)
530 fah=64:gah=5*fah'air heater fin-tube (feet)
540 tcavg=tah-shwreq/(6*gah)'cloudy day tank temp on average day (f)
550 fcd=64:gcd=5*fcd:'cloudy-day fin-tube (feet)
560 closs=(70-tday)*g'cloudy daytime heat loss (btu/h)
570 tcmin=60+closs/gcd'min usable cloudy tank temp (f)
580 cgain=3412*egain+pgain'cloudy-day heat internal gains (btu/day)
590 cdheat=24*(62.5-tavg)*g+shwreq-cgain'cloudy-day heat required (btu/day)
600 tankcap=5*cdheat/(tcavg-tcmin)'cloudyr-ay tank heat cap (btu/f)
610 tankvol=tankcap/8.33'cloudy-day tank volume (gallons)
620 print"1080'tsunspace (f)","tairheat (f)","tcavg (f)","tcmin (f)"
630 print"1090'";ts,tah,tcavg,tcmin
640 print"1100'cloudy-tank volume:",tankvol,"(gallons)"

conductances (btu/h-f):
ceiling  window        wall          air leaks     total
20       24            37.17647      32            119.8946

shelves  depth (in)    tmax (f)      tmin (f)      weight (lb)
4        1.28014       95.51366      71.73969      212.7763

solar water heating:    8012.821     (btu/day)

tsunspace (f)          tairheat (f)  tcavg (f)     tcmin (f)
97.21359               148.1025      143.9291      74.06889

cloudy-tank volume:     618.287      (gallons)

the program below prepares a "winter" file from nrel's rochester tmy2
(typical meteorological year) hourly weather data file, wban 14768.
for other locations, change 14768 in line 50 below to 13739 for phila,
24033 for billings, mt, and so on...

10 pi=4*atn(1):screen 9:key off:cls
20 line (0,0)-(639,349),,b:df=.434
30 for tr= 20 to 70 step 10'temp ref lines
40 line (0,349-5*(tr-10))-(639,349-5*(tr-10)):next
50 f$="14768.tm2"'nrel tmy2 file name (rochester)
60 open f$ for input as #1
70 open "winter" for output as #2
80 line input#1,s$'read header
90 city$=mid$(s$,8,25)
100 lat=val(mid$(s$,40,2))+val(mid$(s$,43,2))/60
110 lon=val(mid$(s$,48,3))+val(mid$(s$,52,2))/60
120 print#2,city$,lat,lon
130 gosub 190
140 pass=1
150 close #1
160 open f$ for input as #1
170 line input#1,s$'ignore header
180 gosub 190
190 for h=1 to 8760'hour of year
200 line input#1,s$
210 month=val(mid$(s$,4,2))'month of year (1-12)
220 if pass=0 and month<10 goto 610
230 if pass=1 and month=2 then end
240 wh=h+8760*pass
250 day=val(mid$(s$,6,2))'day of month
260 hour=val(mid$(s$,8,2))-.5'hour of day
270 n=h/24'day of year (1 to 365)
280 tdb=val(mid$(s$,68,4))*.18+32'dry bulb temp (f)
290 pset(df*(wh-8030),349-5*(tdb-10))
300 igloh=val(mid$(s$,18,4))*.317'global horizontal radiation (btu/ft^2)
310 print#2,month;day;hour;tdb;igloh;
320 idif=val(mid$(s$,30,4))*.317'diffuse horizontal radiation (btu/ft^2)
330 idir=val(mid$(s$,24,4))*.317'direct normal radiation (btu/ft^2)
340 l=pi*lat/180'phila latitude (radians)
350 t=hour'solar time (est)
360 x=-sin(pi*23.45/180)*cos(2*pi*(n+10)/365.25)
370 d=atn(x/sqr(-x*x+1))'sin^-1(x) = declination (radians)
380 w=2*pi*(t-12)/24'hour angle (radians)
390 x=cos(l)*cos(d)*cos(w)+sin(l)*sin(d)
400 thetas=-atn(x/sqr(-x*x+1))+pi/2'cos^-1(x) = sun zenith angle (radians)
410 x=cos(d)*sin(w)/sin(thetas)
420 if x^2 >=1 then phis=-1.570796327#:goto 440
430 phis=atn(x/sqr(-x*x+1))'sin^-1(x) = sun azimuth angle (radians)
440 for phipd=0 to 180 step 90'azimuth angle of plane (degrees)
450 phip=pi*phipd/180
460 x=sin(thetas)*cos(phis-phip)
470 thetai=-atn(x/sqr(-x*x+1))+pi/2'incidence angle to surface (radians)
480 if thetai>=pi/2 then thetai=pi/2
490 rhog=.2'ground reflectance
500 iglop=idir*cos(thetai)+idif/2+igloh*rhog/2'radiation on surface (btu/ft^2)
510 print#2,iglop;
520 'if phipd = 0 then pset(h-8030,349-iglop)
530 next phipd
540 phip=pi*270/180
550 x=sin(thetas)*cos(phis-phip)
560 thetai=-atn(x/sqr(-x*x+1))+pi/2'incidence angle to surface (radians)
570 if thetai>=pi/2 then thetai=pi/2
580 iglop=idir*cos(thetai)+idif/2+igloh*rhog/2'radiation on surface (btu/ft^2)
590 print#2,iglop
600 if hour=.5 then line (df*(wh-8030),349)-(df*(wh-8030),345)
610 next h
620 return

here are 362 worst-case hours ("winter2") from the winter file.
each line has month day hour and sun on the ground and 4 walls.

 1  5  8.5  11.66  24.726  78.53785  9.129583  9.129583  85.46374
 1  5  9.5  12.92  70.374  202.1037  12.74336  12.74336  148.7658
 1  5  10.5  14.36  94.46601  161.2371  30.84408  30.84408  84.11995
 1  5  11.5  15.62  114.437  180.4199  37.27918  37.27918  56.29107
 1  5  12.5  17.06  100.172  96.7497  55.6245  49.32519  49.32519
 1  5  13.5  17.42  97.002  139.1144  78.05207  35.85268  35.85268
 1  5  14.5  17.6  73.544  123.799  96.11642  25.42338  25.42338
 1  5  15.5  17.96  44.063  113.6599  123.81  12.64828  12.64828
 1  5  16.5  17.24  11.729  16.84752  23.54156  6.244897  6.244897
 1  5  17.5  16.7  0  0  0  0  0
 1  5  18.5  15.98  0  0  0  0  0
 1  5  19.5  15.98  0  0  0  0  0
 1  5  20.5  15.98  0  0  0  0  0
 1  5  21.5  15.98  0  0  0  0  0
 1  5  22.5  15.98  0  0  0  0  0
 1  5  23.5  15.98  0  0  0  0  0
 1  6  .5  15.98  0  0  0  0  0
 1  6  1.5  17.06  0  0  0  0  0
 1  6  2.5  17.96  0  0  0  0  0
 1  6  3.5  19.04  0  0  0  0  0
 1  6  4.5  19.76  0  0  0  0  0
 1  6  5.5  20.3  0  0  0  0  0
 1  6  6.5  21.02  0  0  0  0  0
 1  6  7.5  22.28  .9510001  1.066952  .5705999  .5705999  1.381782
 1  6  8.5  23.72  18.703  22.752  10.1123  10.1123  24.04432
 1  6  9.5  24.98  30.432  20.25946  17.7837  17.7837  19.5655
 1  6  10.5  26.24  48.501  31.47705  28.4666  28.4666  29.69874
 1  6  11.5  27.68  57.06  36.34515  33.7605  33.7605  34.10436
 1  6  12.5  28.94  64.351  40.977  38.23839  37.8181  37.8181
 1  6  13.5  29.3  60.23  37.57801  36.44666  35.6625  35.6625
 1  6  14.5  29.66  43.112  26.54027  26.26304  25.5502  25.5502
 1  6  15.5  30.02  30.115  18.27957  18.30125  18.069  18.069
 1  6  16.5  29.66  8.876001  5.3256  5.3256  5.3256  5.3256
 1  6  17.5  29.3  0  0  0  0  0
 1  6  18.5  28.94  0  0  0  0  0
 1  6  19.5  28.94  0  0  0  0  0
 1  6  20.5  28.94  0  0  0  0  0
 1  6  21.5  28.94  0  0  0  0  0
 1  6  22.5  29.3  0  0  0  0  0
 1  6  23.5  29.66  0  0  0  0  0
 1  7  .5  30.02  0  0  0  0  0
 1  7  1.5  29.3  0  0  0  0  0
 1  7  2.5  28.76  0  0  0  0  0
 1  7  3.5  28.04  0  0  0  0  0
 1  7  4.5  26.96  0  0  0  0  0
 1  7  5.5  26.06  0  0  0  0  0
 1  7  6.5  24.98  0  0  0  0  0
 1  7  7.5  24.62  1.902  1.1412  1.1412  1.1412  1.1412
 1  7  8.5  24.44  18.703  11.2218  11.2218  11.2218  11.2218
 1  7  9.5  24.08  44.697  27.06548  26.8182  26.8182  26.99653
 1  7  10.5  24.98  61.498  37.1722  36.8988  36.8988  37.01091
 1  7  11.5  26.06  73.227  44.35154  43.7777  43.7777  43.85418
 1  7  12.5  26.96  82.103  49.2618  49.2618  49.2618  49.2618
 1  7  13.5  26.96  68.472  41.35656  41.19532  41.0832  41.0832
 1  7  14.5  26.96  59.279  35.81461  35.74576  35.5674  35.5674
 1  7  15.5  26.96  35.821  21.4926  21.4926  21.4926  21.4926
 1  7  16.5  26.06  8.876001  5.3256  5.3256  5.3256  5.3256
 1  7  17.5  24.98  0  0  0  0  0
 1  7  18.5  24.08  0  0  0  0  0
 1  7  19.5  23.72  0  0  0  0  0
 1  7  20.5  23.36  0  0  0  0  0
 1  7  21.5  23  0  0  0  0  0
 1  7  22.5  22.64  0  0  0  0  0
 1  7  23.5  22.28  0  0  0  0  0
 1  8  .5  21.92  0  0  0  0  0
 1  8  1.5  21.2  0  0  0  0  0
 1  8  2.5  20.66  0  0  0  0  0
 1  8  3.5  19.94  0  0  0  0  0
 1  8  4.5  19.58  0  0  0  0  0
 1  8  5.5  19.4  0  0  0  0  0
 1  8  6.5  19.04  0  0  0  0  0
 1  8  7.5  19.04  1.585  1.280405  .9509999  .9509999  1.492708
 1  8  8.5  19.04  18.703  18.98735  10.5878  10.5878  19.89165
 1  8  9.5  19.04  48.501  62.43346  22.91909  22.91909  51.47656
 1  8  10.5  19.4  75.763  65.67029  40.5443  40.5443  50.86708
 1  8  11.5  19.58  89.394  83.20006  45.0774  45.0774  50.1676
 1  8  12.5  19.94  97.319  115.6677  51.34473  41.43189  41.43189
 1  8  13.5  19.58  112.218  144.5825  86.19238  45.45778  45.45778
 1  8  14.5  19.4  68.155  76.52848  64.28863  32.33399  32.33399
 1  8  15.5  19.04  33.919  29.62897  30.79042  18.9249  18.9249
 1  8  16.5  17.42  5.389  5.049526  6.326246  3.074899  3.074899
 1  8  17.5  15.62  0  0  0  0  0
 1  8  18.5  14  0  0  0  0  0
 1  8  19.5  13.64  0  0  0  0  0
 1  8  20.5  13.28  0  0  0  0  0
 1  8  21.5  12.92  0  0  0  0  0
 1  8  22.5  14.9  0  0  0  0  0
 1  8  23.5  17.06  0  0  0  0  0
 1  9  .5  19.04  0  0  0  0  0
 1  9  1.5  19.4  0  0  0  0  0
 1  9  2.5  19.58  0  0  0  0  0
 1  9  3.5  19.94  0  0  0  0  0
 1  9  4.5  21.02  0  0  0  0  0
 1  9  5.5  21.92  0  0  0  0  0
 1  9  6.5  23  0  0  0  0  0
 1  9  7.5  24.98  0  0  0  0  0
 1  9  8.5  26.96  8.559  6.234627  4.976899  4.976899  6.373759
 1  9  9.5  28.94  22.19  14.97004  12.997  12.997  14.42619
 1  9  10.5  30.02  73.861  59.99554  40.3541  40.3541  48.44024
 1  9  11.5  30.92  76.714  75.42444  37.62789  37.62789  42.6845
 1  9  12.5  32  86.858  96.09641  46.74127  39.11779  39.11779
 1  9  13.5  29.3  67.204  59.59585  44.99608  34.7749  34.7749
 1  9  14.5  26.78  65.302  75.17882  62.82557  30.30519  30.30519
 1  9  15.5  24.08  33.919  33.42665  35.1322  18.1324  18.1324
 1  9  16.5  22.64  3.487  2.0922  2.0922  2.0922  2.0922
 1  9  17.5  21.38  0  0  0  0  0
 1  9  18.5  19.94  0  0  0  0  0
 1  9  19.5  19.22  0  0  0  0  0
 1  9  20.5  18.68  0  0  0  0  0
 1  9  21.5  17.96  0  0  0  0  0
 1  9  22.5  17.96  0  0  0  0  0
 1  9  23.5  17.96  0  0  0  0  0
 1  10  .5  17.96  0  0  0  0  0
 1  10  1.5  18.32  0  0  0  0  0
 1  10  2.5  18.68  0  0  0  0  0
 1  10  3.5  19.04  0  0  0  0  0
 1  10  4.5  19.4  0  0  0  0  0
 1  10  5.5  19.58  0  0  0  0  0
 1  10  6.5  19.94  0  0  0  0  0
 1  10  7.5  19.94  .317  .1902  .1902  .1902  .1902
 1  10  8.5  19.94  12.046  7.6968  7.069101  7.069101  7.768201
 1  10  9.5  19.94  45.331  29.09809  26.8816  26.8816  28.49098
 1  10  10.5  19.94  45.648  28.16168  27.0718  27.0718  27.52146
 1  10  11.5  19.94  58.328  35.28282  34.9968  34.9968  35.03515
 1  10  12.5  19.94  69.42301  41.93981  41.69215  41.6538  41.6538
 1  10  13.5  20.3  64.351  39.54181  38.90182  38.4521  38.4521
 1  10  14.5  20.66  45.014  27.2546  27.18726  27.0084  27.0084
 1  10  15.5  21.02  27.896  16.99733  17.0453  16.5791  16.5791
 1  10  16.5  20.3  9.51  5.706001  5.706001  5.706001  5.706001
 1  10  17.5  19.76  0  0  0  0  0
 1  10  18.5  19.04  0  0  0  0  0
 1  10  19.5  17.96  0  0  0  0  0
 1  10  20.5  17.06  0  0  0  0  0
 1  10  21.5  15.98  0  0  0  0  0
 1  10  22.5  14.9  0  0  0  0  0
 1  10  23.5  14  0  0  0  0  0
 1  11  .5  12.92  0  0  0  0  0
 1  11  1.5  12.92  0  0  0  0  0
 1  11  2.5  12.92  0  0  0  0  0
 1  11  3.5  12.92  0  0  0  0  0
 1  11  4.5  13.64  0  0  0  0  0
 1  11  5.5  14.36  0  0  0  0  0
 1  11  6.5  15.08  0  0  0  0  0
 1  11  7.5  17.06  1.585  .9510001  .9510001  .9510001  .9510001
 1  11  8.5  19.04  15.85  9.51  9.51  9.51  9.51
 1  11  9.5  21.02  38.357  23.2601  23.0142  23.0142  23.1932
 1  11  10.5  22.1  54.524  43.73963  29.8614  29.8614  35.60037
 1  11  11.5  23  54.841  46.24183  29.1006  29.1006  31.40359
 1  11  12.5  24.08  88.76  81.83446  49.687  44.69699  44.69699
 1  11  13.5  24.08  76.39701  67.79425  51.20267  39.4982  39.4982
 1  11  14.5  24.08  53.89  41.51027  38.37128  29.9565  29.9565
 1  11  15.5  24.08  25.36  16.35991  16.53236  14.899  14.899
 1  11  16.5  23.72  4.438  2.826035  2.934538  2.6628  2.6628
 1  11  17.5  23.36  0  0  0  0  0
 1  11  18.5  23  0  0  0  0  0
 1  11  19.5  23  0  0  0  0  0
 1  11  20.5  23  0  0  0  0  0
 1  11  21.5  23  0  0  0  0  0
 1  11  22.5  22.64  0  0  0  0  0
 1  11  23.5  22.28  0  0  0  0  0
 1  12  .5  21.92  0  0  0  0  0
 1  12  1.5  22.64  0  0  0  0  0
 1  12  2.5  23.36  0  0  0  0  0
 1  12  3.5  24.08  0  0  0  0  0
 1  12  4.5  24.08  0  0  0  0  0
 1  12  5.5  24.08  0  0  0  0  0
 1  12  6.5  24.08  0  0  0  0  0
 1  12  7.5  23.72  1.585  .9510001  .9510001  .9510001  .9510001
 1  12  8.5  23.36  14.582  9.424281  8.5907  8.5907  9.52474
 1  12  9.5  23  36.455  22.94207  21.7145  21.7145  22.61043
 1  12  10.5  23.72  48.501  29.75737  28.9421  28.9421  29.28004
 1  12  11.5  24.26  63.4  38.32534  38.04  38.04  38.07842
 1  12  12.5  24.98  66.887  41.2418  40.00733  39.8152  39.8152
 1  12  13.5  24.98  65.302  39.56612  39.24803  39.0227  39.0227
 1  12  14.5  24.98  58.645  35.27393  35.20773  35.0285  35.0285
 1  12  15.5  24.98  32.968  19.98907  20.01438  19.7808  19.7808
 1  12  16.5  24.26  9.193  5.5158  5.5158  5.5158  5.5158
 1  12  17.5  23.72  0  0  0  0  0
 1  12  18.5  23  0  0  0  0  0
 1  12  19.5  23  0  0  0  0  0
 1  12  20.5  23  0  0  0  0  0
 1  12  21.5  23  0  0  0  0  0
 1  12  22.5  22.64  0  0  0  0  0
 1  12  23.5  22.28  0  0  0  0  0
 1  13  .5  21.92  0  0  0  0  0
 1  13  1.5  21.92  0  0  0  0  0
 1  13  2.5  21.92  0  0  0  0  0
 1  13  3.5  21.92  0  0  0  0  0
 1  13  4.5  22.28  0  0  0  0  0
 1  13  5.5  22.64  0  0  0  0  0
 1  13  6.5  23  0  0  0  0  0
 1  13  7.5  23.72  3.487  2.0922  2.0922  2.0922  2.0922
 1  13  8.5  24.26  23.458  14.28275  14.0748  14.0748  14.30856
 1  13  9.5  24.98  56.109  33.9105  33.6654  33.6654  33.84478
 1  13  10.5  26.06  74.495  45.08125  44.5385  44.5385  44.76404
 1  13  11.5  26.96  85.273  51.73374  51.16381  51.16381  51.24073
 1  13  12.5  28.04  87.809  52.81185  52.56537  52.5269  52.5269
 1  13  13.5  28.04  45.331  31.3571  28.18627  25.9306  25.9306
 1  13  14.5  28.04  71.00801  96.28086  78.177  28.65679  28.65679
 1  13  15.5  28.04  36.772  65.10101  71.50357  13.97969  13.97969
 1  13  16.5  27.32  14.899  19.18325  26.77858  7.988396  7.988396
 1  13  17.5  26.78  0  0  0  0  0
 1  13  18.5  26.06  0  0  0  0  0
 1  13  19.5  26.78  0  0  0  0  0
 1  13  20.5  27.32  0  0  0  0  0
 1  13  21.5  28.04  0  0  0  0  0
 1  13  22.5  28.4  0  0  0  0  0
 1  13  23.5  28.58  0  0  0  0  0
 1  14  .5  28.94  0  0  0  0  0
 1  14  1.5  28.94  0  0  0  0  0
 1  14  2.5  28.94  0  0  0  0  0
 1  14  3.5  28.94  0  0  0  0  0
 1  14  4.5  28.22  0  0  0  0  0
 1  14  5.5  27.68  0  0  0  0  0
 1  14  6.5  26.96  0  0  0  0  0
 1  14  7.5  26.96  .634  .3804  .3804  .3804  .3804
 1  14  8.5  26.96  17.118  10.89322  10.2708  10.2708  10.97287
 1  14  9.5  26.96  39.625  24.83986  23.6165  23.6165  24.5144
 1  14  10.5  27.68  47.55  41.30137  25.043  25.043  31.81666
 1  14  11.5  28.22  86.224  68.57505  46.6624  46.6624  49.62752
 1  14  12.5  28.94  100.172  102.8482  54.37561  46.78919  46.78919
 1  14  13.5  28.94  80.518  62.07062  51.64188  44.1898  44.1898
 1  14  14.5  28.94  56.743  50.72783  45.07649  29.4493  29.4493
 1  14  15.5  28.94  31.383  27.66031  29.02542  17.0863  17.0863
 1  14  16.5  27.68  12.997  23.99868  36.08822  6.371695  6.371695
 1  14  17.5  26.24  0  0  0  0  0
 1  14  18.5  24.98  0  0  0  0  0
 1  14  19.5  23.9  0  0  0  0  0
 1  14  20.5  23  0  0  0  0  0
 1  14  21.5  21.92  0  0  0  0  0
 1  14  22.5  19.94  0  0  0  0  0
 1  14  23.5  17.96  0  0  0  0  0
 1  15  .5  15.98  0  0  0  0  0
 1  15  1.5  18.68  0  0  0  0  0
 1  15  2.5  21.38  0  0  0  0  0
 1  15  3.5  24.08  0  0  0  0  0
 1  15  4.5  24.08  0  0  0  0  0
 1  15  5.5  24.08  0  0  0  0  0
 1  15  6.5  24.08  0  0  0  0  0
 1  15  7.5  23.72  .317  .1902  .1902  .1902  .1902
 1  15  8.5  23.36  12.997  8.419151  7.7982  7.7982  8.501076
 1  15  9.5  23  33.285  20.7894  19.8125  19.8125  20.53165
 1  15  10.5  23  72.593  67.90038  37.0573  37.0573  49.94211
 1  15  11.5  23  87.809  87.40825  42.22439  42.22439  48.35425
 1  15  12.5  23  90.662  87.82219  49.75361  43.7777  43.7777
 1  15  13.5  22.28  104.61  140.5403  81.33593  38.83248  38.83248
 1  15  14.5  21.74  85.907  144.0851  113.1585  26.65968  26.65968
 1  15  15.5  21.02  52.622  117.9352  131.648  14.93068  14.93068
 1  15  16.5  18.32  12.046  32.56771  51.68226  5.008592  5.008592
 1  15  17.5  15.62  0  0  0  0  0
 1  15  18.5  12.92  0  0  0  0  0
 1  15  19.5  14  0  0  0  0  0
 1  15  20.5  14.9  0  0  0  0  0
 1  15  21.5  15.98  0  0  0  0  0
 1  15  22.5  14.36  0  0  0  0  0
 1  15  23.5  12.56  0  0  0  0  0
 1  16  .5  10.94  0  0  0  0  0
 1  16  1.5  11.66  0  0  0  0  0
 1  16  2.5  12.2  0  0  0  0  0
 1  16  3.5  12.92  0  0  0  0  0
 1  16  4.5  11.66  0  0  0  0  0
 1  16  5.5  10.22  0  0  0  0  0
 1  16  6.5  8.959999  0  0  0  0  0
 1  16  7.5  10.58  2.219  4.069766  1.014399  1.014399  6.204218
 1  16  8.5  12.38  23.775  50.92946  10.46099  10.46099  56.43663
 1  16  9.5  14  58.962  100.9546  22.22168  22.22168  80.36186
 1  16  10.5  17.42  101.757  162.0322  34.26768  34.26768  87.79176
 1  16  11.5  20.66  93.832  115.0019  38.38869  38.38869  48.81027
 1  16  12.5  24.08  77.982  62.8388  44.45371  41.5587  41.5587
 1  16  13.5  25.34  66.887  67.64957  46.78406  31.73169  31.73169
 1  16  14.5  26.78  66.887  63.84694  55.95936  33.63369  33.63369
 1  16  15.5  28.04  26.945  21.45337  22.33171  15.0575  15.0575
 1  16  16.5  28.04  10.144  7.854993  9.207159  5.927899  5.927899
 1  16  17.5  28.04  0  0  0  0  0
 1  16  18.5  28.04  0  0  0  0  0
 1  16  19.5  28.94  0  0  0  0  0
 1  16  20.5  30.02  0  0  0  0  0
 1  16  21.5  30.92  0  0  0  0  0
 1  16  22.5  31.28  0  0  0  0  0
 1  16  23.5  31.64  0  0  0  0  0
 1  17  .5  32  0  0  0  0  0
 1  17  1.5  33.98  0  0  0  0  0
 1  17  2.5  35.96  0  0  0  0  0
 1  17  3.5  37.94  0  0  0  0  0
 1  17  4.5  37.94  0  0  0  0  0
 1  17  5.5  37.94  0  0  0  0  0
 1  17  6.5  37.94  0  0  0  0  0
 1  17  7.5  37.58  1.902  1.1412  1.1412  1.1412  1.1412
 1  17  8.5  37.4  15.85  10.17526  9.351499  9.351499  10.29092
 1  17  9.5  37.04  30.749  19.1792  18.4494  18.4494  18.99006
 1  17  10.5  37.76  39.942  24.50452  23.9652  23.9652  24.19179
 1  17  11.5  38.3  61.815  37.62194  36.772  36.772  36.88794
 1  17  12.5  39.02  62.449  37.75269  37.50805  37.4694  37.4694
 1  17  13.5  39.02  61.498  37.3906  36.92174  36.5818  36.5818
 1  17  14.5  39.02  47.55  29.10099  28.9123  28.3715  28.3715
 1  17  15.5  39.02  30.432  18.2592  18.2592  18.2592  18.2592
 1  17  16.5  38.3  11.412  6.8472  6.8472  6.8472  6.8472
 1  17  17.5  37.76  0  0  0  0  0
 1  17  18.5  37.04  0  0  0  0  0
 1  17  19.5  37.04  0  0  0  0  0
 1  17  20.5  37.04  0  0  0  0  0
 1  17  21.5  37.04  0  0  0  0  0
 1  17  22.5  36.68  0  0  0  0  0
 1  17  23.5  36.32  0  0  0  0  0
 1  18  .5  35.96  0  0  0  0  0
 1  18  1.5  36.32  0  0  0  0  0
 1  18  2.5  36.68  0  0  0  0  0
 1  18  3.5  37.04  0  0  0  0  0
 1  18  4.5  38.3  0  0  0  0  0
 1  18  5.5  39.74  0  0  0  0  0
 1  18  6.5  41  0  0  0  0  0
 1  18  7.5  41.72  2.219  1.3314  1.3314  1.3314  1.3314
 1  18  8.5  42.26  15.85  9.715386  9.51  9.51  9.745148
 1  18  9.5  42.98  41.527  25.00046  24.7577  24.7577  24.93815
 1  18  10.5  44.06  62.766  37.6596  37.6596  37.6596  37.6596
 1  18  11.5  44.96  69.42301  41.93665  41.6538  41.6538  41.69249
 1  18  12.5  46.04  73.861  44.59944  44.3553  44.3166  44.3166
 1  18  13.5  44.78  66.887  40.40133  40.24566  40.1322  40.1322
 1  18  14.5  43.34  53.89  32.41815  32.356  32.1755  32.1755
 1  18  15.5  42.08  18.703  11.2218  11.2218  11.2218  11.2218
 1  18  16.5  40.64  11.095  6.657  6.657  6.657  6.657
 1  18  17.5  39.38  0  0  0  0  0
 1  18  18.5  37.94  0  0  0  0  0
 1  18  19.5  37.04  0  0  0  0  0
 1  18  20.5  35.96  0  0  0  0  0
 1  18  21.5  35.06  0  0  0  0  0
 1  18  22.5  35.42  0  0  0  0  0
 1  18  23.5  35.6  0  0  0  0  0
 1  19  .5  35.96  0  0  0  0  0
 1  19  1.5  35.6  0  0  0  0  0
 1  19  2.5  35.42  0  0  0  0  0
 1  19  3.5  35.06  0  0  0  0  0
 1  19  4.5  34.7  0  0  0  0  0
 1  19  5.5  34.34  0  0  0  0  0
 1  19  6.5  33.98  0  0  0  0  0
 1  19  7.5  33.98  2.536  1.5216  1.5216  1.5216  1.5216
 1  19  8.5  33.98  19.337  11.80701  11.6022  11.6022  11.83765
 1  19  9.5  33.98  43.429  26.14113  25.8989  25.8989  26.07958
 1  19  10.5  33.26  65.302  39.1812  39.1812  39.1812  39.1812
 1  19  11.5  32.72  73.544  44.1264  44.1264  44.1264  44.1264
 1  19  12.5  32  79.884  48.05426  47.81065  47.7719  47.7719
 1  19  13.5  33.08  85.907  105.1852  64.25841  34.26769  34.26769
 1  19  14.5  33.98  77.982  122.0816  97.52248  25.23318  25.23318
 1  19  15.5  35.06  58.011  133.1444  151.1249  14.04307  14.04307
 1  19  16.5  33.44  15.216  49.7723  82.01438  5.484086  5.484086
 1  19  17.5  31.64  0  0  0  0  0
 1  19  18.5  30.02  0  0  0  0  0
 1  19  19.5  28.4  0  0  0  0  0
 1  19  20.5  26.6  0  0  0  0  0
 1  19  21.5  24.98  0  0  0  0  0
 1  19  22.5  24.26  0  0  0  0  0
 1  19  23.5  23.72  0  0  0  0  0
 1  20  .5  23  0  0  0  0  0
 1  20  1.5  22.28  0  0  0  0  0
 1  20  2.5  21.74  0  0  0  0  0
 1  20  3.5  21.02  0  0  0  0  0
 1  20  4.5  21.02  0  0  0  0  0
 1  20  5.5  21.02  0  0  0  0  0
 1  20  6.5  21.02  0  0  0  0  0
 1  20  7.5  24.08  3.17  7.126154  1.267998  1.267998  11.42546
 1  20  8.5  26.96  20.605  41.0968  9.034491  9.034491  46.04781
 1  20  9.5  30.02  37.406  27.0514  21.4926  21.4926  25.65353

the simulation below has 3 capacitors: the inherent mass of the 1200 ft^2 
house (2k btu/f); the shelf capacitance, 1200 pounds of water in flat poly
film ducts on 2 stacks of 3 4'x8' shelves under the ceiling; and an 86" diam
x 4' tall 957 gallon stss heat storage tank with a copper coil heat exchanger
for 100% solar space heat and close to 100% solar hot water in december,
with no evacuated tubes nor antifreeze. and no pvs nor milk-based paints
nor bicycle racks, and a low leed score :-)

the frugal electrical use is 227 kwh/month. there's a good greywater
heat exchanger. the shelf and tank caps grew from the sizing calc during
the simulation, which doesn't assume any beneficial tank stratification
nor explicitly use the shelf and fin-tube conductances, but it looks good
and allows margins for those conductances.

there's 32 ft^2 of solar siding and 240 ft^2 of twinwall polycarb passive
sunspace glazing (5 4'x12' thermaglas panels) with 144 ft^2 of one-layer
polycarb higher temp air heater glazing inside the sunspace in a barra-like
system with a kind of thermal cogeneration and 64' of fin-tube inside
the air heater glazing with a small pump and a differential thermostat.

on an average day, sunspace air naturally thermosyphons up through a one-way
plastic film damper into the shelf space under the ceiling and down through
a vertical duct on the inside of the insulated south wall of the house and
back into the lower part of the sunspace. when the house needs heat and
the shelf space is warm, a fan near the floor pulls air out of the duct
and into the room. another small pump moves tank water up through some
fin-tube inside the house on a cloudy day.

the lowest temp in the weather file is -3 f. rochester ny is a hard climate, 
colder than phila, with half the sun, about 500 btu/ft^2 on a south wall on 
an average december day. the min tank temp is 82 f on 1/15 at 9:30 am, when 
the outdoor temp is 23 f. the room temp setpoint is 70 f for 8 hours per day 
and 50 f otherwise, altho it rarely cools to 50 by dawn. the room temp is 
allowed to overshoot to 74 f during the day, which helps. reducing the room 
mass would allow it to cool more quickly at night and reduce the required 
tank and shelf sizes. so would zoning, eg closing bedroom doors during the 
day or automatically turning heat off when rooms are unoccupied. lines 20
and 30 below adjust the graphic display. line 145 sampled the cap temps on
1/5 during a longer simulation, which were set to those temps at the start
of this shorter simulation using the winter2 file.

10 cls:screen 9:line (0,0)-(639,349),,b
20 tmax=100:v=350/tmax'vertical scaling
30 hs=100:he=330'start and end display hours
40 for tr=10 to 50 step 10'plot temp ref lines
50 line (0,349-v*(tr))-(639,349-v*(tr)):next
60 gcloud=120'cloudy day house conductance (btu/h-f)
70 cr=2000!'room capacitance (btu/f)
80 tr=50.58'initial room temp (f)
90 cs=1200'6-shelf capacitance (btu/f)
100 ts=86.85'initial shelf temp (f)
110 ct=957*8.33'tank capacitance (btu/f)
120 tt=115.84:ttmin=1000'initial tank temp (f)
130 open "winter2" for input as #1:line input#1,h$
140 input#1,month,day,hour,ta,sh,ss,sw,sn,se
145 if month=1 and day=5 and hour=9.5 then print tr,ts,tt
150 gsun=120-240/(34+1/.58)'sunny day house conductance (btu/h-f)
160 ssgain=.8*240*ss-(ts-ta)*.58*240'sunspace solar gain (btu)
170 if ssgain<0 then ssgain=0:g=gcloud else g=gsun'sunspace dampers closed
180 sigain=.9*32*ss-(tr-ta)*32'solar siding gain (btu)
190 if sigain<0 then sigain=0'no heat from solar siding
200 if hour<9 or hour>17 then trs=50 else trs=70'room temp setpoint (f)
210 wgain=.5*(24*se+48*ss+24*sw)'solar gain from windows (btu)
220 ir=wgain+sigain+1440-(tr-ta)*g'net heat into room (btu)
230 tr=tr+ir/cr'find new room temp (f)
240 if tr>74 then tr=74'limit upper room temp (f)
250 if tr>trs then is=0:it=0:goto 330'room needs no heat :-)
260 rheat=(trs-tr)*cr'required room heat (btu)
270 if ts


I got ALL of these 85 Solar Panels for FREE and so can you.  Its in our Ebook
Ready for DOWNLOAD NOW.

Site Meter