This study investigates the solubility behavior of Naproxen (NAP) in binary solvent mixtures of 1-propanol (1-PrOH) and 2-propanol (2-PrOH) with ethylene glycol (EG) across a range of temperatures. The solubility of NAP was experimentally determined at five different temperatures (293.15 to 313.15 K), and the data were correlated using various thermodynamic models, including the van't Hoff, Jouyban-Acree, modified Wilson, mixture response surface, Jouyban-Acree-van't Hoff. The results demonstrated that NAP's solubility increases with temperature in both solvent systems. Notably, NAP exhibited higher solubility in mixtures with 1-PrOH compared to 2-PrOH, despite the lower polarity of 2-PrOH. This unexpected trend is attributed to the distinct molecular interactions, including hydrogen bonding, influenced by the structural differences between 1-PrOH and 2-PrOH. The X-ray diffraction analysis confirmed that no polymorphic transformation occurred in NAP during dissolution, maintaining its crystalline structure. The solubility data were well-correlated by the applied models, with overall MRDs% (mean relative deviation percentage) below 6.1.