Next Article in Journal / Special Issue
A Simple Method to Identify the Dominant Fouling Mechanisms during Membrane Filtration Based on Piecewise Multiple Linear Regression
Previous Article in Journal / Special Issue
Diffusion Dialysis for Acid Recovery from Acidic Waste Solutions: Anion Exchange Membranes and Technology Integration
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Operating Conditions Optimization via the Taguchi Method to Remove Colloidal Substances from Recycled Paper and Cardboard Production Wastewater

by
Mayko Rannany S. Sousa
,
Jaime Lora-García
,
María-Fernanda López-Pérez
*,
Asunción Santafé-Moros
and
José M. Gozálvez-Zafrilla
Research Institute for Industrial, Radiophysical and Environmental Safety (ISIRYM) Universitat Politècnica de València (UPV), Plaza Ferrándiz y Carbonell s/n, 03801 Alcoy, Spain
*
Author to whom correspondence should be addressed.
Membranes 2020, 10(8), 170; https://doi.org/10.3390/membranes10080170
Submission received: 19 May 2020 / Revised: 10 July 2020 / Accepted: 22 July 2020 / Published: 29 July 2020
(This article belongs to the Special Issue Membranes for Water and Wastewater Treatment)

Abstract

:
Optimization of the ultrafiltration (UF) process to remove colloidal substances from a paper mill’s treated effluent was investigated in this study. The effects of four operating parameters in a UF system (transmembrane pressure (TMP), cross-flow velocity (CFV), temperature and molecular weight cut-off (MWCO)) on the average permeate flux (Jv), organic matter chemical oxygen demand (COD) rejection rate and the cumulative flux decline (SFD), was investigated by robust experimental design using the Taguchi method. Analysis of variance (ANOVA) for an L9 orthogonal array were used to determine the significance of the individual factors, that is to say, to determine which factor has more and which less influence over the UF response variables. Analysis of the percentage contribution (P%) indicated that the TMP and MWCO have the greatest contribution to the average permeate flux and SFD. In the case of the COD rejection rate, the results showed that MWCO has the highest contribution followed by CFV. The Taguchi method and the utility concept were employed to optimize the multiple response variables. The optimal conditions were found to be 2.0 bar of transmembrane pressure, 1.041 m/s of the cross-flow velocity, 15 °C of the temperature, and 100 kDa MWCO. The validation experiments under the optimal conditions achieved Jv, COD rejection rate and SFD results of 81.15 L·m−2·h−1, 43.90% and 6.01, respectively. Additionally, SST and turbidity decreased by about 99% and 99.5%, respectively, and reduction in particle size from around 458–1281 nm to 12.71–24.36 nm was achieved. The field-emission scanning electron microscopy images under optimal conditions showed that membrane fouling takes place at the highest rate in the first 30 min of UF. The results demonstrate the validity of the approach of using the Taguchi method and utility concept to obtain the optimal membrane conditions for the wastewater treatment using a reduced number of experiments.

1. Introduction

The pulp and paper (P & P) industry is ranked as the world’s third largest consumer of fresh water [1] and an important producer of wastewater with different organic and inorganic contaminants. Depending on the type of processes used in paper manufacture, the integration between production and environmental protection is one of the key topics in the paper industry.
According to the Confederation of European Paper Industries (CEPI) [2], Europe is the second largest producer of paper and paperboard with 22.7% (91.39 million tons) of world production, making it one of the most important industries in the European economic sector. The paper industries hold an important place in Spanish economy, as Spain is one of the European leaders in paper recycling, with 84% of the raw materials used by the paper industry containing recovered paper [2]. However, we cannot forget that water is, also, an essential raw material for manufacturing paper and paperboard, and effluent treatment is a critical part of the process [3]. In order to minimize the amount of freshwater used and the volume of effluent discharged, the European Commission has described the best available techniques to be adopted by the P & P industry [4].
A number of conventional processes have previously been used to treat the different types of paper mill wastewater including coagulation and flocculation [5], adsorption [6,7], advanced oxidation [8,9] and membrane filtration [10,11,12,13]. It is important to mention that paper mills have their own wastewater treatment plants but some water treatment methods typically used for the P & P are not environmentally efficient: e.g.: coagulation/flocculation using inorganic coagulants create disposal problems and conventional aerobic processes have not been efficient in the removal of color or recalcitrant compounds [14,15]. This inadequacy can also make it impossible to reuse water in the papermaking process. Therefore, factories must improve their treatment plants to achieve the pollutant loadings permissible under current regulations and/or to reuse their process water.
Membrane separation technology has been attracting increasing attention as an alternative method for the post-treatment of paper mill wastewater. Some processing methods, such as nanofiltration, ultrafiltration (UF) and reverse osmosis, have recently been used in paper mills to purify secondary and tertiary effluents using external biological treatment [7,16,17,18]. Major advantages of membrane separation processes are their scalability, low installation costs and easiness of operation. However, their technical and economic liability must be carefully assessed for each specific process.
Ultrafiltration is an attractive process for paper mill wastewater treatment and it can be used as an advanced tertiary treatment to remove suspended solids and dissolved and colloidal substances (DCS) during the treatment of paper industry effluent in order to facilitate the reuse of the treated wastewater and reduce fresh water consumption [19,20]. What makes it so attractive is that most of the pollutants consist of high-molecular-weight compounds and these are easily removed by UF [10,21].
However, membrane fouling is still a limiting factor for the adoption and use of UF on a large scale in paper manufacturing applications. This fouling results in a sharp decline in permeate flux and, thus, changes in membranes selectivity [10,22,23]. Membrane fouling also increases the process cost due to repeated plant shutdowns to clean and wash the membranes [24]. Previous studies have shown that the main foulants of the membranes used for paper industry wastewater are DCS including fatty acids, resin acids, lignins and trace amounts of sterols, steryl esters and triglycerides [19,25]. Currently, this treatment technology can only be used to filter paperboard mill effluent that has been pre-treated and that still does not meet discharge standards [26].
Statistical experimental design incorporating design of experiments (DoE) techniques can be used to investigate the effects of all the possible interactions between the factors at one time, while undertaking the fewest possible experiments. A review of the literature revealed that an increasing number of studies are being conducted using DoE approaches in the membrane technology field to optimize operating conditions [10,27,28,29,30,31,32]. The DoE approaches for robust design include the Taguchi method that combines mathematical and statistical techniques to arrive at a special design of experiments with an orthogonal array (OA) to study multiple factors with a small number of experiments. This saves time and money by reducing the number of experiments required in the investigation [16]. It is worth mentioning that this approach is becoming popular because it is easy to adopt and applies an efficient method for optimizing the operating parameters.
This approach also allows studying the influence of each individual factor on the response variables, as well as on the effects of interactions between factors over the response variables, that is to say, all operational conditions varying simultaneously according to design array. This permits the factors that have the greatest and least influence to be determined, along with the optimal level for each factor in an OA [33]. In addition, in many UF approaches it is necessary to consider the application of multiple response optimization, because the process performance is often evaluated using several quality characteristics (responses). In this case, the Taguchi method and utility concept are useful tools for optimizing operating parameters in multiple characteristics responses [34].
The statistical analysis of variance (ANOVA) can be used to provide information on whether the operating parameters (factors) are statistically significant or not, as well as to identify the influence of individual factors and establish the relationships between the factors and operating conditions. In this analysis, the p-value index is used to know which operating parameters have a significant effect on the response variables. This information is complemented with the F-test which helps to identify whether or not each factor is significant at the selected confidence level [35,36]. In this study, ANOVA was also used to analyze the experimental results.
The aim of this work was to determinate the effect of operating conditions such as transmembrane pressure (TMP), cross-flow velocity (CFV), temperature and molecular weight cut-off (MWCO) on the average permeate flux, chemical oxygen demand (COD) rejection rate and cumulative flux decline (SFD), in addition to determining the optimal conditions for the given sets of values and to find the best response variables by using Taguchi experimental design and the utility concept. The results of this study may be used as a guideline when operating UF systems under the best conditions in a wastewater treatment plant (WWTP) in a papermaking factory. The filtration results and analysis of the experimental data presented and discussed in this study were carried out by using ANOVA to find the significance of the controlling factors and optimized using the Taguchi method to find the optimal operating conditions. A standard L9 orthogonal array was selected for experimental planning with four factors and three levels for each factor.

