Determining the chlorine kinetic behavior in surface water using evolutionary metaheuristic algorithms

Farhad Mahmoudi Jalali; Seyyed Roohollah Masoomi; Mostafa Azizi; Reza Aghlmand; Mohammad Gheibi; Zahra Kian; Farhad Mahmoudi Jalali; Seyyed Roohollah Masoomi; Mostafa Azizi; Reza Aghlmand; Mohammad Gheibi; Zahra Kian

ISSN: 2692-4765

Annals of Systems Biology

Research Article Open Access Peer-Reviewed

Determining the chlorine kinetic behavior in surface water using evolutionary metaheuristic algorithms

Farhad Mahmoudi Jalali¹, Seyyed Roohollah Masoomi², Mostafa Azizi³, Reza Aghlmand⁴, Mohammad Gheibi^4* and Zahra Kian⁵

Author and article information

¹Department of Civil Engineering, Faculty of Engineering, Islamic Azad University, Tabriz Branch, Iran
²Department of Environmental Engineering, Civil & Environmental Engineering Faculty, Tarbiat Modares University, Tehran, Iran
³Department of Mine Engineering, Shahid Bahonar University of Kerman, Iran
⁴Department of Civil Engineering, Ferdowsi University of Mashhad, Iran
⁵Department of Chemical Engineering, Amirkabir University of Technology, Tehran, Iran

*Corresponding author: Mohammad Gheibi, Department of Civil Engineering, Ferdowsi University of Mashhad, Iran; E-mail: Mohammadgheibi@ymail.com

DOI: 10.17352/asb.000014

Received: 11 June, 2021| Accepted: 13 July, 2021 | Published: 14 July, 2021

Keywords: Chlorine mass degradation; Ortho-toluidine reagent; Multivariate nonlinear regression; Genetic algorithm

Cite this as

Jalali FM, Masoomi SR, Azizi M, Aghlmand R, Gheibi M, et al. (2021) Determining the chlorine kinetic behavior in surface water using evolutionary metaheuristic algorithms. Ann Syst Biol 4(1): 026-030. DOI: 10.17352/asb.000014

Abstract

The use of a reliable technique in measuring the residual chlorine concentration is of particular importance. Because based on the results of these measurements, the water quality of the network is calculated and evaluated in terms of ensuring the health of the community. The concentration of residual chlorine in the effluent of the treatment plant decreases as water flows through the distribution network. Chlorine decay is divided into two categories: wall and mass decay. Various detectors such as DPD and Ortho-toluidine can be used to measure residual chlorine. This study is the first step in the spectroscopy of chlorine and Ortho-toluidine color complex. The results of this stage of the research indicate that to measure the residual chlorine concentrations, the absorbed values should be read at a peak of 450 nm. In the next step, after measuring the free chlorine concentrations remaining at different times, the degree of chemical reaction was determined using multivariate nonlinear regression methods. In the final step, the reaction coefficient was calculated by simulating the differential equation in the Simulink environment and optimizing the error of the initial values of the reaction constant to the real values using the genetic algorithm (GA). The results of the analysis showed; The mass decay of chlorine in the sample follows the quadratic and the reaction constant is equal to 0.005 $\frac{L . \min}{m g}$ .

Main article text

Introduction

