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.

Characterizing stellar convection in multiple dimensions is a topic at the forefront of stellar astrophysics. Numerical simulations are an essential tool for this task. We present an extension of the existing numerical tool-kit A-MaZe that enables such simulations of stratified flows in a gravitational field. The finite-volume based, cell-centered, and time-explicit hydrodynamics solver of A-MaZe was extended such that the scheme is now well-balanced in both momentum and energy. The algorithm maintains an initially static balance between gravity and pressure to machine precision. Quasi-stationary convection in slab-geometry preserves gas energy (internal plus kinetic) on average despite strong local up- and down-drafts. By contrast, a more standard numerical scheme is demonstrated to result in substantial gains of energy within a short time on purely numerical grounds. The test is further used to point out the role of dimensionality, viscosity, and Rayleigh number for compressible convection. Applications to a young sun in 2D and 3D, covering a part of the inner radiative zone as well as the outer convective zone, demonstrate that the scheme meets its initial design goal. Comparison with results obtained for a physically identical setup with a time-implicit code show qualitative agreement.

We present the first detailed three-dimensional hydrodynamic implicit large eddy simulations of turbulent convection for carbon burning. The simulations start with an initial radial profile mapped from a carbon burning shell within a 15 solar mass stellar evolution model. We considered 4 resolutions from 128^3 to 1024^3 zones. These simulations confirm that convective boundary mixing (CBM) occurs via turbulent entrainment as in the case of oxygen burning. The expansion of the boundary into the surrounding stable region and the entrainment rate are smaller at the bottom boundary because it is stiffer than the upper boundary. The results of this and similar studies call for improved CBM prescriptions in 1D stellar evolution models.

We study lithium depletion in low-mass and solar-like stars as a function of time, using a new diffusion coefficient describing extra-mixing taking place at the bottom of a convective envelope. This new form is motivated by multi-dimensional fully compressible, time implicit hydrodynamic simulations performed with the MUSIC code. Intermittent convective mixing at the convective boundary in a star can be modeled using extreme value theory, a statistical analysis frequently used for finance, meteorology, and environmental science. In this letter, we implement this statistical diffusion coefficient in a one-dimensional stellar evolution code, using parameters calibrated from multi-dimensional hydrodynamic simulations of a young low-mass star. We propose a new scenario that can explain observations of the surface abundance of lithium in the Sun and in clusters covering a wide range of ages, from $\sim$ 50 Myr to $\sim$ 4 Gyr. Because it relies on our physical model of convective penetration, this scenario has a limited number of assumptions. It can explain the observed trend between rotation and depletion, based on a single additional assumption, namely that rotation affects the mixing efficiency at the convective boundary. We suggest the existence of a threshold in stellar rotation rate above which rotation strongly prevents the vertical penetration of plumes and below which rotation has small effects. In addition to providing a possible explanation for the long standing problem of lithium depletion in pre-main sequence and main sequence stars, the strength of our scenario is that its basic assumptions can be tested by future hydrodynamic simulations.

We present the first detailed three-dimensional (3D) hydrodynamic implicit large eddy simulations of turbulent convection of carbon burning in massive stars. Simulations begin with radial profiles mapped from a carbon burning shell within a 15$\,\textrm{M}_\odot$ one-dimensional stellar evolution model. We consider models with $128^3$, $256^3$, $512^3$ and $1024^3$ zones. The turbulent flow properties of these carbon burning simulations are very similar to the oxygen burning case. We performed a mean field analysis of the kinetic energy budgets within the Reynolds-averaged Navier-Stokes framework. For the upper convective boundary region, we find that the numerical dissipation is insensitive to resolution for linear mesh resolutions above 512 grid points. For the stiffer, more stratified lower boundary, our highest resolution model still shows signs of decreasing sub-grid dissipation suggesting it is not yet numerically converged. We find that the widths of the upper and lower boundaries are roughly 30% and 10% of the local pressure scale heights, respectively. The shape of the boundaries is significantly different from those used in stellar evolution models. As in past oxygen-shell burning simulations, we observe entrainment at both boundaries in our carbon-shell burning simulations. In the large P\'eclet number regime found in the advanced phases, the entrainment rate is roughly inversely proportional to the bulk Richardson number, Ri$_{\rm B}$ ($\propto $Ri${\rm_B}^{-\alpha}$, $0.5\lesssim \alpha \lesssim 1.0$). We thus suggest the use of Ri$_{\rm B}$ as a means to take into account the results of 3D hydrodynamics simulations in new 1D prescriptions of convective boundary mixing.

