Numerical calculation of effective permittivity of lossless dielectric mixtures using boundary integral method

The electrostatic environment in a lossless composite structure made of three dielectrics and two dividing closed and disjoined surfaces has been determined, using appropriate surface integrals on these interfaces. These integral equations are transformed to approximate matrix equations and then are solved numerically. The electrostatic potential distribution, dipole moment and the effective permittivity of different composite structures are calculated.>