As in GD2011, our system represents a zoom around an ionisation region. Since we are computing local simulations, the vertical gravity g = −gez and the kinematic viscosity ν are assumed to be constant. Following our purely radiative model of the κ-mechanism (Gastine & Dintrans 2008a,b), the ionisation region is represented by a temperature-dependent radiative conductivity profile that mimics an opacity bump:
(A.2)where Tbump is the position of the hollow in temperature and σ, e, and denote its slope, width, and relative amplitude, respectively. We assume both radiative and hydrostatic equilibria; that is,
where Fbot is the imposed bottom flux. Following GD2011, we chose Lz as the length scale, i.e. [x] = Lz, top density ρtop and top temperature Ttop as density and temperature scales, respectively. The velocity scale is then , while time is given in units of .
Table A.1 then summarises the parameters of the numerical simulations presented in this study in these dimensionless units. The penultimate column of this table contains the value of the frequency ω00 of the fundamental unstable radial mode excited by the κ-mechanism, which lies between 3 and 4 for every DNS. The last column gives the value of the Rayleigh number, which quantifies the strength of the convective motions. It is given by
where Lconv is the width of the convective zone, χ = K0/ρ0cp the radiative diffusivity, and s the entropy.
Dimensionless parameters of the numerical simulations.
With the parallel version of the alternate direction implicit (ADI) solver for the radiative diffusion implemented in the pencil code (see GD2011), we advance the following hydrodynamic equations in time:
where ρ, u, and T denote density, velocity, and temperature, respectively, while K(T) is given by Eq. (A.1). The operator D/Dt = ∂/∂t + u·∇ is the usual total derivative, while S is the (traceless) rate-of-strain tensor given by
Finally, we impose the condition that all fields are periodic in the horizontal direction, while stress-free boundary conditions (i.e., uz = 0 and dux/dz = 0) are assumed for the velocity in the vertical one. Concerning the temperature, a perfect conductor at the bottom (i.e., flux imposed) and a perfect insulator at the top (i.e., temperature imposed) are applied.
To ensure that both the nonlinear saturation and thermal relaxation are achieved, the simulations were computed over very long times, typically t ≳ 3000. As the eigenfrequency of the unstable acoustic mode ω00 ∈ [3 − 4] (see Table A.1), this corresponds approximately to 1800 periods of oscillation.
© ESO, 2011