From 32870ea73f7f12420d9a079dfd4cdc77b2c1b8cf Mon Sep 17 00:00:00 2001
From: Thomas Lichtenegger <thomas.lichtenegger@jku.at>
Date: Wed, 13 Jul 2016 18:06:02 +0200
Subject: [PATCH] Added testing functions. Will be removed again.

---
 .../cfdemSolverRhoPimple.C                    |  4 +-
 .../cfdemSolverRhoPimple/createFields.H       | 57 ++++++++++++++++++-
 2 files changed, 58 insertions(+), 3 deletions(-)

diff --git a/applications/solvers/cfdemSolverRhoPimple/cfdemSolverRhoPimple.C b/applications/solvers/cfdemSolverRhoPimple/cfdemSolverRhoPimple.C
index 64a45dc0..3461c64c 100644
--- a/applications/solvers/cfdemSolverRhoPimple/cfdemSolverRhoPimple.C
+++ b/applications/solvers/cfdemSolverRhoPimple/cfdemSolverRhoPimple.C
@@ -121,7 +121,9 @@ int main(int argc, char *argv[])
                 #include "pEqn.H"
 	      	rhoeps=rho*voidfraction;
             }
-
+//divphi = fvc::div(phi);
+//divphi2 = fvc::div(linearInterpolate(rho*U*voidfraction) & mesh.Sf());
+	gP = fvc::grad(p);
             if (pimple.turbCorr())
             {
                 turbulence->correct();
diff --git a/applications/solvers/cfdemSolverRhoPimple/createFields.H b/applications/solvers/cfdemSolverRhoPimple/createFields.H
index 1c69db2b..6f806fa1 100644
--- a/applications/solvers/cfdemSolverRhoPimple/createFields.H
+++ b/applications/solvers/cfdemSolverRhoPimple/createFields.H
@@ -194,7 +194,60 @@ Info<< "Reading thermophysical properties\n" << endl;
 
 //===============================
 
+volScalarField divphi
+    (
+        IOobject
+        (
+            "divphi",
+            runTime.timeName(),
+            mesh,
+            IOobject::NO_READ,
+            IOobject::AUTO_WRITE
+        ),
+        mesh,
+        dimensionedScalar("0", dimensionSet(1, -3, -1, 0, 0), 0.0)
+    );
+    
+    volScalarField divphi2
+    (
+        IOobject
+        (
+            "divphi2",
+            runTime.timeName(),
+            mesh,
+            IOobject::NO_READ,
+            IOobject::AUTO_WRITE
+        ),
+        mesh,
+        dimensionedScalar("0", dimensionSet(1, -3, -1, 0, 0), 0.0)
+    );
+
+volVectorField gP
+    (
+        IOobject
+        (
+            "gP",
+            runTime.timeName(),
+            mesh,
+            IOobject::NO_READ,
+            IOobject::AUTO_WRITE
+        ),
+        mesh,
+        dimensionedVector("0", dimensionSet(1, -2, -2, 0, 0), vector::zero)
+    );
+//     
+// volVectorField gP2
+//     (
+//         IOobject
+//         (
+//             "gP2",
+//             runTime.timeName(),
+//             mesh,
+//             IOobject::NO_READ,
+//             IOobject::AUTO_WRITE
+//         ),
+//         mesh,
+//         dimensionedVector("0", dimensionSet(1, -2, -2, 0, 0), vector::zero)
+//     );
 
 
- //   singlePhaseTransportModel laminarTransport(U, phi);
-