One advantage of the flotation system is its high applicability as a compact facility. Recently, research for the introduction of additional high-rate dissolved air flotation (DAF) into conventional water treatment plants has been progressing in areas of remodeling water treatment plants, advanced treatment, and combined sewer overflow (CSO) treatment. For successful design and operation of high-rate DAF, laboratory experiments to verify model results and estimate standard values of major parameters by model simulation were performed in this study. With the model simulation having a population balance concept, an operational analysis of bubble-floc aggregates considering a hydraulic situation in the separation zone during high-rate DAF was conducted. Based on the model simulation results, a basic level for high-rate DAF design at a hydraulic loading rate of 30 m/h was determined. On the basis of population balance theory, it can describe the process of bubble-floc aggregate formation at a contact zone and calculate the rate of transfer into the separation zone.

## INTRODUCTION

The key technology of high-rate dissolved air flotation (DAF) is achieved by rapid separation of bubble-floc aggregates. According to the specific gravity of the floc in the separation zone, bubble-floc aggregates have a slower rise velocity than microbubbles. Considering the rise velocity of the aggregates, the downward flow rate in the separation zone of the flotation tank dominates the overall efficiency of the DAF process (Edzwald *et al.* 1999). When designing the high rate DAF, much attention is required because the hydraulic loading rate is faster than the rise velocity of the bubbles. As such, it is highly possible that the treated water is discharged before bubbles and flocs are attached or before the aggregates rise to the water surface. The ideal direction pattern that is maintained in the internal portion of the tank has been suggested as the point that prevents this phenomenon.

The flotation tank consists of a contact zone and separation zone, and the excellence of DAF performance is apparent when the three features of the ideal flow pattern in the separation zone are formed. The three features of the ideal flow pattern include the stratified flow of the surface, the vertical plug flow below the stratified flow, and the horizontal plug flow at the bottom floor. However, this ideal theory for the conventional DAF process with hydraulic load rate (5–15 m/h) does not work in the high-rate DAF process (30 m/h) because the hydraulic load (hydraulic loadings) rate exceeds the rise velocity (Edzwald 2010). The high-speed system can remove free bubbles and the bubble-floc aggregates in the stratified flow area of the separation zone. In particular, research on the hydraulic flow characteristics and stratified flow that occur in the separation zone is required in order to merge the flow pattern in the model. The flow pattern through the separation zone does not follow the simple ideal case of plug flow in the vertical direction. It is influenced by the velocity of water above the baffle (cross-flow velocity), the hydraulic loading, the aspect ratio (length to width ratio), how water is withdrawn at the outlet, and the air bubble suspension in the upper part of the separation zone (Lundh *et al.* 2002).

The basic research carried out over the last 30 years has laid the groundwork for the development of high-rate DAF technology (Wang *et al.* 2005). According to previous studies (Haarhoff & Edzwald 2001; Rubio *et al.* 2007; Oliveira & Rubio 2012), the major trend is to decrease flocculation times, increase DAF loadings, and develop compact (small ‘footprint’ areas) and efficient treatment units. Small aggregates can be easily removed at high DAF (modern design) rates, contradicting the conventional bias that large floc units and bubbles are required for successful separation. In other words, the effective separation of the combined bubble-floc aggregates with small bubbles is determined by the flow rate in the separation zone of the DAF process. While larger separation zones lead to higher separation efficiency of the bubble-floc aggregates, they also create an inefficient situation by increasing the process size. Therefore, technology for the maximum separation of the bubble-floc aggregates through high separation speed in the separation zone while maintaining a minimal physical size is important. Lakghomi *et al.* (2012) reported that increasing the amount of air improved the removal efficiency and generated a useful stratified horizontal flow pattern.

