The problem of phase mixed shear Alfvén waves in the solar corona revisited

The problem of phase mixing of shear Alfvén waves is revisited taking into account dissipative phenomena specific for the solar corona. In regions of space plasmas where the dynamics is controlled by the magnetic field, transport coefficients become anisotropic with transport mechanism having different behavior and magnitude depending on the orientation with respect to the ambient magnetic field. Taking into account realistic values for dissipative coefficients we obtain that the previous results derived in context of torsional Alfvén wave phase mixing are actually heavily underestimated so phase mixing cannot be used to explain the damping of torsional Alfvén waves and heating of open coronal structures. The presented results indicate that in order for phase mixing to still be a viable mechanism to explain heating or wave damping unrealistic assumptions have to be made. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

MHD waves at a spherical interface modelling EIT waves

Eruptive events such as flares and coronal mass ejections (CMEs) are known to generate global waves propagating over distances comparable to the solar radius in different layers of the solar atmosphere. Here we investigate the propagation of coronal EIT waves, modelled as fast magnetoacoustic modes propagating at a spherical interface in the presence of a purely radial magnetic field. Based on a simplified equilibrium we derive the dispersion relation of the waves. The generation and propagation of EIT waves at the spherical interface is studied numerically for different values of spherical degree and preliminary conclusions are reached regarding the properties of EIT waves. (© 2007 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Catastrophe of coronal magnetic rope in partly open multipolar magnetic field

Catastrophe of coronal magnetic rope embedded in a partly open multipolar background magnetic field is studied by using a 2-dimensional, 3-component ideal MHD model in spherical coordinates. The background field is composed of three closed bipolar fields of a coronal streamer and an open field with an equatorial current sheet. The magnetic rope lies below the central bipolar field, and it is characterized by its annular and axial magnetic fluxes. For a given annual flux, there is a critical value of the axial flux, and for a given axial flux, there is a critical value of annual flux such that, below the critical value, the magnetic rope is attached to the solar surface and the system stays in equilibrium, but when the critical value is exceeded, the magnetic rope breaks free and erupts upward. This implies that catastrophe can occur in a coronal magnetic rope embedded in a partly open multipolar background magnetic field. Our computation gives a threshold value of magnetic energy that is about 15% greater than the energy of the partly open magnetic field (the central bipolar field open and the fields on either side closed). The excess energy may serve as source for solar explosions such as coronal mass ejections.     

Solar dynamo models with α ‐effect and turbulent pumping from local 3D convection calculations

Results from kinematic solar dynamo models employing α ‐effect and turbulent pumping from local convection calculations are presented. We estimate the magnitude of these effects to be around 2–3 m s^–1, having scaled the local quantities with the convective velocity at the bottom of the convection zone from a solar mixing‐length model. Rotation profile of the Sun as obtained from helioseismology is applied in the models; we also investigate the effects of the observed surface shear layer on the dynamo solutions. With these choices of the small‐ and large‐scale velocity fields, we obtain estimate of the ratio of the two induction effects, C _α /C _Ω ≈ 10^–3, which we keep fixed in all models. We also include a one‐cell meridional circulation pattern having a magnitude of 10–20 m s^–1 near the surface and 1–2 m s^–1 at the bottom of the convection zone. The model essentially represents a distributed turbulent dynamo, as the α ‐effect is nonzero throughout the convection zone, although it concentrates near the bottom of the convection zone obtaining a maximum around 30° of latitude. Turbulent pumping of the mean fields is predominantly down‐ and equatorward. The anisotropies in the turbulent diffusivity are neglected apart from the fact that the diffusivity is significantly reduced in the overshoot region. We find that, when all these effects are included in the model, it is possible to correctly reproduce many features of the solar activity cycle, namely the correct equatorward migration at low latitudes and the polar branch at high latitudes, and the observed negative sign of B _r B _ϕ . Although the activity clearly shifts towards the equator in comparison to previous models due to the combined action of the α ‐effect peaking at midlatitudes, meridional circulation and latitudinal pumping, most of the activity still occurs at too high latitudes (between 5° … 60°). Other problems include the relatively narrow parameter space within which the preferred solution is dipolar (A0), and the somewhat too short cycle lengths of the solar‐type solutions. The role of the surface shear layer is found to be important only in the case where the α ‐effect has an appreciable magnitude near the surface. (© 2006 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

MHD waves in the solar north polar coronal hole

The effects, hitherto not treated, of the temperature and the number density gradients, both in the parallel and the perpendicular direction to the magnetic field, of O VI ions, on the MHD wave propagation characteristics in the solar North Polar Coronal Hole are investigated. We investigate the magnetosonic wave propagation in a resistive MHD regime where only the thermal conduction is taken into account. Heat conduction across the magnetic field is treated in a non‐classical approach wherein the heat is assumed to be conducted by the plasma waves emitted by ions and absorbed at a distance from the source by other ions. Anisotropic temperature and the number density distributions of O VI ions revealed the chaotic nature of MHD standing wave, especially near the plume/interplume lane borders. Attenuation length scales of the fast mode is shown not to be smoothly varying function of the radial distance from the Sun (© 2010 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

The origin of cutoff frequencies for torsional tube waves propagating in the solar atmosphere

Torsional waves supported by magnetic flux tubes have long been thought to bear a high potential for supplying energy and momentum to the upper solar atmosphere, thereby contributing to its heating and to the driving of dynamic events like spicules. This hope rested on the belief that their propagation is not impeded by cutoff restrictions, unlike longitudinal and kink waves. We point out that this applies only to thin, isothermal tubes. When they widen in the chromosphere, and as a result of temperature gradients, cutoff restrictions arise. We compare them to recent observational reports of such waves and of vortex motions and find that their long period components are already affected by cutoff restrictions. An observational strategy is proposed that should permit the derivation of better information on vortex flows from off‐center observations with next generation telescopes (© 2010 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Modeling a Maunder minimum

We introduce on/off intermittency into a mean field dynamo model by imposing stochastic fluctuations in either the alpha effect or through the inclusion of a fluctuating electromotive force. Sufficiently strong small scale fluctuations with time scales of the order of 0.3–3 years can produce long term variations in the system on time scales of the order of hundreds of years. However, global suppression of magnetic activity in both hemispheres at once was not observed. The variation of the magnetic field does not resemble that of the sunspot number, but is more reminiscent of the ¹⁰Be record. The interpretation of our results focuses attention on the connection between the level of magnetic activity and the sunspot number, an issue that must be elucidated if long term solar effects are to be well understood. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Magneto‐seismology of the solar atmosphere

Magnetohydrodynamic (MHD) waves in solar coronal loops, which were previously only predicted by theory have actually been detected with space‐borne instruments. These observations have given an important and novel tool to measure fundamental parameters in the magnetically embedded solar corona. This paper will illustrate how information about the magnetic and density structure along coronal loops can be determined by measuring the frequency or amplitude profiles of standing fast kink mode oscillations. (© 2007 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Alleviation of catastrophic quenching in solar dynamo model with nonlocal alpha‐effect

The nonlocal alpha‐effect of Babcock‐Leighton type is not prone to the catastrophic quenching due to conservation of magnetic helicity. This is shown with a dynamo model, which jointly applies the nonlocal alpha‐effect, the diamagnetic pumping, and dynamical equation for the magnetic alpha‐effect. The same model shows catastrophic quenching when the alpha‐effect is changed to its local formulation. The nonlocal model shows a preferred excitation of magnetic fields of dipolar symmetry, which oscillate with a period of about ten years and have a toroidal‐to‐polar fields ratio of about a thousand (© 2011 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

3-D non-linear evolution of a magnetic flux tube in a spherical shell: The isentropic case

We present recent 3-D MHD numerical simulations of the non-linear dynamical evolution of magnetic flux tubes in an adiabatically stratified convection zone in spherical geometry, using the anelastic spherical harmonic (ASH) code.We seek to understand the mechanism of emergence of strong toroidal fields from the base of the solar convection zone to the solar surface as active regions. We confirm the results obtained in cartesian geometry that flux tubes that are not twisted split into two counter vortices before reaching the top of the convection zone. Moreover, we find that twisted tubes undergo the poleward-slip instability due to an unbalanced magnetic curvature force which gives the tube a poleward motion both in the non-rotating and in the rotating case. This poleward drift is found to be more pronounced on tubes originally located at high latitudes. Finally, rotation is found to decrease the rise velocity of the flux tubes through the convection zone, especially when the tube is introduced at low latitudes. (© 2007 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Convective dynamos in spherical wedge geometry

Self‐consistent convective dynamo simulations in wedge‐shaped spherical shells are presented. Differential rotation is generated by the interaction of convection with rotation. Equatorward acceleration and dynamo action are obtained only for sufficiently rapid rotation. The angular velocity tends to be constant along cylinders. Oscillatory large‐scale fields are found to migrate in the poleward direction. Comparison with earlier simulations in full spherical shells and Cartesian domains is made (© 2010 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

太阳黑子本影振荡的多波段观测

基于多波段的观测数据, 研究了黑子本影振荡在太阳大气不同高度或温度的振荡特征. 目标黑子发生在2017年9月15日, 位于活动区12680. 多波段的观测数据包括AIA (Atmospheric Imaging Assembly)的极紫外成像及IRIS (Interface Region Imaging Spectrograph)的紫外光谱和成像. 在太阳黑子的本影位置, AIA 1700 Å 光变曲线的振荡周期约为(4.2士0.8) min,而AIA 1600 Å、171 Å和193 Å光变曲线的振荡周期约为(2.8士0.3)min. IRIS的Mg II h和k谱线及远紫外成像(2796 Å)都表现为(3.1士0.5) min的准周期振荡,而谱线Si IV 1393.76 Å的振荡周期约为(2.9+0.4)min.观测结果说明黑子本影的振荡周期随着太阳大气高度的升高而逐渐减小.较短的周期可认为是3 min振荡,很可能与黑子本影上方的慢磁声波有关,它起源于温度极小区并向上传播到日冕.较长的周期可解释为光球的5 min振荡,与太阳的P模振荡有关.     

Dynamo action in simulations of penetrative solar convection with an imposed tachocline

We summarize new and continuing three-dimensional spherical shell simulations of dynamo action by convection allowed to penetrate downward into a tachocline of rotational shear. The inclusion of an imposed tachocline allows us to examine several processes believed to be essential in the operation of the global solar dynamo, including differential rotation, magnetic pumping, and the stretching and organization of fields within the tachocline. In the stably stratified core, our simulations reveal that strong axisymmetric magnetic fields (of ∼ 3000 G strength) can be built, and that those fields generally exhibit a striking antisymmetric parity, with fields in the northern hemisphere largely of opposite polarity to those in the southern hemisphere. In the convection zone above, fluctuating fields dominate over weaker mean fields. New calculations indicate that the tendency toward toroidal fields of antisymmetric parity is relatively insensitive to initial magnetic field configurations; they also reveal that on decade-long timescales, the magnetic fields can briefly enter (and subsequently emerge from) states of symmetric parity.We have not yet observed any overall reversals of the field polarity, nor systematic latitudinal propagation. (© 2007 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Equatorial magnetic helicity flux in simulations with different gauges

We use direct numerical simulations of forced MHD turbulence with a forcing function that produces two different signs of kinetic helicity in the upper and lower parts of the domain. We show that the mean flux of magnetic helicity from the small‐scale field between the two parts of the domain can be described by a Fickian diffusion law with a diffusion coefficient that is approximately independent of the magnetic Reynolds number and about one third of the estimated turbulent magnetic diffusivity. The data suggest that the turbulent diffusive magnetic helicity flux can only be expected to alleviate catastrophic quenching at Reynolds numbers of more than several thousands. We further calculate the magnetic helicity density and its flux in the domain for three different gauges. We consider the Weyl gauge, in which the electrostatic potential vanishes, the pseudo‐Lorenz gauge, where the speed of light is replaced by the sound speed, and the ‘resistive gauge’ in which the Laplacian of the magnetic vector potential acts as a resistive term. We find that, in the statistically steady state, the time‐averaged magnetic helicity density and the magnetic helicity flux are the same in all three gauges (© 2010 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

On global solar dynamo simulations

Global dynamo simulations solving the equations of magnetohydrodynamics (MHD) have been a tool of astrophysicists who try to understand the magnetism of the Sun for several decades now. During recent years many fundamental issues in dynamo theory have been studied in detail by means of local numerical simulations that simplify the problem and allow the study of physical effects in isolation. Global simulations, however, continue to suffer from the age‐old problem of too low spatial resolution, leading to much lower Reynolds numbers and scale separation than in the Sun. Reproducing the internal rotation of the Sun, which plays a crucial role in the dynamo process, has also turned out to be a very difficult problem. In the present paper the current status of global dynamo simulations of the Sun is reviewed. Emphasis is put on efforts to understand how the large‐scale magnetic fields, i.e. whose length scale is greater than the scale of turbulence, are generated in the Sun. Some lessons from mean‐field theory and local simulations are reviewed and their possible implications to the global models are discussed. Possible remedies to some current issues of solar simulations are put forward (© 2011 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

耦合热传导的磁重联的数值模拟

磁重联被认为是太阳耀斑的产生机制,本文数值模拟在日冕中发生在磁重联过程,结果表明耀斑环的表观运动是磁重联的自洽结果;由重联点发出的慢激波对耀斑环的加热有贡献;耀斑环的上升并不意味着重联点的上升。     

Meridional circulation and dynamo waves

The effect of meridional circulation on spherical shell dynamos is considered in the Parker approximation. We demonstrate that the type of the exited solution crucially depend on the intensity of the meridional circulation. If the circulation is equatorward or if it is polarward however do not exceed some critical value, an oscillating solution in form of an equatorward traveling wave is excited. If the polarward meridional circulation becomes too intensive the solution becomes steady growing. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

耀斑大尺度电流片的湍流耗散和能谱分析

磁重联在宇宙的许多动力学现象中都是非常核心的过程.磁流体动力学(MHD)数值模拟是研究磁重联过程以及相应物理图像的一种很有效的手段.通过不同的参数组合,来研究MHD数值模拟中磁雷诺数和空间分辨率对磁重联率、数值耗散和能谱分布的影响.对得到的数据进行分析后,发现磁雷诺数对磁重联率和能谱分布有一定的影响.磁雷诺数越大,磁重联过程进入非线性阶段所需的特征时间越短,磁重联率就越早发生跃升.磁雷诺数R_m对耗散开始发挥作用的Kolmogorov微观尺度l_ko有明显影响:R_m越大,l_ko就越小.研究了磁重联过程中包括数值耗散在内的额外耗散对重联过程的影响.结果表明,撕裂模不稳定性开始之前的额外耗散以纯数值耗散为主,撕裂模不稳定性出现之后,额外耗散出现同步跃升,说明不稳定性导致的湍流明显增强了耗散的效果,相当于在局部湍流区引入了超电阻.能谱分析进一步表明,大尺度电流片的l_ko完全可能出现在宏观的MHD尺度上.