In this work, we investigate the impact of uncertainties due to convective boundary mixing (CBM), commonly called `overshoot', namely the boundary location and the amount of mixing at the convective boundary, on stellar structure and evolution. For this we calculated two grids of stellar evolution models with the MESA code, each with the Ledoux and the Schwarzschild boundary criterion, and vary the amount of CBM. We calculate each grid with the initial masses $15$, $20$ and $25\,\rm{M}_\odot$. We present the stellar structure of the models during the hydrogen and helium burning phases. In the latter, we examine the impact on the nucleosynthesis. We find a broadening of the main-sequence with more CBM, which is more in agreement with observations. Furthermore during the core hydrogen burning phase there is a convergence of the convective boundary location due to CBM. The uncertainties of the intermediate convective zone remove this convergence. The behaviour of this convective zone strongly affects the surface evolution of the model, i.e. how fast it evolves red-wards. The amount of CBM impacts the size of the convective cores and the nucleosynthesis, e.g. the $^{12}$C to $^{16}$O ratio and the weak s-process. Lastly, we determine the uncertainty that the range of parameter values investigated introduce and we find differences of up to $70\%$ for the core masses and the total mass of the star.

We find an asymptotic limit for the integral scale dissipation length in strongly stratified stellar convection; this is not adjustable, and we identify it with the "mixing length parameter", $\alpha \sim 5/3$, a theoretical prediction which agrees to within 10% with the fitted empirical values quoted by eight independent stellar evolutionary groups, as well as several 3D stellar atmosphere groups. For strong stratification, the dissipation length approaches the density scale height; in weak stratification it shrinks to the depth for such thin convective regions. Simulations with different zoning give a monotonic decrease in resolution errors with added zones, down to barely detectable for the most difficult case (lower boundary layer), and an estimate of the dissipation due to the turbulent cascade. No explicit viscosity is required; stellar dissipation occurs from nonlinear fluid effects (Taylor, Onsager). These implicit large eddy simulations of the Euler equations, with simultaneous analysis by Reynolds averaging (RA-ILES), resolve the energy-containing eddies, develop a turbulent cascade down to the grid scale, and produce surfaces of separation at boundary layers. Like experiments (Warhaft), and direct numerical simulations of the Navier-Stokes equations (DNS, Sreenivasan), our simulations develop "anomalous terms" for dissipation in the turbulent cascade (due to intermittency, anisotropy, and interactions between coherent structures). They have more realistic behavior at high order than the K41 theory of Kolmogorov. Mixing length theory (MLT, Boehm-Vitense) overestimates the convective enthalpy flux for strong stratification, which has consequences for the theory of luminous blue variables, for the solar metal abundance, and deep solar convection.

Renzini (1987) wrote an influential critique of mixing-length theory (MLT) as used in stellar evolution codes, and concluded that three-dimensional (3D) fluid dynamical simulations were needed to clarify several important issues. We have critically explored the limitations of the numerical methods and conclude that they are approaching the required accuracy. Implicit large eddy simulations (ILES) automatically connect large scale turbulence to a Kolmogorov cascade below the grid scale, allowing turbulent boundary layers to remove singularities that appear in the theory. Interactions between coherent structures give multi-modal behavior, driving intermittency and fluctuations. Reynolds averaging (RA) allows us to abstract the essential features of this dynamical behavior of boundaries which are appropriate to stellar evolution, and consider how they relate static boundary conditions (Richardson, Schwarzschild or Ledoux). We clarify several questions concerning when and why MLT works, and does not work, using both analytical theory and 3D high resolution numerical simulations. The composition gradients and boundary layer structure which are produced by our simulations suggest a self-consistent approach to boundary layers, removing the need for ad hoc procedures for 'convective overshooting' and `semi-convection'. In a companion paper we quantify the adequacy of our numerical resolution, determine of the length scale of dissipation (the `mixing length') without astronomical calibration, quantify agreement with the four-fifths law of Kolmogorov for weak stratification, and extend MLT to deal with strong stratification.

Convection plays a key role in the evolution of stars due to energy transport and mixing of composition. Despite its importance, this process is still not well understood. One longstanding conundrum in all 1D stellar evolution codes is the treatment of convective boundaries. In this study we compare two convective uncertainties, the boundary location (Ledoux versus Schwarzschild) and the amount of extra mixing, and their impact on the early evolution of massive stars. With increasing convective boundary mixing (CBM), we find a convergence of the two different boundary locations, a decreasing blue to red super giant ratio and a reduced importance of semiconvection.

