The dielectric properties of heterogeneous materials for various condensed-mattersystems are important for several technologies, e.g. impregnated polymers forhigh-density capacitors, polymer carbon black mixtures for automotive tires and currentlimiters in circuit protection. These multiscale systems lead to challenging problemsof connecting microstructural features (shape, spatial arrangement and sizedistribution of inclusions) to macroscopic materials response (permittivity,conductivity). In this paper, we briefly discuss an ab initio computationalelectrostatics approach, based either on the use of the field calculation packageFLUX3D (or FLUX2D) and a conventional finite elements method, or the use of the fieldcalculation package PHI3D and the resolution of boundary integral equations, forcalculating the effective permittivity of two-component dielectric heterostructures.Numerical results concerning inclusions of permittivity ε 1 with variousgeometrical shapes periodically arranged in a host matrix of permittivityε 2 are provided. Next we discuss these results in terms ofphenomenological mixing laws, analytical theory and connectedness. During the pursuitof these activities, several interesting phenomena were discovered that will stimulatefurther investigation.