We examine a penetration layer formed between a central radiative zone and a large convection zone in the deep interior of a young low-mass star. Using the Multidimensional Stellar Implicit Code (MUSIC) to simulate two-dimensional compressible stellar convection in a spherical geometry over long times, we produce statistics that characterize the extent and impact of convective penetration in this layer. We apply extreme value theory to the maximal extent of convective penetration at any time. We compare statistical results from simulations which treat non-local convection, throughout a large portion of the stellar radius, with simulations designed to treat local convection in a small region surrounding the penetration layer. For each of these situations, we compare simulations of different resolution, which have different velocity magnitudes. We also compare statistical results between simulations that radiate energy at a constant rate to those that allow energy to radiate from the stellar surface according to the local surface temperature. Based on the frequency and depth of penetrating convective structures, we observe two distinct layers that form between the convection zone and the stable radiative zone. We show that the probability density function of the maximal depth of convective penetration at any time corresponds closely in space with the radial position where internal waves are excited. We find that the maximal penetration depth can be modeled by a Weibull distribution with a small shape parameter. Using these results, and building on established scalings for diffusion enhanced by large-scale convective motions, we propose a new form for the diffusion coefficient that may be used for one-dimensional stellar evolution calculations in the large P\'eclet number regime. These results should contribute to the 321D link.

We present the results of a numerical benchmark study for the MUlti-dimensional Stellar Implicit Code (MUSIC) based on widely applicable two- and three-dimensional compressible hydrodynamics problems relevant to stellar interiors. MUSIC is an implicit large eddy simulation code that uses implicit time integration, implemented as a Jacobian-free Newton Krylov method. A physics based preconditioning technique which can be adjusted to target varying physics is used to improve the performance of the solver. The problems used for this benchmark study include the Rayleigh-Taylor and Kelvin-Helmholtz instabilities, and the decay of the Taylor-Green vortex. Additionally we show a test of hydrostatic equilibrium, in a stellar environment which is dominated by radiative effects. In this setting the flexibility of the preconditioning technique is demonstrated. This work aims to bridge the gap between the hydrodynamic test problems typically used during development of numerical methods and the complex flows of stellar interiors. A series of multi-dimensional tests are performed and analysed. Each of these test cases is analysed with a simple, scalar diagnostic, with the aim of enabling direct code comparisons. As the tests performed do not have analytic solutions we verify MUSIC by comparing to established codes including ATHENA and the PENCIL code. MUSIC is able to both reproduce behaviour from established and widely-used codes as well as results expected from theoretical predictions. This benchmarking study concludes a series of papers describing the development of the MUSIC code and provides confidence in the future applications.

We present the first 3D simulation of the last minutes of oxygen shell burning in an 18 solar mass supernova progenitor up to the onset of core collapse. A moving inner boundary is used to accurately model the contraction of the silicon and iron core according to a 1D stellar evolution model with a self-consistent treatment of core deleptonization and nuclear quasi-equilibrium. The simulation covers the full solid angle to allow the emergence of large-scale convective modes. Due to core contraction and the concomitant acceleration of nuclear burning, the convective Mach number increases to ~0.1 at collapse, and an l=2 mode emerges shortly before the end of the simulation. Aside from a growth of the oxygen shell from 0.51 to 0.56 solar masses due to entrainment from the carbon shell, the convective flow is reasonably well described by mixing length theory, and the dominant scales are compatible with estimates from linear stability analysis. We deduce that artificial changes in the physics, such as accelerated core contraction, can have precarious consequences for the state of convection at collapse. We argue that scaling laws for the convective velocities and eddy sizes furnish good estimates for the state of shell convection at collapse and develop a simple analytic theory for the impact of convective seed perturbations on shock revival in the ensuing supernova. We predict a reduction of the critical luminosity for explosion by 12--24% due to seed asphericities for our 3D progenitor model relative to the case without large seed perturbations.

