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

Work in progress fixing multidim upwinding.

Browse Source
This commit is contained in:
Atgeirr Flø Rasmussen
2013-08-26 10:43:11 +02:00
parent 4166a0621c
commit a199d78476
2 changed files with 22 additions and 4 deletions
Download Patch File Download Diff File Expand all files Collapse all files

23
opm/core/tof/TofReorder.cpp
View File

@@ -80,6 +80,8 @@ namespace Opm
if (use_multidim_upwind_) { if (use_multidim_upwind_) {
face_tof_.resize(grid_.number_of_faces); face_tof_.resize(grid_.number_of_faces);
std::fill(face_tof_.begin(), face_tof_.end(), 0.0); std::fill(face_tof_.begin(), face_tof_.end(), 0.0);
face_part_tof_.resize(grid_.face_nodepos[grid_.number_of_faces]);
std::fill(face_part_tof_.begin(), face_part_tof_.end(), 0.0);
} }
num_tracers_ = 0; num_tracers_ = 0;
num_multicell_ = 0; num_multicell_ = 0;
@@ -148,6 +150,8 @@ namespace Opm
if (use_multidim_upwind_) { if (use_multidim_upwind_) {
face_tof_.resize(grid_.number_of_faces); face_tof_.resize(grid_.number_of_faces);
std::fill(face_tof_.begin(), face_tof_.end(), 0.0); std::fill(face_tof_.begin(), face_tof_.end(), 0.0);
face_part_tof_.resize(grid_.face_nodepos[grid_.number_of_faces]);
std::fill(face_part_tof_.begin(), face_part_tof_.end(), 0.0);
OPM_THROW(std::runtime_error, "Multidimensional upwind not yet implemented for tracer."); OPM_THROW(std::runtime_error, "Multidimensional upwind not yet implemented for tracer.");
} }
num_multicell_ = 0; num_multicell_ = 0;
@@ -362,9 +366,15 @@ namespace Opm
const double flux_face = std::fabs(darcyflux_[face]); const double flux_face = std::fabs(darcyflux_[face]);
face_term = 0.0; face_term = 0.0;
cell_term_factor = 0.0; cell_term_factor = 0.0;
int num_contrib = 0;
for (int ii = 0; ii < num_adj; ++ii) { for (int ii = 0; ii < num_adj; ++ii) {
const int f = adj_faces_[ii]; const int f = adj_faces_[ii];
const double influx_f = (grid_.face_cells[2*f] == upwind_cell) ? -darcyflux_[f] : darcyflux_[f]; const double influx_f = (grid_.face_cells[2*f] == upwind_cell) ? -darcyflux_[f] : darcyflux_[f];
if (influx_f <= 0.0) {
// We have no contribution from face f, it is an outflow face.
continue;
}
ASSERT(influx_f > 0.0);
const double omega_star = influx_f/flux_face; const double omega_star = influx_f/flux_face;
// SPU // SPU
// const double omega = 0.0; // const double omega = 0.0;
@@ -372,11 +382,18 @@ namespace Opm
// const double omega = omega_star > 0.0 ? std::min(omega_star, 1.0) : 0.0; // const double omega = omega_star > 0.0 ? std::min(omega_star, 1.0) : 0.0;
// SMU // SMU
const double omega = omega_star > 0.0 ? omega_star/(1.0 + omega_star) : 0.0; const double omega = omega_star > 0.0 ? omega_star/(1.0 + omega_star) : 0.0;
face_term += omega*face_tof_[f]; const int* f_nodes_beg = grid_.face_nodes + grid_.face_nodepos[f];
const int* f_nodes_end = grid_.face_nodes + grid_.face_nodepos[f + 1];
for (const int* f_iter = f_nodes_beg; f_iter < f_nodes_end; ++f_iter) {
if (face_part_tof_[*f_iter] > 0.0) {
face_term += omega * face_part_tof_[*f_iter];
cell_term_factor += (1.0 - omega); cell_term_factor += (1.0 - omega);
++num_contrib;
} }
face_term /= double(num_adj); }
cell_term_factor /= double(num_adj); }
face_term /= double(num_contrib);
cell_term_factor /= double(num_contrib);
} }

1
opm/core/tof/TofReorder.hpp
View File

@@ -109,6 +109,7 @@ namespace Opm
// For multidim upwinding: // For multidim upwinding:
bool use_multidim_upwind_; bool use_multidim_upwind_;
std::vector<double> face_tof_; // For multidim upwind face tofs. std::vector<double> face_tof_; // For multidim upwind face tofs.
std::vector<double> face_part_tof_; // For multidim upwind face tofs.
mutable std::vector<int> adj_faces_; // For multidim upwind logic. mutable std::vector<int> adj_faces_; // For multidim upwind logic.
}; };