Determination of moisture surface transfer coefficients under transient conditions

Similar documents
Influence of Air Space on Multi-layered Material Water Vapor Permeability Measurement

WUFI Workshop at NTNU /SINTEF Fundamentals

Liquid water is one of the

Effect of moisture transfer on heat energy storage in simple layer walls

HYGROTHERMAL CHARACTERISTICS OF PUMICE AGGREGATE CONCRETE USED FOR MASONRY WALL BLOCKS

Hygrothermal Properties and Performance of Sea Grass Insulation

MOISTURE PERMEABILITY DATA PRESENTED AS A MATHEMATICAL FUNCTION APPLICABLE TO HEAT AND MOISTURE TRANSPORT MODELS

Chapter 2 Experimental Studies on Hygrothermal Behaviour of ETICS

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

CHAPTER 3. The sun and the seasons. Locating the position of the sun

Drying out capacity and snow melting risk for ventilated wooden roofs - a parameter study

NON STATIONARY DIFFUSION TUBE Functional Sample Diffusion Tube

1D-HAM. Coupled Heat, Air and Moisture Transport in Multi-layered Wall Structures. Manual with brief theory and an example. Version 2.

QIRT th International Conference on Quantitative InfraRed Thermography

123MEAN thermal properties KATEDRA MATERIÁLOVÉHO INŽENÝRSTVÍ A CHEMIE

5. AN INTRODUCTION TO BUILDING PHYSICS

5. Thermal Design. Objective: Control heat flow to: Maintain comfortable indoor conditions

In situ quantitative diagnosis of insulated building walls using passive infrared thermography

Faculty of Mechanical Engineering Institute of Wood and Paper Technology

Evaluation of Comsol as a tool for pinpointing moisture entering locations from inside surface moisture

Daily Hygroscopic Inertia Classes: Application in a Design Method for the Prevention of Mould Growth in Buildings

Durability and hygrothermal performance of building envelope

Regional Climate Simulations with WRF Model

above the land to be warmer, causing convection currents because of the difference in air pressure.

Moisture content in concrete

HYGROTHERMAL BEHAVIOUR OF A HEMP CONCRETE WALL: COMPARISON BETWEEN EXPERIMENTAL AND NUMERICAL RESULTS

Vincent Barraud SOPREMA BASICS OF THERMAL INSULATION

THERMAL DYNAMIC BEHAVIOR OF SINGLE COMPONENT: EXPERIMENTAL ANALYSIS AND NUMERICAL MODELING

Hygrothermal Performance of Ventilated Cold Roofs an Experimental Study

Technical University of Denmark BSIM

Comparison of probabilistic approaches to the mould growth problem

Influence of a Refurbishing on an Old Residential building's Wall in Paris in Summer: Mass and Heat Transfer Approach

The article analyzes the process of vapor transfer through building materials, taking into account their real operation

International Fire Safety Symposium 2015

Determination of water and salt transport parameters of porous materials using methods of inverse modelling

Institut national des sciences appliquées de Strasbourg GENIE CLIMATIQUE ET ENERGETIQUE APPENDICES

Temperature (T) degrees Celsius ( o C) arbitrary scale from 0 o C at melting point of ice to 100 o C at boiling point of water Also (Kelvin, K) = o C

Thermodynamics I Spring 1432/1433H (2011/2012H) Saturday, Wednesday 8:00am - 10:00am & Monday 8:00am - 9:00am MEP 261 Class ZA

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

Assoc. Prof. Dr. Ruta MINIOTAITE

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

Surface Temperatures on Flat Roofs and Hygrothermal Consequences

Section 2 of 6 Applied Psychrometrics

Fundamentals of WUFI-Plus WUFI Workshop NTNU / SINTEF 2008

A Dynamic Method Determining Equilibrium Moisture Content of Porous Materials

CIGRE SC A2 & D1 JOINT COLLOQUIUM 2011, KYOTO JAPAN. Water Saturation Limits and Moisture Equilibrium Curves of Alternative Insulation Systems

IGNITABILITY ANALYSIS USING THE CONE CALORIMETER AND LIFT APPARATUS

Standard Test Method for Measuring the Steady-State Thermal Transmittance of Fenestration Systems Using Hot Box Methods 1

EE-489 Modeling Project: Behavior of the. Temperature Changes in a Room

Projected change in the East Asian summer monsoon from dynamical downscaling

Requirements to perform accurate dielectric material analysis

Documentation of the Solutions to the SFPE Heat Transfer Verification Cases

ON THE CONSERVATION OF MASS AND ENERGY IN HYGROTHERMAL NUMERICAL SIMULATION WITH COMSOL MULTIPHYSICS

