This work is concerned with the development of a numerical model for laminar flow separation from a smooth boundary. The concept of irreversible vorticity generation is used to formulate an algorithm to predict both the location of separation and the shedding rate of vorticity. Results of flow separation from a circular cylinder are presented. The calculated vortex sheets are visualized as streak lines in the wake of the cylinder. Streamline patterns are constructed from these calculations showing location of separation and wake structure.