平板集热器的集热量和热效率计算的Python程序

基本概念

集热器的示意图

集热器有效利用能Q_u

U_L也成为热损失系数

T_pm不易计算，Q_u计算公式可转化为

T_f,i是流体的进口温度

#集热器有效利用能Qu
Q_u = A_c * F_R * (S-U_L*(T_i-T_a))

集热效率η

#集热效率η
Eta = Q_u /(I_T*A_c)

吸热面吸收的太阳辐射量S

集热器性能的预测需要集热板吸收的太阳能信息。倾斜集热器上入射的太阳能可以用第二章的方法求出。这种入射辐射有三种不同的空间分布：直射辐射、散射辐射和地面反射辐射，每种辐射都必须单独处理。计算的细节取决于使用哪个漫射天空模型（第2.14节至第2.16节）。每小时使用各向同性漫射概念，可修改方程式2.15.1，通过将每个项乘以适当的透射率-吸收率积，得出吸收辐射S：

其中（1+cosβ）/2和（1−cosβ）/2分别是从收集器到天空和从收集器到地面的视角因子。下标b、d和g表示直射、散射和地面反射。对于给定的集热器倾斜，图5.4.1给出了散射和地面反射辐射的有效入射角，图4.11.1和5.4可用于确定适当的吸收率和透射率值。方程5.5.1或5.5.2可用于求出（τα）_d和（τα）_g。计算R_b所需的波束辐射角θ可用于求出（τα）_b。或者，可从覆盖层和吸收体的特性中求出（τα）_n，图5.6.1可用于计算每个辐射流的适当入射角测定三种透射率-吸收率乘积。
I_b、I_d、I和R_b的计算见我的博客计算某地的从日出到日落各时刻的水平面直射和散射辐射，倾斜面总辐射的Python程序

#吸热面吸收的太阳辐射量S
S = (Ib*Rb*(Tau_Alpha_s_b) + Id*(Tau_Alpha_s_d)*(1+cos(Beta*pi/180))/2 + Rho_g*I*(Tau_Alpha_s_g)*(1-cos(Beta*pi/180))/2).evalf()

地面反射辐射的入射角θ_g和散射辐射的入射角θ_d

Beta = 60    #集热器倾斜角β
#散射辐射的入射角θd
Thet_d = 59.7 - 0.1388*Beta + 0.001497*Beta**2
#地面反射辐射的入射角θg
Thet_g = 90 - 0.5788*Beta + 0.002693*Beta**2

（τα）的计算

其中θ₁和θ₂是入射角和折射角。 方程5.1.1表示非极化辐射r⊥的垂直分量，和方程5.1.2表示非极化辐射r||的平行分量。方程5.1.3给出了非极化辐射的反射作为两个分量的平均值。

from sympy import *
Thet_1 = 60    #入射角θ1
Thet_2 = 34.58    #折射角θ2
r_per = ((sin((Thet_2-Thet_1)*pi/180))**2 / (sin((Thet_2+Thet_1)*pi/180))**2).evalf()    #r⊥，perpendicular
r_par = ((tan((Thet_2-Thet_1)*pi/180))**2 / (tan((Thet_2+Thet_1)*pi/180))**2).evalf()    #r||，parallel
r = (r_per+r_par)/2

#计算折射角θ2
#空气的折射率n1=1，玻璃的折射率n2=1.526
#假设入射角θ1 = 60°
from sympy import *
n1 = 1
n2 = 1.526
Thet_1 = 60    #入射角θ1 = 60°
Thet_2 = (asin(n1/n2*sin(Thet_1*pi/180))*180/pi).evalf()
print('折射角θ2：',Thet_2)

结果：

折射角θ2： 34.5770069287450

r⊥和r||不相等（正常入射时除外），初始非极化辐射的透射率τ_r为两个部分的平均透射率。其中下标r是提醒只考虑反射损失。

from sympy import *
Thet_1 = 60    #入射角θ1
Thet_2 = 34.58    #折射角θ2
r_per = ((sin((Thet_2-Thet_1)*pi/180))**2 / (sin((Thet_2+Thet_1)*pi/180))**2).evalf()    #r⊥，perpendicular
r_par = ((tan((Thet_2-Thet_1)*pi/180))**2 / (tan((Thet_2+Thet_1)*pi/180))**2).evalf()    #r||，parallel
#非极化辐射的透射率τr
Tau_r = 1/2 * ((1-r_par)/(1+r_par) + (1-r_per)/(1+r_per))

θ₂——折射角
K——消光系数，The extinction coefficient
L——辐射在介质中传播的距离，m。
其中下标a提醒只考虑了吸收损失

from sympy import *
Thet_1 = 60    #入射角θ1
Thet_2 = 34.58    #折射角θ2
#非极化辐射的吸收率τa
K = 32    #消光系数，The extinction coefficient
L = 2.3*10**(-3)    #辐射在介质中传播的距离，m
Tau_a = (exp(-K*L/cos(Thet_2*pi/180))).evalf()

对于极化的垂直分量，覆盖层的透射τ⊥、反射率ρτ⊥和吸收率ατ⊥

from sympy import *
Thet_1 = 60    #入射角θ1
Thet_2 = 34.58    #折射角θ2
r_per = ((sin((Thet_2-Thet_1)*pi/180))**2 / (sin((Thet_2+Thet_1)*pi/180))**2).evalf()    #r⊥，perpendicular
r_par = ((tan((Thet_2-Thet_1)*pi/180))**2 / (tan((Thet_2+Thet_1)*pi/180))**2).evalf()    #r||，parallel
#非极化辐射的吸收率τa
K = 32    #消光系数，The extinction coefficient
L = 2.3*10**(-3)    #辐射在介质中传播的距离，m
Tau_a = (exp(-K*L/cos(Thet_2*pi/180))).evalf()
#透射率τ
Tau = Tau_a * (1-r_per)**2 / (1-(r_per*Tau_a)**2)
#反射率ρ
Rho = r_per * (1+Tau_a*Tau)
#吸收率α
Alpha = (1-Tau_a) * ((1-r_per)/(1-r_per*Tau_a))

（5.3.1）可简化为（5.3.4）

（5.3.3）可简化为（5.3.5）

（5.3.2）可简化为（5.3.6）

与方程5.3.1至5.3.3相比，方程5.3.4至5.3.6的优点在于，偏振在近似方程中通过单个项τr而不是每个单独光学性质的更复杂表达式来解释。

from sympy import *
Thet_1 = 60    #入射角θ1
Thet_2 = 34.58    #折射角θ2
r_per = ((sin((Thet_2-Thet_1)*pi/180))**2 / (sin((Thet_2+Thet_1