In this paper we present a new method for designing a nozzle. In fact the problem is to find the optimal domain for the solution of a linear or nonlinear boundary value PDE, where the boundary condition is defined over an unspecified domain. By an embedding process, the problem is first transformed to a new shape-measure problem, and then this new problem is replaced by another in which we seek to minimize a linear form over a subset of linear equalities. This minimization is global, and the theory allows us to develop a computational method to find the solution by a finite-dimensional linear programming problem.