Re: [DuMuX] Precompiled headers in Dumux?

Hello again, I took a look at ccache but it doesn't provide the functionality I need. As far as I understand, it is useful if you use 'make clean' and 'make' on whole projects but not if you compile a single executable. I am not an expert but for cotire it seems to be the same. One solution to

[DuMuX] Precompiled headers in Dumux?

Hello Dumux, Is there a way to include precompiled headers into Dumux? I would like to include the exprtk parser into my code which allows to use fancy mathematical expressions at runtime. However, if I do so, the compilation will take forever since exprtk.hpp is huge (1.4 MB). Therefore, the

Re: [DuMuX] gravity dependency on global postion

Hi Timo, Thanks for the correction ☺ Georg Von: Dumux [mailto:dumux-boun...@listserv.uni-stuttgart.de] Im Auftrag von Timo Koch Gesendet: Dienstag, 8. Januar 2019 09:22 An: DuMuX User Forum; Felix Feldmann Betreff: Re: [DuMuX] gravity dependency on global postion Hi Felix, hi Georg,

Re: [DuMuX] gravity dependency on global postion

Hi Felix, You could try to run you problem with valgrind. Just type ‘valgrind ./my_problem’ to do that. Then you might see where the error comes from. If the output does not give you an idea try to compile your problem with debug-opts and rerun it with valgrind. Best regards Georg Von: Dumux

Re: [DuMuX] gravity dependency on global postion

Hi Felix, You should try the following: const GravityVector (const GlobalPosition ) const { // initialize a gravity vector with zeros GravityVector gravity(0.0); // modify the entry for the x (???) direction according to your needs if (onLeftArea_( pos)) gravity[0]=100; else // don’t use else

Re: [DuMuX] Property system and definition of static members

Hi Timo, You are right. It works this way: template typename GET_PROP_TYPE(TypeTag, MyType) MyClass::myStaticMember_; I went to the Frühlingsfest yesterday and was therefore not able to understand what my compiler tried to tell me :) Thanks Georg Von: Dumux

Re: [DuMuX] Flux calculation at the scv center (box method)

Hi Bernd, I already implemented the method I described and it seems to work well. Basically, I solve [cid:image002.png@01D218DF.9B80AA40] And [cid:image003.png@01D218DF.9B80AA40] Where i is the current density and N_0 is the molar flux density of H2O. They are defined as

[DuMuX] Flux calculation at the scv center (box method)

Hello Dumuxers, For my fuel cell model I would like to implement the following energy balance in the membrane: [cid:image001.png@01D218A5.6DD790C0] The last term on the right describes the effect of ohmic or joule heating due to charge transport. Here, a dot product of the current density has

Re: [DuMuX] Multiphase balance equations

Hi Bernd, thanks for the information. Now I'll have to understand it... Georg Von: Dumux [mailto:dumux-boun...@listserv.uni-stuttgart.de] Im Auftrag von Bernd Flemisch Gesendet: Donnerstag, 4. August 2016 16:42 An: DuMuX User Forum Betreff: Re: [DuMuX] Multiphase balance equations Maybe from

[DuMuX] Multiphase balance equations

Hello Dumuxers, I am wondering about the advective and diffusive term of the balance equations in the various multi-phase Dumux models. For me, the question arises whether these term should additionally contain the porosity times the phase saturation or not. Starting from the Reynolds

Re: [DuMuX] flux calculation with compute flux

Hi again! For the outflow boundary condition, also computeFlux is used on the boundary. I am not 100% sure but since this is just an approximation it would probably be better to rewrite the evalOutflowSegment-function of the boxlocalresidual according to your explanations? Kind regards Georg

Re: [DuMuX] flux calculation with compute flux

Hi Bernd, Thanks for the information. I was not aware that computeFlux on a boundary face for the box method will be just an approximation. I will have to improve the code. Kind regards Georg Von: Dumux [mailto:dumux-boun...@listserv.uni-stuttgart.de] Im Auftrag von Bernd Flemisch Gesendet:

Re: [DuMuX] flux calculation with compute flux

