We present a numerical study of double-diffusive convection characterized by a stratification unstable to thermal convection while at the same time a mean molecular weight (or solute concentration) difference between top and bottom counteracts this instability. Convective zones can form in this case either by the stratification being locally unstable to the combined action of both temperature and solute gradients or by another process, the oscillatory double-diffusive convective instability, which is triggered by the faster molecular diffusivity of heat in comparison with that one of the solute. We discuss successive layer formation for this problem in the case of an instantaneously heated bottom (plate) which forms a first layer with an interface that becomes temporarily unstable and triggers the formation of further, secondary layers. We consider both the case of a Prandtl number typical for water (oceanographic scenario) and of a low Prandtl number (giant planet scenario). We discuss the impact of a Couette like shear on the flow and in particular on layer formation for different shear rates. Additional layers form due to the oscillatory double-diffusive convective instability, as is observed for some cases. We also test the physical model underlying our numerical experiments by recovering experimental results of layer formation obtained in laboratory setups.

We present results of a large, high resolution 3D hydrodynamical simulation of the surface layers of a DA white dwarf (WD) with $T_{\rm eff}=11800$ K and $\log(g)=8$ using the ANTARES code, the widest and deepest such simulation to date. Our simulations are in good agreement with previous calculations in the Schwarzschild-unstable region and in the overshooting region immediately beneath it. Farther below, in the wave-dominated region, we find that the rms horizontal velocities decay with depth more rapidly than the vertical ones. Since mixing requires both vertical and horizontal displacements, this could have consequences for the size of the region that is well mixed by convection, if this trend is found to hold for deeper layers. We discuss how the size of the mixed region affects the calculated settling times and inferred steady-state accretion rates for WDs with metals observed in their atmospheres.

So far, numerical studies of double-diffusive layering in turbulent convective flows have neglected the effects of rotation. We undertake a first step into that direction by investigating how Coriolis forces affect a double-diffusive layer inside a rotating spherical shell. For this purpose we have run simulations in a parameter regime where these layers are expected to form and successively increased the rate of rotation with the result that fast rotation is found to have a similar stabilising effect on the overall convective flux as an increase of the stability ratio $R_{\rho}$ has in a non-rotating setup. We have also studied to what extent the regimes of rotational constraints suggested by King, Stellmach, and Buffett (2013) for rotation in the case of Rayleigh-B\'enard convection are influenced by double-diffusive convection: their classification could also be applicable to the case of double-diffusive convection in a spherical shell if it is extended to be also a function of the stability ratio $R_{\rho}$. Furthermore, we examined the ratio of saline and thermal Nusselt numbers and compared our results with models of Spruit (2013), Rosenblum et al. (2011) and Wood, Garaud, and Stellmach (2013). We find our data to be fitted best by Spruit's model. Our result that fast rotation further decreases the convective transport, which is already lowered by double-diffusive convection, could play a major role for e.g. the modeling of the interior of some rapidly rotating giant planets, as gaseous giant planets have recently been proposed to be influenced by double-diffusive convection.

We have implemented open boundary conditions into the ANTARES code to increase the realism of our simulations of stellar surface convection. Even though we greatly benefit from the high accuracy of our fifth order numerical scheme (WENO5), the broader stencils needed for the numerical scheme complicate the implementation of boundary conditions. We show that the effective temperature of a numerical simulation cannot be changed by corrections at the lower boundary since the thermal stratification does only change on the Kelvin-Helmholtz time scale. Except for very shallow models, this time scale cannot be covered by multidimensional simulations due to the enormous computational requirements. We demonstrate to what extent numerical simulations of stellar surface convection are sensitive to the initial conditions and the boundary conditions. An ill-conceived choice of parameters for the boundary conditions can have a severe impact. Numerical simulations of stellar surface convection will only be (physically) meaningful and realistic if the initial model, the extent and position of the simulation box, and the parameters from the boundary conditions are chosen adequately.

A grid of numerical simulations of double-diffusive convection is presented for the astrophysical case where viscosity (Prandtl number Pr) and solute diffusivity (Lewis number Le) are much smaller than the thermal diffusivity. As in laboratory and geophysical cases convection takes place in a layered form. The proper translation between subsonic flows in a stellar interior and an incompressible (Boussinesq) fluid is given, and the validity of the Boussinesq approximation for the semiconvection problem is checked by comparison with fully compressible simulations. The predictions of a simplified theory of mixing in semiconvection given in a companion paper are tested against the numerical results, and used to extrapolate these to astrophysical conditions. The predicted effective He-diffusion coefficient is nearly independent of the double-diffusive layering thickness $d$. For a fiducial main sequence model (15 $M_\odot$) the inferred mixing time scale is of the order $10^{10}$ yr. An estimate for the secular increase of $d$ during the semiconvective phase is given. It can potentially reach a significant fraction of a pressure scale height.

A grid of numerical simulations of double-diffusive convection is presented for astrophysical conditions. As in laboratory and geophysical cases convection takes place in a layered form. A translation between the astrophysical fluid mechanics and incompressible (Boussinesq) approximation is given, valid for thin layers. Its validity is checked by comparison of the results of fully compressible and Boussinesq simulations of semiconvection. A fitting formula is given for the superadiabatic gradient as a function of this parameter. The superadiabaticity depends on the thickness $d$ of the double diffusive layers, for which no good theory is available, but the effective He-diffusion coefficient is nearly independent of $d$. For a fiducial main sequence model (15 $M_\odot$) the inferred mixing time scale is of the order $10^{10}$ yr.

We report on modelling in stellar astrophysics with the ANTARES code. First, we describe properties of turbulence in solar granulation as seen in high-resolution calculations. Then, we turn to the first 2D model of pulsation-convection interaction in a cepheid. We discuss properties of the outer and the HEII ionization zone. Thirdly, we report on our work regarding models of semiconvection in the context of stellar physics.