Rackett density Equation or Activity Coeffcients for NaOH - H2O system

[defencedog]

@ https://github.com/reaktoro/reaktoro/issues/275

The above issue is also relevant in this because interaction parameters effect computed densities especially for strong interacting species.

I have mentioned my code here: https://github.com/reaktoro/reaktoro/issues/275#issue-1510423916

Now the density is always reported higher than published data. I am in chloralkali industry which produces caustic. For comparison I am using 32w.t%, 2barA, using both phreeqc & thermofun

code_py.zip

phreeqc

50C 1531.46
70C 1515.08
90C  1507.4

thermofun

50C 1473.27
70C 1454.54
90C  1438.4

Actual (p-13 https://www.griswoldandwagner.com/yabbfiles/Attachments/Caustic_soda_Dow.pdf)

50C 1330
70C 1317
90C  1303