On 08/02/2016 02:18 PM, Vinetou Incucuna wrote:
I have tried to output the solutions from the step one as
well as from the step three. There is some development, but
reasonable,i.e. solution from the first step and the third step look
similar (except the first iteration)
In the attachment you can see the results of my computations
and the partitions.
It appears to me, that the problematic edges in the solution (both
pressure and velocity) correspond the borders between the partitions.
Maybe you could give me a gist, what is going on.
It's hard to tell. The code doesn't look wrong, but without the context
I also can't say that it's right (and I wouldn't have the time to look
at it in sufficient detail even if I had the context). My best guess is
that you may be assembling on the wrong cells. E.g., you may be
assembling also on ghost cells, instead of only on locally owned ones.
Of course I don't know whether that's true.
Are you running on a mesh without hanging nodes? What happens if you run
just on a regular square instead of a complicated domain? My strategy
for debugging is always to make the problem simpler and simpler.
Best
W.
--
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 dealii+unsubscr...@googlegroups.com.
For more options, visit https://groups.google.com/d/optout.
