Hi Jean-Francois,
in addition to the setAll... methods, the boundary types offer methods
(see also http://www.dumux.org/doxygen-stable/html-2.6/a00028.php)
setNeumann(int eqIdx)
setDirichlet(int pvIdx, int eqIdx)
setDirichlet(int eqIdx) (setDirichlet(i) is a shortcut for
setDirichlet(i,i))
That means that Neumann conditions are given with respect to equations
and Dirichlet conditions with respect to primary variables and equations.
In order to know what indices you can give, you can have a look at the
indices class for your model:
http://www.dumux.org/doxygen-stable/html-2.6/a00567.php
Your desired conditions amount to
values.setNeumann(energyEqIdx);
values.setNeumann(contiWEqIdx);
values.setNeumann(contiNEqIdx);
values.setDirichlet(pgIdx, contiGEqIdx);
or, shorter, to
values.setAllNeumann();
values.setDirichlet(pgIdx, contiGEqIdx);
If you do not define the indices like
enum {pgIdx = Indices::pgIdx};
in the beginning of your problem class, you have to use the full names
Indices::pgIdx, ...
Kind regards
Bernd
On 12/04/2014 01:03 PM, [email protected] wrote:
Hello again
I apologize for asking so many questions so quickly but on this one I am
clueless.
I try to understand how to define complex boundary conditions on a 3p3cni
problem.
My question:
How to define on one boundary a Dirichlet pressure conditions and a zero value
Neumann condition on temperature, W and N saturations?
This is a well posed set of condition to the best of my knowledge corresponding to an experimental reality where the pressure on an outlet is controlled.
Looking at the test 3p3cni columnxylolproblem.hh, I note that, as expected, a pressure
dirichlet conditions can be applied but "setALlDirichlet" hides how this
pressure condition is used in the system.
The 4 equations being energy + 3 mass conservations it is unclear to me how
such pressure condition is used.
I tried to go through the source to find it but no luck so far..[ I am not an
expert c++ programmer]
which bring me to reformulate my question:
How do I define a set of BC with Dirichlet on pressure and Neumann on
everything else in the context of a 3p3cni model?
I tried a "trick" :
Assigning all boundary condition to DIrichlet and then changing Temperature,
w and N saturation condition to Neumann but it does not seem to work...
Thanks a lot
JF
_______________________________________________
Dumux mailing list
[email protected]
https://listserv.uni-stuttgart.de/mailman/listinfo/dumux
--
_______________________________________________________________
Bernd Flemisch phone: +49 711 685 69162
IWS, Universität Stuttgart fax: +49 711 685 60430
Pfaffenwaldring 61 email: [email protected]
D-70569 Stuttgart url: www.hydrosys.uni-stuttgart.de
_______________________________________________________________
_______________________________________________
Dumux mailing list
[email protected]
https://listserv.uni-stuttgart.de/mailman/listinfo/dumux