We compare analytic predictions of supernova light curves with recent high quality data from SN2011fe (Ia), from KSN2011b (Ia), and the Palomar Transient Factory (PTF) and the La Silla-QUEST variability survey (LSQ) (Ia). Because of the steady, fast cadence of observations, KSN2011b provides unique new information on SNe Ia: the smoothness of the light curve, which is consistent with significant large-scale mixing during the explosion, possibly due to 3D effects (e.g., Rayleigh-Taylor instabilities), and provides support for a slowly-varying leakage (mean opacity). For a more complex light curve (SN2008D, SNIb), we separate the luminosity due to multiple causes and indicate the possibility of a radioactive plume. The early rise in luminosity is shown to be affected by the opacity (leakage rate) for thermal and non-thermal radiation. A general derivation of Arnett's rule again shows that it depends upon {\em all} processes heating the plasma, not just radioactive ones, so that SNe Ia will differ from SNe Ibc if the latter have multiple heating processes.

Sirius, the seventh-nearest stellar system, is a visual binary containing the metallic-line A1 V star Sirius A, brightest star in the sky, orbited in a 50.13-year period by Sirius B, the brightest and nearest white dwarf (WD). Using images obtained over nearly two decades with the Hubble Space Telescope (HST), along with photographic observations covering almost 20 years, and nearly 2300 historical measurements dating back to the 19th century, we determine precise orbital elements for the visual binary. Combined with the parallax and the motion of the A component, these elements yield dynamical masses of 2.063+/-0.023 Msun and 1.018+/-0.011 Msun for Sirius A and B, respectively. Our precise HST astrometry rules out third bodies orbiting either star in the system, down to masses of ~15-25 Mjup. The location of Sirius B in the H-R diagram is in excellent agreement with theoretical cooling tracks for WDs of its dynamical mass, and implies a cooling age of ~126 Myr. The position of Sirius B in the mass-radius plane is also consistent with WD theory, assuming a carbon-oxygen core. Including the pre-WD evolutionary timescale of the assumed progenitor, the total age of Sirius B is about 228+/-10 Myr. We calculated evolutionary tracks for stars with the dynamical mass of Sirius A, using two independent codes. We find it necessary to assume a slightly sub-solar metallicity, of about 0.85 Zsun, to fit its location in the luminosity-radius plane. The age of Sirius A based on these models is about 237-247 Myr, with uncertainties of +/-15 Myr, consistent with that of the WD companion. We discuss astrophysical puzzles presented by the Sirius system, including the probability that the two stars must have interacted in the past, even though there is no direct evidence for this, and the orbital eccentricity remains high.

