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

Added check for adjustable timestep.

Browse Source
This commit is contained in:
tlichtenegger
2018-06-06 11:25:39 +02:00
parent 09baa0d6d5
commit b2f5221076
1 changed files with 10 additions and 3 deletions
Download Patch File Download Diff File Expand all files Collapse all files

13
src/lagrangian/cfdemParticle/cfdemCloud/cfdemCloud.C
View File

@ -271,6 +271,11 @@ cfdemCloud::cfdemCloud
else
Info << "ignoring ddt(voidfraction)" << endl;
bool adjustTimeStep = mesh_.time().controlDict().lookupOrDefault("adjustTimeStep", false);
if (adjustTimeStep)
FatalError << "CFDEMcoupling does not support adjustable time steps."
<< abort(FatalError);
momCoupleModel_ = new autoPtr<momCoupleModel>[momCoupleModels_.size()];
for (int i=0;i<momCoupleModels_.size();i++)
{
@ -665,10 +670,12 @@ bool cfdemCloud::evolve
//CHECK JUST TIME-INTERPOATE ALREADY SMOOTHENED VOIDFRACTIONNEXT AND UsNEXT FIELD
// IMPLICIT FORCE CONTRIBUTION AND SOLVER USE EXACTLY THE SAME AVERAGED
// QUANTITIES AT THE GRID!
Info << "\n timeStepFraction() = " << dataExchangeM().timeStepFraction() << endl;
if(dataExchangeM().timeStepFraction() > 1.0000001)
scalar timeStepFrac = dataExchangeM().timeStepFraction();
Info << "\n timeStepFraction() = " << timeStepFrac << endl;
if(timeStepFrac > 1.0000001)
{
FatalError << "cfdemCloud::dataExchangeM().timeStepFraction()>1: Do not do this, since dangerous. This might be due to the fact that you used a adjustable CFD time step. Please use a fixed CFD time step." << abort(FatalError);
// FatalError << "cfdemCloud::dataExchangeM().timeStepFraction()>1: Do not do this, since dangerous. This might be due to the fact that you used a adjustable CFD time step. Please use a fixed CFD time step." << abort(FatalError);
Warning << "cfdemCloud::dataExchangeM().timeStepFraction() = " << timeStepFrac << endl;
}
clockM().start(24,"interpolateEulerFields");