Thermal Comfort Zone for Thai People

Similar documents
Designing Air-Distribution Systems To Maximize Comfort

A Field Study of the Thermal Environment in Residential Buildings in Harbin

Thermal Comfort Evaluation of HDB flats

Optimizing Indoor Environments for Occupant Satisfaction. Presented by: Kelli Goldstone April 2016

BUILDING DESIGN FOR HOT AND HUMID CLIMATES IMPLICATIONS ON THERMAL COMFORT AND ENERGY EFFICIENCY. Dr Mirek Piechowski 1, Adrian Rowe 1

Thermal Comfort in Architecture

Adaptive Thermal Comfort in Warm Dry Climate: Economical dwellings in Mexico

IMPLICATIONS OF CLIMATE CHANGE AND OCCUPANT BEHAVIOUR ON FUTURE ENERGY DEMAND IN A ZERO CARBON HOUSE

Thermal Comfort Survey in Office Buildings in Bandung, Indonesia

Field investigation on indoor thermal environment of a high-rise condominium in hot-humid climate of Bangkok, Thailand

Available online at ScienceDirect. Procedia Engineering 145 (2016 ) 18 25

ScienceDirect. Influence of the balcony glazing construction on thermal comfort of apartments in retrofitted large panel buildings

COMPARISON OF THE STANDARDIZED REQUIREMENTS FOR INDOOR CLIMATE IN OFFICE BUILDINGS

Retrofitting Adaptive Comfort Strategies into Conventionally Air Conditioned Commercial Buildings

Lesson 29 Inside And Outside Design Conditions. Version 1 ME, IIT Kharagpur 1

Micheal Okpara University of Agriculture Umudike. Abstract

Thermal Comfort & Green Buildings. Dr. Ommid Saberi. "Green Retrofit Conference. WSP Environment & Energy (Middle East)

Practical Guide Comfort level measurement in the workplace :32. PMV PPD measurement protocol Default Point Graphic Date/Time

Experimental investigation of zero energy office under natural and forced ventilation

A Thermal Comfort Investigation of a Facility Department of a Hospital in Hot-Humid Climate: Correlation between Objective and Subjective Measurements

Assessing the ability of PMV model in predicting thermal sensation in naturally ventilated buildings in UK

ENERGY AND THERMAL COMFORT EVALUATION FOR DIFFERENT PASSIVE SOLUTIONS IN A KINDERGARTEN IN SUMMER CONDITIONS.

125BEPM - Building energy performance modeling. Lecture handouts

Work environment. Microclimate - definition. Main microclimate elements. Core temperature. Microclimate

Integrated BIPV performance assessment for tropical regions: a case study for Bangalore

LOCAL VENTILATION SYSTEMS: SOME INVESTIGATIONS ABOUT COMFORT LEVELS AND ENERGY DEMANDS

Computational and field test analysis of thermal comfort performance of user-controlled thermal chair in an open plan office

POST OCCUPANCY DESIGN INERVENTION TO IMPROVE COMFORT AND ENERGY PERFORMANCE IN A DESERT HOUSE

Assessment of Indoor Climate: Learning from Buildings

THERMAL MASS IMPACT ON ENERGY PERFORMANCE OF A LOW, MEDIUM, AND HEAVY MASS BUILDING IN BELGRADE

THERMAL MASS IMPACT ON ENERGY PERFORMANCE OF A LOW, MEDIUM AND HEAVY MASS BUILDING IN BELGRADE

Bio Climatic Analysis of Bhubaneswar-An Investigation to Arrive Human Comfort through Natural Ventilation

Optimization of air tightness and night ventilation for Passive houses in Italian climates under Fanger and Adaptive comfort models

Energy and Thermal Comfort Management in a Kindergarten School Building in the South of Portugal in Winter Conditions

NUMERICAL SIMULATION OF CAR COCKPIT HEATING DURING WINTER 1

Sensing Occupant Comfort using Wearable Technologies

SICK BUILDING SYNDROME SYMPTOMS AND PERFORMANCE IN A FIELD LABORATORY STUDY AT DIFFERENT LEVELS OF TEMPERATURE AND HUMIDITY

Available online at ScienceDirect. Energy Procedia 78 (2015 )

Nearly Zero Energy Building in Lecco Modern technological building compared to an ancient Villa

