Numerical simulations of laminar premixed flames burning hydrogen and methane in spontaneous ignition mode are performed by harmonically exciting the reactants' temperature at the domain inlet. The results are compared to an analytical model representing the same reactive flow configuration. The model provides a simplified but nevertheless accurate representation of reheat combustion taking place in sequential gas turbine combustors. An analytic expression for autoignition flames transfer functions to entropy waves is derived and used to extend transfer function models from the literature. For validation purposes, results from fully compressible Direct Numerical Simulations (DNS), including a complete representation of the fluctuating acoustic and entropic fields of the reactive flow, are analyzed and compared to incompressible Unsteady Reynolds-Averaged Navier-Stokes (URANS) simulations that only take into account the fluctuating entropic field. Methane flames are found to be more sensitive to entropic forcing than hydrogen flames, featuring nonlinear phenomena even for low excitation amplitudes. In the linear regime, all flames behave as predicted by the analytical model and the URANS simulations are found to correctly predict the fluctuating entropic field. The transition from linear to nonlinear flame response is described in detail and its physical mechanisms are explained. Comparisons with results available in the literature show good prediction capabilities, both in terms of flame describing function and integrated heat release rate. Limitations of the proposed analytical model with respect to real combustion systems are discussed and a simple correction is proposed.