3-D Numerical Modeling of Heat Transport Phenomena in Soil under Climatic Conditions of Southern Thailand



This paper presents a 3-D numerical modeling of heat transport phenomena in soil due to a change of sensible and latent heat, under the ambient conditions of southern Thailand. The vertical soil temperature profile within 3 m was predicted based on energy balance and 3 modes of heat transfer mechanisms, i.e., conduction, convection, and radiation. Mathematical models for estimation of solar radiation intensity, ambient and sky temperatures, relative humidity, and surface wind velocity were used as model inputs. 3-D numerical implicit finite difference schemes, i.e., forward time, and forward, center, and backward spaces were used for discretizing the set of governing, initial, and boundary condition equations. The set of pseudo-linear equations were then solved using the single step Gauss-Seidel iteration method. Computer code was developed by using MATLAB computer software. The soil physical effects; density, thermal conductivity, emissivity, absorptivity, and latent heat on amplitude of soil temperature variation were investigated. Numerical results were validated in comparison to the experimental results. It was found that 3-D numerical modeling could predict the soil temperature to almost the same degree as results that were obtained by experimentation, especially at a depth of 1 m. The root mean square error at ground surface and at depths of 0.5, 1, 1.5, 2, 2.5 and 3 m were 0.169, 0.153, 0.097, 0.116, 0.120, 0.115, and 0.098, respectively. Furthermore, it was found that variation of soil temperature occurred within 0.75 m only.



Finite difference, heat transfer, numerical modeling, soil

Full Text:



G Mihalakakou, M Santamouris and D Asimakoupoulos. Modeling the thermal performance of earth to air heat exchangers. Solar Energ. 1994; 53, 301-7.

M Santamouris, G Mihalakakou, C Balaras, A Argiriou and D Asimakoupoulos. Use of buried pipes for energy conservation in cooling of agricultural greenhouse. Solar Energ. 1995; 55, 111-24.

NH Abu-Hamdeh. Measurement of the thermal conductivity of Sandy Loam and Clay Loam soils using single and dual probes. J. Agri. Energ. Res. 2001; 80, 209-16.

OC Popiel, J Wojtkowiak and B Biernacka. Measurements of temperature distribution in ground. Exp. Therm. Fluid Sci. 2001; 25, 301-9.

E Penrod, MJ Elliot and WK Brown. Soil temperature variation at Lexington, Kentucky. J. Soil Sci. 1960; 90, 275-83.

JE Carson. Analysis of soil and air temperature by Fourier techniques. J. Geophys. Res. 1963; 68, 2217-32.

G Mihalakakou. On estimating soil surface temperature profiles. Energ. Build. 2002; 34, 251-9.

B Givoni and L Katz. Earth temperatures and underground buildings. Energ. Build. 1985; 8, 15-25.

S Krishnaiah and DN Singh. A methodology to determine soil moisture movement due to thermal gradients. Exp. Therm. Fluid Sci. 2003; 27, 715-21.

MMS El-Din. On the heat flow into the ground. Renew. Energ. 1999; 18, 473-90.

G Mihalakakou, M Santamouris and D Asimakoupoulos. Modeling the earth temperature using multiyear measurements. Energ. Build. 1992; 9, 1-9.

P Bloomfield. Fourier Analysis of Time Series: An Introduction. John Wiley & Sons, New York, USA, 1976.

CP Jacovides, G Mihalakakou, M Santamouris and JO Lewis. On the ground temperature profile for passive cooling applications in buildings. Solar Energ. 1996; 57, 167-75.

JP Holman. Heat Transfer. John Wiley & Sons, New York, USA, 1997, p. 1-345.

J Hirunlabh, J Santisirisomboon and P Namprakai. Assessment of solar radiation for Thailand. In: Proceedings of the International Workshop: Calculation Methods for Solar Energy System, Université de Perpignan, France, 1994, p. 1-14.

P Namprakai and J Hirunlabh. A study of hourly to daily radiation ratio correlations at Bangkok, Thailand. KMITT Res. Dev. J. 1989; 12, 40-51.

F Lasnier and WY Juen. The sizing of stand-alone photovoltaic systems using the simulation technique, RERIC. Int. Energ. J. 1990; 12, 22-3.

B Givoni. Review and evaluation of cooling by long-wave radiation. Passive Solar J. 1982; 1, 131-50.

J Waewsak, J Khedari, J Hirunlabh and R Sarachitti. The Study of Sky Emissivity in Thailand. In: Proceedings of the 15th Mechanical Engineering Network, Bangkok, Thailand, 2001, p. 32-6.

J Waewsak, J Khedari, R Sarachitti and J Hirunlabh. Direction and velocity of surface wind in Bangkok. Thammasat Int. J. Sci. Tech. 2002; 7, 56-61.

N Klompong. 2008, The Study of Mathematical Modeling of Heat Transfer in Soil. Master Thesis, Thaksin University, p. 1-102.


  • There are currently no refbacks.


Online ISSN: 2228-835X


Last updated: 17 May 2019