Neumann bcs in mixed space

dolfinx
9

Just to be clear, currently T is a zero vector. Even then this issue occurs? (Of course, if T remains a zero vector, then you don’t need this integral…)

I am updating T in the loop with T.value[0] = n * tval0. I will try to solve it now with the suggested steps using the BlockedNewtonSolver. Thanks again for the help!