Re: [DuMuX] boundary conditions

Hi Elena, 1.I need to set the type of boundary conditions (Neumann or Dirichlet) depending on the fluid pressure at the boundarys , how can I do it? Use a solution dependent Neumann boundary condition: solDependentNeumann() function. See the

Re: [DuMuX] Boundary conditions in decoupled problem

Hi Martin, ok, so i corrected the condition, so now, my code is: void boundaryTypesAtPos(BoundaryTypes , const GlobalPosition& globalPos) const /*@\label{tutorial-decoupled:bctype}@*/ { if ((globalPos[0] > 580- eps_) && (globalPos[1] > 580- eps_) ) {

Re: [DuMuX] Boundary conditions in decoupled problem

no, sorry, there isn't the good files, sorry. the input file is [TimeManager] TEnd = 10 # duration of the simulation [s] DtInitial = 10 # initial time step size [s] [Grid] UpperRightX = 600 # x-coordinate of the upper-right corner of the grid [m] UpperRightY = 600 # y-coordinate of the

Re: [DuMuX] Boundary conditions in decoupled problem

i show that not all the files are joint to my precedent message, so the intersec permeability is: for (int i = 0; i < dim; i++) K_[i][i] = 1e-7; and i use the brook corey law: // residual saturations materialLawParams_.setSwr(0); materialLawParams_.setSnr(0);

Re: [DuMuX] Boundary conditions in decoupled problem

I joint to you the file .input, he tutorialspatialparams_cdecoupled and the file .hh to show all parameters that i use. Thak's for the help. Best regards. 2015-10-21 21:10 GMT+02:00 Martin : > That probably means that the matrix is singular and therefore

Re: [DuMuX] Boundary conditions in decoupled problem

In your input file you have set: 
[Grid] UpperRightX = 300 # x-coordinate of the upper-right corner of the grid [m] UpperRightY = 60 # y-coordinate of the upper-right corner of the grid [m] NumberOfCellsX = 100 # x-resolution of the grid NumberOfCellsY = 1 # y-resolution of the grid 
 that means,

Re: [DuMuX] Boundary conditions in decoupled problem

So i try with TEND = 1000 s, and i have nothing in paraview, there is no flow, you can see this in the picture joint in last time simulation. I don't undertand anything. Have you any idea? Best regards. 2015-10-21 23:37 GMT+02:00 Ait Mahiout Latifa : > I tried with

Re: [DuMuX] Boundary conditions in decoupled problem

I tried with 100 cells but it will run for 25 minutes and has not finished yet, however TEND=1 s only. So what's the problem in your opinion? Best regards. 2015-10-21 21:52 GMT+02:00 Ait Mahiout Latifa : > i try with 100 cells and i will tell you if it gives a

[DuMuX] boundary conditions in coupled problem

Hi, is it possible to implement boundary like this: for example, we have and domain 600m*600m, and we have impos the boundary contitions: if x =20m and y=20m: Sn=0.0 and q_w . n = -e-8, and q_n.n = 0 The method used in decoupled problem to implement in the same times Dirichlet and Neumann

Re: [DuMuX] boundary conditions in coupled problem

Hi, i try to change the boundary conditions in tutorial_coupled problem. The domain is 600m*600m, and the boundary conditions are: if x=580m and y=580m: S_n=1 and P_w=350 bar if x=20m and y=20m:S_n=0, q_w . n = -1.e-8 and q_n.n=0 and nos flux in the rest of domain. So, i try to do like in

Re: [DuMuX] Boundary conditions

Alexander Kissinger Gesendet: Freitag, 31. Juli 2015 16:21 An: DuMuX User Forum Betreff: Re: [DuMuX] Boundary conditions Dear Georg, If I understand correctly this means that whatever is on an outflow boundary is allowed to flow out or into the system. So if we assume pure fickian diffusion

Re: [DuMuX] Boundary conditions

Forum Betreff: Re: [DuMuX] Boundary conditions Dear Georg, I am not quite sure what causes your problem. Could you send the implementation of your boundary conditions, i.e. your boudaryTypes() and your dirichlet() function from your problem? You could also try the following: 1. Just set

Re: [DuMuX] Boundary conditions

; dumux@listserv.uni-stuttgart.de *Betreff:* Re: AW: AW: [DuMuX] Boundary conditions but then grad pc = dpc/dsw*grad sw so setting grad pc*n = 0 should be same as grad s*n =0 if I'm not wrong? Best regards, Timo Am 12.08.2015 um 08:26 schrieb georg.fut...@dlr.de mailto:georg.fut...@dlr.de

Re: [DuMuX] Boundary conditions

An: Futter, Georg; dumux@listserv.uni-stuttgart.de Betreff: Re: AW: [DuMuX] Boundary conditions Hi Georg, If you have grad(S) somewhere in your equation yes. Otherwise no. If you want to have a boundary condition like grad(S)*n = 0, calculate the residual (the whole flux term, i.e. everything

Re: [DuMuX] Boundary conditions

Dear Georg, when you set Dirichlet boundary conditions the equation gets replaced by the Dirichlet boundary condition. I think that's what's clear. When you set Neumann boundary conditions in Dumux you set a value for the entire flux term (advective + diffusive part). So when you set a

Re: [DuMuX] Boundary conditions

? Best regards Georg Von: Timo Koch [mailto:timo.k...@iws.uni-stuttgart.de] Gesendet: Dienstag, 11. August 2015 16:03 An: DuMuX User Forum; Futter, Georg Betreff: Re: [DuMuX] Boundary conditions Dear Georg, when you set Dirichlet boundary conditions the equation gets replaced

Re: [DuMuX] Boundary conditions

von Alexander Kissinger Gesendet: Donnerstag, 30. Juli 2015 08:53 An: DuMuX User Forum Betreff: Re: [DuMuX] Boundary conditions Dear Georg, I am not sure if I got you right, you want to have a fixed Saturation (Dirichlet) for one phase and inject another phase? In Dumux you can choose

Re: [DuMuX] Boundary conditions

*Von:*Dumux [mailto:dumux-boun...@listserv.uni-stuttgart.de] *Im Auftrag von *Alexander Kissinger *Gesendet:* Donnerstag, 30. Juli 2015 08:53 *An:* DuMuX User Forum *Betreff:* Re: [DuMuX] Boundary conditions Dear Georg, I am not sure if I got you right, you want to have a fixed Saturation

Re: [DuMuX] Boundary conditions

Dear Georg, I am not sure if I got you right, you want to have a fixed Saturation (Dirichlet) for one phase and inject another phase? In Dumux you can choose the equation that should be replaced by the Dirichlet condition with the call: setDirichlet(int pvIdx, int eqIdx) The equation you