Почему надо думать, что такие источники будут притягиваться?
Если вместо относительно не удобного трехмерного случая взять пример двумерных источников (в трех измерениях это параллельные нити-источники), то можно воспользоваться красивыми и простыми методами комплексного потенциала, замечательно моделирующими двумерные стационарные задачи течения идеальной жидкости. Комплексный потенциал течения, созданного двумя одинаковыми по обильности источниками будет иметь вид:

,
где

- вещественное число, определяющее обильность источников,

- комплексная переменная,

и

- координаты в комплексных числах, задающие положения обоих источников на комплексной плоскости.
Не сложно показать, что сила взаимодействия такой пары двумерных источников идеальной жидкости направлена в сторону отталкивания. То есть, два одинаковых по знаку точечных источника идеальной жидкости, как и два одинаковых электрических заряда - отталкиваются. А два противоположных по знаку - притягиваются. Тоже самое будет и в трехмерном случае, только тут уже не работает метод комплексного потенциала, но можно рассмотреть обычный скалярный потенциал...
У автора поста, в котором одинаковые по знаку источники - отталкиваются, возможно, рассматривалась не идеальная вязкая жидкость, которая может вести себя не в полной аналогии с электрическими зарядами.