The moulds' filling process may generate flaws consisting of small gas bubbles trapped inside the material, which weaken the solidity of the casted piece. We consider here the inverse problem of determining these small size flaws' locations from velocities boundary measurements. The fluid flow is described by a simplified model based on the Stokes system. A numerical algorithm based on the topological sensitivity analysis applied to an energy-like misfit functional is worked out to that end.