Re: [DuMux] Calculate solution at a given position

Dear Etienne, yes there is. But it’s not direct. You can use evalSolution (from dumux/discretization/evalsolution.hh) but for that you need an element. To find the element that your designated point is in you can use the boundingBoxTree of the gridView for an efficient query. In

Re: [DuMux] forking dumux project on git.iws.uni-stuttgart

Dear Mona, I think due to more restrictive setting on the GitLab server (due to spam problems), new users have a project limit of 0 which first needs to be increased by an admin. I forwarded your request and you should have access now. We are working on improving this situation. Best wishes

Re: [DuMux] Questions about residual saturation

Dear Nicolas,the residual saturation relates to the mobility of the phase (in terms of transport of the phase). But the “immobile” water can still evaporate, so in the compositional models in dumux by phase change the liquid saturation can still go to 0 when the gas phase is mobile. (Only then you

Re: [DuMux] Problem with linear solver residual computation

Hi Leo,did you try parallel vs sequential and is there a difference? This would help narrowing down where a potential problem could come from?BestTimoOn 26 Jun 2023, at 11:13, leopold.stad...@baw.de wrote: Hi Bernd,   Thank you fro your advice. My code does not

Re: [DuMux] richardsnc for more than two components

Dear Mona, yes this is exactly how you would do it. The model deduces the number of components from the fluid system. It’s actually not too difficult to add more components to liquid2c, so you can use that as a base. Set numComponents to 3 (or more). Add a third (fourth…) component whereever

[DuMux] Dune User Meeting 2023 Dresden

Hi Dumux users, the next Dune User Meeting 2023 will take place in Dresden from 18-19 Sep 2023. This is a good chance to present your research with DuMux (which is as a Dune module part of the Dune community)! More information can be found at

Re: [DuMux] [DuMuX installation] Unable to install DuMuX on Ubuntu

Dear Timothé, welcome to the Dumux mailing list and cool that you are considering Dumux. From the logs it appears that something is wrong with your local installation of openmpi. You might have changed your openmpi installation inbetween two installs? I would try to * reinstall openmpi using

Re: [DuMux] Gas Injection - Time

Hi Mohammad, This seems like a C++ problem not a Dumux issue. What have you tried so far? What does it mean “it will abort”? "Method 2” has the wrong logic. “Method 1” looks correct (apart from underscore missing for injectionEnd_). Your attached file neither implements 1 nor 2. Best Timo

Re: [DuMux] Printing Flow Rate or Flux in VTK

Dear Mohammad, please reply to your own question with a follow-up instead of posting the same question twice. Have you looked into the mail archive (https://www.mail-archive.com/dumux@listserv.uni-stuttgart.de/)? A similar question was just answered a few weeks ago. Maybe this solves your

Re: [DuMux] Alugrid - error

ds, Kenza Le sam. 5 nov. 2022 à 12:08, Timo Koch mailto:timok...@math.uio.no>> a écrit : 5. nov. 2022 kl. 12:32 skrev kenza bouznari mailto:bouznari.ke...@gmail.com>>:  Thank you for pointing that out. You are right, I had to change it, but now I get a new error that says : WA

Re: [DuMux] Output face fluxes in VTU

Hi Rafael, this is not implemented in any convenient way. Probably because it’s not clear how to _visualise_ normal face fluxes in a useful way. Vtk doesn’t support face quantities. You can write out the skeleton (only the faces of the grid) if you need something for debugging or

Re: [DuMux] Alugrid - error

egards, Kenza Le ven. 4 nov. 2022 à 19:38, Timo Koch mailto:timok...@math.uio.no>> a écrit : Hi Kenza, the exception tells you exactly what’s wrong. You specified “parameters 1” so it expects one parameter per simplex. But your simplices don’t have parameters associated. In the first line (as

Re: [DuMux] Negative Values -Gas Saturation

On 4 Nov 2022, at 20:54, Mohammad Hodroj (Student) mailto:mn...@mail.aub.edu>> wrote: Dear DuMux Community, I hope this email finds you well. I am writing to kindly ask for your support in the following: a)I am using porous medium flow b) injection problem Issue: I am getting negative gas

Re: [DuMux] Alugrid - error

t file, and it did work. So I believe that the problem is in the *.dgf file and how it's read by the code. Best regards, Kenza Le ven. 4 nov. 2022 à 15:33, Timo Koch mailto:timok...@math.uio.no>> a écrit : Hi Kenza, did you change the AluGrid type in the Grid property from 2D to 3D? This i

Re: [DuMux] Alugrid - error

Hi Kenza, did you change the AluGrid type in the Grid property from 2D to 3D? This is a compile time parameter. Best, Timo 4. nov. 2022 kl. 19:51 skrev kenza bouznari :  Dear DuMux users, I have installed ALUGrid and successfully built DuMux and run test_2p2c_injection_box. I also made all

Re: [DuMux] Computing Boundary Fluxes

Dear Rafael, welcome to the Dumux mailing list! > On 14 Oct 2022, at 23:38, Rafael March wrote: > > Hello, > > I hope this email finds you well. > > I'm trying to use DuMux to perform single and multiphase flow-based upscaling > (Darcy scale). > > This involves running two-phase flow

Re: [DuMux] Printing Permeability in VTK Files

time step? Thanks for your help. Best, Mohammad Get Outlook for Android<https://aka.ms/AAb9ysg> ____ From: Timo Koch mailto:timok...@math.uio.no>> Sent: Wednesday, October 12, 2022, 2:31 PM To: DuMuX User Forum mailto:dumux@listserv.uni-stuttgart.de>>

Re: [DuMux] Printing Permeability in VTK Files

Kindly, find attached the requested files (spatial, problem, and main). Looking forward to your help. Best Regards, Mohammad From: Timo Koch mailto:timok...@math.uio.no>> Sent: Wednesday, October 12, 2022 1:58 PM To: DuMuX User Forum mailto:dumux@listserv.uni-

Re: [DuMux] Printing Permeability in VTK Files

Dear Mohammad, it’s bit difficult to read code snippets as a listing because it matter where things are in the file. If below doesn’t help it would be great if you could attach the spatial params header instead. In (6.) you have "return permeability;” where permeability does not have a

Re: [DuMux] Anisotropic heat conductivity tensor in heat conduction simulation

Hi Helena, you saw correctly that this is not implemented so far. However, this is a relatively straight-forward modification of the EnergyVolumeVariables (https://git.iws.uni-stuttgart.de/dumux-repositories/dumux/-/blob/master/dumux/porousmediumflow/nonisothermal/volumevariables.hh) since

Re: [DuMux] Question - Domain Shape

Dear Mohammad, alternatively, if you have unstrctured grids, use Gmsh to mesh your domain and then specify the Gmsh file [Grid] File = your_grid.msh in the input file. You need a an unstructured grid manager like UGGrid or ALUGrid, and currently only Gmsh version 2 format is supported. (use

Re: [DuMux] GridError: Unknown priority 4

Hi Dmitry, you can follow https://gitlab.dune-project.org/core/dune-grid/-/merge_requests/597 where the problem and hopefully solution will be tracked. Best Timo On 4 Sep 2022, at 23:45, Timo Koch mailto:timok...@math.uio.no>> wrote: Hi Dmitry, thanks for reporting this. I also n

Re: [DuMux] GridError: Unknown priority 4

Hi Dmitry, thanks for reporting this. I also noticed the same issue recently. The problem is that the dune-uggrid gridview is not thread-safe (yet). It works fine in many situation but in combination with MPI there seems to be an unfixed issue. dune-uggrid doesn’t actually claim to be

Re: [DuMux] Problems with DUNE_THROW

Hi Kai, are you sure this is within the try & catch block of the Newton? It be very surprised if the try/catch does not catch the exception. You can use a debugger to find out where this exactly happens. The time step reduction mechanism is only activated if you call newton.solve(x, timeLoop)

Re: [DuMux] DuMuX Simulation-Permeabilityupscaling example

I can see something similar in the test 2p example. > Also, I figured that all the values obtained are in SI units. Is that correct? yes SI units. Best wishes Timo > > Regards, > Rucha > > Von: Timo Koch > Gesendet: Mittwoch, 10. August 2022 15:35 > An: DuMuX Us

Re: [DuMux] DuMuX Simulation-Permeabilityupscaling example

Dear Rucha, the permeability upscaling example solves a stationary equation which has no concept of time. There it is unclear what you are trying to achieve. You also shouldn’t “see” any time steps. Maybe you are referring to two vtu files? Then it’s probably one containing the initial guess

Re: [DuMux] Error while running PoreNetwork Model Example-Weishaupt2019a

egards, Rucha Von: Timo Koch mailto:timok...@math.uio.no>> Gesendet: Dienstag, 19. Juli 2022 13:13 An: Kunte, Rucha mailto:rucha.ku...@dlr.de>> Betreff: Re: [DuMux] Error while running PoreNetwork Model Example-Weishaupt2019a On 19 Jul 2022, at 13:02, rucha.ku...@dlr.de<mailto:ruch

Re: [DuMux] Error while running PoreNetwork Model Example-Weishaupt2019a

to the CMAKE_FLAGS in the *.opts file helps. Remember to also use the modified opts file in the last installation step (dunecontrol). Best wishes Timo On 19 Jul 2022, at 10:31, Timo Koch mailto:timok...@math.uio.no>> wrote: Dear Rucha, from the attached log it looks like you have dune-uggrid p

Re: [DuMux] Doubts for Using DUMUX-Installation and Importing geometries

On 5 Jul 2022, at 12:49, mailto:rucha.ku...@dlr.de>> mailto:rucha.ku...@dlr.de>> wrote: Hello, My name is Rucha Kunte and I am writing my Master Thesis at DLR-Cologne. The objective is to carry out a permeability analysis of a porous structure for which we plan on using Dumux. I have a

Re: [DuMux] Singular matrix and MPI

https://git.iws.uni-stuttgart.de/dumux-repositories/dumux/-/issues/1171 Unfortunately, I can not prepare a proper test at present. Best regards, Dmitry On 6/30/22 13:55, Timo Koch wrote: Hi Dmitry, looks like you are running into a deadlock. Looks like a bug in the exception handler of the Newton (I

Re: [DuMux] Turn parallel on and off

t Timo Thanks a lot Ana Em sex., 17 de jun. de 2022 às 09:00, Timo Koch mailto:timok...@math.uio.no>> escreveu: Dear Ana, you’d have to explain more. What you are trying to achieve is not clear to me. The Newton itself doesn’t really care about parallelism. The assembler assembles on

Re: [DuMux] Singular matrix and MPI

Hi Dmitry, looks like you are running into a deadlock. Looks like a bug in the exception handler of the Newton (I assume https://git.iws.uni-stuttgart.de/dumux-repositories/dumux/-/blob/master/dumux/nonlinear/newtonsolver.hh#L548 is wrong because the exception cannot be assumed to be thrown on

Re: [DuMux] 1D grid and MPI

Dear Dmitry, this is a known current bug (https://git.iws.uni-stuttgart.de/dumux-repositories/dumux/-/issues/864 ) There are some ideas for a fix but due to limited developer resources we are currently focussing on a

Re: [DuMux] DuMux - Help

> On 19 Jun 2022, at 17:52, Mohammad Hodroj (Alumni) wrote: > > Dear DuMux Team, > > I hope my email finds you well. > > My name is Mohammad Hodroj and I am graduate student at the American > University of Beirut. > > I am writing to kindly ask for your help in using DuMux to run my

Re: [DuMux] Turn parallel on and off

Dear Ana, you’d have to explain more. What you are trying to achieve is not clear to me. The Newton itself doesn’t really care about parallelism. The assembler assembles on the gridview it is given (usually the gridview of one process (maybe you read your grid on all processes and don’t

Re: [DuMux] dune-subgrid_FOUND failed

Hi, two things could solve this depending on what the problem is. (1) Remove the build folders (or at least CMakeCache and CMakeFiles) and run dune control again. (2) Make sure you havea compatible Dumux version. The newest Dumux version (3.5) is not compatible with Dune 2.7. You need at

Re: [DuMux] hydrogen injection setup

g (CPE), University of Kansas Email: vuongvanp...@ku.edu<mailto:vuongvanp...@ku.edu> Phone: +1-(785)-979-2664 From: Timo Koch mailto:timok...@math.uio.no>> Sent: 12 Tháng Năm 2022 11:44 SA To: DuMuX User Forum mailto:dumux@listserv.uni-stuttgart.de>> Cc: Pham, Vuong Van ma

Re: [DuMux] hydrogen injection setup

rds, Vuong Van Pham Graduate Research Assistant (GRA) Department of Chemical and Petroleum Engineering (CPE), University of Kansas Email: vuongvanp...@ku.edu<mailto:vuongvanp...@ku.edu> Phone: +1-(785)-979-2664 From: Timo Koch mailto:timok...@math.uio.no>> Sent: 12 Tháng Năm 2022 1

Re: [DuMux] hydrogen injection setup

ku.edu<mailto:vuongvanp...@ku.edu> Phone: +1-(785)-979-2664 From: Timo Koch mailto:koch_t...@hotmail.com>> Sent: 28 Tháng Ba 2022 6:12 SA To: DuMuX User Forum mailto:dumux@listserv.uni-stuttgart.de>> Cc: Pham, Vuong Van mailto:vuongvanp...@ku.edu>> Subject: [DuMux] hydr

Re: [DuMux] info about the software

oad info about State and Private Organisations which have used/are > using DUMUX as research purposes. Also about User Meetings, which is > wonderful. > > > Liebe Grüße, > Sevinj > > > > > > > On 2022-04-20 21:22, Timo Koch wrote: >> Dear Sevinj,

Re: [DuMux] info about the software

Dear Sevinj, welcome to the Dumux mailing list and happy to hear that you are considering Dumux. I’ll try to answer but it would be easier with more detail what you are trying to achieve. 1. Parallelization is possible for many models, both distributed memory parallelism with MPI and shared

Re: [DuMux] Reentrant DuMuX

> On 15. Apr 2022, at 13:07, Dmitry Pavlov wrote: > > Hello, > > I am considering a possibility to use a DuMuX program not as an executable, > but as a shared library. In addition, this library is going to pass the > results to caller via memory, rather than storing .vtu files on disk. Also,

Re: [DuMux] Inheritance of properties

Hi Dmitry, the problem is the line > template > struct TestProperty { > static constexpr bool value = false; }; here you define a default for _all_ TypeTags that don’t have this property specialised. Therefore the property is already defined for TestTag3 and it

[DuMux] hydrogen injection setup

eum Engineering (CPE), University of Kansas > Email: vuongvanp...@ku.edu <mailto:vuongvanp...@ku.edu> > Phone: +1-(785)-979-2664 > > > > From: Timo Koch > Sent: 08 Tháng Ba 2022 3:52 CH > To: DuMuX User Forum > Cc: Pham, Vuong Van > Subject: Re: [DuMux

Re: [DuMux] Installation issues

Dear Vuong, please subscribe to the dumux mailing list (https://listserv.uni-stuttgart.de/mailman/listinfo/dumux), so you can see the answers to your post. Two people already posted answers to your question. Best Timo > On 7. Mar 2022, at 04:56, Pham, Vuong Van wrote: > > To someone who may

Re: [DuMux] Model coupling in DUMUX

Dear Etienne, unfortunately there is currently no example of such a simulation setup with two-phase flow. But there is a simple example for one-phase flow + transport example: https://git.iws.uni-stuttgart.de/dumux-repositories/dumux/-/tree/master/examples/1ptracer Consider two things: (1)

Re: [DuMux] CO2 geological storage

Dear Etienne, we are trying to find out what part of the implementation causes the difference. Could you try using the 3.4 Darcy's law implementation changing line 193 in darcyslaw.hh to const auto& tij = std::abs(fluxVarsCache.advectionTij()); That would be very helpful! One possible reason

Re: [DuMux] CO2 geological storage

Hi Sebastian and Etienne, can any of you two open an issue with the problem description and proposed solution @ https://git.iws.uni-stuttgart.de/dumux-repositories/dumux/-/issues? This looks like something we should definitely look into before the next release. The new version probably has a

Re: [DuMux] DuMuX Installation Problem

Hi, just for documentation purposes: this question was answered on the mailing list here: https://www.mail-archive.com/dumux@listserv.uni-stuttgart.de/msg02594.html Timo > On 15. Feb 2022, at 17:07, Flemisch, Bernd > wrote: > > Hi Farzaneh, > > you can delete the "-march=native" from

Re: [DuMux] Contribution to DuMuX - Suggestion to improve the course-documentation

Dear Axel, welcome to the mailing list. Feedback regarding the beginner material is always welcome. Also improvements to the beginner material and docs are welcome of course. Please consider two things: 1) We want to avoid duplicate resources to reduce maintenance effort regarding docs. The

Re: [DuMux] parallel msh grid

al > mesh format or DGF. > > Kind regards, > Gergely > > Am 12.01.2022 um 17:47 schrieb Timo Koch: >> Dear Gergely, >> >> if the mesh fits into memory, you don’t need to do anything special for >> parallel runs since the default grid manager will call g

Re: [DuMux] parallel msh grid

Dear Gergely, if the mesh fits into memory, you don’t need to do anything special for parallel runs since the default grid manager will call grid.loadBalance() which distributes the grid. The grid data helper also automatically distributes user data like element markers, in case you are using

Re: [DuMux] How to prescribe an uniform background hydraulic head gradients in DuMuX?

> On 15. Dec 2021, at 18:16, Etienne Ahusborde > wrote: > > Dear DuMuX, > > With my colleagues, we are interested to perform a test case of injection of > CO2 into a reservoir with uniform background hydraulic head gradients. > > > > The configuration is represented in the figure above and

Re: [DuMux] Help with Target CPU error

Hi Diya, welcome to the Dumux mailing list. I haven’t tried compiling Dumux on an M1 mac. Thanks for trying. I can give it a try myself in the following weak. Until then, I assume some compiler flag makes it choke. Probably “-march=native” which might not be supported yet on M1. You find the

Re: [DuMux] (no subject)

u.edu> > Phone: +1-(785)-979-2664 > > > > From: Pham, Vuong Van > Sent: Thursday, November 18, 2021 5:02 PM > To: Timo Koch mailto:koch_t...@hotmail.com>>; DuMux > User Mailing List <mailto:dumux@listserv.uni-stuttgart.de>> > Subject: RE: [DuM

Re: [DuMux] Dispersion in radial flow

resent in the model). > > Best regards, > > Dmitry > > > > > >> On 22.11.2021 17:24, Timo Koch wrote: >> Hi Dmitry, >> >> Dumux differentiates between diffusion (molecular diffusion) and dispersion >> (a scale effect). >> >>

Re: [DuMux] DuMux Digest, Vol 128, Issue 12

> On 22. Nov 2021, at 19:09, Ana Carolina Loyola > wrote: > > Dera Timo Koch, > > Yes, this example runs fine ( test_1p_incompressible_tpfa_quad ) here. I > don't know if this issue is related to the Box method, but so far when I run > a

Re: [DuMux] DuMux Digest, Vol 128, Issue 12

> > ___ > > DuMux mailing list > > DuMux@listserv.uni-stuttgart.de <mailto:DuMux@listserv.uni-stuttgart.de> > > https://listserv.uni-stuttgart.de/mailman/listinfo/dumux > > <https://listserv.uni-stuttgart.de/mailman/

Re: [DuMux] Dispersion in radial flow

Hi Dmitry, Dumux differentiates between diffusion (molecular diffusion) and dispersion (a scale effect). (Dispersion used to be implemented in dumux 2 but some problems and had not been reimplemented after 3.0 since no-one had an immediate used case.) However, you are asking in the right time

Re: [DuMux] Step post-processing in MPI mode

Hi Dmitry, you want to only iterate over the internal cells to not sum up things twice on the process boundary, to this end use for (const auto& element : elements(this->gridGeometry().gridView(), Dune::Partitions::interior)) After the loop you actually have to sum up over the process if you

Re: [DuMux] (no subject)

Dear Vuong, it looks like you ran cmake in the source directory. You have to run „make executable“ in the _build_ directory. Per default dumux/build-cmake/… Timo > 18. nov. 2021 kl. 08:19 skrev Kilian Weishaupt > : > >  > Dear Vuong Van, > > did you modify anything in your CMakeLists.txt

Re: [DuMux] (no subject)

Dear Vuong, welcome to the dumux mailing list. The attached output looks like you already ran the script twice and did some things in between. To be able to reproduce any error can you go to a new empty and run the script in there and post the entire terminal output? Timo > 15. nov. 2021

Re: [DuMux] Question to SPE5 model for porous media

Dear Johannes, what is your question? In order to help, we would need a bit more context. What are you trying to solve? What do you mean by it did not work? What is the error message? Best wishes Timo > On 8. Nov 2021, at 16:19, Bauer Johannes Fabian > wrote: > > Dear Dumux-Team, > > I

Re: [DuMux] Error in installing Dumux

Hi Amir, looks like your file system is corrupted somehow or you don’t have permissions to create files on this file system. On which OS are you running this? What does “git clone https://gitlab.dune-project.org/core/dune-common.git” return? Or even simpler: "touch bla.txt”? If you are

Re: [DuMux] [Dumux] Black-oil-model

Hi Johannes, it’s definitely possible. I guess a good starting point for a black-oil type model is 3p(ni), ni: non-isothermal. It would depend on the specific type of model, particular formulation (total velocity, global pressure for example) and solution scheme that you want to use, that

Re: [DuMux] How to make DuMux call VolumeVariables::update less frequently?

But if you stay on a single node and assembly is your main cost then this might indeed be enough and the most performant option. Anyway thanks a lot for your feedback! Timo > Best regards, > > Dmitry > > > > > > On 14.10.2021 12:39, Timo Koch wrote: >> &

Re: [DuMux] How to make DuMux call VolumeVariables::update less frequently?

races and for expensive assembly per element as you certainly have you should get a basically optimal speed-up when using more cores (even on smallish grids). I’d be happy if you give it a try and give feedback if this worked. Best wishes, Timo > On 14. Oct 2021, at 11:39, Timo Ko

Re: [DuMux] How to make DuMux call VolumeVariables::update less frequently?

> On 14. Oct 2021, at 00:12, Dmitry Pavlov wrote: > > Hello, > > I am making a compositional oil-gas simulator in DuMux. An important part of > it, and arguably the most CPU-demanding, is the determination if the given > composition is going to stay in one phase (stability test) and, in

Re: [DuMux] How to stop a gas injection?

Hi Kenza, in your problem you have a “setTime” function. Do you actually call this function in your main.cc ? You have to tell the problem what the current time is at the beginning of each time loop. In an implicit Euler time stepping procedure you probably want to do

Re: [DuMux] Using Embedded Python in DuMux

Big advantage of pybind11 is also that is does automatic type conversion for you between Python and C++, so you can easily convert between C++ std::vector and Python’s list etc... > On 9. Sep 2021, at 15:03, Timo Koch wrote: > > Hi Leo, > > I would definitely look into thi

Re: [DuMux] Using Embedded Python in DuMux

Hi Leo, I would definitely look into this with pybind11 which reads very nice on the C++ side too. There are some instructions at the end here dune-project.org/doc/pythonbindings and in dune-common/dune/python/test/ you also find two small examples for embedded Python. We also have some

Re: [DuMux] MaterialLawParams -> FluidMatrixInteraction

sry, delete ‘static’ from my comment Timo > 5. sep. 2021 kl. 17:17 skrev Timo Koch : > >  should be ok. Creating the new law is the same cost as creating the old > params, because the static functions have zero size. > > Timo > >>> 5. sep. 2021

Re: [DuMux] MaterialLawParams -> FluidMatrixInteraction

K? > > Best regards, > > Dmitry > > > > > >> On 05.09.2021 03:32, Timo Koch wrote: >> Hi Dmitry, >> >> sure this is still possible. If you want you can think of the >> FluidMatrixInteraction as a collection of Mate

Re: [DuMux] MaterialLawParams -> FluidMatrixInteraction

Hi Dmitry, sure this is still possible. If you want you can think of the FluidMatrixInteraction as a collection of MaterialLawParams only that these “params” now additionally have functions like sw(pc). So they are the (parametrized) “law”. Instead of separating parameters and (stateless)

Re: [DuMux] DuMux Digest, Vol 123, Issue 5

Hi Kenza, can you be more specific about where these parameters are? In a quick search I found this https://git.iws.uni-stuttgart.de/dumux-repositories/dumux/-/blob/master/dumux/material/fluidmatrixinteractions/2p/interfacialarea/exponentialcubic.hh

Re: [DuMux] DuMux Digest, Vol 123, Issue 3

I there a reason you have to use such an old version (it’s more than 5 years old)? I will be hard to find someone that can help you with this code (or has time to look into this). We don’t officially support and test versions before Dumux 3 anymore. If you are just starting a new application I

Re: [DuMux] Thermal conduction in Dumux

Dear Mahmoud, as virtually every component in Dumux the thermal conductivity constitutive law is exchangeable. The default for the OnePNI model is “ThermalConduvitivityAverage" (here

Re: [DuMux] Temporal discretization in Dumux

t; > > > Regards, > > > Mahmoud > > > > > > > > > > > > > > > ----- > Il 2021-05-21 16:30 Timo Koch ha scritto: >> You can currently use Forward or Backward differences. But only first >> order one-step methods are implemented. >> Higher

Re: [DuMux] Temporal discretization in Dumux

You can currently use Forward or Backward differences. But only first order one-step methods are implemented. Higher order time one step multi-stage time discretisation methods sort of work on an experimental branch (track https://git.iws.uni-stuttgart.de/dumux-repositories/dumux/-/issues/940

Re: [DuMux] Parallel Multidomain Facet

Hi Ana, unfortunately this is not implemented yet. FoamGrid is not parallel yet which would be needed for parallel assembly. The iterative solvers also currently don’t work in the multidomain setting although there is some ongoing work there (e.g. merge request !1950). Timo > On 18. May 2021,

Re: [DuMux] how to add a kinetic dissolution rate

> On 10. May 2021, at 21:16, Johannes Hommel > wrote: > > Dear Kenza, > > what dissolution do you mean? A dissolution of CO2 into water/brine or a > CO2-caused dissolution of solids? > > A kinetic dissolution of solid could be implemented using the 2pncmin model. > > A kinetic in the

Re: [DuMux] [DuMuX] Numerical derivatives for pressure gradient with the Box method

el.pdf On 28.05.2020 19:07, Timo Koch wrote: Hi Dmitry, glad to hear! We have been discussing the issue in the last developer meeting. The fix introduces a significant overhead in the assembly but is not necessary for most simulations. Also, without the fix although the Jacobian is just an approxi

Re: [DuMux] Setting zero transmissibility along an interface between 2 layers

(), phaseIdx)*tFactor ; >else >return flux_(problem, fluxVarsCache, > fluxVarsCache.advectionPrimaryDataHandle(), phaseIdx)*tFactor ; > > > I would appreciate your feedback and help to fix this. > > > Best Regards, > > Mahmoud > >

Re: [DuMux] Numerical bug in invertCubicPolynomial

Hi Dmitry, thanks for looking into this. It would be very nice if you could provide an MR with the proposed fix together with some unit tests that test corner cases (in test/common/test_math.cc). Timo > On 29. Apr 2021, at 16:05, Dmitry Pavlov wrote: > > Hello, > > I would like to propose

Re: [DuMux] Setting zero transmissibility along an interface between 2 layers

Dear Mahmoud, I understand it in a way that you want to have no mass flow over these boundaries but there can be heat conduction. You could solve this either by using a multidomain model where you only solve heat conduction in some parts of the mesh. This is a bit more evolved. A simpler

Re: [DuMux] Override TwoPNCVolumeVariables::completeFluidState

> On 26. Apr 2021, at 17:07, Dmitry Pavlov wrote: > > Kilian, > > Thank you, I wiped out my GridVariables and GridVolumeVariables as you > suggested. > > It turns out the problem was that I inherited TwoPNCVolumeVariables in > MyVolumeVariables. Once I stopped doing that and just

Re: [DuMux] Question on sourceAtPos function

5);" it gave me > this error: > > > error: ‘liquidEnthalpy’ is not a member of > ‘Dumux::OnePNIProblem::FluidSystem’ > {aka ‘Dumux::FluidSystems::OnePLiquid > >’} > 239 | const auto enthalpy = FluidSystem::liquidEnthalpy(293, > 13.0e5); > &g

Re: [DuMux] Question on sourceAtPos function

e: candidate expects 1 argument, 0 provided > > > I would really appreciate it if you let me know how to correctly accurate for > the 100 m^3/day volumetric rate either using the sourceAtPos function with > elaborating the mistakes or using the pointSource function showing how to > impose the 293 kelvin temperature for the injection well. > I'm ext

Re: [DuMux] Question on sourceAtPos function

Dear Mahmoud, what you are doing seems correct to me. You can get the relevant cell volume with “scv.volume()” so you don’t have to compute that yourself. You can also get the density with "elemVolVars[scv].density(phaseIdx)". But that shouldn’t change anything if you already used the correct

Re: [DuMux] Oil + gas simulation

Hi Dmitry, I think you are misinterpreting the threshold adjustment. If the threshold is higher if the phase state already switched in the last iteration. That means in the last iteration both phases where present and one disappeared. In order to prevent the Newton from switching back and

Re: [DuMux] Importing an Eclipse Grid

Hi Mahmoud, dumux is tested with OPM version: release/2020.04. Can you try this and see if the error persists? Also it might be necessary to add "-DUSE_MPI=ON” to your opts file. You might also need BLAS, LAPACK, Boost, SuiteSparse, Zoltan (see opm (optional) dependencies). This information is

Re: [DuMux] Error building a new module (dumux-trial)

Dear Mahmoud, happy to help. Thank you for reporting your problems. It helps us improving the code and the documentation! Best wishes, Timo > On 8. Mar 2021, at 13:25, Mahmoud Atef Mahmoud Mohamed Aboelseoud S277151 > wrote: > > Dear Timo, > > I really can't thank you enough ! After

Re: [DuMux] Error building a new module (dumux-trial)

Hi Mahmoud, thanks for making us aware that the example in the getting started guide requires Suitesparse. This shouldn’t be the case and is an oversight on our side. We’ll change that in the next days (most likely will be fixed by Monday evening). Otherwise your setup seems to look fine. I

Re: [DuMux] Error building a new module (dumux-trial)

Hi Mahmoud, in oder to provide help we would need the following information: * which OS are you working on? (We only support Linux, macOS should be fine but might have some more complications sometimes) * What commands did you type to produce the error? (Minimal working example) * Which

Re: [DuMux] 1p and symmetry

le are also implemented for Dune::Float128. Timo > > Best regards, > > Dmitry > > > On 11.11.2020 20:16, Timo Koch wrote: >> >>> On 11. Nov 2020, at 17:32, Dmitry Pavlov wrote: >>> >>> Hello, >>> >>> I am runni

Re: [DuMux] 1p and symmetry

> On 11. Nov 2020, at 17:32, Dmitry Pavlov wrote: > > Hello, > > I am running a fairly simple 1p porous flow simulation which is on a 2D > rectangular grid. In a particular case I came at, the initial and boundary > conditions, as well as spatial parameters, do not depend on the "y" >

Re: [DuMux] Fugacity of Brine in gas phase

Hi, if I understand the comment with the TODO correctly, the constraint solver will then not solve the correct phase composition in your case. If your fugacity depends on the mass fraction the system to solve for the phase composition is non-linear and then you need to e.g. use a Newton

Re: [DuMux] Switching linear solver causes simulation running on just a part of a 2D grid

> On 2. Sep 2020, at 18:39, Dmitry Pavlov wrote: > > Timo, > > Thank you. I am exploring my options. Are all three steps (Assemble, Solve, > Update) normally parallelized? yes > > I also noticed that AMG does itself use UMFPack as a coarse solver, it it > (well, SuiteSparse) is installed

Re: [DuMux] Switching linear solver causes simulation running on just a part of a 2D grid

-- _ Timo Koch phone: +49 711 685 64676 IWS, Universität Stuttgart fax: +49 711 685 60430 Pfaffenwaldring 61 email: timo.k...@iws.uni-stuttgart.de D-70569 Stuttgart url

Re: [DuMux] three domain modeling

Hi, I don‘t think you need 3 models for the dumux multidomain in your case. You can just use the default setup with a different grid. As you just have Stokes and Darcy just make a grid that has some hole/channel spared out (or whatever setup you want). For that you can use e.g. the

  1   2   3   >