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

flow reversal tweaks

Browse Source
This commit is contained in:
James E McClure
2019-04-09 07:45:23 -04:00
parent 84e61f7801
commit dd321d3cc6
1 changed files with 2 additions and 13 deletions
Download Patch File Download Diff File Expand all files Collapse all files

15
models/ColorModel.cpp
View File

@@ -642,19 +642,6 @@ void ScaLBL_ColorModel::Run(){
if (rank == 0) printf(" -- adjust force by factor %f \n ",capillary_number / Ca);
Averages->SetParams(rhoA,rhoB,tauA,tauB,Fx,Fy,Fz,alpha,beta);
}
// flow reversal criteria based on fractional flow rate
if (delta_volume_target < 0.0 &&
volA*flow_rate_A/(volA*flow_rate_A+volB*flow_rate_B) < RESIDUAL_ENDPOINT_THRESHOLD){
REVERSE_FLOW_DIRECTION = true;
}
else if (delta_volume_target > 0.0 &&
volB*flow_rate_B/(volA*flow_rate_A+volB*flow_rate_B) < RESIDUAL_ENDPOINT_THRESHOLD){
REVERSE_FLOW_DIRECTION = true;
}
if ( REVERSE_FLOW_DIRECTION ){
delta_volume_target *= (-1.0);
REVERSE_FLOW_DIRECTION = false;
}
}
else{
if (rank==0){
@@ -686,6 +673,8 @@ void ScaLBL_ColorModel::Run(){
// flow direction will reverse after next steady point
MORPH_ADAPT = false;
CURRENT_STEADY_TIMESTEPS=0;
delta_volume_target *= (-1.0);
REVERSE_FLOW_DIRECTION = false;
}
MPI_Barrier(comm);