Numerical simulation of substation grounding system in low‐frequency domain based on the boundary element method with linear basis function