直接式太阳能干燥系统的热性能分析及应用

郝文刚，陆一锋，赖艳华※,，于洪文，吕明新

（1. 山东大学能源与动力工程学院，济南 250061；2. 北京工业大学材料学院，北京 100124；3. 山东大学苏州研究院，苏州 215028；4. 山东桑乐太阳能有限公司，济南 250014）

0 引 言

太阳能是一种取之不尽用之不竭的可再生能源[1]，其在粮食[2-4]、农副产品[5-8]等物料干燥方面应用已经引起很多学者的关注。直接式太阳能热利用技术在干燥领域中的应用最早始于976年Everitt and Stanley提出来具有太阳能透光性的盒子形状干燥器[9]。直接式太阳能干燥（direct solar drying）系统主要由透明盖板、保温层、进风口、出风口等组成的干燥箱[10]。当太阳光通过透明玻璃盖板入射到干燥箱内时，干燥箱内空气温度升高加热干燥物料，同时在热压或风压的作用下室外空气由干燥箱的进口处进入干燥箱内，干燥箱内高温高湿空气通过干燥箱的出口处排到室外。该技术的结构比较简单且制作价格便宜，具有保护干燥物料免受灰尘、雨水、露水等破坏的优点；但其也有一些缺点，例如干燥物料表面过热、干燥物料品质差和干燥能力有限等[11-13]。

为克服直接式太阳能干燥技术的缺点，国外相关学者设计不同形式的直接式太阳能干燥系统研究其热性能、干燥能力等方面问题[14-17]，Ahmad[18]设计搭建了一个覆盖透明保温材料的圆柱形状集热干燥装置，相比开放式太阳能干燥，该装置能够将室内空气温度提升 10 ℃左右。Dissa等[19]研究芒果薄层在直接式太阳能干燥器的干燥特性，其最终的含水率为66.32%，研究结果表明干燥速率和有效性随着干燥时间延长而降低；Eke等[20]学者研制一种新型用于干燥玉米的直接式太阳能器，干燥玉米物料的湿基含水率从29%减少到12%，与开放式太阳能干燥相比，干燥时间节约55%，新研制的干燥器与开放式太阳能干燥方法的干燥效率分别是45.6%和22.7 %。国内关于直接式太阳能干燥技术的研究在文献中还未检索到，但是关于太阳能干燥技术设备的研究有很多[21-22]。李汴生等[23]自行研制太阳能干燥设备，以干湿梅为原料进行自然对流干燥，试验测试结果表明其干燥效率为61.6%，干燥时间为14 h，比开放式太阳能干燥时间缩短72%。综上所述，目前国内外研究大多数集中于干燥系统干燥能力[24-25]、热利用效率[22]、物料干燥特性[26-29]等方面，在预测干燥室内物料表面温度的研究较少，尤其是针对直接式太阳能干燥系统的。因此本文自行设计搭建一种直接式太阳能干燥系统，以红薯为干燥物料进行自然对流干燥试验测试，通过测试相关试验数据，分析该系统的热效率及干燥速率；同时建立该系统的热性能数学模型，用于预测干燥过程中物料表面的温度，以期提高该系统干燥物料的品质。

1 直接式太阳能干燥系统试验

1.1 试验台结构

为建立并验证干燥物料在直接式太阳能干燥系统中表面温度动态预测模型，本文设计搭建一种直接式太阳能干燥系统，图1为直接式太阳能干燥系统结构示意图。

图1 直接式太阳能干燥系统示意图Fig.1 Schematic diagram of indirect solar drying system

该系统主要由空气进口、空气出口、物料托盘、透光玻璃板及不锈钢铁板等组成。干燥箱的框架为铝合金，同时在干燥箱不锈钢板内表面均匀涂抹黑色吸热材料。干燥箱屋顶玻璃的倾斜角度为26.7°，干燥箱的长×宽×高为 1 000 mm×800 mm×800 mm，干燥箱隔热材料采用20 mm厚的聚氨酯保温棉，以防止热量的损失，干燥箱的透光面采用透射率为 90%的普通玻璃。干燥箱内放置干燥物料托盘。物料托盘采用孔隙均匀分布的铁丝网，方便热空气向上流动。

1.2 材料与方法

1）试验原料与测试方案

