Public Member Functions
void	InputParam (const ParamReservoir &rs_param, const ParamOutput &output_param)
	Input parameters from the internal param structure.

void	SetupIsoT ()
	Setup for Isothermal model.

void	SetupT ()
	Setup for thermal model.

void	Setup ()
	Setup the grid information and calculate the properties.

OCP_USI	GetGridNum () const

OCP_INT	GetActIndex (const USI &I, const USI &J, const USI &K) const

void	GetIJKGrid (USI &i, USI &j, USI &k, const OCP_USI &n) const

void	GetIJKBulk (USI &i, USI &j, USI &k, const OCP_USI &n) const

OCP_BOOL	IfUseVtk () const
	return if use vtk format for outputing

void	SetHexaherdronGridOrthogonal ()
	setup polyhedronGrid for orthogonal grid

void	SetHexaherdronGridCorner (const OCP_COORD &mycord)

void	SetupGridTag ()
	Setup grid tag.

void	OutputBaiscInfo () const
	Calculate and return basic informations for grid.

void	CalNumDigutIJK ()
	only used in structured grid

USI	GetNumDigitIJK () const
	Return numDigutIJK.

Protected Member Functions
void	SetupOrthogonalGrid ()
	Setup orthogonal grid.

void	CalDepthVOrthogonalGrid ()
	Calculate the depth and volume for orthogonal grid.

void	SetupNeighborOrthogonalGrid ()
	Setup the neighboring info for an orthogonal grid.

void	SetupCornerGrid ()
	Setup corner-point grid.

void	SetupBasicCornerGrid (const OCP_COORD &CoTmp)
	Setup dx,dy,dz,depth, v for a corner-point grid.

void	SetupNeighborCornerGrid (const OCP_COORD &CoTmp)
	Setup the neighboring info for a corner-point grid.

void	CalActiveGridIsoT (const OCP_DBL &e1, const OCP_DBL &e2)
	Calculate the activity of grid cells. More...

void	CalActiveGridT (const OCP_DBL &e1, const OCP_DBL &e2)
	Calculate the activity of grid cells for ifThermal model.

void	SetupGridLocation ()
	Setup Grid location for Structured grid.

Protected Attributes
USI	gridType
	Orthogonal or Corner grid.

OCP_USI	numGrid
	Number of all cells.

USI	nx
	Number of cells in x-direction.

USI	ny
	Number of cells in y-direction.

USI	nz
	Number of cells in z-direction.

vector< OCP_DBL >	coord
	Lines of a corner-point grid.

vector< OCP_DBL >	zcorn
	ZValues of a corner-point grid.

vector< USI >	gLocation
	Top face, bottom face, side face, numGrid.

vector< OCP_DBL >	tops
	Depth of center of grid cells: numGrid.

vector< OCP_DBL >	dx
	Size of cell in x-direction: numGrid.

vector< OCP_DBL >	dy
	Size of cell in y-direction: numGrid.

vector< OCP_DBL >	dz
	Size of cell in z-direction: numGrid.

vector< OCP_DBL >	v
	Volume of cells: numGrid.

vector< OCP_DBL >	depth
	Depth of center of grid cells: numGrid.

vector< OCP_DBL >	ntg
	Net to gross ratio of cells: numGrid.

vector< OCP_DBL >	poro
	Initial porosity of rock cells: numGrid.

vector< OCP_DBL >	kx
	Absolute permeability in x-direction: numGrid.

vector< OCP_DBL >	ky
	Absolute permeability in y-direction: numGrid.

vector< OCP_DBL >	kz
	Absolute permeability in z-direction: numGrid.

vector< OCP_DBL >	thconr
	Rock if Thermal conductivity: numGrid.

vector< USI >	SATNUM
	Identify SAT region: numGrid.

vector< USI >	PVTNUM
	Identify PVT region for the blackoil model: numGrid.

vector< USI >	ACTNUM

vector< USI >	ROCKNUM
	index of rock table for each grid: numGrid

vector< OCP_DBL >	SwatInit
	Initial water saturation.

vector< vector< GPair > >	gNeighbor
	Neighboring information of grid.