Performance evaluation of an actual building airconditioning

Thermal comfort study of residents in Mlaten housing Semarang, Indonesia

The Response of Human Thermal Perception and Skin Temperature to Step-Changed Activity Level

Busting the Myth: Passive Houses Only Work in Cold Climates. (The Comfort Proposition)

Using passive solutions to improve thermal summer comfort in timber framed houses in South-west France

Natural Ventilation Approached on Circular Courtyard Building in Hot Humid Climate

THERMAL PERFORMANCE USING RADIANT BARRIER AND PHASE CHANGE MATERIAL IN ROOFS

Thermal effect of temperature gradient in a field environment chamber served by displacement ventilation system in the tropics

COMFORT-PRODUCTIVITY Building costs COMFORT-PERFORMANCE. Radiant Heating and Cooling Systems for Better Comfort and Energy Efficiency INDOOR - OUTDOOR

Modelling Analysis of Thermal Performance of Internal Shading Devices for a Commercial Atrium Building in Tropical Climates

Evaluation of Energy Savings of the New Chinese Commercial Building Energy Standard

Methods for Improving Energy Performance of Single-family Buildings in Poland s Climatic Conditions

Outdoor Thermal Comfort in the Hot Arid Climate The effect of socio-economic background and cultural differences

Personalized conditioning with human-in-the-loop control approach

Proceedings: Building Simulation 2007

ANSI/ASHRAE STANDARD , METHODS OF TESTING CHILLED BEAMS

Evaluating Thermal Comfort in a Naturally Conditioned Office in a Temperate Climate Zone

OPERATION AND CONTROL OF THERMALLY ACTIVATED SLAB HEATING AND COOLING SYSTEMS

Thermal Environment of Chinese Residential Buildings with Different Thermal Insulation Materials in Hot-Summer-and-Cold-Winter Zone of China

Adaptive buildings' facades for thermal comfort in hot-humid climates

681: Environmentally Responsive Architecture; Passive Design for School in Southern India

UC Berkeley Indoor Environmental Quality (IEQ)

Thermal Delight in Architecture

Thermal Delight in Architecture

Glazing and Winter Comfort Part 1: An Accessible Web Tool for Early Design Decision-Making

Cool spots in hot climates: a means to achieve pedestrian comfort in hot climates

DISPLACEMENT VENTILATION

A field study of the performance of the Dutch Adaptive Temperature Limits guideline

Articles. Figure 2. Typical set-up of personalized ventilation system.

Energy Consumption Measurement of Energy-saving Buildings and Analysis of Simulation Example

The use of Simulation Software to Evaluate the Thermal Performance of Buildings in Brazil

EVALUATION OF BUILDING ENERGY CONSUMPTION BASED ON FUZZY LOGIC AND NEURAL NETWORKS APPLICATIONS

Cooling Ceiling Panel

teamf [F]oenix Public Library, [F]oenix, AZ agents of change workshop january 9, 2004

Thermal comfort conditions in outdoor spaces

Thermal performance of a closed wet cooling tower for chilled ceilings: measurement and CFD simulation

IMPACT OF COLUMNS AND BEAMS ON THE THERMAL RESISTANCE OF THE BUILDING ENVELOPE

HEATING LOAD PREDICTIONS USING THE STATIC NEURAL NETWORKS METHOD. S.Sholahudin 1, Hwataik Han 2*

Whole Building Energy Performance Evaluation through. Similar Control Strategies

Energy Performance And Long-Term Evaluation Of Internal Thermal Comfort Of An Office Building With Different Kinds Of Glazing Systems And Window Sizes

Bio-climatic Chart for Different Climatic Zones of Northeast India

Assessments of indoor environmental quality on occupant satisfaction and physical parameters in office buildings

ARCHIVES OF CIVIL ENGINEERING, LX, 1, 2014 NUMERICAL INVESTIGATIONS OF EFFECT OF INDOOR AIR QUALITY ON THERMAL COMFORT IN RESIDENTIAL BUILDINGS

The Effects of Set-Points and Dead-Bands of the HVAC System on the Energy Consumption and Occupant Thermal Comfort

DESIGN OF A GREEN DEMO BUILDING IN A HOT AND HUMID CITY IN CHINA. Ana Bacall Sebesta Blomberg & Associates Woburn, MA

