OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2025-02-25 18:55:30 -06:00
Code Issues Packages Projects Releases Wiki Activity

Implements the correct check for aquifer face connections to make sure face is connected to the boundary

Browse Source
This commit is contained in:
kel85uk 2018-04-05 15:10:10 +02:00
parent d8f8646452
commit 0df51a8797
1 changed files with 17 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

18
opm/autodiff/AquiferCarterTracy.hpp
View File

@ -308,6 +308,14 @@ namespace Opm
const auto& grid = eclState.getInputGrid();
cell_idx_ = connection.global_index;
auto globalCellIdx = ugrid.globalCell();
// for (auto globalCells : globalCellIdx){
// std::cout << "global id = " << globalCells << std::endl;
// }
// for (auto cellidx : cell_idx_){
// std::cout << "aqucell id = " << cellidx << std::endl;
// }
assert( cell_idx_ == connection.global_index);
assert( (cell_idx_.size() == connection.influx_coeff.size()) );
@ -318,10 +326,14 @@ namespace Opm
cell_depth_.resize(cell_idx_.size(), d0_);
alphai_.resize(cell_idx_.size(), 1.0);
faceArea_connected_.resize(cell_idx_.size(),0.0);
Scalar faceArea;
auto cell2Faces = Opm::UgGridHelpers::cell2Faces(ugrid);
auto faceCells = Opm::AutoDiffGrid::faceCells(ugrid);
// static_assert(decltype(faceCells)::dummy_error, "DUMP MY TYPE" );
// Translate the C face tag into the enum used by opm-parser's TransMult class
Opm::FaceDir::DirEnum faceDirection;
@ -355,7 +367,11 @@ namespace Opm
if (faceDirection == connection.reservoir_face_dir.at(idx))
{
faceArea_connected_.at(idx) = Opm::UgGridHelpers::faceArea(ugrid, faceIdx);
// Check now if the face is outside of the reservoir, or if it adjoins an inactive cell
// Do not make the connection if the product of the two cellIdx > 0. This is because the
// face is within the reservoir/not connected to boundary. (We still have yet to check for inactive cell adjoining)
faceArea = (faceCells(faceIdx,0)*faceCells(faceIdx,1) > 0)? 0. : Opm::UgGridHelpers::faceArea(ugrid, faceIdx);
faceArea_connected_.at(idx) = faceArea;
denom_face_areas += faceArea_connected_.at(idx);
}
}