OCP_USI	activeGridNum
	Num of active grid.

vector< OCP_USI >	map_Act2All
	Mapping from active grid to all grid: activeGridNum.

vector< GB_Pair >	map_All2Act
	Mapping from grid to active all grid: numGrid.

OCP_USI	fluidGridNum
	Num of fluid grids.

vector< GB_Pair >	map_All2Flu
	Mapping from all grid to fluid grid: numGrid.

OCP_BOOL	useVTK {OCP_FALSE}
	If output in vtk format.

vector< OCPpolyhedron >	polyhedronGrid
	Coordinates of grid points.

vector< USI >	gridTag
	Tag of grid: dead, live(fluid), live(rock)

USI	numDigutIJK
	number of digits of maximum nx,ny,nz

Friends
class	Bulk

class	BulkConn

class	Well

class	ScalePcow

class	Out4RPT

class	Out4VTK

Detailed Description

Definition at line 88 of file Grid.hpp.

Member Function Documentation

◆ CalActiveGridIsoT()

void Grid::CalActiveGridIsoT	(	const OCP_DBL &	e1,
		const OCP_DBL &	e2
	)

protected

Calculate the activity of grid cells.

If porosity or volume of the grid cell is too small, then the cell is inactive.

Definition at line 250 of file Grid.cpp.

 {
     map_Act2All.reserve(numGrid);
     map_All2Act.resize(numGrid);
     OCP_USI count = 0;
     for (OCP_USI n = 0; n < numGrid; n++) {
         if (ACTNUM[n] == 0 || poro[n] * ntg[n] < e1 || v[n] < e2) {
             map_All2Act[n] = GB_Pair(OCP_FALSE, 0);
             ACTNUM[n]      = 0;
             continue;
         }
         map_Act2All.push_back(n);
         map_All2Act[n] = GB_Pair(OCP_TRUE, count);
         count++;
     }
     activeGridNum = count;
     if (numGrid > activeGridNum) {
         cout << "  Number of inactive cells is " << (numGrid - activeGridNum) << " ("
              << (numGrid - activeGridNum) * 100.0 / numGrid << "%)" << endl;
     }
  
     // fluid grid = active grid
     fluidGridNum = activeGridNum;
     map_All2Flu  = map_All2Act;
 }

References ACTNUM, map_Act2All, map_All2Act, ntg, numGrid, poro, and v.

◆ SetHexaherdronGridCorner()

void Grid::SetHexaherdronGridCorner ( const OCP_COORD & mycord )

setup polyhedronGrid for corner grid

Definition at line 392 of file Grid.cpp.

 {
     if (!useVTK) return;
  
     polyhedronGrid.reserve(numGrid);
     OCPpolyhedron tmpP(8);
  
     for (USI k = 0; k < nz; k++) {
         for (USI j = 0; j < ny; j++) {
             for (USI i = 0; i < nx; i++) {
                 tmpP.Points.push_back(mycord.cornerPoints[i][j][k].p4);
                 tmpP.Points.push_back(mycord.cornerPoints[i][j][k].p5);
                 tmpP.Points.push_back(mycord.cornerPoints[i][j][k].p6);
                 tmpP.Points.push_back(mycord.cornerPoints[i][j][k].p7);
                 tmpP.Points.push_back(mycord.cornerPoints[i][j][k].p0);
                 tmpP.Points.push_back(mycord.cornerPoints[i][j][k].p1);
                 tmpP.Points.push_back(mycord.cornerPoints[i][j][k].p2);
                 tmpP.Points.push_back(mycord.cornerPoints[i][j][k].p3);
                 polyhedronGrid.push_back(tmpP);
                 tmpP.Points.clear();
             }
         }
     }
 }

References numGrid, nx, ny, nz, polyhedronGrid, and useVTK.

Member Data Documentation

◆ ACTNUM

vector<USI> Grid::ACTNUM

protected

Indicate activity of grid from input file: numGrid. 0 = inactive, 1 = active.

Definition at line 179 of file Grid.hpp.

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

Public Member Functions

Protected Member Functions