A numerical scheme for free boundary problems arising in jet flows is given by applying the variational inequality formulation.
To find the initially unknown free boundary location, an iterative method is developed and criteria for convergence are given. Numerical results are in good agreement with those by other methods and experimental observations.