From 53f0bf3f0d16bcad9ffabee1d292e662a751cc81 Mon Sep 17 00:00:00 2001
From: henry <henry>
Date: Mon, 25 Jan 2010 12:18:54 +0000
Subject: [PATCH] App settlingFoam: Added under-relaxation to k-equation.

---
 applications/solvers/multiphase/settlingFoam/UEqn.H         | 2 --
 applications/solvers/multiphase/settlingFoam/createFields.H | 3 ++-
 applications/solvers/multiphase/settlingFoam/kEpsilon.H     | 7 +++++--
 tutorials/multiphase/settlingFoam/ras/tank3D/0/U            | 3 ++-
 tutorials/multiphase/settlingFoam/ras/tank3D/0/alpha        | 3 ++-
 tutorials/multiphase/settlingFoam/ras/tank3D/0/epsilon      | 3 ++-
 tutorials/multiphase/settlingFoam/ras/tank3D/0/k            | 3 ++-
 7 files changed, 15 insertions(+), 9 deletions(-)

diff --git a/applications/solvers/multiphase/settlingFoam/UEqn.H b/applications/solvers/multiphase/settlingFoam/UEqn.H
index e903434a..390eeb05 100644
--- a/applications/solvers/multiphase/settlingFoam/UEqn.H
+++ b/applications/solvers/multiphase/settlingFoam/UEqn.H
@@ -22,8 +22,6 @@
           ==
             fvc::reconstruct
             (
-                //fvc::interpolate(rho)*(g & mesh.Sf())
-                //- fvc::snGrad(p)*mesh.magSf()
                 (- ghf*fvc::snGrad(rho) - fvc::snGrad(pmh))*mesh.magSf()
             )
         );
diff --git a/applications/solvers/multiphase/settlingFoam/createFields.H b/applications/solvers/multiphase/settlingFoam/createFields.H
index 54205449..7485fe04 100644
--- a/applications/solvers/multiphase/settlingFoam/createFields.H
+++ b/applications/solvers/multiphase/settlingFoam/createFields.H
@@ -1,4 +1,4 @@
-    Info<< "Reading field p\n" << endl;
+    Info<< "Reading field pmh\n" << endl;
     volScalarField pmh
     (
         IOobject
@@ -379,6 +379,7 @@
         mut + mul
     );
 
+
     Info<< "Calculating field (g.h)f\n" << endl;
     volScalarField gh("gh", g & mesh.C());
     surfaceScalarField ghf = surfaceScalarField("gh", g & mesh.Cf());
diff --git a/applications/solvers/multiphase/settlingFoam/kEpsilon.H b/applications/solvers/multiphase/settlingFoam/kEpsilon.H
index 742fcde7..c95a5408 100644
--- a/applications/solvers/multiphase/settlingFoam/kEpsilon.H
+++ b/applications/solvers/multiphase/settlingFoam/kEpsilon.H
@@ -45,7 +45,7 @@ if(turbulence)
 
     // Turbulent kinetic energy equation
 
-    solve
+    fvScalarMatrix kEqn
     (
         fvm::ddt(rho, k)
       + fvm::div(phi, k)
@@ -60,11 +60,14 @@ if(turbulence)
       - fvm::Sp(rho*epsilon/k, k)
     );
 
+    kEqn.relax();
+    kEqn.solve();
+
     bound(k, k0);
 
 
     //- Re-calculate viscosity
-    mut = rho*Cmu*sqr(k)/(epsilon + epsilon0);
+    mut = rho*Cmu*sqr(k)/epsilon;
 
 #   include "wallViscosity.H"
 }
diff --git a/tutorials/multiphase/settlingFoam/ras/tank3D/0/U b/tutorials/multiphase/settlingFoam/ras/tank3D/0/U
index bac98acc..3e6dfd4b 100644
--- a/tutorials/multiphase/settlingFoam/ras/tank3D/0/U
+++ b/tutorials/multiphase/settlingFoam/ras/tank3D/0/U
@@ -136,7 +136,8 @@ boundaryField
 
     OUTL15
     {
-        type            zeroGradient;
+        type            inletOutlet;
+        inletValue      uniform (0 0 0);
     }
 }
 
diff --git a/tutorials/multiphase/settlingFoam/ras/tank3D/0/alpha b/tutorials/multiphase/settlingFoam/ras/tank3D/0/alpha
index 18e13658..44ad81c6 100644
--- a/tutorials/multiphase/settlingFoam/ras/tank3D/0/alpha
+++ b/tutorials/multiphase/settlingFoam/ras/tank3D/0/alpha
@@ -118,7 +118,8 @@ boundaryField
 
     OUTL15
     {
-        type            zeroGradient;
+        type            inletOutlet;
+        inletValue      uniform 0;
     }
 }
 
diff --git a/tutorials/multiphase/settlingFoam/ras/tank3D/0/epsilon b/tutorials/multiphase/settlingFoam/ras/tank3D/0/epsilon
index 634512f9..ec50c5c5 100644
--- a/tutorials/multiphase/settlingFoam/ras/tank3D/0/epsilon
+++ b/tutorials/multiphase/settlingFoam/ras/tank3D/0/epsilon
@@ -118,7 +118,8 @@ boundaryField
 
     OUTL15
     {
-        type            zeroGradient;
+        type            inletOutlet;
+        inletValue      uniform 1.973e-07;
     }
 }
 
diff --git a/tutorials/multiphase/settlingFoam/ras/tank3D/0/k b/tutorials/multiphase/settlingFoam/ras/tank3D/0/k
index 5bc3ea2d..304c689a 100644
--- a/tutorials/multiphase/settlingFoam/ras/tank3D/0/k
+++ b/tutorials/multiphase/settlingFoam/ras/tank3D/0/k
@@ -118,7 +118,8 @@ boundaryField
 
     OUTL15
     {
-        type            zeroGradient;
+        type            inletOutlet;
+        inletValue      uniform 0.000259;
     }
 }
 
-- 
2.11.4.GIT