1. Introduction
Crystallization is one of the most important processes in industry. In many branches, such as pharmaceutical [
1,
2], agricultural [
3], and food industries [
4], demand for crystalline products is constantly growing. Apart from the quantity, other requirements that need to be met for such products are high purity, specific shape and size of crystals, as well as appropriate properties [
5].
Precipitation of sparingly soluble solids can be either a desired process or an undesired one, depending on the case. Phosphonic acid can be produced in a number of routes [
6] among which there are several ones involving precipitation. The occurrence of such crystalline product in this case will make the purification step more difficult. Sparingly, soluble solids can also deposit on the surface of equipment. In pipes it decreases their diameters by clogging [
7]. When the layer of calcium sulfate is present, the heat exchange is compromised. In membrane processes it can limit the recovery in reverse osmosis systems as well as in forward osmosis systems [
8,
9]. To end up with a well-determined average particle size of the crystals during, e.g., the processing of cane sugar, both the nucleation and growth need to be rigorously controlled [
3].
Calcium sulfate dihydrate is one of the hardly soluble salts and it is a side product, for example, in the desulfurization process of flue gases. Among various methods of SO
2 emission reduction the most widely used is desulfurization. The crystal habit of the produced gypsum is an important aspect of wet limestone scrubbing, since it significantly affects the by-product dewatering properties, and consequently the economics of the whole process [
5]. The crystallization of calcium sulfate has been intensely investigated since it is encountered in a large number of natural and industrial processes. Those include acidic wastewater treatment [
10], phosphoric acid production [
6,
11], wet flue gas desulfurization [
5,
12], reverse osmosis and desalination processes [
8,
9].
Gypsum formation, although it is described by number of authors, is still an issue in industry to which full explanation has not been delivered. The more information is extracted from various perspectives, the closer one is to understanding and properly describing this phenomenon. Influence of temperature and concentration is discussed in a variety of papers, but until now, there was no report on a simple and effective way of determining induction time based on a few points of known values of solution supersaturation. In industrial practice rotational speed is an equally important parameter. Proper determination of the influence of mixing on induction time of a given compound can help to define proper working conditions.
Nucleation is the initial stage in the crystallization process [
13] and is stochastic in its nature [
14]. It is described as the process of random generation of those small particles of the new phase that have the ability to irreversibly overgrow to macroscopic sizes [
15]. During this phenomenon, crystals of organic and inorganic substances are formed from the solution. The nucleation kinetics should be measured for better control of final crystal properties such as size distribution, purity, form, and habit.
The precipitation process is in the range of primary nucleation. Its kinetics can be measured either by adopting a metastable zone width or by induction time method [
16,
17]. Despite the fact that the metastable zone width method is less labor intensive, the induction time method seems to have higher accuracy and is easier to analyze [
17]. The time which elapses between achievement of constant supersaturation and first observed change of phase (appearance of the first crystals) is different for each system and is called the induction time (
tin) [
18]. Therefore, it is “a measure of the ability of a supersaturated solution to remain in the state of metastability” [
19]. More precisely, one can define it as the sum of the time of critical nucleus formation,
tn, and its growth to detectable size—observed first as changing properties of solution,
tg [
13,
16]:
Currently, no exact theory regarding nucleation exists, but more and more studies confirming the presence of two stages of the process could be found (Söhnel and Mullin, 1988) [
20,
21].
Crystal growth to detectable size can be defined in various ways, such as the first “appearance” of crystals, the onset of a change of some physical properties of solution, etc. [
21]. The subjective assessment of the emergence of a new phase in solution is not very exact. Therefore, a method of accurate measurement of the nucleation phenomenon is necessary if precise data about induction period is needed. Detailed knowledge of the crystallization kinetics allows for full control of the process and achievement of the final product with desired properties: size of crystals, humidity, and purity. The detection of the new phase is indicated by the appearance of crystals as an effect of increasing size of crystal nuclei. The determination is possible as long as it exceeds a reasonable period of time. This condition arises from the limitations of many of measuring devices, which tend to fail when
tin is smaller than 5 s [
22,
23]. For successful measurement of short induction time, we need very fast and sensitive devices which can measure changes in solution properties before visual confirmation.
The induction time depends on many factors. It is considerably influenced not only by the initial and process conditions such as supersaturation (concentration) and temperature, but also, among others, by mixing condition (stirring rate), type and composition of solvent [
16,
18,
19,
24].
One of the most important factors influencing induction time are ion collisions which create crystal nuclei [
25]. A larger number of ions (increase in supersaturation level) in the solution increases the probability of nucleation. “The induction period decreases exponentially with supersaturation, which suggests that the nucleation rate increases exponentially” [
26].
Influence of supersaturation was investigated by a number of researchers [
11,
27]. Ziegenheim et al. [
27] investigated precipitation kinetics of gypsum over quite a large range of concentrations: from 0.04 M to 0.2 M. Their research was mainly focused on recognition of kinetics, but they found that over this range of concentrations the induction time decreased by 2 orders of magnitude as concentration increased. Worthwhile to notice is that tested volumes were very small (50 cm
3) and investigated concentrations also were relatively low. Influence of volume on kinetic conditions is also reported in the literature [
14], and the volumes used in the described paper reflect purely laboratory conditions.
Temperature impact was also investigated by number of researchers. The influence of temperature on induction time of calcium sulfate was investigated in the work of Hoang et al. [
7], which confirms that higher temperature favors a significant decrease in induction period. Klepetsanis et al. [
28] investigated precipitation of calcium sulfate in the temperature range of 25 °C to 80 °C but the research was focused on kinetics of reaction. In this case, induction time was only used as a mean to determine precipitation kinetics. In addition, Alimi et al. [
18] investigated dependence of induction time on temperature and supersaturation. From this dependence, it was possible to distinguish between the homogeneous and heterogeneous nucleation mechanisms. It was also observed that the induction period depends highly on the lattice cation/anion molar ratio. It is more important for the lower supersaturations.
In the case of mixing conditions, there is no information on such investigations. Such research is reported in this paper.
1.1. Induction Time Calculation
Mathematical considerations on induction time are based on a classical nucleation theory (CNT). Recently, many studies [
29,
30] have been carried out to measure gypsum induction time and described by CNT. On the other hand, there are papers confirming that gypsum precipitation is a multi-stage process and can be described by “Non-Classical Nucleation Theory” [
31,
32], thus making the CNT unsuitable for describing complex pathways in crystallization process [
33]. There are articles in which research was conducted on the precipitation of gypsum from aqueous solutions, observing the formation of the other two forms of gypsum, namely anhydrous and hemihydrate [
34]. Nevertheless, for assessment of induction time—the first observed change of phase solutions—the CNT is suitable. The induction time is defined as the time which elapsed between achievement of supersaturation and the formation of critical nuclei. Thus, one can conclude that it depends on the nucleation rate [
35]. Induction time is considered inversely proportional to nucleation rate (
J) [
19] and can be expressed with the following formula:
while the nucleation rate (rate of nucleus formation) is given by an Arrhenius equation:
where:
—Boltzmann constant (J/K);
—Pre-exponential factor (-);
—Absolute temperature (K);
—Gibbs free energy for the formation of the crystalline phase (kJ/mol), which can be defined as:
where:
—Supersaturation ratio (-);
—Molar crystal volume (m3/mol);
—Interfacial energy of the nucleus-solution system (J/m2);
—Shape factor (-).
After substitution of Equation (4) to Equation (3) and the result to Equation (1) one can obtain following formula:
This equation describes the curve of the induction time as a function of supersaturation but only at constant temperature. The slope of this straight line (
B) is defined as:
Of course, there are a lot of approaches to describe induction time. One of them is a model proposed by Söhnel [
21]. It is based on interfacial tension between the crystal and aqueous solution as a fundamental parameter in modeling of both nucleation and crystal growth [
18]. Söhnel proposed the following equation:
where
A is constant and
B is slope, given by:
Equation (8) has the same form as the equation of Gibbs free energy. It is created by substituting the equation for Boltzmann constant (Equation (9)) from Equation (7) with changed logarithm base. Boltzmann constant is described as:
where:
R—Universal gas constant (J/mol·K);
NA—Avogadro’s number (1/mol).
Calculation of induction time requires knowledge of certain values, namely: solution supersaturation and constant in the formula proposed by Söhnel [
22]. This equation is only valid for measurements carried out at constant temperature and with varying concentration. As a consequence, all induction time calculations in this paper were carried out based on results obtained at various concentrations.
1.2. Calculation of Solution Supersaturation
For each trial the solution supersaturation was calculated. In order to do this, Equation (10) was used [
5,
18,
23,
28]:
where
and
are the activity (mol/m
3), concentration in molality (mol/kg
solvent) and activity coefficient for ingredient
, respectively, and
is the activity of water, which can be calculated from the osmotic coefficient.
is the solubility product of gypsum.
There are many equations describing activity coefficient, but two of them are the most commonly used. They are Debye–Hückel and Bromley (modified Debye–Hückel) equation.
Both equations are used for solutions of strong electrolytes but there is a small but significant difference. Namely, the first equation can be used only for small values of ionic strength while the second one provides approximation to the activity coefficient variation over the entire range of ionic strength values. It is a linear function described by the following expression [
36]:
where:
—Charge of ion (cation or anion) (-);
—Constants (-);
—Ionic strength (-);
—Individual ions values (-).
The values used in Equation (11) were taken from the work of Bromley [
36] and are presented in
Table 1. These values allow to determine the supersaturation of the system in the case of gypsum precipitation.
The algorithm for supersaturation calculation is as follows. First, calculation of the ions molality in solution is performed using the equation:
where:
—Amount of moles of compound in solution (mole);
—Volume of solution in reactor (0.3 dm3);
—Mass of solvent (kg);
—Concentration of compound (mol/dm3);
Then the ionic strength of solution for various concentrations is calculated using the following formula:
where:
—Number of ions in solution;
—Charge of ion;
—Molality concentration of ion.
Next, calculation of the activity coefficient should be performed. Because of high ionic strength of investigated solutions, calculations were made based on Bromley Equation (11) and values from
Table 1. Moreover, in the formula for the activity coefficient, the absolute values of ion charges and stoichiometric coefficients of the precipitated salt were substituted. Calculations also required knowledge of individual ion values per salt. It was calculated based on the following formula:
After substitution of the values from
Table 1 following result was obtained:
In the next step, the thermodynamic solubility product was calculated using empirical equation [
5,
18,
23]:
The solubility product depends on temperature and thus calculations were made for the whole investigated temperature range. The results are presented in
Table 2. As one can notice, its value is changing with temperature.
Determination of water activity value is a necessary step to calculate supersaturation. If the concentration of a soluble compound rises, the activity coefficient of water decreases below a value of 1. Activity of water should be determined using Equation (16). This correlation contains osmotic coefficient (
ϕ) which should be determined using Bromley theory [
36]:
where:
M—Molar mass (kg/mol);
v—Stoichiometric coefficient of ion of dissolved compounds (-);
c—Concentration of compound (mol/dm3).
Osmotic coefficient was calculated according to following formula [
36]:
where:
5. Conclusions
The presented paper describes precipitation of calcium sulfate dihydrate in a thermostated reactor. Influence of concentration, temperature and rotational speed is reported.
Induction time was measured by means of two methods: turbidimetric and visual. Both methods gave very similar results of induction time. However, the first one proved to be more precise. Moreover, both methods have some limitations. In the case of very short induction times, it would not be possible to perform sufficiently precise measurements using either of them. Classical nucleation theory has limitations regarding gypsum because it forms different structures depending on hydration degree. This work concerns only gypsum production regardless of number of water of molecules.
Measurement of induction time with a turbidimetric probe proved to be an effective and simple solution for investigated range of parameters. The turbidity change, that is also confirmed by visual method, coincides with the values obtained from measurements using the turbidity.
In this work a set of experiments concerning measurement of gypsum induction period for nucleation as a function of temperature, rotational speed and concentration was presented. During measurements a new parameter was tested, which is rotational speed. In each of these three cases, the increase in the value of the tested parameter resulted in faster precipitation of calcium sulfate dihydrate and thus shorter induction time. Better control of each parameter would allow to obtain a more desirable form of calcium sulfate dihydrate and improving its properties.
What is more, the calculation of the calcium sulfate–sodium chloride system in water was performed on the basis of available formulas and knowledge of the concentration of feed solutions.
In this paper, the induction time was also determined based on the empirical equation using a simple, yet previously not reported method. The measured values are close to the calculated values, which proves the correctness of the method and application of the presented algorithm for determination of the induction time for a given temperature and desired supersaturation level without the need to perform experiments. On the other hand, the measurement of turbidity together with the use of the described equation for the time of induction could be used to determine the supersaturation of solutions in which the classic methods are not effective.
The resulting gypsum sediment was filtered and analyzed using a light microscope. In the obtained sediment samples, two different crystalline forms of gypsum could be observed, namely elongated and twining crystals. Both forms are characteristic for calcium sulfate dihydrate. From analysis of obtained photographs one can conclude, that if a specific habit of obtained crystals in desired, then, depending on which process parameter(s) is modified, e.g., concentration, temperature or rotational speed, the specific requirement can be met.