Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userSpherical-shell boundaries for two-dimensional compressible convection in a star
77
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
Context: We study the impact of two-dimensional spherical shells on compressible convection. Realistic profiles for density and temperature from a one-dimensional stellar evolution code are used to produce a model of a large stellar convection zone representative of a young low-mass star. Methods: We perform hydrodynamic implicit large-eddy simulations of compressible convection using the MUltidimensional Stellar Implicit Code (MUSIC). Because MUSIC has been designed to use realistic stellar models produced from one-dimensional stellar evolution calculations, MUSIC simulations are capable of seamlessly modeling a whole star. Simulations in two-dimensional spherical shells that have different radial extents are performed over hundreds of convective turnover times, permitting the collection of well-converged statistics. Results: We evaluate basic statistics of the convective turnover time, the convective velocity, and the overshooting layer. These quantities are selected for their relevance to one-dimensional stellar evolution calculations, so that our results are focused toward the 321D link. The inclusion in the spherical shell of the boundary between the radiative and convection zones decreases the amplitude of convective velocities in the convection zone. The inclusion of near-surface layers in the spherical shell can increase the amplitude of convective velocities, although the radial structure of the velocity profile established by deep convection is unchanged. The impact from including the near-surface layers depends on the speed and structure of small-scale convection in the near-surface layers. Larger convective velocities in the convection zone result in a commensurate increase in the overshooting layer width and decrease in the convective turnover time. These results provide support for non-local aspects of convection.
rate research
Read More
Extending our recent studies of two-dimensional stellar convection to 3D, we compare three-dimensional hydrodynamic simulations to identically set-up two-dimensional simulations, for a realistic pre-main sequence star. We compare statistical quantities related to convective flows including: average velocity, vorticity, local enstrophy, and penetration depth beneath a convection zone. These statistics are produced during stationary, steady-state compressible convection in the stars convection zone. Our simulations with the MUSIC code confirm the common result that two-dimensional simulations of stellar convection have a higher magnitude of velocity on average than three-dimensional simulations. Boundary conditions and the extent of the spherical shell can affect the magnitude and variability of convective velocities. The difference between 2D and 3D velocities is dependent on these background points; in our simulations this can have an effect as large as the difference resulting from the dimensionality of the simulation. Nevertheless, radial velocities near the convective boundary are comparable in our 2D and 3D simulations. The average local enstrophy of the flow is lower for two-dimensional simulations than for three-dimensional simulations, indicating a different shape and structuring of 3D stellar convection. We perform a statistical analysis of the depth of convective penetration below the convection zone, using the model proposed in our recent study (Pratt et al. 2017). Here we analyze the convective penetration in three dimensional simulations, and compare the results to identically set-up 2D simulations. In 3D the penetration depth is as large as the penetration depth calculated from 2D simulations.
We report the finding of an azimuthal dynamo wave of a low-order (m=1) mode in direct numerical simulations (DNS) of turbulent convection in spherical shells. Such waves are predicted by mean field dynamo theory and have been obtained previously in mean-field models. Observational results both from photometry and Doppler imaging have revealed persistent drifts of spots for several rapidly rotating stars, but, although an azimuthal dynamo wave has been proposed as a possible mechanism responsible for this behavior, it has been judged as unlikely, as practical evidence for such waves from DNS has been lacking. The large-scale magnetic field in our DNS, which is due to self-consistent dynamo action, is dominated by a retrograde m=1 mode. Its pattern speed is nearly independent of latitude and does not reflect the speed of the differential rotation at any depth. The extrema of magnetic m=1 structures coincide reasonably with the maxima of m=2 structures of the temperature. These results provide direct support for the observed drifts being due to an azimuthal dynamo wave.
We examine a penetration layer formed between a central radiative zone and a large convection zone in the deep interior of a young low-mass star. Using the Multidimensional Stellar Implicit Code (MUSIC) to simulate two-dimensional compressible stellar convection in a spherical geometry over long times, we produce statistics that characterize the extent and impact of convective penetration in this layer. We apply extreme value theory to the maximal extent of convective penetration at any time. We compare statistical results from simulations which treat non-local convection, throughout a large portion of the stellar radius, with simulations designed to treat local convection in a small region surrounding the penetration layer. For each of these situations, we compare simulations of different resolution, which have different velocity magnitudes. We also compare statistical results between simulations that radiate energy at a constant rate to those that allow energy to radiate from the stellar surface according to the local surface temperature. Based on the frequency and depth of penetrating convective structures, we observe two distinct layers that form between the convection zone and the stable radiative zone. We show that the probability density function of the maximal depth of convective penetration at any time corresponds closely in space with the radial position where internal waves are excited. We find that the maximal penetration depth can be modeled by a Weibull distribution with a small shape parameter. Using these results, and building on established scalings for diffusion enhanced by large-scale convective motions, we propose a new form for the diffusion coefficient that may be used for one-dimensional stellar evolution calculations in the large Peclet number regime. These results should contribute to the 321D link.
(abidged) Context: Stellar convection zones are characterized by vigorous high-Reynolds number turbulence at low Prandtl numbers. Aims: We study the dynamo and differential rotation regimes at varying levels of viscous, thermal, and magnetic diffusion. Methods: We perform three-dimensional simulations of stratified fully compressible magnetohydrodynamic convection in rotating spherical wedges at various thermal and magnetic Prandtl numbers (from 0.25 to 2 and 5, respectively). Results: We find that the rotation profiles for high thermal diffusivity show a monotonically increasing angular velocity from the bottom of the convection zone to the top and from the poles toward the equator. For sufficiently rapid rotation, a region of negative radial shear develops at mid-latitudes as the thermal diffusivity is decreased. This coincides with a change in the dynamo mode from poleward propagating activity belts to equatorward propagating ones. Furthermore, the cyclic solutions disappear at the highest magnetic Reynolds numbers. The total magnetic energy increases with the magnetic Reynolds number in the range studied here ($5-151$), but the energies of the mean magnetic fields level off at high magnetic Reynolds numbers. The differential rotation is strongly affected by the magnetic fields and almost vanishes at the highest magnetic Reynolds numbers. In some of our most turbulent cases we find that two regimes are possible where either differential rotation is strong and mean magnetic fields relatively weak or vice versa. Conclusions: Our simulations indicate a strong non-linear feedback of magnetic fields on differential rotation, leading to qualitative changes in the behaviors of large-scale dynamos at high magnetic Reynolds numbers. Furthermore, we do not find indications of the simulations approaching an asymptotic regime where the results would be independent of diffusion coefficients.
We test the sensitivity of hydrodynamic and magnetohydrodynamic turbulent convection simulations with respect to Mach number, thermal and magnetic boundary conditions, and the centrifugal force. We find that varying the luminosity, which also controls the Mach number, has only a minor effect on the large-scale dynamics. A similar conclusion can also be drawn from the comparison of two formulations of the lower magnetic boundary condition with either vanishing electric field or current density. The centrifugal force has an effect on the solutions, but only if its magnitude with respect to acceleration due to gravity is by two orders of magnitude greater than in the Sun. Finally, we find that the parameterisation of the photospheric physics, either by an explicit cooling term or enhanced radiative diffusion, is more important than the thermal boundary condition. In particular, runs with cooling tend to lead to more anisotropic convection and stronger deviations from the Taylor-Proudman state. In summary, the fully compressible approach taken here with the Pencil Code is found to be valid, while still allowing the disparate timescales to be taken into account.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
