Axial length 3 mm, fine mesh

0.org/AR

@@ -23,12 +23,12 @@ boundaryField
     nzLow
     {
         type            fixedValue;
-        value           uniform 0;
+        value           uniform 0.0;
     }
     nzUpp
     {
         type            fixedValue;
-        value           uniform 0.724242;
+        value           uniform 0.440011;
     }
     coLow
     {

0.org/CH4

@@ -23,7 +23,7 @@ boundaryField
     nzLow
     {
         type            fixedValue;
-        value           uniform 1;
+        value           uniform 0.051354;
     }
     nzUpp
     {

0.org/O2

@@ -23,12 +23,12 @@ boundaryField
     nzLow
     {
         type            fixedValue;
-        value           uniform 0;
+        value           uniform 0.948646;
     }
     nzUpp
     {
         type            fixedValue;
-        value           uniform 0.275758;
+        value           uniform 0.559989;
     }
     coLow
     {

0.org/U

@@ -22,8 +22,14 @@ boundaryField
 {
     nzLow
     {
-        type            fixedValue;
-        value           uniform ( 0.0 0.2 0 );
+        type            uniformFixedValue;
+        uniformValue    table
+            (
+             (0.0 (0 0.2 0))
+             (0.06 (0 0.2 0))
+             (0.2 ( 0.0 2 0 ))
+            );
+        // value           uniform ( 0.0 4.34682536445 0 );
     }
     nzUpp
     {
@@ -32,8 +38,15 @@ boundaryField
     }
     coLow
     {
-        type            fixedValue;
-        value           uniform ( 0.0 0.2 0 );
+        type            uniformFixedValue;
+        uniformValue    table
+            (
+             (0.0 (0 0.2 0))
+             (0.3 (0 0.2 0))
+             (0.5 ( 0.0 4.34682536445 0 ))
+            );
+        // type            fixedValue;
+        // value           uniform ( 0.0 4.34682536445 0 );
     }
     coUpp
     {

constant/mech

@@ -1 +1 @@
-Subproject commit dc422748e3e9d84070ee8d1ca86b8a48ed39a238
+Subproject commit 3bc47be714b0684d143c8aa5524ffb6baf87e1ec

constant/polyMesh/blockMeshDict

@@ -18,36 +18,27 @@ convertToMeters 1e-3; // 1 m / 1000 mm
 
 vertices
 (
-(  0.0000000000000e+00 -5.5000000000000e+00 -0.0000000000000e+00 ) // 0
-(  4.9998096153209e+00 -5.5000000000000e+00 -4.3632677491870e-02 ) // 1
-(  4.9998096153209e+00  5.5000000000000e+00 -4.3632677491870e-02 ) // 2
-(  0.0000000000000e+00  5.5000000000000e+00 -0.0000000000000e+00 ) // 3
-(  4.9998096153209e+00 -5.5000000000000e+00  4.3632677491870e-02 ) // 4
-(  4.9998096153209e+00  5.5000000000000e+00  4.3632677491870e-02 ) // 5
-(  4.9998096153209e+00 -5.5000000000000e+00 -4.3632677491870e-02 ) // 6
-(  9.9996192306417e+00 -5.5000000000000e+00 -8.7265354983739e-02 ) // 7
-(  9.9996192306417e+00  5.5000000000000e+00 -8.7265354983739e-02 ) // 8
-(  4.9998096153209e+00  5.5000000000000e+00 -4.3632677491870e-02 ) // 9
-(  4.9998096153209e+00 -5.5000000000000e+00  4.3632677491870e-02 ) // 10
-(  9.9996192306417e+00 -5.5000000000000e+00  8.7265354983739e-02 ) // 11
-(  9.9996192306417e+00  5.5000000000000e+00  8.7265354983739e-02 ) // 12
-(  4.9998096153209e+00  5.5000000000000e+00  4.3632677491870e-02 ) // 13
-(  9.9996192306417e+00 -5.5000000000000e+00 -8.7265354983739e-02 ) // 14
-(  3.9998476922567e+01 -5.5000000000000e+00 -3.4906141993496e-01 ) // 15
-(  3.9998476922567e+01  5.5000000000000e+00 -3.4906141993496e-01 ) // 16
-(  9.9996192306417e+00  5.5000000000000e+00 -8.7265354983739e-02 ) // 17
-(  9.9996192306417e+00 -5.5000000000000e+00  8.7265354983739e-02 ) // 18
-(  3.9998476922567e+01 -5.5000000000000e+00  3.4906141993496e-01 ) // 19
-(  3.9998476922567e+01  5.5000000000000e+00  3.4906141993496e-01 ) // 20
-(  9.9996192306417e+00  5.5000000000000e+00  8.7265354983739e-02 ) // 21
+ (  0.0000000000000e+00  0.0000000000000e+00 -0.0000000000000e+00 ) // 0
+ (  1.1499562115238e+00  0.0000000000000e+00 -1.0035515823130e-02 ) // 1
+ (  1.1499562115238e+00  3.0000000000000e+00 -1.0035515823130e-02 ) // 2
+ (  0.0000000000000e+00  3.0000000000000e+00 -0.0000000000000e+00 ) // 3
+ (  1.1499562115238e+00  0.0000000000000e+00  1.0035515823130e-02 ) // 4
+ (  1.1499562115238e+00  3.0000000000000e+00  1.0035515823130e-02 ) // 5
+ (  1.1499562115238e+00  0.0000000000000e+00 -1.0035515823130e-02 ) // 6
+ (  5.9997715383850e+00  0.0000000000000e+00 -5.2359212990244e-02 ) // 7
+ (  5.9997715383850e+00  3.0000000000000e+00 -5.2359212990244e-02 ) // 8
+ (  1.1499562115238e+00  3.0000000000000e+00 -1.0035515823130e-02 ) // 9
+ (  1.1499562115238e+00  0.0000000000000e+00  1.0035515823130e-02 ) // 10
+ (  5.9997715383850e+00  0.0000000000000e+00  5.2359212990244e-02 ) // 11
+ (  5.9997715383850e+00  3.0000000000000e+00  5.2359212990244e-02 ) // 12
+ (  1.1499562115238e+00  3.0000000000000e+00  1.0035515823130e-02 ) // 13
 )
 ;
 
 blocks
 (
- hex ( 0 1 2 3 0 4 5 3 ) ( 50 110 1 ) simpleGrading ( 1 1 1 ) // inner
- hex ( 6 7 8 9 10 11 12 13 ) ( 50 110 1 ) simpleGrading ( 1 1 1 ) // coflow
- hex ( 14 15 16 17 18 19 20 21 ) ( 300 110 1 ) simpleGrading ( 1 1 1 ) // outer
+ hex ( 0 1 2 3 0 4 5 3 ) ( 54 90 1 ) simpleGrading ( 1 1 1 ) // inner
+ hex ( 6 7 8 9 10 11 12 13 ) ( 186 90 1 ) simpleGrading ( 1.5 1 1 ) // coflow
 )
 ;
 
@@ -65,52 +56,15 @@ boundary
          (0 1 4 0) // inner down
         );
     }
-    nzUpp
-    {
-        type patch;
-        faces
-        (
-         (2 3 3 5) // inner up
-        );
-    }
-    coLow
-    {
-        type patch;
-        faces
-        (
-         ( 6  7 11 10) // coflow down
-        );
-    }
-    coUpp
-    {
-        type patch;
-        faces
-        (
-         ( 8  9 13 12) // coflow up
-        );
-    }
-    waLow
-    {
-        type patch;
-        faces
-        (
-         (14 15 19 18) // outer down
-        );
-    }
-    waUpp
-    {
-        type patch;
-        faces
-        (
-         (16 17 21 20) // outer up
-        );
-    }
     outlet
     {
         type patch;
         faces
         (
-         (15 16 20 19) // outer right
+         (2 3 3 5) // inner up
+         ( 8  9 13 12) // coflow up
+         ( 7  8 12 11) // coflow right
+         ( 6  7 11 10) // coflow down
         );
     }
     front
@@ -120,7 +74,6 @@ boundary
         (
          (0 4 5 3) // inner front
          (10 11 12 13) // coflow front
-         (18 19 20 21) // outer front
         );
     }
     back
@@ -130,7 +83,6 @@ boundary
         (
          (0 3 2 1) // inner back
          (6 9 8 7) // coflow back
-         (14 17 16 15) // outer back
         );
     }
     axis
@@ -157,28 +109,11 @@ boundary
          ( 9  6 10 13) // coflow left
         );
     }
