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

fix volume bug

Browse Source
This commit is contained in:
James McClure 2020-04-04 09:00:46 -04:00
parent bad52221a8
commit b81d4199a1
1 changed files with 3 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
models/ColorModel.cpp
View File

@ -760,7 +760,7 @@ void ScaLBL_ColorModel::Run(){
double volA = Averages->gnb.V; double volA = Averages->gnb.V;
volA /= Dm->Volume; volA /= Dm->Volume;
volB /= Dm->Volume;; volB /= Dm->Volume;;
initial_volume = volA*Dm->Volume; //initial_volume = volA*Dm->Volume;
double vA_x = Averages->gnb.Px/Averages->gnb.M; double vA_x = Averages->gnb.Px/Averages->gnb.M;
double vA_y = Averages->gnb.Py/Averages->gnb.M; double vA_y = Averages->gnb.Py/Averages->gnb.M;
double vA_z = Averages->gnb.Pz/Averages->gnb.M; double vA_z = Averages->gnb.Pz/Averages->gnb.M;
@ -1221,6 +1221,7 @@ double ScaLBL_ColorModel::SeedPhaseField(const double seed_water_in_oil){
Bq_tmp[n+4*Np] += 0.1111111111111111*random_value; Bq_tmp[n+4*Np] += 0.1111111111111111*random_value;
Bq_tmp[n+5*Np] += 0.1111111111111111*random_value; Bq_tmp[n+5*Np] += 0.1111111111111111*random_value;
Bq_tmp[n+6*Np] += 0.1111111111111111*random_value; Bq_tmp[n+6*Np] += 0.1111111111111111*random_value;
count += 1.0;
} }
mass_loss += random_value*seed_water_in_oil; mass_loss += random_value*seed_water_in_oil;
} }
@ -1248,6 +1249,7 @@ double ScaLBL_ColorModel::SeedPhaseField(const double seed_water_in_oil){
Bq_tmp[n+4*Np] += 0.1111111111111111*random_value; Bq_tmp[n+4*Np] += 0.1111111111111111*random_value;
Bq_tmp[n+5*Np] += 0.1111111111111111*random_value; Bq_tmp[n+5*Np] += 0.1111111111111111*random_value;
Bq_tmp[n+6*Np] += 0.1111111111111111*random_value; Bq_tmp[n+6*Np] += 0.1111111111111111*random_value;
count += 1.0;
} }
mass_loss += random_value*seed_water_in_oil; mass_loss += random_value*seed_water_in_oil;
} }