Performance-Risk Methodology for the Design of High-Performance Affordable Homes

Numerical Modelling of Air Distribution in the Natatorium Supported by the Experiment

Optimisation of an HVAC system for energy saving and thermal comfort in a university classroom

CLIMATE ENERGY INDEX AND BUILDING ENERGY INDEX: NEW INDICES TO ASSESS AND BENCHMARK BUILDING ENERGY PERFORMANCE

A field survey in Calcutta. Architectural issues, thermal comfort and adaptive mechanisms in hot humid climates.

TABS, THERMAL ACTIVATED BUILDING SYSTEM, TO IMPROVE COMFORT IN SCHOOLS

Canada Published online: 28 Feb 2011.

Introduction to basics of energy efficient building design

Available online at ScienceDirect. Energy Procedia 91 (2016 )

INDOOR CLIMATE IN AIR- SUPPORTED STRUCTURE

! #! %%%%&%%% &( &% )) + # +,% #.!. # 0! 1 2! ,%,% 0 # 0! 1 2! ,4&,) 5,%

Simulation of Consumption in District Heating Systems

Theory versus practice of energy and comfort in 4 low energy houses in Belgium

Energy consumption and thermal comfort in a passive house built in Romania

Thermal comfort analysis of public transport passengers in Catania

AHU Models Using Whole Building Cooling and Heating Energy Consumption Data

Transcription:

Engineering, 013, 5, 55-59 http://dx.doi.org/10.436/eng.013.5506 Published Online May 013 (http://www.scirp.org/journal/eng) Thermal Comfort Zone for Thai People Juntakan Taweekun *, Ar-U-Wat Tantiwichien Department of Mechanical Engineering Faculty of Engineering, Prince of Songkla University, Hat Yai, Songkhla, Thailand Email: * juntakan@me.psu.ac.th, aruwat@yahoo.com Received February 1, 013; revised March 3, 013; accepted March 30, 013 Copyright 013 Juntakan Taweekun, Ar-U-Wat Tantiwichien. 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. ABSTRACT Radiant cooling system for thermal comfort of Thai people was developed. Questionnaires for subjective experiments were examined for development of radiant cooling system for thermal comfort. Thermal sensation, humid sensation and air movement sensation were used for thermal comfort assessment. The Predicted Mean Vote (PMV) value was used for thermal comfort evaluation and it was observed that the PMV values are in the comfortable ranges during 18:00 to 10:00 for air velocity 0. m/s and 0.4 m/s. Comfortable periods are extended for higher air velocity. Neutral temperatures are in the ranges of 6.44 C - 33.60 C. The percentage number of dissatisfied (PPD) value for each air velocity was also investigate and the value of PPD is five for zero value of PMV. Simulation and experimental results were evaluated by using three indicators, namely, the coefficient of variation (CV), the mean bias error (MBE) and Chi-Square. The values of CV, MBE and for mean radiant temperature are 14.5%, 1.54% and 3.95 C, respectively. Zones of acceptable thermal environment were developed for Thai People. Results revealed that the comfortable zones of Thai people are in the ranges of relative humidity 50% - 70% and effective temperatures (ET) 4 C - 7 C. Keywords: Thermal Comfort Zone; Radiant Cooling System; Neutral Temperature 1. Introduction Electricity consumption in Thailand is approximately 44%, 4% and 3% of the total electricity use for industrial, commercial and residential sectors, respectively [1]. For residential buildings in urban and suburban areas, air-conditioning is also widely used approximately 70% of the total energy use. Webb (195) pioneered thermal studies in Singapore for humid tropical regions []. Busch (1990) studied on field experiment in naturally ventilated office buildings in Bangkok, Thailand indicated a neutrality of 8.5 C ET, while 4.5 C ET was observed in air-conditioned buildings [3]. Neutral temperature for office workers was 6.7 C in Jakarta, Indonesia [4]. The use of radiant cooling system was suggested to offer quiet comfort and energy efficiency superior to conventional air-conditioning system [5]. Zones of acceptable operative temperature and humidity for people with light activity were developed by ASHRAE. The lower limit of ET for winter and the upper limit of ET for summer are 0 C and 6 C [6]. In the previous studies, * Corresponding author. effects of humidity and air velocity on thermal comfort were seldom discussed. This study concentrates on the use of radiant cooling system considering effects of temperature, humidity, air velocity, metabolic rate, clothing level by subjective experiments. Subjective experiments will be carried out to evaluate thermal comfort under the radiant cooling system of Thai people. The values of PMV and PPD are considered for thermal comfort evaluation. Measured data and data obtained from questionnaire are used for development of radiant cooling system for thermal comfort. The values of CV, MBE and are used to investigate experimental and simulation results.. Methodology.1. Design of Radiant Cooling System In this stage, the radiant cooling panel, control system and relevant equipment were designed. The designed areas of the radiant cooling panels are 3 m which calculated from heat gain into the experimental room divided by the panel cooling capacity. Copyright 013 SciRes.

