The nonlinear dynamics of outflows driven by magnetic explosion on the surface of a compact star is investigated through special relativistic magnetohydrodynamic simulations. We adopt, as the initial equilibrium state, a spherical stellar object embedded in hydrostatic plasma which has a density $\rho(r) \propto r^{- \alpha}$ and is threaded by a dipole magnetic field. The injection of magnetic energy at the surface of compact star breaks the equilibrium and triggers a two-component outflow. At the early evolutionary stage, the magnetic pressure increases rapidly around the stellar surface, initiating a magnetically driven outflow. A strong forward shock driven outflow is then excited. The expansion velocity of the magnetically driven outflow is characterized by the Alfv\'en velocity on the stellar surface, and follows a simple scaling relation $v_{\rm mag} \propto {v_{\rm A}}^{1/2}$. When the initial density profile declines steeply with radius, the strong shock is accelerated self-similarly to relativistic velocity ahead of the magnetically driven component. We find that it evolves according to a self-similar relation $\Gamma_{\rm sh} \propto r_{\rm sh}$, where $\Gamma_{\rm sh}$ is the Lorentz factor of the plasma measured at the shock surface $r_{\rm sh}$. Purely hydrodynamic process would be responsible for the acceleration mechanism of the shock driven outflow. Our two-component outflow model, which is the natural outcome of the magnetic explosion, can provide a better understanding of the magnetic active phenomena on various magnetized compact stars.

We obtained self-similar solutions of relativistically expanding magnetic loops by assuming axisymmetry and a purely radial flow. The stellar rotation and the magnetic fields in the ambient plasma are neglected. We include the Newtonian gravity of the central star. These solutions are extended from those in our previous work (Takahashi, Asano, & Matsumoto 2009) by taking into account discontinuities such as the contact discontinuity and the shock. The global plasma flow consists of three regions, the outflowing region, the post shocked region, and the ambient plasma. They are divided by two discontinuities. The solutions are characterized by the radial velocity, which plays a role of the self-similar parameter in our solutions. The shock Lorentz factor gradually increases with radius. It can be approximately represented by the power of radius with the power law index of 0.25. We also carried out magnetohydrodynamic simulations of the evolution of magnetic loops to study the stability and the generality of our analytical solutions. We used the analytical solutions as the initial condition and the inner boundary conditions. We confirmed that our solutions are stable over the simulation time and that numerical results nicely recover the analytical solutions. We then carried out numerical simulations to study the generality of our solutions by changing the power law index \delta of the ambient plasma density \rho_0 \propto r^{-\delta}. We alter the power law index \delta from 3.5 in the analytical solutions. The analytical solutions are used as the initial conditions inside the shock in all simulations. We observed that the shock Lorentz factor increases with time when \delta is larger than 3, while it decreases with time when \delta is smaller than 3. The shock Lorentz factor is proportional to t^{(\delta-3)/2}. These results are consistent with the analytical studies by Shapiro (1979).

We obtained self-similar solutions of relativistically expanding magnetic loops taking into account the azimuthal magnetic fields. We neglect stellar rotation and assume axisymmetry and a purely radial flow. As the magnetic loops expand, the initial dipole magnetic field is stretched into the radial direction. When the expansion speed approaches the light speed, the displacement current reduces the toroidal current and modifies the distribution of the plasma lifted up from the central star. Since these self-similar solutions describe the free expansion of the magnetic loops, i.e., $Dv/Dt=0$, the equations of motion are similar to those of the static relativistic magnetohydrodynamics. This allows us to estimate the total energy stored in the magnetic loops by applying the virial theorem. This energy is comparable to that of the giant flares observed in magnetars.

We study the magnetic interaction between a neutron star and its disk by solving the time-dependent relativistic force-free equations. At the initial state, we assume that the dipole magnetic field of the neutron star connects the neutron star and its equatorial disk, which deeply enters into the magnetosphere of the neutron star. Magnetic fields are assumed to be frozen to the star and the disk. The rotation of the neutron star and the disk is imposed as boundary conditions. We apply Harten-Lax-van Leer (HLL) method to simulate the evolution of the star-disk system. We carry out simulations for (1) a disk inside the corotation radius, in which the disk rotates faster than the star, and (2) a disk outside the corotation radius, in which the neutron star rotates faster than the disk. Numerical results indicate that for both models, the magnetic field lines connecting the disk and the star inflate as they are twisted by the differential rotation between the disk and the star. When the twist angle exceeds pi radian, the magnetic field lines expand with speed close to the light speed. This mechanism can be the origin of relativistic outflows observed in binaries containing a neutron star.