Heat regenerator occupied by regenerative materials improves thermal efficiency of regenerative combustion system through the recovery of heat of exhaust gaset. By using one-dimensional two-phase fluid dynamics model, the unsteady thermal flow of heat regenerator with spherical particles, was numerically simulated to evaluate the heat transfer and pressure drop and thereby to suggest the parameter for designing heat regenerator. It takes about 7 hours for the steady state of the flow field in regenerator, in which heat absorption of regenerative particle is concurrent with the same magnitude of heat desorption. The regenerative particle experiences small temperature fluctuation below 10 K during the reversing process. The performance of thermal flow in heat regenerator varies with inlet velocity of exhaust gas and air, configuration of regenerator (cross-sectional area and length) and diameter of regenerative particle. As the gas velocity increases, the heat transfer between gas and particle enhances and with the increase the pressure losses. As particle diameter decreases, the air is preheated higher and the exhaust gases are cooled more with the increase of pressure losses.