OilfieldToolsNet/LBPM
4
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

fixed little bug

Browse Source
This commit is contained in:
James E McClure 2015-06-16 17:39:05 -04:00
parent 2ca2203a3c
commit e6e45f6243
1 changed files with 2 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

5
tests/lbpm_permeability_simulator.cpp
View File

@ -1247,9 +1247,8 @@ int main(int argc, char **argv)
double mag_force = sqrt(Fx*Fx+Fy*Fy+Fz*Fz);
double Fo = reference_length*D32*D32*mag_force/viscosity/viscosity;
// .... 1-D flow should be aligned with force ...
double velocity = vawx*Fx/mag_force + vawy*Fy/mag_force
+ vawz*Fz/mag_force;
double err1D = fabs((velocity-sqrt(vawx*vawx+vawy*vawy+vawz*vawz))/velocity;
double velocity = vawx*Fx/mag_force + vawy*Fy/mag_force + vawz*Fz/mag_force;
double err1D = fabs(velocity-sqrt(vawx*vawx+vawy*vawy+vawz*vawz))/velocity;
//.......... Computation of the Reynolds number Re ..............
double Re = D32*velocity/viscosity;
printf("Relative error for 1D representation: %f \n",err1D);