From b93156f14279a6340a90cff0fb87a0167bc9b2e5 Mon Sep 17 00:00:00 2001
From: graham <g.macpherson@opencfd.co.uk>
Date: Mon, 7 Sep 2009 11:45:13 +0100
Subject: [PATCH] Adding momentum and linear kinetic energy monitoring
 functions.

---
 .../InteractingKinematicCloud.C               | 14 +++++++-
 .../InteractingKinematicCloud.H               |  6 ++++
 .../InteractingKinematicCloudI.H              | 34 +++++++++++++++++++
 3 files changed, 53 insertions(+), 1 deletion(-)

diff --git a/src/lagrangian/intermediate/clouds/Templates/InteractingKinematicCloud/InteractingKinematicCloud.C b/src/lagrangian/intermediate/clouds/Templates/InteractingKinematicCloud/InteractingKinematicCloud.C
index 7f8cdecc7c..2ebfb56f80 100644
--- a/src/lagrangian/intermediate/clouds/Templates/InteractingKinematicCloud/InteractingKinematicCloud.C
+++ b/src/lagrangian/intermediate/clouds/Templates/InteractingKinematicCloud/InteractingKinematicCloud.C
@@ -259,6 +259,12 @@ void Foam::InteractingKinematicCloud<ParcelType>::evolve()
 template<class ParcelType>
 void Foam::InteractingKinematicCloud<ParcelType>::info() const
 {
+    vector linearMomentum = linearMomentumOfSystem();
+    reduce(linearMomentum, sumOp<vector>());
+
+    scalar linearKineticEnergy = linearKineticEnergyOfSystem();
+    reduce(linearKineticEnergy, sumOp<scalar>());
+
     Info<< "Cloud: " << this->name() << nl
         << "    Total number of parcels added   = "
         << returnReduce(this->injection().parcelsAddedTotal(), sumOp<label>())
@@ -269,7 +275,13 @@ void Foam::InteractingKinematicCloud<ParcelType>::info() const
         << "    Current number of parcels       = "
         << returnReduce(this->size(), sumOp<label>()) << nl
         << "    Current mass in system          = "
-        << returnReduce(massInSystem(), sumOp<scalar>()) << nl;
+        << returnReduce(massInSystem(), sumOp<scalar>()) << nl
+        << "    Linear momentum                 = "
+        << linearMomentum << nl
+        << "   |Linear momentum|                = "
+        << mag(linearMomentum) << nl
+        << "    Linear kinetic energy           = "
+        << linearKineticEnergy << nl;
 }
 
 
diff --git a/src/lagrangian/intermediate/clouds/Templates/InteractingKinematicCloud/InteractingKinematicCloud.H b/src/lagrangian/intermediate/clouds/Templates/InteractingKinematicCloud/InteractingKinematicCloud.H
index 8ea351540f..e5b3a2e33d 100644
--- a/src/lagrangian/intermediate/clouds/Templates/InteractingKinematicCloud/InteractingKinematicCloud.H
+++ b/src/lagrangian/intermediate/clouds/Templates/InteractingKinematicCloud/InteractingKinematicCloud.H
@@ -361,6 +361,12 @@ public:
             //- Total mass in system
             inline scalar massInSystem() const;
 
+            //- Total linear momentum of the system
+            inline vector linearMomentumOfSystem() const;
+
+            //- Total linear kinetic energy in the system
+            inline scalar linearKineticEnergyOfSystem() const;
+
             //- Print cloud information
             void info() const;
 
diff --git a/src/lagrangian/intermediate/clouds/Templates/InteractingKinematicCloud/InteractingKinematicCloudI.H b/src/lagrangian/intermediate/clouds/Templates/InteractingKinematicCloud/InteractingKinematicCloudI.H
index 4cc79f424c..9b4af061d9 100644
--- a/src/lagrangian/intermediate/clouds/Templates/InteractingKinematicCloud/InteractingKinematicCloudI.H
+++ b/src/lagrangian/intermediate/clouds/Templates/InteractingKinematicCloud/InteractingKinematicCloudI.H
@@ -221,6 +221,40 @@ Foam::InteractingKinematicCloud<ParcelType>::massInSystem() const
 }
 
 
+template<class ParcelType>
+inline Foam::vector
+Foam::InteractingKinematicCloud<ParcelType>::linearMomentumOfSystem() const
+{
+    vector linearMomentum(vector::zero);
+
+    forAllConstIter(typename InteractingKinematicCloud<ParcelType>, *this, iter)
+    {
+        const ParcelType& p = iter();
+
+        linearMomentum += p.mass()*p.U();
+    }
+
+    return linearMomentum;
+}
+
+
+template<class ParcelType>
+inline Foam::scalar
+Foam::InteractingKinematicCloud<ParcelType>::linearKineticEnergyOfSystem() const
+{
+    scalar linearKineticEnergy = 0.0;
+
+    forAllConstIter(typename InteractingKinematicCloud<ParcelType>, *this, iter)
+    {
+        const ParcelType& p = iter();
+
+        linearKineticEnergy += 0.5*p.mass()*(p.U() & p.U());
+    }
+
+    return linearKineticEnergy;
+}
+
+
 template<class ParcelType>
 inline Foam::Random& Foam::InteractingKinematicCloud<ParcelType>::rndGen()
 {