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

reverse flow direction if connected pathway is removed

Browse Source
This commit is contained in:
James E McClure
2019-03-25 10:50:58 -04:00
parent 24f5fff3f3
commit d0a516ecf6
2 changed files with 17 additions and 5 deletions
Download Patch File Download Diff File Expand all files Collapse all files

21
models/ColorModel.cpp
View File

@@ -10,7 +10,7 @@ rank(RANK), nprocs(NP), Restart(0),timestep(0),timestepMax(0),tauA(0),tauB(0),rh
Fx(0),Fy(0),Fz(0),flux(0),din(0),dout(0),inletA(0),inletB(0),outletA(0),outletB(0), Fx(0),Fy(0),Fz(0),flux(0),din(0),dout(0),inletA(0),inletB(0),outletA(0),outletB(0),
Nx(0),Ny(0),Nz(0),N(0),Np(0),nprocx(0),nprocy(0),nprocz(0),BoundaryCondition(0),Lx(0),Ly(0),Lz(0),comm(COMM) Nx(0),Ny(0),Nz(0),N(0),Np(0),nprocx(0),nprocy(0),nprocz(0),BoundaryCondition(0),Lx(0),Ly(0),Lz(0),comm(COMM)
{ {
REVERSE_FLOW_DIRECTION = false;
} }
ScaLBL_ColorModel::~ScaLBL_ColorModel(){ ScaLBL_ColorModel::~ScaLBL_ColorModel(){
@@ -621,11 +621,15 @@ void ScaLBL_ColorModel::Run(){
// flow reversal criteria based on fractional flow rate // flow reversal criteria based on fractional flow rate
if (delta_volume_target < 0.0 && if (delta_volume_target < 0.0 &&
volA*flow_rate_A/(volA*flow_rate_A+volB*flow_rate_B) < RESIDUAL_ENDPOINT_THRESHOLD){ volA*flow_rate_A/(volA*flow_rate_A+volB*flow_rate_B) < RESIDUAL_ENDPOINT_THRESHOLD){
delta_volume_target *= (-1.0); REVERSE_FLOW_DIRECTION = true;
} }
else if (delta_volume_target > 0.0 && else if (delta_volume_target > 0.0 &&
volB*flow_rate_B/(volA*flow_rate_A+volB*flow_rate_B) < RESIDUAL_ENDPOINT_THRESHOLD){ 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); delta_volume_target *= (-1.0);
REVERSE_FLOW_DIRECTION = false;
} }
} }
else{ else{
@@ -651,6 +655,12 @@ void ScaLBL_ColorModel::Run(){
MORPH_ADAPT = false; MORPH_ADAPT = false;
CURRENT_STEADY_TIMESTEPS=0; CURRENT_STEADY_TIMESTEPS=0;
} }
if ( REVERSE_FLOW_DIRECTION ){
// flow direction will reverse after next steady point
MORPH_ADAPT = false;
CURRENT_STEADY_TIMESTEPS=0;
}
MPI_Barrier(comm); MPI_Barrier(comm);
} }
morph_timesteps += analysis_interval; morph_timesteps += analysis_interval;
@@ -760,11 +770,12 @@ double ScaLBL_ColorModel::MorphInit(const double beta, const double target_delta
} }
/* if (volume_connected < 0.025*volume_initial){ if (volume_connected < 0.025*volume_initial){
// if connected volume is less than 2.5% just delete the whole thing // if connected volume is less than 2.5% just delete the whole thing
if (rank==0) printf("Connected region has shrunk to less than 2.5% of total fluid volume (remove the whole thing) \n"); if (rank==0) printf("Connected region has shrunk to less than 2.5% of total fluid volume (remove the whole thing) \n");
REVERSE_FLOW_DIRECTION = true;
} }
else { */ else {
if (rank==0) printf("MorphGrow with target volume fraction change %f \n", target_delta_volume/volume_initial); if (rank==0) printf("MorphGrow with target volume fraction change %f \n", target_delta_volume/volume_initial);
double target_delta_volume_incremental = target_delta_volume; double target_delta_volume_incremental = target_delta_volume;
if (fabs(target_delta_volume) > 0.01*volume_initial) if (fabs(target_delta_volume) > 0.01*volume_initial)
@@ -799,7 +810,7 @@ double ScaLBL_ColorModel::MorphInit(const double beta, const double target_delta
} }
} }
fillDouble.fill(phase); fillDouble.fill(phase);
//} }
count = 0.f; count = 0.f;
for (int k=1; k<Nz-1; k++){ for (int k=1; k<Nz-1; k++){

1
models/ColorModel.h
View File

@@ -32,6 +32,7 @@ public:
void WriteDebug(); void WriteDebug();
bool Restart,pBC; bool Restart,pBC;
bool REVERSE_FLOW_DIRECTION;
int timestep,timestepMax; int timestep,timestepMax;
int BoundaryCondition; int BoundaryCondition;
double tauA,tauB,rhoA,rhoB,alpha,beta; double tauA,tauB,rhoA,rhoB,alpha,beta;