AN EXPERIMENTAL STUDY OF THE FROST FORMATION ON A COLD SURFACE IN FREE CONVECTIVE FLOW

Correlation between Small-Scale Rate of Heat Release and Full-Scale Room Flashover for Surface linings

Mr. P s Science Test!

TRANSSOLAR Energietechnik GmbH Strategies for energy efficient design and thermal comfort in buildings

Design for Airtightness and Moisture Control in New Zealand housing

ARCH 348 BUILDING AND ENVIRONMENTAL SYSTEMS

NFRC THERMAL TEST SUMMARY REPORT Expiration Date: 03/31/06

MODELLING OF FREEZING AND THAWING PROCESSES IN BUILDING MATERIALS

Chapter 2 Inputs for Hygrothermal Simulation Tools

Validation of a coupled heat, vapour and liquid moisture transport model for porous materials implemented in CFD

AquaFlux A New Instrument for Water Vapour Flux Density Measurement

3 rd Annual CHAMPS Developers Workshop

Matter In Our Surroundings

Computational Modelling of the Impact of Solar Irradiance on Chemical Degradation of Painted Wall Hangings in an Historic Interior

AN IMPACT OF AIR INFILTRATION ON HEAT TRANSFER THROUGH LIGHTWEIGHT BUILDING PARTITIONS FILLED WITH LOOSE MINERAL WOOL PHD THESIS ABSTRACT

Computer Evaluation of Results by Room Thermal Stability Testing

Vantage Pro Technical Reference

Specific Heat Measurement of High Temperature Thermal Insulations by Drop Calorimeter Method

/qirt The influence of air humidity on effectiveness of heat sink work. by M. Kopeć*, R. Olbrycht*

CLASSIFIED 2 PRESSURE THERMAL PHYSICS MR. HUSSAM SAMIR

Climate versus Weather

The energy performance of an airflow window

Urea Injection and Preparation in Diesel Applications Multiphase Multicomponent Modeling using Star-CD

Science 1206 Chapter 1 - Inquiring about Weather

Fig 1. Power Tower during Operation

A SIMPLE MODEL FOR THE DYNAMIC COMPUTATION OF BUILDING HEATING AND COOLING DEMAND. Kai Sirén AALTO UNIVERSITY

R13 SET - 1 '' ''' '' ' '''' Code No RT21033

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

Precipitation processes in the Middle East

DYNAMIC INSULATION APPLIED TO THE RESIDENTIAL BUILDING (PART 2) Numerical Evaluation of Thermal Insulation Effect on Air Supply Window System

Name Period 4 th Six Weeks Notes 2013 Weather

IV. Atmospheric Science Section

Influence of water on the total heat transfer in evacuated insulations

MOULD GROWTH DAMAGES DUE TO MOISTURE: COMPARING 1D AND 2D HEAT AND MOISTURE MODELS?

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

SPH3U1 Lesson 03 Energy

Developments at DWD: Integrated water vapour (IWV) from ground-based GPS

Today s Weather 24 Hour Forecaster

Water vapour above Switzerland over the last 12 years

Prediction of Thermal Comfort and Ventilation Efficiency for Small and Large Enclosures by Combined Simulations

Three things necessary for weather are Heat, Air, Moisture (HAM) Weather takes place in the Troposphere (The lower part of the atmosphere).

LECTURE 5 THE CONTENTS OF THIS LECTURE ARE AS FOLLOWS: 1.0 SOME OF THE THEORETICAL CONCEPTS INVOLVED IN HEAT FLOW. 1.1 Sensible Heat Flow

VALIDATION OF REYNOLDS AVERAGED MODEL AND LARGE EDDY SIMULATION IN ACTUAL FLOOR HEATING ROOM. Hiroki Ono 1 and Koji Sakai 1

Part 9: Shrinkage for service and accident conditions

FUNDAMENTALS OF HVAC

STATIC GAS MONITOR Type SGM/DEW Revision A of 20 aprile 2016

INFLUENCE OF HYSTERETIC MOISTURE TRANSFER ON CONCRETE DURABILITY

Transcription:

Determination of moisture surface transfer coefficients under transient conditions T. Bednar & J. Dreyer University of Technology Vienna, Vienna, Austria ABSTRACT: Detailed investigations on surface transfer coefficients for indoor environments under nonstationery conditions, like short term airing, are presented in the paper. The measurement method is carefully investigated to minimize the influence of uneven surface temperature profiles. A small (4x4x2 cm³) insulated specimen with a high liquid conductivity is the best probe for determination of the surface transfer coefficients. From the experiments it can be concluded that under winter conditions the moisture transfer coefficient is increased by 180% from its value with closed windows and that it takes two hours after the windows have been closed till it reaches its old value. 1 INTRODUCTION Determination of Indoor Climate always includes the prediction of thermal and hygral conditions of the indoor air and at the surfaces. Already developed simulation models for the performance of buildings are able to predict those parameters. User behaviour (heat and moisture production, ventilation) and the dynamics of heat and moisture inside materials can be modelled with an acceptable accuracy and input data can be found in literature. Standards recommend different heat transfer coefficients for heat and moisture calculations and neglect surface moisture resistances. Still very rare and incomplete are well documented datasets on the surface transfer coefficients for heat and moisture transport. A documentation of existing models and measurements together with measurements on the dependence of the surface transfer coefficients on the air flow velocity can be found in Worch 2002. To measure realistic surface transfer coefficients a test room has been developed where different heat and moisture production and/or ventilation strategies can be run. During such non-stationary experiments the sorption of moisture at surfaces, the surface transfer coefficients and the dynamics of the room and the outdoor climate are measured. The mass loss of a small specimen (4x4x2 cm³) is determined with a balance. Additionally the surface temperature in the centre of the specimen and the temperature and the relative humidity of the indoor air are measured next to the specimen. The test-room has an area of 13.4 m² and a height of 3.5 m. The east facing exterior wall has one window (1.32 x 1.95 m²). 2 MEASURMENT SETUP The measurement setup for the determination of the surface transfer coefficients is presented in figure 1. Figure 1. Test-Room and measurement setup

3 MATHEMATICAL MODEL The evaporation of water from surfaces can be modelled with the following equation. m& = Aevap β' p (1) where m& = mass loss of specimen; A evap = evaporating area; β ' = moisture transfer coefficient and p = partial vapour pressure difference between room and surface. The density of heat flow across the surface is determined by the convective and radiative heat transfer and the heat which has to be used for evaporation. As the radiative heat transfer is depending on the geometry and the temperature of all surfaces the equation is simplified using an effective heat transfer coefficient. m& q = α ( Troom Tsurf ) hlv (2) Aevap where q =density of heat flow, α =effective heat transfer coefficient, T room =room temperature, T surf =surface temperature and h lv =enthalpy of liquidvapour phase change. To calculate the time-response of a specimen with an insulation ring a three dimensional program has been developed which accounts for the heat flow due to evaporation according to Equation 2. Figure 2 shows the geometry and the network which is used for calculating the heat flows. Not accounted for is the convective heat transport inside the specimen due the liquid water movement or vapour transport. Figure 2. Geometry of model and network for the calculation of the heat flows systematic error of the calculated moisture transfer coefficient. Figure 3 shows the measured temperature profile across the evaporating area. The room temperature has been 22.5 C and the specimen has been insulated with 4cm extruded polystyrene. The specimen had 10 hours to reach equilibrium. The measurement of the surface temperatures has been done with an infrared thermography system SC1000 from Inframetrics calibrated with the temperature at the centre of the specimen measured with thermocouple. As can be seen the temperature profile depends on the material. Temperature C 18.5 18.3 18.1 17.9 17.7 17.5 Profile CaSi AAC 0 0.01 0.02 0.03 0.04 Position m calc. Figure 3. Measured and calculated temperature profile across surface for two different materials. The boundary conditions used for the calculations are α = 12.7W / m² K, 5 β ' = 18 10 kg / m² spa. T room = 22. 5 C With a calcium-silicate fibre board (CaSi, max. moisture content 850 kg/m³), which has a high moisture conductivity, the profile is more uniform comparing to the profile measured with autoclaved aerated concrete (AAC), which has a lower moisture conductivity. The calculated profile agrees very well with the measured profile on the AAC surface. As can bee seen the convective heat transport due to liquid transport flattens the profile and using CaSi as the moisture source is appropriate. 4 SURFACE TEMPERATURE PROFILE As the evaporation of water depends on the partial vapour pressure at the surface, which itself depends on the surface temperature, a uniform temperature across the evaporating surface is important. Usually the surface temperature is measured at one point of the surface. An uneven distribution would result in a 5 STEP RESPONSE AND ACCURACY To estimate the accuracy of the measured transfer coefficients the step response of the system is calculated for two different scenarios. In case 1 the transfer coefficients are varied and in case 2 the room temperature is varied stepwise.