56 J. TAWEEKUN, A. TANTIWICHIEN.. Experimental Set Up and Data Collection The experiment was set up at the low energy house in Prince of Songkla University, Thailand. The radiant cooling panels were constructed from copper tube bonded to aluminium sheet. The cooling tower produces cool water with temperature 5 C and supplies to the wall and ceiling radiant cooling panels with closed system. Temperature sensors installed at surfaces of cooling panels are controlled by PID controller for opening solenoid valve at 0% - 100%. Flow diagram of radiant cooling system is illustrated in Figure 1. Thermal comfort field investigation was referred to Class I as specified in ASHRAE and IS0 standards. Three heights of measurement above floor level are 0.1, 0.6 and 1. m. Experimental conditions for humidity and air movement were controlled by a humidifier and air speed level. Temperature and humidity sensors were used to measure the temperature and relative humidity of interior and exterior air. Surface temperatures of opaque wall and radiant cooling panels were measured by thermocouple type K. Mean radiant temperature was also recorded by globe thermometer. All of these data were continuous recorded to data logger. The interior air speed, electrical consumption and metabolic rate were measured by anemometer, continuous power meter and heart rate meter, respectively..3. Questionnaire Procedure Research objectives, outcomes of the research and how to answer the questionnaire are clarified to the subjects. The first part asks about personal information age, gender, weight, height, education level, type of clothes that the subjects wore on the day of investigation. In the second part, questions were asked about thermal sensation, humid sensation and air movement sensation in response to the room conditions. The subjects are required to make only one choice from the scales as shown in Table 1 for each question. The subjects voted at relative humidity from 50% to 90% and air velocity passing through the body in the ranges of 0. - 1.0 m/s..4. Development of Radiant Cooling System for Thermal Comfort of Thai People The measured data and data obtained from questionnaire are used for development of radiant cooling system for thermal comfort of Thai people. Values of PMV and PPD are investigated for thermal comfort evaluation as given in Equations (1) and (). Value of PMV is used to describe the vote and imbalance condition. Ranges of PMV for thermal comfort conditions are 0.5 to 0.5 (Comfortable); 0.5 to 1.0 (Warm); 1.0 to 0.5 (Cool); over 1.0 (Unacceptably Warm) and under 1.0 (Unacceptably cool) [7]. PMV 0.303exp 0.036M 0.08 L. (1) PPD 100 95exp 0.03353PMV 0.179PMV 4. () where L is imbalance between (M-WR) and rate of heat dissipation Wm, M is metabolic rate Wm and WR is work rate Wm which obtained from equation 5.5-15(M - 0.8). 3. Results and Discussions 3.1. Predicted Mean Vote and Percentage Number of Dissatisfied Evaluation Thermal comfort can be obtained in the wider ranges at the higher air velocity. It was observed that the PMV values are in the range of 1-1 which are in the comfortable ranges during 18:00 to 10:00 at air velocity 0. m/s and 0.4 m/s as illustrated in Figure. Comfortable Figure 1. Flow diagram of radiant cooling system. Copyright 013 SciRes.