Context: We study the impact of two-dimensional spherical shells on compressible convection. Realistic profiles for density and temperature from a one-dimensional stellar evolution code are used to produce a model of a large stellar convection zone representative of a young low-mass star. Methods: We perform hydrodynamic implicit large-eddy simulations of compressible convection using the MUltidimensional Stellar Implicit Code (MUSIC). Because MUSIC has been designed to use realistic stellar models produced from one-dimensional stellar evolution calculations, MUSIC simulations are capable of seamlessly modeling a whole star. Simulations in two-dimensional spherical shells that have different radial extents are performed over hundreds of convective turnover times, permitting the collection of well-converged statistics. Results: We evaluate basic statistics of the convective turnover time, the convective velocity, and the overshooting layer. These quantities are selected for their relevance to one-dimensional stellar evolution calculations, so that our results are focused toward the 321D link. The inclusion in the spherical shell of the boundary between the radiative and convection zones decreases the amplitude of convective velocities in the convection zone. The inclusion of near-surface layers in the spherical shell can increase the amplitude of convective velocities, although the radial structure of the velocity profile established by deep convection is unchanged. The impact from including the near-surface layers depends on the speed and structure of small-scale convection in the near-surface layers. Larger convective velocities in the convection zone result in a commensurate increase in the overshooting layer width and decrease in the convective turnover time. These results provide support for non-local aspects of convection.

Recent numerical simulations suggest that the model by Zahn (1992, A&A, 265, 115) for the turbulent mixing of chemical elements due to differential rotation in stellar radiative zones is valid. We investigate the robustness of this result with respect to the numerical configuration and Reynolds number of the flow. We compare results from simulations performed with two different numerical codes, including one that uses the shearing-box formalism. We also extensively study the dependence of the turbulent diffusion coefficient on the turbulent Reynolds number. The two numerical codes used in this study give consistent results. The turbulent diffusion coefficient is independent of the size of the numerical domain if at least three large turbulent structures fit in the box. Generally, the turbulent diffusion coefficient depends on the turbulent Reynolds number. However, our simulations suggest that an asymptotic regime is obtained when the turbulent Reynolds number is larger than $10^3$. Shear mixing in the regime of small P\'eclet numbers can be investigated numerically both with shearing-box simulations and simulations using explicit forcing. Our results suggest that Zahn's model is valid at large turbulent Reynolds numbers.

This work is the first attempt to describe the multi-dimensional structure of accreting young stars based on fully compressible time implicit multi-dimensional hydrodynamics simulations. One major motivation is to analyse the validity of accretion treatment used in previous 1D stellar evolution studies. We analyse the effect of accretion on the structure of a realistic stellar model of the young Sun. Our work is inspired by the numerical work of Kley \& Lin (1996, ApJ, 461, 933) devoted to the structure of the boundary layer in accretion disks. We analyse the redistribution of accreted material with a range of values of specific entropy relative to the bulk specific entropy of the material in the accreting object's convective envelope. A primary goal is to understand whether and how accreted energy deposited onto a stellar surface is redistributed in the interior. This study focusses on the high accretion rates characteristic of FU Ori systems. We find that the highest entropy cases produce a distinctive behaviour in the mass redistribution, rms velocities, and enthalpy flux in the convective envelope. This change in behaviour is characterised by the formation of a hot layer on the surface of the accreting object, which tends to suppress convection in the envelope. We analyse the long-term effect of such a hot buffer zone on the structure and evolution of the accreting object with 1D stellar evolution calculations. We study the relevance of the assumption of redistribution of accreted energy into the stellar interior used in the literature. One conclusion is that, for a given amount of accreted energy transferred to the accreting object, a treatment assuming accretion energy redistribution throughout the stellar interior could significantly overestimate the effects on the stellar structure, in particular, on the resulting expansion.

Recently Pasetto et al. have proposed a new method to derive a convection theory appropriate for the implementation in stellar evolution codes. Their approach is based on the simple physical picture of spherical bubbles moving within a potential flow in dynamically unstable regions, and a detailed computation of the bubble dynamics. Based on this approach the authors derive a new theory of convection which is claimed to be parameter free, non-local and time-dependent. This is a very strong claim, as such a theory is the holy grail of stellar physics. Unfortunately we have identified several distinct problems in the derivation which ultimately render their theory inapplicable to any physical regime. In addition we show that the framework of spherical bubbles in potential flows is unable to capture the essence of stellar convection, even when equations are derived correctly.

