Optimization of Calcium Sulfate Precipitated in the Laminar Flow Pipe through Response Surface Modeling of Temperature, Ca2+ Concentration and Citric acid Additives

This paper presents results of investigation of scaling of calcium sulfate (CaSO4) on metallic pipes. In this study, the optimizing variables; namely temperature (50-600C), concentration of Ca2+ (2000-3000 ppm), citric acid concentration (10-20 ppm) were set-up to provide the optimum yield of the mass scales. The mineral scale detected using XRD is mainly gypsum, and the scale has a plate like morphology under SEM examination. SRM (surface response methodology) prediction provided that the temperature is a significant factor, while the Ca2+ concentration and the citrate are insignificant variables determining the optimal condition of the mass scale yields. An optimized mass scale response of 119.99 mg was obtained at a temperature of 56.36 OC, concentration Ca2+ of 2649.21 ppm and citric acid concentration of 12.11 ppm, respectively. The addition of citric acid did not modify the crystal morphology, but may control the size of crystals.


INTRODUCTION
Calcium sulfate (i.e., gypsum, anhydrite and hemi-hydrate), is commonly formed as scale minerals in pipes of industries ranging from oil and gas to desalination.Scaling of calcium sulfate occurs when the surfaces of metallic pipes contact with the supersaturated water, which can create in technological problems such as the diminution of heat transfer, increase of energy consumption and unscheduled equipment shutdown [1][2] .Additionally, this existing scale cannot be removed by acid, but only it can be taken away by mechanical method.Nowadays, much effort has been made in meditating on the precipitation of calcium sulfate in view of their applications in a routine of industrial and environmental engineering 2 .Written report on this mineral scale formation have been previously focused on examining precipitation formed in the bulk solution by using laboratory beaker or bulk jar tests 3 .
Here two mechanisms of crystallization of calcium sulfate are commonly experienced in the water system.Firstly, homogenous crystallization frequently occurs in the aqueous solution in the absence of any foreign ions.Secondly, the calcium sulfate particles are mainly formed in heterogeneous processes by the presence of a foreign solid phase [4][5] .
Further study on the aspect of calcium sulfate scaling has resulted in numerous studies for the methods to detect and assess the scale formation on metal surfaces.Calcium sulfate has been recognized as a mineral with relative insolubility, thus it can be easy to be precipitated where calcium and sulfate present in the aqueous solution and subsequently lead to the scale formation even at low pH.Theoretical solubility model of the calcium sulfate compounds in water and in multi-component aqueous solutions, has been discussed previously 2 .The solubility of calcium sulfate hydrates may be expressed as the following equation: where n = 0, 0.5 and 2 corresponding to anhydrite, bassanite and gypsum, respectively.Here, gypsum is considered as the most stable phase found in lower temperatures, while anhydrite can be commonly formed at higher temperatures 6 .Additionally the scale formation of calcium sulfate is influenced by many factors such as concentration of solution, pH value, temperature, pressures, and ionic strength 7 .Correspondingly, sulfate scaling of Ca-ions is of particular interest because this precipitation has shown inverse solubility, i.e. its solubility reduces with increasing temperature.Because calcium sulfate polymorphism is a complicating factor, the detailed mineralogical study of deposits is required to understand the mechanism tendency scaling of the surface [6][7][8] .
Mode of scale deposit for mation is reported to be classified as two groups 2 : (i) surface crystallization in which mineral precipitated more or less selectively onto the surfaces of the equipment in contact with the aqueous fluids (usually at elevated temperature) and (ii) bulk crystallization which involves the accumulation of precipitates as a result of sedimentation or transport by fluid flow.In this second mode, the scale can be created spontaneously in the bulk, once the supersaturated solution increases, or corrosion by-products are formed at a second stage sediment out.Moreover, scaling in pipe system may be imputed by a combination of these two pathways and also kept in line by the process condition in pipes.Correspondingly, the mitigation of mineral scaling may be performed by the utilization of chemicals and antiscalants for inhibiting nucleation, crystal growth or both.For instance, polyacrylic acid, polyacrylamide, polymaleic anhydride, and polyphosphates are commonly employed [9][10] .At that place are likewise a large number of commercial antiscalants with new formulations being continuously upgraded and tried out for a form of scaling species 7 .Among the various systems of using chemical additives, the prevention of scaling depends on their chemistry and the nature of the solids forming 7 .Here, the use of poly (citric acid) has been demonstrated to suppress the formation of calcium sulfate scale at low concentrations 11 .The poly (citric acid) has the absorption ability of the surface of calcium sulfate scale crystals, thus distorting their scale crystal polymorphs.Therefore, it is potential as a calcium sulfate scale inhibitor, which can be prepared from monomer citric acid.Additionally, the citric acid is usually utilized as additive in the transformation of bassanite to gypsum 12 .
Instead, citrate has widely known as scale inhibitor for calcium carbonate (CaCO 3 ) 13-14-15-16 crystallization, because it can absorb on the calcium carbonate crystal for prevention of crystal particles growing.Nevertheless, a problem of using citric acid in the water system may be linked to the long retention time of crystal growth.Therefore, in that location is still considerable debate and doubt over the mechanism of the effects of variable concentrations of a carboxylic acid and/or various carboxylic acid used in the prevention of crystal growth [17][18] .To come up to this gap results may need a detailed research on the effects of the carboxylic acids (citric, maleic and tartaric) and variable additive concentrations on the kinetics and phase morphologies that develop during calcium sulfate formation reactions and it can subsequently gain a more mechanistic understanding of the crystallization 19 .
The present study was undertaken for examining the calcium sulfate scale precipitation in pipes through controlling the independent operating variables.The variables investigated were: temperature, Ca 2+ concentration and citric acid additives.Material characterization, including SEM for morphological analysis and XRPD for phase composition was applied in the study.In this work, the central composite design (CCD) for calculating the surface response methodology (SRM), involved temperature and Ca 2+ concentration, the citric acid additives to maximize calcium sulfate precipitation.A second order transferred polynomial model (as an inverse model) was determined for the yield of mass scale as a function of these variables.

