The nature of carbon burning flames in Type Ia supernovae is explored as they interact with Kolmogorov turbulence. One-dimensional calculations using the Linear Eddy Model of Kerstein (1991) elucidate three regimes of turbulent burning. In the simplest case, large scale turbulence folds and deforms thin laminar flamelets to produce a flame brush with a total burning rate given approximately by the speed of turbulent fluctuations on the integral scale, U_L. This is the regime where the supernova explosion begins and where most of its pre-detonation burning occurs. As the density declines, turbulence starts to tear the individual flamelets, making broader structures that move faster. For a brief time, these turbulent flamelets are still narrow compared to their spacing and the concept of a flame brush moving with an overall speed of U_L remains valid. However, the typical width of the individual flamelets, which is given by the condition that their turnover time equals their burning time, continues to increase as the density declines. Eventually, mixed regions almost as large as the integral scale itself are transiently formed. At that point, a transition to detonation can occur. The conditions for such a transition are explored numerically and it is estimated that the transition will occur for densities near 1 x 10**7 g/cm**3, provided the turbulent speed on the integral scale exceeds about 15% sonic. An example calculation shows the details of a detonation actually developing.

We present a statistical model which shows the influence of turbulence on a thermonuclear flame propagating in C+O white dwarf matter. Based on a Monte Carlo description of turbulence, it provides a method for investigating the physics in the so-called distributed burning regime. Using this method we perform numerical simulations of turbulent flames and show that in this particular regime the flamelet model for the turbulent flame velocity loses its validity. In fact, at high turbulent intensities burning in the distributed regime can lead to a deceleration of the turbulent flame and thus induces a competing process to turbulent effects that cause a higher flame speed. It is also shown that in dense C+O matter turbulent heat transport is described adequately by the Peclet number, Pe, rather than by the Reynolds number, which means that flame propagation is decoupled from small-scale turbulence. Finally, at the onset of our results we argue that the available turbulent energy in an exploding C+O white dwarf is probably too low in order to make a deflagration to detonation transition possible.

Gibson scaling and related properties of flame-surface geometry in turbulent premixed combustion are demonstrated using a novel computational model, Deterministic Turbulent Mixing (DTM). In DTM, turbulent advection is represented by a sequence of maps applied to the computational domain. The structure of the mapping sequence incorporates pertinent scaling properties of the turbulent cascade. Here, combustion in Kolmogorov turbulence (kinetic-energy cascade) and in Bolgiano-Obukhov convective turbulence (potential-energy cascade) is simulated. Implications with regard to chemical flames and astrophysical (thermonuclear) flames are noted.