Next Article in Journal
High-Speed and High-Power 940 nm Flip-Chip VCSEL Array for LiDAR Application
Previous Article in Journal
Strain Localization Modes within Single Crystals Using Finite Deformation Strain Gradient Crystal Plasticity
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Influence of Supersaturation, Temperature and Rotational Speed on Induction Time of Calcium Sulfate Crystallization

1
Department of Chemical Engineering and Process Design, Faculty of Chemistry, Silesian University of Technology, ul. Strzody 7, 44-100 Gliwice, Poland
2
Łukasiewicz Research Network, Institute of Non-Ferrous Metals, Centre of Hydroelectrometallurgy, ul. Sowińskiego 5, 44-100 Gliwice, Poland
*
Author to whom correspondence should be addressed.
Crystals 2021, 11(10), 1236; https://doi.org/10.3390/cryst11101236
Submission received: 7 September 2021 / Revised: 23 September 2021 / Accepted: 12 October 2021 / Published: 14 October 2021
(This article belongs to the Section Crystal Engineering)

Abstract

:
Calcium sulfate is a very important product of precipitation reactions in various branches of industry. The most common applications include building materials, impression materials in dentistry, immobilizing casts and an inactive ingredient of tablet excipients. It is also used as a drying agent and color glaze. In this paper, influence of various conditions, i.e., substrates concentration (supersaturation), temperature and rotational speed on induction time was investigated. The range of investigated parameters was 0.15–0.35 M for concentration with step of 0.05 M, 298–323 K with step of 5 K for temperature, and 1.67–8.33 s−1 for rotational speed with change every 1.67 s−1. The experiments were performed using visual and turbidimetric methods. It was proven that increase in all three investigated parameters resulted in shorter induction time. The turbidimetric method was found to be more precise, but still insufficient for concentrations below 0.15 M. Moreover, as a result of data analysis, a simple method of induction time calculation was proposed when arbitrary measurements in various concentrations are known.

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 SO2 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]:
t i n = t n + t g  
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 cm3) 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:
t i n 1 J
while the nucleation rate (rate of nucleus formation) is given by an Arrhenius equation:
J = A · exp ( Δ G k T )
where:
k —Boltzmann constant (J/K);
A —Pre-exponential factor (-);
T —Absolute temperature (K);
Δ G —Gibbs free energy for the formation of the crystalline phase (kJ/mol), which can be defined as:
Δ G = β δ 3 V m 2 k 2 T 2 ( ln S ) 2
where:
S —Supersaturation ratio (-);
V m —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:
l n ( t i n ) = K + β δ 3 V m 2 k 3 T 3 ( ln S ) 2
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:
B = β δ 3 V m 2 k 3 T 3
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:
log ( t i n ) = A + B ( l o g S ) 2
where A is constant and B is slope, given by:
B = β δ 3 V m 2 N A ( 2.3 R T ) 3  
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:
k = R N A
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]:
S = a C a 2 + ·   a S O 4 2 a H 2 O 2 K s p = m C a 2 + ·   m S O 4 2 · γ C a 2 + ·   γ S O 4 2 a H 2 O 2 K s p
where a i ,   m i and γ i are the activity (mol/m3), concentration in molality (mol/kgsolvent) and activity coefficient for ingredient i , respectively, and a H 2 O 2 is the activity of water, which can be calculated from the osmotic coefficient. K s p 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]:
log γ i = A γ · | Z + Z | · I 0.5 1 + ρ I 1 / 2 + ( 0.06 + 0.06 B ) | Z + Z | I ( 1 + 1.5 | Z + Z | I ) 2 + B I + C I 2
where:
z + , z —Charge of ion (cation or anion) (-);
A γ ,   ρ ,   C —Constants (-);
I —Ionic strength (-);
B —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:
m i = n i m s o l v e n t = V · C i m s o l v e n t
where:
n i —Amount of moles of compound in solution (mole);
V —Volume of solution in reactor (0.3 dm3);
m s o l v e n t —Mass of solvent (kg);
C i —Concentration of compound (mol/dm3);
Then the ionic strength of solution for various concentrations is calculated using the following formula:
I = 1 2 i = 1 n m i · z i 2
where:
n —Number of ions in solution;
z i —Charge of ion;
m i —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:
B = B C a 2 + + B S O 4 2 + σ C a 2 + · σ S O 4 2
After substitution of the values from Table 1 following result was obtained:
B = 0.0347 + 0 + 0.119 · ( 0.4 ) = 0.0102
In the next step, the thermodynamic solubility product was calculated using empirical equation [5,18,23]:
log K s p = 390.9619 152.6246 log T 12545.62 T + 0.0818493 T
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]:
( a H 2 O ) = ϕ M H 2 O v c i
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]:
1 ϕ = 2.303 A | z + z | I 0.5 3 σ ( ρ I ) 2.303 ( 0.06 + 0.6 · B ) I 2 ψ ( a I ) 2.303 B I 2
where:
σ ( ρ I ) = 3 ( ρ I ) 3 [ 1 + ρ I 1 1 + ρ I 2 l n ( 1 + ρ I ) ]
ψ ( a I ) = 2 a I [ 1 + 2 a I ( 1 + a I ) 2 l n ( 1 + a I ) a I ]
a = 1.5 | Z + Z |