2. Materials and Methods

2.1. Paperboard Mill Treated Effluent Feedstock

The paperboard mill treated effluent (PMTE) used in this work came from a secondary clarifier effluent from a WWTP in a papermaking factory located in the south of the Valencian autonomous region in Spain. In order to prevent early membrane fouling, remove large suspended solids, and reduce initial turbidity and COD in the PMTE, the raw feed solution was pre-filtered by conventional filtration (low-pressure pump at around 1 bar) with a Cintropur®® NW 50 filter element, and centrifugal propeller and filter cloths with a 5 μm nominal pore size. The significant characteristics of the PMTE samples are listed in Table 1.

2.2. Membranes and Experimental Setup

This study used polyethersulfone (PES) membranes provided by Synder Filtration™ (Vacaville, CA, USA) with MWCO of 10, 50 and 100 kDa, denoted 10-ST, 50-MQ and 100-LY, respectively.
The experiments were performed in a typical UF pilot plant with a flat-sheet membrane module (Rhône-Poulenc, Lyon, France), that allowed working with two membranes with similar or different MWCO (depending on the experiment being carried out). The effective area for each membrane in the module was 154.8 cm2. The details of the experimental set-up have been described previously by Sousa et al. [21]. The pilot plant had a data acquisition system (temperature, module input and output pressure) from LabVIEW® (National Instruments, Austin, TX, USA). The permeate was collected during the filtration in a beaker placed on an electronic balance connected to a computer in order to continuously register the weighting data. This data was then automatically logged every 30 s and subsequently used to calculate permeate flux through the membranes.

2.3. Analytical Methods

The PMTE used as the feed solution and the UF permeate samples were analyzed according to the methods described below. The suspended solids analyses were carried out in accordance with the standard methods [37]. Turbidity was measured using a Dinko 112 turbidimeter (ASTM D1889, Barcelona, Spain). Conductivity was measured using a WTW level 3 conductivity device (ASTM D1125–14, Weilheim in Oberbayern, Germany). COD and total nitrogen in the effluent were analyzed using a Merck photometer (Merck KGaA, Darmstadt, Germany) and a Merck TR-300 thermoreactor (Merck KGaA, Darmstadt, Germany) in accordance with the standard methods.

2.4. Field-Emission Scanning Electron Microscopy (FESEM)

Field-emission scanning electron microscopy (FESEM) measurements were used to provide information on the fouling that formed on the membranes. The surface and cross-section morphologies of the fresh and fouled membranes were observed by field-emission scanning electron microscope, (ZEISS ULTRA 55, Oxford Instruments, Berkshire, UK), operated with a voltage of 200 kV and an accelerating voltage of 0.02–5 kV. Before analysis, the dried membrane samples were attached to double-sided adhesive carbon tape on an aluminum holder, and subsequently coated with a thin layer of gold prior to analysis.

2.5. Experimental Procedure

2.5.1. Membrane Characterization

Before the UF runs, permeability experiments were carried out to determine the intrinsic membrane resistance (Rm). Distilled water was used as the feed solution and measurements were taken for 1.0, 1.5, 2.0, and 3.0 bar of transmembrane pressure (TMP) at 1.041 m/s and 22.5 °C, in total recirculation mode to generate a quasi-steady state. The characterization process was undertaken for an operation time of 2 h to stabilize the flux through the membrane during this time. Previously each membrane was worked under compaction conditions with pure water at 5 bar for 1 h, in order to obtain a stable membrane structure. Rm values were calculated using the resistance model (Equation (1)), where under this condition, there was no membrane fouling resistance (Rf) [38]:
J P = T M P μ   R t = T M P μ   ( R m + R f )
where JP is the permeate flux, μ is the viscosity of the permeate stream and Rt is the total membrane resistance. Rf can be understood as the result of the sum of the three main fouling mechanisms: (i) pore blockage resistance when colloids block the membrane pores, (ii) adsorption resistance as a result of foulant adsorption inside or over the membrane and, (iii) cake layer resistance as a consequence of the accumulation of particles on the membrane [35,39,40].

2.5.2. Ultrafiltration Experiments

The UF experiments were performed in crossflow filtration mode. The studied parameters were varied in the following ranges: TMP (1.0–3.0 bar), CFV (0.463–1.041 m/s), MWCO (10–100 kDa) and temperature (15–30 °C). These values were selected based on the operational limits of the experimental setup, industrial scale-up and economic considerations [21].
The evolution of permeate flux was gravimetrically measured at different time intervals according to Equation (2),
J P = 1 A m   ρ d m p d t
where Am is the effective membrane area, mp is the total mass of permeate, ρ is the water density, and t is the filtration time.
In order to keep the feed concentration constant, both the permeate and the retained streams were continuously recirculated to the feed tank.
To evaluate the UF performance in terms of permeability the average permeate flux was calculated by the following equation [41]:
J ¯ P = 1 t M t 1 t M J P ( t ) · d t
where JP(t) is the permeate flux evolution over time determined by regression analysis on the experimental data, t1 is the initial time operation (first data collected), and tM is 2 h.
To analyze the effect of the operating conditions on UF resistance, the cumulative flux decline (SFD) [41] was calculated from the following equation:
S F D = i = 1 M J P ( 0 ) J P ( i ) J P ( 0 )
where M is the number of experimental data collected and, JP (0) is the initial permeate flux measured at t 1 .
This parameter gives information about how the flux declines over the duration of the experiment (not just the difference between the initial and final permeate flux). Therefore, the higher the SFD value, the faster and more noticeable is the flux decline, indicating that the membrane fouling is more severe.
The average flux decline index ( F D ¯ ) provides information on the decrease of feed permeate flux throughout the experiment estimated as follows:
F D ¯ = [ 1 t M t 1 t M J P ( 0 ) J P ( t ) J P ( 0 ) d t ] × 100
To evaluate the UF efficiency in removing organic matter, the COD rejection rate was chosen as the response variable calculated as:
R ( % ) = ( 1 C p C f ) × 100
where: Cp and Cf are the COD concentration in the permeate and feed, respectively.

2.6. Experimental Design Based on the Taguchi Method

An experimental design based on the Taguchi method was used to design the experiments. The Taguchi method applies fractional experimental designs, called orthogonal arrays (OA), to reduce the number of experiments required to determine the optimum conditions based on the results [29,30,42]. One of the important steps in the Taguchi approach is the appropriate selection of OAs, which depends on the number of control factors and their levels. The minimum number of experimental trails required in an OA is given by Nmin = (L − 1)·F + 1, where F and N are the number of factors and levels respectively [36,43].
As mentioned above, the four factors (parameters) chosen were the transmembrane pressure, the cross-flow velocity, the temperature and the MWCO of the membrane; and three response variables were analyzed: the average permeate flux, the COD rejection and the cumulative flux decline. The selected factors, their designated symbols and levels are presented in Table 2.
For an experimental design with four factors and three levels for each factor, an L9 (34) orthogonal array was selected. In this case, 27 runs were conducted (three repetitions at each trial condition). The design of the experiments planning matrix for the L9 array [43] is shown in Table 3. The experiments were carried out in a randomized order. New membrane samples were preconditioned and used in the experimental runs (27 membrane samples). In this way, the lurking effect of possible irreversible fouling was avoided and the intrinsic variability of the membrane material is included in the replications.
The aim of this DoE was to determine the operating parameters (factors) under which the average permeate flux and COD rejection rate achieve their maximum values, and the SFD achieved its minimum value. The Taguchi method uses a statistical measure of the process performance, called signal-to-noise ratio (S/N), which depends on the criterion for the response variable to be optimized. The S/N ratios are divided into three different categories and data sets, the larger-the better, the smaller-the-better and the nominal-the-better [16,29]. In this study, the system was optimized when the average permeate flux and COD rejection rate were as large as possible (Equation (7)), and the SFD was as small as possible (Equation (8)):
The   larger the better   ( S / N ) = 10   log ( 1 n i = 1 n 1 Y i 2 )
The   smaller the better   ( S / N ) = 10   log ( 1 n i = 1 n Y i 2 )
where Yi the response variable at each experiment and n is the total number of repetitions in a trial.
The sequence of steps to be followed using the Taguchi method to optimize the UF process is shown in Figure 1.
Minitab Statistical and Statgraphics Centurion XVII Software were used to analyze the Taguchi experiments and optimize the operating conditions.

