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

Bug fixed in normal velocity of wn interface

Browse Source
This commit is contained in:
James McClure
2013-12-31 17:55:18 -05:00
parent 68f0610716
commit 32cacd8afd
1 changed files with 12 additions and 4 deletions
Download Patch File Download Diff File Expand all files Collapse all files

16
include/pmmc.h
View File

@@ -3957,10 +3957,18 @@ inline void pmmc_InterfaceSpeed(DoubleArray &dPdt, DoubleArray &P_x, DoubleArray
z = SurfaceVector(2*npts+p);
norm = sqrt(x*x+y*y+z*z);
// Save the surface values and normal vector
SurfaceValues(p) = -zeta/norm;
SurfaceVector(p) = x/norm;
SurfaceVector(npts+p) = y/norm;
SurfaceVector(2*npts+p) = z/norm;
if (norm > 0.0){
SurfaceValues(p) = -zeta/norm;
SurfaceVector(p) = x/norm;
SurfaceVector(npts+p) = y/norm;
SurfaceVector(2*npts+p) = z/norm;
}
else{
SurfaceValues(p) = 0.0;
SurfaceVector(p) = 0.0;
SurfaceVector(npts+p) = 0.0;
SurfaceVector(2*npts+p) = 0.0;
}
}
//.............................................................................
// Compute the average speed of the interface