Re: [DuMuX] BC apparently not respected

2019-07-16 Thread Timo Koch
On 16.07.19 12:13, lc wrote: Hello, Maybe this helps to understand what I’m talking about. https://en.m.wikipedia.org/wiki/Minimal_working_example In particular, “The important feature of a minimal working example is that it is as small and as simple as possible, such that it is just

Re: [DuMuX] BC apparently not respected

2019-07-16 Thread Timo Koch
Am 16.07.2019 um 10:24 schrieb lc mailto:lorenzo.camp...@uniroma1.it>>: Hello On 16.07.2019 10:22, Timo Koch wrote: please post a minimal working example to reproduce the behaviour. Also, we need to know what you are computing, what is the outcome, what is the expected outcome and why you

Re: [DuMuX] BC apparently not respected

2019-07-16 Thread Timo Koch
> Am 16.07.2019 um 08:55 schrieb lc : > > Hello Timo, > > thanks for your answer. > > Since, I'm working with Dumux 2.12 I cannot say anything about 1) or 2). > > So, for the moment, I attach al the material necessary to replicate the > behaviour that I get on Dumux 2.12, at least to

Re: [DuMuX] BC apparently not respected

2019-07-16 Thread lc
Hello Timo, thanks for your answer. Since, I'm working with Dumux 2.12 I cannot say anything about 1) or 2). So, for the moment, I attach al the material necessary to replicate the behaviour that I get on Dumux 2.12, at least to understand if it is a bug, my implementation mistake or should

Re: [DuMuX] BC apparently not respected

2019-07-15 Thread Timo Koch
Hi Lorenzo, I think the status is as follows (Dumux 3.0): 1) For the box scheme Dirichlet boundaries are always strongly enforced so nothing special to do here. 2) For cell-centered schemes there is a new feature called "internal Dirichlet constraints" that allows you to enforce Dirichlet

Re: [DuMuX] BC apparently not respected

2019-07-15 Thread lc
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

Re: [DuMuX] BC apparently not respected

2019-03-29 Thread Timo Koch
On 29.03.19 13:46, lc wrote: 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?

Re: [DuMuX] BC apparently not respected

2019-03-29 Thread lc
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

Re: [DuMuX] BC apparently not respected

2019-03-29 Thread lc
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

Re: [DuMuX] BC apparently not respected

2019-03-29 Thread lc
Good afternoon Beatrix, I tried to apply your suggestion, imposing setAllDirichlet at the inlet (left) and  everywhere else setAllNeumann, as was also done in the tutorial handbook. Unfortunately, the behaviour at the inlet regarding the "degradation" of the BC is still present even if I

Re: [DuMuX] BC apparently not respected

2019-03-27 Thread Beatrix Becker
Dear Lorenzo, in my opinion it makes no sense to set a Dirichlet and a Neumann condition at the same time. Either you know the flux over your boundary or you know the values of your primary variables on the boundary. If you change your code to setAllNeumann on the left boundary of your

Re: [DuMuX] BC apparently not respected

2019-03-27 Thread Timo Koch
Hi Lorenzo, you still have large capillary pressure gradients at the left boundary at t=600 thus the saturation changes. You also have them in the initial conditions, so they probably never go away. Why this is so, depends on your setup. Does it also occur when your grid is fine from the

Re: [DuMuX] BC apparently not respected

2019-03-27 Thread Timo Koch
Dear Lorenzo, generally for cell-centered discretizations in Dumux, Dirichlet boundary conditions are weakly enforced on the boundary. Your first degree of freedom is half a cell away from the boundary, so it doesn't necessarily have to be the same as the Dirichlet value. In order to see if

[DuMuX] BC apparently not respected

2019-03-27 Thread lc
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