An Enviro-Economic Comparison of Two Different CO 2 Avoidance Processes for NGCC-Derived Electricity

Similar documents
Optimal Design Technologies for Integration of Combined Cycle Gas Turbine Power Plant with CO 2 Capture

Automated Optimisation of Multi Stage Refrigeration Systems within a Multi-Objective Optimisation Framework

Optimisation of Pressure Swing Adsorption (PSA) Process for Producing High Purity CO 2 for Sequestration Purposes

The Novel Design of an IGCC System with Zero Carbon Emissions

Heat Integration of an Oxy-Combustion Process for Coal- Fired Power Plants with CO 2 Capture by Pinch Analysis

Novel Approach for Integrated Biomass Supply Chain

Cost Effective Heat Exchangers Network of Total Site Heat Integration

Waste Management Pinch Analysis (WAMPA) with Economic Assessment

Energy-Efficient Co-production of Hydrogen and Power from Brown Coal Employing Direct Chemical Looping

Assessment of Hydrogen and Power Co-Generation Based on Biomass Direct Chemical Looping Systems

CFD Study of Biomass Cooking Stove using Autodesk Simulation CFD to Improve Energy Efficiency and Emission Characteristics

Utility-Heat Exchanger Grid Diagram: A Tool for Designing the Total Site Heat Exchanger Network

Maximizing the Efficiency of a Heat Recovery Steam Generator for Solid Oxide Fuel Cell-Based Trigeneration Systems

Pathways for deploying CCS at Australian power plants

Optimization on Scheduling for Cleaning Heat Exchangers in the Heat Exchanger Networks

Description and evaluation of flowsheet modifications and their interaction for an efficient monoethanolamine based post-combustion CO 2 capture

Targeting and Design of Evacuated-Tube Solar Collector Networks

Performance and Costs of CO 2 Capture at Gas Fired Power Plants

A Linear Mathematical Model to Determine the Minimum Utility Targets for a Batch Process

Assessment of Carbon Capture Options for Super-Critical Coal-based Power Plants

Participant Plants and Streams Selection for Interplant Heat Integration among Three Plants

Evaluating the Potential of a Process Site for Waste Heat Recovery

Utilizing Solar Energy in Order to Reduce Energy Loads of Building

Modelling and Simulation of MSF Desalination Plant: The Effect of Venting System Design for Non-Condensable Gases

Techno-Economic Analysis of Low Temperature Waste Heat Recovery and Utilization at an Integrated Steel Plant in Sweden

Comparison Between a Detailed Pinch Analysis and the Heat Load Model for Pulp and Paper Case Study for a Swedish Thermo-Mechanical Pulp and Paper Mill

Roman Hackl*, Simon Harvey VOL. 29, Introduction

Water Network Optimisation with Consideration of Network Complexity

Resilient decision making in steam network investments

Minimization of Embodied Carbon Footprint from Housing Sector of Malaysia

David J. Kukulka a, *, Rick Smith b, Wei Li c VOL. 45, Introduction

Optimal Selection of Steam Mains in Total Site Utility Systems

Evaluating the Economic Efficiency of the Technologies for Greenhouse Gas Footprint Reduction

Thermo-Economic Modelling and Process Integration of CO 2 -Mitigation Options on Oil and Gas Platforms

Targeting Minimum Heat Transfer Area for Heat Recovery on Total Sites

Design of Integrated Solar Thermal Energy System for Multi- Period Process Heat Demand

H 2 S Oxidative Decomposition for the Simultaneous Production of Sulphur and Hydrogen

PERFORMANCE EVALUATION OF NGCC AND COAL-FIRED STEAM POWER PLANTS WITH INTEGRATED CCS AND ORC SYSTEMS

Integration of Steam Power Plant with Process Utility System

Available online at GHGT-9

A Modified Energy Transfer Diagram for Heat Exchanger Network Retrofit Bridge Analysis

Investigation of Heat Exchanger Network Flexibility of Distillation Unit for Processing Different Types of Crude Oil

Data Processing of Thermal Power Plants Based on Dynamic Data Reconciliation

Exploitation of Low-Grade Heat in Site Utility Systems

DISPATCHABLE CCS Minimizing the Energy Penalty Costs

Numerical Investigation of the Air Flow in a Novel PV-Trombe Wall Based on CFD Method