本试验采用的干燥物料为新鲜红薯，购于京东商城，原产地为云南玉溪。试验前将红薯切成大小、厚度均匀一致的圆片，用天平将其分成质量相等的 2份，一份用于室内干燥试验组，一份用于室外干燥对比组。试验测试时间为2017年5月7日上午8:30至下午17:00，试验组和对比组干燥物料的初始质量均为240 g。对比组采用开放式太阳能干燥方式，试验组和对比组均在相同的室外环境下进行干燥，对比组设计的目的是为与试验组进行物料干燥速率的对比。

2）试验仪器

太阳能辐射强度采用TBQ-2太阳能总辐射表（锦州阳光气象科技有限公司，量程：0～2 000 W/m2；精度：2%）进行测试，室内外空气温湿度由温湿度自计议174H，（德图仪器国际贸易（上海）有限公司，量程：温度（-20～70℃）、相对湿度（0%～100%）；精度：相对湿度（±3%）、温度（±0.5℃）），红薯表面温度采用T型热电偶（美国Omega公司，量程：-200～350℃；精度：±0.5℃）进行测试，温度传感器采集的数据通过数据采集器ADAM4117（研华科技（中国）有限公司）进行记录，记录时间步长为10 min，红薯质量测试仪器为电子天平（上海花潮电器有限公司，量程：0～10 kg；精度：±0.1 g）。

1.3 试验误差分析

为说明试验测试数据的准确性，本文对上述所用试验测试仪器的测试参数进行不确定度分析，假设x1，x2，x3,…, xn是n个相互独立的变量，其代表影响试验测试参数准确性的因素，R=f (x1, x2,…, xn)，代表试验测试参数准确性的函数关系式。则R的不确定度WR[30]见公式（1）。

利用公式（1）计算不同测试参数的不确定度值，如表1所示。

表1 测试参数的不确定度分析Table 1 Uncertainty analysis of testing parameters

利用公式（1）对本试验测试的总不确定度进行计算，计算其值为5.2%，该试验测试的不确定度在试验可接受的范围之内[16]。

2 数学模型及求解方法

2.1 假设条件

在保证模拟结果的准确性和方便求解的前提下，对直接式太阳能干燥箱的模型进行简化，作出以下假设[4]。1）干燥箱的铝合金、玻璃等材料的比热容量被忽略；2）干燥箱的各表面由于长、宽均远远大于其厚度，因此被按照一维热传导进行处理；

3）各表面的温度均按照均匀温度进行计算；4）干燥箱内整个换热过程是准稳态过程。

2.2 热平衡方程

基于上述假设条件，为准确建立直接式太阳能干燥箱与干燥物料之间热交换方程，本文基于热网络法建立直接式太阳能干燥系统的热阻及热能量流动图，如图 2所示。

2.2.1 干燥箱的热平衡方程

干燥箱内空气能量变化主要受进出干燥箱空气的能量损失、干燥箱内空气与物料之间热交换、通过干燥箱壁面向室外散失的热量的影响，因此干燥箱的热平衡方程如公式（2）所示[4]。

式中 hcr为干燥物料与干燥箱内空气综合换热系数，W/(m2·K)；Tc为干燥物料表面温度，K；Td为干燥箱内空气温度，K；Ta为室外空气温度，K；Ac为干燥物料表面积，m2；Mf为进出干燥箱空气的质量流量，kg/s；Cf为干燥箱内空气比热容，J/(kg·K)；U为不同方向干燥箱体的传热系数，W/(m2·K)；A为不同方向干燥箱体的面积，m2。

其中不同方位的传热系数和箱体面积可以由公式（3）计算所得。

式中下标S、R、W、E、N、F分别为南向、坡屋顶、西向、东向、北向、底部方位。

式中 hv,go-a为玻璃外表面与室外空气的对流换热系数，W/(m2·K)；dg为玻璃厚度，m；λg为玻璃的导热系数，W/(m2·K)；hv,gi-d为玻璃内表面与室内空气的对流换热系数，W/(m2·K)。

图2 直接式太阳能干燥系统热阻及热能量流动图Fig.2 Thermal resistance and thermal energy circuit diagram of direct solar drying system

式中 hv,wo-a为非透明干燥箱体与室外空气的对流换热系数，W/(m2·K)；dw为干燥箱箱体不同组成的厚度，m；λw为干燥箱箱体不同组成的导热系数，W/(m2·K)；hv,wi-d为非透明干燥箱体与室内空气的对流换热系数，W/(m2·K)。

