FortMin=70;tPrimeMin=tMin;ac=Input["Please input the area of the collector in meters squared."];ϵL=.7;mDotCpMin=3400*31*12*3600;cL=ϵLmDotCpMin;d=15000*31*12*3600;ms=100;cp=4190;h=-4.002;k=4.702;tenv=21;uas=6;deltat=31*12*3600;b=3.85;c=-0.15;dd=-1.959;frPrimeUc=3.5;hBar=19.16*1000000;frPrimeταBar=0.72*0.94;n=31;kT=0.480;capitalA=7.10-20.00kT+12.08;capitalB=-8.02+18.16kT-10.68;capitalC=-1.02+4.10kT-1.96;hdbyh=1.317-3.023kT+3.372-1.769;frPrimeταn=0.72;ρg=0.2;s=39.8;ll=39.8;tambPrime=16.8;z=dcL;δ=23.45Sin(284+135)Degree;a=0.015;g=0.2136;f1=;x=acfrPrimeUc;qu=qMax-a(Exp[bf1]-1)(1-Exp[cx])(Exp[ddz])*d;hs=ArcCos[-Tan[llDegree]Tan[δDegree]];smallrdn=;smallrn=smallrdn(1.07+.025Sin[hsDegree-60Degree]);rB=(Cos[(ll-s)Degree]*Cos[δDegree]*Cos[hsDegree]+Sin[(ll-s)Degree]*Sin[δDegree])(Cos[llDegree]Cos[δDegree]*Cos[hsDegree]+Sin[llDegree]Sin[δDegree]);r=(1-hdbyh)rB+hdbyh+ρg;rBn=(Cos[(ll-s)Degree]*Cos[δDegree]+Sin[(ll-s)Degree]*Sin[δDegree])(Cos[llDegree]Cos[δDegree]+Sin[llDegree]Sin[δDegree]);rn=1-*hdbyhrBn+*hdbyh+ρg;xcMin=;ϕMax=ExpcapitalA+capitalB(xcMin+capitalC);htBar=rhBar;qMax=acfrPrimeταBarhtBarnϕMax;f2=(qu-uas((tPrimeMin+g(Exp[kf1]-1)Exp[hd/cL])-tenv)deltat),Abs[f1-f2]>0.001,tPrimeMin=tPrimeMin+.001;f1=;x=acfrPrimeUc;qu=qMax-a(Exp[bf1]-1)(1-Exp[cx])(Exp[ddz])*d;hs=ArcCos[-Tan[llDegree]Tan[δDegree]];smallrdn=;smallrn=smallrdn(1.07+.025Sin[hsDegree-60Degree]);rB=(Cos[(ll-s)Degree]*Cos[δDegree]*Cos[hsDegree]+Sin[(ll-s)Degree]*Sin[δDegree])(Cos[llDegree]Cos[δDegree]*Cos[hsDegree]+Sin[llDegree]Sin[δDegree]);r=(1-hdbyh)rB+hdbyh+ρg;rBn=(Cos[(ll-s)Degree]*Cos[δDegree]+Sin[(ll-s)Degree]*Sin[δDegree])(Cos[llDegree]Cos[δDegree]+Sin[llDegree]Sin[δDegree]);rn=1-*hdbyhrBn+*hdbyh+ρg;xcMin=;ϕMax=ExpcapitalA+capitalB(xcMin+capitalC);htBar=rhBar;qMax=acfrPrimeταBarhtBarnϕMax;f2=(qu-uas((tPrimeMin+g(Exp[kf1]-1)Exp[hd/cL])-tenv)deltat),=f2;
2
kT
2
kT
2
kT
2
kT
3
kT
360
365
-0.76
ms*
cp
1000
350
-0.704
ms*
cp
1000
350
(ϵL(mDotCpMin))(tPrimeMin-tMin)
d
deltat
d
Pi
24
1-Cos[hsDegree]
Sin[hsDegree]-hsDegreeCos[hsDegree]
Pi
180
1+Cos[sDegree]
2
1-Cos[sDegree]
2
smallrdn
smallrn
smallrdn
smallrn
1+Cos[sDegree]
2
1-Cos[sDegree]
2
1
smallrnrnhBar
frPrimeUc(tPrimeMin-tambPrime)
frPrimeταBar
rn
r
2
xcMin
1
d
(ϵL(mDotCpMin))(tPrimeMin-tMin)
d
deltat
d
Pi
24
1-Cos[hsDegree]
Sin[hsDegree]-hsDegreeCos[hsDegree]
Pi
180
1+Cos[sDegree]
2
1-Cos[sDegree]
2
smallrdn
smallrn
smallrdn
smallrn
1+Cos[sDegree]
2
1-Cos[sDegree]
2
1
smallrnrnhBar
frPrimeUc(tPrimeMin-tambPrime)
frPrimeταBar
rn
r
2
xcMin
1
d
c
ac
Forac=10,ac≤100,FortMin=70;tPrimeMin=tMin;ϵL=.7;mDotCpMin=3400*31*12*3600;cL=ϵLmDotCpMin;d=15000*31*12*3600;ms=100;cp=4190;h=-4.002;k=4.702;tenv=21;uas=6;deltat=31*12*3600;b=3.85;c=-0.15;dd=-1.959;frPrimeUc=3.5;hBar=19.16*1000000;frPrimeταBar=0.72*0.94;n=31;kT=0.480;capitalA=7.10-20.00kT+12.08;capitalB=-8.02+18.16kT-10.68;capitalC=-1.02+4.10kT-1.96;hdbyh=1.317-3.023kT+3.372-1.769;frPrimeταn=0.72;ρg=0.2;s=39.8;ll=39.8;tambPrime=16.8;z=dcL;δ=23.45Sin(284+135)Degree;a=0.015;g=0.2136;f1=;x=acfrPrimeUc;qu=qMax-a(Exp[bf1]-1)(1-Exp[cx])(Exp[ddz])*d;hs=ArcCos[-Tan[llDegree]Tan[δDegree]];smallrdn=;smallrn=smallrdn(1.07+.025Sin[hsDegree-60Degree]);rB=(Cos[(ll-s)Degree]*Cos[δDegree]*Cos[hsDegree]+Sin[(ll-s)Degree]*Sin[δDegree])(Cos[llDegree]Cos[δDegree]*Cos[hsDegree]+Sin[llDegree]Sin[δDegree]);r=(1-hdbyh)rB+hdbyh+ρg;rBn=(Cos[(ll-s)Degree]*Cos[δDegree]+Sin[(ll-s)Degree]*Sin[δDegree])(Cos[llDegree]Cos[δDegree]+Sin[llDegree]Sin[δDegree]);rn=1-*hdbyhrBn+*hdbyh+ρg;xcMin=;ϕMax=ExpcapitalA+capitalB(xcMin+capitalC);htBar=rhBar;qMax=acfrPrimeταBarhtBarnϕMax;f2=(qu-uas((tPrimeMin+g(Exp[kf1]-1)Exp[hd/cL])-tenv)deltat),Abs[f1-f2]>0.001,tPrimeMin=tPrimeMin+.001;f1=;x=acfrPrimeUc;qu=qMax-a(Exp[bf1]-1)(1-Exp[cx])(Exp[ddz])*d;hs=ArcCos[-Tan[llDegree]Tan[δDegree]];smallrdn=;smallrn=smallrdn(1.07+.025Sin[hsDegree-60Degree]);rB=(Cos[(ll-s)Degree]*Cos[δDegree]*Cos[hsDegree]+Sin[(ll-s)Degree]*Sin[δDegree])(Cos[llDegree]Cos[δDegree]*Cos[hsDegree]+Sin[llDegree]Sin[δDegree]);r=(1-hdbyh)rB+hdbyh+ρg;rBn=(Cos[(ll-s)Degree]*Cos[δDegree]+Sin[(ll-s)Degree]*Sin[δDegree])(Cos[llDegree]Cos[δDegree]+Sin[llDegree]Sin[δDegree]);rn=1-*hdbyhrBn+*hdbyh+ρg;xcMin=;ϕMax=ExpcapitalA+capitalB(xcMin+capitalC);htBar=rhBar;qMax=acfrPrimeταBarhtBarnϕMax;f2=(qu-uas((tPrimeMin+g(Exp[kf1]-1)Exp[hd/cL])-tenv)deltat),=f2;;ac=ac+10,{};ListLinePlot[{{10,},{20,},{30,},{40,},{50,},{60,},{70,},{80,},{90,},{100,}}]
2
kT
2
kT
2
kT
2
kT
3
kT
360
365
-0.76
ms*
cp
1000
350
-0.704
ms*
cp
1000
350
(ϵL(mDotCpMin))(tPrimeMin-tMin)
d
deltat
d
Pi
24
1-Cos[hsDegree]
Sin[hsDegree]-hsDegreeCos[hsDegree]
Pi
180
1+Cos[sDegree]
2
1-Cos[sDegree]
2
smallrdn
smallrn
smallrdn
smallrn
1+Cos[sDegree]
2
1-Cos[sDegree]
2
1
smallrnrnhBar
frPrimeUc(tPrimeMin-tambPrime)
frPrimeταBar
rn
r
2
xcMin
1
d
(ϵL(mDotCpMin))(tPrimeMin-tMin)
d
deltat
d
Pi
24
1-Cos[hsDegree]
Sin[hsDegree]-hsDegreeCos[hsDegree]
Pi
180
1+Cos[sDegree]
2
1-Cos[sDegree]
2
smallrdn
smallrn
smallrdn
smallrn
1+Cos[sDegree]
2
1-Cos[sDegree]
2
1
smallrnrnhBar
frPrimeUc(tPrimeMin-tambPrime)
frPrimeταBar
rn
r
2
xcMin
1
d
c
ac
c
10
c
20
c
30
c
40
c
50
c
60
c
70
c
80
c
90
c
100