Total Site Utility Systems Optimisation for Milk Powder Production

Using Experimental Measurements of the Concentrations of Carbon Dioxide for Determining the Intensity of Ventilation in the Rooms

LNG Port Interconnection with the Natural Gas Distribution Network in Finland

Heat Integration Across Plants Considering Distance Factor

Sustainable Milk Powder Production using Enhanced Process Integration and 100 % Renewable Energy

Available online at Energy Procedia 37 (2013 ) GHGT-11

Calin-Cristian Cormos* VOL. 37, Introduction

Numerical Simulation of Spray Cooling Gas with High Temperature and Velocity in Compressor Inter-Stage

Thermo-hydraulic Design of Solar Collector Networks for Industrial Applications

CHEMICAL ENGINEERING TRANSACTIONS

Simulation of CO 2 Absorption Using the System K 2 CO 3 - Piperazine

A Novel Synergistic 4-column Methanol Distillation Process

Thermal Dynamic Model and Analysis of Residential Buildings

Secondary Combustion Chamber with Inbuilt Heat Transfer Area Thermal Model for Improved Waste-to- Energy Systems Modelling

Techno-economic Assessment of the Fermentative Hydrogen Production from Sugar Beet

Optimisation economics for CO 2 Capture and Storage in Central Queensland (Australia)

Life Cycle Assessment of Electricity Production from Natural Gas Combined Cycle

Optimal High-Speed Passenger Transport Network with Minimal Integrated Construction and Operation Energy Consumption Based on Superstructure Model

On the Production of Liquid Synthetic Motor Fuels

Available online at ScienceDirect. Energy Procedia 63 (2014 ) GHGT-12

Overview of the CO2CRC Capture Program

MODERN COAL-FIRED OXYFUEL POWER PLANTS WITH CO 2 CAPTURE ENERGETIC AND ECONOMIC EVALUATION

Industrial Pasteurisation Process Improvement Possibilities

Sustainable Design and Operation of Shale Gas Supply Chain Networks with Life Cycle Economic and Environmental Optimisation

Techno-Economic Analysis of Seawater Freezing Desalination using Liquefied Natural Gas

Optimal Design of Solvent Based Post Combustion CO 2 Capture Processes in Quicklime Plants

Tuning Parameters of PID Controllers for the Operation of Heat Exchangers under Fouling Conditions

Biomass-Photovoltaic Hybrid Plant for Hydrogen Production via Steam Electrolysis

Comparison of Heat Exchanger Designs Using Vipertex 1EHT Enhanced Heat Transfer Tubes

Carbon Sources Diagram - A Tool for Carbon-Constrained Energy Sector Planning

Process Integration for Bromine Plant

Efficient Power Generation from Low Temperature Heat Source

Improving the energy efficiency of industrial processes using mathematical modelling

Carbon Emission Pinch Analysis for Sustainable Landfill

P-graph for Optimising Industrial Symbiotic Networks

Systematic Analysis of Proton Electrolyte Membrane Fuel Cell Systems Integrated with Biogas Reforming Process

Power Grid Simulation Model for Long Term Operation Planning

Optimisation of Fin Selection and Thermal Design of Plate- Fin Heat Exchangers

Practical Approach for Engineers to Optimise Industrial Ovens for Energy Saving

Comparison of Sweepgas and Vacuum Membrane Distillation as In-Situ Separation of Ethanol from Aqueous Solutions

Thermal Desalination Process Based on Self-Heat Recuperation

Exergetic Modelling of a 30-kW Gas Microturbine and Cogeneration System by Artificial Neural Networks

Environmental Index for Palm Oil Mill

Optimization of an Existing Coal-fired Power Plant with CO 2 Capture

Hybrid-Heating-Systems for Optimized Integration of Low-Temperature-Heat and Renewable Energy

A New Optimisation-based Design Methodology for Energyefficient Crude Oil Distillation Systems with Preflash Units

MULTI-OBJECTIVE OPTIMIZATION OF A FLUIDIZED BED CATALYTIC CRACKER UNIT TO MINIMIZE CO 2 EMISSIONS

DYNAMIC MODELS AND CONTROL STRATEGIES FOR ABSORPTION- BASED CARBON CAPTURE PROCESSES

CO2CRC Capture & Storage update

Integrating Renewable Energy to Power, Heat and Water Systems