2. Materials and Methods

2.1. Solutions Preparation

Two reagents were used to prepare appropriate solutions: calcium chloride dihydrate (VWR International) and sodium sulfate (VWR International). Both of them were of analytical purity, kept under anhydrous conditions. The calculated amount of reagents was weighed on an analytical balance (accuracy +/− 0.0001 g) in a balance room and transferred quantitatively to flasks. All solutions were prepared in 0.5 dm3 flasks by dissolving salts in deionized water (Hydrolab, Poland). Prepared solutions were left in sealed glass flasks and stored in a dark and dry container. Both solutions were prepared in a number of concentrations equal to 0.15, 0.20, 0.25, 0.30 and 0.35 mol/dm3.
The preparation procedure for both visual and turbidimetric measurement was the same, although measurement was different. Equal volumes of 150 cm3 of both solutions were placed in separate vessels and agitated using magnetic stirrer at constant temperature for about 30 min to achieve thermal equilibrium. The solutions were equimolar. The mixing of two solutions (Na2SO4 and CaCl2) of a given concentration took place directly in the jacked reactor (working volume: 0.3 dm3) kept at a constant set temperature using thermostat (Huber Ministat).

2.2. Measurement Procedure

Induction time was measured using two methods: turbidimetric and visual. Each measurement was repeated three times and the obtained data were averaged.
Turbidity measurements were performed using a special equipment consisting of: digital meter (MULTI 3510 IDS) for digital IDS sensor (VISO TURBO 900-P) connected to a computer via dedicated software for collecting data (MULTI LAB USER). Such set up allowed to gather results of measurement directly in a spreadsheet file in very short time intervals (1 s). Moreover, during each measurement the turbidity probe was always immersed in the solution at the same height and at the same angle. This was necessary to obtain repeatable results.
Turbidimetric method is based on measuring the critical nuclei formation with detectable size. Appearance of visible crystals (with size larger than critical size) causes changes in intensity of transmitted light along the detector [27].
In the case of visual assessment of the induction time videos were recorded by digital camera (LUMIX Quad Indicator) illustrating the turbidity of the solution over time. Based on the recorded film, the change in the solution turbidity could be determined.
During tests the influence of the following parameters on induction time were examined: concentration, temperature, and rotational speed. The following ranges were selected:
  • Concentration: 0.15–0.35 M (0.05 M interval).
  • Temperature: 298–323 K (5 K interval).
  • Rotational speed: 1.67–8.33 s−1 (1.67 s−1 interval).
Each of above-mentioned parameters were examined under constant conditions. When the determining the influence of temperature the supersaturation and rotational speed were at constant level. The investigated ranges were selected in such a way that the resulting induction time was long enough to be observable and short enough to be reasonable.

3. Results and Discussion

3.1. Influence of Concentration

The first investigated parameter was concentration. All measurements were made at constant temperature 298 K and with mixing rate of 1.67 s−1.
The exemplary turbidity measurement results are presented in Figure 1. As it could be observed, the higher the concentration the shorter induction time. In the case of concentration of 0.15 M the induction time was found to be equal to 160 s, for 0.20 M it was 92 s, for 0.25 M—69 s, for 0.30 M—49 s, and for 0.35 M—28 s. Change in concentration between two neighboring values was constant and equal to 0.05 M, but difference between induction time for subsequent values of concentration becomes smaller when concentration increases. What is also quite interesting, the point at which precipitation starts is also more clearly marked with increasing concentration.
In Table 3 the results of visual observations are presented. For an easier comparison between investigated concentrations a time-lapse from exactly the same time points are presented. As one can observe, with increasing concentration the turbidity of solution becomes visible faster. The pink color is a result of using colored cooling agent in the thermostat. The results obtained using both methods are similar.
In Figure 2 induction time is presented as a function of concentration. The effect is clearly visible: time decreases with increasing concentration and so the slope is negative. This result is as expected. The reason is that more concentrated solutions have a higher supersaturation and thus a higher driving force for nucleation to occur. In most research obtained points lay along a straight line [37,38]. In the case of the presented measurements, the result for lowest concentration is visibly out of range. The reason can be that the smallest concentration and thus supersaturation results in smaller concentration of ions and, as a consequence, the number of successful collisions is smaller. Therefore, determination of the exact moment at which turbidity increases can be hard. Excluding this one point from the chart one can obtain an almost exactly straight line (Figure 3) with Pearson’s correlation coefficient almost equal to 1. It means that the relationship between concentration (supersaturation) in the range between 0.2 and 0.35 M and induction time is linear.
After analyzing the results, the concentration of 0.25 M was chosen for further investigations.

3.2. Influence of Temperature

