Space-borne missions CoRoT and Kepler have provided a rich harvest of high-quality photometric data for solar-like pulsators. It is now possible to measure damping rates for hundreds of main-sequence and thousands of red-giant. However, among the seismic parameters, mode damping rates remain poorly understood and thus barely used for inferring the physical properties of stars. Previous approaches to model mode damping rates were based on mixing-length theory or a Reynolds-stress approach to model turbulent convection. While able to grasp the main physics of the problem, those approaches are of little help to provide quantitative estimates as well as a definitive answer on the relative contribution of each physical mechanism. Our aim is thus to assess the ability of 3D hydrodynamical simulations to infer the physical mechanisms responsible for damping of solar-like oscillations. To this end, a solar high-spatial resolution and long-duration hydrodynamical 3D simulation computed with the ANTARES code allows probing the coupling between turbulent convection and the normal modes of the simulated box. Indeed, normal modes of the simulation experience realistic driving and damping in the super-adiabatic layers of the simulation. Therefore, investigating the properties of the normal modes in the simulation provides a unique insight into the mode physics. We demonstrate that such an approach provides constraints on the solar damping rates and is able to disentangle the relative contribution related to the perturbation of the turbulent pressure, the gas pressure, the radiative flux, and the convective flux contributions. Finally, we conclude that using the normal modes of a 3D numerical simulation is possible and is potentially able to unveil the respective role of the different physical mechanisms responsible for mode damping provided the time-duration of the simulation is long enough.

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.

We present a new fourth-order, finite-volume hydrodynamics code named Apsara. The code employs a high-order, finite-volume method for mapped coordinates with extensions for nonlinear hyperbolic conservation laws. Apsara can handle arbitrary structured curvilinear meshes in three spatial dimensions. The code has successfully passed several hydrodynamic test problems, including the advection of a Gaussian density profile and a nonlinear vortex and the propagation of linear acoustic waves. For these test problems, Apsara produces fourth-order accurate results in case of smooth grid mappings. The order of accuracy is reduced to first-order when using the nonsmooth circular grid mapping. When applying the high-order method to simulations of low-Mach number flows, for example, the Gresho vortex and the Taylor-Green vortex, we discover that Apsara delivers superior results to codes based on the dimensionally split, piecewise parabolic method (PPM) widely used in astrophysics. Hence, Apsara is a suitable tool for simulating highly subsonic flows in astrophysics. In the first astrophysical application, we perform implicit large eddy simulations (ILES) of anisotropic turbulence in the context of core collapse supernova (CCSN) and obtain results similar to those previously reported.

Recent results from high-resolution solar granulation observations indicate the existence of a population of small granular cells that are smaller than 600 km in diameter. These small convective cells strongly contribute to the total area of granules and are located in the intergranular lanes, where they form clusters and chains. We study high-resolution radiation hydrodynamics simulations of the upper convection zone and photosphere to detect small granular cells, define their spatial alignment, and analyze their physical properties. We developed an automated image-segmentation algorithm specifically adapted to high-resolution simulations to identify granules. The resulting segmentation masks were applied to physical quantities, such as intensity and vertical velocity profiles, provided by the simulation. A new clustering algorithm was developed to study the alignment of small granular cells. This study shows that small granules make a distinct contribution to the total area of granules and form clusters of chain-like alignments. The simulation profiles demonstrate a different nature for small granular cells because they exhibit on average lower intensities, lower horizontal velocities, and are located deeper inside of convective layers than regular granules. Their intensity distribution deviates from a normal distribution as known for larger granules, and follows a Weibull distribution.

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 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.

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 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.