release on 2013-03-26_15-27-17

README

@@ -1,81 +0,0 @@
-/*---------------------------------------------------------------------------*\
-    CFDEMcoupling - Open Source CFD-DEM coupling
-
-    CFDEMcoupling is part of the CFDEMproject
-    www.cfdem.com
-                                Christoph Goniva, christoph.goniva@cfdem.com
-                                Copyright 2009-2012 JKU Linz
-                                Copyright 2012-     DCS Computing GmbH, Linz
--------------------------------------------------------------------------------
-License
-    This file is part of CFDEMcoupling.
-
-    CFDEMcoupling is free software; you can redistribute it and/or modify it
-    under the terms of the GNU General Public License as published by the
-    Free Software Foundation; either version 3 of the License, or (at your
-    option) any later version.
-
-    CFDEMcoupling is distributed in the hope that it will be useful, but WITHOUT
-    ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
-    FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
-    for more details.
-
-    You should have received a copy of the GNU General Public License
-    along with CFDEMcoupling; if not, write to the Free Software Foundation,
-    Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
-
-Description
-    This code is designed to realize coupled CFD-DEM simulations using LIGGGHTS
-    and OpenFOAM. Note: this code is not part of OpenFOAM (see DISCLAIMER).
-\*---------------------------------------------------------------------------*/
-
-
-CFDEM coupling provides an open source parallel coupled CFD-DEM framework 
-combining the strengths of LIGGGHTS DEM code and the Open Source 
-CFD package OpenFOAM(R)(*). The CFDEMcoupling toolbox allows to expand 
-standard CFD solvers of OpenFOAM(R)(*) to include a coupling to the DEM 
-code LIGGGHTS. In this toolbox the particle representation within the 
-CFD solver is organized by "cloud" classes. Key functionalities are organised 
-in sub-models (e.g. force models, data exchange models, etc.) which can easily 
-be selected and combined by dictionary settings.
-
-The coupled solvers run fully parallel on distributed-memory clusters. 
-
-Features are:
-
-- its modular approach allows users to easily implement new models
-- its MPI parallelization enables to use it for large scale problems
-- the "forum"_lws on CFD-DEM gives the possibility to exchange with other 
-  users / developers
-- the use of GIT allows to easily update to the latest version
-- basic documentation is provided
-
-The file structure:
-
-- "src" directory including the source files of the coupling toolbox and models
-- "applications" directory including the solver files for coupled CFD-DEM simulations
-- "doc" directory including the documentation of CFDEMcoupling
-- "tutorials" directory including basic tutorial cases showing the functionality
-
-
-
-Details on installation are given on the "www.cfdem.com"
-
-The functionality of this CFD-DEM framwork is described via "tutorial cases" showing 
-how to use different solvers and models.
-
-CFDEMcoupling stands for Computational Fluid Dynamics (CFD) -
-Discrete Element Method (DEM) coupling.
-
-CFDEMcoupling is an open-source code, distributed freely under the terms of the 
-GNU Public License (GPL).
-
-Core development of CFDEMcoupling is done by 
-Christoph Goniva and Christoph Kloss, both at DCS Computing GmbH, 2012
-
-
-\*---------------------------------------------------------------------------*/
-(*) "OpenFOAM(R)"_of is a registered trade mark of the ESI Group.
-This offering is not affiliated, approved or endorsed by ESI Group, 
-the producer of the OpenFOAM® software and owner of the OpenFOAM® trade mark.
-\*---------------------------------------------------------------------------*/

applications/solvers/cfdemSolverPimple/Make/files

@@ -0,0 +1,3 @@
+cfdemSolverPimple.C
+
+EXE = $(CFDEM_APP_DIR)/cfdemSolverPimple

applications/solvers/cfdemSolverPimple/Make/options

@@ -0,0 +1,17 @@
+EXE_INC = \
+    -I$(LIB_SRC)/turbulenceModels/incompressible/turbulenceModel \
+    -I$(LIB_SRC)/transportModels \
+    -I$(LIB_SRC)/transportModels/incompressible/singlePhaseTransportModel \
+    -I$(LIB_SRC)/finiteVolume/lnInclude \
+    -I$(LIB_SRC)/meshTools/lnInclude\
+    -I$(CFDEM_SRC_DIR)/lnInclude \
+
+EXE_LIBS = \
+    -L$(FOAM_USER_LIBBIN)\
+    -lincompressibleTransportModels \
+    -lincompressibleTurbulenceModel \
+    -lincompressibleRASModels \
+    -lincompressibleLESModels \
+    -lfiniteVolume \
+    -lmeshTools \
+    -l$(CFDEM_LIB_NAME)

applications/solvers/cfdemSolverPimple/cfdemSolverPimple.C

@@ -0,0 +1,116 @@
+/*---------------------------------------------------------------------------*\
+    CFDEMcoupling - Open Source CFD-DEM coupling
+
+    CFDEMcoupling is part of the CFDEMproject
+    www.cfdem.com
+                                Christoph Goniva, christoph.goniva@cfdem.com
+                                Copyright (C) 2011      OpenFOAM Foundation
+                                Copyright (C) 2009-2012 JKU, Linz
+                                Copyright (C) 2012-     DCS Computing GmbH,Linz
+-------------------------------------------------------------------------------
+License
+    This file is part of CFDEMcoupling.
+
+    CFDEMcoupling is free software: you can redistribute it and/or modify it
+    under the terms of the GNU General Public License as published by
+    the Free Software Foundation, either version 3 of the License, or
+    (at your option) any later version.
+
+    CFDEMcoupling is distributed in the hope that it will be useful, but WITHOUT
+    ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
+    FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
+    for more details.
+
+    You should have received a copy of the GNU General Public License
+    along with CFDEMcoupling.  If not, see <http://www.gnu.org/licenses/>.
+
+Application
+    cfdemSolverPimple
+
+Description
+    Large time-step transient solver for incompressible, flow using the PIMPLE
+    (merged PISO-SIMPLE) algorithm. Turbulence modelling is generic,
+    i.e. laminar, RAS or LES may be selected.
+    The code is an evolution of the solver pimpleFoam in OpenFOAM(R) 2.1.x., 
+    where additional functionality for CFD-DEM coupling is added.
+\*---------------------------------------------------------------------------*/
+
+#include "fvCFD.H"
+#include "singlePhaseTransportModel.H"
+#include "turbulenceModel.H"
+#include "pimpleControl.H"
+#include "IObasicSourceList.H"
+
+#include "cfdemCloud.H"
+#include "implicitCouple.H"
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+int main(int argc, char *argv[])
+{
+    #include "setRootCase.H"
+    #include "createTime.H"
+    #include "createMesh.H"
+    #include "createFields.H"
+    #include "initContinuityErrs.H"
+
+    pimpleControl pimple(mesh);
+
+    // create cfdemCloud
+    #include "readGravitationalAcceleration.H"
+    cfdemCloud particleCloud(mesh);
+    #include "checkModelType.H"
+
+    // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+    Info<< "\nStarting time loop\n" << endl;
+
+    while (runTime.run())
+    {
+        #include "readTimeControls.H"
+        #include "CourantNo.H"
+        #include "setDeltaT.H"
+
+        runTime++;
+
+        Info<< "Time = " << runTime.timeName() << nl << endl;
+
+        // do particle stuff
+        Info << "- evolve()" << endl;
+        particleCloud.evolve(voidfraction,Us,U);
+
+        Info << "update Ksl.internalField()" << endl;
+        Ksl.internalField() = particleCloud.momCoupleM(0).impMomSource();
+        Ksl.correctBoundaryConditions();
+
+        // --- Pressure-velocity PIMPLE corrector loop
+        while (pimple.loop())
+        {
+            #include "UEqn.H"
+
+            // --- Pressure corrector loop
+            while (pimple.correct())
+            {
+                #include "pEqn.H"
+            }
+
+            if (pimple.turbCorr())
+            {
+                turbulence->correct();
+            }
+        }
+
+        runTime.write();
+
+        Info<< "ExecutionTime = " << runTime.elapsedCpuTime() << " s"
+            << "  ClockTime = " << runTime.elapsedClockTime() << " s"
+            << nl << endl;
+    }
+
+    Info<< "End\n" << endl;
+
+    return 0;
+}
+
+
+// ************************************************************************* //

