Investigation of Effect of Rock Storage System Parameters on Thermal Cooling Performance
Alex Xavery Matofali^{1}, Estomih S. Massawe^{2}
^{1}Faculty of Science, Department of Biometry and Mathematics, Sokoine University of Agriculture, Morogoro, Tanzania
^{2}Mathematics Department, University of Dar es Salaam, Dar es Salaam, Tanzania
Email address:
To cite this article:
Alex Xavery Matofali, Estomih S. Massawe. Investigation of Effect of Rock Storage System Parameters on Thermal Cooling Performance. Applied and Computational Mathematics. Vol. 5, No. 1, 2016, pp. 1017. doi: 10.11648/j.acm.20160501.12
Abstract: We investigate the effects of key parameters of the rock bed system on thermal cooling performance of the system after a fixed time of operation. The method of solving the mathematical model uses a semidiscretization finite difference approximation for discretizing space in solid problem domain. A finite element approximation is used in the fluid problem domain. Graphical results on the effects of parameter variation on damping and time delay on the peaking of the outlet air temperature through the bed are presented and discussed.
Keywords: Rock Beds, Effects, Thermal Storage Systems, Parameters
1. Introduction
Passive cooling systems using building fabric storage have been used and researched over years, as means of providing building cooling using much less energy than electricity powered individual airconditioning unit. The need to use sustainable energy technology is of great interest. Using renewable energy and energy efficiency in buildings is technoeconomically feasible, and significantly cheaper than building new power stations. A longterm solution in the form of new lowcost technology that relies only on renewable sources of energy and on locally available materials is needed. Sustainable energy sources researches are continuing to address the concern over climate changes, pollution and nonrenewable sources (Anderson et al., 2015).
The use of thermal energy storage (TES) plays an important role in building energy conservation initiatives, which are greatly enhanced by incorporation of latent heat storage in building products. Thermal energy storage deals with selecting of media/devices that absorb and store heat/cold during peak power operation and releases the same during reduced power operation. Building fabric or packed beds (rock pile, pebble bed and rock bed) materials are some of the thermal storage devices (Ravikumar and Srinivasan, 2005). However, storage and recovery of thermal energy must be done efficiently to achieve high capacity factors and low levelized cost of electricity. As described in the review of Kuravi et al., (2013) thermal energy technology should meet several requirements; high energy density, good heat transfer between the heat transfer fluid and solid storage media, stability of the storage medium, low thermal losses, low cost and reversibility through many charging and discharging cycles.
The thermal energy storage capacity of buildings if increased can increase human comfort by decreasing the frequency of internal air temperature swings so that indoor air temperature is closer to the desired temperature for a longer period of time (Dincer et al., 1997). TES systems provide the potential to attain energy savings, which in turn reduce the environmental impact related to nonrenewable energy use. In fact, these systems can reduce time or rate mismatch that is often found between energy supply and energy demand (Ataer, 2006; Pasupathy and Velraj, 2006; Manohar and Adeyanju, 2009).
The thermal mass is that it can be used in conjunction with night ventilation of a building to provide passive cooling. Outside air is circulated through the building where it comes into contact with and cools the building fabric. The cooling that is stored in the building fabric is then available to offset the heat gains of the following day and keep temperatures within comfort limits (Barnard, 2006). A comparison of thermal energy storage design is given by Li et al. (2012). TES can be done with sensible heat storage systems or latent heat storage. The present study explores sensible thermal energy storage and the solid arrangement studied here is to store the heat/cool in a packed bed.
According to Coutier and Farber (1982) packed bed (rock bed) represents the most sensible energy storage unit for air based systems. A packed bed is a volume of porous media obtained by packing particles of selected materials into a container. A packed bed storage system consists of loosely packed solid materials through which the heat transport fluid is circulated. Heated fluid (usually air) flows from solar collectors into a bed of graded particles from top to bottom in which thermal energy is transferred during the charging phase (Balaras, C. A. 1996; Singh et al, 2010; Beasley and Clack, 1984).
Aly and ElSharkawy (1990) studied the effects of storage material properties on the thermal behaviour of packed beds during charging. A transient onedimensional twophase mathematical model was used to describe the temperature fields in the air and solid media constituting the bed. They found out that one of the main sets of parameters affecting the design of solid packed beds is the physical properties of the solid phase used as storage material thus choice of material is particularly important.
The aim of this study is to investigate the effect of key parameters of the rock bed system on the amount of heat stored after a fixed time of operation. The rock bed acts as regenerative heat exchanger. Cool air is used to lower the temperature of the rocks at night so that rocks absorb the heat from the higher temperature outside air during the day.
2. The Mathematical Model
The problem considered is a mathematical representation of heat transfer in the packed bed unit. It predicts the response of the unit to the time dependent fluid inlet temperature. The rock bed thermal storage unit under investigation is shown diagrammatically in Figure 1. The rocks are contained in a rectangular shaped wire cage set against the sides of a brick underground. The bed is used as a thermal storage to regulate temperature in the room space environment for a typical building. The energy transporting fluid is air, and energy transporting materials are rocks.
In the model, the air flow is assumed to be evenly spread in crosssection and rocks at a crosssection are assumed to have the same temperature distribution. For the sake of ease of analysing the temperature distribution within the rock, the shape is assumed to be spherical. Thus, for the mathematical model the bed is assumed to be evenly packed with spherical rocks of equal radius (that is, we assume the bed contains the same mass and size of rock with uniform crosssection area), with a fluid (air) spread uniformly in crosssection, flowing through the bed starting from the initial time.
The temperature of the rocks satisfies the heat diffusion equation. The energy equation for the fluid contains diffusion and heat transport terms. Coupling of the solid (rock) and fluid phases is through heat convection between the rock surface and the fluid, which is represented by a term in the equation of each phase. The boundary conditions for the fluid are that the initial temperature is known, and at the outlet end the fluid is assumed not to lose heat significantly to the duct. An additional term in the fluid equation allows for heat loss to the environment. This takes the form of convection to the surrounding earth or bricks which are at a known constant temperature. In the development of the model, the following assumptions are taken into consideration. The bed material temperature, satisfies the onedimensional parabolic equation (1) for.
(1)
with boundary conditions:
(2)
(3)
and initial condition:
(4)
also the fluid temperature, satisfies,
(5)
subject to the boundary conditions:
(6)
(7)
and initial condition:
(8)
where, is the area of convective heat transfer, is the specific heat capacity, is the bed void fraction (which is defined as the volume of gas per unit volume of rockbed), is the convective heat transfer coefficient, is the thermal conductivity, is the bed length, is the bed perimeter, is the radial distance, is the density, is the radius of the average solid particles (rocks), is the surface frontal area, is the temperature, is time, is the heat loss coefficient from bed to surrounding, is the distance in the direction of fluid flow, and is the mass flow rate of heat transfer fluid. For subscripts: stands for fluid, is fluid inlet, is bed material and stands for surface of a rock.
The problem variable is defined for any in the range, although the spheres occur in random discrete positions in the bed. Thus represents the temperature in an average sphere at position as it varies with time. Similarly is the fluid temperature at position varying with time. The term on the left hand side of equation (5) represents an energy rate associated with the motion of the fluid through the bed. On the right hand side, the first term represents fluid heat dispersion and, the second and third terms represent convective heat transfer on the fluidsolid interface and, thermal interaction with the environment, respectively.
3. Numerical Procedures
The mathematical model in equations (1)(4) is solved numerically using semidiscretization finite difference approximation. For equations (5)(8) finite element approximation is used. In order to develop a numerical method which uses either finite differences or finite elements, nodes are introduced for the space dimension along the bed and along the radius of the spheres. We define the following discretization in space for fixed time,, where and are node spacing and are the number of node points used in discretization of the sphere and bed respectively. Consider approximating the space derivatives with finite differences. Equation (1) is equivalent to
(9)
where and satisfies the boundary condition . A consideration of the symmetry of the solution shows that at the origin. Keeping constant, we treat as a function of only and expand the equation on the right hand side of equation (1) using a Taylor series around to obtain
(10)
which transforms equation (1) to
(11)
Using central difference quotient, at the first node, equation (11) becomes
(12)
since which follows from boundary condition (6). For the middle nodes the discrete form of equation (1) is
(13)
and at node is
(14)
Suppose an artificial node is introduced, from the surface boundary condition equation (2) we use central difference quotient for the first derivative for the surface boundary condition equation (5) to obtain,
(15)
By substituting this result into equation (14) we get
(16)
where is the interface rock temperature.
Equations (12) through (16) form a system of equations for the thermal energy material (solid) and may be written in matrix form as
(17)
where,
,
,
(The constant vector from the convection boundary condition) and is an matrix shown below:
(18)
where, , .
The system of equations (17) is consistent and consists of equations in unknowns for the solid.
4. Finite Element Solution of the Fluid Equation
Finite element method is used with the fluid temperature approximated by linear variation on each element. This has been chosen as the temperature along the bed does not vary sharply. Using isoparametric formulation of the element stiffness matrix and Garlekin's approach (Logan, 2007) we obtain the following matrix form.
(19)
where , , while and are the nodal temperatures to be determined. Applying the element equation to equation (19), we see that the element stiffness matrix is now composed of three parts:
(20)
where,
(21)
and the element force matrix and unknown temperature matrix are andrespectively. From equation (21) we observe that the mass transport stiffness matrix is asymmetric and, hence is asymmetric. Thus, the equations for fluid without boundary conditions becomes,
(22)
The stiffness matrix comes from the terms involving the fluid temperatures, and the element nodal force matrix from the convection and environmental terms. The condition that has to be imposed, that is, we set and zero for the rest of the rows and. If equation (22) is multiplied by and the derivative boundary condition is incorporated the resulting equation has the form:
(23)
where.
5. Semianalytic Method for Solution in Time
The rock and fluid temperatures are assumed to be known at a particular time, and we needed to find the rock temperatures at. The fluid temperatures are known at from the given inlet temperature. We assume that the fluid temperaturevaries linearly between times and, that is, the fluid temperature responds very quickly to a change in input temperature because of the flow speed. This contrasts with the change of temperature in the rocks. That is,
(24)
Then equation (17) becomes,
(25)
for with known and . Setting and, equation (25) becomes
(26)
Let, where the subscript stands for "additional" and. Then, equation (26) changes to
(27)
Suppose is a matrix that diagonalizes, then (the diagonal matrix) and since then equation (27) becomes
(28)
Letting and equation (28) can be written as
(29)
This is a set of uncoupled equations for the temperature at each along the rock cage. The general equation of uncoupled system is
(30)
Equation (30) is a first order differential equation and can be solved using the usual methods. For each row of the system (30) the solution is, at, we have
In general, we have
(31)
and,
6. Results and Discussions
The results of mathematical simulation of the system have been discussed in the system. The effect of the system parameters (convective heat transfer coefficient, mass flow rate, the length of the bed, and rock sizes (stones dimension)) on the thermal performance of the system has been investigated. It is thought that the temperature distribution in the bed is the most aspect, it is important to discuss the effect of system and operating variable on temperature distribution. Mathematical simulation has been used to yield the temperature of element of storage system during charging at a given instant of time.
Rock  Air  Other parameters 










 