Disinfection is one of the main missions of water treatment systems. This process in water treatment plants usually involves primary and secondary disinfection. According to USEPA guidelines, the primary purpose of disinfection is to control biological growth in treatment ponds, prevent the formation of biological masses on filters, improve coagulation and flocculation, and reduce the taste and odor of raw water entering the treatment plant [1]. According to the definition of Standard 1053 of the Standards and Industrial Research Organization of Iran, drinking water is water that has physical, chemical, and biological properties to such an extent that its consumption for drinking in the short and long term does not cause human health [2]. In water treatment plants, depending on the type of source and its quality, different treatment units can be used. But the disinfection unit in all water treatment plants must be used to ensure microbial health. Currently, the use of chlorine for disinfection is the most common disinfection method in Iran due to its cheapness and antimicrobial power, and relatively good stability. In the disinfection process, five general mechanisms lead to the disruption of microorganisms and damage to their cellular structure. These mechanisms include damage to the cell wall, change in cell permeability, a change in the colloidal property of the protoplasm, a change in DNA or RNA, and an inhibition of enzymatic activity. Oxidizing chemicals such as chlorine inhibit the enzymatic activity of microorganisms [3]. Chlorine may come in many forms, compounds, and phases, such as; Sodium hypochlorite, calcium hypochlorite, chlorine dioxide, chloramine compounds, and chlorine gas should be used [4]. Ammar, et al. (2014) have investigated the degradation of chlorine dioxide in the desalinated water mass. In this study, for both primary and secondary disinfection, the effect of concentrations of 0.6 to 1.4 mg/l chlorine dioxide on mass decay rate was evaluated. Also, due to the effect of temperature on the degradation rates of chlorinated compounds, the experiments were performed at four temperatures of 20, 27, 35, and 45 degrees Celsius. The results of these experiments were computationally evaluated using T-test statistical methods [5]. Liu, et al. (2014) have investigated the effect of pH, water temperature, and the ambient temperature around pipelines in relation to chlorine mass degradation. This study was performed on case samples of two separate regions, which led to the presentation of models to express the relationship between chlorine degradation and the mentioned parameter [6]. Hua, et al. (2015) have conducted a study on chlorine decay and its reaction with organic matter (DOM) in water. In this study, models with quadratic reactions were presented to calculate the residual chlorine concentration [7]. Monteiro, et al. (2014) have investigated chlorine degradation in the urban water network using the EPANET MSX multiple simulation model. In this study, based on the reported results of other studies, chlorine mass degradation was considered as a quadratic type, and simulations were performed in the EPANET MSX environment. The results of this study showed; EPANET MSX environment in chlorine degradation simulation has similar results to the EPANET model, but due to the graphical features of the EPANET model than the EPANET MSX model, users are more interested in using the EPANET model [8]. In the present study, the intention is to first determine the optimal wavelength for measuring chlorine concentration, then using the synthetic experiments and computational nonlinear multivariate regression method, the degree of reaction is determined. Finally, using the combination of MATLAB software Simulink environment, and genetic algorithm, the values of the constant coefficient of the traditional reaction of chlorine mass decay are calculated.

Materials and methods

Chlorine mass decay

Chlorine reacts with organic and mineral molecules in the water matrix during movement along with the water supply network and its concentration decreases. This type of reduction in chlorine concentration is called chloride mass decay. Over the past few decades, scientists have conducted extensive research on the mass (volumetric) decay of chlorine. The results of the research have led to the development of several models for determining the mass degradation coefficient of chlorine (Table 1). Haas and Karra (1984) first proposed five practical models for chlorine mass degradation. In their research, all substances that can react with free chlorine were assumed to be a composite (A). The general form of the reaction and the differential equation used in the study are presented in equations 1 and 2, respectively [9].

Table 1: Models of chlorine mass degradation in various research results.
Proposed equations	Reference
$\begin{array}{l} D_{t} {= D}_{1} t^{n} \\ D_{1} \to Chlorine Concentration in first hour \end{array}$	[10]
$\frac{d C}{d t} = k_{1} C_{t} F_{1} + k_{2} C_{t} F_{2}$	[11]
$\frac{d C}{d t} = - k_{1} C_{1} - k_{2} C_{2}$	[9]
$\frac{d C}{d t} = - k C C_{R}$	[12]
$\frac{d C}{d t} = - k C$	[13-17]
$\begin{array}{l} C_{t} = C_{0} + k_{1} \ln (\frac{C_{0}}{C_{t}}) - k_{2} t \\ for treated water \to k_{1} = - 0.8147 C_{0}, k_{2} = 2.2808 (^{\frac{C_{0}}{[U V A]}) - 1.2971} [T O C] \\ for untreated water \to k_{1} = - 0.8404 C_{0}, k_{2} = 0.404 (^{\frac{C_{0}}{[U V A]}) - 0.9108} [T O C] \\ T O C \to Total Organic Carbon \\ UVA \to Ultra Violet A (315-400 nm) \end{array}$	[18]
$\frac{d C}{d t} = - k_{f} C C_{f} - k_{s} C C_{s}$	[19]

$Cl+A \to Product$ (1)