Because many previous studies report a threshold for high-rate DAF with a hydraulic load rate over 0.4167 m/min (or 25 m/h), modeling and experimental studies are required for water treatment through various methods of flotation separation and measuring the rise velocity of the bubble-floc aggregates for the proper operation and design of high-rate DAF. Thus, the purpose of this study is to understand the characteristics of operating and designing parameters for the DAF process and to evaluate the scope and applicability of the proposed threshold for high-rate DAF with a hydraulic load rate of 0.5 m/min (or 30 m/h) for treating combined sewer overflows (CSOs). In order to successfully design and operate high-rate DAF, standard values of major parameters were determined through model simulation, and laboratory tests were performed to verify the estimated model values in this study. Considering the hydraulic situation in the separation zone and the operational analysis of the bubble-floc aggregates, a basic level of high-rate DAF with a hydraulic loading rate of 30 m/h for design has been determined by simulating the model introducing the population balance concept. The model simulation can describe the formation process of bubble-floc aggregates in the contact zone and calculate the transfer rate of the formed aggregates that transfer into the separation zone. Focusing on the change of the bubble volume fraction in the separation zone, depending upon the hydraulic loading rate, the model simulation was conducted using the population balance turbulence (PBT) model. Laboratory tests were conducted to measure the rise velocity of the bubbles and the bubble-floc aggregates as a function of the run-time change and the depth change of the milky water layer.

## MATERIALS AND METHODS

### Theoretical approach

*i*-bubbles

*n*at mixing time

_{i}*t*). On the basis of the assumptions reported by Fukushi

*et al.*(1995) and modified by Leppinen

*et al.*(2001), the number of flocs attached,

*i*-bubbles

*n*, is represented by ordinary differential equations as follows: where

_{i}*k*is the turbulent collision rate constant,

*n*is the number of bubbles per unit volume, and

_{bubbles}*α*is the adhesion efficiency.

_{i}*V*,

_{b}*V*, and

_{f}*V*). Considering the bubble rise velocity,

_{fb}*V*, for a water column under no-loading condition with stagnating water flow, Stokes' law describes the rise velocity of bubbles,

_{b}*V*, for Reynolds number,

_{b}*Re*≤ 1, of bubble diameter

*d*≤ 125 μm at 20 °C according to Equation (3): where

_{b}*V*is the terminal rise velocity of the bubbles, cm/s;

_{b}*g*is the gravitational constant, 980 cm/s

^{2};

*d*is the effective diameter of the bubbles, cm;

_{b}*ρ*is the density of the bubbles, g/cm

_{b}^{3};

*ρ*is the density of the aqueous phase, g/cm

_{w}^{3}; and

*μ*is the viscosity of the aqueous phase,

_{w}*cp*.

*et al.*2011). Therefore, it is necessary to compare the bubble layer formed by the 120 μm bubbles with the results in the time domain forming the stratified flow according to the bubble layer. Haarhoff & Edzwald (2004) proposed the theoretical rise velocity of the bubble-floc aggregates for

*Re*≤ 1 as shown in Equation (4). If

*Re*is in the transition region, 1 <

*Re*< 50, Equation (5) can be used to calculate the rise velocity. For solving Equation (5), it is required to first calculate Equation (6) for the diameter of the aggregates, and Equation (7) for the density of the aggregates, to be dependent on the attached bubble number: where

*V*is the terminal rise velocity of the agglomerate, cm/s;

_{fb}*ρ*is the density of the agglomerate, g/cm

_{fb}^{3};

*ρ*is the density of the aqueous phase, g/cm

_{w}^{3};

*μ*is the viscosity of the aqueous phase,

_{w}*cp*; and

*N*is the number of bubbles per unit volume.

_{b}*K*in Equation (8) is described in the form of agglomerates, and the effects of

*K*will be described for the resistance increase in drag. When large bubbles that are greater than 100 μm are attached on a small floc of 40 μm, the aggregates are substantially spherical and

*K*is 24. If the floc is significantly larger than 100 μm, the aggregates approach the form of floc and

*K*is 45. Since

*K*is 24 of the floc in 40 μm, and

*K*is 45 of the floc in 170 μm, it is assumed that

*K*becomes gradually smaller:

*et al.*(1986) formulized Equation (9) to calculate the maximum number of bubbles attached on the floc:

Many aggregates formed in the contact zone rise to the top layer. However, a portion of them is moved into the separation zone and rises into the top layer to be separated if the rise velocity in the laminar state is greater than the hydraulic loading rate. Design and performance of the footprint area for removal of free bubbles and of the bubble-floc aggregates applying the Hazen theory of the sedimentation tank have been considered the idealized DAF tank. The vertical plug flow is assumed to occur over the main part of the separation zone (i.e., the clarification section). The horizontal flow below the bubble-floc aggregate layer and above the baffle for separating the contact zone and the separation zone is ignored when evaluating the clarification performance. The free bubbles and the bubble-floc aggregates are removed if the rise velocities, *V _{b}* and

