) ANALYTICAL AND NUMERICAL STUDIES OF HEAT AND MOISTURE TRANSFER THROUGH POROUS INSULATION By BI-FENG ZHENG B. Sc. (Thermal Science and Engineering) Zhejiang University, China M. Sc. (Thermal Science and Engineering) Zhejiang University, China A THESIS SUBMITTED IN PARTIAL FULFILLMENT OF THE REQUIREMENTS FOR THE DEGREE OF MASTER OF APPLIED SCIENCE in THE FACULTY OF GRADUATE STUDIES DEPARTMENT OF MECHANICAL ENGINEERING We accept this thesis as conforming to the required standard THE UNIVERSITY OF BRITISH COLUMBIA October 1993 © BI-FENG ZHENG, 1993 In presenting this thesis in partial fulfilment of the requirements for an advanced degree at the University of British Columbia, I agree that the Library shall make it freely available for reference and study. I further agree that permission for extensive copying of this thesis for scholarly purposes may be granted by the head of my department or by his or her representatives. It is understood that copying or publication of this thesis for financial gain shall not be allowed without my written permission. Department of Mechanical Engineering The University of British Columbia 2324 Main Mall Vancouver, Canada V6T 1Z4 Date:^ 0 cto bex^9) Abstract This work contains both analytical and numerical studies of heat and moisture transport through a porous insulation in the presence of condensation, with impermeable, adiabatic vertical boundaries, and with one horizontal boundary facing a warm humid ambient and the other facing a cold impermeable surface. The analytical model is developed for heat and water vapor transfer in flat-slab and round-pipe thermal insulations. The model is validated by comparing its predictions with available experimental data. The effective thermal conductivity of the insulation in the presence of condensation depends on seven design and operating variables. The effect of these variables is determined by a parametric study. For practical operating conditions, the effective thermal conductivity varies from about 1.5 to 15 times the drystate value. The computed data are presented in the form of design curves which may be used to estimate the effective thermal conductivity for fiat-slab and round-pipe insulation systems. The analysis in this work quantifies the process of energy and mass transport in a porous insulation. A rigorous and fundamental formulation of heat and mass transfer in the insulation system is presented. The problem is modeled as one-dimensional, transient, multiphase flow with variable properties. Four stages in the energy and moisture transport process are identified, and they are formulated by a system of transient intercoupled equations and several thermodynamic relations using a local volume averaging technique. The numerical results are compared with experimental data for five different operating conditions and for times up to 600 hours. The model predicts the temperature distribution, heat transfer rate, the total moisture gain successfully. The predicted liquid distributions agree well with measured data for a period of up to 70 hours. The interesting effects of pertinent parameters on the energy and moisture transfer in the porous insulation are investigated. The present study, which for the first time presents a full simulation of the problem considering the mobile condensate, can be applied to other classes of problems on heat and mass transfer with phase change through a porous medium. iii Table of Contents Abstract List of Tables^ viii List of Figures^ ix Nomenclature^ xiii Acknowledgement^ 1 2 3 xvii INTRODUCTION 1 1.1 Background ^ 1 1.2 Physical Process Description ^ 2 1.3 Geometrical Characteristic of Fibrous Insulation ^ 5 1.4 Motivation of the Present Study ^ 5 A BRIEF REVIEW OF LITERATURE 8 2.1 Background ^ 8 2.2 Experimental Studies ^ 10 2.3 Theoretical Studies ^ 12 2.4 Scope and Objectives of the Present Study ^ 14 MATHEMATICAL DESCRIPTION 16 3.1 Introduction ^ 16 3.2 Continuum Approach ^ 16 iv 4 3.3 Local Volume Average Technique ^ 17 3.4 Basic Assumptions ^ 18 3.5 Governing Equations ^ 19 3.6 Constitutive Correlations ^ 20 QUASI-STEADY ANALYTICAL MODEL 22 4.1 Background ^ 22 4.2 Analysis ^ 22 4.3 Nondimensional Form ^ 24 4.4 Solution Procedure ^ 26 4.4.1^Dry Region ^ 26 4.4.2^Wet Region ^ 27 4.4.3^Wet-Dry Interface ^ 28 4.4.4^Trial and Error Procedure ^ 29 Heat Transfer ^ 30 4.5.1^Heat Flux at Cold Plate ^ 30 4.5.2^Effective Thermal Conductivity ^ 30 4.6 Round-Pipe Insulation ^ 31 4.7 Concluding Remarks ^ 33 4.5 5 ANALYTICAL RESULTS AND DISCUSSIONS ^ 34 5.1 Introduction ^ 34 5.2 Comparison with Experimental Results ^ 34 5.3 Comparison with Literature Reports ^ 36 5.4 Effects of Condensation on Heat Transfer ^ 38 5.5 Concluding Remarks ^ 44 v 6 TRANSIENT NUMERICAL MODEL ^ 46 6.1 Introduction ^ 46 6.2 Analysis and Formulation ^ 47 6.2.1 General Formulation ^ 6.2.2 Stage 1: Initial Process 47 ^50 6.2.3 Stage 2: Immobile Liquid Accumulation ^ 51 6.2.4 Stage 3: Mobile Liquid Diffusion ^ 53 6.2.5 Stage 4: Long Term Liquid Accumulation and Diffusion ^ 56 6.3 Solution Methodology ^ 56 6.3.1 Discretized Formulations ^ 57 6.3.2 The Solution Algorithm ^ 63 6.4 Numerical Considerations ^ 7 NUMERICAL RESULTS AND DISCUSSIONS ^ 67 69 7.1 Introduction ^ 69 7.2 Physical Data ^ 69 7.3 Diffusion and Condensation Processes ^ 70 7.3.1 Initial Process ^70 7.3.2 Quasi-steady State Period ^ 73 7.3.3 Liquid Accumulation and Diffusion ^ 73 7.4 Comparison with Experimental Results ^ 75 7.4.1 Temperature and Heat Transfer ^ 76 7.4.2 The Moisture Gain and Liquid Transport ^ 76 7.5 Thermal and Transport Performances ^ 85 7.5.1 The Variation of Heat Flux ^ 85 7.5.2 Effect of Humidity Levels ^ 87 vi 7.5.3 Effect of Convective Heat Transfer Boundary Conditions ^ 88 7.5.4 Effect of Slab Thickness ^ 7.5.5 Effect of Porosity 88 ^93 7.5.6 Comparison of Different Thermal Conductivity Models ^ 93 8 CONCLUSIONS AND RECOMMENDATIONS^ 95 8.1 Analytical Study ^ 95 8.2 Numerical Study ^ 96 8.3 Recommendations ^ 97 Bibliography^ 98 Appendices^ 102 A Water Vapor-Air Mixture Diffusion^ 102 A.1 Vapor Diffusion Coefficient in Porous Media ^ 102 A.2 Mass Transfer Coefficient in Ambient ^ 103 A.3 Calculation of Vapor Concentration in Ambient ^ 104 B Thermal Conductivity of Moist Insulation ^ vii 105 List of Tables 5.1 Comparison of Predicted and Measured [23] Heat Flux ^ 35 5.2 Comparison of Predicted and Measured [21,22] Effective Thermal Conductivity ^38 7.1 Physical Data. ^ viii 70 List of Figures 1.1 Schematic Diagram of Physical Problem ^3 1.2 Four Stages in Energy and Moisture Transport Process ^4 3.1 Macroscopic Average Volume Element ^ 18 4.1 Schematic Diagram of Round-Pipe Insulation. ^ 32 5.1 Temperature distribution in slab; Comparison with experimental Data [23]. 35 5.2 Vapor Flux Distribution in slab, Comparison with experimental data [23]. 36 5.3 Comparison of analytical and numerical results, Conditions: L = 0.099 m, T. = 20°C, kdry = 0.037 W/m K, h* = 12 W/m2 K ^ 37 5.4 Variation of heat flux with temperature difference, T. = 40°C. ^ 39 5.5 Variation of heat flux with temperature difference, T. = 30°C. ^ 39 5.6 Variation of heat flux with temperature difference, Ta = 20°C. ^ 40 5.7 Variation of effective thermal conductivity ratio with temperature difference, T. = 40°C ^ 40 5.8 Variation of effective thermal conductivity ratio with temperature difference, T. = 30°C ^ 41 5.9 Variation of effective thermal conductivity ratio with temperature difference, Ta = 20°C ^ 42 5.10 Effect of h on the variation of the effective thermal conductivity ratio. . ^ 43 5.11 Effect of K on the variation of the effective thermal conductivity ratio. ^ 43 5.12 Effect of e on the variation of the effective thermal conductivity ratio. . ^ 44 ix 6.1 Grid-point Cluster for One Dimensional Transport Problem: (a). Control Volume for the Internal Points, (b). Control Volume for the Boundary 57 Points ^ 6.2 Grid-point Cluster for the Interface Location of Wet-dry Regions ^ 61 7.1 Temperature Distribution in the Initial Stage; p* = 53 kg/m3, L = 66 mm, RH . 96%, h* = 12 W/m2 K ^ 71 7.2 The Distribution of Vapor Density in the Initial Stage; p* = 53 kg/m3, L . 66 mm, RH = 96%, h* = 12 W/m2 K ^ 71 7.3 The Time Variation of Condensation Rate in the Initial Stage; p* = 53 kg/m', L = 66 mm, RH = 96%, h* = 12 W/m2 K. ^ 72 7.4 The Time Variation of Volumetric Liquid Fraction in the Initial Stage; p* = 53 kg/m3, L = 66 mm, RH = 96%, h* = 12 W/m2 K ^ 72 7.5 Temperature Distribution in the Quasi-steady State Period, Comparison of Numerical, Analytical and Experimental Results [22]; p* = 53 kg/m2, L = 66 mm, RH = 96%, h* = 12 W/m2 K ^ 74 7.6 The Time Variation of Liquid Fraction over a Long Term Period; p* = 53 kg/m2, L = 66 mm, RH = 96%, 17,* = 12 W/m2 K, 71, = 33°C , T2 = 6.8°C. 74 7.7 The Temperature Distributions in Run 1, Comparison of Numerical Results with Experiment Data [24]; p* = 53 kg/m2, L = 62.02 mm, RH = ^97%, h* = 12 W/m2 K 77 7.8 The Temperature Distributions in Run 2, Comparison of Numerical Results with Experiment Data [24]; p* = 53 kg/m', L = 68.65 mm, RH = ^96.5%, h* = 12 W/m2 K x 78 7.9 The Temperature Distributions in Run 3, Comparison of Numerical Results with Experiment Data [24]; p* = 53 kg/m', L = 70 mm, RH = 96%, h* = 12 W/m2 K ^ 79 7.10 The Temperature Distributions in Run 5, Comparison of Numerical Results with Experiment Data [24]; p* = 53 kg/m3, L = 76.13 mm, RH = 96%, h* = 12 W/m2 K 80 7.11 The Heat Flux at the Cold Plate, Comparison of Numerical Results with Experiment Data for Four Runs [24]; p* = 53 kg/m3, h* = 12 W/m2 K. . 81 7.12 The Total Moisture Gain, Comparison of Numerical Results with Experiment Data for Four Runs [24]; p* = 53 kg/m3, h* = 12 W/m2 K ^ 83 7.13 The Liquid Distribution in Slab at Different Times, Comparison of Numerical Results with Experiment Data [24]; Run 1. p* = 53 kg/m3, h* = 12 W/m2 K, L = 62.02 mm ^ 84 7.14 The Variation of Heat Flux at Cold Boundary with Times; p* = 53 kg/m', h* = 12 W/m2 K, Ta = 33°C, T2 = 6.8°C, RH = 90%. ^ 86 7.15 The Variation of Equivalent Thermal Conductivity along the Insulation Slab; p* = 53 kg/m3, h* = 12 W/m2 K, T = 33°C, T2 = 6.8°C, RH = 90%. 86 7.16 The Variation of keff/kdry at Cold Boundary and Total Liquid Concentration under Different Humidity Levels; p* = 53 kg/m3 ^ 89 7.17 The Variation of Iceffikh.v at Cold Plate and Total Liquid Concentration under Different Convective Heat Transfer Boundary Conditions; p* = 53 kg/m', T. = 33°C, T2 = 6.8°C, L = 66 mm, RH = 90%. ^ 90 7.18 The Variation of Iceff/kdrv at Cold Boundary and Total Liquid Concentration under Different Slab Thickness; p* = 53 kg/m3, T. = 33°C, T2 = 6.8°C, h* = 12 W/m2 K, RH = 90% xi 91 7.19 The Variation of Iceff/kdry at Cold Boundary and Total Liquid Concentra- tion under Different Porousities; p* = 53 kg/m', T. = 33°C, L = 66 mm, h* = 12 W/m2 K, RH = 90% 7.20 The Variation of keff/kdry 53 kg/m3, T. = 33°C, T2 = 6.8°C, 92 at Cold Surface Using Different k Models; p* = T2 = 6.8°C, L = 66 mm, h* = 12W/m2 K, RH = 90%. 94 xii Nomenclature a^constant in Clausius-Clapeyron relation, Eq. (4.26) a^coefficient in linear discretized equations b^constant in Clausius-Clapeyron relation, Eq. (4.26) bi^constant in Equation (4.16) b2^constant in Equation (4.16) B Biot number, h*LIk*ff B,„^mass transfer Biot number, km * L/Dff cp^dimensionless heat capacity at constant pressure C(z)^dimensionless concentration of vapour C.^dimensionless concentration of vapour in ambient Cc^dimensionless concentration of vapour at dry-wet interface Dt,^dimensionless vapor diffusion coefficient Di^dimensionless liquid diffusion coefficient defined in Eq. (6.26) Dfi^dimensionless liquid diffusion coefficient defined in Eq. (3.4) D12^binary diffusion coefficient of Fick's law [m/s2] f(9)^dimensionless function for concentration l(0)^df/d8 Fo t. Fourier number, L2/ao.,eff g(T)^Clausius - Clapeyron relation [kg/kg(dry air)] G dimensionless condensation rate defined in analytical study G*^dimensionless condensation rate defined in numerical study h^dimensionless heat transfer coefficient in ambient h19^ enthalpy of vaporization [kJ/kg] hm^ dimensionless mass transfer coefficient in ambient dimensionless thermal conductivity kdr y^ dry thermal conductivity [W/m ke f f^ effective thermal conductivity defined in Eq. (4.29) [W/m K] keg^ equivalent thermal conductivity defined in Eq. (7.1) [W/m K] K(w)^hydaulic conductivity length of slab [m] Le^ Lewis number, a* ID: Leg^ equivalent slab length [m] m(x)^vapor concentration [kg/kg (dry air)] ma^vapor concentration in ambient [kg/kg (dry air)] m2^saturation concentration at temperature Pa^ T2 [kg/kg(dry air)] dimensionless partial pressure of air Pal^dimensionless total pressure of gas phases Pv^ dimensionless partial pressure of vapor — P5^dimensionless parameters defined in Eq. (6.10) Ra^ air gas constant [N m/kg K] vapor gas constant [N m/kg K] dimensionless heat flux at cold surface 4^heat flux at cold surface [W/m2} dimensionless heat flux Q'^dimensionless heat flux in Reference [38] radius of pipe insulation [m] To^outer radius of pipe insulation [m] 3dv ri^inner radius of pipe insulation [m] T(x)^temperature distribution [K] To^reference temperature in Clausius-Clapeyron relation [K] Ta^ ambient temperature [K] T2^ temperature of cold surface [K] S^liquid saturation ^ volume of porous insulation [m3 V„,^main stream velocity [m/s] w^weight liquid content x^distance along slab [m] y^transformed variable [m] z^dimensionless distance zd^dimensionless dry length a0,eff *^effective thermal diffusivity, kcjil4C0 P^dimensionless group E^ volumetric fraction €o^porosity (void fraction) T^vapour flux [kg/m2 s] (1)^dimensionless vapour flux Pa^ density of air [kg/m3] I'^condensation rate per unit volume [kg/m3 s] T^ tortuosity factor B(z)^dimensionless temperature distribution 19e^dimensionless temperature at the dry-wet interface XV Subscripts 0^initial value 2^value at cold surface a^ambient a^air phase dry east side control volume face eff^effective east side grid point i-th phase central grid point ref^reference vapor phase west east side control volume face west side grid point liquid phase 7^gas phase o-^solid phase Superscripts dimensional value value in previous time step xvi Acknowledgement I wish to express sincere gratitude to Dr. M. Iqbal and Dr. E. G. Hauptmann for their supervision and continuous encouragement through all stages of this work. Financial assistance from the National Science and Engineering Research Council of Canada is gratefully appreciated. My thanks go to Dr. Wijeysundera from The National University of Singapore who generously offered the experimental data. His valuable suggestions, comments, and discussions during his sabbatical leave at the Department of Mechanical Engineering, UBC were also very helpful. This thesis is dedicated to my wife, Wei-Hung, her understanding and support made the completion of this work possible. xvii Chapter 1 INTRODUCTION 1.1 Background Heat and mass transfer in porous media accompanied by phase change is a phenomenon which occurs frequently in nature and in many engineering applications. These applications are in the area of chemical , environmental, mechanical and petroleum engineering, geology and others. In recent years, heat and moisture transfer through fibrous insulation has drawn considerable attention among many investigators because of applications in the energy management of buildings, and heated and refrigerated building envelopes. Fibrous insulation is used within building sections, such as partition walls, ceilings, roofs, floors etc, to separate the heated space from the unheated space. In air conditioning systems, fibrous insulation is commonly used to insulate the chilled water pipes and cold air ducts. In these applications, vapor barriers are applied to prevent water vapor ingress from the ambient. However, because of factors such as poor installation, cracks develop with age, particularly near joints, clamps and supports. The water vapor will migrate into the insulation. Water vapor condensation can take place anywhere in a porous insulation when the vapor density is greater than the saturation vapor density which corresponds to the local temperature at that point. The condensation effect is especially prominent when the insulating material is exposed to large temperature differences and high humidity environments. As condensation occurs, the liquid phase resulting from condensation will cause a significant increase in heat transfer across the insulation and 1 Chapter 1. INTRODUCTION^ 2 hence it affects the thermal performance of the insulation. The long term effects of condensation are corrosion of the metallic components of the insulated system and the eventual destruction of insulation. 1.2 Physical Process Description The heat and moisture transfer in porous insulation is generally a complex multidimensional problem including vapor transport by diffusion and convection, flow of liquid due to gravity and capillary action, and condensation or freezing accompanied by release or consumption of latent heat. The solution of this entire problem is not attempted here. Instead, the problem is reduced to a one-dimensional configuration which still retains most of the important parameters of the original problem. The one-dimensional configuration investigated here is shown in Figure 1.1 An initial dry fibrous insulation has two boundaries which encounter two different environments: one of the boundary faces a cold, impermeable plate with temperature T2, the other boundary is exposed to a warm humid ambient air, with temperature T., humidity RH, vapor concentration C., convective heat transfer coefficient h., and mass transfer coefficient hm. The water vapor will migrate into the insulation slab due to the vapor concentration gradient and the thermal gradient, and it is expected to exhibit four stages of transport processes as shown in Figure 1.2. The first stage is a relatively short initial transient stage in which the temperature and vapor concentration fields are developing within the insulation slab. During this process, a very small quantity of liquid water accumulates in the porous medium. Condensation is defined here as the accumulation of liquid beyond the adsorption process in which the vapor reaches the saturation concentration at a given temperature. The region in which Chapter 1. INTRODUCTION^ 3 Warm, Humid Ambient Air Ta, Ma, h, hm Figure 1.1: Schematic Diagram of Physical Problem. the air-vapor mixture is saturated is referred to as the wet region where condensation takes place. The wet region is established , and the wet and dry region is separated by the wet-dry interface after this stage. In the second stage, the heat and vapor transfer processes reach a quasi-steady state, and the temperature and vapor concentration fields are invariable with time. Liquid is accumulated in the wet region, however the amount of liquid content is still low and does not have a significant effect on transport properties. When liquid accumulation exceeds a value high enough as to affect the properties significantly, the liquid starts to move due to generated liquid pressure and will flow towards the wet-dry interface and towards the drier, warmer surface by capillary action. The wet-dry interface is moving as the wet region expands due to the liquid outflow. In the last stage, the wet-dry interface eventually reaches the exposed surface, The Chapter 1. INTRODUCTION 4.1 ^ //// / ./ ..••^ It- Condensate Dry Insulation Ta + +^ 4 T2 ^ P.Wet-dry Interface Ta ++44+4+444 ^ (b) Quasi-steady Stage: (a) Initial Stage: Condensation ^ at Cold Plate. Liquid Accumulation in the Wet Region. T2 Dry Region Ta # # # # # # # ^ (d) Long Term Stage: (c) Liquid Diffusion Stage: ^ Liquid Accumulation and Liquid Diffusion Starting ^ Diffusion in Entire Slab. from Cold Surface. Dry Region Wet Region Liquid Diffusion Region Figure 1.2: Four Stages in Energy and Moisture Transport Process. Chapter 1. INTRODUCTION^ 5 effects of liquid accumulation and liquid flow become dominant. The accumulation of liquid in the insulation is believed to have two important effects: (1). decrease in the local vapor diffusion coefficient, (2) increase in local thermal conductivity of the slab. These factors will both influence the heat transfer and condensation rate in the slab. 1.3 Geometrical Characteristic of Fibrous Insulation The Fiberglass batts used as thermal insulation consist of layers of Fiberglass filaments that are held together with a phenolic binder. The batts usually have a very high void fraction (95-99 percent). The Fiberglass filaments have varying diameters with an average value of 3.56(10)-6m. Each Fiberglass layer differs from its adjacent layers and has a different void fraction. The Fiberglass filaments are dispersed randomly in each layer. However, the number of fibers along the batt is slightly larger than that across the batt, and the number of fibers that are perpendicular to the batt and run from one layer to the other is very small. Generally speaking, fibrous insulation is an inhomogeneous assemblage of Fiberglass filaments and phenolic binder. The insulation is characterized by an anisotropic fiber density and a nonuniform void fraction. 1.4 Motivation of the Present Study A detailed study for the entire transport process in fibrous insulation is complicated by the following factors: • Heat and mass transport in porous insulation is a problem with multiphase flow. In general a wet porous insulation consists of three phases: the solid matrix, the liquid water, and a binary gas phase composed of air and water vapor. Chapter 1. INTRODUCTION^ 6 • The transport mechanisms involved in the process are quite complicated. Energy transport in such a medium occurs by conduction in all of the phases as well as by convection with those phases which are able to move. In addition, there is heat transfer caused by phase change at the interface between liquid and gas. Mass transfer occurs within the voids of the medium. In an unsaturated state, these voids are partially filled with a liquid, whereas the rest of the voids contains air and water vapor. In the gas phase, there is vapor diffusion due to the vapor concentration gradients, bulk convection due to the density variation induced by temperature gradients, and air infiltration due to the small difference in gas pressure across the insulation. The flow of liquid is caused by internal forces, such as capillary, intermolecular and osmotic forces, and external forces, such as imposed pressure difference, and gravity. • The transport properties involved vary strongly with structure of the porous medium, moisture content and temperature. • It is almost impossible to study quantitatively the transport processes for the irregular void configurations which exist in general in porous insulation. It is still difficult even for a regularly shaped matrix. Thus there is indeed a need for a rigorous and extensive investigation of the heat and mass transport processes in porous insulation. This may lead to better estimation of insulation properties for building design purposes. In the present work, the focus is on the analysis of the process of moisture accumulation and transport in the fore-mentioned insulation system, and as well as the effect of moisture on thermal performance of insulation. The entire process of energy and mass transport in the porous insulation system is quantified and the significant transport mechanisms are identified. A simplified formulation for four stages of the transport process is Chapter 1. INTRODUCTION^ 7 presented to simulate the transient and spatial variation of the pertinent variables. An analytical model is also presented to perform a parametric study on the thermal performance of porous insulation in the presence of condensation. Both the analytical and transient numerical model are validated by comparing with experimental data. Finally the interesting effects of variation of pertinent parameters on the energy and moisture transfer are investigated. Chapter 2 A BRIEF REVIEW OF LITERATURE 2.1 Background The importance of heat and mass transport in a porous medium is well appreciated. The early treatments of liquid flow and heat transfer in a porous medium were more empirical rather than rigorously theoretical. The initial macroscopic treatments on water flow through a porous medium were made in the 1880's [2]. The porous medium, which is a heterogeneous system made of a solid matrix with its void filled with gases and liquids, can be treated as a continuum by properly accounting for the role of transport through each phase in this system of phases. The attempt at microscopic (pore-level) transport studies began in the early 1900's [2]. The pore-level analysis offered a better understanding of heat and mass transport, however the description and solution of a transport problem at the microscopic level is impractical and, perhaps impossible because the geometry of the surface that bounds the phase is not observable and is too complex to be described. The macroscopic level approach, at which continuous and differentiable quantities may be determined, is needed to model and predict the heat and mass transport in a porous medium. Progress was achieved in the second half of this century with the more rigorous approaches of local volume averaging techniques [3]. Early studies on transport in porous insulation were mainly on the motion of liquid through unsaturated porous media. In early papers, heat transfer was not given sufficient attention until Krischer (1940) [4]. Probably one of the most significant contributors to 8 Chapter 2. A BRIEF REVIEW OF LITERATURE^ 9 heat and mass transfer in porous media, Krischer first considered seriously the intimate role that transport of energy may play in moisture transfer in a porous medium. Prior to Philip and DeVries (1957) [5], moisture movement was explained by a simple diffusion theory, the flux being given by Fides law. Philip and DeVries (1957) [5] and DeVries (1958) [6] extented the theory and generalized the effects of capillary flow, vapor transport and energy transfer, and represented the transport due to capillary forces in terms of gradients of moisture content and temperature. Their work consolidated most of the previous knowledge on the influence of temperature gradients on moisture movement and laid the groundwork for a better understanding of simultaneous heat and water movement in soil. Extensive research on heat and mass transfer in porous media has been done by the Soviet scientist, Luikov (1966, 1975)[7, 8]. He described heat and mass transfer as caused both by the temperature and the liquid concentration gradient, and developed a mechanistic approach for the representation of simultaneous heat and moisture transfer in dying processes. A set of linear mass and heat transport equations were derived. However, the effect of phase change was not considered and the dependence of transport coefficients on temperature and concentration were neglected. Whitaker (1977) [3] presented general formulations of simultaneous heat, mass and momentum transfer in porous media. The intrinsic phase averaged quantities were defined at the appropriate level and the continuum point equations which are averaged were described. Although simultaneous heat and moisture transfer in a porous medium has been studied by many other researchers [9, 10, 11], most of the work has been directed towards studies of soils and sands. Very little attention has been given to fibrous insulating materials until the last fifteen years. Chapter 2. A BRIEF REVIEW OF LITERATURE^ 10 2.2 Experimental Studies The experimental investigation of porous insulation has focused mainly on two aspects: 1). the effect of moisture gain by porous insulation on thermal conductivity of and heat transfer through the insulation; 2). transport of liquid condensate within insulation. Jespersen (1953) [12] measured the thermal conductivity of a 62 kg/m3 glass fiber, using a steady heat flow method, but with a very small gradient of temperature across the sample in order to get a uniform moisture distribution. Joy (1957) [13] worked on 96 kg/m3 fibrous insulation, using a transient heat flow method. Their results show a significant increase of thermal conductivity even for a very low moisture contents. Similar measurements to Jespersen and Joy's work have been performed by Langlais (1982) [14], using a classic steady-state method. The thermal conductivity of a high density mineral fiberboard was determined as a function of moisture content. However, their results show that thermal conductivity is far less affected for low moisture content. They explained that their differences with Joy's results are due to the non-uniform moisture distribution under the influence of the temperature gradient. Bomberg and Shirtliffe (1978)[15] also demonstrated the influence of this redistribution on thermal conductivity. The moisture redistribution process was also studied experimentally by Kumaran (1987, 1988) [16, 17] in an insulation slab under the influence of a temperature gradient. Langlas et al. (1983) [18] and Langlais and Klarsfeld (1984) [19] further studied the effects of moisture transfer through an insulation slab under two sets of boundary conditions and pointed out that the rapid variations in thermal conductivity values were due to moisture phase change and diffusion. Thomas et al. (1983) [20] measured the moisture and temperature distribution , as well as the thermal conductivity versus time. Modi and Benner (1985, 1986)[21, 22] measured the moisture gain of spray-applied fiber-glass and cellulose insulation slabs from the surroundings and also the effect of the moisture Chapter 2. A BRIEF REVIEW OF LITERATURE^ 11 gain on the thermal conductivity. Wijeysundera et al. (1989) [23, 24] conducted two series of experiments to study the diffusion and condensation of water vapor in fiber-glass. The first set of experiments were tests on fibrous insulation. In the second series of experiments, one of the faces of the insulation slab was exposed to a warm and humid ambient and the other one was maintained at a low temperature. The temperature and liquid distribution, the total moisture gain, the heat flux were measured for a range of experimental conditions. Further experimental studies on effects of moisture gain by a fiber glass insulation slab in a second set experimental conditions were reported by Wijeysundera (1992) [25]. The major heat and moisture transfer parameters were measured over a period of up to 600 hours. The transport of the liquid condensate within fibrous insulation has received less attention. Motakef and El-Masri (1985) [26] presented a model for the isothermal liquid diffusion in fibrous insulation in the absence of gravitational forces. The model related the liquid diffusion to the characteristics of the insulation, such as void fraction, spatial distribution, directionally index, tortuousity factor etc. Experimental data was also reported to verify the model. Timusk and Tenende (1988) [27] studied the capillary rise of water in fiber glass slabs. Cid and Crausse (1990) [28] measured the coefficient of liquid diffusion which was correlated to moisture content and the nominal density of the solid phase, based on the principle of attenuation of two radioactive emissions. Although the thermal effects were not added to the measurements, so doing would makes it possible to predict the behavior of the materials when they are submitted to moisture. Extensive field studies of heat transfer through wet insulation have been reported by Heldlin (1988) [29], during all seasons of the year. Chapter 2. A BRIEF REVIEW OF LITERATURE^ 12 2.3 Theoretical Studies As mentioned in the beginning of this Chapter, simultaneous heat and mass transfer has been extensively studied for various systems [3-11]. However, these studies have been only recently extended to heat and moisture transfer in insulation . Dinulescu and Eckert (1980) [30] analyzed the moisture migration in a slab of an unsaturated porous medium between two impermeable surfaces. The problem of condensation in insulation was first studied by Ogniewicz and Tien (1980) [31]. The condensation process in insulation is characterized in terms of three regimes, and a quasi-steady model was presented to studied the condensation effects in the second regime, where coupling between temperature and concentration of the condensing vapor was taken into account. Motakef and El-Masri (1986) [32] presented a one-dimensional analytical model for heat and mass transport with phase change in a porous slab similar to Oganiewcz and Tien [31], and the closed-form approximate solutions for two limiting regimes of condensate diffusivity were reported. The same problem was also modeled by Vafai and Sarkar(1986) [33] as a transient, multiphase flow, with variable properties. The formulation was based on the local volume-averaging technique for each phase to come up with the governing equations for the condensation process in fibrous insulation. However, the simplifying assumption that liquid accumulation is small was made in arriving at the solution. Shapiro and Motekef (1989) [34] extended their analysis to unsteady transport processes and reduced the unsteady process to that of quasi-steady fields in time-dependent domains corresponding to mobile and immobile condensates. The model underpredicts the energy transfer in the wet zone and the range of validity of the solutions is limited, especially for the mobile condensate. Chapter 2. A BRIEF REVIEW OF LITERATURE ^ 13 Vafai and Whitaker (1986) [35] first reported a two-dimensional unsteady numerical work with simplifying assumptions. Vafai and Tien (1988, 1989) [36, 37] performed a more thorough numerical simulation of the two-dimensional multiphase transport process in a porous medium. Information was presented on the variation and the intercoupling effects of the important field variables, the effects of some parameters such as humidity levels, the porosity etc. on the condensation rate, liquid accumulation and energy transfer. However, No attempt was made to verify numerical solutions by comparison with experimental data. The computational results presented in the paper only show the transport process within the first 500 seconds. In such a short initial period, the absorption process of bounded liquid is still in the initial stage, the effects of temperature gradient and condensation are not significant, and the contribution to the heat and mass transfer in insulation is small considering the long term impact of accumulation and transport of the liquid condensate through the insulation. In addition to, the treatment used in their scheme to determine when condensation should be taken into account lacks either experimental validity or theoretical support. Wijeysundera et al. [23, 25] used a semi-empirical model to interpret some of their experimental results. The analytical model developed for transport process during the quasi-steady phase has given good predictions of heat and vapor flow through insulation for low level of moisture transfer. Tao et al. (1991) [38] analyzed the moisture and frost accumulation in a glass-fiber slab. An approach similar to the study of Vafai and Sarker [32], who used the local volume average technique, was applied to the case with temperature below the triple point of water. A numerical simulation was performed for one-dimensional, transient, vapor diffusion with phase changes and variable properties. Chapter 2. A BRIEF REVIEW OF LITERATURE^ 14 2.4 Scope and Objectives of the Present Study It is clear from the literature survey that the transfer processes in insulation can be quantified by a). field and laboratory studies, b). analytical/numerical studies. Generally speaking, the field and laboratory studies are time consuming due to the slow process of moisture transfer. Furthermore, it is difficult to cover the effects of such a wide range of operating and design parameters. The available analytical studies mainly dealt with the quasi-steady phase of the transfer processes. Although several investigators have reported good predictions with their models, the performance of insulation in the presence of condensation has not been sufficiently analyzed. From a design stand-point, it is important to study the effects of moisture on thermal performance of insulation under practical operating conditions. Less work has been done to simulate numerically the full transport processes which the insulation may undergo for a long time period in the presence of condensation. The long term effects of condensation on the performance of insulation has not been studied theoretically. The present study focus on the analytical and numerical investigation of the transfer processes through fibrous insulation in the presence of condensation. The objectives of this study are 1). To perform a detailed and parametric study on effects of condensation on the thermal performance of insulation in the quasi-steady phase. The analytical models based on the quasi-steady models by Wijeysundera et al. [22] will be extended to cover the main design and operating parameters on which the effective thermal conductivity depends, with the hope that the modified model will be able to to predict the most important heat and mass transfer performance for engineering design and application purposes. Chapter 2. A BRIEF REVIEW OF LITERATURE^ 15 2). To perform a numerical simulation study that covers all the phases mentioned in Chapter one, and also uses measured physical property variations available in the literature. The numerical results will be compared with the experimental data by Wijeysundera et al. {22,24 Accordingly, the analytical and numerical models will be derived, and the computer codes developed. Chapter 3 MATHEMATICAL DESCRIPTION 3.1 Introduction This chapter presents the background for the general mathematical models that describe transport of heat and mass in porous insulation at a macroscopic level. The general formulation of this problem is based on the local volume-averaging technique [3] for mass, momentum, and energy equations for each phase which make up the governing equations for heat and mass transport process in insulation. The general governing equations for the transport of energy, water vapor and liquid, together with constitutive relations are described. The limitations and restrictions of the governing equations are stated and the validity of the assumptions are discussed. 3.2 Continuum Approach Porous insulation is a multiphase system, where the various phases are separated from each other by an abrupt interface. It is difficult to formulate the transport problem including the interface boundaries between each phase. The conventional treatment is to replace the real system, where the multiphases together occupy disjoint subdomains within a porous medium, with a model in which each phase is assumed to behave as a continuum over the whole domain. With the continuum approach, the average values of phase variable are taken over elementary volumes, centered at the point. By assigning average values to every point, 16 17 Chapter 3. MATHEMATICAL DESCRIPTION^ the variables which are differentiable functions of the space coordinates can be obtained. The continuum approach is thought reasonable since it circumvents the needs to specify the exact configuration of the interface boundaries, and it also is useful in solving field problems of practical interest due to its measurable quantities. The information regarding interface boundaries and the actual variation of quantities within each phase still remain in the form of coefficients. The numerical values of these coefficients must be determined experimentally for a specific porous media, in the laboratory, or in the field. 3.3 Local Volume Average Technique The local volume technique is applied by associating with every point in the porous medium an averaging volume V that is bounded by a closed spatial surface A, as shown in Fig. 3.1. In general, the volume V is composed of three phases. These are: the solid phase V,, the liquid phase Vfi(t), and the gas phase V.,f(t). At any location in space z, a quantity y is said to be spatially averaged when it is defined as 1 < y > (z) = v(z) fv(z)ydV. (3.1) A quantity in phase a is said to be intrinsic phase averaged when it is defined as 1 <y >a= — v (t) iva(t)yadV,^ (3.2) where ya is zero in phases other than the a phase, and <y >a is the mean value of y in Va. In this study, the averaging symbols are omitted in order to simplify the notation. The derivation of the general governing equations for heat and mass transfer in a porous insulation is based on Whitaker's work [3]. Several assumptions are made to arrive at the governing equations. Chapter 3. MATHEMATICAL DESCRIPTION^ 18 Liquid Phase (13) Figure 3.1: Macroscopic Average Volume Element. 3.4 Basic Assumptions (1). The fibrous insulation is homogeneous and isotropic. This assumption is the common simplying procedure to rationally tackle the problems for heat and mass transfer in porous materials. The structure of the porous media is difficult to characterize. Motakef et al. [26] has been able to correlate diffusion of liquid using the void fraction, pore size distribution, tortuosity and so on. It seems clear that describing the parameters characterizing the structure will be a most difficult task for porous insulation, therefore the present theoretical development is restricted to the case where the structure parameters are independent of the spatial coordinates. (2). The solid-liquid-gas system is in local thermal equilibrium. Since the heat and mass transport process investigated here, as well as most of drying processes, are characterized by no or relatively low convective transport rates, under these circumstances one could therefore assume that conductive transport is sufficient to eliminate significant temperature differences between the separated phases. Chapter 3. MATHEMATICAL DESCRIPTION^ (3). 19 No convective gas phase flow occurs in the insulation matrix; any moisture accumulation is caused by vapor diffusion only. In the present physical model, the effects of air infiltration are not included. This restriction might be removed in further studies. (4). The total gas phase pressure in the insulation matrix is constant. This assumption is usually justified for drying process and building envelopes. 5). The gravitational effect on the liquid motion is negligible. Under the assumptions (3) and (4), the gas and liquid pressure distribution can be considered to be hydrostatic. Thus the forces exerted on the liquid phase are of two kinds: surface tension and gravity forces. For porous insulation with very small pores, the effect of gravity can be neglected due to the large surface tension forces. The experimental observation by Wijeysundera et al. [25] showed that after long term testing up to more than 300 hours, no liquid was observed to drip from a totally wet insulation slab. 3.5 Governing Equations The derivation of governing equations requires considerable algebraic manipulation. As a result, they are taken directly from Whitaker's work [3]. The governing equations for water vapor, liquid, and energy transport are given as the following: Vapor diffusion equation: a^op:r _ a(c-r4) -[D:;(11)-a-;1 - - at* (3.3) Liquid transport equation: a^*^act;*aeg, PA-ai[DA(w)Til+ r PA at** Energy equation: (3.4) ^ ^p: Chapter 3. MATHEMATICAL DESCRIPTION ^ a (371 aT — [e(w)— j + hfgr = p`c* . Ox^ax^P at. 20 (3.5) 3.6 Constitutive Correlations Volumetric Constrain: The volume fractions for the three phases are defined as e„ = 11,7711, ^Eg = vo ( t ) / v,^el, = 14/V(t).^(3.6) Therefore the sum of these fraction is one, fa-Fe9+ey=1.^ (3.7) Thermodynamic Relations: Some thermodynamic relations are needed to connect the thermal energy equation with the vapor diffusion equation. By treating the gas phase as ideal gas, the air and vapor density can be determined by = Rep:T,^ P: = Rv P:71 . ^ (3.8) (3.9) The vapor pressure-temperature relation for the vaporizing gas can be represented by the Clausius-Clapeyron equation: h191^1 P: = P:oexP[—(— Re T —^)iTref (3.10) Chapter 3. MATHEMATICAL DESCRIPTION^ 21 The limitation of the Clausius-Clapeyron relation is that the gas-liquid interface is assumed flat, that is surface tension effects are not significant. When the effects of curvature and surface tension are important, the Kelvin equation may be used to represent the vapor pressure-temperature relation [45]. Assumption (3) gives the other relation for gas phase pressures: P: = Ptotal P:. (3.11) The above three governing equations and five correlations can be used to solve the following eight unknown: pa, pa, egg, c, T, r, p„, pa. The effort required to solve the coupled transport equation is significant, although the numerical methods to solve these equations are well known. What appears to be extremely difficult is that special models need to be constructed so that the various transport stages in the problem can be studied properly, and comparison of the theoretical model with experiments in order to determine the parameters that appear in the transport equations. The general mathematical description of transport process needs to be extended to a series of specified theories that adequately describe the whole heat and mass transport process. Also we need some theoretical and experimental basis for estimating the values of the parameters which appear in the transport equations. Chapter 4 QUASI-STEADY ANALYTICAL MODEL 4.1 Background Experimental studies [23, 24] have shown that after a short initial transient period, the energy and vapor transfer processes reach a quasi-steady state. The numerical simulations by Tao et al. [38] have shown that quasi-steady conditions are reached within about one hour for practical situations. During this quasi-steady phase, heat flux and therefore the effective thermal conductivity attain high values. In this section a detailed analysis is developed for energy transfer during the quasisteady phase. Expressions for heat flux and effective thermal conductivity are derived for the flat-slab system. These are easily extended to a round-pipe system by a simple transformation of coordinates. 4.2 Analysis A schematic diagram of the slab is shown in Fig. 1.1. The outer surface of the insulation is exposed to warm humid ambient air from which heat and water vapor flow towards the cold surface. These flows are caused by the temperature and vapor pressure gradients within the insulation. When the local vapor concentration at a point exceeds the saturation concentration corresponding to the local temperature, condensation occurs. The wet region is established and separated from dry region by the wet-dry interface. The vapor that reaches the impermeable surface condenses at this surface. 22 Chapter 4. QUASI-STEADY ANALYTICAL MODEL ^ 23 In addition to the basic assumptions made in Chapter Three, the following assumptions are made in the analytical model. (1) The liquid phase is immobile and its effect on properties such as the vapor diffusion coefficient and the thermal conductivity is negligible. Since the quasi-steady state is reached in a short time (about an hour) for practical situations [38], the quantity of liquid condensed is relatively small, and it will therefore have little effect on the properties. (2) All properties are independent of temperature. This assumption is justified because the temperature range of the present study is limited to about 0-40°C. The properties of the slab are therefore, evaluated at the mean temperature. (3) The initial adsorption process is not taken into account because of its relatively small contribution to the total liquid fraction; only the liquid accumulation due to condensation in the wet-region is considered. Subject to the above assumptions, the governing equations of energy and mass transfer for the quasi-steady state may be written as: d2T le— + 'lig = 0, dx2 (4.1) and d2m dx2 In the dry region the condensation rate r = O. r. (4.2) 0 and Equations (4.1) and (4.2) are uncoupled. In the wet region the vapor concentration m is a function of the local temperature. The boundary conditions at the exposed surface, x = 0, of the slab are: d – k*— T = h* [T„ – T(0)], dx (4.3) Chapter 4. QUASI-STEADY ANALYTICAL MODEL^ 24 and dm - p:D:,— = h:n[rna — m(0)]. dx (4.4) At the impermeable surface, x = L, the boundary conditions are: T(L) = T2 and m2 = g(T2).^ (4.5) 4.3 Nondimensional Form In order to generalize the equations, the following nondimensional variables are introduced. Chapter 4. QUASI STEADY ANALYTICAL MODEL^ - = Dimensionless distance, z Dimensionless temperature, 9 Dimensionless concentration, _ a _ Heat transfer Biot number, B Mass transfer Biot number, X L ^T —T2 — Ta—T2 171- m 2 ma—m2 h* L = k* (4.6) Bm = 11L AID: Dimensionless condensation rate, 25 — nr, P:DZ (In a—m2)I L _ _AL_ — lc* (Ta—T2) Dimensionless heat flux, Dimensionless vapou flux, - Ratio of latent heat to conductive heat, # T.L p:D.(rna—m2) p:D:h f kna—m2) — lc* (Ta—T2) The meaning of the symbols is defined in the Nomenclature. The governing equations in terms of the dimensionless variables are: d2B + GO = 0, dz2 (4.7) and d2C dz2 +G=0. (4.8) Chapter 4. QUASI-STEADY ANALYTICAL MODEL ^ 26 The boundary conditions become: at z = 0, d9 — — = B[1— 0(0)j,^ dz (4.9) and 1 dC — — = B„,[1^C(0)]; — dz at z = 1, 0(1)^=^0, (4.10) (4.11) C(1)^=^0. 4.4 Solution Procedure 4.4.1 Dry Region The solution of Equations (4.7) and (4.8) in the dry-region are obtained by setting the condensation rate G = 0 and using the following conditions at the wet-dry interface. At Z = zd, ^0(zd) = Oc, 1^ (4.12) O(zd) = C. The expressions for the temperature and concentration in the dry region are as follows: 0(z) — —[1 — Oc]z + Pc/ B + zd] ^ (4.13) [zd + 1/13] and —41 — Cclz + [Cc/ Bm + zd] C(z) —^ + 1/13,7,] ^• ^[zd (4.14) Chapter 4. QUASI-STEADY ANALYTICAL MODEL ^ 27 4.4.2 Wet Region Eliminating the condensation rate G between Equations (4.7) and (4.8), the total energy equation can be obtained as: d20 ad2C dz2 + ''' dz2 = °' (4.16) where i3 is a nondimensional group which may be interpreted as the ratio of latent heat transfer due to vapor flow to that of sensible heat transfer. Since in the wet region the vapor is saturated, the vapor concentration is a unique function of local temperature, C = 1(0); the form of the function f(9) will be given later in this section. Therefore, Equation (4-15) can be directly integrated to obtain the temperature distribution in the wet region as 0 + fiC(0) = biz + b2,^ (4.16) where the concentration C is a function of the local temperature, bi and b2 are constants of integration. The following boundary and interface conditions are used to obtain b1 and b2. Z = zd, 0 = 0c, and C = f (0 c); 1 (4.17) z =1, 8 = 0, and C = 0. This gives the temperature distribution in the wet-region as: 0 + Of(0) = [ec + 0 f (9 c)1( 1— z 1 — Xcl ) (1 < z < zd).^(4.18) Chapter 4. QUASI-STEADY ANALYTICAL MODEL ^ 28 4.4.3 Wet - Dry Interface Both temperature distributions given by Equation (4.13) for the dry region and Equation (4.18) for the wet region contain two unknowns, Oc and zd. This requires a careful study of the wet-dry interface. At the dry-wet interface the concentration, C = AO), therefore we have: dC^de dz I.=^ dz I. • (4.19) Substituting in Equation (4.19) for^and ficz- from Equations (4.13) and (4.14), the following expression is obtained for the dry region length zd [1— f(0c)]/B —[f(0c)(1-0c)]/B,„, zd^[P(0,)(1 — 8) — (1 — 1(8))]^• (4.20) By applying the following interface condition: dO^de , (4.21) 6-1.; Itvet= (7E; Wry)^ this gives the equation [9.+ fif(9.)1 ^(1 —Or) (4.22) (1 — zd)[1 ger(Oc)]^zd + 1 / B^ • Therefore 1 — zd _ ^Oc fii(Oc) zd +1/B — [1 + f(Oc)](1 — Oci (4.23) Substituting for zd from Equation (4.20) in Equation (4.23), the following nonlinear equation is obtained for the interface temperature Oc. (1 --0c).= (1 +1/13) [1 —f(Oc)]^1^1^Oc-f-Af(Oc) ( (1 +1/B,n) P(Oc) -F B Brn 1+10P(9c) )( ].^ (4.24) Chapter 4. QUASI-STEADY ANALYTICAL MODEL^ 29 The function f(0) in Equation (4.24) is the temperature dependence of the saturation concentration. In dimensional form this is given by 1^1 g(T) = a exp[b( — )],^ (4.25) which is the Clausius - Clapeyron relation. The values of the constants a, b, and To for the temperature range 0 58°C are [24 a = 0.0232,^b = 5271.2,^and To = 302K. In nondimensional terms, Equation (4.25) has the form aexp[b(A ^ ^e(Ta-T2 )4-T2 )1^M2 f(19) = ma — 7712^ ma — 7712 (4.26) Note that at the impermeable surface, 8 = 0, and AO) = 0. An examination of Equation (4.24) shows that Bc depends on three nondimensional numbers B, Bfl, and fl. However, Equation (4.26) reveals that f(9c) depends on the operating conditions Ta, T2 and ma. When these variables are specified, the solution of Equation (4.24) for the unknown Oc is easily accomplished by a trial and error procedure. Once Oc is found, all the important heat transfer parameters can be obtained directly. 4.4.4 Trial and Error Procedure The trial and error procedure for solving the interface temperature 9c is as following: 1). Assuming an initial value for the dry-wet interface temperature 9, the local saturated vapor concentration and its first derivation are obtained from the ClausiusClapeyron relationship. These values are substituted into Equation (4.20) to obtain the dry region length zd. 2). Applying the boundary condition 9 = 83 at the cold plate, one can obtained the wet region length from Equation (4.18). Chapter 4. QUASI-STEADY ANALYTICAL MODEL^ 30 3). The total slab length L is the sum of dry and wet region lengths. The computed slab length L is compared with the actual slab length. 4). If the computed slab length does not match the actual slab length, the Oc value is adjusted and the procedure from step 1 to 3 is repeated until the computed slab length matches the actual slab length. 4.5 Heat Transfer 4.5.1 Heat Flux at Cold Plate The heat flux at the cold surface is given by dT^dm qL.—[k*--Fhigp:D: dx (4.27) It should be noted that the heat flux includes both sensible heat transfer and the latent heat released at the impermeable surface due to condensation. Using Equations (4.18) and (4.20), the dimensionless heat flux is obtained as: [19c + Pf(9.)][(1 — f(9c))— f'(9)(1 — 9c)] 41, = [(1 + 1/B,„),P(0,)(1 — Oc)— (1 + 1/B)(1 f(9c))). (4.28) — 4.5.2 Effective Thermal Conductivity The equivalent thermal conductivity of the slab in the presence of condensation is defined by the same expression that is used for a dry slab. This is given by _ 1^L 41,(1-17 +^f.) = (Ta — T2). In the nondimensional form (4.29) 31 Chapter 4. QUASI-STEADY ANALYTICAL MODEL^ 91{1+ ( kdry )1 = 1. B^ff ke Substituting for qL in Equation (4.30) gives the following expression for (4.30) (keff/kdry) in terms of Oc: [9c + flf(ec)][1— f(e9c)— f(8c)(1— 9c)} k^[Oc + 13 f (9c) — (1 + 11 Bm)jfi(e9c)(1 — Oc) — [Oc + 13 f (I) — (1 + 1/B)(1 — f(Oc))]• (4.31) Iceff^ 4.6 Round-Pipe Insulation For the round-pipe insulation system shown in Fig. 4.1 , moisture diffusion is assumed to occur uniformly across the circumference of the cross-section. Making the same assumptions as for the flat-slab, the governing equations of heat and vapor transfer can be written as: d dT — [rk*— ]-1- higrr = 0, dr^dr (4.32) dm dr " dr (4.33) and d — [rp* D* — ] — rr = 0. Consider the transformation: y = ro inPi,^ r (4.34) where 7.0 is the outer radius of the insulation. Applying the transformation (4.34), Equations (4.32) and (4.33) become (4.35) Chapter 4. QUASI-STEADY ANALYTICAL MODEL ^ 32 Wet-dry Interface Ta, ma, h, hm Figure 4.1: Schematic Diagram of Round-Pipe Insulation. and dm c-1-[*D* — ]- r(-7-7 . dp y " dy^ro .o (4.36) The boundary conditions at the outer surface in the transformed system are: at r = 7-0, i.e y=0 dT k* — = /eV'. — T(0)], dy (4.37) and p:D: dm dy = h:i[mc, — m(0)].^ It can be seen that in the dry region where (4.38) r . o, Equations (4.35) and (4.36) are the same as for a flat-slab as are the boundary conditions (4.37) and (4.38). In the wet-region the combined energy equation obtained by eliminating II(,'; )2 between Equations (4.35) and (4.36) is the same as Equation (4.15). Therefore it is clear that the final expressions Chapter 4. QUASI STEADY ANALYTICAL MODEL^ - 33 derived for the heat flux and effective thermal conductivity for flat slab can be applied to round-pipe insulations by taking the equivalent flat-slab thickness as ro Leg = ro ln( —) (4.39) where ri is the inner radius of insulation. 4.7 Concluding Remarks The analytical expressions for temperature distribution, heat flux and effective thermal conductivity in the presence of moisture gain have been derived. The main design and operating parameters on which the effective thermal conductivity depends are identified. The results of the parametric study using the above expressions and the verification of the model will be presented in the next Chapter. Chapter 5 ANALYTICAL RESULTS AND DISCUSSIONS 5.1 Introduction The quasi-steady model presented in the previous chapter was used to perform a detailed parametric study. It is useful to examine carefully the variables on which the heat flux and the effective thermal conductivity depend. The main dimensionless groups involved are B, B,, and 0. The Lewis number for heat and mass transfer at outer surface may be assumed to be approximately unity [39]. With this assumption the dimensionless groups depend on three basic characteristics of the insulation viz, the thermal conductivity, the porosity and the thickness. The other variable is the outside heat transfer coefficient. The parameter (3 and the function 1(8), depend on the ambient temperature, the ambient relative humidity and the cold surface temperature. In view of the above seven independent design and operating variables involved, it is difficult to represent the results in the most general form. Therefore for ease of application the results are presented in dimensional form, covering ranges of parameters that are of practical interest. 5.2 Comparison with Experimental Results The predictions made by the analytical model derived in previous chapter for the heat flux and temperature distribution were compared with the experimental data presented in Ref. [23] for a flat-slab. The comparison of the measured [23] and computed temperature 34 Chapter 5. ANALYTICAL RESULTS AND DISCUSSIONS^ 35 1.00 Wet Region Line^tAnalYlleal results Symbols :Experimental data 0.75 o Dry Wet Region 0.50 0.25 ^ RH: 96% - RH: 80% RH: 70% Dry Wet Region 0.00 0.00 , 0.25 ^ 0.50 ^ 0.75 • 'N ^ 1.00 Nondimensional Distance along the Slab Z/L Figure 5.1: Temperature distribution in slab; Comparison with experimental Data [23]. T. (°C) (°C) 33.0 32.0 33.1 33.2 6.8 6.8 6.8 6.8 Tcold Relative Humidity (%) 96 90 80 70 Heat Flux Measured (W/m2) 35.0 32.0 30.2 27.4 Heat Flux Predicted (W /m2) 38.46 34 32.95 29.56 Table 5.1: Comparison of Predicted and Measured [23} Heat Flux distribution is shown in Fig. 5.1. The conditions are as following: T. = 30°C, T2 = 6.8°C, dry density of insulation slab p. 53 kg/m2, L = 66 mm, and the relative humidities are 96%, 80% and 70% respectively. The computed temperature distributions agree well with the measurements. The computed wet-dry interface is also indicated in the figure. The comparison of the computed and measured heat flux at the cold plate is shown in Table 5.1. There is satisfactory agreement between the calculated and measured heat flux. Chapter 5. ANALYTICAL RESULTS AND DISCUSSIONS ^ 36 1.00 e 0.90 /". 0.80 .2 0.70 0.60 0.50 ^ 0.00^0.25^0.50^0.75 1.00 Nondimensional Distance along the Slab ZIL Figure 5.2: Vapor Flux Distribution in slab, Comparison with experimental data [23]. The comparison of predicted and measured vapor flux in the slab is shown in Fig. 5.2. The conditions are the same as in Fig. 5.1. 5.3 Comparison with Literature Reports Figure 5.3 shows a comparison of the predicted heat flux using the present analytical model with the numerical results of Tao et al. [38]. The nondimensional heat flux was defined by Tao et al. [38] as the ratio of sensible heat flux with condensation in the slab to the dry-state heat flux with the same temperatures across the slab. This heat flux is shown as Q' in Fig 5.4. The slight discrepancy between the results can be attributed to the physical properties used. The present analytical model assumes constant mean physical properties while the numerical model [38] uses variable properties that depend on the liquid content. In the present study the heat flux at the cold plate was defined by Equation (4.29) . Chapter 5. ANALYTICAL RESULTS AND DISCUSSIONS ^ 37 2.5 1 = 0.5 Present analytical Model Numerical Model [38] 0.0 0 5^10^15 Cold Plate Temperature T2 (C) 20 Figure 5.3: Comparison of analytical and numerical results, Conditions: L = 0.099 m, T. = 20°C, kdry = 0.037 W/m K, h* = 12 W/m2 K. This includes the sensible heat flux and the latent heat released by the vapor condensing at the cold plate. The variation of the nondimensional total heat flux is indicated as Q in Fig. 5.3. It is seen that the magnitude of Q is more than twice that of Q'. The variation of the total heat flux with cold surface temperature will be discussed later in this section. Modi and Benner [21, 22] measured the effective thermal conductivity of flat-slab insulations for different relative humidities and two values of the cold surface temperature. There is considerable scatter in their experimental data. However, there is an identifiable quasi-steady regime in the graphs. The quasi-steady values of the measured, and computed values of (Iceff/kdry) are given in Table 5.2. There is satisfactory agreement when the cold surface temperature is 3°C, but for a cold surface temperature of 10°C there is some disagreement. Chapter 5. ANALYTICAL RESULTS AND DISCUSSIONS ^ T. Td (°C) (°C) 30 30 30 30 30 30 3 3 3 10 10 10 Relative Humidity (%) 70 50 30 70 50 30 ke f f I kdry Measured 2.22 1.75 1.375 1.72 1.4 1.0 38 Ice f f I kdry Predicted 2.37 1.807 1.285 2.50 1.72 1.02 Table 5.2: Comparison of Predicted and Measured [21,22] Effective Thermal Conductivity 5.4 Effects of Condensation on Heat Transfer The variation of the heat flux with the temperature difference is shown in Figs. 5.4-5.6, for different values of the ambient relative humidity and slab thickness. The variation of the (icaffikdry) with the temperature difference is shown in Figs. 5.7 - 5.9, for different values of the ambient relative humidity and slab thickness. It is instructive to explain in physical terms the shape of these curves. Consider a practical situation where the ambient temperature and the ambient humidity are maintained constant and the cold surface temperature is decreased progressively starting from a value close to the ambient temperature. When the cold surface temperature T2 is such that the saturation concentration, m(T2) > ma, the vapor concentration in the ambient, no vapour diffusion will occur. The value of (ken /kdry) is then equal to unity and the heat flux is only due to sensible heat transfer. When m(T2) < ma, vapour will diffuse through the dry slab and condense at the cold surface. This small quantity of condensed liquid may be assumed to have negligible influence on the transport processes in the porous slab. The heat flux leaving the impermeable cold surface will include both the sensible heat flux and the latent heat flux due to condensation. The heat flux at the cold surface increases with the temperature difference due to the increased sensible heat Chapter 5. ANALYTICAL RESULTS AND DISCUSSIONS^ 39 150 Ta ■ 40°C k ■ 0.035 W/m X h ■ 8 W/m2X E ■ 98% 0 0 5 10 15 20 25 30 35 40 Temperature Difference Between Warm And Cold Sides Figure 5.4: Variation of heat flux with temperature difference, T. = 40°C. 75 0 Ta ■ 30°C k ■ 0.035 W/m X h ■ 8 W/m2IC E ■■ 98% 0 5 10 15 RH■90% L■25 mm • 20 , 40%, 25 30 Temperature Difference Between Warm And Cold Sides Figure 5.5: Variation of heat flux with temperature difference, T. = 30°C. ^ Chapter 5. ANALYTICAL RESULTS AND DISCUSSIONS ^ 40 50 40 ri g 30 Ta ■ 20°C k ■ 0.035 W/m K h 8 W/m2 K E ■ 98% g 20 ^..' 10 L■25 mm-% . ,- ....-^ • „...• •'^ 50%,...". . .,^.-..^ .^ .^.... . .. .. •...- ,^..^L■75 mm ,^ .. ... . . .^.^ RH-90% ^. . ...•^,..-^....., .^ ..^,-^...-....^. 7 0 0^5^10^15^20 Temperature Difference Between Warm And Cold Sides Figure 5.6: Variation of heat flux with temperature difference, Ta = 20°C. 0^5^10^15^20^25^30^35^40 Temperature Difference Between Warm And Cold Sides Figure 5.7: Variation of effective thermal conductivity ratio with temperature difference, = 40°C. Chapter 5. ANALYTICAL RESULTS AND DISCUSSIONS ^ 41 5.0 -r1 V 0 0.0 o a .14 $4 Ei .22.0 1 .0 0 0.0 0^5^10^15^20^25^30 Temperature Difference Between Warm And Cold Sides Figure 5.8: Variation of effective thermal conductivity ratio with temperature difference, Ta = 30°C. transfer and the increased vapour flux. The ratio (Iceffilcdry) also increases due to the same reason. At the turning point P in Fig. 5.7 and the corresponding point Q in Fig. 5.4, the vapour condensation is just beginning to occur within the insulation slab. At this point the wet-dry interface is just at the cold surface and the dry length zd given by Equation (4.20) is equal to one. When the temperature difference is increased further, the wet-dry interface moves towards the exposed surface i.e. zd decreases progressively. Figures. 5.45.6 show that the heat flux increases continuously with increasing temperature difference. However, the rate of increase of the total heat flux is reduced because of the decrease in the function I'M with temperature. Due to the different rates of variation of and (T. — T2), the value of (ke f f kdry) increases up to P and then decreases slowly. The definition of ke f f is based on Equation (4.29), which is the design energy equation Chapter 5. ANALYTICAL RESULTS AND DISCUSSIONS^ 42 25 0 0 V a M2.0 .-1 0 g 1•5 0 . L■50 mm ___ /^------ _ _ _ L■75^_______________ • g g2 to R11■50% Ei o Ta • 20°C k • 0.035 W/m X h 8 W/siX E ■ 98% OZ 0.0 0^5^10^15^20 Temperature Difference Between Warm And Cold Sides Figure 5.9: Variation of effective thermal conductivity ratio with temperature difference, Ta = 20°C. for dry insulation systems. This definition facilitates the inclusion of the effects of condensation in more broader energy conservation studies by the use of a modified thermal conductivity value. This definition, however, leads to the somewhat unusual variation of keff with temperature difference as discussed earlier. The heat flux and the effective thermal conductivity increase with decreasing slab thickness. As expected higher ambient temperatures and humidities lead to higher effective thermal conductivities. The sensitivity of the effective thermal conductivity to the heat transfer coefficient, the thermal conductivity and porosity of the insulation is shown in Figs. 5.10, 5.11 and 5.12 respectively. Its magnitude increases with decreasing ha and k and increasing porosity. Design data for round-pipe insulations may be easily deduced from Figs. 5.7-5.9 using • Chapter 5. ANALYTICAL RESULTS AND DISCUSSIONS ^ >4 4.0 ko 'CI 3.0 "• ▪ 11^2.0 0 14 0 O -443 *-1 4, 1) g 1.0 0.0 0^5^10^15^20^25^30 Temperature Difference Between Warm and Cold Sides Figure 5.10: Effect of h on the variation of the effective thermal conductivity ratio. ,_ 1 / ,.. _. _ - ___— •^---^"---- -- - 14 /„P 1r 14 Ta ■ 30°C L ■ 67 mm h ■ 8 W/m2 X 6 ■ 98% RH ■ 96% 0.0 ^ K=0.038 W/M K ^ K=0.035 W/M K K=0.033 W/M — — - K=0.030 W /m K I^1^I^I^1^I^I^I^I^I^, 0 5^10 ^ 15^20^25 ^ 30 Temperature Difference Between Warm and Cold Sides Figure 5.11: Effect of K on the variation of the effective thermal conductivity ratio. 43 Chapter 5. ANALYTICAL RESULTS AND DISCUSSIONS ^ Ta ■ 30°C h ■ 8 W/m2 X L ■ 67 mm k ■ 0.035 W/m RE ■ 96% 44 ^ EPSIW98% ^ EPSIV=96% EPSIV=94% 0.0 0^5^10^15^20^25^30 Temperature Difference Between Warm and Cold Sides Figure 5.12: Effect of c on the variation of the effective thermal conductivity ratio. the equivalent flat-slab length given by Equation (4.39). For example, a 35 mm thick pipe insulation applied on a pipe of diameter 50 mm has an equivalent fiat-slab thickness of about 52.5mm. The curves in Figs. 5.7-5.9 may now be used to find the effective thermal conductivity for L = 52.5 mm. 5.5 Concluding Remarks An analytical model is developed for the quasi-steady phase of the heat and water vapor transfer through flat and pipe insulations. The computational effort involves the solution of the nonlinear equation (4.18), by a trial and error method, to determine the dimensionless temperature 8, at the wet-dry interface. With a knowledge of Oc, the heat flux and the effective thermal conductivity are obtained by directed substitution in Equation (4.28) and (4.31) respectively. The same procedure is appliable for pipe insulations with Chapter 5. ANALYTICAL RESULTS AND DISCUSSIONS ^ 45 the use of the equivalent flat slab thickness given by Equation (4.39). The effective thermal conductivity and the heat flux in the presence of condensation depend on seven independent design and operating variables of system. The effective thermal conductivity increases with increasing temperature difference across the insulation and reaches a maximum when condensation begins to occur in the slab. Thereafter, it decreases slightly as the temperature difference is increased further. For practical operating conditions the effective thermal conductivity varies from about 1.5 to 15 times the dry-state value. The results of the parametric study are presented as curves which may be used to estimate the effective thermal conductivity over practical conditions for design purposes. The same curves can be used for pipe insulations with the use of equivalent flat-slab thickness. Chapter 6 TRANSIENT NUMERICAL MODEL 6.1 Introduction The moisture transport processes in this physical model have four main stages which have been stated in Chapter One. There is an initial transient period which may last up to about an hour for most practical situations before a quasi-steady phase is attained. The quasi-steady state may last anywhere up to several days before the accumulation and transport of liquid water begins to significantly affect the heat transfer due to reevaporation and changes in physical properties of the insulation. The quasi-steady analytical model presented in Chapter Four has given reasonable prediction of the heat and mass transfer parameters for the quasi-steady phase. However, many simplifying assumptions have been made in order to obtain the quasi-steady solutions. For numerical studies on transient moisture transfer through porous insulation, only work on the effect of condensation during the early time period has been reported in the literature [35-37]. In this Chapter, the long-term transient moisture transport processes are analyzed, and the significant transport mechanisms are identified. Accordingly, the rigorous and detailed formulation of heat and moisture transfer for all the stages is presented. The long term effects of condensation are investigated numerically, and the mobile condensate is taken into account. 46 Chapter 6. TRANSIENT NUMERICAL MODEL^ 47 6.2 Analysis and Formulation 6.2.1 General Formulation The assumptions made in the transient formulation are the same as that stated in Chapter Three. The formulation is based on the local volume-averaged technique. In order to generalize the mathematical formulation, the following new set of nondimensional variables have been introduced: Dimensionless distance, Dimensionless temperature, 0 Dimensionless density, Pi Ta-T2 P? Pij t. Dimensionless time, L2/a,I01, Dimensionless condensation rate, G* Dimensionless heat capacity, cp Dimensionless thermal conductivity, = = Dimensionless vapor diffusivity, 1)(0) Dimensionless liquid diffusivity, DAM The symbols are also defined in the Nomenclature. = • Poao,ff/L2 cp„ — )4 P421 a 0,ef = 0,eff (6.1) Chapter 6. TRANSIENT NUMERICAL MODEL^ 48 The governing equations (3.3)-(3.5) in terms of the dimensionless variables are: Vapor diffusion equation; any a a(p) _ EDv(0)---] —^ G* = ax az at (6.2) Liquid transport equation; a^ae,^c* aft, Pt3— az[Di3(w)— az1+ ^= at' (6.3) a^a()^UT PCP at • az^) az 1 (6.4) Energy equation; —[k(w^ - + P2G* = The corresponding constitutive correlations in terms of the dimensionless variables are: Volumetric Constraint; ear + eo 4- E.7 ^1, (6.5) Thermodynamic Relations; ^1 Pv = P3Pve, (6.6) Pa = P4P09 ) (6.7) ^1 Pv = ExP[P5(w — eref)1, (6.8) and Matta = Pa + Pv•^ (6.9) Chapter 6. TRANSIENT NUMERICAL MODEL^ 49 The dimensionless parameters are defined as the following: - • AT cp0 F2 — D ? *. =^7_2, 3 Pt, , 0 (6.10) AT p: hf .0 g P5^R,, AT The spatial average density of porous insulation is defined as P =^efiP)9 67P7.^ (6.11) The mass fraction weighted average quantity of heat capacity is defined as cp = EcrPacna €13P0cA15 c7P7cw (6.12) The problem is modeled as a porous insulation with impermeable and adiabatic vertical boundaries. The upper horizontal boundary is impermeable and subjected to a cold temperature, and the lower one is exposed to moist ambient at a warm and humid condition. The insulation slab is assumed to be fully dry and initially has a uniform temperature. It encounters a temperature drop at upper surface and is suddenly exposed to the ambient air at the lower surface. Therefore, the initial condition for the model is Chapter 6. TRANSIENT NUMERICAL MODEL ^ 50 pv(z,0) = p„,, (6.13) ei3(z,0) = 0 11(z,0) = Ga. The boundary conditions to complete the formulation are ae 5; 10,0= B[90 — 19(0,t)], aPv Oz ko,o= Bm[pv, - pi,(0,t)b 0(1,0 = 02,^ pu(1,t)==!7(02), ^ (6.14) (6.15) (6.16) (6.17) where g(0) is saturated vapor density which can be represented by the Clausius-Clapeyron relation 1^1 g(0) = exp[P5(— n — n )j/P382. (6.18) U2^Uref 6.2.2 Stage 1: Initial Process Initially, it is assumed that there is no water vapor condensate inside the insulation slab. The temperature and vapor density fields are governed by: Vapor diffusion equation; Chapter 6. TRANSIENT NUMERICAL MODEL^ 51 apt?,^apv 1[1),,(61)1 -= 67^ , az (6.19) a r,^ a19,^ae — azilgAu)-a–;1= PCP-a-t-, (6.20) Energy equation; where the dry thermal conductivity is the only thermal property which is a function of temperature. Although there might be some initial adsorbed water inside the slab in practical situations, it is assumed that the insulation slab is totally dry initially, and the liquid fraction is zero, - 0.^ (6.21) 6.2.3 Stage 2: Immobile Liquid Accumulation Stage 2 is identified as that when condensation occurs within the insulation slab. The air-vapor mixture is saturated when the actual vapor density in a region is larger than or equal to the local saturated vapor density which can be represented by equation (6.18). The new two-phase regions (dry-wet) are then determined and the wet-dry boundary is located where the actual vapor density is just equal to the local saturation vapor density. It should be noted physically, due to the inhomogeneities of actual fibrous insulation, the wet-dry boundary is not a single line but has a finite volume. However according to the local average technique, this fuzzy boundary volume can be averaged and represented by a line. • Dry Region In the dry region, the condensation rate G* is set to zero, as well as the liquid fraction for the above mentioned reasons. Equations (6.20)-(6.19) can be still used Chapter 6. TRANSIENT NUMERICAL MODEL ^ 52 for solving the temperature and vapor density fields. • Wet Region During this stage, the liquid which is accumulated due to condensation is small and has an insignificant effect on the properties. Therefore, the liquid is practically immobile, and is accumulated at a constant rate. The temperature and vapor density fields are governed by Equations (6.2) and (6.3). Since the vapor is saturated in the wet region, vapor density is the only property which is a function of the local temperature, and it can be represented by the Clausius- Clapeyronrelation (6.18). Differentiating of Equation (6.18) against t and z respectively, one obtains: and apt,^ao g'(e)— LIZ az' (6.22) apt,^ae wet g1(8)— at^at. (6.23) Eliminating the condensation rate G* between Equations (6.2) and (6.3), and making use of relations (6.22) and (6.23), the vapor diffusion equation (6.2) and energy equation (6.3) can be combined as: ao aE a ae P2g(0)—1 = -[(k(W) P2D„(0)V(6))— ]. [pc., + P2E-d(o)]— at + a t az az (6.24) Together with boundary conditions on the cold side (6.16) and on the wet-dry interface, which will be stated later, Equation (6.24) can be used to solve for the temperature field in the wet region without the need to deal with an unknown source term. Once the temperature distribution in the wet region is known, the Chapter 6. TRANSIENT NUMERICAL MODEL^ 53 saturated vapor density profile in the wet region can be obtained easily from Eq. (6.18), and the condensation rate can be solved either from the vapor diffusion equation (6.2) or the energy equation (6.3). Since the liquid is in a pendulous state, the liquid diffusion equation (6.4) is reduced to aEa G* at = pi (6.25) • 6.2.4 Stage 3: Mobile Liquid Diffusion When the amount of liquid has been collected to some extent, the liquid can no longer be trapped in the pores. The liquid will diffuse in the direction of liquid content gradient. Since all the vapor reaching the impermeable cold plate is considered to have condensed at the cold plate, the liquid fraction is highest near the cold plate at the beginning, and diffuses towards the warmer, and drier direction. The critical liquid fraction eo denotes the liquid fraction below which the liquid is immobile. It is a complex function of the structure of fibers, liquid pressure, surface tension etc.. According to the experimental results of Cid and Crausse [28] for fibrous insulation with dry state density of 4. 53 kg/m3, the liquid diffusivity abruptly becomes significant (about 10 m2/s) when the liquid fraction is between 0.035-0.075. After the liquid starts to flow, there are two situations to be distinguished. In the first case, the liquid diffuses inside the wet region. There is no apparent liquid front, the liquid diffusion equation is applied where the liquid fraction exceeds the critical value. In the second case, the liquid flows out the wet region. It enlarges the wet region and makes the wet-dry interface move due to the mobile liquid front. The temperature and vapor density fields in the wet region are formulated in the same Chapter 6. TRANSIENT NUMERICAL MODEL^ 54 way as that in stage 2. However, treatment of the moving wet-dry interface is needed for this stage, and will be stated in the section of solution methodology. The Equation (6.4) was developed to govern the liquid diffusion. However, the liquid diffusivity is a phenomenological quantity. Without any available experimental data, the liquid diffusion coefficient is not a known quantity. Cid and Crausse [28] measured the liquid diffusion coefficient for several commonly used fibrous insulation. Their liquid diffusion coefficient Di(w) is defined as in the following liquid transport equation: aa * Ow + e^ 11 w =^ —[.D1(w)-e; 4 3 l-K (w )J + —.^(6.26) at* Ox^ PO where w is weight liquid content, K(w) is hydraulic conductivity and is the function of the liquid saturation S. The relations between w, S and liquid fraction ep are (6.27) = ^ 43E0 = COS = (1 — ca)S,^ (6.28) (6.29) where 60 is the void fraction of porous insulation. Note that some important relations exist between the differentiations of w, S, K(w) and liquid fraction co: OK (w) 0.1C(w) Ow Ox — Ow ax' OK(w) OK(w) as — eo as aw _ 1 p;,` OK(w) 4 as (6.30) (6.31) Chapter 6. TRANSIENT NUMERICAL MODEL^ ow4 ô€13 _, — ox^(4; ax' aw _.= 4 ae,, at*^pc*, at 55 (6.32) (6.33) Substituting the above relations into the Equation (6.26), the liquid diffusion equation with experimentally determined liquid diffusion coefficient, in terms of the liquid fraction, is obtained: ae,3^a^ae„,^1 OK(w) act,^r at* = a—x[Dr(Ei9)— axi+ Eo as ax + --7 pf3. (6.34) Applying the same dimensionless variables as in Eq. (6.1), the liquid diffusion equation in nondimensional form is ae,^a^ae,3^(9E, G* at = -a-i[D1(")-ai]+ F(S)-; + --P1' (6.35) where Pi is defined in (6.10), and F(S) can be viewed as a coefficient due to hydraulic conductivity, and is defined as: 1 OK(w) L F(S)= ^.^ €0 as cq; ,ef f (6.36) The liquid diffusion coefficient and hydraulic conductivity in terms of the liquid saturation are regressed from Cid and Crausse's experimental data [28]. The following polynomials have less than 5% discrepancy from the original experimental data: Dr(S) = 1.29(10)-8 + 9.913(10)-7S + 1.6018(10)-6S2 — 7.8408(10)-6V +1.053(10)-5S4;^0.45 > S > 0.05, (6.37) Chapter 6. TRANSIENT NUMERICAL MODEL ^ K(S) = 56 4.502(10)-i — 4.658(10)-6S + 1.05486(10)-3V — 2.56279(10)-3S-3 —1.7555(10)'S4 + 3.956(10)-3S5; ^0.45 > S > 0.05. (6.38) 6.2.5 Stage 4: Long Term Liquid Accumulation and Diffusion The liquid front eventually reaches the exposed surface, the entire insulation slab becomes wet. The liquid is accumulated continuously because the water vapor from the ambient air keeps condensing inside the slab and on the cold plate. The liquid will move under the liquid fraction gradient, and will be redistributed continuously. The formulation for this stages is quite the same as that for the wet region in previous stage except that convective boundary conditions are introduced into the formulation for the wet region. The treatment of these boundary conditions on the exposed surface will be stated in detail in the next section. 6.3 Solution Methodology The above introduced formulation can not be solved analytically. The numerical scheme is based on the finite difference form of the above-stated formulation. The solution domain is divided into a finite number of control volumes or cells. A grid-point is placed at the geometric center of each control volume. This arrangement has the following advantages: • The value of the general variable which is available at the center of the control volume represents the volume averaged value over the control volume. • The physical properties and the source terms can be calculated at the center of the control volume. Chapter 6. TRANSIENT NUMERICAL MODEL ^ 57 k^(6 z) e ..1 w w ^ P^ E •^ ______.Z A z^ kis^ interface^interface E ..,^Z A Z interface^interface (a) ell^ (b) Figure 6.1: Grid-point Cluster for One Dimensional Transport Problem: (a). Control Volume for the Internal Points, (b). Control Volume for the Boundary Points. • Discontinuities at the boundaries can be conveniently handled by locating boundary cells where the discontinuities occur. 6.3.1 Discretized Formulations The discretized system of the governing equations is formed with a grid-point cluster shown in Fig. 6.1. The dashed line show the faces of the control volume. For the onedimensional problem under consideration, we assume an unit thickness in the directions other than z. Modification of these control volumes near the boundaries is straightforward. For the numerical purpose, the governing equations described in the previous section can be represented by the general transport equation a^a^as:k -&(4)) = a-i[B(Tz.)] + S4^(6.39) where (I) can be replaced for different equations, B is a general diffusion coefficient, and St Chapter 6. TRANSIENT NUMERICAL MODEL^ 58 a general source term. The discretized form of Equation (6.39) is obtained by integrating over the respective control volumes, i.e., ft+At ft+at al,^ti-At it+At a^al, 1 it — at dtdz = ft^4^ az (B(-5---z ) + Ndzdt.^(6.40) [— For detailed derivation method of the discretized equation, one could refer to Patankar's work [42]. The fully implicit scheme is chosen for the requirements of stability, simplity and physically satisfactory behavior. In general, the central difference form is used for internal nodes and the backward or forward difference used for the boundary nodes. The backward difference form is used for the time derivative. The discretized forms of the formulations for this problem are given as the following. Stage 1: Initial Process The discretized equation for energy transfer is: at the internal nodes, apOp = aEOE + awew + b, (6.41) ke aE = ^ (6x)e' (6.42) kw aw = ^ (64,' (6.43) ,,o^Pcpllz At ' (6.44) b= 4,4, (6.45) where -'1' - ap = aE + aw + 4. (6.46) For the dry slab, the conductivity variation only depends on temperature distribution. The interface conductivities Ice, kw in terms of these grid-point values are Chapter 6. TRANSIENT NUMERICAL MODEL^ 7._^2kpkE kp + kE' 2kpkw kw — • Ne = ^ kp + kw 59 (6.47) (6.48) At z = 0, half cell near boundary is introduced, the boundary condition is pcP Az At 2 (9 p - OCIDI ) = ke (6z)e (9E — Op) — h(OP — O.), (6.49) where h is dimensionless convective heat transfer coefficient in the ambient, and € 7 is constant for a dry slab. The vapor diffusion equation has a similar discretized form to that of the energy equation. The equation system (6.41)-(6.46) can be used for vapor diffusion, if the corresponding variable and coefficients are replaced. The boundary condition for the vapor diffusion equation at z = 0 is 677 ^,^Dye — kPv P — P° p) = At 2^'^ v' 6z^' ^Pv,p) — hni(pv,p — p,,,a), (6.50) where kr, is dimensionless mass transfer coefficient in the ambient. At z = 0, the temperature at the cold plate is given; no problems are introduced in the discretized energy equation. The vapor density at this boundary is equal to the saturated vapor density at T2. Stage 2: Immobile Liquid Accumulation In this stage, the discretization equations for the dry region are the same as that in stage one. For the wet region, the combined energy equation (6.24) can not be solved directly because the gas phase volume fraction Gy is not known. An iterative procedure is needed to solve for the temperature field in the wet region. Chapter 6. TRANSIENT NUMERICAL MODEL^ 60 The first derivative term of c,, is omitted for the first iteration, the combined energy equation becomes ae a^ae [pcp + P2e.yi(0)]-5--t- = --a—z-[(k(w) + P2D„(0)i(0))--a-z], which has the following discretized form for the internal grids, [p cp + p2e_ygt( 9)]p rt ( op 07,) = i( k(tV)+ P2aD:(9)g ' (9)1, (I 9 E .. 0 p) PM+ P2 f: (8 )g '(9)118 (6.51) 1 (6.52) (Op ■.■ 0 W) . In the dry-wet interface, the boundary condition for the combined energy equation can be derived from the energy balance at the interface. Assuming that the grid size can be taken fine enough so that the dry-wet interface is approximately located at the interface of control volumes, as shown in Fig. 6.2, the energy balance of the grid where the dry-wet interface locates can be written as [PCP + P2e-ygi (0)]13 11(9 p — Op) = P(w)-1-P2 ap:(9),91 PA. (0 E _ O) —[11(OP — Ow) + (P2f:)w (Pp — Pw)]. 1 (6.53) In the above discretized equation, the vapor densities in grid P and W are unknown at the present time step. As a numerical approximation, the values at the previous time step are used instead. Thus the vapor density term is treated as a source term, and an iterative procedure is performed to minimize the error. The vapor density is calculated directly from the Clausius- Clapeyron equation (6.18) in the wet region, and is computed from the vapor diffusion equation in the dry region. The condensation rate in the wet region can be computed explicitly from the vapor diffusion equation or from energy equation. The vapor diffusion equation was used to calculate the condensation rate, Chapter 6. TRANSIENT NUMERICAL MODEL ^ 61 Dry Region^wet Region ;P ; Wet-dry Interface Figure 6.2: Grid-point Cluster for the Interface Location of Wet-dry Regions Gp^ a 7z 2 le lPv,E f^ - - —- I^Di, Pv,P) ^ (Az)2 lev(P.,P P.N) (67Pv)P — ■EyPv t (6.54) The following boundary conditions are used for solving the condensation rate: at the wet-dry interface, z = zd, = 0,^ (6.55) at the cold plate, z 1, ,, 2D„ e = Az2 lPv,P P.,w), (6.56) which indicates that all vapor which reaches the cold plate will condense there. The liquid fraction is computed from GaP 0 EI3'P = + E". (6.57) Once the gas phase volumetric fraction El. is obtained from Eq.(6.5), a source term 0 P2g(0)61'PZ"-P is added to the right hand side of Equation (6.52), and the discretized energy equation for the second iteration becomes 62 Chapter 6. TRANSIENT NUMERICAL MODEL^ [pcp + p2e7g1(95111(0 p _ 90p) :____ j(k(w)-FP2AD;(8)9i(eyie [(k(w)+P2D„(8)91 _p2g (0 Az ) (9E _ 9 p) (Ow (9 P Cry ,P -9,67)-€137,P . 1 (6.58) The other variables can be obtained from the corresponding formula for the first iteration. Satisfactory results can be obtained after 4-5 time iteration. Stage 3: Mobile Liquid Diffusion The discretized formulations and solution procedure for solving the temperature are similar to that of stage 2. The discretized form of the liquid transport equation is Az^DI ler^, Di ,^ N G* A k€0,E 60,P) ^ keis p— €13,w)-1- F(S) IP (6.13,E — etLp)+ (€0,P — e°0,P) =^ Az^ Az ' At Pi (6.59) The liquid starts to diffuse from the grid close to the cold boundary. Before the liquid front reaches the quasi-steady wet-dry boundary, no special boundary formula is needed for the liquid front. The liquid front is traced where the liquid volumetric fraction is equal to, or larger than the critical liquid fraction cox,. = 0.5€0 (taken from Cid and Crausse's experimental data). When the liquid flows out the quasi-steady wet-dry interface, the new wet-dry interface is located at the interface of two grids where the liquid fraction in one grid exceeds the critical value in the last time step. that is assuming at the i time step, the dry-wet interface is located at the interface between grids j-1 and j, and the value of liquid fraction in grid j reaches the critical value, for the i+1 time step, the new wet-dry interface will move to interface of grid j-1, and j-2. This approximation for tracing the moving boundary is reasonable as long as the grid size is taken fine enough. Chapter 6. TRANSIENT NUMERICAL MODEL^ 63 At z = 1, the boundary condition is Az ,^0 \^I,„^ G* Az Lut kfo,P — Ef3,13) = Az (EON Efi,P)+ F(S) Ip (Ei3,P EON) + (6.60) the moving wet-dry boundary, At the Az ,^DI le, G* Az 0\ .^(6.61) kEl3 P E ) ^ kefi E Ef3 P) F(S) IP (€13,p — Es,w) — " Pi 2 Az^ Stage 4: Long Term Liquid Accumulation and Diffusion In this stage, the wet-dry interface reaches the exposed surface. The entire insulation slab is wet. The solution procedure for temperature is the same as mentioned for the wet region. The discretized equations (6.54) and (6.58) for the wet region are applicable for the internal grids. The boundary condition on the exposed surface is introduced as {PCP + p2e7g/(8)} ip tist.(ep _ eop) ,_ t(k(w)-FP2AD,;(91g1(9)11. (0E _ op) —11.(8p — 0.) — P2h,„(pp — pa)]. 1 (6.62) Similar boundary conditions on the exposed surface are also introduced for computing the condensation rate and liquid fraction. 6.3.2 The Solution Algorithm A Fortran 77 computer code has been developed to solve the above-stated system of the coupled nonlinear equations. Starting with the initial values of the temperature and vapor density, the solutions of the above model for the four stages involve several formats, as well as several iteration procedures. These formats are: Chapter 6. TRANSIENT NUMERICAL MODEL^ 64 • Format 1. This format was followed for any time and location along slab for which no condensation occurs in the insulation slab. 1. Initial values of 0 and pt, are given from Equation (6.13), (step 1 is applied only at t = 0). 2. The condensation rate is set to zero, 0 is obtained from Equation (6.20) and A, is obtained from Equation(6.19). 3. The saturated vapor density p„,, is obtained from Eq. (6.18) based on the computed temperature profile. 4. If, at any location, the pi, obtained in step 2 is greater than or equal to pi,,, obtained from step 3, then Format 2 is adopted for the next time step. Otherwise the procedure starting from step 2 is repeated for the next time step. • Format 2. This format was followed for any time and location for which the wet region is established and co < co,,.. 1. The temperature in the wet region is solved from Equation (6.52) with boundary conditions (6.53) and (6.16) as the first iteration. 2. The vapor density in the wet region is obtained from Equation (6.18). 3. The temperature and vapor density in the dry region are obtained from Equation (6.20) and (6.19). 4. The condensation rate is obtained from Equations (6.54) -(6.56). 5. The liquid phase fraction co is obtained from Equation (6.57), and the gas phase fraction ey is obtained from the volumetric constraint Equation (6.5). 6. The second iteration value of the temperature in the wet region is obtained from Equation (6.58). Chapter 6. TRANSIENT NUMERICAL MODEL ^ 65 7. The same procedures are followed from step 2 to step 5 to obtain the second iteration values of vapor density, condensation rate and the liquid fraction. 8. If, at any location, the liquid fraction obtained from step 7 is greater than or equal to the critical liquid fraction, Format 3 is adopted for the next time step. Otherwise, the procedure starting from step 1 is repeated for the next time step. • Format 3. This format was followed for any time and location for which the liquid volumetric fraction exceeds the critical liquid volumetric fraction. 1. The temperature field in the wet region is solved from Equation (6.52) with boundary conditions (6.53) and (6.16) as the first iteration. 2. The vapor density in the wet region is obtained from Equation (6.18). 3. The temperature and vapor density in the dry region are obtained from Equation (6.20) and (6.19). 4. The condensation rate is obtained from Equations (6.54) -(6.56). 5. The liquid phase fraction co is obtained from Equation (6.59), with the boundary conditions (6.60), (6.61). The gas phase fraction e-y is obtained from the volumetric constraint Equation (6.5). 6. The second iteration value of the temperature in the wet region is obtained from Equation (6.58). 7. The same procedures are followed from step 2 to step 5 to obtained the second iteration values of vapor density, condensation rate and the liquid fraction. 8. If zd=0, Format 4 is adopted for the next time step. Otherwise, the procedure starting from step 1 is repeated for the next time step. Chapter 6. TRANSIENT NUMERICAL MODEL^ 66 • Format 4 This format was followed for any time and location for which all the locations are wet. 1. The first iteration value of temperature is obtained from Equation (6.52) with boundary conditions (6.16) and (6.6.62). 2. The vapor density in the wet region is obtained from Equation (6.18). 4. The condensation rate is obtained from Equations (6.54),(6.56) and (6.63). 5. The liquid phase fraction efi is obtained from Equation (6.59), with the boundary conditions (6.60), (6.61). The gas phase fraction e.,,, is obtained from the volumetric constraint Equation (6.5). 6. The second iteration value of the temperature in the wet region is obtained from Equation (6.58). 7. The same procedures are followed from step 2 to step 5 to obtained the second iteration values of vapor density, condensation rate and the liquid fraction. The discretization equations for the above one-dimensional formulation can be solved by a delightly convenient algorithm, called TriDiagonal -Matrix Algorithm (TDMA). In principle, all the above-mentioned formulations have a form similar to the equation system (6.41)-(6.46). The nonzero coefficient of the set of equations defined by Equation (6.41) form a tri-diagonal matrix. The convergence of TDMA is fast because the boundary condition information is transmitted at once to the nodal points lying inside the solution domain. It should be pointed out that although there are two groups of formulation for the wet and dry regions respectively in the stage 2 and 3, these two groups of discretization equations should be solved in one matrix. Numerical errors may occur when the solutions of field variables are attempted separately from the dry and wet regions. There are two Chapter 6. TRANSIENT NUMERICAL MODEL ^ 67 reasons for this. First, the dry region is being reduced due to the moving interface, and the TDMA will lose effectiveness when the dimension of the matrix is less than 3. The second reason is that nodal points inside the solution domain are blind to the information of another domain. 6.4 Numerical Considerations • Grid System and Time Step Two checks were made in order to investigate the accuracy of the scheme. First, the three different grid systems, given by 50, 100 and 150 grids, were used to test the effects of the grid size on heat transfer rate, the wet-dry interface location, and the distributions of the field variables. It was found that an excellent agreement exists between the results from 100 and 150 grid systems. Next, the time step size was checked by fixing the grid size with 100 grids while varying the time step size as 30s, 60s, 120s, 180s, 240s, 360s etc.. As the results show, decreasing the time step size beyond 180 seconds did not have an effect on the numerical results. The accuracy of the wet-dry interface position depends on the selection of the grid size. For this reason, the grid size should be chosen as fine as possible. However, the required time step size becomes much smaller as the grid is refined. Considering the balance between accuracy and economy of computing time, a grid system with 100 grids and a time step of 120 seconds were chosen for the numerical computation. • Iterative Scheme In computing the field variables such as temperature 0 in the wet region, several iteration schemes were needed. In order to handle nonlinearities in the iteration solutions of the discretization equations, it was necessary to underrelax the iteration Chapter 6. TRANSIENT NUMERICAL MODEL^ 68 process. The solution was considered to be converged when the deviation of any variable from the last iterated value was within 0.01%. • Numerical Errors and Instabilities Numerical errors and instabilities can be caused by inappropriate numerical pa- rameters such as coarse grids, too large a time step and underrelaxation factors, as well as unrealistic initial values. The correct evaluation of numerical parameters can be achieved from the numerical experiments. Chapter 7 NUMERICAL RESULTS AND DISCUSSIONS 7.1 Introduction The numerical model presented in the previous Chapter was used to investigate transient heat and mass transport accounting for phase change in a fibrous insulation slab. One of the main objectives of this work was to simulate the thermal behavior of porous insulation including dynamic response during the diffusion and condensation process. The numerical formulation and its solution methodology are validated by the comparison of numerical results with the experimental data of Wijeysundera et al.. The variation and intercoupling effects of important field variables such as temperature, vapor density, condensation rate and liquid fraction are presented. Heat transfer rate through the insulation is quantified, and the moisture accumulation and distribution are numerically simulated. The long term effects of the moisture transport on thermal performance of the insulation were studied, and the interesting effects of variation of physical characteristics on the moisture and energy transfer are discussed. 7.2 Physical Data In order to examine the validity of the numerical model, the physical data used in the numerical computation are based on the fibrous insulation which was tested as an experimental sample, and also based on the real experimental conditions [23]-[25]. The physical data are summarized in Table 7.1. 69 Chapter 7. NUMERICAL RESULTS AND DISCUSSIONS ^ p(*) cp*,0 1c; c4,eff p: c;. k:, EQ. 14 53 kg/m3 840 J/kg K 0.037 W/m K 8.84(10)-7 m2/s 2600 kg/m3 836.8 J/kg K 0.762 W/m K 0.02 461.89 J/kg K 4 c;` k; p: c„ k: c; p,0 hfg 70 999.87 kg/m' 4200 J/kg K 0.57 W/m K 1.15 kg/m3 1005 J/kg K 0.57 W/m K 1882 J/kg K 2337 bar 2450000 J/kg Table 7.1: Physical Data. 7.3 Diffusion and Condensation Processes 7.3.1 Initial Process The dynamic response of porous insulation subjected to the specified boundary conditions has been investigated by case studies. Figures 7.1-7.4 illustrate the distribution of the temperature, vapor density, condensation rate, and liquid fraction in the initial stage. As can be seen in Figure 7.1, when a porous insulation slab is suddenly exposed to the ambient air on one side, and is subjected to a temperature drop on the other side, this temperature drop propagates from the cold side into the insulation slab with time. Under the given boundary conditions, the temperature field becomes quasi-steady after about half an hour (2040 seconds). The propagating behavior is also observed for vapor density as shown in Fig. 7.2, and the vapor density field reaches the quasi-steady state after about the same period. It is apparent that the drop of vapor density propagates faster than that of temperature. It can be expected as the Lewis number which is a measure of the relative importance of heat transport to the vapor transport, is less than one. The distribution of condensation rate with time is shown in Fig. 7.3. The condensation occurs first at the cold plate, and then takes place within the insulation slab after ^ Chapter 7. NUMERICAL RESULTS AND DISCUSSIONS^ 71 35 g 30 E-1 25 e i 20 tt., E 15 E-( 10 5 ^ 0.00^0.25^0.50^0.75 1.00 Dimensionless Distance Z Figure 7.1: Temperature Distribution in the Initial Stage; p* = 53 kg/m', L = 66 mm, RH = 96%, h* = 12 W/m2 K. A 0.0350 ok 0' 0300 -.V ca: 0.0250 Ta = 33°C 0.0200 ...^ ^....^T2 = 6.8°C --- A 0.0150 - g. ett0' 0100 ›. 0.0050 0.00 .^I^-^I 0.25^0.50^0.75 1.00 Dimensionless Distance Z Figure 7.2: The Distribution of Vapor Density in the Initial Stage; p* = 53 kg/m3, L = 66 mm, RH = 96%, h* = 12 W/m2 K. ^ Chapter 7. NUMERICAL RESULTS AND DISCUSSIONS^ 72 e'E 0.00200 0.00150 0.00100 2 0.00050 (.5 0.00000 I^,^,^I 1 ^ ^ 0.25^0.50^0.75 1.00 0.00 ^- Dimensionless Distance Z Figure 7.3: The Time Variation of Condensation Rate in the Initial Stage; p* 53 kg/m3, L 66 mm, RH .96%, h* = 12 W/m2 K. `c" 0.00300 — CA)^Ta = 33°C T2 = 6.8°C = 0.00250 —^ 0 trc fro P S cr 0.00200 — 0.00150 — .41 0.00100 — c. • ^ 3600s^2040s^960s^480s N ^( 0.00050 5 "O °Loom ^ ^,^,^1^, ^0.00^0.25^0.50^0.75 1.00 Dimensionless Distance Z Figure 7.4: The Time Variation of Volumetric Liquid Fraction in the Initial Stage; p* . 53 kg/m3, L 66 mm, RH = 96%, h* = 12 W/m2 K. Chapter 7. NUMERICAL RESULTS AND DISCUSSIONS ^ 73 about 480 seconds. The condensation occurs within the slab when the local saturation density is decreased below the actual vapor density. With the propagation of temperature drop in the slab, the location of the condensation front moves towards the warmer side, and stops at some place when the temperature field becomes quasi-steady. The condensation rate is highest at the cold plate because all the vapor flux is condensed there. The liquid fraction is accumulated with time due to condensation as shown in Fig. 7.4. In the initial stage, it attains its highest value at the cold plate. 7.3.2 Quasi-steady State Period Figures 7.1 and 7.2 show that after a period of time the distributions of temperature and vapor density remain constant. Fig. 7.3 shows that after about the same period, the location of the condensation front is fixed at some place which indicates the quasi-steady wet-dry interface. The quasi-steady behavior of temperature and vapor density with respect to time indicates reduction of transient thermal, mass and diffusion processes. The liquid fraction et, still increases with time due to a constant condensation rate. For the illustrative case shown in Figs. 7.1-7.4, the quasi-steady state is reached in 2040 seconds. Fig. 7.5 shows the temperature distributions along the insulation slab after a quasi-steady state has been established; the numerical results are compared with the analytical results in Chapter Five. The experimental results under the same conditions are also plotted in Fig. 7.5, which shows good agreement. '7.3.3 Liquid Accumulation and Diffusion The profiles of liquid volumetric fraction over a period of 120 hours are shown in Figure 7.6. Initially, vapor condensation occurs mainly on the cold surface. The liquid generated is accumulated at the impermeable cold plate and the nearby insulation layer. With the establishment of a quasi-steady state, liquid is accumulated in the wet region at a constant ▪▪ Chapter 7. NUMERICAL RESULTS AND DISCUSSIONS^ 74 1.00 ^ 0.75 0.50 E 0.25 Num .-96%RH Ana.-96%.RH Exp.-96%RH Num.-70%RH Ana.-70%RH Exp. 70%RH 0.00 0.00^0.25^0.50^0.75 1.00 Dimensionless Distance Z Figure 7.5: Temperature Distribution in the Quasi-steady State Period, Comparison of Numerical, Analytical and Experimental Results [22]; p* = 53 kg/m3, L = 66 mm, RH = 96%, h* = 12 W/m2 K. 0.) .C$ 120 h 0.075 100 h ^0.050 *C^0.025 0.000 0.00^0.25^0.50^0.75 1.00 Dimensionless Distance Z Figure 7.6: The Time Variation of Liquid Fraction over a Long Term Period; p* = 53 kg/m3, L = 66 mm, RH = 96%, h* = 12 W/m2 K, T. = 33°C , T2 = 6.8°C. Chapter 7. NUMERICAL RESULTS AND DISCUSSIONS^ 75 rate. The curve at 1 hour indicates that the liquid fraction at the impermeable cold boundary is higher than the critical value, and the liquid diffuses from the cold surface into the insulation slab. From the curves from 5 hours to 60 hours, we can see clearly that the liquid diffusion front moves towards the dry side. Before the liquid diffusion front arrives, the liquid fraction accumulates slowly with time. It increases sharply when the liquid diffusion front arrives. The liquid front can flow out the quasi-steady boundary as indicated by the curve at 70 hours, and eventually reaches the exposed surface, as shown by the curves from 80 hours to 120 hours. In the later stage, the liquid fraction is higher than the critical value in the entire insulation slab. The liquid diffuses from one layer into the adjoining layer with the resultant increase in liquid content and a corresponding leveling of the liquid content in the upstream layers. In addition to this transport mechanism, the local liquid content still increases due to the local condensation. Since the liquid amount produced by condensation is still highest on the impermeable cold plate, the liquid will diffuse under the liquid fraction gradient, i.e. from the impermeable surface towards the exposed surface. 7.4 Comparison with Experimental Results Wijeysundera et al. [23, 24, 25]] had measured the temperature distribution, heat flux, total moisture gain, and the liquid distribution for a range of experimental conditions through five runs. The testing times ranged form 300 to 600 hours. In the present work, case studies have been performed by comparing the computational results with the experimental data. The comparisons of the temperature distributions, heat transfer rates at the impermeable plate, the total moisture gain and liquid distribution with the experimental data are presented respectively in the following sections. Chapter 7. NUMERICAL RESULTS AND DISCUSSIONS^ 76 7.4.1 Temperature and Heat Transfer Sixty six cases of measured temperature distribution tn different times for five operating conditions have been simulated. For brevity, only the comparisons with the experimental data for 16 cases are presented in Figs. 7.7 - 7.10. There is very good agreement between the measured and computed temperature distributions. It should be pointed out that the uncertainty of the thermocouples location in the insulation slab may be responsible for some fluctuations in the measured data. The comparisons of computed and measured heat flux at the impermeable cold plate for four runs are shown in Fig. 7.11. The computed heat flux at the cold plate is summed up from two components. The first component is conductive heat flux which depends on the thermal conductivity of the moisture laden insulation; so called 'solid thermal conductivity'. The second is latent heat due to the condensation of the vapor reached the cold plate. Fig. 7.11 shows that the numerical model underpredicts the heat flux at the cold plate. The reason for the discrepancy may be due to the accuracy of moisture laden thermal conductivity models. Several possible thermal conductivity models are listed in Appendix B. Among these, the bead arrangement model is chosen for representative computation after a comparison of these models which will be discussed later. The fluctuations in the heat flux curves are due to the nonuniform temperature difference between the ambient air and the cold plate in the experiments. 7.4.2 The Moisture Gain and Liquid Transport The total moisture gain per unit volume of the slab is obtained by integrating numerically the liquid content over all the finite control volumes. The comparison of computed total moisture gain with measured data under four different operating conditions is shown in Fig. 7.12. The numerical results show good agreement with measurements for the first Chapter 7. NUMERICAL RESULTS AND DISCUSSIONS ^ 77 45 Ta=40.5°C T2=20.8°C Time: 24 hours 4 U ....e 35 E. E 4s3° 4 K25 E cl) E* 20 0.000 0.025 0.050 15 0.000^0.025^0.050 Dimensionless Distance Z Dimensionless Distance Z 45 Ta=39.9°C T2=22.2°C Time: 214.3 hours 0.000 0.025 0.050 Dimensionless Distance Z Ta=39.9°C T2=21.1°C Time: 311 hours 15 0.000^0.025^0.050 Dimensionless Distance Z Figure 7.7: The Temperature Distributions in Run 1, Comparison of Numerical Results with Experiment Data [24]; p* = 53 kg/m3, L = 62.02 mm, RH = 97%, h* = 12 W/rn3 K. Chapter 7. NUMERICAL RESULTS AND DISCUSSIONS ^ 78 45 40 P 35 Ta=40.7°C T2=14.4°C Time: 168.8 hours = Ta=41.4°C T2=11.9°C Time: 48 hours 30 25 125 et.5 20 15 10 0.000^0.025^0.050 0.000 Dimensionless Distance Z 0.025^0.050 Dimensionless Distance Z 45 4 ;3.35 Ta=41.9°C T2=1 6.3°C Time: 337 hours E.+ 30 I4 eu j 25 25 elo 520 a) E 20 a) E-1 E-1 15 15 10 0.000^0.025^0.050 10 0.000^0.025 Dimensionless Distance Z 0.050 Dimensionless Distance Z Figure 7.8: The Temperature Distributions in Run 2, Comparison of Numerical Results with Experiment Data [24]; p* = 53 kg/m3, L = 68.65 mm, RH = 96.5%, h* = 12 W/m2 K. Chapter 7. NUMERICAL RESULTS AND DISCUSSIONS ^ 79 40 35 5...) 3° Ta=36.5°C T2= 8.1°C Time: 94.8 hours E.., 25 20 a415 E4 10 5 0.000 0.025 0.050 Dimensionless Distance Z Dimensionless Distance Z 40 35° c. ) 30 Ta=35.9°C T2= 9.9°C Time: 190.8 hours 1 ,••••••■ 25 25 20 [0 Q 15 E 15 10 E 10 5 5 a 5 CLo 0 0.000^0.025^0.050 Dimensionless Distance Z Dimensionless Distance Z Figure 7.9: The Temperature Distributions in Run 3, Comparison of Numerical Results with Experiment Data [24]; p* = 53 kg/m', L = 70 mm, RH = 96%, h* = 12 W/m2 K. Chapter 7. NUMERICAL RESULTS AND DISCUSSIONS 80 45 45 Ta=40.0°C T2=20.2°C Time: 139.8 hours -_ a 0.... 35 E-4 Ta=39.9°C T2=20.4°C Time: 455.5 hours 4 35 e a eti 25 to25 5 4,) &..4 E „CD 20 20 15 ^15 0.000^0.025^0.050^0.075^0.000^0.025^0.050^0.075 Dimensionless Distance Z^Dimensionless Distance Z Ta=39.8°C Ta=39.9°C T2=20.5°C Time: 551.2 hours 0.000 0.025 0.050 Dimensionless Distance Z 0.075 T2=20.5°C Time: 620.8 hours 0.000 0.025 0.050 0.075 Dimensionless Distance Z Figure 7.10: The Temperature Distributions in Run 5, Comparison of Numerical Results with Experiment Data [24]; p* = 53 kg/m3, L = 76.13 mm, RH = 96%, h* = 12 Win? K. Chapter 7. NUMERICAL RESULTS AND DISCUSSIONS^ - 60 0 0 0 0 0 0 0 0 0 - 50 • • 'w 40 ,...... .5 30 10 _ U •^ •• • • • • •, M^: • Numerical O Experimental [24] 20 300 350 0 0: - L = 68.65 mm -- - RH = 96.5% CT4 50^100^150 200^250 0 81 I^.^. 0 50^100^150 200 250 300 350 Time (Hams) Time (Hams) (b) Run 2 (a) Run 1 50 - 45 40 - ,..... 35 a - . 0 0 0 0 0 • ." . •^o • • • 0 II - Cri 20 10 _ 30 _ • • • m m m mm • 0 0 0^0 . . Numerical 13 Experimental [24] 50^100^150 200^250 Time (Hums) (c) Run 3 • - .5 20 - L = 76.13 mm RH = 96% Cr 10 • Numerical 0 Experimental [24] I 0 • • • /...... .5 25 - L=70 mm RH = 96% - 00 00 0 0 0 0 0 0 30 15 40 300 350 1 250^500 Time (Hums) (d) Run 5 Figure 7.11: The Heat Flux at the Cold Plate, Comparison of Numerical Results with Experiment Data for Four Runs [24]; p 53 kg/m3, h* = 12 W/m2 K. Chapter 7. NUMERICAL RESULTS AND DISCUSSIONS ^ 82 three runs, especially for run 3 shown in Fig. 7.12(c). Fig. 7.12(d) shows that in the longer times, the predicted total moisture gain is higher than the measured data. This may be due to the decrease in vapor diffusion coefficient with liquid concentration which has not been taken into account in this model. The measured spatial distribution of the average liquid concentration was obtained by measuring the increase in mass of the different insulation layers. The comparison of the computed and measured spatial distribution of the average liquid concentration is shown in Fig. 7.13 for run 1. Fig. 13 (a) shows good agreement between the computed results with the measured data for a period up to 70 hours. The liquid distributions in this period display similar trends. The liquid produced by condensation is deposited mostly in the insulation layer adjacent to the cold plate, and is forced towards the next layer with time. Fig. 13(b) shows the comparison of computed results with experimental data at 120 and 144.5 hours respectively. Two apparent discrepancies can be observed from this plot. First the predicted liquid movement is faster than the measured. In the measurement, the liquid is trapped in the layer adjacent to the cold plate untill the layer attains a very high liquid content. This indicates that the liquid diffusivity becomes significant only when liquid fraction reaches a very high value (approximately co > 0.2 for this case). This phenomenon did not appear in Cid and Crausse's experimental work [28]. The second discrepancy is on the average liquid concentration of the layer with the exposed surface. The experimental data shows a rapid rise of average liquid concentration of that layer after a certain period while the numerical model predicts a far slower growth of average liquid concentration in the layer. The present model can not explain why the layer with the exposed surface has higher liquid content than in the inner layers. A physical explanation for this behavior given by Wijeysundera et al. [25] is the tendency for the layer with the exposed surface to accumulate liquid due to the surface tension forces between water and the last layer of fibers on the exposed 83 Chapter 7. NUMERICAL RESULTS AND DISCUSSIONS^ 350 300 : L= 62.02 mm - RH = 97% ....... V) 250 .1 • cA 0.1 3o E-1 0 • - 2 '—" 200 0 .a 0 150 0 li 100 - _ • tfl .0" 100 M • Numerical ° Experimental [241 3 0 -D o o 50^100^150 200 250 300 - oo 350 ^ to-. in . • O e le o o - 300 0 2 250 ...." o _ • 200 0 i 150 o . • Numerical 0 Experimental [241- 50^100^150 200 250 S • Numerical 0^° Experimental [24] 2 30 E-1 ^0^50^100^150^200 L= 76.13 mm E RH = 96% 300 - ^ ^ 50 _ - o - o ^ 250 300 ^ Time (Hours) ^ (c) Run 3 • s - .4 o 100 0 ^oo^ U . . 350 350 L=70 mm LW) 150 - RH = 96% • 0 o 0. (b) Run 2 (a) Run 1 4 0 • Time (Hums) Time (Horns) 200 • 0 0 - 50 F-1 0 • _ ,S 150 . O o^• - .8 200 0 GI -^o 50 o o 0 • _ 1 L = 68.65 mm RH = 96.5% ^ o^• • o s . o- o - o 0 _ 0 0 . Numerical 0 Experimental [24]250^500 Time (Hums) (d) Run 5 Figure 7.12: The Total Moisture Gain, Comparison of Numerical Results with Experiment Data for Four Runs [24]; p* = 53 kg/m3, li* = 12 W/m2 K. • Chapter 7. NUMERICAL RESULTS AND DISCUSSIONS^ 84 250 200 I 150 Predicted, 24h • o Measured [25], 24h --IV— Predicted, 48h 0^Measured [25], 48h ---A-- Predicted, 69.5h • Measured [25], 69.5h a 100 • 50 0 0.000 0.010 ^ 0.020 ^ ^ ^ ^ 0.040 0.030 0.050 0.060 Location, X (mm) (a) 300 Predicted, 120h E 250 0^Measured [25], 120h "--•^Predicted, 144.5h 0 Measured [25], 144.5h 200 - o 150 100 1-4 50 0 1^1^ 0.000 0.010 0.020^0.030 0.040 1^1^1^1^1^1^- 0.050 0.060 Location, X (mm) (b) Figure 7.13: The Liquid Distribution in Slab at Different Times, Comparison of Numerical Results with Experiment Data [24]; Run 1. p* = 53 kg/m3, h* = 12 W/m2 K, L = 62.02 mm. Chapter 7. NUMERICAL RESULTS AND DISCUSSIONS^ 85 surface. However, the reason for liquid penetrating the inner layers without becoming trapped there is unknown. The effect of gravity may be an explanation for this. Since the mechanism of moisture movement and the liquid path through the fibers in fibrous insulation are not known, theoretical simulation of liquid transport process in fibrous materials requires a more rigorous model of liquid transport coefficients. 7.5 Thermal and Transport Performances 7.5.1 The Variation of Heat Flux The variation of heat flux at cold boundary with time under fixed temperature difference between the ambient and cold plate, as well as the variation of the two components which contribute to the total heat flux, are shown in Fig. 7.14. The conductive component shows an increasing trend with times because of the increasing 'solid' thermal conductivity (which is a result of increased liquid content). The latent heat component decreases with time due to the reduction of vapor flux from ambient. However, the total heat flux shows an increasing trend with time, which suggests that the apparent thermal conductivity due to laden moisture has been increased adequately to bring this about. A dramatic increase in heat flux occurs when the liquid flows out the wet-dry interface. It attains its peak value when liquid front reaches the exposed surface, then decreases to normal values and rates. This phenomenon is complicated by the combination of evaporation and condensation which occur when liquid enters an unsaturated region. Liquid flows into the dry region forced by capillary action. When liquid encounters the warmer, unsaturated water vapor-air mixture, evaporation will occur at the liquid front. Some vapor is produced in the local space. Since the liquid amount which is forced from high liquid content region to the dry region is large, the result of liquid mixed with unsaturated water vapor is a saturated liquid-gas mixture which reaches a new Chapter 7. NUMERICAL RESULTS AND DISCUSSIONS^ I^111111111114 , 40 - N, - 86 The moment when the liquid _reaches the exposed surface ...----= CI cond+q latent^,-----\ The moment when the liquid : flows out the wet-dry interface i q cond qlatent ^■ 10 0 ^ 0 1 25^50^75 ^ 100 Fourier Number Fo Figure 7.14: The Variation of Heat Flux at Cold Boundary with Times; p* = 53 kg/m', h* = 12 W/m2 K. T = 33°C. T. = 6.8°C. Rif = 90%. - 3.0 - - 1.0 0.0 - i i Location of wet-dry interface 0.25^0.50^0.75^1.00 Dimensionless Distance Z Figure 7.15: The Variation of Equivalent Thermal Conductivity along the Insulation Slab; p* = 53 kg/m3, h* = 12 W/m2 K, T. = 33°C, T2 = 6.8°C, RH = 90%. Chapter 7. NUMERICAL RESULTS AND DISCUSSIONS^ 87 equilibrium state. The region which the liquid front occupies becomes saturated. The outflowing liquid empties its space in the wet region, which enables an increasing vapor flux into the insulation slab. The condensation rate and the condensation amount increase in the wet region. Subsequently, both the moisture laden conductive heat and latent heat increase. After the liquid front reaches the exposed surface, the entire insulation slab is saturated. The temperature and vapor density fields regain their quasi-steady behavior, and the condensation rate returns to a constant rate. The equivalent thermal conductivity at any cross section is defined by the following expression: dT qi = —Ice — 1x, . g dx (7.1) The ratio of equivalent thermal conductivity to the dry state thermal conductivity represents the effect of condensation on the thermal performance of insulation. The variation of equivalent thermal conductivity ratio along the slab at selected moments is shown in Fig. 7.15. As can be expected, the equivalent thermal conductivity ratio is unity in the dry region. It increases significantly in the wet region and attains a maximum value near the wet-dry interface. The curves of equivalent thermal conductivity decline gradually in the direction of the cold side due to the reduction of vapor flux in that direction. 7.5.2 Effect of Humidity Levels The effective thermal conductivity ratio and the total moisture gain under different humidity levels in the ambient are presented in Fig. 7.16. The term Iceffilcary in Fig. 7.16(a) represents the ratio of conventional effective thermal conductivity at cold boundary to the dry state thermal conductivity of insulation (for brevity, all thermal conductivity Chapter 7. NUMERICAL RESULTS AND DISCUSSIONS ^ 88 ratios which appeare in the following refer to the ratio at the cold boundary). Humidity levels on the exterior boundary have a significant influence on the heat flux and the total liquid concentration. This is because increasing the humidity level enhances the vapor transport, and the enhancement in vapor transfer will in turn cause an increase in heat transfer. It is also noted that a higher humidity level on the exterior boundary reduces the time to reach a quasi-steady state. As can be observed in Fig. 7.16(a), decreasing the humidity level reduces the quasi-steady wet region, thus reducing the time for liquid to flow out of the wet region. 7.5.3 Effect of Convective Heat Transfer Boundary Conditions Interesting results are obtained through the examination of the effect of the ambient convective heat transfer coefficient. Fig. 7.17 shows the variations of keffikdry and the moisture gain with time. It is noted that keff/kdry decreases with increasing heat transfer coefficient. Since increasing convective heat transfer at the exterior boundary enhances the temperature penetration inside the slab, and the saturation vapor density is mainly dependent on the temperature distribution, the wet region in the slab is reduced as a result. The reduction of the wet region leads to decreases of condensation rate. Subsequently, the moisture laden thermal conductivity is decreased, as well as the latent heat flux. Therefore the total heat flux through the insulation slab decreases with an increasing of ambient convective heat transfer coefficient. The above-mentioned results indicate the complex intercoupled nature of heat and mass transport in these types of problems. 7.5.4 Effect of Slab Thickness The effects of insulation slab thickness on thermal performance and liquid accumulation in an insulation slab are shown in Fig. 7.18. As can be expected, with a decrease of slab Chapter 7. NUMERICAL RESULTS AND DISCUSSIONS^ 89 100 (4—'s ,... 5 °JD 75 -V 0 'M 0 F0 w .-, o 30 =1 25 0 0^25^50^75^100 Fourier Number Fo (b) Figure 7.16: The Variation of Ice f f /1cdry at Cold Boundary and Total Liquid Concentration under Different Humidity Levels; p* = 53 kg/m3. Chapter 7. NUMERICAL RESULTS AND DISCUSSIONS ^ 90 100 ^ h= 8 Wm-2K ^ h=12 Wrii2 K -^h=16 Wm"2K 0 0^25^50^75^100 Fourier Number Fo (b) Figure 7.17: The Variation of Iceff/kdry at Cold Plate and Total Liquid Concentration under Different Convective Heat Transfer Boundary Conditions; p* = 53 kg/m3, Ta = 33°C, T2 = 6.8°C, L = 66 mm, RH = 90%. Chapter 7. NUMERICAL RESULTS AND DISCUSSIONS^ 0 0^25^50^75^100^125 91 150 Fourier Number Fo (b) Figure 7.18: The Variation of keff/kdry at Cold Boundary and Total Liquid Concentration under Different Slab Thickness; p* = 53 kg/m2, T. = 33°C, 7'2 = 6.8°C, h* = 12 W/m2 K, RH = 90%. Chapter 7. NUMERICAL RESULTS AND DISCUSSIONS^ g 92 2.75 4. 1::=^2.50 e 2.25 C.7 stt 2.00 25^50^75 ^ 100 Fourier Number Fo (a) 0 ^ 25^50^75 ^ 100 Fourier Number Fo (b) Figure 7.19: The Variation of ke f f kdry at Cold Boundary and Total Liquid Concentration under Different Porousities; p* = 53 kg/m', T. = 33°C, T2 = 6.8°C, L = 66 mm, h* = 12 W/m2 K, RH = 90%. Chapter 7. NUMERICAL RESULTS AND DISCUSSIONS ^ 93 thickness, both keffilcdry and the average liquid concentration increase. This is because temperature diffusion is fast in a thinner slab, and vapor saturation will occur earlier with a thinner slab. Within the same period, the liquid accumulated in a slab increases with decreasing slab thickness; thus the moisture laden conductive heat and latent heat also increase. 7.5.5 Effect of Porosity Fig. 7.19 shows that the variation in porosity does not have significant effect on the thermal performance and liquid accumulation in a porous insulation. For practical use, the porosities for most of fibrous insulation fall into a narrow range due to their highly porous characteristics. However, examining the effect of porosity can still give us the trend of energy and mass transfer under the influence of porosity. The ken ficdry and average liquid concentration increase with increasing porosity. 7.5.6 Comparison of Different Thermal Conductivity Models The apparent thermal conductivity of porous insulation in the presence of liquid is greatly dependent on the liquid distribution in an insulation. Appendix B gives several possible models for k available in literature [2], [43]. These models include bead arrangement, series arrangement, parallel arrangement and form arrangement, based on the manner of liquid location in the pores of insulation. The variation of keff/kdri, with time using different k models is shown in Fig. 7.20. The parallel arrangement model gives maximum rates of heat transfer. The bead arrangement model gives better trend-wise agreement with the measured heat flux at cold plate for all runs; it is recommended for computational use. The actual distribution of liquid in the insulation is not completely known. However, the above models give an indication of the sensitivity of the predicted quantities to the variation in these important Chapter 7. NUMERICAL RESULTS AND DISCUSSIONS ^ 25^50^75 ^ 94 100 Fourier Number Fo Figure 7.20: The Variation of keff/kary at Cold Surface Using Different k Models; p* = 53 kg/m', T. = 33°C, T2 = 6.8°C, L = 66 mm, h* = 12W/m2 K, RH = 90%. physical properties. Chapter 8 CONCLUSIONS AND RECOMMENDATIONS The following conclusions can be drawn from the analytical and numerical studies on heat and moisture transfer through a porous insulation given in the previous chapters. 8.1 Analytical Study • The analytical model developed for the heat and water vapor transfer through flatslab and round pipe insulation is valid for describing the problem of condensation in porous insulation due to thermal and vapor diffusion. • The effective thermal conductivity and other heat transfer parameters are strongly affected by condensation. The effective thermal conductivity is a maximum when condensation first occurs in the slab. For practical operating conditions its magnitude varies from about 1.5 to 15 times the dry-state value. • The effective thermal conductivity and other heat transfer parameters depend on seven independent design and operating variables. A parametric study shows that the effective thermal conductivity increases with the increasing ambient humidity level and porosity, and decrease with increasing ambient heat transfer coefficient, slab thickness and dry-state thermal conductivity of insulation. • The design curves are presented for flat-slabs. They may be used for round-pipe insulation by computing the equivalent flat-slab thickness of the insulation. 95 Chapter 8. CONCLUSIONS AND RECOMMENDATIONS^ 96 8.2 Numerical Study • The analysis in the present work has identified the four stages of heat and moisture transfer in a porous insulation. A wave-like propagation phenomenon has been observed for the important field variables such as the temperature and vapor density in the initial stage. The dynamic response and transient behavior with condensation and liquid diffusion processes are clearly shown. The numerical prediction of quasisteady behavior has verified the analytical results, and the long term behavior of energy and moisture transfer through a porous insulation has been presented. • The numerical model simulates the heat and moisture transport process through a porous insulation successfully. The predicted temperature distributions, heat transfer rates and total moisture gains have good agreements with the experimental results over various operating conditions for a time period up to 600 hours. • The present model gives reasonable prediction of liquid distribution for the situations in which liquid concentration in an insulation is low. For a complete interpretation of liquid distribution, a more accurate liquid diffusion coefficient model based on experimental work is required. • The effects of important parameters on heat and mass transfer have been investigated. The heat transfer is strongly affected by moisture distribution. The transient effective thermal conductivity increases with increasing exterior humidity level and porosity, and also increases with decreasing convective heat transfer coefficient and slab thickness. Chapter 8. CONCLUSIONS AND RECOMMENDATIONS ^ 97 8.3 Recommendations It is clear from previous comparisons and discussions that a complete simulation of heat and moisture transfer through a porous insulation will require more rigorous theoretical work, as well as more thorough and accurate experimental work. The following recommendations are suggested for a future work. • Due to lack of information on liquid diffusivity and hydraulic conductivity of fibrous insulations, more accurate experimental work on physical properties and transport performance of fibrous insulation is greatly needed. • A two-dimensional numerical model is needed in order to interpret a tendency of liquid movement along the fibers. This may be important for the cases where the pressure gradient in an insulation is low. • The present numerical model can be improved by removing some simplifications. The effects of air infiltration, gravitation, and varied pressure in an insulation would be taken into account for practical applications. • It would be interesting to extend the subject of the present study to a similar problem with two permeable boundaries which represents the common situation in building envelopes. Bibliography [1] Lotz, W. A., "Moisture problems in buildings in hot humid climates", ASHRAE J., April 1989. [2] Kaviany, M., "Principles of heat transfer in porous media", Springer-Verlag, New York, 1991. [3] Whitaker, S., "Simultaneous heat, mass and momentum transfer in porous mediva: A theory of drying", Advance in Heat Transfer, Vol. 13, 119-203, 1977. [4] Krischer, 0., " The heat, moisture, and vapor movement during drying porous materials", VDIZ., Beih. 1, 17-24, 1990. [5] Philip, J.R., and De Vries, D.A, "Moisture movement in porous materials under temperature gradients", Trans. Am. Geo. Union, 38:222-232, 1957. [6] De Vries, D.A., "Simultaneous transfer of heat and moisture in porous media", Trans. Am. Geo. Union, 39: 909-916, 1958. [7] Luikov, A.V., "Heat and mass transfer in capillary-porous bodies", Pergamon, Oxford, 1966. [8] Luikov, A.V., "System of differential equations of heat and mass transfer in capillary porous bodies (review)", Int. J. Heat Mass Transfer 18: 1-14, 1975. [9] Cary, J.W., and Taylor, S.A., "The interaction of the simultaneous diffusion of heat and water vapor", Soil Sci. Soc. Am. Proc. 26:413:416, 1962. [10] Huang., C.L.D., "Multi-phase moisture transfer in porous media subjected to temperature gradient", Int. J. Heat Mass Transfer, 22:1295-1307, 1979. [11] Eckert, E.R.G., and Faghri, M., "A general analysis of moisture migration caused by temperature differences in an unsaturated porous medium", Int. J. Heat Mass Transfer, 23: 1613-1623, 1980. [12] Jespersen, H.B., "Thermal conductivity of moisture materials and its measurements", J. Inst. Heat Vent. Engrs 1, pp. 216-222, 1953. [13] Joy, F.A., "Symposium on thermal conductivity measurements and applications of thermal insulations", ASTM STP 217 65-80, 1957. 98 Bibliography^ 99 [14] Langlais, C., Hyrien,M., and Klarsfeld,S., "Moisture migration in buildings", ASTM STP 779 192-206, 1982. [15] Bomberg, M. and Shirtliffe, C. J., "Influence of moisture gradients on heat transfer through porous building materials, thermal transmission measurement of insulation", ASTM STP 660, R.P. Tye, ed., American Society for Testing and Materials, 211-233, 1978. [16] Kumaran, M. K.," Moisture transport through glass-fibre insulation in the presence of a thermal gradient", J. of Thermal Insulation, 10:243-255, 1987. [17] Kumaran, M. K., "Comparison of simultaneous heat and moisture transport through glass-fiber and spray-cellulose insulations", J. of Thermal Insulation, 12 :6-16, 1988. [18] Langlais, L. C., Hydrien, M. and Klarsfeld, S.," Influence of moisture on heat transfer through fibrous insulating materials, thermal insulation, materials and systems for energy conservation in the 80's." ASTM STP 789, Govan, F. A., Greason, D. M. and McAllister, J. D., eds., American Society for testing and materials, 563-581,1983. [19] Langlais, C. and Klarsfeld, S., " Heat and mass transfer in fibrous insulation", J. Thermal Insulation, 8:49-80, 1984. [20] Thomas, W. C., Bal, G. P. and Onega, R. J.," Heat and mass transfer in glass fibre roof insulating materials", ASTM STP American Society for Testing and Materials, pp. 582-601, 1984. [21] Modi, D. K., and Benner, S. M.," Moisture gain of spray applied insulations and its effect on effective thermal conductivity - Part I", J. of Thermal Insulation, 8:259-277, 1985. [22] Benner, S. M. and Modi, D. K., "Moisture gain of spray applied insulations and its effect on effective thermal conductivity - Part II", J. of Thermal Insulation, 9:211223, 1986. [23] Wijeysundera, N. E., Hawlader, M. N. A. and Tan, Y. T., "Water vapour diffusion and condensation in fibrous insulations", Int. J. Heat Mass Transfer, 32(10): 18651878, 1989. [24] Wijeysundera, N. E., Hawlader, M. N. A. and Lian, S. C.," An experimental study of condensation in fiberglass insulations", Paper FE89-32 ASHRAE Far East Conference on Air Condition in Hot Climates, Kuala Lumpur, Malaysia, October 25-28, 1989. Bibliography^ 100 [25] Wijeysundera, N.E., and Hawlader, M.N.A., "Effect of condensation and liquid transport on the thermal performance of fibrous insulations", Int. J. Heat Mass Transfer, 35: 2605-2616, 1992. [26] Motakef, S. and El-Masri, M. A., " Liquid diffusion in fibrous insulation", ASME, J. Heat Transfer, 107:229-306, 1985. [27] Timusk, J., and Tenender, L.M., "Mechanism of drainage and capillary rise in glass fibre insulation", J. Thermal Insulation, 11: 231-241, 1988. [28] Cid, J., and Crausse, P., "Influence of the structural characteristic of fibrous heat insulators upon their properties of moisture transfer", J. Thermal Insulation, 14: 123-134, 1990. [29] Hedlin, C. P., "Heat transfer in a wet porous thermal insulation in a flat roof', J. of Thermal Insulation, 11:165-188, 1988. [30] Dinulescu, H.A., and Eckert, E.R.G., "Analysis of the one-dimensional moisture migration caused by temperature gradients in porous medium", Int. J. Heat Mass Transfer, 23: 1069-1078, 1980. [31] Ogniewiez, Y. and Tien, C. L., "Analysis of consensation in porous insulation" , Int. J. Heat Mass Transfer, 24:421-429, 1981. [32] Motakef, S. and El-Masri, M. A., "Simultaneous heat and mass transfer with phase change in a porous slab", Int. J. Heat Mass transfer, 29(10): 1503- 1512, 1986. [33] Vafai, K., and Sarkar, S., "Condensation effects in a fibrous insulation slab", ASME J. Heat Transfer, 108:667-675, 1986. [34] Shapiro, A. P. and Motakef. S., " Unsteady heat and mass transfer with phase change in porous slab: analytical solutions and experimental results", Int. J. Heat Mass Transfer, 33(1): 163-173, 1990. [35] Vafai, K and Witaker, S., " Simultaneous heat and mass transfer accompanied by phase change in porous insulation", Trans. ASME, J. Heat Transfer, 108:132-140, 1986. [36] Vafai, K and Tien, H. C., " A numerical investigation of phase change effects in porous materials", Int. J. Heat and Mass Transfer, 23(7): 1261-1277, 1989. [37] Tien, H.C., and Vafai, K., "A synthesis of infiltration effects on an insulation matrix", Int. J. Heat Mass Transfer, 33:1263-1280, 1990. Bibliography^ 101 [38] Tao, Y. -X., Besant, R. W. and Rezkallah, "Unsteady heat and mass transfer with phase changes in an insulation slab: frosting effects", Int. J. Heat Mass Transfer, 34(7): 1593-1603, 1991. [39] Edwards D. K., Denny, D. E. and Mills, A. F., "Transfer process: a introduction to diffusion, convection and radiation", New York, McGraw-Hill, 1979. [40] Woodside, W., "Probe for thermal conductivity measurements of dry and moist materials", Heat pipe. Air Cond. 30, pp. 163-170, 1958. [41] Crausse, P., Bacon, C. and Langlais C., " Experimental and theoretical study of simultaneous heat and moisture transfer in a fibrous insulation", J. of Thermal Insulation, 9:46-67, 1985. [42] Patankar, S.V., "Numerical heat transfer and fluid flow", New York, McGraw-Hill, 1980. [43] Batty, W.J., O'Callaghan P.W. and Probert, S.D, "Apparent thermal conductivity of glass-fiber insulant: effects of compression and moisture content", Applied Energy, 9:55:76, 1981. [44] Pallady, P.II. and Handley, P.J., "Evaluating moist air properties", Chemical Engineering, 10:, 1984. [45] Defay, R., Prigogine, I., and Bellemans, A., "Surface tension and adsorption", Wiley, New York, 1966 Appendix A Water Vapor-Air Mixture Diffusion A.1 Vapor Diffusion Coefficient in Porous Media The mass transfer of a water vapor-air mixture within the pores of a porous medium occurs by two mechanisms: ordinary diffusion and Knudsen diffusion [39]. The Knudsen diffusion refers to the case that the pores are small or the gas density is very low, therefore the molecules collide with the pore walls more frequently, and diffusion of molecules along the pore wall is described by the equations for free molecule or Knudsen flow. Porous insulation usually has a very high porosity. The vapor-air mixture is relatively dense. The vapor diffusion in porous insulation is dominated by ordinary diffusion which can be described by Fick's law. The vapor-air mixture diffusion coefficient in porous insulation, therefore, can be taken to be [39]: Dv* = 5-9-D12) T (A.1) where D12 is the binary diffusion coefficient of Pick's law; eo is the porosity; T is the tortuousity, accounting for the increased diffusion length due to tortuous paths of real pores, and for the effects of constrictions. Diffusion coefficients of gases at low pressure are almost composition independent, and increase with temperature. According to reference [23], the binary diffusion coefficient for water vapor-air mixture is given by: 102 Appendix A. Water Vapor-Air Mixture Diffusion^ 7' )1.685^ D12 = 1.97(10)-5( 255.2 ' 103 (A.2) A.2 Mass Transfer Coefficient in Ambient To determine the mass transfer coefficient between the ambient air and the insulation slab, an analysis of a slug flow over a slab surface is conducted. According to reference [39], the solution for the average heat transfer coefficient over a length of the flat surface is: 2* k^ V h* = 70/2 atL Y/2 (A.3) A similar solution can be derived for the average mass transfer coefficient over a length L: 2p*12D* m V )1/2 h* = ^a m^71/2 DI2L ( '1/4 (A.4) where V,n is the main stream velocity. Comparing the above two expressions, one could obtain: 14,^P:1)12( a* )1/2.^ h*^k D12 (A.5) Note that (a*/D:) is the Lewis number, which characterizes the ratio of thermal and mass diffusion. For quasi-steady phase, the Lewis number may be assumed to be unity. Wijeysundera et al. [23] gave a modified expression for mass transfer coefficient: *^h*^p:D12 m1047) k* ). (A.6) Appendix A. Water Vapor-Air Mixture Diffusion^ 104 A.3 Calculation of Vapor Concentration in Ambient The vapor concentration (or vapor density) of the ambient air needs to be calculated. This is done based on the saturated vapor pressure and relative humidity in the section. The equations for the determination of the saturated vapor pressure were obtained by using a mathematical fit of the pyschometric chart by Pallady and Hanley [44]. KD = —8.83(10)-13TR3 + 3.07(10)-6TR2 — 3.47(10)-3TR+ 4.4,^(A.7) F = KD(1 — 1165.7/TR) + 5, ^ (A.8) PD = 1.668(10)F,^ (A.9) where, PD is the saturated vapor partial pressure in mm of mercury, and TR is temperature in °R. The vapor partial pressure in bar is: p: = RH PD 760 . (A.10) The vapor density in ambient can be obtained from gas state equation: P: = P:RvT. (A.11) The concentration of vapor is: P: mv =^ (A.12) The air density can be obtained from gas state equation by assuming the total pressure is 1 atmospheres. Appendix B Thermal Conductivity of Moist Insulation The presence of moisture has a significant effect on apparent thermal conductivity of porous insulation. It is believed that the apparent thermal conductivity of moist insulation depends on the manner that the liquid distributes inside the insulation [43]. Four different physical models for liquid distribution in a porous medium are considered in the present study. (a) Bead arrangement. In this model liquid is located as small beads throughout the insulant. The apparent thermal conductivity can be represent by: k* (€13) = Ic:1 (1 eon (1 — esY) (B.13) where y^ (Pij — (leh 21c:1) (B.14) (b). Series arrangement. In this model the liquid is distributed in layers perpendicular to the direction of heat flow. The liquid would then have the maximum effect in inhibiting the heat flow. For this arrangement, 41q; k* (€0) = IC13 A 105 (B.15) Appendix B. Thermal Conductivity of Moist Insulation ^ 106 where Ce(ko *^a X= — — kd). €7 (B.16) (c). Parallel arrangement. In this model the liquid is located in a continuous shape parallel to the direction of heat flow. For this arrangement, k*(ep) = lc: + X, where X =— " (Ivo* — kd*). el. (d) Form arrangement. In this model a liquid film surrounds each insulant fiber or particle; the films thereby forming a honeycomb structure. The same formula as that for bead arrangement can be used for the form arrangement.
- Library Home /
- Search Collections /
- Open Collections /
- Browse Collections /
- UBC Theses and Dissertations /
- Analytical and numerical studies of heat and moisture...
Open Collections
UBC Theses and Dissertations
Featured Collection
UBC Theses and Dissertations
Analytical and numerical studies of heat and moisture transfer through porous insulation Zheng, Bi-Feng 1993
pdf
Page Metadata
Item Metadata
Title | Analytical and numerical studies of heat and moisture transfer through porous insulation |
Creator |
Zheng, Bi-Feng |
Date | 1993 |
Date Issued | 2008-08-28T22:39:39Z |
Description | This work contains both analytical and numerical studies of heat and moisture transport through a porous insulation in the presence of condensation, with impermeable, adiabatic vertical boundaries, and with one horizontal boundary facing a warm humid ambient and the other facing a cold impermeable surface. The analytical model is developed for heat and water vapor transfer in flat-slab and round-pipe thermal insulations. The model is validated by comparing its predictions with available experimental data. The effective thermal conductivity of the insulation in the presence of condensation depends on seven design and operating variables. The effect of these variables is determined by a parametric study. For practical operating conditions, the effective thermal conductivity varies from about 1.5 to 15 times the dry-state value. The computed data are presented in the form of design curves which may be used to estimate the effective thermal conductivity for flat-slab and round-pipe insulation systems. The analysis in this work quantifies the process of energy and mass transport in a porous insulation. A rigorous and fundamental formulation of heat and mass transfer in the insulation system is presented. The problem is modeled as one-dimensional, transient, multiphase flow with variable properties. Four stages in the energy and moisture transport process are identified, and they are formulated by a system of transient intercoupled equations and several thermodynamic relations using a local volume averaging technique. The numerical results are compared with experimental data for five different operating conditions and for times up to 600 hours. The model predicts the temperature distribution, heat transfer rate, the total moisture gain successfully. The predicted liquid distributions agree well with measured data for a period of up to 70 hours. The interesting effects of pertinent parameters on the energy and moisture transfer in the porous insulation are investigated. The present study, which for the first time presents a full simulation of the problem considering the mobile condensate, can be applied to other classes of problems on heat and mass transfer with phase change through a porous medium. |
Extent | 4378792 bytes |
Genre |
Thesis/Dissertation |
Type |
Text |
File Format | application/pdf |
Language | eng |
Collection |
Retrospective Theses and Dissertations, 1919-2007 |
Series | UBC Retrospective Theses Digitization Project |
Date Available | 2008-08-28 |
Provider | Vancouver : University of British Columbia Library |
Rights | For non-commercial purposes only, such as research, private study and education. Additional conditions apply, see Terms of Use https://open.library.ubc.ca/terms_of_use. |
DOI | 10.14288/1.0080848 |
URI | http://hdl.handle.net/2429/1573 |
Degree |
Master of Applied Science - MASc |
Program |
Mechanical Engineering |
Affiliation |
Applied Science, Faculty of Mechanical Engineering, Department of |
Degree Grantor | University of British Columbia |
Graduation Date | 1993-11 |
Campus |
UBCV |
Scholarly Level | Graduate |
Aggregated Source Repository | DSpace |
Download
- Media
- [if-you-see-this-DO-NOT-CLICK]
- ubc_1993_fall_zheng_bi_feng.pdf [ 4.18MB ]
- [if-you-see-this-DO-NOT-CLICK]
- Metadata
- JSON: 1.0080848.json
- JSON-LD: 1.0080848+ld.json
- RDF/XML (Pretty): 1.0080848.xml
- RDF/JSON: 1.0080848+rdf.json
- Turtle: 1.0080848+rdf-turtle.txt
- N-Triples: 1.0080848+rdf-ntriples.txt
- Original Record: 1.0080848 +original-record.json
- Full Text
- 1.0080848.txt
- Citation
- 1.0080848.ris
Full Text
Cite
Citation Scheme:
Usage Statistics
Country | Views | Downloads |
---|---|---|
China | 25 | 6 |
United States | 11 | 0 |
United Kingdom | 7 | 3 |
Japan | 7 | 0 |
India | 3 | 0 |
Chile | 3 | 0 |
Russia | 2 | 0 |
Poland | 2 | 0 |
Canada | 2 | 0 |
Philippines | 1 | 0 |
Kazakhstan | 1 | 0 |
Australia | 1 | 0 |
Mexico | 1 | 0 |
City | Views | Downloads |
---|---|---|
Beijing | 21 | 0 |
Unknown | 9 | 6 |
Tokyo | 7 | 0 |
Ashburn | 6 | 0 |
London | 4 | 3 |
Fort Worth | 2 | 0 |
Saint Petersburg | 2 | 0 |
Shenzhen | 2 | 5 |
University Park | 2 | 0 |
Uxbridge | 2 | 0 |
Xi'an | 1 | 0 |
Glen Waverley | 1 | 0 |
Toronto | 1 | 0 |
{[{ mDataHeader[type] }]} | {[{ month[type] }]} | {[{ tData[type] }]} |
Share
Embed
Customize your widget with the following options, then copy and paste the code below into the HTML
of your page to embed this item in your website.
<div id="ubcOpenCollectionsWidgetDisplay">
<script id="ubcOpenCollectionsWidget"
src="{[{embed.src}]}"
data-item="{[{embed.item}]}"
data-collection="{[{embed.collection}]}"
data-metadata="{[{embed.showMetadata}]}"
data-width="{[{embed.width}]}"
async >
</script>
</div>
Our image viewer uses the IIIF 2.0 standard.
To load this item in other compatible viewers, use this url:
http://iiif.library.ubc.ca/presentation/dsp.831.1-0080848/manifest