Issue 
Emergent Scientist
Volume 1, 2017
International Physicists' Tournament 2016



Article Number  5  
Number of page(s)  7  
DOI  https://doi.org/10.1051/emsci/2017005  
Published online  22 September 2017 
Research Article
Cooling jug physics
^{1}
School of Physics and Technology, V. N. Karazin Kharkiv National University,
Svobody Sq. 4,
61022
Kharkiv, Ukraine
^{2}
Department of Aerospace Engineering Sciences, University of Colorado Boulder,
429 UCB,
Boulder,
CO
80309, USA
^{*} email: lunya4ek@gmail.com
Received:
5
August
2017
Accepted:
7
August
2017
We discuss the physics of the potinpot cooler. By balancing temperature decrease due to evaporation and temperature increase due to heat exchange, we find the equilibrium temperature of the pot. In this simplified model, the cooling jug acts as a psychrometer, and the theoretical prediction of our model is in a good agreement with psychrometric tables. Next, we study dynamics of the jug cooling. The cooling rate is limited by water vapour diffusion through air, heat conduction through air, and heat conduction through the body of the jug. The derived rate of temperature decrease is in general agreement with the result of our experiment. In the end, we discuss some additional factors, such as capillary effects in the raw clay, water viscosity in the capillaries, and impact of complex shape of the jug.
Key words: phase transitions / physical kinetics
© O. Luniachek et al., Published by EDP Sciences, 2017
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1 Introduction
Keeping food cool has always been an effective way of extending its storage life, and different devices have been used for this sake throughout history. Most of them consume energy to compensate for the energy exchange of the food with the environment owing to the zeroth law of thermodynamics, so that currently the energy spent for food preservation is about 0.5 kWh per person per day [1].
Evaporative cooling constitutes one of the most effective and widely spread ways of lowering the temperature, with its manifestations ranging from cooling hot food by blowing on it to lowering our own body’s temperatures by sweating, and from evaporating fluorocarbon refrigerants for retaining low temperature in common fridges to attaining the coldest temperature ever reached in ultra cold atom experiments [2].
In this article, we focus on the potinpot cooler, which harnesses evaporative cooling by water to keep foods fresh with zero energy consumption. It is a system made of two clay pots and a wet sand as a separator between them. The jug cools down due to the evaporative cooling of water from pores on the surface of the outer jug, to which the water gets from the wet sand. Evaporation is caused by heating of the outer surface by heat fluxes from both the atmosphere and the inner pot. At first we explain working principle of such a system in the static case, i.e. when the temperature of the jug surface and the external air are in equilibrium. In this model, we consider only the contribution of heat flow from the atmosphere to water evaporation. Then we make our model dynamic by taking into account heat fluxes from the inner pot.
Our consideration of this cooling system involves many simplifications. We explore only one particular case of an evaporative cooling system. More general approaches are presented in [3,4].
2 Methods
2.1 Static equilibrium temperature of the jug
We will find equilibrium temperature T_{1} of a wet jug when the air is of temperature T_{air} and relative humidity ϕ_{air}. Our model is shown in Figure 1. When T_{air} > T_{1}, the surface of the jug gets heat via heat conduction from the exterior air, and uses the heat to evaporate water from the surface. The water vapour is then transported away from the surface by diffusion. In equilibrium, the mass of water evaporated due to the heat coming in by heat conduction, equals mass of water vapour transported away from the surface by diffusion, and it is the equilibrium between heat conduction and diffusion that determines the equilibrium surface temperature T_{1}. We assume that both heat and diffusive fluxes propagate across a flat boundary layer of thickness l, in which gradients of both temperature and absolute humidity are uniform. Beyond this layer (at the distance l from the surface) air is assumed to be well mixed, and temperature and humidity are assumed to be the same as in the bulk of air. We also assume that the jug itself is in thermal equilibrium with its surface, so that the surface exchanges heat only with air, and that water covers all the surface, so that the pressure of water vapour at the surface is the same as the pressure of saturated vapour at temperature T_{1}. We compare air temperature at the distances equal or bigger than boundary layer to the temperature of the dry thermometer (air is well mixed there), while temperature of saturated water vapour at the surface of the jug is compared to the wet thermometer. Later we are going to analyze these assumptions more scrupulously, and to go beyond them.
Following Newton’s law [5], the heat flux to the surface is (1) where κ is the heat conductivity of air. This heat evaporates the following mass of water from unit surface per unit time, (2)where L is the latent heat of vapourization of water.
The rate at which the evaporated water diffuses away from the surface, is governed by Fick’s law, (3) where D is the diffusion coefficient of water vapour in air, and ρ(T_{1}) is density of the saturated water vapour as function of temperature.
In equilibrium, equation (2) equals to equation (3), resulting into the following equation for T_{1}, (4)
Importantly, l cancels from this equation, and we do not have to make any assumptions about the boundary layer to determine the equilibrium temperature.
To solve this equation, we need to substitute the functional form of ρ(T_{1}). Although such approximations for the pressure dependence p(T_{1}) are well developed [6,7], approximations of ρ(T_{1}) are scarce. For this reason, we had to construct our own approximation. We chose a functional form similar to the August equation or, more generally, to Clausius–Clapeyron equation, (5) where the temperature T_{1} is given in degrees Celsius, while A and B are free constants, chosen to best fit the experimental data in the temperature range 10–50 °C. In Figure 3, 50 °C is the theoretical extrapolation, our approach can not define humidity precisely at such temperatures. The best fit is provided by A = 369 000 kg m^{−3}, B = 4943 K. It is shown in Figure 3, overplotted with the experimental data [8].
When we substitute equation (5) into equation (4), we get a transcendental equation, that cannot be solved using elementary functions. Still, we notice that ΔT = T_{air} − T_{1} is usually only a few degrees Celsius, so that the difference between ρ(T_{air}) and ρ(T_{1}) is small. That allows us to decompose ρ(T_{1}) = ρ(T_{air} + ΔT) into a Taylor series in terms of ΔT, truncating the series after the linear term. The general form of the solution to the resulting linearized equation is (6)
To verify our method, we reverse equation (6) using equation (5) to obtain humidity, (7)
Then in Figure 4 we overplot this theoretical expression for humidity with data from a psychrometric table^{1} [9].
We see a decent agreement between the data and the theoretical prediction, which justifies the methods we are using. Still, for big temperature differences ΔT we see an increasing discrepancy between the theory and the data. This should be expected, as the accuracy of the Taylor approximation decreases with ΔT.
Next we use equation (6) to predict the ΔT for the cooling jug, and plot the results in Figure 5. As expected, the cooling ΔT is larger when the ambient temperature is high, and the humidity low. According to our theory, at T_{air} = 50 °C the temperature difference can reach up to 20 °C, although as we see from Figure 4 the theory tends to give bigger errors at big ΔT.
Although our approach is less precise than the available phenomenological formulae [10], it gives a clearer understanding of the underlying physics, and allows for simple generalizations to nonstationary cases, as we demonstrate in the following subsection.
Fig. 1 Geometry of the wall of potinpot with boundary layer and the distribution of the heat fluxes at static equilibrium. 
Fig. 2 Geometry of potinpot cooler and the distribution of the heat fluxes at dynamic equilibrium. 
Fig. 3 Experimental relation between the temperature and the density of the saturated water vapour [8], approximated by equation (5). 
Fig. 4 Data from the psychrometric table (circles) overplotted with the theoretical formula – equation (7) (solid lines). The table shows humidity ϕ_{air} versus temperature T_{air} of the exterior air (dry thermometer) . Different temperature differences ΔT between the dry and the wet thermometer are plotted in different colours. 
Fig. 5 The temperature difference ΔT between the cooling jug and the exterior air as given by equation (6), plotted against the temperature of the exterior air T_{air}. Situations corresponding to different humidities ϕ_{air} are shown in different colours. 
2.2 Dynamics of jug cooling
We want to consider how the finite heat conductivity of the jug walls hinder its cooling speed. We assume the jug to be filled with a well mixed liquid, whose temperature T_{in} is uniform throughout its volume, and whose specific heat capacity is c_{liq}. The wall consists of three layers, as shown in Figure 2. Two layers are composed of clay and one is composed of sand.
To focus on the evaporation itself, we consider a simplified geometry where the wall thickness of both jugs, the distance between the jugs and the thickness of the air boundary layer are all much smaller than the inner radius of the inner jug. This allows us to treat the heat conduction and diffusion problems as onedimensional. The theory required to go beyond this assumption will be briefly reviewed in Section 4.
So far we have only considered the influence of the outside heat flux on water evaporation from the surface. The inner heat flux should be also taken into account. The heat fluxes in the different layers q_{1}, q_{2}, q_{3} and q_{air} (see Fig. 2) are determined by the heat conductivity equation: (8) (9) (10) (11)
The heat flow from inside the jug q_{1} contributes to water evaporation from the surface. The sum of the outer and the inner heat flows to the surface is the reason for water evaporation, as it is represented by the heat balance equation, (12)
The same heat flow q_{1} causes the jug to cool down in accordance with the energy balance equation, (13) Here S is again the total surface area of the jug, and C is the total heat capacity of the jug, (14)with m and c with different subscripts denoting masses and specific heat capacities of liquid, sand and clay.
These equations allow us to calculate cooling rate of the jug. We assume that the thermal inertia of the liquid inside the jug is much bigger than the thermal inertia of the jug itself. Then approximately q_{1} = q_{2} = q_{3}, which allows us to eliminate T_{2} and T_{3} from equations (9) to (11), and to express q_{1} in terms of T_{1} and T_{in}. Then we substitute this equation together with equation (8) into equation (12), express dm/(dtdS) from equation (3), and thus get an equation for T_{1}. Again applying the Taylor decomposition to ρ(T), we linearize this equation, and easily find its solution T_{1}. From the resulting T_{1} we calculate q_{1}, and substitute it into equation (13). Thus we arrive at a linear differential equation for T_{in}, whose solution is (15) with the relaxation time τ determined as (16)
The boundary layer for a laminar flow can be expressed in the following way [11]: (17) Here η is the dynamic viscosity, R is the characteristic length of the pot, ρ is the density of the oncoming flow, and v is the velocity of this flow, and α is a constant of the order of unity, which is commonly set at α = 5 [11]. Then the thickness l corresponds to the distance from the surface, where the air speed reaches 99% of the speed of the wind [11]. Assumption that the air is well mixed just outside this area and poorly mixed inside is somewhat arbitrary, and setting different speed limits instead of 99% would correspond to a different α. Moreover, the boundary layer has different thickness at different parts of the jug, thus making equation (17) a mere estimate. To account for these uncertainties, we take α = 5 ± 2.
Other parameters relevant for our model are listed in Table 1 together with their uncertainties. L, D, c_{liq} , κ_{air} and η are given at room temperature, and the uncertainty is computed by assuming temperature variations of 10 °C. Uncertainties of c_{clay}, c_{sand} , κ_{clay} strongly depend on the kind of sand and clay used, and are estimated from the range of values cited for these substances (e.g. [12,13]). Uncertainty in κ_{sand} is predominantly dictated by the uncertain wetness of sand [14]. For the speed of wind we took typical speeds of air in a room.
Other important limitations of the proposed theory include different surface areas of different boundaries, as well as inequality of heat fluxes at these boundaries, because some heat is stored in heat capacity of the material. Heat exchange through the top and the bottom of the jug was also neglected.
It is interesting though, that despite these shortcomings, our model can give reasonable estimates for the cooling rate, as we will see in the following section.
Values of the parameters used in our models.
3 Results
3.1 Experiment
For the experiment we use the two clay pots of different size. The space between the pots is filled with sand. Although we give theory for the case when the inner pot is filled with water, our experiment is carried out with wet sand. Such replacement is appropriate as we only need to change specific heat capacity and heat conductivity in our calculations. In order for refrigeration to occur, we need moisture, so we pour water into the sand between the pots. The sand becomes wet. Water enters the micropores of unbaked clay, due to the capillary effect. Shortly after the start of the experiment, we observe the formation of droplets on the outside of the pot. Drops evaporate from the entire surface of the pot and the system is cooled.
If the outer pot is made of baked clay or covered with paint, the micropores on the surface are closed, and water is not able to appear on the outer surface for evaporation.
In order to neglect the evaporation of water directly from the space between pots and to reduce heat loss to the environment the system is sealed with foil on top.
The humidity of the ambient air is ϕ_{air} = 55%. The thermometer is put inside of the cup. The room temperature is 18.0 °C. After 5 h the system temperature is equal to 13.5 °C. Further, the temperature slowly decreases, approximately following an exponential law. In order for the system to continue cooling, we need to add water, because otherwise the sand becomes drier. In the beginning of the experiment we poured about 0.5 L of water into the sand. Then we add water every half an hour, resulting in the total of 300 mL of water added. As the heat capacity of the added water is much smaller than the heat capacity of sand, adding the water could not much influence the cooling rate of the system, although certainly should have been accounted for in a more precise theory.
Our setup is shown in Figure 6. Noteworthy, in our experiment the wall thickness of the jug is not small as compared to the radius, in direct contradiction with assumptions of the theory. Therefore, our theory can describe the experiment only approximately.
Temperature versus time is shown in Figure 7. Substitution of the relevant physical parameters and their uncertainties into equation (15) with ΔT = 5.6 °C and h. This theoretical uncertainty is shown in the figure as lightgreen shaded area.
Note, that one experiment is not sufficient to prove a theory with so many parameters, but merely to illustrate it. Moreover, giving substantial uncertainties in the parameters involved in the theory (especially the coefficient in Eq. (17)), such a good agreement of the theory and the experiment should be coincidental to some extent. Still, it allows us to argue, that our theory well reproduces the general temperature trend and gives at least a reasonable orderofmagnitude estimate for the cooling rate.
Fig. 6 View of the experimental system. Its properties are as follows. Inner pot: height – 15 cm, diameter – 9 cm, thickness – 6 mm. Outer pot: height – 40 cm, diameter – 18 cm, thickness – 6 mm. We used 10 kg of sand and 0.8 L of water. 
Fig. 7 Results of the experiment (red dots) overplotted with the leastsquare exponential fit (blue line) and the theoretical predictions by equation (15) (green line). The area shadowed in light green shows the theoretical uncertainty. 
4 Discussion
Some physical factors omitted by our theory could in principle be incorporated into the model, although their inclusion would make the model overly complicated analytically, while their practical significance is hard to estimate. Here we only briefly discuss the underlying physics, and its effect on our model.
4.1 Geometric shape
Our previous calculations were done under the assumption of a jug with a thin shell and an interior of a constant temperature. To understand the cooling rate of the shell or the temperature inhomogeneity in the interior, we must do a rigorous 3D simulation of heat diffusion in the jug. This means going from an ordinary differential equation like equation (13) to a partial differential equation.
Matter and energy balance on the surface results into the boundary condition for this PDE. To derive it, we equate equations (2) and (3). By transforming and simplifying this expression, we get an inhomogeneous boundary condition of the third kind: (18) with the coefficients (19) (20) (21)
To get the temperature distribution, the heat conductivity equation has to be solved: (22) Here Δ denotes the Laplace operator.
For a cylindrical pot, the Laplace operator can be rewritten in cylindrical coordinates, and the solution will be obtained as a series in terms of the Bessel functions, although the coefficients of this series can be found only from a transcendental algebraic equation. For a spherical pot, the problem can be simplified in spherical coordinates, and the solution can be obtained in a similar manner as a series in terms of spherical functions. For a pot of an arbitrary shape, there are no general methods to solve the equation, and it should be simulated numerically.
4.2 Viscosity and surface tension
One more limitation to water evaporation from the surface is set by the viscosity of water. On its way through the clay, water has to overcome viscous friction in capillaries, which in the case of a cylindric capillary is governed by Hagen–Poiseuille equation, (23) Here μ is the dynamic viscosity of water, ρ is its density, L is the length of the capillary, R is its radius, Δp is the pressure difference between the two ends of the capillary, and dm/dt is the mass flow rate of water.
If the sand is not soaking wet and can hold water inside, water pressure inside it is not greater than the atmospheric pressure. Then to move water through the capillary, the water pressure on the outer side of the jug must be less than the atmospheric pressure by some amount Δp, which could only be attributed to the Laplace pressure on the outer side of the capillary, (24) Here σ is water surface tension, and the radius of curvature R is of the same order of magnitude as the radius of the capillary.
Now we can see the nontrivial interplay between these phenomena. If water experiences high friction in capillaries, the water pressure on the outer surface of the jug falls. Now water is sucked inside the thinnest capillaries, which have the biggest Laplace pressure. But the same Laplace pressure makes it harder for water to evaporate, and diminishes the pressure of saturated vapour on the surface of the jug. By this mechanism viscosity decreases the cooling rate of the jug.
4.3 Deadends
The experiments, with which we started modelling this problem and which we presented at the tournament, were much cruder than the experiments described above. The major problem was that we did not have pots made of raw clay. So we used materials at hand to construct a setup shown in Figure 8.
A paper cup was covered with a cloth rag from inside and outside. Another cup was placed inside it. Water was put in the first cup, as well as to the space between the cups. Water rose by the cloth due to the capillary effect, and evaporated, cooling the setup. The thermometer was put inside the cup. The room temperature was 16.2 °C, and in less than an hour the water temperature inside the cup fell down to 11.5 °C. Figure 9 shows the temperature as a function of time. The general trend is similar to Figure 7, yet the exponential law now provides a worse fit to the data.
Although representing the same physics, this model significantly deviated from the problem statement. Moreover, water oozing through a folded cloth of complicated geometry presents a problem, whose precise physical description is not easier than the original problem of potinpot refrigerator.
Fig. 8 The experimental setup used for our fist experiments. 
Fig. 9 Results of our first experiment (red dots) ant the leastsquare exponential fit (blue line). 
5 Conclusions
The cooling of the pot is caused by water evaporation from its surface. The minimal temperature attainable by the pot is determined from the mass and energy conservation laws at the surface. The rate at which this temperature will be achieved is limited by the heat conductivity through the pot, the heat conductivity through the boundary layer, and the diffusion of water vapour through the boundary layer.
In the simplest model we get a linearized algebraic equation for the minimal temperature, and a linear differential equation for the process of cooling to this temperature. The solution of the former is in reasonably good agreement with the psychrometric table. The solution of the latter as a function of time is an exponential function, which qualitatively agrees with our experiments.
To make the model more precise, one can include additional factors, such as the complex geometric shape of the jug, the influence of resistance to water flow in the capillaries, the curvature of the water surface in the capillaries, and air convection as a mixing agent that supplements the exterior wind. But this added precision will come on an intolerable cost – tremendously increased complexity.
References
 D.J. MacKay, Sustainable energy – without the hot air, 1st edn. (UIT Cambridge, Cambridge, 2009), pp. 50–53 (In the text)
 M. Weidemüller, C. Zimmermann, Interactions in ultracold gases: from atoms to molecules (WileyVCH Verlag GmbH & Co. KGaA, Weinheim, 2003), pp. 46–54 (In the text)
 B. Halasz, A general mathematical model of evaporative cooling devices, Rev. Gén. Therm. 37, 245–255 (1998) [CrossRef] (In the text)
 J.M. Wu, X. Huang, H. Zhang, Theoretical analysis on heat and mass transfer in a direct evaporative cooler, Appl. Therm. Eng. 29, 980–984 (2009) [CrossRef] (In the text)
 D.V. Sivukhin, General physics course: thermodynamics and molecular physics (Nauka, Moscow, 1990), Vol. II, pp. 347–349 (In the text)
 G.W. Thomson, The Antoine equation for vapour pressure data, Chem. Rev. 38, 1–39 (1946) [CrossRef] (In the text)
 O.A. Alduchov, R.E. Eskridge, Improved magnus form approximation of saturation vapor pressure, J. Appl. Meteorol. 35, 601–609 (1996) [CrossRef] (In the text)
 L. Haar, J.S. Gallagher, G.S. Kell, NBS/NRC steam tables: thermodynamic and transport properties and computer programs for vapor and liquid states of water in SI units, in Hemisphere (1984), p. 8 (In the text)
 N.I. Koshkin, M.I. Schirkevitch, Handbook of elementary physics (Nauka, Moscow, 1975), p. 74 (In the text)
 R. Stull, Wetbulb temperature from relative humidity and air temperature, J. Appl. Meteorol. Climatol. 50(11), 2267–2269 (2011) [CrossRef] (In the text)
 H. Schlichting, et al., Boundarylayer theory (McGrawHill, New York, 1960), Vol. 7, pp. 31–140 (In the text)
 Engineering ToolBox, http://www.engineeringtoolbox.com (In the text)
 H.M. AbuelNaga, D.T. Bergado, A. Bouazza, Thermal conductivity evolution of saturated clay under consolidation process, Int. J. Geomech. 8(2), 114–122 (2008) [CrossRef] (In the text)
 K.M. Smits, T. Sakaki, A. Limsuwat, T.H. Illangasekare, Determination of the thermal conductivity of sands under varying moisture, drainage/wetting, and porosity conditionsapplications in nearsurface soil moisture distribution analysis, in Hydrology Days (2009), pp. 57–65 (In the text)
Psychrometric table is used to determine relative humidity based on measurement results of a psychrometer, i.e. a device consisting of two thermometers, a dry one and a wet one. The dry thermometer shows the temperature of the ambient air, while the thermometer wrapped in wet fabric shows the temperature of evaporating water. The readings of the thermometers can then be translated into the relative humidity of the air.
Cite this article as: Oleh Luniachek, Ruslan Timchenko, Oleksiy Golubov, Cooling jug physics, Emergent Scientist 1, 5 (2017)
All Tables
All Figures
Fig. 1 Geometry of the wall of potinpot with boundary layer and the distribution of the heat fluxes at static equilibrium. 

In the text 
Fig. 2 Geometry of potinpot cooler and the distribution of the heat fluxes at dynamic equilibrium. 

In the text 
Fig. 3 Experimental relation between the temperature and the density of the saturated water vapour [8], approximated by equation (5). 

In the text 
Fig. 4 Data from the psychrometric table (circles) overplotted with the theoretical formula – equation (7) (solid lines). The table shows humidity ϕ_{air} versus temperature T_{air} of the exterior air (dry thermometer) . Different temperature differences ΔT between the dry and the wet thermometer are plotted in different colours. 

In the text 
Fig. 5 The temperature difference ΔT between the cooling jug and the exterior air as given by equation (6), plotted against the temperature of the exterior air T_{air}. Situations corresponding to different humidities ϕ_{air} are shown in different colours. 

In the text 
Fig. 6 View of the experimental system. Its properties are as follows. Inner pot: height – 15 cm, diameter – 9 cm, thickness – 6 mm. Outer pot: height – 40 cm, diameter – 18 cm, thickness – 6 mm. We used 10 kg of sand and 0.8 L of water. 

In the text 
Fig. 7 Results of the experiment (red dots) overplotted with the leastsquare exponential fit (blue line) and the theoretical predictions by equation (15) (green line). The area shadowed in light green shows the theoretical uncertainty. 

In the text 
Fig. 8 The experimental setup used for our fist experiments. 

In the text 
Fig. 9 Results of our first experiment (red dots) ant the leastsquare exponential fit (blue line). 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.