The methodology presented in Part I of this work is applied to a large number of pressure-ternperature flash calculations, and to the automated construction of constant temperature pressure-composition phase diagrams, and constant pressure temperature-composition phase diagrams for binary mixtures modeled with an augmented van der Waals equation of state. An automated prototype implementation of the algorithm is developed for this purpose. We follow the classification of Scott and van Konynenburg [R.L Scott, P.H. van Konynenburg, Discuss. Faraday Soc. 49 (1970) 871 and present phase diagrams corresponding to non-azeotropic mixtures of the five main types of fluid phase behavior (I-V), studying in detail representative diagrams at constant pressure and constant temperature. Special attention is given to the solution of numerically problematic equilibrium regions, such as those close to three-phase equilibria where metastable and unstable critical points can also be found. Of the order of 104 flash calculations at varying temperatures and pressures, and for different intermolecular parameters of the components in the mixture, have been cat tied out. The algorithm provides the correct stable equilibrium state for all of the points considered. Despite the fact that our implementation is not optimised for performance, we find that the algorithm identifies the stable Solution in difficult regions of the phase space without any penalty in terms of computational time, when compared to simpler regions. (C) 2008 Elsevier B.V. All rights reserved.