Relativistic effects in the oscillatory damping of magnetic disturbances near two-dimensional X-points are investigated. By taking into account displacement current, we study new features of extremely magnetized systems, in which the Alfv\'en velocity is almost the speed of light. The frequencies of the least-damped mode are calculated using linearized relativistic MHD equations for wide ranges of the Lundquist number S and the magnetization parameter $\sigma$. These timescales approach constant values in the large resistive limit: the oscillation time becomes a few times the light crossing time, irrespective of $\sigma$, and the decay time is proportional to $\sigma$ and therefore is longer for a highly magnetized system.

We present a numerical model in which a cold pair plasma is ejected with relativistic speed through a polar cap region and flows almost radially outside the light cylinder. Stationary axisymmetric structures of electromagnetic fields and plasma flows are self-consistently calculated. In our model, motions of positively and negatively charged particles are assumed to be determined by electromagnetic forces and inertial terms, without pair creation and annihilation or radiation loss. The global electromagnetic fields are calculated by the Maxwell's equations for the plasma density and velocity, without using ideal MHD condition. Numerical result demonstrates the acceleration and deceleration of plasma due to parallel component of the electric fields. Numerical model is successfully constructed for weak magnetic fields or highly relativistic fluid velocity, i.e, kinetic energy dominated outflow. It is found that appropriate choices of boundary conditions and plasma injection model at the polar cap should be explored in order to extend present method to more realistic pulsar magnetosphere, in which the Poynting flux is dominated.