CN108664720B - 一种非均匀辐照度下的聚光光伏光热系统性能计算方法 - Google Patents
一种非均匀辐照度下的聚光光伏光热系统性能计算方法 Download PDFInfo
- Publication number
- CN108664720B CN108664720B CN201810409878.9A CN201810409878A CN108664720B CN 108664720 B CN108664720 B CN 108664720B CN 201810409878 A CN201810409878 A CN 201810409878A CN 108664720 B CN108664720 B CN 108664720B
- Authority
- CN
- China
- Prior art keywords
- photovoltaic
- thermal
- photo
- photovoltaic cell
- micro
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000004364 calculation method Methods 0.000 title abstract description 9
- 230000008878 coupling Effects 0.000 claims abstract description 25
- 238000010168 coupling process Methods 0.000 claims abstract description 25
- 238000005859 coupling reaction Methods 0.000 claims abstract description 25
- 238000009826 distribution Methods 0.000 claims abstract description 24
- 230000000694 effects Effects 0.000 claims abstract description 9
- 239000006185 dispersion Substances 0.000 claims abstract description 5
- 239000012530 fluid Substances 0.000 claims description 62
- 239000011521 glass Substances 0.000 claims description 37
- 238000012360 testing method Methods 0.000 claims description 21
- 238000000034 method Methods 0.000 claims description 18
- 238000001816 cooling Methods 0.000 claims description 11
- 238000006243 chemical reaction Methods 0.000 claims description 8
- 230000004907 flux Effects 0.000 claims description 7
- 239000000463 material Substances 0.000 claims description 7
- 238000004134 energy conservation Methods 0.000 claims description 5
- 230000005855 radiation Effects 0.000 claims description 5
- 238000012546 transfer Methods 0.000 claims description 4
- 238000003491 array Methods 0.000 claims description 3
- 238000010276 construction Methods 0.000 claims description 3
- 230000008569 process Effects 0.000 claims description 3
- 238000002834 transmittance Methods 0.000 claims description 3
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 2
- 238000005206 flow analysis Methods 0.000 claims description 2
- 230000005684 electric field Effects 0.000 abstract description 3
- 210000004027 cell Anatomy 0.000 description 119
- 238000011161 development Methods 0.000 description 4
- 230000008646 thermal stress Effects 0.000 description 3
- 230000005494 condensation Effects 0.000 description 2
- 238000009833 condensation Methods 0.000 description 2
- 230000003247 decreasing effect Effects 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 229910021420 polycrystalline silicon Inorganic materials 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 150000001875 compounds Chemical class 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000011438 discrete method Methods 0.000 description 1
- 230000001678 irradiating effect Effects 0.000 description 1
- 238000004321 preservation Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Photovoltaic Devices (AREA)
Abstract
一种非均匀辐照度下的聚光光伏光热系统性能计算方法，通过对非均匀辐照下的光伏光热组件进行二维离散，建立光伏光热微元的光‑热‑电耦合模型。运用二维离散光‑热‑电耦合模型可以得到离散的光伏光热微元的热电数据，如电场分布和温度分布等。最后根据光伏光热微元的热电数据来确定聚光光伏光热系统总体的热电性能参数。与传统的集中参数模型相比，本发明提出的二维离散光‑热‑电耦合模型具有计算准确性更高、应用范围更广以及可以获得系统局部的细节信息等优势。另外，二维离散光‑热‑电耦合模型也可以用来评估辐照度非均匀性对聚光光伏光热系统总体性能参数的影响。
Description
技术领域
本发明属于太阳能聚光光伏光热研究领域，涉及一种非均匀辐照度下的聚光光伏光热系统性能计算方法。
背景技术
在太阳能聚光光伏光热研究领域，由于聚光器自身的聚光原理与特性，或者是不完全对准的跟踪策略等客观因素，非均匀的辐照度总是存在于光伏光热组件表面的。这种非均匀的辐照度直接会导致光伏电池的电流分布和温度分布的非均匀性。局部过大的电流会使局部的欧姆损失增大，从而造成光伏电池的开路电压以及填充因子下降，最终导致光伏电池的光电转换效率下降。局部过高的电池温度还会导致光伏电池表面的热斑和热应力集中，从而造成光伏电池的失效。在现有的用于评估聚光光伏光热系统热电输出性能的模型中，绝大部分采用的是集中参数模型，即采用光伏光热组件表面的平均辐照度以及平均温度进行计算。然而在非均匀辐照度和电池温度工况下，集中参数模型显然无法准确获得系统的热电输出性能，而且也无法获得系统局部的细节信息。因此，对于非均匀辐照度下的聚光光伏光热系统，需要用一种离散的光-热-电耦合模型来预测电池表面的温度分布以及系统总体的热电输出性能。
发明内容
针对非均匀辐照度下的聚光光伏光热系统，本发明的目的在于提供一种非均匀辐照度下的聚光光伏光热系统性能计算方法，利用离散的光-热-电耦合模型准确预测在非均匀辐照度下光伏电池的温度分布以及系统的热电输出性能，同时能够评估聚光器辐照非均匀度对光伏电池的温度分布以及系统热电性能的影响。
为实现上述目的，本发明采用的技术方案如下：
一种非均匀辐照度下的聚光光伏光热系统性能计算方法，包括以下步骤：
1)计算区域离散：在光伏光热组件的宽度方向，将离散单元个数设置为n；在光伏光热组件的长度方向，假设串联的单体光伏电池片个数是m，因此离散单元个数设置为m；光伏光热组件被离散的单元个数为m×n，每个离散的单元称之为光伏光热微元；每一个光伏光热微元都包括一个光伏电池微元和一个流体微元；
2)光伏光热微元上接收到的辐照度的确定：光伏光热组件上接收到的辐照度分布用蒙特卡洛光线追踪方法或者用光学模拟软件计算得到；光伏光热微元上接收到的辐照度由微元上接收到的辐射通量除以微元的受光面积得到；
3)光伏光热微元的能量流动分析，获得关于光伏光热微元的能量守恒方程：以任意一个光伏光热微元(i,j)作为控制容积，对其能量流动进行分析，获得关于光伏光热微元的能量守恒方程：
Qpv,abs(i,j)＝Pm(i,j)+Qf(i,j)+Qloss,conv(i,j)+Qloss,rad(i,j) (1)
其中，Qpv,abs(i,j)是光伏电池微元吸收的辐射通量，Pm(i,j)是光伏电池微元输出的最大电功率，Qf(i,j)是流体微元吸收的热能，Qloss,conv(i,j)是光伏光热微元通过玻璃盖板与环境的对流热损失，Qloss,rad(i,j)是光伏光热微元通过玻璃盖板与环境的辐射热损失；
4)根据光伏光热微元玻璃盖板表面的热平衡方程得到光伏电池微元与玻璃盖板表面的导热量，即光伏电池微元与玻璃盖板表面的导热量等于玻璃盖板与环境之间的辐射换热量和对流换热量之和：
其中，Tpv(i,j)是光伏电池微元的温度，Rpv-c是光伏电池和玻璃盖板之间的热阻，由光伏光热组件的结构和材料的导热系数以及厚度计算；
5)根据流体微元的热平衡方程，得到流体微元获得的热能：
其中，Rpv-f是光伏电池和流体之间的热阻，由光伏电池与冷却通道壁面之间的导热热阻以及冷却通道壁面与流体之间的对流热阻构成；
步骤3)-5)构成了光伏光热微元的热模型；
下面给出光伏光热微元的电模型：
6)计算光伏电池微元的最大功率点电流，由下面两式得到：
其中，Im(i,j)是光伏电池微元的最大功率点电流，Gref是光伏电池标准测试条件下的辐照度，Tref是光伏电池标准测试条件下的温度，αI是光伏电池的电流温度系数；
光伏电池微元在标准测试条件下的最大功率点电流Im,ref(i,j)根据单体光伏电池片在标准测试条件下的最大功率点电流Im,ref来计算：
Im,ref(i,j)＝Im,ref/n (9)
7)计算光伏电池微元的最大功率点电压：
Vm(i,j)＝Vm,ref(i,j){1+βV[Tpv(i,j)-Tref]}-Rs(i,j)[Im(i,j)-Im,ref(i,j)](10)
其中，Vm(i,j)是光伏电池微元的最大功率点电压，βV是光伏电池的电压温度系数；光伏电池微元在标准测试条件下的最大功率点电压Vm,ref(i,j)等于单体光伏电池片在标准测试条件下的最大功率点电压Vm,ref：
Vm,ref(i,j)＝Vm,ref (11)
8)光伏电池微元输出的最大电功率，等于光伏电池微元最大功率点电流和电压之积：
Pm(i,j)＝Im(i,j)·Vm(i,j) (14)
步骤6)-8)完成了光伏光热微元电模型的构建，步骤2)-8)构成了光伏光热微元的光-热-电耦合模型；
根据光伏光热微元的光-热-电耦合模型能够得到光伏电池表面的温度分布，同时能够获得聚光光伏光热系统总体的热电性能参数。
本发明进一步的改进在于，
Qpv,abs(i,j)＝τcαpvGpv(i,j)ΔA (2)
其中，τc是光伏光热组件玻璃盖板的透过率，αpv是光伏电池的吸收率，Gpv(i,j)是光伏光热微元上接收到的辐照度，ΔA是光伏光热微元的受光面积。
本发明进一步的改进在于，
其中，
是流体微元的质量流量，Cp是流体的比热容，Tf(i,j)是流体微元的出口温度，Tf(i,j-1)是流体微元的进口温度。
本发明进一步的改进在于，
Qloss,conv(i,j)＝hwΔA[Tc(i,j)-Tamb] (4)
其中，hw是玻璃盖板和环境的对流换热系数，由风速决定；Tc(i,j)是光伏光热微元玻璃盖板的温度，Tamb是环境温度。
本发明进一步的改进在于，
其中，σ是斯忒藩-玻尔兹曼常数，εc是玻璃盖板的发射率，Tsky是有效天空温度。
本发明进一步的改进在于，公式(10)中，光伏电池微元的串联电阻Rs(i,j)根据离散前后光伏电池阵列的I2R相等的原则进行计算：
nm·[I2(i,j)·Rs(i,j)]＝m[nI(i,j)]2Rs (12)
其中，I(i,j)是光伏电池微元的电流，Rs是单体光伏电池片的串联电阻，公式(12)化简为：
Rs(i,j)＝nRs (13)
本发明进一步的改进在于，根据光伏光热微元的光-热-电耦合模型能够得到光伏电池表面的温度分布，同时能够获得聚光光伏光热系统总体的热电性能参数的具体过程如下：
1)光伏电池阵列最大功率点电流的确定：根据串联效应，光伏电池阵列最大功率点电流由串联的单体光伏电池片最大功率点电流的最小值确定；每个单体光伏电池片最大功率点电流由宽度方向所有光伏电池微元的最大功率点电流叠加得到：
2)光伏电池阵列最大功率点电压的确定：根据并联效应，光伏电池阵列最大功率点电压由长度方向所有单体光伏电池片最大功率点电压叠加得到；每个单体光伏电池片最大功率点电压采用宽度方向所有光伏电池微元的最大功率点电压的平均值：
3)光伏电池阵列输出的最大电功率，等于光伏电池阵列最大功率点电流和电压之积：
Pm＝ImVm (17)
4)聚光光伏光热系统出口流体温度的确定：聚光光伏光热系统的出口流体温度由光伏光热组件冷却流道出口处宽度方向所有流体微元的出口温度的平均值确定：
5)根据步骤3)和4)，获得了光伏电池阵列输出的最大电功率以及聚光光伏光热系统出口流体温度之后，进一步得到聚光光伏光热系统的光电转换效率和光热转换效率：
其中，Gt是入射太阳总辐照度，Ain是聚光器入射光孔面积，
是冷却通道中流体的流量，Tf,in是流体进口温度。
与现有技术相比，本发明具有如下有益效果：通过对非均匀辐照下的光伏光热组件进行二维离散，建立光伏光热微元的光-热-电耦合模型。运用二维离散光-热-电耦合模型可以得到离散的光伏光热微元的热电数据，如电场分布和温度分布等。最后根据光伏光热微元的热电数据来确定聚光光伏光热系统总体的热电性能参数。与传统的集中参数模型相比，本发明提出的二维离散光-热-电耦合模型具有计算准确性更高、应用范围更广以及可以获得系统局部的细节信息等优势。本发明提出的光-热-电多场耦合模型可以获得光伏光热组件的温度场，这对于判断光伏电池的局部热斑以及计算光伏电池的热应力提供了条件。另外，二维离散光-热-电耦合模型也可以用来评估辐照度非均匀性对聚光光伏光热系统总体性能参数的影响。
附图说明
图1为光伏光热组件二维离散示意图。
图2为二维离散光-热-电耦合模型求解整体思路流程图。
图3为离散光伏光热微元热电耦合模型求解思路流程图。
图4为HEMR与LEMR聚光时光伏光热组件宽度方向的辐照度分布曲线。
图5为HEMR与LEMR聚光时光伏电池宽度方向温度分布曲线。
具体实施方式
下面结合附图对本发明作进一步的阐述。
一种非均匀辐照度下的聚光光伏光热系统性能计算方法，包括以下步骤：
1)计算区域离散：在光伏光热组件的宽度方向，即光伏电池的细栅线方向，将离散单元个数设置为n；在光伏光热组件的长度方向，即光伏电池的主栅线方向，假设串联的单体光伏电池片个数是m，因此离散单元个数设置为m；因此光伏光热组件被离散的单元个数为m×n，每个离散的单元称之为光伏光热微元。需要注意的是，每一个光伏光热微元都包括一个光伏电池微元和一个流体微元。为了体现光伏光热微元的二维离散性质，用坐标(i,j)进行标记；
2)光伏光热微元上接收到的辐照度的确定：光伏光热组件上接收到的辐照度分布用蒙特卡洛光线追踪方法或者用光学模拟软件计算得到；光伏光热微元上接收到的辐照度由微元上接收到的辐射通量除以微元的受光面积得到；
3)光伏光热微元的能量流动分析；以任意一个光伏光热微元(i,j)作为控制容积，对其能量流动进行分析，获得关于光伏光热微元的能量守恒方程：
Qpv,abs(i,j)＝Pm(i,j)+Qf(i,j)+Qloss,conv(i,j)+Qloss,rad(i,j) (1)
其中，Qpv,abs(i,j)是光伏电池微元吸收的辐射通量，Pm(i,j)是光伏电池微元输出的最大电功率，Qf(i,j)是流体微元吸收的热能，Qloss,conv(i,j)是光伏光热微元通过玻璃盖板与环境的对流热损失，Qloss,rad(i,j)是光伏光热微元通过玻璃盖板与环境的辐射热损失。
公式(1)中各部分能量分别由下面的公式给出：
Qpv,abs(i,j)＝τcαpvGpv(i,j)ΔA (2)
其中，τc是光伏光热组件玻璃盖板的透过率，αpv是光伏电池的吸收率，Gpv(i,j)是光伏光热微元上接收到的辐照度，ΔA是光伏光热微元的受光面积；
其中，
是流体微元的质量流量，Cp是流体的比热容，Tf(i,j)是流体微元的出口温度，Tf(i,j-1)是流体微元的进口温度。
Qloss,conv(i,j)＝hwΔA[Tc(i,j)-Tamb] (4)
其中，hw是玻璃盖板和环境的对流换热系数，由风速决定；Tc(i,j)是光伏光热微元玻璃盖板的温度，Tamb是环境温度。
其中，σ是斯忒藩-玻尔兹曼常数，εc是玻璃盖板的发射率，Tsky是有效天空温度。
4)光伏光热微元玻璃盖板表面的热平衡方程，即光伏电池微元与玻璃盖板表面的导热量等于玻璃盖板与环境之间的辐射换热量和对流换热量之和：
其中，Tpv(i,j)是光伏电池微元的温度，Rpv-c是光伏电池和玻璃盖板之间的热阻，由光伏光热组件的结构和材料的导热系数以及厚度计算。
5)流体微元的热平衡方程，此处认为光伏光热组件的保温效果非常好，忽略流体与环境之间的热损失，因此光伏电池微元与流体微元的换热量即为流体微元获得的热能：
其中，Rpv-f是光伏电池和流体之间的热阻，由光伏电池与冷却通道壁面之间的导热热阻以及冷却通道壁面与流体之间的对流热阻构成。
至此，步骤3)-5)就构成了光伏光热微元的热模型。从公式(1)可以看出，光伏电池微元输出的最大电功率Pm(i,j)会影响微元的能量流动，进而影响光伏电池微元的温度；而反过来光伏电池微元的温度又会影响其输出的电功率，因此光伏光热微元的热电性能是耦合在一起的。
下面给出光伏光热微元的电模型。
6)光伏电池微元的最大功率点电流。综合考虑辐照度和温度的影响，光伏电池微元的最大功率点电流可由下面两式得到：
其中，Im(i,j)是光伏电池微元的最大功率点电流，Gref是光伏电池标准测试条件下的辐照度(1000W/m2)，Tref是光伏电池标准测试条件下的温度(298K)，αI是光伏电池的电流温度系数。
光伏电池微元在标准测试条件下的最大功率点电流，即Im,ref(i,j)根据单体光伏电池片在标准测试条件下的最大功率点电流Im,ref来计算：
Im,ref(i,j)＝Im,ref/n (9)
7)光伏电池微元的最大功率点电压。对于光伏电池微元的最大功率点电压，一般不仅需要考虑温度的影响，而且还需额外考虑串联电阻的影响：
Vm(i,j)＝Vm,ref(i,j){1+βV[Tpv(i,j)-Tref]}-Rs(i,j)[Im(i,j)-Im,ref(i,j)](10)
其中，Vm(i,j)是光伏电池微元的最大功率点电压，βV是光伏电池的电压温度系数；光伏电池微元在标准测试条件下的最大功率点电压，即Vm,ref(i,j)等于单体光伏电池片在标准测试条件下的最大功率点电压Vm,ref：
Vm,ref(i,j)＝Vm,ref (11)
公式(10)中，光伏电池微元的串联电阻Rs(i,j)根据离散前后光伏电池阵列的I2R相等的原则进行计算：
nm·[I2(i,j)·Rs(i,j)]＝m[nI(i,j)]2Rs (12)
其中，I(i,j)是光伏电池微元的电流，Rs是单体光伏电池片的串联电阻，公式(12)可化简为：
Rs(i,j)＝nRs (13)
8)光伏电池微元输出的最大电功率，等于光伏电池微元最大功率点电流和电压之积：
Pm(i,j)＝Im(i,j)·Vm(i,j) (14)
至此，步骤6)-8)就完成了光伏光热微元电模型的构建，而步骤2)-8)则是构成了光伏光热微元完整的光-热-电耦合模型。基于上述光伏光热微元的光-热-电耦合模型可以得到光伏电池表面的温度分布，同时也可以获得聚光光伏光热系统总体的热电性能参数。
9)光伏电池阵列最大功率点电流的确定。根据串联效应，光伏电池阵列最大功率点电流由串联的单体光伏电池片最大功率点电流的最小值确定；而每个单体光伏电池片最大功率点电流由宽度方向所有光伏电池微元的最大功率点电流叠加得到：
10)光伏电池阵列最大功率点电压的确定。根据并联效应，光伏电池阵列最大功率点电压由长度方向所有单体光伏电池片最大功率点电压叠加得到；而每个单体光伏电池片最大功率点电压采用宽度方向所有光伏电池微元的最大功率点电压的平均值：
11)光伏电池阵列输出的最大电功率，等于光伏电池阵列最大功率点电流和电压之积：
Pm＝ImVm (17)
12)聚光光伏光热系统出口流体温度的确定。聚光光伏光热系统的出口流体温度由光伏光热组件冷却流道出口处宽度方向所有流体微元的出口温度的平均值确定：
13)根据步骤11)和12)，获得了光伏电池阵列输出的最大电功率以及聚光光伏光热系统出口流体温度之后，进一步可得到聚光光伏光热系统的光电转换效率和光热转换效率：
其中，Gt是入射太阳总辐照度，Ain是聚光器入射光孔面积，
是冷却通道中流体的流量，Tf,in是流体进口温度。
下面通过一个具体实施例进行详细说明。
1)参见图1，首先确定光伏光热组件宽度和长度方向要离散的单元数量n和m，然后对光伏光热组件进行单元划分，每一个单元称为光伏光热微元，并给定光伏光热微元的编号(i,j)，其中i表示宽度方向的离散单元序号，j表示长度方向的离散单元序号；
2)运用蒙特卡洛光线追踪方法或者光学模拟软件计算太阳光经聚光器后照射在光伏光热组件上的辐照度分布，同时计算任意光伏光热微元上的辐照度Gpv(i,j)，即由光伏光热微元上接收到的辐照通量除以微元的受光面积得到；
3)首先计算光伏光热微元(1,1)的热电数据，具体过程叙述如下；
4)参见图3，首先输入该微元上接收到的辐照度Gpv(i,j)以及风速、环境温度、进口流体温度和流量等初始条件，然后假设出口流体温度初值等于进口流体温度；
5)根据进出口流体温度和流量求出流体微元获得的热能Qf(i,j)；根据公式(7)计算出光伏电池微元的温度Tpv(i,j)，再由公式(4)-(10)求出光伏电池微元输出的最大电功率Pm(i,j)；
6)假设光伏光热微元玻璃盖板温度等于光伏电池微元的温度，通过牛顿迭代求解关于玻璃盖板温度的四次方程，即公式(6)；
7)判断流进和流出光伏光热微元的能量是否守恒，即公式(1)等号两边是否在给定误差范围内成立；
8)若公式(1)在给定误差范围内不成立则更新出口流体温度，并重复步骤5)-7)直到公式(1)在给定误差范围内成立为止；参见图3，这里更新出口流体温度时采用的是0.01K的步长往上叠加；
9)参见图2，若公式(1)在给定误差范围内成立则说明进出光伏光热微元的能量是守恒的，则保存该光伏光热微元的热电数据，并在宽度方向计算下一个光伏光热微元的热电数据，即重复步骤4)-8)，直到计算完宽度方向所有光伏光热微元的热电数据；
10)将步骤9)中保存的宽度方向所有流体微元的出口温度赋值给进口温度，然后在长度方向推进重复步骤4)-9)，直到计算完长度方向所有光伏光热微元的热电数据，最后根据公式(15)-(20)求得聚光光伏光热系统总体的热电性能参数。
该实施例的初始气象数据如下，总辐照度为956.8W/m2，散射辐射比为8.5％，环境温度为32.1℃，风速为0.1m/s。本算例中，聚光光伏光热系统采用的聚光器有两种，即XieHuling等在学术论文Design and performance study of truncated CPC byeliminating multiple reflections of solar radiation in hybrid CPV/T system:Highest and lowest truncation position中设计的两种截短式复合抛物面聚光器HEMR和LEMR。系统所采用的光伏光热组件材料参数如表1所示，光伏电池阵列中采用的单体多晶硅太阳能电池片在标准测试条件下的电性能参数如表2所示。HEMR和LEMR聚光器的聚光比均设置为6倍，在光伏光热组件的宽度方向和长度方向分别离散成26和60个单元。运用本发明提出的二维离散光-热-电耦合模型，根据上述步骤所述进行实施，可以得到光伏电池在宽度方向的温度分布以及系统总体的性能参数。
参见图4，在HEMR和LEMR聚光下，光伏光热组件宽度方向的辐照度分布是不均匀的，而且LEMR聚光时的不均匀性更明显。与之相对应的，图5中光伏电池在两种聚光器聚光时宽度方向的温度分布也是不均匀的，且LEMR聚光时的不均匀性更明显。图5中也可以看出光伏电池宽度方向的温度分布存在两个对称的峰，这两处可能会形成热斑。
参见表3，与假设的均匀辐照条件相比，在HEMR和LEMR聚光条件下，最大功率点电压、填充因子、光电转换效率及其温度系数都有不同程度的降低，而且LEMR聚光时下降的程度更大。这表明在非均匀辐照度条件下，以平均辐照度和光伏电池温度来计算系统的性能参数时，总是会存在偏差的，辐照度的非均匀性越大，偏差越大。
以上实施算例结果表明，在非均匀辐照度条件下，与采用平均辐照度和平均光伏电池温度进行计算的集中参数模型相比，二维离散的光-热-电多场耦合模型计算准确性更高，而且还可以获得光伏电池的温度分布等细节信息，因此本发明提出的计算方法是具有先进性的。
表1光伏光热组件的材料参数
表2标准测试条件下单体多晶硅太阳能电池片的电性能参数
表3聚光光伏光热系统总体性能参数
本发明利用二维离散的方法，将光场、温度场以及电场耦合起来进行求解，可以准确地预测不均匀辐照度下聚光光伏光热系统总体性能，从而有效评估辐照度非均匀性对系统性能的影响。集中参数模型只能获得光伏光热组件的平均温度，而本发明提出的光-热-电多场耦合模型可以获得光伏光热组件的温度场，这对于判断光伏电池的局部热斑以及计算光伏电池的热应力提供了条件。另外，相比于集中参数模型，本发明提出的光-热-电多场耦合模型也具有更广的应用范围。
Claims (7)
1.一种非均匀辐照度下的聚光光伏光热系统性能计算方法，其特征在于，包括以下步骤：
1)计算区域离散：在光伏光热组件的宽度方向，将离散单元个数设置为n；在光伏光热组件的长度方向，假设串联的单体光伏电池片个数是m，因此离散单元个数设置为m；光伏光热组件被离散的单元个数为m×n，每个离散的单元称之为光伏光热微元；每一个光伏光热微元都包括一个光伏电池微元和一个流体微元；
2)光伏光热微元上接收到的辐照度的确定：光伏光热组件上接收到的辐照度分布用蒙特卡洛光线追踪方法或者用光学模拟软件计算得到；光伏光热微元上接收到的辐照度由微元上接收到的辐射通量除以微元的受光面积得到；
3)光伏光热微元的能量流动分析，获得关于光伏光热微元的能量守恒方程：以任意一个光伏光热微元(i,j)作为控制容积，对其能量流动进行分析，获得关于光伏光热微元的能量守恒方程：
Qpv,abs(i,j)＝Pm(i,j)+Qf(i,j)+Qloss,conv(i,j)+Qloss,rad(i,j) (1)
其中，Qpv,abs(i,j)是光伏电池微元吸收的辐射通量，Pm(i,j)是光伏电池微元输出的最大电功率，Qf(i,j)是流体微元吸收的热能，Qloss,conv(i,j)是光伏光热微元通过玻璃盖板与环境的对流热损失，Qloss,rad(i,j)是光伏光热微元通过玻璃盖板与环境的辐射热损失；
4)根据光伏光热微元玻璃盖板表面的热平衡方程得到光伏电池微元与玻璃盖板表面的导热量，即光伏电池微元与玻璃盖板表面的导热量等于玻璃盖板与环境之间的辐射换热量和对流换热量之和：
其中，Tpv(i,j)是光伏电池微元的温度，Rpv-c是光伏电池和玻璃盖板之间的热阻，由光伏光热组件的结构和材料的导热系数以及厚度计算；Tc(i,j)是光伏光热微元玻璃盖板的温度，εc是玻璃盖板的发射率，σ是斯忒藩-玻尔兹曼常数，Tsky是天空温度，hw是玻璃盖板与环境之间的对流换热系数，Tamb是环境温度；
5)根据流体微元的热平衡方程，得到流体微元获得的热能：
其中，Rpv-f是光伏电池和流体之间的热阻，由光伏电池与冷却通道壁面之间的导热热阻以及冷却通道壁面与流体之间的对流热阻构成；Tf(i,j)是流体微元的出口温度，Tf(i,j-1)是流体微元的进口温度，
是流体微元的质量流量，Cp是流体的比热容，ΔA是光伏光热微元的受光面积；
步骤3)-5)构成了光伏光热微元的热模型；
下面给出光伏光热微元的电模型：
6)计算光伏电池微元的最大功率点电流，由下面两式得到：
其中，Im(i,j)是光伏电池微元的最大功率点电流，Gref是光伏电池标准测试条件下的辐照度，Tref是光伏电池标准测试条件下的温度，αI是光伏电池的电流温度系数；Tref是光伏电池标准测试条件下的温度，Rs(i,j)是光伏电池微元的串联电阻；
光伏电池微元在标准测试条件下的最大功率点电流Im,ref(i,j)根据单体光伏电池片在标准测试条件下的最大功率点电流Im,ref来计算：
Im,ref(i,j)＝Im,ref/n (9)
7)计算光伏电池微元的最大功率点电压：
Vm(i,j)＝Vm,ref(i,j){1+βV[Tpv(i,j)-Tref]}-Rs(i,j)[Im(i,j)-Im,ref(i,j)] (10)
其中，Vm(i,j)是光伏电池微元的最大功率点电压，βV是光伏电池的电压温度系数；光伏电池微元在标准测试条件下的最大功率点电压Vm,ref(i,j)等于单体光伏电池片在标准测试条件下的最大功率点电压Vm,ref：
Vm,ref(i,j)＝Vm,ref (11)
8)光伏电池微元输出的最大电功率，等于光伏电池微元最大功率点电流和电压之积：
Pm(i,j)＝Im(i,j)·Vm(i,j) (14)
步骤6)-8)完成了光伏光热微元电模型的构建，步骤2)-8)构成了光伏光热微元的光-热-电耦合模型；
根据光伏光热微元的光-热-电耦合模型能够得到光伏电池表面的温度分布，同时能够获得聚光光伏光热系统总体的热电性能参数。
2.根据权利要求1所述的一种非均匀辐照度下的聚光光伏光热系统性能计算方法，其特征在于，
Qpv,abs(i,j)＝τcαpvGpv(i,j)ΔA (2)
其中，τc是光伏光热组件玻璃盖板的透过率，αpv是光伏电池的吸收率，Gpv(i,j)是光伏光热微元上接收到的辐照度，ΔA是光伏光热微元的受光面积。
3.根据权利要求1所述的一种非均匀辐照度下的聚光光伏光热系统性能计算方法，其特征在于，
其中，
是流体微元的质量流量，Cp是流体的比热容，Tf(i,j)是流体微元的出口温度，Tf(i,j-1)是流体微元的进口温度。
4.根据权利要求1所述的一种非均匀辐照度下的聚光光伏光热系统性能计算方法，其特征在于，
Qloss,conv(i,j)＝hwΔA[Tc(i,j)-Tamb] (4)
其中，hw是玻璃盖板和环境的对流换热系数，由风速决定；Tc(i,j)是光伏光热微元玻璃盖板的温度，Tamb是环境温度，ΔA是光伏光热微元的受光面积。
5.根据权利要求1所述的一种非均匀辐照度下的聚光光伏光热系统性能计算方法，其特征在于，
其中，σ是斯忒藩-玻尔兹曼常数，εc是玻璃盖板的发射率，Tsky是有效天空温度，ΔA是光伏光热微元的受光面积。
6.根据权利要求1所述的一种非均匀辐照度下的聚光光伏光热系统性能计算方法，其特征在于，公式(10)中，光伏电池微元的串联电阻Rs(i,j)根据离散前后光伏电池阵列的I2R相等的原则进行计算：
nm·[I2(i,j)·Rs(i,j)]＝m[nI(i,j)]2Rs (12)
其中，I(i,j)是光伏电池微元的电流，Rs是单体光伏电池片的串联电阻，公式(12)化简为：
Rs(i,j)＝nRs (13)。
7.根据权利要求1所述的一种非均匀辐照度下的聚光光伏光热系统性能计算方法，其特征在于，根据光伏光热微元的光-热-电耦合模型能够得到光伏电池表面的温度分布，同时能够获得聚光光伏光热系统总体的热电性能参数的具体过程如下：
1)光伏电池阵列最大功率点电流的确定：根据串联效应，光伏电池阵列最大功率点电流由串联的单体光伏电池片最大功率点电流的最小值确定；每个单体光伏电池片最大功率点电流由宽度方向所有光伏电池微元的最大功率点电流叠加得到：
2)光伏电池阵列最大功率点电压的确定：根据并联效应，光伏电池阵列最大功率点电压由长度方向所有单体光伏电池片最大功率点电压叠加得到；每个单体光伏电池片最大功率点电压采用宽度方向所有光伏电池微元的最大功率点电压的平均值：
3)光伏电池阵列输出的最大电功率，等于光伏电池阵列最大功率点电流和电压之积：
Pm＝ImVm (17)
4)聚光光伏光热系统出口流体温度的确定：聚光光伏光热系统的出口流体温度由光伏光热组件冷却流道出口处宽度方向所有流体微元的出口温度的平均值确定：
5)根据步骤3)和4)，获得了光伏电池阵列输出的最大电功率以及聚光光伏光热系统出口流体温度之后，进一步得到聚光光伏光热系统的光电转换效率和光热转换效率：
其中，Gt是入射太阳总辐照度，Ain是聚光器入射光孔面积，
是冷却通道中流体的流量，Tf,in是流体进口温度。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810409878.9A CN108664720B (zh) | 2018-05-02 | 2018-05-02 | 一种非均匀辐照度下的聚光光伏光热系统性能计算方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810409878.9A CN108664720B (zh) | 2018-05-02 | 2018-05-02 | 一种非均匀辐照度下的聚光光伏光热系统性能计算方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108664720A CN108664720A (zh) | 2018-10-16 |
| CN108664720B true CN108664720B (zh) | 2020-04-28 |
Family
ID=63780430
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810409878.9A Active CN108664720B (zh) | 2018-05-02 | 2018-05-02 | 一种非均匀辐照度下的聚光光伏光热系统性能计算方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108664720B (zh) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111177884B (zh) * | 2019-12-06 | 2024-07-16 | 清华大学合肥公共安全研究院 | 基于温度场相似的辐照-高温-强风耦合分部实验方法 |
| CN112944702A (zh) * | 2021-02-07 | 2021-06-11 | 中国科学院重庆绿色智能技术研究院 | 太阳能热电联供系统投入产出比计算方法、装置及介质 |
| CN112944698B (zh) * | 2021-02-07 | 2023-01-24 | 中国科学院重庆绿色智能技术研究院 | 一种太阳能热电联供组件的瞬态热电输出计算方法及系统 |
| CN113446656B (zh) * | 2021-06-29 | 2021-11-09 | 天津滨电电力工程有限公司 | 功率-负荷匹配的光伏光热pv/t热电冷联供系统调控方法 |
| CN115615153B (zh) * | 2022-08-09 | 2024-06-07 | 珠海格力电器股份有限公司 | 干燥系统控制方法、装置及干燥系统 |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104598716A (zh) * | 2014-11-28 | 2015-05-06 | 楚雄师范学院 | 一种基于模型分析的聚光光伏/热热水太阳能系统设计方法 |
| CN104767483A (zh) * | 2015-03-31 | 2015-07-08 | 国家电网公司 | 一种光热特性耦合光伏电池温度预测方法 |
| CN105485939A (zh) * | 2015-12-14 | 2016-04-13 | 西安交通大学 | 一种太阳能聚光光伏光热联产系统热电输出性能的测量与计算方法 |
| CN105718422A (zh) * | 2016-03-07 | 2016-06-29 | 华东理工大学 | 一种中空光伏玻璃幕墙热性能计算方法 |
| CN106372350A (zh) * | 2016-09-12 | 2017-02-01 | 西安交通大学 | 评估二维聚光器对聚光光伏系统电性能影响的计算方法 |
| CN106651627A (zh) * | 2016-09-12 | 2017-05-10 | 西安交通大学 | 一种低倍聚光光伏光热系统的热电性能评估方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP6573809B2 (ja) * | 2014-09-16 | 2019-09-11 | アイメック・ヴェーゼットウェーＩｍｅｃ Ｖｚｗ | 太陽電池モジュールのシミュレーション |
-
2018
- 2018-05-02 CN CN201810409878.9A patent/CN108664720B/zh active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104598716A (zh) * | 2014-11-28 | 2015-05-06 | 楚雄师范学院 | 一种基于模型分析的聚光光伏/热热水太阳能系统设计方法 |
| CN104767483A (zh) * | 2015-03-31 | 2015-07-08 | 国家电网公司 | 一种光热特性耦合光伏电池温度预测方法 |
| CN105485939A (zh) * | 2015-12-14 | 2016-04-13 | 西安交通大学 | 一种太阳能聚光光伏光热联产系统热电输出性能的测量与计算方法 |
| CN105718422A (zh) * | 2016-03-07 | 2016-06-29 | 华东理工大学 | 一种中空光伏玻璃幕墙热性能计算方法 |
| CN106372350A (zh) * | 2016-09-12 | 2017-02-01 | 西安交通大学 | 评估二维聚光器对聚光光伏系统电性能影响的计算方法 |
| CN106651627A (zh) * | 2016-09-12 | 2017-05-10 | 西安交通大学 | 一种低倍聚光光伏光热系统的热电性能评估方法 |
Non-Patent Citations (2)
| Title |
|---|
| "Design and performance study of truncated CPC by eliminating multiple reflections of solar radiation in hybrid CPV/T system: Highest and lowest truncation position ";Huling Xie.et al;《Solar Energy》;20161015;第136卷;Page217-225 * |
| "低倍聚光光伏光热系统热性能分析";高阳等;《工程热物理学报》;20160331(第3期);第624-628页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108664720A (zh) | 2018-10-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108664720B (zh) | 一种非均匀辐照度下的聚光光伏光热系统性能计算方法 | |
| Yousef et al. | Performance investigation of low–Concentration photovoltaic systems under hot and arid conditions: Experimental and numerical results | |
| Slimani et al. | A detailed thermal-electrical model of three photovoltaic/thermal (PV/T) hybrid air collectors and photovoltaic (PV) module: Comparative study under Algiers climatic conditions | |
| Radwan et al. | The influence of microchannel heat sink configurations on the performance of low concentrator photovoltaic systems | |
| Corbin et al. | Experimental and numerical investigation on thermal and electrical performance of a building integrated photovoltaic–thermal collector system | |
| Coventry | Performance of a concentrating photovoltaic/thermal solar collector | |
| Zhang et al. | Thermal and electrical performance of low-concentrating PV/T and flat-plate PV/T systems: A comparative study | |
| Chemisana et al. | Characterization of a photovoltaic-thermal module for Fresnel linear concentrator | |
| Abd El-Samie et al. | Three-dimensional numerical investigation of a hybrid low concentrated photovoltaic/thermal system | |
| Gorouh et al. | Thermal modelling and experimental evaluation of a novel concentrating photovoltaic thermal collector (CPVT) with parabolic concentrator | |
| Simonetti et al. | Development and validation of a comprehensive dynamic mathematical model for hybrid PV/T solar collectors | |
| Assoa et al. | Dynamic study of a new concept of photovoltaic–thermal hybrid collector | |
| Wang et al. | Design and performance study on a large-scale hybrid CPV/T system based on unsteady-state thermal model | |
| Haiping et al. | Numerical simulation and experimental analysis of an LCPV/T system under real operating conditions | |
| Kazem et al. | Comparison and evaluation of solar photovoltaic thermal system with hybrid collector: An experimental study | |
| Gholami et al. | Evaluation of optimal dual axis concentrated photovoltaic thermal system with active ventilation using Frog Leap algorithm | |
| Li et al. | A coupled optical-thermal-electrical model to predict the performance of hybrid PV/T-CCPC roof-top systems | |
| Khanalizadeh et al. | Experimental investigation of a PV/T system containing a TEG section between water-based heat exchanger and air-based heat sink | |
| Slimani et al. | Analysis of thermal and electrical performance of a solar PV/T air collector: Energetic study for two configurations | |
| Chen et al. | Comparative investigation on photo-thermal performance of both compound parabolic concentrator and ordinary all-glass evacuated tube absorbers: an incorporated experimental and theoretical study | |
| Liu et al. | Experimental and numerical analyses of parameter optimization of photovoltaic cooling system | |
| Dey et al. | Numerical modeling and performance assessment of elongated compound parabolic concentrator based LCPVT system | |
| Hussain et al. | Parametric performance analysis of a concentrated photovoltaic co-generation system equipped with a thermal storage tank | |
| Baig et al. | Modelling and experimental analysis of a seasonally tracked V-trough PV/T system in India | |
| Lv et al. | Preliminary experiment and performance evaluation of a terrestrial solar thermoelectric generators under fluctuant solar radiation |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |