Search

Explicit radiation hydrodynamic simulations of the atmospheres of massive stars and of convection in accretion discs around white dwarfs suffer from prohibitively short time steps due to radiation. This constraint is related to the cooling time rather than the radiative pressure, which also becomes important in hot stars and discs. We show that the radiative time step constraint is governed by the minimum of the sum of the optically thick and thin contributions rather than the smaller one of the two. In simulations with the Pencil Code, their weighting fractions are found empirically. In three-dimensional convective accretion disc simulations, the Deardorff term is found to be the main contributor to the enthalpy flux rather than the superadiabatic gradient. We conclude with a discussion of how the radiative time step problem could be mitigated in certain types of investigations. [ABSTRACT FROM AUTHOR]

Hydromagnetic turbulence produced during phase transitions in the early universe can be a powerful source of stochastic gravitational waves (GWs). GWs can be modelled by the linearised spatial part of the Einstein equations sourced by the Reynolds and Maxwell stresses. We have implemented two different GW solvers into the Pencil Code – a code which uses a third order timestep and sixth order finite differences. Using direct numerical integration of the GW equations, we study the appearance of a numerical degradation of the GW amplitude at the highest wavenumbers, which depends on the length of the timestep – even when the Courant–Friedrichs–Lewy condition is ten times below the stability limit. This degradation leads to a numerical error, which is found to scale with the third power of the timestep. A similar degradation is not seen in the magnetic and velocity fields. To mitigate numerical degradation effects, we alternatively use the exact solution of the GW equations under the assumption that the source is constant between subsequent timesteps. This allows us to use a much longer timestep, which cuts the computational cost by a factor of about ten. [ABSTRACT FROM AUTHOR]

We present a high-resolution convergence study of detonation initiated by a temperature gradient in a stoichiometric hydrogen–oxygen mixture using the PENCIL CODE and compare with a code that employs a fifth order weighted essentially non-oscillating (WENO) scheme. With Mach numbers reaching 10–30, a certain amount of shock viscosity is needed in the PENCIL CODE to remove or reduce numerical pressure oscillations on the grid scale at the position of the shock. Detonation is found to occur for intermediate values of the shock viscosity parameter. At fixed values of this parameter, the numerical error associated with those small wiggles in the pressure profile is found to decrease with decreasing mesh width δ x like δ x − 1.4 down to δ x = 0.2 μ m . With the WENO scheme, solutions are smooth at δ x = 10 μ m , but no detonation is obtained for δ x = 5 μ m . This is argued to be an artifact of a decoupling between pressure and reaction fronts. [ABSTRACT FROM AUTHOR]

Recent numerical work in helioseismology has shown that a periodically varying subsurface magnetic field leads to a fanning of the f-mode, which emerges from a density jump at the surface. In an attempt to model a more realistic situation, we now modulate this periodic variation with an envelope, giving thus more emphasis on localised bipolar magnetic structures in the middle of the domain. Some notable findings are: (i) compared to the purely hydrodynamic case, the strength of the f-mode is significantly larger at high horizontal wavenumbers k, but the fanning is weaker for the localised subsurface magnetic field concentrations investigated here than the periodic ones studied earlier; (ii) when the strength of the magnetic field is enhanced at a fixed depth below the surface, the fanning of the f-mode in the k ω diagram increases proportionally in such a way that the normalised f-mode strengths remain nearly the same in different such cases; (iii) the unstable Bloch modes reported previously in case of harmonically varying magnetic fields are now completely absent when more realistic localised magnetic field concentrations are imposed beneath the surface, thus suggesting that the Bloch modes are unlikely to be supported during most phases of the solar cycle; (iv) the f-mode strength appears to depend also on the depth of magnetic field concentrations such that it shows a relative decrement when the maximum of the magnetic field is moved to a deeper layer. We argue that detections of f-mode perturbations such as those being explored here could be effective tracers of solar magnetic fields below the photosphere before these are directly detectable as visible manifestations in terms of active regions or sunspots. [ABSTRACT FROM AUTHOR]

We present the results of simulations of forced turbulence in a slab where the mean kinetic helicity has a maximum near the mid-plane, generating gradients of magnetic helicity of both large and small-scale fields. We also study systems that have poorly conducting buffer zones away from the midplane in order to assess the effects of boundaries. The dynamical α quenching phenomenology requires that the magnetic helicity in the small-scale fields approaches a nearly static, gauge independent state. To stress-test this steady state condition we choose a system with a uniform sign of kinetic helicity, so that the total magnetic helicity can reach a steady state value only through fluxes through the boundary, which are themselves suppressed by the velocity boundary conditions. Even with such a set up, the small-scale magnetic helicity is found to reach a steady state. In agreement with the earlier work, the magnetic helicity fluxes of small-scale fields are found to be turbulently diffusive. By comparing results with and without halos, we show that artificial constraints on magnetic helicity at the boundary do not have a significant impact on the evolution of the magnetic helicity, except that 'softer' (halo) boundary conditions give a lower energy of the saturated mean magnetic field. [ABSTRACT FROM AUTHOR]