Re: [DuMuX] Problem with singular matrix in 1D for thermalnonequilibrium

2019-03-25 Thread Katharina Heck
Hi Max, The heat source in that test is applied directly at the boundary of the porous medium, so the value of the solid thermal conductivity at that face influences the test result a lot. In the dumux version 2.12 we applied a arithmetic mean for the calculation of the solid thermal

Re: [DuMuX] Problem with singular matrix in 1D for thermalnonequilibrium

2019-03-11 Thread Katharina Heck
Hi Max, I am not sure if I get your question correctly. You are using the mpnc-thermalnonequilibrium test and you change the solid thermal conductivity in the outflow? From what I understand about that test it would seem logical to me that the solid thermal conductivity changes the results

Re: [DuMuX] Problem with singular matrix in 1D for thermalnonequilibrium

2019-03-10 Thread Flemisch, Bernd
Hi Max, > I'm working with version 3.0 and have made the following changes. > > Old: > > using LinearSolver = Dumux::AMGBackend; > auto linearSolver = std::make_shared(leafGridView, > fvGridGeometry->dofMapper()); > > New: > > using LinearSolver = Dumux::UMFPackBackend; > auto

Re: [DuMuX] Problem with singular matrix in 1D for thermalnonequilibrium

2019-03-07 Thread Flemisch, Bernd
Hi Max, you should first check if the problem is the linear solver by switching to a direct solver. If you use Dumux 2.x, that's setting a SET_TYPE_PROP(YourTypeTag, LinearSolver, UMFPackBackend); in the problem file, for 3.0 you have to change the corresponding line in the main file. Let us

[DuMuX] Problem with singular matrix in 1D for thermalnonequilibrium

2019-03-07 Thread Maximilian Johannes Lueftner
Dear DuMuX experts, i have adapted the simulation from the test folder „dumux/test/porousmediumflow/mpnc/implicit/thermalnonequilibrium" . I used the YASPGrid for 1D with grading: struct Grid { using type = Dune::YaspGrid<1, Dune::TensorProductCoordinates >; }; With 800 Elements (400 in