Next Article in Journal
Questions and Concerns Re: Blue Water Footprints Reported in “Water Footprint of Meat Analogs: Selected Indicators According to Life Cycle Assessment”
Previous Article in Journal
Hydrologic Validation of MERGE Precipitation Products over Anthropogenic Watersheds
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Micellar-Enhanced Ultrafiltration to Remove Nickel Ions: A Response Surface Method and Artificial Neural Network Optimization

Faculty of Engineering and Applied Science, Memorial University, St. John’s, NL A1B 3X5, Canada
*
Author to whom correspondence should be addressed.
Water 2020, 12(5), 1269; https://doi.org/10.3390/w12051269
Submission received: 29 March 2020 / Revised: 20 April 2020 / Accepted: 27 April 2020 / Published: 30 April 2020
(This article belongs to the Section Wastewater Treatment and Reuse)

Abstract

:
Nickel ions from aqueous solutions were removed by micellar-enhanced ultrafiltration (MEUF), using the surfactant sodium dodecyl sulfate (SDS) as a chelating agent. Process variables and indicators were modeled and optimized by a response surface methodology (RSM), using the Box–Behnken design (BBD). The generated quadratic models described the relationship between a performance indicator (nickel rejection rate or permeate flux) and process variables (pressure, nickel concentration, SDS concentration, and molecular weight cut-off (MWCO)). The analysis of variance (ANOVA) showed that both models are statistically significant. To remove 1 mM of nickel ions, the optimal condition for maximum nickel removal and flux were: pressure = 30 psi, CSDS = 10.05 mM, and MWCO = 10 kDa, resulting in a rejection rate of 98.16% and a flux of 119.20 L/h∙m2. Experimental verification indicates that the RSM model could adequately describe the performance indicators within the examined ranges of the process variables. An artificial neural network (ANN) modelling followed to predict the MEUF performance and validate the RSM results. The obtained ANN models showed good fitness to the experimental data.

1. Introduction

Nickel is a common heavy metal generated from various industrial activities such as mining, electroplating, batteries manufacturing, metal finishing, and forging. It is carcinogenic, non-biodegradable, and could accumulate and persist in the nature and living organisms [1,2]. Even at low concentrations, nickel can be toxic to the environment and humans. Conventional methods (e.g., chemical precipitation, adsorption, ion exchange, electrodialysis) when treating large-volume aqueous solution containing low concentration of heavy metals (e.g., nickel) can be challenged by secondary pollution of deposition, high cost, low selectivity, and difficulties of recycling metals [3,4]. Although membrane technologies such as reverse osmosis have been successfully used to remove metal ions from aqueous solution with high removal efficiency, their high operational and maintenance costs hinder their wider application. To overcome these drawbacks, micellar-enhanced ultrafiltration (MEUF) provides an alternative for heavy metal removal as it can achieve high removal rate and high permeate flux under mild conditions with lower energy costs [5]. An MEUF integrates a surfactant, which with sufficient dose self-assembles and forms micelles. The micelles then bind metal ions through electrostatic interactions and can be retained by an ultrafiltration membrane [6].
Most MEUF studies for nickel removal used the conventional one-factor-at-a-time method, namely, to examine one operational variable while fixing the others. For example, Karate and Marathe [7] examined the MEUF removal of nickel by testing a series of factors individually: flow rate, surfactant to metal (S/M) ratio, pH, feed metal ion concentration, pressure, and presence of electrolyte. Tanhaei, et al. [8] investigated the MEUF removal of nickel using single and mixed surfactants. Similarly, they determined the optimum SDS and nickel concentrations by separately examining the effect of SDS concentration, nickel concentration, pressure, and pH. Danis and Aydiner [9] examined the MEUF performance in four stages, by changing surfactant concentrations, nickel concentrations, transmembrane pressure, and electrolyte content separately.
System optimization is important in engineering applications because it is directly related to costs. In MEUF studies, it is crucial to find an optimal operating condition that yields high rejection and high permeate flux simultaneously, with minimal dosage of surfactant and power consumption. Though easy to conduct, the one-factor-at-a-time method tends to involve much labor and resources (many experimental runs) for a multi-variable system and does not provide adequate information on factor interactions or estimate the effects [10]. Besides, the method is difficult to find a true optimal condition with a reasonable number of experimental runs. These limitations can be avoided by using a more systematic experimental method, such as a response surface methodology (RSM). RSM is an experiment- and statistic-based technique that involves multiple factors and their interactions to optimize a process [11]. It has been increasingly used in environmental studies, such as to optimize the process condition for wastewater treatment [12,13,14,15,16,17,18]. However, only a few attempts of MEUF have been made to remove heavy metals using the RSM method, mostly focusing on copper, cadmium, and zinc [3,19]. Reports on RSM-based nickel removal using MEUF were rare.
Further, computer modeling can be integrated to describe a complex input-output relationship of a given system. Such approaches are suitable for uncertain or approximate reasoning when the systems are complex to describe with a mathematical model. In recent years, artificial neural network (ANN) has been developed to understand non-linear multi-variable systems [20]. It is a statistical modeling tool to approximate complex functions, especially nonlinear ones between system inputs and outputs, with few requirements for data. A typical ANN consists an input layer, one or more hidden layer with multiple hidden neurons, and an output layer. These layers are mathematically linked by weights and biases. ANN has become an attractive tool for non-linear modelling. It has been used in many fields of science and engineering [21,22,23] but extremely limited in MEUF studies. Table 1 summarizes the MEUF studies integrating RSM and other optimization models. The integration of ANN to RSM can provide additional information of the process behavior [24], but research efforts in MEUF are scarce (Table 1). These studies examined the removal of lead and zinc from cross-flow ultrafiltration systems, mostly conducted by the same researcher. Though the cross-flow operation could better scale-up to industrial application, most laboratory MEUF studies were carried out under batch operation. The removal of nickel ions from the common dead-end ultrafiltration system is desired.
This study examines the process of MEUF to remove nickel ions from dilute aqueous streams. The objectives are to (1) optimize MEUF process conditions using RSM, (2) predict the maximum nickel removal and flux rate under optimal conditions, and (3) verify RSM results using ANN modeling.

2. Materials and Methods

2.1. Materials

All chemicals were of analytical grade and were used as received. The anionic surfactant sodium dodecyl sulfate (SDS, 20% in H2O) was purchased from Sigma-Aldrich, Canada. Its properties are listed in Table 2. Nickel sulfate hexahydrate (NiSO4∙6H2O, J.T. Baker) were used as sources of metal ions. The pH of feed solutions was adjusted to 8 ± 0.1. Nickel reference standard solutions (1000 ppm ± 1%/certified) for Flame Atomic Absorption (FAA) tests were purchased from Fisher Scientific and diluted as needed. Distilled water was used in all experimental procedures. Permeate samples were collected and stored using sorption-free materials.

2.2. Dead-End Ultrafiltration Experiments

Batch experiments were conducted in a stirred ultrafiltration Cell (Amicon Model 8400, EMD Millipore) with a maximum volume uptake of 400 mL. Regenerated cellulose membrane (EMD Millipore, Canada) was used, with 3, 5, and 10 kDa MWCO (diameter of 76 mm and effective area of 0.00418 m2). An initial 250-mL feed solution was filled and continuously stirred (at a constant rate to get effective agitation and prevent membrane fouling) in each experimental run. All experiments were conducted at room temperature (23 ± 1 °C). The applied transmembrane pressure was controlled by pressurized nitrogen gas.
For each ultrafiltration run, 250 mL of feed solution was prepared with designated nickel and surfactant concentrations. When an ultrafiltration run starts, the first 10 mL sample was discarded, then every 20 mL permeate was sampled. The run was terminated when successive five samples were collected and timed. Nickel concentrations of the permeate samples were measured, and their permeate fluxes and rejection rates were determined. For both rejection rate and permeate flux, the average values of five permeate samples for each experimental run were calculated and used as inputs for RSM and ANN modeling. The membrane was cleaned after each run to recover its permeability (indicated by the flux rate of distilled water measured at 40 psi) and can be repeatedly used if over 90% of the original water flux (i.e., flux of distilled water passing the clean membrane at 40 psi) was recovered. Pretreatment of sampling apparatus, storage of samples, and recovery of membranes followed the procedures described by Lin et al. [29]. The experimental scheme and mechanism of MEUF are illustrated in Figure 1.

2.3. Sample and Data Analysis

The nickel concentration in permeate samples (CP) were measured using a Varian Model 55B SpectrAA FAA Spectrophotometer at 232.0 nm. The mean values of triplicate measurements for each permeate sample were calculated (%RSD ≤ 1.3%). FAA calibration curves were made before each set of measurement (R2 > 0.999).
To evaluate the removal efficiency of nickel ions using MEUF, the nickel rejection rate (R) and permeate flux (J) of the metal were calculated as follows:
R ( % ) = ( 1 C p C r ) × 100
where Cp and Cr denote the nickel concentration in the permeate and retentate, respectively. Cr was calculated using material balance.
J ( L / h · m 2 ) = V p t × A m
where V p is the volume of the permeate sample; t is the sampling time; and A m is the effective area of the membrane.

2.4. Response Surface Modeling

The RSM modeling and optimization consist the following steps: (1) statistical design of experiment, where all process variables vary simultaneously over experimental runs; (2) define coefficients of variables (and their interactions) in the mathematical model based on experimental results; (3) check the adequacy of the regressed model; and (4) predict the optimal experimental condition and responses using the model.

2.4.1. Design of Experiments

In this study, an RSM model based on Box–Behnken design (BBD) was used to optimize the four independent variables (factors) and to observe their effect on MEUF performance in terms of rejection rate and permeate flux. A BBD entails factors at high (+1), basic (0), and low (−1) levels. The center points (coded level 0 or the basic level), which were the midpoints between the high and low levels, were repeated multiple times. Table 3 presents the factors and levels set by the BBD. The design consists of 29 experimental runs, including 5 replicates of the central experiments to check the analysis repeatability and to estimate the experimental error. The responses (rejection rate and permeate flux) were determined experimentally according to designed runs. Design-Expert (version 11.1) was used for RSM modeling.

2.4.2. RSM Modeling

