I am currently working through the tutorial pages and have a problem understanding the boundary conditions. I am currently stuck at the tutorial dealing with the Navier-Stokes equations for Cylinder Flow (https://fenicsproject.org/pub/tutorial/html/._ftut1009.html). Why is it that the inflow nodes defined by
at first do not have the values of the inflow profile? What I noticed is that it that the inflow-node values start at 0 (of course, this is the default initialization value) and then quickly converge to the inflow-profile after a few time steps. Why does the employed DirichletBoundary conditions not hard-set these values such that they are enforced within the first time step?
There is very little value in setting these, as all those dofs at the boundary are strongly enforced in the first tentative solve. You would only get a slight difference in the region between the boundary and the first element inside the boundary…
Thank you for the quick response! To check my understanding: The boundary condition is only enforced in the intermediate solution but not at the final solution? Therefore, it is not a bug but rather a consequence of the technique that was employed?