Modeling and Numerical Simulation of Heat Transfers in a Metallic Pressure Cooker Isolated with Kapok Wool

In this work, a numerical study of heat transfers in a metallic pressure cooker isolated with kapok wool was carried out. This equipment works like a thermos, allowing finishing cooking meals only thanks to the heat stored at the beginning of cooking, which generates energy savings. Cooked meals are also kept hot for long hours. In our previous work, we have highlighted the performances of the pressure cooker when making common dishes in Burkina Faso. Also, the parameters (thickness and density) of the insulating matrix allowing having such performances as well as the influence of the climatic conditions on the pressure cooker operation were analyzed in detail in this present work. The numerical methodology is based on the nodal method and the transfer equations obtained by making an energy balance on each node have been discretized using an implicit scheme with finite differences and resolved by the Gauss algorithm. Numerical results validated experimentally show that the thickness of the kapok wool as well as its density play an important role in the pressure cooker operation. In addition, equipment performances are very little influenced by the weather conditions of the city of Ouagadougou (Burkina Faso).


Introduction
Thermal insulation is one of the most used solutions to save energy in different How to cite this paper: Ouedraogo, D., Igo, S.W., Sawadogo, G.L., Compaore, A., Zeghmati, B. and Chesneau, X. (2020) Modeling and Numerical Simulation of Heat Transfers in a Metallic Pressure Cooker Isolated with Kapok Wool. Modeling and Numerical Simulation of Material Science, 10, sectors (building, catering, industry, etc.). In many sectors, natural insulation is more and more preferred over synthetic insulation because of the requirements related to sustainable development. Natural insulators can be vegetable, animal or mineral origin. Among these natural insulators that are found in abundance in tropical countries like Burkina Faso, is kapok wool. Numerous studies have demonstrated the ability of this fiber to compete with synthetic insulation.
Experimental studies conducted by Voumbo et al. [1] on kapok wool, using the plate method, show very interesting thermophysical properties. Indeed, they have shown that the thermal conductivity of kapok wool varies between 0.03 and 0.04 W/m•K for a density between 5 and 40 kg/m 3 and an average thermal diffusivity of 17.1 × 10 −7 m 2 /s. J. C. Damfeu et al. [2] determined by an experimental measurement using the radial flux method the thermo-physical properties at low density of the following natural fibers: fibers kapok; peanut shell fibers; rattan fiber and coconut fiber.
The results presented show that the thermal conductivity of kapok wool (λ = 0.045 W/m•K) is in agreement with that reported by Voumbo et al. [1].
A. Wereme et al. [3] compared the storage capacity of ice in a container insulated with kapok wool and in the shell of a refrigerator insulated with polyurethane. They found that the insulated container with kapok wool could keep the ice for 6 days while the shell of the refrigerator could do it in 7 days, which shows that kapok wool can compete with polyurethane.
In the same way aming Wang, [4] experimentally analyzed the heat transfer in the kapok using a sensor (HIH-3610). It appears that the use of kapok as a thermal insulating material in the fabric helps minimize heat loss from the human body. The influence of temperature and wind speed on heat transfer through the kapok was analyzed and compared to that of cotton. It has been shown that with kapok the heat losses of the human body are lower than those obtained by cotton.
In the building sector, a study focused on a comparative analysis of the thermal performance of a non-insulated house with that of a house whose roof is insulated with kapok fiber shows that kapok fiber is suitable for insulating the roof in residential buildings in a hot climate [5].
Other work has shown that the combination of kapok wool and other materials can lead to very interesting materials in terms of thermal insulation. Adulkareem S et al. [6] show through an experimental study that composite materials resulting from the combination of kapok wool and sugarcane bagasse or coconut fiber at percentages of 50% lead to very good insulating materials, have shown that adding a small amount of kapok wool to the plaster leads to an improvement in the thermophysical properties of the composite material [1]. Because of its very interesting thermophysical properties, we have recently used kapok wool in the thermal insulation of a pressure cooker intended for cooking and meals conservation. The experimental study by Drissa Ouedraogo et al. [7] which followed showed that the device makes it possible to achieve enormous energy savings (between 30% and 75%) in the cooking of local dishes, and its constant time is about 60 hours. However, even if these results are mainly due to the thermophysical properties of kapok wool, the influence of other parameters such as the thickness and density of the insulating matrix on the performance of the device was not elucidated during this study. However, these parameters are essential when reproducing or resizing the device. In addition, the behavior of the pressure cooker in a typical climate of Burkina Faso, throughout the year is not known, which could limit its application in a given period. Consequently, the objective of this work is to carry out a numerical simulation in order on the one hand to understand the influence of the parameters (thickness, density) of the insulating matrix on the operation of the device and on the other hand to know the influence of the climate on his behavior during the year.

