2 if(particleCloud.verbose())
4 Info <<
"=============================================================================" << endl;
5 Info <<
"totaldragforceEuler calculus" << endl;
6 vector totaldragforceEuler(0,0,0);
9 totaldragforceEuler += Ksl[cellI]*(Us[cellI]-U[cellI])/rho[cellI] * Ksl.mesh().V()[cellI];
11 Pout <<
"totaldragforceEuler = "<< mag(totaldragforceEuler) << endl;
12 Pout <<
"dv/dt =" << sum(fvc::ddt(voidfraction)) << endl;
13 Info <<
"=============================================================================" << endl;