Tutorial case for ore reducing BF.

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/CFD/0/Ksl

@@ -0,0 +1,54 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  1.6                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volScalarField;
+    location    "0";
+    object      Ksl;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [1 -3 -1 0 0 0 0];
+
+internalField   uniform 0;
+
+boundaryField
+{
+    top
+    {
+        type            zeroGradient;
+    }
+
+    bottom
+    {
+        type            zeroGradient;
+    }
+
+    wall
+    {
+        type            zeroGradient;
+    }
+
+    inner
+    {
+        type            zeroGradient;
+    }
+    inletLeft
+    {
+        type            zeroGradient;
+    }
+    inletRight
+    {
+        type            zeroGradient;
+    }
+}
+
+
+// ************************************************************************* //

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/CFD/0/KslNext

@@ -0,0 +1,53 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  1.6                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volScalarField;
+    location    "0";
+    object      KslNext;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [1 -3 -1 0 0 0 0];
+
+internalField   uniform 0;
+
+boundaryField
+{
+    top
+    {
+        type            zeroGradient;
+    }
+
+    bottom
+    {
+        type            zeroGradient;
+    }
+
+    wall
+    {
+        type            zeroGradient;
+    }
+
+    inner
+    {
+        type            zeroGradient;
+    }
+    inletLeft
+    {
+        type            zeroGradient;
+    }
+    inletRight
+    {
+        type            zeroGradient;
+    }
+}
+
+// ************************************************************************* //

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/CFD/0/KslPrev

@@ -0,0 +1,54 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  1.6                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volScalarField;
+    location    "0";
+    object      KslPrev;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [1 -3 -1 0 0 0 0];
+
+internalField   uniform 0;
+
+boundaryField
+{
+    top
+    {
+        type            zeroGradient;
+    }
+
+    bottom
+    {
+        type            zeroGradient;
+    }
+
+    wall
+    {
+        type            zeroGradient;
+    }
+
+    inner
+    {
+        type            zeroGradient;
+    }
+    inletLeft
+    {
+        type            zeroGradient;
+    }
+    inletRight
+    {
+        type            zeroGradient;
+    }
+}
+
+
+// ************************************************************************* //

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/CFD/0/T

@@ -0,0 +1,56 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  1.6                                   |
+|   \\  /    A nd           | Web:      http://www.OpenFOAM.org               |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volScalarField;
+    object      T;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 0 0 1 0 0 0];
+
+internalField   uniform 293;
+
+boundaryField
+{
+    top
+    {
+        type            zeroGradient;
+    }
+
+    bottom
+    {
+        type            fixedValue;
+        value		uniform 293;
+    }
+
+    wall
+    {
+        type            fixedValue;
+        value		uniform 293;
+    }
+
+    inner
+    {
+        type            zeroGradient;
+    }
+    inletLeft
+    {
+        type            fixedValue;
+        value		uniform 2600;
+    }
+    inletRight
+    {
+        type            fixedValue;
+        value		uniform 2600;
+    }
+}
+
+// ************************************************************************* //

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/CFD/0/U

@@ -0,0 +1,66 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  2.3.0                                 |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volVectorField;
+    object      U;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 1 -1 0 0 0 0];
+
+internalField   uniform (0 0 0);
+
+boundaryField
+{
+    bottom
+    {
+        type            uniformFixedValue;
+        uniformValue    constant (0 0 0);
+    }
+
+   
+    top
+    {
+        type            zeroGradient;
+    }
+
+    wall
+    {
+        type            uniformFixedValue;
+        uniformValue    constant (0 0 0);
+    }
+
+    inner
+    {
+        type            uniformFixedValue;
+        uniformValue    constant (0 0 0);
+    }
+
+    defaultFaces
+    {
+        type            empty;
+    }
+
+    inletLeft
+    {
+        type            uniformFixedValue;
+        uniformValue    constant (200 0 0);
+    }
+    inletRight
+    {
+        type            uniformFixedValue;
+        uniformValue    constant (-200 0 0);
+    }
+
+
+}
+
+// ************************************************************************* //

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/CFD/0/Us

@@ -0,0 +1,54 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  1.6                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volVectorField;
+    location    "0";
+    object      Us;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 1 -1 0 0 0 0];
+
+internalField   uniform (0 0 0);
+
+boundaryField
+{
+    top
+    {
+        type            zeroGradient;
+    }
+
+    bottom
+    {
+        type            zeroGradient;
+    }
+
+    wall
+    {
+        type            zeroGradient;
+    }
+
+    inner
+    {
+        type            zeroGradient;
+    }
+    inletLeft
+    {
+        type            zeroGradient;
+    }
+    inletRight
+    {
+        type            zeroGradient;
+    }
+}
+
+
+// ************************************************************************* //

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/CFD/0/UsNext

@@ -0,0 +1,55 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  1.6                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volVectorField;
+    location    "0";
+    object      UsNext;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 1 -1 0 0 0 0];
+
+internalField   uniform (0 0 0);
+
+boundaryField
+{
+    top
+    {
+        type            zeroGradient;
+    }
+
+    bottom
+    {
+        type            zeroGradient;
+    }
+
+    wall
+    {
+        type            zeroGradient;
+    }
+
+    inner
+    {
+        type            zeroGradient;
+    }
+    inletLeft
+    {
+        type            zeroGradient;
+    }
+    inletRight
+    {
+        type            zeroGradient;
+    }
+
+}
+
+
+// ************************************************************************* //

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/CFD/0/UsPrev

@@ -0,0 +1,55 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  1.6                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volVectorField;
+    location    "0";
+    object      UsPrev;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 1 -1 0 0 0 0];
+
+internalField   uniform (0 0 0);
+
+boundaryField
+{
+    top
+    {
+        type            zeroGradient;
+    }
+
+    bottom
+    {
+        type            zeroGradient;
+    }
+
+    wall
+    {
+        type            zeroGradient;
+    }
+
+    inner
+    {
+        type            zeroGradient;
+    }
+    inletLeft
+    {
+        type            zeroGradient;
+    }
+    inletRight
+    {
+        type            zeroGradient;
+    }
+
+
+}
+
+// ************************************************************************* //

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/CFD/0/_wallQFactor

@@ -0,0 +1,63 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  1.6                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volScalarField;
+    location    "0";
+    object      voidfraction;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 0 0 0 0 0 0];
+
+internalField   uniform 1;
+
+boundaryField
+{
+   bottom
+    {
+        type            zeroGradient;
+    }
+
+    
+    top
+    {
+        type            zeroGradient;
+
+    }
+
+    wall
+    {
+        type            fixedValue;
+        value           uniform 10.0;
+    }
+
+    inner
+    {
+        type            zeroGradient;
+    }
+
+    defaultFaces
+    {
+        type            empty;
+    }
+
+    inletLeft
+    {
+        type            zeroGradient;
+    }
+    inletRight
+    {
+        type            zeroGradient;
+    }
+
+}
+
+// ************************************************************************* //

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/CFD/0/addSource

@@ -0,0 +1,53 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  1.6                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volScalarField;
+    object      addSource;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [1 -1 -3 0 0 0 0]; 
+
+internalField   uniform 0;  
+
+boundaryField
+{
+    top
+    {
+        type            zeroGradient;
+    }
+
+    bottom
+    {
+        type            zeroGradient;
+    }
+
+    wall
+    {
+        type            zeroGradient;
+    }
+
+    inner
+    {
+        type            zeroGradient;
+    }
+    inletLeft
+    {
+        type            zeroGradient;
+    }
+    inletRight
+    {
+        type            zeroGradient;
+    }
+}
+
+
+// ************************************************************************* //

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/CFD/0/orig.0/Ksl

@@ -0,0 +1,54 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  1.6                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volScalarField;
+    location    "0";
+    object      Ksl;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [1 -3 -1 0 0 0 0];
+
+internalField   uniform 0;
+
+boundaryField
+{
+    top
+    {
+        type            zeroGradient;
+    }
+
+    bottom
+    {
+        type            zeroGradient;
+    }
+
+    wall
+    {
+        type            zeroGradient;
+    }
+
+    inner
+    {
+        type            zeroGradient;
+    }
+    inletLeft
+    {
+        type            zeroGradient;
+    }
+    inletRight
+    {
+        type            zeroGradient;
+    }
+}
+
+
+// ************************************************************************* //

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/CFD/0/orig.0/KslNext

@@ -0,0 +1,53 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  1.6                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volScalarField;
+    location    "0";
+    object      KslNext;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [1 -3 -1 0 0 0 0];
+
+internalField   uniform 0;
+
+boundaryField
+{
+    top
+    {
+        type            zeroGradient;
+    }
+
+    bottom
+    {
+        type            zeroGradient;
+    }
+
+    wall
+    {
+        type            zeroGradient;
+    }
+
+    inner
+    {
+        type            zeroGradient;
+    }
+    inletLeft
+    {
+        type            zeroGradient;
+    }
+    inletRight
+    {
+        type            zeroGradient;
+    }
+}
+
+// ************************************************************************* //

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/CFD/0/orig.0/KslPrev

@@ -0,0 +1,54 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  1.6                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volScalarField;
+    location    "0";
+    object      KslPrev;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [1 -3 -1 0 0 0 0];
+
+internalField   uniform 0;
+
+boundaryField
+{
+    top
+    {
+        type            zeroGradient;
+    }
+
+    bottom
+    {
+        type            zeroGradient;
+    }
+
+    wall
+    {
+        type            zeroGradient;
+    }
+
+    inner
+    {
+        type            zeroGradient;
+    }
+    inletLeft
+    {
+        type            zeroGradient;
+    }
+    inletRight
+    {
+        type            zeroGradient;
+    }
+}
+
+
+// ************************************************************************* //

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/CFD/0/orig.0/T

@@ -0,0 +1,56 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  1.6                                   |
+|   \\  /    A nd           | Web:      http://www.OpenFOAM.org               |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volScalarField;
+    object      T;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 0 0 1 0 0 0];
+
+internalField   uniform 293;
+
+boundaryField
+{
+    top
+    {
+        type            zeroGradient;
+    }
+
+    bottom
+    {
+        type            fixedValue;
+        value		uniform 293;
+    }
+
+    wall
+    {
+        type            fixedValue;
+        value		uniform 293;
+    }
+
+    inner
+    {
+        type            zeroGradient;
+    }
+    inletLeft
+    {
+        type            fixedValue;
+        value		uniform 2600;
+    }
+    inletRight
+    {
+        type            fixedValue;
+        value		uniform 2600;
+    }
+}
+
+// ************************************************************************* //

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/CFD/0/orig.0/U

@@ -0,0 +1,66 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  2.3.0                                 |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volVectorField;
+    object      U;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 1 -1 0 0 0 0];
+
+internalField   uniform (0 0 0);
+
+boundaryField
+{
+    bottom
+    {
+        type            uniformFixedValue;
+        uniformValue    constant (0 0 0);
+    }
+
+   
+    top
+    {
+        type            zeroGradient;
+    }
+
+    wall
+    {
+        type            uniformFixedValue;
+        uniformValue    constant (0 0 0);
+    }
+
+    inner
+    {
+        type            uniformFixedValue;
+        uniformValue    constant (0 0 0);
+    }
+
+    defaultFaces
+    {
+        type            empty;
+    }
+
+    inletLeft
+    {
+        type            uniformFixedValue;
+        uniformValue    constant (200 0 0);
+    }
+    inletRight
+    {
+        type            uniformFixedValue;
+        uniformValue    constant (-200 0 0);
+    }
+
+
+}
+
+// ************************************************************************* //

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/CFD/0/orig.0/Us

@@ -0,0 +1,54 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  1.6                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volVectorField;
+    location    "0";
+    object      Us;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 1 -1 0 0 0 0];
+
+internalField   uniform (0 0 0);
+
+boundaryField
+{
+    top
+    {
+        type            zeroGradient;
+    }
+
+    bottom
+    {
+        type            zeroGradient;
+    }
+
+    wall
+    {
+        type            zeroGradient;
+    }
+
+    inner
+    {
+        type            zeroGradient;
+    }
+    inletLeft
+    {
+        type            zeroGradient;
+    }
+    inletRight
+    {
+        type            zeroGradient;
+    }
+}
+
+
+// ************************************************************************* //

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/CFD/0/orig.0/UsNext

@@ -0,0 +1,55 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  1.6                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volVectorField;
+    location    "0";
+    object      UsNext;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 1 -1 0 0 0 0];
+
+internalField   uniform (0 0 0);
+
+boundaryField
+{
+    top
+    {
+        type            zeroGradient;
+    }
+
+    bottom
+    {
+        type            zeroGradient;
+    }
+
+    wall
+    {
+        type            zeroGradient;
+    }
+
+    inner
+    {
+        type            zeroGradient;
+    }
+    inletLeft
+    {
+        type            zeroGradient;
+    }
+    inletRight
+    {
+        type            zeroGradient;
+    }
+
+}
+
+
+// ************************************************************************* //

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/CFD/0/orig.0/UsPrev