To determine the mathematical relationship between the responses and factors the following second-order polynomial equation was used to fit the experimental data obtained from the BBD experimental design. The response surface model includes the main, quadratic, and interactions terms:
Y = b 0 + i = 1 n b i X i + i = 1 n b i i X i 2 + i = 1 n 1 j = i + 1 n b i j X i X j
where Y is the predicted response; b0 the constant coefficient, bi the linear coefficients, bii the quadratic coefficients, bij the interaction coefficients; n the number of design variables, and Xi, Xj the coded levels of design variables.
Stepwise regression procedure was performed using backward elimination method to excludes non-significant terms (p-values > 0.05) from the initial response surface model. The regression coefficients of the reduced model are computed by multiple linear regression (MLR) method to minimize the sum of square of the residuals. The validity of the empirical model was tested using analysis of variance (ANOVA) at 95% confidence level. The fitted model was assessed by the R-squared (R2), the adjusted R-squared (R2-adj), and the predicted R-squared (R2-pre). The R2 value increases with the number of model terms, even when non-significant terms are added to the model. Therefore, the R2 value of a refined model is usually smaller than that of the full model. The R2-adjusted coefficient is used to adjust to the number of model terms, where the addition of non-significant terms usually decreases the R2-adjusted value. The predicted R-squared shows how well a model predicts responses for new observations. Based on the obtained response surface models, optimal conditions were determined by maximizing the nickel rejection and the permeate flux.

2.5. Artificial Neural Network (ANN) Modeling

The BBD design (factors and levels) and the corresponding responses were used to develop the ANN model, using the neural network toolbox for MATLAB 2016b. An ANN model with one hidden layer was trained to simulate nickel removal by MEUF. The inputs were pressure, feed nickel concentration, feed SDS concentration, and MWCO, whereas the outputs were rejection rate and permeate flux. All inputs and outputs were normalized into the range of (0, 1) to avoid putting too much weight on variables with a large variance. Twelve neurons in the hidden layer were optimized for the ANN model following Jing, et al. [30]. Datasets were randomly divided into training (70%), validation (15%), and testing (15%) subsets. The model was trained by minimizing the mean squared error (MSE) while maximizing the correlation coefficients (R) between the experimental and modeling outputs. The two outputs, i.e., rejection fate and permeate flux, were given equal weightings when calculating R2 for the ANN model. For comparison purposes, an inverse range scaling was performed on all modeling outputs to transfer them from (0, 1) to their original scales.

3. Results and Discussion

3.1. Ultrafiltration Experimental Results

Experimental results (i.e., nickel rejection rates and permeate flux) of the BBD design are reported in Table 4. The maximum rejection rate of nickel is 98.70% (flux = 23.03 L/h∙m2) in run 17, with a transmembrane pressure of 30 psi, nickel concentration of 1.25 mM, SDS concentration of 16.6 mM, and MWCO of 3 kDa. The maximum flux 178.28 L/h∙m2 (R = 91.83%) was found in run 26 with 50 psi pressure, 1.25 mM nickel, and 16.6 mM SDS using membrane MWCO of 10 kDa. It can be seen that higher rejection (or flux) tends to compromise on lower flux (or rejection), yet in practice high values of both rejection (indicates MEUF effectiveness) and flux (indicates efficiency) are desired. As such, an operating condition generating high rejection and flux is needed.

3.2. RSM Models

