Dear Alex, Thanks a lot. This is very helpful
I have done a few test [ nothing systematic and definitive at that point] and find reasonable results. Jean-Francois ________________________________________ De : Alexander Kissinger [[email protected]] Envoyé : mercredi 10 décembre 2014 08:09 À : DuMuX User Forum; LEON Jean-Francois Objet : Re: [DuMuX] Fugacity of NAPL in 3p3c system 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
