Commit d716f535 authored by Katharina Heck's avatar Katharina Heck Committed by Timo Koch
Browse files

[fix][richards] add extrusion factor in storage derivative

parent 8d375bab
......@@ -196,7 +196,7 @@ public:
static_assert(!FluidSystem::isCompressible(0),
"richards/localresidual.hh: Analytic Jacobian only supports incompressible fluids!");
const auto poreVolume = Extrusion::volume(scv)*curVolVars.porosity();
const auto poreVolume = Extrusion::volume(scv)*curVolVars.porosity()*curVolVars.extrusionFactor();
static const auto rho = curVolVars.density(0);
// partial derivative of storage term w.r.t. p_w
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment