Solution converged

TODO

@@ -0,0 +1,2 @@
+* fvSolution: cellLimited
+* collapseDict: collapseEdges

run.sh

@@ -1,24 +0,0 @@
-#!/usr/bin/env bash
-
-path=.
-#build/simple-cubic/0.1
-
-# python src/genmesh.py
-# python src/prefoam.py
-
-rm -r postProcessing processor* *.log logs constant/polymesh
-
-ideasUnvFoam -case $path mesh.unv
-polyDualMesh 70 -overwrite
-checkMesh -case $path -allGeometry -allTopology
-foamDictionary -case $path constant/polyMesh/boundary -entry entry0.wall.type -set wall
-
-
-potentialFoam -case $path
-
-decomposePar -case $path
-
-foamDictionary -case $path 0/U -entry boundary.inlet.type -set pressureInletVelocity
-foamDictionary -case $path 0/U -entry boundary.inlet.value -set 'uniform (0 0 0)'
-mpirun -np 4 --oversubscribe simpleFoam -parallel -case $path > $path/simpleFoam.log
-

src/simple-cubic/system/fvSolution

@@ -57,7 +57,7 @@ solvers
 
 potentialFlow
 {
-    nNonOrthogonalCorrectors 10;
+    nNonOrthogonalCorrectors 20;
     PhiRefCell  0;
     PhiRefPoint 0;
     PhiRefValue 0;

src/tools/changeInitial.sh

@@ -0,0 +1,7 @@
+#!/usr/bin/env bash
+
+for n in {0..3}; do
+
+    foamDictionary "processor${n}/0/U" -entry boundaryField.inlet.type -set pressureInletVelocity
+    foamDictionary "processor${n}/0/U" -entry boundaryField.inlet.value -set 'uniform (0 0 0)'
+done

src/tools/run.sh

@@ -0,0 +1,36 @@
+#!/usr/bin/env bash
+
+path=.
+#build/simple-cubic/0.1
+
+# python src/genmesh.py
+# python src/prefoam.py
+
+# Clean case directory
+rm -rf postProcessing processor* *.log logs *.obj constant/polyMesh
+
+# Export and transform mesh
+ideasUnvFoam -case $path mesh.unv
+transformPoints -scale '(0.001 0.001 0.001)'
+#polyDualMesh 70 -overwrite
+checkMesh -case $path -allGeometry -allTopology > checkMesh.log
+
+# Change boundary type for mesh
+foamDictionary -case $path constant/polyMesh/boundary -entry entry0.wall.type -set wall
+
+# Case decomposition
+decomposePar -case $path
+
+# Initial approximation
+potentialFoam -case $path -parallel
+
+# Change boundary type for simpleFoam
+for n in {0..3}; do
+
+    foamDictionary "processor${n}/0/U" -entry boundaryField.inlet.type -set pressureInletVelocity
+    foamDictionary "processor${n}/0/U" -entry boundaryField.inlet.value -set 'uniform (0 0 0)'
+done
+
+# Main calculation
+mpirun -np 4 --oversubscribe simpleFoam -parallel -case $path > $path/simpleFoam.log
+