Dear Niek,
This does indeed look strange. I cannot give an immediate answer
unfortunately. It looks like the problem appears in a single element
layer close to the boundary and does not seem to propagate into the
domain. Can you confirm this? Since I would expect a wrong Neumann
boundary condition to extend further into the domain (through the
diffusive term), I wonder whether you really look at a Neumann boundary
or whether it actually is a Dirichlet boundary. (If I remember
correctly, the problem setting in step-51 specified Dirichlet conditions
on that particular boundary.) In case of a Dirichlet boundary, my guess
is that you do not correctly take the time-dependent conditions into
account. Essentially, you would need to modify the constraint matrix
each time with data from project_boundary_conditions and make sure that
your solution vectors get the correct values set.
Best,
Martin
--
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.