Adapt wdfac to wpimult

This commit is contained in:
Tor Harald Sandve
2023-11-23 15:22:51 +01:00
parent 8f6b2b99e7
commit e1aa233026

View File

@@ -1332,7 +1332,7 @@ namespace Opm
double Q = std::abs(ws.perf_data.phase_rates[perf*pu.num_phases + pu.phase_pos[Gas]]);
const auto& connection = this->well_ecl_.getConnections()[ws.perf_data.ecl_index[perf]];
double Kh = connection.Kh();
double scaling = 3.141592653589 * Kh;
double scaling = 3.141592653589 * Kh * connection.wpimult();
const unsigned gas_comp_idx = Indices::canonicalToActiveComponentIndex(FluidSystem::gasCompIdx);
wi[gas_comp_idx] = 1.0/(1.0/(trans_mult * this->well_index_[perf]) + (Q/2 * d / scaling));
return wi;