Hello everyone,
I am new to the Dumux community and therefore, also after completing the
tutorials, still quite unexperienced.
Since over one week I am focusing on material laws and capillary pressure
especially, but there are two questions
for which I cant find any answers. The versions I am working with right now
is dune 2.3 and dumux 2.7
I am working with the Brooks-Corey material law and call it via the
TwoPAdapter ("2padapter.hh"), but there are two things I dont understand so
far:
1. In function static void capillaryPressures(): Contrary to the
commentary lines, the nonwetting phase gets nothing added and the capillary
pressure is actually
subtracted from the pressure of the wetting phase, what is the purpose of
that? As far as I know, the capillary pressure is defined as the difference
between
the pressures from nonwetting and wetting phase i.e. pc = pn - pw
The passage of code I am talking about is lines 61-65 in
dumux/material/fluidmatrixinteractions/mp/2padapter.hh:
"
l. 61: // non-wetting phase gets the capillary pressure added
l. 62: values[nPhaseIdx] = 0;
l. 63:
l. 64: // wetting phase does not get anything added
l. 65: values[wPhaseIdx] = - TwoPLaw::pc(params,
state.saturation(wPhaseIdx));
"
2. The description in file 2padater.hh states, that the TwoPAdapter can
be used to apply Brooks-Corey to M-phase fluid systems, but I dont get it
how that can be accomplished
since the number of phases is defined inside class TwoPAdapter as 2 (l. 51
"enum { numPhases = 2};"). It is also referred to another file called
"MpBrooksCoreyMaterialParams" (l. 42)
but I cant find the respective file anywhere. On the other hand I found in
the same folder "mplinearmaterial.hh" and "mplinearmaterialparams.hh" which
make the linear material law
applicable for multiple phases. Is it possible that multiphase Brooks-Corey
was not fully implemented yet in dumux 2.7? If yes, what are the
alternatives? I was thinking about defining
temporary fluidsystem, each containing two phases with decreasing
wettability and then calling the TwoPAdapter for those systems, would that
be possible?
Thanks already for any support, I appreciate any advice!
Best Regards
Patrick
_______________________________________________
Dumux mailing list
[email protected]
https://listserv.uni-stuttgart.de/mailman/listinfo/dumux
