Abstract

Colemanite is a preferred boron mineral in industry, such as boric acid production, fabrication of heat resistant glass, and cleaning agents. Dissolution of the mineral is one of the most important processes for these industries. In this study, dissolution of colemanite was examined in water saturated with carbon dioxide solutions. Also, prediction of dissolution rate was determined using artificial neural networks (ANNs) which are based on the multilayered perceptron. Reaction temperature, total pressure, stirring speed, solid/liquid ratio, particle size, and reaction time were selected as input parameters to predict the dissolution rate. Experimental dataset was used to train multilayer perceptron (MLP) networks to allow for prediction of dissolution kinetics. Developing ANNs has provided highly accurate predictions in comparison with an obtained mathematical model used through regression method. We conclude that ANNs may be a preferred alternative approach instead of conventional statistical methods for prediction of boron minerals.

1. Introduction

Natural resources play an immense role in creating the country’s wealth and spurring economic growth, and they have a significant impact on the human-social development of countries [1, 2]. In addition, they play an effective role in the developed countries’ technologies and level of prosperity, provide employment, satisfy energy needs and services, encourage subindustry and manufacturing, give regional development prominence, and provide foreign exchange [3]. Turkey has very diverse natural mineral deposits because of the loam and strategic position of its geographical features. Except for the petroleum and coal, there are 53 exploitable minerals and metals and 4,500 mineral deposits in Turkey [4]. Among these minerals, undoubtedly boron is the most important mineral in terms of mineral reserve and production capacity in Turkey, which has approximately 72% of the world’s boron reserves and by the end of 2007 has risen to the first place in the boron market [5].

Boron is one of the most important elements in the world due to its strategic and industrial value [6]. It is frequently used in nuclear engineering, high-quality steels, production of heat resistant polymers, cosmetic, leather, ceramics, rubber, paint, textile, agricultural catalysts, and in other industries [7, 8]. It is never found free in nature but invariably occurs as the B2O3, in combination with the oxides of other elements, forming borates of greater or lesser complexity [9]. It is known that there are more than 230 boron types in the world. Among these minerals are commercially important ones: tincal (Na2B4O7·10H2O), colemanite (Ca2B6O11·5H2O), and ulexite (NaCaB5O9·8H2O) [10, 11].

Boron mineral products, which contain impurities, are comprised mostly of clays and other constituent minerals that tend to be associated with clays [12]. Dissolution processes are applied to overcome these negative effects or to obtain product in industry.

Dissolution of the mineral is related to typical industrial processes, such as those in hydrometallurgy, medicine, oceanography, crystallography, ceramics, and desalination. On the other hand, it is also used for biological and environmental precipitation processes [13].

Colemanite, which is one of the most important boron minerals, is currently used on a large industrial scale. The upswing in the demands for minerals will continue in the years ahead as they are directly linked to rapid development of the technology. The increasing demand and new industrial use of boron compounds have increased their importance, and these compounds have been used as raw material in various areas of industry. There are numerous studies in the literature about dissolution of boron in order to meet the industry demands. Some of these studies are shown in Table 1. Statistical and numerical methods were usually used to determine dissolution kinetics in these studies.

Some obtained mathematical models are given as follows: Çopur et al. studied dissolution kinetics of colemanite in water saturated with carbon dioxide solutions. They have obtained a mathematical equation as follows by using statistical method [14]: where , , , , and indicate total pressure, solid/liquid ratio, particle size, temperature, and reaction time, respectively. Guliyev et al. studied colemanite in potassium hydrogen sulphate solutions. They used heterogeneous reaction models to determine the correlation between the dissolution rate and the parameters. Model equation was determined using both numerical and analytical methods as below [8]:

Kuşlu et al. investigated ulexite in borax pentahydrate solutions. In order to determine a mathematical model, they used integrated rate equations for the unreacted shrinking core model applying nonlinear regression analysis. The mathematical model was described as follows [10]:

In (2) and (3) , , , , and indicate concentration, stirring speed, mean particle size, solid/liquid ratio, and reaction time, respectively. was designated as dissolution rate for all equations. However, there are some disadvantages of using numerical and statistical methods:(i)a large data set is necessary in order to obtain reliable results,(ii)prediction data cannot accord with experimental data,(iii)complex calculating is needed,(iv)being time consuming.

Neural networks are one of the artificial intelligence techniques. It is based on present understanding of the biological nervous system and its ability to learn through example [15, 16]. Neural networks are used to solve the problems which cannot be modeled, in particular. A neural network can learn, adapt, predict, and classify. Prediction of parameters capacity of neural networks is very high. It provides more accurate results than the conventional statistical methods for prediction. Therefore, it has been used in different engineering applications [1619].

There is nearly no article in the scientific literature specifically devoted to a study of an integrated prediction of dissolution kinetics of boron minerals based on ANNs (see Table 1). In general, in order to determine prediction of dissolution kinetics boron mineral, conventional statistical methods were used by authors [14, 20, 21].

