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.