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

tweak to morph

Browse Source
This commit is contained in:
James E McClure 2018-11-02 14:17:17 -04:00
parent 9f7b530587
commit f9d5c2965a
1 changed files with 5 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
models/ColorModel.cpp
View File

@ -568,7 +568,9 @@ void ScaLBL_ColorModel::Run(){
double volB = Averages->wet_morph->V();
double volA = Averages->nonwet_morph->V();
double delta_volume = MorphInit(beta,morph_delta);
if ((1.f+delta_volume)*volA/(volA + volB) > 1.f - TARGET_SATURATION){
volA += delta_volume;
volB -= delta_volume;
if (volB/(volA + volB) > TARGET_SATURATION){
MORPH_ADAPT = false;
}
MPI_Barrier(comm);
@ -750,8 +752,8 @@ double ScaLBL_ColorModel::MorphInit(const double beta, const double morph_delta)
MPI_Allreduce(&count,&count_global,1,MPI_DOUBLE,MPI_SUM,comm);
volume_final=count_global;
double delta_volume = (volume_final-volume_initial)/volume_initial;
if (rank == 0) printf("MorphInit: change fluid volume fraction by %f \n", delta_volume);
double delta_volume = (volume_final-volume_initial);
if (rank == 0) printf("MorphInit: change fluid volume fraction by %f \n", delta_volume/volume_initial);
// 6. copy back to the device
//if (rank==0) printf("MorphInit: copy data back to device\n");