ARTICLE IN PRESS Building and Environment 42 (2007) 1071 1080 www.elsevier.com/locate/buildenv A simple method for estimating transient heat transfer in slab-on-ground floors Zhipeng Zhong, James E. Braun Ray W. Herrick Laboratories, School of Mechanical Engineering, Purdue University, 140 S. Intramural Drive, West Lafayette, IN 47907-2031, USA Received 19 November 2004; received in revised form 3 January 2005; accepted 20 January 2005 Abstract The problem of calculating transient heat transfer in concrete floor slabs is complicated due to ground coupling, which can require the numerical solution of two or three-dimensional transient conduction equations. This paper presents a simplified method for calculating transient slab-on-ground heat transfer that can be incorporated within hourly simulation programs. The method assumes that there are two primary one-dimensional paths for heat transfer from a ground-coupled floor slab: (1) one-dimensional heat transfer from the perimeter of the slab to the ambient and (2) one-dimensional heat transfer between the slab interior surface and a portion of the soil beneath the slab. The perimeter heat transfer is assumed to occur at quasi-steady state and is characterized in terms of a perimeter heat loss factor (F p ). Transient heat transfer within the slab and ground are modeled using a simple thermal circuit employing three nodes with an adiabatic boundary condition at a specified depth within the soil underneath the slab. Although some simulation models consider this type of two-path model, there appears to be no validation of this approach and there is no guidance for specifying perimeter heat loss factors and underfloor soil depths and node locations for the thermal circuit. In the current paper, results from detailed two-dimensional finite-element models for typical floor constructions and soil properties were used to identify (1) locations for nodes within the slab and soil, (2) correlations for soil depth as a function of soil properties associated with the underfloor adiabatic boundary condition, and (3) correlations for perimeter heat loss factor as a function of soil properties and edge insulation levels for different constructions. Transient heat transfer results from the simple model compared well with results from the finite-element program for different floor constructions, edge insulation, soil properties, locations, and times of year. r 2005 Elsevier Ltd. All rights reserved. Keywords: Transient heat transfer; Numerical method; Ground-coupled heat transfer; Perimeter heat loss 1. Introduction Ground-coupled heat transfer through concrete floor slabs is typically a significant component of the total load for heating or cooling in low-rise buildings like residential buildings, etc. It was estimated that earthcontact heat loss accounted for about 10% of the annual heat loss for average homes in the US during the early 1970s [1]. Since that time, buildings have been constructed with more insulation, better windows, and more air tight envelopes. As a result, buildings are generally more energy efficient today and ground-coupled heat Corresponding author. loss may account for 30 50% of the total heat loss for a well-built house [1 3]. Furthermore, transients associated with floor slabs can be very significant and important to consider in estimating both peak loads for sizing of equipment and total energy requirements for economic analyses. Recently, there has been considerable interest in utilizing energy storage within floor slabs and other building materials to achieve load shifting and reduction in on-peak cooling requirements and costs (e.g., [4]). The load shifting is accomplished simply by manipulating zone setpoint temperatures over time within the comfort zone. The motivation for the work presented in this paper was the need for a simple model for predicting transient heat transfer in ground-coupled 0360-1323/$ - see front matter r 2005 Elsevier Ltd. All rights reserved. doi:10.1016/j.buildenv.2005.01.030
1072 ARTICLE IN PRESS Z. Zhong, J.E. Braun / Building and Environment 42 (2007) 1071 1080 Nomenclatures a absorptance for solar radiation A f total floor surface area, m 2 A s amplitude of the ground surface temperature wave c p specific thermal capacity, J/kg-K C node thermal capacitance, J/K; D* ratio of the soil depth to the soil depth associated with the baseline (45 cm) e transfer function coefficient for previous outputs, W/m 2 F p perimeter heat loss factor, W/m-K h convective heat transfer coefficient, W/m 2 -K k thermal conductivity, w/m-k k* ratio of the thermal conductivity to the baseline thermal conductivity of 1.73 w/m-k k v vegetation shade factor L length of the internal zone, m P external perimeter of the slab, m R node thermal resistance, K/W R in surface unit thermal resistance, K/W s transfer function coefficient for previous temperatures, W/m 2 -K q 00 f heat flux, W/m 2 _Q T z T z T sol air T sol air T g T m y Subscripts total heat transfer rate, W internal zone air temperature, 1C 24-h average for T z, 1C sol-air temperature, 1C 24-h average for T sol air, 1C far depth soil temperature, 1C annual average temperature, 1C soil depth, m t time in hourly series 1,2,3,4 node series 4 HWC 4-in heavy weight concrete layer 6HWC 6-in heavy weight concrete layer gravel gravel layer soil soil layer FE finite-element method Greek letters a thermal diffusivity, m 2 /s j phase angle r density, kg/m 3 t time, s floor slabs that could be integrated within an hourly simulation program and used to assess the potential for load shifting opportunities associated with building energy storage in small commercial buildings. There are well-established models for steady-state, ground-coupled heat transfer from basements and slab floors. For instance, ASHRAE [5] presents a steadystate model for heat loss from ground-coupled floor slabs that are typically used for modeling in residential applications. The model assumes that the predominate path for heat loss from the slab occurs between the slab perimeter and the ambient. The floor heat transfer is characterized in terms of a perimeter heat loss factor (F p ). ASHRAE gives typical values of F p for a limited number of foundation insulation configurations and climatic conditions. The method and heat loss factors were developed from measurements conducted for specific foundation sizes and soil properties under steady-state heat transfer conditions. No general approach has been presented for estimating perimeter heat loss factors from material and soil properties. EN ISO 13370 [6] presents a similar method for estimating steady-state ground floor heat flux. Medved and Cerne [7] used this general formulation and developed simplified method to calculate parameters of the model for different sizes and shapes of buildings with and without basements with a variety of properties. Krarti et al. [8] developed the ITPE (Interzone Temperature Profile Estimation) method that combines analytical and numerical techniques to obtain two- and three-dimensional solutions of the heat conduction equation for slab-on-grade floors and basements [9]. The ITPE method provides a general solution for the temperature field within soil under steady conditions for a slab-on-grade and slab-on-ground floor. The ITPE formalism has been applied to develop a simplified method for seasonal foundation heat loss [10], evaluate the thermal bridge effects of floor foundation [11], optimize insulation R-value by life-cycle cost analysis [12,13] and develop simplified tool for the calculation of foundation heat gain [12,13]. For analyzing the use of building thermal mass, it is very important to consider transient heat transfer within the floor slab and underfloor materials [14]. One approach involves solution of two- or three-dimensional transient conduction equations using finite-difference or finite-element methods. The solution generally requires the numerical solution of a large number of coupled equations. Because of the computation required, studies have been performed to develop fast computational schemes or modeling simplifications for this problem. For instance, Seem et al. [15] developed a method for calculating transfer functions in multi-dimensional conduction heat transfer problems that is useful for
ARTICLE IN PRESS Z. Zhong, J.E. Braun / Building and Environment 42 (2007) 1071 1080 1073 modeling ground-coupled floor slabs. The approach utilizes the following steps: (1) convert a finite-difference representation of the partial differential equations to a state-space formation, (2) utilize the analytical solution to the state-space formation with a ramp or step input to determine a transfer function representation, (3) utilize model-order-reduction techniques to obtain a simplified transfer function representation. The transfer function representation is determined prior to or at the beginning of a simulation and is then used at each timestep to calculate energy flows at the boundaries of the domain (e.g., at the floor-interior interface). For annual simulations, the transfer function representation requires much less computation than that associated with numerical solutions to the partial differential equations at each timestep. Adjali et al. [16] performed numerical simulations for a ground-coupled floor slab and compared predictions with measurements. Two- and three-dimensional conductive heat transfer models based on a finite-volume method for the floor and soil were added to APACHE, a whole building thermal simulation program widely used in Europe. In-situ measurements were made at the Cardiff School of Engineering, University of Wales, over a period of 18 months. Temperatures (internal and external air, at different depths in the slab and the soil down to 4 m beneath the floor surface), heat fluxes through the concrete floor slab, the thermal conductivity of the ground and the water table and soil moisture content were measured. The results showed that the purely conductive earth-contact model accurately predicted the thermal behavior beneath the building. A two-dimensional analysis was found to be adequate. Three-dimensional effects were only important near the corners of the slab. Even with simplifications, the computational requirements for solving two or three-dimensional heat transfer within floor slabs and ground soil is very significant. Simulation tools available for estimating cooling and heating loads in commercial and residential buildings typically do not employ detailed multi-dimensional models for ground-coupled heat transfer. For instance, TRNSYS [17] considers one-dimensional transient conduction in the floor and materials immediately beneath the floor. A portion of the soil can be considered with a specified underfloor boundary condition. However, no guidance is given for the amount of soil to be included and the temperature of underfloor boundary condition to employ. Some simulation models (e.g., [18]) consider two paths for heat transfer from a ground-coupled slab floor: (1) steady-state heat loss from the perimeter to the ambient and (2) one-dimensional transient conduction in the floor and materials immediately beneath the floor with an adiabatic underfloor boundary condition. However, no literature could be found that validates this modeling approach. Furthermore, it s not obvious how to specify perimeter heat loss factors and underfloor soil depths and the number and location of nodes for the simplified thermal circuit. The current paper starts with the general two-path modeling approach for transient slab heat transfer and develops heuristics and correlations for general application of the method. In particular, results from detailed two-dimensional finite-element models for typical floor constructions and soil properties were used to identify (1) locations for nodes within the slab and soil, (2) correlations for soil depth as a function of soil properties associated with the underfloor adiabatic boundary condition, and (3) correlations for perimeter heat loss factor as a function of soil properties and edge insulation levels for different constructions. 2. Methodology 2.1. Case study descriptions A number of ground-coupled floor slabs were considered in the model development and evaluation study. The basic geometry and some of the parameters are shown in Fig. 1. The finite-element modeling was applied to a two-dimensional slice including the slab floor and surrounding soil. The construction and dimensions are typical of a small commercial building. Many of the results were determined using a slab length Fig. 1. Schematic of the slab-on-ground geometry.
1074 ARTICLE IN PRESS Z. Zhong, J.E. Braun / Building and Environment 42 (2007) 1071 1080 of 20 m. It was determined that the spatially averaged heat flux at the indoor floor surface is insensitive to slab length for values greater than about 10 m. The dimensions of the soil (depth of 10 m and width of 50 m) were chosen so that the bottom and side boundaries were nearly adiabatic corresponding to the far field assumption [16]. Due to symmetry, only half of the floor and soil shown in Fig. 1 was modeled. The internal zone boundary condition was convection to a zone temperature that varies within a relatively narrow range. The boundary condition at the soil surface considered both convection and radiation using an equivalent sol air temperature [5]. Typical convective heat transfer coefficients and a solar surface absorptance for summer conditions were specified for the upper surfaces of the floor and ground. Two widely used constructions were selected for floor slabs from the Builder s Foundation Handbook [19]. Configuration A used a 15 cm (6 in) heavy weight concrete slab laying directly on ground soil. Configuration B had a 10 cm (4 in) heavy weight concrete slab with 4 in of sand and pea gravel above the ground soil. A heavy weight concrete foundation (20 cm in width and 1 m in depth) with various foundation insulation levels was considered. For a given floor construction and set of boundary conditions, the primary factors that affect slab heat flux are the level of insulation on the foundation and the soil properties. If insulation is employed on the foundation, then typically about 1 in (2.54 cm) to 2 in (5.1 cm) of styrofoam is attached to the outer foundation wall below ground. In this study, outer foundation insulation was assumed to extend from the top surface of the ground to a 1 m depth. Base case results were determined for insulation having a thermal resistance corresponding to 1 in (2.54 cm) of styrofoam insulation (R ¼ 0:95 K m 2 =WorR ¼ 5:4ft 2 21F h=btu). A range of different insulation levels were considered (R ¼ 1:425 K m 2 =W; R ¼ 0:95 K m 2 =W; R ¼ 0:475 K m 2 =Wand R ¼ 0K m 2 =W) in developing correlations for perimeter heat loss factors. Table 1 gives thermal properties of floor, underfloor materials, and soil used for base case results. Soil properties can vary considerably according to soil type and moisture content. It is extremely difficult to consider the time dependence of soil properties due to changes in moisture content. Typically, average properties are assumed for a particular site. The base case soil properties in Table 1 are representative of saturated Table 1 Base case thermal properties of floor and ground materials k (w/m-k) r (kg/m 3 ) c p (J/kg-K) a (m 2 /s) Heavy weight concrete 1.80 2245 837 9.58E-7 Pea gravel 2.52 1800 1000 1.40E-6 Ground soil 1.73 1844 1090 8.61E-7 soil of medium sand with fine gravel [20]. In addition, soil properties were varied over a wide range to develop correlations for perimeter heat loss and soil depth for the simplified model. For steady-state conduction, the only soil property of interest is the thermal conductivity. For transient conduction problems, thermal diffusivity (ratio of thermal conductivity to the product of the density and specific heat) is the relevant parameter. Thermal diffusivity, a; is a measure of the ability of the material to diffuse thermal energy through conduction. The higher the thermal diffusivity the faster a material will respond to changes in the boundary conditions. In addition to the base case, soil thermal conductivity was varied in the range of 0.6 3.5 W/m-K. This corresponds to thermal diffusivities in the range of 2.99 E-6 1.74E-6 m 2 /s. For the simulation cases considered in this study, the zone and ambient sol air temperature boundary conditions were varied using 24-h periodic functions of the following form. TðtÞ ¼A 0 þ A 1 sin p 12 t þ j 1 þ A 3 sin p 3 t þ j 3 þ A 2 sin p 6 t þ j 2. ð1þ For the ambient sol air temperatures, different sets of coefficients were determined through regression using hourly weather data from different locations and different times of year. Hourly ambient temperatures and horizontal radiation were used along with the solar absorptance and convection heat transfer coefficient given in Fig. 1 to compute sol-air temperatures for the different cases. For instance, Fig. 2a shows a sample 24- h variation in sol-air temperature for Arcata, CA (California climate zone 1) on August 1st based on TMY2 data. The zone temperature variation was chosen based upon typical results associated with conventional night setup control in small commercial buildings [14], as shown in Fig. 2b. In this study, Fig. 2a and b were used as the base case to build the correlations only. Ambient sol air temperatures and zone temperatures for small commercial buildings from other locations in different seasons under various controls were used in validating the model. 2.2. Finite-element model A commercially available finite-element program [21] was employed to solve the two-dimensional transient heat conduction problem for the geometry shown in Fig. 1. The slab floor slab and surrounding soil, footer and concrete wall structure were modeled using approximately 3000 nodes. Fig. 3 depicts the nonuniform mesh arrangement employed. The element sizes were smaller in areas of larger temperature gradients, such as near the top surface and at the edge of the floor
ARTICLE IN PRESS Z. Zhong, J.E. Braun / Building and Environment 42 (2007) 1071 1080 1075 slab. Adiabatic boundary conditions were specified at 15 m from the edge of the slab and at a depth of 10 m from the top surface. The internal zone (slab and wall) and ground surfaces were subjected to convective heat transfer using the boundary conditions described in the previous section. For the transient modeling, the initial temperatures of the floor slab and soil directly beneath the slab were chosen as the average of the zone air temperatures over Tsol-air (C) (a) Tz (C) 60.0 50.0 40.0 30.0 20.0 10.0 Tsol-air for Arcata, CA Aug-1st 0.0 0:00 6:00 12:00 18:00 0:00 24.5 24.0 23.5 23.0 22.5 22.0 Tz for Arcata, CA Aug-1st 21.5 0:00 6:00 12:00 18:00 0:00 (b) Fig. 2. Example variations in ambient sol air and zone temperatures. (a) Ambient sol air temperature (b) Zone temperature. the 24-h period. The initial temperatures of the surrounding soil were determined using the model of Baggs given in Eq. (2) [22], which was developed for naturally occurring ground temperatures away from building structures. Tðy; tþ ¼T m 1:07 k v A s Expð 0:00031552ya 0:5 Þ cos 2p 365 ðt t 0 þ 0:018335ya 0:5 Þ, ð2þ where the vegetation shade factor, k v ; was taken to be 1.0 for bare ground like a parking lot or side walk outside a retail store. A 1-h time step was used for all simulations. The model was run for many identical days in a row until the solution approached a steady-periodic condition. This required simulation lengths of up to 12 000 h. 2.3. Simplified ground-coupled floor modeling For estimating zone cooling and heating loads, the quantity of interest for modeling a ground-coupled floor is the heat transfer rate at the top surface of the floor. The simplified modeling approach considers two primary one-dimensional paths for heat transfer from a groundcoupled floor slab: (1) one-dimensional steady-state heat transfer from the perimeter of the slab to the ambient and (2) one-dimensional transient heat transfer between the slab interior surface and a portion of the soil beneath the slab. At any time t; the total heat transfer rate from the zone air to the floor is determined as _Q t ¼ F p Pð T z T sol air ÞþA f q 00 f;t (3) where F p is the slab perimeter heat loss factor, P is the slab perimeter, T z and T sol air are 24-h average zone and sol air temperatures (updated continuously), A f is the total floor surface area, and q 00 f is the heat flux (heat transfer rate per unit area) at the floor surface for a onedimensional transient analysis. The one-dimensional heat flux is essentially the heat flux that occurs at Fig. 3. Schematic of the two-dimensional finite-element modeling mesh.
1076 ARTICLE IN PRESS Z. Zhong, J.E. Braun / Building and Environment 42 (2007) 1071 1080 the centerline of the slab (i.e., the left-most boundary of Fig. 3). The one-dimensional transient analysis uses a simple thermal circuit employing 3 nodes with an adiabatic boundary condition at a specified depth within the soil underneath the slab as depicted in Fig. 4. The number and location of for the nodes was determined through analysis of and comparisons with finite-element modeling results. The resulting model can be represented with three continuous, linear, time-invariant differential equations as: C 1 dt 1 dg ¼ T z T 1 R 1 þ T 2 T 1 R 2 C 2 dt 2 dg ¼ T 1 T 2 R 2 þ T 3 T 2 R 3 dt 3 C 3 dg ¼ T 2 T 3 þ T g T 3 ð4þ R 3 R 4 The three differential equations can be put in matrix form as 2 6 4 dt 1 dt dt 2 dt dt 3 dt 3 2 7 5 ¼ 6 4 1 R 1 C 1 1 R 2 C 1 1 R 2 C 1 0 2 6 4 3 7 1 R 2 C 2 1 R 2 C 2 1 1 R 3 C 2 R 3 C 2 7 5 1 0 R 3 C 3 1 R 3 C 3 3 2 32 3 1 R 1 C 1 0 0 T z 7 5 þ 6 4 0 0 076 0 7 54 5. ð5þ 0 0 0 T g T 1 T 2 T 3 Eq. (5) is a state-space representation for the simplified one-dimensional conduction heat transfer problem. It is solved analytically and converted to a transfer function using the method of Seem et al. [15]. The transfer function solution gives the current one-dimensional heat flux at the floor surface in terms of current and past hourly values of the zone temperature and past hourly values of the heat flux. q 00 f;t ¼ s 0T z;t þ s 1 T z;t 1 þ s 2 T z;t 2 þ s 3 T z;t 3 þ e 1 q 00 f;t 1 þ e 2 q 00 f;t 2 þ e 3q 00 f;t 3. The s and e transfer function coefficients are determined from the solution to the state-space problem with linear changes in zone temperature over each hour as outlined by Seem et al. [15]. The transfer function representation is determined prior to a simulation and is then used at each hourly timestep to calculate the floor heat transfer. The initial temperatures of the floor slab and soil directly beneath the slab were chosen as the average of the zone air temperatures over the 24-h period. A 1-h time step was used for all cases and the model was run to a steady-periodic condition for each day type considered. 2.4. Process for correlating perimeter heat loss factor and soil depth ð6þ In order to apply the simple model described in the last section, it is necessary to specify a perimeter heat loss factor and depth of soil under the floor. The perimeter heat loss depends on the floor construction, amount of foundation insulation, and soil properties. The appropriate depth of soil to include in the analysis depends on the floor materials and soil properties. Results from finite-element modeling were used to develop correlations for these two parameters. The following steps were followed to develop and evaluate the correlations and simplified model. Fig. 4. Simplified 3-node models for a one-dimensional groundcoupled floors. (a) Floor Configuration A. (b) Floor Configuration B. 1. The finite-element model was used to generate onedimensional heat fluxes at the adiabatic center line of the floor slab (i.e., the left-most boundary in Fig. 3) for floor configurations A and B with different soil properties and the baseline ambient and zone temperature variations. 2. A simplified one-dimensional transient model was also developed for each configuration and soil type in terms of the unknown soil depth. 3. For a given floor configuration and soil properties, errors in predictions of hourly centerline floor surface heat fluxes were minimized with respect to soil depth
ARTICLE IN PRESS Z. Zhong, J.E. Braun / Building and Environment 42 (2007) 1071 1080 1077 for the simplified model. The problem involved minimizing the following cost function: vffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi ux 24 2 J ¼ t (7) t¼1 q 00 f;t q00 FE;t with respect to the soil depth, D; where q 00 f;t was determined from the simplified model and q 00 FE;t was from the finite-element model. 4. Correlations of optimized soil depth D in terms of soil conductivity, k; were developed for estimating centerline floor heat transfer. 5. The two-dimensional finite-element model was used to generate transient heat transfer results over the entire domain depicted in Fig. 3 for floor configurations A and B with different soil properties and the baseline ambient and zone temperature variations. The total heat transfer rates at the top floor surface were determined for each hour from the twodimensional results. 6. For a given floor configuration, foundation insulation, and soil properties, errors in predictions of hourly floor surface heat transfer rates were minimized with respect to the perimeter heat loss factor for the simplified model. The problem involved minimizing the following cost function: vffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi ux 24 J ¼ t _Q f;t Q _ 2 FE;t (8) t¼1 with respect to perimeter heat loss factor, F p. 7. Correlations of optimized perimeter heat loss factor in terms of soil conductivity and foundation insulation R value were developed for the different floor configurations. 8. Predictions from the resulting simplified model were compared with two-dimensional finite-element modeling results for a range of ambient conditions (winter and summer in different locations), indoor temperature variations, soil conductivities, and foundation insulation levels. 3. Results and discussion 3.1. Floor configuration A: 4 heavy-weight concrete on soil Fig. 5 gives an example of the effect of soil depth on the performance of the simplified one-dimensional heat flux model (Fig. 4a) for a soil thermal conductivity of 1.73 w/m-k. The optimal soil depth to use in this simplified model is about 45 cm for this case. Fig. 6 compares the centerline heat flux determined with the one-dimensional finite-element model and simplified model with the optimized soil depth. Mean Square Error (W/m 2 ) 0.30 0.28 0.26 0.24 0.22 0.20 0.18 10 20 30 40 50 60 70 Depth (cm) Fig. 5. Example optimization of soil depth for centerline heat flux. Heat Flux (W/m 2 ) 12.0 8.0 4.0-8.0-12.0 k*=1, R*=0 Simplified (CTF) 0.0 0:00 6:00 12:00 18:00 0:00-4.0 Fig. 6. Comparison of centerline heat flux for finite-element and simplified models (Configuration A, base case). Optimized soil depths for a range of soil thermal conductivities were determined and used to develop a correlation. For floors without gravel below the slab (Configuration A), the resulting correlation is D ¼ 0:0489 k 2 þ 0:5692 k þ 0:4709, (9) where k is the ratio of the thermal conductivity to the baseline thermal conductivity of 1.73 w/m-k and D is the ratio of the soil depth to the soil depth associated with the baseline conductivity (45 cm). Fig. 7 shows comparisons between the curve fit and optimized soil depths. The overall simplified model utilizes the simplified centerline heat flux model with the soil depth correlation and an edge heat loss model as given in Eq. (3). Errors between model predictions from Eq. (3) and twodimensional finite-element results were minimized to obtain perimeter heat loss factors for various edge insulation levels and soil conductivities. Fig. 8 gives an example of hourly comparisons for average floor surface heat flux obtained using an optimized soil depth and perimeter heat loss factor for baseline parameters and conditions. The simplified model is able to accurately characterize the transient heat transfer at the floor surface.
1078 1.60 ARTICLE IN PRESS Z. Zhong, J.E. Braun / Building and Environment 42 (2007) 1071 1080 4.00 D* 1.40 1.20 1.00 0.80 0.60 Curve Fit 0.40 0.00 0.50 1.00 1.50 2.00 k* Fig. 7. Curve fit for optimized soil depth for centerline heat flux (Configuration A). Fp (W/m-K) 3.50 3.00 2.50 2.00 1.50 1.00 L=10m L=7.5m L=5.0m L=2.5m 0.50 0.00 0.50 1.00 1.50 2.00 k* Fig. 9. Effect of soil conductivity and half-floor length on optimized perimeter heat loss factor (Configuration A, no edge insulation). Average Heat Flux(W/m2) 10.0 8.0 6.0 4.0 2.0-6.0-8.0-10.0 k*=1, R*=0 Simplified (CTF+Fp) 0.0 0:00 6:00 12:00 18:00 0:00-2.0-4.0 Fig. 8. Comparison of average floor surface heat flow flux for finiteelement and simplified models (Configuration A, base case). Fp (W/m-K) 4.00 3.50 3.00 2.50 2.00 1.50 1.00 0.50 R*=0.0 R*=0.5 R*=1.0 R*=1.5 0.00 0.00 0.50 1.00 1.50 2.00 k* Fig. 10. Curve-fit results for edge heat loss factor (Configuration A). The perimeter heat loss factor depends strongly on the edge insulation and soil conductivity and weakly on the dimensions of the zone. Fig. 9 shows the effect of dimensionless soil conductivity and half-floor length (distance from centerline to edge) on the optimized perimeter heat loss factor for baseline parameters and conditions, except with no foundation insulation. Without edge insulation, the impact of floor length is negligible for overall floor lengths greater than about 10 m. Adding edge insulation leads to even lower impact of the zone dimension on the edge heat loss. Therefore, it is reasonable to neglect the effect of zone dimensions on edge heat loss for most buildings. Fig. 9 indicates that perimeter heat loss factor is linear with soil conductivity. However, it also depends upon the edge insulation. The following correlation was developed for floor configuration A using two-dimensional finite-element results for the baseline geometry and conditions over a range of soil conductivities and edge insulation levels. F p ¼ ak þ b, (10) where the slope a and intercept b are correlated as a ¼ 0:1935R 3 þ 0:4594R 2 0:5029R þ 1:4911, b ¼ 0:6539R 3 þ 1:5876R 2 1:4163R þ 0:5628, and where R is the ratio of the edge insulation R-value to the base case edge insulation R ¼ 0.95 W/m 2 -K. Fig. 10 shows comparisons between the optimized perimeter heat loss factors and values determined from the correlation of Eq. (10). 3.2. Floor configuration B: 10 cm (4 in) HWC slab with 4 inches of gravel Optimized soil depths for a range of soil thermal conductivities were determined for configuration B (Fig. 5b) and used to develop the correlation in Eq. (11). Fig. 11 shows comparisons between the curve fit and optimized soil depths. D ¼ 0:1106 k 2 þ 0:7305 k þ 0:3809. (11) For floor configuration B, the linear form given in Eq. (10) was used to correlate perimeter heat loss factor
ARTICLE IN PRESS Z. Zhong, J.E. Braun / Building and Environment 42 (2007) 1071 1080 1079 D* 1.60 1.40 1.20 1.00 0.80 0.60 Curve Fit 0.40 0.00 0.50 1.00 1.50 2.00 k* Fig. 11. Curve fit for optimized soil depth for centerline heat flux (Configuration B). Average Heat Flux (W/m 2 ) 6.0 5.0 4.0 3.0 2.0 1.0-2.0 Simplified (CTF+Fp) k*=1.15, R*=0.70 0.0 0:00 6:00 12:00 18:00 0:00-1.0 Fig. 13. Example simplified model accuracy for New York City in summer (Configuration A). Fp (W/m-K) 4.00 3.50 3.00 2.50 2.00 1.50 1.00 0.50 0.00 0.00 0.50 1.00 1.50 2.00 k* R*=0.0 R*=0.5 R*=1.0 R*=1.5 Fig. 12. Curve-fit results for edge heat loss factor (Configuration B). Average Heat Flux (W/m 2 ) 10.0 8.0 6.0 4.0-4.0-6.0-8.0-10.0 Simplified (CTF+Fp) k*=0.85, R*=1.20 2.0 0.0 0:00-2.0 6:00 12:00 18:00 0:00 Fig. 14. Example simplified model accuracy for Chicago in winter (Configuration B). with soil conductivity and edge insulation. The resulting expressions for the slope a and intercept b in Eq. (10) are a ¼ 0:118R 3 þ 0:3928R 2 0:4419R þ 1:4463, b ¼ 0:2937R 3 þ 1:0164R 2 1:2702R þ 0:5897. Fig. 12 shows comparisons between the optimized perimeter heat loss factors and values determined from the correlation of Eq. (10). 3.3. Validation The correlations were developed using a single profile for ambient temperature variations with an associated deep ground temperature that were representative of summer in California and a single profile for zone temperature variation that was associated with night ventilation precooling. However, the resulting model was tested over a wide range of ambient and zone temperature conditions, including summer and winter conditions in a variety of locations. Fig. 13 shows example results in New York City in summer (August 1st) under night-ventilation control for floor configuration A. Fig. 14 gives results in Chicago for winter (February 1st) with conventional night setup control for floor configuration B. These results are typical of the accuracy of the simplified model in predicting slab heat transfer rate as compared with the finite-element results. Part of the explanation for the good agreement between the simplified model and the numerical simulations is the fact that the correlations for the simplified model were developed using the same software and assumptions used for the validation. Agreement with measurements would not be expected to be as good because of un-modeled effects, such as time variation in soil moisture content. 4. Conclusions A simple model was developed for ground-coupled concrete floor slabs that can predict spatially averaged transient heat flux at the interior surface and could be integrated with existing building load simulation tools. The model incorporates correlations for the effects of ground soil properties and edge insulation that were
1080 ARTICLE IN PRESS Z. Zhong, J.E. Braun / Building and Environment 42 (2007) 1071 1080 developed for two common floor configurations. Hourly heat flux predictions compared very well with predictions from a two-dimensional finite-element program for these two geometries. References [1] Labs K, Carmody J, Sterling R, Shen L, Huang Y, Parker D. Building foundation design handbook. ORNL Report Sub/86-72143/1, Oak Ridge National Laboratory. 1988. [2] Deru M. A model for ground-coupled heat and moisture transfer from buildings, NREL/TP-550-33954, National Renewable Energy Laboratory Technical Report. 2003. [3] Claesson J, Hagentoft CE. Heat loss to the ground from building I. General theory. Building and Environment, 1991; 26(2):195 208. [4] Braun JE. Load Control using building thermal mass. ASME Journal of Solar Energy Engineering 2003;125:292 301. [5] ASHRAE, ASHRAE handbook fundamentals. Atlanta: American Society of Heating, Refrigerating and Air-Conditioning Engineers, Inc; 2001. pp. 28.10 28.13. 2001. [6] European Committee for Standardization. Thermal Performance of Buildings Heat Transfer Via the Ground Calculation Methods. European Standard EN ISO 13370. 1998. [7] Medved S, Cerne B. A simplified method for calculating heat losses to the ground according to the EN ISO 13370 standard. 2002; 34: 523 528. [8] Krarti M, Claridge DE, Kreider JF. ITPE method applications to time-varying two-dimensional ground-coupling problems. Internatinal Journal of Heat Mass Transfer 1988;31:1899 911. [9] Claridge DE, Krarti M, Kreider JF. Energy calculations for basements, slabs, and crawl Spaces, Final Report, ASHRAE 666- TR, Atlanta: American Society of Heating, Refrigerating and Air- Conditioning Engineers, Inc. 1993. [10] Krarti M, Choi S. Simplified method for foundation heat loss calculation. ASHRAE Transaction 1996;102 part 1. [11] Al-Anzi A, Krarti M. Evaluation of the thermal bridging effects on the thermal performance of slab-on-grade floor foundation. 1997. vol. 103. part 1. [12] Chuangchid P, Krarti M. Economic analysis of foundation heat gain for refrigerated warehouses. ASHRAE Transaction 2000;106 part 2. [13] Chuangchid P, Krarti M. Parametric analysis and development of a design tool for foundation heat gain for coolers. ASHRAE Transaction 2000;106 part 2. [14] Braun JE, Zhong Z. Development and evaluation of a night ventilation precooling algorithm. To appear in International Journal of Heating, Ventilating, Air-Conditioning and Refrigeration Research 2005. [15] Seem JE, Klein SA, Beckman WA, Mitchell JW. Transfer function for efficient calculation of multidimensional transient heat transfer. ASME Journal of Heat Transfer 1989;111:5 12. [16] Adjali MH, Davies M, Rees SW, Litter J. Temperature in and under a slab-on-ground floor: two- and three-dimensional numerical simulations and comparison with experiment data. Building and Environment 2000;35:655 62. [17] Solar Energy Laboratory. TRNSYS: a transient system simulation program (reference manual). Madison, WI: University of Wisconsin/Madison Press; 2000. [18] Balcomb JD, Crowder RS. (NREL), ENERGY-10, A design tool computer program for buildings. Proceedings of Solar 95 (20th Passive Solar Conference), American Solar Energy Society, July 15 20, 1995, Minneapolis, MN. 1995. [19] Carmody J, Christian J, Labs K. Builder s foundation hand book. ORNL/CON-295, ORNL Report Sub/86-72143/1, Oak Ridge National Laboratory. 1991. [20] Henk JL, van Gelder GJ. In Situ measurement of ground thermal conductivity: a Dutch perspective. ASHRAE Transactions 2002; 108:0 part 1. [21] Klein SA, Beckman WA, Myers GE. FEHT A finite element analysis program. F-Chart Software. 1995. [22] Popiel CO, Wojtkowiak J, Biernacka B. Measurement of temperature distribution in ground. Experimental thermal and fluid science 2001;25:301 9.