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

fix a few bugs

Browse Source
This commit is contained in:
James McClure 2021-06-11 17:29:07 -04:00
parent 2dd652d495
commit 9d968ee2ca
1 changed files with 5 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

10
tests/lbpm_color_simulator.cpp
View File

@ -74,7 +74,7 @@ int main( int argc, char **argv )
/* flow adaptor keys to control */ /* flow adaptor keys to control */
int SKIP_TIMESTEPS = 0; int SKIP_TIMESTEPS = 0;
int MAX_STEADY_TIME = 1000000; int MAX_STEADY_TIME = 1000000;
double FRACTIONAL_FLOW_INCREMENT = 0.05 double FRACTIONAL_FLOW_INCREMENT = 0.05;
if (ColorModel.db->keyExists( "FlowAdaptor" )){ if (ColorModel.db->keyExists( "FlowAdaptor" )){
auto flow_db = ColorModel.db->getDatabase( "FlowAdaptor" ); auto flow_db = ColorModel.db->getDatabase( "FlowAdaptor" );
MAX_STEADY_TIME = flow_db->getWithDefault<int>( "max_steady_timesteps", 1000000 ); MAX_STEADY_TIME = flow_db->getWithDefault<int>( "max_steady_timesteps", 1000000 );
@ -98,15 +98,15 @@ int main( int argc, char **argv )
int skip_time = 0; int skip_time = 0;
timestep = ColorModel.timestep; timestep = ColorModel.timestep;
double SaturationChange = 0.0; double SaturationChange = 0.0;
double volB = M.Averages->gwb.V; double volB = ColorModel.Averages->gwb.V;
double volA = M.Averages->gnb.V; double volA = ColorModel.Averages->gnb.V;
double initialSaturation = volB/(volA + volB); double initialSaturation = volB/(volA + volB);
while (skip_time < SKIP_TIMESTEPS && fabs(SaturationChange) < fabs(FRACTIONAL_FLOW_INCREMENT) ){ while (skip_time < SKIP_TIMESTEPS && fabs(SaturationChange) < fabs(FRACTIONAL_FLOW_INCREMENT) ){
timestep += ANALYSIS_INTERVAL; timestep += ANALYSIS_INTERVAL;
Adapt.UpdateFractionalFlow(ColorModel); Adapt.UpdateFractionalFlow(ColorModel);
MLUPS = ColorModel.Run(timestep); MLUPS = ColorModel.Run(timestep);
double volB = M.Averages->gwb.V; double volB = ColorModel.Averages->gwb.V;
double volA = M.Averages->gnb.V; double volA = ColorModel.Averages->gnb.V;
SaturationChange = volB/(volA + volB) - initialSaturation; SaturationChange = volB/(volA + volB) - initialSaturation;
skip_time += ANALYSIS_INTERVAL; skip_time += ANALYSIS_INTERVAL;
} }