Reagents
Materials needed in the preparation of the supersaturated solutions were calcium chloride (CaCl 2 ) and sodium sulfate (Na 2 SO 4 ) powder with the analytical grade (Merck™).The citric acid (C 6 H 8 O 7 ) with analytical grade (Merck™) was also used as additives.The crystal-forming solutions were made using distilled water.

Procedures and measurements
Experiments were carried out using a laboratory equipment for calcium sulfate crystallization is shown in Fig. 1.Crystallization experiments were performed similar to the method proposed by Muryanto et al (2014) 20 , by dissolving each of CaCl 2 and Na 2 SO 4 crystals in the distilled water (500 ml) for providing Ca 2+ concentration of 2000, 2500 and 3000 ppm.Citric acid (10, 15 and 20 ppm) was added by dissolving the crystals C 6 H 8 O 7 in a vessel containing 500 ml of solution CaCl 2 , then stirred until well mixed.For the experiment, one liter containing each of equimolar CaCl 2 and Na 2 SO 4 was prepared.The solutions at predetermined concentrations were separately placed in the two vessels and equilibrated until the designated temperature (50-60 O C) was reached.The conductivity of solutions, leaving the test pipe section was continuously checked for up to 2 hours.The scale deposited on the surface of the coupons within the housing sample was carefully removed and dried in an oven overnight at 60 O C. The dried coupons were taken out to be weighed.The weight difference of the coupons obtained before and after the experiments was identified as the scale mass.The scale deposited on the surface of the coupons was carefully removed and kept in a plastic container for subsequent characterization.

Materials characterization
All powder samples were carbon-coated prior to investigation by scanning electron microscopy (SEM) (FEI Inspect S50) with energy dispersive spectroscopy (EDS) system fitted with a field emission source and operating at an accelerating voltage of 15 KV.Phase identification of scale mineral was conducted by XRPD (X-ray powder diffraction) analysis.XRPD data were acquired by a conventional Bragg-Brentano (BB) diffractometry with parafocusing geometry and Cu-Ká monochromated radiation.The scan parameters (5-90 2è O , 0.020 steps, 15 s/step) were taken as required for observation.A PC-based search match program, the Philips X'Pert Software (Philips Electronics N.V) was to identify possible crystalline phases in the samples.In this method, the peak positions and peak heights were verified against the entries in the ICDD-PDF (International Centre for Diffraction Data-Powder Diffraction File).The identified mineral phases were subsequently adjusted by Rietveld method using Fullprof-2k, software, program version 3.30 [21][22] .The crystal structure model for the Rietveld refinement was obtained in the AMCSD (American Mineralogist crystal structure database).The detailed discussion of the method is provided elsewhere [23][24] .