@@ -0,0 +1,55 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  1.6                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volVectorField;
+    location    "0";
+    object      UsPrev;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 1 -1 0 0 0 0];
+
+internalField   uniform (0 0 0);
+
+boundaryField
+{
+    top
+    {
+        type            zeroGradient;
+    }
+
+    bottom
+    {
+        type            zeroGradient;
+    }
+
+    wall
+    {
+        type            zeroGradient;
+    }
+
+    inner
+    {
+        type            zeroGradient;
+    }
+    inletLeft
+    {
+        type            zeroGradient;
+    }
+    inletRight
+    {
+        type            zeroGradient;
+    }
+
+
+}
+
+// ************************************************************************* //

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/CFD/0/orig.0/_wallQFactor

@@ -0,0 +1,63 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  1.6                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volScalarField;
+    location    "0";
+    object      voidfraction;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 0 0 0 0 0 0];
+
+internalField   uniform 1;
+
+boundaryField
+{
+   bottom
+    {
+        type            zeroGradient;
+    }
+
+    
+    top
+    {
+        type            zeroGradient;
+
+    }
+
+    wall
+    {
+        type            fixedValue;
+        value           uniform 10.0;
+    }
+
+    inner
+    {
+        type            zeroGradient;
+    }
+
+    defaultFaces
+    {
+        type            empty;
+    }
+
+    inletLeft
+    {
+        type            zeroGradient;
+    }
+    inletRight
+    {
+        type            zeroGradient;
+    }
+
+}
+
+// ************************************************************************* //

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/CFD/0/orig.0/addSource

@@ -0,0 +1,53 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  1.6                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volScalarField;
+    object      addSource;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [1 -1 -3 0 0 0 0]; 
+
+internalField   uniform 0;  
+
+boundaryField
+{
+    top
+    {
+        type            zeroGradient;
+    }
+
+    bottom
+    {
+        type            zeroGradient;
+    }
+
+    wall
+    {
+        type            zeroGradient;
+    }
+
+    inner
+    {
+        type            zeroGradient;
+    }
+    inletLeft
+    {
+        type            zeroGradient;
+    }
+    inletRight
+    {
+        type            zeroGradient;
+    }
+}
+
+
+// ************************************************************************* //

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/CFD/0/orig.0/p

@@ -0,0 +1,60 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  2.3.0                                 |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volScalarField;
+    object      p;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [1 -1 -2 0 0 0 0];
+
+internalField   uniform 1.0e5;
+
+boundaryField
+{
+    bottom
+    {
+        type             zeroGradient;
+    }
+
+    
+    top
+    {
+        type            fixedValue;
+        value           uniform 1.0e5;
+    }
+
+    wall
+    {
+        type            zeroGradient;
+    }
+
+    inner
+    {
+        type            zeroGradient;
+    }
+
+    defaultFaces
+    {
+        type            empty;
+    }
+
+    inletLeft
+    {
+        type            zeroGradient;
+    }
+    inletRight
+    {
+        type            zeroGradient;
+    }
+}
+
+// ************************************************************************* //

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/CFD/0/orig.0/partEffThermCondField

@@ -0,0 +1,54 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  1.6                                   |
+|   \\  /    A nd           | Web:      http://www.OpenFOAM.org               |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volScalarField;
+    object      partEffThermCondField;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [1 1 -3 -1 0 0 0];
+
+internalField   uniform 0;
+
+boundaryField
+{
+    top
+    {
+        type            zeroGradient;
+    }
+
+    bottom
+    {
+        type            zeroGradient;
+    }
+
+    wall
+    {
+        // fixed heat loss, put all in gas phase
+        type            fixedValue;
+        value		uniform 0.0;
+    }
+
+    inner
+    {
+        type            zeroGradient;
+    }
+    inletLeft
+    {
+        type            zeroGradient;
+    }
+    inletRight
+    {
+        type            zeroGradient;
+    }
+}
+
+// ************************************************************************* //

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/CFD/0/orig.0/particleForces

@@ -0,0 +1,53 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  1.6                                   |
+|   \\  /    A nd           | Web:      http://www.OpenFOAM.org               |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volVectorField;
+    object      particleForces;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [2 -5 -2 0 0 0 0];
+
+internalField   uniform (0 0 0);
+
+boundaryField
+{
+    top
+    {
+        type            zeroGradient;
+    }
+
+    bottom
+    {
+        type            zeroGradient;
+    }
+
+    wall
+    {
+        type            zeroGradient;
+    }
+
+    inner
+    {
+        type            zeroGradient;
+    }
+    inletLeft
+    {
+        type            zeroGradient;
+    }
+    inletRight
+    {
+        type            zeroGradient;
+    }
+
+}
+
+// ************************************************************************* //

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/CFD/0/orig.0/rho

@@ -0,0 +1,56 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  1.6                                   |
+|   \\  /    A nd           | Web:      http://www.OpenFOAM.org               |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volScalarField;
+    object      rho;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [1 -3 0 0 0 0 0];
+
+internalField   uniform 1.1885;
+
+boundaryField
+{
+    top
+    {
+        type            zeroGradient;
+    }
+
+    bottom
+    {
+        type            zeroGradient;
+    }
+
+    wall
+    {
+        type            zeroGradient;
+    }
+
+    inner
+    {
+        type            zeroGradient;
+    }
+    inletLeft
+    {
+        type            fixedValue;
+        value		uniform 1.5;
+    }
+    inletRight
+    {
+        type            fixedValue;
+        value		uniform 1.5;
+    }
+
+}
+
+
+// ************************************************************************* //

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/CFD/0/orig.0/voidfraction

@@ -0,0 +1,57 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  1.6                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volScalarField;
+    location    "0";
+    object      voidfraction;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 0 0 0 0 0 0];
+
+internalField   uniform 1;
+
+boundaryField
+{
+    top
+    {
+        type            zeroGradient;
+    }
+
+    bottom
+    {
+        type            zeroGradient;
+    }
+
+    wall
+    {
+        type            zeroGradient;
+    }
+
+    inner
+    {
+        type            zeroGradient;
+    }
+    inletLeft
+    {
+        type            fixedValue;
+        value           uniform 0.4;
+    }
+    inletRight
+    {
+        type            fixedValue;
+        value           uniform 0.4;
+    }
+
+}
+
+
+// ************************************************************************* //

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/CFD/0/orig.0/voidfractionNext

@@ -0,0 +1,57 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  1.6                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volScalarField;
+    location    "0";
+    object      voidfractionNext;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 0 0 0 0 0 0];
+
+internalField   uniform 1;
+
+boundaryField
+{
+    top
+    {
+        type            zeroGradient;
+    }
+
+    bottom
+    {
+        type            zeroGradient;
+    }
+
+    wall
+    {
+        type            zeroGradient;
+    }
+
+    inner
+    {
+        type            zeroGradient;
+    }
+    inletLeft
+    {
+        type            fixedValue;
+        value           uniform 1;
+    }
+    inletRight
+    {
+        type            fixedValue;
+        value           uniform 1;
+    }
+
+}
+
+
+// ************************************************************************* //

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/CFD/0/orig.0/voidfractionPrev

@@ -0,0 +1,57 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  1.6                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volScalarField;
+    location    "0";
+    object      voidfractionPrev;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 0 0 0 0 0 0];
+
+internalField   uniform 1;
+
+boundaryField
+{
+    top
+    {
+        type            zeroGradient;
+    }
+
+    bottom
+    {
+        type            zeroGradient;
+    }
+
+    wall
+    {
+        type            zeroGradient;
+    }
+
+    inner
+    {
+        type            zeroGradient;
+    }
+    inletLeft
+    {
+        type            fixedValue;
+        value           uniform 1;
+    }
+    inletRight
+    {
+        type            fixedValue;
+        value           uniform 1;
+    }
+
+}
+
+
+// ************************************************************************* //

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/CFD/0/orig.0/wallHeatLoss

@@ -0,0 +1,65 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  1.6                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volScalarField;
+    location    "0";
+    object      wallHeatLoss;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [1 0 -3 -1 0 0 0];
+
+internalField   uniform 1;
+
+boundaryField
+{
+   bottom
+    {
+        type            fixedValue;
+        value           uniform 20.0;
+    }
+
+    
+    top
+    {
+        type            zeroGradient;
+
+    }
+
+    wall
+    {
+        // value according to Zhang et al. (2016)
+        type            fixedValue;
+        value           uniform 20.0;
+    }
+
+    inner
+    {
+        type            zeroGradient;
+    }
+
+    defaultFaces
+    {
+        type            empty;
+    }
+
+    inletLeft
+    {
+        type            zeroGradient;
+    }
+    inletRight
+    {
+        type            zeroGradient;
+    }
+
+}
+
+// ************************************************************************* //

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/CFD/0/p

@@ -0,0 +1,60 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  2.3.0                                 |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volScalarField;
+    object      p;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [1 -1 -2 0 0 0 0];
+
+internalField   uniform 1.0e5;
+
+boundaryField
+{
+    bottom
+    {
+        type             zeroGradient;
+    }
+
+    
+    top
+    {
+        type            fixedValue;
+        value           uniform 1.0e5;
+    }
+
+    wall
+    {
+        type            zeroGradient;
+    }
+
+    inner
+    {
+        type            zeroGradient;
+    }
+
+    defaultFaces
+    {
+        type            empty;
+    }
+
+    inletLeft
+    {
+        type            zeroGradient;
+    }
+    inletRight
+    {
+        type            zeroGradient;
+    }
+}
+
+// ************************************************************************* //

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/CFD/0/partEffThermCondField

@@ -0,0 +1,54 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  1.6                                   |
+|   \\  /    A nd           | Web:      http://www.OpenFOAM.org               |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volScalarField;
+    object      partEffThermCondField;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [1 1 -3 -1 0 0 0];
+
+internalField   uniform 0;
+
+boundaryField
+{
+    top
+    {
+        type            zeroGradient;
+    }
+
+    bottom
+    {
+        type            zeroGradient;
+    }
+
+    wall
+    {
+        // fixed heat loss, put all in gas phase
+        type            fixedValue;
+        value		uniform 0.0;
+    }
+
+    inner
+    {
+        type            zeroGradient;
+    }
+    inletLeft
+    {
+        type            zeroGradient;
+    }
+    inletRight
+    {
+        type            zeroGradient;
+    }
+}
+
+// ************************************************************************* //

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/CFD/0/particleForces

@@ -0,0 +1,53 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  1.6                                   |
+|   \\  /    A nd           | Web:      http://www.OpenFOAM.org               |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volVectorField;
+    object      particleForces;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [2 -5 -2 0 0 0 0];
+
+internalField   uniform (0 0 0);
+
+boundaryField
+{
+    top
+    {
+        type            zeroGradient;
+    }
+
+    bottom
+    {
+        type            zeroGradient;
+    }
+
+    wall
+    {
+        type            zeroGradient;
+    }
+
+    inner
+    {
+        type            zeroGradient;
+    }
+    inletLeft
+    {
+        type            zeroGradient;
+    }
+    inletRight
+    {
+        type            zeroGradient;
+    }
+
+}
+
+// ************************************************************************* //

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/CFD/0/rho

@@ -0,0 +1,56 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  1.6                                   |
+|   \\  /    A nd           | Web:      http://www.OpenFOAM.org               |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volScalarField;
+    object      rho;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [1 -3 0 0 0 0 0];
+
+internalField   uniform 1.1885;
+
+boundaryField
+{
+    top
+    {
+        type            zeroGradient;
+    }
+
+    bottom
+    {
+        type            zeroGradient;
+    }
+
+    wall
+    {
+        type            zeroGradient;
+    }
+
+    inner
+    {
+        type            zeroGradient;
+    }
+    inletLeft
+    {
+        type            fixedValue;
+        value		uniform 1.5;
+    }
+    inletRight
+    {
+        type            fixedValue;
+        value		uniform 1.5;
+    }
+
+}
+
+
+// ************************************************************************* //

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/CFD/0/voidfraction

@@ -0,0 +1,57 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  1.6                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volScalarField;
+    location    "0";
+    object      voidfraction;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 0 0 0 0 0 0];
+
+internalField   uniform 1;
+
+boundaryField
+{
+    top
+    {
+        type            zeroGradient;
+    }
+
+    bottom
+    {
+        type            zeroGradient;
+    }
+
+    wall
+    {
+        type            zeroGradient;
+    }
+
+    inner
+    {
+        type            zeroGradient;
+    }
+    inletLeft
+    {
+        type            fixedValue;
+        value           uniform 0.4;
+    }
+    inletRight
+    {
+        type            fixedValue;
+        value           uniform 0.4;
+    }
+
+}
+
+
+// ************************************************************************* //

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/CFD/0/voidfractionNext

@@ -0,0 +1,57 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  1.6                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volScalarField;
+    location    "0";
+    object      voidfractionNext;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 0 0 0 0 0 0];
+
+internalField   uniform 1;
+
+boundaryField
+{
+    top
+    {
+        type            zeroGradient;
+    }
+
+    bottom
+    {
+        type            zeroGradient;
+    }
+
+    wall
+    {
+        type            zeroGradient;
+    }
+
+    inner
+    {
+        type            zeroGradient;
+    }
+    inletLeft
+    {
+        type            fixedValue;
+        value           uniform 1;
+    }
+    inletRight
+    {
+        type            fixedValue;
+        value           uniform 1;
+    }
+
+}
+
+
+// ************************************************************************* //

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/CFD/0/voidfractionPrev

