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

Merge pull request #63 from thomaram/adjust_perm

Browse Source 
Add adjusted perms
This commit is contained in:
JamesEMcClure
2022-03-21 14:16:13 -04:00
committed by GitHub
parent fe2496ebdd 766dfc299a
commit 70b12830cb
2 changed files with 18 additions and 10 deletions
Download Patch File Download Diff File Expand all files Collapse all files

19
models/ColorModel.cpp
View File

@@ -758,7 +758,7 @@ double ScaLBL_ColorModel::Run(int returntime) {
double volA = Averages->gnb.V;
volA /= Dm->Volume;
volB /= Dm->Volume;
;
//initial_volume = volA*Dm->Volume;
double vA_x = Averages->gnb.Px / Averages->gnb.M;
double vA_y = Averages->gnb.Py / Averages->gnb.M;
@@ -1015,14 +1015,19 @@ double ScaLBL_ColorModel::Run(int returntime) {
"Ca "
"eff.pressure\n");
}
// Adjust the effperms with porosity term to make the nominal
// values closer to measured data
fprintf(scal_log_file, "%i %.5g ", CURRENT_TIMESTEP,
current_saturation);
fprintf(scal_log_file, "%.5g %.5g ", kAeff_low, kBeff_low);
fprintf(scal_log_file, "%.5g %.5g ", kAeff_connected_low,
kBeff_connected_low);
fprintf(scal_log_file, "%.5g %.5g ", kAeff_disconnected,
kBeff_disconnected);
fprintf(scal_log_file, "%.5g %.5g ", krnf_low, krwf_low);
fprintf(scal_log_file, "%.5g %.5g ", kAeff_low * Mask->Porosity(),
kBeff_low * Mask->Porosity());
fprintf(scal_log_file, "%.5g %.5g ", kAeff_connected_low * Mask->Porosity(),
kBeff_connected_low * Mask->Porosity());
fprintf(scal_log_file, "%.5g %.5g ", kAeff_disconnected * Mask->Porosity(),
kBeff_disconnected * Mask->Porosity());
fprintf(scal_log_file, "%.5g %.5g ", krnf_low * Mask->Porosity(),
krwf_low * Mask->Porosity());
fprintf(scal_log_file, "%.5g %.5g %.5g ", pAB,
pAB_connected, Ca);
fprintf(scal_log_file, "%.5g\n", eff_pres);

9
models/MRTModel.cpp
View File

@@ -254,7 +254,7 @@ void ScaLBL_MRTModel::Run() {
if (WriteHeader) {
log_file = fopen("Permeability.csv", "a+");
fprintf(log_file, "time Fx Fy Fz mu Vs As Js Xs vx vy vz k\n");
fprintf(log_file, "time Fx Fy Fz mu Vs As Js Xs vx vy vz absperm absperm_adjusted\n");
fclose(log_file);
}
}
@@ -384,14 +384,17 @@ void ScaLBL_MRTModel::Run() {
double h = Dm->voxel_length;
double absperm =
h * h * mu * Mask->Porosity() * flow_rate / force_mag;
double absperm_adjusted =
h * h * mu * Mask->Porosity() * Mask->Porosity() * flow_rate / force_mag;
if (rank == 0) {
printf(" %f\n", absperm);
FILE *log_file = fopen("Permeability.csv", "a");
fprintf(log_file,
"%i %.8g %.8g %.8g %.8g %.8g %.8g %.8g %.8g %.8g %.8g "
"%.8g %.8g\n",
"%.8g %.8g %.8g\n",
timestep, Fx, Fy, Fz, mu, h * h * h * Vs, h * h * As,
h * Hs, Xs, vax, vay, vaz, absperm);
h * Hs, Xs, vax, vay, vaz, absperm, absperm_adjusted);
fclose(log_file);
}
}