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

increase precision of time step fraction output [ci skip]

Browse Source
This commit is contained in:
danielque
2019-07-26 16:11:09 +02:00
parent b462327d2d
commit d3d7d07cdb
1 changed files with 6 additions and 1 deletions
Download Patch File Download Diff File Expand all files Collapse all files

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

@ -45,6 +45,7 @@ Description
#include "smoothingModel.H" #include "smoothingModel.H"
#include "liggghtsCommandModel.H" #include "liggghtsCommandModel.H"
#include "otherForceModel.H" #include "otherForceModel.H"
#include "IOmanip.H"
namespace Foam namespace Foam
{ {
@ -668,12 +669,16 @@ bool cfdemCloud::evolve
//CHECK JUST TIME-INTERPOLATE ALREADY SMOOTHENED VOIDFRACTIONNEXT AND UsNEXT FIELD //CHECK JUST TIME-INTERPOLATE ALREADY SMOOTHENED VOIDFRACTIONNEXT AND UsNEXT FIELD
// IMPLICIT FORCE CONTRIBUTION AND SOLVER USE EXACTLY THE SAME AVERAGED // IMPLICIT FORCE CONTRIBUTION AND SOLVER USE EXACTLY THE SAME AVERAGED
// QUANTITIES AT THE GRID! // QUANTITIES AT THE GRID!
scalar timeStepFrac = dataExchangeM().timeStepFraction(); const scalar timeStepFrac = dataExchangeM().timeStepFraction();
int old_precision = Info().precision(10);
Info << "\n timeStepFraction() = " << timeStepFrac << endl; Info << "\n timeStepFraction() = " << timeStepFrac << endl;
Info().precision(old_precision);
if(timeStepFrac > 1.0000001) 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);
old_precision = Warning().precision(10);
Warning << "cfdemCloud::dataExchangeM().timeStepFraction() = " << timeStepFrac << endl; Warning << "cfdemCloud::dataExchangeM().timeStepFraction() = " << timeStepFrac << endl;
Warning().precision(old_precision);
} }
clockM().start(24,"interpolateEulerFields"); clockM().start(24,"interpolateEulerFields");