For the numerical evaluation of the thermoacoustic stability of rocket engines often hybrid methods are applied, which separate the computation of wave propagation in the combustor from the analysis of the flame response to acoustic perturbations. Closure requires a thermoacoustic feedback model which provides the heat release fluctuation in the source term of the employed wave transport equations. The influence of the acoustic fluctuations in the combustion chamber on the heat release fluctuations from the modulation of the atomization of the propellants in a hypergolic upper stage rocket engine is studied. Numerical modeling of a single injector provides the time mean reacting flow field. A network of transfer functions representing all aspects relevant for the feedback model is presented. Analytical models for the injector admittances and for the atomization transfer functions are provided. The dynamics of evaporation and combustion are studied numerically and the numerical results are analyzed. An analytical approximation of the computed flame transfer function is combined with the analytical models for the injector and the atomization quality to derive the feedback model for the wave propagation code. The evaluation of this model on the basis of the Rayleigh index reveals the thermoacoustic driving potential originating from the fluctuating spray quality.