Electrostatic levitations of anisotropic dielectric oblate and prolate spheroids

The paper presents analytical approach to electrostatic field for dielectric spheroids problem. For both oblate and prolate spheroids fields distributions are derived. The variable separation method is applied. Electrostatic force is evaluated with the help of the Maxwell stress tensor generalized method, material force density, coenergy and equivalent dipole. Levitation forces for both oblate and prolate dielectric spheroids versus axis permittivities and spheroids height are presented.


List of symbols
Imposed electric field strength along zaxis (constant item) E η , E θ ,E ϕ Electric field strength components F 1 Imposed electric field gradient at z 0 axis V m -2 g The acceleration of gravity of free space or its changed value due to acting upward buoyant force in fluids h Height of spheroid, semi-axis parallel to the z-axis i u Versor for uth coordinates u η, θ , ϕ Lame coefficients for spheroidal coordinates: oblate L η L θ c e ch 2 (η) − sin 2  The Maxwell stress tensor component for u−w axes

Introduction
Levitation of samples in electrostatic field is one of possible method for arising objects in gravitational field by means of electrostatic forces (Table 1).Nowadays, technology is often focused on electrostatic separations in medicament treatments, levitations of liquid drops and particle, alloys and some laboratory investigations [3,4,6,7,14].There are also technologies that apply static levitations arising without capacitances [8].The paper considers electrostatic levitation driven only by material force acting on dielectric samples.Nevertheless, the electrostatic levitation force may act together with another electromagnetic forces [10,11].
The material forces of electric nature act at the boundary of spheroid (Table 1).While material forces are acting on dielectric spheroid, they could raise it in gravitational field [1,5,9].The spheroid levitation is feasible if electrostatic force is greater than gravitation force mg.
In this paper are carried out the analyses of electrostatic levitation forces for oblate and prolate dielectric spheroids, that extends the analysis for the ball [13].The analytical solutions are derived by means of variable separation method for electric potential.The distributions of electric field are given by Legendre functions of first and second kinds.The Legendre functions are given by means of hypergeometric series [3].
The analytical solutions are clearly revealing the influence of spheroids parameters (the most interesting is anisotropy) on the field distributions and electrostatic forces.The anisotropy can appear as a result of crystalline feature, rapid solidification, layering, ball crushing or 3D printing of produced spheroid.The permittivity may be adjusted due to saturation or nonlinearity for particular work point, too.The analytical solution presented could be regarded as a default solution.Moreover, the analytical solution is a benchmark tests for numerical algorithms.It may constitute start points for multistep design and can be also a part of hybrid, i.e. analytical-numerical algorithms [2,5].
The approaches present analytical solutions of electrostatic levitations for diagonally anisotropic dielectric spheroids.The approaches cover analysis of electric field distribution inside and outside the spheroids and multiple levitation force evaluations.The levitations forces are evaluated by the Maxwell stress tensor generalized method, material force, coenergy and equivalent dipole formula (valid only for isotopic spheroids).
The novelty of the approaches lies in the fact that are analytically described: • The electric field distributions inside and outside of both oblate and prolate anisotropic dielectric spheroids, and • The force arising at the surface of anisotropic dielectric both oblate and prolate spheroids by four methods.
The three semi-axes of spheroid are denoted as follows for either oblate (R, R, h < R) or prolate (R, R, h > R) spheroids, respectively.Both spheroids are geometrically different thus appropriate coordinate systems used are different, too.Field distributions for oblate and prolate spheroids are given by different series and will be separately analysed.