Stress Study on Different Designs of Ceramic High Temperature Heat Exchangers

Available online at Energy Procedia 4 (2011) Energy Procedia 00 (2010) GHGT-10

Sensitivity Analysis of Industrial Heat Exchanger Network Design

Case Study: Optimisation of Federal Land Development Authority Jengka Supply Chain Management Using P-Graph Approach

Transcription:

289 A publication of CHEMICAL ENGINEERING TRANSACTIONS VOL. 45, 2015 Guest Editors: Petar Sabev Varbanov, Jiří Jaromír Klemeš, Sharifah Rafidah Wan Alwi, Jun Yow Yong, Xia Liu Copyright 2015, AIDIC Servizi S.r.l., ISBN 978-88-95608-36-5; ISSN 2283-9216 The Italian Association of Chemical Engineering www.aidic.it/cet DOI: 10.3303/CET1545049 An Enviro-Economic Comparison of Two Different CO 2 Avoidance Processes for NGCC-Derived Electricity Tristan Lambert a, Andrew Hoadley* a, Barry Hooper b a Monash University, Dept. Chemical Engineering Monash University Clayton Campus, Melbourne 3800, Australia b UNO Technology, 368 St Kilda Rd Melbourne 3004, Australia Andrew.Hoadley@monash.edu In order to combat the rising carbon dioxide emission level several technologies have been proposed to deal with the challenge. This paper compares the economic and emission performance of two gas turbine based technologies: integration with solar thermal energy (STE) and post combustion carbon capture and storage (CCS). Both technologies are integrated into a natural gas combined cycle (NGCC) power plant. The comparison is presented using a multi-objective optimisation (MOO) framework to provide the tradeoffs between the increased plant CO 2 equivalent (CO 2e) emission reduction and economic performance for both cases. This has been performed using a dynamic optimisation methodology enabling the calculation of the net present value of the power plants including real electricity price data along with the variable solar radiation. The economic impacts of CO 2 pricing and government co-investment are also analysed to determine the required assistance to achieve the economic breakeven point. 1. Introduction In order to address the growing problem of climate change, it will be insufficient to swap coal to gas fired power to meet the 2 C temperature rise scenario (IEA, 2012). In order to reduce the emissions of NGCC different process modifications have been suggested. The two processes studied here are carbon capture and storage (CCS) and the solar thermal energy (STE). Both of these processes are integrated into an NGCC power plant. The addition of CCS is detrimental to the plant s economic performance. Alongside the capital investment and associated capture plant operating cost, CCS reduces the electrical output from the power plant. The most commercially ready process for CCS is solvent absorption where CO 2 is absorbed from the flue gas into a solvent. The solvent is then regenerated, requiring energy input. One method of defining this energy reduction is with the energy penalty (ΔE), defined as: Pw / ocap Pwcap E % (1) P w / ocap where P w/o cap and P cap are the power plant outputs without and with capture. Sipöcz et al. (2011) shows a ΔE of 14.5 % at a capture rate of 90 %. Exhaust gas recirculation (EGR), where a portion of the flue gas is recycled to the gas turbine air inlet to increase the CO 2, which reduces the energy penalty to 12.0 % (Sipöcz et al., 2011). ElKady et al. (2009) showed the maximum EGR possible without gas turbine modifications was 35 %. The most commonly used solvent for post-combustion capture (PCC), monoethanolamine (MEA), has been shown by Supap et al. (2006) to degrade in the presence of oxygen. As a result potassium carbonate (K 2CO 3), similar to that used by Pandit et al (2014), is used for this study. Solar thermal energy integration has also been studied as a method of reducing the emissions from gas turbine power plants. The SOLGATE project (European Commission, 2005) looked at the use of solar thermal energy to preheat the combustion air prior to entering the combustion chamber. This process has been further developed by the SOLUGAS project (Krozynietz et al, 2012). Spelling et al (2012) have investigated the trade-offs between levelised cost of electricity (LCOE) and investment costs of solar Please cite this article as: Lambert T., Hoadley A., Hooper B., 2015, An enviro-economic comparison of two different co2 avoidance processes for ngcc-derived electricity, Chemical Engineering Transactions, 45, 289-294 DOI:10.3303/CET1545049

290 assisted NGCC power plants. As with carbon capture the addition of solar thermal energy requires significant capital investment, and additional operating costs. With two different processes, and the competing objectives of economic and environmental performance, multi-objective optimisation (MOO) is used for the optimisation. MOO allows for the optimisation of two or more variables at the same time allowing for the comparison of different cases on both economic and environmental benefits. The MOO results in a Pareto-optimal front show the optimised performance of the power plants over all allowable operating conditions. Similar to Kang et al. (2014), real electricity price data is used to determine the power plant performance over the whole year. 2. Methodology The power plant and capture process were modelled within the commercial process modelling software Aspen Plus V8.4. The Heat Integration and steam cycle process optimisation was calculated from the Grand Composite Curve as in work by Harkin et al. (2012). The gas turbine used as the basis for both processes is the Alstom GT26 (Alstom, 2007). The gas turbine power plant performance is given in Table 1. The focus in this paper is the comparison between the CO 2 reducing process of CCS and STE integration. In the case of CCS the solvent absorption process is modelled with steam extraction from the steam turbine used to provide the energy for solvent regeneration. The captured CO 2 is then compressed to 100 bar(a) for storage. For the solar thermal assisted case the solar tower is modelled using a heater to provide the energy to the compressed air. The amount of energy provided to the heat is given by solar field size and efficiencies. The two processes which are optimised are shown in Figure 1, and performance parameters are given in Table 2. The Aspen Plus models were used to generate surrogate models of the processes in order to improve computation time. The surrogate model was created using the SUMO software toolbox (Gorissen et al. 2010) using artificial neural networks. The SUMO generates a set of empirical formulas that can then be used to calculate the process performance and this reduces the total CPU time required for the multiple simulations required in the MOO framework. A comparison of the surrogate model performance to the Aspen Plus simulation for the NGCC CCS process is shown in Figure 2. The MOO used in the optimisation is the Genetic Algorithm based code NSGA-II developed by Deb et al. (2002), which Sharma et al. (2012) have adapted to operate in VBA through Microsoft Excel. NSGA-II has been applied to CO 2 reduction by Li Yuen Fong et al (2014) and Sharma et al. (2014). The competing values of economic performance, defined by the Net Present Value (NPV), and the environmental performance, defined by the reduction in CO 2 emissions, will be optimised using the MOO. First defined by Ho et al. (2009) the economic parameters are given as part of Table 2. In line with Pandit et al. (2014), economic performance is calculated in 2011 Australian dollars ($AU). The South Australian electricity prices are taken from the 2013-14 financial year (AEMO, 2014). During this economic period Australia had a carbon price of 24.15 $AU/tCO 2e and this has been included in the economic analysis. 2011-12 financial year data, inflated to 2013-14, has also been used to compare the effect of the use of a carbon price. Three carbon prices are analysed 0, 24.15 and 50 $AU/tCO 2 along with a government coinvestment of 50 % on the carbon emission reduction equipment. This co-investment affects the carbon capture plant and solar thermal equipment, but not the base NGCC plant. Table 1: Alstom GT26 combined cycle power plant performance (Alstom, 2007) Natural Gas Combined Cycle Gas Turbine (GT) ISO Rating (MW) 288 GT Exhaust Gas Flow (kg/s) 650 GT Exhaust Gas Temperature ( C) 614 Steam Turbine Power Output (MW) 157 Steam Turbine HP Pressure/Temperature (bar/ C) 135.3/565 Steam Turbine MP Pressure/Temperature (bar/ C) 28/565 Steam Turbine LP Pressure/Temperature (bar/ C) 4.74/287

291 (a) NGCC with Post-Combustion Carbon Capture (b) Solar thermal assisted NGCC Figure 1: Basic process diagrams for the two processes studied, (a) NGCC power plant with solvent absorption post-combustion carbon capture, and (b) solar thermal assisted NGCC power plant The dynamics of the solar irradiance and electricity prices are handled using a statistical approach. Rather than model each 8,760 h individually a frequency distribution of the solar and electricity data is used to determine the power plant performance as well as costs and revenue over the whole year. Real solar and electricity data is used where available. The calculations are performed by integrating across the price and solar ranges to determine overall plant performance. In the solar assisted NGCC power plant, for each solar irradiance level the power output of the plant is calculated, along with the expected value of the electricity produced. Table 2: Process and economic performance parameters. Solar efficiency data from Xu et al. (2011) Solar Plant Data Location Port Augusta, South Type Australia Solar Tower Maximum Temperature ( C) 950 Pressure drop (kpa) 30 Field Efficiency (%) 75 Receiver Efficiency (%) 90 Economic Data Discount Rate (%) 7 Load Factor (%) 90 Project Lifetime (y) 25 Natural Gas Price ($AU/GJ) 4 Build Time (y) 2 Carbon Prices ($AU/tCO 2e) 0, 24.15, 50 CO 2 Storage Cost ($AU/tCO 2e) 6.03 Government Co-investment (%) 0, 50

