Dear Bernd, Thanks a lot for your comments and help, they helped me a lot to solved my problem. >From each of your comments, I learned something about DUMUX and now I have a >clearer idea about using classes and template in DUMUX. Thanks.
In the near future, I will try to use DUMUX for 3d 2pni and 2p2cni. I wanted to ask you if you can suggest me some papers and reports that can be useful. For instance, I started reading: CO2 sequestration simulations: A comparison study between DuMuX and Eclipse, "Walter Vogel". Best regards Ali De : [email protected] [mailto:[email protected]] De la part de Bernd Flemisch Envoyé : mardi 21 août 2012 10:20 À : DuMuX User Forum Objet : Re: [DuMuX] RE : 2pni problem Dear Ali, you are right about the naming of the indices. This is an inconsistency that hopefully will be removed in future releases. Having a look at the error message or at 2pimmisciblefluidsystem.hh reveals that this fluid system only has an init method that accepts no arguments at all. So you should change the corresponding lines in injectionproblem2pni.hh (213-218 in the unmodified file) to FluidSystem::init(); Kind regards Bernd On 08/20/2012 08:52 PM, [email protected]<mailto:[email protected]> wrote: Dear Bernd, Thank you so much for your comment. I don't know why i didn't use this optim.opts in the first installation. I reinstall DUMUX this time using the optim.opts and I don't have the warnings (great). Now I believe that I have a correct error message. If i understood correctly my problem is coming from injectionspatailparameters.hh where the lphaseIdx is defined (test_2pni-error-1). I tried to change it wPhaseIdx = FluidSystem::wPhaseIdx as i think lPhaseIdx = FluidSystem::lPhaseIdx is only valid when we use fluidsystem for example h2on2. But, I obtained another error about : no matching function for call to 'Dumux::TwoPImmiscibleFluidSystem<Dumux::Properties::TTag::InjectionProblem2PNI>::init(double, double, int, double, double, int). Itriedto solve this problem by adding : SET_TYPE_PROP(InjectionProblem2PNI, FluidSystem, Dumux::TwoPImmiscibleFluidSystem<TypeTag>); but it didn't work. can you help me please with this problem. bests, Ali ________________________________ De : [email protected]<mailto:[email protected]> [[email protected]<mailto:[email protected]>] de la part de Bernd Flemisch [[email protected]<mailto:[email protected]> Date d'envoi : lundi 20 août 2012 17:54 À : DuMuX User Forum Objet : Re: [DuMuX] 2pni problem Hey Ali, this is not the full output of the build command, there is no error message in there. Most likely, you have a limit in the number of lines stored in your terminal window which you can enlarge. You also should consider to use the compiler options from the files debug.opts or optim.opts in the dumux directory (which you should pass to dunecontrol) such that you don't get the comparison warnings. I also suggest to change the two lines declaring the phase types to typedef Dumux::LiquidPhase<Scalar, Dumux::Brine<Scalar, Dumux::H2O<Scalar> > > type; typedef Dumux::GasPhase<Scalar, Dumux::SimpleCO2<Scalar> > type; Kind regards Bernd On 08/20/2012 05:37 PM, [email protected]<mailto:[email protected]> wrote: Dear DUMUX, I want to use DUMUX to simulate for Non-isothermal two-phase flow. I started first working on isothermal two-phase flow and I successfully modified tutorial_coupled (in DUMUX /tutorial). Now, I am working on test/boxmodels/ 2pni and I tried simply using brine and simpleco2 as fluids. In this case, I believe I can use this test without using a fluidsystem (like h2on2 used as an example in the original code). Therefore, I made this simple modification in injectionproblem2pni.hh: #if 1 // Set the wetting phase SET_PROP(InjectionProblem2PNI, WettingPhase) { private: typedef typename GET_PROP_TYPE(TypeTag, Scalar) Scalar; public: typedef Dumux::LiquidPhase<Scalar, Dumux::Brine<Scalar,H2O<Scalar>>> type; }; // Set the non-wetting phase SET_PROP(InjectionProblem2PNI, NonwettingPhase) { private: typedef typename GET_PROP_TYPE(TypeTag, Scalar) Scalar; public: typedef Dumux::GasPhase<Scalar, Dumux::SimpleCO2<Scalar>> type; }; #else SET_TYPE_PROP(InjectionProblem2PNI, FluidSystem, FluidSystems::H2OAir<typename GET_PROP_TYPE(TypeTag, Scalar)>); #endif However it didn't work and I received an error. I hope that you can help me with this problem and let me know please if I have to make further modification in this test to use brine and co2 instead of h2on2fluidsystem. I made some changes (ex. I thought it can come from how the indices are defined but I don't still find know how I can modify them) Thanks a lot for your help and I look forward to hearing from you soon. Please let me know if you need more information. **I send you the injectionproblem2pni.hh where I made my modifications and also the error that I received in terminal for your perusal. In .hh file , I commented my modifications by "my modification". Bests Ali _______________________________________________ Dumux mailing list [email protected]<mailto:[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]<mailto:[email protected]> D-70569 Stuttgart url: www.hydrosys.uni-stuttgart.de<http://www.hydrosys.uni-stuttgart.de> _____________________________________________________________________ _______________________________________________ Dumux mailing list [email protected]<mailto:[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]<mailto:[email protected]> D-70569 Stuttgart url: www.hydrosys.uni-stuttgart.de<http://www.hydrosys.uni-stuttgart.de> _____________________________________________________________________
_______________________________________________ Dumux mailing list [email protected] https://listserv.uni-stuttgart.de/mailman/listinfo/dumux