2.6.1. Utility Concept

The implementation of the utility concept in the Taguchi method helps to obtain the best combination of operating parameters to optimize multiple response S/N ratios (MRSN) simultaneously by differentiating the relative importance (weights) of various responses [34,45,46]. In this work, it is assumed that the overall utility is the sum of the responses of the individual utilities and it can be written as [34,47]:
U ( x 1 , x 2 , x n ) = f   [ U 1 ( x 1 ) , U 2 ( x 2 ) , U n ( x n ) ]
where U(x1,x2,…xn) is the overall utility of n response variables and Ui(xi) is the utility index of ith response.
The response variables can be attributed priorities depending upon the process goals to be achieved. The priorities can be adjusted by providing a weight to the individual utility index. Therefore, by assigning weights to the response variables, the overall utility function can be expressed as:
U ( x 1 , x 2 , x n ) = i = 1 n W i U i ( x i )
where Wi is the weight assigned to the ith response variable.
It is worth noting that the assignment of weights is a purely subjective (empirical) step and depends on each experiment or process that will be carried out [48]. Therefore, in this paper, the analytic hierarchy process (AHP) method, developed by [49] was used to determine the associate weight criteria for each response variable in the multiple optimization required to calculate the overall utility index. The relative normalized weight Wi of each criterion is calculated using the AHP geometric mean method GMi on the rows in the pairwise comparison matrix, A║aij║and it can be calculated from the follow equation [49]:
W i = G M i i = 1 M G M i
where i,j = 1, 2, …, M and M is the number of factors in judgement matrix A.
aij = 1 for i = j,
aij = 1/aij for i ≠ j.
In addition, the total sum of the weight for all the responses must be assigned to hold the following condition:
i = 1 n W i = 1
For this study, as stated above, the objective was to maximize permeate flux and COD rejection rate, and minimize the SFD, simultaneously. From the utility concept, the MRSN of the overall utility value is given by Equation (13):
μ M R S N = W J ¯ P μ 1 + W C O D _ r e j e c t i o n   μ 2 + W S F D μ 3
where:
μ 1 = 10 log ( 1 J P ¯ 2 )
μ 2 = 10 log ( 1 C O D r e j e c t i o n 2 )
μ 3 = 10 log ( S F D 2 )
WJP, WCOD_rejection, and WSFD are the weights assigned to the permeate flux, COD rejection rate and SFD.
It is worth mentioning that the utility function is of the higher-the-better type. If the composite measure (the overall utility) is maximized, the quality characteristics considered for the evaluation of utility will automatically be optimized (maximized or minimized) [44].

2.6.2. Optimal Performance Prediction

Once the optimal level of the operating conditions has been selected, it is possible to predict and verify the utility responses using the optimal parameters. The predicted response values under optimal conditions Yopt can be calculated from Equation (17) [50,51]:
Y opt = m + j = 1 p [ ( m i , j ) m ]
where m is the overall mean value of μ MRSN over nine trials; mi,j is the mean value of the quality response under optimal conditions; and p is the number of significant operating parameters that affect the UF process.
The 95% confidence interval for the confirmation experiments (CICE) must be evaluated at the selected error level according to the following expression [43,50,51]:
CI CE = ± F α ( 1 , f e ) · M S e ( 1 n e f f + 1 R )
where Fα (1,fe) is the F-ratio at a confidence level of (1 − α) against a DOF equal to one and an error degree of freedom fe and, neff is the effective sample size calculated as:
n e f f = N 1 + ( DOF   of   all   factors   used   to   estimate   the   mean )
where N is the number of experiments, R is the number of repetitions, and MSe is the error variance.

2.7. Analysis of Variance (ANOVA)

In order to determine the relative importance of the factors, ANOVA was employed by calculating the sum of squares (SS), degrees of freedom (DOF), mean of square (MS), associated F-test of significance (F) and percentage contribution (P%) [52].

3. Results and Discussions

3.1. Experimental Results

As previously mentioned, permeability tests were carried out prior to each Taguchi experiment with pure water as the feed, in order to determine the intrinsic resistance of the membrane for each membrane used, as illustrated in Figure 2.
The specific resistance values obtained from the permeability test for the membranes of 10, 50 and 100 kDa were 3.46 × 1012, 4.56 × 1012 and 9.88 × 1012 m−1, respectively.
The values for the average permeate flux, COD rejection rate and the cumulative flux decline caused by membrane fouling for each trial experiment according to the Taguchi design are shown in Table 4. The highest average permeate flux was obtained in Trial 8 ( J ¯ P = 95.16 L·m−2·h−1) and the lowest value was obtained in Trial 1 ( J ¯ P = 15.23 L−1·h−1). The corresponding average flux decline indices were respectively F D ¯ = 44.87% and a F D ¯ = 14.0%.
For some trials, significant differences on the results can be observed between replicates. The reason can be found in the inhomogeneous behavior between samples taken from the same sheet. The replication used in the experimental method aims to diminish the effect of the membrane variability on the response.
Figure 3 shows the evolution of the permeate average flux for all the trials carried out according the Taguchi design of Table 3. It can be seen that for the same MWCO, the permeate flux decreased with increasing TMP in all trials because of the membrane fouling. An increase in the TMP leads to higher accumulation of colloidal substances on the membrane surface and pore blocking. At the beginning of the process, flux declined very quickly, mainly for the trials corresponding to higher TMP and MWCO (Figure. 4c), possibly because of the membrane pores becoming blocked more rapidly by adsorption and the accumulation of colloidal substances. Afterwards, the permeate flux continued to decline due to the growth of a cake layer on the membrane surface, until the permeate flux reached a quasi-stationary state [53,54,55,56].

3.2. Taguchi Results

The corresponding S/N ratio (in dB) calculated for the response variables at each trial are listed in Table 5.
In order to analyze the influence of each factor on the response variable, the S/N ratio for a single factor can be determined by averaging the S/N ratios at their levels. The range of the effect for each factor (Δs) is calculated as the difference between the two readings, the higher the range, the stronger the effect of the factor, in other words, it shows which parameter has the greatest effect on the response.
The mean S/N ratio curves for each factor are shown in Figure 4. It is worth mentioning that the peak points in these plots correspond to the optimal condition.
As can be seen in Figure 4, the variations (Δs) around the mean S/N value were different for the different factors. TMP and membrane MWCO had the greatest effect on the average permeate flux as they have the steepest slope, Δs = 7.10 and 5.95 respectively. CFV was the next one with a Δs = 1.42, and temperature had the lowest variation around the mean S/N value, with Δs = 0.78. In addition, from Table 5, the overall mean value was calculated as 33.59 (dB) from all the trial experiment results. It can be observed that the increase in J ¯ P was stronger when the TMP changes from 1.0 to 2.0 bar than when it changes over the range from 2.0 to 3.0 bar, this could be due to the effects of polarization and cake compaction on the membrane surface. For CFV, the slope of the line between the different levels is not the same (0.463–0.752 m/s is higher than at 0.752–1.041 m/s), but with a small variation around the J ¯ P value. Also, it can be seen for MWCO and temperature that the slopes from 10 to 100 kDa and 15 °C to 30 °C (respectively) are almost the same. Therefore, the maximum average permeate flux can be obtained for 3.0 bar, 100 kDa, 1.041 m/s CFV and high temperature (30 °C).
Under optimal COD rejection conditions, a positive and larger value of S/N is desired. In Figure 4b when comparing the S/N between different factors, it was shown that the most significant variation around the mean S/N ratios is observed for MWCO and CFV (Δs = 2.07 and 1.32 respectively). Also, it can be seen that the S/N ratio increased with TMP and CFV and decreased with MWCO. Hence maximum COD removal occurred at higher TMP and CFV (3.0 bar and 1.041 m/s), and 10 kDa. It is worth mentioning that the DCS found in the PMTE are a mixture of high and low molecular weight organic and inorganic compounds, thus the contribution of the smaller particles gives lower rejection during high MWCO UF in membranes with 50 and 100 kDa MWCO.
Figure 4c shows that an increase in TMP, temperature and MWCO caused a decrease in the S/N ratio for SFD, that is to say, these factors intensified the membrane fouling effects. On the other hand, an increase in CFV induced an increase in the S/N ratio, this resulted in a decrease in the fouling effect. The highest variations around the mean S/N ratio were found for MWCO and TMP (Δs = 5.21 and 4.81). Generally, the permeate flux increased with increasing MWCO and TMP. However, under these operating conditions, DCS in PMTE can easily pass through the membrane and blocking can be observed within the pores and on the membrane surface. In addition, the highest S/N ratio for the SFD factor (−8.98 ± 0.28) was achieved in Trial 1, whereas the lowest S/N ratio (−19.43 ± 0.29) was obtained in Trial 8. The optimal conditions that minimized the SFD (lowest level of fouling) were obtained at the lowest TMP (1.0 bar), highest CFV (1.041 m/s), at temperature 15 °C and at the smallest MWCO (10 kDa).