Update 2: I forgot to change a similar function for the other subdomain. Now everything is fine. I guess the code is correct. Von: Dumux [mailto:dumux-boun...@listserv.uni-stuttgart.de] Im Auftrag von georg.fut...@dlr.de Gesendet: Dienstag, 15. März 2016 15:13 An:

Re: [DuMuX] flux calculation with compute flux

Update: I just ran the code and got a read error so maybe something is still wrong with the code but I think that I go into the right direction. Georg Von: Dumux [mailto:dumux-boun...@listserv.uni-stuttgart.de] Im Auftrag von georg.fut...@dlr.de Gesendet: Dienstag, 15. März 2016 15:06 An:

[DuMuX] flux calculation with compute flux

Hello Dumux, I tried to calculate the fluxes over a certain part of the boundary of the model domain. For this purpose I used an ElementIterator to iterate over the elements of the grid. In each element I iterate over the intersections and in each intersection over the vertices. With

Re: [DuMuX] Multidomain in parallel

Hello Markus, hello Bernd, > Does that mean you used SuperLU or SuperLU-Dist? In this case No just the standard SuperLU. > Anyway, increased walltime sounds fishy. How much does it increase? On two cores the walltime increased by 27%, on four by 114%. >A good test would be to setup the AMG

Re: [DuMuX] Multidomain in parallel

Hi Bernd, I tried running my multidomain model with the AMG-Solver (not yet in parallel). First error that popped up was about AmgTraits not being set. As a first try I used the AmgTraits defined in the el2ppropertydefaults. However, it does not converge. For larger time steps I get warnings

Re: [DuMuX] Multidomain in parallel

Hi Christoph, We ran our multidomain application in parallel using SuperLU as solver. However, this lead to an increase in wall time. Since, for now, we do not solve Stokes flow it might be possible to use an iterative solver for our fuel cell models and speed up the simulation through

[DuMuX] Multidomain in parallel

Hello Dumuxers, This question may be a little naïve but would it be possible to run a multidomain application in parallel? What would I have to do? Is there a test case in Dumux which can be run in parallel? I found the Thesis "Electrodiffusion Models of Axon and Extracellular Space Using the

[DuMuX] Setting Properties

