Numerical simulations of the kappa-mechanism with convection

A strong coupling between convection and pulsations is known to play a major role in the disappearance of unstable modes close to the red edge of the classical Cepheid instability strip. As mean-field models of time-dependent convection rely on weakly-constrained parameters, we tackle this problem by the means of 2-D Direct Numerical Simulations (DNS) of kappa-mechanism with convection. Using a linear stability analysis, we first determine the physical conditions favourable to the kappa-mechanism to occur inside a purely-radiative layer. Both the instability strips and the nonlinear saturation of unstable modes are then confirmed by the corresponding DNS. We next present the new simulations with convection, where a convective zone and the driving region overlap. The coupling between the convective motions and acoustic modes is then addressed by using projections onto an acoustic subspace.Comment: 5 pages, 6 figures, accepted for publication in Astrophysics and Space Science, HELAS workshop (Rome june 2009

A New Look at Mode Conversion in a Stratified Isothermal Atmosphere

Recent numerical investigations of wave propagation near coronal magnetic null points (McLaughlin and Hood: Astron. Astrophys. 459, 641,2006) have indicated how a fast MHD wave partially converts into a slow MHD wave as the disturbance passes from a low-beta plasma to a high-beta plasma. This is a complex process and a clear understanding of the conversion mechanism requires the detailed investigation of a simpler model. An investigation of mode conversion in a stratified, isothermal atmosphere, with a uniform, vertical magnetic field is carried out, both numerically and analytically. In contrast to previous investigations of upward-propagating waves (Zhugzhda and Dzhalilov: Astron. Astrophys. 112, 16, 1982a; Cally: Astrophys. J. 548, 473, 2001), this paper studies the downward propagation of waves from a low-beta to high-beta environment. A simple expression for the amplitude of the transmitted wave is compared with the numerical solution.Comment: 14 pages, 6 figure

Properties of high-frequency wave power halos around active regions: an analysis of multi-height data from HMI and AIA onboard SDO

We study properties of waves of frequencies above the photospheric acoustic cut-off of $\approx$5.3 mHz, around four active regions, through spatial maps of their power estimated using data from Helioseismic and Magnetic Imager (HMI) and Atmospheric Imaging Assembly (AIA) onboard Solar Dynamics Observatory (SDO). The wavelength channels 1600 {\AA} and 1700 {\AA} from AIA are now known to capture clear oscillation signals due to helioseismic p modes as well as waves propagating up through to the chromosphere. Here we study in detail, in comparison with HMI Doppler data, properties of the power maps, especially the so called 'acoustic halos' seen around active regions, as a function of wave frequencies, inclination and strength of magnetic field (derived from the vector field observations by HMI) and observation height. We infer possible signatures of (magneto-)acoustic wave refraction from the observation height dependent changes, and hence due to changing magnetic strength and geometry, in the dependences of power maps on the photospheric magnetic quantities. We discuss the implications for theories of p mode absorption and mode conversions by the magnetic field.Comment: 22 pages, 12 figures, Accepted by journal Solar Physic

The Ruled Vertex and Nontoric del Pezzo Surfaces

We construct the topological partition function of local nontoric del Pezzo surfaces using the ruled vertex formalism.Comment: 16 pages, 4 figure

Three Dimensional MHD Wave Propagation and Conversion to Alfven Waves near the Solar Surface. I. Direct Numerical Solution

The efficacy of fast/slow MHD mode conversion in the surface layers of sunspots has been demonstrated over recent years using a number of modelling techniques, including ray theory, perturbation theory, differential eigensystem analysis, and direct numerical simulation. These show that significant energy may be transferred between the fast and slow modes in the neighbourhood of the equipartition layer where the Alfven and sound speeds coincide. However, most of the models so far have been two dimensional. In three dimensions the Alfven wave may couple to the magneto-acoustic waves with important implications for energy loss from helioseismic modes and for oscillations in the atmosphere above the spot. In this paper, we carry out a numerical ``scattering experiment'', placing an acoustic driver 4 Mm below the solar surface and monitoring the acoustic and Alfvenic wave energy flux high in an isothermal atmosphere placed above it. These calculations indeed show that energy conversion to upward travelling Alfven waves can be substantial, in many cases exceeding loss to slow (acoustic) waves. Typically, at penumbral magnetic field strengths, the strongest Alfven fluxes are produced when the field is inclined 30-40 degrees from the vertical, with the vertical plane of wave propagation offset from the vertical plane containing field lines by some 60-80 degrees.Comment: Accepted for the HELAS II/ SOHO 19/ GONG 2007 Topical Issue of Solar Physic

Non-linear numerical simulations of magneto-acoustic wave propagation in small-scale flux tubes

We present results of non-linear, 2D, numerical simulations of magneto-acoustic wave propagation in the photosphere and chromosphere of small-scale flux tubes with internal structure. Waves with realistic periods of three to five minutes are studied, after applying horizontal and vertical oscillatory perturbations to the equilibrium model. Spurious reflections of shock waves from the upper boundary are minimized thanks to a special boundary condition. This has allowed us to increase the duration of the simulations and to make it long enough to perform a statistical analysis of oscillations. The simulations show that deep horizontal motions of the flux tube generate a slow (magnetic) mode and a surface mode. These modes are efficiently transformed into a slow (acoustic) mode in the vA < cS atmosphere. The slow (acoustic) mode propagates vertically along the field lines, forms shocks and remains always within the flux tube. It might deposit effectively the energy of the driver into the chromosphere. When the driver oscillates with a high frequency, above the cut-off, non-linear wave propagation occurs with the same dominant driver period at all heights. At low frequencies, below the cut-off, the dominant period of oscillations changes with height from that of the driver in the photosphere to its first harmonic (half period) in the chromosphere. Depending on the period and on the type of the driver, different shock patterns are observed.Comment: 22 pages 6 color figures, submitted to Solar Physics, proceeding of SOHO 19/ GONG 2007 meeting, Melbourne, Australi

Oscillations and waves in solar spicules

Since their discovery, spicules have attracted increased attention as energy/mass bridges between the dense and dynamic photosphere and the tenuous hot solar corona. Mechanical energy of photospheric random and coherent motions can be guided by magnetic field lines, spanning from the interior to the upper parts of the solar atmosphere, in the form of waves and oscillations. Since spicules are one of the most pronounced features of the chromosphere, the energy transport they participate in can be traced by the observations of their oscillatory motions. Oscillations in spicules have been observed for a long time. However the recent high-resolutions and high-cadence space and ground based facilities with superb spatial, temporal and spectral capacities brought new aspects in the research of spicule dynamics. Here we review the progress made in imaging and spectroscopic observations of waves and oscillations in spicules. The observations are accompanied by a discussion on theoretical modelling and interpretations of these oscillations. Finally, we embark on the recent developments made on the presence and role of Alfven and kink waves in spicules. We also address the extensive debate made on the Alfven versus kink waves in the context of the explanation of the observed transverse oscillations of spicule axes

Numerical Simulations of Magnetoacoustic-Gravity Waves in the Solar Atmosphere

We investigate the excitation of magnetoacoustic-gravity waves generated from localized pulses in the gas pressure as well as in vertical component of velocity. These pulses are initially launched at the top of the solar photosphere that is permeated by a weak magnetic field. We investigate three different configurations of the background magnetic field lines: horizontal, vertical and oblique to the gravitational force. We numerically model magnetoacoustic-gravity waves by implementing a realistic (VAL-C) model of solar temperature. We solve two-dimensional ideal magnetohydrodynamic equations numerically with the use of the FLASH code to simulate the dynamics of the lower solar atmosphere. The initial pulses result in shocks at higher altitudes. Our numerical simulations reveal that a small-amplitude initial pulse can produce magnetoacoustic-gravity waves, which are later reflected from the transition region due to the large temperature gradient. The atmospheric cavities in the lower solar atmosphere are found to be the ideal places that may act as a resonator for various oscillations, including their trapping and leakage into the higher atmosphere. Our numerical simulations successfully model the excitation of such wave modes, their reflection and trapping, as well as the associated plasma dynamics

Alfv\'en Reflection and Reverberation in the Solar Atmosphere

Magneto-atmospheres with Alfv\'en speed [a] that increases monotonically with height are often used to model the solar atmosphere, at least out to several solar radii. A common example involves uniform vertical or inclined magnetic field in an isothermal atmosphere, for which the Alfv\'en speed is exponential. We address the issue of internal reflection in such atmospheres, both for time-harmonic and for transient waves. It is found that a mathematical boundary condition may be devised that corresponds to perfect absorption at infinity, and, using this, that many atmospheres where a(x) is analytic and unbounded present no internal reflection of harmonic Alfv\'en waves. However, except for certain special cases, such solutions are accompanied by a wake, which may be thought of as a kind of reflection. For the initial-value problem where a harmonic source is suddenly switched on (and optionally off), there is also an associated transient that normally decays with time as O(t-1) or O(t-1 ln t), depending on the phase of the driver. Unlike the steady-state harmonic solutions, the transient does reflect weakly. Alfv\'en waves in the solar corona driven by a finite-duration train of p-modes are expected to leave such transients.Comment: Accepted by Solar Physic

Structural Invariance of Sunspot Umbrae Over the Solar Cycle: 1993-2004

Measurements of maximum magnetic flux, minimum intensity, and size are presented for 12 967 sunspot umbrae detected on the NASA/NSO spectromagnetograms between 1993 and 2004 to study umbral structure and strength during the solar cycle. The umbrae are selected using an automated thresholding technique. Measured umbral intensities are first corrected for a confirming observation of umbral limb-darkening. Log-normal fits to the observed size distribution confirm that the size spectrum shape does not vary with time. The intensity-magnetic flux relationship is found to be steady over the solar cycle. The dependence of umbral size on the magnetic flux and minimum intensity are also independent of cycle phase and give linear and quadratic relations, respectively. While the large sample size does show a low amplitude oscillation in the mean minimum intensity and maximum magnetic flux correlated with the solar cycle, this can be explained in terms of variations in the mean umbral size. These size variations, however, are small and do not substantiate a meaningful change in the size spectrum of the umbrae generated by the Sun. Thus, in contrast to previous reports, the observations suggest the equilibrium structure, as testified by the invariant size-magnetic field relationship, as well as the mean size (i.e. strength) of sunspot umbrae do not significantly depend on solar cycle phase.Comment: 17 pages, 6 figures. Published in Solar Physic