From 598f705e35ba7b39d6dd537d9ce10f49b3528eaf Mon Sep 17 00:00:00 2001
From: graham <g.macpherson@opencfd.co.uk>
Date: Wed, 9 Sep 2009 17:09:43 +0100
Subject: [PATCH] Adding omega (angular velocity) helper function

---
 .../InteractingKinematicParcel.C                          | 2 ++
 .../InteractingKinematicParcel.H                          | 6 ++++++
 .../InteractingKinematicParcelI.H                         | 8 ++++++++
 3 files changed, 16 insertions(+)

diff --git a/src/lagrangian/intermediate/parcels/Templates/InteractingKinematicParcel/InteractingKinematicParcel.C b/src/lagrangian/intermediate/parcels/Templates/InteractingKinematicParcel/InteractingKinematicParcel.C
index 9f3394c28d..3773321506 100644
--- a/src/lagrangian/intermediate/parcels/Templates/InteractingKinematicParcel/InteractingKinematicParcel.C
+++ b/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;
         }
 
diff --git a/src/lagrangian/intermediate/parcels/Templates/InteractingKinematicParcel/InteractingKinematicParcel.H b/src/lagrangian/intermediate/parcels/Templates/InteractingKinematicParcel/InteractingKinematicParcel.H
index 475c2db163..48ae97ed79 100644
--- a/src/lagrangian/intermediate/parcels/Templates/InteractingKinematicParcel/InteractingKinematicParcel.H
+++ b/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;
 
diff --git a/src/lagrangian/intermediate/parcels/Templates/InteractingKinematicParcel/InteractingKinematicParcelI.H b/src/lagrangian/intermediate/parcels/Templates/InteractingKinematicParcel/InteractingKinematicParcelI.H
index 0f95bceb25..ee2df1df80 100644
--- a/src/lagrangian/intermediate/parcels/Templates/InteractingKinematicParcel/InteractingKinematicParcelI.H
+++ b/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
 {