Nickel based superalloys have superior high temperature mechanical strength, corrosion and creep resistance in harsh environments and found applications in the hot sections as turbine blades and turbine discs in jet engines and gas generator turbines in the aerospace and energy industries. The efficiency of these turbine engines depends on the turbine inlet temperature, which is determined by the high temperature strength and behavior of these superalloys. The microstructure of nickel based superalloys usually contains coherently precipitated gamma prime (?) Ni3Al phase within the random solid solution of the gamma () matrix, with the ? phase being the strengthening phase of the superalloys. How the alloying elements partition into the and ? phases and especially in the site occupancy behaviors in the strengthening ? phases play a critical role in their high temperature mechanical behaviors. The goal of this dissertation is to study the site substitution behavior of the major alloying elements including Cr, Co and Ti through first principles based calculations. Site substitution energies have been calculated using the anti-site formation, the standard defect formation formalism, and the vacancy formation based formalism. Elements such as Cr and Ti were found to show strong preference for Al sublattice, whereas Co …