292 Figure 2: Comparison of the NGCC CCS plant surrogate model to the Aspen Plus model. Number of generations and simulation CPU time is shown The expected value, EV, is defined as: EV (2) where x is the electricity price and p(x) is the probability of that price. Real solar and electricity data is used where available, and numerical integration is used. The solar irradiance distribution for Port Augusta, South Australia is shown in Figure 3, along with the electricity price data for South Australia. The electricity price data shows three price distributions based on the carbon price. Real 2013-14 financial year data, with a carbon price of 24.15 $AU/tCO 2e is shown along with zero carbon price data from the 2011-12 data inflated to 2013-14 prices and a 50 $AU/tCO 2e hypothetical data set modelled with a log-normal distribution. Comparing the real data from 2013-14 to the inflated data from 2011-12 shows an increase in average electricity price of 29.92 $AU/MWh (61.71 $AU/MWh in 2013-14 to 31.79 $AU/MWh for the inflated data). This gives an average electricity price increase of 1.24 $AU/MWh per 1 $AU/tCO 2e. This is then used to create the hypothetical carbon price cases. For a carbon price of 50 $AU/tCO 2e this results in an estimated average electricity price of 93.74 $AU/MWh. 3. Results Results of the different scenarios are given in Figure 4. Figure 4 (a) shows the effect of changing the carbon price of 0 $AU/tCO 2, 24.15 $AU/tCO 2 and 50 $AU/tCO 2. Figure 4 (b) shows the effect of a 50 % government co-investment on the CO 2 reduction equipment, i.e. the carbon capture and solar thermal plants. Solar assisted NGCC power plants cannot achieve more than around 15 % CO 2 reduction as the amount of potential energy from the solar field that can be harnessed by the NGCC plant is reduced. The cost of CO 2 captured is shown in Figure 5. (a) (b) Figure 3: Frequency plots used with (a) solar irradiance (bin size = 1 W/m 2 ), and (b) electricity price data (bin size = 1 $AU/MWh), including three carbon price cases

293 (a) (b) Figure 4: Net Present Value (NPV) for all cases with (a) comparing the effect of carbon price, and (b) government co-investment. Base NGCC performance is shown at 0 % CO 2 reduction 4. Discussion Increasing the carbon price increases the NPV of all cases, as the emission intensity of all plants is lower than the increase in electricity price compared to the carbon price increase. The increase in average electricity price corresponds to an increase comparable to brown coal power plant emission intensities of around 1.24 tco 2/MWh. As a result the increase in the electricity price, means that the electricity produced is more valuable than reducing CO 2 emissions for NGCC power plants. The base NGCC has an emission intensity of 0.34 tco 2/MWh. A government co-investment also improves the NPV performance of both cases as the capital expenditure is reduced. The effect of the reduced capital requirement is greater at higher capture rates where the amount of additional capital has a higher cost, and hence reducing the cost has a greater impact. Above solar field areas of around 600,000 m 2 the amount of extra energy usable in the GT is reduced and the extra field size offers reduced impact as most of the energy cannot be utilised. As the energy cannot be used the cost of maintaining the equipment exceeds the savings in natural gas and carbon price, increasing the cost of the solar assisted cases further. This leads to the drop in NPV at a 15 % CO 2 reduction. (a) (b) Figure 5: Cost of CO 2 captured for each of the cases: (a) carbon price cases, and (b) government coinvestment. Cost of CO 2 captured calculated by ΔNPV/Σ Discounted CO 2 Reduced reduction level.

