We use 2.5-D electromagnetic particle-in-cell simulation code to investigate the acceleration of electrons in solar coronal holes through the interaction of Alfv\'en waves with an interplume region. The interplume is modeled by a cavity density gradients that are perpendicular to the background magnetic field. The aim is to contribute to explain the observation of suprathermal electrons under relatively quiet sun. Simulations show that Alfv\'en waves in interaction with the interplume region gives rise to a strong local electric field that accelerates electrons in the direction parallel to the background magnetic field. Suprathermal electron beams and small-scale coherent structures are observed within interplume of strong density gradients. These features result from non linear evolution of the electron beam plasma instability.

We use numerical simulations to investigate the interaction of an $f$-mode wave packet with small and large models of a sunspot in a stratified atmosphere. While a loose cluster model has been largely studied before, we focus in this study on the scattering from an ensemble of tightly compact tubes. We showed that the small compact cluster produces a slight distorted scattered wave field in the transverse direction, which can be attributed to the simultaneous oscillations of the pairs of tubes within the cluster aligned in a perpendicular direction to the incoming wave. However, no signature of a multiple-scattering regime has been observed from this model, while it has been clearly observable for the large compact cluster model. Furthermore, we pointed out the importance of the geometrical shape of the monolithic model on the interaction of $f$-mode waves with a sunspot in a high frequency range ($\nu =$ 5 mHz). These results are a contribution to the observational effort to distinguish seismically between different configurations of magnetic flux tubes within sunspots and plage.

Motivated by the problem of local solar subsurface magnetic structure, we have used numerical simulation to investigate the propagation of waves through monolithic magnetic flux tubes of different size. A cluster model can be a good approximation to simulate sunspots as well as solar plage regions which are composed of an ensemble of compactly packed thin flux tubes. Simulations of this type is a powerful tool to probe the structure and the dynamic of various solar features which are related directly to solar magnetic field activity.

Motivated by the question of how to distinguish seismically between monolithic and cluster models of sunspots, we have simulated the propagation of an $f$-mode wave packet through two identical small magnetic flux tubes (R=200 km), embedded in a stratified atmosphere. We want to study the effect of separation $d$ and incidence angle $\chi$ on the scattered wave. We have demonstrated that the horizontal compact pair of tubes ($d/R=2$, $\chi=0$) oscillate as a single tube when the incident wave is propagating, which gives a scattered wave amplitude of about twice that from a single tube. The scattered amplitude decreases with increasing $d$ when $d$ is about $\lambda/2\pi$ where $\lambda$ is the wavelength of the incident wave packet. In this case the individual tubes start to oscillate separately in the manner of near-field scattering. When $d$ is about twice of $\lambda/2\pi$, scattering from individual tubes reaches the far-field regime, giving rise to coherent scattering with an amplitude similar to the case of the compact pair of tubes. For perpendicular incidence ($\chi=\pi/2$), the tubes oscillate simultaneously with the incident wave packet. Moreover, simulations show that a compact cluster oscillates almost as a single individual small tube and acts more like a scattering object, while a loose cluster shows multiple-scattering in the near-field and the absorption is largest when $d$ within the cluster is about $\lambda/2\pi$. This is the first step to understand the seismic response of a bundle of magnetic flux tubes in the context of sunspot and plage helioseismology.

Three-dimensional numerical simulations have been used to study the scattering of a surface-gravity wave packet by vertical magnetic flux tubes, with radii from 200 km to 3 Mm, embedded in stratified polytropic atmosphere. The scattered wave was found to consist primarily of m=0 (axisymmetric) and m=1 modes. It was found that the ratio of the amplitude of these two modes is strongly dependant on the radius of the flux tube: The kink mode is the dominant mode excited in tubes with a small radius while the sausage mode is dominant for large tubes. Simulations of this type provide a simple, efficient and robust way to start understanding the seismic signature of flux tubes, which have recently began to be observed.

In this paper we describe the semi-spectral linear MHD (SLiM) code which we have written to follow the interaction of linear waves through an inhomogeneous three-dimensional solar atmosphere. The background model allows almost arbitrary perturbations of density, temperature, sound speed as well as magnetic and velocity fields. We give details of several of the tests we have used to check the code. The code will be useful in understanding the helioseismic signatures of various solar features, including sunspots.