Dear developers,
I'm a postdoc student and I was asked to study the viscous finger
phenomenon due to water/oil interface.
After literature review I'm convinced that dumux is a good choice. So,
before jump into the problem, I'd like to ask you few questions to
understand if it could be the
Dear developers and users,
looking at equations 3.1 and 3.2 of the handbook 2.12 it seems that no
capillarity effect is taken into account.
Indeed, I would expect something like the equations reported in the
attachment. Nevertheless,
values for capillarity pressure are loaded in the
Dear developers and users,
I have some other questions.
1) I successfully ran the tutorial_implicit example from the handbook
2.12. Now, I need to modify such test case in order to be consistent
with my requirements.
In particular, I need to use a different relative permeability and
Good afternoon,
I added "my laws" for the capillarity pressure and relative
permeability. I would like, at least at first, to use them without any
regularization so I call them directly "as they are" for example as a
modified copy of the Brooks & Corey.
So now my question is: if I don't use
sh file. You can
tell ALUGrid to use this file through the input file
[Grid]
File = mygrid.msh
For more informations on grid managers and how to use them, please
have a look at the Dumux handbook
(http://www.dumux.org/documents/dumux-handbook-2.12.pdf).
Best wishes
Timo
On 24.10.18
,
Lorenzo
On 17/10/2018 16:53, Timo Koch wrote:
Hi Lorenzo,
On 17.10.2018 15:21, lc wrote:
Dear developers and users,
I have some other questions.
1) I successfully ran the tutorial_implicit example from the
handbook 2.12. Now, I need to modify such test case in order to be
consistent
Hello,
I'd like to ask if CUDA or any graphical processing is supported in DUMUX.
Kind regards,
Lorenzo
___
Dumux mailing list
Dumux@listserv.uni-stuttgart.de
https://listserv.uni-stuttgart.de/mailman/listinfo/dumux
Good morning,
I'd like to ask how to impose an initial condition on the water
saturation and not on the oil saturation.
Actually, what I implemented is the following:
const auto pos = fvGeometry.subContVol[scvIdx].global;
double x_0 = 0.;
double smin = 0.2121;
setup?
Thanks a lot,
Lorenzo
On 02/11/2018 13:56, Timo Koch wrote:
Hi Lorenzo,
it all depends on your setup. Probably something wrong with your
custom material law.
What have you tried so far to fix your issue?
Timo
On 2. Nov 2018, at 11:53, lc <mailto:lorenzo.camp...@uniroma1.it>&
Good morning,
I'd like to ask which types of parallelism (mpi/openmp/cuda/...) is
implemented in dumux and how is it possible to run in parallel
(compilation settings)?
Kind regards,
Lorenzo
___
Dumux mailing list
Dumux@listserv.uni-stuttgart.de
Dear developers and users,
thanks for your fast and appropriate.
Once again I have a doubt. How can I apply an initial condition (for
example on saturation profile) which depends on the x and y spatial
coordinates of the domain?
//! Evaluates the initial value for a control volume. For this
Dear all,
I would like to ask few more questions about Dumux.
1) About *Restarting*:
if I set (in the problem.hh file):
bool shouldWriteRestartFile() const
{
return *true*;
}
then I get, for example:
*tutorial_sequential_time=3000_rank=0.drs*
Is it correct that, in order
Good morning,
I'm facing two types of issues:
1) On DuMux 2.12, trying to run any test which calls set_singular_limit
results in the following error.
/home/lorenzo/DUMUX/dumux/dumux/material/constraintsolvers/compositionfromfugacities.hh:110:59:
error: ‘set_singular_limit’ is not a member
Hello,
which Dune version are you using? DuMux 2.12 is only compatible with
the 2.6 release not master.
I'm using DuMux 2.12 and using compatible branch.
If I use Dumux 3.0.0 I do not encouter such problem, but I would like to
fix it on 2.12 also.
I assume you follow the instructions in
Good morning,
I have one general question:
is it possible to run 1 dimensional test case, for example, buckely
leverett?
Kind regards,
Lorenzo
On 19.12.2018 13:26, Timo Koch wrote:
Hi Lorenzo,
does it get better if you lower the CFL factor?
Timo
On 19.12.18 11:10, lc wrote:
Hello
Hello,
yes, I saw the example in the dumux-lecture.
Is it the dimension induced by the Grid definition in the input file:
Cells 100 1 ?
And consequently the system of equations solved is 1d?
Thank you,
Lorenzo
On 21.12.2018 11:52, Timo Koch wrote:
Am 21.12.2018 um 08:39 schrieb lc
Dear all,
I catch the occasion to greeting you the best wishes for Christmas and
holidays!
I understand this may not be the best moment for questions, but maybe
not for everybody is Xmas ;)
Since I got interesting results with DuMux 2.12 in the study of viscous
fingers, I was asked to
Good morning,
I have a general question.
What is the highest formal order of accuracy in space and time in DuMux?
As I read, using central FV and Box model and FE and explicit/implicit
Euler the maximum should be 2. Is it correct?
Which is the combination which allows to get the most
Hello,
On 02.11.2018 10:45, Dennis Gläser wrote:
per default, the formulation for the 2p model is pw-sn. That means
your primary variables are the water pressure and the non-wetting
phase saturation (in your case oil I assume). Therefore,
Indices::swIdx does not exists as it is not part of
ly you can open the files and read
through the code.
I do not understand your question about the alugrid.cfg file... I
never did anything in particular with it.
Best regards,
Beatrix
Am 05.12.18 um 10:56 vorm. schrieb lc:
Hello,
thanks for the previous answer, it was helpful.
Hello,
I have the following questions:
1) Is it possible to run IMPES algorithm in sequential (not implicit)
mode without considering capillarity forces (pressure)?
2) Is it possible to do some sort of grid adaptation, for example
locally refine the (unstructrured) mesh where the water
to realize and personally I would try to avoid them and
restrict myself to Neumann and/or Dirichlet BCs where you actually
know what exactly happens.
Best wishes,
Dennis
*Von:* Dumux im Auftrag von
lc
*Gesendet:* Montag
Hello,
As I understand from previous emails, for the porous 2d implicit box
model, the primary variables are wetting phase pressure and non-wetting
phase saturation.
If so, when imposing BCs and ICs, it can be done on non-primary
variables? Or it is always better to define the problem in
enzo,
yes, in the new version, mixed boundary conditions (e.g., setting
Dirichlet for p_w and Neumann for q_n) are only possible for the box
scheme. Also, the outflow BC was removed, so you have to implemented
that yourself, e.g., using an Neumann condition.
Best wishes
Kilian
On 11/23/2018 10:00 AM
Hello,
decreasing the CFL worked.
Thanks,
Lorenzo
On 19.12.2018 13:26, Timo Koch wrote:
Hi Lorenzo,
does it get better if you lower the CFL factor?
Timo
On 19.12.18 11:10, lc wrote:
Hello,
I observed an unexpected behaviour for which I'd like to ask your
help. In the attached figure
Hello Dumux team,
We are, finally, in the post-processing phase, aiming at publishing soon
... but we observed a scaring issue concerning the verification of the
BC at the inlet.
I'm using Dumux 2.12 with the 2p sequential model for a rectangular
domain discretized with a cartesian
ue with Dumux.
Best wishes
Timo
On 27.03.19 14:53, lc wrote:
Hello Timo,
I enclose 2 solutions, t = 0 and t = 600 days.
Indeed, I have a saturation profile as IC as you can see at t = 0.h
Best regards,
Lorenzo
On 27.03.2019 16:41, Timo Koch wrote:
Dear Lorenzo,
generally for cell-cent
Hello Timo,
Does it also occur when your grid is fine from the beginning and you
don't have hanging nodes?
Yes, I tried different variants, strctured, adapted, unstructured and
the behaviour is the same.
What are your criterions for Newton convergence?
At the beginning, I just copied
Hello Timo,
Does it also occur when your grid is fine from the beginning and you
don't have hanging nodes?
Yes, I tried different variants, strctured, adapted, unstructured and
the behaviour is the same.
What are your criterions for Newton convergence?
At the beginning, I just copied
Hello,
I'm still fighting to fix the parallel execution.
Now, I re-install mpich/openmpi and configured Dumux 2.12 (and also 3.0.0).
I'm running: porousmediumflow/2p/sequential/test_impes testcase with 1
and 4 core which should be ok for the features of my machine.
This time seems that the
y parallel solver is currently the AMGBackend. There
is also a test called "test_impeswithamg". Please try that one in
parallel.
Timo
On 06.02.19 11:58, lc wrote:
Hello,
I'm still fighting to fix the parallel execution.
Now, I re-install mpich/openmpi and configured Dumux 2.12 (and
Here is the log,
Thank you,
Lorenzo
On 30.01.2019 14:53, Bernd Flemisch wrote:
Are you sure that dunecontrol finds MPI? Can you post the output of
dunecontrol?
Bernd
On 01/30/2019 12:46 PM, lc wrote:
Good morning,
I updated the virtual machine as you may see from the lscpu_new file
___
*Von:* Dumux im Auftrag von
lc
*Gesendet:* Donnerstag, 24. Januar 2019 10:44:20
*An:* dumux@listserv.uni-stuttgart.de
*Betreff:* [DuMuX] about grid adaptation
Good morning,
is it possible to obtain
Good morning,
but so, how can I impose the flux of the non-wetting phase (oil) to be
constant, dsn/dn = 0, at the right boundary?
Kind regards,
Lorenzo
On 31.01.2019 0:43, Flemisch, Bernd wrote:
1) It seems that you are setting no-flow for the saturation on the
right boundary. I'm not
Hello,
yes, sorry, so what I actually want to impose is that the normal
derivative of saturation at the right boundary is = 0 and as I
understand since, sw and sn are related, it doesn't matter which
saturation.
Regards,
Lorenzo
On 01.02.2019 12:12, Timo Koch wrote:
It’s unclear what
Good afternoon,
after having searched on internet, I didn't find the info which I need.
Could you please explain the details of the implementation of the
Outflow (even if it has been removed in the new release) and Neuman BC
in DuMux and in what they differ, please? In particular, in the
Good evening,
I'd like to ask about the chemical reactions handling capability in
DuMux. Is it able to handle the kind of processes which occur typically
with surfactants and reactive polimers?
Best regards,
Lorenzo
___
Dumux mailing list
Good morning,
I'd like to ask whether it would be possible to simulate a typical
*quarter five spot problem*?
How is it managed such configuration given that we would like to
consider *cylindrical wells*?
Kind regards,
Lorenzo Campoli
--
___
*Il
Hello,
I'd like to ask your opinion about a behaviour of the wetting phase
saturation at inlet.
Here I simulate an horizontal well [150 x 100] m, without gravity but
with capillarity (Pc) with 2p (pwsw formulation) on Dumux 2.12.
Neuman on pressure and Dirichlet on saturation is imposed at
Dear Dumux community,
I'm simulating a 2d horizontal well with water/oil with 2p model.
At the inlet (left) boundary I impose Neuman condition on pressure,
dp/dx = const. and so pressure will not stay constant.
Is it possible to impose a "variable" pressure in time (as explained in
the
Hello,
can you explain how can I strongly impose the Dirichet condition in the
first cell (or the closest one) on the inlet (left ) boundary in order
to achieve what you said, please?
I need to use 2p, sequential algorithm.
Kind regards,
Lorenzo
On 29.03.2019 16:31, Timo Koch wrote:
I
Dear community,
I'd like to report you some question/observation:
1) The links in https://dumux.org/projects don't work.
2) I'm going to deal with surfactants, which in my case means, 2p3c or
3p3c: water and oil phases, and polymer/surfactant agent as third
component.
Now, as I understant,
Dear Dumux developers,
how can I assign a random permeability value for each cell?
I asa thinking something like this but I need to loop over all cell ...
Scalar RandPerm = 0.;
for (int i = 0; i < dim; i++) {
RandPerm = rand() % 5 + 98;
RandPerm *=
Hello dear Dumux community,
I'd like to ask 2 question:
1) how it is discretized an equation like: du/dt = - du/dx in Dumux. Of
course, I'm not asking for the all course of numerical discretization of
PDF, but the specific solutions which are implemented in Dumux since I
noticed that even if
Dear Dumux community,
I'd like to ask how can I simulate a typical polymer slug injection
scenario in an horizontal well, WAP (water alternate polymer).
In particular, for example, 100 days of water injection, followed by
other 100 days of polymer (not important specifically what it is),
Hello Dear Dumux developers,
I'd like to ask:
1) Are the Langmuir and Todd-Longstaff models for the description of
polymer adsorption implemented?
If yes, is there any example?
If no, some hint about where to add.
I only found a params.KdNAPL.
Secondly, is there an example of the
Hello,
is it possible to simulate adsorption/desorption in Dumux in 1d and 2d,
with some model?
If yes, is there any example?
Thank you,
Kind regards,
Lorenzo
--
___
*Il tuo 5 diventa 1000*
Fai
crescere la tua università
Dona il 5 per mille alla
al schemes don't get very much attention right now from the
core developers.
Sorry I couldn't help more, open-source research software is always a
work in progress,
Timo
On 15.07.19 13:15, lc wrote:
Hello,
can you explain how can I strongly impose the Dirichet condition in
the first cell (or
Hello,
thank you.
On 28.01.2020 15:15, Kilian Weishaupt wrote:
Do you consider a water/polymer mixture at a fixed concentration?
Well, polymer concentration is "the big question". Indeed, the task
could be to find the optimal value or variation curve to optimize oil
output. So, now another
Good morning,
is it possible to retrieve somewhere the testcase reported in Sect. 5.2
of https://doi.org/10.1016/S0309-1708(02)00014-3 ?
Kind regards,
Lorenzo
--
Le informazioni
contenute in questo messaggio di posta elettronica
Hello,
2) You have to implement your custom "FluidSystem" class that can,
e.g., inherits from whatever fluid system you used before and overload
the viscosity function. You can use the surfactant concentration
passed in through the fluidState from the VolumeVariables class.
I'm trying to do
*Scalar SolidDensity = getParam("Solid.Density");*
in the *computeStorage *function, which is quite costly.
So you should define it as
*static const Scalar SolidDensity = getParam("Solid.Density");*
Regards,
Martin
Am 19.02.20 um 14:09 schrieb lc:
Hello,
On 19.02.2020 14
Dear Dumux community,
speaking about adsorption in Dumux, I'd like to ask some clarification.
With respect to the following answer and considering a 1p2c model, in
equilibrium.
I need to implement an adsorption term such as d(c+A(c))/dt ... where c
is the generic concentration (mass or
Good morning,
I did as suggested and I observe a remarkable slow down down of the
execution, keeping everything the same.
Is this due to the fact that we add a new term in the time derivative
and the system should be re-assembled?
Is it an expected behaviour?
I checked that this happens
:* Dumux im Auftrag von
lc
*Gesendet:* Mittwoch, 19. Februar 2020 11:34:27
*An:* dumux@listserv.uni-stuttgart.de; Koch, Timo
*Betreff:* Re: [DuMuX] adsorption/desorption term
Good morning,
I did as suggested and I observe a remarkable slow down down of the
execution, keeping everything th
55 matches
Mail list logo