We discuss equilibrium conditions for heterogeneous substances subject to electrostatic or magnetostatic effects. We demonstrate that the force-like aleph tensor N-ij and the energy-like beth tensor beth(ij) for polarizable deformable substances are divergence-free: del N-i(ij) = 0 and del(i)beth(ij)=0. We introduce two additional tensors: the divergence-free energy-like gimel tensor beth(ij) for rigid dielectrics and the general electrostatic gamma tensor Gamma(ij) which is not divergence-free. Our approach is based on a logically consistent extension of the Gibbs energy principle that takes into account polarization effects. While the model is mathematically rigorous, we caution against the assumption that it can reliably predict physical phenomena. On the contrary, clear models often lead to conclusions that are at odds with experiment and therefore should be treated as physical paradoxes that deserve the attention of the scientific community.