J. TAWEEKUN, A. TANTIWICHIEN 57 periods are extended in the ranges of 15:00 to 1:00 at higher air velocity. The relationships between PMV and PPD values for each air velocity are also examined. The value of PPD is five for zero value of PMV. The number reaches 10 when PMV deviates from zero by 0.5 and 5 when PMV deviates by 1.0. 3.. Thermal Comfort Assessment All subjects are in seat and quiet (met 1.0) with short shirt and long sleeve. The subjects voted in the range of relative humidity 50% - 90% with air velocity 0. - 1.0 m/s. 3..1. Humid Sensation Analysis Results implied that the subjects feel mostly just right and slightly humid at air velocity 0., 0.4 and 0.6 m/s. Moreover, the subjects prefer mainly just right and also feel slightly dry at air velocity 0.8 and 1.0 m/s. The humid sensations are mostly voted just right in the ranges of relative humidity 50% to 60% which is in the ranges refer to ASHRAE Standard. The humid acceptability for Thai subjects decreases at the higher relative humidity and are below 80% for all cases of relative humidity above 70% as illustrated in Figure 3. Table 1. Thermal, humid and air movement sensation scales. Value Thermal Sensation +3 hot Humid Sensation Air Movement Sensation + warm humid too breezy +1 slightly warm slightly humid breezy 0 neutral just right just right 1 slightly cool slightly dry too still cool dry 3 cold Table. Thermal sensation voted (TSV) equations and neutral temperature. Air Velocity (m/s) TSV Equation R Neutral Temperature ( C) 0. TSV 0.193T i 5.10 0.8966 6.44 0.9 0.4 TSV 0.17T i 4.573 0.919 6.90 0.9 0.6 TSV 0.06T i 5.919 0.7969 8.73 0.9 0.8 TSV 0.75T i 8.165 0.918 9.69 0.9 1.0 TSV 0.19T i 6.59 0.81 3.60 0.9 3... Air Movement Sensation Analysis The subjects feel mostly just right at air velocity 0. m/s and 0.4 m/s. The votes for breezy and too breezy are obtained at the higher air velocity 0.8 m/s and 1 m/s. The highest air movement acceptability can be obtained in the ranges of 5 C and 9 C for air velocity 0. m/s and in the ranges of 6 C and 30 C for air velocity 0.4 m/s. Figure. Comparisons of PMV values. Figure 3. Humid acceptability for air velocity 0. - 1.0 m/s. 3..3. Thermal Sensation Analysis The relationships between thermal sensation vote and interior air temperature T i in the ranges of relative humidity 50% to 90% at air velocity 0.4 and 0.6 m/s are given in Figures 4 and 5, respectively. The corresponding neutral temperature is the temperature at which the regression line crosses x-axis. Results showed that neutral temperatures are in the ranges of 6.44 C - 3.60 C as given in Table. Subjective experimental results for air velocity 0.4 m/s and various relative humidity are shown in Table 3. The relations between the percentage of the subjects voting 0 and the interior temperature are shown in Figure 6 for air velocity 0.4 m/s and Figure 7 for air velocity 0.6 m/s. Temperature and air velocity have significant effect to thermal sensation vote. The higher percent of vote 0 were obtained in case of air velocity 0. m/s and 0.4 m/s. The comfort temperature can be extended to 3.60 C in case of air velocity 1.0 m/s. Never- Copyright 013 SciRes.

58 J. TAWEEKUN, A. TANTIWICHIEN Temperature ( C) Table 3. Distribution of thermal sensation votes at air velocity 0.4 m/s. RH Number of Thermal Sensation Vote (%) 3 1 0 1 3 Frequency Number Percent of 0 Vote Mean Vote 4 + 0.9 80-90 4 8 0 3 6.50 0.50 5 + 0.9 50-60 6 8 78.57 0.1 70-80 6 8 78.57 0.1 6 + 0.9 50-60 6 30 86.67 0.0 60-70 6 8 78.57 0.1 70-80 100 0.00 80-90 45 45 100 0.00 7 + 0.9 50-60 33 6 39 84.6 0.15 60-70 13 5 38 65.79 0.34 70-80 13 78 6 97 80.41 0.07 80-90 6 3 6 6 50 64.00 0.4 8 + 0.9 60-70 58 58 100 0.00 9 + 0.9 60-70 39 6 45 86.67 0.13 70-80 5 1 8 89.9 0.14 30 + 0.9 60-70 5 8 13 6 19.3 1.31 70-80 3 6 38 84.1 0.16 31 + 0.9 60-70 9 15 3 7 33.33 0.78 3 + 0.9 70-80 8 15 4 7 9.63 0.85 Figure 4. Relationships between mean thermal sensation vote and interior temperature at air velocity 0.4 m/s. Figure 6. Thermal acceptability for air velocity 0.4 m/s. Figure 5. Relationships between mean thermal sensation vote and interior temperature at air velocity 0.6 m/s. Figure 7. Thermal acceptability for air velocity 0.6 m/s. Copyright 013 SciRes.

