Hi Dennis,
Thanks a lot, it works now.
Best regards,
Nikolai
From: Dumux on behalf of Dennis
Gläser
Sent: Tuesday, August 21, 2018 2:03:46 PM
To: dumux@listserv.uni-stuttgart.de
Subject: Re: [DuMuX] Check global conservation error
Hi Nikolai,
you have
Looking forward for your feedback.
Thanks,
Nikolai
*From:* Dumux on behalf of
Dennis Gläser
*Sent:* Tuesday, August 21, 2018 11:39:34 AM
*To:* dumux@listserv.uni-stuttgart.de
*Subject:* Re: [DuMuX] Check global co
for your feedback.
Thanks,
Nikolai
From: Dumux on behalf of Dennis
Gläser
Sent: Tuesday, August 21, 2018 11:39:34 AM
To: dumux@listserv.uni-stuttgart.de
Subject: Re: [DuMuX] Check global conservation error
Hi Nikolai,
I think your implementation is almost cor
Hi Nikolai,
I think your implementation is almost correct. If you want to balance
masses, you have to add the density to your upwind terms though, i.e.:
auto upwindTermN = [](const auto& volVars) { return
volVars.mobility(nPhaseIdx)*volVars.density(nPhaseIdx); };
...
Also, you have to
Dear DuMuX experts,
Please advise how can I implement a global conservation error check, that is,
evaluate the residual ( mass(t^{n+1}) - mass(t^n) ) / dt + sum of fluxes over
the domain boundaries.
So far I can evaluate the mass in the domain by looping through the elements as
in void