Thermodynamic Modeling of Binary, Ternary and Quaternary Aqueous Solutions Containing Ba2+, Na+, Cl−, and SO42− By eNRTL Model