Ensure tracer flow for cross-flowing injector connections is included

This commit is contained in:
Vegard Kippe 2024-12-06 15:26:19 +01:00
parent 5d4f578319
commit 2b59222abc

View File

@ -515,6 +515,10 @@ protected:
}
else if (rate_f < 0) {
for (int tIdx = 0; tIdx < tr.numTracer(); ++tIdx) {
// Store _injector_ tracer rate for cross-flowing well connections (can be done here since WTRACER is constant)
this->wellTracerRate_.at(std::make_pair(eclWell.name(),this->name(tr.idx_[tIdx]))) += rate_f*wtracer[tIdx];
this->wellFreeTracerRate_.at(std::make_pair(eclWell.name(),this->wellfname(tr.idx_[tIdx]))) += rate_f*wtracer[tIdx];
// Production of free tracer
tr.residual_[tIdx][I][0] -= rate_f * tr.concentration_[tIdx][I][0];