Determination of Hydrodynamics in Municipal Waste Water by a Lagoon System with Screens

Stabilization lagoons are economic systems that are built for treatment of municipal, industrial and agricultural wastewater; these systems are used in rural communities. Objective of this article is to present the hydrodynamics flow in lagoon system considering 6 screens with 7 channels containing curved forms with slopes suitable to stabilize the flow along each channel, and approach a piston flow. Hydrodynamics of this system with IBER software was analyzed, where was considered the velocity and hydraulic gradient, using Froude number. Also transport of total suspended solids was modelled. Efficiency in the treatment with this design was evaluated, using parameters such as, pH, conductivity, alkalinity, hardness, total solids, dissolved oxygen, redox potential and Chemical Oxygen Demand (COD). Through the results, a homogeneous transport was observed, mainly dissolved oxygen which was concordant with redox potential and COD, also through the curves, short circuits were minimized, avoiding dead zones and making treatment more efficient, finally were possible to comply with regulations of Mexico NOM-001-SEMARNAT-1996 of discharges and the NOM-003-SEMARNAT-1997 for water of agricultural use.


Introduction
Natural treatment systems, such as stabilization lagoon systems and wetlands either natural or artificial, represent a very useful alternative to being environmen-tally friendly systems, as well as being an economic option due to low operating and maintenance costs [1].They are ideal to be implemented in places where there is not a large population or in rural areas, which makes the construction of conventional wastewater treatment plants difficult [2] [3] [4].Stabilization lagoons aim to treat wastewater mainly from domestic, municipal and industrial sources [5] and also meet quantitative parameters, set by regulations in each country that allows their discharge into the receiving environment without causing environmental problems [6].Stabilization lagoons are classified according to nature of biological activity that is carried out: aerobic, anaerobic and facultative [7] [8] [9].
Lagoons systems are usually constructed as the final stage of treatment of previous depuration systems, such as activated sludge, methods that replace chlorination, being the most common disinfection method in these systems [8] [10].
A particular case of stabilization lagoons is maturation lagoons; main functions are: 1) nitrification of ammoniac nitrogen, 2) nutrient removal and 3) clarification of effluent [11] [12].The case study in this journal is a maturation lagoon system with screens which are constructed in Autonomous Metropolitan University-Azcapotzalco in Mexico City, like a pilot system to be used in mainly rural zones, where they design screens to approximate the flow of water to a piston model [13], improving the hydraulic retention times, on the other hand, to achieve an optimum treatment efficiency.
To avoid short circuits and dead zones, curved shapes with slope were designed in every channel change, which also generated more homogeneous flow.In order to achieve this, a study of hydrodynamic was carried out with IBER Software considering two periods of sampling: winter and spring, taking into account hydraulic flow velocities in each channel, and Froude number for ratio of inertia forces and gravity forces acting on water flow.In order to analyze the behavior in curves shapes with slope, a monitoring parameter of total solids in each channel was used.So as to determine efficiency of screens and curves in lagoon system: temperature, pH, conductivity, dissolved oxygen (DO), chemical oxygen demand (COD), biochemical oxygen demand (BOD), total solids suspended (TSS), hardness and alkalinity were monitored.

Materials and Methods
In order to complete municipal wastewater treatment of a horizontal flow subsurface wetland for rural areas [14] a pilot lagoon system in Autonomous Metropolitan University-Azcapotzalco located in Mexico City was built.System measure 12.15 m long, 8 m wide, and 0.70 m high, divided into 7 channels with 6 screens built with light porous block, whose purpose is to fix bacterial film on the walls [15].The main purpose of wetland-maturation lagoon system is to obtain quality water at least for agricultural irrigation.Figure 1 shows a diagram of wastewater treatment system which is fed by a primary treatment plant near the Campus and from wastewater of University.Figure 1  Hydraulic previous consisting of an anaerobic pre-treatment system, Figure 1(a), where wastewater passes to two grids allowing retain sand and/or coarse settle able solids; after this water passes a regulation tank where it is accumulates and flow rate is controlled for system rest, next step is to bring water to two serially septic tanks where anaerobic digestion is performed to remove a portion of organic matter; for finally the wastewater enters to wetland through 12 water shedders in triangular form for influent and a design flow Q of 2 L/s. Figure 1(a) corresponds to the wetland, where it was entered with 1728 ± 0.048 mg/L of TSS, 364.88 ± 14.99 mg/L of BOD and COD 568.35 mg/L ± 18.58 mg/L in winter and 1528 ± 0.055 mg/L of TSS, 560.38 ± 20.57mg/L BOD and 987, 112 ± 51.77 mg/L COD in spring.In Figure 1(c) the lagoon system with screens is presented, which completes the study treatment presented in this article.

