Study on Performance of Anisotropic Materials of Thermal Conductivity

Similar documents
Double-Skin Facade in Low-Latitude: Study on the Absorptance, Reflectance, and Transmittance of Direct Solar Radiation

Experimental Performance and Numerical Simulation of Double Glass Wall Thana Ananacha

2014 International Conference on Computer Science and Electronic Technology (ICCSET 2014)

Δ q = ( ψ L) HDH (1) here, Δq is the additional heat transfer caused by the thermal bridge, Ψ and L are the linear thermal transmittance and length of

Simulation on the Temperature Drop Rule of Hot Oil Pipeline

STUDY OF FLAT PLATE COLLECTOR SOLAR WATER HEATER

Available online at ScienceDirect. Procedia Engineering 121 (2015 ) 19 26

Hourly Solar Radiation Analysis of Buildings within 48 Facings in FuZhou, China

The energy performance of an airflow window

STUDY OF A PASSIVE SOLAR WINTER HEATING SYSTEM BASED ON TROMBE WALL

Study on Failure Plane of Bulk Solid in Squat Silos

Evaluation of Overall Stability for Garage Foundation Excavation of High Buildings

AR/IA 241 LN 231 Lecture 4: Fundamental of Energy

Orientation of Building

TREES Training for Renovated Energy Efficient Social housing

Aalborg Universitet. Publication date: Document Version Publisher's PDF, also known as Version of record

ADVANCED ROOF COATINGS: MATERIALS AND THEIR APPLICATIONS

Preliminary Experimental Study on Heat Transfer Characteristics of Wall with Automatic Adjustment of Heat Transfer Coefficient

STUDY ON THE THERMAL PERFORMANCE AND AIR DISTRIBUTION OF A DISPLACEMENT VENTILATION SYSTEM FOR LARGE SPACE APPLICATION

Study on the Pressure and Temperature Distribution of Solid-Plug Conveying Element on Centrifugal Extruder

Cooling of Electronics Lecture 2

P = ρ{ g a } + µ 2 V II. FLUID STATICS

Available online at ScienceDirect. Procedia Engineering 121 (2015 )

Test Study on Strength and Permeability Properties of Lime-Fly Ash Loess under Freeze-Thaw Cycles

THERMAL COMFORT IN HIGHLY GLAZED BUILDINGS DETERMINED FOR WEATHER YEARS ON ACCOUNT OF SOLAR RADIATION. Dominika Knera 1 and Dariusz Heim 1

CAE 463/524 Building Enclosure Design Fall 2012

ISO INTERNATIONAL STANDARD

2-D CFD analysis of passenger compartment for thermal comfort and ventilation

Comparison of Numerical Method for Forward. and Backward Time Centered Space for. Long - Term Simulation of Shoreline Evolution

Numerical Simulation for Flow and Heat Transfer Characteristics of L- Type Chaotic Channel

Here Comes the Sun. Rachel Wagner. Duluth Energy Design Conference 26 February

CAE 331/513 Building Science Fall 2016

Natural Convection in a Trapezoidal Enclosure Heated from the Side with a Baffle Mounted on Its Upper Inclined Surface

ARCH 348 BUILDING AND ENVIRONMENTAL SYSTEMS

Engineering Thermodynamics. Chapter 3. Energy Transport by Heat, Work and Mass

Thermal Comfort; Operative Temperature in the Sun

Solar Radiation Protections on Façades: A Case Study in a Hot Semi-Humid Climate

Climatically Adapted Piloti Arrangement and Ratio of Residential Blocks in a Subtropical Climate City

Experimental Evaluation of Natural Heat Transfer in Façade Integrated Triangular Enclosures

Modeling and Simulation of an Air-cooling Condenser under Transient Conditions

INFLUENCE OF FRICTION PENDULUM SYSTEM ON THE RESPONSE OF BASE ISOLATED STRUCTURES

Influence of Chord Lengths of Splitter Blades on Performance of Small Axial Flow Fan

INVESTIGATING GLAZING SYSTEM SIMULATED RESULTS WITH REAL MEASUREMENTS