-    dummy3
-    {
-        type empty;
-        faces
-        (
-         ( 7  8 12 11) // coflow right
-        );
-    }
-    dummy4
-    {
-        type empty;
-        faces
-        (
-         (17 14 18 21) // outer left
-        );
-    }
 );
 
 mergePatchPairs
 (
   (dummy1 dummy2)
-  (dummy3 dummy4)
 );
 
 // ************************************************************************* //

constant/thermophysicalProperties

@@ -26,12 +26,14 @@ thermoType
     specie          specie;
 }
 
-inertSpecie N2;
+inertSpecie AR;
 
 chemistryReader foamChemistryReader;
 
-foamChemistryFile       "$FOAM_CASE/constant/mech/ethylene-air-2step-BFER/reactions";
+foamChemistryFile       "$FOAM_CASE/constant/mech/GRI30/grimech30.foam";
 
-foamChemistryThermoFile "$FOAM_CASE/constant/mech/ethylene-air-2step-BFER/thermos";
+foamChemistryThermoFile "$FOAM_CASE/constant/mech/GRI30/thermo30.foam";
+
+foamTransportFile       "$FOAM_CASE/constant/mech/GRI30/transport30.foam";
 
 // ************************************************************************* //

job.mpi

@@ -1,12 +1,12 @@
 #!/bin/bash
 
