Hi Ali,
it is nice to hear that you already use Dumux for advanced problems.
What kind of linear solver do you use? From our experience, the standard
linear solvers will not work nicely for more complex problems in parallel.
Since Dumux 2.3, we have a backend for the ISTL AMG solver. We are using
that for our parallel CO2 computations. You can have a look at how to
use it from test/implicit/1p/test_cc1pwithamg and the problem file in
there. You first need to install the DUNE module dune-pdelab,
http://www.dune-project.org/pdelab/, version 1.0 or 1.1 and to apply
patches to dune-istl and dune-pdelab. See dumux/patches/README for
details. You then need to rerun dunecontrol.
For heterogeneous parameters, it can also be benefitial to use the
cell-centered version instead of box (vertex-centered). You can have a
look at the co2 test problem file to see how this is done.
Currently, the salinity is a public static variable of the Brine
component and not intended to be varying over space. If you use the
BrineCO2FluidSystem in the standard way, the density and viscosity
values of the brine will be tabulated depending on this salinity. You
can use an untabulated version, see the documentation of
BrineCO2FluidSystem at line 745ff. Beware that this can slow down your
computation.
Once you use an untabulated Brine component, you could manipulate the
static salinity value. In a quick and dirty way, you could do this at
the very beginning of the update function in co2volumevariables.hh, line
122, for example
if (element.geometry().center()[2] < 1000.0)
FluidSystem::Brine::salinity = 0.1;
else
FluidSystem::Brine::salinity = 0.05;
Then all the secondary variables like density would be calculated
depending on this salinity value.
Please let us know how this works out.
Kind regards
Bernd
On 10/22/2013 01:59 AM, Ali NOWAMOOZ wrote:
Dear Dumux,
First, I would like to congratulate all DuMux team for this new version.
I have a question about the new version of DuMux :
I would like to know If any progress has been made on the parallel computing in
this version comparing to version 2.3. We could install dumux on Clumeq (server
of laval universityhttps://www.calculquebec.ca/index.php/en/) and run a CO2
2p2c model on 96 processors (12 nodes * 8 processors). However, we recognized a
problem with the linear solvers of DUNE that couldn't converge (actually, it
converged but only with time steps lower than 500 s). I found some explanations
on the dumux website at this
addresshttp://www.mail-archive.com/[email protected]/msg00006.html.
According to this mail, the parallel version of dumux was experimental at that
moment. Please let me know if it is also the case in the new version. If you
need any more information, please let me know. I can also run a model and send
you the results.
I have also a question about CO2 model. I would like to know if it is possible to define different salinities in the domain. I don't want to simulate the salt transport. I only need to distinct fresh water from brine in my domain.
Best regards
Ali NOWAMOOZ
_______________________________________________
Dumux mailing list
[email protected]
https://listserv.uni-stuttgart.de/mailman/listinfo/dumux
--
_______________________________________________________________
!!!! CMWR 2014: 10th - 13th June 2014 in Stuttgart !!!!
Please visit www.cmwr14.de
_______________________________________________________________
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
