Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userA two-dimensional mixing length theory of convective transport
423
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
The helioseismic observations of the internal rotation profile of the Sun raise questions about the two-dimensional (2D) nature of the transport of angular momentum in stars. Here we derive a convective prescription for axisymmetric (2D) stellar evolution models. We describe the small scale motions by a spectrum of unstable linear modes in a Boussinesq fluid. Our saturation prescription makes use of the angular dependence of the linear dispersion relation to estimate the anisotropy of convective velocities. We are then able to provide closed form expressions for the thermal and angular momentum fluxes with only one free parameter, the mixing length. We illustrate our prescription for slow rotation, to first order in the rotation rate. In this limit, the thermodynamical variables are spherically symetric, while the angular momentum depends both on radius and latitude. We obtain a closed set of equations for stellar evolution, with a self-consistent description for the transport of angular momentum in convective regions. We derive the linear coefficients which link the angular momentum flux to the rotation rate ($Lambda$- effect) and its gradient ($alpha$-effect). We compare our results to former relevant numerical work.
rate research
Read More
Stellar convection is customarily described by Mixing-Length Theory, which makes use of the mixing-length scale to express the convective flux, velocity, and temperature gradients of the convective elements and stellar medium. The mixing-length scale is taken to be proportional to the local pressure scale height, and the proportionality factor (the mixing-length parameter) must be determined by comparing the stellar models to some calibrator, usually the Sun. No strong arguments exist to suggest that the mixing-length parameter is the same in all stars and at all evolutionary phases. The aim of this study is to present a new theory of stellar convection that does not require the mixing length parameter. We present a self-consistent analytical formulation of stellar convection that determines the properties of stellar convection as a function of the physical behaviour of the convective elements themselves and of the surrounding medium. This new theory is formulated starting from a conventional solution of the Navier-Stokes/Euler equations, i.e. the Bernoulli equation for a perfect fluid, but expressed in a non-inertial reference frame co-moving with the convective elements. In our formalism the motion of stellar convective cells inside convectively-unstable layers is fully determined by a new system of equations for convection in a non-local and time-dependent formalism. We obtain an analytical, non-local, time-dependent sub-sonic solution for the convective energy transport that does not depend on any free parameter. The theory is suitable for the outer convective zones of solar type stars and stars of all mass on the main sequence band. The predictions of the new theory are compared with those from the standard mixing-length paradigm for the most accurate calibrator, the Sun, with very satisfactory results.
Turbulent convection is certainly one of the most important and thorny issues in stellar physics. Our deficient knowledge of this crucial physical process introduces a fairly large uncertainty concerning the internal structure and evolution of stars. A striking example is overshoot at the edge of convective cores. Indeed, nearly all stellar evolutionary codes treat the overshooting zones in a very approximative way that considers both its extent and the profile of the temperature gradient as free parameters. There are only a few sophisticated theories of stellar convection such as Reynolds stress approaches, but they also require the adjustment of a non-negligible number of free parameters. We present here a theory, based on the plume theory as well as on the mean-field equations, but without relying on the usual Taylors closure hypothesis. It leads us to a set of eight differential equations plus a few algebraic ones. Our theory is essentially a non-mixing length theory. It enables us to compute the temperature gradient in a shrinking convective core and its overshooting zone. The case of an expanding convective core is also discussed, though more briefly. Numerical simulations have quickly improved during recent years and enabling us to foresee that they will probably soon provide a model of convection adapted to the computation of 1D stellar models.
Classical novae are thermonuclear explosions that take place in the envelopes of accreting white dwarfs in binary systems. The material piles up under degenerate conditions, driving a thermonuclear runaway. The energy released by the suite of nuclear processes operating at the envelope heats the material up to peak temperatures about 100 - 400 MK. During these events, about 10-3 - 10-7 Msun, enriched in CNO and, sometimes, other intermediate-mass elements (e.g., Ne, Na, Mg, Al) are ejected into the interstellar medium. To account for the gross observational properties of classical novae (in particular, the large concentrations of metals spectroscopically inferred in the ejecta), models require mixing between the (solar-like) material transferred from the secondary and the outermost layers (CO- or ONe-rich) of the underlying white dwarf. Recent multidimensional simulations have demonstrated that Kelvin-Helmholtz instabilities can naturally produce self-enrichment of the accreted envelope with material from the underlying white dwarf at levels that agree with observations. However, the feasibility of this mechanism has been explored in the framework of CO white dwarfs, while mixing with different substrates still needs to be properly addressed. Three-dimensional simulations of mixing at the core-envelope interface during nova outbursts have been performed with the multidimensional code FLASH, for two types of substrates: CO- and ONe-rich. We show that the presence of an ONe-rich substrate, as in neon novae, yields larger metallicity enhancements in the ejecta, compared to CO,rich substrates (i.e., non-neon novae). A number of requirements and constraints for such 3-D simulations (e.g., minimum resolution, size of the computational domain) are also outlined.
The $s$-process nucleosynthesis in Asymptotic Giant Branch (AGB) stars depends on the modeling of convective boundaries. We present models and s-process simulations that adopt a treatment of convective boundaries based on the results of hydrodynamic simulations and on the theory of mixing due to gravity waves in the vicinity of convective boundaries. Hydrodynamics simulations suggest the presence of convective boundary mixing (CBM) at the bottom of the thermal pulse-driven convective zone. Similarly, convection-induced mixing processes are proposed for the mixing below the convective envelope during third dredge-up where the 13C pocket for the s process in AGB stars forms. In this work we apply a CBM model motivated by simulations and theory to models with initial mass $M = 2$ and $M = 3M_odot$, and with initial metal content Z = 0.01 and Z = 0.02. As reported previously, the He-intershell abundance of 12C and 16O are increased by CBM at the bottom of pulse-driven convection zone. This mixing is affecting the $^{22}Ne(alpha,n)^{25}Mg$ activation and the s-process effciency in the 13C-pocket. In our model CBM at the bottom of the convective envelope during the third dredgeup represents gravity wave mixing. We take further into account that hydrodynamic simulations indicate a declining mixing efficiency already about a pressure scale height from the convective boundaries, compared to mixing-length theory. We obtain the formation of the 13C-pocket with a mass of $approx 10^{-4}M_odot$. The final $s$-process abundances are characterized by 0.36 < [s=Fe] < 0.78 and the heavy-to-light s-process ratio is 0.23 < [hs=ls] < 0.45. Finally, we compare our results with stellar observations, pre-solar grain measurements and previous work.
The treatment of mixing processes is still one of the major uncertainties in 1D stellar evolution models. This is mostly due to the need to parametrize and approximate aspects of hydrodynamics in hydrostatic codes. In particular, the effect of hydrodynamic instabilities in rotating stars, for example, dynamical shear instability, evades consistent description. We intend to study the accuracy of the diffusion approximation to dynamical shear in hydrostatic stellar evolution models by comparing 1D models to a first-principle hydrodynamics simulation starting from the same initial conditions. We chose an initial model calculated with the stellar evolution code GENEC that is just at the onset of a dynamical shear instability but does not show any other instabilities (e.g., convection). This was mapped to the hydrodynamics code SLH to perform a 2D simulation in the equatorial plane. We compare the resulting profiles in the two codes and compute an effective diffusion coefficient for the hydro simulation. Shear instabilities develop in the 2D simulation in the regions predicted by linear theory to become unstable in the 1D model. Angular velocity and chemical composition is redistributed in the unstable region, thereby creating new unstable regions. Eventually the 2D simulation settles in a symmetric, steady state, which is Richardson stable everywhere, whereas the instability remains for longer in the 1D model due to current limitations in the 1D code. A spatially resolved diffusion coefficient is extracted by comparing the initial and final profiles of mean atomic mass. The presented simulation gives a first insight on hydrodynamics of shear instabilities in a real stellar environment and even allows us to directly extract an effective diffusion coefficient. We see evidence for a critical Richardson number of 0.25 as regions above this value remain stable for the course of the simulation.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