Effect of Installing a Curved Venetian Blind to the Glass Window on Heat Transmission

Estimation of Hottest Spot Temperature in Power Transformer Windings with NDOF and DOF Cooling

Different Materials with High Thermal Mass and its Influence on a Buildings Heat Loss An Analysis based on the Theory of Dynamic Thermal Networks

cen29305_ch08.qxd 11/30/05 3:05 PM Page 451 INTERNAL FORCED CONVECTION CHAPTER 8 Liquid or gas flow through pipes or ducts is commonly used in heating

A Study on Quantitative Analysis Method for Project Bidding Decision

COMMERCIAL HEAT GAIN - HEAT LOSS CALCULATION LOAD FORM

Performance of Fixed Horizontal Shading Devices in South Facing Residential Buildings in Dhaka

TOPICAL PROBLEMS OF FLUID MECHANICS 17 ONE-DIMENSIONAL TEMPERATURE DISTRIBUTION OF CONDENSING ANNULAR FINS OF DIFFERENT PROFILES

Professor Haig Gulvanessian CBE Civil Engineering and Eurocodes Consultant, Visiting Professor, Imperial College London

5. AN INTRODUCTION TO BUILDING PHYSICS

Thermal mass vs. thermal response factors: determining optimal geometrical properties and envelope assemblies of building materials

Open Access Measuring Method for Diameter of Bearings Based on the Edge Detection Using Zernike Moment

Research Article Analysis Bending Solutions of Clamped Rectangular Thick Plate

CZECH TECHNICAL UNIVERSITY IN PRAGUE FACULTY OF MECHANICAL ENGINEERING DEPARTMENT OF ENVIRONMENTAL ENGINEERING

COVENANT UNIVERSITY NIGERIA TUTORIAL KIT OMEGA SEMESTER PROGRAMME: MECHANICAL ENGINEERING

Solar Control in Traditional Architecture, Potentials for Passive Design in Hot and Arid Climate

Simulation of Temperature Stratification in Atriums: Validation of Basic Model Features. Susan Bajracharya

NUMERICAL ANALYSIS OF HEAT STORAGE AND HEAT CONDUCTIVITY IN THE CONCRETE HOLLOW CORE DECK ELEMENT

The influence of solar radiation on the distribution of temperatures in historic masonry

Heat Loss from Cavity Receiver for Solar Micro- Concentrating Collector

Natural Ventilation. CFD modelling of a double-skin facade. Huw Birch. Supervisor: Abigail Hathway

Hygrothermal Performance of Ventilated Cold Roofs an Experimental Study

Applied Thermodynamics HEAT TRANSFER. Introduction What and How?

The Wind-Induced Vibration Response for Tower Crane Based on Virtual Excitation Method

Synthetic asymptote formula for surface-printed resistor

FUNDAMENTALS OF HVAC

Project 2. Introduction: 10/23/2016. Josh Rodriguez and Becca Behrens

Analyzing Sun Angles of Solar Time Dial for Design of Building Envelope Components for the Region of Andhra Pradesh, India

Chapter 2 HEAT CONDUCTION EQUATION

Chapter 2 HEAT CONDUCTION EQUATION

Investigations on effect of the orientation on thermal comfort in terraced housing in Malaysia

Study on Stack Effect of Stairwell by Numerical Model of Leakage Flow through Gap of Door

Open Access. Suman Chakraborty* Q T + S gen = 1 S 1 S 2. Department of Mechanical Engineering, Indian Institute of Technology, Kharagpur , India

Investigation of different wall profiles on energy consumption and baking time in domestic ovens

Chapter 1 INTRODUCTION AND BASIC CONCEPTS

Application of Time Sequence Model Based on Excluded Seasonality in Daily Runoff Prediction

A Numerical Analysis of Indoor Thermal Environment and Human Thermophysiological Responses under Natural Ventilation S. Iizuka 1,*, T. Sakoi 2, T. Sai

A Numerical Investigation on Active Chilled Beams for Indoor Air Conditioning

Open Access Experimental Research and Analysis of Vortex Excited Vibration Suppression of Spiral Stripe Strake