3.3. ANOVA Results

A statistical analysis of variance (ANOVA) was carried out to quantitatively determine the effect of each factor on the UF process indicators, with the aim of estimating whether the process parameters are statistically significant or not on the results responses. The ANOVA results are shown in Table 6.
In order to determine the qualitative significance of each factor on the responses, Fisher’s test (F-value) was employed in the ANOVA analysis. An F-value is defined as the ratio of variance due to the effect of a factor on the variance due to the inherent error in the system [57]. The F-value was compared to the critical F-value (Fcr) [52]. A calculated F-value lower than the Fcr-value means that the effect of that factor is not significant at the selected confidence level or/and it is not important in comparison with the error term. In this study, with four factors, three levels for each factor and three repetitions at each trial condition, the DOF for each factor is 2 and the DOF for the error is 18, so the Fcr-value at a confidence level of 95% is equal to 3.55. In accordance with the ANOVA table, for average permeation flux, the F-value for TMP and MWCO (114.1 and 81.52, respectively) are greater than the Fcr-value. This means that the variance of these factors is significant compared with the variance of error and they have a significant effect on the response. On the other hand, temperature and CFV had no meaningful qualitative effect on J ¯ P , as their F-values were less than the Fcr-value. Furthermore, COD rejection rate and SFD presented F-values for all factors greater than the Fcr-value which means that the effect of these factors is significant at the 95% confidence level and they have a meaningful qualitative effect on responses.
Another statistical tool that is helpful for qualitative evaluation in ANOVA is the p-value, which is used to indicate which factors had a significant effect on the responses. The smaller the p-value at an α level of significance, the more significant is the corresponding factor [58,59]. In this study, based on p-values at the 95% confidence level ( α = 0.05), all the factors had a statistically significant (p-value < 0.05) effect on the COD rejection rate and SFD. For J ¯ P , CFV and temperature had a p-value higher than 0.05, thus the effect may be regarded as insignificant and it can be ignored.
The use of the percentage contribution (P%) in ANOVA analysis is helpful for the quantitative evaluation of the factorial effects of the performance indicators. The percentage contributions P% of all factors on average permeate flux, COD rejection rate, and SFD are shown in Figure 5. TMP (P% = 54.62) was the most important factor on average permeate flux, as higher pressure resulted in higher permeate flux, according to Darcy’s law. MWCO (P% = 39.06) was the second most important factor, followed by temperature and CFV. For the COD rejection rate, the order of importance for the factors is as follows MWCO > CFV > temperature > TMP. In addition, MWCO and TMP (46.57% and 40.89%, respectively) were the most significant parameters on membrane fouling resistance, followed by the CFV and temperature. TMP and MWCO were the most important factors for responses. Higher TMP and MWCO resulted in higher permeate flux. However, more intensive flux decline, due to membrane fouling, occurred at higher permeate flux.
It is important to mention that the values reported due to error resulting from uncontrollable noises should be below 50% for the results to be reliable [28,29,60]. Therefore, it can be seen in Figure 5 that for average permeate flux, COD rejection rate and SFD, the error values are 4.31%, 3.07%, and 3.81%, respectively. This means that the error values for the experiment are not significant for the UF process.

3.4. Optimal Results Obtained from the Taguchi Method and Utility Concept

The aim of optimizing the process was to find the operating conditions that led to a maximum average permeate flux and COD rejection rate based on the levels that gave the highest S/N ratios for the factors (desirable values) and to minimize the SFD, that is, the levels that gave the smallest S/N ratios (adverse values).

3.4.1. Analysis of Individual Response Optimization

After identifying the optimal operating conditions, the optimal responses were predicted individually using the Taguchi method and ANOVA. Table 7 shows the Taguchi prediction results for the optimal conditions for average permeate flux, COD rejection rate and SFD.
According to the Taguchi predictions, the average permeate flux at TMP 3.0 bar, CFV 0.752 m/s, at 22.5 °C and with a 100 kDa MWCO, achieves 81.20 L·m−2·h−1. The COD rejection rate predicted under optimal conditions indicates a 57.92% rejection, higher than any value obtained in the DoE combinations. For the SFD under optimal conditions estimated by the Taguchi method, the minimum SFD predicted is approximately 1.80, equivalent to a F D ¯ of 8.65%. Therefore, we can see that the values of the three response variables combined are far from those values obtained experimentally (see Table 4).

3.4.2. Analysis of Multi-Response Optimization

As mentioned above, in order to determine the weight for each response variable, a pairwise comparison matrix was compiled using the AHP method as presented in Table 8.
Thus, the weights assigned to response variables were WJP = 0.568, WCOD Rejection = 0.098 and WSFD = 0.334. The consistency ratio index (CR) is used to evaluate the consistency of AHP estimates. In this case, it was calculated as 0.021, which should be less than the allowed value of CR = 0.1, this means that the pairwise comparison matrix was considered acceptable.
The overall utility index for the μMRSN was calculated using Equation (13) with values associated with the weights of each response, using the lager the better (S/N) and the results are presented in Table 9.
ANOVA analysis was also performed for the multiple response variables using the utility concept. From Table 10 it is clear that when F-value is compared with Fcr (3.55), TMP, MWCO and CFV had a qualitatively significant effect (at a confidence level of 95%) on MRSN. The percentage contributions extracted from the ANOVA table were also used to determine the significance of each operating parameter in the process. The P% values were arranged as follows: TMP > MWCO > CFV > Temperature. Therefore, according to the results, TMP and MWCO were the most important factors in optimizing the multi–response UF system.
The optimal operating conditions for the simultaneous response were obtained based on the criteria that both J ¯ P and COD rejection rate must be maximized and SFD should be minimized. The variation in the overall utility for the operating parameters at different levels is presented in Figure 6.
It is clear from Figure 6 that the optimal combination of operating conditions (maximum value of the overall utility) was found at the second level of transmembrane pressure (2.0 bar), the second level of cross-flow velocity (1.041 m/s), the first level of temperature (15 °C), and third level of MWCO (100 kDa).
Once the optimal levels had been selected the next step was to estimate the multi-response S/N ratio and predict the optimal values for the simultaneous optimization response, calculated using Equation (17) and presented in Table 11.

3.5. Confirmation Experiment under Optimal Conditions

