There is a wide consensus in the astrophysics community that the mechanism underlying the observed Classical Nova eruptions is a surface thermonuclear runaway. We start this short review with the main observational facts that lead to the theoretical model of a thermonuclear runaway that takes place in an accreted hydrogen rich envelope placed on top of a cool degenerate core of a white dwarf. According to the theory, the accreted envelope becomes unstable to convection days to weeks prior to the runaway. During the extreme stages of the runaway itself, when the burning is most efficient, the envelope is fully convective. Therefore, the elements processed under such extreme conditions are lifted to the outermost regions of the star. A significant fraction of the envelope is ejected during the outburst. The complicated combination of hydrodynamic instabilities and explosive hydrogen burning, close to the surface of the star, gives us a unique opportunity to study this complex reactive flow. The range of core masses, core temperatures and accretion rates introduce a whole range of burning temperatures and densities. Following the description of the "standard" cases, we then focus on rare, but still possible, portions of the relevant parameter space, in which "breakout" of the traditional CNO cycle can occur and lead to heavy element enrichment patterns caused only by breakout burning. We conclude our review with the main challenges that nova theorists face today, with special emphasis on problems related to the nucleosynthesis issues.

For very slow white dwarf accretors in CV's Townsley and Bildsten (2004) found a relation between the accretion rate and the central temperature of the white dwarf Tc. According to this relation for accretion rates less than 10^-10 solar masses per year Tc is much lower than 10^7 K. Motivated by this study we follow the thermonuclear runaway on massive white dwarfs (M_WD=1.25 - 1.40 solar masses) with Tc lower than 10^7 K, accreting matter of solar composition. We demonstrate that in that range of the relevant parameter space (Tc,M_WD and accretion rate) the slope of the relation between the peak temperatures achieved during the runaway and Tc becomes much steeper than its value for Tc above 10^7 K. The peak temperatures we derive can lead to nuclear breakout from the conventional "hot carbon-nitrogen-oxygen" cycle. When breakout conditions are achieved the heavy element abundances can show a much wider variety than what is possible with the common enrichment mechanisms.

Multidimensional reactive flow models of accreted hydrogen rich envelopes on top of degenerate cold white dwarfs are very effective tools for the study of critical, non spherically symmetric, behaviors during the early stages of nova outbursts. Such models can shed light both on the mechanism responsible for the heavy element enrichment observed to characterize nova envelope matter and on the role of perturbations during the early stages of ignition of the runaway. The complexity of convective reactive flow in multi-dimensions makes the computational model itself complex and sensitive to the details of the numerics. In this study, we demonstrate that the imposed outer boundary condition can have a dramatic effect on the solution. Several commonly used choices for the outer boundary conditions are examined. It is shown that the solutions obtained from Lagrangian simulations, where the envelope is allowed to expand and mass is being conserved, are consistent with spherically symmetric solutions. In Eulerian schemes which utilize an outer boundary condition of free outflow, the outburst can be artificially quenched.