A 2D axisymmetric numerical analysis was performed to study the characteristics of charge process inside solar thermal storage tank. The interfacial area density and inertial resistance of filler material are selected as simulation parameters. The interfacial area density is varied as 800, 2000, and 4000 1/m. The inertial resistance is varied as 1, 3, and 5 1/m. When the interfacial area density increases from 800 to 4000 1/m, the discrepancy of the temperature distributions between the filler and heat transfer fluid decreases. As inertial resistance increases from 1 to 5, both of the temperature and fluid flow pattern changes considerably.