A numerical method to determine the electromagnetic field of the pulsar magnetosphere with inclined magnetic moment

A numerical method to determine the electromagnetic field of a steadily rotating magnetosphere with an inclined magnetic moment under a given boundary condition on an arbitrary shaped boundary surface is presented. The region may include the light cylinder. The present method, together with a companion method giving particle motion and creation, makes an iterative scheme to obtain a global model of the pulsar magnetosphere. A key problem for explaining the particle acceleration in pulsars is to solve field-aligned electric field in an accelerating region bounded by an ideal-MHD region. The present method is fit to connect a solution for the non-ideal-MHD region with another solution for the ideal-MHD region on a boundary surface whose location should also be solved (i.e., a floating boundary). The integration scheme is based on the boundary element method and it has great advantage as compared with other methods like the finite difference method and the Fourier transformation method.