ENH: Added calculation for Eotvos number

src/lagrangian/intermediate/parcels/Templates/KinematicParcel/KinematicParcel.H

@@ -572,6 +572,14 @@ public:
                 const scalar sigma      // particle surface tension
             ) const;
 
+            //- Eotvos number
+            inline scalar Eo
+            (
+                const vector& a,        // acceleration
+                const scalar d,         // particle diameter
+                const scalar sigma      // particle surface tension
+            ) const;
+
 
         // Main calculation loop
 

src/lagrangian/intermediate/parcels/Templates/KinematicParcel/KinematicParcelI.H

@@ -567,4 +567,17 @@ inline Foam::scalar Foam::KinematicParcel<ParcelType>::We
 }
 
 
+template<class ParcelType>
+inline Foam::scalar Foam::KinematicParcel<ParcelType>::Eo
+(
+    const vector& a,
+    const scalar d,
+    const scalar sigma
+) const
+{
+    vector dir = U_/(mag(U_) + ROOTVSMALL);
+    return mag(a & dir)*(rho_ - rhoc_)*sqr(d)/(sigma + ROOTVSMALL);
+}
+
+
 // ************************************************************************* //