added flow reversal criterion at endpoint

This commit is contained in:
James McClure 2019-03-06 06:33:01 -05:00
parent 33ac071b1b
commit 1588356f0d

View File

@ -606,6 +606,10 @@ void ScaLBL_ColorModel::Run(){
fprintf(kr_log_file,"%i %.5g %.5g %.5g %.5g %.5g %.5g ",timestep-analysis_interval+20,muA,muB,5.796*alpha,Fx,Fy,Fz);
fprintf(kr_log_file,"%.5g %.5g %.5g %.5g %.5g %.5g %.5g %.5g\n",volA,volB,vA_x,vA_y,vA_z,vB_x,vB_y,vB_z);
fclose(kr_log_file);
// Add exit criteria based on relative permeability ratio
if (volA*sqrt(vA_x*vA_x + vA_y*vA_y + vA_z*vA_z)/(volB*sqrt(vB_x*vB_x + vB_y*vB_y + vB_z*vB_z)) < 0.05){
timestep = timestepMax;
}
printf(" Measured capillary number %f \n ",Ca);
}