The ANTARES radiation hydrodynamics code is capable of simulating the solar granulation in detail unequaled by direct observation. We introduce a state-of-the-art numerical tool to the solar physics community and demonstrate its applicability to model the solar granulation. The code is based on the weighted essentially non-oscillatory finite volume method and by its implementation of local mesh refinement is also capable of simulating turbulent fluids. While the ANTARES code already provides promising insights into small-scale dynamical processes occurring in the quiet-Sun photosphere, it will soon be capable of modeling the latter in the scope of radiation magnetohydrodynamics. In this first preliminary study we focus on the vertical photospheric stratification by examining a 3-D model photosphere with an evolution time much larger than the dynamical timescales of the solar granulation and of particular large horizontal extent corresponding to $25\!" \!\! \times \, 25\!"$ on the solar surface to smooth out horizontal spatial inhomogeneities separately for up- and downflows. The highly resolved Cartesian grid thereby covers $\sim 4~\mathrm{Mm}$ of the upper convection zone and the adjacent photosphere. Correlation analysis, both local and two-point, provides a suitable means to probe the photospheric structure and thereby to identify several layers of characteristic dynamics: The thermal convection zone is found to reach some ten kilometers above the solar surface, while convectively overshooting gas penetrates even higher into the low photosphere. An $\approx 145\,\mathrm{km}$ wide transition layer separates the convective from the oscillatory layers in the higher photosphere.

Previous studies have discovered a population of small granules with diameters less than 800 km showing differing physical properties. High resolution simulations and observations of the solar granulation, in combination with automated segmentation and tracking algorithms, allow us to study the evolution of the structural and physical properties of these granules and surrounding vortex motions with high temporal and spatial accuracy. We focus on the dynamics of granules (lifetime, fragmentation, size, position, intensity, vertical velocity) over time and the influence of strong vortex motions. Of special interest are the dynamics of small granules compared to regular-sized granules. We developed a temporal tracking algorithm based on our developed segmentation algorithm for solar granulation. This was applied to radiation hydrodynamics simulations and high resolution observations of the quiet Sun by SUNRISE/IMaX. The dynamics of small granules differ in regard to their diameter, intensity and depth evolution compared to regular granules. The tracked granules in the simulation and observations reveal similar dynamics (lifetime, evolution of size, vertical velocity and intensity). The fragmentation analysis shows that the majority of granules in simulations do not fragment, while the opposite was found in observations. Strong vortex motions were detected at the location of small granules. Regions of strong vertical vorticity show high intensities and downflow velocities, and live up to several minutes. The analysis of granules separated according to their diameter in different groups reveals strongly differing behaviors. The largest discrepancies can be found within the groups of small, medium-sized and large granules and have to be analyzed independently. The predominant location of vortex motions on and close to small granules indicates a strong influence on the dynamics of granules.

After an overview of general aspects of modelling the pulsation- convection interaction we present reasons why such simulations (in multidimensions) are needed but, at the same time, pose a considerable challenge. We then discuss, for several topics, what insights multidimensional simulations have either already provided or can be expected to yield in the future. We finally discuss properties of our ANTARES code. Many of these features can be expected to be characteristic of other codes which may possibly be applied to these physical questions in the foreseeable future.

Non-local, time-dependent convection models have been used to explain the location of double-mode pulsations in Cepheids in the HR diagram as well as the existence and location of the red edge of the instability strip. These properties are highly sensitive to model parameters. We use 2D radiation hydrodynamical simulations with realistic microphysics and grey radiative-transfer to model a short period Cepheid. The simulations show that the strength of the convection zone varies significantly over the pulsation period and exhibits a phase shift relative to the variations in radius. We evaluate the convective flux and the work integral as predicted by the most common convection models. It turns out that over one pulsation cycle the model parameter $\alpha_{\rm c}$, has to be varied by up to a factor of beyond 2 to match the convective flux obtained from the simulations. To bring convective fluxes integrated over the He II convection zone and the overshoot zone below into agreement, this parameter has to be varied by a factor of up to $\sim 7.5$ (Kuhfu{\ss}). We then present results on the energetics of the convection and overshoot zone by radially symmetric and fluctuating quantities. To successfully model this scenario by a static, one dimensional or even by a simple time-dependent model appears extremely challenging. We conclude that significant improvements are needed to make predictions based on 1D models more robust and to improve the reliability of conclusions on the convection-pulsation coupling drawn from them. Multidimensional simulations can provide guidelines for developing descriptions of convection then applied in traditional 1D modelling.

We investigate the achievable efficiency of both the time and the space discretisation methods used in Antares for mixed parabolic-hyperbolic problems. We show that the fifth order variant of WENO combined with a second order Runge-Kutta scheme is not only more accurate than standard first and second order schemes, but also more efficient taking the computation time into account. Then, we calculate the error decay rates of WENO with several explicit Runge-Kutta schemes for advective and diffusive problems with smooth and non-smooth initial conditions. With this data, we estimate the computational costs of three-dimensional simulations of stellar surface convection and show that SSP RK(3,2) is the most efficient scheme considered in this comparison.

We investigate the applicability of curvilinear grids in the context of astrophysical simulations and WENO schemes. With the non-smooth mapping functions from Calhoun et al. (2008), we can tackle many astrophysical problems which were out of scope with the standard grids in numerical astrophysics. We describe the difficulties occurring when implementing curvilinear coordinates into our WENO code, and how we overcome them. We illustrate the theoretical results with numerical data. The WENO finite difference scheme works only for high Mach number flows and smooth mapping functions whereas the finite volume scheme gives accurate results even for low Mach number flows and on non-smooth grids.

We have extended the ANTARES code to simulate the coupling of pulsation with convection in Cepheid-like variables in an increasingly realistic way, in particular in multidimensions, 2D at this stage. Present days models of radially pulsating stars assume radial symmetry and have the pulsation-convection interaction included via model equations containing ad hoc closures and moreover parameters whose values are barely known. We intend to construct ever more realistic multidimensional models of Cepheids. In the present paper, the first of a series, we describe the basic numerical approach and how it is motivated by physical properties of these objects which are sometimes more, sometimes less obvious. - For the construction of appropriate models a polar grid co-moving with the mean radial velocity has been introduced to optimize radial resolution throughout the different pulsation phases. The grid is radially stretched to account for the change of spatial scales due to vertical stratification and a new grid refinement scheme is introduced to resolve the upper, hydrogen ionisation zone where the gradient of temperature is steepest. We demonstrate that the simulations are not conservative when the original weighted essentially non-oscillatory method implemented in ANTARES is used and derive a new scheme which allows a conservative time evolution. The numerical approximation of diffusion follows the same principles. Moreover, the radiative transfer solver has been modified to improve the efficiency of calculations on parallel computers. We show that with these improvements the ANTARES code can be used for realistic simulations of the convection-pulsation interaction in Cepheids. We discuss the properties of several models which include the upper 42% of a Cepheid along its radial coordinate, assume different opening angles, and are suitable for an in-depth study of convection and pulsation.

The ANTARES code has been designed for simulation of astrophysical flows in a variety of situations, in particular in the context of stellar physics. Here, we describe extensions as necessary to model the interaction of pulsation and convection in classical pulsating stars. These extensions encomprise the introduction of a spherical grid, movable in the radial direction, specific forms of grid-refinement and considerations regarding radiative transfer. We then present the basic parameters of the cepheid we study more closely. For that star we provide a short discussion of patterns of the H+HeI and the HeII convection zones and the interaction with pulsation seen in the pdV work or atmospheric structures.

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.

Context. Magnetic Bright Points (MBPs) are small-scale magnetic features in the solar photosphere. They may be a possible source of coronal heating by rapid footpoint motions that cause magnetohydrodynamical waves. The number and size distribution are of vital importance in estimating the small scale-magnetic-field energy. Aims. The size distribution of MBPs is derived for G-band images acquired by the Hinode/SOT instrument. Methods. For identification purposes, a new automated segmentation and identification algorithm was developed. Results. For a sampling of 0.108 arcsec/pixel, we derived a mean diameter of (218 +- 48) km for the MBPs. For the full resolved data set with a sampling of 0.054 arcsec/pixel, the size distribution shifted to a mean diameter of (166 +- 31) km. The determined diameters are consistent with earlier published values. The shift is most probably due to the different spatial sampling. Conclusions. We conclude that the smallest magnetic elements in the solar photosphere cannot yet be resolved by G-band observations. The influence of discretisation effects (sampling) has also not yet been investigated sufficiently.