Experimental measurements and analysis of the indoor conditions in Italian museum storerooms: a study case

Analysis and Calculation of Double Circular Arc Gear Meshing Impact Model

Numerical Study of Heat Propagation in Living Tissue Subjected to Instantaneous Heating

Diurnal Change of Wall Flows and Energy Balance on the South Facing Wall When Background Wind is Absent

ARCHITECTURE IN THE DAYLIGHT

Characteristics of flow over the free overfall of triangular channel

THE EFFECTS OF CALORIMETER TILT ON THE INWARD-FLOWING FRACTION OF ABSORBED SOLAR RADIATION IN A VENETIAN BLIND

TRANSPARENT INNOVATIVE MATERIALS: ENERGETIC AND LIGHTING PERFORMANCES EVALUATION

On the Independence of the Formal System L *

Buckling Behavior of Long Symmetrically Laminated Plates Subjected to Shear and Linearly Varying Axial Edge Loads

Solar radiation and architectural design in Barcelona

Experiment Study on Rheological Model of Soft Clay

THERMAL TRANSMITTANCE OF MULTI-LAYER GLAZING WITH ULTRATHIN INTERNAL PARTITIONS. Agnieszka A. Lechowska 1, Jacek A. Schnotale 1

SIMULATION OF FRAME CAVITY HEAT TRANSFER USING BISCO v10w

Ceiling Radiant Cooling Panels Employing Heat-Conducting Rails: Deriving the Governing Heat Transfer Equations

Optimization of a Multi-pendulum Wave Energy Converter

Numerical simulation of high pressure gas quenching of H13 steel

Effect of Periodic Variation of Sol-air Temperature on the Performance of Integrated Solar Collector Storage System

Transcription:

68 The Open Civil Engineering Journal, 0, 5, 68-7 Open Access Study on Performance of Anisotropic Materials of Thermal Conductivity Shunyu Su *, Jian Chen and Chunzhi Zhang College of Uran Construction, Wuhan University of Science and Technology, Wuhan 430070, China Astract: The advantage and shortcoming of anisotropic material of thermal conductivity using as uilding envelope were analyzed in this paper. The thermal performance of anisotropic material of thermal conductivity was presented y rotating the original rectangular axes so as to determine the principal axes and make the cross-derivative terms disappear. The results indicated that the heat flux vector is commonly not perpendicular to the isothermal surface in anisotropic material. Doule skin facade has een successfully applied in many uilding designs. The application of laminated composite with anisotropic thermal conductivity in doule skin envelope was proposed to avoid its disadvantage. The interior envelope in doule skin system may e made of laminated glaze or other laminated materials. Basing on the comination of anisotropic material and doule skin envelope, the indoor cooling and heating load decrease in summer and winter respectively while the anisotropic material was used as interior uilding envelope. Especially in summer, the effect of energy saving is ovious since the partial magnitude of conductive heat in the interior envelope could e rought out from the cavity y natural ventilation through it. Keywords: Anisotropic material, doule skin façade, laminated composite, thermal conductivity.. INTRODUCTION Thermal conductivity of a material k has great influence on its thermal performance. It is a measure of how fast heat will flow in that material. In heat transfer analysis, a material is normally assumed to e isotropic. It has uniform properties in all directions. Therefore, it is common practice to evaluate the thermal conductivity k at the average temperature and treat it as a constant in calculations []. This assumption is realistic for most materials, ut some materials, such as laminated composite materials and wood, they are anisotropic and exhiit different structural characteristics in different directions. For example, the thermal conductivity of wood across the grain is different than that parallel to the grain. In such cases the thermal conductivity k should e expressed as a tensor uantity to account for the variation with direction. Shiah and Tan [-3] developed an exact transformation approach that can e applied to numerically treat the two-dimensional and even three-dimensional anisotropic prolems in field theory y the oundary element method (BEM). The thermal conductive theory of anisotropic materials plays an important role in many areas of science and technology. But works aout it are still limited and relatively few reported. Doule skin envelope or doule skin facade (DSF) is widely applied in uilding and architecture area. Experimental and analytical studies have een carried out on the thermal characteristics of DSF [4-6]. Pasuay [7] presented the advantages and disadvantages of DSF on the asis of actual measurements of high-rise uildings. Ding et al. [8] finished *Address correspondence to this author at the College of Uran Construction, Wuhan University of Science and Technology, No. 5 mail ox, No. 99 Xiongchu Avenue, Wuhan 430070, China; Tel: +86-07-6889356; Fax: +86-07-6889364; E-mail: shunyusu@63.com reduced-scale model experiments and a computational fluid dynamics analysis to evaluate the natural ventilation performance in the prototype of an eight-storey uilding with DSF in the south. Wong et al. [9] reported that DSF used on the eastern and western walls of uilding can effectively conserve cooling energy y shielding uildings from solar radiation and reducing the heat y natural ventilation. Nowadays studies aout DSF are undertaken in different climatic conditions such as hot and humid; moderate; cold and hot arid climates [0-3]. Moreover, detailed information on the thermal distriution in DSF is necessary for uilding envelope so as to not only provide comfortale indoor environment ut also save energy. The purpose of this paper is to analyze the thermal conductive characteristics of anisotropic materials, and to discuss the application of anisotropic material of thermal conductivity in doule skin envelope system so as to make use of its advantage and avoid its disadvantage.. CONDUCTIVE MECHANISM OF ANISOTROPIC MATERIALS OF THERMAL CONDUCTIVITY Comparing with isotropic material, the conductive process of anisotropic material has two important differences. One is that the thermal conductivity of anisotropic material is different in different directions. The other is that the magnitude of heat transfer rate of anisotropic material in one direction is not only relevant to the temperature gradient in this direction, ut also relevant to the temperature gradient which is normal to the direction [4]. In rectangular coordinates (x, ) for two-dimensional conduction, the magnitudes of heat flux x and in x and directions can e determined from Fourier s law as 874-495/ 0 Bentham Open