The second analyzed parameter was temperature. In this case, concentration was constant and equal to 0.25 M. Additionally, mixing rate was constant and equal 1.67 s−1. Figure 4 presents the results of turbidimetric measurements. As one can notice, a change in induction time with change of temperature can be observed. The longest induction time (49 s) is for temperature 298 K, then it decreases to 41 s for 303 K, 35 s for 308 K, 26 s for 313 K, 18 s for 318 K down to 15 s for 323 K.
Results from the visual method confirm the correctness of the turbidimetric method. Again, for comparison, time-lapse photos taken at the same time point are presented in Table 4. The influence of temperature on induction time is again indisputable. However, it is smoother than in the case of concentration. What is interesting, the result for measurements taken at 298 K is slightly different from the measurements taken previously. There could be many reasons for such discrepancy. It could be caused by some mistakes during preparation of an initial solution, despite keeping extreme caution. In addition, some microimpurities could be present in the reactor. However, as crystallization is a challenging process such differences in obtained results can happen due to many reasons.
As in the previous case, both methods were comparable. A noticeable difference in the turbidity of solution in Table 5 for the 30 s time-lapse exist between 303 K and 308 K. A very similar result (35 s) was also obtained from measurements using the turbidity probe.
In Figure 5 induction time as a function of temperature is presented. In comparison to Figure 2 where induction time is a function of concentration, this change is more uniform, which is also visible in the presented trend line and its equation. The slope is again negative, which gives increasing influence of temperature as its value is higher. The highest investigated temperature was 323 K as above this value the product of precipitation is no longer calcium sulfate dihydrate (gypsum) only, but calcium sulfate semihydrate is also obtained. Its presence could distort the results.
After analyzing the obtained results, a temperature of 298 K was chosen for further investigations.

3.3. Influence of Rotational Speed

The third examined parameter was rotational speed. In this case, concentration and temperature were constant and equal to 0.25 M and 298 K, respectively. In Figure 6 results of turbidimetric measurements are presented. Again, the highest value of induction time is reported for the lowest rotational speed (67 s). For rotational speeds of 3.33 s−1, 5.00 s−1, 6.67 s−1 and 8.33 s−1 the resulting induction times are 56 s, 53 s, 45 s and 41 s, respectively. For the highest rotational speed, the measuring points are very scattered which results from high turbulence.
The averaged values are illustrated in Figure 7. As expected, an increase in agitation speed causes a decrease in the induction time of calcium sulfate dihydrate crystals. However, the differences characterizing various rotational speeds are much smaller than that observed in the case of concentration and temperature dependences. The slope of the function is, as in previous cases, negative and the overall fit is at good level. The higher the rotational speed is the better the dispersion of molecules of the two species and thus the reaction can take place more effectively.
However, time-lapse photos in Table 5 taken after 60 s indicate a small change in the turbidity of the solution. Almost in all cases the solution looked similar except for 1.67 s−1 where change is very subtle and 8.33 s−1 where the solution is visibly more turbid. This indicates that the visual method is not precise enough to determine the induction time in this case. Small changes in the measured times for the corresponding speeds make it difficult to subjectively assess the moment at which crystallization occurs.
As by definition induction time is inversely proportional to nucleation rate, one can deduce that nucleation rate increases with decreasing induction time.

3.4. Crystal Shape Analysis

The dried gypsum precipitate was subjected to microscopic analysis. The CaSO4 crystals appear to look like a fine powder. In Figure 8 microscopic photographs of calcium sulfate dihydrate for varying concentration are presented. For lowest concentration (Figure 8a) crystals are the largest. Additionally, both crystal forms reported for gypsum, namely needle-like and dove tails, can be noticed. As the concentration increases (Figure 8b–e) the observed crystals become smaller. The reason for such situation is rapid nucleation. For low values of concentration and thus supersaturation, the induction time is relatively longer, and the ions can arrange in a more orderly manner. In the case of higher concentrations, the structure of precipitate changed. This does not mean that the structure of gypsum was changed, but it means that at high concentrations (high supersaturation), spontaneous precipitation occurred, which resulted in the very fast precipitation out of the solution with very short induction time and may have resulted in inhibition of the crystals’ growth over time [39]. This resulted in the formation of a large number of crystals, but of small size.
A long induction period permits the growth of larger crystals (Figure 8a) while a shorter induction period favors creation of smaller crystals (Figure 8e).
A quite different relationship can be noticed by analyzing photographs of crystals synthesized at different temperatures presented in Figure 9. As one can notice, an increase in temperature caused the precipitation of larger and larger gypsum crystals. Moreover, it can be seen that at higher temperatures (shorter induction times) a large number of polycrystals is formed. All photos show two crystalline forms of calcium sulfate dihydrate. The speed at which crystals precipitate affects their size.
For all investigated rotational speeds both forms of gypsum can also be noticed, as presented in Figure 10a–e. However, the higher the rotational speed the more aggregated crystals appear. In addition, the ratio of dovetail to needle-like shaped crystals is shifted toward dovetail shaped crystals. Additionally, a slight increase in size of particles can be observed but this is not as clear of an increase as in previous cases.
An obvious conclusion comes to mind, which is that if a specific habit of obtained crystals is desired, then depending on which process parameter(s) is modified, e.g., concentration, temperature or rotational speed, the specific requirement can be met.

4. Calculations

4.1. Supersaturation Calculation

