Adapting molar data (without density) for molal models

Theoretical geochemical models for electrolyte solutions based on classical thermodynamic principles rely largely upon molal concentrations as input because molality (wt/wt) is independent of temperature and pressure. On the other hand, there are countless studies in the literature where concentrations are expressed as molarity (wt/vol) because these units are more easily measured. To convert from molarity to molality requires an estimate of solution density. Unfortunately, in many, if not most, cases where molarity is the concentration of choice, solution densities are not measured. For concentrated brines such as seawater or even more dense brines, the difference between molarity and molality is significant. Without knowledge of density, these brinish, molar-based studies are closed to theoretical electrolyte solution models. The objective of this paper is to present an algorithm that can accurately calculate the density of molar-based solutions, and, as a consequence, molality.The algorithm consist of molar inputs into a molal-based model that can calculate density (FREZCHEM). The algorithm uses an iterative process for calculating absolute salinity (SA)(SA), density (ρ)(ρ), and the conversion factor (CF)(CF) for molarity to molality. Three cases were examined ranging in density from 1.023 to 1.203 kg(soln.)/l. In all three cases, the SASA, ρρ, and CFCF values converged to within 1ppm by nine iterations. In all three cases, the calculated densities agreed with experimental measurements to within ±0.1%±0.1%. This algorithm opens a large literature based on molar concentrations to exploration with theoretical models based on molal concentrations and classical thermodynamic principles.