294 5. Conclusions Using a statistical approach, dynamic MOO was performed on two different technologies for reducing NGCC based power plant CO 2 emissions. Using real electricity price and solar irradiance data carbon capture and storage and solar thermal energy were assessed based on economic and emission reduction potential. At most CO 2 reduction level the CCS cases outperformed the solar assisted cases. However, even with co-investment, there is still a significant cost associated with reducing CO 2 emissions. Acknowledgements The authors would like to acknowledge the funding provided by the Australian Government through its CRC Program to support this CO2CRC research project. References AEMO (Australian Energy Market Operator), 2014, <www.aemo.com.au/electricity/data/price-and- Demand/Aggregated-Price-and-Demand-Data-Files/Aggregated-Price-and-Demand-2011-to-2015>, accessed 15.11.2014. Alstom, 2007, Transforming Tallawarra: GT26 based plant for booming New South Wales, Australia, Modern Power Systems, TBC 04 03E. Deb K., Pratap A., Agarwal S., Meyariven T., 2002, A fast and elitist multi-objective genetic algorithm: NSGA-II, Evolutionary Computing, IEEE Transactions, 6(2), 182-197. ElKady A., Evulet A., Brand A., Ursin T., Lynghjem A., 2009, application of exhaust gas recirculation in a DLN F-Class combustion system for post-combustion carbon capture, Journal of Engineering for Gas Turbines and Power, 131, 034505-1-6. European Commission, 2005, SOLGATE: solar hybrid gas turbine electric power system, Final Publishable Report, Directorate-General for Research, Brussels, No. KI-NA-21615\_EN-C. Gorissen D., Crombecq K., Couckuyt I., Dhaene T., Demeester P., 2010, A surrogate modelling and adaptive sampling toolbox for computer based design. Journal of Machine Learning Research, 11, 2051-2055. Harkin T., Hoadley A., Hooper, B., 2012, Optimisation of power stations with carbon capture plants - tradeoffs between costs and net power, Journal of Cleaner Production, 34, 98-109. Ho M., Allinson G., Wiley D., 2009, Factors affecting the cost of capture for Australian lignite coal fired power plants, Energy Procedia, 1, 763-770. IEA, 2012, Energy Technology Perspectives - How to secure a clean energy future, International Energy Agency (IEA), Paris, France. Kang C.A., Brandt A.R., Durlofsky L.J., 2014, Optimizing heat integration in a flexible coal-natural gas power station with CO 2 capture, International Journal of Greenhouse Gas Control, 31, 138-152. Korzynietz R., Quero M., Uhlig R., 2012, SOLUGAS - future solar hybrid technology, International SolarPACES Conference, Marrakesh, <cms.solarpaces2012.org/proceedings/paper/ 7ee7e32ece8f2f8e0984d5ebff9d77b> accessed 15.11.2014. Li Yuen Fong J.C., Anderson, C., Hooper, B., Xiao, G., Webley, P., Hoadley, A., 2014, Multi-objective optimisation of hybrid CO 2 capture processes using exergy analysis, Chemical Engineering Transactions, 29, 1501-1510. Pandit J., Hatkin T., Anderson C., Ho M., Wiley D., Hooper B., 2014, CO 2 emission reduction from natural gas power stations using a precipitating solvent absorption process, International Journal of Greenhouse Gas Control, 28, 234-247. Sharma I., Hoadley A., Mahajani S., Ganesh A., 2014, Automated optimisation of multi stage refrigeration within a multi-objective optimisation framework, Chemical Engineering Transactions, 29, 25-30. Sharma S., Rangaiah G., Cheah K., 2012, Multi-objective optimization using MS Excel with an application to design of a falling film evaporator system, Food and Bioproducts Processing, 90, 2, 123-134. Sipöcz N., Tobiesen A., Assadi M., 2011, Integrated modelling and simulation of a 400 MW power plant with CO 2 capture, Energy Procedia, 4, 1941-1948. Spelling J., Favrat D., Martin A., Augsburger G., 2012, Thermoeconomic optimization of a combined-cycle solar tower power plant, Energy, 41, 113-120. Supap T., Idem R., Tontiwachwuthikul P., Saiwan C., 2006, The roles of O 2 and SO 2 in the degradation of monoethanolamine during CO 2 absorption from industrial flue gas streams, EIC Climate Change Technology, IEEE, 1-6. Xu C., Wang Z., Li X., Sun F., 2011, Energy and exergy analysis of solar power tower plants, Applied Thermal Engineering, 31, 3904-3913.