*V*, exceed the hydraulic loading rate (

_{fb}*V*) in the separation zone. As mentioned above, the standard value of rise velocity of the bubbles and the bubble-floc aggregates is approximately 20 m/h. Therefore, these bubbles and the aggregates are removed when the hydraulic loading rate is less than 20 m/h in the separation zone. Edzwald (2010) summarizes that the hydraulic loading rate is conservative enough in order to evaluate the performance of the conventional DAF system. The typical hydraulic loading rate of the conventional DAF system is 5–15 m/h. When considering the footprint area (

_{sz−hl}*A*) of only the separation zone and having a recycle rate of 10%, the hydraulic loading rate in the separation zone is still 20 m/h or less. Thus, even though the plug flow does not occur, the hydraulic loading rate is used to design the separation zone of the conventional DAF system.

_{sz}*V*acting on the separation zone is as follows: If considering the recycle ratio: If considering the dead space m, which exists in the real separation zone: The rise velocity of the aggregates reaches the maximum when the occurrence of such dead space is suppressed (Behin & Bahrami 2012). The separation area can be calculated according to changing the recycle flow rate as functions of the varying recycle ratio

_{sz−hl}*R*and the fixed hydraulic loading rate of 30 m/h (0.5 m/min). In addition, the rise velocity

*V*of the aggregates may be used to calculate the hydraulic loading rate

_{fb}*V*and the dead space

_{sz−hl}*m*.

For theoretical experiments, the PBT model written by MATLAB software was used for simulation of size distribution, rising velocity, time, and removal efficiency of bubbles, flocs, and floc-bubble aggregates. The simulation was performed using the following information: floc sizes between 0.1 and 210 μm, bubble size of 100 μm, initial density of flocs 1,100 kg/m^{3} and 500 kg/m^{3} (Santana *et al.* 2012), and hydraulic loading rate *V _{sz−hl}* of 0.5 m/min (or 30 m/h) in order to examine the change in the rise velocity of one aggregate formed by one bubble attached on one floc; the initial number concentration of flocs,

*N*, 2,000, the number of bubbles

_{p0}*N*from 10

_{bi}^{5}bubbles to 10

^{6}for estimating the removal efficiency of bubble-floc aggregates; inflow rate of 1.62 m

^{3}/h and recycle rate of 10%, bubble diameter

*d*of 120 μm, floc diameter

_{b}*d*of 250 μm, and an initial density of 900 kg/m

_{f}^{3}in order to investigate critical values of parameters.

### Labs Test

The Labs Test for rise velocity verification of bubble and bubble-floc aggregates was performed in order to measure the rise velocity of the bubbles and the bubble-floc aggregates. In the experiments, two types of flocs were used: the coagulation floc and the flocculation floc. The coagulation floc (Floc Type 1) was formed by rapid mixing at 60 revolutions per minute (rpm) for 30 s, while the flocculation floc (Floc Type 2) was prepared by slow mixing at 30 rpm for 2 min in order to make the floc larger than after the rapid mixing.

## RESULTS AND DISCUSSION

### Rise velocity of bubble and bubble-floc aggregates

According to the hydraulic loading rate increase, the mean velocity of the separation tank is increased and the microbubbles may move to the deep portion of the separation tank. These microbubbles determine the rise velocity as a function of the bubble volume. While many aggregates are separated by rising to the top layer, a portion of them are separated by rising to the top layer when their rise velocities are greater than the hydraulic loading rate in the laminar flow state after moving to the separation zone. If the microbubbles move into the separation zone and do not increase in size, the rise rate is low. However, the bubbles in the separation zone are bigger than those in the contact zone (Haarhoff & Edzwald 2004). The bubbles that are released from the nozzle are distributed at a size of 40 ∼ 60 μm in the contact zone, but they are distributed at a size of 50–150 μm in the separation zone through growth, coalescence, and clustering phenomena (Leppinen & Dalziel 2004).

*Re*) and the rise velocity of the microbubbles. The rise velocity of the bubbles smaller than 150 μm is 0.83 cm/s (0.5 m/min or 30 m/h) for

*Re*< 1.

^{3}(Santana

*et al.*2012). The simulation results are shown in Figure 4 compared with the hydraulic loading rate

*V*of 0.5 m/min (or 30 m/h). Figure 4 indicates the rise velocity of the aggregates for one bubble of 100 μm attached on the various floc sizes with an initial density of 1,100 and 500 kg/m

