Solar magnetic field is believed to play a central role in solar activities and flares, filament eruptions as well as CMEs are due to the magnetic field re-organization and the interaction between the plasma and the field. At present the reliable magnetic field measurements are still confined to a few lower levels like in photosphere and chromosphere. Although IR technique may be applied to observe the coronal field but the technique is not well-established yet. Radio techniques may be applied to diagnose the coronal field but assumptions on radiation mechanisms and propagations are needed. Therefore extrapolation from photospheric data upwards is still the primary method to reconstruction coronal field. Potential field has minimum energy content and a force-free field can provide the required excess energy for energy release like flares, etc. Linear models have undesirable properties and it is expected to consider non-constant-alpha force-free field model. As the recent result indicates that the plasma beta is sandwich-ed distributed above the solar surface (Gary, 2001), care must be taken in modeling the coronal field correctly. As the reconstruction of solar coronal magnetic fields is an open boundary problem, it is desired to apply some technique that can incorporate this property. The boundary element method is a well-established numerical techniques that has been applied to many fields including open-space problems. It has also been applied to solar magnetic field problems for potential, linear force-free field and non-constant-alpha force-free field problems. It may also be extended to consider the non-force-free field problem. Here we introduce the procedure of the boundary element method and show its applications in reconstruction of solar magnetic field problems.
