zhyang-dev/CFDEMcoupling-PFM
1
0
Fork 0
mirror of https://github.com/ParticulateFlow/CFDEMcoupling-PFM.git synced 2025-12-08 06:37:44 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

add interpolation

Browse Source
This commit is contained in:
tmjnijssen
2022-04-14 15:50:22 +02:00
committed by Tim MJ Nijssen
parent 4ddf691936
commit d6dab59bfd
2 changed files with 27 additions and 7 deletions
Download Patch File Download Diff File Expand all files Collapse all files

31
src/lagrangian/cfdemParticle/subModels/forceModel/transferFluidProperties/transferFluidProperties.C
View File

@ -52,14 +52,17 @@ transferFluidProperties::transferFluidProperties
) )
: :
forceModel(dict,sm), forceModel(dict,sm),
propsDict_(dict.subDict(typeName + "Props")), propsDict_(dict.subDict(typeName + "Props"))
verbose_(propsDict_.lookupOrDefault<bool>("verbose",false))
{ {
particleCloud_.registerParticleProperty<double**>("fluidDensity",1); particleCloud_.registerParticleProperty<double**>("fluidDensity",1);
particleCloud_.registerParticleProperty<double**>("fluidViscosity",1); particleCloud_.registerParticleProperty<double**>("fluidViscosity",1);
// init force sub model // init force sub model
setForceSubModels(propsDict_); 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();
} }
@ -82,22 +85,40 @@ void transferFluidProperties::setForce() const
const volScalarField& rhoField = forceSubM(0).rhoField(); const volScalarField& rhoField = forceSubM(0).rhoField();
const volScalarField& nufField = forceSubM(0).nuField(); const volScalarField& nufField = forceSubM(0).nuField();
interpolationCellPoint<scalar> rhoInterpolator_(rhoField);
interpolationCellPoint<scalar> nufInterpolator_(nufField);
label cellI = 0; label cellI = 0;
double rho = 0.;
double nuf = 0.;
vector position(0,0,0);
for(int index = 0; index < particleCloud_.numberOfParticles(); ++index) for(int index = 0; index < particleCloud_.numberOfParticles(); ++index)
{ {
cellI = particleCloud_.cellIDs()[index][0]; cellI = particleCloud_.cellIDs()[index][0];
if (cellI >= 0) if (cellI >= 0)
{ {
fluidDensity_[index][0] = rhoField[cellI]; if(forceSubM(0).interpolation())
fluidViscosity_[index][0] = nufField[cellI] * rhoField[cellI]; {
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("fluidDensity","scalar-atom",fluidDensity_);
particleCloud_.dataExchangeM().giveData("fluidViscosity","scalar-atom",fluidViscosity_); particleCloud_.dataExchangeM().giveData("fluidViscosity","scalar-atom",fluidViscosity_);
if (verbose_) Info << "give data done" << endl; if (forceSubM(0).verbose()) Info << "give data done" << endl;
} }

3
src/lagrangian/cfdemParticle/subModels/forceModel/transferFluidProperties/transferFluidProperties.H
View File

@ -24,6 +24,7 @@ SourceFiles
#define transferFluidProperties_H #define transferFluidProperties_H
#include "forceModel.H" #include "forceModel.H"
#include "interpolationCellPoint.H"
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
@ -42,8 +43,6 @@ private:
dictionary propsDict_; dictionary propsDict_;
bool verbose_;
public: public:
//- Runtime type information //- Runtime type information