From dee99d968350b65cad39a763fa349f3b5b40dd43 Mon Sep 17 00:00:00 2001 From: Melanie Lipp <melanie.lipp@iws.uni-stuttgart.de> Date: Tue, 7 Sep 2021 11:11:50 +0200 Subject: [PATCH] [cvdstaggeredrefinedlocalresdidual] Fix factor for storage term at refinement boundary. --- dumux/assembly/cvdstaggeredrefinedlocalresidual.hh | 13 +++++++++++-- 1 file changed, 11 insertions(+), 2 deletions(-) diff --git a/dumux/assembly/cvdstaggeredrefinedlocalresidual.hh b/dumux/assembly/cvdstaggeredrefinedlocalresidual.hh index c17846e..386d94e 100644 --- a/dumux/assembly/cvdstaggeredrefinedlocalresidual.hh +++ b/dumux/assembly/cvdstaggeredrefinedlocalresidual.hh @@ -363,9 +363,18 @@ public: // multiply by 0.5 because we only consider half of a staggered control volume here storage *= 0.5*scv.volume()*extrusionFactor; - if (scvf.neighbor() && scvf.outside().level() > scvf.inside().level()) + // was used for coarse side to divide the volume to 1/4 + // if (scvf.neighbor() && scvf.outside().level() > scvf.inside().level()) + // { + // storage *= 0.5; + // } + + //we take the small volume (which is already half of what we want) and multiply by 0.5 + //when adding to the residual also by 0.5 + //to counteract 1/4 volume, we multiply by 2 + if(!scvf.isSimple() && scvf.isInsideFiner()) { - storage *= 0.5; + storage *= 2.0; } storage /= timeLoop_->timeStepSize(); -- GitLab