Stellar evolution models of massive stars are important for many areas of astrophysics, for example nucleosynthesis yields, supernova progenitor models and understanding physics under extreme conditions. Turbulence occurs in stars primarily due to nuclear burning at different mass coordinates within the star. The understanding and correct treatment of turbulence and turbulent mixing at convective boundaries in stellar models has been studied for decades but still lacks a definitive solution. This paper presents initial results of a study on convective boundary mixing (CBM) in massive stars. The 'stiffness' of a convective boundary can be quantified using the bulk Richardson number ($\textrm{Ri}_B$), the ratio of the potential energy for restoration of the boundary to the kinetic energy of turbulent eddies. A 'stiff' boundary ($\textrm{Ri}_B \sim 10^4$) will suppress CBM, whereas in the opposite case a 'soft' boundary ($\textrm{Ri}_B \sim 10$) will be more susceptible to CBM. One of the key results obtained so far is that lower convective boundaries (closer to the centre) of nuclear burning shells are 'stiffer' than the corresponding upper boundaries, implying limited CBM at lower shell boundaries. This is in agreement with 3D hydrodynamic simulations carried out by Meakin and Arnett [The Astrophysical Journal 667:448-475, 2007]. This result also has implications for new CBM prescriptions in massive stars as well as for nuclear burning flame front propagation in Super-Asymptotic Giant Branch stars and also the onset of novae.

This work is a continuation of our efforts to develop an efficient implicit solver for multidimensional hydrodynamics for the purpose of studying important physical processes in stellar interiors, such as turbulent convection and overshooting. We present an implicit solver that results from the combination of a Jacobian-Free Newton-Krylov method and a preconditioning technique tailored to the inviscid, compressible equations of stellar hydrodynamics. We assess the accuracy and performance of the solver for both 2D and 3D problems for Mach numbers down to $10^{-6}$. Although our applications concern flows in stellar interiors, the method can be applied to general advection and/or diffusion-dominated flows. The method presented in this paper opens up new avenues in 3D modeling of realistic stellar interiors allowing the study of important problems in stellar structure and evolution.

Several parametrizations for overshooting in 1D stellar evolution calculations coexist in the literature. These parametrizations are used somewhat arbitrarily in stellar evolution codes, based on what works best for a given problem, or even for historical reasons related to the development of each code. We bring attention to the fact that these different parametrizations correspond to different physical regimes of overshooting, depending whether the effects of radiation are dominant, marginal, or negligible. Our analysis is based on previously published theoretical results, as well as multidimensional hydrodynamical simulations of stellar convection where the interaction between the convective region and a stably-stratified region is observed. Although the underlying hydrodynamical processes are the same, the outcome of the overshooting process is profoundly affected by radiative effects. Using a simple picture of the scales involved in the overshooting process, we show how three regimes are obtained, depending on the importance of radiative effects. These three regimes correspond to the different behaviors observed in hydrodynamical simulations so far, and to the three types of parametrizations used in 1D codes. We suggest that the existing parametrizations for overshooting should coexist in 1D stellar evolution codes, and should be applied consistently at convective boundaries depending on the local physical conditions.

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.

In this contribution we present initial results of a study on convective boundary mixing (CBM) in massive stellar models using the GENEVA stellar evolution code. Before undertaking costly 3D hydrodynamic simulations, it is important to study the general properties of convective boundaries, such as the: composition jump; pressure gradient; and `stiffness'. Models for a 15Mo star were computed. We found that for convective shells above the core, the lower (in radius or mass) boundaries are `stiffer' according to the bulk Richardson number than the relative upper (Schwarzschild) boundaries. Thus, we expect reduced CBM at the lower boundaries in comparison to the upper. This has implications on flame front propagation and the onset of novae.

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.

We present a statistical analysis of turbulent convection in stars within our Reynolds-Averaged Navier Stokes (RANS) framework in spherical geometry which we derived from first principles. The primary results reported in this document include: (1) an extensive set of mean-field equations for compressible, multi-species hydrodynamics, and (2) corresponding mean-field data computed from various simulation models. Some supplementary scale analysis data is also presented. The simulation data which is presented includes: (1) shell convection during oxygen burning in a 23 solar mass supernova progenitor, (2) envelope convection in a 5 solar mass red giant, (3) shell convection during the helium flash, and (4) a hydrogen injection flash in a 1.25 solar mass star. These simulations have been partially described previously in Meakin [2006], Meakin and Arnett [2007a,b, 2010], Arnett et al. [2009, 2010], Viallet et al. [2011, 2013a,b] and Mocak et al. [2009, 2011]. New data is also included in this document with several new domain and resolution configurations as well as some variations in the physical model such as convection zone depth and driving source term. The long term goal of this work is to aid in the development of more sophisticated models for treating hydrodynamic phenomena (e.g., turbulent convection) in the field of stellar evolution by providing a direct link between 3D simulation data and the mean fields which are modeled by 1D stellar evolution codes. As such, this data can be used to test previously proposed turbulence models found in the literature and sometimes used in stellar modeling. This data can also serve to test basic physical principles for model building and inspire new prescriptions for use in 1D evolution codes.

Time-implicit schemes are attractive since they allow numerical time steps that are much larger than those permitted by the Courant-Friedrich-Lewy criterion characterizing time-explicit methods. This advantage comes, however, with a cost: the solution of a system of nonlinear equations is required at each time step. In this work, the nonlinear system results from the discretization of the hydrodynamical equations with the Crank-Nicholson scheme. We compare the cost of different methods, based on Newton-Raphson iterations, to solve this nonlinear system, and benchmark their performances against time-explicit schemes. Since our general scientific objective is to model stellar interiors, we use as test cases two realistic models for the convective envelope of a red giant and a young Sun. Focusing on 2D simulations, we show that the best performances are obtained with the quasi-Newton method proposed by Broyden. Another important concern is the accuracy of implicit calculations. Based on the study of an idealized problem, namely the advection of a single vortex by a uniform flow, we show that there are two aspects: i) the nonlinear solver has to be accurate enough to resolve the truncation error of the numerical discretization, and ii) the time step has be small enough to resolve the advection of eddies. We show that with these two conditions fulfilled, our implicit methods exhibit similar accuracy to time-explicit schemes, which have lower values for the time step and higher computational costs. Finally, we discuss in the conclusion the applicability of these methods to fully implicit 3D calculations.

We present 3D implicit large eddy simulations (ILES) of the turbulent convection in the envelope of a 5 Msun red giant star and in the oxygen-burning shell of a 23 Msun supernova progenitor. The numerical models are analyzed in the framework of 1D Reynolds-Averaged Navier-Stokes (RANS) equations. The effects of pressure fluctuations are more important in the red giant model, owing to larger stratification of the convective zone. We show how this impacts different terms in the mean-field equations. We clarify the driving sources of kinetic energy, and show that the rate of turbulent dissipation is comparable to the convective luminosity. Although our flows have low Mach number and are nearly adiabatic, our analysis is general and can be applied to photospheric convection as well. The robustness of our analysis of turbulent convection is supported by the insensitivity of the mean-field balances to linear mesh resolution. We find robust results for the turbulent convection zone and the stable layers in the oxygen-burning shell model, and robust results everywhere in the red giant model, but the mean fields are not well converged in the narrow boundary regions (which contain steep gradients) in the oxygen-burning shell model. This last result illustrates the importance of unresolved physics at the convective boundary, which governs the mixing there.

Studies of beryllium abundance in large samples of solar-type stars show a small fraction of extremely beryllium-deficient stars, which challenges our current understanding of light element depletion in these stars. We suggest two possible scenarios that may explain this high level of Be depletion: early accretion and rotational mixing. We show that in both cases, the conditions required to reach the observed level of Be depletion are quite extreme, which explains the very small fraction of detected Be outliers. We suggest that substantial Be depletion can be obtained in stars if they were fast rotators in the past, with high initial rotational velocities and short disc lifetimes. Our analysis suggests that rotational mixing may not be efficient enough to deplete Be in less than 10 Myr. Consequently, the detection of strongly Be-deficient stars in clusters younger than $\sim$ 10 Myr may provide a genuine signature of accretion process and the proof that some protostars may undergo many extreme bursts of accretion during their embedded phases of evolution.

This paper describes the first steps of development of a new multidimensional time implicit code devoted to the study of hydrodynamical processes in stellar interiors. The code solves the hydrodynamical equations in spherical geometry and is based on the finite volume method. Radiation transport is taken into account within the diffusion approximation. Realistic equation of state and opacities are implemented, allowing the study of a wide range of problems characteristic of stellar interiors. We describe in details the numerical method and various standard tests performed to validate the method. We present preliminary results devoted to the description of stellar convection. We first perform a local simulation of convection in the surface layers of a A-type star model. This simulation is used to test the ability of the code to address stellar conditions and to validate our results, since they can be compared to similar previous simulations based on explicit codes. We then present a global simulation of turbulent convective motions in a cold giant envelope, covering 80% in radius of the stellar structure. Although our implicit scheme is unconditionally stable, we show that in practice there is a limitation on the time step which prevent the flow to move over several cells during a time step. Nevertheless, in the cold giant model we reach a hydro CFL number of 100. We also show that we are able to address flows with a wide range of Mach numbers (10^-3 < Ms< 0.5), which is impossible with an anelastic approach. Our first developments are meant to demonstrate that the use of an implicit scheme applied to a stellar evolution context is perfectly thinkable and to provide useful guidelines to optimise the development of an implicit multi-D hydrodynamical code.

