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

disable steady point rejection

Browse Source
This commit is contained in:
James E McClure 2020-05-18 11:12:27 -04:00
parent 02643365f9
commit ea614489df
1 changed files with 5 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

7
models/ColorModel.cpp
View File

@ -815,7 +815,7 @@ void ScaLBL_ColorModel::Run(){
color_db->putVector<double>("F",{Fx,Fy,Fz});
}
if ( isSteady ){
if (SET_CAPILLARY_NUMBER && fabs(capillary_number - Ca) / capillary_number > 2.0){
/*if (SET_CAPILLARY_NUMBER && fabs(capillary_number - Ca) / capillary_number > 2.0){
// reject steady points if they don't match the Ca well enough
double RESCALE_FORCE_FACTOR = capillary_number / Ca;
if (RESCALE_FORCE_FACTOR > 2.0) RESCALE_FORCE_FACTOR = 2.0;
@ -836,6 +836,7 @@ void ScaLBL_ColorModel::Run(){
CURRENT_STEADY_TIMESTEPS=0;
}
else {
*/
MORPH_ADAPT = true;
CURRENT_MORPH_TIMESTEPS=0;
delta_volume_target = Dm->Volume*volA *morph_delta; // set target volume change
@ -956,7 +957,9 @@ void ScaLBL_ColorModel::Run(){
}
morph_timesteps=0;
Ca_previous = Ca;
}
/*}
THIS IS WHERE YOU DISABLE STEADY POINT REJECTION
*/
}
if (MORPH_ADAPT ){