Merge experiment branch

Author: phresher

.gitignore

@@ -1,3 +1,4 @@
+*.grf
 *.gz
 log.*
 postProcessing
@@ -13,3 +14,4 @@ processor*
 polyMesh/
 *.OpenFOAM
 dynamicCode
+sequencedVTK/

README.md

@@ -15,10 +15,6 @@ Here's a brief description:
 
   Abrupt contraction of a rectangular channel where cavitation occures. The *blockMeshDict* is written with variables. *interPhaseChangeFoam*, *pisoFoam* and  *pimpleFoam* are used.
 
-* **abruptContractionPar**
-
-  Parallel edition of abruptContraction.
-
 * **waterChannel**
 
   Two-phase open channel flow. It's a fork  of the original interFoam tutorial *waterChannel*. The *surfaces* function is configured to generate water surfaces (namely, the isosurfaces where *alpha.water* is equal to 0.5). 

abruptContraction/interPhaseChangeFoam/0/U

@@ -27,14 +27,12 @@ boundaryField
         referenceField  uniform (5 0 0);
         fluctuationScale (0.02 0.01 0.01);
         value           uniform (5 0 0);
- 
     }
     outlet
     {
         type            inletOutlet;
         inletValue      uniform (0 0 0);
         value           uniform (0 0 0);
- 
     }
     side
     {
@@ -46,7 +44,7 @@ boundaryField
     }
     symmetry
     {
-        type            noSlip;
+        type            symmetry;
     }
 }
 

abruptContraction/interPhaseChangeFoam/0/alpha.water

@@ -41,7 +41,7 @@ boundaryField
     }
     symmetry
     {
-        type            zeroGradient;
+        type            symmetry;
     }
 }
 

abruptContraction/interPhaseChangeFoam/0/k

@@ -47,7 +47,7 @@ boundaryField
 
     symmetry
     {
-        type            fixedValue;
+        type            symmetry;
         value           uniform 0;
     }
 }

abruptContraction/interPhaseChangeFoam/0/nut

@@ -42,7 +42,7 @@ boundaryField
 
     symmetry
     {
-        type            zeroGradient;
+        type            symmetry;
     }
 }
 

abruptContraction/interPhaseChangeFoam/0/p_rgh

@@ -42,7 +42,7 @@ boundaryField
     }
     symmetry
     {
-        type            zeroGradient;
+        type            symmetry;
     }
 }
 

abruptContraction/interPhaseChangeFoam/0/phi

@@ -43,7 +43,7 @@ boundaryField
     }
     symmetry
     {
-        type            calculated;
+        type            symmetry;
         value           uniform 0;
     }
 }

abruptContraction/interPhaseChangeFoam/Allmesh

@@ -0,0 +1,29 @@
+#!/bin/sh
+cd ${0%/*} || exit 1    # Run from this directory
+
+# Source tutorial run functions
+. $WM_PROJECT_DIR/bin/tools/RunFunctions
+
+refineMeshByCellSet()
+{
+   while [ $# -ge 1 ]
+   do
+      if [ ! -e log.refineMesh.$1 ]
+      then
+          echo "Creating cell set for primary zone - $1"
+          topoSet -dict system/topoSetDict.$1 > log.topoSet.$1 2>&1
+          echo "Refining primary zone - $1"
+          refineMesh -dict system/refineMeshDict.$1 -overwrite \
+              > log.refineMesh.$1 2>&1
+      fi
+      shift
+   done
+}
+
+runApplication blockMesh
+refineMeshByCellSet 1
+refineWallLayer -overwrite '(walls)' 0.2 >> log.refineWallLayer.1
+
+
+
+#------------------------------------------------------------------------------

abruptContraction/interPhaseChangeFoam/Allrun

@@ -4,15 +4,10 @@ cd ${0%/*} || exit 1    # Run from this directory
 # Source tutorial run functions
 . $WM_PROJECT_DIR/bin/tools/RunFunctions
 
-# Generate the base block mesh
-runApplication blockMesh
+# Generate the base block mesh and refine
+./Allmesh
 
-# cp -r 0.orig 0
-
-# Initialise with potentialFoam solution
-# runApplication potentialFoam -pName p_rgh
-
-# Run the solver
+# Run the solver in parallel
 runApplication `getApplication`
 
 # ----------------------------------------------------------------- end-of-file