On 3/19/25 00:51, Jean-Paul Pelteret wrote:
You might want to try either reformulating your boundary conditions as updates
to the BCs applied at the previous timestep (i.e. constrain the increment) or
you should refactor the code (assembly, constraints creation + something
before the Newton iterations start) such that you apply the constraints
directly to the solution vector itself.
There is also a discussion on boundary values for Newton methods in step-15
that you might want to read through. In short, you likely only want to apply
Dirichlet boundary conditions to the first Newton iteration, and then use zero
boundary conditions after that.
Best
W.
--
------------------------------------------------------------------------
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
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.
To view this discussion visit
https://groups.google.com/d/msgid/dealii/d229e447-a652-4877-8b03-5659c14b83f5%40colostate.edu.
