Hi Nikolai,
maybe the easiest way is to also implement an priVarFormulation() function*
*in your TwoP volVars, as it is for example done for the 2p2c volVars.
Then you can access the formulation using
ElementVolumeVariables::VolumeVariables::priVarFormulation()
in your darcyslaw.
I guess this is related to the previous discussion where you mentioned
that you had to change
darcyslaw as follows:
* if (!scvf.boundary()) {**
****
** // Inside the domain calculate the phase flux using
phase pressures (with pn = pw + pc)**
** pInside = insideVolVars.pressure(phaseIdx);**
** pOutside = outsideVolVars.pressure(phaseIdx);**
** }**
** else {**
** // Do not include the contributions from capillary
forces for boundary fluxes**
****//static_assert( (VolumeVariables::priVarFormulation() ==
TwoPFormulation::p0s1), "darcyslaw.hh: BC implemented for p0s1
formulation only!"); // uncommenting this yields an error**
** // WARNING: implemented for p1s0 formulation only!**
** pInside = insideVolVars.pressure(1);**
** pOutside = outsideVolVars.pressure(1);**
** }*
Which means that you do not want to have a capillary flux over your
domain boundary.
I think the easier alternative would be to use a solution-dependent
Neumann condition in your problem
* NumEqVector neumann(const Element& element,**
** const FVElementGeometry& fvGeometry,**
** const ElementVolumeVariables& elemVolVars,**
** const SubControlVolumeFace& scvf) const
*where you can get the nonwetting phase pressure of the boundary cell
via the *elemVolVars*,
and the fluxes can then be calculated in this function similar to
darcyslaw but now you can specify
that pc=0 for your boundary fluxes (i.e. pw = pn).
Kind regards,
Martin
On 08/24/2018 09:34 AM, Nikolai Andrianov wrote:
Hi Martin,
I would like to know from within
dumux/discretization/cellcentered/tpfa/darcyslaw.hh if
TwoPFormulation::p0s1 or TwoPFormulation::p1s0 is used.
Thanks,
Nikolai
------------------------------------------------------------------------
*From:* Martin Schneider <[email protected]>
*Sent:* Thursday, August 23, 2018 4:26:20 PM
*To:* DuMuX User Forum; Nikolai Andrianov
*Subject:* Re: [DuMuX] Formulation used in darcyslaw.hh
Dear Nikolai,
could you please specify what you mean with "formulation".
Do you mean which implementation of CCTpfaDarcysLaw?
Thanks,
Martin
On 08/23/2018 03:39 PM, Nikolai Andrianov wrote:
Dear DuMuX experts,
Please advise how can I identify, which formulation is used in
dumux/discretization/cellcentered/tpfa/darcyslaw.hh.
Thanks,
Nikolai
_______________________________________________
Dumux mailing list
[email protected]
https://listserv.uni-stuttgart.de/mailman/listinfo/dumux
--
M.Sc. Martin Schneider
University of Stuttgart
Institute for Modelling Hydraulic and Environmental Systems
Department of Hydromechanics and Modelling of Hydrosystems
Pfaffenwaldring 61
D-70569 Stuttgart
Tel: (+49) 0711/ 685-69159
Fax: (+49) 0711/ 685-60430
E-Mail:[email protected]
--
M.Sc. Martin Schneider
University of Stuttgart
Institute for Modelling Hydraulic and Environmental Systems
Department of Hydromechanics and Modelling of Hydrosystems
Pfaffenwaldring 61
D-70569 Stuttgart
Tel: (+49) 0711/ 685-69159
Fax: (+49) 0711/ 685-60430
E-Mail: [email protected]
_______________________________________________
Dumux mailing list
[email protected]
https://listserv.uni-stuttgart.de/mailman/listinfo/dumux
