The equilibrium solubility of sulfadiazine (SD, 3) in several {acetonitrile (MeCN) + ethanol (EtOH)} mixtures at nine temperatures from T/K = (278.15 K to 318.15) has been determined by following the shake flask method. SD solubility increased with temperature-arising as well as with the MeCN proportion-increasing in the mixtures. Thus, x3 increased from 7.74 × 10−5 in neat EtOH to 6.20 × 10−4 in neat MeCN at T/K = 298.15. SD solubility was adequately correlated with a second-order multivariate equation as function of both mixtures composition and temperature. Moreover, two models including the Jouyban–Acree and Jouyban–Acree–van’t Hoff models were applied to mathematical SD solubility data modeling in solvent mixtures. The accuracy of each model is investigated by the mean relative deviations (MRD%) of the back-calculated solubility. A full predictive model was provided by training the Jouyban–Acree–van’t Hoff model with only seven experimental solubility data which provided excellent predictions with the MRD% of 3.7 %. All used models show a low MRD% values (< 4.0 %) for the calculated data indicating a good correlation of SD solubility data with the given mathematical models. By means of the van’t Hoff and Gibbs equations, the apparent thermodynamic quantities relative to SD dissolution and mixing processes, namely Gibbs energies, enthalpies, and entropies, were calculated and reported. Apparent dissolution quantities of SD were positive in all cases indicating endothermic and entropy-driven behaviors. A non-linear enthalpy–entropy relationship was observed for SD in the plot of SD dissolution enthalpy vs. Gibbs energy. Observed trend exhibits negative slope in the composition from neat EtOH to the mixture of 0.05 in mass fraction of MeCN indicating entropy-driving mechanism for this SD transfer process. Moreover, variant but positive slopes were found in the composition interval of 0.05 < w1 < 1.00 indicating enthalpy-driving mechanism for these SD transfer processes. Furthermore, the preferential solvation of SD by MeCN or EtOH was analyzed by using the inverse Kirkwood–Buff integrals. Thus, SD is preferentially solvated by EtOH molecules in EtOH-rich mixtures but preferentially solvated by MeCN in MeCN-rich mixtures. In this way, this research expands the literature investigations about the solubility of SD in some non-aqueous cosolvent mixtures conformed by MeCN and other alcohols.