Study on Performance of Anisotropic Materials of Thermal Conductivity The Open Civil Engineering Journal, 0, Volume 5 69 x = k k () = k k () where t represents the temperature in the anisotropic material, and k ij, (i, j=, ) are the thermal conductivities, where from thermodynamic principles and Onsagar s reciprocity relation, the coefficients must satisfy the following relations: k 0, k 0, k k k 0, k =k (3) To solve the anisotropic prolem is evidently a more complicated process since k ij are invariant coefficients with respect to x and axes in rectangular coordinate system. A commonly adopted approach to treat the anisotropic prolem is to determine the principal axes (, ) y rotating the original rectangular axes so that the cross-derivative terms will disappear. In other words, in rectangular coordinates (, ) for two-dimensional conduction, the magnitude of heat flux in the i (i=, ) direction is relevant to the temperature gradient in this direction, ut is not relevant to the temperature gradient in orthogonal axis. The magnitudes of heat flux can also e expressed from Fourier s law as = k (4) Since temperature gradient and heat flux are not relevant to coordinate system in the same oject, the following may e otained: x = cos sin (6) = sin + cos (7) = cos sin (8) x = sin + cos (9) From E. (6) to E. (9), the following may e deduced: x = ( k cos + k sin ) ( k x k )sin cos (0) = ( k k )sin cos ( k x sin + k cos ) () If the two surfaces of the laminated composite maintain constant temperatures of t w and t w. There will e one- dimensional temperature field in the, and x = 0 () = k (5) = t t w w (3) where k and k represent the principal thermal conductivities. Fig. () shows the two-dimensional steady-state conduction of an anisotropic material. There is a laminated composite with anisotropic conductivity. The inclined lines in the figure express laminated surfaces. The directions along the laminated surface and its normal are thermal conductivity principal axes and. The thermal conductivities along the principal axes and are k and k. The x axis represents the direction of isothermal surface. Rectangular coordinates (x, ) are inconsistent with the principal axes and. represents the anisotropic angle in the figure. x x x x where represents thickness of the laminated composite. So that x = ( k k )sin cos = ( k k )sin cos t t w w = ( k sin + k cos ) = ( k sin + k cos ) t t w w (4) (5) Euations (4) and (5) indicate, for anisotropic material, the temperature gradient in one direction develops magnitude of heat transfer rate in orthogonal direction. 3. ADVANTAGE OF ANISOTROPIC MATERIALS OF THERMAL CONDUCTIVITY O x Fig. (). Two-dimensional steady-state conduction of the anisotropic material. x Basing on aove analysis, it can e concluded that the heat flux vector is commonly not perpendicular to the isothermal surface in anisotropic material. For thermal conduction in a as Fig. () shows, the conductive heat flux is more concerned. The apparent thermal conductivity is defined

