Dear Jean,
I will try to give you some hints:
1- No gas dissolution is allowed in NAPL (xna = 1e-10 while it could be
something like
xna= xgg * pg_ / (FluidSystem::fugacityCoefficient(fluidState_,
wPhaseIdx,gCompIdx) * pn_); )
You are right, in the 3p3cvolumevariables.hh no dissolution in the NAPL
phase is allowed, if the useConstraintSolver property is set to false.
My suggestion is setting that property to true and defining the
fugacityCoefficient appropriately in the FluidSystem you are developing,
that should take care of dissolution. Set the property in your problem file:
SET_BOOL_PROP(YourProblemTypeTag, UseConstraintSolver, true);
2....itt SEEMS pretty easy to add ( to easy?)
Putting the proper value for xna in 3p3cvolume variable.hh should do the job
Yes, I think the best way is using the constraint solver, since this
functionality is also used in the other compositional models 2p2c, MPNC.
What are the reasons (if any) that this feature is explicitely not available?
I don't really know. I guess the kind of problems that were solved with
this model never required dissolution in the NAPL phase until now. So
thoroughly check your results to see that you have the correct
dissolution behavior.
best regards
Alex
On 12/09/2014 10:45 PM, [email protected] wrote:
Hello there
I was wondering if I could get some feedback and pointers on what work and dont
in the 3p3c solver regarding gas dissolution in napl phase. [ and gas diffusion
as well].
Looking at the code I notice that:
1- No gas dissolution is allowed in NAPL (xna = 1e-10 while it could be
something like
xna= xgg * pg_ / (FluidSystem::fugacityCoefficient(fluidState_,
wPhaseIdx,gCompIdx) * pn_); )
2....itt SEEMS pretty easy to add ( to easy?)
Putting the proper value for xna in 3p3cvolume variable.hh should do the job
or use the multicomposotional solver
The fact that these "apparently easy solution" are disabled puzzle me. So my
question:
What are the reasons (if any) that this feature is explicitely not available?
(The fact that this is mentionned in the comments means that somebody gave it
some thoughts and decided not to do it...)
On the other hand, if it is just a matter of spending some time and doing some
testing I can do it and let you know what I find....
I will just appreciate to have some pointers : what is known to work ( or dont
work) in order to avoid wandering in dead end that I could avoid....
Thanks a lot in advance...
Jean-Francois
_______________________________________________
Dumux mailing list
[email protected]
https://listserv.uni-stuttgart.de/mailman/listinfo/dumux
--
Alexander Kissinger
Institut für Wasser- und Umweltsystemmodellierung
Lehrstuhl für Hydromechanik und Hydrosystemmodellierung
Pfaffenwaldring 61
D-70569 Stuttgart
Telefon: +49 (0) 711 685-64729
E-Mail: [email protected]
_______________________________________________
Dumux mailing list
[email protected]
https://listserv.uni-stuttgart.de/mailman/listinfo/dumux