Changes in supersaturation, ionic strength, and activity coefficient as a function of concentration and temperature are presented in Table 6 and Table 7. As one can notice, ionic strength in the case of both varying concentration and temperature is above the applicability range of classic Debye–Hückel equation, hence the use of its extended version proposed by Bromley is used. The investigated range of concentrations results in supersaturation level of 4.264 to 12.505. Usually for sparingly soluble substances crystallized by means of chemical reaction, supersaturation level is very high (above 100). However, this is true for low ionic strength of solutions. In the case of high ionic strength, such as for calcium sulfate, such high supersaturation cannot be obtained, and much lower values are accepted. Change in supersaturation for constant temperature and varying concentration is three-fold, whereas for constant concentration and varying temperature it is very small. In the case of changing rotational speed both concentration and temperature are constant and therefore supersaturation level will be also constant.

4.2. Calculation of Induction Time

As mentioned earlier, Equation (7) can be used to calculate induction time but only when its influence on concentration is considered. In the case of temperature impact, supersaturation as well as variable B depends on this parameter, while for the influence of mixing the temperature and supersaturation are constant. The above limitations make it impossible to calculate the induction times in all presented cases—apart from varying concentration one is reliant on measurements.
In the literature there are very few attempts made to calculate induction time. Usually, they involve many complicated steps or are not given in a straightforward way. In this paper authors present a simple way to obtain induction time of calcium sulfate at a given temperature for chosen concentration within selected range. Calculations were carried out according to the following algorithm.
Based on Equation (8) and values in Table 8 variable B was calculated. As a temperature-dependent parameter, it greatly influences the slope on the chart describing induction time as a function of supersaturation.
After substituting values from Table 8 and the selected temperature into the Equation (7), one can obtain value for parameter B. Results for the whole investigated temperature range are presented in Table 9.
The curve of induction time as a function of supersaturation is described by Equation (7), in which the coefficient B comes from Table 9 while the constant A was determined using the least squares method according to the algorithm presented in Figure 11.
After plotting log(tin) vs. 1/(logS)2 for various concentrations and determining the fitting curve, one can also read the sought values of the constants A and B. In this way one can determine the induction time for various concentrations at a given temperature. In Figure 11 two cases are presented: one with all investigated values of supersaturation and one without the smallest value of investigated supersaturations, as this value has an influence on the equation for direct calculation of induction time. As one can notice, in Figure 11a the Pearson’s coefficient is quite high, but the one presented in Figure 11b is closer to unity and thus the fitting is better. Apart from this paper there exist many data on gypsum, thus if there is a need to know the induction time in a given temperature for an arbitrary concentration, then using the described method is simple and definitely less labor intensive than performing experiments. The values read from the chart are in good agreement to values obtained using the least square method.
Calculated induction times were compared with the ones obtained from the measurements in Table 10. Moreover, a standard deviation was calculated for each case. As one can see the deviations are of reasonable value—only for one concentration deviation is high. The used method of the least squares is a simple one giving comparable results as well as the equations used proved to be correct. The reason for values of relative error is probably due to high sensitivity of induction time measurement to external conditions. Despite the fact that authors made every effort to perform each trial under exactly the same conditions including use of the same initial solutions, for each concentration one result was significantly different from the other two. This can cause a slight variation of received averaged values and some grater fluctuations in final calculations.

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.

Author Contributions

Conceptualization, A.M.; methodology, A.M. and S.O.; validation, A.M. and S.O.; formal analysis, A.M., S.O.; investigation, S.O.; resources, A.M.; writing—original draft preparation, S.O.; writing—review and editing, A.M. and M.D.; visualization, S.O.; supervision, A.M. and M.D.; project administration, A.M.; funding acquisition, A.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the Polish National Science Centre (NCN) under Grant No. 2016/21/D/ST8/01714.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

Symbols
aactivity (mol/dm3)
Aconstant (-)
Ayconstant (-)
Bconstant (-) individual ion values (-)
Cconstant (-)
Cmmolar concentration (mol/dm3)
fanephelometric signal (FNU)
Ggrowth rate (s)
ΔGGibbs free energy (J)
Iionic strength (mol/dm3)
Jnucleation rate (1/m3s)
kBoltzmann constant
kaarea shape factor (-)
K S P solubility product (-)
nnumber of moles (-)
NAAvogadro’s number (1/mol)
mmass (kg) concentration in molality (mol/kgsolvent)
Mmolar mass (kg/kmol)
Runiversal gas constant (J/K·mol)
Srelative supersaturation (-)
tindinduction time (s)
Ttemperature (K)
Vvolume (m3)
Vmcrystal molar volume (m3/mol)
zion charge (-)
Greek letters:
α degree of conversion (-)
β shape factor (-)
γ activity coefficient (-)
γ s shear stress (Pa)
δ constant (-) surface energy for gypsum (J/m2)
ϑ crystal grow rate (m/s)
μ potential (J/mol)
ρ constant (-)
σ a relative degree of supersaturation (-)
σ c dimensionless coefficient being a measure of interaction between neighboring
τ time (s)
Φ osmotic coefficient (-)
φ volume fraction of solid in liquid (-)
ν stoichiometric coefficient of ion of dissolved compound
Subscripts:
avaverage
ii-compound
ngnano gram
s−1revolution per minute
solvsolvent
+positive ion charge
negative ion charge
H2Owater