Small-scale magnetic fields in the solar photosphere can be identified in high-resolution magnetograms or in the G-band as magnetic bright points (MBPs). Rapid motions of these fields can cause magneto-hydrodynamical waves and can also lead to nanoflares by magnetic field braiding and twisting. The MBP velocity distribution is a crucial parameter for estimating the amplitudes of those waves and the amount of energy they can contribute to coronal heating. The velocity and lifetime distributions of MBPs are derived from solar G-band images of a quiet sun region acquired by the Hinode/SOT instrument with different temporal and spatial sampling rates. We developed an automatic segmentation, identification and tracking algorithm to analyse G-Band image sequences to obtain the lifetime and velocity distributions of MBPs. The influence of temporal/spatial sampling rates on these distributions is studied and used to correct the obtained lifetimes and velocity distributions for these digitalisation effects. After the correction of algorithm effects, we obtained a mean MBP lifetime of (2.50 +- 0.05) min and mean MBP velocities, depending on smoothing processes, in the range of (1 - 2) km/s. Corrected for temporal sampling effects, we obtained for the effective velocity distribution a Rayleigh function with a coefficient of (1.62 +- 0.05) km/s. The x- and y- components of the velocity distributions are Gaussians. The lifetime distribution can be fitted by an exponential function.

We discuss the general design of the ANTARES code which is intended for simulations in stellar hydrodynamics with radiative transfer and realistic microphysics in 1D, 2D and 3D. We then compare the quality of various numerical methods. We have applied ANTARES in order to obtain high resolution simulations of solar granulation which we describe and analyze. In order to obtain high resolution, we apply grid refinement to a region predominantly occupied by an exploding granule. Strong, rapidly rotating vortex tubes of small diameter (~100 km) generated by the downdrafts and ascending into the photosphere near the granule boundaries evolve, often entering the photosphere from below in an arclike fashion. They essentially contribute to the turbulent velocity field near the granule boundaries.

We present results from 2D radiation-hydrodynamical simulations of fully compressible convection for the surface layers of A-type stars with the ANTARES code. Spectroscopic indicators for photospheric convective velocity fields show a maximum of velocities near Teff ~8000 K. In that range the largest values are measured for the subgroup of Am stars. Thus far, no prognostic model, neither theoretical nor numerical, is able to exactly reproduce the line profiles of sharp line A and Am stars in that temperature range. In general, the helium abundance of A stars is not known from observations. Hence, we have considered two extreme cases for our simulations: a solar helium abundance as an upper limit and zero helium abundance as a lower limit. The simulation for the helium free case is found to differ from the case with solar helium abundance by larger velocities, larger flow structures, and by a sign reversal of the flux of kinetic energy inside the hydrogen ionisation zone. Both simulations show extended shock fronts emerging from the optical surface, as well as mixing far below the region of partial ionisation of hydrogen, and vertical oscillations emerging after initial perturbations have been damped. We discuss problems related to the rapid radiative cooling at the surface of A-type stars such as resolution and efficient relaxation. The present work is considered as a step towards a systematic study of convection in A- to F-type stars, encouraged by the new data becoming available for these objects from both asteroseismological missions and from high resolution spectroscopy.

Using grid refinement, we have simulated solar granulation in 2D. The refined region measures 1.97*2.58 Mm (vertical*horizontal). Grid spacing there is 1.82*2.84 km. The downflows exhibit strong Kelvin-Helmholtz instabilities. Below the photosphere, acoustic pulses are generated. They proceed laterally (in some cases distances of at least the size of our refined domain) and may be enhanced when transversing downflows) as well as upwards where, in the photosphere they contribute significantly to 'turbulence' (velocity gradients, etc.) The acoustic pulses are ubiquitous in that at any time several of them are seen in our high-resolution domain. Their possible contributions to p-mode excitation or heating of the chromosphere needs to be investigated.