CFDEMcoupling  2.4
 All Classes
solverDebugInfo.H
1 
2 if(particleCloud.verbose())
3 {
4  Info <<"=============================================================================" << endl;
5  Info << "totaldragforceEuler calculus" << endl;
6  vector totaldragforceEuler(0,0,0);
7  forAll(Ksl,cellI)
8  {
9  totaldragforceEuler += Ksl[cellI]*(Us[cellI]-U[cellI])/rho[cellI] * Ksl.mesh().V()[cellI];
10  }
11  Pout <<"totaldragforceEuler = "<< mag(totaldragforceEuler) << endl;
12  Pout << "dv/dt =" << sum(fvc::ddt(voidfraction)) << endl;
13  Info <<"=============================================================================" << endl;
14 }