Description of the Pressure Cooker
The pressure cooker has a parallelepipedic shape (L × W × H = 65 cm × 60 cm × 50 cm) box made of 2 mm thick steel sheet which encloses the insulation system consisting of kapok wool. The thickness of kapok wool varies from 1 to 20 cm while its density varies from 10 to 50 kg/m 3 . The kapok wool is protected from the outer wall by a 15 mm thick wooden sheet. The pot containing hot water is placed in the middle of the pressure cooker and surrounded by kapok wool (Figure 1).

Simplifying Hypotheses
We adopt the following simplifying hypotheses: • The hot water temperature inside the pressure cooker is uniform; • The thermo-physical properties of air and materials are constant; • The external temperature is taken equal to the ambient temperature; • Heat exchanges through the cover wall are negligible; • Materials are assimilated to gray bodies; • The celestial vault is considered to be a black body; • The atmospheric diffuse radiation is isotropic.
According to the above hypotheses, the model of heat transfers in the pressure cooker is represented in the Figure 2.

Heat Transfers Equations in the Pressure Cooker
In general, modeling a thermal system using the nodal method amounts to setting up a network of thermal capacities ( i , heat sources and thermal coefficients. Heat transfers are represented by a thermal resistance network. By applying Kirchhoff's law [8], to the network of thermal resistances schematizing the different modes of heat transfer, we give the general heat balance equation by expression (1). This equation is based on nodal analysis by Boyer et al., [9] and Thomas Nganyaa et al. [10].
j index of the medium for which T j is a potential connected to the potential T i , Q i heat source (W), G ijx heat transfer coefficient between the media i and j (W/m 2 K) according to the heat transfer mode x (conduction, convection, radiation) and S, surface of the section considered (m 2 ). The thermal heat transfer coefficients G ijx are defined by the formal expression of the energy flow Ф ij , exchanged between two adjacent nodes i and j ( ) Figure 3 show the heat transfer model.

Expressions of Heat Transfer Coefficients
The heat exchanges in the pressure cooker take place according to the three modes (convective, radiative and by conduction).
The convective exchange coefficient between the outer wall of the pressure cooker and the ambient air depends on the action of the wind speed. This coefficient h cv is deduced from the relation of Mac Adams [11] [12].
• The heat transfer coefficient between the underside of the horizontal wall of the pot and the ambient air is deduced from the relation of Incroprera & De Witt [15].
• The coefficient of radiative transfer between the external wall of the pressure cooker and the celestial vault is determined by the Expression (11): In the case of a vertical wall, the radiative form factor F sky verifies the expression determineted by Ivanova et al. [16] and J. Ramirez et al. [17].
( ) b: represents a parameter which is a function of the anisotropy of the sky.
For an isotropic sky (b = 0) the radiative from factor is equal to 0.5. In the literature, there are several correlations for determining the temperature of the sky. We selected the one proposed by Swinbank [18].

Discretization of Equations
The Equations (3)-(7) are discretized by an implicit method with finite differences. This method is based on a Taylor series development which transforms partial differential equations into a system of algebraic equations. So, the discretized equations are given by the Expressions (13) • At the wood sheet • Within the pot The algebraic Equations (14)-(23) can be written in the form of a system of equations as follows

Solving the System of Equations
At t 0 + Δt, we assigned an arbitrary value to the temperatures on components of the pressure cooker. The heat transfer coefficients (convective, radiative) are thus calculated. The resolution of the system of algebraic Equations (3)-(7) leads to new values of the temperature of the components of the pres-sure cooker which are compared with the arbitrary value. If the difference between these two temperatures is greater than the desired precision, the values of the calculated temperatures replace the arbitrary value and the procedure described below is repeated until the convergence is obtained. The convergence was obtained when the following criterion was satisfied:

Model Validation
Using the methodology described in Section 2.4, we obtain the numerical curve of hot water cooling under the conditions: Hot water initial temperature: 90˚C, ambient temperature: 25˚C, kapok wool density: 50 kg/m 3 , kapok wool thickness: 15 cm. This numerical curve is then compared to the hot water cooling curve obtained experimentally under the same conditions ( Figure 4). We note a good qualitative agreement between these two results. The discrepancies do not exceed 5% and are due to the various empirical correlations used for the calculation of the heat transfer coefficients.

Influence of Kapok Wool Thickness and Density on the
Pressure Cooker Operation Figure 5 shows the influence of the thickness e k (for 1 to 20 cm) of kapok wool on the cooling curve of the hot water contained in the pressure cooker. The initial temperature inside the pot is taken at 90˚C and the outside temperature is 25˚C. We note that the temporal evolution of the water temperature decreases rapidly to tend asymptotically towards a constant value for small thicknesses of kapok wool (1, 2, 3, 4 and 5 cm). This leads to rapid cooling of hot water for low values of the thickness of the insulation. For large thicknesses of kapok wool, particularly from 15 cm, we note that the cooling curve flows to straight lines, which indicates a low cooling of the hot water over time. This result is explained by the fact that at low thickness, the thermal resistance of the insulating material is low, which increases transfers to the outside environment. Figure 6 shows the influence of the density (for 5 to 50 kg/m 3 ) of kapok wool on the cooling curve of the hot water contained in the pressure cooker. The initial temperature inside the pot is taken at 90˚C and the outside temperature is 25˚C. Figure 6. Influence of the density of the kapok wool on the cooling curve of the hot water. Modeling and Numerical Simulation of Material Science As in the previous case, it can be seen that the increase of the kapok wool density leads to poor cooling of the hot water over time. On the other hand, for low density value, a very rapid cooling of the hot water is observed. This result is due to the fact that at low density, the thermal diffusivity of the insulating material is high, which increases transfers to the outside environment.

Simulation of the Pressure Cooker Operation with
Meteorological Data from the City of Ouagadougou Figure 7 shows the monthly change in ambient air temperature following a seasonal cycle at Ouagadougou. We note that the highest temperatures are observed in March April and May. The maximum temperature was recorded during the month of April and reached the value of 40˚C. The lowest temperatures are observed during the months of December, January and August. January is the coldest month with a maximum and minimum temperature of 32˚C and 16˚C respectively. Figure 8 shows the influence of the hottest and coldest months on the operation of the pressure cooker. The initial temperature inside the pot is taken at 90˚C.

Influence of the Hottest and Coldest Months on the Operation of the Pressure Cooker
We notice that there is a similar change in the cooling curve for the two months. The difference between the two curves is more pronounced for the high times because of the temperature gradient between the two months. However, the transfers are more important for the month of January which present lower temperatures. The GAP in final temperatures between the two months is approximately 4˚C, which is relatively small compared to the initial temperature of the hot water.  Figure 9 shows the influence of the hottest month on the wall of the pressure cooker. The initial temperature inside the pot is taken at 90˚C. We note that the temperature of the wall follows that of the ambient air which reaches a maximum between 12 h and 15 h because of the sunshine. Despite the fluctuations in the wall temperature, it can be seen that the cooling curve of the hot water is not disturbed, which shows the effectiveness of the insulation. Figure 10 shows the influence of the 12 typical day of year on the operation of the pressure cooker. The initial temperature inside the pot is taken at 90˚C.

Influence of the 12 Typical Days of the Year on the Operation of the Pressure Cooker
The observation in Figure 10 shows that the evolution over time of the cooling curve for the typical 12 days is similar. We also note that the temperature gradient between the different cooling curves is very low from one month to the next. These results show that the use of the pressure cooker is not conditioned by climatic by climatic variations during the year.

Conclusions
A metallic pressure cooker insulated with kapok wool has been numerically investigated in the present work. The numerical methodology is based on the nodal method and the transfer equations obtained by making an energy balance on each node have been discretized using an implicit scheme with finite differences and resolved by the Gauss algorithm. The main results are summarized as follows: -The pressure cooker performances are linked to the good thermophysical properties of kapok wool, -The increase of the thickness and density of kapok wool leads to optimal operation of the pressure cooker, -The pressure cooker operation is very little influenced by the climatic conditions of the city of Ouagadougou. Consequently, it can be used at any time of the year. Its permanent use will considerably reduce energy consumption in the catering sector, as we have shown in our previous work.