No Arabic abstract
We study the combined effects of convection and radiative diffusion on the evolution of thin magnetic flux tubes in the solar interior. Radiative diffusion is the primary supplier of heat to convective motions in the lower convection zone, and it results in a heat input per unit volume of magnetic flux tubes that has been ignored by many previous thin flux tube studies. We use a thin flux tube model subject to convection taken from a rotating spherical shell of turbulent, solar-like convection as described by Weber, Fan, and Miesch (2011, Astrophys. J., 741, 11; 2013, Solar Phys., 287, 239), now taking into account the influence of radiative heating on flux tubes of large-scale active regions. Our simulations show that flux tubes of less than or equal to 60 kG subject to solar-like convective flows do not anchor in the overshoot region, but rather drift upward due to the increased buoyancy of the flux tube earlier in its evolution as a result of the inclusion of radiative diffusion. Flux tubes of magnetic field strengths ranging from 15 kG to 100 kG have rise times of less than or equal to 0.2 years, and exhibit a Joys Law tilt-angle trend. Our results suggest that radiative heating is an effective mechanism by which flux tubes can escape from the stably stratified overshoot region, and that flux tubes do not necessarily need to be anchored in the overshoot region to produce emergence properties similar to those of active regions on the Sun.
Investigation of the turbulent properties of solar convection is extremely important for understanding the multi-scale dynamics observed on the solar surface. In particular, recent high-resolution observations have revealed ubiquitous vortical structures, and numerical simulations have demonstrated links between vortex tube dynamics and magnetic field organization and have shown the importance of vortex tube interactions in the mechanisms of acoustic wave excitation on the Sun. In this paper we investigate the mechanisms of the formation of vortex tubes in highly-turbulent convective flows near the solar surface by using realistic radiative hydrodynamic LES simulations. Analysis of data from the simulations indicates two basic processes of vortex tube formation: 1) development of small-scale convective instability inside convective granules, and 2) a Kelvin-Helmholtz type instability of shearing flows in intergranular lanes. Our analysis shows that vortex stretching during these processes is a primary source of generation of small-scale vorticity on the Sun.
We use a thin flux tube model in a rotating spherical shell of turbulent convective flows to study how active region scale flux tubes rise buoyantly from the bottom of the convection zone to near the solar surface. We investigate toroidal flux tubes at the base of the convection zone with field strengths ranging from 15 kG to 100 kG at initial latitudes ranging from 1 degree to 40 degrees with a total flux of 10^22 Mx. We find that the dynamic evolution of the flux tube changes from convection dominated to magnetic buoyancy dominated as the initial field strength increases from 15 kG to 100 kG. At 100 kG, the development of Omega-shaped rising loops is mainly controlled by the growth of the magnetic buoyancy instability. However, at low field strengths of 15 kG, the development of rising Omega-shaped loops is largely controlled by convective flows, and properties of the emerging loops are significantly changed compared to previous results in the absence of convection. With convection, rise times are drastically reduced (from years to a few months), loops are able to emerge at low latitudes, and tilt angles of emerging loops are consistent with Joys Law for initial field strengths of greater than or equal to 40 kG. We also examine other asymmetries that develop between the leading and following legs of the emerging loops. Taking all the results together, we find that mid-range field strengths of approximately 40 - 50 kG produce emerging loops that best match the observed properties of solar active regions.
(abridged) Context: The mechanisms that cause the formation of sunspots are still unclear. Aims: We study the self-organisation of initially uniform sub-equipartition magnetic fields by highly stratified turbulent convection. Methods: We perform simulations of magnetoconvection in Cartesian domains that are $8.5$-$24$ Mm deep and $34$-$96$ Mm wide. We impose either a vertical or a horizontal uniform magnetic field in a convection-driven turbulent flow. Results: We find that super-equipartition magnetic flux concentrations are formed near the surface with domain depths of $12.5$ and $24$ Mm. The size of the concentrations increases as the box size increases and the largest structures ($20$ Mm horizontally) are obtained in the 24 Mm deep models. The field strength in the concentrations is in the range of $3$-$5$ kG. The concentrations grow approximately linearly in time. The effective magnetic pressure measured in the simulations is positive near the surface and negative in the bulk of the convection zone. Its derivative with respect to the mean magnetic field, however, is positive in the majority of the domain, which is unfavourable for the negative effective magnetic pressure instability (NEMPI). Furthermore, we find that magnetic flux is concentrated in regions of converging flow corresponding to large-scale supergranulation convection pattern. Conclusions: The linear growth of large-scale flux concentrations implies that their dominant formation process is tangling of the large-scale field rather than an instability. One plausible mechanism explaining both the linear growth and the concentrate on of the flux in the regions of converging flow pattern is flux expulsion. Possible reasons for the absence of NEMPI are that the derivative of the effective magnetic pressure with respect to the mean magnetic field has an unfavourable sign and that there may not be sufficient scale separation.
Most 1d hydrodynamic models of plasma confined to magnetic flux tubes assume circular cross-section of these tubes. We use potential field models to show that flux tubes in circumstances relevant to the solar corona do not in general maintain the same cross-sectional shape through their length and therefore the assumption of a circular cross-section is rarely true. We support our hypothesis with mathematical reasoning and numeric experiments. We demonstrate that lifting this assumption in realistic non-circular loops make apparent expansion of magnetic flux tubes consistent with that of observed coronal loops. We propose that in a bundle of ribbon-like loops those that are viewed along the wide direction would stand out against those that are viewed across the wide direction, due to the difference in their column depths. That would impose a bias towards selecting loops that appear not to be expanding seen projected in the plane of sky. An implication of this selection bias is that the preferentially selected non-circular loops would appear to have increased pressure scale height even if they are resolved by current instruments.
Magnetic flux tubes in the solar wind can be twisted as they are transported from the solar surface, where the tubes are twisted owing to photospheric motions. It is suggested that the twisted magnetic tubes can be detected as the variation of total (thermal+magnetic) pressure during their passage through observing satellite. We show that the total pressure of several observed twisted tubes resembles the theoretically expected profile. The twist of isolated magnetic tube may explain the observed abrupt changes of magnetic field direction at tube walls. We have also found some evidence that the flux tube walls can be associated with local heating of the plasma and elevated proton and electron temperatures. For the tubes aligned with the Parker spiral, the twist angle can be estimated from the change of magnetic field direction. Stability analysis of twisted tubes shows that the critical twist angle of the tube with a homogeneous twist is 70$^0$, but the angle can further decrease owing to the motion of the tube with regards to the solar wind stream. The tubes with a stronger twist are unstable to the kink instability, therefore they probably can not reach 1 AU.