This paper deals with approximate integral solutions to the one-dimensional model describing the charging process of stratified thermal storage tanks. Temperature is assumed to be the form of Fermi-Dirac distribution function, which can be separated to two sets of cubic polynomials for each hot and cold side of thermal boundary layers. Proposed approximate integral solutions are compared to the previous works of the approximate analytic solutions and show reasonable agreement. The approach, however, has benefits in mathematical difficulties, complicated solution form and unstable convergence of series solution founded in the previous analytic solutions. Solutions for a semi-infinite region, which have simple closed form solutions, give close agreement to those for a finite region. Thermocline thickness is obtained in closed form and shows proportional behavior to the square root of time and inverse proportional behavior to the square root of flow rate.