_{sz−hl}^{3}. For flocs that are 10 μm or less, the rise velocity of the aggregates is approximately 0.33 m/min (or 20 m/h), and it is the same with the rise velocity for the same size bubbles only. In the case of an initial floc density of 1,100 kg/m

^{3}and a floc of 200 μm, the rise velocity of the aggregates approached zero. In the case of an initial floc density of 500 kg/m

^{3}, the rise velocity of the aggregates increased when one bubble was attached on the floc size of 100 μm or more, and it was over 0.5 m/min (or 30 m/h) when the floc size was 210 μm.

### Rise velocity of bubble-floc aggregates in terms of bubble size and number

While the rise velocity of the large floc achieved a high rate, the attachment of many bubbles was required. The rise velocity *V _{fbi}* was defined as the proper rise velocity of the proper bubble-floc aggregates when the rise velocity

*V*became greater than the hydraulic loading rate

_{fbi}*V*of 0.5 m/min (or 30 m/h) and the rise time

_{sz−hl}*t*was faster than the retention time

_{r}*τ*. In order to check the effect of the number of bubbles attached on the floc, the rise velocity of the aggregates was calculated according to the bubble number attached on a 210 μm floc having an initial floc density of 500 kg/m

_{s}^{3}. The calculated rise velocity is summarized in Table 1. When one bubble is attached on one floc,

*N*is defined as the number of bubbles,

_{b1}*V*is indicated as the rise velocity when

_{fb1}*V*≥

_{fb1}*V*and

_{sz−hl},*t*≤

_{r}*τ*for meeting flotation condition in the separation zone. Therefore, in the case of a bubble diameter of 100 μm, even though the rise velocity is faster than the hydraulic loading rate, two bubbles are needed to produce a rise time faster than the retention time, as shown in Table 2. When the rise velocity

_{s}*V*is faster than the hydraulic loading rate

_{fbi}*V*and the rise time

_{sz−hl}*t*is faster than the retention time

_{r}*τ*,

_{s}*N*is defined as the proper number of bubbles in order to form the bubble-floc aggregates.

_{bi}Bubble diameter (μm) | 10 | 20 | 30 | 40 | 50 | 60 | 70 | 80 | 90 | 100 |
---|---|---|---|---|---|---|---|---|---|---|

N _{bmax} | 1,385 | 346 | 154 | 87 | 55 | 38 | 28 | 22 | 17 | 14 |

N _{b1} | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 |

V m/min _{fb1} | 0.38 | 0.38 | 0.39 | 0.39 | 0.40 | 0.42 | 0.43 | 0.45 | 0.48 | 0.51 |

N _{bi} | 398 | 83 | 32 | 16 | 10 | 6 | 4 | 3 | 3 | 2 |

V m/min _{fbi} | 0.59 | 0.59 | 0.59 | 0.59 | 0.61 | 0.59 | 0.59 | 0.60 | 0.68 | 0.64 |

V m/min _{sz−hl} | 0.5 | 0.5 | 0.5 | 0.5 | 0.5 | 0.5 | 0.5 | 0.5 | 0.5 | 0.5 |

Bubble diameter (μm) | 10 | 20 | 30 | 40 | 50 | 60 | 70 | 80 | 90 | 100 |
---|---|---|---|---|---|---|---|---|---|---|

N _{bmax} | 1,385 | 346 | 154 | 87 | 55 | 38 | 28 | 22 | 17 | 14 |

N _{b1} | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 |

V m/min _{fb1} | 0.38 | 0.38 | 0.39 | 0.39 | 0.40 | 0.42 | 0.43 | 0.45 | 0.48 | 0.51 |

N _{bi} | 398 | 83 | 32 | 16 | 10 | 6 | 4 | 3 | 3 | 2 |

V m/min _{fbi} | 0.59 | 0.59 | 0.59 | 0.59 | 0.61 | 0.59 | 0.59 | 0.60 | 0.68 | 0.64 |

V m/min _{sz−hl} | 0.5 | 0.5 | 0.5 | 0.5 | 0.5 | 0.5 | 0.5 | 0.5 | 0.5 | 0.5 |

Bubble diameter (μm) | Attached bubbles (number) | Rise velocity, V (m/min) _{fb} | Reynolds number, Re (number) | Dead space, m (number) | Rise time, t (min) _{r} |
---|---|---|---|---|---|