Hello Dumuxers, I have a question concerning the property system. Is it possible to set a property depending on the value of another property? I have something in mind like this: SET_PROP(MyProblem, MyProperty) { private: enum { boolValue = GET_PROP_VALUE(TypeTag, MyBoolValue) }; public:

Re: [DuMuX] time-dependent boundary conditions

Thanks for your help! Von: Dumux [mailto:dumux-boun...@listserv.uni-stuttgart.de] Im Auftrag von Bernd Flemisch Gesendet: Freitag, 22. Januar 2016 09:48 An: DuMuX User Forum Betreff: Re: [DuMuX] time-dependent boundary conditions Hi Georg, to answer your questions: 1. There are several tests.

[DuMuX] time-dependent boundary conditions

Hello Dumuxers, Is there an example with time-dependent boundary conditions? I am wondering which time to use in the dirichletAtPos and solDependentNeumann functions. Is it timeManager().time() or timeManager().time() + timeManager().timeStepSize()? Best regards Georg Futter

Re: [DuMuX] A fundamental question concerning diffusion in Dumux

Hallo Rainer, Let's meet next week any day is fine for me. Best regards Georg Von: Dumux [mailto:dumux-boun...@listserv.uni-stuttgart.de] Im Auftrag von Rainer Helmig Gesendet: Freitag, 4. Dezember 2015 15:56 An: DuMuX User Forum Betreff: Re: [DuMuX] A fundamental question concerning diffusion

[DuMuX] Multidomain problem in 3D

Hello Dumux, Does anyone know where the dimensionality of a multidomain problem is defined? I am unable to find it and I would like to test it 3D. Thanks for your help. Best regards Georg Futter -- German Aerospace Center (DLR) Institute of Engineering Thermodynamics

[DuMuX] A fundamental question concerning diffusion in Dumux

Hello Dumuxers, I was wondering why the diffusive fluxes in Dumux are defined as D*rho*grad(mol-or-massfraction). Typically one would use D*grad(c) (e.g. Ficks law) where c=rho*mol-or-massfraction. Using the Dumux equation means that local differences in the density are neglected for diffusive

[DuMuX] UMFPack

Hello Dumuxers, I would like to test UMFPack on my application. However, I failed to compile it with this solver. I got the latest version of SuitSparse (4.4.5), compiled it and successfully tested UMFPack. When I run dunecontrol it tells me that UMFPack is found. Now, if I compile my

Re: [DuMuX] UMFPack

Hi Christoph, I get the same compilation errors when compiling umfpacktest of dune-istl and a dumux test case. How do I reconfigure dune-istl? Bye Georg -Ursprüngliche Nachricht- Von: Dumux [mailto:dumux-boun...@listserv.uni-stuttgart.de] Im Auftrag von Christoph Grüninger Gesendet:

Re: [DuMuX] Boundary conditions

Hi Alex, I tried to do the same thing as in the 1p2coutflowproblem only with a 2p3c Fluidsystem and the mpnc model. Initially, only the gas phase is present in the whole domain (Sg=1) and I set Dirichlet boundary conditions for all primary variables on the inlet (pg = 1.6 bar, Sg=1 and the

Re: [DuMuX] Boundary conditions

Hi Alex, I attached the relevant pieces of code in the file. Just to make sure we talk about the same thing: in the mpnc model there are N+M equations/primary variables for an isothermal system. This gives 5 equations for the 2p3c case and the primary variables are: Sg, pg, and three

Re: [DuMuX] Boundary conditions

Dear Alex, Thanks for your reply! 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)

[DuMuX] GET_RUNTIME_PARAM macro

Hello Dumux, I was wondering whether it is possible to read arrays from the input file using the GET_RUNTIME_PARAM macro. Or is it possible to read values in a loop? I mean something like: For (int i=0; inumOfStuff; i++) Value[i]=GET_RUNTIME_PARAM(TypeTag, Scalar,

[DuMuX] creation of multidomain grids

Hello DuMuXers, The grid geometry I want to create is attached as .jpg-file. I want to create a sub domain Omega1 as depicted. Here is my code from the coupled problem: void initializeGrid() { MDGrid mdGrid = this-mdGrid(); mdGrid.startSubDomainMarking(); //

Re: [DuMuX] Check for phase presence in the mpnc model

Hi Bernd, I already use SuperLU. Btw what about the PARDISO solver? Can it be used for multidomain applications and could this be an improvement? Kind regards Georg Von: Dumux [mailto:dumux-boun...@listserv.uni-stuttgart.de] Im Auftrag von Bernd Flemisch Gesendet: Mittwoch, 27. Mai 2015 12:36

Re: [DuMuX] Check for phase presence in the mpnc model

Hi Bernd, I checked for the phase presence in the evalCoupling function of my localoperator-file. Here I couple the water component balance of the mpnc model to my membrane model by setting CouplingOutflow boundary conditions for the membrane and CouplingInflow for the mpnc domain. As I

Re: [DuMuX] Check for phase presence in the mpnc model

Hi Bernd, That's a good idea. I will try it and see what happens. Kind regards Georg Von: Dumux [mailto:dumux-boun...@listserv.uni-stuttgart.de] Im Auftrag von Bernd Flemisch Gesendet: Dienstag, 26. Mai 2015 15:44 An: DuMuX User Forum Betreff: Re: [DuMuX] Check for phase presence in the mpnc

Re: [DuMuX] Check for phase presence in the mpnc model

By the solver I mean the linear solver. At least the values are correct. Maybe the numerical error from the coupling just gets too high if I use the old solution? Georg Von: Dumux [mailto:dumux-boun...@listserv.uni-stuttgart.de] Im Auftrag von Bernd Flemisch Gesendet: Dienstag, 26. Mai 2015

Re: [DuMuX] Check for phase presence in the mpnc model

No, this does not work at all. The solver does not converge. Von: Dumux [mailto:dumux-boun...@listserv.uni-stuttgart.de] Im Auftrag von georg.fut...@dlr.de Gesendet: Dienstag, 26. Mai 2015 15:47 An: dumux@listserv.uni-stuttgart.de Betreff: Re: [DuMuX] Check for phase presence in the mpnc model

[DuMuX] Check for phase presence in the mpnc model

Hello DuMuX, I started working with the mpnc model and I wonder how I can check for the phase presence. Just accessing the phase saturation does not seem to be a safe way since during newton iterations the value may differ from zero even though the phase is not present (am I right here?).

[DuMuX] Cauchy boundary conditions

Hello DuMuX users, Is there a possibility to set a Cauchy boundary condition in DuMuX. If so, how does it work? Best regards Georg Futter -- German Aerospace Center (DLR) Institute of Engineering Thermodynamics | Computational Electrochemistry | Pfaffenwaldring 38-40

Re: [DuMuX] strange artifacts

Hi Bernd, Thanks a lot for the quick help. This solved the problem. Best regards Georg Von: Dumux [mailto:dumux-boun...@listserv.uni-stuttgart.de] Im Auftrag von Bernd Flemisch Gesendet: Freitag, 17. April 2015 11:06 An: DuMuX User Forum Betreff: Re: [DuMuX] strange artifacts Hi Georg, it

Re: [DuMuX] strange artifacts

Hi Christoph, Here it is. Georg -Ursprüngliche Nachricht- Von: Dumux [mailto:dumux-boun...@listserv.uni-stuttgart.de] Im Auftrag von Christoph Grüninger Gesendet: Freitag, 17. April 2015 10:52 An: dumux@listserv.uni-stuttgart.de Betreff: Re: [DuMuX] strange artifacts Hi Georg, can you

Re: [DuMuX] usage of restart file

Hi Thomas, Thanks a lot! :-) Best regards Georg -Ursprüngliche Nachricht- Von: Dumux [mailto:dumux-boun...@listserv.uni-stuttgart.de] Im Auftrag von Thomas Fetzer Gesendet: Freitag, 30. Januar 2015 14:30 An: dumux@listserv.uni-stuttgart.de Betreff: Re: [DuMuX] usage of restart file Hi

[DuMuX] usage of restart file

Hello Dumux, I tried to use a restart file for the first time but it didn't work. So I guess I am missing something. My input file looks like this: ### # Simulation restart # # DuMux simulations can be restarted from *.drs files # Set

[DuMuX] Where to set ISOTHERMAL?

Hi Dumux, Just a quick and probably easy question: The fuel cell model can be run isothermal or non-isothermal. This depends on the ISOTHERMAL variable. For example: #if ISOTHERMAL #include ... #else #include... #endif Where is it set? Best regards Georg Futter

Re: [DuMuX] Where to set ISOTHERMAL?

Hi Alex, Thanks for your help. For the fuelcell_lisaproblem ISOTHERMAL is defined in the application file. Have a nice weekend Georg Von: Dumux [mailto:dumux-boun...@listserv.uni-stuttgart.de] Im Auftrag von Alexander Kissinger Gesendet: Freitag, 17. Oktober 2014 14:22 An: DuMuX User Forum

[DuMuX] Coupled problems

Hello Dumux experts, I am currently working on the coupling between the (Acosta) fuel cell model and my model for the polymer membrane. For this purpose I adapted the 2cstokes2p2c-test of dumux stable. The coupled problem compiles and runs now but I'm getting numerical problems. To check what

Re: [DuMuX] Eigenvalues of the jacobian matrix

Hi Christoph, Thank you for your quick answer. I was wondering whether I could use the functions in eigenvalues.hh. Is this not a good idea? Regards Georg -Ursprüngliche Nachricht- Von: dumux-boun...@listserv.uni-stuttgart.de [mailto:dumux-boun...@listserv.uni-stuttgart.de] Im Auftrag

[DuMuX] Eigenvalues of the jacobian matrix

Hello Everyone, To investigate some instabilities in my model I would like to have a look at the eigenvalues of the jacobian matrix. Is there an example how this can be done? Thanks Georg -- Deutsches Zentrum für Luft- und Raumfahrt e.V. (DLR) Institut für Technische

[DuMuX] Multiple simulations in a row

Hello Dumux, as some of you might already know I am working on a dumux fuel cell model here at DLR. For my project I need to run lots of quick simulations (10^6) with varying boundary conditions. Is there a possibility to do this, and if not, what would be the way to do it? Best regards Georg