References

  1. Gao, Z.; Rohani, S.; Gong, J.; Wang, J. Recent Developments in the Crystallization Process: Toward the Pharmaceutical Industry. Engineering 2017, 3, 343–353. [Google Scholar] [CrossRef]
  2. Shekunov, B.Y.; York, P. Crystallization Processes in Pharmaceutical Technology and Drug Delivery Design. J. Cryst. Growth 2000, 211, 122–136. [Google Scholar] [CrossRef]
  3. Kikuchi, Y.; Suginobe, R.; Kanzaki, Y.; Ohara, S.; Hirao, M. Integrated Modeling of Sugarcane Cultivation and Ethanol Fermentation from Agriculture and Engineering Perspectives. In Computer Aided Chemical Engineering; Elsevier: Amsterdam, The Netherlands, 2012; Volume 30, pp. 182–186. [Google Scholar]
  4. Deora, N.S.; Misra, N.N.; Deswal, A.; Mishra, H.N.; Cullen, P.J.; Tiwari, B.K. Ultrasound for Improved Crystallisation in Food Processing. Food Eng. Rev. 2013, 5, 36–44. [Google Scholar] [CrossRef] [Green Version]
  5. Lancia, A.; Musmarra, D.; Prisciandaro, M. Measuring Induction Period for Calcium Sulfate Dihydrate Precipitation. AIChE J. 1999, 45, 390–397. [Google Scholar] [CrossRef]
  6. Sevrain, C.M.; Berchel, M.; Couthon, H.; Jaffrès, P.-A. Phosphonic Acid: Preparation and Applications. Beilstein J. Org. Chem. 2017, 13, 2186–2213. [Google Scholar] [CrossRef] [PubMed]
  7. Hoang, T.A.; Ang, H.M.; Rohl, A.L. Effects of Temperature on the Scaling of Calcium Sulphate in Pipes. Powder Technol. 2007, 179, 31–37. [Google Scholar] [CrossRef]
  8. Mi, B.; Elimelech, M. Gypsum Scaling and Cleaning in Forward Osmosis: Measurements and Mechanisms. Environ. Sci. Technol. 2010, 44, 2022–2028. [Google Scholar] [CrossRef]
  9. Pomerantz, N.; Ladizhansky, Y.; Korin, E.; Waisman, M.; Daltrophe, N.; Gilron, J. Prevention of Scaling of Reverse Osmosis Membranes by “Zeroing” the Elapsed Nucleation Time. Part I. Calcium Sulfate. Ind. Eng. Chem. Res. 2006, 45, 2008–2016. [Google Scholar] [CrossRef]
  10. Deng, L.; Zhang, Y.; Chen, F.; Cao, S.; You, S.; Liu, Y.; Zhang, Y. Reactive Crystallization of Calcium Sulfate Dihydrate from Acidic Wastewater and Lime. Chin. J. Chem. Eng. 2013, 21, 1303–1312. [Google Scholar] [CrossRef]
  11. Abdel-Aal, E.A.; Rashad, M.M.; El-Shall, H. Crystallization of Calcium Sulfate Dihydrate at Different Supersaturation Ratios and Different Free Sulfate Concentrations. Cryst. Res. Technol. 2004, 39, 313–321. [Google Scholar] [CrossRef]
  12. Cheng, T.; Zhou, X.; Yang, L.; Sun, Z.; Wu, H. Emission Characteristics of Soluble Ions in Fine Particulates in Limestone–Gypsum Wet Flue Gas Desulfurization System. Energy Fuels 2020, 34, 3836–3842. [Google Scholar] [CrossRef]
  13. Mullin, J. Crystallization; Butterworths: New York, NY, USA, 1961. [Google Scholar]
  14. Kubota, N. Effect of Sample Volume on Metastable Zone Width and Induction Time. J. Cryst. Growth 2012, 345, 27–33. [Google Scholar] [CrossRef]
  15. Knezic, D.; Zaccaro, J.; Myerson, A.S. Nucleation Induction Time in Levitated Droplets. J. Phys. Chem. B 2004, 108, 10672–10677. [Google Scholar] [CrossRef]
  16. Karthika, S.; Radhakrishnan, T.K.; Kalaichelvi, P. Measurement of Nucleation Rate of Ibuprofen in Ionic Liquid Using Induction Time Method. J. Cryst. Growth 2019, 521, 55–59. [Google Scholar] [CrossRef]
  17. Kulkarni, S.A.; Kadam, S.S.; Meekes, H.; Stankiewicz, A.I.; ter Horst, J.H. Crystal Nucleation Kinetics from Induction Times and Metastable Zone Widths. Cryst. Growth Des. 2013, 13, 2435–2440. [Google Scholar] [CrossRef]
  18. Alimi, F.; Elfil, H.; Gadri, A. Kinetics of the Precipitation of Calcium Sulfate Dihydrate in a Desalination Unit. Desalination 2003, 158, 9–16. [Google Scholar] [CrossRef]
  19. Lenka, M.; Sarkar, D. Determination of Metastable Zone Width, Induction Period and Primary Nucleation Kinetics for Cooling Crystallization of L-Asparaginenohydrate. J. Cryst. Growth 2014, 408, 85–90. [Google Scholar] [CrossRef]
  20. Van Driessche, A.E.S.; Stawski, T.M.; Kellermeier, M. Calcium Sulfate Precipitation Pathways in Natural and Engineered Environments. Chem. Geol. 2019, 530, 119274. [Google Scholar] [CrossRef]
  21. Söhnel, O.; Mullin, J.W. Interpretation of Crystallization Induction Periods. J. Colloid Interface Sci. 1988, 123, 43–50. [Google Scholar] [CrossRef]
  22. Söhnel, O.; Mullin, J.W. A Method for the Determination of Precipitation Induction Periods. J. Cryst. Growth 1978, 44, 377–382. [Google Scholar] [CrossRef]
  23. Alimi, F.; Gadri, A. Kinetics and Morphology of Formed Gypsum. Desalination 2004, 166, 427–434. [Google Scholar]
  24. Nemdili, L.; Koutchoukali, O.; Mameri, F.; Gouaou, I.; Koutchoukali, M.S.; Ulrich, J. Crystallization Study of Potassium Sulfate-Water System, Metastable Zone Width and Induction Time Measurements Using Ultrasonic, Turbidity and 3D-ORM Techniques. J. Cryst. Growth 2018, 500, 44–51. [Google Scholar] [CrossRef]
  25. Choi, J.Y.; Lee, T.; Cheng, Y.; Cohen, Y. Observed Crystallization Induction Time in Seeded Gypsum Crystallization. Ind. Eng. Chem. Res. 2019, 58, 23359–23365. [Google Scholar] [CrossRef]
  26. Selvaraju, K.; Valluvan, R.; Kumararaman, S. Experimental Determination of Metastable Zonewidth, Induction Period, Interfacial Energy and Growth of Non-Linear Optical l-Glutamic Acid Hydrochloride Single Crystals. Mater. Lett. 2006, 60, 1565–1569. [Google Scholar] [CrossRef]
  27. Ziegenheim, S.; Peintler, G.; Pálinkó, I.; Sipos, P. The Kinetics of the Precipitation of Gypsum, CaSO4·2H2O, over a Wide Range of Reactant Concentrations. React. Kinet. Mech. Catal. 2020, 131, 75–88. [Google Scholar] [CrossRef]
  28. Klepetsanis, P.G.; Dalas, E.; Koutsoukos, P.G. Role of Temperature in the Spontaneous Precipitation of Calcium Sulfate Dihydrate. Langmuir 1999, 15, 1534–1540. [Google Scholar] [CrossRef]
  29. Van Driessche, A.E.S.; Stawski, T.M.; Benning, L.G.; Kellermeier, M. Calcium Sulfate Precipitation Throughout Its Phase Diagram. In New Perspectives on Mineral Nucleation and Growth; Van Driessche, A.E.S., Kellermeier, M., Benning, L.G., Gebauer, D., Eds.; Springer International Publishing: Cham, Switzerland, 2017; pp. 227–256. [Google Scholar]
  30. Conley, R.F.; Bundy, W.M. Mechanism of Gypsyfication. Geochim. Cosmochim. Acta 1958, 15, 57–72. [Google Scholar] [CrossRef]
  31. Liu, S.-T.; Nancollas, G.H. Linear Crystallization and Induction-Period Studies of the Growth of Calcium Sulphate Dihydrate Crystals. Talanta 1973, 20, 211–216. [Google Scholar] [CrossRef]
  32. Fan, C.; Kan, A.T. Quantitative Evaluation of Calcium Sulfate Precipitation Kinetics in the Presence and Absence of Scale Inhibitors. SPE J. 2010, 15, 977–988. [Google Scholar] [CrossRef]
  33. Nicoleau, L.; Van Driessche, A.E.S.; Kellermeier, M. A Kinetic Analysis of the Role of Polymers in Mineral Nucleation. The Example of Gypsum. Cem. Concr. Res. 2019, 124, 105837. [Google Scholar] [CrossRef]
  34. Wang, Y.-W.; Kim, Y.-Y.; Christenson, H.K.; Meldrum, F.C. A New Precipitation Pathway for Calcium Sulfate Dihydrate (Gypsum) via Amorphous and Hemihydrate Intermediates. Chem. Commun. 2012, 48, 504–506. [Google Scholar] [CrossRef]
  35. Shiau, L.-D.; Lu, T.-S. A Model for Determination of the Interfacial Energy from the Induction Time or Metastable Zone Width Data Based on Turbidity Measurements. Cryst. Eng. Comm. 2014, 16, 9743–9752. [Google Scholar] [CrossRef]
  36. Bromley, L.A. Thermodynamic Properties of Strong Electrolytes in Aqueous Solutions. AIChE J. 1973, 19, 313–320. [Google Scholar] [CrossRef]
  37. Prisciandaro, M.; Olivieri, E.; Lancia, A.; Musmarra, D. Gypsum Scale Control by Nitrilotrimethylenephosphonic Acid. Ind. Eng. Chem. Res. 2009, 48, 10877–10883. [Google Scholar] [CrossRef]
  38. Prisciandaro, M.; Olivieri, E.; Lancia, A.; Musmarra, D. PBTC as an Antiscalant for Gypsum Precipitation: Interfacial Tension and Activation Energy Estimation. Ind. Eng. Chem. Res. 2012, 51, 12844–12851. [Google Scholar] [CrossRef]
  39. Sheikholeslami, R.; Ong, H.W.K. Kinetics and Thermodynamics of Calcium Carbonate and Calcium Sulfate at Salinities up to 1.5 M. Desalination 2003, 157, 217–234. [Google Scholar] [CrossRef]
