wellinfo.cpp

Go to the documentation of this file.

#include "wellinfo.h"
#include "numericmethods.h"
WellInfo::WellInfo(Point3D p1,Point3D p2,double value,ConditionType bcType) 
  :P(p1),Q(p2),m_bcType(bcType),m_value(value)
{
  assert(p1[0] < p2[0] &&
         p1[1] < p2[1] &&
         p1[2] < p2[2]);
}

bool WellInfo::isPointInWell(const Point3D &p1) const 
{
  return (p1[0] >= P[0] && p1[0] <= Q[0] &&
          p1[1] >= P[1] && p1[1] <= Q[1] &&
          p1[2] >= P[2] && p1[2] <= Q[2]);
          
}
bool WellInfo::isPointInWell(const VecDouble &p1) const 
{
  assert(p1.size() ==3);
  return (p1(0) >= P[0] && p1(0) <= Q[0] &&
          p1(1) >= P[1] && p1(1) <= Q[1] &&
          p1(2) >= P[2] && p1(2) <= Q[2]);
          
}


double WellInfo::volume() 
{
  return fabs( (Q[0]-P[0])*(Q[1]-P[1])*(Q[2]-P[2]) );
}




double WellInfo::area() 
{
  Point3D DX = Q;
  DX -=P;

  return 2*(DX[0]*DX[1] + DX[0]*DX[2] + DX[1]*DX[2]);
  
}


void WellInfo::adjustBoundaryWithGrid(const Point3D &DX) 
{
  double tol = NumericMethods::min(DX[0],DX[1],DX[2])/10.0;
  P[0] = floor(P[0]/DX[0])*DX[0];
  P[1] = floor(P[1]/DX[1])*DX[1];
  P[2] = floor(P[2]/DX[2])*DX[2];

  Q[0] = ceil(Q[0]/DX[0])*DX[0];
  Q[1] = ceil(Q[1]/DX[1])*DX[1];
  Q[2] = ceil(Q[2]/DX[2])*DX[2];

  //check if the well has at least one cell. If the user specifies a small
  //well region, than the well could have no cells, this is considered a
  //mistake, after wall if the user specified a well he wants that well
  //in the mesh, right ? So in that case we adjust the well
  //so it can contain at least one cell of the mesh
  if (NumericMethods::a_equal(P[0],Q[0],tol))
    Q[0] = P[0]+DX[0];
  if (NumericMethods::a_equal(P[1],Q[1],tol))
    Q[1] = P[1]+DX[1];
  if (NumericMethods::a_equal(P[2],Q[2],tol))
    Q[2] = P[2]+DX[2];

}


bool WellInfo::isInnerCell(Point3D &cellCenter,Point3D &DX) 
{
  if (!isPointInWell(cellCenter))
      return false;

  Point3D p = P;
  p+=DX;

  Point3D q = Q;
  q-=DX;

  if (p[0] >= q[0] || p[1] >= q[1] ||p[2] >= q[2]) 
  {
    return false;
  }


  return NumericMethods::isInCube(p,q,cellCenter);
}

---