Source author record

Friedrich Kupka

Friedrich Kupka appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

astro-ph.SR physics.comp-ph math.NA physics.flu-dyn astro-ph.EP physics.ao-ph physics.geo-ph

Catalog footprint

What is connected

10works

7topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2020arXiv

Thermal Convection in Stars and in Their Atmosphere

Thermal convection is one of the main mechanisms of heat transport and mixing in stars in general and also in the photospheric layers which emit the radiation that we observe with astronomical instruments. The present lecture notes first introduce the role of convection in astrophysics and explain the basic physics of convection. This is followed by an overview on the modelling of convection. Challenges and pitfalls in numerical simulation based modelling are discussed subsequently. Finally, a particular application for the previously introduced concepts is described in more detail: the study of convective overshooting into stably stratified layers around convection zones in stars.

preprint2016arXiv

Multidimensional modelling of classical pulsating stars

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.

preprint2015arXiv

Multidimensional realistic modelling of Cepheid-like variables-II: Analysis of a Cepheid model

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 $α_{\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ß). 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.

preprint2014arXiv

Achievable Efficiency of Numerical Methods for Simulations of Solar Surface Convection

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.

preprint2014arXiv

The effects of rotation on a double-diffusive layer in a rotating spherical shell

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_ρ$ 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énard 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_ρ$. 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.

preprint2013arXiv

Curvilinear Grids for WENO Methods in Astrophysical Simulations

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.

preprint2012arXiv

A Low Mach Number Solver: Enhancing Applicability

In astrophysics and meteorology there exist numerous situations where flows exhibit small velocities compared to the sound speed. To overcome the stringent timestep restrictions posed by the predominantly used explicit methods for integration in time the Euler (or Navier-Stokes) equations are usually replaced by modified versions. In astrophysics this is nearly exclusively the anelastic approximation. Kwatra et al. have proposed a method with favourable time-step properties integrating the original equations (and thus allowing, for example, also the treatment of shocks). We describe the extension of the method to the Navier-Stokes and two-component equations. - However, when applying the extended method to problems in convection and double diffusive convection (semiconvection) we ran into numerical difficulties. We describe our procedure for stabilizing the method. We also investigate the behaviour of Kwatra et al.'s method for very low Mach numbers (down to Ma = 0.001) and point out its very favourable properties in this realm for situations where the explicit counterpart of this method returns absolutely unusable results. Furthermore, we show that the method strongly scales over 3 orders of magnitude of processor cores and is limited only by the specific network structure of the high performance computer we use.

preprint2012arXiv

Multidimensional realistic modelling of Cepheid-like variables. I: Extensions of the ANTARES code

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.

preprint2011arXiv

Total-Variation-Diminishing Implicit-Explicit Runge-Kutta Methods for the Simulation of Double-Diffusive Convection in Astrophysics

We put forward the use of total-variation-diminishing (or more generally, strong stability preserving) implicit-explicit Runge-Kutta methods for the time integration of the equations of motion associated with the semiconvection problem in the simulation of stellar convection. The fully compressible Navier-Stokes equation, augmented by continuity and total energy equations, and an equation of state describing the relation between the thermodynamic quantities, is semi-discretized in space by essentially non-oscillatory schemes and dissipative finite difference methods. It is subsequently integrated in time by Runge-Kutta methods which are constructed such as to preserve the total variation diminishing (or strong stability) property satisfied by the spatial discretization coupled with the forward Euler method. We analyse the stability, accuracy and dissipativity of the time integrators and demonstrate that the most successful methods yield a substantial gain in computational efficiency as compared to classical explicit Runge-Kutta methods.

preprint2010arXiv

Simulations of stellar convection, pulsation and semiconvection

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.