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

communicate pressure with DEM

Browse Source
This commit is contained in:
ekinaci
2018-02-08 17:42:37 +01:00
parent 3b8db2aa53
commit 34897bad97
2 changed files with 17 additions and 1 deletions
Download Patch File Download Diff File Expand all files Collapse all files

12
src/lagrangian/cfdemParticle/subModels/chemistryModel/diffusionCoefficients/diffusionCoefficients.C
View File

@ -71,6 +71,8 @@ diffusionCoefficient::diffusionCoefficient
tempField_(sm.mesh().lookupObject<volScalarField> (tempFieldName_)),
pressureFieldName_(propsDict_.lookupOrDefault<word>("pressureFieldName","p")),
P_(sm.mesh().lookupObject<volScalarField>(pressureFieldName_)),
partPressureName_(propsDict_.lookupOrDefault<word>("partPressureName","partP")),
partPressure_(NULL),
densityFieldName_(propsDict_.lookupOrDefault<word>("densityFieldName","rho")),
rho_(sm.mesh().lookupObject<volScalarField> (densityFieldName_)),
molarConcFieldName_(propsDict_.lookupOrDefault<word>("totalMoleFieldName","molarConc")),
@ -101,6 +103,8 @@ diffusionCoefficient::~diffusionCoefficient()
coeffs.clearStorage();
molWeight.clearStorage();
particleCloud_.dataExchangeM().destroy(partPressure_,1);
int nP_ = particleCloud_.numberOfParticles();
for (int i=0; i<diffusantGasNames_.size(); i++)
{
@ -121,6 +125,8 @@ diffusionCoefficient::~diffusionCoefficient()
double initVal=0.0;
if (particleCloud_.dataExchangeM().maxNumberOfParticles() > 0)
{
particleCloud_.dataExchangeM().allocateArray(partPressure_,initVal,1,"nparticles");
for (int i=0; i<diffusantGasNames_.size(); i++)
{
particleCloud_.dataExchangeM().allocateArray(diffusionCoefficients_[i],initVal,1,"nparticles");
@ -133,6 +139,7 @@ void diffusionCoefficient::reAllocMyArrays() const
if (particleCloud_.numberOfParticlesChanged())
{
double initVal=0.0;
particleCloud_.dataExchangeM().allocateArray(partPressure_,initVal,1,"nparticles");
for (int i=0; i<diffusantGasNames_.size(); i++)
{
@ -239,6 +246,9 @@ void diffusionCoefficient::execute()
Pfluid = Pfluid/101325;
Info << "pressure field" << Pfluid << nl << endl;
partPressure_[index][0] = Pfluid;
Info << "partPressure_[index][0] = " << partPressure_[index][0] << endl;
// Texp = pow(Tfluid,1.75);
Texp = Tfluid*sqrt(sqrt(Tfluid*Tfluid*Tfluid));
Info << "T - exponent calculation" << Texp << nl << endl;
@ -324,6 +334,8 @@ void diffusionCoefficient::execute()
}
}
particleCloud_.dataExchangeM().giveData(partPressureName_, "scalar-atom",partPressure_);
for (int i=0; i<diffusantGasNames_.size();i++)
{
word pushName = diffusantGasNames_[i] + "_diffCoeff";

6
src/lagrangian/cfdemParticle/subModels/chemistryModel/diffusionCoefficients/diffusionCoefficients.H
View File

@ -70,7 +70,11 @@ private:
// gas pressure at particle location
word pressureFieldName_;
const volScalarField& P_;
const volScalarField& P_;
word partPressureName_;
mutable double **partPressure_;
word densityFieldName_;