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

added possibility for doing well contributions separately, still error

Browse Source
This commit is contained in:
hnil 2022-06-20 14:45:55 +02:00 committed by Atgeirr Flø Rasmussen
parent cfd727b012
commit e245f0ae50
1 changed files with 17 additions and 11 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

28
opm/models/discretization/common/linearizertpfa.hh
View File

@ -432,8 +432,8 @@ private:
}
}
void setResAndJacobi(VectorBlock& res,MatrixBlock& bMat,const ADVectorBlock& resid){
public:
void setResAndJacobi(VectorBlock& res,MatrixBlock& bMat,const ADVectorBlock& resid) const{
for (unsigned eqIdx = 0; eqIdx < numEq; eqIdx++)
res[eqIdx] = resid[eqIdx].value();
@ -447,9 +447,10 @@ private:
}
}
}
private:
void linearizeGlobalTPFA_()
{
const bool well_local = true;
resetSystem_();
unsigned numCells = model_().numTotalDof();
for(unsigned globI = 0; globI < numCells; globI++){
@ -476,14 +477,16 @@ private:
residual_[globI] += res;
jacobian_->addToBlock(globI, globI, bMat);
// wells sources for now (should be moved out)
res = 0.0;
bMat = 0.0;
adres = 0.0;
LocalResidual::computeSource(adres, problem_(), globI, 0);
adres *= -volume;
setResAndJacobi(res, bMat, adres);
residual_[globI] += res;
jacobian_->addToBlock(globI, globI, bMat);
if(well_local){
res = 0.0;
bMat = 0.0;
adres = 0.0;
LocalResidual::computeSource(adres, problem_(), globI, 0);
adres *= -volume;
setResAndJacobi(res, bMat, adres);
residual_[globI] += res;
jacobian_->addToBlock(globI, globI, bMat);
}
for(const auto& globJ: neighbours){
assert(globJ != globI);
res = 0.0;
@ -508,6 +511,9 @@ private:
bMat *= -1.0;
jacobian_->addToBlock(globJ, globI, bMat);
}
if(not(well_local)){
problem_().wellModel().addReseroirSourceTerms(residual_,*jacobian_);
}
}
// before the first iteration of each time step, we need to update the