式中At为干燥箱透光面积，m2。

式中Ant为干燥箱非透光面积，m2。

由公式（2）得到公式（8）

2.2.2 物料的热平衡方程

干燥物料的能量变化主要受入射到干燥箱内的太阳能、干燥物料与干燥箱内空气的热交换的影响，因此干燥物料的能量平衡方程如公式（9）所示。

式中Mc为干燥物料的质量，kg；Cc为干燥物料的比热容，J/kg·K)；I(i)为太阳能逐时辐射强度，W/m2；αc为干燥物料的太阳能吸收率；τg为玻璃的透光率。

将式（8）代入式（9）得到式（10）

2.3 模型参数的确定

2.3.1 综合换热系数

干燥物料与干燥箱内空气综合换热系数见式（11）。

式中 hv,d-c为干燥物料与干燥箱内空气对流换热系数，W/(m2·K)；hew为干燥物料水分的蒸发换热系数，W/(m2·K)。

干燥物料与干燥箱内空气对流换热系数见式12[4]。

式中vd为干燥箱内空气流速，m/s。

干燥物料水分蒸发换热系数计算如式（13）所示[31]。

不同温度下的水蒸气分压力计算如式（14）、式（15）所示。

非透明干燥箱体与室外空气的对流换热系数为[32]

式中va为干燥箱外空气流速，m/s。

玻璃与室外空气的对流换热系数为[4]非透明干燥箱体与室内空气的对流换热系数为[4]

玻璃与室内空气的对流换热系数为[4]

2.3.2 初始值和物性参数

本文数学模型求解所需的初始值和物性参数[3,33]如表2所示。

表2 模型中初始和物性参数计算值Table 2 Initial and physical parameter values in model

2.4 模型求解方法

本模型建立的目的是根据监测室外气象参数、干燥箱内空气温湿度进而预测干燥物料表面的温度。数学模型式（10）属于显性常微分方程，其中式（11）～（19）为其参数求解边界条件，通过欧拉差分法将式（10）进行离散并利用 MATLAB2014a编制模型求解程序。室外气象参数、干燥箱内空气温湿度、红薯的表面温度初始值均采用试验数据。

大鼠进行脑缺血后瞳孔散大，眼睛由红变白，缺血前期呼吸加快变浅，缺血后期呼吸减慢变深；再灌注期间瞳孔由散大转为缩小，眼睛由苍白转为红润，且呼吸加深加快慢慢转为正常平稳呼吸。

2.5 热性能评价指标

1）系统所获得的热能量。

式中Qu,th为干燥箱所获得的能量，J；i为逐时时刻。

2）太阳能利用效率。

式中 Id(i)为入射到干燥箱内太阳能辐射强度，W/m2；η

为干燥箱的太阳能利用效率，%。

3 模型验证

为了验证模拟结果的准确性，本文采用决定系数 R2和均方根误差RMSE作为其判断的依据，其计算公式如式（22）、（23）所示。

式中 Xi为试验值；Yi为模拟值；X为试验平均值；Y为模拟平均值。

直接式太阳能干燥系统中红薯表面温度模拟值与试验值变化曲线如图 3所示，干燥过程中红薯的表面温度变化范围为，其峰值出现在13:00时刻。试验值与模拟值之间的决定系数R2为0.98，均方根误差为1，说明该模型能够较好的根据相关参数预测出红薯表面温度，可以有效的防止干燥物料表面过热，提高干燥品质。

图3 物料表面温度模拟值与试验值对比Fig.3 Comparison curves between simulated and experimental drying materials surface temperature

4 系统性能测试

试验测试的室内外太阳能辐射强度及干燥物料质量逐时变化曲线如图4所示。

图4 太阳辐射强度及干燥物料质量逐时变化曲线Fig.4 Curve of hourly variation of solar radiation and drying material mass

