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

Added wetting phase properties to blobstats.tcat from tests/BlobAnalysis.cpp

Browse Source
This commit is contained in:
James E McClure
2015-05-18 12:12:36 -04:00
parent 49753252c8
commit a6b8fe00a8
1 changed files with 2 additions and 2 deletions
Download Patch File Download Diff File Expand all files Collapse all files

4
tests/BlobAnalysis.cpp
View File

@@ -954,7 +954,7 @@ int main(int argc, char **argv)
printf("Sauter mean diamter = %f \n",D); printf("Sauter mean diamter = %f \n",D);
printf("WARNING: lazy coder hard-coded the surface tension as 0.058 \n"); printf("WARNING: lazy coder hard-coded the surface tension as 0.058 \n");
FILE *BLOBSTATES= fopen("blobstates.tcat","a"); FILE *BLOBSTATES= fopen("blobstates.tcat","a");
fprintf(BLOBSTATES,"%.5g %.5g %.5g\n",vol_w,pw,aws);
// Compute the averages over the entire non-wetting phsae // Compute the averages over the entire non-wetting phsae
for (a=0; a<nblobs; a++){ for (a=0; a<nblobs; a++){
nwp_volume += BlobAverages(0,a); nwp_volume += BlobAverages(0,a);
@@ -1030,7 +1030,7 @@ int main(int argc, char **argv)
Gws(4) += BlobAverages(24,a)*BlobAverages(3,a); Gws(4) += BlobAverages(24,a)*BlobAverages(3,a);
Gws(5) += BlobAverages(25,a)*BlobAverages(3,a); Gws(5) += BlobAverages(25,a)*BlobAverages(3,a);
if (fabs(1.0 - nwp_volume*iVol/porosity - sw) > 0.0025 || a == 0){ if (fabs(1.0 - nwp_volume*iVol/porosity - sw) > 0.0025 || a == 1){
sw = 1.0 - nwp_volume*iVol/porosity; sw = 1.0 - nwp_volume*iVol/porosity;
JwnD = -Jwn*D/awn; JwnD = -Jwn*D/awn;