In this study, an artificial neural network was developed to predict dissolution kinetics of colemanite which was determined using the feedforward backpropagation neural network algorithm. ANNs performance was determined by using different learning methods, activation function, hidden layer, and neuron numbers to determine prediction of dissolution rate. ANNs were trained using experimental data [14]. Afterwards, prediction data was obtained from ANNs in comparison with previous mathematical model [14] which was formed using conventional prediction technique. ANNs model has given more accurate results than mathematical model.

This paper is organized as follows. In Section 2, the experimental study and used method are presented. Designing of an artificial neural network for prediction of dissolution rate and simulations method are considered in Section 3. Result of the simulations is addressed in Section 4. Finally, conclusions are presented in Section 5.

2. Material and Methods

2.1. Material

Dissolution processes were carried out by using colemanite mineral which was obtained from real boric acid plant (ETI Mine, Bandırma, Turkey). The samples were crushed by a jaw crusher and sieved using ASTM Standard sieves to obtain the following size fractions; 137.5, 213.5, 446, and 563.5 μm. The chemical characteristics of colemanite minerals, chosen parameters, and their ranges used in experiments are shown in Table 2.

2.2. Methods

An artificial neural network was developed to predict dissolution rate of colemanite in this study. Developed neural network results were compared with a regression based mathematical model which was obtained in a previous study. This mathematical model was formed using a series of experiments. Neural network was trained with data obtained from these experiments.

Experiments were carried out in a reactor (Parr 4848 reactor controller—Parr pressure reactor—Figure 1) which provides the temperature, pressure, stirring speed, and pH control. The schematic diagrams of experimental setup are shown in Figure 1. The solid prepared in accordance with the certain solid/liquid ratio was put into the reactor and 200 mL distilled water was added onto it. The system was set into the desired conditions and after the temperature of the reactor content reached to the determined value, CO2 gas was passed through the reactor until the air inside went to the outside. Later the gas outlet valve was closed and after the pressure value was adjusted to the desired level, the experiments were started.

At the end of the reaction, solutions were filtered by using a blue band filter and B2O3, Na+, and Ca2+ analyses were measured in permeate by using volumetric method and flame photometry, respectively. The mole fractions of B2O3, Na2O, and CaO passing into the solution from ulexite mineral were calculated by formulas given in (4), where designates quantity dissolution and shows the soluble compound:

3. Designing of an Artificial Neural Network for Prediction of Dissolution Rate

Artificial neural networks are computational systems that simulate the microstructure (neurons) of a biological nervous system. The most basic components of ANNs are modeled after the structure of the brain, and therefore even the terminology is borrowed from neuroscience [27]. ANNs consist of a large number of processing elements with their interconnections. They are essentially parallel computing systems similar to biological neural networks [28] called neurons, with each layer being fully connected to the proceeding layer by interconnection fully connected to the proceeding layer by interconnection strengths or weights [16].

ANNs are widely used to approximate complex systems that are difficult to model using conventional modeling techniques, such as mathematical modeling. There is no a certain method for selection of proper ANNs structure and training algorithm. The best solution is obtained by trial and error. On the other hand, the neural networks have a high prediction capability. Therefore a neural network, which has multilayer feedforward structure with two hidden layers, was designed (see Figure 2). Pressure, temperature, particle size of colemanite, solid liquid ratio, reaction time, and stirring speed are inputs of ANNs structure. The ANN was trained by using multilayer perceptron (MLP) networks. Backpropagation (BP) algorithm is the typical means of adjusting the weights and biases to obtain minimum the mean square error between the target and the network output by using a gradient descent algorithm [17, 29, 30]. Backpropagation gradient descent to minimize the target error, which is approximated in the vector space created by the weights and biases [31, 32].

In order to implement learning algorithm, each iteration of training includes the following procedures.

(1) Set the initial values of weights and .

(2) Compute the outputs for all neurons and layer, starting with the input layer as shown below: where is the weight between the input neurons and the hidden neurons, is the weight between the hidden and the output nodes, is the value of the input which consists of pressure, temperature, particle size of colemanite, solid liquid ratio, reaction time, and stirring speed, is the number of inputs of neuron in the hidden layer, is the value of the output for hidden nodes, is the number of neurons of the hidden layer, is the number of inputs of neuron in the output layer, is the output signals (dissolution of colemanite), and is the number of neurons of the output layer. In order to convert the input signals to the output signals, sigmoid transfer function can be used in ANN. Sigmoid transfer function formula is given below:

(3) Compute the error. In order to determine ANN performance, root mean square errors (RMSE), mean absolute errors (MAE), and coefficient of correlation () parameters were used. These are defined as where, is the mean of the measured data . RMSE measures residual errors that give a global idea of the difference between the observed and modeled values. And provides the variability measure of the data reproduced in the model. is the average of .

(4) The number of iterations is determined according to minimum value of the error function. When the error function reaches to sufficient value, iteration is stopped.

(5) Learning error computed for every neuron in all layers: where represents the total number of patterns, the desired outputs (experimental data) and the actual outputs.