We report on the development of an implicit multi-D hydrodynamic code for stellar evolution. We present two test-cases relevant for the first scientific goal of the code: the simulation of convection in pulsating stars. First results on a realistic stellar model are also presented.

Accretion disks in AGN should be subject to the same type of instability as in cataclysmic variables (CVs) or in low-mass X-ray binaries (LMXBs), which leads to dwarf nova and soft X-ray transient outbursts. It has been suggested that this thermal/viscous instability can account for the long term variability of AGNs. We test this assertion by presenting a systematic study of the application of the disk instability model (DIM) to AGNs. We are using the adaptative grid numerical code we have developed in the context of CVs, enabling us to fully resolve the radial structure of the disk. We show that, because in AGN disks the Mach numbers are very large, the heating and cooling fronts are so narrow that they cannot be resolved by the numerical codes that have been used until now. In addition, these fronts propagate on time scales much shorter than the viscous time. As a result, a sequence of heating and cooling fronts propagate back and forth in the disk, leading only to small variations of the accretion rate onto the black hole, with short quiescent states occurring for very low mass transfer rates only. Truncation of the inner part of the disk by e.g. an ADAF does not alter this result, but enables longer quiescent states. Finally we discuss the effects of irradiation by the central X-ray source, and show that, even for extremely high irradiation efficiencies, outbursts are not a natural outcome of the model.

We discuss two mechanisms that could result in an enhancement of the mass transfer rate during outbursts of dwarf novae and soft X-ray transients: the hot outer disc rim itself could heat the L1 point and scattered radiation by optically thin outflowing matter could also heat L1 significantly. We determine quantitatively the increase of the mass transfer rate resulting from an extra heating. During outbursts, the disc edge heats up the upper layer of the secondary with a flux of the order of the intrinsic stellar flux. This probably has no large effect on the mass transfer rate. In soft X-ray transients, the environing medium of the disc (corona+wind) could back-scatter a certain fraction of the accretion luminosity toward L1. Since soft X-ray transients reach high luminosities, even a low efficiency of this effect could yield a significant heating of L1, whereas we show that in dwarf novae this effect is negligible. Initially the incoming radiation does not penetrate below the photosphere of the secondary. Depending on the heating efficiency, which has to be determined, the mass transfer rate could be significantly increased (~x100).

We investigate numerically the surface flow on the secondary star during outbursts. We use a simple model for the irradiation and the geometry of the secondary star: the irradiation temperature is treated as a free parameter and the secondary is replaced by a spherical star with a space-dependent Coriolis force that mimics the effect of the Roche geometry. The Euler equations are solved in spherical coordinates with the TVD-MacCormack scheme. We show that the Coriolis force leads to the formation of a circulation flow from high latitude region to the close vicinity of the $L_1$ point. However no heat can be efficiently transported to the $L_1$ region due to the rapid radiative cooling of the hot material as it enters the equatorial belt shadowed from irradiation. Under the assumption of hydrostatic equilibrium, the Coriolis force could lead to a moderate increase of the mass transfer rate by pushing the gas in the vertical direction at the $L_1$ point, but only during the initial phases of the outburst (about 15 -- 20 orbital periods). We conclude that the Coriolis force does not prevent a flow from the heated regions of the secondary towards the $L_1$ region, at least during the initial phase of an outburst, but the resulting increase of the mass transfer rate is moderate, and it is unlikely to be able to account for the duration of long outbursts.

The long term magnetohydrodynamic stability of magnetized transonic/supersonic jets is numerically investigated using a spatial approach. We focus on two-dimensional linearly-unstable slab configurations where the jet is embedded in a flow-aligned uniform magnetic field of weak amplitude. We compare our results with previous studies using a temporal approach where longitudinally periodic domains were adopted. The finite-volume based versatile advection code (VAC) is used to solve the full set of ideal compressible MHD equations. We follow the development of Kelvin-Helmholtz modes that are driven by a white noise perturbation continuously introduced at the jet inlet. No noticeable difference is observed in spatial simulations versus analogous temporal ones during the linear and early non-linear evolution of the configuration. However, in the case of transonic flows, a different long-term scenario occurs in our spatial runs. Indeed, after the large-scale disruption of the flow, a sheath region of enhanced magnetic field encompassing the jet core forms along the whole flow. This provides a partial stabilization mechanism leading to enhanced stability for later times, which is almost independent of the initial magnitude of the magnetic field. The implication of this mechanism for the stability of astrophysical jets is discussed.