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

added forcing floor for morphLBM

Browse Source
This commit is contained in:
James E McClure 2018-10-27 18:12:47 -04:00
parent 23dd734a40
commit 757805628c
1 changed files with 1 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

1
models/ColorModel.cpp
View File

@ -552,6 +552,7 @@ void ScaLBL_ColorModel::Run(){
if (lead_timesteps > 5000){ if (lead_timesteps > 5000){
Fz *= capillary_number / Ca; Fz *= capillary_number / Ca;
if (Fz > 1e-3) Fz = 1e-3; // impose ceiling for stability if (Fz > 1e-3) Fz = 1e-3; // impose ceiling for stability
if (Fz < 1e-6) Fz = 1e-6; // impose floor so we don't do something super dumb
tolerance = fabs(capillary_number - Ca) / capillary_number ; tolerance = fabs(capillary_number - Ca) / capillary_number ;
if (rank == 0) printf(" -- adjust force by %f \n ",tolerance); if (rank == 0) printf(" -- adjust force by %f \n ",tolerance);
} }