Topology optimization of systems governed by external Bernoulli free boundary problems is studied. The level set approach combined with radial basis functions permit to convert the problem into a parametric optimization problem, solved by a gradient -based method.