2006/11/09 by Casey Meakin, Casey A. Meakin, David Arnett · 2 citations
Physics and Astronomy · #Astro and Planetary Science #Astrophysics and Star Formation Studies #Stellar, planetary, and galactic studies #astro-ph
paper · pdf · doi:10.1086/520318
published as Astrophys.J.667:448-475,2007 · 32 pages(emulateapj), 26 figures, 6 tables. Submitted to ApJ. Comments will be warmly welcomed. For version with full resolution figures and related pre-prints goto http://celsus.as.arizona.edu/preprints/
arxiv created 2006/11/09 · openalex publication_date 2007/09/17 · arxiv updated 2011/02/11 · openalex created_date 2025/10/10 · openalex updated_date 2026/07/31
We describe the results of 3D numerical simulations of oxygen shell burning and hydrogen core burning in a 23 M ☉ stellar model. A detailed comparison is made to stellar mixing-length theory (MLT) for the shell-burning model. Simulations in 2D are significantly different from 3D, in terms of both flow morphology and velocity amplitude. Convective mixing regions are better predicted using a dynamic boundary condition based on the bulk Richardson number than by purely local, static criteria like Schwarzschild or Ledoux. MLT gives a good description of the velocity scale and temperature gradient for shell convection; however, there are other important effects that it does not capture, mostly related to the dynamical motion of the boundaries between convective and nonconvective regions. There is asymmetry between upflows and downflows, so the net kinetic energy flux is not zero. The motion of convective boundaries is a source of gravity waves; this is a necessary consequence of the deceleration of convective plumes. Convective "overshooting" is best described as an elastic response by the convective boundary, rather than ballistic penetration of the stable layers by turbulent eddies. The convective boundaries are rife with internal and interfacial wave motions, and a variety of instabilities arise that induce mixing through a process best described as turbulent entrainment. We find that the rate at which material entrainment proceeds at the boundaries is consistent with analogous laboratory experiments and simulation and observation of terrestrial atmospheric mixing. In particular, the normalized entrainment rate E = u E /σ H is well described by a power-law dependence on the bulk Richardson number Ri B = Δ bL /σ for the conditions studied, 20 ≲ Ri B ≲ 420. We find E = A Ri , with best-fit values log A = 0.027 ± 0.38 and n = 1.05 ± 0.21. We discuss the applicability of these results to stellar evolution calculations.