OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2024-12-02 13:59:12 -06:00
Code Issues Packages Projects Releases Wiki Activity

Added initFacePressure() function.

Browse Source
This commit is contained in:
Atgeirr Flø Rasmussen 2012-05-21 10:10:52 +02:00
parent 2ad94fc573
commit cd72a50977
1 changed files with 35 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

35
opm/core/utility/initState_impl.hpp
View File

@ -262,6 +262,41 @@ namespace Opm
}
}
// Initialize face pressures to distance-weighted average of adjacent cell pressures.
template <class State>
void initFacePressure(const UnstructuredGrid& grid,
State& state)
{
const int dim = grid.dimensions;
const std::vector<double>& cp = state.pressure();
std::vector<double>& fp = state.facepressure();
for (int f = 0; f < grid.number_of_faces; ++f) {
double dist[2] = { 0.0, 0.0 };
double press[2] = { 0.0, 0.0 };
int bdy_idx = -1;
for (int j = 0; j < 2; ++j) {
const int c = grid.face_cells[2*f + j];
if (c >= 0) {
dist[j] = 0.0;
for (int dd = 0; dd < dim; ++dd) {
double diff = grid.face_centroids[dim*f + dd] - grid.cell_centroids[dim*c + dd];
dist[j] += diff*diff;
}
dist[j] = std::sqrt(dist[j]);
press[j] = cp[c];
} else {
bdy_idx = j;
}
}
if (bdy_idx == -1) {
fp[f] = press[0]*(dist[1]/(dist[0] + dist[1])) + press[1]*(dist[0]/(dist[0] + dist[1]));
} else {
fp[f] = press[(bdy_idx + 1) % 2];
}
}
}
} // anonymous namespace