Acoustically forced lean premixed turbulent bluff-body stabilized flames are investigated using turbulent combustion CFD. The calculations simulate aspects of the experimental investigation by Balachandran et al. [R. Balachandran, B. Ayoola, C. Kaminski, A. Dowling, E. Mastorakos, Combust. Flame 143 (2005) 37-55] and focus on the amplitude dependence of the flame response. For the frequencies of interest in this investigation an unsteady Reynolds-averaged Navier-Stokes (URANS) approach is appropriate. The combustion is represented using a modified laminar flamelet approach with an algebraic representation of the flame surface density. The predictions are compared with flame surface density (FSD) and OH* chemiluminescence measuremems. In the experiments the response of the flame has been quantified by means of a number of single-frequency, amplitude-dependent transfer functions. The predicted flame shape and position are in good agreement with the experiment. The dynamic response of the flame to inlet velocity forcing is also well captured by the calculations. At moderate frequencies nonlinear behavior of the transfer functions is observed as the forcing amplitude is increased. In the experiments this nonlinearity was attributed in part to the rollup of the reacting shear layer into vortices and in part to the collision of the inner and outer flame sheets. This transition to nonlinearity is also observed in the transfer functions obtained from the predictions. Furthermore, the vortex shedding and flame-sheet collapse may be seen in snapshots of the predicted flow field taken throughout the forcing cycle. The URANS methodology successfully predicts the behavior of the forced premixed turbulent flames and captures the effects of saturation in the transfer function of the response of the heat release to velocity fluctuations. (c) 2006 The Combustion Institute. Published by Elsevier Inc. All rights reserved.