$\frac{d C}{d t} = - k^{'} C^{n} R^{m}$ (2)

As can be seen in the above relation, the parameters k’, C, R, m, and n are related to the mass synthesis coefficient is the chlorine concentration (mg/L), the concentration of chlorinated oxidants (mg/L), the degree of reaction of the oxidants, and the degree of reaction of chlorine. Since the concentration of oxidizing agents is much higher than the concentration of disinfectant (chlorine), the concentration of oxidizing agents is assumed to be constant over time. For this reason, equation 3 was proposed by Haas and Karra (1984) [9], as follows.

$\begin{array}{l} k = k^{'} R^{m} \\ \frac{d C}{d t} = - k C^{n} (3) \end{array}$

Finally, the researchers performed complementary modeling to provide a relationship for limit concentrations. It is assumed in some models; Over time, the chlorine concentration does not fall below known values (CL). The differential equation is partially described by equation 4.

$\frac{d C}{d t} = - k (C - C_{l}) C^{n - 1}$ (4)

As shown in Table 1; Many researchers divide the chlorine decay process into two parts: fast reaction (k1) and slow reaction (k2). Al-Jasser, et al. (2007) [17] argue for two phases (fast and slow) in considering the chlorine decay process; The fast phase of the reactions takes place in the treatment plant or inside the tanks, in places where a load of oxidizing substances (impurities in water) is high, and also the slow phase takes place along the lines of the water supply network. It is worth noting; Where re-chlorination occurs; Both fast and slow reactions occur.

In this study, to determine the rate of chlorine mass degradation, water outlet samples of the filtration unit of a treatment plant located in Khorasan Razavi province were used. The treatment plant is fed by surface water sources. In these samples, the chlorine degradation rate was evaluated over a period of 4 hours. An Agilent 8453 spectrophotometer and spectrocell were used to measure the chlorine concentration. According to Beer-Lambert law, there is a direct and linear relationship between the amount of light absorbed by the samples and the concentration of the substance in them [20]. The detector used in the experiments was Ortho-toluidine, which produces a distinct yellow color in the presence of chlorine. For this purpose, at the beginning of the measurements, it was necessary to perform a series of spectroscopic studies to determine the maximum absorption wavelength by the yellow color of the chlorine and Ortho-toluidine complex. Therefore, standard samples with different concentrations of chlorine were prepared and spectroscopically examined. It should be noted that to prepare standard solutions, bleach-containing 8% NaOCl was used. To determine the percentage of NaOCl content of purchased bleach, the iodometric method was used [21]. In the titration steps of the iodometric method, three chemicals of potassium dichromate, sodium thiosulfate, and starch adhesive were used. In order to increase the accuracy of the test and to prevent possible errors, the three mentioned substances should be titrated separately and the purity of each of them should be calculated and evaluated. In the next step, the values of the standard samples (containing concentrations of 0.25, 0.5, and 1 mg / l of bleach water) were prepared using the above-mentioned bleach water solution according to the iodometric test outputs. A portion of Ortho-toluidine was added and each sample was spectroscopically measured in four replications in the visible wavelength range (380-600 nm) and their average absorption was calculated. According to the results of spectroscopy, the wavelength of 450 nm (Lambda = 450 nm) was obtained as the maximum absorption index. The spectroscopic results of standard solutions are presented in Figure 1, which well shows the maximum wavelength. After determining the maximum absorption wavelength, the necessary measurements were made to determine the residual chlorine concentrations in the desired water samples. Finally, after performing synthetic experiments on chlorine mass degradation, the values of reaction degree and reaction constant were calculated.

Calculation of reaction degree and constant coefficient of chlorine mass degradation using statistical techniques and differential equations.

In this part of the study, a combination of Simulink environment and genetic algorithm (GA) was used to calculate the synthetic coefficients of chlorine mass degradation. The method of solving this equation (Haas and Karra model [9]) is based on the ODE 45 method of MATLAB software in the analysis of differential equations (Figure 2). As can be seen in equation 5, the equation contains unknowns including n (degree of reaction) and k (synthetic coefficient of constant reaction). To calculate the degree of chemical reaction from regression analysis, changes in chlorine concentration over time were performed as shown in Table 2. The error optimization equation (Error function in Figure 2) and the regulatory parameters of the genetic algorithm are also expressed in Equation 6 [22].

$\frac{d c}{d t} = K C^{n}$ (5)

$\begin{array}{l} {Error Function = min (K}_{I n i t i a l} - K_{Re a l}) \\ M u t a t i o n Rate = 0 .01 \\ Cross Over Probability = 0 .8 \\ Initial Population = 50 \\ Number of Iteration = 500 \\ K_{I n i t i a l} \to Random Value of K \\ K_{Re a l} \to K Value Calculate by ODE45 method (6) \end{array}$

Results and discussion

The synthetic experiment on chlorine mass degradation was performed at different time steps. The results are summarized in Table 3. The table contains the results of calculations for determining the degree of reactions zero, one, and two. Also, the data in Table 4 for zero, one, and two-degree reactions are shown in Figures 3, 4, and 5, respectively.

The application of nonlinear multivariate regression analysis to determine the degree of reaction of chlorine mass degradation led to the results summarized in Table 3.

From the analysis of the results of the above table, it is clear; The quadratic reaction creates the best compatibility due to the maximum amount of R-Square and the lowest amount of error (SE). In other words, the second degree can be chosen as the most appropriate value of the degree of reaction. A review of the research background reveals; Due to the small difference between the R-Square and SE statistical indices in the first and second degree, the reaction can be declared as the first degree. Given that the qualitative-hydraulic modeling processes are performed on first-grade enamel and there is little difference, the existing claim can be stated as the result of measurements. Therefore, the engineering judgment of the researcher determines the degree used in qualitative-hydraulic modeling. In order to calculate the constant mass degradation coefficient k, a combination of the mentioned model in Simulink environment and GA algorithm was used. The results of the integrated model showed that the value of the constant mass decay coefficient of chlorine is equal to 0.005 $\frac{L . \min}{m g}$ , which in the linearization mode gives a value equal to 0.0045 $\frac{L . \min}{m g}$ .

However, determining the value of the constant reaction coefficient by the combined method is more accurate and takes precedence over the linearization method (straight-line slope) which resulted in an error of 10%.

Conclusion

Measuring and calculating chlorine deterioration in the treatment plant and water supply network is of particular importance. There are various solutions and methods for the synthetic study of chlorine mass decay. One of these methods is the use of Ortho-toluidine reagent, which showed the results of spectrophotometry using spectrophotometry (UV-Visible); The readings should be 450 nm in length. After measuring chlorine degradation, the degree of chemical reaction, as well as the reaction constant, should be determined. This study has modeled the general shape of the equation in the Simulink environment of MATLAB software to determine the degree of reaction from multivariate nonlinear regression analysis and to determine the reaction constant. Then, with the help of the GA algorithm, it optimizes the error of the initial values of the reaction until it reaches the real value. The results showed that the case reaction of this study followed the second degree and its reaction constant was equal to 0.005 $\frac{L . \min}{m g}$ .

References

Copyright

© 2021 Jalali FM, et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

View Similar Articles

Order for reprints

Article Alerts

Subscribe to our articles alerts and stay tuned.

This work is licensed under a Creative Commons Attribution 4.0 International License.

Quick Enquiry

Table 2: Nonlinear regression analysis of the kinetic reduction of chlorine concentration in fluid mass to determine the degree of the chemical reaction.
Describe the steps	Steps of statistical analysis of data
In this section, regression analyzes and the R2 coefficient for graph (c – t) are examined.	Zero degrees control
In this section, regression analyzes and coefficient R2 for the graph (Ln (c) - t) are examined.	First-degree control
In this section, regression analyzes and the R2 coefficient for graph (1/c - t) are examined.	Second-degree control
It is determined based on the highest value of R2 and the lowest degree of reaction.	Determine the degree of reaction

Table 3: Multivariate nonlinear regression analysis to determine the degree of chloride mass degradation reaction.
SE	R-Square	Degree of reaction
0.09	0.92	0
0.058	0.972	1
0.04	0.986	2

Table 4: Average results of chlorine mass degradation experiments and calculations related to the degree of reaction.
Time(min)	[C] - ppm	Ln[C]	[C]-1
0	1.43	0.35	0.69
10	1.40	0.34	0.71
20	1.20	0.18	0.82
30	1.16	0.15	0.85
50	1.13	0.12	0.88
70	1.00	0.008	0.99
100	0.87	-0.14	1.14
130	0.79	-0.22	1.25
160	0.67	-0.39	1.48
220	0.62	-0.47	1.60
254	0.53	-0.63	1.87

Annals of Systems Biology

Determining the chlorine kinetic behavior in surface water using evolutionary metaheuristic algorithms