(6) Update weights along negative gradient of error: where is the learning rate, is the momentum, and is the learning cycle.

(7) These procedures are repeated until the desired value of error.

All calculations were performed with software. The neural network structure is shown in Figure 2. In this figure, , , , and denote nodes input layer, hidden layer, and output layer, respectively. Weight of the nodes is referred to as . Subscripts specify the connections between the nodes. For example, is the weight between nodes and . The data were randomized and divided into two parts, training and testing. After the randomizing process, 65 data were used for training and 20 data were used for testing. Before applying the ANNs to the data, the training input and output values were normalized using the equation: where and symbolize the minimum and maximum of the data. Different values can be assigned for the scaling factors and . There are no fixed rules as to which standardization approach should be used in particular circumstances. This range increases the extrapolation ability of the ANN models. Therefore, these factors were assigned as 0.6 and 0.2, respectively [16].

In order to determine the best BP training algorithm, Bayesian regulation, BFGS quasi-Newton methods, and Levenberg-Marquardt (LM) algorithm with one and two hidden layers were trained and validated. Additionally, logsig, tansig, and radbas were used to find the best ANN transfer function. Different combinations of ANN structures with one or two hidden layers are tested in terms of iterations.

The following factors are employed during the training design: constant learning rate, constant moment, and learning cycles. Weights and bias values have been iteratively renewed using the various algorithms to minimize the RMSE between the network output and the target output. These parameters should be preferred to be as small as possible to obtain a best performance from ANN models, since the algorithm goes unstable [3134].

The ANN networks training was stopped after 625 learning cycles (epochs) since the variation of error was too small after this epoch. MSE is shown in Figure 3 according to epoch. The training was realized using the following parameters:(i)constant learning rate = 0.01,(ii)constant momentum factor = 0.9,(iii)mean-square error target = 10−6.

4. Results and Discussion

The performance comparison of BP algorithms in dissolution kinetics of colemanite mineral in water saturated with CO2 is given in Table 3. The different structures of ANN models and their outputs are given in Table 3. The performance of ANN outputs for each model was evaluated using three parameters: mean absolute error (MAE), the coefficient of correlation (), and root mean square errors (RMSE). As can be seen from Table 3, LM backpropagation and logsig with two hidden layers are more suitable as training and transfer function than others.

The test RMSE statistics of the ANN models are given in Table 4. The best result for minimum RMSE was selected. In this table, and represent the number of neurons in the first and second hidden layers, respectively. As can be seen from this table, the ANN model which has 7 neurons in the first hidden layer and 4 neurons in the second hidden layer has the lowest RMSE (0.0073), MAE (0.0061), and the highest (0.9975) value in test period. According to this result optimum ANN model was determined and structure of it is shown in Figure 4.

The performance of neural networks for training is shown in Figure 5. As can be seen from this figure ANNs results perfectly follow experimental results. The network was evaluated by comparing its predicted output values with experimental data which was shown in Figure 6. As it is seen in the figure the experimental data were found to be compatible with ANNs output. The maximum residual between the experimental data and ann output was determined to be around 0.0017 for training dataset. The residual characteristics have a decreasing trend. After training the neural network, test performance was checked. The performance of test was shown in Figure 7. Figure 7 also shows an analysis between the network outputs (estimations) and the corresponding targets (observed data) for the test data set. It is obvious that the predicted values from the trained neural network outputs catch the targets well. Residuals between testing and experimental data were shown in Figure 8. The maximum residual between estimations and observed data is determined about 0.002 for testing dataset.

Çopur et al. used heterogeneous reaction models in order to determine dissolution kinetics. They obtained a mathematical model by using numerical methods based on regression method. This model is given in (1). The dissolution rate was calculated using ANN test inputs with (1). Experimental, mathematical, and ANN results were compared. Results were shown in Figure 9. According to Figure 9 developed ANN results are closer to experimental data than the mathematical model results obtained from numerical methods. As is obvious from Figure 9 ANN gives better results than numerical methods.

5. Conclusion

The work presented here has demonstrated that ANN can be successfully employed to predict dissolution kinetics of colemanite mineral. In order to develop ANN models, total pressure, reaction temperature, particle size, solid/liquid ratio, and stirring speed parameters were used as input parameters and dissolution rate as the output. Experimental dataset was used to train multilayer perceptron (MLP) networks to allow for prediction of dissolution kinetics. So as to obtain most suitable prediction data, different learning methods, activation function, hidden layer, and neuron numbers were used. Levenberg-Marquardt backpropagation algorithm and Log-sigmoid (logsig) with two hidden layers were determined as training and transfer function. Also, ANN structure is comprised of 6 input neurons, 7 first hidden, 4 second hidden, and one output layers. This structure has the lowest RMSE (0.0073) and the highest (0.9975) values.

Developed ANN has given highly accurate predictions in comparison with an obtained mathematical model used through regression method. We conclude that ANN may be preferred as an alternative approach instead of conventional statistical methods for prediction of boron minerals. The prediction of dissolution kinetics can be obtained using with ANN quickly and accurately.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.