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

fix laplacian

Browse Source
This commit is contained in:
James McClure 2021-01-19 11:38:23 -05:00
parent d442f890f6
commit e21774d021
1 changed files with 2 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
analysis/SubPhase.cpp
View File

@ -157,7 +157,7 @@ void SubPhase::SetParams(double rhoA, double rhoB, double tauA, double tauB, dou
} }
void SubPhase::Basic(){ void SubPhase::Basic(){
int i,j,k,n,imin,jmin,kmin,kmax, nq; int i,j,k,n,imin,jmin,kmin,kmax;
// If external boundary conditions are set, do not average over the inlet // If external boundary conditions are set, do not average over the inlet
kmin=1; kmax=Nz-1; kmin=1; kmax=Nz-1;
@ -182,7 +182,7 @@ void SubPhase::Basic(){
double fx = (Phi(i+1,j,k) - 2.0*Phi(i,j,k) + Phi(i-1,j,k)); double fx = (Phi(i+1,j,k) - 2.0*Phi(i,j,k) + Phi(i-1,j,k));
double fy = (Phi(i,j+1,k) - 2.0*Phi(i,j,k) + Phi(i,j-1,k)); double fy = (Phi(i,j+1,k) - 2.0*Phi(i,j,k) + Phi(i,j-1,k));
double fz = (Phi(i,j,k+1) - 2.0*Phi(i,j,k) + Phi(i,j,k-1)); double fz = (Phi(i,j,k+1) - 2.0*Phi(i,j,k) + Phi(i,j,k-1));
Laplacian(i,j,k) = sqrt(fx*fx+fy*fy+fz*fz); Laplacian(i,j,k) = 0.5*(fx + fy + fz);
} }
} }
} }