Data-driven models were constructed for the mechanical properties of multi-component Ni-based superalloys, based on systematically planned, limited experimental data using a number of evolutionary approaches. Novel alloy design was carried out by optimizing two conflicting requirements of maximizing tensile stress and time-to-rupture using a genetic algorithm-based multi-objective optimization method. The procedure resulted in a number of optimized alloys having superior properties. The results were corroborated by a rigorous thermodynamic analysis and the alloys found were further classified in terms of their expected levels of hardenabilty, creep, and corrosion resistances along with the two original objectives that were optimized. A number of hitherto unknown alloys with potential superior properties in terms of all the attributes ultimately emerged through these analyses. This work is focused on providing the experimentalists with linear correlations among the design variables and between the design variables and the desired properties, non-linear correlations (qualitative) between the design variables and the desired properties, and a quantitative measure of the effect of design variables on the desired properties. Pareto-optimized predictions obtained from various data-driven approaches were screened for thermodynamic equilibrium. The results were further classified for additional properties.