5 scalar totalParticleWeights=0;
6 vector totalForce_array(0,0,0);
7 vector totalForce_field(0,0,0);
8 scalar particleVolume_radius=0;
9 scalar particleVolume_field=0;
10 scalar particleVolume_array=0;
12 vector meanUs_array(0,0,0);
13 vector meanUs_field(0,0,0);
14 scalar meanAlpha_field=0;
15 vector meanU_field(0,0,0);
19 for(
int index = 0;index < numberOfParticles(); ++index)
22 totalForce_array[i] += impForces_[index][i];
23 meanUs_array[i] += velocities_[index][i];
25 meanR_array += radii_[index][0];
26 particleVolume_radius += radii_[index][0]*radii_[index][0]*radii_[index][0]*4./3.*3.1415;
29 for(
int subCell=0;subCell<voidFractionM().cellsPerParticle()[index][0];subCell++)
32 totalParticleWeights += particleWeights_[index][subCell];
33 particleVolume_array += particleVolumes_[index][subCell];
41 bool particleInside=
false;
42 for(
int index = 0;index < numberOfParticles(); ++index){
43 for(
int subCell=0;subCell<voidFractionM().cellsPerParticle()[index][0];subCell++){
44 if(cellIDs_[index][subCell] == cellI){
54 meanAlpha_field += alpha[cellI];
55 meanU_field += U[cellI];
56 meanUs_field += Us[cellI];
57 particleVolume_field += (1-alpha[cellI])*alpha.mesh().V()[cellI];
58 totalForce_field += forceM(0).expParticleForces()[cellI]+forceM(0).impParticleForces()[cellI];
65 meanAlpha_field /= countCell;
66 meanU_field /= countCell;
67 meanUs_field /= countCell;
72 meanU_field = vector(0,0,0);
73 meanUs_field = vector(0,0,0);
75 meanUs_array /= numberOfParticles()+SMALL;
76 meanR_array /= numberOfParticles()+SMALL;
78 Info <<
"=============================================================================" << endl;
79 Info <<
"Debug Info, only serial and not tested!" << endl;
80 Info <<
" numberOfParticles_ = "<< numberOfParticles() <<
" != " << endl;
81 Info <<
"totalParticleWeights = "<< totalParticleWeights << endl;
82 Info <<
" points= "<< points << endl;
83 Info <<
"countCell= "<< countCell << endl;
84 Info <<
" totalForce_array = "<< mag(totalForce_array)<<
" != " << endl;
85 Info <<
" totalForce_field = "<< mag(totalForce_field) << endl;
86 Info <<
" particleVolume_field = "<< particleVolume_field <<
" != " << endl;
87 Info <<
" particleVolume_array = "<< particleVolume_array <<
" != " << endl;
88 Info <<
" particleVolume_radius = "<< particleVolume_radius << endl;
89 Info <<
"meanUs_field = "<< mag(meanUs_field) <<
" ~= " << endl;
90 Info <<
"meanUs_array = "<< mag(meanUs_array) << endl;
91 Info <<
"meanU_field = "<< mag(meanU_field) << endl;
92 Info <<
"meanAlpha_field = "<< meanAlpha_field << endl;
93 Info <<
"meanR_array = "<< meanR_array << endl;
94 Info <<
"=============================================================================" << endl;