The dark strips could just be incorrect normals produced by the filter. Recall I mentioned that the filter is a little unstable - geometric algorithms are hard to get right! You could try to generate new normals by adding a vtkPolyDataNormals
filter to the output of the boolean operations filter.
You’ll need to create additional boolean operation filters, setting one input to the accumulating running union object and the other input to the additional tube.