Implementation of the electrohydrodynamics’ perfect dielectric model in OpenFOAM®


  • Stefan A. Bošković Innovation Center of the Faculty of Technology and Metallurgy in Belgrade Ltd., Karnegijeva 4, 11120 Belgrade, Serbia
  • Aleksandar Karač Polytechnic Faculty, University of Zenica, Fakultetska 1, 72000 Zenica, Bosnia and Herzegovina
  • Slobodan B. Vrhovac University of Belgrade, Institute of Physics Belgrade, Pregrevica 118, 11080 Belgrade, Serbia
  • Aleksandar Belić University of Belgrade, Institute of Physics Belgrade, Pregrevica 118, 11080 Belgrade, Serbia
  • Ranko Bugarski University of Belgrade, Faculty of Technology and Metallurgy, Karnegijeva 4, 11120 Belgrade, Serbia



The electrohydrodynamics’ (EHD) perfect dielectric model was added into computational fluid dynamics (CFD) software OpenFOAM® in order to improve its usability for the EHD field and specifically for the mentioned model. Based on the investigated literature, it can be said that this is the most complete implementatiton of the said model. Two sets of numerical simulations with two different fluids are presented and analyzed. One set is one-dimensional. The other set is with a drop of one fluid surrounded by other fluid. Oscillations can be observed with certain expressions or calculation strategies for the electrostrictive force, and used for disregarding them. Results that are closer to analytical predictions can be obtained by using appropriate expression for the dielectric force. The electrostrictive force was implemented not only for nonpolar, but also for polar fluids, and it is shown that it might significantly influence the drop deformation. Calculated and analytically predicted drop deformations were close or comparable even up to around 0.25, what is significantly higher and different from a previous study made by other authors. Different expressions for the electric permittivity and usage of limiters for volume fractions were investigated. Conclusions from this paper can be transferred to more complicated models.