from pylab import *
from scipy.optimize import fsolve
|
|
# Part a
G = 700 # Solar Irradiation (W/m-K)
eps_g = 0.9 # Emissivity of the glass
sigma = 5.67E-8 # Stefan-Boltzman Constant (W/m2-K4)
Tsurr = 25+273.15 # Surrounding Temperature (K)
def T1(Ts):
eta = 0.28-0.001*(Ts-273.15) # Describes efficiency of semiconductor
lg = 0.003 # length of glass (m)
kg = 1.4 # thermal conductivity of glass (W/m-K)
la = 0.001 # length of Adhesive (m)
ka = 145 # thermal conductivity of adhesive (W/m-K)
T1 = Ts - 0.83*G*(1-eta)*(lg/kg + la/ka)
return T1
def f(Ts):
eta = 0.28-0.001*(Ts-273.15)
hrad = eps_g*sigma*(T1(Ts)+Tsurr)*(T1(Ts)**2+Tsurr**2)
hconv = 16.55
return 0.83*G*(1-eta) + 0.1*G - (T1(Ts) - Tsurr)*(hrad + hconv)
Tsi = fsolve(f,350)
eta = 0.28-0.001*(Tsi-273.15)
Pelec = 0.83*G*0.1*1*eta
T1_a = T1(Tsi)
print 'Tsi = ',(Tsi-273.15), 'degrees C'
print 'eta = ',eta
print 'Pelec = ', Pelec, 'W'
print 'T1 = ', T1_a
Tsi = 48.8077827068 degrees C eta = 0.231192217293 Pelec = 13.4322678247 W T1 = 320.997536484 Tsi = 48.8077827068 degrees C eta = 0.231192217293 Pelec = 13.4322678247 W T1 = 320.997536484 |
# Part b
def f(Ts):
eta = 0.28-0.001*(Ts-273.15)
hrad = eps_g*sigma*(T1(Ts)+Tsurr)*(T1(Ts)**2+Tsurr**2)
hconv = 26.24
return 0.83*G*(1-eta) + 0.1*G - (T1(Ts) - Tsurr)*(hrad + hconv)
Tsi = fsolve(f,350)
eta = 0.28-0.001*(Tsi-273.15)
Pelec = 0.83*G*0.1*1*eta
T1_b = T1(Tsi)
print 'Tsi = ',(Tsi-273.15), 'degrees C'
print 'eta = ',eta
print 'Pelec = ', Pelec, 'W'
print 'T1 = ', T1_b
Tsi = 41.9227309533 degrees C eta = 0.238077269047 Pelec = 13.8322893316 W T1 = 314.121084208 Tsi = 41.9227309533 degrees C eta = 0.238077269047 Pelec = 13.8322893316 W T1 = 314.121084208 |
Tsi_plot = []
Pelec_plot = []
um = linspace(0,10,22)
for i in range(22):
ReL = um[i]*0.1/20.92E-6
hconv = 30.0E-3/0.1*(0.037*ReL**(4./5.)*(0.7)**(1./3.))
def f(Ts):
eta = 0.28-0.001*(Ts-273.15)
hrad = eps_g*sigma*(T1(Ts)+Tsurr)*(T1(Ts)**2+Tsurr**2)
return 0.83*G*(1-eta) + 0.1*G - (T1(Ts) - Tsurr)*(hrad + hconv)
Tsi = fsolve(f,350)-273.15
eta = 0.28-0.001*(Tsi)
Pelec = 0.83*G*0.1*1*eta
Tsi_plot.append(Tsi)
Pelec_plot.append(Pelec)
close('all')
subplot(211)
plot(um,Tsi_plot)
ylabel('Silicon Temperature (C)')
xlabel('velocity (m/s)')
subplot(212)
plot(um,Pelec_plot,'g')
xlabel('velocity (m/s)')
ylabel('Power (W)')
savefig('.')
|
|
|