Water Sampling in the Lagoon System
Figure 2(a) shows an overview of the lagoon system with screens and curved shapes.Figure 2(b) shows a schematic with 21 sampling sites, where site 1 is located at the entrance in lagoon system and site 21 is located at output in lagoon system.
Lagoon system is built with porous block, as it aims to build such facilities in rural areas for economic reasons and on the other hand it is to take advantage of porosity of the material for fixing bacterial film.
Sampling in each sites selected was monitored for 2 periods: winter and spring 2016, at the same time from 10:00 a.m. to 4:00 p.m. with the purpose of ensure microbial activity.Sampling was carried out at a regular depth of 25 cm with kind bailer samplers at 13:00 hours.

Measures in Situ
In order to determine the behavior of lagoon system as a function of parameters: temperature, pH, conductivity, redox potential (ORP), dissolved oxygen (DO), total solids suspended (TSS) sites in seven mentioned channels were defined.pH was determined according to the NMX-AA-008-SCFI-2000 [16] standard, while temperature was determined according to NMX-AA-007-SCFI-2000 [17].For each parameter a Hanna model HI 1618 sensor was used, with a Hanna HI98150 interface.Finally, for conductivity and ORP measurements a LabQuest Vernier Interface system with its respective sensors was used.

Ex Situ Measures
Alkalinity [18] and Hardness was determined [19].The COD was determined by Hach method with a Hanna brand reactor model HI-839800 and a UV/visible spectrophotometer DR/2400 brand Hach.To obtain total solids (TSS) the technique of standard methods [20] was followed.The volume collected in all cases was one litter in Nalgen containers, kept under refrigeration at 4˚C.In situ measurements, for maximum dissolved oxygen concentration, corresponding to a maximum microbiological activity was performed in all samples at 1:00 p.m. [21].
All determinations were performed on an average of two samplings per week per season (winter and spring).

Design of Landfills for Incoming Water to the Lagoon
It was important to consider optimization and homogenization of water flow in the input of lagoon system, for which landfills were designed, using Thomson Equation (1) for triangular landfills [22]. where: Flow rate that enters in lagoon system from wetland was Q = 2 L/s, and for design of landfill was considered from Thomson Table [22] a coefficient of 1.4, obtaining a total load H equal to 0.867 cm, which was divided among 12 landfills, so that H/landfill of 0.072 cm was obtained (Figure 3).

Modeling Using IBER Software
For hydrodynamic study in lagoon system, IBER software was used [23] which is a numerical model of simulation of free-flow in non-permanent regime and of environmental processes in fluvial hydraulics.Application of IBER covers fluvial hydrodynamics, flood zones evaluation, transport of suspended material and aquatic flows calculation.Hydrodynamic module for two dimensions was used for this study, where conservation equations of mass and momentum were solved.

Results
The variation of parameters: temperature, pH, conductivity, hardness and alka-linity during two sampling periods are presented in Table 1, corresponding to spring and winter.It is observed that values of Temperature did not vary significantly and ranged from 17˚C to 24˚C, adequate values for development of photosynthetic activity [24].In terms of pH, a decrease was observed in both periods.Conductivity determination is related to the presence of ions in water [15] values along the system are presented in Table 1, a decrease in values as you move from one channel to another is observed.However as in COD when entering water flow to curves shapes with slope an increase in its concentration is noticed, due to a reduction of velocity generated by slope, thus allowing ions concentration and therefore conductivity tends to increase.At system output with an increase in velocity, a homogeneous distribution of ions is observed.In addition a dilution and/or adsorption of ions are able to results in site 21 lower concentrations than first channels.As for the Hardness of water in each channel during two sampling seasons, an increase mainly in calcium ions at curves input is observed, behavior that coincides with conductivity variation.Same effect on conductivity variation again presents, water in all sampling sites until system output, resulted with values of 75 -150, reason why it is considered Moderately Hard (MD) according criteria of international hardness [25].Table 1 also shows alkalinity variation along the system, pH values fluctuated in two sampling periods between 7.9 -7.5 [25] region that covers the predominance of bicarbonates, therefore alkalinity belongs to 3 HCO − .Behavior in terms of distribution in different channels is consistent with conductivity and hardness.

Total Solids Suspended (TSS) Variation
In Figure 4, it is observed a higher concentration of solids (mg/L) in every channel change (3 -4, 6 -7, 9 -10, 12 -13, 15 -16 and 18 -19), due to the entrance of water flow in each curve with slope, phenomenon that occurred in both periods.
Overall, a tendency to decrease in concentration of TSS was observed, being scarce on departure at exit of lagoon system, which indicates a homogeneous sedimentation process along each channel.

Dissolved Oxygen (DO) Variation
Dissolved Oxygen (DO) variation in lagoon system is observed in Figure 5.It is Figure 4. Total Solids Suspended (TSS) variation in lagoon system in two seasons.
Figure 5. Dissolved Oxygen variation in lagoon system in two seasons.
important to mention that water which entering in system, comes from a treatment sub-surface wetland of horizontal flow [14] later enters in lagoon through 12 landfills, where flow is distributed and conducted to first channel in lagoon system, a jump of 50 cm height allowing aeration of water is produced.DO oscillate in values 3.04 ± 0.144 mg/L at the entrance of first channel during the winter and 3.00 ± 0.132 in spring.DO concentration increases as water flow along each channel is advanced, later it is reduce to input of each curve with slope, and is indicated by vertical dotted lines.Finally, the concentration improves at the exit of each curve and during straight sections.DO behavior is consistent with the redox potential (ORP) (Figure 6) and COD (Figure 7).Oxygenation was better in winter due to presence of algae and also due to lower evapotranspiration.In spring, oxygenation was lower despite presence of algae and an increase in temperature was observed.Due to action of wind and rain season during April and May respectively, there was a significant fall of leaves, reason for which degradation of organic matter could have consumed part of oxygen.Overall a decrease of DO is observed, it was not so abrupt in first channels, and in spite of this the behavior profiles were retained.

Redox Potential Variation (ORP)
In Figure 6, it was observed that ORP approached increasingly positive values as the flow advanced from channel to channel, behavior that coincided with DO variation.In winter, the ORP presented values more negative than in the spring, but from site 5 (corresponding to channel 2) acquired values slightly more positive than those of spring, behavior that lasted until the exit of the system.
It is also observed that at site 1 (entrance to the system) is the site with the most negative value in both stations, as the flow advances along the channels, an improvement in conductivity is observed, behavior that coincides with the Improvement of DO and COD.It can also be observed that at the beginning of each ORP curve it acquired negative values.

Biochemical Oxygen Demand (BOD)
Figure 8 shows results using the two graphs according two stations, a very marked trend of decreasing in BOD was observed.At entrances of curves shapes with Figure 6.ORP variation in lagoon system in two seasons.slope, BOD increased, due to the fact that there was more accumulation of biodegradable material, the degradation was increasing until its decrease in the output of lagoon system.

Analysis of Chemical Oxygen Demand (COD)
COD is a representative measure of organic contamination in an effluent, parameter which will be controlled within the different discharge regulations.Figure 7 presents results of COD as a parameter for monitoring the degradation of organic material, both biodegradable and non-biodegradable in the different channels along the lagoon system.There was a general trend of decreasing COD along the channels at both stations.At the beginning of each banked curve a higher concentration of COD was presented, where more oxygen was required for oxidation of organic matter and it is related with ORP values because these parameters are closely related.

Modelling Results with Software IBER
In Figure 9(a), the result of modeling with IBER software is presented, observing a velocity distribution in different sections of lagoon system, identifying the variation of water flow with different colors: light blue means a uniform flow and changes to dark blue in external curve section, where the flow decreases its velocity.In internal curve section there is an increase that impels it to exit of the curve, stabilizing this flow along of channel.The phenomenon is repeated in each curve until last channel, where flow acquires a greater velocity until output of lagoon, which is represented with yellow-green color (Figure 9

Transport of Total Suspended Solids with IBER
This software was also used to model TSS transport in the lagoon system (Figure 10).The concentration of TSS in water is corresponding to field measured values which were introduced in the software.The results indicate that the highest concentration occurred in the first channel (red color), although it seems that before the entrances to the curves the TSS tend to concentrate, this can be due to reduction velocity.In the other channels also there is a drag of TSS, always concentrating at the entrance of each curve, until concentration decreases.It is important to note with this model that in zones of the curves, blue color is homogeneous, which indicates that there are no short circuits, therefore there are no "dead zones" and the hydraulic flow improves significantly.up in each channel.Regarding BOD in the first and seventh channels, the best efficiency was observed in the sixth and seventh channels for COD.

System Efficiency
Finally in Table 3, efficiency of all system is shown, considering Figure 1.In pre-treatment, practically the same efficiency was achieved in two seasons with a little higher value in winter for the three parameters.In the wetland, TSS concentration was slightly higher in winter, however BOD and COD were better in spring, possibly the rise of temperature improved the ability of microorganisms to degrade MO.
In case of lagoon system, the treatment for TSS was better in winter, because fewer solids were introduced in wetland.The BOD was slightly better in winter, where the flow velocity can be influenced by the curved forms of each channel that avoid short circuits and therefore dead zones, same can be said of COD.

Conclusions
Design of 7 screens in lagoon helped to approximate the total water flow to a piston flow, which was optimized by the slopes in curves that induced to cross to the flow in sections where tangential forces were distributed evenly, avoiding thus presence of dead zones, short circuits and dispersed flows.
In addition, a good distribution of the hydraulic flow was achieved, which was demonstrated through behavior of DO, OPR, TSS, BOD and COD, as shown in graphs.With IBER software, it was possible to obtain velocity distribution and Froude number in different sections, also modelling the transport of TSS, which allowed knowing areas susceptible to TSS saturation.The slope change in curve reduces the bottom drag of sediments, producing a slight decrease in flow velocity, causing the total solids to be retained in pre-curve sections.
When analyzing removal rates of both the lagoon and entire treatment system, it was observed that lagoon system allowed completing the treatment in compliance with Mexican norm NOM-001-SEMARNAT-1996 of discharges and NOM-003-SEMARNAT-1997 for agricultural water.

Figure 1 .
Figure 1.Diagram of natural treatment system (a) Pre-treatment system, (b) Wetland system, (c) Maturation lagoon system.

Figure 2 .
Figure 2. (a) Lagoon system, (b) Location of sampling sites in maturation lagoon.

Figure 3 .
Figure 3. Landfills for entry of water from the wetland.
: h = depth; U x , U y = Horizontal velocities averaged in depth; G = Acceleration of gravity; Z s = Free surface elevation; τ s = Friction in free surface due to the friction produced by wind; ρ = Density of water; Ω = Angular velocity of rotation of the earth; Λ = Latitude of the point considered; e xx τ , e xy τ , e yy τ = Horizontal effective tangential tensions; M s , M x , M y = Source/sink sum of mass and moment, through which the modelling is performed.

Figure 7 .
Figure 7. COD variation in lagoon system in two seasons.
(a)).Note the almost yellow color where flow is virtually free.Froude number was determinate, which relates forces of inertia and gravity, implicitly considering velocity and hydraulic gradient at each point of path the water flow, evaluating type of regime.
Figure 9(b) shows increasing Froude number as the water flow increased in each channels, ending with a tendency toward a critical regime in the seventh channel.It can also be observed, as in velocity case

Figure 8 .
Figure 8. BOD variation in lagoon system in two seasons.

Table 1 .
Physicochemical parameters variation in lagoon system with screens.

Table 2 summarizes
BOD and COD values as parameters for efficiency follow

Table 2 .
Efficiency of treatment by the lagoon system with screens and curve shapes with slope.

Table 3 .
Efficiency at each treatment stage for season and efficiency in all systems.