Orthocircles

bb = 3; ff = 0.075;
i = -2; j = 2;
ContourPlot3D[
((x^2 + y^2 - 1)^2 + z^2)*((y^2 + z^2 - 1)^2 + x^2)*((z^2 + x^2 - 1)^2 + y^2) - ff^2*(1 + bb*(x^2 + y^2 + z^2)) == 0,
{x, i, j}, {y, i, j}, {z, i, j},
PlotRange -> All, PlotPoints -> 100, Axes -> False,
 Boxed -> False,
 ContourStyle -> Directive[RandomColor[], Opacity[0.7], Specularity[White, 30]]]

The equation is derived by placing 3 circles perpendicular to each other.