The paper deals with a fictitious domain approach applied to the Neumann boundary value problem by using boundary lagrange multipliers. Convergence results are established.The fictitious domain solver is used in optimal shape design problem.