Scaling of Na2SO4, Na2CO3 and their double salts has long been recognized as a problem in black liquor evaporators. We have used published data for Na-CO3-SO4-H2O system to develop a model for predicting the solubility and solid phase composition within the double salt precipitation region. The model is based on a thermodynamic definition of a solubility constant with Pitzer ionic activity coefficients in the solution. It is capable of calculating the solubility and the equilibrium solid composition in Na2SO4-Na2CO3 aqueous solutions in the temperature range of 100-115°C and 150°C and within the precipitation region of a solid solution consisting of burkeite (2Na2SO4Na2CO3) and sodium carbonate salts. The model construction procedure and its possible extension to other Na2SO4-Na2CO3 precipitates will be addressed.