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

built passed

Browse Source
This commit is contained in:
Rex Zhe Li
2021-02-08 00:00:26 -05:00
parent 98491ccd80
commit e34170d232
3 changed files with 9 additions and 10 deletions
Download Patch File Download Diff File Expand all files Collapse all files

9
cpu/FreeLee.cpp
View File

@@ -1738,12 +1738,11 @@ extern "C" void ScaLBL_D3Q19_AAeven_FreeLeeModel_SingleFluid_BGK(double *dist, d
m18 = dist[17*Np+n];
//compute fluid velocity
ux = 3.0/rho0*(m1-m2+m7-m8+m9-m10+m11-m12+m13-m14+0.5*(chem*nx+Fx));
uy = 3.0/rho0*(m3-m4+m7-m8-m9+m10+m15-m16+m17-m18+0.5*(chem*ny+Fy));
uz = 3.0/rho0*(m5-m6+m11-m12-m13+m14+m15-m16-m17+m18+0.5*(chem*nz+Fz));
ux = 3.0/rho0*(m1-m2+m7-m8+m9-m10+m11-m12+m13-m14+0.5*(Fx));
uy = 3.0/rho0*(m3-m4+m7-m8-m9+m10+m15-m16+m17-m18+0.5*(Fy));
uz = 3.0/rho0*(m5-m6+m11-m12-m13+m14+m15-m16-m17+m18+0.5*(Fz));
//compute pressure
p = (m0+m2+m1+m4+m3+m6+m5+m8+m7+m10+m9+m12+m11+m14+m13+m16+m15+m18+m17)
+0.5*(rhoA-rhoB)/2.0/3.0*(ux*nx+uy*ny+uz*nz);
p = (m0+m2+m1+m4+m3+m6+m5+m8+m7+m10+m9+m12+m11+m14+m13+m16+m15+m18+m17);
//------------------------------------------------- BCK collison ------------------------------------------------------------//
// q=0

2
models/FreeLeeModel.cpp
View File

@@ -1006,7 +1006,7 @@ void ScaLBL_FreeLeeModel::WriteDebug_TwoFluid(){
void ScaLBL_FreeLeeModel::WriteDebug_SingleFluid(){
DoubleArray PhaseData(Nxh,Nyh,Nzh);
DoubleArray PhaseField(Nx,Ny,Nz);
// Copy back final phase indicator field and convert to regular layout
ScaLBL_Comm->RegularLayout(Map,Pressure,PhaseField);

8
tests/lbpm_freelee_SingleFluidBGK_simulator.cpp
View File

@@ -51,10 +51,10 @@ int main( int argc, char **argv )
LeeModel.ReadParams( filename );
LeeModel.SetDomain();
LeeModel.ReadInput();
LeeModel.Create_SingleFluid()();
LeeModel.Initialize_SingleFluid()();
LeeModel.Run_SingleFluid()();
LeeModel.WriteDebug_SingleFluid()();
LeeModel.Create_SingleFluid();
LeeModel.Initialize_SingleFluid();
LeeModel.Run_SingleFluid();
LeeModel.WriteDebug_SingleFluid();
PROFILE_STOP("Main");
auto file = db->getWithDefault<std::string>( "TimerFile", "lbpm_freelee_SingleFluidBGK_simulator" );