在试验测试期间，室内外太阳能辐射强度的变化范围分别为 398.91～979.9 W/m2、192～864.1 W/m2，室内太阳能辐射强度最高值出现在 12:30，比室外晚 1 h，干燥箱的透射比范围为51.7%～89.6%，其最大透射比值也出现在12:30，主要是因为干燥箱透光面位于坡屋顶及南面方位，太阳光能够较多的直射进干燥箱内。从图 4中可以得出，直接式太阳能干燥比开放式太阳能干燥（对照组）的平均干燥速率高 7.7 g/h，红薯干燥的最终质量低7 g，说明直接式太阳能干燥的干燥能力比开放式太阳能干燥能力高。15:00后直接式太阳能干燥速率突然增大的原因是由于干燥箱作为一种类似的被动式建筑，具有一定的蓄热能力，对室外环境的影响有一定的抗延迟性，导致室内干燥物料质量下降幅度相对于室外突然增大。

室内外空气温湿度逐时变化曲线如图 5所示，干燥箱内空气温度最高值出现在13:00，与图4中太阳能辐射强度峰值出现时刻晚半个小时，说明干燥箱体对室外环境变化具有一定的延迟性和抗干扰性。干燥箱内空气温度范围在50～70℃的分布时间段占总干燥时间的80%，可以说明此搭建的直接式太阳能干燥系统能够为大多数干燥物料提供较适宜的干燥温度。

图5 室内外空气温湿度逐时变化曲线Fig.5 Curve of hourly variation of indoor and outdoor temperature and humidity

直接式太阳能干燥系统所获得的热能及热效率逐时变化曲线如图 6所示。干燥系统所获得的总热能量为3.92 kW·h，其平均太阳能热利用效率为21.23%。

图6 太阳能利用效率及所获得的热能变化曲线Fig.6 Curve of solar energy utilization efficiency and thermal energy obtained of drying system

5 讨 论

1）在建立物料表面温度预测模型时，忽略了干燥箱的铝合金、玻璃等材料的比热容量，因为物料在干燥的过程中主要受室内热空气及太阳能辐射两方面的影响。本文的目的是预测干燥物料表面温度。因此在建立热网络方程时，为方便计算忽略干燥箱各非透光面吸收太阳辐射的得热量。干燥箱体吸收的太阳能辐射得热量最终也是转化为室内空气温升和散失到室外。而本计算采用的监测的室内空气逐时温度作为已知条件，可以说包括部分干燥箱体的太阳能辐射得热量。

2）本文中试验组和对比组的2条干燥曲线差异较少，但是二者间平均干燥速率相差7.7g/h，造成这种现象的原因可能是试验过程中采用的干燥红薯质量较少。

3）在对直接式太阳能干燥系统的应用效果进行试验分析时，发现直接式太阳能干燥系统是否应该增加风机来调节干燥箱内空气质量流量应值得深入探讨与研究。在众多前人研究中发现干燥箱内空气质量流量是影响物料干燥速率的重要因素[6]，但通常在间接式或混合式太阳能干燥系统应用比较多，目的是为提高集热器集热效率。本文设计的直接式太阳能干燥系统是纯依靠被动式太阳能技术营造物料适宜的干燥环境，如果加入风机来控制干燥箱内空气质量流量可以在中午室内高温时短时间运行，目的是将室内温度降低在适宜的物料干燥温度范围内，但是如果长期运行会造成直接式太阳能干燥系统的热利用效率下降，同时也会增加该系统的初投资和运行费用。因此直接式太阳能干燥系统中是否加入风机可以从全生命周期评价方面作进一步研究。在本文中不作考虑。

6 结 论

通过对搭建的直接式太阳能干燥系统进行实验测试及理论分析，得出以下结论：

1）干燥箱的透射比范围为51.7%～89.6%，干燥箱内空气温度最高值出现的时刻比太阳能辐射强度峰值出现时刻晚半个小时，说明干燥箱体对室外环境变化具有一定的延迟性和抗干扰性。直接式太阳能干燥比开放式太阳能干燥的平均干燥速率高 7.7 g/h，干燥红薯最终质量低7 g，说明直接式太阳能干燥的干燥能力比开放式太阳能干燥强。

2）干燥箱内空气温度范围在50～70℃的分布时间段占总干燥时间的 80%，可以说明此搭建的直接式太阳能干燥系统能够为大多数干燥物料提供较适宜的干燥温度。干燥系统所获得的总热能量为3.92 kW·h，其平均太阳能热利用效率为21.23%。

