SF flame bc

0.org/C2H4

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

0.org/N2

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

0.org/O2

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

0.org/U

@@ -22,14 +22,8 @@ boundaryField
 {
     nzLow
     {
-        type            uniformFixedValue;
+        type            fixedValue;
-        uniformValue    table
+        value           uniform ( 0.0 0.2 0 );
-            (
-             (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
     {
@@ -38,15 +32,8 @@ boundaryField
     }
     coLow
     {
-        type            uniformFixedValue;
+        type            fixedValue;
-        uniformValue    table
+        value           uniform ( 0.0 0.2 0 );
-            (
-             (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 3bc47be714b0684d143c8aa5524ffb6baf87e1ec
+Subproject commit dc422748e3e9d84070ee8d1ca86b8a48ed39a238

constant/polyMesh/blockMeshDict

@@ -18,27 +18,36 @@ convertToMeters 1e-3; // 1 m / 1000 mm
 
 vertices
 (
- (  0.0000000000000e+00  0.0000000000000e+00 -0.0000000000000e+00 ) // 0
+(  0.0000000000000e+00 -5.5000000000000e+00 -0.0000000000000e+00 ) // 0
- (  1.1499562115238e+00  0.0000000000000e+00 -1.0035515823130e-02 ) // 1
+(  4.9998096153209e+00 -5.5000000000000e+00 -4.3632677491870e-02 ) // 1
- (  1.1499562115238e+00  3.0000000000000e+00 -1.0035515823130e-02 ) // 2
+(  4.9998096153209e+00  5.5000000000000e+00 -4.3632677491870e-02 ) // 2
- (  0.0000000000000e+00  3.0000000000000e+00 -0.0000000000000e+00 ) // 3
+(  0.0000000000000e+00  5.5000000000000e+00 -0.0000000000000e+00 ) // 3
- (  1.1499562115238e+00  0.0000000000000e+00  1.0035515823130e-02 ) // 4
+(  4.9998096153209e+00 -5.5000000000000e+00  4.3632677491870e-02 ) // 4
- (  1.1499562115238e+00  3.0000000000000e+00  1.0035515823130e-02 ) // 5
+(  4.9998096153209e+00  5.5000000000000e+00  4.3632677491870e-02 ) // 5
- (  1.1499562115238e+00  0.0000000000000e+00 -1.0035515823130e-02 ) // 6
+(  4.9998096153209e+00 -5.5000000000000e+00 -4.3632677491870e-02 ) // 6
- (  5.9997715383850e+00  0.0000000000000e+00 -5.2359212990244e-02 ) // 7
+(  9.9996192306417e+00 -5.5000000000000e+00 -8.7265354983739e-02 ) // 7
- (  5.9997715383850e+00  3.0000000000000e+00 -5.2359212990244e-02 ) // 8
+(  9.9996192306417e+00  5.5000000000000e+00 -8.7265354983739e-02 ) // 8
- (  1.1499562115238e+00  3.0000000000000e+00 -1.0035515823130e-02 ) // 9
+(  4.9998096153209e+00  5.5000000000000e+00 -4.3632677491870e-02 ) // 9
- (  1.1499562115238e+00  0.0000000000000e+00  1.0035515823130e-02 ) // 10
+(  4.9998096153209e+00 -5.5000000000000e+00  4.3632677491870e-02 ) // 10
- (  5.9997715383850e+00  0.0000000000000e+00  5.2359212990244e-02 ) // 11
+(  9.9996192306417e+00 -5.5000000000000e+00  8.7265354983739e-02 ) // 11
- (  5.9997715383850e+00  3.0000000000000e+00  5.2359212990244e-02 ) // 12
+(  9.9996192306417e+00  5.5000000000000e+00  8.7265354983739e-02 ) // 12
- (  1.1499562115238e+00  3.0000000000000e+00  1.0035515823130e-02 ) // 13
+(  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
 )
 ;
 
 blocks
 (
- hex ( 0 1 2 3 0 4 5 3 ) ( 54 90 1 ) simpleGrading ( 1 1 1 ) // inner
+ 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 ) ( 186 90 1 ) simpleGrading ( 1.5 1 1 ) // coflow
+ 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
 )
 ;
 
@@ -56,17 +65,54 @@ boundary
          (0 1 4 0) // inner down
         );
     }
-    outlet
+    nzUpp
     {
         type patch;
         faces
         (
          (2 3 3 5) // inner up
-         ( 8  9 13 12) // coflow up
+        );
-         ( 7  8 12 11) // coflow right
+    }
+    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
+        );
+    }
     front
     {
         type wedge;
@@ -74,6 +120,7 @@ boundary
         (
          (0 4 5 3) // inner front
          (10 11 12 13) // coflow front
+         (18 19 20 21) // outer front
         );
     }
     back
@@ -83,6 +130,7 @@ boundary
         (
          (0 3 2 1) // inner back
          (6 9 8 7) // coflow back
+         (14 17 16 15) // outer back
         );
     }
     axis
@@ -109,11 +157,28 @@ 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,14 +26,12 @@ thermoType
     specie          specie;
 }
 
-inertSpecie AR;
+inertSpecie N2;
 
 chemistryReader foamChemistryReader;
 
-foamChemistryFile       "$FOAM_CASE/constant/mech/GRI30/grimech30.foam";
+foamChemistryFile       "$FOAM_CASE/constant/mech/ethylene-air-2step-BFER/reactions";
 
-foamChemistryThermoFile "$FOAM_CASE/constant/mech/GRI30/thermo30.foam";
+foamChemistryThermoFile "$FOAM_CASE/constant/mech/ethylene-air-2step-BFER/thermos";
-
-foamTransportFile       "$FOAM_CASE/constant/mech/GRI30/transport30.foam";
 
 // ************************************************************************* //

job.mpi

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

system/controlDict

@@ -90,16 +90,6 @@ 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 240;
+numberOfSubdomains 100;
 
 // method         manual;
 method         simple;
@@ -24,7 +24,7 @@ method         simple;
 simpleCoeffs
 {
     // n               ( 4 5 1 );
-    n               ( 60 4 1 );
+    n               ( 100 1 1 );
     delta           0.001;
 }
 

system/fvOptions

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

system/fvSchemes

@@ -36,9 +36,6 @@ 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

@@ -1,23 +0,0 @@
-/*--------------------------------*- 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.0005 -1) (0.0005 0.001 1);
+            box (-0.001 -0.001 -1) (0.001 0.001 1);
         }
     }