src/lagrangian/cfdemParticle/Make/files

@@ -6,7 +6,6 @@ voidFractionModels = subModels/voidFractionModel
 locateModels = subModels/locateModel
 meshMotionModels = subModels/meshMotionModel
 momCoupleModels = subModels/momCoupleModel
-regionModels = subModels/regionModel
 dataExchangeModels = subModels/dataExchangeModel
 averagingModels = subModels/averagingModel
 clockModels = subModels/clockModel
@@ -21,18 +20,33 @@ $(forceModels)/forceModel/newForceModel.C
 $(forceModels)/noDrag/noDrag.C
 $(forceModels)/checkCouplingInterval/checkCouplingInterval.C
 $(forceModels)/DiFeliceDrag/DiFeliceDrag.C
+$(forceModels)/DiFeliceDragNLift/DiFeliceDragNLift.C
 $(forceModels)/GidaspowDrag/GidaspowDrag.C
 $(forceModels)/SchillerNaumannDrag/SchillerNaumannDrag.C
 $(forceModels)/Archimedes/Archimedes.C
 $(forceModels)/ArchimedesIB/ArchimedesIB.C
 $(forceModels)/interface/interface.C
 $(forceModels)/ShirgaonkarIB/ShirgaonkarIB.C
+$(forceModels)/fieldTimeAverage/fieldTimeAverage.C
+$(forceModels)/fieldBound/fieldBound.C
+$(forceModels)/volWeightedAverage/volWeightedAverage.C
+$(forceModels)/totalMomentumExchange/totalMomentumExchange.C
 $(forceModels)/KochHillDrag/KochHillDrag.C
+$(forceModels)/BeetstraDrag/multiphaseFlowBasic/multiphaseFlowBasic.C
+$(forceModels)/BeetstraDrag/BeetstraDrag.C
+$(forceModels)/LaEuScalarLiquid/LaEuScalarLiquid.C
 $(forceModels)/LaEuScalarTemp/LaEuScalarTemp.C
+$(forceModels)/LaEuScalarDust/LaEuScalarDust.C
 $(forceModels)/virtualMassForce/virtualMassForce.C
 $(forceModels)/gradPForce/gradPForce.C
+$(forceModels)/gradULiftForce/gradULiftForce.C
 $(forceModels)/viscForce/viscForce.C
 $(forceModels)/MeiLift/MeiLift.C
+$(forceModels)/KochHillDragNLift/KochHillDragNLift.C
+$(forceModels)/solidsPressureForce/solidsPressureForce.C
+$(forceModels)/periodicPressure/periodicPressure.C
+$(forceModels)/periodicPressureControl/periodicPressureControl.C
+$(forceModels)/averageSlipVel/averageSlipVel.C
 
 $(forceModelsMS)/forceModelMS/forceModelMS.C
 $(forceModelsMS)/forceModelMS/newForceModelMS.C
@@ -62,21 +76,18 @@ $(locateModels)/turboEngineSearch/turboEngineSearch.C
 $(locateModels)/turboEngineSearchM2M/turboEngineSearchM2M.C
 $(locateModels)/engineSearchIB/engineSearchIB.C
 
-
 $(meshMotionModels)/meshMotionModel/meshMotionModel.C
 $(meshMotionModels)/meshMotionModel/newMeshMotionModel.C
 $(meshMotionModels)/noMeshMotion/noMeshMotion.C
+$(meshMotionModels)/DEMdrivenMeshMotion/DEMdrivenMeshMotion.C
 
 $(momCoupleModels)/momCoupleModel/momCoupleModel.C
 $(momCoupleModels)/momCoupleModel/newMomCoupleModel.C
 $(momCoupleModels)/explicitCouple/explicitCouple.C
+$(momCoupleModels)/explicitCoupleSource/explicitCoupleSource.C
 $(momCoupleModels)/implicitCouple/implicitCouple.C
 $(momCoupleModels)/noCouple/noCouple.C
 
-$(regionModels)/regionModel/regionModel.C
-$(regionModels)/regionModel/newRegionModel.C
-$(regionModels)/allRegion/allRegion.C
-
 $(dataExchangeModels)/dataExchangeModel/dataExchangeModel.C
 $(dataExchangeModels)/dataExchangeModel/newDataExchangeModel.C
 $(dataExchangeModels)/oneWayVTK/oneWayVTK.C
@@ -102,4 +113,4 @@ $(liggghtsCommandModels)/runLiggghts/runLiggghts.C
 $(liggghtsCommandModels)/writeLiggghts/writeLiggghts.C
 $(liggghtsCommandModels)/readLiggghtsData/readLiggghtsData.C
 
-LIB = $(FOAM_USER_LIBBIN)/lib$(CFDEM_LIB_NAME)
+LIB = $(CFDEM_LIB_DIR)/lib$(CFDEM_LIB_NAME)

tutorials/cfdemSolverPimple/ErgunTestMPI/Allrun.sh