Electrostatic levitation of oblate spheroid
Let imposed static electric field be along to z-axis oriented upwards-Fig.1.In the field is set non-charged, nonconductive, homogeneous and diagonally anisotropic dielectric oblate spheroid.
The imposed electric field strength (far from the spheroid) parallel component to axis z is as given below where E 0 and F k are given constants (z 0 for the centre of spheroid).For N > 2 the imposed field gradient is spatially variable.If N 2, thus the field gradient is constant.
In case E 0 0 and F k 0 (k 1, 2, 3, …) the field is symmetrical above and below the spheroid, thus force does not appear.In such a case, the presented below solutions constitute solutions for the academic problem of spheroid in a uniform electrostatic field.The most appropriate orthogonal coordinate system for oblate spheroid electric field analysis is oblate spheroidal coordinates (Fig. 2).
For c e > 0 the oblate spheroidal coordinates η, θ are as follows where d 1 (ρ − c e ) 2 + z 2 , d 2 (ρ + c e ) 2 + z 2 are distances of the point P(ρ, z) from the two focuses (+ c e , 0) and (c e , 0), respectively.The coordinate η max acosh(R/c e ) describes the ellipse (d 1 + d 2 2R) at which is located a certain point P(η max ,θ ), and θ is the angle between hyperbole asymptote and axis z (Fig. 2).These coordinates must be always supplemented by the linear eccentricity c e .
Inversely, cylindrical coordinates are as follows and z(η, θ ) c e sh(η) cos(θ ). ( The assumption c e > 0 implicates that analytical solutions for neither oblate nor prolate spheroids can be applied for ball (c e ≡ 0).However, the solutions for both spheroids, in limit c e → 0, converge to the solution of the ball [13].
Axially symmetrical electrostatic field problem describes electric scalar potential that is independent from longitudinal angle ϕ Hence, according to relation all electric field components are determined, e.g.
For diagonal (normal, principal axes) anisotropy of medium, the electric displacement field in oblate spheroidal coordinates is as follows Gauss's law for divergence-free electric field is as follows From ( 10) results partial differential equation as follows For variable separation in the form of is satisfied The Legendre polynomials P n (cos(θ )), for integer n > 0, are solutions (θ ) of the Lagrange equation in the form of The second linearly independent solution constitutes Legendre functions denoted by Q n (cos(θ )) ( [3], pp.958-960).Subsequently, it can be written equation where While the sign minus in the front of square root in ( 16) is being rejected, then ν is always positive for natural n.
The first solutions of (15) take the form of and are built of associated Legendre functions P ν μ (z) ( [3] Eq. 8.820.1).P ν μ (z) are convergent for |1 -z|< 2. While μ 0-electrostatic problem described by (13), and it is denoting for simplicity P ν 0 (z) P ν (z).The second linearly independent solutions are the functions Electric scalar potential outside the oblate spheroid (isotropic region, ν n) takes the form of and inside oblate spheroid is as follows where d n 0 due to singularities of Q ν (z) and its derivatives ( [3], p. 975).
The real parts of complex solutions are approved as physical results and further are applied for force calculations.
Basing on electric scalar potential the electric field strength (7) components can be easily derived.Exemplary, inside the spheroid and The electric field strength distribution must be described by (1) far from the spheroid.Thus, the constants a n result directly from the imposed field (1) and depend on polynomial coefficients as follows and for n 2, …, N where constants k n 1/3, 2/15, 2/35, …, g max j•sh(η max ), and S n ε ηη ε out P ν (g max )P n (g max ) − P ν (g max )P n (g max ) The electric field must obey two boundary conditions, i.e. continuities of both tangential electric strength E θ and normal electric displacement D η at the spheroid surface η η max (outside the spheroid ε θθ ε ηη ε out ).The continuity conditions constitute two equations for unknown constants b n , c n .Hence, the constants b n , c n (n 1, 2, …, N) are as follows The electric scalar potential leads to electric field strength by means of gradient formula (7) and electric displacement field by means of (9), subsequently.The electric field strength and electric displacement determine the force acting on the spheroid, i.e. electrostatic levitation force.Levitation force physically acts at the boundary of the spheroid, because only there appears permittivity step change in electrostatic field.
For oblate spheroidal orthogonal coordinates electric force density components are as follows and where ("Appendix I") It should be underlined that the above equations are derived from the Maxwell electromagnetic field equations and Lorentz force formula ( [10,12], Table 1).
Force density along z-axis is given as below where α n is angle between z-axis and normal to the oblate spheroid surface.It should be underline that the angle α n is not equal to the coordinate θ (Fig. 2).The angle α n is described by the trigonometric functions as follows .
Equivalently, the force density formula can be derived by means of coordinate transformation as it is given below (33) .
Hence, the force density along z axis (31) describes the following relation where vectors are built of the Maxwell stress tensor components [1,5,11,12] as follows The main Eq.(34) leads to the Maxwell stress tensor generalized method where only the Maxwell stress tensor is integrated over the spheroid surface at η η max+ as follows The additional term of the Maxwell stress tensor generalized method F z is a volume integral.For oblate spheroidal coordinates value of F z results directly from (34) and takes the following extended forms given by ( 38) or (39) or equivalently Physically, the electrostatic levitation force, i.e. material force F Nz is acting only at surface of the spheroid where permittivity changes in step way.The material force density is equal to zero inside the homogeneous spheroid.Mathematically, the material force describes inhomogeneity component [5,10] as follows where repeated indices are implicitly summed over.For diagonally (normally) anisotropy of spheroid permittivity (9) force can be presented as follows where appear only continuous field components D η and E θ at the boundary η η max .
In order to confirm the value of the electrostatic force coenergy method [1,5] is applied as follows hence Formulas ( 41) and ( 42) are mathematically equivalent.
For isotropic spheroid, the force can be also calculated by means of equivalent electric dipole model [1,12,13] in constant gradient field (1), i.e. for N 2. Namely, for a effective dipole moment as follows the force acting on the dipole given by the formula is equal to electrostatic force.The oblate spheroid shape is similar to a ball for h → R (c e → 0, η max → ∞).For h → 0 (c e → R, η max → 0) spheroid is similar to a thin disc (is coin-like).In this second case the force (41) decreases to zero almost linearly for h → 0, and tangent line has got a slope of The tangent for thin spheroid is drawn in Fig. 5 (and Fig. 11).
The forces evaluated for electrostatic field by means of the Maxwell stress tensor generalized method F z + F z , material force F Nz (inhomogeneity component), coenergy F Cz and equivalent dipole F Dz (only for isotropic dielectric spheroid) are equal one another In Figs. 3 and 4, the electrostatic levitation forces are presented vs. relative axis permittivity change + /-20% of permittivity for the isotropic case.The greater axis permittivity the greater force.The increase is greater for latitudinal permittivity ε θθ rises than for hyperbolic permittivity ε ηη rises.It should be also pointed out that the residual component F z (the Maxwell stress tensor generalized method volume integral [12,13]) vanishes for isotropic dielectric spheroid (central abscissa point of each curve).
Figure 5 presents the electrostatic force vs. oblate spheroid height h (0, R].The tangent for thin disc h → 0 (coin-like) is drawn due to (46).The star * in the right-upper corner of Fig. 5 is incorporated basing on analytical solution obtained independently in [13] for the ball, i.e. h ≡ R (analogously, the star * is set in Fig. 9).It can be seen that the solution for spheroid at limit h → R is convergent to that for the ball h ≡ R.
The force evaluated by both the Maxwell stress tensor generalized and material force methods always leads to equal results Fig. 5 Electrostatic levitation force for anisotropic dielectric oblate evaluated by Maxwell stress tensor generalized method (points), material force (dash-dot line) and residual component F z (line near abscissa) versus height h For checking, the global average relative error is defined as follows Exemplary, for forces in Fig. 5 (n 20 intervals per curve) the average relative error is equal to 4.414E−003.
If F zi were neglected for anisotropic region, thus maximal error would be equal to 5.1E−002 (Fig. 5).

Electrostatic levitation of prolate spheroid
Let us consider consecutively prolate spheroid at the same conditions of electric field excitation (Fig. 6).
The partial differential Eq. ( 11) for prolate spheroid takes slightly different form, due to different Lame coefficients, than for oblate spheroid problem.Subsequently, the separation Eq. ( 13) leads to the following relations The first solution for H() function in prolate spheroidal coordinates is as follows The second linearly independent solutions constitute Legendre functions of second kind Q ν μ (ch(η)) (for μ 0 are shortly denoted as Q ν (ch(η)) [3]).The solutions for angular function (θ ) are Legendre polynomials P n (cos(θ )).Electric scalar potential outside the prolate spheroid (isotropic region) takes the form of and inside prolate spheroid is as follows below where d n 0 due to irregularity of function Q ν () Q ν 0 ().Basing on solutions of electric scalar potential the electric field strength is derived by gradient operator for prolate spheroidal coordinates.Subsequently, the electric field force can be evaluated by analogous formulas as those for the oblate spheroid.The mathematical differences in oblate and prolate formulas lay in different relations for Lame coefficients and normal angle β n (Fig. 7).Exemplary, for prolate spheroid the additional term for the Maxwell stress tensor generalized method takes the form as follows The relations for material force, stress tensor surface integral, coenergy method formula and equivalent dipole moment are analogous.
In Figs. 8 and 9, electrostatic forces for prolate spheroid vs. relative axis permittivity change ± 20% of permittivity for the isotropic case.Force increase is greater for hyperbolic permittivity ε ηη changes than for latitudinal permittivity ε θθ changes.In Fig. 10 the force vs. prolate spheroid height is shown.

Oblate and prolate spheroids
Electrostatic forces, that may cause electrostatic levitation in gravitation field, appears at the boundary of spheroids.The forces are of material nature and are appearing when the inside permittivity is different from the outside permittivity.
Among others, the forces depend on spheroid axes length R and h.The relative differences between these length enable to present the geometry influence on forces for each spheroid height h.Relative axes difference (a-the greater semi-axis R, b-the smaller semi-axis h) for the oblate spheroid is defined as given below and converges to minus unity for infinitely thin disc (coinlike).
For prolate spheroid relative axes difference (a-the greater semi-axis h, b-the smaller semi-axis R) is defined as follows and converges to plus unity for the infinitely high prolate spheroid (stick-like).Relative axes difference is equal to zero for the ball.While relative axes difference is converging to zero (h → R), the shape of each spheroid is converging to the shape of ball (h ≡ R).
The definitions (55) and (56) enable to present forces for any spheroid h 0 ÷ ∞ at one common curve shown in Fig. 11.
Obviously, the curve presented in Fig. 11 does not include points h 0 and h ∞.
For h → 0 the slope of force curve is equal to tangent given by the limit relation (46).
Figure 11 proves that the forces for both oblate (h → R-) and prolate (h → R +) spheroids are converging to the force of the ball (h ≡ R).The electrostatic force for the ball is the greatest one (for given imposed electrostatic field).
It should be pointed out that the analytical solutions given by series (especially by functions Q ν μ () [3], pp.958-960) are weakly converging for prolate spheroid at limits both h → R + and h → ∞.

Conclusions
The analytical solution of oblate and prolate spheroids in electrostatic field has been presented.The presence of electrostatic field may result in levitation of dielectric spheroids either oblate or prolate.The considered spheroids are nonconductive and diagonal (normal) anisotropic dielectric in oblate and prolate spheroidal coordinates The electric fields distributions for oblate and prolate spheroids are obtained by the variable separation method in adequate co-ordinate system.
It should be pointed out that the differential equations for oblate and prolate spheroids are different due to different formulas for Lame coefficients.Hence, mathematical formulas for field components and forces are slightly different.As a consequence, the forces for oblate and prolate spheroids at the same ratio of axes length are different (Fig. 11).
The correctness of the methodology and calculations are confirmed among others by multiple force calculations for both spheroids.
Particularly, four methods for each spheroid are applied.There are the Maxwell stress tensor generalized, material force, coenergy and equivalent dipole methods.The first three methods lead always to the same values of forces.The fourth method, i.e. equivalent dipole method, is valid only for isotropic dielectric spheroid.
The analytical solutions obtained for dielectric spheroids are valid over a large range of parameters, e.g.permittivities, height, radius and imposed field parameters that is not so easy to obtain by numerical analysis.
The presented solutions for oblate and prolate dielectric spheroids can be applied for rapid design, testing numerical algorithms and prototyping.
The presented solution for electrostatic problem can be applied also for no-current magnetostatic problem.The introducing relations ( 6) and ( 7) can be applied directly for magnetic scalar potential, too.However, such a solution for scalar magnetic potential cannot be applied when current appears [13].
a n , b n , c n , d n Constants for electric scalar potential series e E u D u /2 Electric field energy density bilinear form; u η, θ , ϕ e e η + e θ + e ϕ Electric field energy density bilinear form E 0

Fig. 1
Fig. 1 Anisotropic dielectric and nonconductive oblate spheroid in imposed electrostatic field

Fig. 6
Fig. 6 Anisotropic dielectric and nonconductive prolate spheroid in imposed electrostatic field

Table 1
Forces densities in electromagnetic field