Skip to content
Snippets Groups Projects
Commit 04821e5d authored by Thomas Fetzer's avatar Thomas Fetzer
Browse files

[multidomain] Simplify pressure correction term

parent 4c4df4ae
No related branches found
No related tags found
1 merge request!252Fix/multidomain normal momentum
......@@ -147,10 +147,8 @@ public:
// if p.n comes from the pm
if (bcTypes.isCouplingNeumann(momentumYIdx) || bcTypes.isCouplingMortar(momentumYIdx))
{
DimVector pressureCorrection(intersection.centerUnitOuterNormal());
pressureCorrection *= volVars.pressure();
this->residual_[scvIdx][momentumYIdx] += pressureCorrection[momentumYIdx]
* boundaryVars.face().area;
this->residual_[scvIdx][momentumYIdx] += volVars.pressure()
* boundaryVars.face().normal[momentumYIdx];
Valgrind::CheckDefined(this->residual_[scvIdx][momentumYIdx]);
}
......
......@@ -157,10 +157,8 @@ public:
// if p.n comes from the pm
if (bcTypes.isCouplingNeumann(momentumYIdx) || bcTypes.isCouplingMortar(momentumYIdx))
{
DimVector pressureCorrection(intersection.centerUnitOuterNormal());
pressureCorrection *= volVars.pressure();
this->residual_[scvIdx][momentumYIdx] += pressureCorrection[momentumYIdx]
* boundaryVars.face().area;
this->residual_[scvIdx][momentumYIdx] += volVars.pressure()
* boundaryVars.face().normal[momentumYIdx];
Valgrind::CheckDefined(this->residual_[scvIdx][momentumYIdx]);
}
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment