Hi Daniel, I found one problem with inp file - in line "1 0 quad 1 2 3 4". I guess the second number should be nonzero here. I have try to set 1 and get better results. Now I can see velocity changes around the input-output bounds. I think with cell face iterator I have to use cell quad iterator to set all bounds when I create my inp, am I right?
For input bound I use the same function as in the 35 tutorial. I tried to set ConstantFunction() with value about 1 or -1 for the output bound. But I stiil got the same error on the 3th step. I can see the pressure was about 10^18 before that. I get a good looking flow between input and output bounds when I use ConstantFunction() with the same value for both input and output bounds. But only for the about 5 steps. And then I can see the pressure grows dramatically. How I can chooze correct boundary conditions when I have variable input boundary? Thank you! Evgeny четверг, 28 февраля 2019 г., 18:20:48 UTC+5 пользователь Daniel Arndt написал: > > Eugeniy, > > my first guess would be that you are setting boundary conditions > inadequately and that this gives you a singular matrix. > What happens if you apply homogeneous Dirichlet boundary conditions for > the velocity everywhere? > > Best, > Daniel > -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to [email protected]. For more options, visit https://groups.google.com/d/optout.
