We can augment the overall system to implement dirichlet boundary condition in implicit case, but how to do this in multi-predictor-corrector method?
The question is too broad, and is not FEniCS-specific. We are happy to help, but we can’t replace you doing a literature review or asking suggestions to your advisor.
Closing.