research article a novel parabolic trough concentrating...
TRANSCRIPT
-
Research ArticleA Novel Parabolic Trough Concentrating Solar Heating forCut Tobacco Drying System
Jiang Tao Liu, Ming Li, Qiong Fen Yu, and De Li Ling
Solar Energy Research Institute, Yunnan Normal University, Kunming 650092, China
Correspondence should be addressed to Ming Li; [email protected]
Received 31 October 2013; Accepted 4 February 2014; Published 26 March 2014
Academic Editor: Adel M. Sharaf
Copyright © 2014 Jiang Tao Liu et al. This is an open access article distributed under the Creative Commons Attribution License,which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
A novel parabolic trough concentrating solar heating for cut tobacco drying system was established. The opening width effect ofV type metal cavity absorber was investigated. A cut tobacco drying mathematical model calculated by fourth-order Runge-Kuttanumerical solution method was used to simulate the cut tobacco drying process. And finally the orthogonal test method was usedto optimize the parameters of cut tobacco drying process. The result shows that the heating rate, acquisition factor, and collectorsystem efficiency increase with increasing the opening width of the absorber. The simulation results are in good agreement withexperimental data for cut tobacco drying process. The relative errors between simulated and experimental values are less than 8%,indicating that thismathematical model is accurate for the cut tobacco airflow drying process.The optimumpreparation conditionsare an inlet airflow velocity of 15m/s, an initial cut tobacco moisture content of 26%, and an inlet airflow temperature of 200∘C.The thermal efficiency of the dryer and the final cut tobacco moisture content are 66.32% and 14.15%, respectively.The result showsthat this parabolic trough concentrating solar heating will be one of the heat recourse candidates for cut tobacco drying system.
1. Introduction
Solar energy currently represents the most abundant inex-haustible, nonpolluting, and free energy resources that couldhave a positive meaning in alleviating the global energyshortage and environmental pollution. There are abundantsolar energy resources and good economic tobacco industryin Yunnan Province, located in the southwest of China.However, drying cut tobacco of the production process incigarette factory needs a large amount of heat energy (150∘C∼300∘C). If abundant solar energy resources are used for dryingcut tobacco of the production process in cigarette factory, wecan not only save a large amount of fossil energy but alsowiden the field of utilization of solar energy. In this paper,solar energy drying cut tobacco process will be investigated.
However, the solar energy is intermittent in nature andthat received on earth is of small flux density due toatmospheric scattering and abortion, making it necessaryto use large surfaces to collect solar energy for drying cuttobacco process utilization. The major technologies usedfor solar energy conversion to heat are thermal processes
comprising of solar collectors. There are two basic types ofsolar collectors, the flat-plate and the concentrating solarcollectors. The flat-plate collector has advantage of absorbingboth beam and diffused radiation and, therefore, still func-tions when beam radiation is cut off by the cloud. The areaabsorbing solar radiation is the same as the area interceptingsolar radiation. However, flat-plate collectors are designedfor applications requiring energy delivered at temperaturesquite lower than 100∘C, indicating that they are not suitablefor drying cut tobacco process. The concentrating collectorutilizes optical systems like reflectors, refractors, and soforth to increase the intensity of solar radiation incidentson energy-absorbing surfaces. It has the main advantage ofgenerating high temperatures andmay be used for drying cuttobacco process.
Parabolic trough solar collector (PTC) is one type ofthe intermediate/high concentrating collectors, being widelyused in the field of electricity generation, air-conditioning,heating, desalination, and so on. Conventional receiver usedin PTC is the vacuum tube. Because of different expandingcoefficients between the metallic tube and the glass tube, the
Hindawi Publishing CorporationInternational Journal of PhotoenergyVolume 2014, Article ID 209028, 10 pageshttp://dx.doi.org/10.1155/2014/209028
-
2 International Journal of Photoenergy
cost of vacuum tube receiver is relatively high, and its usefullife is short.Therefore,many researches focus on investigationof cavity absorbers. Boyd et al. [1] firstly developed a receiverthat consists of an annular cylindrical tube with an apertureparallel to the axis of the cylinder.The aperture is illuminatedby means of a focusing concentrator, such as a lens, and theradiation energy is scattered inside the cavity, absorbed by thewalls and transmitted, as heat, to the working fluid flowingaxially inside the annulus. To reduce radiative losses the tubeis surrounded by a layer of a thermal insulator. Principaladvantages of this design lie in exclusive use of commonlyavailable materials and current fabrication technology. Inorder to minimize the heat loss, this annular cylindricaltube was improved by Barra and Franceschi [2]. Eight smallcylindrical tubes were scattered inside the big cylindrical tubeand used for fluid transferring. Qiaoli et al. [3] proposedthe scheme of tightly connecting cluster and the inner wallof pipe. Moreover, Zhang et al. [4] established a series ofabsorbers with triangular, circular, hemispherical and square,and simulated the optical performance and thermal prop-erty using the software programs TracePro and Fluent. Theresults showed that the triangular cavity absorber structureowned the best optical performance and thermal property,which could gain collecting efficiency above 40%, under thecondition of collecting temperature 150∘C. Based on theseabove studies, a novel cavity absorber was investigated bythis paper. PTC system with V type metal cavity absorbersimilar to triangle structure would be used for drying cuttobacco process and the thermal loss is constrained bythe cavity radiation. In order to lower the absorber cost,aluminum alloy was selected as the material of absorb-er.
Because of the coupling effect among the heat transfer,mass transfer, and energy, air drying process became com-plex. In order to analyze air drying cut tobacco process, sim-ulation model was discussed. Fukuchi et al. [5] regarded thecut tobacco as equal volume sphere and simulated themotioncharacteristics, indicating the simulation results agreed wellwith the experimental results. The model equivalent spherewas established to simulate drying process of superheatedsteam by Pakowski et al. [6]. The result showed the mois-ture content of shredded tobacco simulated and the outlettemperature simulated agreed well with the experimentalvalues. Based on these simulation results, simulation modelof air drying cut tobacco process was studied and the processparameters of process of drying cut tobacco were optimizedin this research.
In this paper, parabolic trough concentrating systemwith V type metal cavity absorber was used for providingheat for cut tobacco drying system. In order to reduce theincident light energy loss and prevent cavity absorber fromdeformation, the width of V type metal cavity absorberwas calculated, ensuring that stable heat energy providedby the PTC system is enough for cut tobacco dryingprocess. Moreover, a mathematical model of cut tobaccodrying process was calculated by using fourth-order Runge-Kutta numerical solution method and the accuracy of themodel compared with the experimental results was ana-lyzed.
Table 1: Main parameters of cut tobacco drying system.
Parabolic trough concentrator mirrorReflectivity 0.9Opening width 3.2mLighting area 30m2
Length 10mMetallic V cavity absorber
Opening width 4 cm, 7 cmInternal surface coating Aluminum anodic oxide filmAbsorptivity 0.85∼0.9Emissivity 0.1∼0.2
Variable diameter drying pipePart 1 Height: 1.2m, diameter: 200mmPart 2 Height: 2m, diameter: 300mmTotal height 3.25m
2. Materials and Methods
2.1. Materials. Before solar energy drying process, the cuttobacco samples were pretreated by humidifier and equili-brated water content for 24 h in conditions of temperature20±2
∘Cand relative humidity 60±2%.Moisture content of cuttobacco samples was determined by YC/T31-1996 standardtest methods [7].
2.2. Experimental Equipment
2.2.1. Parabolic Trough Concentrating Solar Heating for CutTobacco Drying System. Diagram of parabolic trough con-centrating solar heating for cut tobacco drying system wasshown in Figure 1(a). The main parameters of the design aregiven in Table 1. The parabolic trough concentrating solarheating for cut tobacco drying process was illustrated asfollows. Firstly, when sunlight arrived at the parabolic troughconcentrator mirror, the solar tracker was adjusted to makethemetal V cavity body coincidewith focal length of the PTC.Because of photothermal conversion, heat transfer oil wasused to exchange heat with the cavity body and then heatedto the predetermined temperature. Secondly, the circulatingfan was opened and the hot air exchanged with heat transferoil was transferred into the drying pipe. When temperatureof the hot air reached target temperature in the drying pipe,cut tobacco pretreated was put into the drying pipe for dryingprocess. Thirdly, when finishing the drying process, productof cut tobacco dried was followed with airflow into thewhirlwind separator and separated from the airflow. Fourthly,part of airflow was recycled into next drying process and theleft airflow was supplied by the valve S2. The electric heaterwas used to keep the output temperature of heat conductionoil stable when solar irradiance changed.
2.2.2. Test Schematic of Energy Flux Density Distribution.In order to analyze the effect of cavity opening width onenergy flux density distribution, the test schematicwas shownin Figure 2. As shown in Figure 2, the width and the focallength of parabolic trough concentrator mirror were 3.2m
-
International Journal of Photoenergy 3
M
P
Vent line
Fresh steam
Product in
Product out
1
2
3
4
5
6
78
9
S1
S2
Sunlight
(1) V-shaped metal cavity(2) Cyclone(3) Circulating fan(4) Air flow meter(5) Variable diameter drying pipe
(6) Heat exchanger(7) Electric heater(8) Oil pump(9) Collecting mirror
(a) Diagram of parabolic trough concentrating solar heating for cut tobacco dryingsystem
(b) Photograph of parabolic trough concen-trating solar heating for cut tobacco dryingsystem
Figure 1: Parabolic trough concentrating solar heating for cut tobacco drying system.
and 117 cm, respectively.There were a width of 10 cm Lamberttarget, neutral density attenuator (ND-filter), CCD industrialcameras, computers and solar tracking control system, andother parts. Test method was detailed by our previous work[8].
3. Cut Tobacco Drying Mathematical Model
3.1. Assumption. In this work, a relatively simple one-dimensional model [9, 10] was adopted assuming no radialor axial dispersion of flow. The shape of cut tobacco particleswas considered as equal volume sphere and had a uniformsize. Diameter and density variations of cut tobacco particleswere ignored in the drying process. The wall of drying pipehas good insulation and there was no heat exchange with theatmospheric environment.
3.2. Mathematical Model of the Drying Process. Themomen-tum equation of cut tobacco particle was as follows:
𝑑V𝑝
𝑑𝑧
=
3
4
𝜉
𝜌
𝑔(V𝑔− V𝑝)
2
V𝑝𝑑
𝑝𝜌
𝑝
− (1 −
𝜌
𝑔
𝜌
𝑝
)
𝑔
V𝑝
−
0.1
𝐷
.
(1)
The relationship between 𝜉 and Re in (1) was shown inTable 2 [11].
PTC
ComputerData line
ND-filter
Sunlight
x
Lamberttarget
z
CCD camera
O
Sunlight
Figure 2: The measurement schematic of focal line energy fluxdensity distribution.
The movement equation of air stream was as follows:
𝑑V𝑔
𝑑𝑧
= −
𝑑𝑃
𝜌
𝑔V𝑔𝑑𝑧
−
𝑔
V𝑔
−
3𝑁𝜉(V𝑔− V𝑝)
2
4V𝑔𝑑
𝑝
−
𝜆V𝑔
2𝐷
.
(2)
The parameter 𝑁 of (2) was the number of particles perunit volume in the drying pipe and identified according to theexperimental conditions.
-
4 International Journal of Photoenergy
Table 2: The relationship between 𝜉 and Re.
Parameter Numerical valueRe 0-1 1–500 500–15000𝜉 24/Re 24/Re0.5 0.44
The differential equations of cut tobacco particles andairflow humidity distribution were as follows:
𝐺
𝑝𝑑𝑋 + 𝐺
𝑔𝑑𝑌 = 0 (3)
𝐺
𝑃𝑑𝑋 +𝑊
𝑑𝑑𝐴 = 0. (4)
The equation of dry mass transfer rate was as follows:
𝑊
𝑑= 𝐾
𝑦(𝑌eq − 𝑌) . (5)
The differential equation of cut tobacco particles humiditywas as follows:
𝑑𝑋
𝑑𝑧
= −
6𝐾
𝑦(𝑌eq − 𝑌) (1 + 𝑋)
𝜌
𝑝𝑑
𝑝V𝑝
.(6)
The mass transfer coefficient 𝐾𝑦, the Schmidt number
Sc, and the Prandtl number Pr in (6) were associated andcalculated according to the correlation proposed by literature[12]:
ℎ
𝐾
𝑦
= 𝐶
𝑔(1 + 𝑌) (
ScPr
)
0.56
Pr =(𝐶
𝑔]𝑔𝜌
𝑔)
𝜆
𝑔
,
Sc =]𝑔
𝐷V𝑎.
(7)
Combining (3) and (6), the differential equation of airflowand humidity could be obtained as
𝑑𝑌
𝑑𝑧
= −
6𝐺
𝑝𝐾
𝑦(𝑌eq − 𝑌) (1 + 𝑋)
𝐺
𝑔𝜌
𝑝𝑑
𝑝V𝑝
.(8)
The heat balance equation of cut tobacco particles was
𝑑𝑡
𝑝
𝑑𝑧
= −
𝐶
𝑤𝑡
𝑝
𝐶
𝑠+ 𝐶
𝑤𝑋
𝑑𝑋
𝑑𝑧
+
6 (1 + 𝑋) [ℎ (𝑡
𝑔− 𝑡
𝑝) − 𝐾
𝑦(𝑌eq − 𝑌) (𝑟0 + 𝐶𝑤𝑡𝑝) ]
𝜌
𝑝V𝑝𝑑
𝑝(𝐶
𝑠+ 𝐶
𝑤𝑋)
.
(9)
The convective heat transfer coefficient ℎ between the airstream and cut tobacco particles in (9) was associated withRussell number Nu and also calculated by literature [12].
The equation of gas heat balance was as follows:
𝑑𝑡
𝑔
𝑑𝑧
= −
𝑟
0
𝐶
𝑔
𝑑𝑌
𝑑𝑧
−
6𝐺
𝑝(1 + 𝑋) [ℎ (𝑡
𝑔− 𝑡
𝑝) − 𝐾
𝑦(𝑌eq − 𝑌) (𝑟0 + 𝐶𝑤𝑡𝑝)]
𝐺
𝑔V𝑝𝑑
𝑝𝐶
𝑔
.
(10)
The differential equation of pipe diameter was as follows:
𝑑𝐷
𝑑𝑧
= 2 cot 𝜃. (11)
There, 𝜃 was the expansion angle of the adjustable pipe.According to the equations from (1) to (2) and (6) to
(11), the tapered tubular air drying cut tobacco mathematicalmodel was formed. Then this mathematical model was cal-culated using fourth-order Runge-Kutta method and Matlabsoftware program.
4. Results and Discussion
4.1. The heating Experiments of Parabolic Trough Concentrat-ing System. Experiments were studied to analyze the effect ofcavity opening width on collector temperatures as shown inFigure 3.
As can be seen from Figure 3, when the solar irradiationintensity is low, the cavity with opening width of 7 cm hasfaster heating rate because the cavity with opening widthof 4 cm has relatively larger heat loss. When the run timeincreases to about 1 month, the cavity opening width of 4 cmhas larger deformation. As for the cavity with opening widthof 7 cm, it was almost not deformed except in enclosure seamsof insulated enclosure.These changes before and after heatingabout the cavity shape are shown in Figure 4.
In order to analyze the effect of cavity opening width onenergy flux density distribution, the test results were shownin Figure 5. During test process, the direct solar radiationduring measurement was 568W/m2. As can be seen fromFigure 5, the 95% of energy is concentrated in the focal lineof concentration at the 0–7.5 cm. When the cavity absorberwith opening width of 4 cm is used, there is part of theenergy outside the cavity focused on the insulation enclosure.This may be because the insulation housing material and theabsorber material were different. The cavity absorber withopening width of 4 cm has severe deformation when it isstressed. The cavity absorber with opening width of 7 cm canabsorb most of the energy. The insulation enclosure absorbsless energy and only causes deformation of shell commissure.
When the impact of tracking accuracy is ignored, theoptical efficiency of this collector system is equal to theproduct of mirror reflectivity, acquisition factor, and thecavity absorptivity. The acquisition factors of the cavityabsorber with the opening width of 4 cm and 7 cm are
-
International Journal of Photoenergy 5
13:20 13:25 13:30 13:35 13:40 13:45 13:50 13:5540
60
80
100
120
140
160
180
200
220
240
Collector outlet temperatureCollector inlet temperatureDirect solar radiation
Time
Tem
pera
ture
(∘C)
400
500
600
700
800
900
1000
Dire
ct so
lar r
adia
tion
(W/m
2)
(a) 7 cm
Collector outlet temperatureCollector inlet temperatureDirect solar radiation
11:50 12:10 12:30 12:50 13:10 13:3080
100
120
140
160
180
200
220
240
400
500
600
700
800
900
1000
1100
1200
1300
Time
Tem
pera
ture
(∘C)
Dire
ct so
lar r
adia
tion
(W/m
2)
(b) 4 cm
Figure 3: Effect of cavity opening width on collector temperature(experimental conditions: mass of heat conduction oil = 42.65 kgand flow rate = 0.15 kg/s).
0.81 and 0.96, respectively. The collector system efficiency iscalculated using the following formula:
𝜂sys =𝑚𝐶oil (𝑡out − 𝑡in)
𝐴
𝑎𝐻
𝑏
. (12)
When the heat conduction oil temperature is 230∘C,the collector system efficiency of the cavity absorber withthe opening width of 4 cm and 7 cm is 10.9% and 25.2%,respectively. The result shows that the heating rate, acqui-sition factor, and collector system efficiency increase withincreasing the opening width of the absorber and can alsoprevent the cavity from deformation. It naturally follows thatstability of heat conduction oil output temperature can beimproved.
(a) 4 cm
(b) 7 cm
Figure 4: The changes about the cavity shape before and afterheating.
0.0 1.0 2.0 3.0 4.0 5.0 6.0 7.0
0
7
14
21
28
35
42
49
56
×103
Flux
(W/m
2)
x0 (cm)
x1 = 7 cm
x2 = 4 cm
Energy flux density distribution
Figure 5: The test results of focal line energy flux density distribu-tion.
4.2. Design and Construction of the VariableDiameter Drying Pipe
4.2.1. Approximation of Variable Diameter Drying Pipe. Thispaper uses the Fedorov method on the drying tube diameterand height of approximate calculation.
The number𝐾𝑖is calculated by the following equation:
𝐾
𝑖= 𝑑
𝑝×
3√
4𝑔 (𝜌
𝑝− 𝜌
𝑔)
3𝜇
2
𝑔
𝜌
𝑔
.
(13)
Through the Figure 6 found the relationship betweenRe𝑝and 𝐾
𝑖and through the Figure 7 found the relationship
between Re𝑝and Nu, the heat transfer coefficient of 𝐾
between air and material particles is as follows:
𝐾 =
Nu ⋅ 𝜆𝑔
𝑑
𝑝
. (14)
-
6 International Journal of Photoenergy
1 10 102 1030.1
1
10
102
103
104
Rep
Ki
2345
2345
2345
23
23
45
0.20.30.4
2 3 4 5 6789 2 3 4 5 2 3 4 5
Figure 6: The relationship between Re𝑝
and 𝐾𝑖
.
0 50 1002
5
8
Nu
Rep
Figure 7: The relationship between Re𝑝
and Nu.
Suspension velocity of particles is as follows:
𝜐
𝑎=
Re𝑝𝜇
𝑔
𝑑
𝑝
. (15)
The total surface area of the particle is as follows:
𝐴
𝑠=
6𝐺
𝑝
𝑑
𝑝𝜌
𝑝
. (16)
The mean temperature difference between the material andthe air is
Δ𝑡 =
𝑡
1+ 𝑡
2
2
− 𝑡
𝑠. (17)
The time of the drying process is
𝜏 =
𝑄
𝐾𝐴
𝑠Δ𝑡
. (18)
The length 𝑧 of the drying pipe is as follows:
𝑧 = (V𝑔− V𝑎) 𝜏. (19)
(a) 𝜃 = 50∘
(b) 𝜃 = 52∘
Figure 8:The velocity of flow field in variable diameter drying pipe.
The diameter 𝐷 of the drying pipe is
𝐷 = √
4𝐺
𝑔𝜌
𝑠
3600𝜋𝜐
𝑔
. (20)
The design parameters are substituted into the aboveformula of cut tobacco drying. The parameters are shown inTable 3.
4.2.2. The Simulation Speed Flow of Variable Diameter DryingPipe. Reducing type dry pipe consists of two different diam-eters of straight pipe and an expansion pipe connection. Inthe actual environment, the improper value of 𝜃 will lead tosediment in the drying pipe on drying process. According tothe actualmodel usingANSYS Fluent for the fluid flowmodelthe sediment will make the flow channel blockage and have agreat influence on the drying effect even an accident.
Figure 8 shows the entrance of hot air temperature is200∘C, velocity of 15m/s; the variable angle 𝜃 is 50∘ and52
∘, respectively. We get the velocity of flow field in variablediameter drying pipe through ANSYS Fluent.
Figure 8 shows that, when the 𝜃 is 50∘, there is refluxarea in variable diameter of variable diameter drying pipe.The reflux area will make sediment in the drying pipe ondrying process. When the 𝜃 is 52∘, there is no reflux area invariable diameter drying pipe. In practical design, the heightand heat loss of variable diameter drying pipe are considered;
-
International Journal of Photoenergy 7
1200mm
250mm
300mm
2000mm
150mm
82∘360mm
93.8mm
23∘ 30∘
200
mm
45∘
60∘50mm
200mm
(a) The parameters of the drying pipe (b) The photo of variable diameter drying pipe
Figure 9: The diagram of variable diameter drying pipe.
Table 3: The design parameters of cut tobacco drying.
Parameters Density(kg/m3)Specific heat(kJ/kg⋅∘C)
Initial velocity(m/s)
Initial moisturecontent (kg/kg)
Initial temperature(∘C)
The water contentafter drying (kg/kg)
Diameter(m)
Cut tobacco 800 1.467 0 0.25 20 0.115∼0.14 1.2 × 10−3
Air 0.75 1.025 15 0.025 190 — —
the variable angle 𝜃 of 60∘ can well meet the requirements ofcut tobacco drying process.
4.2.3. Production of the Variable Diameter Drying Pipe.According to the above calculation and analysis, the variablediameter drying pipe is designed and ismade of stainless steel304 with the thickness of 2mm. It is fixed through the bracketto ensure the drying tube vertically. The parameters of thedrying pipe and the photo are shown in Figure 9.
4.3. Cut Tobacco Drying Process
4.3.1. Simulation of Cut Tobacco Drying Process. The distri-bution simulation of airflow and cut tobacco parameters isshown in Table 4. The velocities shown in Figure 10 varywith the height of drying process. In the accelerating section(range from 0m to 1.2m), the airflow velocity decreasesslightly and that of cut tobacco increases quickly withincreasing the height of drying pipe. In the transition section(from 1.2m to 1.25m), the airflow velocity decreases sharplyand that of cut tobacco increases slowly with increasing the
height of drying pipe. But in the constant section (from 1.25mto 3.25m), the velocities of the airflow and the cut tobaccodecrease gradually with the height increasing.
As shown in Figure 10(b), the temperature of air flowdecreases with the increases of drying pipe height. The tem-perature of cut tobacco increases rapidly in the acceleratingsection and then increases slowly with the increases of dryingpipe height. This is the cut tobacco is put into the dryingpipe and separated quickly by high-speed airflow. Then thecontact area between the airflow and the cut tobacco isincreased, indicating that both heat and mass transfer ratesare improved. In the accelerating section, because the relativespeed of cut tobacco particles is large, both heat transfer areasfor gas phase together with solid phase and the volumetricheat transfer coefficient are improved. It naturally follows thatcut tobacco temperature increases rapidly in the acceleratingsection. But in the constant section, the relative speed ofairflow and cut tobacco particles basically keeps unchanged,so the speed of cut tobacco particles no longer increasesand the time of cut tobacco in drying pipe is prolong-ed.
-
8 International Journal of Photoenergy
Table 4: Parameters of cut tobacco drying (cut tobacco).
Intensity/kg⋅m−3 Specificheat/J⋅(kg⋅∘C)−1Initial moisture
content/% Diameter/mSolids-gas
ratioAirflow
temperature/∘CAirflow
speed/ms−1Initial temperature
/∘C255 1467 25.05 1.2×10-3 0.1 200 15 20
0.0 0.4 0.8 1.2 1.6 2.0 2.4 2.8 3.20
2
4
6
8
10
12
14
16
Cut tobacco velocityAir velocityCut tobacco moisture content
Drying pipe height z (m)
Moi
sture
cont
ent (%
d.b.
)
Velo
city
(m/s
)12
14
16
18
20
22
24
26
(a)
Cut tobacco temperatureAir temperatureCut tobacco moisture content
0.0 0.4 0.8 1.2 1.6 2.0 2.4 2.8 3.212
14
16
18
20
22
24
26
30
60
90
120
150
180
210
Tem
pera
ture
(∘C)
Drying pipe height z (m)M
oistu
re co
nten
t (%
d.b.
)
(b)
Figure 10: Distribution simulation of airflow and cut tobacco parameters (simulation experimental conditions: inlet air temperature = 200∘C,flow rate = 15m/s, and initial moisture content of cut tobacco = 25%).
A1 B1 C1 C1 C2
11.2
11.6
12.0
12.4
12.8
13.2
13.6
14.0
14.4
14.8
Moi
sture
cont
ent (%
)
Run number
MeasuredSimulated
Figure 11: Final cut tobacco moisture contents of simulation andexperiment results.
4.3.2. Validation of Mathematical Model. The mathematicalmodel of the cut tobacco particles drying process was cal-culated using fourth-order Runge-Kutta numerical solutionmethod. Experimental conditions [12] and the results of both
simulation and experiment are shown in Table 5 and Figures11, 12, and 13, respectively.
Figures 11, 12, and 13 show the final cut tobacco moisturecontents, temperatures, and the final airflow temperaturesalong with the simulation results for the same. It can be seenfrom it that the simulation results are in good agreement withexperimental data for cut tobacco drying process.The relativeerrors between simulated and experimental values are lessthan 8%, indicating that this mathematical model is accuratefor the cut tobacco airflow drying process.
4.4. Optimization Parameters of Cut Tobacco Drying Process.The effects of initial cut tobacco moisture content, inlet air-flow velocity, and temperature on cut tobacco drying processare studied by the orthogonal experiment in three factors andthree levels. Criteria for determining the optimum parameterconditions are the final cut tobacco moisture content andthermal efficiency of the dryer. The optimization results areshown in Table 6.
According to the orthogonal test result, it can be seen thatthe thermal efficiency increases with increasing the airflowtemperature and the inlet cut tobacco moisture content butdecreases with the decrease of airflow velocity. The resultshows the optimum preparation conditions are an inletairflow velocity of 15m/s, an initial cut tobacco moisturecontent of 26%, and an inlet airflow temperature of 200∘C.The thermal efficiency of the dryer and the final cut tobaccomoisture content are 66.32% and 14.15%, respectively. The
-
International Journal of Photoenergy 9
Table 5: Summaries of experimental conditions.
Run number 𝑡𝑔in V𝑔in V𝑔out 𝑡𝑝out 𝑋in 𝑋out
A1 199.5 15 163.4 53.2 24.03 12.37B1 199.8 15 156.9 58.7 26.02 14.15C1 200.4 15 157.7 59.3 25.05 13.39C1 200.4 20 152.6 61.2 25.05 11.23C2 190.3 15 148.5 52.8 25.05 14.33
Table 6: Result of orthogonal experiment.
Run number V𝑔in (m/s) 𝑋in (%) 𝑡𝑔in (
∘C) 𝜂 (%)1 15 24 180 62.462 15 25 190 63.543 15 26 200 66.324 18 24 190 59.465 18 25 200 61.346 18 26 180 60.647 20 24 200 54.568 20 25 180 52.359 20 26 190 57.49K1 64.107 58.827 58.483 —K2 60.480 59.077 60.163 —K3 54.800 61.483 60.740 —Range 9.307 2.656 2.257 —
MeasuredSimulated
A1 B1 C1 C2Run number
20
25
30
35
40
45
50
55
60
65
70
Tem
pera
ture
(∘C)
C1
Figure 12: Final cut tobacco temperatures of simulation andexperiment results.
final cut tobacco moisture content meets the requirementsof cut tobacco drying process. The result shows that thisparabolic trough concentrating solar heating will be a poten-tial heat recourse for cut tobacco drying system.
5. Conclusion
A novel parabolic trough concentrating solar heating forcut tobacco drying system was established. The result shows
100
110
120
130
140
150
160
170
180
MeasuredSimulated
A1 B1 C1 C2Run number
Tem
pera
ture
(∘C)
C1
Figure 13: Final airflow temperatures of simulation and experimentresults.
that the heating rate, acquisition factor, and collector systemefficiency increase with increasing the opening width of theabsorber and can also prevent the cavity from deformation.The simulation results using the cut tobacco drying math-ematical model were in good agreement with experimentaldata for cut tobacco drying process. The optimum prepara-tion conditions were an inlet airflow velocity of 15m/s, aninitial cut tobacco moisture content of 26%, and an inlet
-
10 International Journal of Photoenergy
airflow temperature of 200∘C. The thermal efficiency of thedryer and the final cut tobaccomoisture content were 66.32%and 14.15%, respectively.The result showed that this parabolictrough concentrating solar heating would be one of the heatrecourse candidates for cut tobacco drying system.
Nomenclature
𝐴: Drying pipe cross-sectional area, m2𝐶
𝑠: Specific heat capacity of dry material,
kJ/(kg⋅∘C)𝐶
𝑤: Specific heat capacity of water, kJ/(kg⋅∘C)
𝐶
𝑔: Specific heat capacity of moist air,kJ/(kg⋅∘C)
𝐶oil: Specific heat capacity of heat transfer oil,kJ/(kg⋅∘C)
𝐷: Drying pipe diameter, m𝑑
𝑝: Particle size, m
𝐺
𝑔: Oven dry air flow, kg/s
𝐺
𝑝: Oven dry particle size flow, kg/s
ℎ: Convective heat transfer coefficients,W/(m2 ⋅ ∘C)
𝐾
𝑦: Mass transfer coefficient, kg/(m2⋅s)
𝑁: The number of particles per unit volume,1/m3
Nu: Nusselt number𝑃sat: Saturated vapor pressure, PaPr: Prandtl number𝑟
0: 0∘C latent heat of vaporization of water,
kJ/(kg⋅∘C)Re: Reynolds number𝑡
𝑔: Air temperature, ∘C
𝑡
𝑝: Particle temperature, ∘C
V𝑔: Steam velocity, m/s
V𝑝: Particle velocity, m/s
𝑊
𝑑: Mass transfer rate, W/(m2 ⋅ ∘C⋅s)
𝑋: The moisture content of the material, %d.b.
𝑌: Absolute humidity of air, % d.b.𝑌eq: Moisture content of air saturation𝑧: Drying pipe length, m.
Greek
𝜉: Drag coefficient𝜆
𝑔: Thermal conductivity of air, J/(m⋅∘C ⋅s)
𝜌
𝑝: Wet density of the material, kg/m3
𝜌
𝑔: Wet density of the air, kg/m3
𝜌
𝑠: Oven dry density of the material, kg/m3
𝜃: Variable diameter angle𝜂sys: Collector system efficiency𝜂: Drying pipe thermal efficiency.
Subscript
In: inletOut: outlet.
Conflict of Interests
The authors declare that there is no conflict of interestsregarding the publication of this paper.
Acknowledgments
The present study was supported by the National NaturalScience Foundation, China (Grant no. U1137605), and theApplied Basic Research Programs of Science and TechnologyDepartment Foundation of Yunnan Province, China (Grantno. 2011DFA60460).
References
[1] D. A. Boyd, R. Gajewski, and R. Swift, “A cylindrical blackbodysolar energy receiver,” Solar Energy, vol. 18, no. 5, pp. 395–401,1976.
[2] O. A. Barra and L. Franceschi, “The parabolic trough plantsusing black body receivers: experimental and theoretical analy-ses,” Solar Energy, vol. 28, no. 2, pp. 163–171, 1982.
[3] C. Qiaoli, G. Xinshi, C. Shuxin et al., “Numerical analysis ofthermal characteristics of solar cavity receiver with absorber ofa bundle of pipes,” Acta Energiae Solaris Sinica, vol. 16, no. 1, pp.21–28, 1995.
[4] L. Zhang, H. Zhai, Y. Dai, and R.-Z. Wang, “Thermal analysisof cavity receiver for solar energy heat collector,” Journal ofEngineering Thermophysics, vol. 29, no. 9, pp. 1453–1457, 2008.
[5] J. Fukuchi, Y. Ohtaka, and C. Hanaoka, “A numerical fluidsimulation for a pneumatics conveying dryer of cut tobacco,”in Proceedings of the CORESTA Congress, pp. 201–212, Lisbon,Portugal, 2000.
[6] Z. Pakowski, A. Druzdzel, and J. Drwiega, “Validation of amodel of an expanding superheated steam flash dryer for cuttobacco based on processing data,” Drying Technology, vol. 22,no. 1-2, pp. 45–57, 2004.
[7] L. Huimin and M. Ming, YC/T31-1996 Tobacco and TobaccoProducts—Preparation of Test Specimens and Determinationof Moisture Content—Oven Method, China Standard Press,Beijing, China, 1996.
[8] X. Chengmu, L. Ming, J. Xu et al., “Frequency statistics analysisfor energy-flux-density distribution on focal plane of parabolictrough solar concentrators,” Acta Optica Sinica, vol. 33, no. 4,Article ID 040800, pp. 1–7, 2013.
[9] R. Legros, C. A. Millington, and R. Clift, “Drying of tobaccoparticles in a mobilized bed,” Drying Technology, vol. 12, no. 3,pp. 517–543, 1994.
[10] R. Blasco, R. Vega, and P. I. Alvarez, “Pneumatic drying withsuperheated steam: bi-dimensional model for high solid con-centration,” Drying Technology, vol. 19, no. 8, pp. 2047–2061,2001.
[11] P. Yongkang andW. Zhong,Modern Drying Technology, Chem-ical Industry Press, Beijing, China, 1998.
[12] W. E. Ranz andW.R.Marshall, “Evaporation fromdrops, part I,”Chemical Engineering Progress, vol. 48, no. 3, pp. 141–146, 1952.
-
Submit your manuscripts athttp://www.hindawi.com
Hindawi Publishing Corporationhttp://www.hindawi.com Volume 2014
Inorganic ChemistryInternational Journal of
Hindawi Publishing Corporation http://www.hindawi.com Volume 2014
International Journal ofPhotoenergy
Hindawi Publishing Corporationhttp://www.hindawi.com Volume 2014
Carbohydrate Chemistry
International Journal of
Hindawi Publishing Corporationhttp://www.hindawi.com Volume 2014
Journal of
Chemistry
Hindawi Publishing Corporationhttp://www.hindawi.com Volume 2014
Advances in
Physical Chemistry
Hindawi Publishing Corporationhttp://www.hindawi.com
Analytical Methods in Chemistry
Journal of
Volume 2014
Bioinorganic Chemistry and ApplicationsHindawi Publishing Corporationhttp://www.hindawi.com Volume 2014
SpectroscopyInternational Journal of
Hindawi Publishing Corporationhttp://www.hindawi.com Volume 2014
The Scientific World JournalHindawi Publishing Corporation http://www.hindawi.com Volume 2014
Medicinal ChemistryInternational Journal of
Hindawi Publishing Corporationhttp://www.hindawi.com Volume 2014
Chromatography Research International
Hindawi Publishing Corporationhttp://www.hindawi.com Volume 2014
Applied ChemistryJournal of
Hindawi Publishing Corporationhttp://www.hindawi.com Volume 2014
Hindawi Publishing Corporationhttp://www.hindawi.com Volume 2014
Theoretical ChemistryJournal of
Hindawi Publishing Corporationhttp://www.hindawi.com Volume 2014
Journal of
Spectroscopy
Analytical ChemistryInternational Journal of
Hindawi Publishing Corporationhttp://www.hindawi.com Volume 2014
Journal of
Hindawi Publishing Corporationhttp://www.hindawi.com Volume 2014
Quantum Chemistry
Hindawi Publishing Corporationhttp://www.hindawi.com Volume 2014
Organic Chemistry International
ElectrochemistryInternational Journal of
Hindawi Publishing Corporation http://www.hindawi.com Volume 2014
Hindawi Publishing Corporationhttp://www.hindawi.com Volume 2014
CatalystsJournal of