Numerical Simulations of Penetration and Overshoot in the Sun Tamara M. Rogers Astronomy and Astrophysics Department, University of California, Santa Cruz, CA 95064
[email protected] Gary A. Glatzmaier Earth Sciences Department, University of California, Santa Cruz, CA 95064 C.A. Jones Department of Applied Math, University of Leeds, Leeds LS2 9JT, UK ABSTRACT We present numerical simulations of convective overshoot in a two- dimensional model of the solar equatorial plane. The model equations are solved in the anelastic approximation with enhanced thermal conductivity and viscosity for numerical stability. The simulated domain extends from 0.001 R to 0.93 R, spanning both convective and radiative regions. We show that convective pen- etration leads to a slightly extended, mildly subadiabatic temperature gradient beneath the convection zone, spanning approximately 0.05 Hp,belowwhichthere is a rapid transition to a strongly subadiabatic region. A slightly higher tempera- ture is maintained in the overshoot region by adiabatic heating from overshooting plumes. This enhanced temperature may partially account for the sound speed discrepancy between the standard solar model and helioseismology. Simulations conducted with tracer particles suggest that a fully mixed region exists down to at least 0.684 R. Subject headings: convection: overshoot,mixing 1. Introduction One of the main unsolved problems in all of stellar evolution theory is the treatment of convective-radiative boundaries. In the Sun, understanding convective overshoot is cru- cial in determining properties of the solar tachocline and the solar dynamo. Helioseismic –2– observations have shown that the differential rotation observed at the solar surface persists throughout the convection zone (Thompson et al.