WellInfo Class Reference

#include <wellinfo.h>

Collaboration diagram for WellInfo:

[legend]

List of all members.

Public Types
enum	ConditionType { FLUX, PRESSURE, SOURCE_INJ_RATE, SOURCE_FIXED_PRESSURE }
Public Member Functions
	WellInfo (Point3D p1, Point3D p2, double value, ConditionType bcType)
	~WellInfo ()
bool	isPointInWell (const VecDouble &p1) const
bool	isPointInWell (const Point3D &p1) const
double	getInjRate ()
double	getPressureBC ()
double	getSourceInjRate ()
double	getSourceFixedPressure ()
void	adjustBoundaryWithGrid (const Point3D &DX)
bool	isInnerCell (Point3D &p, Point3D &DX)
ConditionType	getBCType ()
bool	isSourceWell () const
const VecDouble &	getTransportBC () const
void	setTransportBC (const VecDouble &v)
double	volume ()
double	area ()
Public Attributes
Point3D	P
Point3D	Q
ConditionType	m_bcType
double	m_value
VecDouble	m_transBC

---

Detailed Description

Store well information for a particular well. The wells specify boundary conditions that can be of three type

FLUX:= The wells are represented as holes in the mesh where we specify the Neumman conditions for the pressure system

PRESSURE:= The wells are represented as holes where we specify the Dirichlet boundary conditions for pressure For FLUX and PRESSURE, the wells have boundary conditions for the transport variables (Saturations, Concentrations, etc) SOURCE:= The well is represented as a source term inside, the region of the well. There are no holes in this case The conditions for the transport is specified using Fixed Values Conditions in all the region ocupied by the wells

Definition at line 24 of file wellinfo.h.

---

Member Enumeration Documentation

enum WellInfo::ConditionType

Enumerator:

FLUX
PRESSURE
SOURCE_INJ_RATE
SOURCE_FIXED_PRESSURE

Definition at line 27 of file wellinfo.h.

{FLUX,PRESSURE,SOURCE_INJ_RATE,SOURCE_FIXED_PRESSURE};

---

Constructor & Destructor Documentation

WellInfo::WellInfo	(	Point3D	p1,
		Point3D	p2,
		double	value,
		ConditionType	bcType
	)

Definition at line 3 of file wellinfo.cpp.

  :P(p1),Q(p2),m_bcType(bcType),m_value(value)
{
  assert(p1[0] < p2[0] &&
         p1[1] < p2[1] &&
         p1[2] < p2[2]);
}

WellInfo::~WellInfo

(

)

[inline]

Definition at line 35 of file wellinfo.h.

{}

---

Member Function Documentation

void WellInfo::adjustBoundaryWithGrid

(

const Point3D &

DX

)

Definition at line 46 of file wellinfo.cpp.

{
  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];

}

double WellInfo::area

(

)

Definition at line 36 of file wellinfo.cpp.

{
  Point3D DX = Q;
  DX -=P;

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

ConditionType WellInfo::getBCType

(

)

[inline]

Definition at line 44 of file wellinfo.h.

{return m_bcType;}

double WellInfo::getInjRate

(

)

[inline]

Definition at line 38 of file wellinfo.h.

{assert(m_bcType==FLUX);return m_value;}

double WellInfo::getPressureBC

(

)

[inline]

Definition at line 39 of file wellinfo.h.

{assert(m_bcType==PRESSURE);return m_value;}

double WellInfo::getSourceFixedPressure

(

)

[inline]

Definition at line 41 of file wellinfo.h.

{assert(m_bcType==SOURCE_FIXED_PRESSURE);return m_value;}

double WellInfo::getSourceInjRate

(

)

[inline]

Definition at line 40 of file wellinfo.h.

{assert(m_bcType==SOURCE_INJ_RATE);return m_value;}

const VecDouble& WellInfo::getTransportBC

(

)

const [inline]

Definition at line 46 of file wellinfo.h.

{assert(isSourceWell());return m_transBC;}

bool WellInfo::isInnerCell	(	Point3D &	p,
		Point3D &	DX
	)

Definition at line 72 of file wellinfo.cpp.

{
  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);
}

bool WellInfo::isPointInWell

(

const Point3D &

p1

)

const

Definition at line 11 of file wellinfo.cpp.

{
  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

Definition at line 18 of file wellinfo.cpp.

{
  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]);
          
}

bool WellInfo::isSourceWell

(

)

const [inline]

Definition at line 45 of file wellinfo.h.

{return m_bcType == SOURCE_INJ_RATE;}

void WellInfo::setTransportBC

(

const VecDouble &

v

)

[inline]

Definition at line 47 of file wellinfo.h.

{m_transBC=v;}

double WellInfo::volume

(

)

Definition at line 28 of file wellinfo.cpp.

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

---

Member Data Documentation

ConditionType WellInfo::m_bcType

Definition at line 29 of file wellinfo.h.

VecDouble WellInfo::m_transBC

Definition at line 31 of file wellinfo.h.

double WellInfo::m_value

Definition at line 30 of file wellinfo.h.

Point3D WellInfo::P

Definition at line 28 of file wellinfo.h.

Point3D WellInfo::Q

Definition at line 28 of file wellinfo.h.

---

The documentation for this class was generated from the following files:

• wellinfo.h
• wellinfo.cpp