Near-polytropic simulations with a radiative surface
by
Atefeh Barekat
→
Europe/Stockholm
122:022
122:022
Description
Studies of solar and stellar convection often employ simple polytropic setups using the diffusion approximation instead of solving the proper radiative transfer equation. This allows one to control separately the polytropic index of the hydrostatic reference solution, the temperature contrast between top and bottom, and the Rayleigh and Peclet numbers. Here we extend such studies by including radiative transfer in the gray approximation using a Kramers-like opacity with freely adjustable coefficients. We study the properties of such models and compare with results from the diffusion approximation. We use the Pencil Code, which is a high-order finite difference code where radiation is treated using the method of long characteristics. The source function is given by the Planck function. The opacity is written as kappa=kappa_0 rho^a T^b. We consider sets of one-dimensional models and perform a comparison with the diffusion approximation in a two-dimensional model. Except for the case with b=5, we find one dimensional hydrostatic equilibria with a nearly polytropic stratification and a polytropic index close to n=(3-b)/(1+a), covering both convectively stable (n>3/2) and unstable (n<3/2) cases. For b=3 and a=-1, the value of n is undefined, but the final equilibrium stratification turns out to be nearly isentropic. For large values of kappa_0, the thermal adjustment time becomes long, the Peclet and Rayleigh numbers become large, and the temperature contrast increases and is thus no longer an independent input parameter, unless the Stefan--Boltzmann constant is considered adjustable. Proper radiative transfer with Kramers-like opacities provides a useful tool for studying stratified layers with a radiative surface in ways that are more physical than what is possible with polytropic models using the diffusion approximation.