10 | 398 | 0.588 | 2.08 | 0.0003 | 0.397276 |

20 | 83 | 0.589 | 2.10 | 0.0020 | 0.396595 |

30 | 32 | 0.590 | 2.12 | 0.0022 | 0.396502 |

40 | 16 | 0.589 | 2.13 | 0.0020 | 0.396594 |

50 | 10 | 0.606 | 2.20 | 0.0289 | 0.385897 |

60 | 6 | 0.593 | 2.16 | 0.0083 | 0.394101 |

70 | 4 | 0.589 | 2.15 | 0.0016 | 0.396761 |

80 | 3 | 0.600 | 2.20 | 0.0197 | 0.389543 |

90 | 3 | 0.679 | 2.54 | 0.1335 | 0.34432 |

100 | 2 | 0.641 | 2.38 | 0.0819 | 0.364826 |

Bubble diameter (μm) | Attached bubbles (number) | Rise velocity, V (m/min) _{fb} | Reynolds number, Re (number) | Dead space, m (number) | Rise time, t (min) _{r} |
---|---|---|---|---|---|

10 | 398 | 0.588 | 2.08 | 0.0003 | 0.397276 |

20 | 83 | 0.589 | 2.10 | 0.0020 | 0.396595 |

30 | 32 | 0.590 | 2.12 | 0.0022 | 0.396502 |

40 | 16 | 0.589 | 2.13 | 0.0020 | 0.396594 |

50 | 10 | 0.606 | 2.20 | 0.0289 | 0.385897 |

60 | 6 | 0.593 | 2.16 | 0.0083 | 0.394101 |

70 | 4 | 0.589 | 2.15 | 0.0016 | 0.396761 |

80 | 3 | 0.600 | 2.20 | 0.0197 | 0.389543 |

90 | 3 | 0.679 | 2.54 | 0.1335 | 0.34432 |

100 | 2 | 0.641 | 2.38 | 0.0819 | 0.364826 |

The maximum number of bubbles attached on the 210 μm floc is 1,385 for the 10 μm bubbles and 14 for 100 μm bubbles, that is, the smaller bubbles lead to more bubble attachments (Han 2002). In order to examine the rise velocity of the aggregates according to the attached number of bubbles, the rise velocity when only one bubble is attached needs to be calculated. When the calculated rise velocity was compared with the hydraulic loading rate, the rise velocity of the aggregates attached by a bubble of only 100 μm was greater than the hydraulic loading rate of 0.5 m/min (30 m/h). However, it did not meet the flotation condition in the separation zone. The flotation condition in the separation zone is met when the rise velocity of the aggregates is greater than the hydraulic loading rate (i.e., *V _{fb}* ≥

*V*) and the rise time

_{sz−hl}*t*is faster than the retention time

_{r}*τ*To reach

_{s}.*t*≤

_{r}*τ*, the attached number of bubbles,

_{s}*N*, was 398 bubbles for the bubble size 10 μm and two bubbles for the bubble size 100 μm. Ranging between 0.59 and 0.68 m/min, the calculated rise velocity of the aggregates was faster than the hydraulic loading rate of 0.5 m/min (or 30 m/h), but the Reynolds number

_{bi}*Re*was in the transition region between 1 and 50.

As shown in Table 1, the rise velocity changed according to the size of bubbles. As an increase in the hydraulic loading rate *V _{sz−hl}* of 0.5 m/min (or 30 m/h) indicates an increase of inflow rate, the recycling flow rate was increased if the recycle ratio was fixed. If the flow rate of the circulating water is not increased, the density of the bubble layer in the separating tank is reduced, and the quality of the treated water worsens. If the flow rate in the contact zone is too rapid, the possibility of collision between bubble and floc becomes smaller and the contact zone will not be sufficiently used.

### Removal efficiency of bubble-floc aggregates in terms of attached bubble number on a floc

*X*when the bubbles injected into the contact zone collide and attach on a floc during the contact time. The floc is flocculated particles which form floc-bubble aggregates when they meet with bubbles in the contact zone during the contact time. The initial number of the concentration of flocs,

*N*, is 2,000 in this study. By colliding and attaching during the contact time, most of the formed aggregates were removed by floating to the top bubble layer as fast as the removal efficiency

