Several schemes of energy minimization of drying process including intermittent drying have been attempted. Intermittent drying is conducted by applying different heat inputs in each drying period. An effective and physically meaningful drying model is useful for process design and product technology. The lumped reaction engineering approach (L-REA) has been shown previously to be accurate to model the intermittent drying In L-REA, the REA (reaction engineering approach) is used to describe the global drying rate. In this study, the REA is used to model the local evaporation/condensation rate and combined with the mechanistic drying models to yield the spatial reaction engineering approach (S-REA), a non-equilibrium multiphase drying model. The accuracy of the S-REA to model the intermittent drying under time-varying drying air temperature is evaluated here. In order to incorporate the effect of time-varying drying air temperature, the equilibrium activation energy and boundary condition of heat balance implement the corresponding drying settings in each drying period. The results of modeling using the S-REA match well with the experimental data. The S-REA can yield the spatial profiles of moisture content, concentration of water vapor, temperature and local evaporation/condensation rate so that better understanding of transport phenomena of intermittent drying can be obtained. It is argued here that the REA can describe the local evaporation rate under time-varying external conditions well. The S-REA is an effective non-equilibrium multiphase approach for modeling of intermittent drying process. (C) 2013 Elsevier B.V. All rights reserved.