Skip to content
Snippets Groups Projects
Commit fc3df5de authored by Yue Wang's avatar Yue Wang
Browse files

[poroelastic] fix the average density in localresidual

parent 3e150ea1
No related branches found
No related tags found
1 merge request!3459[poroelastic] fix the average density in localresidual
......@@ -92,7 +92,7 @@ public:
const auto& vv = elemVolVars[scv];
const auto phi = vv.porosity();
const auto rhoFluid = problem.spatialParams().effectiveFluidDensity(element, scv);
const auto rhoAverage = phi*rhoFluid + (1.0 - phi*vv.solidDensity());
const auto rhoAverage = phi*rhoFluid + (1.0 - phi)*vv.solidDensity();
// add body force
const auto& g = problem.spatialParams().gravity(scv.center());
......
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