@@ -0,0 +1,57 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  1.6                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volScalarField;
+    location    "0";
+    object      voidfractionPrev;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 0 0 0 0 0 0];
+
+internalField   uniform 1;
+
+boundaryField
+{
+    top
+    {
+        type            zeroGradient;
+    }
+
+    bottom
+    {
+        type            zeroGradient;
+    }
+
+    wall
+    {
+        type            zeroGradient;
+    }
+
+    inner
+    {
+        type            zeroGradient;
+    }
+    inletLeft
+    {
+        type            fixedValue;
+        value           uniform 1;
+    }
+    inletRight
+    {
+        type            fixedValue;
+        value           uniform 1;
+    }
+
+}
+
+
+// ************************************************************************* //

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/CFD/0/wallHeatLoss

@@ -0,0 +1,65 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  1.6                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volScalarField;
+    location    "0";
+    object      wallHeatLoss;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [1 0 -3 -1 0 0 0];
+
+internalField   uniform 1;
+
+boundaryField
+{
+   bottom
+    {
+        type            fixedValue;
+        value           uniform 20.0;
+    }
+
+    
+    top
+    {
+        type            zeroGradient;
+
+    }
+
+    wall
+    {
+        // value according to Zhang et al. (2016)
+        type            fixedValue;
+        value           uniform 20.0;
+    }
+
+    inner
+    {
+        type            zeroGradient;
+    }
+
+    defaultFaces
+    {
+        type            empty;
+    }
+
+    inletLeft
+    {
+        type            zeroGradient;
+    }
+    inletRight
+    {
+        type            zeroGradient;
+    }
+
+}
+
+// ************************************************************************* //

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/CFD/Allclean.sh

@@ -0,0 +1,3 @@
+rm -r postProcessing
+rm -r proc*
+rm -r dynamicCode

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/CFD/clockData/0/timeEvalFull.txt

@@ -0,0 +1,19 @@
+Name	avgdeltaT	maxdeltaT	nOfRuns	level	parentNr	parentName
+Global	9.315003e+04	9.316150e+04	88300	0	0	
+Coupling	7.612455e+04	7.613103e+04	88300	1	1	Global
+LIGGGHTS	3.098020e+04	3.098410e+04	88300	2	2	Coupling
+LIGGGHTS_reallocArrays	6.624312e+02	6.864300e+02	88300	2	2	Coupling
+resetVolFields	6.060625e+01	6.236000e+01	88300	2	2	Coupling
+getDEMdata	2.395094e+03	2.421820e+03	88300	2	2	Coupling
+findCell	4.393824e+03	4.823750e+03	88300	2	2	Coupling
+setvoidFraction	1.877871e+04	2.019341e+04	88300	2	2	Coupling
+setAverages	1.967899e+03	3.430440e+03	88300	2	2	Coupling
+setForce	9.589802e+03	9.621400e+03	88300	2	2	Coupling
+setParticleForceField	2.630923e+03	2.674020e+03	88300	2	2	Coupling
+giveDEMdata	7.310400e+02	7.625200e+02	88300	2	2	Coupling
+interpolateEulerFields	1.636500e+02	4.082800e+02	88300	2	2	Coupling
+dumpDEMdata	8.875000e-02	1.200000e-01	88300	2	2	Coupling
+Flow	1.702489e+04	1.702987e+04	88300	1	1	Global
+calcEnergyContributions	3.632901e+03	3.654920e+03	88300	2	2	Coupling
+speciesExecute	5.375000e-02	1.300000e-01	88300	2	2	Coupling
+postFlow	6.246600e+02	6.281200e+02	264900	2	15	Flow

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/CFD/clockData/0/timeEval_0.txt

@@ -0,0 +1,20 @@
+Time Evaluation for Processor Nr.0
+Name	deltaT	nOfRuns	level	parentNr	parentName
+Global	9.312911e+04	88300	0	0	
+Coupling	7.610723e+04	88300	1	1	Global
+LIGGGHTS	3.096977e+04	88300	2	2	Coupling
+LIGGGHTS_reallocArrays	6.376200e+02	88300	2	2	Coupling
+resetVolFields	5.900000e+01	88300	2	2	Coupling
+getDEMdata	2.421820e+03	88300	2	2	Coupling
+findCell	3.464180e+03	88300	2	2	Coupling
+setvoidFraction	2.009698e+04	88300	2	2	Coupling
+setAverages	1.627500e+03	88300	2	2	Coupling
+setForce	9.551710e+03	88300	2	2	Coupling
+setParticleForceField	2.670290e+03	88300	2	2	Coupling
+giveDEMdata	7.172200e+02	88300	2	2	Coupling
+interpolateEulerFields	1.164400e+02	88300	2	2	Coupling
+dumpDEMdata	1.100000e-01	88300	2	2	Coupling
+Flow	1.702118e+04	88300	1	1	Global
+calcEnergyContributions	3.639770e+03	88300	2	2	Coupling
+speciesExecute	1.300000e-01	88300	2	2	Coupling
+postFlow	6.203700e+02	264900	2	15	Flow

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/CFD/clockData/0/timeEval_1.txt

@@ -0,0 +1,20 @@
+Time Evaluation for Processor Nr.1
+Name	deltaT	nOfRuns	level	parentNr	parentName
+Global	9.315669e+04	88300	0	0	
+Coupling	7.613023e+04	88300	1	1	Global
+LIGGGHTS	3.098080e+04	88300	2	2	Coupling
+LIGGGHTS_reallocArrays	6.402300e+02	88300	2	2	Coupling
+resetVolFields	5.931000e+01	88300	2	2	Coupling
+getDEMdata	2.419350e+03	88300	2	2	Coupling
+findCell	3.510330e+03	88300	2	2	Coupling
+setvoidFraction	2.019341e+04	88300	2	2	Coupling
+setAverages	1.504670e+03	88300	2	2	Coupling
+setForce	9.555540e+03	88300	2	2	Coupling
+setParticleForceField	2.674020e+03	88300	2	2	Coupling
+giveDEMdata	7.245400e+02	88300	2	2	Coupling
+interpolateEulerFields	1.009100e+02	88300	2	2	Coupling
+dumpDEMdata	8.000000e-02	88300	2	2	Coupling
+Flow	1.702593e+04	88300	1	1	Global
+calcEnergyContributions	3.628730e+03	88300	2	2	Coupling
+speciesExecute	6.000000e-02	88300	2	2	Coupling
+postFlow	6.259300e+02	264900	2	15	Flow

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/CFD/clockData/0/timeEval_2.txt

@@ -0,0 +1,20 @@
+Time Evaluation for Processor Nr.2
+Name	deltaT	nOfRuns	level	parentNr	parentName
+Global	9.315974e+04	88300	0	0	
+Coupling	7.612933e+04	88300	1	1	Global
+LIGGGHTS	3.098410e+04	88300	2	2	Coupling
+LIGGGHTS_reallocArrays	6.387900e+02	88300	2	2	Coupling
+resetVolFields	5.832000e+01	88300	2	2	Coupling
+getDEMdata	2.418620e+03	88300	2	2	Coupling
+findCell	4.723250e+03	88300	2	2	Coupling
+setvoidFraction	1.910096e+04	88300	2	2	Coupling
+setAverages	1.384060e+03	88300	2	2	Coupling
+setForce	9.577370e+03	88300	2	2	Coupling
+setParticleForceField	2.649620e+03	88300	2	2	Coupling
+giveDEMdata	7.446100e+02	88300	2	2	Coupling
+interpolateEulerFields	1.032900e+02	88300	2	2	Coupling
+dumpDEMdata	9.000000e-02	88300	2	2	Coupling
+Flow	1.702973e+04	88300	1	1	Global
+calcEnergyContributions	3.609680e+03	88300	2	2	Coupling
+speciesExecute	5.000000e-02	88300	2	2	Coupling
+postFlow	6.281200e+02	264900	2	15	Flow

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/CFD/clockData/0/timeEval_3.txt

@@ -0,0 +1,20 @@
+Time Evaluation for Processor Nr.3
+Name	deltaT	nOfRuns	level	parentNr	parentName
+Global	9.316150e+04	88300	0	0	
+Coupling	7.613103e+04	88300	1	1	Global
+LIGGGHTS	3.098362e+04	88300	2	2	Coupling
+LIGGGHTS_reallocArrays	6.423200e+02	88300	2	2	Coupling
+resetVolFields	6.202000e+01	88300	2	2	Coupling
+getDEMdata	2.415760e+03	88300	2	2	Coupling
+findCell	4.790430e+03	88300	2	2	Coupling
+setvoidFraction	1.902303e+04	88300	2	2	Coupling
+setAverages	1.394970e+03	88300	2	2	Coupling
+setForce	9.577260e+03	88300	2	2	Coupling
+setParticleForceField	2.655070e+03	88300	2	2	Coupling
+giveDEMdata	7.347000e+02	88300	2	2	Coupling
+interpolateEulerFields	1.005100e+02	88300	2	2	Coupling
+dumpDEMdata	7.000000e-02	88300	2	2	Coupling
+Flow	1.702987e+04	88300	1	1	Global
+calcEnergyContributions	3.613100e+03	88300	2	2	Coupling
+speciesExecute	6.000000e-02	88300	2	2	Coupling
+postFlow	6.272700e+02	264900	2	15	Flow

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/CFD/clockData/0/timeEval_4.txt

@@ -0,0 +1,20 @@
+Time Evaluation for Processor Nr.4
+Name	deltaT	nOfRuns	level	parentNr	parentName
+Global	9.314680e+04	88300	0	0	
+Coupling	7.612486e+04	88300	1	1	Global
+LIGGGHTS	3.098179e+04	88300	2	2	Coupling
+LIGGGHTS_reallocArrays	6.842300e+02	88300	2	2	Coupling
+resetVolFields	6.200000e+01	88300	2	2	Coupling
+getDEMdata	2.370950e+03	88300	2	2	Coupling
+findCell	4.683930e+03	88300	2	2	Coupling
+setvoidFraction	1.887134e+04	88300	2	2	Coupling
+setAverages	1.634340e+03	88300	2	2	Coupling
+setForce	9.604220e+03	88300	2	2	Coupling
+setParticleForceField	2.619940e+03	88300	2	2	Coupling
+giveDEMdata	7.129900e+02	88300	2	2	Coupling
+interpolateEulerFields	1.188100e+02	88300	2	2	Coupling
+dumpDEMdata	7.000000e-02	88300	2	2	Coupling
+Flow	1.702148e+04	88300	1	1	Global
+calcEnergyContributions	3.643300e+03	88300	2	2	Coupling
+speciesExecute	6.000000e-02	88300	2	2	Coupling
+postFlow	6.236400e+02	264900	2	15	Flow

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/CFD/clockData/0/timeEval_5.txt

@@ -0,0 +1,20 @@
+Time Evaluation for Processor Nr.5
+Name	deltaT	nOfRuns	level	parentNr	parentName
+Global	9.314844e+04	88300	0	0	
+Coupling	7.612639e+04	88300	1	1	Global
+LIGGGHTS	3.097941e+04	88300	2	2	Coupling
+LIGGGHTS_reallocArrays	6.864300e+02	88300	2	2	Coupling
+resetVolFields	6.236000e+01	88300	2	2	Coupling
+getDEMdata	2.371250e+03	88300	2	2	Coupling
+findCell	4.823750e+03	88300	2	2	Coupling
+setvoidFraction	1.887920e+04	88300	2	2	Coupling
+setAverages	1.502390e+03	88300	2	2	Coupling
+setForce	9.621400e+03	88300	2	2	Coupling
+setParticleForceField	2.604740e+03	88300	2	2	Coupling
+giveDEMdata	7.001900e+02	88300	2	2	Coupling
+interpolateEulerFields	1.029400e+02	88300	2	2	Coupling
+dumpDEMdata	7.000000e-02	88300	2	2	Coupling
+Flow	1.702153e+04	88300	1	1	Global
+calcEnergyContributions	3.654920e+03	88300	2	2	Coupling
+speciesExecute	4.000000e-02	88300	2	2	Coupling
+postFlow	6.215200e+02	264900	2	15	Flow

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/CFD/clockData/0/timeEval_6.txt

@@ -0,0 +1,20 @@
+Time Evaluation for Processor Nr.6
+Name	deltaT	nOfRuns	level	parentNr	parentName
+Global	9.315006e+04	88300	0	0	
+Coupling	7.612536e+04	88300	1	1	Global
+LIGGGHTS	3.098098e+04	88300	2	2	Coupling
+LIGGGHTS_reallocArrays	6.853600e+02	88300	2	2	Coupling
+resetVolFields	6.005000e+01	88300	2	2	Coupling
+getDEMdata	2.373050e+03	88300	2	2	Coupling
+findCell	4.533870e+03	88300	2	2	Coupling
+setvoidFraction	1.691624e+04	88300	2	2	Coupling
+setAverages	3.430440e+03	88300	2	2	Coupling
+setForce	9.619720e+03	88300	2	2	Coupling
+setParticleForceField	2.578760e+03	88300	2	2	Coupling
+giveDEMdata	7.625200e+02	88300	2	2	Coupling
+interpolateEulerFields	4.082800e+02	88300	2	2	Coupling
+dumpDEMdata	1.000000e-01	88300	2	2	Coupling
+Flow	1.702413e+04	88300	1	1	Global
+calcEnergyContributions	3.636580e+03	88300	2	2	Coupling
+speciesExecute	1.000000e-02	88300	2	2	Coupling
+postFlow	6.241300e+02	264900	2	15	Flow

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/CFD/clockData/0/timeEval_7.txt