70 The Open Civil Engineering Journal, 0, Volume 5 Su et al. k = k sin + k cos (6) Therefore, = k t w t w (7) For thermal conduction in the made of isotropic material, the conductive heat flux across it is expressed as ' = k t w t w (8) k k If k =k, then the apparent thermal conductivity =. The thermal conductive material ecomes isotropic. Ordinarily, 0k k. Tale shows the apparent thermal conductivity k varies with the anisotropic angle. For special case, while k =0, k = k sin k ; while =0, k = k. Since the value of k is more lower than that of k, the thermal conductive material ecomes a poor heat conductor or insulator; while =90, k = k. The thermal conductive performance is very good for this material, or this material is a good heat conductor, ut the insulating effect is conseuentially worse for it. The results of Tale show that generally the apparent thermal conductivity k ' < k. There fore,. In other words, if the thermal conductivity of isotropic material in thickness direction is eual to that of anisotropic material in laminated direction, the magnitude of heat flux across the made of anisotropic material is less than that made of isotropic material. This advantage of anisotropic material can e made use in uilding envelope. Fig. () shows the conduction of uilding envelope made of anisotropic material in summer while the outdoor temperature is higher than that of indoor. The heat transfer rate vector is not perpendicular to the uilding wall. It will develop partial magnitude paralleling to the wall. Then the real heat flux across the uilding envelope from outdoor to indoor decreases when the wall temperatures of two sides maintain constant respectively. If 0k k, it can easily come true y laying out laminated composite material in reasonale direction, the partial magnitude of heat transfer rate will flow upwards. Finally, it will spread to the amient on the roof of the single-storey uilding y thermal convection with the around air. Therefore, the indoor cooling load in summer decreases while the anisotropic material is used as uilding envelope. Fig. (3) shows the conduction of uilding envelope made of anisotropic material in winter while the indoor temperature is higher than outdoor. The real heat flux across the uilding envelope from indoor to outdoor decreases when the wall temperatures of two sides maintain constant respectively. The heat transfer rate vector will develop partial magnitude paralleling to the wall. If 0k k, the partial magnitude of heat transfer rate will flow downwards. It can partially spread to the indoor air y thermal convection. Then, the indoor heating load in winter may decrease a little while the anisotropic material is used as uilding envelope. Fig. (). Conduction of the anisotropic material in summer. x x x Outdoor O x x Outdoor x Fig. (3). Conduction of the anisotropic material in winter. 4. DOUBLE SKIN ENVELOPE SYSTEM AND LAMI- NATED COMPOSITE ANISOTROPIC MATERIALS Although the anisotropic material of thermal conductivity made of laminated composite has its advantage used in x x x x x x Tale. The Apparent Thermal Conductivities Varying with the Anisotropic Angles 0 30 45 60 90 k k k + 3k 4 k + k 3k + k 4 k