J. TAWEEKUN, A. TANTIWICHIEN 59 theless, thermal acceptability for Thai subjects at air velocity 1.0 m/s for percent of vote 0 are only 37% which supports recommendation by ASHRAE as air velocity limited to 0.8 m/s to avoid a sensation of draft, which occurs when there is non-uniform heat transfer between different parts of the body. 3.3. Statistic Evaluation This study uses three indicators to evaluate the experi- mental and simulation results, namely CV, MBE and. EnergyPlus Program was used in this study. These indicators have been extensively used for estimation of the accuracy of neural networks and other black-box type models predicting electrical consumption [8-10]. The coefficient of variation is the root mean square error divided by the mean value of the actual output. The resultant values of CV, MBE and for mean radiant temperature are 14.5%, 1.54% and 3.95 C, respectively. 3.4. Thermal Comfort Zone As mentioned earlier, the lower limit and upper limit of effective temperature is 4 C, and 6 C, respectively and the upper limit of relative humidity is 60% developed by ASHRAE. In this study, experimental results with the values of clothing level, type of activity and air speed are used to develop zone of acceptable thermal environment of Thai People. Experimental results revealed that the comfortable zones of Thai people are in the ranges of relative humidity 50% - 70% and effective temperatures 4 C - 7 C for air velocity 0. m/s. 4. Conclusion Experimental results implied that radiant cooling system can be applied for buildings under tropical climate and cooling tower can provide cool water for radiant cooling panels. At the higher air velocity, thermal comfort range can be obtained at the higher temperature. Thermal comfort zone for Thai people has wider ranges than ASH- RAE comfort zone. This is due to high relative humidity throughout the year in Thailand and the people can adapt in the wider ranges. Policy and Planning Office (EPPO), Ministry of Energy, Thailand and Prince of Songkla University for providing the research scholarships. REFERENCES [1] http://www.eppo.go.th [] C. G. Webb, On Some Observation of Indoor Climate in Malaya, Journal of the Institution of Heating and Ventilating Engineers, Vol. 0, No. 04, 1959, pp. 189-195. [3] J. F. Busch, Thermal Responses to the Thai Office Environment, ASHRAE Transactions, Vol. 96, No. 1, 1990, pp. 859-87. [4] T. H. Karyono, Higher PMV Causes Higher Energy Consumption in Air-Conditioned Buildings: A Case Study in Jakarta, Indonesia, In: Standards for Thermal Comfort: Indoor Air Temperature Standards for the 1st Century, Chapman & Hall, New York, 1995, pp. 19-6. [5] T. Imarori, T. Omori and K. Bogaki, Thermal Comfort and Energy Consumption of the Radiant Ceiling Panel System: Comparison with All-Air System, Energy and Buildings, Vol. 30, No., 1999, pp. 167-175. doi:10.1016/s0378-7788(98)00084-x [6] American Society of Heating, Refrigerating and Air Conditioning Engineers (ASHRAE), Thermal Environmental Conditions for Human Occupancy, ANSI/ASHRAE Standard 55-199, ASHRAE, Atlanta, 199. [7] P. O. Fanger, Thermal Comfort Analysis and Applications in Environmental Engineering, Mcgraw-Hill, New York, 1970. [8] W. J. Stevenson, Using Artificial Neural Nets to Predict Building Energy Parameters, ASHRAE Transactions, Vol. 100, No., 1994, pp. 1081-1087. [9] J. F. Kreider, D. Claridge, P. S. Curtiss, R. Dodier, J. Haberl and M. Krarti, Recurrent Neural Networks for Building Energy Use Prediction and System Identification: A Progress Report, ASME Journal of Solar Energy Engineering, Vol. 117, No. 3, 1995, pp. 161-166. doi:10.1115/1.847757 [10] A. Dhar, T. A. Reddy and D. E. Claridge, Modeling Hourly Energy Use in Commercial Buildings with Fourier Series Functional Forms, ASME Journal of Solar Energy Engineering, Vol. 10, No. 3, 1998, pp. 17-3. doi:10.1115/1.88807 5. Acknowledgements This research was financially supported by the Energy Copyright 013 SciRes.