3）通过对搭建的直接式太阳能干燥系统建立热性能动态模型并求解，试验值与模拟值之间的决定系数R2为0.98，均方根误差为1，说明该模型能够较好的根据相关参数预测出干燥物料表面温度。同时也说明建立的热性能动态数学模型具有一定普遍适用性。

[1] Panwar N L, Siva R V, Ranjan K R, et al. Sustainable development with renewable energy resources: A review[J].World Review of Science, Technology and Sustainable Development, 2013, 10(4): 163－184.

[2] Dissa A O, Bathiebo J, Kam S, et al. Modelling and experimental validation of thin layer indirect solar drying of mango slices[J]. Renewable Energy, 2009, 34(4): 1000－1008.

[3] Ali I, Abdelkader L, El H B, et al. Solar convective drying in thin layers and modeling of municipal waste at three temperatures[J]. Applied Thermal Engineering, 2016, 108:41－47.

[4] Tiwari S, Tiwari G N, Al-Helal I M. Performance analysis of photovoltaic-thermal (PVT) mixed mode greenhouse solar dryer[J]. Solar Energy, 2016, 133: 421－428.

[5] Orbegoso E M, Saavedra R, Marcelo D, et al. Numerical characterisation of one-step and three-step solar air heating collectors used for cocoa bean solar drying[J]. Journal of Environmental Management, 2017, 203(1): 1080－1094.

[6] Dissa A O, Bathiebo D J, Desmorieux H, et al. Experimental characterisation and modelling of thin layer direct solar drying of Amelie and Brooks mangoes[J]. Energy, 2011,36(5): 2517－2527.

[7] Kareem M W, Habib Khairul, Sopian K, et al. Multi-pass solar air heating collector system for drying of screw-pine leaf (Pandanus tectorius)[J]. Renewable Energy, 2017, 112:413－424.

[8] Condorí M, Duran G, Echazú R, et al. Semi-industrial drying of vegetables using an array of large solar air collectors[J].Energy for Sustainable Development, 2017, 37: 1－9.

[9] Sharma A, Chen C R, Vu L N. Solar-energy drying systems:A review[J]. Renewable and Sustainable Energy Reviews,2009, 13(6-7): 1185－1210.

[10] Vijayan S, Arjunan T V, Kumar A. Mathematical modeling and performance analysis of thin layer drying of bitter gourd in sensible storage based indirect solar dryer[J]. Innovative Food Science & Emerging Technologies, 2016, 36: 59－67.

[11] Kumar M, Sansaniwal S K, Khatak P. Progress in solar dryers for drying various commodities[J]. Renewable and Sustainable Energy Reviews, 2016, 55: 346－360.

[12] Mekhilef S, Faramarzi S Z, Saidur R, et al. The application of solar technologies for sustainable development of agricultural sector[J]. Renewable and Sustainable Energy Reviews, 2013,18: 583－594.

[13] Pirasteh G, Saidur R, Rahman S M A, et al. A review on development of solar drying applications[J]. Renewable and Sustainable Energy Reviews, 2014, 31: 133－148.

[14] Banout J, Ehl P, Havlik J, et al. Design and performance evaluation of a Double-pass solar drier for drying of red chilli(Capsicum annum L.)[J]. Solar Energy, 2011, 85(3): 506－515.

[15] Dissa A O, Desmorieux H, Bathiebo J, et al. Convective drying characteristics of Amelie mango (Mangifera Indica L.cv. ‘Amelie’) with correction for shrinkage[J]. Journal of Food Engineering, 2008, 88(4): 429－437.

[16] ELkhadraoui A, Kooli S, Hamdi I, et al. Experimental investigation and economic evaluation of a new mixed-mode solar greenhouse dryer for drying of red pepper and grape[J].Renewable Energy, 2015, 77: 1－8.

[17] Fudholi A, Sopian K, Othman M Y, et al. Energy and exergy analyses of solar drying system of red seaweed[J]. Energy and Buildings, 2014, 68: 121－129.

[18] Ahmad N T. Agricultural solar air collector made from low-cost plastic packing film[J]. 2001, 23(3): 663－671.

[19] Dissa A O, Bathiebo D J, Desmorieux H, et al. Experimental characterisation and modelling of thin layer direct solar drying of Amelie and Brooks mangoes[J]. Energy, 2011,36(5): 2517－2527.

[20] Eke A, Arinze E. Natural convection mud type solar dryers for rural farmers[J]. Niger J Technol Dev, 2011, 8(2): 92－106.