@@ -0,0 +1,65 @@
+#!/bin/bash
+
+#===================================================================#
+# allrun script for testcase as part of test routine 
+# run settlingTest
+# Christoph Goniva - March 2012
+#===================================================================#
+
+#- define variables
+casePath="$(dirname "$(readlink -f ${BASH_SOURCE[0]})")"
+pizzaPath="$CFDEM_PIZZA_DIR"
+postproc="false"
+nrPostProcProcessors=2
+
+# check if mesh was built
+if [ -d "$casePath/CFD/constant/polyMesh/boundary" ]; then
+    echo "mesh was built before - using old mesh"
+else
+    echo "mesh needs to be built"
+    cd $casePath/CFD
+    blockMesh
+fi
+
+
+
+#- run parallel CFD-DEM in new terminal
+gnome-terminal --title='cfdemSolverPimple ErgunTestMPI CFD'  -e "bash $casePath/parCFDDEMrun.sh" 
+
+if [ $postproc == "true" ]
+  then
+
+    #- keep terminal open (if started in new terminal)
+    echo "simulation finisehd? ...press enter to proceed"
+    read
+
+    #- get VTK data from liggghts dump file
+    cd $casePath/DEM
+    python $pizzaPath/pizza.py -f pizzaScriptRestart
+
+    #- get VTK data from CFD sim
+    cd $casePath/CFD
+    #foamToVTK                                                   #- serial run of foamToVTK
+    source $CFDEM_SRC_DIR/etc/functions.sh                       #- include functions
+    pseudoParallelRun "foamToVTK" $nrPostProcProcessors          #- pseudo parallel run of foamToVTK
+
+    #- start paraview
+    paraview
+
+    #- keep terminal open (if started in new terminal)
+    echo "...press enter to clean up case"
+    echo "press Ctr+C to keep data"
+    read
+
+    #- clean up case
+    echo "deleting data at: $casePath"
+    rm -r $casePath/CFD/0.*
+    rm -r $casePath/CFD/VTK
+    rm -r $casePath/CFD/couplingFiles/*
+    rm -r $casePath/CFD/particles/*
+    rm -r $casePath/DEM/post/*
+    rm -r $casePath/DEM/log.*
+    rm -r $casePath/log_*
+    echo "done"
+
+fi

tutorials/cfdemSolverPimple/ErgunTestMPI/CFD/0/Ksl

@@ -0,0 +1,41 @@
+/*--------------------------------*- 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
+{
+    wall
+    {
+        type            zeroGradient;
+        //type            fixedValue;
+        //value           uniform 0;
+    }
+    inlet
+    {
+        type            zeroGradient;
+    }
+    outlet
+    {
+        type            zeroGradient;
+    }
+}
+
+
+// ************************************************************************* //

tutorials/cfdemSolverPimple/ErgunTestMPI/CFD/0/U

@@ -0,0 +1,81 @@
+/*--------------------------------*- 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      U;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 1 -1 0 0 0 0];
+
+internalField   uniform (0 0 0);
+
+boundaryField
+{
+
+    wall
+    {
+        //type            fixedValue;
+        //value           uniform (0 0 0);
+        type            slip;
+    }
+
+    inlet
+    {
+        /*type            flowRateInletVelocity;
+        flowRate        0.001;
+        value           uniform (0 0 0);*/
+
+        /*type            fixedValue;
+        value           uniform (0 0 0.0001);*/
+
+        //type            zeroGradient;
+
+/*        type            groovyBC;
+        variables       "Uend=vector(0,0,0.02);tEnd=0.1;";
+        valueExpression "((time() < tEnd) ? Uend/tEnd*time():Uend)";
+        value           uniform (0 0 0);*/
+
+        type            uniformFixedValue;
+        uniformValue     table
+        (
+            (0.000 (0 0 0.002))
+            (0.010 (0 0 0.002))
+            (0.011 (0 0 0.004))
+            (0.020 (0 0 0.004))
+            (0.021 (0 0 0.006))
+            (0.030 (0 0 0.006))
+            (0.031 (0 0 0.008))
+            (0.040 (0 0 0.008))
+            (0.041 (0 0 0.010))
+            (0.050 (0 0 0.010))
+            (0.051 (0 0 0.012))
+            (0.060 (0 0 0.012))
+            (0.061 (0 0 0.014))
+            (0.070 (0 0 0.014))
+            (0.071 (0 0 0.016))
+            (0.080 (0 0 0.016))
+            (0.081 (0 0 0.018))
+            (0.090 (0 0 0.018))
+            (0.091 (0 0 0.020))
+            (0.100 (0 0 0.020))
+        ); 
+    }
+    outlet
+    {
+        /*type            fluxCorrectedVelocity; //inletOutlet;
+        value           uniform (0 0 0);
+        inletValue      uniform (0 0 0);*/
+        type            zeroGradient;
+    }
+}
+
+// ************************************************************************* //

tutorials/cfdemSolverPimple/ErgunTestMPI/CFD/0/Us

@@ -0,0 +1,45 @@
+/*--------------------------------*- 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
+{
+    wall
+    {
+        type            zeroGradient;
+        //type            fixedValue;
+        //value           uniform (1 0 0);
+    }
+    inlet
+    {
+        type            zeroGradient;
+        //type            fixedValue;
+        //value           uniform (1 0 0);;
+    }
+    outlet
+    {
+        type            zeroGradient;
+        //type            fixedValue;
+        //value           uniform (1 0 0);
+    }
+}
+
+
+// ************************************************************************* //

tutorials/cfdemSolverPimple/ErgunTestMPI/CFD/0/epsilon

@@ -0,0 +1,48 @@
+/*--------------------------------*- 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      epsilon;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 2 -3 0 0 0 0];
+
+internalField   uniform 0.003;
+
+boundaryField
+{
+    wall
+    {
+        type            epsilonWallFunction;
+        Cmu             0.09;
+        kappa           0.41;
+        E               9.8;
+        value           uniform 0.003;
+    }
+    outlet
+    {
+        type            inletOutlet;
+        inletValue      uniform 0.003;
+        value           uniform 0.003;
+    }
+    inlet
+    {
+        /*type            turbulentMixingLengthDissipationRateInlet;
+        mixingLength    0.005;
+        value           uniform 0.003;*/
+        type            zeroGradient;
+    }
+}
+
+
+// ************************************************************************* //

tutorials/cfdemSolverPimple/ErgunTestMPI/CFD/0/k

@@ -0,0 +1,45 @@
+/*--------------------------------*- 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      k;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 2 -2 0 0 0 0];
+
+internalField   uniform 0.001;
+
+boundaryField
+{
+    wall
+    {
+        type            kqRWallFunction;
+        value           uniform 0.001;
+    }
+    outlet
+    {
+        type            inletOutlet;
+        inletValue      uniform 0.001;
+        value           uniform 0.001;
+    }
+    inlet
+    {
+        /*type            turbulentIntensityKineticEnergyInlet;
+        intensity       0.05;
+        value           uniform 0.001;*/
+        type            zeroGradient;
+    }
+}
+
+
+// ************************************************************************* //

tutorials/cfdemSolverPimple/ErgunTestMPI/CFD/0/nut

@@ -0,0 +1,45 @@
+/*--------------------------------*- 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      nut;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 2 -1 0 0 0 0];
+
+internalField   uniform 0;
+
+boundaryField
+{
+    wall
+    {
+        type            nutWallFunction;
+        Cmu             0.09;
+        kappa           0.41;
+        E               9.8;
+        value           uniform 0;
+    }
+    outlet
+    {
+        type            calculated;
+        value           uniform 0;
+    }
+    inlet
+    {
+        type            calculated;
+        value           uniform 0;
+    }
+}
+
+
+// ************************************************************************* //

tutorials/cfdemSolverPimple/ErgunTestMPI/CFD/0/p

@@ -0,0 +1,43 @@
+/*--------------------------------*- 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      p;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 2 -2 0 0 0 0];
+
+internalField   uniform 1.0e5;
+
+boundaryField
+{
+
+    wall
+    {
+        type            zeroGradient;
+    }
+
+    inlet
+    {
+        type            zeroGradient;
+        //type            fixedValue;
+        //value           uniform 100000;
+    }
+    outlet
+    {
+        //type            zeroGradient;
+        type            fixedValue;
+        value           uniform 100000;
+    }
+}
+
+// ************************************************************************* //

tutorials/cfdemSolverPimple/ErgunTestMPI/CFD/0/rho

@@ -0,0 +1,39 @@
+/*--------------------------------*- 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 10;
+
+boundaryField
+{
+
+    wall
+    {
+        type            zeroGradient;
+    }
+
+    inlet
+    {
+        type            zeroGradient;
+    }
+    outlet
+    {
+        type            zeroGradient;
+    }
+}
+
+// ************************************************************************* //

tutorials/cfdemSolverPimple/ErgunTestMPI/CFD/0/voidfraction

@@ -0,0 +1,41 @@
+/*--------------------------------*- 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
+{
+    wall
+    {
+        type            zeroGradient;
+    }
+    outlet
+    {
+        type            fixedValue;
+        value           uniform 1;
+    }
+    inlet
+    {
+        type            fixedValue;
+        value           uniform 1;
+    }
+}
+
+
+// ************************************************************************* //

tutorials/cfdemSolverPimple/ErgunTestMPI/CFD/constant/RASProperties

@@ -0,0 +1,25 @@
+/*--------------------------------*- 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    "constant";
+    object      RASProperties;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+RASModel        laminar;
+
+turbulence      off;
+
+printCoeffs     on;
+
+
+// ************************************************************************* //

tutorials/cfdemSolverPimple/ErgunTestMPI/CFD/constant/couplingProperties

@@ -0,0 +1,189 @@
+/*---------------------------------------------------------------------------*\
+| =========                 |                                                 |
+| \\      /  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 B; // A or B
+
+couplingInterval 100;
+
+voidFractionModel divided; // might crash with centre;
+
+locateModel engine;
+
+meshMotionModel noMeshMotion;
+
+regionModel allRegion;
+
+IOModel off;
+
+dataExchangeModel twoWayMPI;//twoWayFiles;//oneWayVTK;//
+
+averagingModel dense;//dilute;//
+
+clockModel off;
+
+forceModels
+(
+    //GidaspowDrag
+    DiFeliceDrag
+    //KochHillDrag
+    //gradPForce
+    //viscForce
+    Archimedes
+    //volWeightedAverage
+    //totalMomentumExchange
+);
+
+momCoupleModels
+(
+    implicitCouple
+);
+
+turbulenceModelType RASProperties;//LESProperties;//
+
+//===========================================================================//
+// sub-model properties
+
+implicitCoupleProps
+{
+    velFieldName "U";
+    granVelFieldName "Us";
+    voidfractionFieldName "voidfraction";
+}
+
+ArchimedesProps
+{
+    densityFieldName "rho";
+    gravityFieldName "g";
+}
+gradPForceProps
+{
+    pFieldName "p";
+    densityFieldName "rho";
+    voidfractionFieldName "voidfraction";
+    velocityFieldName "U";
+    interpolation;
+}
+
+viscForceProps
+{
+    velocityFieldName "U";
+    densityFieldName "rho";
+    interpolation;
+}
+volWeightedAverageProps
+{
+    scalarFieldNames
+    (
+        voidfraction
+    );
+    vectorFieldNames
+    (
+    );
+    upperThreshold 0.999;
+    lowerThreshold 0;
+    verbose;
+}
+totalMomentumExchangeProps
+{
+    implicitMomExFieldName "Ksl";
+    explicitMomExFieldName "none";
+    fluidVelFieldName "U";
+    granVelFieldName "Us";
+    densityFieldName "rho";
+}
+GidaspowDragProps
+{
+    velFieldName "U";
+    densityFieldName "rho";
+}
+DiFeliceDragProps
+{
+    velFieldName "U";
+    densityFieldName "rho";
+    voidfractionFieldName "voidfraction";
+}
+
+KochHillDragProps
+{
+    velFieldName "U";
+    densityFieldName "rho";
+}
+
+virtualMassForceProps
+{
+    velFieldName "U";
+    densityFieldName "rho";
+}
+
+oneWayVTKProps
+{
+    couplingFilename "vtk_out%4.4d.vtk";
+    maxNumberOfParticles 30000;
+}
+
+twoWayFilesProps
+{
+    maxNumberOfParticles 10100;
+}
+
+centreProps
+{
+    alphaMin 0.10;
+}
+
+dividedProps
+{
+    alphaMin 0.1;
+    scaleUpVol 1.0;
+}
+
+bigParticleProps
+{
+    alphaMin 0.05;
+    scaleUpVol 10.0;
+    maxCellsPerParticle 1000;
+}
+
+GaussProps
+{
+    maxCellsPerParticle 5000;
+    alphaMin 0.05;
+    scaleUpVol 10.0;
+}
+
+engineProps
+{
+    treeSearch true;
+}
+
+twoWayMPIProps
+{
+    maxNumberOfParticles 10100;
+    liggghtsPath "../DEM/in.liggghts_resume";
+}
+// ************************************************************************* //

tutorials/cfdemSolverPimple/ErgunTestMPI/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/cfdemSolverPimple/ErgunTestMPI/CFD/constant/liggghtsCommands

@@ -0,0 +1,31 @@
+/*---------------------------------------------------------------------------*\
+| =========                 |                                                 |
+| \\      /  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
+);
+
+// ************************************************************************* //

tutorials/cfdemSolverPimple/ErgunTestMPI/CFD/constant/polyMesh/blockMeshDict

@@ -0,0 +1,154 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  2.0                                   |
+|   \\  /    A nd           | Web:      http://www.OpenFOAM.org               |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+version  2.0;
+format   ascii;
+class       dictionary;
+object      blockMeshDict;
+}
+// ************************************
+
+
+
+
+
+
+   meshGenApp blockMesh;
+   convertToMeters 1;
+
+    //64 mm column diameter
+    //15 cm length
+   
+   
+   
+    //Width of middle square section
+   
+   
+   
+   
+    //how many cells in the square section
+    //how many cells from square section to perimeter
+    // how many cells from top to bottom
+
+   vertices
+   (
+    ( 0.0069 0.0069 0.0  ) // Vertex fiveoclocksqb = 0 
+    (-0.0069 0.0069 0.0  ) // Vertex sevenoclocksqb = 1 
+    (-0.0069 -0.0069 0.0 ) // Vertex elevenoclocksqb = 2 
+    ( 0.0069 -0.0069 0.0 ) // Vertex oneoclocksqb = 3 
+   
+    ( 0.00975807358913172 0.00975807357161699 0.0 ) // Vertex fiveoclockcb = 4 
+    (-0.00975807358913172 0.00975807357161699 0.0) // Vertex sevenoclockcb = 5 
+    (-0.00975807358913172 -0.00975807357161699 0.0) // Vertex elevenoclockcb = 6 
+    ( 0.00975807358913172 -0.00975807357161699 0.0) // Vertex oneoclockcb = 7 
+
+    ( 0.0069 0.0069 0.0553) // Vertex fiveoclocksqt = 8 
+    (-0.0069 0.0069 0.0553) // Vertex sevenoclocksqt = 9 
+    (-0.0069 -0.0069 0.0553) // Vertex elevenoclocksqt = 10 
+    ( 0.0069 -0.0069 0.0553) // Vertex oneoclocksqt = 11 
+   
+    ( 0.00975807358913172 0.00975807357161699 0.0553) // Vertex fiveoclockct = 12 
+    (-0.00975807358913172 0.00975807357161699 0.0553) // Vertex sevenoclockct = 13 
+    (-0.00975807358913172 -0.00975807357161699 0.0553) // Vertex elevenoclockct = 14 
+    ( 0.00975807358913172 -0.00975807357161699 0.0553) // Vertex oneoclockct = 15 
+   );				
+
+   blocks
+   (
+    //square block
+    hex (
+       2 3 0 1   
+       10 11 8 9   
+       )
+    (8 8 24)
+    simpleGrading (1 1 1)
+
+    //slice1
+    hex (
+       1 0 4 5
+       9 8 12 13
+       )
+    (8 4 24)
+    simpleGrading (1 1 1)
+
+    //slice2
+    hex (
+       6 2 1  5 
+       14 10 9   13 
+       )
+   ( 4 8 24)
+simpleGrading (1 1 1)
+
+   //slice3
+   hex (
+         6 7 3 2   
+         14 15 11 10   
+       )
+   (8 4 24)
+simpleGrading (1 1 1)
+
+   //slice4
+   hex (
+         3 7 4 0  
+         11 15 12 8  
+       )
+   (4 8 24)
+simpleGrading (1 1 1)
+
+   );
+
+
+   //create the quarter circles
+   edges
+   (
+    arc 4 5 (0.0 0.0138 0.0 )
+    arc 5 6 (-0.0138 0.0 0.0)
+    arc 6 7 (0.0 -0.0138 0.0 )
+    arc 7 4 (0.0138 0.0 0.0)
+
+    arc 12 13 (0.0 0.0138 0.0553 )
+    arc 13 14 (-0.0138 0.0 0.0553 )
+    arc 14 15 (0.0 -0.0138 0.0553 )
+    arc 15 12 (0.0138 0.0 0.0553 )
+
+   );
+
+   patches
+   (
+    patch inlet
+    (
+     (0 3 2 1)
+     (0 4 7 3)
+     (4 0 1 5)
+     (1 2 6 5)
+     (3 7 6 2)
+    )
+
+    patch outlet
+    (
+     (8 11 10 9)
+     (8 12 15 11)
+     (12 8 9 13)
+     (9 10 14 13)
+     (11 15 14 10)
+    )
+
+    wall wall
+    (
+     (5 4 12 13)
+     (5 13 14 6)
+     (6 14 15 7)
+     (7 15 12 4)
+    )
+
+);
+
+mergePatchPairs
+(
+);

tutorials/cfdemSolverPimple/ErgunTestMPI/CFD/constant/transportProperties

@@ -0,0 +1,37 @@
+/*--------------------------------*- 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    "constant";
+    object      transportProperties;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+transportModel  Newtonian;
+
+nu              nu [ 0 2 -1 0 0 0 0 ] 1.5e-04;
+
+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/cfdemSolverPimple/ErgunTestMPI/CFD/constant/turbulenceProperties

@@ -0,0 +1,21 @@
+/*--------------------------------*- 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    "constant";
+    object      turbulenceProperties;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+simulationType  RASModel;
+
+
+// ************************************************************************* //

tutorials/cfdemSolverPimple/ErgunTestMPI/CFD/octave/loaddata.m

@@ -0,0 +1,69 @@
+function [data,colname]=loaddata(filename,columns,headerlines)
+
+% data = loaddata(filename,columns,headerlines)
+% data = loaddata(filename,0,headerlines): automatc detection of 
+%           columns by word count in first headerline or first data line
+%           (if headerline==0), separated by blanks or tabs
+%           columns=-1: write colum assignment code lines
+%              works only if headerline>0 and colnames separated by 
+%              not more than 1 blank!
+%           columns=-2: columns separated by 2 blanks (fluent/scheme auswertung-instat.scm)
+% [data, colname] = loaddata(filename,0,headerlines) 
+%           colname = cell string array of column names from headerline
+
+fprintf(1,'loading %s ...  ', filename);
+f=fopen(filename,'r');
+if f==-1
+  fprintf(1,'\n*** error: could not open "%s" ...\n', filename);
+  data=[];
+else
+  for i=1:headerlines
+    if i==1 s=fgets(f); else fgets(f); end
+  end
+  if headerlines==0
+    s=fgets(f); frewind(f);
+  end
+  if columns<=0 % & headerlines>=1
+    fprintf(1,'\n');
+    fprintf(1,'  %s',s);
+    pos = findstr(sprintf('\t'),s); % trennzeichenpositionen in zeile suchen
+    if length(pos)>0
+      endpos=length(s);
+      %if isspace(s(endpos-1)) endpos=endpos-1; end
+      pos = [0 pos(1,:) endpos]; % anfangs- und endposition hinzuf<EFBFBD>gen
+    elseif columns==-2 % spalteneintr<EFBFBD>ge durch 2 blanks getrennt (fluent/scheme)
+      pos=[0];
+      for i=1:length(s)-1
+        if isspace(s(i)) & isspace(s(i+1))
+          pos(end+1)=i+1;
+        end
+      end
+      pos(end+1)=length(s);      
+    else % spalteneintr<EFBFBD>ge durch blanks getrennt     
+      %pos = findstr(' ',s); % trennzeichenpositionen in zeile suchen
+      pos=[];
+      word=0;
+      for i=1:length(s)
+        if isspace(s(i))
+          word=0;
+        elseif ~word
+          word=1;
+          pos(end+1)=i-1;
+        end
+      end
+      pos(end+1)=length(s);
+    end
+    if headerlines>=1
+      for i=1:length(pos)-1 % alle spalten
+        colname{i}=s(pos(i)+1:pos(i+1)-1);
+        if columns==-1 fprintf(' = data(strmatch(''%s'',colname),:); %% column %d\n', colname{i},i); end
+        if columns==-2 fprintf('%% column %d: %s\n', i, colname{i}); end
+      end   
+    end
+    columns = length(pos)-1;
+    fprintf(1,'  total: %d columns ',columns);
+  end
+  data=fscanf(f,'%f',[columns,inf]);
+  fclose(f);
+  fprintf(1,'done.\n');
+end

tutorials/cfdemSolverPimple/ErgunTestMPI/CFD/octave/totalPressureDrop.m

@@ -0,0 +1,107 @@
+close all;
+clear;
+clc;
+
+%====================================%
+% simulation data 1
+%====================================%
+%path = '../probes/0/p';
+path = '../probes/0/p';
+columns=22;
+headerlines=4;
+data = loaddata(path,columns,headerlines);
+data=transpose(data);
+[x,y]=size(data)
+dp_sim = (data(:,2)-data(:,y))/10000;
+t_sim = data(:,1);
+%fprintf('final pressureDrop of sim = %f Pa\n',dp_sim(length(dp_sim)) )
+
+%====================================%
+% analytical calculation
+%====================================%
+
+%===================
+% Ergun Equation
+%===================
+fprintf('\ncalc Ergun eqn:\n')
+dp = 0.001			% particle diameter
+phip = 1			% sphericity
+epsilon = 0.451335              % void fraction
+Ustart = 0.002
+Uend = 0.02
+timeStepSize = 0.001;            % time interval of pressure data
+Tstart = 0;
+Tend = t_sim(length(t_sim));
+deltaU=(Uend-Ustart)/((Tend-Tstart)/timeStepSize);
+U = Ustart+deltaU:deltaU:Uend;  % velocity over time
+Ua = U / epsilon;		% physical velocity
+L = 0.0156			% length of bed
+rhoG = 10			% density in kg/m3
+nuG = 1.5*10^-4			% kinemat Visk in m2/s
+muG = nuG*rhoG			% dynam visc in Pa s
+
+dpErgun= L * (
+                150*((1-epsilon)^2/epsilon^3)*((muG.*U)/(phip*dp)^2) 
+              +1.75*((1-epsilon)/epsilon^3)*((rhoG.*U.^2)/(phip*dp))
+        )/10000/rhoG;
+
+fprintf('NOTE: this pressure is divided by density (according to CFD solver)\n')
+fprintf('so the result does not depend on density\n')
+
+%fprintf('final pressure drop (Ergun eqn)= %f Pa\n',dpErgun)
+
+%==================================
+% min fluidization velocity in m/s
+%==================================
+rhoP = 2000                      % particle density in kg/m3
+g = 9.81                        % gravity m/s2
+
+Umf = dp^2*(rhoP-rhoG)*g/(150*muG)*(epsilon^3*phip^2)/(1-epsilon);
+ReMF = Umf*dp*rhoG/muG;
+if(ReMF<20)
+    fprintf('applying eqn1 for Umf.\n')
+elseif(ReMF>20 && ReMF<1000)
+    fprintf('applying eqn1 for Umf.\n')
+elseif (ReMF>=1000)
+    fprintf('applying eqn2 for Umf.\n')
+    Umf = sqrt(dp*(rhoP-rhoG)*g/(1.75*rhoG)*epsilon^3*phip);
+    ReMF = Umf*dp*rhoG/muG;
+end
+Umf
+ReMF
+
+dpUmf= L * (
+                150*((1-epsilon)^2/epsilon^3)*((muG.*Umf)/(phip*dp)^2) 
+              +1.75*((1-epsilon)/epsilon^3)*((rhoG.*Umf.^2)/(phip*dp))
+        )/10000/rhoG;
+
+%====================================%
+% plot data
+%====================================%
+length(U)
+length(dp_sim)
+figure(2)
+plot(U,dp_sim)
+title("Ergun pressure drop vs. simulation")
+a=strcat("analytical (Ergun), Umf=",num2str(Umf),", dpUmf=",num2str(dpUmf));
+legend(a,"simulation")
+xlabel("velocity in [m/s]")
+ylabel("pressure drop [bar]")
+axis([0,Uend,0,dpErgun(length(dpErgun))])
+
+figure(1)
+plot(U,dpErgun,U,dp_sim,[Umf,Uend],dpUmf*ones(1,2))
+title("Ergun pressure drop vs. simulation")
+a=strcat("analytical (Ergun), Umf=",num2str(Umf),", dpUmf=",num2str(dpUmf));
+legend(a,"simulation")
+xlabel("velocity in [m/s]")
+ylabel("pressure drop [bar]")
+axis([0,Uend,0,dpErgun(length(dpErgun))])
+
+%print('cfdemSolverPiso_settlingTest.eps','-deps2')
+print -color "cfdemSolverPimple_ErgunTestMPI.eps"
+replot;
+
+
+
+

tutorials/cfdemSolverPimple/ErgunTestMPI/CFD/steps_0p1s

@@ -0,0 +1,22 @@
+(
+ (0.000 (0 0 0.002))
+ (0.010 (0 0 0.002))
+ (0.011 (0 0 0.004))
+ (0.020 (0 0 0.004))
+ (0.021 (0 0 0.006))
+ (0.030 (0 0 0.006))
+ (0.031 (0 0 0.008))
+ (0.040 (0 0 0.008))
+ (0.041 (0 0 0.010))
+ (0.050 (0 0 0.010))
+ (0.051 (0 0 0.012))
+ (0.060 (0 0 0.012))
+ (0.061 (0 0 0.014))
+ (0.070 (0 0 0.014))
+ (0.071 (0 0 0.016))
+ (0.080 (0 0 0.016))
+ (0.081 (0 0 0.018))
+ (0.090 (0 0 0.018))
+ (0.091 (0 0 0.020))
+ (0.100 (0 0 0.020))
+)

tutorials/cfdemSolverPimple/ErgunTestMPI/CFD/system/controlDict

@@ -0,0 +1,117 @@
+/*--------------------------------*- 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      controlDict;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+application     pisoFoam;
+
+startFrom       startTime;
+
+startTime       0;
+
+stopAt          endTime;
+
+endTime         0.1;//0.01;
+
+deltaT          0.001;
+
+writeControl    adjustableRunTime;
+
+writeInterval   0.001;//0.01;
+
+purgeWrite      0;
+
+writeFormat     ascii;
+
+writePrecision  6;
+
+writeCompression uncompressed;
+
+timeFormat      general;
+
+timePrecision   6;
+
+runTimeModifiable yes;
+
+adjustTimeStep  no;
+
+maxCo 0.1;
+
+libs ("libOpenFOAM.so" "libgroovyBC.so");
+
+functions
+(
+
+    probes
+    {
+        type        probes;
+        // Where to load it from
+        functionObjectLibs ( "libsampling.so" );  
+        // Name of the directory for probe data
+        name        probes;
+        probeLocations
+        (
+            (0 0 0.0001)
+            (0 0 0.0026)
+            (0 0 0.0051)
+            (0 0 0.0076)
+            (0 0 0.0101)
+            (0 0 0.0126)
+            (0 0 0.0151)
+            (0 0 0.0176)
+            (0 0 0.0201)
+            (0 0 0.0226)
+            (0 0 0.0251)
+            (0 0 0.0276)
+            (0 0 0.0301)
+            (0 0 0.0326)
+            (0 0 0.0351)
+            (0 0 0.0375)
+            (0 0 0.0401)
+            (0 0 0.0426)
+            (0 0 0.0451)
+            (0 0 0.0476)
+            (0 0 0.0529)
+        );
+
+        // Fields to be probed
+        fields ( p U voidfraction volAverage_voidfraction);
+
+        // Write at same frequency as fields
+        outputControl   timeStep;//outputTime;
+        outputInterval  1;
+    }
+
+   /*pressureDrop
+   {
+     type patchAverage;
+     functionObjectLibs
+     (
+         "libsimpleFunctionObjects.so"
+     );
+     verbose true;
+     patches
+     (
+         inlet
+         outlet
+     );
+     fields
+     (
+         p
+     );
+     factor  1;
+   }*/
+);
+// ************************************************************************* //

tutorials/cfdemSolverPimple/ErgunTestMPI/CFD/system/controlDict.foam

@@ -0,0 +1,117 @@
+/*--------------------------------*- 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      controlDict;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+application     pisoFoam;
+
+startFrom       startTime;
+
+startTime       0;
+
+stopAt          endTime;
+
+endTime         0.1;//0.01;
+
+deltaT          0.001;
+
+writeControl    adjustableRunTime;
+
+writeInterval   0.001;//0.01;
+
+purgeWrite      0;
+
+writeFormat     ascii;
+
+writePrecision  6;
+
+writeCompression uncompressed;
+
+timeFormat      general;
+
+timePrecision   6;
+
+runTimeModifiable yes;
+
+adjustTimeStep  no;
+
+maxCo 0.1;
+
+libs ("libOpenFOAM.so" "libgroovyBC.so");
+
+functions
+(
+
+    probes
+    {
+        type        probes;
+        // Where to load it from
+        functionObjectLibs ( "libsampling.so" );  
+        // Name of the directory for probe data
+        name        probes;
+        probeLocations
+        (
+            (0 0 0.0001)
+            (0 0 0.0026)
+            (0 0 0.0051)
+            (0 0 0.0076)
+            (0 0 0.0101)
+            (0 0 0.0126)
+            (0 0 0.0151)
+            (0 0 0.0176)
+            (0 0 0.0201)
+            (0 0 0.0226)
+            (0 0 0.0251)
+            (0 0 0.0276)
+            (0 0 0.0301)
+            (0 0 0.0326)
+            (0 0 0.0351)
+            (0 0 0.0375)
+            (0 0 0.0401)
+            (0 0 0.0426)
+            (0 0 0.0451)
+            (0 0 0.0476)
+            (0 0 0.0529)
+        );
+
+        // Fields to be probed
+        fields ( p U voidfraction volAverage_voidfraction);
+
+        // Write at same frequency as fields
+        outputControl   timeStep;//outputTime;
+        outputInterval  1;
+    }
+
+   /*pressureDrop
+   {
+     type patchAverage;
+     functionObjectLibs
+     (
+         "libsimpleFunctionObjects.so"
+     );
+     verbose true;
+     patches
+     (
+         inlet
+         outlet
+     );
+     fields
+     (
+         p
+     );
+     factor  1;
+   }*/
+);
+// ************************************************************************* //

tutorials/cfdemSolverPimple/ErgunTestMPI/CFD/system/decomposeParDict

@@ -0,0 +1,86 @@
+/*-------------------------------*- 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  2;
+
+//- Keep owner and neighbour on same processor for faces in zones:
+// preserveFaceZones (heater solid1 solid3);
+
+//method          scotch;
+// method          hierarchical;
+ method          simple;
+// method          metis;
+// method          manual;
+
+simpleCoeffs
+{
+    n           (1 1 2);
+    delta       0.001;
+}
+
+hierarchicalCoeffs
+{
+    n           (2 2 1);
+    delta       0.001;
+    order       xyz;
+}
+
+metisCoeffs
+{
+ /*
+    processorWeights
+    (
+        1
+        1
+        1
+        1
+    );
+  */
+}
+
+scotchCoeffs
+{
+    //processorWeights
+    //(
+    //    1
+    //    1
+    //    1
+    //    1
+    //);
+    //writeGraph  true;
+    //strategy "b";
+}
+
+manualCoeffs
+{
+    dataFile    "decompositionData";
+}
+
+
+//// Is the case distributed
+//distributed     yes;
+//// Per slave (so nProcs-1 entries) the directory above the case.
+//roots           
+//(
+//    "/tmp"
+//    "/tmp"
+//);
+
+
+// ************************************************************************* //

tutorials/cfdemSolverPimple/ErgunTestMPI/CFD/system/funkySetFieldsDict

@@ -0,0 +1,33 @@
+FoamFile
+{
+    version         2.0;
+    format          ascii;
+
+    root            "/.automount/werner/Werner/bgschaid/bgschaid-foamStuff/DepositionTests";
+    case            "neutralTest";
+    instance        "system";
+    local           "";
+
+    class           dictionary;
+    object          funkySetFieldsDict;
+}
+
+expressions
+(
+	init
+	{
+		field voidfraction;
+		expression "1";
+                valuePatches 2(inlet outlet); // set value 1 at those patches
+	}
+	setGamma
+	{
+		field voidfraction;
+		expression "0.43";
+		//condition "pow(pow(pos().x-0.2,2) + pow(pos().y-0,2)+ pow(pos().z-0,2),0.5) < 0.04";
+                condition "pos().z < 0.01467 && pos().z > 0.0";
+                keepPatches true;
+                //valuePatches 1(inlet); // set value 1 at those patches
+	}
+	
+);

tutorials/cfdemSolverPimple/ErgunTestMPI/CFD/system/fvSchemes

@@ -0,0 +1,76 @@
+/*--------------------------------*- 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)         Gauss linear;
+    grad(U)         Gauss linear;
+}
+
+divSchemes
+{
+    default         Gauss linear;
+    div(phi,U)      Gauss limitedLinearV 1;
+    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;
+}
+
+interpolationSchemes
+{
+    default         linear;
+    interpolate(U)  linear;
+}
+
+snGradSchemes
+{
+    default         corrected;
+}
+
+fluxRequired
+{
+    default         no;
+    p               ;
+}
+
+
+// ************************************************************************* //

tutorials/cfdemSolverPimple/ErgunTestMPI/CFD/system/fvSolution

@@ -0,0 +1,75 @@
+/*--------------------------------*- 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      fvSolution;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+solvers
+{
+    p
+    {
+        solver          GAMG;
+        tolerance       1e-06;
+        relTol          0.01;
+        smoother        GaussSeidel;
+        cacheAgglomeration true;
+        nCellsInCoarsestLevel 10;
+        agglomerator    faceAreaPair;
+        mergeLevels     1;
+    }
+
+    pFinal
+    {
+        solver          GAMG;
+        tolerance       1e-06;
+        relTol          0;
+        smoother        GaussSeidel;
+        cacheAgglomeration true;
+        nCellsInCoarsestLevel 10;
+        agglomerator    faceAreaPair;
+        mergeLevels     1;
+    }
+
+    "(U|k|epsilon)"
+    {
+        solver          PBiCG;
+        preconditioner  DILU;
+        tolerance       1e-05;
+        relTol          0.1;
+    }
+
+    "(U|k|epsilon)Final"
+    {
+        $U;
+        tolerance       1e-05;
+        relTol          0;
+    }
+}
+
+PIMPLE
+{
+    nOuterCorrectors 2;
+    nCorrectors     2;
+    nNonOrthogonalCorrectors 0;
+    pRefCell        0;
+    pRefValue       0;
+}
+
+relaxationFactors
+{
+    "U.*"           1;
+    "k.*"           1;
+    "epsilon.*"     1;
+}
+
+// ************************************************************************* //

tutorials/cfdemSolverPimple/ErgunTestMPI/DEM/in.liggghts_init

@@ -0,0 +1,64 @@
+# Pour granular particles into chute container, then induce flow
+echo both 
+atom_style	granular
+atom_modify	map array
+communicate	single vel yes
+
+boundary	m m m
+newton		off
+
+units		si
+
+region		reg block -0.015 0.015 -0.015 0.015 -0.001 0.0554 units box
+create_box	1 reg
+
+neighbor	0.001 bin
+neigh_modify	delay 0
+
+
+#Material properties required for new pair styles
+
+fix 		m1 all property/global youngsModulus peratomtype 5.e6
+fix 		m2 all property/global poissonsRatio peratomtype 0.45
+fix 		m3 all property/global coefficientRestitution peratomtypepair 1 0.3
+fix 		m4 all property/global coefficientFriction peratomtypepair 1 0.5
+
+#pair style
+pair_style 	gran/hertz/history #Hertzian without cohesion
+pair_coeff	* *
+
+#timestep, gravity
+timestep	0.00001
+fix		gravi all gravity 9.81 vector 0.0 0.0 -1.0
+
+#walls
+fix			zwalls1 all wall/gran/hertz/history   primitive type 1  zplane 0.0
+fix			zwalls2 all wall/gran/hertz/history   primitive type 1  zplane 0.0553
+fix         cylwalls all wall/gran/hertz/history primitive type 1 zcylinder 0.01385 0. 0. tangential_damping off
+
+#particle distributions and insertion
+region		bc cylinder z 0.0 0.0 0.012 0. 0.055 units box
+fix		    pts1 all particletemplate/sphere 1 atom_type 1 density constant 200 radius constant 0.0005
+fix		    pdd1 all particledistribution/discrete 1.  1 pts1 1.0
+
+fix		ins all insert/pack seed 100001 distributiontemplate pdd1 vel constant 0. 0. -1. insert_every once overlapcheck yes all_in yes particles_in_region 10000 region bc
+
+#apply nve integration to all particles that are inserted as single particles
+fix		integr all nve/sphere
+
+#screen output
+compute		1 all erotate/sphere
+thermo_style	custom step atoms ke c_1 vol
+thermo		1000
+thermo_modify	lost ignore norm no
+compute_modify	thermo_temp dynamic yes
+
+#insert the first particles so that dump is not empty
+run		1
+dump		dmp all custom 5000 post/dump.liggghts_init id type type x y z ix iy iz vx vy vz fx fy fz omegax omegay omegaz radius 
+#force : f_couple_cfd[0] f_couple_cfd[1] f_couple_cfd[2]
+#node : f_couple_cfd[6]
+#cell id : f_couple_cfd[7]
+
+run		10000 upto
+write_restart 	liggghts.restart

tutorials/cfdemSolverPimple/ErgunTestMPI/DEM/in.liggghts_resume

@@ -0,0 +1,71 @@
+# Pour granular particles into chute container, then induce flow
+
+atom_style	granular
+atom_modify	map array
+communicate	single vel yes
+
+boundary	m m m
+newton		off
+
+units		si
+processors	2 1 1
+
+#read the restart file
+read_restart 	../DEM/liggghts.restart
+
+#do not do this here, the simulation box is in the restart file!
+#region		reg block -0.015 0.015 -0.015 0.015 -0.001 0.0554 units box
+#create_box	1 reg
+
+neighbor	0.0005 bin
+neigh_modify	delay 0
+
+#Material properties required for new pair styles
+
+fix 		m1 all property/global youngsModulus peratomtype 5.e6
+fix 		m2 all property/global poissonsRatio peratomtype 0.45
+fix 		m3 all property/global coefficientRestitution peratomtypepair 1 0.3
+fix 		m4 all property/global coefficientFriction peratomtypepair 1 0.5
+
+#pair style
+pair_style 	gran/hertz/history #Hertzian without cohesion
+pair_coeff	* *
+
+#timestep, gravity
+timestep	0.00001
+fix		gravi all gravity 9.81 vector 0.0 0.0 -1.0
+
+fix			zwalls1 all wall/gran/hertz/history   primitive type 1  zplane 0.0
+fix			zwalls2 all wall/gran/hertz/history   primitive type 1  zplane 0.0553
+fix         cylwalls all wall/gran/hertz/history primitive type 1 zcylinder 0.01385 0. 0. tangential_damping off
+
+# change the particles density
+set             group all density 2000
+
+#cfd coupling
+fix		cfd all couple/cfd couple_every 100 mpi
+fix		cfd2 all couple/cfd/force
+
+#apply nve integration to all particles that are inserted as single particles
+fix		integr all nve/sphere
+
+#center of mass
+compute         centerOfMass all com
+
+#compute total dragforce
+compute		dragtotal all reduce sum f_dragforce[1] f_dragforce[2] f_dragforce[3]
+
+#screen output
+compute		1 all erotate/sphere
+thermo_style	custom step atoms ke c_1 vol c_centerOfMass[3] c_dragtotal[1] c_dragtotal[2] c_dragtotal[3]
+thermo		10
+thermo_modify	lost ignore norm no
+compute_modify	thermo_temp dynamic yes
+
+#insert the first particles so that dump is not empty
+dump 		myDump all stl 1 post/dump_*.stl
+#run		1
+dump		dmp all custom 5000 ../DEM/post/dump*.liggghts_restart id type type x y z vx vy vz fx fy fz f_dragforce[1] f_dragforce[2] f_dragforce[3] radius 
+undump		myDump
+
+run		1

tutorials/cfdemSolverPimple/ErgunTestMPI/DEM/post/dummy

@@ -0,0 +1,2 @@
+dummyfile
+dummyfile

tutorials/cfdemSolverPimple/ErgunTestMPI/parCFDDEMrun.sh

@@ -0,0 +1,64 @@
+#!/bin/bash
+
+#===================================================================#
+# allrun script for testcase as part of test routine 
+# run settlingTest CFD part
+# Christoph Goniva - March 2012
+#===================================================================#
+
+#- source CFDEM env vars
+. ~/.bashrc
+
+#- include functions
+source $CFDEM_SRC_DIR/etc/functions.sh
+
+#--------------------------------------------------------------------------------#
+#- define variables
+casePath="$(dirname "$(readlink -f ${BASH_SOURCE[0]})")"
+logpath=$casePath
+headerText="run_parallel_cfdemSolverPimple_ErgunTestMPI_CFDDEM"
+logfileName="log_$headerText"
+solverName="cfdemSolverPimple"
+nrProcs="2"
+machineFileName="none"   # yourMachinefileName | none
+debugMode="off"          # on | off
+testHarnessPath="$CFDEM_TEST_HARNESS_PATH"
+#--------------------------------------------------------------------------------#
+
+#- call function to run a parallel CFD-DEM case
+parCFDDEMrun $logpath $logfileName $casePath $headerText $solverName $nrProcs $machineFileName $debugMode
+
+#------------------------------#
+#  octave
+
+#- change path
+cd octave
+
+#- rmove old graph
+rm cfdemSolverPiso_ErgunTestMPI.eps
+
+#- run octave
+octave totalPressureDrop.m
+
+#- show plot 
+evince cfdemSolverPimple_ErgunTestMPI.eps
+#------------------------------#
+
+#- copy log file to test harness
+cp ../../$logfileName $testHarnessPath
+cp cfdemSolverPimple_ErgunTestMPI.eps $testHarnessPath
+
+#- clean up case
+cd ..
+rm -rf 0.*
+rm -rf processor*
+rm -rf patchAverage_pressureDrop
+rm -rf probes
+rm log.liggghts
+rm ../DEM/post/dump*
+rm -rf particles
+
+#- preserve post directory
+echo "dummyfile" >> $casePath/DEM/post/dummy
+
+

tutorials/cfdemSolverPiso/settlingTestMPI/CFD/constant/couplingProperties

@@ -38,7 +38,7 @@ meshMotionModel noMeshMotion;
 
 regionModel allRegion;
 
-IOModel trackIO; //basicIO; // 
+IOModel basicIO; //trackIO; // 
 
 dataExchangeModel twoWayMPI;//twoWayFiles;//oneWayVTK;//
 

tutorials/cfdemSolverPiso/settlingTestMPI/parCFDDEMrun.sh

@@ -88,7 +88,7 @@ rm -r $casePath/CFD/couplingFiles/*
 rm -r $casePath/DEM/post/*
 rm -r $casePath/DEM/log.*
 rm -r $casePath/CFD/probes
-rm -r $casePath/CFD/particles
+rm -r $casePath/CFD/lagrangian
 rm -r $casePath/CFD/clockData
 echo "done"
 

tutorials/cfdemSolverPisoMS/ErgunTestMPI/DEM/post/dummy

@@ -2,3 +2,4 @@ dummyfile
 dummyfile
 dummyfile
 dummyfile
+dummyfile