I think you have used correct boundary conditions. You need to find the error somewhere else.nIs your model 2D? The Swiss-roll recuperator is tubes with circular cross-sections, but if your model is a 2D and the pictures are showing your domain, note that the solver takes your model as a rectangular channel with theoretically infinite depth (perpendicular to the screen) so that gradients in the z-direction are zero.n