  We performed 2.5D magnetohydrodynamic (MHD) simulations showing the
propagation of fast-mode MHD waves of different initial amplitudes and their
interaction with a coronal hole (CH), using our newly developed numerical code.
We find that this interaction results in, first, the formation of reflected,
traversing and transmitted waves (collectively, secondary waves) and, second,
in the appearance of stationary features at the CH boundary. Moreover, we
observe a density depletion that is moving in the opposite direction to the
incoming wave. We find a correlation between the initial amplitude of the
incoming wave and the amplitudes of the secondary waves as well as the peak
values of the stationary features. Additionally, we compare the phase speed of
the secondary waves and the lifetime of the stationary features to
observations. Both effects obtained in the simulation, the evolution of
secondary waves, as well as the formation of stationary fronts at the CH
boundary, strongly support the theory that coronal waves are fast-mode MHD
waves.
