Convection is one of the most problematic components of almost every model of a pulsating star. RR Lyrae stars, in particular, have a narrow convection zone at the outer layers, and a model without convection cannot get the red edge of the instability strip (see for example Kollath et al., 2000).

In recent years, the common approach is to use non-local time-dependent extensions to the Mixing Length Theory of convection (Bono et al., 1997; Feuchtinger, 1999). These models have many free parameters that need to be calibrated. In this paper we report on a preliminary work in which a different approach is used: 2D hydro simulations similar to the red giant model in Asida (2000). This approach was used by Deupree (1975, 1977) but we used a much better numerical resolution and focused on comparison with MLT results.