Experimental design
In this study, the optimization of the three variables (temperature, Ca 2+ concentration, and citric acid) to yield the optimum mass scale was performed using SRM within the CCD (Table 1).SRM calculation was conducted by the statistical v. 6 software packages (StatSoft, Tulsa, OK, USA).Using this method, the proper response value and mathematical model fitted to the measured data was acquired from the experiments, and the independent variables of optimal conditions.

Properties of solid precipitated crystals
The corresponding solid crystals were subjected to XRPD as shown in Figures 2a and 2b.The quality of Rietveld refinements of the crystalline scale may be gauged from the diffraction plot (Fig. 2a), where the assigned phases clearly stand out in the difference curve of the calculated and the measured diffraction profile.The Rietveld analysis of crystals proved that only the gypsum polymorph present in scale obtained from the solution in the absence of additives at a temperature of 50 O C. In this experiment, metastable crystalline phases of bassanite and amorphous calcium sulfate were not identified.Depending on the amount of citric acid additives, gypsum crystals were predominantly precipitated on the surface of all coupon samples (Fig. 2b).Nevertheless, in the presence of 10 ppm citric acid additive, it is shown that the probability amorphous phase of calcium sulfate, which have peaks with low intensity, could be identified by their XRPD patterns.
Further, gypsum crystal habit, morphology and size did not vary throughout the sample depth depending on the additive concentration.The surface of gypsum crystals was composed of micrometersized monoclinic crystals with a plate-like morphology (Fig. 3).This gypsum crystals are same as the finding in the literature 18 .These crystals formed a layer at the surface within a size of around 10 µm (Fig. 3a-c).In the presence of 20 ppm citric acid, the shape, size or habit of the end-product gypsum were similar to the CaSO4 experiment without additives, although several small size of the gypsum needles was observed.The citric acid additive may control the growth of the scale mineral, however, there was no
Based on multiple regression analysis of the experimental data, the optimization resulted in the following second-order polynomial equation in term of code values: Y = -452.9011+ 15.24867 X 1 -0.03793X 1 2 + 0.0868 X 2 + 0.000002 X 2 2 -21.21898X 3 + 0.3642 X 3 2 -0.00143X 1 X 2 + 0.495 X 1 X 3 -0.00116X 2 X 3 where Y is the yield of mass scale (mg), and X 1 , X 2 and X 3 are the coded variables for temperature, Ca 2+ concentration and citric acid concentration, respectively.Analysis of variance (ANOVA) for the statistical testing of the model is shown in Table 3.
Here, the ANOVA fitted the quadratic polynomial model of mass scale yield.The influence of the significance of a factor can be seen from F-value and p-value.The quadratic regression model showed the value of determination coefficient (R 2 ) of 0.965 with no significant lack of fit at p > 0.05 which means that the calculated model fitted 96.5 % of the result and only 3.5% of the total variation did not fit to the model.The p-value less than 0.05 with an accuracy of 95% indicates that the relationship between the response and the independent variables can be fitted using the model.The significance of the model was also judged by F-test, where F -value is defined as the ratio between MSF (mean squares of factor) of  the MSE (mean squares of error).A factor can be said to have a significant effect when the F -value is greater than F-table.
The F -value of model compared with the F-table shown in Table 3, provided that the F -value is greater than the F-table.This demonstrates a substantial outcome on the response variable of the mass scales.The same matter can be seen in chart Pareto, as depicted in Figure 4.
It shows that p-value is less than 0.05 providing that the independent variables have insignificant effect.Apparently, the optimization variables for the response of the scale mass (mg) may be related to the temperature (X 1 ).In contrast, the influence of linear concentration (X 2 ), linear of citric acid (X 2 ), the interaction between temperature and citric acid (X 1 X 3 ), quadratic concentration (X 2 2 ), quadratic of citric acid (X 3 2 ), the interaction between temperature and concentration ( X 1 X 2 ), the interaction and the concentration of citric acid (X 1 X 2 ), quadratic temperature (X 1 2) can be ignored because they may not provide a significant effect on the response of the mass scales.

