Compressible magnetoconvection in three dimensions: pattern formation in a strongly stratified layer

The interaction between magnetic fields and convection is interesting both because of its astrophysical importance and because the nonlinear Lorentz force leads to an especially rich variety of behaviour. We present several sets of computational results for magnetoconvection in a square box, with periodic lateral boundary conditions, that show transitions from steady convection with an ordered planform through a regime with intermittent bursts to complicated spatiotemporal behaviour. The constraints imposed by the square lattice are relaxed as the aspect ratio is increased. In wide boxes we find a new regime, in which regions with strong fields are separated from regions with vigorous convection. We show also how considerations of symmetry and associated group theory can be used to explain the nature of these transitions and the sequence in which the relevant bifurcations occur

The abrupt development of penumbrae in sunspots

A sunspot is distinguished from a pore by having a filamentary penumbra, corresponding to convective motions that carry energy into the spot from the surrounding field-free plasma. A simplified model of energy transport in sunspots is developed in order to model the transition from pores to spots as the magnetic flux is varied. The observed overlap between the radii of large pores and small spots implies that the filamentary convective mode sets in suddenly and rapidly, as in the idealized case where pore solutions lose stability at a bifurcation

Modelling photospheric magnetoconvection

The increasing power of computers makes it possible to model the nonlinear interaction between magnetic fields and convection at the surfaces of solar-type stars in ever greater detail. We present the results of idealized numerical experiments on two-dimensional magnetoconvection in a fully compressible perfect gas. We first vary the aspect ratio lambda of the computational box and show that the system runs through a sequence of convective patterns, and that it is only for a sufficiently wide box (lambda>=6) that the flow becomes insensitive to further increases in lambda. Next, setting lambda=6, we decrease the field strength from a value strong enough to halt convection and find transitions to small-scale steady convection, next to spatially modulated oscillations (first periodic, then chaotic) and then to a new regime of flux separation, with regions of strong field (where convection is almost completely suppressed) separated by broad convective plumes. We also explore the effects of altering the boundary conditions and show that this sequence of transitions is robust. Finally, we relate these model calculations to recent high-resolution observations of solar magnetoconvection, in plage regions as well as in light bridges and the umbrae of sunspots

Pulsating waves in nonlinear magnetoconvection

Numerical experiments on compressible magnetoconvection reveal a new type of periodic oscillation, associated with alternating streaming motion. Analogous behaviour in a Boussinesq fluid is constrained by extra symmetry. A low-order model confirms that these pulsating waves appear via a pitchfork-Hopf-gluing bifurcation sequence from the steady state

The three-dimensional development of the shearing instability of convection

Two-dimensional convection can become unstable to a mean shear flow. In three dimensions, with periodic boundary conditions in the two horizontal directions, this instability can cause the alignment of convection rolls to alternate between the x and y axes. Rolls with their axes in the y-direction become unstable to a shear flow in the x-direction that tilts and suppresses the rolls, but this flow does not affect rolls whose axes are aligned with it. New rolls, orthogonal to the original rolls, can grow, until they in turn become unstable to the shear flow instability. This behaviour is illustrated both through numerical simulations and through low-order models, and the sequence of local and global bifurcations is determined

The butterfly diagram in the 18th century

Digitized images of the drawings by J.C. Staudacher were used to determine sunspot positions for the period of 1749-1796. From the entire set of drawings, 6285 sunspot positions were obtained for a total of 999 days. Various methods have been applied to find the orientation of the solar disk which is not given for the vast majority of the drawings by Staudacher. Heliographic latitudes and longitudes in the Carrington rotation frame were determined. The resulting butterfly diagram shows a highly populated equator during the first two cycles (Cycles 0 and 1 in the usual counting since 1749). An intermediate period is Cycle 2, whereas Cycles 3 and 4 show a typical butterfly shape. A tentative explanation may be the transient dominance of a quadrupolar magnetic field during the first two cycles.Comment: Accepted for publication in Solar Physics, 1 table, 2 figure

Small-scale magnetic buoyancy and magnetic pumping effects in a turbulent convection

