A numerical model for thin liquid film (<100 nm) drainage in the presence of an external electric field is developed. Long-wave theory is applied to approximate and simplify the governing equations. A spatiotemporal film morphology evolution equation thus obtained is then solved using a combination of finite difference to resolve the spatial dimensions and an adaptive time step ODE solver for the temporal propagation. The effect of fluid properties, namely, viscosity and surface tension, on the film drainage time is observed for a homogeneous electric field, which leads to random dewetting spots. Electrically heterogeneous fields, achieved by modeling electrodes with various periodic patterns, are explored to identify their effect on the drainage time and behavior. Finally, the chemical heterogeneity of the substrate is coupled with the periodic electric heterogeneity to understand the implications of combined heterogeneity. It is observed that the introduction of any heterogeneity results in faster drainage of the film when compared to that of the homogeneous field. In all cases, the thin film is drained, leaving submicrometer-scale structures at the interface. Well controlled surface patterns are found on the application of periodic heterogeneity. This study effectively demonstrates the immense potential of electrically induced thin film drainage as a means for faster de-emulsification and for the creation of ordered submicrometer-scale surface patterns on soft materials.