From 624b228f2d1131f1ec15523d29a01ff5ca38975a Mon Sep 17 00:00:00 2001
From: Henry <Henry>
Date: Tue, 28 Feb 2012 15:55:07 +0000
Subject: [PATCH] twoPhaseEulerFoam: Add a residualSlip velocity magnitude for
 the drag

---
 .../multiphase/twoPhaseEulerFoam/createFields.H       | 11 ++++++++++-
 .../multiphase/twoPhaseEulerFoam/liftDragCoeffs.H     |  2 +-
 2 files changed, 11 insertions(+), 2 deletions(-)

diff --git a/applications/solvers/multiphase/twoPhaseEulerFoam/createFields.H b/applications/solvers/multiphase/twoPhaseEulerFoam/createFields.H
index 95395d97d4..137c63827b 100644
--- a/applications/solvers/multiphase/twoPhaseEulerFoam/createFields.H
+++ b/applications/solvers/multiphase/twoPhaseEulerFoam/createFields.H
@@ -112,7 +112,9 @@
         (
             "phi",
             runTime.timeName(),
-            mesh
+            mesh,
+            IOobject::NO_READ,
+            IOobject::AUTO_WRITE
         ),
         fvc::interpolate(alpha)*phia + fvc::interpolate(beta)*phib
     );
@@ -194,6 +196,13 @@
         }
     }
 
+    dimensionedScalar residualSlip
+    (
+        "residualSlip",
+        dimVelocity,
+        interfacialProperties.lookup("residualSlip")
+    );
+
     Info << "dragPhase is " << dragPhase << endl;
     kineticTheoryModel kineticTheory
     (
diff --git a/applications/solvers/multiphase/twoPhaseEulerFoam/liftDragCoeffs.H b/applications/solvers/multiphase/twoPhaseEulerFoam/liftDragCoeffs.H
index b614a4dfbe..27519eaf2e 100644
--- a/applications/solvers/multiphase/twoPhaseEulerFoam/liftDragCoeffs.H
+++ b/applications/solvers/multiphase/twoPhaseEulerFoam/liftDragCoeffs.H
@@ -1,5 +1,5 @@
     volVectorField Ur(Ua - Ub);
-    volScalarField magUr(mag(Ur));
+    volScalarField magUr(mag(Ur) + residualSlip);
 
     volScalarField Ka(draga->K(magUr));
     volScalarField K(Ka);