Study on Performance of Anisotropic Materials of Thermal Conductivity The Open Civil Engineering Journal, 0, Volume 5 7 uilding envelope, it is still imperfect when it used in multistorey uilding or high-rise uilding. Since partial magnitude of heat transfer rate will flow upwards in laminated composite material in summer, for single-storey uilding, this part of heat will spread to the amient on the roof of the uilding y thermal convection with the around air; ut for multi-storey uilding or high-rise uilding, it may spread to the indoor of upper storey y thermal convection with the indoor air when it flows upwards in the uilding envelop. In order to avoid this shortcoming, it could e related to doule skin facade in uilding. The techniue of doule skin facade has een successfully applied in many uilding designs. Doule skin facade using as the uilding envelope consists of an outer glazed skin and an inner glazed or partially glazed one. The two skins together with top and ottom dividing s form a cavity. This air layer works as a thermal uffer zone. And this kind of uilding envelope is proposed not only to reduce noise from outdoors, ut most importantly, to decrease indoor heating load in winter y increasing outside temperature of the inner facade, using solar incidence, and decrease indoor cooling load in summer y lowering the temperature via ventilation in the cavity. Fig. (4) shows schematic diagram of this doule skin envelope in summer. The top and ottom dividing s are open in summer, and outdoor air is induced into the cavity and naturally ventilated through it y uoyancy arising from the outside and inside temperature difference of the cavity. Heat in the cavity and from the two skins spread out in this convection process. Fig. (5) shows schematic diagram of this doule skin envelope in winter. The top and ottom dividing s are close. Heat from sunshine stores in the cavity and then transfers to indoor through the interior envelope. This uilding envelope system is considered to e an environment-friendly design since it utilizes natural ventilation and sunshine so that comfort can e achieved with less energy consumption [5]. Doule skin envelope used in summer could avoid the disadvantage of anisotropic material of thermal conductivity y ventilation in its cavity. Ventilation in the cavity of this doule skin uilding envelope system can e natural, mechanical or a comination of oth [6]. Outside Exterior glass Interior envelope Floor Air outlet Fig. (5). Schematic diagram of doule skin envelope in winter. Outside Exterior glass Air inlet Floor Interior envelope Fig. (4). Schematic diagram of doule skin envelope in summer. 5. APPLICATION OF LAMINATED COMPOSITE WITH ANISOTROPIC THERMAL CONDUCTIVITY IN DOUBLE SKIN ENVELOPE The selection of uilding material has important effect to the design of uilding thermal performance and energy saving. According to aove analysis, laminated composite with anisotropic thermal conductivity could e applied in doule skin envelope. As Fig. () shows, the partial magnitude of heat transfer rate in the uilding envelope with anisotropic material will flow upwards in summer. It will spread to outside on the roof of the single-storey uilding y thermal convection with the around air. But in multi-storey or high-rise uilding, this partial magnitude of heat transfer rate will flow upstairs to the interior envelope of higher floor. It may transfer to indoor of higher floor y convection etween the interior envelope surface and indoor air of the higher floor. Thus the indoor cooling load of the higher floor will increase. To avoid this disadvantage, the laminated composite with anisotropic thermal conductivity could e used as the mate-