-#SBATCH -J testLaminar        # Job name
-#SBATCH -o job.%j.out         # Name of stdout output file (%j expands to jobId)
-#SBATCH -n 100                # Total number of mpi tasks requested
-#SBATCH --ntasks-per-node=20
+#SBATCH -J iCoflow        # Job name
+#SBATCH -o log.%j.out         # Name of stdout output file (%j expands to jobId)
+#SBATCH -n 240                # Total number of mpi tasks requested
+#SBATCH --ntasks-per-node=24
 #SBATCH --mail-user=ignis@postech.ac.kr
 #SBATCH --mail-type=ALL
 
 # Launch MPI-based executable
 
-prun reactingFoam -parallel
+prun laminarReactingFoam -parallel

system/controlDict

@@ -90,6 +90,16 @@ functions
              Phi
             );
     }
+
+    writeObjects2
+    {
+        type        writeObjects;
+        libs        ("libutilityFunctionObjects.so");
+        exclusiveWriting     false;
+        writeControl   clockTime;
+        writeInterval  900;
+        objects     (U T OH CH);
+    }
 }
 
 // ************************************************************************* //

system/decomposeParDict

@@ -15,7 +15,7 @@ FoamFile
 }
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
-numberOfSubdomains 100;
+numberOfSubdomains 240;
 
 // method         manual;
 method         simple;
@@ -24,7 +24,7 @@ method         simple;
 simpleCoeffs
 {
     // n               ( 4 5 1 );
-    n               ( 100 1 1 );
+    n               ( 60 4 1 );
     delta           0.001;
 }
 

system/fvOptions

@@ -23,7 +23,7 @@ source1
     fixedTemperatureConstraintCoeffs
     {
         timeStart       0.0;
-        duration        0.1;
+        duration        0.05;
         selectionMode   cellSet;
         cellSet         ignitionCells;
 
@@ -31,8 +31,8 @@ source1
         temperature     table
             2
             (
-             (0 300)
-             (0.05 2000)
+             (0 1800)
+             (0.0001 2240)
             )
             ;
     }

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

system/mapFieldsDict

@@ -0,0 +1,23 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  4.x                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       dictionary;
+    location    "system";
+    object      mapFieldsDict;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+patchMap        (nzLow nzLow);
+
+cuttingPatches  ();
+
+
+// ************************************************************************* //

system/topoSetDict

@@ -416,7 +416,7 @@ actions
         source boxToCell;
         sourceInfo
         {
-            box (-0.001 -0.001 -1) (0.001 0.001 1);
+            box (-0.001 0.0005 -1) (0.0005 0.001 1);
         }
     }