_{p0}*X*. Free bubbles and flocs that did not collide/attach and the remainder of the formed aggregates with a rise velocity slower than the hydraulic loading rate would be removed by increasing the rise velocity and causing the attachment of free bubbles in the separation zone or their release according to the water flow. By increasing the number of bubbles from 10

^{5}bubbles to 10

^{6}, the removal efficiency

*X*approached near 100% within the first few seconds, as shown in Figure 5. If the removal efficiency

*X*increases, it is considered to increase the separation efficiency of the untreated flocs and aggregates in the separation zone of the high-rate DAF.

The rise velocity of the aggregates as a function of the recycle ratio was examined in order to check the degree of dead space to be considered when designing the separation tank. Considering several parameters, including a hydraulic loading rate of 30 m/h, an inflow rate of 1.62 m^{3}/h, and a recycle rate of 10%, the results of the retention time, the rise time, the rise velocity of the aggregates, the Reynolds number, and the dead space are summarized in Table 2.

*d*of 120 μm and the floc diameter

_{b}*d*of 250 μm, and having an initial density of 900 kg/m

_{f}^{3}, Figure 6 presents the results of several parameters including aggregate diameter (

*D*), critical floc diameter (

_{fb}*Critical D*), and critical aggregate diameter (

_{f}*Critical D*). The aggregate diameter

_{fb}*D*increased as the number of bubbles attached increased; it increased from 125 μm when one bubble attached to 160 μm when ten bubbles attached, and the critical aggregate diameter,

_{fb}*Critical D*, decreased from 60 to 20 μm by increasing the number of bubbles attached. If the critical aggregate diameter is smaller than the diameters of the microbubbles, flocs, and aggregates, the fate of the aggregates in the separation zone was influenced by the difference between the drag force due to the surrounding stream flow and the buoyancy due to the bubble volume. The smaller microbubbles, especially those with diameters less than 10 μm, were discharged easily with low flow velocity.

_{fb}### Labs Test for rise velocity verification of bubble and bubble-floc aggregates

^{3}.

## CONCLUSIONS

For flotation in the separation zone during laminar flow, problems arise with a hydraulic loading rate over 0.4167 m/min (25 m/h) because it causes turbulent flow. Thus, a change in the design of the flotation tank is required. When the hydraulic loading rate in the separation zone is 0.4167 m/min to 0.667 m/min (25 to 40 m/h), obtaining sufficient flotation space or maintaining a thicker microbubble cloud layer is required in the separation zone; that is, a thicker bubble layer leads to a higher removal efficiency. The bubble layer can be made much thicker to raise capacity when using the high-rate DAF. In this study, a simulation using a population balance in turbulence (PBT) model concept was conducted and simulation results were verified through laboratory experiments. The change in bubble volume fraction in the high-rate DAF tank using a hydraulic loading rate of 30 m/h was examined through simulation and laboratory experiments.

For the high-rate DAF processes such as CSO treatment, using a hydraulic load rate of 30 m/h to be feasible, these conditions must be satisfied: (1) the flow pattern in the separation zone and stratified flow must be maintained and (2) when dispersion of flocs is very rare, the probability of contact between flocs and bubbles is increased, as many bubbles are supplied. Therefore, the rise velocity of the aggregates is increased when these conditions are satisfied. Especially, the following two things must be considered carefully: (1) if the bubble diameter is less, the collision efficiency between the floc and bubble is increased but the rise velocity of the bubble-floc aggregate becomes slow; and (2) along the flow of the effluent, bubble-floc aggregates smaller than the critical bubble-floc aggregate are discharged with flocs that are not combined with bubbles.

In this study, the effect of the bubble volume fraction in the separation zone of the high-rate DAF on hydraulic loading rate was confirmed. In addition, the hydraulic loading rate may influence the change of flow. The change of flow in the high-rate DAF tank as a function of the hydraulic loading rate may change the tendency of rise along the sloped baffle after inflow water is mixed with recycled water. Computational fluid dynamics (CFD) performance is required to estimate a more accurate change of flow. Therefore, future studies may focus on CFD simulation using information obtained from the results of the current study.

## ACKNOWLEDGEMENTS

This study is supported by the Korea Ministry of Environment (MOE) as ‘Eco-Innovation (EI) Program’ (E314-00015-0412-1) and the National Research Foundation of Korea (NRF) with grants from the Ministry of Education in Korea, 2014 (NRF-2012R1A1A4A01010342).