From c8078e73e1de8fe67c0af342552a0f1c5078ec62 Mon Sep 17 00:00:00 2001
From: Henry <Henry>
Date: Fri, 18 Mar 2011 10:25:25 +0000
Subject: [PATCH 1/2] engineFoam: Corrected mesh construction and tutorial

---
 src/engine/include/createEngineMesh.H                     | 3 ++-
 tutorials/combustion/engineFoam/kivaTest/system/fvSchemes | 2 +-
 2 files changed, 3 insertions(+), 2 deletions(-)

diff --git a/src/engine/include/createEngineMesh.H b/src/engine/include/createEngineMesh.H
index 88a9e06686..263ebed5ff 100644
--- a/src/engine/include/createEngineMesh.H
+++ b/src/engine/include/createEngineMesh.H
@@ -9,7 +9,8 @@
             (
                 engineMesh::defaultRegion,
                 runTime.timeName(),
-                runTime
+                runTime,
+                Foam::IOobject::MUST_READ
             )
         )
     );
diff --git a/tutorials/combustion/engineFoam/kivaTest/system/fvSchemes b/tutorials/combustion/engineFoam/kivaTest/system/fvSchemes
index 7d9e15bae0..0a28cdbe58 100644
--- a/tutorials/combustion/engineFoam/kivaTest/system/fvSchemes
+++ b/tutorials/combustion/engineFoam/kivaTest/system/fvSchemes
@@ -60,7 +60,7 @@ laplacianSchemes
     laplacian(DepsilonEff,epsilon) Gauss linear limited 0.5;
     laplacian(DREff,R) Gauss linear limited 0.5;
     laplacian((rho*(1|A(U))),p) Gauss linear limited 0.5;
-    laplacian(muEff,b) Gauss linear limited 0.5;
+    laplacian(alphaEff,b) Gauss linear limited 0.5;
     laplacian(muEff,ft) Gauss linear limited 0.5;
     laplacian(alphaEff,h) Gauss linear limited 0.5;
     laplacian(alphaEff,hu) Gauss linear limited 0.5;

From 172c026e67140e44fcfe0be4fa84971ecff863fb Mon Sep 17 00:00:00 2001
From: Henry <Henry>
Date: Fri, 18 Mar 2011 10:25:56 +0000
Subject: [PATCH 2/2] porousSimpleFoam: Completed relocation

---
 .../incompressible/simpleFoam/Allwmake        |  8 ++++++
 .../porousSimpleFoam/createPorousZones.H      | 26 +++++++++++++++++++
 2 files changed, 34 insertions(+)
 create mode 100755 applications/solvers/incompressible/simpleFoam/Allwmake
 create mode 100644 applications/solvers/incompressible/simpleFoam/porousSimpleFoam/createPorousZones.H

diff --git a/applications/solvers/incompressible/simpleFoam/Allwmake b/applications/solvers/incompressible/simpleFoam/Allwmake
new file mode 100755
index 0000000000..577d11018f
--- /dev/null
+++ b/applications/solvers/incompressible/simpleFoam/Allwmake
@@ -0,0 +1,8 @@
+#!/bin/sh
+cd ${0%/*} || exit 1    # run from this directory
+set -x
+
+wmake
+wmake porousSimpleFoam
+
+# ----------------------------------------------------------------- end-of-file
diff --git a/applications/solvers/incompressible/simpleFoam/porousSimpleFoam/createPorousZones.H b/applications/solvers/incompressible/simpleFoam/porousSimpleFoam/createPorousZones.H
new file mode 100644
index 0000000000..d14adbb493
--- /dev/null
+++ b/applications/solvers/incompressible/simpleFoam/porousSimpleFoam/createPorousZones.H
@@ -0,0 +1,26 @@
+    porousZones pZones(mesh);
+    Switch pressureImplicitPorosity(false);
+
+    // nUCorrectors used for pressureImplicitPorosity
+    int nUCorr = 0;
+    if (pZones.size())
+    {
+        // nUCorrectors for pressureImplicitPorosity
+        if (mesh.solutionDict().subDict("SIMPLE").found("nUCorrectors"))
+        {
+            nUCorr = readInt
+            (
+                mesh.solutionDict().subDict("SIMPLE").lookup("nUCorrectors")
+            );
+        }
+
+        if (nUCorr > 0)
+        {
+            pressureImplicitPorosity = true;
+            Info<< "Using pressure implicit porosity" << endl;
+        }
+        else
+        {
+            Info<< "Using pressure explicit porosity" << endl;
+        }
+    }