We explore the effects of an outer stably stratified coronal envelope on rotating turbulent convection, differential rotation, and large-scale dynamo action in spherical wedge models of the Sun. We solve the compressible magnetohydrodynamic equations
in a two-layer model with unstable stratification below the surface, representing the convection zone, and a stably stratified coronal envelope above. The interface represents a free surface. We compare our model to models that have no coronal envelope. The presence of a coronal envelope is found to modify the Reynolds stress and the $Lambda$ effect resulting in a weaker and non-cylindrical differential rotation. This is related to the reduced latitudinal temperature variations that are caused by and dependent on the angular velocity. Some simulations develop a near-surface shear layer that we can relate to a sign change in the meridional Reynolds stress term in the thermal wind balance equation. Furthermore, the presence of a free surface changes the magnetic field evolution since the toroidal field is concentrated closer to the surface. In all simulations, however, the migration direction of the mean magnetic field can be explained by the Parker--Yoshimura rule, which is consistent with earlier findings. A realistic treatment of the upper boundary in spherical dynamo simulations is crucial for the dynamics of the flow and magnetic field evolution.
We present results from four convectively-driven stellar dynamo simulations in spherical wedge geometry. All of these simulations produce cyclic and migrating mean magnetic fields. Through detailed comparisons we show that the migration direction can
be explained by an $alphaOmega$ dynamo wave following the Parker--Yoshimura rule. We conclude that the equatorward migration in this and previous work is due to a positive (negative) $alpha$ effect in the northern (southern) hemisphere and a negative radial gradient of $Omega$ outside the inner tangent cylinder of these models. This idea is supported by a strong correlation between negative radial shear and toroidal field strength in the region of equatorward propagation.