7 The Open Civil Engineering Journal, 0, Volume 5 Su et al. rial of interior envelope in doule skin envelope. The interior envelope may e made of laminated glaze or other laminated material as Fig. () or Fig. (3) shows. The top and ottom dividing s of the doule skin envelope will e open in summer. The heat in the surface of interior envelope will e rought out y flowing air through the cavity. And in winter, the dividing s will e close. In this case, the heat in the surface of interior envelope can only transfer to the air in the close cavity. But this heat transfer is limited. In fact, the cavity is warmer than outside. And with solar incidence on the facade of exterior envelope, its temperature increases. Comparing with outside temperature in winter, the higher cavity temperature has a significant role in reducing heat loss from the uilding. It can e concluded that this comined doule skin envelope system, while the interior envelope is made of laminated composite with anisotropic thermal conductivity, will lower indoor cooling load in summer as well as indoor heating load in winter. The effect of energy saving is ovious in summer since the partial magnitude of conductive heat in the interior envelope could e rought out from the cavity y natural ventilation through it. 6. CONCLUSIONS The thermal performance of conductivity anisotropic material was analyzed in this paper. The characteristics in anisotropic materials are different from those in isotropic materials. The heat flux vector is commonly not perpendicular to the isothermal surface in anisotropic material. The magnitude of heat transfer rate of anisotropic material in one direction is not only relevant to the temperature gradient in this direction, ut also relevant to the temperature gradient which is normal to the direction. The doule skin facade in winter preheats the interior space of the envelope and decreases the heating load. And in summer, the ventilation of the cavity can effectively lower the cooling load. The interior envelope in doule skin system may e made of laminated glaze or other laminated material with anisotropy. Basing on the comination of anisotropic material and doule skin envelope, the indoor cooling and heating load decrease in summer and winter respectively while the anisotropic material was used as interior uilding envelope. CONFLICTS OF INTEREST None Declared 7. ACKNOWLEDGMENTS The authors acknowledge the financial support from Doctoral Science Research Foundation of Wuhan University of Science and Technology. The authors also would like to thank the Organizing Committee of the nd International Conference on Manufacturing Science and Engineering (ICMSE 0). The paper was revised according to their comments. REFERENCES [] Y. A. Cengel, Heat transfer: a practical approach, nd ed., New York: McGraw Hill, 003. [] Y. C. Shiah and C. L. Tan, "BEM treatment of two-dimensional anisotropic field prolems y direct domain mapping", Engng. Anal. Bounda. Elem., vol. 0, pp. 347-35, 997. [3] Y. C. Shiah and C. L. Tan, "BEM treatment of three-dimensional anisotropic field prolems y direct domain mapping", Engng. Anal. Bounda. Elem., vol. 8, pp. 43-5, 004. [4] H. Manz, "Total solar energy transmittance of glass doule facades with free convection", Energy Build., vol. 36, pp. 7-36, 004. [5] E. Gratia, A. D. Herde, "The most efficient position of shading devices in a doule-skin facade", Energy Build., vol. 39, pp. 364-373, 007. [6] T. Hiroaki, O. Masaya, T. Hideki, Y. Y. Gyu, W. Kenji, "Thermal characteristics of a doule-glazed external wall system with roll screen in cooling season", Build. Environ., vol. 44, pp. 509-56, 009. [7] T. Pasuay, "Natural ventilation in high-rise uildings with doule facades, saving or waste of energy", Energy Build., vol. 36, pp. 38-389, 004. [8] W. T. Ding, Y. Hasemi, T. Yamada, "Natural ventilation performance of a doule-skin facade with a solar chimney", Energy Build., vol. 37, pp. 4-48, 005. [9] N. H. Wong, L. Wang, N. C. Aida, R. P. Anupama, X. Wei, "Effects of doule glazed facade on energy consumption, thermal comfort and condensation for a typical office uilding in Singapore", Energy Build., vol. 37, pp. 563-57, 005. [0] P. C. Wong, D. Parsad, M. Behnia, "A new type of doule skin facade configuration for the hot and humid climate", Energy Build., vol. 40, pp. 94-945, 008. [] E. Gratia, A. D. Herde, "Optimal operation of a south doule-skin facade", Energy Build., vol. 36, pp. 4-60, 004. [] R. Hossegen, B. J. Wachenfeldt, S. O. Hanssen, "Building simulation as an assisting tool in decision making, case study: with or without a doule-skin facade?", Energy Build., vol. 40, pp. 8-87, 008. [3] N. Hashemi, R. Fayaz, M. Sarshar, "Thermal ehaviour of a ventilated doule skin facade in hot arid climate", Energy Build., vol. 4, pp. 83-83, 00. [4] Q. S. Yang and B. Y. Pu, Advanced heat transfer, nd ed., Shanghai: Shanghai Jiao Tong University Press, 00. (In Chinese) [5] Y. M. Kim, S. Y. Kim, S. W. Shin and J. Y. Sohn, "Contriution of natural ventilation in a doule skin envelope to heating load reduction in winter", Build. Environ., vol. 44, pp. 36-44, 009. [6] D. Saelens, S. Roels, H. Hens, "Strategies to improve the energy performance of multiple-skin facades", Build. Environ., vol. 43, pp. 638-650, 008. Received: Octoer 0, 00 Revised: Feruary 4, 0 Accepted: July 30, 0 Su et al.; Licensee Bentham Open. This is an open access article licensed under the terms of the Creative Commons Attriution Non-Commercial License (http://creativecommons.org/licenses/ y-nc/3.0/) which permits unrestricted, non-commercial use, distriution and reproduction in any medium, provided the work is properly cited.