From b2efc2094db91871fc834362ebb66b7c6fed0c8f Mon Sep 17 00:00:00 2001
From: ignis <ignis@postech.ac.kr>
Date: Tue, 8 Jan 2019 17:37:10 +0900
Subject: [PATCH] Axial length 3.8 cm, mixture-avg diffusivity

---
 constant/polyMesh/blockMeshDict   | 18 +++++++++---------
 constant/thermophysicalProperties |  6 ++++--
 job.mpi                           |  2 +-
 system/fvOptions                  |  8 ++++----
 system/fvSchemes                  |  3 +++
 system/topoSetDict                |  2 +-
 6 files changed, 22 insertions(+), 17 deletions(-)

diff --git a/constant/polyMesh/blockMeshDict b/constant/polyMesh/blockMeshDict
index 2f5b4ab..62fc7d4 100644
--- a/constant/polyMesh/blockMeshDict
+++ b/constant/polyMesh/blockMeshDict
@@ -20,25 +20,25 @@ vertices
 (
 (  0.0000000000000e+00  0.0000000000000e+00 -0.0000000000000e+00 ) // 0
 (  1.1499562115238e+00  0.0000000000000e+00 -1.0035515823130e-02 ) // 1
-(  1.1499562115238e+00  2.0000000000000e+02 -1.0035515823130e-02 ) // 2
-(  0.0000000000000e+00  2.0000000000000e+02 -0.0000000000000e+00 ) // 3
+(  1.1499562115238e+00  3.2000000000000e+01 -1.0035515823130e-02 ) // 2
+(  0.0000000000000e+00  3.2000000000000e+01 -0.0000000000000e+00 ) // 3
 (  1.1499562115238e+00  0.0000000000000e+00  1.0035515823130e-02 ) // 4
-(  1.1499562115238e+00  2.0000000000000e+02  1.0035515823130e-02 ) // 5
+(  1.1499562115238e+00  3.2000000000000e+01  1.0035515823130e-02 ) // 5
 (  1.1499562115238e+00  0.0000000000000e+00 -1.0035515823130e-02 ) // 6
 (  1.4999428845963e+01  0.0000000000000e+00 -1.3089803247561e-01 ) // 7
-(  1.4999428845963e+01  2.0000000000000e+02 -1.3089803247561e-01 ) // 8
-(  1.1499562115238e+00  2.0000000000000e+02 -1.0035515823130e-02 ) // 9
+(  1.4999428845963e+01  3.2000000000000e+01 -1.3089803247561e-01 ) // 8
+(  1.1499562115238e+00  3.2000000000000e+01 -1.0035515823130e-02 ) // 9
 (  1.1499562115238e+00  0.0000000000000e+00  1.0035515823130e-02 ) // 10
 (  1.4999428845963e+01  0.0000000000000e+00  1.3089803247561e-01 ) // 11
-(  1.4999428845963e+01  2.0000000000000e+02  1.3089803247561e-01 ) // 12
-(  1.1499562115238e+00  2.0000000000000e+02  1.0035515823130e-02 ) // 13
+(  1.4999428845963e+01  3.2000000000000e+01  1.3089803247561e-01 ) // 12
+(  1.1499562115238e+00  3.2000000000000e+01  1.0035515823130e-02 ) // 13
 )
 ;
 
 blocks
 (
- hex ( 0 1 2 3 0 4 5 3 ) ( 18 2000 1 ) simpleGrading ( 1 1 1 ) // inner
- hex ( 6 7 8 9 10 11 12 13 ) ( 102 2000 1 ) simpleGrading ( 3.8 1 1 ) // coflow
+ hex ( 0 1 2 3 0 4 5 3 ) ( 18 320 1 ) simpleGrading ( 1 1 1 ) // inner
+ hex ( 6 7 8 9 10 11 12 13 ) ( 102 320 1 ) simpleGrading ( 3.8 1 1 ) // coflow
 )
 ;
 
diff --git a/constant/thermophysicalProperties b/constant/thermophysicalProperties
index 9565965..00fadf1 100644
--- a/constant/thermophysicalProperties
+++ b/constant/thermophysicalProperties
@@ -30,8 +30,10 @@ inertSpecie AR;
 
 chemistryReader foamChemistryReader;
 
-foamChemistryFile       "$FOAM_CASE/constant/mech/GRI30/reactions";
+foamChemistryFile       "$FOAM_CASE/constant/mech/GRI30/grimech30.foam";
 
-foamChemistryThermoFile "$FOAM_CASE/constant/mech/GRI30/thermos";
+foamChemistryThermoFile "$FOAM_CASE/constant/mech/GRI30/thermo30.foam";
+
+foamTransportFile       "$FOAM_CASE/constant/mech/GRI30/transport30.foam";
 
 // ************************************************************************* //
diff --git a/job.mpi b/job.mpi
index 552b705..ee384f1 100644
--- a/job.mpi
+++ b/job.mpi
@@ -9,4 +9,4 @@
 
 # Launch MPI-based executable
 
-prun reactingFoam -parallel
+prun laminarReactingFoam -parallel
diff --git a/system/fvOptions b/system/fvOptions
index 71ae9e3..aead1cd 100644
--- a/system/fvOptions
+++ b/system/fvOptions
@@ -18,12 +18,12 @@ FoamFile
 source1
 {
     type            fixedTemperatureConstraint;
-    active          false;
+    active          true;
 
     fixedTemperatureConstraintCoeffs
     {
         timeStart       0.0;
-        duration        0.03;
+        duration        0.05;
         selectionMode   cellSet;
         cellSet         ignitionCells;
 
@@ -31,8 +31,8 @@ source1
         temperature     table
             2
             (
-             (0 300)
-             (0.015 2240)
+             (0 1800)
+             (0.0001 2240)
             )
             ;
     }
diff --git a/system/fvSchemes b/system/fvSchemes
index 5370070..8f9f7a5 100644
--- a/system/fvSchemes
+++ b/system/fvSchemes
@@ -36,6 +36,9 @@ divSchemes
     div(phi,epsilon) Gauss limitedLinear 1;
     div(phi,k) Gauss limitedLinear 1;
     div(((rho*nuEff)*dev2(T(grad(U))))) Gauss linear;
+    div((mu*dev2(T(grad(U))))) Gauss linear;
+    div((mu*(U&(dev2(T(grad(U)))+grad(U))))) Gauss linear;
+    div(rho*Sum(hYV)) Gauss linear;
 }
 
 laplacianSchemes
diff --git a/system/topoSetDict b/system/topoSetDict
index 8500d75..035988b 100644
--- a/system/topoSetDict
+++ b/system/topoSetDict
@@ -416,7 +416,7 @@ actions
         source boxToCell;
         sourceInfo
         {
-            box (-0.001 0.001 -1) (0.001 0.002 1);
+            box (-0.001 0.0005 -1) (0.0005 0.001 1);
         }
     }