SUMMARYThis paper describes triangular surface charge methods (SCM) using a curved shape function for calculating electric fields in composite dielectrics. In SCM, the boundary condition on dielectric surface is conventionally represented by the point matching of a normal component of electric flux density (D n ), which we call the D p method. We have proposed a method based on the continuity equation of electric flux, which is termed the Φ s method. The electric flux is calculated by numerically integrating D n on each corresponding partial area with proper weight functions. When the continuity of the electric potential (V) is also needed, the point matching of V is adopted conventionally (V p method). We have applied the same integral procedure as the Φ s method to the continuity condition of V by replacing D n with V for the kernel of the integral equation, and we obtained a new representation of potential condition (V s method).We have computed the electric field for a spherical dielectric under a uniform field. The calculated results show that the Φ s method improves the accuracy of the field by about two orders compared with the D p method, and that the V s method gives accuracy comparable to the V p method. The proposed techniques attain the accuracy of the electric field at the spherical center almost equal to the accuracy of the total surface area of the simulating sphere.