Goodness-of-fit of the regression model is evaluated using ANOVA by testing the significance of the regression model, significance of individual model coefficients, and lack-of-fit. For both rejection and flux models the assumptions for ANOVA are met, e.g., the residuals are normally and randomly distributed (figures not shown). Table 5 and Table 6 summarize the ANOVA analysis for rejection and flux, respectively, showing the goodness-of-fit of the quadratic models. The regression models for nickel rejection and permeate flux are statistically significant at the 95% confidence level in the studied range.
The significance of the model on rejection rate was determined by Fisher test, indicated by the F-value. The model F-value of 12.61 indicates that the model is significant, with a 0.01% chance that an F-value could occur due to noise. The lack-of-fit F-value of 3.46 indicates that there is a 12% chance that an F-value could occur due to noise. The calculated R2 (0.8486) and adjusted R2 (0.7813) was reasonably close to 1, showing good fitness of the regressed model. The difference between predicted R2 and the adjusted R2 is over 0.02. This may be due to the close values of the response (which can be sensitive to experimental and measurement errors).
The flux model shows great fitness. The model F-value of 2171.32 indicates that the model is highly significant, with only 0.01% chance that the value could occur due to noise. Non-significant lack-of-fit (p = 0.4387) also indicates good fitness of the model. Both R2 (0.9972) and adjusted R2 (0.9968) show good fitness of the regressed model. High predicted R2 (0.9958) indicates that the model can well predict response for new observations.
Regression model for nickel rejection:
Significant model terms (p < 0.05) are coded factors A, B, C, D, AD, BC, CD, and C². The reduced regression model (coded factors) for nickel rejection was determined as:
Rejection = 96.48 − 1.15 A − 0.99 B + 1.77 C − 1.17 D − 1.65 AD + 2.02 BC − 1.84 CD − 2.17 C²
where coded factor subject to the level of (−1,1).
The regressed model in terms of actual factors is:
Rejection = 80.35 + 0.19 Pressure − 6.69 CNi + 1.26 CSDS + 2.60 MWCO − 0.05 (Pressure)(MWCO) + 0.32 CNi CSDS − 0.06 (CSDS)(MWCO) − 0.03 (CSDS
where factors subjected to: 30 ≤ pressure ≤ 50 psi, 0.5 ≤ CNi ≤ 2 mM, 8.3 ≤ CSDS ≤ 24.9 mM, 3 ≤ MWCO ≤ 10 kDa. Equations (4) and (5) can be used to predict the nickel rejection for given levels of each factor.
The coefficients of coded factors indicate that the importance of the factor is in the order: BC > CD > C> AD > D A > B, i.e., interaction of nickel concentration and SDS concentration > interaction of SDS concentration and MWCO > SDS concentration > interaction of pressure and MWCO > MWCO pressure > nickel concentration.
Regression model for permeate flux:
Table 6 shows that A, D, AD, A², D² are significant model terms. The final equation in terms of coded factors is:
ln(Flux) = 3.89 + 0.22 A − 0.02 C + 0.80 D + 0.30 D²
where coded factor subject to the level of (−1,1).
Final equation in terms of actual factors is:
ln(Flux) = 2.59 + 0.02 Pressure − 0.002 CSDS − 0.09 MWCO + 0.024 (MWCO)²
where factors subjected to: 30 ≤ pressure ≤ 50 psi, 0.5 ≤ CNi ≤ 2 mM, 8.3 ≤ CSDS ≤ 24.9 mM, 3 ≤ MWCO ≤ 10 kDa.
The importance of factors is: D > A > C, i.e., MWCO > pressure > CSDS.

3.3. Effect of Factors on Rejection Rate and Permeate Flux

The response surface plots show the effect of pressure, nickel concentration, SDS concentration, and MWCO on rejection rate and permeate flux. The response surface and contour plot enable visualization of parameter interaction. Based on the ANOVA results, three interaction effect (i.e., pressure and MWCO, feed nickel and SDS concentration, feed SDS concentration and MWCO) on rejection rate and three individual effect (i.e., pressure, feed SDS concentration, and MWCO) on flux will be discussed.

3.3.1. Effect of Factors on Rejection

ANOVA results indicate significant interaction effect between pressure and MWCO on rejection. Figure 2 shows the effect of pressure and MWCO on rejection, when feed nickel and SDS concentrations were fixed at their central levels. Pressure seems to affect the rejection rate more at higher MWCO than the lower end. The rejection rate was relatively stable with the increase in pressure at MWCO of 3 kDa but considerably dropped at MWCO of 10 kDa (Figure 2a). Previous one-factor-at-a-time studies showed that pressure alone had a small effect on the rejection rate. For example, Huang, et al. [31] examined the rejection under a transmembrane pressure of 40 to 800 kPa and found that pressure did not significantly change the rejection rate. Mulligan et al. (2011) reported similar conclusions for a pressure range from 30 to 140 kPa. This observation can be explained, because the pressure does not affect the interaction between metal ions and the surfactant but mainly provides a driving force for mass transport across the membrane.
In terms of MWCO, smaller MWCO tend to generate higher rejection (Figure 2 and Figure 4). The observation is in agreement with previous findings. For example, Baek and Yang [32] reported higher chromate rejection (>99%) using membrane MWCO of 3 kDa than that of 10 kDa (98%). Bade and Lee [33] reported 98% rejection of chromate using the surfactant cetylpyridinium chloride (CPC) with 100 kDa membrane and 97% with 300 kDa membrane.
Figure 3 shows the predicted response under different metal and surfactant concentrations, when the pressure and MWCO were fixed as their central values. It can be seen that, at a lower SDS concentration (1 CMC), low metal concentrations result in high rejection rate. At higher nickel concentration, the decrease in rejection might be attributed to a lack of available binding sites. To sum up, MEUF is more efficient to treat dilute (i.e., low concentration) nickel streams, showing an advantage to traditional techniques (e.g., precipitation) that are inefficient at dilute streams. Alternatively, MEUF could be used as a secondary treatment method.
In the examined concentration range, higher SDS concentration resulted in higher rejection of nickel ions. When SDS concentration increased to approximately 20 mM no further increase in rejection is observed. Therefore, increasing the SDS feed concentration enhances the rejection of heavy metals until certain limits. The maximum rejection might be due to the competition between the surfactant sodium ions and nickel ions. The electrostatic interaction between the anionic micellar surface and nickel cations depends on the ion charge and concentration. At first, when increasing the SDS feed concentration, a higher fraction of surfactants will be in the micellar form. This will increase the surface charge; hence, more divalent nickel ions will be adsorbed on the micellar surface displacing the sodium ions. This ion exchange will consequently enhance heavy metal rejection. However, at low heavy metal feed concentration, when SDS concentration is further increased to concentration up to 20 mM, the sodium counter ion concentration might increase to an extent that the adsorption of sodium counter ions is favored. Therefore, no further increase in nickel rejection is achieved, as shown in Figure 3a and Figure 4a.

3.3.2. Effect of Factors on Flux

Pressure, SDS concentration, and MWCO significantly contribute to the flux rate. When the pressure was increased from 30 to 50 psi, the permeate flux rate increased. The permeate flux follows the Darcy’s law [34], i.e.,   J = L p × P where the membrane permeability   L p = 1 / ( η · R M ) , where η is the viscosity of the solution and R M is the membrane resistance. If the permeate flux linearly increases with pressure, the separation process is under the pressure controlled region, where the concentration polarization is negligible [4]. This linear relationship was observed in the present study (figure not shown), indicating that concentration polarization was not obvious, and that the membranes performed well.
SDS concentration was found negatively and linearly related to the flux rate. Increasing amount of SDS forms more SDS micelles which are retained by the membrane. The retained micelles may concentrate on the membrane surface or in its pores, hence reducing the permeate flux. The decrease in flux rate with the increase in surfactant concentration has been reported in literature [34]. In addition, nickel concentration seems to have little effect on flux. This can be explained by the small size of nickel ions which can easily pass the ultrafiltration membrane.
A higher MWCO (i.e., bigger pore size) of the membrane increases the permeate flux. Nonlinear relationship between MWCO and flux was observed. Flux rate gradually increases with MWCO in its lower ranges (3–7 kDa) and quickly increase at higher ranges (7–10 kDa).

3.4. RSM Optimization

The economic operation of the membrane processes draws attention to achieve lower costs in practice. As such, the MEUF process desires to use lower transmembrane pressures (minimize pressure) and lower dosages of surfactant (minimize CSDS) to treat large volumes of water (maximize flux rate), as well as obtaining a high efficiency in removing nickel ions (maximize rejection).
The optimal conditions of the MEUF of nickel were obtained using the desirability function approach in Design Expert. The condition was found by maximizing rejection and flux (defined by equations 5 and 7, respectively) when setting minimum pressure, CNi =1 mM, minimum CSDS, and 3 ≤ MWCO ≤ 10 kDa. The predicted maximum rejection rate (major response) and flux (secondary response) are 98.16% and 119.20 L/h∙m2, respectively, where pressure = 30 psi, CNi = 1.0 mM, CSDS = 10.05 mM, and MWCO = 10 kDa.

3.5. ANN Modeling

To predict the values of rejection rate and permeate flux using the ANN model, 75% of the data were randomly used for training purpose. The remainders were categorized as testing and validation data. In order to evaluate the ANN model, the model was presented with new values of rejection and flux that were not used during the training. The rejection and flux values estimated by ANN models were then compared with their corresponding actual values. The scatter regression plots of the ANN model shows the predicted values of rejection and flux against their experimental values (Figure 5 and Figure 6, respectively). Due to the inverse rescaling, two outputs, i.e., rejection rate and permeate flux, were first converted from (0, 1) and then plotted together within their orginal ranges. The best linear fit equations for the training, validation, tesing, and overall subsets mostly had a slope between 0.99 and 1, and the values of R2 were all higher than 0.99 (except for the testing values for rejection model, R2 = 0.719), indicating a close match between the experimental and modeling results. Therefore, the trained ANN model was able to accurately simulate the rejection rate and permeate flux for nickel removal process.

4. Conclusions

This paper demonstrated the feasibility of using BBD (an RSM design) to study the effect of process variables (pressure, nickel concentration, SDS concentration, and MWCO) on MEUF performance to remove nickel ions from aqueous solutions. RSM generated quadratic models for rejection and flux, respectively. The results showed that all factors significantly contribute to the rejection rate, namely the effectiveness of a MEUF process (importance follows: SDS concentration > MWCO pressure > nickel concentration). Pressure and MWCO are significant factors contributing to the permeate flux (importance: MWCO > pressure). Among the range of factors in the study, the optimal conditions to remove 1 mM nickel ions from aqueous solutions while obtaining the highest rejection (98.16%) and flux (119.20 L/h∙m2) are: pressure = 30 psi, CSDS = 10.05 mM, and MWCO = 10 kDa. Verification experiments showed that the quadratic models could adequately predict the MEUF performance. Furthermore, ANN modeling showed good model fitness to the experimental data. This study shows that RSM and ANN models could be used and provide information for the MEUF treatment of nickel-contaminated water. In future works, a cross-flow MEUF system will be used to better reflect the industrial practice. Furthermore, the recycle and reuse of SDS will be attempted to further reduce the capital cost.

Author Contributions

W.L. contributed to project conceptualization, experimental design, data collection, analysis and interpretation, modeling, validation, and drafted the manuscript; L.J. contributed to project conceptualization, experimental design, modeling, and manuscript review; B.Z. contributed to project conceptualization, project administration, funding acquisition, supervision, and manuscript revision. All authors have read and agreed to the published version of the manuscript.

Funding

This research was partially funded by Natural Sciences and Engineering Research Council of Canada (NSERC) and Canada Foundation for Innovation (CFI)

Acknowledgments

The authors would like to thank the Core Research Equipment and Instrument Training (CREAIT) Network of Memorial University for equipment support.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Fu, F.; Wang, Q. Removal of heavy metal ions from wastewaters: A review. J. Environ. Manag. 2011, 92, 407–418. [Google Scholar] [CrossRef]
  2. Tchounwou, P.B.; Yedjou, C.G.; Patlolla, A.K.; Sutton, D.J. Heavy Metal Toxicity and the Environment; Springer: Basel, Switzerland, 2012; pp. 133–164. [Google Scholar] [CrossRef] [Green Version]
  3. Xiarchos, I.; Jaworska, A.; Zakrzewska-Trznadel, G. Response surface methodology for the modelling of copper removal from aqueous solutions using micellar-enhanced ultrafiltration. J. Membr. Sci. 2008, 321, 222–231. [Google Scholar] [CrossRef]
  4. Landaburu-Aguirre, J.; Pongrácz, E.; Perämäki, P.; Keiski, R.L. Micellar-enhanced ultrafiltration for the removal of cadmium and zinc: Use of response surface methodology to improve understanding of process performance and optimisation. J. Hazard. Mater. 2010, 180, 524–534. [Google Scholar] [CrossRef] [PubMed]
  5. Tung, C.-C.; Yang, Y.-M.; Chang, C.-H.; Maa, J.-R. Removal of copper ions and dissolved phenol from water using micellar-enhanced ultrafiltration with mixed surfactants. Waste Manag. 2002, 22, 695–701. [Google Scholar] [CrossRef]
  6. Huang, J.; Yuan, F.; Zeng, G.; Li, X.; Gu, Y.; Shi, L.; Liu, W.; Shi, Y. Influence of pH on heavy metal speciation and removal from wastewater using micellar-enhanced ultrafiltration. Chemosphere 2017, 173, 199–206. [Google Scholar] [CrossRef] [PubMed]
  7. Karate, V.D.; Marathe, K. Simultaneous removal of nickel and cobalt from aqueous stream by cross flow micellar enhanced ultrafiltration. J. Hazard. Mater. 2008, 157, 464–471. [Google Scholar] [CrossRef]
  8. Tanhaei, B.; Pourafshari Chenar, M.; Saghatoleslami, N.; Hesampour, M.; Kallioinen, M.; Sillanpää, M.; Mänttäri, M. Removal of nickel ions from aqueous solution by micellar-enhanced ultrafiltration, using mixed anionic–non-ionic surfactants. Sep. Purif. Technol. 2014, 138, 169–176. [Google Scholar] [CrossRef]
  9. Danis, U.; Aydiner, C. Investigation of process performance and fouling mechanisms in micellar-enhanced ultrafiltration of nickel-contaminated waters. J. Hazard. Mater. 2009, 162, 577–587. [Google Scholar] [CrossRef]
  10. Czitrom, V. One-factor-at-a-time versus designed experiments. Am. Stat. 1999, 53, 126–131. [Google Scholar] [CrossRef]
  11. Montgomery, D.C. Design and Analysis of Experiments; John wiley & sons: Hoboken, NJ, USA, 2017. [Google Scholar]
  12. Ahmadi, M.; Vahabzadeh, F.; Bonakdarpour, B.; Mofarrah, E.; Mehranian, M. Application of the central composite design and response surface methodology to the advanced treatment of olive oil processing wastewater using Fenton’s peroxidation. J. Hazard. Mater. 2005, 123, 187–195. [Google Scholar] [CrossRef]
  13. Chavalparit, O.; Ongwandee, M. Optimizing electrocoagulation process for the treatment of biodiesel wastewater using response surface methodology. J. Environ. Sci. 2009, 21, 1491–1496. [Google Scholar] [CrossRef]
  14. Kiran, B.; Kaushik, A.; Kaushik, C.P. Response surface methodological approach for optimizing removal of Cr (VI) from aqueous solution using immobilized cyanobacterium. Chem. Eng. J. 2007, 126, 147–153. [Google Scholar] [CrossRef]
  15. Körbahti, B.K.; Aktaş, N.; Tanyolaç, A. Optimization of electrochemical treatment of industrial paint wastewater with response surface methodology. J. Hazard. Mater. 2007, 148, 83–90. [Google Scholar] [CrossRef] [PubMed]
  16. Sadri Moghaddam, S.; Alavi Moghaddam, M.R.; Arami, M. Coagulation/flocculation process for dye removal using sludge from water treatment plant: Optimization through response surface methodology. J. Hazard. Mater. 2010, 175, 651–657. [Google Scholar] [CrossRef]
  17. Wang, J.-P.; Chen, Y.-Z.; Ge, X.-W.; Yu, H.-Q. Optimization of coagulation–flocculation process for a paper-recycling wastewater treatment using response surface methodology. Colloids Surf. A Physicochem. Eng. Asp. 2007, 302, 204–210. [Google Scholar] [CrossRef]
  18. Zhu, X.; Tian, J.; Liu, R.; Chen, L. Optimization of Fenton and electro-Fenton oxidation of biologically treated coking wastewater using response surface methodology. Sep. Purif. Technol. 2011, 81, 444–450. [Google Scholar] [CrossRef]
  19. Landaburu-Aguirre, J.; García, V.; Pongrácz, E.; Keiski, R.L. The removal of zinc from synthetic wastewaters by micellar-enhanced ultrafiltration: Statistical design of experiments. Desalination 2009, 240, 262–269. [Google Scholar] [CrossRef]
  20. Desai, K.M.; Vaidya, B.K.; Singhal, R.S.; Bhagwat, S.S. Use of an artificial neural network in modeling yeast biomass and yield of β-glucan. Process Biochem. 2005, 40, 1617–1626. [Google Scholar] [CrossRef]
  21. Kasiri, M.; Aleboyeh, H.; Aleboyeh, A. Modeling and optimization of heterogeneous photo-fenton process with response surface methodology and artificial neural networks. Environ. Sci. Technol. 2008, 42, 7970–7975. [Google Scholar] [CrossRef]
  22. Yi, Q.-X.; Huang, J.-F.; Wang, F.-M.; Wang, X.-Z.; Liu, Z.-Y. Monitoring rice nitrogen status using hyperspectral reflectance and artificial neural network. Environ. Sci. Technol. 2007, 41, 6770–6775. [Google Scholar] [CrossRef]
  23. Zhang, K.; Zhang, B.; Chen, B.; Jing, L.; Zhu, Z.; Kazemi, K. Modeling and optimization of Newfoundland shrimp waste hydrolysis for microbial growth using response surface methodology and artificial neural networks. Mar. Pollut. Bull. 2016, 109, 245–252. [Google Scholar] [CrossRef] [PubMed]
  24. Balkin, S.D.; Lin, D.K. A neural network approach to response surface methodology. Commun. Stat. Theory Methods 2000, 29, 2215–2227. [Google Scholar] [CrossRef]
  25. Rahmanian, B.; Pakizeh, M.; Mansoori, S.A.A.; Esfandyari, M.; Jafari, D.; Maddah, H.; Maskooki, A. Prediction of MEUF process performance using artificial neural networks and ANFIS approaches. J. Taiwan Inst. Chem. Eng. 2012, 43, 558–565. [Google Scholar] [CrossRef]
  26. Rahmanian, B.; Pakizeh, M.; Mansoori, S.A.A.; Abedini, R. Application of experimental design approach and artificial neural network (ANN) for the determination of potential micellar-enhanced ultrafiltration process. J. Hazard. Mater. 2011, 187, 67–74. [Google Scholar] [CrossRef] [PubMed]
  27. Rahmanian, B.; Pakizeh, M.; Esfandyari, M.; Heshmatnezhad, F.; Maskooki, A. Fuzzy modeling and simulation for lead removal using micellar-enhanced ultrafiltration (MEUF). J. Hazard. Mater. 2011, 192, 585–592. [Google Scholar] [CrossRef] [PubMed]
  28. Jana, D.K.; Roy, K.; Dey, S. Comparative assessment on lead removal using micellar-enhanced ultrafiltration (MEUF) based on a type-2 fuzzy logic and response surface methodology. Sep. Purif. Technol. 2018, 207, 28–41. [Google Scholar] [CrossRef]
  29. Lin, W.; Jing, L.; Zhu, Z.; Cai, Q.; Zhang, B. Removal of heavy metals from mining wastewater by micellar-enhanced ultrafiltration (MEUF): Experimental investigation and monte carlo-based artificial neural network modeling. Water Air Soil Pollut. 2017, 228, 206. [Google Scholar] [CrossRef]
  30. Jing, L.; Chen, B.; Zhang, B. Modeling of UV-induced photodegradation of naphthalene in marine oily wastewater by artificial neural networks. Water Air Soil Pollut. 2014, 225, 1–14. [Google Scholar] [CrossRef]
  31. Huang, Y.; Wu, D.; Wang, X.; Huang, W.; Lawless, D.; Feng, X. Removal of heavy metals from water using polyvinylamine by polymer-enhanced ultrafiltration and flocculation. Sep. Purif. Technol. 2016, 158, 124–136. [Google Scholar] [CrossRef]
  32. Baek, K.; Yang, T.-W. Competitive bind of anionic metals with cetylpyridinium chloride micelle in micellar-enhanced ultrafiltration. Desalination 2004, 167, 101–110. [Google Scholar] [CrossRef]
  33. Bade, R.; Lee, S.-H. Chromate removal from wastewater using micellar enhanced ultrafiltration and activated carbon fibre processes; validation of experiment with mathematical equations. Environ. Eng. Res. 2008, 13, 98–104. [Google Scholar] [CrossRef]
  34. Kamble, S.; Marathe, K. Membrane characteristics and fouling study in MEUF for the removal of chromate anions from aqueous streams. Sep. Sci. Technol. 2005, 40, 3051–3070. [Google Scholar] [CrossRef]
Figure 1. A Schematic diagram of (a) SDS monomer, (b) SDS micelle (when SDS concentration > CMC) (c) micellar-enhance ultrafiltration (MEUF) setup, and (d) mechanism of MEUF removal of metal ions.
Figure 1. A Schematic diagram of (a) SDS monomer, (b) SDS micelle (when SDS concentration > CMC) (c) micellar-enhance ultrafiltration (MEUF) setup, and (d) mechanism of MEUF removal of metal ions.
Water 12 01269 g001
Figure 2. Response surface (a) and contour (b) showing the effect of pressure and MWCO on rejection rate. CNi = 1.25mM, CSDS = 16.6mM.
Figure 2. Response surface (a) and contour (b) showing the effect of pressure and MWCO on rejection rate. CNi = 1.25mM, CSDS = 16.6mM.
Water 12 01269 g002
Figure 3. Response surface (a) and contour (b) showing the effect of nickel and SDS concentrations on rejection rate. Pressure = 40 psi, MWCO = 5 kDa.
Figure 3. Response surface (a) and contour (b) showing the effect of nickel and SDS concentrations on rejection rate. Pressure = 40 psi, MWCO = 5 kDa.
Water 12 01269 g003
Figure 4. Response surface (a) and contour (b) showing the effect of SDS concentration and MWCO on rejection rate. Pressure = 40 psi, CNi = 1.25 mM.
Figure 4. Response surface (a) and contour (b) showing the effect of SDS concentration and MWCO on rejection rate. Pressure = 40 psi, CNi = 1.25 mM.
Water 12 01269 g004
Figure 5. The scatter plots of ANN model predicted values (rejection rate of nickel ions) versus experimental values for (a) training, (b) validation, (c) testing, and (d) all data sets.
Figure 5. The scatter plots of ANN model predicted values (rejection rate of nickel ions) versus experimental values for (a) training, (b) validation, (c) testing, and (d) all data sets.
Water 12 01269 g005
Figure 6. The scatter plots of ANN model predicted values (permeate flux) versus experimental values for (a) training, (b) validation, (c) testing, and (d) all data sets.
Figure 6. The scatter plots of ANN model predicted values (permeate flux) versus experimental values for (a) training, (b) validation, (c) testing, and (d) all data sets.
Water 12 01269 g006
Table 1. Summary of MEUF studies integrating RSM and other optimization models.
Table 1. Summary of MEUF studies integrating RSM and other optimization models.
SoluteUF System (Surfactant and Flow)RSM DesignIndependent VariablesOptimization ModelReferences
Pb2+SDS, cross-flowBBD (3 factors and 3 levels, 17 runs)CSDS (2–6 mM), S/M (5–15), pH (2–12)ANN and ANFIS[25]
Zn2+SDS and Brij-35, cross-flowFFD (7 factors, 22 runs)Pressure, pH, CSDS, S/M, L/M, CNaCl, Brij35/SDS ratioANN, R2 > 0.91[26]
Pb2+SDS, cross-flowBBD (3 factors, 3 levels)CSDS (2–6 mM), S/M (5–15), pH (2–12)Fuzzy logic models, R > 0.91[27]
Pb2+CTAB, cross-flowBBD (3 factors, 3 levels)CSDS (1.61–6.43 mM), S/M (5.64–13.8), pH (2.34–12.1)Interval type-2 fuzzy logic[28]
Surfactant: SDS = surfactant sodium dodecyl sulfate, CTAB = cetyltrimethylammonium bromide. Models: ANFIS = adaptive neuro-fuzzy inference system; ANN = artificial neural network; BBD = Box–Behnken Design; FFD = full factorial design.
Table 2. Properties of the surfactant used in this study.
Table 2. Properties of the surfactant used in this study.
PropertiesSpecifications
NameSodium dodecyl sulfate (SDS)
Chemical structure Water 12 01269 i001
Ionic typeAnionic
Molecular weight288.38 g/mol
Critical micellar concentration (CMC)8.2-8.3 mM
Table 3. Factors and levels set by Box–Behnken design (BBD).
Table 3. Factors and levels set by Box–Behnken design (BBD).
Factors Levels
Minimum (−1)Center (0)Maximum (+1)
(A) Pressure (psi)304050
(B) Ni concentration (mM)0.51.252
(C) SDS concentration (mM)8.316.624.9
(D) Molecular weight cut-off, or MWCO (kDa)35*10
* The center point 5 kDa MWCO was used instead of 6.5k Da due to the size of commercial membrane.
Table 4. Design layout and experimental results of BBD design.
Table 4. Design layout and experimental results of BBD design.
Std.RunFactor Input VariablesResponse Variable
Factor A
Pressure (psi)
Factor B
Ni conc. (Mm)
Factor C
SDS conc. (Mm)
Factor D
MWCO (kDa)
Rejection a (%)Flux a
(L/h∙m2)
131400.58.3594.8637.93
182501.258.3592.9845.15
253401.2516.6598.1336.83
74401.258.31094.30158.67
295401.2516.6597.0937.43
206501.2524.9598.1343.31
67401.2524.9397.1529.96
198301.2524.9598.1728.74
22940216.6397.9831.03
2310400.516.61097.76148.64
14114028.3588.0637.41
1012501.2516.6398.6738.25
31330216.6596.1529.27
2814401.2516.6596.5939.51
1115301.2516.61097.84115.56
2716401.2516.6596.3237.78
917301.2516.6398.7023.03
2618401.2516.6596.4736.45
819401.2524.91080.53 b149.23
42050216.6595.7045.36
221500.516.6595.0846.16
1722301.258.3591.3128.78
162340224.9598.2035.00
2124400.516.6398.4029.10
125300.516.6590.43 b30.19
1226501.2516.61091.83178.28
1527400.524.9596.9435.67
242840216.61093.53138.17
529401.258.3392.6128.96
a Rejection/flux values of an ultrafiltration run are the mean values of rejection/flux of all permeate samples (n=5) in that run. b Observed outliers; eliminated from analysis.
Table 5. ANOVA for reduced quadratic model (response: rejection).
Table 5. ANOVA for reduced quadratic model (response: rejection).
SourceSum of SquaresdfMean SquareF-Valuep-Value
Model158.79819.8512.61<0.0001significant
A-Pressure12.24112.247.780.0121
B-C-Ni10.51110.516.680.0187
C-C-SDS17.18117.1810.920.0039
D-MWCO13.18113.188.370.0097
AD12.10112.107.690.0125
BC16.26116.2610.330.0048
CD7.3317.334.660.0447
27.39127.3917.400.0006
Residual28.33181.57
Lack of Fit26.16141.873.460.1200not significant
Pure Error2.1640.5409
Cor Total187.1226
Fit statistics: R² = 0.8486, Adjusted R² =0.7813, Predicted R² = 0.4481. df = degree of freedom.
Table 6. ANOVA for reduced quadratic model (response: flux).
Table 6. ANOVA for reduced quadratic model (response: flux).
SourceSum of SquaresdfMean SquareF-Valuep-Value
Model10.6842.672173.32<0.0001significant
A-Pressure0.591410.5914481.41<0.0001
C-C-SDS0.003310.00332.670.1151
D-MWCO7.6517.656229.95<0.0001
0.383210.3832311.97<0.0001
Residual0.0295240.0012
Lack of Fit0.0256200.00131.310.4387not significant
Pure Error0.003940.0010
Cor Total10.7128
Data were transformed into natural log. R² = 0.9972, Adjusted R² = 0.9968, Predicted R² = 0.9958.

Share and Cite

MDPI and ACS Style

Lin, W.; Jing, L.; Zhang, B. Micellar-Enhanced Ultrafiltration to Remove Nickel Ions: A Response Surface Method and Artificial Neural Network Optimization. Water 2020, 12, 1269. https://doi.org/10.3390/w12051269

AMA Style

Lin W, Jing L, Zhang B. Micellar-Enhanced Ultrafiltration to Remove Nickel Ions: A Response Surface Method and Artificial Neural Network Optimization. Water. 2020; 12(5):1269. https://doi.org/10.3390/w12051269

Chicago/Turabian Style

Lin, Weiyun, Liang Jing, and Baiyu Zhang. 2020. "Micellar-Enhanced Ultrafiltration to Remove Nickel Ions: A Response Surface Method and Artificial Neural Network Optimization" Water 12, no. 5: 1269. https://doi.org/10.3390/w12051269

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