@@ -0,0 +1,20 @@
+Time Evaluation for Processor Nr.7
+Name	deltaT	nOfRuns	level	parentNr	parentName
+Global	9.314787e+04	88300	0	0	
+Coupling	7.612196e+04	88300	1	1	Global
+LIGGGHTS	3.098116e+04	88300	2	2	Coupling
+LIGGGHTS_reallocArrays	6.844700e+02	88300	2	2	Coupling
+resetVolFields	6.179000e+01	88300	2	2	Coupling
+getDEMdata	2.369950e+03	88300	2	2	Coupling
+findCell	4.620850e+03	88300	2	2	Coupling
+setvoidFraction	1.714848e+04	88300	2	2	Coupling
+setAverages	3.264820e+03	88300	2	2	Coupling
+setForce	9.611200e+03	88300	2	2	Coupling
+setParticleForceField	2.594940e+03	88300	2	2	Coupling
+giveDEMdata	7.515500e+02	88300	2	2	Coupling
+interpolateEulerFields	2.580200e+02	88300	2	2	Coupling
+dumpDEMdata	1.200000e-01	88300	2	2	Coupling
+Flow	1.702528e+04	88300	1	1	Global
+calcEnergyContributions	3.637130e+03	88300	2	2	Coupling
+speciesExecute	2.000000e-02	88300	2	2	Coupling
+postFlow	6.263000e+02	264900	2	15	Flow

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/CFD/constant/RASProperties

@@ -0,0 +1,25 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  2.3.0                                 |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       dictionary;
+    location    "constant";
+    object      RASProperties;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+RASModel        laminar;//kEpsilon;
+
+turbulence      off;
+
+printCoeffs     off;
+
+
+// ************************************************************************* //

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/CFD/constant/couplingProperties

@@ -0,0 +1,176 @@
+/*---------------------------------------------------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  1.4                                   |
+|   \\  /    A nd           | Web:      http://www.openfoam.org               |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+
+
+FoamFile
+{
+    version         2.0;
+    format          ascii;
+
+    root            "";
+    case            "";
+    instance        "";
+    local           "";
+
+    class           dictionary;
+    object          couplingProperties;
+}
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+//===========================================================================//
+// sub-models & settings
+
+modelType "A";
+
+couplingInterval 10;
+
+voidFractionModel divided;
+
+locateModel engine;
+
+meshMotionModel noMeshMotion;
+
+regionModel allRegion;
+
+IOModel "off";
+
+dataExchangeModel twoWayMPI;
+
+averagingModel dense;
+
+clockModel standardClock;
+
+smoothingModel off;
+
+probeModel off;
+
+forceModels
+(
+    particleDeformation
+    dSauter
+    BeetstraDragPoly
+    gradPForce
+    viscForce
+);
+
+momCoupleModels
+(
+    implicitCouple
+);
+
+energyModels
+(
+    heatTransferRanzMarshall
+);
+
+thermCondModel SyamlalThermCond;
+
+chemistryModels ();
+
+turbulenceModelType "turbulenceProperties";
+
+getParticleTypes true;
+getParticleDensities false;
+getParticleEffVolFactors true;
+//===========================================================================//
+// sub-model properties
+
+implicitCoupleProps
+{
+    velFieldName "U";
+    granVelFieldName "Us";
+    voidfractionFieldName "voidfraction";
+}
+
+gradPForceProps
+{
+    pFieldName "p";
+    densityFieldName "rho";
+    voidfractionFieldName "voidfraction";
+    velocityFieldName "U";
+}
+
+viscForceProps
+{
+    velocityFieldName "U";
+    densityFieldName "rho";
+}
+
+
+BeetstraDragProps
+{
+    velFieldName "U";
+    granVelFieldName "Us";
+    densityFieldName "rho";
+    voidfractionFieldName "voidfraction";
+    phi 1;
+    coarseGrainingFactors (4.0 6.0 10.0 20.0);
+}
+
+dSauterProps
+{
+    coarseGrainingFactors (4.0 6.0 10.0 20.0);
+}
+
+particleDeformationProps
+{
+    refFieldName "partTempRef";
+    partTypes (3 4);
+    defaultDeformCellsName "noOreReg";
+// same as Huo et al. (2017)
+    lowerBounds (1475 1475);
+    upperBounds (1675 1675);
+    backgroundVoidage 0.05;
+}
+
+heatTransferRanzMarshallProps
+{
+    partTempName "Temp";
+    partHeatFluxName "convectiveHeatFlux";
+    calcPartTempField false;
+    coarseGrainingFactors (4.0 6.0 10.0 20.0);
+    partRefTemp 293;
+    interpolation false;
+    NusseltScalingFactor 0.1;
+    initPartTemp true;
+    Tmin 293;
+    Tmax 3000;
+}
+
+SyamlalThermCondProps
+{
+}
+
+
+centreProps
+{
+    alphaMin 0.10;
+}
+
+engineProps
+{
+    treeSearch true;
+}
+
+dividedProps
+{
+    alphaMin 0.1;
+    scaleUpVol 1.0;
+    weight 1.0;
+    verbose;
+}
+
+twoWayMPIProps
+{
+    liggghtsPath "../DEM/in.liggghts_CFDDEM";
+    liggghtsEndOfRunPath "../DEM/in.final";
+}
+
+
+// ************************************************************************* //

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/CFD/constant/g

@@ -0,0 +1,22 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  1.6                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       uniformDimensionedVectorField;
+    location    "constant";
+    object      g;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 1 -2 0 0 0 0];
+value           ( 0 0 -9.81 );
+
+
+// ************************************************************************* //

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/CFD/constant/liggghtsCommands

@@ -0,0 +1,42 @@
+/*---------------------------------------------------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  1.4                                   |
+|   \\  /    A nd           | Web:      http://www.openfoam.org               |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+
+
+FoamFile
+{
+    version         2.0;
+    format          ascii;
+
+    root            "";
+    case            "";
+    instance        "";
+    local           "";
+
+    class           dictionary;
+    object          liggghtsCommands;
+}
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+liggghtsCommandModels
+(
+   runLiggghts
+   writeLiggghts
+);
+// ************************************************************************* //
+
+/*runLiggghtsProps
+{
+    preNo false;
+}
+
+writeLiggghtsProps
+{
+    writeLast off;
+    overwrite off;
+}*/

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/CFD/constant/symmetrizerPropertiesPost

@@ -0,0 +1,24 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  2.3.0                                 |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       dictionary;
+    location    "constant";
+    object      symmetrizerProperties;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+scalarFields (voidfractionMean);
+vectorFields (UsMean UMean);
+
+refPoint (0 0 0);
+refDirection (1 0 0);
+
+// ************************************************************************* //

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/CFD/constant/symmetrizerPropertiesPre

@@ -0,0 +1,24 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  2.3.0                                 |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       dictionary;
+    location    "constant";
+    object      symmetrizerProperties;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+scalarFields (partTempRef);
+vectorFields ();
+
+refPoint (0 0 0);
+refDirection (1 0 0);
+
+// ************************************************************************* //

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/CFD/constant/thermophysicalProperties

@@ -0,0 +1,49 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  2.3.0                                 |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       dictionary;
+    location    "constant";
+    object      thermophysicalProperties;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+thermoType
+{
+    type            hePsiThermo;
+    mixture         pureMixture;
+    transport       const;
+    thermo          eConst;
+    equationOfState perfectGas;
+    specie          specie;
+    energy          sensibleInternalEnergy;
+}
+
+mixture
+{
+    specie
+    {
+        nMoles      1;
+        molWeight   28.9;
+    }
+    thermodynamics
+    {
+        Cv          718;
+        Hf          0;
+    }
+    transport
+    {
+        mu          2e-05;
+        Pr          0.8;
+    }
+}
+
+
+// ************************************************************************* //

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/CFD/constant/transportProperties

@@ -0,0 +1,43 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  2.3.0                                 |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       dictionary;
+    location    "constant";
+    object      transportProperties;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+transportModel  Newtonian;
+
+nu              nu [ 0 2 -1 0 0 0 0 ] 1e-05;
+
+kf		kf [ 1 1 -3 -1 0 0 0 ] 0.026;
+
+Cp		Cp [ 0 2 -2 -1 0 0 0 ] 1005;
+
+CrossPowerLawCoeffs
+{
+    nu0             nu0 [ 0 2 -1 0 0 0 0 ] 1e-06;
+    nuInf           nuInf [ 0 2 -1 0 0 0 0 ] 1e-06;
+    m               m [ 0 0 1 0 0 0 0 ] 1;
+    n               n [ 0 0 0 0 0 0 0 ] 1;
+}
+
+BirdCarreauCoeffs
+{
+    nu0             nu0 [ 0 2 -1 0 0 0 0 ] 1e-06;
+    nuInf           nuInf [ 0 2 -1 0 0 0 0 ] 1e-06;
+    k               k [ 0 0 1 0 0 0 0 ] 0;
+    n               n [ 0 0 0 0 0 0 0 ] 1;
+}
+
+
+// ************************************************************************* //

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/CFD/constant/turbulenceProperties

@@ -0,0 +1,21 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  2.3.0                                 |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       dictionary;
+    location    "constant";
+    object      turbulenceProperties;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+simulationType  laminar;
+
+
+// ************************************************************************* //

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/CFD/mesh_slot.sh

@@ -0,0 +1,13 @@
+./Allclean.sh
+cd system
+m4 -P blockMeshDict_slot.m4 > blockMeshDict
+cd ..
+blockMesh
+mirrorMesh -overwrite
+
+topoSet -dict system/topoSetDict
+createPatch -overwrite
+#rm -r 0
+#cp -r orig.0/ 0
+
+setsToZones -noFlipMap

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/CFD/orig.0/Ksl

@@ -0,0 +1,54 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  1.6                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volScalarField;
+    location    "0";
+    object      Ksl;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [1 -3 -1 0 0 0 0];
+
+internalField   uniform 0;
+
+boundaryField
+{
+    top
+    {
+        type            zeroGradient;
+    }
+
+    bottom
+    {
+        type            zeroGradient;
+    }
+
+    wall
+    {
+        type            zeroGradient;
+    }
+
+    inner
+    {
+        type            zeroGradient;
+    }
+    inletLeft
+    {
+        type            zeroGradient;
+    }
+    inletRight
+    {
+        type            zeroGradient;
+    }
+}
+
+
+// ************************************************************************* //

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/CFD/orig.0/KslNext

@@ -0,0 +1,53 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  1.6                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volScalarField;
+    location    "0";
+    object      KslNext;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [1 -3 -1 0 0 0 0];
+
+internalField   uniform 0;
+
+boundaryField
+{
+    top
+    {
+        type            zeroGradient;
+    }
+
+    bottom
+    {
+        type            zeroGradient;
+    }
+
+    wall
+    {
+        type            zeroGradient;
+    }
+
+    inner
+    {
+        type            zeroGradient;
+    }
+    inletLeft
+    {
+        type            zeroGradient;
+    }
+    inletRight
+    {
+        type            zeroGradient;
+    }
+}
+
+// ************************************************************************* //

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/CFD/orig.0/KslPrev

@@ -0,0 +1,54 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  1.6                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volScalarField;
+    location    "0";
+    object      KslPrev;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [1 -3 -1 0 0 0 0];
+
+internalField   uniform 0;
+
+boundaryField
+{
+    top
+    {
+        type            zeroGradient;
+    }
+
+    bottom
+    {
+        type            zeroGradient;
+    }
+
+    wall
+    {
+        type            zeroGradient;
+    }
+
+    inner
+    {
+        type            zeroGradient;
+    }
+    inletLeft
+    {
+        type            zeroGradient;
+    }
+    inletRight
+    {
+        type            zeroGradient;
+    }
+}
+
+
+// ************************************************************************* //

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/CFD/orig.0/T

@@ -0,0 +1,56 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  1.6                                   |
+|   \\  /    A nd           | Web:      http://www.OpenFOAM.org               |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volScalarField;
+    object      T;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 0 0 1 0 0 0];
+
+internalField   uniform 293;
+
+boundaryField
+{
+    top
+    {
+        type            zeroGradient;
+    }
+
+    bottom
+    {
+        type            fixedValue;
+        value		uniform 293;
+    }
+
+    wall
+    {
+        type            fixedValue;
+        value		uniform 293;
+    }
+
+    inner
+    {
+        type            zeroGradient;
+    }
+    inletLeft
+    {
+        type            fixedValue;
+        value		uniform 2600;
+    }
+    inletRight
+    {
+        type            fixedValue;
+        value		uniform 2600;
+    }
+}
+
+// ************************************************************************* //

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/CFD/orig.0/U

@@ -0,0 +1,66 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  2.3.0                                 |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volVectorField;
+    object      U;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 1 -1 0 0 0 0];
+
+internalField   uniform (0 0 0);
+
+boundaryField
+{
+    bottom
+    {
+        type            uniformFixedValue;
+        uniformValue    constant (0 0 0);
+    }
+
+   
+    top
+    {
+        type            zeroGradient;
+    }
+
+    wall
+    {
+        type            uniformFixedValue;
+        uniformValue    constant (0 0 0);
+    }
+
+    inner
+    {
+        type            uniformFixedValue;
+        uniformValue    constant (0 0 0);
+    }
+
+    defaultFaces
+    {
+        type            empty;
+    }
+
+    inletLeft
+    {
+        type            uniformFixedValue;
+        uniformValue    constant (200 0 0);
+    }
+    inletRight
+    {
+        type            uniformFixedValue;
+        uniformValue    constant (-200 0 0);
+    }
+
+
+}
+
+// ************************************************************************* //

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/CFD/orig.0/Us

@@ -0,0 +1,54 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  1.6                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volVectorField;
+    location    "0";
+    object      Us;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 1 -1 0 0 0 0];
+
+internalField   uniform (0 0 0);
+
+boundaryField
+{
+    top
+    {
+        type            zeroGradient;
+    }
+
+    bottom
+    {
+        type            zeroGradient;
+    }
+
+    wall
+    {
+        type            zeroGradient;
+    }
+
+    inner
+    {
+        type            zeroGradient;
+    }
+    inletLeft
+    {
+        type            zeroGradient;
+    }
+    inletRight
+    {
+        type            zeroGradient;
+    }
+}
+
+
+// ************************************************************************* //

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/CFD/orig.0/UsNext

@@ -0,0 +1,55 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  1.6                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volVectorField;
+    location    "0";
+    object      UsNext;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 1 -1 0 0 0 0];
+
+internalField   uniform (0 0 0);
+
+boundaryField
+{
+    top
+    {
+        type            zeroGradient;
+    }
+
+    bottom
+    {
+        type            zeroGradient;
+    }
+
+    wall
+    {
+        type            zeroGradient;
+    }
+
+    inner
+    {
+        type            zeroGradient;
+    }
+    inletLeft
+    {
+        type            zeroGradient;
+    }
+    inletRight
+    {
+        type            zeroGradient;
+    }
+
+}
+
+
+// ************************************************************************* //

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/CFD/orig.0/UsPrev

@@ -0,0 +1,55 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  1.6                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volVectorField;
+    location    "0";
+    object      UsPrev;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 1 -1 0 0 0 0];
+
+internalField   uniform (0 0 0);
+
+boundaryField
+{
+    top
+    {
+        type            zeroGradient;
+    }
+
+    bottom
+    {
+        type            zeroGradient;
+    }
+
+    wall
+    {
+        type            zeroGradient;
+    }
+
+    inner
+    {
+        type            zeroGradient;
+    }
+    inletLeft
+    {
+        type            zeroGradient;
+    }
+    inletRight
+    {
+        type            zeroGradient;
+    }
+
+
+}
+
+// ************************************************************************* //

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/CFD/orig.0/_wallQFactor

@@ -0,0 +1,63 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  1.6                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volScalarField;
+    location    "0";
+    object      voidfraction;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 0 0 0 0 0 0];
+
+internalField   uniform 1;
+
+boundaryField
+{
+   bottom
+    {
+        type            zeroGradient;
+    }
+
+    
+    top
+    {
+        type            zeroGradient;
+
+    }
+
+    wall
+    {
+        type            fixedValue;
+        value           uniform 10.0;
+    }
+
+    inner
+    {
+        type            zeroGradient;
+    }
+
+    defaultFaces
+    {
+        type            empty;
+    }
+
+    inletLeft
+    {
+        type            zeroGradient;
+    }
+    inletRight
+    {
+        type            zeroGradient;
+    }
+
+}
+
+// ************************************************************************* //

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/CFD/orig.0/addSource

@@ -0,0 +1,53 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  1.6                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volScalarField;
+    object      addSource;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [1 -1 -3 0 0 0 0]; 
+
+internalField   uniform 0;  
+
+boundaryField
+{
+    top
+    {
+        type            zeroGradient;
+    }
+
+    bottom
+    {
+        type            zeroGradient;
+    }
+
+    wall
+    {
+        type            zeroGradient;
+    }
+
+    inner
+    {
+        type            zeroGradient;
+    }
+    inletLeft
+    {
+        type            zeroGradient;
+    }
+    inletRight
+    {
+        type            zeroGradient;
+    }
+}
+
+
+// ************************************************************************* //

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/CFD/orig.0/p

@@ -0,0 +1,60 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  2.3.0                                 |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volScalarField;
+    object      p;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [1 -1 -2 0 0 0 0];
+
+internalField   uniform 1.0e5;
+
+boundaryField
+{
+    bottom
+    {
+        type             zeroGradient;
+    }
+
+    
+    top
+    {
+        type            fixedValue;
+        value           uniform 1.0e5;
+    }
+
+    wall
+    {
+        type            zeroGradient;
+    }
+
+    inner
+    {
+        type            zeroGradient;
+    }
+
+    defaultFaces
+    {
+        type            empty;
+    }
+
+    inletLeft
+    {
+        type            zeroGradient;
+    }
+    inletRight
+    {
+        type            zeroGradient;
+    }
+}
+
+// ************************************************************************* //

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/CFD/orig.0/partEffThermCondField

@@ -0,0 +1,54 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  1.6                                   |
+|   \\  /    A nd           | Web:      http://www.OpenFOAM.org               |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volScalarField;
+    object      partEffThermCondField;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [1 1 -3 -1 0 0 0];
+
+internalField   uniform 0;
+
+boundaryField
+{
+    top
+    {
+        type            zeroGradient;
+    }
+
+    bottom
+    {
+        type            zeroGradient;
+    }
+
+    wall
+    {
+        // fixed heat loss, put all in gas phase
+        type            fixedValue;
+        value		uniform 0.0;
+    }
+
+    inner
+    {
+        type            zeroGradient;
+    }
+    inletLeft
+    {
+        type            zeroGradient;
+    }
+    inletRight
+    {
+        type            zeroGradient;
+    }
+}
+
+// ************************************************************************* //

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/CFD/orig.0/particleForces

@@ -0,0 +1,53 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  1.6                                   |
+|   \\  /    A nd           | Web:      http://www.OpenFOAM.org               |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volVectorField;
+    object      particleForces;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [2 -5 -2 0 0 0 0];
+
+internalField   uniform (0 0 0);
+
+boundaryField
+{
+    top
+    {
+        type            zeroGradient;
+    }
+
+    bottom
+    {
+        type            zeroGradient;
+    }
+
+    wall
+    {
+        type            zeroGradient;
+    }
+
+    inner
+    {
+        type            zeroGradient;
+    }
+    inletLeft
+    {
+        type            zeroGradient;
+    }
+    inletRight
+    {
+        type            zeroGradient;
+    }
+
+}
+
+// ************************************************************************* //

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/CFD/orig.0/rho

@@ -0,0 +1,56 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  1.6                                   |
+|   \\  /    A nd           | Web:      http://www.OpenFOAM.org               |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volScalarField;
+    object      rho;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [1 -3 0 0 0 0 0];
+
+internalField   uniform 1.1885;
+
+boundaryField
+{
+    top
+    {
+        type            zeroGradient;
+    }
+
+    bottom
+    {
+        type            zeroGradient;
+    }
+
+    wall
+    {
+        type            zeroGradient;
+    }
+
+    inner
+    {
+        type            zeroGradient;
+    }
+    inletLeft
+    {
+        type            fixedValue;
+        value		uniform 1.5;
+    }
+    inletRight
+    {
+        type            fixedValue;
+        value		uniform 1.5;
+    }
+
+}
+
+
+// ************************************************************************* //

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/CFD/orig.0/voidfraction

@@ -0,0 +1,57 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  1.6                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volScalarField;
+    location    "0";
+    object      voidfraction;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 0 0 0 0 0 0];
+
+internalField   uniform 1;
+
+boundaryField
+{
+    top
+    {
+        type            zeroGradient;
+    }
+
+    bottom
+    {
+        type            zeroGradient;
+    }
+
+    wall
+    {
+        type            zeroGradient;
+    }
+
+    inner
+    {
+        type            zeroGradient;
+    }
+    inletLeft
+    {
+        type            fixedValue;
+        value           uniform 0.4;
+    }
+    inletRight
+    {
+        type            fixedValue;
+        value           uniform 0.4;
+    }
+
+}
+
+
+// ************************************************************************* //

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/CFD/orig.0/voidfractionNext

@@ -0,0 +1,57 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  1.6                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volScalarField;
+    location    "0";
+    object      voidfractionNext;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 0 0 0 0 0 0];
+
+internalField   uniform 1;
+
+boundaryField
+{
+    top
+    {
+        type            zeroGradient;
+    }
+
+    bottom
+    {
+        type            zeroGradient;
+    }
+
+    wall
+    {
+        type            zeroGradient;
+    }
+
+    inner
+    {
+        type            zeroGradient;
+    }
+    inletLeft
+    {
+        type            fixedValue;
+        value           uniform 1;
+    }
+    inletRight
+    {
+        type            fixedValue;
+        value           uniform 1;
+    }
+
+}
+
+
+// ************************************************************************* //

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/CFD/orig.0/voidfractionPrev

@@ -0,0 +1,57 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  1.6                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volScalarField;
+    location    "0";
+    object      voidfractionPrev;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 0 0 0 0 0 0];
+
+internalField   uniform 1;
+
+boundaryField
+{
+    top
+    {
+        type            zeroGradient;
+    }
+
+    bottom
+    {
+        type            zeroGradient;
+    }
+
+    wall
+    {
+        type            zeroGradient;
+    }
+
+    inner
+    {
+        type            zeroGradient;
+    }
+    inletLeft
+    {
+        type            fixedValue;
+        value           uniform 1;
+    }
+    inletRight
+    {
+        type            fixedValue;
+        value           uniform 1;
+    }
+
+}
+
+
+// ************************************************************************* //

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/CFD/orig.0/wallHeatLoss

@@ -0,0 +1,65 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  1.6                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volScalarField;
+    location    "0";
+    object      wallHeatLoss;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [1 0 -3 -1 0 0 0];
+
+internalField   uniform 1;
+
+boundaryField
+{
+   bottom
+    {
+        type            fixedValue;
+        value           uniform 20.0;
+    }
+
+    
+    top
+    {
+        type            zeroGradient;
+
+    }
+
+    wall
+    {
+        // value according to Zhang et al. (2016)
+        type            fixedValue;
+        value           uniform 20.0;
+    }
+
+    inner
+    {
+        type            zeroGradient;
+    }
+
+    defaultFaces
+    {
+        type            empty;
+    }
+
+    inletLeft
+    {
+        type            zeroGradient;
+    }
+    inletRight
+    {
+        type            zeroGradient;
+    }
+
+}
+
+// ************************************************************************* //

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/CFD/system/blockMeshDict

@@ -0,0 +1,195 @@
+// Parametrized test case for a BF geometry
+
+
+//Run using:
+//m4 -P blockMeshDict.m4 > blockMeshDict
+
+//m4 definitions:
+
+
+
+
+
+//Mathematical constants:
+
+
+//Geometry
+
+// width of wedge
+
+
+
+// height levels
+
+
+
+
+
+
+
+// xlevels
+
+
+
+
+
+
+
+
+//Grid points (integers!):
+
+
+
+
+
+
+
+
+
+//m4_define(rGrading, 0.5)
+
+
+
+
+/*---------------------------------------------------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  1.4.1                                 |
+|   \\  /    A nd           | Web:      http://www.openfoam.org               |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+
+FoamFile
+{
+    version         2.0;
+    format          ascii;
+
+    root            "";
+    case            "";
+    instance        "";
+    local           "";
+
+    class           dictionary;
+    object          blockMeshDict;
+}
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+convertToMeters 1;
+
+vertices
+(
+(-6 0.6 -2) // Vertex V0 = 0 
+(0 0.6 -2) // Vertex V1 = 1 
+(0 0.6 0.75) // Vertex V2 = 2 
+(-6 0.6 0.75) // Vertex V3 = 3 
+
+(-7.4 0.6 3.9) // Vertex V4 = 4 
+(0 0.6 3.9) // Vertex V5 = 5 
+(0 0.6 6.4) // Vertex V6 = 6 
+(-7.4 0.6 6.4) // Vertex V7 = 7 
+
+(-4.7 0.6 24) // Vertex V8 = 8 
+(0 0.6 24) // Vertex V9 = 9 
+(0 0.6 26) // Vertex V10 = 10 
+(-4.7 0.6 26) // Vertex V11 = 11 
+
+// neg. y values
+(-6 -0.6 -2) // Vertex V12 = 12 
+(0 -0.6 -2) // Vertex V13 = 13 
+(0 -0.6 0.75) // Vertex V14 = 14 
+(-6 -0.6 0.75) // Vertex V15 = 15 
+
+(-7.4 -0.6 3.9) // Vertex V16 = 16 
+(0 -0.6 3.9) // Vertex V17 = 17 
+(0 -0.6 6.4) // Vertex V18 = 18 
+(-7.4 -0.6 6.4) // Vertex V19 = 19 
+
+(-4.7 -0.6 24) // Vertex V20 = 20 
+(0 -0.6 24) // Vertex V21 = 21 
+(0 -0.6 26) // Vertex V22 = 22 
+(-4.7 -0.6 26) // Vertex V23 = 23 
+
+
+);
+
+// Defining blocks:
+blocks
+(
+    hex ( 0 1 2 3 12 13 14 15 ) AB
+    (32 25 11)
+    simpleGrading (1.0 1 1)
+
+    hex ( 3 2 5 4 15 14 17 16) BC
+    (32 15 11)
+    simpleGrading (1.0 1 1)
+
+    hex ( 4 5 6 7 16 17 18 19) CD
+    (32 12 11)
+    simpleGrading (1.0 1 1)
+
+    hex ( 7 6 9 8 19 18 21 20 ) EF
+    (32 90 11)
+    simpleGrading (1.0 1 1)
+
+    hex ( 8 9 10 11 20 21 22 23 ) GH
+    (32 10 11)
+    simpleGrading (1.0 1 1)
+);
+
+// Defining patches:
+boundary
+(
+    bottom
+    {
+        type patch;
+        faces
+        (
+         (0 1 13 12)
+        );
+    }
+    wall
+    {
+        type wall;
+        faces
+        (
+         (0 3 15 12)
+         (3 4 16 15)
+         (4 7 19 16)
+         (7 8 20 19)
+         (8 11 23 20)
+        );
+    }
+    top
+    {
+        type patch;
+        faces
+        (
+         (11 10 22 23)
+        );
+    }
+    inner
+    { 
+        type patch;
+        faces
+        (
+         (0 1 2 3)
+         (3 2 5 4)
+         (4 5 6 7)
+         (7 6 9 8)
+         (8 9 10 11)
+
+         (12 13 14 15)
+         (15 14 17 16)
+         (16 17 18 19)
+         (19 18 21 20)
+         (20 21 22 23)
+        );
+    }
+);
+
+mergePatchPairs 
+(
+);
+
+// ************************************************************************* //

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/CFD/system/blockMeshDict_slot.m4

@@ -0,0 +1,195 @@
+// Parametrized test case for a BF geometry
+
+
+//Run using:
+//m4 -P blockMeshDict.m4 > blockMeshDict
+
+//m4 definitions:
+m4_changecom(//)m4_changequote([,])
+m4_define(calc, [m4_esyscmd(perl -e 'use Math::Trig; printf ($1)')])
+m4_define(VCOUNT, 0)
+m4_define(vlabel, [[// ]Vertex $1 = VCOUNT m4_define($1, VCOUNT)m4_define([VCOUNT], m4_incr(VCOUNT))])
+
+//Mathematical constants:
+m4_define(pi, 3.1415926536)
+
+//Geometry
+
+// width of wedge
+m4_define(y0, 0.6)
+m4_define(y1, -0.6)
+
+// height levels
+m4_define(z0, -2)
+m4_define(z1, 0.75)
+m4_define(z2, 3.9)
+m4_define(z3, 6.4)
+m4_define(z4, 24)
+m4_define(z5, 26)
+
+// xlevels
+m4_define(x0, -6)
+m4_define(x1, -6)
+m4_define(x2, -7.4)
+m4_define(x3, -7.4)
+m4_define(x4, -4.7)
+m4_define(x5, -4.7)
+
+
+//Grid points (integers!):
+
+m4_define(xNumberOfCells, 32)
+m4_define(yNumberOfCells, 11)
+m4_define(zNumberOfCells0, 25)
+m4_define(zNumberOfCells1, 15)
+m4_define(zNumberOfCells2, 12)
+m4_define(zNumberOfCells3, 90)
+m4_define(zNumberOfCells4, 10)
+m4_define(rGrading, 1.0)
+//m4_define(rGrading, 0.5)
+
+
+
+
+/*---------------------------------------------------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  1.4.1                                 |
+|   \\  /    A nd           | Web:      http://www.openfoam.org               |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+
+FoamFile
+{
+    version         2.0;
+    format          ascii;
+
+    root            "";
+    case            "";
+    instance        "";
+    local           "";
+
+    class           dictionary;
+    object          blockMeshDict;
+}
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+convertToMeters 1;
+
+vertices
+(
+(x0 y0 z0) vlabel(V0)
+(0 y0 z0) vlabel(V1)
+(0 y0 z1) vlabel(V2)
+(x1 y0 z1) vlabel(V3)
+
+(x2 y0 z2) vlabel(V4)
+(0 y0 z2) vlabel(V5)
+(0 y0 z3) vlabel(V6)
+(x3 y0 z3) vlabel(V7)
+
+(x4 y0 z4) vlabel(V8)
+(0 y0 z4) vlabel(V9)
+(0 y0 z5) vlabel(V10)
+(x5 y0 z5) vlabel(V11)
+
+// neg. y values
+(x0 y1 z0) vlabel(V12)
+(0 y1 z0) vlabel(V13)
+(0 y1 z1) vlabel(V14)
+(x1 y1 z1) vlabel(V15)
+
+(x2 y1 z2) vlabel(V16)
+(0 y1 z2) vlabel(V17)
+(0 y1 z3) vlabel(V18)
+(x3 y1 z3) vlabel(V19)
+
+(x4 y1 z4) vlabel(V20)
+(0 y1 z4) vlabel(V21)
+(0 y1 z5) vlabel(V22)
+(x5 y1 z5) vlabel(V23)
+
+
+);
+
+// Defining blocks:
+blocks
+(
+    hex ( V0 V1 V2 V3 V12 V13 V14 V15 ) AB
+    (xNumberOfCells zNumberOfCells0 yNumberOfCells)
+    simpleGrading (rGrading 1 1)
+
+    hex ( V3 V2 V5 V4 V15 V14 V17 V16) BC
+    (xNumberOfCells zNumberOfCells1 yNumberOfCells)
+    simpleGrading (rGrading 1 1)
+
+    hex ( V4 V5 V6 V7 V16 V17 V18 V19) CD
+    (xNumberOfCells zNumberOfCells2 yNumberOfCells)
+    simpleGrading (rGrading 1 1)
+
+    hex ( V7 V6 V9 V8 V19 V18 V21 V20 ) EF
+    (xNumberOfCells zNumberOfCells3 yNumberOfCells)
+    simpleGrading (rGrading 1 1)
+
+    hex ( V8 V9 V10 V11 V20 V21 V22 V23 ) GH
+    (xNumberOfCells zNumberOfCells4 yNumberOfCells)
+    simpleGrading (rGrading 1 1)
+);
+
+// Defining patches:
+boundary
+(
+    bottom
+    {
+        type patch;
+        faces
+        (
+         (V0 V1 V13 V12)
+        );
+    }
+    wall
+    {
+        type wall;
+        faces
+        (
+         (V0 V3 V15 V12)
+         (V3 V4 V16 V15)
+         (V4 V7 V19 V16)
+         (V7 V8 V20 V19)
+         (V8 V11 V23 V20)
+        );
+    }
+    top
+    {
+        type patch;
+        faces
+        (
+         (V11 V10 V22 V23)
+        );
+    }
+    inner
+    { 
+        type patch;
+        faces
+        (
+         (V0 V1 V2 V3)
+         (V3 V2 V5 V4)
+         (V4 V5 V6 V7)
+         (V7 V6 V9 V8)
+         (V8 V9 V10 V11)
+
+         (V12 V13 V14 V15)
+         (V15 V14 V17 V16)
+         (V16 V17 V18 V19)
+         (V19 V18 V21 V20)
+         (V20 V21 V22 V23)
+        );
+    }
+);
+
+mergePatchPairs 
+(
+);
+
+// ************************************************************************* //

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/CFD/system/controlDict

@@ -0,0 +1,132 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  2.3.0                                 |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       dictionary;
+    location    "system";
+    object      controlDict;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+application     potentialFoam;
+
+startFrom       startTime;
+
+startTime       0;
+
+stopAt          endTime;
+
+endTime         50.0;
+
+deltaT          0.00025;
+
+writeControl    timeStep;
+
+writeInterval   10000;
+
+purgeWrite      0;
+
+writeFormat     ascii;
+
+writePrecision  6;
+
+writeCompression off;
+
+timeFormat      general;
+
+timePrecision   6;
+
+runTimeModifiable true;
+
+adjustTimeStep no;
+
+functions
+{
+    fieldAverage1
+    {
+        type            fieldAverage;
+        libs            ("libfieldFunctionObjects.so");
+        writeControl    writeTime;
+
+        fields
+        (
+            U
+            {
+                mean        on;
+                prime2Mean  off;
+                base        time;
+                window	    10.0001;
+            }
+	    Us
+            {
+                mean        on;
+                prime2Mean  off;
+                base        time;
+                window	    10.0001;
+            }
+	    voidfraction
+            {
+                mean        on;
+                prime2Mean  off;
+                base        time;
+                window	    10.0001;
+            }
+        );
+    }
+
+    fieldOutput
+    {
+        type writeObjects;
+        functionObjectLibs ( "libutilityFunctionObjects.so" );
+        exclusiveWriting     true;
+        objects ("Us" "UsMean" "voidfraction" "voidfractionMean" "partType" "dSauter" "dSauterMean");
+        writeControl timeStep;
+        writeInterval       10000;
+    }
+
+    probes1
+    {
+        type probes;
+
+	functionObjectLibs ("libsampling.so");
+
+	#include "probesDict";
+    }
+
+    cellObj1
+    {
+        type            volFieldValue;
+        libs ("libfieldFunctionObjects.so");
+        enabled         true;
+        writeControl    timeStep;
+        writeInterval   400;
+        log             true;
+        writeFields     false;
+        regionType      cellZone;
+        name            topReg;
+        operation       volAverage;
+
+        fields
+        (
+            UsMean
+        );
+    }
+
+    conditionalStop
+    {
+        type abort;
+        libs    ("libutilityFunctionObjects.so");
+        action writeNow;
+        file "ABORT";
+    }
+
+}
+
+// ************************************************************************* //

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/CFD/system/createPatchDict

@@ -0,0 +1,72 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  2.3.0                                 |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       dictionary;
+    object      createPatchDict;
+}
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+// This application/dictionary controls:
+// - optional: create new patches from boundary faces (either given as
+//   a set of patches or as a faceSet)
+// - always: order faces on coupled patches such that they are opposite. This
+//   is done for all coupled faces, not just for any patches created.
+// - optional: synchronise points on coupled patches.
+
+// 1. Create cyclic:
+// - specify where the faces should come from
+// - specify the type of cyclic. If a rotational specify the rotationAxis
+//   and centre to make matching easier
+// - always create both halves in one invocation with correct 'neighbourPatch'
+//   setting.
+// - optionally pointSync true to guarantee points to line up.
+
+// 2. Correct incorrect cyclic:
+// This will usually fail upon loading:
+//  "face 0 area does not match neighbour 2 by 0.0100005%"
+//  " -- possible face ordering problem."
+// - in polyMesh/boundary file:
+//      - loosen matchTolerance of all cyclics to get case to load
+//      - or change patch type from 'cyclic' to 'patch'
+//        and regenerate cyclic as above
+
+
+// Do a synchronisation of coupled points after creation of any patches.
+// Note: this does not work with points that are on multiple coupled patches
+//       with transformations (i.e. cyclics).
+pointSync false;
+
+// Patches to create.
+patches
+(
+   {
+        name inletLeft;
+        patchInfo
+        {
+            type patch;
+        }
+        constructFrom set;
+        set inletLeft;
+    }
+
+{
+        name inletRight;
+        patchInfo
+        {
+            type patch;
+        }
+        constructFrom set;
+        set inletRight;
+    }
+);
+
+// ************************************************************************* //

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/CFD/system/decomposeParDict

@@ -0,0 +1,32 @@
+/*-------------------------------*- C++ -*---------------------------------*\
+|    =========                                                              |
+|    \\      /     OpenFOAM                                                 |
+|     \\    /                                                               |
+|      \\  /       The Open Source CFD Toolbox                              |
+|       \\/                                        http://www.OpenFOAM.org  |
+\*-------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       dictionary;
+    note        "mesh decomposition control dictionary";
+    location    "system";
+    object      decomposeParDict;
+}
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+numberOfSubdomains  8;
+
+method          simple;
+
+simpleCoeffs
+{
+    n		(2 1 4); 
+  // n           (2 2 8);
+    delta       0.001;
+}
+
+// ************************************************************************* //
+

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/CFD/system/fvOptions

@@ -0,0 +1,30 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  2.3.0                                 |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       dictionary;
+    location    "system";
+    object      fvOptions;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+source1
+{
+    type        limitTemperature;
+    active          yes;
+    limitTemperatureCoeffs
+    {
+        active          yes;
+        selectionMode   all;
+        min     273;
+        max     313;
+    }
+}
+

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/CFD/system/fvSchemes

@@ -0,0 +1,82 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  1.6                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       dictionary;
+    location    "system";
+    object      fvSchemes;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+ddtSchemes
+{
+    default         Euler;
+}
+
+gradSchemes
+{
+    default         Gauss linear;
+    grad(p)         cellLimited Gauss linear 1;
+    // grad(U)         Gauss linear;
+    grad(U)         cellLimited Gauss linear 1;
+}
+
+divSchemes
+{
+    default         Gauss linear;
+     div(phi,U)      Gauss upwind;
+   // div(phi,U)      bounded Gauss linearUpwindV grad(U);
+    div(phi,k)      Gauss limitedLinear 1;
+    div(phi,epsilon) Gauss limitedLinear 1;
+    div(phi,R)      Gauss limitedLinear 1;
+    div(R)          Gauss linear;
+    div(phi,nuTilda) Gauss limitedLinear 1;
+    div((viscousTerm*dev(grad(U).T()))) Gauss linear;
+    div((nu*dev(grad(U).T()))) Gauss linear;
+    div((nuEff*dev(grad(U).T()))) Gauss linear;
+}
+
+laplacianSchemes
+{
+    /*
+    default         Gauss linear corrected;
+    laplacian(viscousTerm,U) Gauss linear corrected;
+    laplacian(nu,U) Gauss linear corrected;
+    laplacian(nuEff,U) Gauss linear corrected;
+    laplacian((1|A(U)),p) Gauss linear corrected;
+    laplacian((voidfraction2|A(U)),p) Gauss linear corrected;
+    laplacian(DkEff,k) Gauss linear corrected;
+    laplacian(DepsilonEff,epsilon) Gauss linear corrected;
+    laplacian(DREff,R) Gauss linear corrected;
+    laplacian(DnuTildaEff,nuTilda) Gauss linear corrected;
+    */
+    default         Gauss linear limited 0.333;
+}
+
+interpolationSchemes
+{
+    default         linear;
+    interpolate(U)  linear;
+}
+
+snGradSchemes
+{
+    // default         corrected;
+    default         limited 0.333;
+}
+
+fluxRequired
+{
+    default         no;
+    p               ;
+}
+
+
+// ************************************************************************* //

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/CFD/system/fvSolution

@@ -0,0 +1,115 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  2.3.0                                 |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       dictionary;
+    location    "system";
+    object      fvSolution;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+solvers
+{
+    "(p|rho)"
+    {
+        solver          PCG;
+        preconditioner  DIC;
+        tolerance       1e-8;//1e-6;
+        relTol          0.01;
+    }
+
+    "(p|rho)Final"
+    {
+        $p;
+        relTol          0;
+    }
+
+    "(U|k|e|epsilon|R|nuTilda|c)"
+    {
+        solver          smoothSolver;
+        smoother        GaussSeidel;
+        tolerance       1e-7;//1e-5;
+        relTol          0;
+    }
+
+    "(U|e|k|nuTilda)Final"
+    {
+        $U;
+        relTol          0;
+    }
+
+    T
+    {
+        solver          PBiCG;
+        preconditioner  DILU;
+        tolerance       1e-8;
+        relTol          0;
+    }
+
+    alphaSt
+    {
+        //solver          smoothSolver;
+        //smoother        GaussSeidel;
+	solver          PBiCG;
+        preconditioner  DILU;
+        tolerance       1e-05;
+        relTol          0;
+    }
+
+    alphaDyn
+    {
+        //solver          smoothSolver;
+        //smoother        GaussSeidel;
+	solver          PBiCG;
+        preconditioner  DILU;
+        tolerance       1e-05;
+        relTol          0;
+    }
+
+    "(voidfraction|Ksl|UsNext|voidfractionNext|sSmoothField|vSmoothField|fSmooth)"
+    {
+        solver          PCG;
+        preconditioner  DIC;
+        tolerance       1e-09;
+        relTol          1e-06;
+    }
+}
+
+PIMPLE
+{
+    momentumPredictor yes;
+    nOuterCorrectors 2;
+    nCorrectors     1;
+    nNonOrthogonalCorrectors 0;
+    rhoMin          rhoMin [ 1 -3 0 0 0 ] 0.5;
+    rhoMax          rhoMax [ 1 -3 0 0 0 ] 2.0;
+}
+
+PISO
+{
+    nCorrectors     4;
+    nNonOrthogonalCorrectors 1;
+    pRefCell        0;
+    pRefValue       0;
+}
+
+relaxationFactors
+{
+    fields
+    {
+    }
+    equations
+    {
+        ".*"  1;
+    }
+}
+
+
+// ************************************************************************* //

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/CFD/system/mirrorMeshDict

@@ -0,0 +1,27 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  2.3.0                                 |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       dictionary;
+    object      mirrorMeshDict;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+planeType           pointAndNormal;
+
+pointAndNormalDict
+{
+    basePoint       (0 0 0);
+    normalVector    (-1 0 0);
+}
+
+planeTolerance      1e-3;
+
+// ************************************************************************* //

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/CFD/system/probesDict

@@ -0,0 +1,63 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \    /   O peration     | Version:  2.1.x                                 |
+|   \  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+	version		2.0;
+	format		ascii;
+	class		dictionary;
+	object		probesDict;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+
+
+fields
+(
+    voidfraction
+    voidfractionMean
+    rho
+    p
+    Us
+    UsMean
+    U
+    T
+    massFluxField
+    massFluxFieldMean
+);
+
+writeControl   timeStep;
+writeInterval  400;
+
+
+// Locations to be probed.
+probeLocations
+(
+	(-5.9 0.0 0.0)
+	(-5.5 0.0 0.0)
+	(-5.0 0.0 0.0)
+	(-4.5 0.0 0.0)
+	(-4.0 0.0 0.0)
+	(-3.5 0.0 0.0)
+	(5.9 0.0 0.0)
+	(5.5 0.0 0.0)
+	(5.0 0.0 0.0)
+	(4.5 0.0 0.0)
+	(4.0 0.0 0.0)
+	(3.5 0.0 0.0)
+	(0 0.0 -1.0)
+	(0 0.0 0.0)
+        (0 0.0 1.0)
+        (0 0.0 2.0)
+	(0 0.0 5.0)
+        (0 0.0 10.0)
+        (0 0.0 15.0)
+        (0 0.0 20.0)
+
+);
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/CFD/system/setFieldsDict

@@ -0,0 +1,41 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     | Website:  https://openfoam.org
+    \\  /    A nd           | Version:  6
+     \\/     M anipulation  |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       dictionary;
+    location    "system";
+    object      setFieldsDict;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+defaultFieldValues
+(
+    volScalarFieldValue partTempRef 300
+    volScalarFieldValue partTemp 300
+);
+
+regions
+(
+    rotatedBoxToCell
+    {
+       origin   (0.0 -0.6 8);
+       i        (20 0 -10);
+       j        (-20 0 -10);
+       k        (0 1.2 0);
+        fieldValues
+        (
+            volScalarFieldValue partTempRef 2700
+            volScalarFieldValue partTemp 2700
+        );
+    }
+);
+
+
+// ************************************************************************* //

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/CFD/system/topoSetDict

@@ -0,0 +1,134 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  2.3.0                                 |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       dictionary;
+    object      topoSetDict;
+}
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+tuydia 0.25;
+actions
+(
+  {
+    name    c_left;
+    type    cellSet;
+    action  new;
+    source  boxToCell;
+    sourceInfo
+    {
+       box (-6.5 #calc "-0.5*$tuydia" #calc "-0.25*$tuydia")(-5.5 #calc "0.5*$tuydia" #calc "0.75*$tuydia");
+    }
+  }
+
+ {
+    name    inletLeft;
+    type    faceSet;
+    action  new;
+    source  patchToFace;
+    sourceInfo
+    {
+       name "wall";
+    }
+  }
+
+ {
+    name    inletLeft;
+    type    faceSet;
+    action  subset;
+    source  cellToFace;
+    sourceInfo
+    {
+       set c_left;
+       option all;
+    }
+  }
+
+  {
+    name    c_right;
+    type    cellSet;
+    action  new;
+    source  boxToCell;
+    sourceInfo
+    {
+       box (5.5 #calc "-0.5*$tuydia" #calc "-0.25*$tuydia")(6.5 #calc "0.5*$tuydia" #calc "0.75*$tuydia");
+    }
+  }
+
+ {
+    name    inletRight;
+    type    faceSet;
+    action  new;
+    source  patchToFace;
+    sourceInfo
+    {
+       name "wall";
+    }
+  }
+
+ {
+    name    inletRight;
+    type    faceSet;
+    action  subset;
+    source  cellToFace;
+    sourceInfo
+    {
+       set c_right;
+       option all;
+    }
+  }
+
+{
+    name    rw_left;
+    type    cellSet;
+    action  new;
+    source  boxToCell;
+    sourceInfo
+    {
+       box (-6.0 -0.5 -0.5)(-4.5 0.5 0.5);
+    }
+  }
+
+{
+    name    rw_right;
+    type    cellSet;
+    action  new;
+    source  boxToCell;
+    sourceInfo
+    {
+       box (4.5 -0.5 -0.5)(6.0 0.5 0.5);
+    }
+  }
+
+{
+    name    topReg;
+    type    cellSet;
+    action  new;
+    source  boxToCell;
+    sourceInfo
+    {
+       box (-5.0 -1.0 25)(5.0 1.0 26);
+    }
+}
+
+{
+    name    noOreReg;
+    type    cellSet;
+    action  new;
+    source  boxToCell;
+    sourceInfo
+    {
+       box (-8.0 -1.0 -2.0)(8.0 1.0 3.0);
+    }
+}
+
+);
+
+// ************************************************************************* //

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/DEM/DEMparameters

@@ -0,0 +1,94 @@
+variable dt equal 0.000025
+variable skin equal 0.001
+variable Ndump equal 100000
+variable NsetDM equal 1000
+
+# DEM-parameters COKE
+variable yM_C equal 5.e7          # Young's modulus
+variable pR_C equal 0.25          # Poisson ratio
+variable rest_C equal 0.5         # coefficient of restitution
+variable fric_C equal 0.9         # coefficient of friction
+variable roll_C equal 0.5         # coefficient of rolling friction
+
+# values according to information from VAS
+# mean diameter coke = 0.044 m
+# mean density coke = 950
+# bulk density = 500
+
+variable dens_C equal 950	  # density
+variable rad_C equal 0.088	  # c.g. radius (cg = 4)
+variable poly_C equal 0.81	  # polydisp factor
+
+
+# DEM-parameters FINE COKE
+variable yM_fC equal 5.e7          # Young's modulus
+variable pR_fC equal 0.25          # Poisson ratio
+variable rest_fC equal 0.5         # coefficient of restitution
+variable fric_fC equal 0.9         # coefficient of friction
+variable roll_fC equal 0.5         # coefficient of rolling friction
+
+# mean diameter fine coke = 0.02 m
+# mean density fine coke = 950
+# bulk density fine coke = 570
+
+variable dens_fC equal 950	  # density
+variable rad_fC equal 0.06	  # c.g. radius (cg = 6)
+variable poly_fC equal 0.92	  # polydisp factor
+
+
+# DEM-parameters ORE (represents the mixture of sinter, pellets and actual ore)
+variable yM_O equal 5.e7          # Young's modulus
+variable pR_O equal 0.25          # Poisson ratio
+variable rest_O equal 0.5         # coefficient of restitution
+variable fric_O equal 0.8         # coefficient of friction
+variable roll_O equal 0.4        # coefficient of rolling friction
+
+# values according to Modern Blast Furnace Ironmaking
+# mean diameter ore mixture = 0.015 m
+# mean density ore mixture = 3300
+# bulk density = 2200
+
+variable dens_O equal 3300	  # density
+variable rad_O equal 0.075	  # c.g. radius (cg = 10)
+variable poly_O equal 1.03	  # polydisp factor
+
+
+# DEM-parameters FINE ORE (represents the mixture of sinter, pellets and actual ore)
+variable yM_fO equal 5.e7          # Young's modulus
+variable pR_fO equal 0.25          # Poisson ratio
+variable rest_fO equal 0.5         # coefficient of restitution
+variable fric_fO equal 0.8         # coefficient of friction
+variable roll_fO equal 0.4        # coefficient of rolling friction
+
+# values according to G. Holzinger's literature study
+# mean diameter fine ore mixture = 0.006 m
+# mean density fine ore mixture = 3300
+# bulk density = 2200
+
+variable dens_fO equal 3300	  # density
+variable rad_fO equal 0.06	  # c.g. radius (cg = 20)
+variable poly_fO equal 1.03	  # polydisp factor
+
+
+
+
+# DEM-parameters wall
+variable yM_W equal 5.e6          # Young's modulus
+variable pR_W equal 0.25          # Poisson ratio
+variable rest_W equal 0.5         # coefficient of restitution
+variable fric_W equal 0.85         # coefficient of friction
+variable roll_W equal 0.45        # coefficient of rolling friction
+
+# DEM-parameters mixed
+variable rest_CO equal 0.2         # coefficient of restitution
+variable rest_WO equal 0.2         # coefficient of restitution
+variable rest_CW equal 0.5         # coefficient of restitution
+variable fric_CO equal 0.8         # coefficient of friction
+variable fric_WO equal 0.2         # coefficient of friction
+variable fric_CW equal 0.5         # coefficient of friction
+variable roll_CO equal 0.02         # for particles and walls with friction
+variable roll_WO equal 0.02         # for particles and walls with friction
+variable roll_CW equal 0.02         # for particles and walls with friction
+
+
+variable cV equal 0.1           # characteristic velocity

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/DEM/in.DM

@@ -0,0 +1,13 @@
+# adapt particle parameters in deadman region
+
+group DMcoke region DM
+set group DMcoke property/atom/polydispparcel effvolfactor 1.35
+
+group hearthcoke region hearth
+
+group nonDMcoke subtract hearthcoke DMcoke
+set group nonDMcoke property/atom/polydispparcel effvolfactor ${poly_C}
+
+group DMcoke clear
+group hearthcoke clear
+group nonDMcoke clear

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/DEM/in.RW

@@ -0,0 +1,8 @@
+#########
+# removal
+#########
+
+# raceways
+variable rCMass equal 200
+fix         delRW1 all remove nevery 1000 massrate ${rCMass} style delete seed 10007143 region RW1 verbose no restart_read no
+fix         delRW2 all remove nevery 1000 massrate ${rCMass} style delete seed 10007134 region RW2 verbose no restart_read no

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/DEM/in.abort

@@ -0,0 +1 @@
+print " " file ../CFD/ABORT

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/DEM/in.coupling_integration

@@ -0,0 +1,22 @@
+##############
+# cfd coupling
+##############
+
+variable Ncouple equal 10
+variable Tin equal 293
+
+fix		cfd1 all couple/cfd couple_every ${Ncouple} mpi
+fix		cfd2 all couple/cfd/force/implicit transfer_density no transfer_type yes transfer_property name effvolfactor type scalar-atom
+fix		cfd3 all couple/cfd/convection T0 293
+fix		cfd4 all couple/cfd/parttempfield T0 293
+fix		cfd5 all couple/cfd/deform rmin 0.01
+
+# make sure particle outside CFD domain are initialized with reasonable temperature
+# this is executed at the beginning of each run 
+region outsideCFD block INF INF INF INF 26.0 INF units box
+variable NFixTemp equal ${Ncouple}+2
+fix eeTemp all execute ${NFixTemp} " set region outsideCFD property/atom Temp ${Tin} " once
+
+fix		integr all nve/sphere
+fix		limV all limit/vel 10.0
+fix		limV RW limit/vel 0.01

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/DEM/in.final

@@ -0,0 +1,27 @@
+#remove out-of-CFD domain particles and uppermost, partial layer
+#write out uppermost complete layers
+
+delete_atoms group highestCoke
+delete_atoms group highestOre
+delete_atoms group highestFineMix
+run 0
+
+write_dump secHighestOre custom ../DEM/initOre id type radius mass density x y z
+write_dump secHighestCoke custom ../DEM/initCoke id type radius mass density x y z
+write_dump secHighestFineMix custom ../DEM/initFineMix id type radius mass density x y z
+
+write_dump secHighestOre custom/vtk ../DEM/initOre.vtk id type radius mass density x y z
+write_dump secHighestCoke custom/vtk ../DEM/initCoke.vtk id type radius mass density x y z
+write_dump secHighestFineMix custom/vtk ../DEM/initFineMix.vtk id type radius mass density x y z
+
+write_restart  	../DEM/liggghts.restartCFDEM
+
+# exclude near-wall artifacts
+region centerReg block -4.0 4.0 -0.59 0.59 INF INF units box
+group center region centerReg
+group remainder subtract center secHighestOre secHighestCoke secHighestFineMix
+write_dump remainder custom/vtk ../DEM/remainder.vtk id type radius mass density x y z
+variable zInsLower equal bound(remainder,zmax)
+
+print 'variable zInsLower equal ${zInsLower}' screen yes
+print "variable zInsLower equal ${zInsLower}" file ../DEM/zInsLower

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/DEM/in.insert

@@ -0,0 +1,26 @@
+print 'Performing insertion number ${counter}.' screen yes
+
+group secHighestCoke clear
+group secHighestOre clear
+group secHighestFineMix clear
+group secHighestCoke union secHighestCoke highestCoke
+group secHighestOre union secHighestOre highestOre
+group secHighestFineMix union secHighestFineMix highestFineMix
+
+group highestCoke clear
+group highestOre clear
+group highestFineMix clear
+
+if "${counter} > 1" then "unfix ins${counter}C"
+if "${counter} > 1" then "unfix ins${counter}O"
+if "${counter} > 1" then "unfix ins${counter}O1"
+if "${counter} > 1" then "unfix ins${counter}O2"
+
+next counter
+
+fix ins${counter}C highestCoke insert/pack/dense seed 125 region insRegTop2 distributiontemplate pdd1 volumefraction_region 0.56
+fix ins${counter}O highestOre insert/pack/dense seed 124 region insRegTop1b distributiontemplate pdd2 volumefraction_region 0.56
+fix ins${counter}O1 highestFineMix insert/pack/dense seed 123 region insRegTop1a distributiontemplate pdd3 volumefraction_region 0.56
+fix ins${counter}O2 highestFineMix insert/pack/dense seed 122 region insRegTop1c distributiontemplate pdd3 volumefraction_region 0.56
+
+run 0

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/DEM/in.insertion_setup

@@ -0,0 +1,38 @@
+#################
+# insertion setup
+#################
+
+variable outerRing equal 4.61
+variable innerRing equal 4.1
+region insRegTop1a block -${outerRing} -${innerRing} -0.59 0.59 26.125 26.95 units box
+region insRegTop1b block -${innerRing} ${innerRing} -0.59 0.59 26.125 26.95 units box
+region insRegTop1c block ${innerRing} ${outerRing} -0.59 0.59 26.125 26.95 units box
+region insRegTop2 block -${outerRing} ${outerRing} -0.59 0.59 26.95 27.57 units box
+
+region regTop union 4 insRegTop1a insRegTop1b insRegTop1c insRegTop2
+variable nTop equal count(all,regTop)
+
+group secHighestCoke empty
+group secHighestOre empty
+group secHighestFineMix empty
+group highestFineMix empty
+group highestCoke empty
+group highestOre empty
+
+# ensure reasonable initialization
+group secHighestCoke union all secHighestCoke
+group highestCoke union all highestCoke
+group highestFineMix union all highestFineMix
+variable zMaxSecHighCoke equal bound(secHighestCoke,zmax)
+
+variable Nins equal 1000
+variable counter loop 1000
+variable insertedCoke equal mass(highestCoke)
+variable insertedOre equal mass(highestOre)
+variable insertedFineMix equal mass(highestFineMix)
+fix eeIns all execute ${Nins} "../../CFDDEM/DEM/in.insert" file conditional nTop 0 0.1
+
+variable time equal step*dt
+variable np equal count(all)
+fix printN all print 4000 "${time} ${np}" file ../DEM/monitor/Npart_CFDDEMflow.txt title "#time   n_p"
+fix printM all print 4000 "${time} ${insertedCoke} ${insertedOre} ${insertedFineMix} " file ../DEM/monitor/Mpart_CFDDEMflow.txt title "#time   m_C m_O m_FM"

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/DEM/in.liggghts_CFDDEM

@@ -0,0 +1,42 @@
+# material properties
+include ../DEM/DEMparameters
+
+
+###############################################
+atom_style		granular
+atom_modify		map array
+communicate	single vel yes
+
+boundary		f p f
+newton			off
+units			si
+processors      1 1 8
+
+read_restart 		../DEM/liggghts.restart
+reset_timestep 0
+
+neighbor	${skin} bin
+neigh_modify	delay 0
+
+timestep    ${dt}
+
+include ../DEM/in.material_settings
+include ../DEM/in.regions
+include ../DEM/in.RW
+include ../DEM/in.coupling_integration
+include ../DEM/in.output
+
+run	1
+
+include ../DEM/in.set_polydisp
+include ../DEM/in.insertion_setup
+
+# conditionally stop simulation
+variable tmin equal 20.0
+variable zmaxCFD equal 26.0
+variable zminCFD equal 25.8
+variable abort equal (v_time>=${tmin})&&(v_zMaxSecHighCoke<=${zmaxCFD})&&(v_zMaxSecHighCoke>=${zminCFD})
+fix printAbort all print 4000 "${time} ${zMaxSecHighCoke} ${abort}" file ../DEM/monitor/Nabort.txt title "#time   z_max_2 abort"
+fix eeAbort all execute 1000 "../DEM/in.abort" file conditional abort 1 0.1
+
+restart 200000 ../DEM/restart.A ../DEM/restart.B

tutorials/rcfdemSolverRhoSteadyPimpleChem/oreReducingBlastFurnaceSlot/CFDDEM/DEM/in.material_settings

@@ -0,0 +1,48 @@
+##################################################
+# material properties required for new pair styles
+##################################################
+
+fix 		m1 all property/global youngsModulus peratomtype ${yM_C} ${yM_C} ${yM_O} ${yM_O} ${yM_W}
+fix 		m2 all property/global poissonsRatio peratomtype ${pR_C} ${pR_C} ${pR_O} ${pR_O} ${pR_W}
+fix 		m3 all property/global coefficientRestitution peratomtypepair 5 ${rest_C} ${rest_C} ${rest_CO} ${rest_CO} ${rest_CW} ${rest_C} ${rest_C} ${rest_CO} ${rest_CO} ${rest_CW} ${rest_CO} ${rest_CO} ${rest_O} ${rest_O} ${rest_WO}  ${rest_CO} ${rest_CO} ${rest_O} ${rest_O} ${rest_WO} ${rest_CW} ${rest_CW} ${rest_WO} ${rest_WO} ${rest_W}
+fix 		m4 all property/global coefficientFriction peratomtypepair 5 ${fric_C} ${fric_C} ${fric_CO} ${fric_CO} ${fric_CW} ${fric_C} ${fric_C} ${fric_CO} ${fric_CO} ${fric_CW} ${fric_CO} ${fric_CO} ${fric_O} ${fric_O} ${fric_WO} ${fric_CO} ${fric_CO} ${fric_O} ${fric_O} ${fric_WO} ${fric_CW} ${fric_CW}  ${fric_WO} ${fric_WO} ${fric_W}
+fix         	m5 all property/global coefficientRollingFriction peratomtypepair 5 ${roll_C} ${roll_C} ${roll_CO} ${roll_CO} ${roll_CW} ${roll_C} ${roll_C} ${roll_CO} ${roll_CO} ${roll_CW} ${roll_CO} ${roll_CO} ${roll_O} ${roll_O} ${roll_WO} ${roll_CO} ${roll_CO} ${roll_O} ${roll_O} ${roll_WO} ${roll_CW} ${roll_CW} ${roll_WO} ${roll_WO} ${roll_W}
+fix         	m6 all property/global characteristicVelocity scalar ${cV}
+
+fix 		ftco all property/global thermalConductivity peratomtype 1.7 1.7 80 80 1.0
+fix 		ftca all property/global thermalCapacity peratomtype 850 850 600 600 1.0 
+fix             heattransfer all heat/gran initial_temperature 293
+
+# pair style
+pair_style gran model hertz tangential history rolling_friction cdt
+pair_coeff  * *
+
+# gravity
+fix         gravi all gravity 9.81 vector 0.0 0.0 -1.0
+
+fix         bottom all mesh/surface file ../../CFDDEM/DEM/stl_files/surface_bottom.stl type 5
+fix         wall all mesh/surface file ../../CFDDEM/DEM/stl_files/surface_wall.stl type 5
+
+fix granwalls all wall/gran model hertz tangential history mesh n_meshes 2 meshes bottom wall
+
+
+#################################
+# particle types and distribution
+#################################
+
+group coke type 1
+group finecoke type 2
+group ore type 3
+group fineore type 4
+
+fix         pts1 all particletemplate/sphere 2 atom_type 1 density constant ${dens_C} radius constant ${rad_C}
+fix	    pdd1 all particledistribution/discrete 3 1 pts1 1.0
+
+fix         pts2 all particletemplate/sphere 5 atom_type 3 density constant ${dens_O} radius constant ${rad_O}
+fix	    pdd2 all particledistribution/discrete 7 1 pts2 1.0
+
+fix         pts3a all particletemplate/sphere 11 atom_type 2 density constant ${dens_fC} radius constant ${rad_fC}
+fix         pts3b all particletemplate/sphere 13 atom_type 4 density constant ${dens_fO} radius constant ${rad_fO}
+fix	    pdd3 all particledistribution/discrete 17 2 pts3a 0.2 pts3b 0.8 
+
+fix	    fp all property/atom/polydispparcel effvolfactor ${poly_C} ${poly_fC} ${poly_O} ${poly_fO} 1.0