Research Article | Open Access
Analysis of the Hybrid Power System for High-Altitude Unmanned Aircraft
The application of single solar array on high-altitude unmanned aircraft will waste energy because of its low conversion efficiency. Furthermore, since its energy utilization is limited, the surface temperature of solar array will rise to 70°C due to the waste solar energy, thus reducing the electrical performance of the solar array. In order to reuse the energy converted into heat by solar array, a hybrid power system is presented in this paper. In the hybrid power system, a new electricity-generating method is adopted to spread the photovoltaic cell on the wing surface and arrange photothermal power in the wing box section. Because the temperature on the back of photovoltaic cell is high, it can be used as the high-temperature heat source. The lower wing surface can be a low-temperature cold source. A high-altitude unmanned aircraft was used to analyze the performances of pure solar-powered aircraft and hybrid powered aircraft. The analysis result showed that the hybrid system could reduce the area of wing by 19% and that high-altitude unmanned aircraft with a 35 m or less wingspan could raise the utilization rate of solar energy per unit area after adopting the hybrid power system.
The use of electric power for flight vehicle propulsion is not new . The first one was the hydrogen-filled dirigible France in year 1884 that won a 10 km race around Villacoublay and Meudon . However, in order to realize “eternal flight,” solar energy is the best option due to its characteristic of recycling.
Since the first flight of a solar-powered aircraft (Sunrise I, designed by R. J. Boucher, Astro Flight Inc.) was succeeded on the dry lake at Camp Irwin, California, on November 4, 1974 , more and more solar-powered aircraft had been designed and validated (see Figure 1).
The most important solar-powered aircraft are introduced as follows. The first one is the Pathfinder series program developed by AeroVironment Inc. and supported by the US government funding [4–6]. During the implementation of this program, four prototypes were established: Pathfinder, Pathfinder Plus, Centurion, and Helios. The Helios aimed to demonstrate sustained flight at an altitude near 30,480 m (100,000 ft) and to fly without stop for at least 24 h. Unfortunately, it never reached the second objective as it was destroyed when it fell into the Pacific Ocean on June 26, 2003, due to the too large structure and structural failures . The second one is Zephyr series program developed by QinetiQ in the UK with the UK Ministry of Defense [7, 8]. Zephyr was a lightweight UAV (e.g., the latest Zephyr 7 was only 53 kg) designed to fly at the altitude up to 20,000 m for months at a time. However, its payload was only 2.5 kg showing the narrow application prospect. At present, a genuinely useful solar-powered aircraft is not available.
Currently, the main technology of photovoltaic cell is to generate electricity from photoelectric characteristics , but only the shorter wavelength light is utilized. Taking the silicon-based photovoltaic cells as an example, silicon solar battery only absorbs the light with the shorter wavelength and the light energy is fully converted into heat energy when the wavelength of the light is greater than 1.2 μm. The heat will not only make the temperature of battery rise to 70°C or above in certain circumstance [10, 11] but also lead to the sharp decline of cell efficiency.
Therefore, due to the low conversion efficiency and low specific power of photovoltaic cells, the current high-altitude unmanned aircraft design size is large, while the load capacity is limited [12, 13]. On the other hand, in order to realize “eternal flight,” we also need to use the fuel cell or lithium battery or both of the combinations  to cooperate with the solar cell to achieve the energy recycling.
Zephyr (shown in Figure 2) is categorized as the lightweight solar-powered UAV (unmanned aerial vehicle).
The aircraft is made of carbon-fiber material which can utilize solar energy through the solar array (from UNI-SOLAR) to drive two propellers and charge the installed lithium sulphur batteries (from Sion Power). Zephyr 7 holds the official aircraft-endurance world record for the UAV. The record has been verified by the FAI (Fédération Aéronautique Internationale) in the duration (14 days and 21 min in Yuma, 2010) and altitude (the maximum altitude is 24,560 m) [8, 15]. Based on Zephyr, a new solar-powered UAV has been constructed at Mälardalen University . In this paper, taking the Zephyr 7 solar-powered UAV as a model, we analyzed the performances of both pure photovoltaic cell and hybrid power used in the solar aircraft.
2. Traditional Design of Pure Photovoltaic Cell
The main components of the aircraft include fuselage, motor, solar array, storage battery, and straight wing with the high lift coefficient. Relevant data of the wing are provided as shown in Table 1 .
Because the wing is a straight wing manufactured with carbon-fiber composite materials, the structural weight of wing has an approximately linear correlation with wingspan. Therefore, the structural weight of wing per wingspan is about kg/m.
The Sun Power CSS125-125, which was one of the most efficient solar panels commercially available nowadays, was selected to constitute the solar panel. The solar panel with 96 photovoltaic cells has the conversion efficiency of 19.6%. The related parameters of the solar panel are described as follows: max power 320 W, temperature range between −40°C and 85°C, the power per wingspan ( W/m), and the weight of photovoltaic system per wingspan ( kg/m).
Welgard Motor C1823-20 is the most efficient power system with the minimal requirement of motor power input of 33.2 W and the thrust of 192 g. The engine power has a linear relationship with thrust. The formula is , where is the thrust of engine; is power of engine; is the scaling factor.
Important data of the weights of other components are described as follows: storage battery (including charger) ( kg), motor including propeller ( kg), and other systems ( kg). The total weight of solar-powered aircraft with only photovoltaic cells is
The specific power-to-weight ratio for the aircraft is
The lift can be calculated according to (2) :where is lift; is the area of wing; is the lift coefficient. The equation of drag is expressed as  where is the drag; is the area of wing; is the drag coefficient. In addition, the lift is equal to the aircraft weight , and the thrust is equal to the aircraft drag in straight and level flight.
3. A Photovoltaic and Thermal Module
In order to improve the solar-powered aircraft flight capability, many scholars have proposed a variety of ideas. For example, Brandt and Gilliam proposed storing energy by gravitational potential and believed that the plane could use plenty of solar energy to climb up in the day and be an unpowered glider at night to reduce battery mass . However, Brandt only gave the initial idea. Sachs and Costa studied the solar-powered aircraft flight trajectory with minimum energy storage and found that when climbing in the day and gliding at night, theoretically, solar-powered aircraft could realize the flight in day and night without any storage battery. However, the maximum climbing height in the day is more than 20 km and the minimum gliding height at night is only about 1 km. The difference between the maximum climbing height and the minimum gliding height does not meet the design requirements because of large flight altitude span . Gao et al. from National University of Defense Technology studied the high-altitude long-endurance flight based on energy extraction from wind shear and put forward the idea of dynamics simulation. However, the idea requires the infrastructure for monitoring wind field [7, 21].
In addition, solar radiation in the normal direction of the wing plane can be increased by adjusting the aircraft’s flight attitude. In order to maximize remaining energy, Klesh and Kabamba established the analysis models of flight dynamics, solar radiation intensity, and energy consumption on flight attitude and studied the relationship between the solar-powered aircraft flight trajectory of hovering at a given height and remaining energy . The analysis results showed that remaining energy obtained by optimized hovering trajectory was greater than the horizontal straight flight. Spangelo et al. expanded the flight space further to study the relationship between flight trajectory and the remaining energy in three-dimensional cylindrical space, presented the design method for flight trajectory optimization, and finally obtained the same results to Spangelo et al.’s . However, in these methods, the flying platform should meet the established mission profile to some extent, thus limiting the wide application of solar aircraft.
These studies could not enhance the actual power output capability of flying platform essentially. Therefore, the long-endurance flight capacity of solar-powered aircraft cannot be effectively utilized according to the mission requirements. Based on the structural characteristics of the aircraft, the analysis of solar-powered aircraft flight environment, and the temperature change on the surface of the photovoltaic cell, in the paper, we innovatively designed a photovoltaic-thermoelectric (photothermal cell) hybrid power system with the wing structure, which could improve the flight performance of solar-powered aircraft by increasing the transmission capacity of the energy system. In the system, the photovoltaic cells are spread on the wing surface and thermoelectric batteries are arranged in the wing box section because the temperature on the back of photovoltaic battery is high. The back of photovoltaic cells can be used as high-temperature heat source and the lower wing surface can be a low-temperature cold source, as shown in Figures 3 and 4. This design can reuse the solar energy converted into heat. Therefore, the efficiency of solar energy per area is increased. Moreover, the ambient temperature is also decreased.
4. Impact of Thermal Conductivity on Thermoelectric Battery
In the schematic diagram of thermoelectric battery structure (Figure 5), , , , , , and , respectively, represent the temperatures of heat source, cold source, heat plate, cold plate, the contact layer of hot side, and the contact layer of cold side. The unit is K. , , , and , respectively, represent the thermal conductivities of heat plate, cold plate, the contact layer of hot side, and the contact layer of cold side. The unit is W/K.
We made the following four assumptions on the model. Firstly, air convection and radiation between the thermoelectric cold side and the hot side and air convection and radiation between P type and N type semiconductor can be negligible. Secondly, heat exchange does not occur between the environment and other components (heat plate, cold plate, P type, and N type semiconductor). Thirdly, heat conduction occurs only along the length direction of the electric arm. Fourthly, the heat conduction coefficient does not vary with the arm length.
Supposing that the heat of the generator absorbed from the contact layer of hot side per unit time is and that the heat released from cold side is , when and , the following equation can be obtained:
Supposing that the heat of the flow generator device absorbed from the heat plate per unit time is and that the heat released from cold side is , the heat generated by Peltier effect can be expressed aswhere is Seebeck coefficient of thermocouple usually expressed in (V/K); is the current generated inside the thermoelectric generator.
When the generator works, according to the principle of heat conduction, conduction heat caused by the temperature difference between both sides of the semiconductor component can be expressed aswhere , the thermal conductivity of components, can be expressed aswhere , , and and , , and are the length, conductivity, and cross-sectional area of P-type and N-type legs, respectively.
Joule heat flux , generated by one pair of semiconductor legs when the current goes through the generator, is
Let ; thus the heat flow Joule heat transferring to the hot side and cold side is and , respectively. Therefore, the heat of the generator absorbed from the contact layer of hot side can be expressed as
At the same time, the heat released from cold side is
Therefore, the output power of the thermoelectric battery is
When , , , and , (12) can be converted into
The power generation efficiency of the thermoelectric battery can be expressed as . Then we get
The effects of thermal conduction on the output power and efficiency of the hybrid system designed in this paper are described as follows. The cold side has an important influence on the efficiency. Moreover, the heat dissipation efficiencies of both P-N junction and the cold side have a decisive role in heat dissipation of the system. Therefore, in the paper, we analyzed the influences of the thermal conductivity of cold side and P-N junction on the output power and efficiency.
Supposing that the back of photovoltaic cell is the heat plate, because the back temperature of photovoltaic cells can reach 70°C, the temperature of heat source is assumed to be 100°C. Moreover, the temperature of the cold plate is the same to ambient temperature. Let be 5 × 10−4 Ω. The contact layer of the hot side is made of ceramic material of the high thermal conductivity. Both and are 1 W/K. Then, the curves of the variations of output power and efficiency with the thermal conductivities of cold side () and P-N junction () can be obtained, as shown in Figures 6–9.
When the thermal conductivities of components vary from 0.01 W/K to 0.02 W/K and 0.03 W/K, the maximum output power is reduced from 0.0159 W to 0.0116 W and 0.0089 W. The reduction ratios are 37.1% and 30.3%, respectively. The efficiency decreases from 2.7% to 1.15% and 0.67%. The reduction ratios are 13.4% and 71.5%, respectively. The output power and efficiency of thermoelectric generator decrease with the increase of thermal conductivity of components and the decreasing rate is growing. The change can be interpreted as follows. When thermal conductivities of components are decreased, conduction heat going through the components is decreased, thus greatly increasing the temperature difference between both sides of the generator. In this way, both the output power and efficiency are increased. From the magnitude of the trend, we can see that the power and efficiency of generator are very sensitive to thermal conductivities of components. Therefore, in the selection of the materials for the hybrid system, the priority should be given to the material properties of the generator.
When the thermal conductivity of the contact layer of cold side varies from 0.075 W/K to 0.1 W/K, the maximum output power is increased from 0.0122 W to 0.0137 W and the increase ratio is 12.3%. When it varies from 0.1 W/K to 0.125 W/K, the maximum output power is increased from 0.0137 W to 0.0146 W and the increase ratio is 6.6%. The output power and efficiency of semiconductor components increase with the increase in the thermal conductivity of the contact layer of cold side, but the increasing rate is becoming smaller. The change can be interpreted as follows. When the temperatures of heat source and cold source are constant, the higher heat conductivity coefficient of cold side indicates the better heat conduction. In this way, the flow of the heat across the components is increased and the thermoelectric energy conversion of the battery and its electromotive force are improved. Finally, the output power and efficiency are increased. With the further increase in the thermal conductivity of the cold side, the heat conduction of components as well as the total heat flow is limited. This is why the increasing rate becomes smaller. For the solar-powered aircraft, it is impossible to obtain a high heat transfer coefficient through the complex cold side. Moreover, the complex cold side does not bring out a huge increase in the efficiency.
5. The Improved Design of the Hybrid Power System
Because the size of PN junction has an important influence on the efficiency of thermoelectric battery, it is necessary to study the proper number of thermoelectric batteries in the box and find the optimal length of PN junction coming firstly. Considering the specific characteristics of solar-powered aircraft , in this paper, we used the power-to-weight ratio as the best indicator to measure the quality of the design of photothermal cell.
Photothermal cell power generation system involves the heat collecting process, heat-electricity conversion process, and cooling process. The heat collection is realized by photovoltaic cells. Thermoelectric generator is composed of the hot deflector, the PN junction, the cold deflector, and heat dissipation through convection between the lower wing surface and air . According to Newton’s law of thermodynamics, the temperature difference between two-arm ends of thermopile leads to thermal current flows. The heat obtained by the couple arm from the heat end is the sum of Pult heat, Joule heat, and conduction heat . The formulas are expressed as follows:where is heat absorbed by the generator from heat source; is heat absorbed from cold side; is the loop current; is load resistance; is couple arm resistance; is the Seebeck coefficient; is the thermal conductivity of the device; , , , , , and are, respectively, the temperatures of heat source and cold side, thermal conductivity, length, and the cross-sectional area, and electrical resistivity of the P type and N type semiconductor; is the number of generators.
Therefore, the expression of electromotive force , electrical current , and the output power can be given as follows:
The basic expression of the power-to-weight ratio is
With above formulas, the power-to-weight ratio can be given as
P-N junction area does not affect the output power per unit area of the wing . Therefore, we studied the influence of single PN junction (P-end has the same length with N-end). The cross-sectional area is 1 cm × 1 cm. The parameters are provided as follows. P junction has the same parameters with N junction: cm2, = 4.83 × 10−4 V/K, = 0.003 Ω, and W/(m·K) . The environment temperature is considered as the temperature of the cold side. Figure 10 shows the variation of power output with temperature under the influences of PN junction arm length. Figure 11 shows the variation of power-to-weight ratio with temperature under the influence of length of P-N junction arm.
The result shows that when the arm length of P-N junction is 1.75 cm, the output power reaches the peak. However, when the ratio of power-to-weight reaches the peak, the arm length of P-N junction should be 0.5 cm. Therefore, the size of PN junction should be set as 1 cm × 1 cm × 0.5 cm in the design of this paper. The thickness of wing is about 43.5 mm. Considering the thickness of the fin and the heat-transferring ceramic, we set to be 5, indicating that 5 layers of thermoelectric batteries are arranged in the wing box. At room temperature, the rated power is set to be 7.6 W/m2; the thermoelectric power per wingspan is W/m. The weight of the thermoelectric system per wingspan is kg/m. Therefore, the gross weight of the hybrid solar system-powered aircraft is expressed as
In order to compare the performances of those two aircraft, we take the same value for and . Therefore, in the design of pure PV, while the plane is under the cruise state, we get the following equations:
In the design of the hybrid system, while the plane is under cruise state, we get the following equations:
The next equation can be obtained from the above four equations:where is 22.5 m; then can be figured out as 18.2 m according to (22). The result shows that the hybrid system can reduce the wingspan by 4.3 m (19%) under the condition of Zephyr 7. Therefore, the wing area can be reduced by 19%.
The power-to-weight ratio for the aircraft can be given as follows:
After comparing the power-to-weight ratio of hybrid power system with that of photovoltaic system under different wingspan states, we got the following equation:
Figure 12 can be obtained with (24). As shown in Figure 12, when wingspan is less than 35 m, the power-to-weight ratio of hybrid solar system-powered aircraft is higher than that of photovoltaic solar system-powered aircraft. An opposite result will be obtained when the wingspan is higher than 35 m.
Based on the photovoltaic system, the addition of the thermoelectric system can increase the output power per unit area of photovoltaic cell. When the output power meets the requirement of engine power, the area and weight of wing can be reduced. Therefore, the aircraft requires a lower lift and reduces the requirement of engine power. However, at the same time, the weight per unit area of solar system is increased. Therefore, the weight per unit area of the wing is also increased. The overall lift is decreased because the area of wing is reduced. Due to the above influences, the aircraft requires the higher engine power in return. The solar-powered aircraft with hybrid power system has the excellent performance when its wingspan is less than 35 m.
Combined with the wing configuration, we presented a hybrid energy design method with the higher efficiency. Taking Zephyr 7 as an example, we analyzed the weight and power of two aircraft designs which, respectively, used photovoltaic system and hybrid system. Comparative analysis shows that the hybrid system can effectively reduce the wing area by 19%.
Moreover, because power-to-weight ratio is an important performance indicator, we compared the two designs in this paper. Aircraft can adopt a photovoltaic-thermoelectric hybrid power system when the wingspan is less than 35 m. To some extent, the hybrid system can increase the utilization of solar energy and reduce the weight of the whole structure.
Conflict of Interests
The authors declare that they do not have any commercial or associative interests that represent a conflict of interests in connection with the work submitted.
This work was supported by the National Natural Science Foundation of China under Grant no. 51307004 and the Fundamental Research Funds for the Central Universities (YWF -15-GJSYS-025).
- M. Hoffborn, A historical survey of solar powered airplanes and evaluation of it's potential market [Bachelor thesis], Malardalen University, Västerås, Sweden, 2009.
- A. Noth, “History of solar flight,” in Proceedings of the 26th Congress of International Council of the Aeronautical Sciences (ICAS '08), Anchorage, Alaska, USA, September 2008.
- A. Noth, Design of solar powered airplanes for continuous flight [Ph.D. dissertation], Ecole Polytechnique Fédérale de Lausanne, ETH Zürich, Zürich, Switzerland, 2008.
- T. E. Noll, J. M. Brown, M. E. Perez-Davis, S. D. Ishmael, and G. C. Tiffany, Investigation of the Helios Prototype Aircraft Mishap: Volume I Mishap Report, NASA Langley Research Center, Hampton, Va, USA, 2004.
- N. J. Colella and G. S. Wenneker, Pathfinder and the Development of Solar Rechargeable Aircraft, 1994, https://str.llnl.gov/etr/pdfs/07_94.1.pdf.
- P. G. Carey, R. C. Aceves, N. J. Colella, K. A. Williams, R. A. Sinton, and G. S. Glenn, “A solar module fabrication process for hale solar electric UAVS,” in Proceedings of the 24th IEEE Photovoltaic Specialists Conference, Waikoloa, Hawaii, USA, December 1994.
- X.-Z. Gao, Z.-X. Hou, Z. Guo, R.-F. Fan, and X.-Q. Chen, “The equivalence of gravitational potential and rechargeable battery for high-altitude long-endurance solar-powered aircraft on energy storage,” Energy Conversion and Management, vol. 76, pp. 986–995, 2013.
- X. F. Zhu, Z. Guo, and Z. X. Hou, “Solar-powered airplanes: a historical perspective and future challenges,” Progress in Aerospace Sciences, vol. 71, pp. 36–53, 2014.
- A. E. Shahat, “PV module optimum operation modeling,” Journal of Power Technologies, vol. 94, no. 1, pp. 50–66, 2014.
- K. Khouzam, K. K. Chen, and Y. N. Poo, “Simulation and real-time modelling of space photovoltaic systems,” in Proceedings of the 24th IEEE Photovoltaic Specialists Conference. Part 2 (of 2), pp. 38–41, 1994.
- K. Sun, M. Zhu, and Q. Liu, “Membrane material-based rigid solar array design and thermal simulation for stratospheric airships,” Advances in Materials Science and Engineering, vol. 2014, Article ID 192707, 7 pages, 2014.
- E. Koutroulis, K. Kalaitzakis, and N. C. Voulgaris, “Development of a microcontroller-based, photovoltaic maximum power point tracking control system,” IEEE Transactions on Power Electronics, vol. 16, no. 1, pp. 46–54, 2001.
- C. L. Nickol, M. D. Guynn, L. L. Kohout, and T. A. Ozoroski, High Altitude Long Endurance Air Vehicle Analysis of Alternatives and Technology Requirements Development., American Institute of Aeronautics and Astronautics, 2007.
- M. Dudek, P. Tomczyk, P. Wygonik, M. Korkosz, P. Bogusz, and B. Lis, “Hybrid fuel cell—battery system as a main power unit for small unmanned aerial vehicles (UAV),” International Journal of Electrochemical Science, vol. 8, pp. 8442–8463, 2013.
- A. Rapinett, Zephyr. A high altitude long endurance unmanned air vehicle [M.S. thesis], University of Surrey, 2009, http://personal.ph.surrey.ac.uk/~phs1pr/mphys-dissertations/2009/Rapinett-MPhys09.pdf.
- A. Sultan, Solaris project: the design of a solar powered UAV [Bachelor thesis], School of Innovation, Design and Engineering, Mälardalen University, Västerås, Sweden, 2011.
- C. Y. Yin, P. Y. Gao, J. Guo, Z. Guo, and Z. X. Hou, “Mechanical performance analysis on wing's of Zephyr7 UAV,” Structure & Environment Engineering, vol. 39, no. 3, pp. 19–25, 2012.
- A. Noth, R. Siegwart, and W. Engel, “Design of solar powered airplanes for continuous flight,” in Aircraft and Spacecraft Systems: Design, Modeling and Control, ETHZ Lecture, 2007.
- S. A. Brandt and F. T. Gilliam, “Design analysis methodology for solar-powered aircraft,” Journal of Aircraft, vol. 32, no. 4, pp. 703–709, 1995.
- G. Sachs and O. Costa, “Optimization of dynamic soaring at ridges,” in Proceedings of the AIAA Atmospheric Flight Mechanics Conference and Exhibit, Austin, Tex, USA, 2003.
- X.-Z. Gao, Z.-X. Hou, Z. Guo, P. Wang, and J.-T. Zhang, “Research on characteristics of gravitational gliding for high-altitude solar-powered unmanned aerial vehicles,” Proceedings of the Institution of Mechanical Engineers, Part G: Journal of Aerospace Engineering, vol. 227, no. 12, pp. 1911–1923, 2013.
- A. T. Klesh and P. T. Kabamba, “Energy-optimal path planning for solar-powered aircraft in level flight,” in Proceedings of the AIAA Guidance, Navigation and Control Conference and Exhibit, Hilton Head, SC, USA, August 2007.
- S. C. Spangelo, E. G. Gilberty, A. T. Klesh, P. T. Kabamba, and A. R. Girard, “Periodic energy-optimal path planning for solar-powered aircraft,” in Proceedings of the AIAA Guidance, Navigation, and Control Conference, Chicago, Ill, USA, August 2009.
- F. Sánchez Lasheras, J. A. Vilán Vilán, P. J. García Nieto, and J. J. del Coz Díaz, “The use of design of experiments to improve a neural network model in order to predict the thickness of the chromium layer in a hard chromium plating process,” Mathematical and Computer Modelling, vol. 52, no. 7-8, pp. 1169–1176, 2010.
- J. J. del Coz Diaz, P. J. Garcia-Nieto, F. P. Alvarez-Rabanall, M. Alonso-Martínez, J. Dominguez-Hernandez, and J. M. Perez-Bella, “The use of response surface methodology to improve the thermal transmittance of lightweight concrete hollow bricks by FEM,” Construction and Building Materials, vol. 52, pp. 331–344, 2014.
- Y. Vorobiev, J. González-Hernández, P. Vorobiev, and L. Bulat, “Thermal-photovoltaic solar hybrid system for efficient solar energy conversion,” Solar Energy, vol. 80, no. 2, pp. 170–176, 2006.
- T. Zhang, Z. M. Zhou, X. Z. Huang, G. W. Liang, and B. Xu, “Thermal electric coupling analysis of matched load of semiconductor thermoelectric generator,” Chinese Journal of Power Sources, vol. 10, pp. 12–13, 2010.
- G. Nakamoto, T. Souma, M. Yamaba, and M. Kurisu, “Thermoelectric properties of ()4Sb3 below room temperature,” Journal of Alloys and Compounds, vol. 377, no. 1-2, pp. 59–65, 2004.
- J. C. Chen, B. H. Lin, H. J. Wang, and G. Lin, “Optimal design of a multi-couple thermoelectric generator,” Semiconductor Science and Technology, vol. 15, no. 2, pp. 184–188, 2000.
- M. Kubo, M. Shinoda, T. Furuhata, and K. Kitagawa, “Optimization of the incision size and cold-end temperature of a thermoelectric device,” Energy, vol. 30, no. 11-12, pp. 2156–2170, 2005.
Copyright © 2015 Kangwen Sun 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.