Merge branch 'develop' of github.com:ParticulateFlow/CFDEMcoupling into develop

applications/solvers/cfdemSolverMultiphaseScalar/Allwclean

@@ -2,7 +2,7 @@
 cd ${0%/*} || exit 1    # Run from this directory
 set -x
 
-wclean libso multiphaseMixture
+wclean libso multiphaseMixtureScalar
 wclean
 
 #------------------------------------------------------------------------------

applications/solvers/cfdemSolverMultiphaseScalar/Allwmake

@@ -6,7 +6,7 @@ targetType=libso
 . $WM_PROJECT_DIR/wmake/scripts/AllwmakeParseArguments
 set -x
 
-wmake $targetType multiphaseMixture
+wmake $targetType multiphaseMixtureScalar
 wmake
 
 #------------------------------------------------------------------------------

applications/solvers/cfdemSolverMultiphaseScalar/Make/options

@@ -6,7 +6,7 @@ include $(CFDEM_ADD_LIBS_DIR)/additionalLibs
 EXE_INC = \
      $(PFLAGS) \
     -I$(CFDEM_OFVERSION_DIR) \
-    -ImultiphaseMixture/lnInclude \
+    -ImultiphaseMixtureScalar/lnInclude \
     -I$(LIB_SRC)/transportModels \
     -I$(LIB_SRC)/transportModels/incompressible/lnInclude \
     -I$(LIB_SRC)/transportModels/interfaceProperties/lnInclude \

applications/solvers/cfdemSolverMultiphaseScalar/cfdemSolverMultiphaseScalar.C

@@ -30,7 +30,7 @@ Description
 \*---------------------------------------------------------------------------*/
 
 #include "fvCFD.H"
-#include "multiphaseMixture.H"
+#include "multiphaseMixtureScalar.H"
 #include "turbulentTransportModel.H"
 #include "pimpleControl.H"
 #include "fvOptions.H"

applications/solvers/cfdemSolverMultiphaseScalar/createFields.H

@@ -88,7 +88,7 @@ surfaceScalarField phi
      linearInterpolate(U*voidfraction) & mesh.Sf()
  );
 
-multiphaseMixture mixture(U, phi, voidfraction);
+multiphaseMixtureScalar mixture(U, phi, voidfraction);
 
 // Need to store rho for ddt(rho, U)
 volScalarField rho

applications/solvers/cfdemSolverMultiphaseScalar/multiphaseMixtureScalar/Make/files

@@ -1,5 +1,5 @@
 phase/phase.C
 alphaContactAngle/alphaContactAngleFvPatchScalarField.C
-multiphaseMixture.C
+multiphaseMixtureScalar.C
 
 LIB = $(CFDEM_LIB_DIR)/libcfdemMultiphaseInterFoamScalar

applications/solvers/cfdemSolverMultiphaseScalar/multiphaseMixtureScalar/alphaContactAngle/alphaContactAngleFvPatchScalarField.H

@@ -26,7 +26,7 @@ Class
 
 Description
     Contact-angle boundary condition for multi-phase interface-capturing
-    simulations.  Used in conjuction with multiphaseMixture.
+    simulations.  Used in conjuction with multiphaseMixtureScalar.
 
 SourceFiles
     alphaContactAngleFvPatchScalarField.C
@@ -37,7 +37,7 @@ SourceFiles
 #define alphaContactAngleFvPatchScalarField_H
 
 #include "zeroGradientFvPatchFields.H"
-#include "multiphaseMixture.H"
+#include "multiphaseMixtureScalar.H"
 
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
@@ -117,8 +117,8 @@ public:
     typedef HashTable
     <
         interfaceThetaProps,
-        multiphaseMixture::interfacePair,
-        multiphaseMixture::interfacePair::hash
+        multiphaseMixtureScalar::interfacePair,
+        multiphaseMixtureScalar::interfacePair::hash
     > thetaPropsTable;
 
 

applications/solvers/cfdemSolverMultiphaseScalar/multiphaseMixtureScalar/multiphaseMixtureScalar.C

@@ -18,7 +18,7 @@ License
 
 \*---------------------------------------------------------------------------*/
 
-#include "multiphaseMixture.H"
+#include "multiphaseMixtureScalar.H"
 #include "alphaContactAngleFvPatchScalarField.H"
 #include "Time.H"
 #include "subCycle.H"
@@ -31,13 +31,13 @@ License
 
 // * * * * * * * * * * * * * * * Static Member Data  * * * * * * * * * * * * //
 
-const Foam::scalar Foam::multiphaseMixture::convertToRad =
+const Foam::scalar Foam::multiphaseMixtureScalar::convertToRad =
     Foam::constant::mathematical::pi/180.0;
 
 
 // * * * * * * * * * * * * * Private Member Functions  * * * * * * * * * * * //
 
-void Foam::multiphaseMixture::calcAlphas()
+void Foam::multiphaseMixtureScalar::calcAlphas()
 {
     scalar level = 0.0;
     alphas_ == 0.0;
@@ -51,7 +51,7 @@ void Foam::multiphaseMixture::calcAlphas()
 
 
 Foam::tmp<Foam::volScalarField>
-Foam::multiphaseMixture::calcNu() const
+Foam::multiphaseMixtureScalar::calcNu() const
 {
     PtrDictionary<phase>::const_iterator iter = phases_.begin();
 
@@ -74,7 +74,7 @@ Foam::multiphaseMixture::calcNu() const
 }
 
 Foam::tmp<Foam::surfaceScalarField>
-Foam::multiphaseMixture::calcStf() const
+Foam::multiphaseMixtureScalar::calcStf() const
 {
     tmp<surfaceScalarField> tstf
     (
@@ -134,7 +134,7 @@ Foam::multiphaseMixture::calcStf() const
 
 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
 
-Foam::multiphaseMixture::multiphaseMixture
+Foam::multiphaseMixtureScalar::multiphaseMixtureScalar
 (
     const volVectorField& U,
     const surfaceScalarField& phi,
@@ -230,7 +230,7 @@ Foam::multiphaseMixture::multiphaseMixture
 // * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * * //
 
 Foam::tmp<Foam::volScalarField>
-Foam::multiphaseMixture::rho() const
+Foam::multiphaseMixtureScalar::rho() const
 {
     PtrDictionary<phase>::const_iterator iter = phases_.begin();
 
@@ -247,7 +247,7 @@ Foam::multiphaseMixture::rho() const
 
 
 Foam::tmp<Foam::scalarField>
-Foam::multiphaseMixture::rho(const label patchi) const
+Foam::multiphaseMixtureScalar::rho(const label patchi) const
 {
     PtrDictionary<phase>::const_iterator iter = phases_.begin();
 
@@ -264,9 +264,9 @@ Foam::multiphaseMixture::rho(const label patchi) const
 
 
 Foam::tmp<Foam::volScalarField>
-Foam::multiphaseMixture::mu() const
+Foam::multiphaseMixtureScalar::mu() const
 {
-    Info << "In multiphasemixture mu()" << endl;
+    Info << "In multiphaseMixtureScalar mu()" << endl;
     return rho()*nu();
 //    PtrDictionary<phase>::const_iterator iter = phases_.begin();
 
@@ -283,7 +283,7 @@ Foam::multiphaseMixture::mu() const
 
 
 Foam::tmp<Foam::scalarField>
-Foam::multiphaseMixture::mu(const label patchi) const
+Foam::multiphaseMixtureScalar::mu(const label patchi) const
 {
     PtrDictionary<phase>::const_iterator iter = phases_.begin();
 
@@ -306,7 +306,7 @@ Foam::multiphaseMixture::mu(const label patchi) const
 
 
 Foam::tmp<Foam::surfaceScalarField>
-Foam::multiphaseMixture::muf() const
+Foam::multiphaseMixtureScalar::muf() const
 {
 
     return nuf()*fvc::interpolate(rho());
@@ -327,13 +327,13 @@ Foam::multiphaseMixture::muf() const
 
 
 Foam::tmp<Foam::volScalarField>
-Foam::multiphaseMixture::nu() const
+Foam::multiphaseMixtureScalar::nu() const
 {
     return nu_;
 }
 
 Foam::tmp<Foam::scalarField>
-Foam::multiphaseMixture::nu(const label patchi) const
+Foam::multiphaseMixtureScalar::nu(const label patchi) const
 {
     //return nu_.boundaryField()[patchi];
     PtrDictionary<phase>::const_iterator iter = phases_.begin();
@@ -355,7 +355,7 @@ Foam::multiphaseMixture::nu(const label patchi) const
 
 
 Foam::tmp<Foam::surfaceScalarField>
-Foam::multiphaseMixture::nuf() const
+Foam::multiphaseMixtureScalar::nuf() const
 {
     //return muf()/fvc::interpolate(rho());
     PtrDictionary<phase>::const_iterator iter = phases_.begin();
@@ -374,7 +374,7 @@ Foam::multiphaseMixture::nuf() const
 }
 
 Foam::tmp<Foam::volScalarField>
-Foam::multiphaseMixture::Cp() const
+Foam::multiphaseMixtureScalar::Cp() const
 {
     PtrDictionary<phase>::const_iterator iter = phases_.begin();
 
@@ -395,7 +395,7 @@ Foam::multiphaseMixture::Cp() const
 }
 
 Foam::tmp<Foam::volScalarField>
-Foam::multiphaseMixture::kf() const
+Foam::multiphaseMixtureScalar::kf() const
 {
     PtrDictionary<phase>::const_iterator iter = phases_.begin();
 
@@ -417,7 +417,7 @@ Foam::multiphaseMixture::kf() const
 }
 
 Foam::tmp<Foam::volScalarField>
-Foam::multiphaseMixture::D() const
+Foam::multiphaseMixtureScalar::D() const
 {
     PtrDictionary<phase>::const_iterator iter = phases_.begin();
 
@@ -439,7 +439,7 @@ Foam::multiphaseMixture::D() const
 }
 
 Foam::tmp<Foam::volScalarField>
-Foam::multiphaseMixture::Cs() const
+Foam::multiphaseMixtureScalar::Cs() const
 {
     PtrDictionary<phase>::const_iterator iter = phases_.begin();
 
@@ -456,7 +456,7 @@ Foam::multiphaseMixture::Cs() const
 }
 
 Foam::tmp<Foam::surfaceScalarField>
-Foam::multiphaseMixture::diffusionCorrection() const
+Foam::multiphaseMixtureScalar::diffusionCorrection() const
 {
 
     surfaceScalarField numerator
@@ -517,7 +517,7 @@ Foam::multiphaseMixture::diffusionCorrection() const
     return correction;
 }
 
-void Foam::multiphaseMixture::solve()
+void Foam::multiphaseMixtureScalar::solve()
 {
     correct();
 
@@ -570,7 +570,7 @@ void Foam::multiphaseMixture::solve()
 }
 
 
-void Foam::multiphaseMixture::correct()
+void Foam::multiphaseMixtureScalar::correct()
 {
     forAllIter(PtrDictionary<phase>, phases_, iter)
     {
@@ -579,7 +579,7 @@ void Foam::multiphaseMixture::correct()
 }
 
 
-Foam::tmp<Foam::surfaceVectorField> Foam::multiphaseMixture::nHatfv
+Foam::tmp<Foam::surfaceVectorField> Foam::multiphaseMixtureScalar::nHatfv
 (
     const volScalarField& alpha1,
     const volScalarField& alpha2
@@ -605,7 +605,7 @@ Foam::tmp<Foam::surfaceVectorField> Foam::multiphaseMixture::nHatfv
 }
 
 
-Foam::tmp<Foam::surfaceScalarField> Foam::multiphaseMixture::nHatf
+Foam::tmp<Foam::surfaceScalarField> Foam::multiphaseMixtureScalar::nHatf
 (
     const volScalarField& alpha1,
     const volScalarField& alpha2
@@ -622,7 +622,7 @@ Foam::tmp<Foam::surfaceScalarField> Foam::multiphaseMixture::nHatf
 // The dynamic contact angle is calculated from the component of the
 // velocity on the direction of the interface, parallel to the wall.
 
-void Foam::multiphaseMixture::correctContactAngle
+void Foam::multiphaseMixtureScalar::correctContactAngle
 (
     const phase& alpha1,
     const phase& alpha2,
@@ -726,7 +726,7 @@ void Foam::multiphaseMixture::correctContactAngle
 }
 
 
-Foam::tmp<Foam::volScalarField> Foam::multiphaseMixture::K
+Foam::tmp<Foam::volScalarField> Foam::multiphaseMixtureScalar::K
 (
     const phase& alpha1,
     const phase& alpha2
@@ -742,7 +742,7 @@ Foam::tmp<Foam::volScalarField> Foam::multiphaseMixture::K
 
 
 Foam::tmp<Foam::volScalarField>
-Foam::multiphaseMixture::nearInterface() const
+Foam::multiphaseMixtureScalar::nearInterface() const
 {
     tmp<volScalarField> tnearInt
     (
@@ -768,7 +768,7 @@ Foam::multiphaseMixture::nearInterface() const
 }
 
 
-void Foam::multiphaseMixture::solveAlphas
+void Foam::multiphaseMixtureScalar::solveAlphas
 (
     const scalar cAlpha
 )
@@ -901,7 +901,7 @@ void Foam::multiphaseMixture::solveAlphas
 }
 
 
-bool Foam::multiphaseMixture::read()
+bool Foam::multiphaseMixtureScalar::read()
 {
     if (transportModel::read())
     {

applications/solvers/cfdemSolverMultiphaseScalar/multiphaseMixtureScalar/multiphaseMixtureScalar.H

@@ -17,10 +17,10 @@ License
     Copyright (C) 2018- Mathias Vångö, JKU Linz, Austria
 
 Class
-    multiphaseMixture
+    multiphaseMixtureScalar
 
 Description
-    This class is based on the OpenFOAM(R) Foam::multiphaseMixture class,
+    This class is based on the OpenFOAM(R) Foam::multiphaseMixtureScalar class,
     which is an incompressible multi-phase mixture with built in solution
     for the phase fractions with interface compression for interface-capturing.
     It has been extended to include the void fraction in the volume fraction
@@ -33,11 +33,11 @@ Description
     between each phase-pair.
 
 SourceFiles
-    multiphaseMixture.C
+    multiphaseMixtureScalar.C
 \*---------------------------------------------------------------------------*/
 
-#ifndef multiphaseMixture_H
-#define multiphaseMixture_H
+#ifndef multiphaseMixtureScalar_H
+#define multiphaseMixtureScalar_H
 
 #include "incompressible/transportModel/transportModel.H"
 #include "IOdictionary.H"
@@ -52,10 +52,10 @@ namespace Foam
 {
 
 /*---------------------------------------------------------------------------*\
-                      Class multiphaseMixture Declaration
+                      Class multiphaseMixtureScalar Declaration
 \*---------------------------------------------------------------------------*/
 
-class multiphaseMixture
+class multiphaseMixtureScalar
 :
     public IOdictionary,
     public transportModel
@@ -191,7 +191,7 @@ public:
     // Constructors
 
         //- Construct from components
-        multiphaseMixture
+        multiphaseMixtureScalar
         (
             const volVectorField& U,
             const surfaceScalarField& phi,
@@ -200,7 +200,7 @@ public:
 
 
     //- Destructor
-    virtual ~multiphaseMixture()
+    virtual ~multiphaseMixtureScalar()
     {}
 
 

applications/solvers/cfdemSolverMultiphaseScalar/multiphaseMixtureScalar/phase/phase.H

@@ -26,7 +26,7 @@ Class
 
 Description
     Single incompressible phase derived from the phase-fraction.
-    Used as part of the multiPhaseMixture for interface-capturing multi-phase
+    Used as part of the multiphaseMixtureScalar for interface-capturing multi-phase
     simulations.
 
 SourceFiles

doc/CFDEMcoupling_models.txt

@@ -124,6 +124,7 @@ particleDeformation,
 potentialRelaxation,
 "surfaceTensionForce"_forceModel_surfaceTensionForce.html,
 terminalVelocity,
+"transferFluidProperties"_forceModel_transferFluidProperties.html,
 turbulentDispersion,
 turbulentVelocityFluctuations,
 "virtualMassForce"_forceModel_virtualMassForce.html,

doc/forceModel_transferFluidProperties.txt

@@ -0,0 +1,42 @@
+"CFDEMproject Website"_lws - "Main Page"_main :c
+
+:link(lws,http://www.cfdem.com)
+:link(main,CFDEMcoupling_Manual.html)
+
+:line
+
+forceModel transferFluidProperties command :h3
+
+[Syntax:]
+
+Defined in "couplingProperties"_CFDEMcoupling_dicts.html#couplingProperties
+dictionary.
+
+forceModels
+(
+    transferFluidProperties
+);
+transferFluidPropertiesProps
+\{
+    verbose       switch1;
+    interpolation switch2;
+\} :pre
+
+{switch1} = (optional, default false) sub model switch, see "forceSubModel"_forceSubModel.html for details :ulb,l
+{switch2} = (optional, default false) sub model switch, see "forceSubModel"_forceSubModel.html for details :l
+:ule
+
+[Description:]
+
+This "force model" does not influence the particles or the flow - it transfer to fluid density and (dynamic)
+viscosity from OpenFOAM to LIGGGHTS.
+
+
+[Restrictions:]
+
+This model requires {fix cfd/coupling/fluidproperties} to work.
+
+[Related commands:]
+
+"forceModel"_forceModel.html
+

etc/library-list.txt

@@ -3,4 +3,4 @@ lagrangian/cfdemParticleComp/dir
 recurrence/dir
 finiteVolume/dir
 ../applications/solvers/cfdemSolverMultiphase/multiphaseMixture/dir
-../applications/solvers/cfdemSolverMultiphaseScalar/multiphaseMixture/dir
+../applications/solvers/cfdemSolverMultiphaseScalar/multiphaseMixtureScalar/dir

src/lagrangian/cfdemParticle/Make/files

@@ -96,6 +96,7 @@ $(forceModels)/potentialRelaxation/potentialRelaxation.C
 $(forceModels)/BeetstraDrag/BeetstraDrag.C
 $(forceModels)/BeetstraDragPoly/BeetstraDragPoly.C
 $(forceModels)/dSauter/dSauter.C
+$(forceModels)/transferFluidProperties/transferFluidProperties.C
 $(forceModels)/Fines/Fines.C
 $(forceModels)/Fines/FinesFields.C
 $(forceModels)/Fines/FanningDynFines.C

src/lagrangian/cfdemParticle/subModels/forceModel/transferFluidProperties/transferFluidProperties.C

@@ -0,0 +1,129 @@
+/*---------------------------------------------------------------------------*\
+License
+    This 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.
+    This code 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 this code.  If not, see <http://www.gnu.org/licenses/>.
+
+    Copyright (C) 2015- Thomas Lichtenegger, JKU Linz, Austria
+
+Description
+    transfer fluid properties to LIGGGHTS
+
+SourceFiles
+    transferFluidProperties.C
+\*---------------------------------------------------------------------------*/
+
+#include "error.H"
+
+#include "transferFluidProperties.H"
+#include "addToRunTimeSelectionTable.H"
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+namespace Foam
+{
+
+// * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
+
+defineTypeNameAndDebug(transferFluidProperties, 0);
+
+addToRunTimeSelectionTable
+(
+    forceModel,
+    transferFluidProperties,
+    dictionary
+);
+
+
+// * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
+
+// Construct from components
+transferFluidProperties::transferFluidProperties
+(
+    const dictionary& dict,
+    cfdemCloud& sm
+)
+:
+    forceModel(dict,sm),
+    propsDict_(dict.subDict(typeName + "Props"))
+{
+    particleCloud_.registerParticleProperty<double**>("fluidDensity",1);
+    particleCloud_.registerParticleProperty<double**>("fluidViscosity",1);
+
+    // init force sub model
+    setForceSubModels(propsDict_);
+    // define switches which can be read from dict
+    forceSubM(0).setSwitchesList(SW_VERBOSE,true); // activate search for verbose switch
+    forceSubM(0).setSwitchesList(SW_INTERPOLATION,true); // activate search for interpolate switch
+    forceSubM(0).readSwitches();
+}
+
+
+// * * * * * * * * * * * * * * * * Destructor  * * * * * * * * * * * * * * * //
+
+transferFluidProperties::~transferFluidProperties()
+{
+}
+
+// * * * * * * * * * * * * * * * private Member Functions  * * * * * * * * * * * * * //
+
+
+// * * * * * * * * * * * * * * * public Member Functions  * * * * * * * * * * * * * //
+
+void transferFluidProperties::setForce() const
+{
+    double**& fluidDensity_   = particleCloud_.getParticlePropertyRef<double**>("fluidDensity");
+    double**& fluidViscosity_ = particleCloud_.getParticlePropertyRef<double**>("fluidViscosity");
+
+    const volScalarField& rhoField = forceSubM(0).rhoField();
+    const volScalarField& nufField = forceSubM(0).nuField();
+
+    interpolationCellPoint<scalar> rhoInterpolator_(rhoField);
+    interpolationCellPoint<scalar> nufInterpolator_(nufField);
+
+    label cellI = 0;
+    double rho = 0.;
+    double nuf = 0.;
+    vector position(0,0,0);
+
+    for(int index = 0; index < particleCloud_.numberOfParticles(); ++index)
+    {
+        cellI = particleCloud_.cellIDs()[index][0];
+        if (cellI >= 0)
+        {
+            if(forceSubM(0).interpolation())
+            {
+                position = particleCloud_.position(index);
+                rho = rhoInterpolator_.interpolate(position,cellI);
+                nuf = nufInterpolator_.interpolate(position,cellI);
+            }
+            else
+            {
+                rho = rhoField[cellI];
+                nuf = nufField[cellI];
+            }
+
+            fluidDensity_[index][0]   = rho;
+            fluidViscosity_[index][0] = nuf*rho;
+        }
+    }
+
+    particleCloud_.dataExchangeM().giveData("fluidDensity","scalar-atom",fluidDensity_);
+    particleCloud_.dataExchangeM().giveData("fluidViscosity","scalar-atom",fluidViscosity_);
+    
+    if (forceSubM(0).verbose()) Info << "give data done" << endl;
+}
+
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+} // End namespace Foam
+
+// ************************************************************************* //

src/lagrangian/cfdemParticle/subModels/forceModel/transferFluidProperties/transferFluidProperties.H

@@ -0,0 +1,79 @@
+/*---------------------------------------------------------------------------*\
+License
+    This 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.
+    This code 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 this code.  If not, see <http://www.gnu.org/licenses/>.
+
+    Copyright (C) 2015- Thomas Lichtenegger, JKU Linz, Austria
+
+Description
+    transfer fluid properties to LIGGGHTS
+
+SourceFiles
+    transferFluidProperties.C
+\*---------------------------------------------------------------------------*/
+
+#ifndef transferFluidProperties_H
+#define transferFluidProperties_H
+
+#include "forceModel.H"
+#include "interpolationCellPoint.H"
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+namespace Foam
+{
+
+/*---------------------------------------------------------------------------*\
+                    Class transferFluidProperties Declaration
+\*---------------------------------------------------------------------------*/
+
+class transferFluidProperties
+:
+    public forceModel
+{
+private:
+
+    dictionary propsDict_;
+
+public:
+
+    //- Runtime type information
+    TypeName("transferFluidProperties");
+
+    // Constructors
+
+        //- Construct from components
+        transferFluidProperties
+        (
+            const dictionary& dict,
+            cfdemCloud& sm
+        );
+
+    // Destructor
+
+        ~transferFluidProperties();
+
+
+    // Member Functions
+        void setForce() const;
+
+};
+
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+} // End namespace Foam
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+#endif
+
+// ************************************************************************* //

src/lagrangian/cfdemParticleComp/Make/files

@@ -90,6 +90,7 @@ $(forceModels)/directedDiffusiveRelaxation/directedDiffusiveRelaxation.C
 $(forceModels)/BeetstraDrag/BeetstraDrag.C
 $(forceModels)/BeetstraDragPoly/BeetstraDragPoly.C
 $(forceModels)/dSauter/dSauter.C
+$(forceModels)/transferFluidProperties/transferFluidProperties.C
 $(forceModels)/Fines/Fines.C
 $(forceModels)/Fines/FinesFields.C
 $(forceModels)/Fines/FanningDynFines.C