Santanu Sarkara,
Ratul Chowdhurya,
Ranjana Dasa,
Sudip Chakrabortya,
Heechul Choib and
Chiranjib Bhattacharjee*a
aDepartment of Chemical Engineering, Jadavpur University, Kolkata, India. E-mail: c.bhatta@gmail.com; cbhattacharyya@chemical.jdvu.ac.in; Fax: +91 33 2414 6203; Tel: +91 98364 02118
bSchool of Environmental Science and Engineering, Gwangju Institute of Science and Technology, Gwangju, Korea
First published on 26th March 2014
The degradation of chlorhexidine digluconate (CHD) has been investigated using TiO2 suspension under UV irradiation with variation of pH, catalyst loading and substrate concentration within the range of 4.0–10.5, 0.1–0.4 g L−1 and 0.5–1.5 g L−1, respectively. The molecules of CHD have been adsorbed on TiO2 nano particles and degraded under UV irradiation. The influence of reaction parameters has been studied and reasons behind those have been correspondingly established. The molecular weights of by-products have been identified after photocatalytic degradation using mass spectroscopy. During the optimization of the process parameters, subtractive clustering fuzzy inference system has been employed to generate the optimum fuzzy rule base sets. A stable estimation technique has been attained by applying hybrid learning algorithm. Moreover, the performance and accuracy of the developed model has been tested using distinct test data set randomly selected from the experimental domain. The effects of different inherent parameters of subtractive clustering have also been studied and among them only squash factor is found to be the most important one for the present study. ANFIS has thus proved to be a congenial tool for the optimization of process parameters of the photocatalytic reaction discussed here.
Advanced oxidation processes (AOPs) have been adopted by several groups of researchers to mineralise different wastes before their discharge. Heterogeneous photocatalysis has now become a promising technique to eliminate organic and inorganic substances from waste streams in the presence of semiconductor oxides such as TiO2, MnO2 etc.10 Efficiency of this process has become well established in several research publications.11–14 Photocatalytic oxidation process depends on several process parameters like initial substrate concentration, catalyst loading, pH of the solution matrix, UV intensity etc. So, assessment of the process efficiency necessitates optimisation of the process conditions. For chemical reactions with multiple process variables, detection of optimum conditions is tedious with experimental data; however, availability of several optimization techniques has made it easier.11–15 The mentioned research groups have used response surface methodology (RSM) to optimize initial substrate concentration, catalyst loading and pH of the medium during photocatalytic degradation of pharmaceutical drugs and has found good agreement between experimental data with theoretical prediction.16 However, artificial neural network (ANN) is more advantageous than RSM considering the fact that ANN does not require any specified fitting function and it is very efficient to approximate all non-linear equations of different order. ANN approach has now been attempted by several research groups for optimisation of the chemical processes with multiple process variables with minimised error conditions.17 ANN is a computer programme based predictive model to simulate a system with the working principle of the human brain. It has the ability to make a system virtually stable and can be applied to establish the non linear relationship between process controlling parameters. The estimated outputs of the system using ANN were quite appreciable as the model was trained using experimental data. In 2008, Calza et al.18 adopted ANN and checked the predictability of the model. According to them the percentage error between the predictive and experimental values was very low. In general, the accuracy level of the ANN is far better than RSM as before prediction it is trained with experimental data.17 Basically, ANN can describe a non linear system through mapping between reasons and responses. Moreover, it is a simple simulation tool, and cannot be utilized to optimize the process parameters for a larger range of those parameters without the help of experimental data. Thus, such mechanistic approaches make the system more complicated. Very recently, to enhance the applicability and the performance of ANN, fuzzy logic was incorporated with it and it was called an artificial neuro-fuzzy inference system (ANFIS). The main advantage of the ANFIS is that it approximates the non linear system by setting IF-THEN rules19 and thus it is considered as a universal approximator.20 Moreover, in the research article by Sen et al.35 using ANN and KBHNN (knowledge based hybrid neural networks) on membrane filtration systems, they have reported that the final principle of working of the system was entirely dependent on the assumptions taken into account while formulating the problem. Again Sen et al.36 has reported that ANN needed rigorous theoretical modeling which has been avoided to make things simpler, but in the present study no assumptions have been incorporated as ANFIS has its learning algorithm suited better to analyze unknown dynamics.
The present study is focusing on the photocatalytic degradation of chlorhexidine digluconate (CHD) using TiO2 suspension and optimization of the process parameters using ANFIS. CHD, a class of antimicrobial drugs is used along with regular tooth brushing/flossing to treat gingivitis and to decrease the formation of mouth sores. It is also used in topical antiseptics and disinfectants as a preservative. In addition chlorhexidine is used as a pharmaceutical preservative.21,22 Although, it may cause damage or create irritation to the skin and eye on direct contact.23 Since it directly enters the environment after use, the gradual growth in its concentration is observed and thus it hampers aquatic life.24 Moreover, it has an adverse effect on aquatic and sewage microorganisms.24 Therefore, it is recommended that it should not enter directly to any type of water resources.23 To control the hostile effect of CHD as a pollutant, the present experimental study is undertaken to establish proper methodology for photocatalysis of CHD, and thereafter development of a suitable ANFIS structure for the process optimization.
Before the start of the reaction a continuous stirring was provided in the presence of TiO2 only in the reaction broth for thirty minutes to form a homogeneous suspension of nano particles. Then the reactor was placed under UV irradiation with the addition of CHD. At certain intervals samples had been collected and the nano particles were separated from the reaction mixture using cold centrifuge. The supernatant fluid was analyzed to measure the remaining CHD concentration after photocatalytic oxidation. The initial substrate concentration, catalyst loading and pH of the solution matrix had been chosen as process parameter, and the rate of degradation of CHD was observed with variation of one of them when the others were kept constant. Though intensity of UV source was varied, it was not chosen as the variable parameter to reduce the complexity in optimization process. Trial degradation procedure ensured that maximum removal of antibiotic was possible within one hour irradiation time and UV intensity of 80 μW cm−2. The effect of intensity of UV irradiation has already been investigated by Das et al.34
The removal of percentage of CHD was calculated using the following equation.
![]() | (1) |
:
40), operating at a flow rate of 1 mL min−1.
After the photocatalytic degradation, the molecular weights of the by-products were identified using mass spectroscopy of Quadrupole-TOF Micromass spectrometer (Waters Co., USA). The details of the mass spectroscopy and the HPLC chromatogram of depredated by-products have already been published by the same research group.34
The pH of the reaction medium during the experiment was measured with Sartorius AG Gottingen (PT-10P18602675) pH Meter.
A two rules based Sugeno type ANFIS can be written using the following equations.
| If x is A1 and y is B1 then f1 = p1x + q1y + r1 | (2) |
| If x is A2 and y is B2 then f2 = p2x + q2y + r2 | (3) |
The forward pass and backward pass are required to train the network. The forward pass propagates the input vector through the network layer by layer. In the backward pass, the error is sent back through the network in a similar manner to back propagation. In the layer 1 the output of the each node is calculated using eqn (4) and (5). O1,i(X) is the membership for x and y.
| O1,i = μAi(x), for i = 1, 2 | (4) |
| O1,i = μBi−2(y), for i = 3, 4 | (5) |
In layer 2 every nodes are fixed and t-Norm is used to ‘AND’ the membership grades. That can be explained with the help of eqn (6).
| O2,i = wi = μAi(x)μBi(y), i = 1, 2 | (6) |
Next, layer 3 contains fixed nodes which calculate the ratio of the firing strengths of the rules and can be represented as below.
![]() | (7) |
The nodes of layer 4 are adaptive in nature and it calculates the consequent of the parameters of ANFIS architecture i.e. On,i (where n = 4 and i = 1, 2,...), and can be written as eqn (8). The output of the system is calculated though single node of the layer 5 and that can be mathematically represented as eqn (9).
O4,i = ifi = i(pix + qiy + ri)
| (8) |
![]() | (9) |
Therefore the input vector is fed to the network through layer by layer and ANFIS learns the premise and consequent parameters for the membership functions and the rules. It includes both the advantages of ANN and fuzzy-logic model as well as eliminating some disadvantages during their lonely use. Moreover, the ambiguity of ANFIS is less than ANN. The training and convergence of ANFIS is faster compared to ANN based model.20 Overall, optimization of any process with minimum error is possible using the mentioned computational simulation method.
In the present study the initial CHD concentration, TiO2 loading and pH of the reaction medium had been selected as input parameters of ANFIS, and removal percentage of CHD was the corresponding output due to photocatalytic oxidation. For the present study, three types of values for initial concentration of CHDG, three types of TiO2 values and three types of pH values were chosen from the experimental data. Therefore, it would have A1, A2 and A3 for CHD initial concentration, B1, B2 and B3 for TiO2 and C1, C2 and C3 for pH values. Hence, the ANFIS architecture was schemed as described above.
Finally, it has been clarified that concentration of CHD is a model variable as well as being a part of the training data since, while training the system, the model had been fed with input parameter data and also the output data (concentration of CHD) so that when any unknown set of input data is entered, it would generate the CHD concentration as the simulated output. This would help to get the CHD concentration without performing the real experiment.
In the present context, the ANFIS Editor graphical user interface (GUI) is inbuilt in Fuzzy Logic Toolbox within the framework of MATLAB®V7.8.0 (R2009a) (The MathWorks, Inc., USA, R14) software, which was used for modelling and simulation purposes. The optimized sets of rules were generated using grid partition and subtractive clustering method, which contains eight different types of membership functions. Several numbers of rules may be produced by the system, and it is very much essential to optimize them either manually or automatically. It is assumed that each data point is a potential cluster centre in the subtractive clustering, which calculates the probability of each data point to be converted into a potential cluster centre, considering the availability of data points surrounding to it. The point with the highest potential is to be considered as the first cluster centre and calculate the potentials of other data points by subtracting the amount of potential as a function of its distance from the first cluster centre. Thereafter, the highest potential point is identified among the remaining point as the second cluster centre. Again the potentials of other data points with respect to the previous one should be revised. This process continues until the potential of the last centre should be less than the first one, thus almost all the data points come under consideration. The four principle parameters under consideration for subtractive clustering method are: Range of Influence, Accept Ratio, Reject Ratio and Squash Factor. Out of these, squash factor is the most significant as it helps in controlling the cluster size under consideration which indirectly affects the % error associated in ANFIS training of a given system.
![]() | ||
| Fig. 3 Effect of initial concentration of CHD on removal when TiO2 concentration and pH are constant. | ||
Overall, initial concentration of the substrate is a controlling parameter of the AOP and lower initial concentration should be maintained to achieve better degradation. Fig. 3 also demonstrates the effects of catalyst loading and pH on CHD conversion, which have been discussed in more detail in subsequent sections.
![]() | ||
| Fig. 4 Effect of TiO2 concentration on removal of CHD when initial concentration of CHD and pH are constant. | ||
Therefore, a proper ratio of catalyst to substrate should be maintained during photocatalytic degradation. Effects of initial substrate concentration and pH have also been demonstrated in Fig. 4; the explanation for the former has been given in a previous section and the effect of pH has been discussed in the following section.
After training the system with nearly 70% of the existing data set, we tried to simulate the remaining 30% data with minimum error (tolerance kept to 3% deviation). Our main challenge was to select the correct set of 70% data which would be able to mimic the system dynamics as good as possible so that the remaining 30% can be simulated without much error. This 70% data has to be selected using a heuristic method with a view to incorporate output values of both higher and lower magnitudes from the existing data set. Furthermore, 70% of the entire data set of 108 data points (all were experimental data) was used for training the system in silico. 70% of the entire data set was chosen as the training set as it had been seen that taking less data than that could not characterize the in silico system well enough hence the training error was quite high. At the same time, taking more data than the 70% data set as chosen would be all redundant as it would use up more computational memory and yet not improve the system training further. So 70% was in a way an optimal training set. The trained set is shown in Fig. 7.
Using only higher values to train the fis (fuzzy inference system) would necessarily incorporate a lot of error in predicting the lower range values and vice versa. Fig. 8 shows the prediction accuracy achieved in our case and that the % error (2.694) laid below the set tolerance limit of 3%.
Once the system has been tested to have replicated the “testing data” within tolerance limit, it can be used effectively as a dry-lab alternative of the experiment itself and hence can be used for further use with actually performing the experiment in wet-lab. However, it should be kept in mind that for successful prediction of an output corresponding to a given set of inputs, the input parameters should lie within the respective ranges as set within the ANFIS system.
We have shown a variation of ANFIS training error with increasing squash factor values ranging from 0.01 to 3 in Fig. 9. Our objective was to keep the % error below 3% (tolerance), so while increasing the squash factor values it was found that over the range 1 to 1.25 (encircled in Fig. 9) the error % was minimum (2.694%). So we chose the value of squash factor as 1.25 as a higher squash factor creates a bigger cluster which conforms to our objective of analysis.
This very applicability makes the ANFIS network all the more robust in its performance and serves as a suitable and more effective tool to predict the optimum conditions for a given multi-parameter reaction, than performing the experiment itself all over again.
The principal advantage of ANFIS over other approaches is that we just have to measure an influencing parameter and feed it to train the network to find its influence on the output; it need not have to provide its kinetics. Thus it saves a lot of computational load by avoiding complex differential equations which would otherwise have existed in the case of kinetic models. Hence ANFIS converges faster even for high degrees of non linearity.
Overall the present research article has been a portrait of the fate of heterogeneous photocatalysis of pharmaceutical waste water treatment and process optimization through ANFIS. The details of architecture of ANFIS have also been studied here.
| This journal is © The Royal Society of Chemistry 2014 |