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

tweaking morph adapt

Browse Source
This commit is contained in:
James E McClure
2018-11-04 15:03:19 -05:00
parent 11f9fe00cd
commit 37702347cc
1 changed files with 6 additions and 5 deletions
Download Patch File Download Diff File Expand all files Collapse all files

11
models/ColorModel.cpp
View File

@@ -612,16 +612,17 @@ void ScaLBL_ColorModel::Run(){
Fy *= 1e-3/force_magnitude;
Fz *= 1e-3/force_magnitude;
}
if (force_magnitude < 1e-7){
Fx *= 1e-7/force_magnitude; // impose floor
Fy *= 1e-7/force_magnitude;
Fz *= 1e-7/force_magnitude;
if (force_magnitude < 1e-6){
Fx *= 1e-6/force_magnitude; // impose floor
Fy *= 1e-6/force_magnitude;
Fz *= 1e-6/force_magnitude;
}
if (Fx*vA_x + Fy*vA_y + Fz*vA_z < 0.f || Fx*vB_x + Fy*vB_y +Fz*vB_z < 0.f){
/* if (Fx*vA_x + Fy*vA_y + Fz*vA_z < 0.f || Fx*vB_x + Fy*vB_y +Fz*vB_z < 0.f){
Fx *= 2.f; // bigger forces are needed if flow rate is "too small"
Fy *= 2.f;
Fz *= 2.f;
}
*/
tolerance = fabs(1.f - (0.95 + 0.05*capillary_number / Ca));
if (rank == 0) printf(" -- adjust force by %f \n ",tolerance);
}