The purpose of this work is to develop a procedure to obtain the normal modes of a coronal loop from time-dependent numerical simulations with the aim of better understanding observed transverse loop oscillations. To achieve this goal, in this paper we present a new method and test its performance with a problem for which the normal modes can be computed analytically. In a follow-up paper, the application to the simulations of \citet{rial2013} is tackled. The method proceeds iteratively and at each step consists of (i) a time-dependent numerical simulation followed by (ii) the Complex Empirical Orthogonal Function (CEOF) analysis of the simulation results. The CEOF analysis provides an approximation to the normal mode eigenfunctions that can be used to set up the initial conditions for the numerical simulation of the following iteration, in which an improved normal mode approximation is obtained. The iterative process is stopped once the global difference between successive approximate eigenfunctions is below a prescribed threshold. The equilibrium used in this paper contains material discontinuities that result in one eigenfunction with a jump across these discontinuities and two eigenfunctions whose normal derivatives are discontinuous there. After 6 iterations, the approximation to the frequency and eigenfunctions are accurate to $\lesssim 0.7$\% except for the eigenfunction with discontinuities, which displays a much larger error at these positions.

We investigate the temporal evolution of impulsively generated perturbations in a potential coronal arcade with an embedded loop. As the initial configuration we consider a coronal loop, represented by a density enhancement, which is unbounded in the ignorable direction of the arcade. The linearized time-dependent magnetohydrodynamic equations have been numerically solved in field-aligned coordinates and the time evolution of the initial perturbations has been studied in the zero-beta approximation. For propagation constrained to the plane of the arcade, the considered initial perturbations do not excite trapped modes of the system. This weakness of the model is overcome by the inclusion of wave propagation in the ignorable direction. The inclusion of perpendicular propagation produces two main results. First, damping by wave leakage is less efficient because the loop is able to act as a wave trap of vertical oscillations. Second, the consideration of an inhomogeneous corona enables the resonant damping of vertical oscillations and the energy transfer from the interior of the loop to the external coronal medium.

We numerically investigate the excitation and temporal evolution of oscillations in a two-dimensional coronal arcade by including the three-dimensional propagation of perturbations. The time evolution of impulsively generated perturbations is studied by solving the linear, ideal magnetohydrodynamic (MHD) equations in the zero-beta approximation. As we neglect gas pressure the slow mode is absent and therefore only coupled MHD fast and Alfven modes remain. Two types of numerical experiments are performed. First, the resonant wave energy transfer between a fast normal mode of the system and local Alfven waves is analyzed. It is seen how, because of resonant coupling, the fast wave with global character transfers its energy to Alfvenic oscillations localized around a particular magnetic surface within the arcade, thus producing the damping of the initial fast MHD mode. Second, the time evolution of a localized impulsive excitation, trying to mimic a nearby coronal disturbance, is considered. In this case, the generated fast wavefront leaves its energy on several magnetic surfaces within the arcade. The system is therefore able to trap energy in the form of Alfvenic oscillations, even in the absence of a density enhancement such as that of a coronal loop. These local oscillations are subsequently phase-mixed to smaller spatial scales. The amount of wave energy trapped by the system via wave energy conversion strongly depends on the wavelength of perturbations in the perpendicular direction, but is almost independent from the ratio of the magnetic to density scale heights.