After determining the optimal operating conditions for the overall utility value and the significance of factors, validation experiments (for multi-responses) were carried out at the optimal levels in order to validate the predicted UF responses suggested using the Taguchi method with utility concept [34,51].
The observed permeate flux results as a function of time under optimized conditions during the UF of PMTE are plotted in Figure 7. As described previously, the flux decline was mainly the result of two phenomena, pore blocking and cake layer formation, which mostly occurred in the first hour of the process [21,29]. During the first 30 min of the UF, the flux decreased by 22.63%. Furthermore, immediately after pore blockage, the permeate flux continued to decline due to the formation and growth of a cake layer until the system approached the quasi-steady state. At the end of the process (after 2 h), the final permeate flux was about 67.0 L·m−2·h−1 and flux decline was around 39.72%, which confirms that the membrane fouling took place with a higher rate in the first 30 min and at a slower rate when the system had achieved a steady state. Therefore, the observed experimental values of average permeate flux and cumulative flux decline were about 81.15 L·m−2·h−1 and 6.01 (SFD equivalent to a F D ¯ value of 28.96%).
The total resistance ( R t ) at the end of the 2-h experiment under optimal conditions was 1.13 × 1013 m−1, which is sum of the intrinsic membrane resistance (Rm = 3.40 × 1012 m−1) and fouling resistance (Rf = 7.91 × 1012 m−1).
In addition, the membrane surface morphologies were observed by FESEM. Figure 8 shows the images of the membrane (PES 100 kDa) before and after the UF experiments were carried out. As can clearly be seen, before UF there is no blocking on the pores and no cake layer on the membrane surface. Figure 8b,c shows the surface of the membrane fouling after 30 min and after 2 h. In both cases, the images show the existence of pore blocking due to DCS adsorption within the membrane pores and sediments deposited on the surface (cake layer) acting to resist the UF [61]. Figure 8d shows the morphologies of the fouling sediments on the membrane.
Furthermore, it can be seen that the membrane was indeed fouled after 30 min filtration. However, the FESEM images of the membranes after filtration, at 2 h, were highly similar to the membranes after 30 min filtration. Therefore, it may be concluded that the permeate flux decline might result from the pore blocking as opposed to the formation and growth of the cake layer on the membrane.
Additionally, to verify the permeate applicability for paper mill reuse, the physical and chemical properties of the treated effluent obtained under optimal operational conditions (PES 100 kDa membrane at TMP = 2.0 bar, CFV = 0.752 m/s, and T = 15 °C) were compared with treated paper mil effluent used in this study as feed solution (Table 1). The results obtained for the physical-chemical parameters are given in Table 12. From the results obtained, all properties showed high retention efficiencies and proved the effectiveness of the UF under optimal conditions.
The optimum predicted results at the 95% confidence interval calculated using Equation (18) and the observed experimental results for the response variables are given in Table 13.
The observed multi-response of the overall utility falls within the 95% confidence interval for the optimal range of the response variables. In addition, it is clearly observed that the deviation between predicted and experimental results is very small, which confirms that the Taguchi method and utility concept can be used to predict the multi–response UF for any parametric combination, while individual optimization don’t got good predictions.

4. Conclusions

In this study, the Taguchi method, utility concept and ANOVA analysis were used as statistical tools to investigate the effects and significance of four operating parameters and to optimize the UF process with respect to average permeate flux, COD rejection rate and cumulative flux decline.
ANOVA was used to determine the most significant factors affecting the response variables. From the percentage contribution, the order of importance of each factor in maximum J ¯ P was TMP > MWCO > T > CFV; for maximum COD rejection rate it was MWCO > CFV > T > TMP; and to achieve the minimum SFD: MWCO > TMP > CFV > T.
The optimal UF operating parameters, based on the Taguchi method and utility concept, were found at TMP (2.0 bar), CFV (1.041 m/s), temperature (15 °C) and MWCO (100 kDa). Under these optimal conditions, J ¯ P , COD rejection rate and SFD resistance of 81.15 L·m−2·h−1, 43.90% and 6.01 (around and F D ¯ value of 28.96 %), respectively, were obtained and they were within of the predicted range at the 95% confidence interval.
Measurements of turbidity, COD and particle size in the permeate showed a significant decrease 3.21 to 0.0002 NTU, 146 mg/L to 81.8 mg/L and 458–1281 nm to 12.71–24.36 nm, respectively, which confirms a substantial reduction in colloidal compounds. Therefore, it can be said that UF is suitable for removing dissolved and colloidal substances from wastewater effluents from recycled paperboard manufacturing.
Finally, we can say that the Taguchi method and utility allow membrane conditions for the P & P wastewater treatment to be optimized using a reduced number of experiments. The methodology used in this study could be used as a guideline for operating UF systems applied as a tertiary treatment for paperboard mill treated effluents under optimal conditions.

Author Contributions

Conceptualization, J.L.-G., M.-F.L.-P., M.R.S.S., A.S.-M. and J.M.G.-Z.; methodology, M.R.S.S., A.S.-M. and J.M.G.-Z.; software, M.R.S.S., A.S.-M. and J.M.G.-Z.; validation, J.L.-G., M.-F.L.-P., M.R.S.S., A.S.-M. and J.M.G.-Z.; formal analysis, J.L.-G., M.-F.L.-P., M.R.S.S., A.S.-M.; investigation, M.R.S.S. resources, M.R.S.S.; data curation, J.L.-G., M.-F.L.-P., M.R.S.S., A.S.-M. and J.M.G.-Z.; writing—original draft preparation, M.R.S.S.; writing—review and editing, M.-F.L.-P., M.R.S.S., A.S.-M. and J.M.G.-Z.; visualization, M.R.S.S., A.S.-M. and J.M.G.-Z.; supervision, J.L.-G. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

The authors of this work wish to gratefully acknowledge the assistance received from the UPV through its research program, and the Research Institute for Industrial, Radiophysical and Environmental Safety (ISIRYM).

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Sevimli, M. Post-Treatment of Pulp and Paper Industry Wastewater by Advanced Oxidation Processes. Ozone Sci. Eng. Ozone Sci. Eng. 2005, 27, 37–43. [Google Scholar] [CrossRef]
  2. CEPI. Key Statistics Report 2017|CEPI-CONFEDERATION OF EUROPEAN PAPER INDUSTRIES. Available online: http://www.cepi.org/keystatistics2017 (accessed on 20 September 2018).
  3. Rajkumar, K. An Evaluation of Biological Approach for the Effluent Treatment of Paper Boards Industry-An Economic Perspective. J. Bioremediation Biodegrad. 2016, 7. [Google Scholar] [CrossRef] [Green Version]
  4. Suhr, M.; Klein, G.; Kourti, I.; Gonzalo, M.R.; Santonja, G.G.; Roudier, S.; Sancho, L.D. Best Available Techniques (BAT) Reference Document for the Production of Pulp, Paper and Board. Eur. Comm. 2015, 906. [Google Scholar] [CrossRef]
  5. Ahmad, A.; Wong, S.; Teng, T.; Zuhairi, A. Improvement of Alum and PACl Coagulation by Polyacrylamides (PAMs) for the Treatment of Pulp and Paper Mill Wastewater. Chem. Eng. J. 2008, 137, 510–517. [Google Scholar] [CrossRef]
  6. Temmink, H.; Grolle, K. Tertiary Activated Carbon Treatment of Paper and Board Industry Wastewater. Bioresour. Technol. 2005, 96, 1683–1689. [Google Scholar] [CrossRef] [PubMed]
  7. Zhang, Q.; Chuang, K.T. Adsorption of Organic Pollutants from Effluents of a Kraft Pulp Mill on Activated Carbon and Polymer Resin. Adv. Environ. Res. 2001, 5, 251–258. [Google Scholar] [CrossRef]
  8. Catalkaya, E.C.; Kargi, F. Advanced Oxidation Treatment of Pulp Mill Effluent for TOC and Toxicity Removals. J. Environ. Manag. 2008, 87, 396–404. [Google Scholar] [CrossRef]
  9. Pérez, M.; Torrades, F.; García-Hortal, J.A.; Domènech, X.; Peral, J. Removal of Organic Contaminants in Paper Pulp Treatment Effluents under Fenton and Photo-Fenton Conditions. Appl. Catal. B Environ. 2002, 36, 63–74. [Google Scholar] [CrossRef]
  10. Gönder, Z.B.; Arayici, S.; Barlas, H. Treatment of Pulp and Paper Mill Wastewater Using Utrafiltration Process: Optimization of the Fouling and Rejections. Ind. Eng. Chem. Res. 2012, 51, 6184–6195. [Google Scholar] [CrossRef]
  11. Liua, G.; Liu, Y.; Shib, H.; Qianb, Y. Treatability of Krafi Spent Liquor by Microfiltration and Ultrafiltration. Desalination 2004, 11, 131–141. [Google Scholar] [CrossRef]
  12. Nuortila-Jokinen, J.; Mänttäri, M.; Huuhilo, T.; Kallioinen, M.; Nyström, M. Water Circuit Closure with Membrane Technology in the Pulp and Paper Industry. Water Sci. Technol. A J. Int. Assoc. Water Pollut. Res. 2004, 50, 217–227. [Google Scholar] [CrossRef]
  13. Zaidi, A.; Buisson, H.; Sourirajan, S.; Wood, H. Ultra- and Nano-Filtration in Advanced Effluent Treatment Schemes for Pollution Control in the Pulp and Paper Industry. Water Sci. Technol. 1992, 25, 263–276. [Google Scholar] [CrossRef]
  14. Toczyłowska-Mamińska, R. Limits and perspectives of pulp and paper industry wastewater treatment—A review. Renew. Sustain. Energy Rev. 2017, 78, 764–772. [Google Scholar] [CrossRef]
  15. Kamali, M.; Khodaparast, Z. Review on recent developments on pulp and paper mill wastewater treatment. Ecotoxicol. Environ. Saf. 2015, 114, 326–342. [Google Scholar] [CrossRef]
  16. Beril Gönder, Z.; Arayici, S.; Barlas, H. Advanced Treatment of Pulp and Paper Mill Wastewater by Nanofiltration Process: Effects of Operating Conditions on Membrane Fouling. Sep. Purif. Technol. 2011, 76, 292–302. [Google Scholar] [CrossRef]
  17. Neves, L.C.; Souza, J.B.D. Pulp and paper mill effluent post-treatment using microfiltration and ultrafiltration membranes. Cell. Chem. Technol. 2017, 51, 579–588. [Google Scholar]
  18. Shukla, S.K.; Kumar, V.; Van Doan, T.; Yoo, K.; Kim, Y.; Park, J. Combining activated sludge process with membrane separation to obtain recyclable quality water from paper mill effluent. Clean Technol. Environ. Policy 2015, 17, 781–788. [Google Scholar] [CrossRef]
  19. Chen, C.; Mao, S.; Wang, J.; Bao, J.; Xu, H.; Su, W.; Dai, H. Application of Ultrafiltration in a Paper Mill: Process Water Reuse and Membrane Fouling Analysis. BioResources 2015, 10, 2376–2391. [Google Scholar] [CrossRef] [Green Version]
  20. Krawczyk, H.; Oinonen, P.; Jönsson, A.-S. Combined Membrane Filtration and Enzymatic Treatment for Recovery of High Molecular Mass Hemicelluloses from Chemithermomechanical Pulp Process Water. Chem. Eng. J. 2013, 225, 292–299. [Google Scholar] [CrossRef]
  21. Sousa, M.R.S.; Lora-Garcia, J.; López-Pérez, M.-F. Modelling Approach to an Ultrafiltration Process for the Removal of Dissolved and Colloidal Substances from Treated Wastewater for Reuse in Recycled Paper Manufacturing. J. Water Process Eng. 2018, 21, 96–106. [Google Scholar] [CrossRef]
  22. Karthik, M.; Dhodapkar, R.; Manekar, P.; Aswale, P.; Nandy, T. Closing water loop in a paper mill section for water conservation and reuse. Desalination 2011, 281, 172–178. [Google Scholar] [CrossRef]
  23. Mänttäri, M.; Nuortila-Jokinen, J.; Nyström, M. Evaluation of Nanofiltration Membranes for Filtration of Paper Mill Total Effluent. Filtr. Sep. 1997, 34, 275–280. [Google Scholar] [CrossRef]
  24. Cassano, A.; Conidi, C.; Drioli, E. Comparison of the Performance of UF Membranes in Olive Mill Wastewaters Treatment. Water Res. 2011, 45, 3197–3204. [Google Scholar] [CrossRef] [PubMed]
  25. Puro, L.; Tanninen, J.; Nyström, M. Analyses of Organic Foulants in Membranes Fouled by Pulp and Paper Mill Effluent Using Solid-Liquid Extraction. Desalination 2002, 143, 1–9. [Google Scholar] [CrossRef]
  26. Puro, L.; Kallioinen, M.; Mänttäri, M.; Nyström, M. Evaluation of Behavior and Fouling Potential of Wood Extractives in Ultrafiltration of Pulp and Paper Mill Process Water. J. Membr. Sci. 2011, 368, 150–158. [Google Scholar] [CrossRef]
  27. Hesampour, M.; Krzyzaniak, A.; Nyström, M. The Influence of Different Factors on the Stability and Ultrafiltration of Emulsified Oil in Water. J. Membr. Sci. 2008, 325, 199–208. [Google Scholar] [CrossRef]
  28. Pourjafar, S.; Jahanshahi, M.; Rahimpour, A. Optimization of TiO2 Modified Poly(Vinyl Alcohol) Thin Film Composite Nanofiltration Membranes Using Taguchi Method. Desalination 2013, 315, 107–114. [Google Scholar] [CrossRef]
  29. Reyhani, A.; Sepehrinia, K.; Shahabadi, S.M.S.; Rekabdar, F.; Gheshlaghi, A. Optimization of Operating Conditions in Ultrafiltration Process for Produced Water Treatment via Taguchi Methodology. Desalin. Water Treat. 2015, 54, 2669–2680. [Google Scholar] [CrossRef]
  30. Rezvanpour, A.; Roostaazad, R.; Hesampour, M.; Nyström, M.; Ghotbi, C. Effective Factors in the Treatment of Kerosene–Water Emulsion by Using UF Membranes. J. Hazard. Mater. 2009, 161, 1216–1224. [Google Scholar] [CrossRef]
  31. Salahi, A.; Abbasi, M.; Mohammadi, T. Permeate Flux Decline during UF of Oily Wastewater: Experimental and Modeling. Desalination 2010, 251, 153–160. [Google Scholar] [CrossRef]
  32. Salahi, A.; Mohammadi, T. Oily Wastewater Treatment by Ultrafiltration Using Taguchi Experimental Design. Water Sci. Technol. Lond. 2011, 63, 1476–1484. [Google Scholar] [CrossRef] [PubMed]
  33. Ezzati, A.; Gorouhi, E.; Mohammadi, T. Separation of Water in Oil Emulsions Using Microfiltration. Desalination 2005, 185, 371–382. [Google Scholar] [CrossRef]
  34. Kaladhar, M.; Subbaiah, K.V.; Rao, C.S.; Rao, K.N. Application of Taguchi Approach and Utility Concept in Solving the Multi-Objective Problem When Turning AISI 202 Austenitic Stainless Steel. J. Eng. Sci. Technol. Rev. 2011, 4, 55–61. [Google Scholar] [CrossRef]
  35. Mohammadi, T.; Safavi, M.A. Application of Taguchi Method in Optimization of Desalination by Vacuum Membrane Distillation. Desalination 2009, 249, 83–89. [Google Scholar] [CrossRef]
  36. Park, S. Robust Design and Analysis for Quality Engineering; Springer US: New York, NY, USA, 1996. [Google Scholar]
  37. APHA. Standard Methods for the Examination of Water and Wastewater, 20th ed.; American Public Health Association: Washington, DC, USA, 1998. [Google Scholar]
  38. Cheryan, M. Ultrafiltration and Microfiltration Handbook; CRC Press: Washington, DC, USA, 1998. [Google Scholar]
  39. Khan, M.T.; Takizawa, S.; Lewandowski, Z.; Jones, W.L.; Camper, A.K.; Katayama, H.; Kurisu, F.; Ohgaki, S. Membrane Fouling Due to Dynamic Particle Size Changes in the Aerated Hybrid PAC–MF System. J. Membr. Sci. 2011, 371, 99–107. [Google Scholar] [CrossRef] [Green Version]
  40. Rezaei, H.; Ashtiani, F.Z.; Fouladitajar, A. Fouling Behavior and Performance of Microfiltration Membranes for Whey Treatment in Steady and Unsteady-State Conditions. Braz. J. Chem. Eng. 2014, 31, 503–518. [Google Scholar] [CrossRef] [Green Version]
  41. Cojocaru, C.; Zakrzewska-Trznadel, G. Response Surface Modeling and Optimization of Copper Removal from Aqua Solutions Using Polymer Assisted Ultrafiltration. J. Membr. Sci. 2007, 298, 56–70. [Google Scholar] [CrossRef]
  42. Idris, A.; Ismail, A.F.; Noordin, M.Y.; Shilton, S.J. Optimization of Cellulose Acetate Hollow Fiber Reverse Osmosis Membrane Production Using Taguchi Method. J. Membr. Sci. 2002, 205, 223–237. [Google Scholar] [CrossRef] [Green Version]
  43. Phadke, M.S. Quality Engineering Using Robust Design; Prentice Hall: Englewood Cliffs, NJ, USA, 1989. [Google Scholar]
  44. Roy, R.K. A Primer on the Taguchi Method; Society of Manufacturing Engineers: Southfield, MI, USA, 1990. [Google Scholar]
  45. Barua, P.B.; Kumar, P.; Gaindhar, J.L. Optimal Setting of Process Parameters for Multi-Characteristic Products Using Taguchi Design and Utility-Conceptan Approach. Proc. ICAMIE Univ. Roorkee (India) 1997, 1997, 839–842. [Google Scholar]
  46. Goyal, S.; Singh Kandra, V.; Yadav, P. Experimental study of turning operation and optimization of MMR and surface roughness using Taguchi method. Int. J. Innov. Res. Adv. Eng. 2016. [Google Scholar] [CrossRef]
  47. Bunn, D.W. Analysis for Optimal Decisions; John Wiley & Sons: New York, NY, USA, 1982. [Google Scholar]
  48. Kumar, Y.; Singh, H. Multi-Response Optimization in Dry Turning Process Using Taguchi’s Approach and Utility Concept. Procedia Mater. Sci. 2014, 5, 2142–2151. [Google Scholar] [CrossRef] [Green Version]
  49. Saaty, T.L. Analytic Hierarchy Process; McGraw Hill Higher Education: New York, NY, USA, 1980. [Google Scholar]
  50. Kansal, H.K.; Singh, S.; Kumar, P. Performance Parameters Optimization (Multi-Characteristics) of Powder Mixed Electric Discharge Machining (PMEDM) through Taguchi’s Method and Utility Concept. Indian J. Eng. Mater. Sci. 2006, 29, 81–94. [Google Scholar]
  51. Ross, P.J. Taguchi Techniques for Quality Engineering: Loss Function, Orthogonal Experiments, Parameter and Tolerance Design; McGraw Hill Professional: Washington, DC, USA, 1996. [Google Scholar]
  52. Montgomery, D.C. Design and Analysis of Experiments; John Wiley & Sons: Hoboken, NJ, USA, 2008. [Google Scholar]
  53. Song, L. Flux Decline in Crossflow Microfiltration and Ultrafiltration: Mechanisms and Modeling of Membrane Fouling. J. Membr. Sci. 1998, 139, 183–200. [Google Scholar] [CrossRef]
  54. Mulder, J. Basic Principles of Membrane Technology, 2nd ed.; Springer: Dordrecht, The Netherlands, 1996. [Google Scholar]
  55. Xu, J.; Chang, C.-Y.; Gao, C. Performance of a Ceramic Ultrafiltration Membrane System in Pretreatment to Seawater Desalination. Sep. Purif. Technol. 2010, 75, 165–173. [Google Scholar] [CrossRef]
  56. Seyed Shahabadi, S.M.; Reyhani, A. Optimization of Operating Conditions in Ultrafiltration Process for Produced Water Treatment via the Full Factorial Design Methodology. Sep. Purif. Technol. 2014, 132, 50–61. [Google Scholar] [CrossRef]
  57. Ennil Köse, T. Agricultural Residue Anion Exchanger for Removal of Dyestuff from Wastewater Using Full Factorial Design. Desalination 2008, 222, 323–330. [Google Scholar] [CrossRef]
  58. Cochran, W.G. Experimental Designs, 2nd ed.; WILEY: New York, NY, USA, 1992. [Google Scholar]
  59. Myers, R.H.; Montgomery, D.C.; Anderson-Cook, C.M. Response Surface Methodology: Process and Product Optimization Using Designed Experiments, 4th ed.; Wiley: Hoboken, NJ, USA, 2016. [Google Scholar]
  60. Gönder, Z.B.; Kaya, Y.; Vergili, I.; Barlas, H. Optimization of Filtration Conditions for CIP Wastewater Treatment by Nanofiltration Process Using Taguchi Approach. Sep. Purif. Technol. 2010, 70, 265–273. [Google Scholar] [CrossRef]
  61. Reyhani, A.; Hemmati, M. Wastewater Treatment by Ultrafiltration System, Considering the Effects of Operating Conditions: Experimental and Modeling. Desalin. Water Treat. 2014, 52, 6282–6294. [Google Scholar] [CrossRef]