We determine the nonlinear drift velocities of the mean magnetic field and nonlinear turbulent magnetic diffusion in a turbulent convection. We show that the nonlinear drift velocities are caused by the three kinds of the inhomogeneities, i.e., inhomogeneous turbulence; the nonuniform fluid density and the nonuniform turbulent heat flux. The inhomogeneous turbulence results in the well-known turbulent diamagnetic and paramagnetic velocities. The nonlinear drift velocities of the mean magnetic field cause the small-scale magnetic buoyancy and magnetic pumping effects in the turbulent convection. These phenomena are different from the large-scale magnetic buoyancy and magnetic pumping effects which are due to the effect of the mean magnetic field on the large-scale density stratified fluid flow. The small-scale magnetic buoyancy and magnetic pumping can be stronger than these large-scale effects when the mean magnetic field is smaller than the equipartition field. We discuss the small-scale magnetic buoyancy and magnetic pumping effects in the context of the solar and stellar turbulent convection. We demonstrate also that the nonlinear turbulent magnetic diffusion in the turbulent convection is anisotropic even for a weak mean magnetic field. In particular, it is enhanced in the radial direction. The magnetic fluctuations due to the small-scale dynamo increase the turbulent magnetic diffusion of the toroidal component of the mean magnetic field, while they do not affect the turbulent magnetic diffusion of the poloidal field.Comment: 13 pages, 4 figure, REVTEX4, Geophysical and Astrophysical Fluid Dynamics, in pres

Stokes Diagnostis of 2D MHD-simulated Solar Magnetogranulation

We study the properties of solar magnetic fields on scales less than the spatial resolution of solar telescopes. A synthetic infrared spectropolarimetric diagnostics based on a 2D MHD simulation of magnetoconvection is used for this. We analyze two time sequences of snapshots that likely represent two regions of the network fields with their immediate surrounding on the solar surface with the unsigned magnetic flux density of 300 and 140 G. In the first region we find from probability density functions of the magnetic field strength that the most probable field strength at logtau_5=0 is equal to 250 G. Weak fields (B < 500 G) occupy about 70% of the surface, while stronger fields (B 1000 G) occupy only 9.7% of the surface. The magnetic flux is -28 G and its imbalance is -0.04. In the second region, these parameters are correspondingly equal to 150 G, 93.3 %, 0.3 %, -40 G, and -0.10. We estimate the distribution of line-of-sight velocities on the surface of log tau_5=-1. The mean velocity is equal to 0.4 km/s in the first simulated region. The averaged velocity in the granules is -1.2 km/s and in the intergranules is 2.5 km/s. In the second region, the corresponding values of the mean velocities are equal to 0, -1.8, 1.5 km/s. In addition we analyze the asymmetry of synthetic Stokes-V profiles of the Fe I 1564.8 nm line. The mean values of the amplitude and area asymmetry do not exceed 1%. The spatially smoothed amplitude asymmetry is increased to 10% while the area asymmetry is only slightly varied.Comment: 24 pages, 12 figure

Physical Properties of Wave Motion in Inclined Magnetic Fields Within Sunspot Penumbrae

At the surface of the Sun, acoustic waves appear to be affected by the presence of strong magnetic fields in active regions. We explore the possibility that the inclined magnetic field in sunspot penumbrae may convert primarily vertically propagating acoustic waves into elliptical motion. We use helioseismic holography to measure the modulus and phase of the correlation between incoming acoustic waves and the local surface motion within two sunspots. These correlations are modeled assuming the surface motion is elliptical, and we explore the properties of the elliptical motion on the magnetic field inclination. We also demonstrate that the phase shift of the outward propagating waves is opposite to the phase shift of the inward propagating waves in stronger, more vertical fields, but similar to the inward phase shifts in weaker, more inclined fields.Comment: 22 pages, 13 figure

Modeling the Subsurface Structure of Sunspots

While sunspots are easily observed at the solar surface, determining their subsurface structure is not trivial. There are two main hypotheses for the subsurface structure of sunspots: the monolithic model and the cluster model. Local helioseismology is the only means by which we can investigate subphotospheric structure. However, as current linear inversion techniques do not yet allow helioseismology to probe the internal structure with sufficient confidence to distinguish between the monolith and cluster models, the development of physically realistic sunspot models are a priority for helioseismologists. This is because they are not only important indicators of the variety of physical effects that may influence helioseismic inferences in active regions, but they also enable detailed assessments of the validity of helioseismic interpretations through numerical forward modeling. In this paper, we provide a critical review of the existing sunspot models and an overview of numerical methods employed to model wave propagation through model sunspots. We then carry out an helioseismic analysis of the sunspot in Active Region 9787 and address the serious inconsistencies uncovered by \citeauthor{gizonetal2009}~(\citeyear{gizonetal2009,gizonetal2009a}). We find that this sunspot is most probably associated with a shallow, positive wave-speed perturbation (unlike the traditional two-layer model) and that travel-time measurements are consistent with a horizontal outflow in the surrounding moat.Comment: 73 pages, 19 figures, accepted by Solar Physic