Parameter  Range 
Coefficient of heat transfer coefficient, ()  6 to 12 
Length of rock bed ()  2 to 3 
Air flow rate ()  0.4 to 0.8 
Radius ()  0.05 to 0.2 
6.1. Effects of Convective Heat Transfer Coefficients
Figure 2 depicts the transient response of the packed unit for different peak values of the convective heat transfer coefficient. Calculation of the fluid temperature leaving the rock bed for various values of convective heat transfer coefficient () shows that this parameter has significant influence on the thermal performance of the system. A period of 18 h is chosen to allow closer look at the behaviour rather than show the whole fortnight. This covers two periods of fan operation from 12.00 noon to 15: 30 hrs, and 22.00 hrs to 6.00 hrs.
From the figure it can be seen that, the effect of increasing the convective heat transfer coefficient is to dampen the response of the fluid outlet temperature. This is because the interaction between the rocks and the air is enhanced. The bigger the convective heat transfer coefficient the greater the gain in the fluid temperature and this brings the predicted temperature closer to the measured value. In practice, mechanical ventilation is required to achieve reasonable values of convective heat transfer coefficient. However the actual values of convective heat transfer coefficient depend not only on the air flow rate but also on the level of turbulence created in the air flow in the channels, which is a function of surface roughness and obstructions in the air flow path (IsanskaCwiek, 2005). For each of the three convective heat transfer coefficients, the flow rate of is used for forced convection and the other parameters are as given in Table 1.
6.2. Effects of the Air Flow Rates
This section investigates the effect of the air flow rateby simulating low rates of, and to observe the agreement between predicted and measured fluid outlet temperature. Figure 3 shows the transient response of the morning unit for different peak values of denoted by in the figure. For each of the three air flow rates, the convective heat transfer coefficient is assumed to be constant at and the other parameters are as given in Table 1. Since air flow rate has significant impact on the value of the heat transfer coefficient value, the results are comparable to results given in section 6.1. This shows that the air flow rate parameter has a strong dependence on the thermal performance of the system.
From Figure 3 it can be seen that using the lower flow rate of is to dampen the fluid outlet temperature. The greater temperature gain achieved with a lower flow rate may be due to the air being given more time to be in contact with the rocks. The situation is complicated if it is assumed that convection is linked to air flow rate; then lowering the air flow rate in itself increases the peak damping, but the related result of decreasing the convective heat transfer coefficient reverses the effect. Also, the measured average value of recommended for use in the modelling is not a correct value. The fluid must be well mixed before taking a measurement otherwise the reading is not a true representative of the fluid flow rate.
6.3. Effects of Length of the Packed Bed
Figure 4 compares the outlet air temperature from rock beds of the different lengths. The figure indicates that the effect of increasing the bed length is to increase the peak damping of the fluid outlet temperature. That is, the temperature difference is greater as the bed length increases.
6.4. Effects of Rock Size
In Figure 5 the effects of radius of rocks (rock sizes) on the outlet air temperature from packed bed is displayed. It is found from the figure that, lowering the radius of the rocks increases the peak damping of the fluid outlet temperature. Hence, for efficient operation, smaller rocks allow for more efficient storage of energy than large rocks.
7. Conclusions
The combined effect of convective heat transfer coefficient, air flow rates, length of the packed bed and the rock size on the transient responses of the packed bed unit to the time dependent fluid inlet temperature was investigated. The partial differential equations are solved numerically using a semidiscretization finite difference method together with finite element method. The results are summarized as follows:
• The greater convective heat transfer coefficients tend to produces the greater gain in the outlet fluid temperature leaving the rock bed.
• The greater temperature gain can be achieved with lower flow rate and this is due to the air being given more time to be in contact with the rocks.
• For the efficient operation it is found also that, long bed length and smaller rocks allow more efficient storage of energy than small bed length and large rocks respectively.
References