mirror of
https://github.com/OPM/opm-simulators.git
synced 2025-02-25 18:55:30 -06:00
remove derivative in thermal flux for injectors
This commit is contained in:
parent
d8c0057601
commit
7553e22726
@ -681,10 +681,13 @@ namespace Opm
|
|||||||
fs.setDensity(phaseIdx, rho);
|
fs.setDensity(phaseIdx, rho);
|
||||||
const auto& h = FluidSystem::enthalpy(fs, paramCache, phaseIdx);
|
const auto& h = FluidSystem::enthalpy(fs, paramCache, phaseIdx);
|
||||||
fs.setEnthalpy(phaseIdx, h);
|
fs.setEnthalpy(phaseIdx, h);
|
||||||
|
cq_r_thermal *= this->extendEval(fs.enthalpy(phaseIdx)) * this->extendEval(fs.density(phaseIdx));
|
||||||
|
connectionRates[perf][Indices::contiEnergyEqIdx] += getValue(cq_r_thermal);
|
||||||
|
} else {
|
||||||
|
// compute the thermal flux
|
||||||
|
cq_r_thermal *= this->extendEval(fs.enthalpy(phaseIdx)) * this->extendEval(fs.density(phaseIdx));
|
||||||
|
connectionRates[perf][Indices::contiEnergyEqIdx] += Base::restrictEval(cq_r_thermal);
|
||||||
}
|
}
|
||||||
// compute the thermal flux
|
|
||||||
cq_r_thermal *= this->extendEval(fs.enthalpy(phaseIdx)) * this->extendEval(fs.density(phaseIdx));
|
|
||||||
connectionRates[perf][Indices::contiEnergyEqIdx] += Base::restrictEval(cq_r_thermal);
|
|
||||||
}
|
}
|
||||||
}
|
}
|
||||||
|
|
||||||
|
Loading…
Reference in New Issue
Block a user