FEM simulation of microwave absorption by a water drop

The FEM electrical model of a cubic drop (edge=32 10⁻¹⁰ m) of water containing 1536 H20 molecules; i.e. the exact density, is solved by the association Comsol-Matlab. Each water molecule is reduced to a permanent dipole µ=1.8546 D; the O^--atoms are supposed to be fixed and each H⁺⁺ hydrogen is randomly oriented around its oxygen (S. Lefeuvre, O. Gomonova, chap 16, Computational Finite Element Methods in Nanotechnology, CRC Press, 2013). For each run, a Matlab script locates the 3072 points and their electrical charges, following the Comsol order, and then gives the hand to Comsol to solve three times the Poisson equation, corresponding to three boundary conditions: isolation on the six faces (said 0), isolation on 4 faces+ground and V on 2 opposite ones (said 1), isolation on 4+ ground and -V (said 2). The script collects the results, namely the coordinates of each dipole D and the electric field E corresponding to the state´s 0,1,2 and makes the post-processing for instance (E1-E0).D to get the dipole energy. Further, it computes the mean value of the scalar products in a period of time to evaluate the mean microwave power in the drop.