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

fix bug in permeability tolerance

Browse Source
This commit is contained in:
JamesEMcclure 2019-06-05 14:57:06 -04:00
parent f360eb1ecb
commit 8ad659ff17
1 changed files with 3 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
models/MRTModel.cpp
View File

@ -23,7 +23,7 @@ void ScaLBL_MRTModel::ReadParams(string filename){
tau = 1.0; tau = 1.0;
timestepMax = 100000; timestepMax = 100000;
tolerance = 0.01; tolerance = 1.0e-8;
Fx = Fy = 0.0; Fx = Fy = 0.0;
Fz = 1.0e-5; Fz = 1.0e-5;
@ -203,12 +203,12 @@ void ScaLBL_MRTModel::Run(){
double starttime,stoptime,cputime; double starttime,stoptime,cputime;
ScaLBL_DeviceBarrier(); MPI_Barrier(comm); ScaLBL_DeviceBarrier(); MPI_Barrier(comm);
starttime = MPI_Wtime(); starttime = MPI_Wtime();
if (rank==0) printf("Beginning AA timesteps...\n"); if (rank==0) printf("Beginning AA timesteps, timestepMax = %i \n", timestepMax);
if (rank==0) printf("********************************************************\n"); if (rank==0) printf("********************************************************\n");
timestep=0; timestep=0;
double error = 1.0; double error = 1.0;
double flow_rate_previous = 0.0; double flow_rate_previous = 0.0;
while (timestep < timestepMax && error < tolerance) { while (timestep < timestepMax && error > tolerance) {
//************************************************************************/ //************************************************************************/
timestep++; timestep++;
ScaLBL_Comm->SendD3Q19AA(fq); //READ FROM NORMAL ScaLBL_Comm->SendD3Q19AA(fq); //READ FROM NORMAL