Adding omega (angular velocity) helper function

src/lagrangian/intermediate/parcels/Templates/InteractingKinematicParcel/InteractingKinematicParcel.C

@@ -241,6 +241,8 @@ bool Foam::InteractingKinematicParcel<ParcelType>::move(TrackData& td)
 
             p.U() += 0.5*deltaT*p.f()/p.mass();
 
+            pi_ += 0.5*deltaT*tau_;
+
             break;
         }
 

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

@@ -34,6 +34,9 @@ Description
     - turbulent dispersion
     - wall interactions
 
+    Includes many-body collisions and rotational motion.  Only
+    spherical particles currently supported.
+
 SourceFiles
     InteractingKinematicParcelI.H
     InteractingKinematicParcel.C
@@ -432,6 +435,9 @@ public:
             //- Particle moment of inertia around diameter axis
             inline scalar momentOfInertia() const;
 
+            //- Particle angular velocity
+            inline vector omega() const;
+
             //- Particle volume
             inline scalar volume() const;
 

src/lagrangian/intermediate/parcels/Templates/InteractingKinematicParcel/InteractingKinematicParcelI.H

@@ -437,6 +437,14 @@ Foam::InteractingKinematicParcel<ParcelType>::momentOfInertia() const
 }
 
 
+template <class ParcelType>
+inline Foam::vector
+Foam::InteractingKinematicParcel<ParcelType>::omega() const
+{
+    return pi_/momentOfInertia();
+}
+
+
 template <class ParcelType>
 inline Foam::scalar Foam::InteractingKinematicParcel<ParcelType>::volume() const
 {