[21] 张璧光，伊松林. 我国太阳能干燥现状及在热利用中应有的地位[J]. 干燥技术与设备，2012，10(1)：32－37.Zhang Biguang, Yi Songlin. The present state of the solar drying in china and the status which it should have in thermal application [J]. Drying Technology& Equipment, 2012, 10(1):32－37. (in Chinese with English abstract)

[22] 王云峰，李明，王六玲，等. 太阳能干燥装置性能及三七干燥效果[J]. 农业工程学报，2010，26(10)：377－383.Wang Yunfeng, Li Ming, Wang Liuling, et al. Performance of solar dryer and drying effect for panax notoginseng[J].Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2010, 26(10): 377－383. (in Chinese with English abstract)

[23] 李汴生，申晓曦，刘伟涛. 太阳能干燥设备的干燥实验及数学建模[J]. 华南理工大学学报(自然科学版)，2011，39(4)：115－120.Li Biansheng, Shen Xiaoxi, Liu Weitao. Experimental investigation and mathematical modeling of solar drying equipment[J]. Journal of South China University of Technology (Nature Science Edition), 2011, 39(4): 115－120.(in Chinese with English abstract)

[24] 王伟华，王海，何思鲁，等. 南美白对虾太阳能干燥能耗参数优化及中试[J]. 农业工程学报，2016，32(18)：271－278.Wang Weihua, Wang Hai, He Silu, et al. Parameter optimization for energy consumption of solar drying of Penaeus vannamei and pilot scale test[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2016, 32(18): 271－278. (in Chinese with English abstract)

[25] 尹丽洁，王中慧，陈德珍. 污泥冬季太阳能干燥实验和过程预测[J]. 太阳能学报，2013，34(12)：2088－2094.Yin Lijie, Wang Zhonghui, Chen Dezhen. Experiment and prediction of solar drying process of sewage sludge in winter[J]. Journal of solar energy, 2013, 34(12): 2088－2094. (in Chinese with English abstract)

[26] 李菁，萧夏，蒲晓璐，等. 紫薯热风干燥特性及数学模型[J]. 食品科学，2012，33(15)：90－94.Li Jing, Xiao Xia, Pu Xiaolu, et al. Characteristics and mathematical model of hot-air drying for purple sweet potato[J]. Food Science, 2012, 33(15): 90－94. (in Chinese with English abstract)

[27] 凌德力，李明，罗熙，等. 基于槽式聚光太阳能供热的烟丝干燥特性研究[J]. 太阳能学报，2015,36(2)：460－466.Ling Deli, Li Ming, Luo Xi, et al. Study on drying characteristics of cut tobacco based on trough concentrating solar heating[J]. Journal of solar energy, 2015, 36(2): 460－466. (in Chinese with English abstract)

[28] 马煜，李明，魏生贤，等. 三七薄层干燥特性研究[J]. 太阳能学报，2012,33(6)：937－943.Ma Yu, Li Ming, Wei Shengxian, et al. Lamina drying investigation of notoginseng [J]. Journal of solar energy,2012, 33(6): 937－943. (in Chinese with English abstract)

[29] 孟岳成，王君，房升，等. 熟化红薯热风干燥特性及数学模型适用性[J]. 农业工程学报，2011，27(7)：387－392.Meng Yuecheng, Wang Jun, Fang Sheng, et al. Drying characteristics and mathematical modeling of hot air drying of cooked sweet potatoes [J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2011, 27(7): 387－392. (in Chinese with English abstract)

[30] Fudholi A, Sopian K, Yazdi M H, et al. Performance analysis of solar drying system for red chili[J]. Solar Energy, 2014, 99:47－54.

[31] Putranto A, Chen X, Xiao Z, et al. Mathematical modeling of intermittent and convective drying of rice and coffee using the reaction engineering approach (REA)[J]. Journal of Food Engineering, 2011, 105(4): 638－646.

[32] Reddy K S, Aravindhan S, Mallick T K. Investigation of performance and emission characteristics of a biogas fuelled electric generator integrated with solar concentrated photovoltaic system[J]. Renewable Energy, 2016, 92: 233－243.

[33] Jain D, Tiwari G N. Thermal aspects of open sun drying of various crops[J]. Energy, 2003, 28(1): 37－54.