working on BC switch

This commit is contained in:
James E McClure
2019-05-16 14:44:31 -04:00
parent 46f04e6ab6
commit f92cd802fd

View File

@@ -741,6 +741,9 @@ void ScaLBL_ColorModel::Run(){
BoundaryCondition = 4;
flux = morph_delta*5.796*alpha*capillary_number/(fabs(morph_delta)*muA*Dm->Lz);
Fx = Fy = Fz = 0.0;
ScaLBL_Comm->BoundaryCondition = 4;
ScaLBL_Comm_Regular->BoundaryCondition = 4;
if (rank==0) printf("***Direct simulation with flux(A) %f, volume change %f / %f ***\n",flux, delta_volume, delta_volume_target);
}
else if (USE_SEED){
@@ -770,6 +773,8 @@ void ScaLBL_ColorModel::Run(){
Fy = capillary_number*dir_y*force_mag / Ca;
Fz = capillary_number*dir_z*force_mag / Ca;
}
ScaLBL_Comm->BoundaryCondition = 0;
ScaLBL_Comm_Regular->BoundaryCondition = 0;
}
else if (!(USE_DIRECT) && CURRENT_MORPH_TIMESTEPS > MAX_MORPH_TIMESTEPS) {
MORPH_ADAPT = false;