mirror of
https://github.com/OPM/opm-simulators.git
synced 2025-02-25 18:55:30 -06:00
Merge pull request #4282 from totto82/fixInjTemp2
remove derivative in thermal flux for injectors
This commit is contained in:
commit
fa1ea7f997
@ -678,10 +678,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);
|
|
||||||
}
|
}
|
||||||
}
|
}
|
||||||
|
|
||||||
|
@ -74,7 +74,7 @@ add_test_compare_parallel_simulation(CASENAME spe1_thermal
|
|||||||
FILENAME SPE1CASE2_THERMAL
|
FILENAME SPE1CASE2_THERMAL
|
||||||
SIMULATOR flow
|
SIMULATOR flow
|
||||||
ABS_TOL ${abs_tol}
|
ABS_TOL ${abs_tol}
|
||||||
REL_TOL ${rel_tol}
|
REL_TOL ${coarse_rel_tol_parallel}}
|
||||||
DIR spe1
|
DIR spe1
|
||||||
TEST_ARGS --linear-solver-reduction=1e-7 --tolerance-cnv=5e-6 --tolerance-mb=1e-8)
|
TEST_ARGS --linear-solver-reduction=1e-7 --tolerance-cnv=5e-6 --tolerance-mb=1e-8)
|
||||||
|
|
||||||
|
Loading…
Reference in New Issue
Block a user