Optimization of independent variables for the response of mass scale
The interaction of independent variables in response of mass scale yield was depicted by 3D  response surface and 2D contour plots created by the model (Figs.5-7).Different shapes of the contour plots showed different interactions between the variables, where the significant interactions between the variables was represented as an elliptical contour plot.In contrast, a circular contour plot presents the insignificant relationships between independent and dependent variables.Fig. 5 presents the interaction between temperature (X 1 ) and Ca 2+ concentration (X 2 ) on the yield of mass scale.When the temperature below 50 O C, there was a small amount in the response, and the temperature below 50 O C did not show any obvious effect on extraction yield.However, the increase of temperature from 50 to 66 O C make improve the mass scale yield.The possible explanation could be that the increase in temperature can increase reaction speed and collisions between molecules occurring more quickly, so that the mass scale formation increases 25 .
Fig. 6 describes the effect of Ca 2+ concentration (X 2 ) and additive concentration (X 3 ) on the yield of mass scale.It can be observed that when Ca 2+ concentration (X 2 ) was fixed at 0 levels, additive concentration (X 3 ) displayed a quadratic effect on the response yield.Varying Ca 2+ concentration from 1500 to 3000 ppm with an increase of citric acid additives from 8 to 16 ppm, the target compounds' mass scale yield was increasing with increase of Ca 2+ concentration.Fig. 7 illustrates that an increase in temperature (X 1 ) leads to increased production of mass scale.Maximal mass scale production was obtained at a citric acid concentration (X 3 ) of 12.11 ppm.Here, the optimal range for mass scale production was from 50 to 66 °C and citrate of 10 to 20 ppm.
To determine the optimum mass scales was done by entering values into the equation, which include the optimum ratio of the mass response for variable optimization.Optimum mass scales are shown in Table 4.

Validation of the predicted value for the optimal variable
The results of SRM optimization were then verified using the experimental data.Verification was conducted by comparing the results of SRM optimization with the results of laboratory experiments.The calculation generated to the optimum yields of mass scale was based on temperatures (56.36 O C), Ca 2+ concentration (2649.21ppm) and citric acid additives (12.11 ppm), respectively.As shown in Table 4, three parallel experiments were carried out under the optimal conditions.The software predicted that the precipitated yield of total mass scale was 119.99 mg.The comparison was performed to acquire the % error of difference in optimizing response.Table 5 present results of experiments conducted in conformity with the optimum ratio of SRM for analysis yielding the mass scale of 121.34 mg.From the results of these experiments, the % error was calculated.Here the calculated % error for the mass scale response was 1.11 %.This means for the mass scale response has an accuracy of 99.99 %.So the results of the SRM optimization analysis were in agreement with the data of experimental mass scales.

CONCLUSIONS
SRM optimization of variable operating conditions, provided that the most influence on the mass of the scale is temperature.The optimum result with the mass response scale (119.99 mg) was related to temperature of 56.36 O C; concentration Ca 2+ of 2649.21 ppm; citric acid of 12.11ppm, respectively.The crystalline phase of the scale was found to be mostly gypsum as shown by the XRD, although the amorphous phase may be developed on the scale in the solution with 10 ppm additive and a temperature of 50 O C. The SEM analysis also showed that the scale has a plate like morphology.The addition of citric acid seemed to not change the crystal morphology, only the possibility of citric acid additives may control crystal growth on the surface of the particles.

Fig. 1 :Fig. 2 :
Fig. 1: Schematics of equipment used for calcium sulfate scale formation in pipes

Fig. 3 :
Fig. 3: Morphology of the calcium sulfate precipitated in the solution of 2000 ppm Ca 2+ concentration with; a) 0 ppm; b) 10 ppm and c) 20 ppm citric acid additives and temperature of 50 O C

Fig. 4 :
Fig. 4: Pareto chart of optimization independent variables on the response of the mass scales (mg)