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

Change of averaging procedure.

Browse Source
This commit is contained in:
Thomas Lichtenegger
2016-05-17 08:21:27 +02:00
parent d4ccd008ba
commit a731d14e72
1 changed files with 7 additions and 5 deletions
Download Patch File Download Diff File Expand all files Collapse all files

12
src/lagrangian/cfdemParticle/subModels/energyModel/heatTransferGunn/heatTransferGunn.C
View File

@ -110,6 +110,8 @@ heatTransferGunn::heatTransferGunn
partRefTemp_.value()=readScalar(propsDict_.lookup ("partRefTemp"));
partTempField_.writeOpt() = IOobject::AUTO_WRITE;
partRelTempField_.writeOpt() = IOobject::AUTO_WRITE;
partTempField_.write();
partRelTempField_.write();
Info << "Particle temperature field activated." << endl;
}
}
@ -146,11 +148,7 @@ void heatTransferGunn::calcEnergyContribution()
particleCloud_.dataExchangeM().getData(partTempName_,"scalar-atom",partTemp_);
if(calcPartTempField_)
{
dimensionedScalar aveTemp("aveTemp", dimensionSet(0,0,0,1,0,0,0), 0.0);
for(int index = 0;index < particleCloud_.numberOfParticles(); ++index)
aveTemp.value() += partTemp_[index][0];
aveTemp.value() /= particleCloud_.numberOfParticles();
{
partTempField_.internalField() = 0.0;
particleCloud_.averagingM().resetWeightFields();
particleCloud_.averagingM().setScalarAverage
@ -161,6 +159,10 @@ void heatTransferGunn::calcEnergyContribution()
particleCloud_.averagingM().UsWeightField(),
NULL
);
volScalarField volP (1 - voidfraction_);
volScalarField weigthedTp (volP * partTempField_);
// average per cell-value, not per volume * cell-value
dimensionedScalar aveTemp = weigthedTp.average() / volP.average();
partRelTempField_ = (partTempField_ - aveTemp) / (aveTemp - partRefTemp_);
}