Similarly,
Fs
X
(1−ǫ)2nFn,1−t
with
Fn,1−t=q(wt+w1)2+(2n+1)2h2
t−qw2
1+h2
t−2nht2w1.
Also
Fs
1−1=ρsF1(t1)−1+(ρs)3F1(t1t1)−1+. . . =∞
X
n=1
(1 −ǫ)2n−1Fs
n,1−1
Fs
n,1−1=qw2
1+(nht)2−nhtw1.
And, finally,
Fs
t−1=A1
At
Fs
1−t=1
2Fs
1−t.
Because of the infinite series for these view factors, and because we would like to see the variation of q2with
variable shield position (subject to hb+ht=h=const), this is best solved with little computer program given
below. Of course, the minimum heat loss is obtained when hb→0, since only then all radiation is blocked
from directly going to outer space.
PROGRAM POBLEM6.19
IMPLICIT NONE
INTEGER :: i,II,n
REAL :: sigma=5.670e-8,w1,w2,hh,T2,eps,rho,hb,ht,F2b,F22,Ftt,F1t,Ft1,F11,add,denom,q2
! Dimensions
w1=1.
w2=2.
T2=600.
eps=.1
rho=1.-eps
II=20
! Loop over different hb
Ftt=0.
DO n=1,100
add=rho**(2*n-1)*(sqrt(w2**2+(n*ht)**2)-sqrt(w1**+ht**2)-(n-1)*ht)/w2
Ftt=Ftt+add
IF(add<1.e-6) EXIT