Figure 1. Flow diagram of Taguchi method steps to optimize a UF process to remove dissolved and colloidal substances (DCS) from paperboard mill treated effluent [43,44].
Figure 1. Flow diagram of Taguchi method steps to optimize a UF process to remove dissolved and colloidal substances (DCS) from paperboard mill treated effluent [43,44].
Membranes 10 00170 g001
Figure 2. Volumetric flux as a function of transmembrane pressure for polyethersulfone (PES) membranes of different MWCO (T = 22.5 °C).
Figure 2. Volumetric flux as a function of transmembrane pressure for polyethersulfone (PES) membranes of different MWCO (T = 22.5 °C).
Membranes 10 00170 g002
Figure 3. Evolution of the permeate flux through the operating time for each MWCO: (a) 10 kDa, (b) 50 kDa, and (c) 100 kDa (values are the average of the three replicates).
Figure 3. Evolution of the permeate flux through the operating time for each MWCO: (a) 10 kDa, (b) 50 kDa, and (c) 100 kDa (values are the average of the three replicates).
Membranes 10 00170 g003
Figure 4. Mean effect curves for S/N ratios for (a) the average permeate flux, (b) COD rejection rate, and (c) the cumulative flux decline (SFD).
Figure 4. Mean effect curves for S/N ratios for (a) the average permeate flux, (b) COD rejection rate, and (c) the cumulative flux decline (SFD).
Membranes 10 00170 g004aMembranes 10 00170 g004b
Figure 5. ANOVA results for the percentage contribution of each factor to the response processes.
Figure 5. ANOVA results for the percentage contribution of each factor to the response processes.
Membranes 10 00170 g005
Figure 6. Effect of operating parameters on mean utility value ( J ¯ P COD rejection rate, SFD).
Figure 6. Effect of operating parameters on mean utility value ( J ¯ P COD rejection rate, SFD).
Membranes 10 00170 g006
Figure 7. Permeate flux as a function of time under optimized conditions during UF of paper mill wastewater: PES 100 kDa membrane at TMP = 2.0 bar, CFV = 0.752 m/s and T = 15 °C.
Figure 7. Permeate flux as a function of time under optimized conditions during UF of paper mill wastewater: PES 100 kDa membrane at TMP = 2.0 bar, CFV = 0.752 m/s and T = 15 °C.
Membranes 10 00170 g007
Figure 8. Field-emission scanning electron microscope (FESEM) image of fresh and fouled membranes (PES 100 kDa) at different operating times, (a) clean membrane surface, (b) membrane surface fouled after 30 min filtration, (c) at the end of the UF (2 h) with pore blocking and cake layer (d) membrane foulant sediments.
Figure 8. Field-emission scanning electron microscope (FESEM) image of fresh and fouled membranes (PES 100 kDa) at different operating times, (a) clean membrane surface, (b) membrane surface fouled after 30 min filtration, (c) at the end of the UF (2 h) with pore blocking and cake layer (d) membrane foulant sediments.
Membranes 10 00170 g008
Table 1. Average compositions of the ultrafiltration (UF) feed solution (prefiltered solution from a secondary clarifier effluent from a wastewater treatment plant (WWTP)).
Table 1. Average compositions of the ultrafiltration (UF) feed solution (prefiltered solution from a secondary clarifier effluent from a wastewater treatment plant (WWTP)).
ParameterUnitsValue
Suspended solids (TSS)g·L−10.046 ± 0.01
TurbidityNTU3.21 ± 0.5
ConductivitymS·cm−14.20 ± 1.0
CODmg·L−1146 ± 5.0
Total nitrogenmg·L−10.8 ± 0.01
pH8.30 ± 0.5
Pre-filtered by conventional filtration (pretreatment).
Table 2. Operating parameters and their levels.
Table 2. Operating parameters and their levels.
ParametersLabelsLevels
L1L2L3
TMP (bar)A1.02.03.0
CFV (m/s)B0.4630.7521.041
Temperature (°C)C15.022.530.0
MWCO (kDa)D1050100
MWCO = Molecular weight cut-off.
Table 3. Experimental layout using L9 (34) orthogonal array in accordance with the Taguchi method.
Table 3. Experimental layout using L9 (34) orthogonal array in accordance with the Taguchi method.
Experimental Trial (*)Levels
ABCD
11111
21222
31333
42123
52231
62312
73132
83213
93321
* All experiments were carried out in a randomized run.
Table 4. Experimental responses of the Taguchi orthogonal array L9 (34) for the three repetitions of each trial (R1, R2 and R3).
Table 4. Experimental responses of the Taguchi orthogonal array L9 (34) for the three repetitions of each trial (R1, R2 and R3).
Trial n°Factors (Parameters)Responses
TMP
(bar)
CFV
(m/s)
T
(°C)
MWCO
(kDa)
J ¯ P  
(L·m−2·h−1)
COD
Rejection (%)
SFD
R1R2R3R1R2R3R1R2R3
11.00.46315.01015.2325.7415.9746.2548.946.852.922.732.80
21.00.75222.55029.4732.8529.1838.8535.8334.234.904.593.99
31.01.04130.010045.1542.0747.7941.8341.7742.856.195.986.18
42.00.46322.510072.0469.6774.3036.9234.2333.758.077.687.06
52.00.75230.01049.5439.0241.4350.4250.9651.964.295.973.99
62.01.04115.05057.7967.8153.6947.6946.2546.853.774.224.56
73.00.46330.05059.1459.5973.5740.4641.3740.688.297.948.46
83.00.75215.010082.5090.5595.1643.7543.8542.059.739.169.21
93.01.04122.51048.1150.3249.3052.553.4655.385.654.815.46
Table 5. Signal-to-noise results (mean ± standard deviation for the three repetitions at each trial).
Table 5. Signal-to-noise results (mean ± standard deviation for the three repetitions at each trial).
Trial n°S/N Ratio
J ¯ P   COD RejectionSFD
124.88 ± 2.5233.50 ± 0.25−8.98 ± 0.28
229.65 ± 0.5731.16 ± 0.56−13.13 ± 0.94
333.03 ± 0.5532.49 ± 0.12−15.73 ± 0.16
437.14 ± 0.2830.85 ± 0.42−17.64 ± 0.59
532.60 ± 1.0934.17 ± 0.13−13.67 ± 1.86
635.41 ± 1.0433.43 ± 0.14−12.46 ± 0.82
736.01 ± 1.0832.22 ± 0.10−18.31 ± 0.28
839.01 ± 0.6732.71 ± 0.20−19.43 ± 0.29
933.84 ± 0.1934.61 ± 0.24−14.52 ± 0.74
Table 6. Analysis of variance (ANOVA) results for average permeate flux, COD rejection, and SFD for each factor.
Table 6. Analysis of variance (ANOVA) results for average permeate flux, COD rejection, and SFD for each factor.
ResponsesFactorsDOFSSMSF-Valuep-Value
TMP (bar)26357.683178.84114.010.000
CFV (m/s)254.7627.380.980.394
Average permeate fluxTemperature (°C)2179.1589.583.210.064
MWCO (kDa)24546.302273.1581.520.000
Error/others18501.9027.88
Total2611,639.80
TMP (bar)273.5536.7821.700.000
CFV (m/s)2198.2199.1158.470.00
COD rejectionTemperature (°C)282.5841.2924.360.00
MWCO (kDa)2607.95303.9817.340.00
Error/others1830.511.70
Total26992.79
TMP (bar)246.4823.2496.570.000
CFV (m/s)26.123.0612.720.000
SFDTemperature (°C)23.801.907.890.003
MWCO (kDa)252.9526.47110.000.000
Error/others184.330.24
Total26113.68
Table 7. Individual Taguchi predictions for average permeate flux, COD rejection rate and SFD.
Table 7. Individual Taguchi predictions for average permeate flux, COD rejection rate and SFD.
Response VariablesOptimum Operating ConditionsSignificant FactorsPredicted Optimal Responses
S/N Ratio (dB)Value
J ¯ P A3, B3, C3, D3A, D38.1981.20 L·m−2·h−1
COD RejectionA3, B3, C1, D1A, B, C, D35.2657.92%
SFDA1, B3, C1, D1A, B, C, D−5.101.80
Table 8. Pairwise comparison matrix.
Table 8. Pairwise comparison matrix.
Response μ 1 μ 2 μ 3
μ 1 1.05.02.0
μ 2 1/51.01/4
μ 3 1/24.01.0
Table 9. Utility value based on UF responses ( J ¯ P , COD rejection, SFD) for each repetition of the experiment.
Table 9. Utility value based on UF responses ( J ¯ P , COD rejection, SFD) for each repetition of the experiment.
Trial n°Utility Value
R1R2R3
113.6016.4213.96
215.1915.8115.64
316.6916.4316.99
418.1118.0318.58
518.3716.2417.72
619.4519.8918.53
717.1417.3218.17
818.3919.0219.27
917.4618.1617.32
R1, R2 and R3 is the number of repetitions of the experiment.
Table 10. ANOVA analysis results for multi-response UF (overall utility function).
Table 10. ANOVA analysis results for multi-response UF (overall utility function).
ResponsesFactorsDOFSSMSF-valuep-valueP (%)
TMP (bar)238.8519.4335.500.00060.93
CFV (m/s)25.842.925.340.0159.16
Utility concept Temperature (°C)20.940.470.860.4401.47
MWCO (kDa)28.284.147.570.00412.99
Error/others189.850.55 15.45
Total2663.76 100.00
Table 11. Optimal conditions for multi–response UF predicted using the utility concept.
Table 11. Optimal conditions for multi–response UF predicted using the utility concept.
MethodResponseOptimal Conditions Optimal Values
Multi-response MRSN (dB)A2 B3 C1 D319.82
J ¯ P (L·m−2·h−1)77.22
COD (%)45.69
SFD/ F D ¯ 6.24/30%
Table 12. Permeate quality (process performance) under optimum conditions, at the end of 2 h operating.
Table 12. Permeate quality (process performance) under optimum conditions, at the end of 2 h operating.
ParameterUF Permeate QualityPercent Removal (%)
TSS (g/L)0.000299.57
Turbidity (NTU)0.0897.51
COD (mg/L)81.843.90
Total nitrogen (mg/L)0.5333.75
particle size (nm)12.76–24.36
Table 13. Summary and comparison of experimental and predicted optimal conditions for PMTE.
Table 13. Summary and comparison of experimental and predicted optimal conditions for PMTE.
MethodResponseOptimal Conditions Optimal Values95% CICEConfirmation Experiments% Deviation
Multi-response MRSN (dB)A2 * B3 * C1 D3 *19.7818.54 ≤ MRSN ≤ 21.0819.780.27
J ¯ P (L·m−2·h−1)77.2268.16 ≤ J ¯ P ≤ 86.2881.154.84
COD (%)45.6943.45 ≤ COD ≤ 47.9243.904.07
SFD6.245.39 ≤ SFD ≤ 7.086.013.75
* Significant at the 95% confidence interval.