Turbulent mixing of chemical elements by convection has fundamental effects on the evolution of stars. The standard algorithm at present, mixing-length theory (MLT), is intrinsically local, and must be supplemented by extensions with adjustable parameters. As a step toward reducing this arbitrariness, we compare asteroseismically inferred internal structures of two Kepler slowly pulsating B stars (SPB's; $M\sim 3.25 M_\odot$) to predictions of 321D turbulence theory, based upon well-resolved, truly turbulent three-dimensional simulations (Arnett , et al. 2015, Christini, et al. 2016) which include boundary physics absent from MLT. We find promising agreement between the steepness and shapes of the theoretically-predicted composition profile outside the convective region in 3D simulations and in asteroseismically constrained composition profiles in the best 1D models of the two SPBs. The structure and motion of the boundary layer, and the generation of waves, are discussed.

Asteroseismology of non-radial pulsations in Hot B Subdwarfs (sdB stars) offers a unique view into the interior of core-helium-burning stars. Ground-based and space-borne high precision light curves allow for the analysis of pressure and gravity mode pulsations to probe the structure of sdB stars deep into the convective core. As such asteroseismological analysis provides an excellent opportunity to test our understanding of stellar evolution. In light of the newest constraints from asteroseismology of sdB and red clump stars, standard approaches of convective mixing in 1D stellar evolution models are called into question. The problem lies in the current treatment of overshooting and the entrainment at the convective boundary. Unfortunately no consistent algorithm of convective mixing exists to solve the problem, introducing uncertainties to the estimates of stellar ages. Three dimensional simulations of stellar convection show the natural development of an overshooting region and a boundary layer. In search for a consistent prescription of convection in one dimensional stellar evolution models, guidance from three dimensional simulations and asteroseismological results is indispensable.

The numerical simulation of turbulence in stars has led to a rich set of possibilities regarding stellar pulsations, asteroseismology, thermonuclear yields, and formation of neutron stars and black holes. The breaking of symmetry by turbulent flow grows in amplitude as collapse is approached, which insures that the conditions at the onset of collapse are not spherical. This lack of spherical symmetry has important implications for the mechanism of explosion and ejected nucleosynthesis products. Numerical resolution of several different types of three--dimensional (3D) stellar simulations are compared; it is suggested that core collapse simulations may be under-resolved. New physical effects which appear in 3D are summarized. Connections between simulations of progenitor explosion and observations of supernova remnants (SNR) are discussed. Present treatment of boundaries, for mixing regions during He--burning, requires revision.

We explore the effects of rotation on convective carbon, oxygen, and silicon shell burning during the late stages of evolution in a 20Msun star. Using the Modules for Experiments in Stellar Astrophysics (MESA) we construct 1D stellar models both with no rotation and with an initial rigid rotation of 50% of critical. At different points during the evolution, we map the 1D models into 2D and follow the multidimensional evolution using the FLASH compressible hydrodynamics code for many convective turnover times until a quasi-steady state is reached. We characterize the strength and scale of convective motions via decomposition of the momentum density into vector spherical harmonics. We find that rotation influences the total power in solenoidal modes, with a slightly larger impact for carbon and oxygen shell burning than for silicon shell burning. Including rotation in one-dimensional (1D) stellar evolution models alters the structure of the star in a manner that has a significant impact on the character of multidimensional convection. Adding modest amounts of rotation to a stellar model that ignores rotation during the evolutionary stage, however, has little impact on the character of resulting convection. Since the spatial scale and strength of convection present at the point of core collapse directly influence the supernova mechanism, our results suggest that rotation could play an important role in setting the stage for massive stellar explosions.

Stellar models provide a vital basis for many aspects of astronomy and astrophysics. Recent advances in observational astronomy -- through asteroseismology, precision photometry, high-resolution spectroscopy, and large-scale surveys -- are placing stellar models under greater quantitative scrutiny than ever. The model limitations are being exposed and the next generation of stellar models is needed as soon as possible. The current uncertainties in the models propagate to the later phases of stellar evolution, hindering our understanding of stellar populations and chemical evolution. Here we give a brief overview of the evolution, importance, and substantial uncertainties of core helium burning stars in particular and then briefly discuss a range of methods, both theoretical and observational, that we are using to advance the modelling.

We present the first three dimensional (3D) simulation of the final minutes of iron core growth in a massive star, up to and including the point of core gravitational instability and collapse. We self-consistently capture the development of strong convection driven by violent Si burning in the shell surrounding the iron core. This convective burning builds the iron core to its critical (Chandrasekhar) mass and collapse ensues, driven by electron capture and photodisintegration. The non-spherical structure and motion (turbulent fluctuations) generated by 3D convection is substantial at the point of collapse. We examine the impact of such physically-realistic 3D initial conditions on the core-collapse supernova mechanism using 3D simulations including multispecies neutrino leakage. We conclude that non-spherical progenitor structure should not be ignored, and has a significant and favorable impact on the likelihood for neutrino-driven explosions.

We examine the physical basis for algorithms to replace mixing-length theory (MLT) in stellar evolutionary computations. Our 321D procedure is based on numerical solutions of the Navier-Stokes equations. These implicit large eddy simulations (ILES) are three-dimensional (3D), time-dependent, and turbulent, including the Kolmogorov cascade. We use the Reynolds-averaged Navier-Stokes (RANS) formulation to make concise the 3D simulation data, and use the 3D simulations to give closure for the RANS equations. We further analyze this data set with a simple analytical model, which is non-local and time-dependent, and which contains both MLT and the Lorenz convective roll as particular subsets of solutions. A characteristic length (the damping length) again emerges in the simulations; it is determined by an observed balance between (1) the large-scale driving, and (2) small-scale damping. The nature of mixing and convective boundaries is analyzed, including dynamic, thermal and compositional effects, and compared to a simple model. We find that (1) braking regions (boundary layers in which mixing occurs) automatically appear {\it beyond} the edges of convection as defined by the Schwarzschild criterion, (2) dynamic (non-local) terms imply a non-zero turbulent kinetic energy flux (unlike MLT), (3) the effects of composition gradients on flow can be comparable to thermal effects, and (4) convective boundaries in neutrino-cooled stages differ in nature from those in photon-cooled stages (different P\'eclet numbers). The algorithms are based upon ILES solutions to the Navier-Stokes equations, so that, unlike MLT, they do not require any calibration to astronomical systems in order to predict stellar properties. Implications for solar abundances, helioseismology, asteroseismology, nucleosynthesis yields, supernova progenitors and core collapse are indicated.

Stellar evolution calculations have had great success reproducing the observed atmospheric properties of different classes of stars. Recent detections of g-mode pulsations in evolved He burning stars allow a rare comparison of their internal structure with stellar models. Asteroseismology of subdwarf B stars suggests convective cores of $0.22-0.28\,M_\odot$, $\gtrsim 45\,\%$ of the total stellar mass. Previous studies found significantly smaller convective core masses ($\lesssim 0.19\,M_\odot$) at a comparable evolutionary stage. We evolved stellar models with MESA (Modules for Experiments in Stellar Astrophysics) to explore how well the interior structure inferred from asteroseismology can be reproduced by standard algorithms. Our qualitative evolutionary paths, position in the $\log g-T_{\rm{eff}}$ diagram and model timescales are consistent with previous results. SdB masses from our full evolutionary sequences fall within the range of the empirical sdB mass distribution, but are nearly always lower than the median. Using standard MLT with atomic diffusion we find convective core masses of $\sim 0.17-0.18 M_\odot$, averaged over the entire sdB lifetime. We can increase the convective core sizes to be as large as those inferred from asteroseismology, but only for extreme values of the overshoot parameter (overshoot gives numerically unstable and physically unrealistic behavior at the boundary). High resolution three-dimensional (3D) simulations of turbulent convection in stars suggest that the Schwarzschild criterion for convective mixing sytematically {\it underestimates the actual extent of mixing because a boundary layer forms.} Accounting for this would decrease the errors in both sdB total and convective core masses.

This is a sketch of a 321D approximation for stellar convection which is nonlocal, and thus has nonzero fluxes of KE (to be published in more detail elsewhere). Boundary conditions are discussed in a fluid dynamics context (i.e., predictions for overshoot, semiconvection and entrainment are analyzed). We plan to add this as an option to MESA. Inclusion of KE fluxes seems to help resolve the solar abundance problem (Asplund 2009). Smaller cores may ease the explosion problems with core collapse supernova simulations.

Three-dimensional (3D), time dependent numerical simulations, of flow of matter in stars, now have sufficient resolution to be fully turbulent. The late stages of the evolution of massive stars, leading up to core collapse to a neutron star (or black hole), and often to supernova explosion and nucleosynthesis, are strongly convective because of vigorous neutrino cooling and nuclear heating. Unlike models based on current stellar evolutionary practice, these simulations show a chaotic dynamics characteristic of highly turbulent flow. Theoretical analysis of this flow, both in the Reynolds-averaged Navier-Stokes (RANS) framework and by simple dynamic models, show an encouraging consistency with the numerical results. It may now be possible to develop physically realistic and robust procedures for convection and mixing which (unlike 3D numerical simulation) may be applied throughout the long life times of stars. In addition, a new picture of the presupernova stages is emerging which is more dynamic and interesting (i.e., predictive of new and newly observed phenomena) than our previous one.

Symbiotic binaries are systems containing white dwarfs (WDs) and red giants. Symbiotic novae are those systems in which thermonuclear eruptions occur on the WD components. These are to be distinguished from events driven by accretion disk instabilities analogous to dwarf novae eruptions in cataclysmic variable outbursts. Another class of symbiotic systems are those in which the WD is extremely luminous and it seems likely that quiescent nuclear burning is ongoing on the accreting WD. A fundamental question is the secular evolution of the WD. Do the repeated outbursts or quiescent burning in these accreting systems cause the WD to gain or lose mass? If it is gaining mass, can it eventually reach the Chandrasekhar Limit and become a supernova (a SN Ia if it can hide the hydrogen and helium in the system)? In order to better understand these systems, we have begun a new study of the evolution of Thermonuclear Runaways (TNRs) in the accreted envelopes of WDs using a variety of initial WD masses, luminosities and mass accretion rates. We use our 1-D hydro code, NOVA, which includes the new convective algorithm of Arnett, Meakin and Young, the Hix and Thielemann nuclear reaction solver, the Iliadis reaction rate library, the Timmes equation of state, and the OPAL opacities. We assume a solar composition (Lodders abundance distribution) and do not allow any mixing of accreted material with core material. This assumption strongly influences our results. We report here (1) that the WD grows in mass for all simulations so that canonical `steady burning' does not occur, and (2) that only a small fraction of the accreted matter is ejected in some (but not all) simulations. We also find that the accreting systems, before thermonuclear runaway, are too cool to be seen in X-ray searches for SN Ia progenitors.

SBW1 is a B supergiant surrounded by a ring nebula that is a nearby twin of SN 1987A. We present images and spectra of SBW1 obtained with HST, Spitzer, and Gemini South. HST images of SBW1 do not exhibit long Rayleigh-Taylor fingers, which are presumed to cause the hotspots in the SN1987A ring, but instead show a geometrically thin clumpy ring. The radial mass distribution and size scales of inhomogeneities in SBW1's ring closely resemble those in the SN1987A ring, but the more complete disk expected to reside at the base of the RT fingers is absent. This structure may explain why portions of the SN1987A ring between the hotspots have not yet brightened after more than 15 years. The model we suggest does not require a fast wind colliding with a previous red supergiant wind. More surprisingly, images of SBW1 also reveal diffuse emission filling the interior of the ring in H-alpha and thermal-IR emission; 190K dust dominates the 8-20 micron luminosity (but contains only 1e-5Msun). Cooler (85K dust resides in the equatorial ring itself (dust mass of 5e-3Msun). Diffuse emission extends inward to 1 arcsecond from the central star, where a paucity of emission suggests an inner hole excavated by the wind. We propose that diffuse emission inside the ring arises from an ionized flow of material photoevaporated from the dense ring, and it prevents the supergiant wind from advancing in the equator. This inner emission could correspond to a structure hypothesized to reside around SN1987A that was never directly detected. We suggest that photoionization can play an important dynamical role in shaping the ring nebula, and we speculate that this might help explain the origin of the polar rings around SN1987A. abridged.

We present a brief summary of the Single Degenerate Scenario for the progenitors of Type Ia Supernovae in which it is assumed that a low mass carbon-oxygen white dwarf is growing in mass as a result of accretion from a secondary star in a close binary system. Recent hydrodynamic simulations of accretion of solar material onto white dwarfs without mixing always produce a thermonuclear runaway and steady burning does not occur. For a broad range in WD mass (0.4 Solar masses to 1.35 Solar Masses), the maximum ejected material occurs for the 1.25 Solar Mass sequences and then decreases as the white dwarf mass decreases. Therefore, the white dwarfs are growing in mass as a consequence of the accretion of solar material and as long as there is no mixing of accreted material with core material. In contrast, a thermonuclear runaway in the accreted hydrogen-rich layers on the low luminosity WDs in close binary systems where mixing of core matter with accreted material has occurred is the outburst mechanism for Classical, Recurrent, and Symbiotic novae. The differences in characteristics of these systems is likely the WD mass and mass accretion rate. The high levels of enrichment of CN ejecta in elements ranging from carbon to sulfur confirm that there is dredge-up of matter from the core of the WD and enable them to contribute to the chemical enrichment of the interstellar medium. Therefore, studies of CNe can lead to an improved understanding of Galactic nucleosynthesis, some sources of pre-solar grains, and the Extragalactic distance scale. The characteristics of the outburst depend on the white dwarf mass, luminosity, mass accretion rate, and the chemical composition of both the accreting material and WD material. The properties of the outburst also depends on when, how, and if the accreted layers are mixed with the WD core and the mixing mechanism is still unknown.

Three-dimensional (3D) hydrodynamic simulations of shell oxygen burning (Meakin and Arnett, 2007b) exhibit bursty, recurrent fluctuations in turbulent kinetic energy. These are shown to be due to a general instability of the convective cell, requiring only a localized source of heating or cooling. Such fluctuations are shown to be suppressed in simulations of stellar evolution which use mixing-length theory (MLT). Quantitatively similar behavior occurs in the model of a convective roll (cell) of Lorenz (1963), which is known to have a strange attractor that gives rise to chaotic fluctuations in time of velocity and, as we show, luminosity. Study of simulations suggests that the behavior of a Lorenz convective roll may resemble that of a cell in convective flow. We examine some implications of this simplest approximation, and suggest paths for improvement. Using the Lorenz model as representative of a convective cell, a multiple-cell model of a convective layer gives total luminosity fluctuations which are suggestive of irregular variables (red giants and supergiants (Schwarzschild 1975)), and of the long secondary period feature in semi-regular AGB variables (Stothers 2010, Wood, Olivier and Kawaler 2004). This "tau-mechanism" is a new source for stellar variability, which is inherently non-linear (unseen in linear stability analysis), and one closely related to intermittency in turbulence. It was already implicit in the 3D global simulations of Woodward, Porter and Jacobs (2003). This fluctuating behavior is seen in extended 2D simulations of CNeOSi burning shells (Arnett and Meakin 2011b), and may cause instability which leads to eruptions in progenitors of core collapse supernovae PRIOR to collapse.

Two-dimensional (2D) hydrodynamical simulations of progenitor evolution of a 23 solar mass star, close to core collapse (about 1 hour, in 1D), with simultaneously active C, Ne, O, and Si burning shells, are presented and contrasted to existing 1D models (which are forced to be quasi-static). Pronounced asymmetries, and strong dynamical interactions between shells are seen in 2D. Although instigated by turbulence, the dynamic behavior proceeds to sufficiently large amplitudes that it couples to the nuclear burning. Dramatic growth of low order modes is seen, as well as large deviations from spherical symmetry in the burning shells. The vigorous dynamics is more violent than that seen in earlier burning stages in the 3D simulations of a single cell in the oxygen burning shell, or in 2D simulations not including an active Si shell. Linear perturbative analysis does not capture the chaotic behavior of turbulence (e.g., strange attractors such as that discovered by Lorenz), and therefore badly underestimates the vigor of the instability. The limitations of 1D and 2D models are discussed in detail. The 2D models, although flawed geometrically, represent a more realistic treatment of the relevant dynamics than existing 1D models, and present a dramatically different view of the stages of evolution prior to collapse. Implications for interpretation of SN1987A, abundances in young supernova remnants, pre-collapse outbursts, progenitor structure, neutron star kicks, and fallback are outlined. While 2D simulations provide new qualitative insight, fully 3D simulations are needed for a quantitative understanding of this stage of stellar evolution. The necessary properties of such simulations are delineated.

Three-dimensional (3D) hydrodynamic simulations of shell oxygen burning (Meakin and Arnett 2007) exhibit bursty, recurrent fluctuations in turbulent kinetic energy. These are shown to be due to a global instability in the convective region, which has been suppressed in calculations of stellar evolution which use mixing-length theory (MLT). Quantitatively similar behavior occurs in the model of a convective roll (cell) of Lorenz (1963), which is known to have a strange attractor that gives rise to random fluctuations in time.An extension of the Lorenz model, which includes Kolmogorov damping and nuclear burning, is shown to exhibit bursty, recurrent fluctuations like those seen in the 3D simulations. A simple model of a convective layer (composed of multiple Lorenz cells) gives luminosity fluctuations which are suggestive of irregular variables (red giants and supergiants, Schwarzschild 1975). Apparent inconsistencies between Arnett, Meakin, and Young (2009) and Nordlund, Stein, and Asplund (2009) on the nature of convective driving have been resolved, and are discussed.

A program is outlined, and first results described, in which fully three-dimensional, time dependent simulations of hydrodynamic turbulence are used as a basis for theoretical investigation of the physics of turbulence in stars. The inadequacy of the treatment of turbulent convection as a diffusive process is discussed. A generalization to rotation and magnetohydrodynamics is indicated, as are connection to simulations of 3D stellar atmospheres.

As a step toward a complete theoretical integration of 3D compressible hydrodynamic simulations into stellar evolution, convection at the surface and sub-surface layers of the Sun is re-examined, from a restricted point of view, in the language of mixing-length theory (MLT) . Requiring that MLT use a hydrodynamically realistic dissipation length gives a new constraint on solar models. While the stellar structure which results is similar to that obtained by YREC and Garching models, the theoretical picture differs. A new quantitative connection is made between macro-turbulence, micro-turbulence, and the convective velocity scale at the photosphere, which has finite values. The "geometric parameter" in MLT is found to correspond more reasonably with the size of the strong downward plumes which drive convection (Stein and Nordlund 1998), and thus has a physical interpretation even in MLT. Use of 3D simulations of both adiabatic convection and stellar atmospheres will allow the determination of the dissipation length and the geometric parameter (i.e., the entropy jump), with no astronomical calibration. A physically realistic treatment of convection in stellar evolution will require additional modifications beyond MLT, including effects of kinetic energy flux, entrainment (the most dramatic difference from MLT found by Meakin and Arnett 2007), rotation, and magnetic fields (Balbus 2009}.