Mean Massflowrate 10-5 kg/m² s 4.0 3.5 3.0 2.5 2.0 1.5 4x4x 2.5 cm³ specimen 5cm insulation one-dimensional drying 1.0 Temperature C 23 21 19 17 15 T room T surf 11.5 11.5 11.0 10.5 10.0 9.5 "measured" transfer coefficient boundary condition 11.0 10.5 10.0 9.5 "measured" transfer coefficient boundary condition 9.0 9.0 Figure 4. Step-response of the system for case 1. The transfer coefficient is varied stepwise between 10 10-5 and 11 10-5 kg/m² h Pa. Initial temperature: 23 C. The temperature of the surface is measured in the centre of the drying area. Figure 5. Step-response of the system for case 2. The room temperature is varied stepwise between 22.5 C and 18 C. Initial temperature: 19 C In Figure 4 the results for the simulated measurement are presented for case 1. The uneven temperature field leads to a systematic deviation of 3%. As the temperature in the centre of the drying area is a slightly lower than the mean value the vapour pressure difference is underestimated and therefore the measured transfer coefficient is overestimated. Additionally it can be seen that if the specimen starts with an initial temperature equal to the room temperature it takes about 4 hours until equilibrium is reached. Due to the geometry, the additional heat flow across the insulated surfaces increases the surface temperature and the mass flow rate. It has been shown that, increasing the area to 20x20cm², results in a more uneven temperature distribution across the drying surface (Bednar 2002)and therefore the systematic deviation of the measured transfer coefficient increases too. In Figure 5 the results of case 2 are presented. The stepwise variation of the room temperature results in variation of the surface temperature. As different parts of the surface react in a different way, the relation between the temperature in the centre of the area and the mean temperature changes. If the room temperature increases rapidly the deviation of the measured transfer coefficient increases to 5%. From the calculations above one can conclude that if the mean surface temperature is estimated from the temperature in the centre of the drying area the systematic error is less than 5%, especially if the material has high liquid water conductivity.

6 RESULTS FOR SHORT TERM AIRING 35 The measurement setup has been used to measure the moisture transfer coefficient during short term airing of the test room. The initial temperature of the test room was 22.5 C. The outdoor conditions during the measuring period were 8 C, 95% rel. humidity, wind: W, 1km/h, rainfall. The window was opened to an angle of 90 during 30 minutes. The recorded mass loss of the specimen and the measured difference of the partial vapour pressure are shown in Figure 6. 0.0 1200 30 25 20 15 10 5 0 open Window 0 1 2 3 4 Mass loss kg/m² -0.1-0.2-0.3-0.4-0.5-0.6 open Window 0 1 2 3 4 1000 800 600 400 200 0 p Pa Figure 6. Measured mass loss of specimen and partial vapour pressure difference during test period. Outdoor conditions 8 C and 94% relative humidity. The initial moisture content was 700 kg/m³ to assure that the surface stays wet during the measurement. Due to the temperature difference between room and outdoor environment the air flow during the period with the open window influenced the balance in a way that the determination of the time-dependent mass loss during that period was not possible. Therefore only the mass loss over 30 minute intervals was used to determine the mean mass flow rate. The mean partial vapour pressure difference was calculated out of the measured surface temperature and room conditions. As shown in Figure 7 the moisture transfer coefficient is nearly doubled during intensive airing of a room. After the window has been closed the air flow in the room is still higher and the transfer coefficient is increased. This could be explained with the disturbed temperature fields in the room which induce air movement until equilibrium is reached again. Already published results for the transfer coefficients (Dreyer et.al. 2002) with closed, tilted and Figure 7. 30 minutes average of the moisture transfer coefficient during the test period. open windows showed lower transfer coefficients for closed (12 10-5 kg/m² h Pa) and constant open windows (16 10-5 kg/m² h Pa) under spring weather conditions (outdoor temperature 18 C). This could be explained with the lower temperature differences between indoor air, window surface and outdoor air resulting in lower air flows. 7 CONCLUSIONS The measurement of the surface transfer coefficients can be done with a small specimen with high liquid moisture conductivity. The systematic error due the uneven surface temperatures is less than 5%. With closed windows and still air the moisture transfer coefficient is around 18 10-5 kg/m² h Pa. During short term airing under winter conditions the moisture transfer coefficient is increased to 180% and it takes about two hours till it reaches the old value. REFERENCES Bednar, T., 2002. Approximation of liquid moisture transport coefficient of porous building materials by suction and drying experiments. Demands on determination of drying curve. In Building Physics in the Nordic Countries; Proc. intern. symp., Trondheim, 2002. Dreyer, J. & Bednar, T. & Jung, M., 2002. Simulation der Feuchtedynamik in Räumen Einfluß der Übergangswiderstände unter instationären Bedingungen. 11. Bauklimatisches Symposium; Proc. Intern. Symp. Dresden 2002. Worch, A., 2002. Überlegungen und Versuche zur Erfassung des Wasserdampfüberganges an Bauteiloberflächen. Dissertation, Universität Dortmund.