Share and Cite

MDPI and ACS Style

Sousa, M.R.S.; Lora-García, J.; López-Pérez, M.-F.; Santafé-Moros, A.; Gozálvez-Zafrilla, J.M. Operating Conditions Optimization via the Taguchi Method to Remove Colloidal Substances from Recycled Paper and Cardboard Production Wastewater. Membranes 2020, 10, 170. https://doi.org/10.3390/membranes10080170

AMA Style

Sousa MRS, Lora-García J, López-Pérez M-F, Santafé-Moros A, Gozálvez-Zafrilla JM. Operating Conditions Optimization via the Taguchi Method to Remove Colloidal Substances from Recycled Paper and Cardboard Production Wastewater. Membranes. 2020; 10(8):170. https://doi.org/10.3390/membranes10080170

Chicago/Turabian Style

Sousa, Mayko Rannany S., Jaime Lora-García, María-Fernanda López-Pérez, Asunción Santafé-Moros, and José M. Gozálvez-Zafrilla. 2020. "Operating Conditions Optimization via the Taguchi Method to Remove Colloidal Substances from Recycled Paper and Cardboard Production Wastewater" Membranes 10, no. 8: 170. https://doi.org/10.3390/membranes10080170

APA Style

Sousa, M. R. S., Lora-García, J., López-Pérez, M. -F., Santafé-Moros, A., & Gozálvez-Zafrilla, J. M. (2020). Operating Conditions Optimization via the Taguchi Method to Remove Colloidal Substances from Recycled Paper and Cardboard Production Wastewater. Membranes, 10(8), 170. https://doi.org/10.3390/membranes10080170

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