Figure 1. Influence of concentration on induction time: • 0.15 M; ■ 0.20 M; ♦ 0.25 M; ▲ 0.30 M; × 0.35 M.
Figure 1. Influence of concentration on induction time: • 0.15 M; ■ 0.20 M; ♦ 0.25 M; ▲ 0.30 M; × 0.35 M.
Crystals 11 01236 g001
Figure 2. Induction time as a function of concentration.
Figure 2. Induction time as a function of concentration.
Crystals 11 01236 g002
Figure 3. Altered induction time as a function of concentration.
Figure 3. Altered induction time as a function of concentration.
Crystals 11 01236 g003
Figure 4. Influence of temperature on induction time: • 298 K; + 303 K; ■ 308 K; ♦ 313 K; ▲ 318 K; × 323 K.
Figure 4. Influence of temperature on induction time: • 298 K; + 303 K; ■ 308 K; ♦ 313 K; ▲ 318 K; × 323 K.
Crystals 11 01236 g004
Figure 5. Induction time as a function of temperature.
Figure 5. Induction time as a function of temperature.
Crystals 11 01236 g005
Figure 6. Influence of rotational speed on induction time: • 1.67 (s−1); ■ 3.33 (s−1); ▲ 5.00 (s−1); × 6.67 (s−1); ♦ 8.33 (s−1).
Figure 6. Influence of rotational speed on induction time: • 1.67 (s−1); ■ 3.33 (s−1); ▲ 5.00 (s−1); × 6.67 (s−1); ♦ 8.33 (s−1).
Crystals 11 01236 g006
Figure 7. Induction time versus rotational speed.
Figure 7. Induction time versus rotational speed.
Crystals 11 01236 g007
Figure 8. Microscopic photographs (magnification 50×) of obtained calcium sulfate dihydrate for T = 298 K, n = 1.67 s−1 and various concentration: (a) C = 0.15 M; (b) C = 0.20 M; (c) C = 0.25 M; (d) C = 0.30 M; (e) C = 035 M.
Figure 8. Microscopic photographs (magnification 50×) of obtained calcium sulfate dihydrate for T = 298 K, n = 1.67 s−1 and various concentration: (a) C = 0.15 M; (b) C = 0.20 M; (c) C = 0.25 M; (d) C = 0.30 M; (e) C = 035 M.
Crystals 11 01236 g008
Figure 9. Microscopic photographs (magnification 50×) of obtained calcium sulfate dihydrate for C = 0.25 M, n = 1.67 s−1 and various temperatures: (a) T = 298 K; (b) T = 303 K; (c) T = 308 K; (d) T = 313 K; (e) T = 318 K; (f) T = 323 K.
Figure 9. Microscopic photographs (magnification 50×) of obtained calcium sulfate dihydrate for C = 0.25 M, n = 1.67 s−1 and various temperatures: (a) T = 298 K; (b) T = 303 K; (c) T = 308 K; (d) T = 313 K; (e) T = 318 K; (f) T = 323 K.
Crystals 11 01236 g009
Figure 10. Microscopic photographs (magnification 50×) of obtained calcium sulfate dihydrate for C = 0.25 M, T = 298 K and various speeds: (a) n = 1.67 s−1; (b) n = 3.33 s−1; (c) n = 5.00 s−1; (d) n = 6.67 s−1; (e) n = 8.33 s−1.
Figure 10. Microscopic photographs (magnification 50×) of obtained calcium sulfate dihydrate for C = 0.25 M, T = 298 K and various speeds: (a) n = 1.67 s−1; (b) n = 3.33 s−1; (c) n = 5.00 s−1; (d) n = 6.67 s−1; (e) n = 8.33 s−1.
Crystals 11 01236 g010
Figure 11. Induction time versus supersaturation (a) with all investigated values (b) without the smallest value.
Figure 11. Induction time versus supersaturation (a) with all investigated values (b) without the smallest value.
Crystals 11 01236 g011
Table 1. Values used in Bromley equation [36].
Table 1. Values used in Bromley equation [36].
ConstantValueUnit
A γ0.509kg1/2/mol1/2
ρ 1kg1/2/mol1/2
C0.05-
B C a 2 + 0.0374-
δ C a 2 + 0.119-
B S O 4 2 0-
δ S O 4 2 −0.4-
Table 2. Thermodynamic solubility product.
Table 2. Thermodynamic solubility product.
T (K)log KspKsp (105)
298.15−4.372574.24
303.15−4.371684.25
308.15−4.375274.21
313.15−4.382964.14
318.15−4.393984.04
323.15−4.408203.91
Table 3. Comparison of turbidity of solution during precipitation for various concentrations.
Table 3. Comparison of turbidity of solution during precipitation for various concentrations.
Concentration (mol/dm3)t = 0 (s)t = 30 (s)t = 60 (s)
0.15 Crystals 11 01236 i001 Crystals 11 01236 i002 Crystals 11 01236 i003
0.20 Crystals 11 01236 i004 Crystals 11 01236 i005 Crystals 11 01236 i006
0.25 Crystals 11 01236 i007 Crystals 11 01236 i008 Crystals 11 01236 i009
0.30 Crystals 11 01236 i010 Crystals 11 01236 i011 Crystals 11 01236 i012
0.35 Crystals 11 01236 i013 Crystals 11 01236 i014 Crystals 11 01236 i015
Table 4. Comparison turbidity of solution during precipitation for various temperature.
Table 4. Comparison turbidity of solution during precipitation for various temperature.
Temperature (K)t = 0 (s)t = 30 (s)t = 60 (s)
298 Crystals 11 01236 i016 Crystals 11 01236 i017 Crystals 11 01236 i018
303 Crystals 11 01236 i019 Crystals 11 01236 i020 Crystals 11 01236 i021
308 Crystals 11 01236 i022 Crystals 11 01236 i023 Crystals 11 01236 i024
313 Crystals 11 01236 i025 Crystals 11 01236 i026 Crystals 11 01236 i027
318 Crystals 11 01236 i028 Crystals 11 01236 i029 Crystals 11 01236 i030
323 Crystals 11 01236 i031 Crystals 11 01236 i032 Crystals 11 01236 i033
Table 5. Comparison of turbidity of solution during precipitation for various rotational speed.
Table 5. Comparison of turbidity of solution during precipitation for various rotational speed.
Rotational Speed (s−1)t = 0 (s)t = 30 (s)t = 60 (s)
1.66 Crystals 11 01236 i034 Crystals 11 01236 i035 Crystals 11 01236 i036
3.33 Crystals 11 01236 i037 Crystals 11 01236 i038 Crystals 11 01236 i039
5.00 Crystals 11 01236 i040 Crystals 11 01236 i041 Crystals 11 01236 i042
6.67 Crystals 11 01236 i043 Crystals 11 01236 i044 Crystals 11 01236 i045
8.33 Crystals 11 01236 i046 Crystals 11 01236 i047 Crystals 11 01236 i048
Table 6. Supersaturation, molality, ionic strength, and activity of water for various concentration and constant temperature of 25 °C.
Table 6. Supersaturation, molality, ionic strength, and activity of water for various concentration and constant temperature of 25 °C.
Concentration (mol/dm3) m S O 4 2   / m C a 2 +
(mol/kgsolvent)
m N a +   / m C l
(mol/kgsolvent)
I a H 2 O γCa2+ γSO42−S
0.150.0750.1500.4500.9780.1834.264
0.200.1000.2000.6000.9720.1645.992
0.250.1250.2500.7500.9670.1517.906
0.300.1500.3000.9000.9620.14310.055
0.350.1750.3501.0500.9570.13812.505
Table 7. Supersaturation, molality, ionic strength and activity of water for various temperature and constant concentration of 0.25 (mol/dm3).
Table 7. Supersaturation, molality, ionic strength and activity of water for various temperature and constant concentration of 0.25 (mol/dm3).
Temperature (°C) m S O 4 2   / m C a 2 +
(mol/kgsolvent)
m N a +   / m C l
(mol/kgsolvent)
I a H 2 O   γCa2+ γSO42−S
250.1250.2500.7500.9670.1517.906
307.890
357.955
408.095
458.305
508.582
Table 8. Constants in equation on values B [5].
Table 8. Constants in equation on values B [5].
ConstantValuesUnit
N A 6.023·10231/mol
δ 0.012J/m2
R8.3145J/mol·K
V m 7.49·10−5m3/mol
β 16.755-
Table 9. Changes of value B as function of temperature.
Table 9. Changes of value B as function of temperature.
t (°C)T (K)B
25298.1515.494
30303.1514.740
35308.1514.034
40313.1513.372
45318.1512.752
50323.1512.169
Table 10. Comparison of the calculated induction times with measured ones.
Table 10. Comparison of the calculated induction times with measured ones.
Concentration (mol/dm3)SupersaturationAverage Measured Induction Time (s)Standard DeviationCalculated Induction Time (s)
0.154.2641604.81195
0.205.9929216.3481
0.257.906696.7456
0.3010.055491.3545
0.3512.505281.4739
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Małysiak, A.; Orda, S.; Drzazga, M. Influence of Supersaturation, Temperature and Rotational Speed on Induction Time of Calcium Sulfate Crystallization. Crystals 2021, 11, 1236. https://doi.org/10.3390/cryst11101236

AMA Style

Małysiak A, Orda S, Drzazga M. Influence of Supersaturation, Temperature and Rotational Speed on Induction Time of Calcium Sulfate Crystallization. Crystals. 2021; 11(10):1236. https://doi.org/10.3390/cryst11101236

Chicago/Turabian Style

Małysiak, Agata, Szymon Orda, and Michał Drzazga. 2021. "Influence of Supersaturation, Temperature and Rotational Speed on Induction Time of Calcium Sulfate Crystallization" Crystals 11, no. 10: 1236. https://doi.org/10.3390/cryst11101236

APA Style

Małysiak, A., Orda, S., & Drzazga, M. (2021). Influence of Supersaturation, Temperature and Rotational Speed on Induction Time of Calcium Sulfate Crystallization. Crystals, 11(10), 1236. https://doi.org/10.3390/cryst11101236

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop