Performance of Closed Loop Venturi Aspirated Aeration System: Experimental Study and Numerical Analysis with Discrete Bubble Model

In wastewater treatment plants, aeration plays a significant role as it encourages aerobic respiration of microbes, which are necessary to break down carbonaceous matter in the waste stream. This process can account for the majority of energy use in wastewater treatment plants. The aeration process is also necessary in odor control in lagoons and in the aquaculture industry. Generally, the aeration process is accomplished with compressors or blowers which can be of low efficiency due to ideal gas laws. This study introduces the idea of increasing aeration efficiency by looping water from a reservoir through a piping network which includes a venturi aspirator at its inlet. For this purpose, an experimental study has been conducted in a laboratory setup with a pump which pulls water from a tank, passes it through a Venturi aspirator and a helical piping network intended to increase bubble residence time, before depositing it back into the bulk fluid tank. This same process is modeled computationally, using a discrete bubble method (DBM), with good agreement with experiments. The overall purpose here is to determine the optimal configurations for standard aeration efficiency (SAE) and the standard oxygen transfer rate (SOTR). A parametric study has been implemented using the DBM based on different hydraulic and flow parameters. The model is also used to predict the SAE of a hypothetical aeration system. Results indicate that it is possible to achieve SAE values in the range of surface aerators or submersed jet aerators using the proposed aeration system with less complex components, thereby decreasing overall costs.


Introduction
The ecological quality of a water system depends on the availability of the dissolved oxygen (DO) in the system. Although water molecules contain an oxygen atom, this is not readily available to cellular life for respiration. In fact, it is a measure of the amount of diatomic oxygen (O 2 ) freely available in the water. Therefore, the dissolved oxygen level in the water is an important indicator of water quality as aquatic organisms and fish live on the dissolved oxygen in the water. If the dissolved oxygen level is too low, the aquatic life may not survive and the system may have odor problems as well. Additionally, excessive dissolved nitrogen can be harmful to fish and zooplankton. Naturally, oxygen enters the stream from the atmosphere where groundwater discharges into the stream. Therefore, rapidly moving water contains more dissolved oxygen than stagnant water. Bacteria in water depending on this dissolved oxygen consume organic matter. As the DO level decreases due to temperature variation throughout the day or in different seasonal periods, the organic matter can grow to excessive levels and create undesirable bacterial growth. Entraining air bubbles can help oxygen enter the water. As a result, many physical processes have been used in industrial and environmental water streams to aerate liquid by the entrainment of air bubbles.
The goal of this study is to explore the effect of different parameters in venturi aeration systems to find an optimized system that will have higher standard aeration efficiency (SAE). An analytical model is developed based on the discrete bubble model premise. One of the objectives of this study is to investigate whether this analytical model can be applied to model the oxygen transfer for an accompanying experiment. With model validation, the DBM technique is then used to determine the effects of pipe diameter, pipe length and number of aerators on resulting SAE. This paper is completely unique in that it is the first study conducted to consider the possibility of an accompanying "closed loop tube" to increase bubble residence time and hence, aeration efficiency, when compared to other venturi aeration systems.

Clean Water Tests
To compare different aeration systems and configurations, the clean water testing procedure described in ASCE 2007 [22] is followed. This reference describes the procedure to calculate the standard oxygen transfer efficiency (SOTE, %), the standard oxygen transfer rate (SOTR, kgO 2 /h), and the standard aeration efficiency (SAE, kgO 2 /kWh). To calculate the SAE, the power measurement is important and can be expressed in different ways. Generally, power measured directly from the electrical cabinet is usually desirable as it includes blower, coupling, gearbox, and motor inefficiencies. To measure SAE, the standard conditions are 20 • C temperature, zero DO, mean atmospheric pressure, and zero impact of water salinity or other contaminants (e.g., α factor = 1.0, β factor = 1.0). Clean water test results can be used as a warranty of performance between different systems.

Experimental Procedure
Before starting each experiment (see Figure 1), a large aeration tank was filled with water. The aeration tank maximum capacity is 300 gallons. The height of the tank is 6.75 ft with 2.92 ft in diameter. The tank is filled with 250 gallons of tap water and the height of the water level from the bottom of the tank is approximately 5 ft. In industrial boiler systems, it is very common to use anhydrous sodium sulfite (Na 2 SO 3 ) in the presence of a catalyst (CoCl 2 ) to scavenge DO in water [23]. Sodium Sulfite (Na 2 SO 3 ) was added to the water which produced Sodium Sulfate Na 2 SO 4 . The following Equation shows the chemical reaction achieving oxygen scavenging: Water 2020, 12, x FOR PEER REVIEW 3 of 23 model is developed based on the discrete bubble model premise. One of the objectives of this study is to investigate whether this analytical model can be applied to model the oxygen transfer for an accompanying experiment. With model validation, the DBM technique is then used to determine the effects of pipe diameter, pipe length and number of aerators on resulting SAE. This paper is completely unique in that it is the first study conducted to consider the possibility of an accompanying "closed loop tube" to increase bubble residence time and hence, aeration efficiency, when compared to other venturi aeration systems.

Clean Water Tests
To compare different aeration systems and configurations, the clean water testing procedure described in ASCE 2007 [22] is followed. This reference describes the procedure to calculate the standard oxygen transfer efficiency (SOTE, %), the standard oxygen transfer rate (SOTR, kgO2/h), and the standard aeration efficiency (SAE, kgO2/kWh). To calculate the SAE, the power measurement is important and can be expressed in different ways. Generally, power measured directly from the electrical cabinet is usually desirable as it includes blower, coupling, gearbox, and motor inefficiencies. To measure SAE, the standard conditions are 20 °C temperature, zero DO, mean atmospheric pressure, and zero impact of water salinity or other contaminants (e.g., α factor = 1.0, β factor = 1.0). Clean water test results can be used as a warranty of performance between different systems.

Experimental Procedure
Before starting each experiment (see Figure 1), a large aeration tank was filled with water. The aeration tank maximum capacity is 300 gallons. The height of the tank is 6.75 ft with 2.92 ft in diameter. The tank is filled with 250 gallons of tap water and the height of the water level from the bottom of the tank is approximately 5 ft. In industrial boiler systems, it is very common to use anhydrous sodium sulfite (Na2SO3) in the presence of a catalyst (CoCl2) to scavenge DO in water [23]. Sodium Sulfite ( ) was added to the water which produced Sodium Sulfate . The following Equation shows the chemical reaction achieving oxygen scavenging:  According to Equation (1), the amount of Na 2 SO 3 required can be determined by knowing the water volume in the system, the present DO concentration, and the required level of DO concentration. CoCl 2 was used to increase the reaction speed and the suggested amount was 0.1% of Na 2 SO 3 by weight [23,24]. Theoretically, 7.88 mg/L of Sodium Sulfite is required to scavenge 1 mg/L of dissolved oxygenation concentrations.
The aeration system testing has two steps: (1) deoxygenate the testing water and then (2) reaerate the water while recording the DO concentration. For this purpose, two Vernier optical DO probes (Model Number: ODO-BTA) are used to measure the instantaneous dissolved oxygen concentrations in the aeration tank at two suspended heights. The Vernier DO probe comes with its own data acquisition system where data are collected at a frequency of 1 Hz. The DO probe has the capability of self-calibration based on atmospheric pressure and temperature.
A small scale experimental system was built to see the performance of the venturi injectors with a coiled pipe of varied lengths. The system is shown in Figure 2. The list of equipment required to do the test is listed in Table 1. The setup with a single venturi injector is shown in Figure 1. The list of equipment required to do the test is listed in Table 1.
Water 2020, 12, x FOR PEER REVIEW 4 of 23 According to Equation (1), the amount of Na2SO3 required can be determined by knowing the water volume in the system, the present DO concentration, and the required level of DO concentration. CoCl2 was used to increase the reaction speed and the suggested amount was 0.1% of Na2SO3 by weight [23,24]. Theoretically, 7.88 mg/L of Sodium Sulfite is required to scavenge 1 mg/L of dissolved oxygenation concentrations.
The aeration system testing has two steps: (1) deoxygenate the testing water and then (2) reaerate the water while recording the DO concentration. For this purpose, two Vernier optical DO probes (Model Number: ODO-BTA) are used to measure the instantaneous dissolved oxygen concentrations in the aeration tank at two suspended heights. The Vernier DO probe comes with its own data acquisition system where data are collected at a frequency of 1 Hz. The DO probe has the capability of self-calibration based on atmospheric pressure and temperature.
A small scale experimental system was built to see the performance of the venturi injectors with a coiled pipe of varied lengths. The system is shown in Figure 2. The list of equipment required to do the test is listed in Table 1. The setup with a single venturi injector is shown in Figure 1. The list of equipment required to do the test is listed in Table 1.  The pumps and aeration tanks are on at floor level in the aeration testing lab of The University of Alabama. During tests, the aeration tank is closed and therefore, surface oxygen enrichment is neglected. The system creates a positive pressure head in the pump's suction side. The pump is selfpriming and has its own six-gallon pressurized tank attached to it. The suction side water flow rate  The pumps and aeration tanks are on at floor level in the aeration testing lab of The University of Alabama. During tests, the aeration tank is closed and therefore, surface oxygen enrichment is neglected. The system creates a positive pressure head in the pump's suction side. The pump is self-priming and has its own six-gallon pressurized tank attached to it. The suction side water flow rate is controlled with a ball valve. The discharge side of the pump is attached to the flowmeter and pressure gauge to measure the flow rate and pressure in the discharge side of the pump. Instead of adding the venturi injector directly to the main line, it has been added to the side line and whenever the operation venturi injector is required, the mainline valve is shut and the side line valve is opened. From the discharge side of the injector, the coiled, cross-linked polyethylene (PEX) pipe is attached. Since the PEX pipe is coiled and can be kept hanging over the wall with a mounting, it does not hold much space and the length of the coiled pipe is also variable. Finally, the coiled pipe discharges over the aeration tank, making the outlet thermodynamic pressure to be atmospheric pressure. The pressurized water, when traveling through the convergence part of the injector, experiences a pressure drop to below atmospheric pressure due to the venture effect in the throat, which results in the suction of air through the venturi suction port and a large number of small air bubbles are produced. The suction diameter of the venturi injector is 0.75 inches.
Three individual tests have been conducted with the experimental setup shown in Figure 2. In the first two experiments, the pressure differential between the injector inlet to the outlet is varied from 6.6 psi to 12.4 psi. In the last experiment, the flow rate of the moving fluid was varied with a ball valve from 1.26 cfm to 1.13 cfm along with a reduction in the injector inlet pressure. In all experiments, two oxygen probes are used in different depths of the tank to observe the dissolved oxygen level. The oxygen probes are able to collect one sample of data per second. The variation of dissolved oxygen data from the two probes was insignificant, confirming that the water in the aeration tank had been well mixed. Since, in all cases, the pressure differential between the injector inlet and outlet has been varied, the suction air through the injector hole is also different in each case. Table 2 lists the important experimental parameters that varied throughout the tests.

Theory
This section provides an overview of the theory that is being used in calculating aeration efficiency for the experimental setup. Moreover, this section discusses the fundamentals behind the discrete bubble modeling and the algorithm developed for numerical analysis.

Two Film Theory
Mass transfer in multiphase flow (gas-liquid-solid) processes is a significant phenomenon occurring in chemical, petrochemical, and biological engineering problems. The transport of species among phases is accomplished by two means: diffusion (physical) and chemical reactions. There are several parameters that govern this mass transfer process, for example, system pressure, temperature, concentration gradients, reaction kinetics, conductance of mass transfer, and activation energy.
There are different theories available for dealing with mass transfer among the phases, for example, two-film theory, penetration theory, and the surface renewal theory. Among them, the two-film theory developed by Lewis and Whitman [26] is one of the oldest for gas transfer to a particular liquid. The theory assumes that there exists a stagnant gas and a stagnant liquid film with thickness (δ) which are separated by an interface. The thin films are covered by gas bulk and liquid bulk on their respective sides. To dissolve a gas molecule into the liquid, it has to pass these four distinct regions. The reversed path applies for a gas molecule leaving the liquid. The assumptions made in this theory are: (1) linear concentration profile in stagnant films, (2) mass transfer through the films occurs as steady-state conditions, (3) instantaneous equilibrium, and (4) dilute solutions and Henry's Law is applicable. The theory states that the boundary (or the interface) of gas-liquid phases is two distinct films. Figure 3 shows an illustrative diagram for the two-film theory.
Water 2020, 12, x FOR PEER REVIEW 6 of 23 film theory developed by Lewis and Whitman [26] is one of the oldest for gas transfer to a particular liquid. The theory assumes that there exists a stagnant gas and a stagnant liquid film with thickness (δ) which are separated by an interface. The thin films are covered by gas bulk and liquid bulk on their respective sides. To dissolve a gas molecule into the liquid, it has to pass these four distinct regions. The reversed path applies for a gas molecule leaving the liquid. The assumptions made in this theory are: (1) linear concentration profile in stagnant films, (2) mass transfer through the films occurs as steady-state conditions, (3) instantaneous equilibrium, and (4) dilute solutions and Henry's Law is applicable. The theory states that the boundary (or the interface) of gas-liquid phases is two distinct films. Figure 3 shows an illustrative diagram for the two-film theory.

Experimental Results Evaluation Technique
Bubble aeration is a mass transfer process. The rate of gas transfer is generally proportional to the difference between the existing concentration and the equilibrium concentration of the gas in the solution. This rate of gas transfer is basically governed by the liquid phase mass transfer coefficient, . However, it is non-trivial to determine the interfacial area for mass transfer ( ) per unit volume (V) experimentally. In that case, the overall mass transfer coefficient is used without obtaining the factors and ( ). In short, the overall mass transfer coefficient is written as follows: In Equation (2), is an overall volumetric oxygen mass transfer coefficient in clean water, h −1 , V is water volume in the tank, m 3 and is an interfacial area of mass transfer, m 2 . Moreover, = is interfacial area per unit volume m 2 /m 3 . To determine the oxygen transfer coefficient nonsteady methods were employed [26]. According to two-film theory and Fick's law [27][28][29], the nonsteady oxygen transfer differential equation is: After integration, Equation (3) becomes, where, is the saturation concentration of oxygen under test conditions, is the concentration of oxygen at time, = 0 and is oxygen concentration at time . If the initial concentration in the liquid is = 0, then Equation (4) becomes:

Experimental Results Evaluation Technique
Bubble aeration is a mass transfer process. The rate of gas transfer is generally proportional to the difference between the existing concentration and the equilibrium concentration of the gas in the solution. This rate of gas transfer is basically governed by the liquid phase mass transfer coefficient, K L . However, it is non-trivial to determine the interfacial area for mass transfer (A t ) per unit volume (V) experimentally. In that case, the overall mass transfer coefficient K L a is used without obtaining the factors K L and ( A t V ). In short, the overall mass transfer coefficient is written as follows: In Equation (2), K L a is an overall volumetric oxygen mass transfer coefficient in clean water, h −1 , V is water volume in the tank, m 3 and A t is an interfacial area of mass transfer, m 2 . Moreover, a = A t V is interfacial area per unit volume m 2 /m 3 . To determine the oxygen transfer coefficient K L a non-steady methods were employed [26]. According to two-film theory and Fick's law [27][28][29], the non-steady oxygen transfer differential equation is: After integration, Equation (3) becomes, where, C s is the saturation concentration of oxygen under test conditions, C 0 is the concentration of oxygen at time, t = 0 and C t is oxygen concentration at time t. If the initial concentration in the liquid is C 0 = 0, then Equation (4) becomes: Finally, the total mass transfer coefficient K L a with respect to the total volume of liquid in the system can be calculated in the following way [30]: Water 2020, 12, 1637 7 of 22 Here, K L a is the mass transfer coefficient of the system measured in (h −1 ) at water temperature T. Applying non-linear regression analysis of the logarithmic function of concentration deficit (C s − C t ) with aeration time t, it is possible to determine K L a. K L a is then converted to a standard reference temperature of 20 • C (K L a 20 ) followed by Equation (6): In Equation (7), T represents the water temperature ( • C) and θ is the temperature correction factor which stands for 1.024 for pure water. For any aeration system to compare with the existing system, it is very important to discuss the standard oxygen transfer rate (SOTR). The SOTR is defined as the mass of oxygen transferred to the water volume per unit time when zero initial dissolved oxygen and standard conditions are considered (i.e., water temperature at 20 • C and pressure at 1 atm). SOTR is therefore defined as follows: The efficiency of the aeration system is determined by computing standard aeration efficiency, which is stated as rate of oxygen transfer per unit power input and expressed in the following way: In this Equation, P in stands for the total power drawn from the outlet of the compressor or pump. It is expressed in Kg O 2 /kWh or lbs O 2 /hp-h.

Discrete Bubble Model (Analytical Technique)
The analytical model used here is a discrete bubble model. The discrete bubble model was first used in literature [31] in the phenomena of plug flow through a tank of well-mixed water. Bubbles affect the flow significantly when they are experiencing large pressure gradients and velocity variations. This model tracks every bubble motion as discrete entities and also the bubble deformation due to mass transfer across the bubble surface. The main characteristics of this model are listed below: 1.
The model is based on the Euler-Euler (E-E) approach.

2.
This model is capable of solving homogeneous multiphase flow.

3.
This E-E model employs the volume averaged mass and momentum conservation equation to describe time dependent motion of both phases. 4.
The bubble number contained in a computational cell is represented by a volume fraction.

5.
The bubble size information is obtained by incorporating population balance equations with break-up and coalesce of bubbles as well as growth or shrinkage of bubbles due to mass transfer across the bubble surface.
This model is being extensively used in modeling gas-liquid dispersed flow in bubble column reactors to investigate complex phenomena like hydrodynamics, mass transfer, and chemical reactions [32][33][34]. The benefit of using the Euler-Euler approach is that the number of bubbles is not limited and storage requirement and computational power demand are not dependent on bubble numbers. However, this approach comes with diffusion errors and is sensitive to bubble diameter. The distribution of bubble diameters is either estimated by population balance along with coalescence and break up events [35] or directly obtained from the experiment [36,37]. DBM is applied successfully in airlift aerators [38,39], bubble plume [40], and diffused bubble and wastewater ozonation systems [41,42]. In this present paper, the bubble and fluid is considered as a continuous interpenetrating fluid which gives the opportunity to use the DBM method to predict the rate of oxygen transfer during bubble aeration. The bubbles are considered spherical in shape, which simplifies the forces exerted on the bubbles. In DBM studies, no turbulence of continuous phase is considered. The specific interface area is significant since the mass transfer among the two phases is dependent on it. In fact, the bubble size distribution is time-dependent, as processes like bubble coalesce and break up affect this phenomenon. Consequently, this affects the interface area and mass transfer. In this study, the size of the bubble is a function of the pressure gradient and mass transfer only, and the effects of bubbles coalescing and breaking up is overlooked.
The bubbles are assumed to be moving with the fluid with an average mixture velocity. For each experiment, the number of bubbles in the pipe at any given time is assumed to be invariable considering that no input of the system has been changed. The bubbles are simply experiencing pressure forces inside the pipe which varies along the pipe length. Surface tension and liquid head pressure is negligible in this case. For simplifying the model, it is assumed that bubble break-up and coalescing are not occurring. The flow is incompressible with the volume-average mixture properties (e.g., viscosity, density) being used. The distribution of the bubble size is represented by a single Sauter mean diameter [43]. The water and air temperature are assumed to be equal and constant. The mass transfer flux across the surface of a bubble is: where, K L is liquid side mass transfer coefficient, C S is equilibrium concentration at the gas/water interface, C is bulk aqueous-phase concentration, and J is mass transfer flux. The gas side mass transfer resistance has been neglected. The equilibrium concentration can be expressed as: where, H is Henry's constant and P i is partial pressure of the gas. With this, Equation (10) can be written as follows: Henry's constant (H) (mol.m −3 .bar −1 ) and mass transfer coefficient (ms −1 ) are determined in the following way [44]: and The rate of mass transfer across the surface area of the bubble of radius (r) can be obtained from the mass transfer flux Equation: The discrete bubble is moving with the water with mixture velocity v m . At any time, the velocity of the bubble inside the pipe along the horizontal axis can be related by the following equation: Since the bulk aqueous-phase concentration does not change significantly during the travel of bubbles inside the pipe, the pseudo-steady state assumption may be considered. From this equation, it is possible to determine the mass transfer of gaseous species per bubble per unit length of pipe: Water 2020, 12, 1637 9 of 22 The number of bubbles in the pipe at any instant is constant. This number of bubbles is calculated by knowing the volumetric flow of air through the venturi injector, Q 0 and the initial bubble volume, V 0 .
Multiplying the mass transfer equation for a single bubble with the total number of bubbles (N) will give the total mass transfer per unit length at per unit time.
The above equation is a one-dimensional ordinary differential equation. On the right hand side, the variables are H, which is function of temperature only and K L , which is a function of bubble radius.
Bubble radius will change with the change in pressure inside the pipe and also with the mass transfer of oxygen and nitrogen gas across the bubble and liquid interface. As said before, the ideal gas law is used in this case to measure the change of bubble radius with time.
In this experiment, a coiled pipe is used to transfer water from the source to the aeration tank. To include the effect of the coiled pipe in the calculation, the friction loss of this type of pipe is being calculated using the Mishra and Gupta correction factor with Blasius correlations [45]. The friction factor for curved or helically coiled tubes, taking into consideration the equivalent curve radius, R C , is as follows: where, D P , R, H, and Re represent pipe diameter, curve radius, helicoidal pitch, and Reynolds Number respectively. In the current analysis, the two phase flow is modeled as the homogeneous flow model which is the simplest technique for analyzing multiphase flows. In the homogeneous model, both liquid and gas phases move at the same velocity (slip ratio = 1). It is also known as the zero slip model. The homogeneous model considers the two-phase flow as a single-phase flow having the average fluid properties in which the properties depend upon the mixture quality. For this purpose, one parameter which is known as mass quality is defined as follows: where, Q water and Q air stand for volumetric flow rate of water (m 3 /s) and volumetric flow rate of air (m 3 /s) respectively. The definition of two phase viscosity (µ m ) based on the mass averaged value of the reciprocals of the gas viscosity (µ g ) and liquid viscosity (µ l ) is defined as follows: The mixture density (ρ m ) also obtained in the mass averaged process of gas density (ρ g ) and liquid density (ρ l ) is calculated as follows: Water 2020, 12, 1637 10 of 22 The velocity and the Reynolds number of the mixture is calculated in the following way: The Reynolds number of the mixture can be found in this way: where, A pipe is pipe area which is equal to π 4 D p 2 (m 2 ).
Finally, the frictional head loss (h f ) across the length of the coiled pipe (l) is calculated by using mixture density, velocity, and Reynolds number which is written as follows: Gordiychuk [46] have performed an experiment to calculate the Sauter mean diameter of the bubble (d 32 ) generated in venturi-type bubble generators. They have derived the relation to determine the Sauter mean diameter as a function of water flow rate, air flow rate, and air inlet diameter. The relation is given in the following way (see Table 2 of [46] for ML-BGS algorithm): In this relation, α stands for air to water fraction. Re w and Re air are Reynolds number for water and air respectively. The Reynolds number of air is evaluated based on suction diameter (d inlet,(suction) ), suction velocity, and air density. On the other hand, for water, Reynolds number is calculated based on pipe diameter, water velocity and water density. Now that all necessary equations to calculate mass transfer in the discrete bubble model have been documented, an overview of the algorithm structure used for this model is shown in Figure 4: The Reynolds number of the mixture can be found in this way: where, is pipe area which is equal to (m 2 ).
Finally, the frictional head loss (hf) across the length of the coiled pipe (l) is calculated by using mixture density, velocity, and Reynolds number which is written as follows: Gordiychuk [46] have performed an experiment to calculate the Sauter mean diameter of the bubble ( ) generated in venturi-type bubble generators. They have derived the relation to determine the Sauter mean diameter as a function of water flow rate, air flow rate, and air inlet diameter. The relation is given in the following way (see Table 2 of [46] for ML-BGS algorithm): In this relation, stands for air to water fraction. and are Reynolds number for water and air respectively. The Reynolds number of air is evaluated based on suction diameter ( ,( ) ), suction velocity, and air density. On the other hand, for water, Reynolds number is calculated based on pipe diameter, water velocity and water density. Now that all necessary equations to calculate mass transfer in the discrete bubble model have been documented, an overview of the algorithm structure used for this model is shown in Figure 4:

Discussion
In this section, three different experimental results are shown. These results are compared with the results obtained from the discrete bubble model to determine the model efficiency in predicting the oxygen transfer efficiency of the discussed system. In the end, the discrete bubble model is used to predict a proposed aeration system for an aquaculture farm to obtain an SAE over 1.5. The following section is divided in three subsections: the first subsection will validate the analytical model, the second section discusses the effects of different parameters on oxygen transfer efficiency using the discrete bubble model, and the third will cover a proposed aeration system suitable for aquaculture.

Validation of the Numerical Model
In all experiments, the air suction through the venturi hole was natural. The air suction is dependent on pressure differential and water flow rates for a specific injector. Mazzei documented the air suction rate through the injector hole for different scenarios in their product catalogue for the specific injector model. Since then, the pressure differentials and water flow rate conditions for different experiments in this study are close to what is mentioned in the product catalogue, the air suction rate is not measured but rather is directly used from the catalogue. This also gives the independence of testing different scenarios in numerical studies as it subsides the requirement of measuring suction flow rate. While doing numerical analysis, in each case, the water flow rate is known and the code is also capable of calculating the injector inlet pressure head from different hydraulic parameters. As a result of this, varieties of pressure of differentials can be tested along with different water flow conditions.
The analytical results were obtained using the same physical parameters as the experiments. Data from the two probes are averaged over the time period. Figure 5 shows the average dissolved oxygen vs. time results for three different experiments. The analytical results obtained through discrete bubble models for each case are plotted in the same graph to show their conformity with the experimental results.
Water 2020, 12, x FOR PEER REVIEW 11 of 23 In this section, three different experimental results are shown. These results are compared with the results obtained from the discrete bubble model to determine the model efficiency in predicting the oxygen transfer efficiency of the discussed system. In the end, the discrete bubble model is used to predict a proposed aeration system for an aquaculture farm to obtain an SAE over 1.5. The following section is divided in three subsections: the first subsection will validate the analytical model, the second section discusses the effects of different parameters on oxygen transfer efficiency using the discrete bubble model, and the third will cover a proposed aeration system suitable for aquaculture.

Validation of the Numerical Model
In all experiments, the air suction through the venturi hole was natural. The air suction is dependent on pressure differential and water flow rates for a specific injector. Mazzei documented the air suction rate through the injector hole for different scenarios in their product catalogue for the specific injector model. Since then, the pressure differentials and water flow rate conditions for different experiments in this study are close to what is mentioned in the product catalogue, the air suction rate is not measured but rather is directly used from the catalogue. This also gives the independence of testing different scenarios in numerical studies as it subsides the requirement of measuring suction flow rate. While doing numerical analysis, in each case, the water flow rate is known and the code is also capable of calculating the injector inlet pressure head from different hydraulic parameters. As a result of this, varieties of pressure of differentials can be tested along with different water flow conditions.
The analytical results were obtained using the same physical parameters as the experiments. Data from the two probes are averaged over the time period. Figure 5 shows the average dissolved oxygen vs. time results for three different experiments. The analytical results obtained through discrete bubble models for each case are plotted in the same graph to show their conformity with the experimental results.    Table 2: (a) Test-1; (b) Test-2; (c) Test-3. Figure 5 shows the increase in dissolved oxygen in the aeration tank with time. From the experimental results, it can be observed that the first test reaches saturation more quickly than the other two. The second test takes the longest time to reach saturation. This is due to the fact that in the second test, less air is drawn into the aerator. The analytical results are reasonably able to predict the measured DO. In most cases, the analytical model results fall within 10-15% of the measured data, although it underestimates the results. In the beginning of the test, the analytical results divert mostly from the experimental results, but with time the measured and predicted values started to conform more with each other. This difference in predicted results with measured values could happen for several reasons. One of the main factors is the degree of variability of the bubble size in real measurement compared to the analytical model. In the analytical model, all bubbles are considered spherical in shape which is not true in real phenomena. Bubble radius also changes due to coalesce and breakage. Moreover, the effect of turbulent mixing inside the pipe after the injector is not included in the model. The velocity of the bubble and water is assumed to be the same with having a mixture velocity which might be another source of error. This average velocity is added in the model without considering the time-dependent velocity change due to the turbulent flow inside the pipe.
Based on K L a, the standard aeration efficiency (SAE) and the standard oxygen transfer rate (SOTR) are determined. All values of K L a, SOTR and SAE are presented in Table 3 for both the experiment and analytical results. By comparing the experimental results, the highest SAE was achieved for the third test, while the maximum oxygen transfer rate was observed for test-1. This is also valid for the analytical results. According to Equation (6), a semi log-based graph is plotted with respect to time in Figure 6. The slope of the curve represents the oxygen transfer coefficient (K L a) of the aerator module which is an important parameter in determining the efficiency of the aerator. Oxygen transfer coefficient expresses the relationship between the rate of mass transfer and the concentration deficit which is the driving force to dissolve oxygen in the water. In this case, the K L a values obtained through experiment under standard conditions for test-1, test-2, and test-3 are 0.067 min −1 , 0.048 min −1 , and 0.053 min −1 respectively. The highest oxygen transfer coefficient is achieved through test-1. In the analytical model, a similar log deficit plot is generated to determine the K L a. The analytical model K L a are −0.056 min −1 , 0.042 min −1 , and 0.047 min −1 for test-1,-2, and -3 respectively. In all cases, the analytical model K L a is less than the experimental K L a values and the difference in them are between 10-16%.
Water 2020, 12, x FOR PEER REVIEW 12 of 23 Figure 5 shows the increase in dissolved oxygen in the aeration tank with time. From the experimental results, it can be observed that the first test reaches saturation more quickly than the other two. The second test takes the longest time to reach saturation. This is due to the fact that in the second test, less air is drawn into the aerator. The analytical results are reasonably able to predict the measured DO. In most cases, the analytical model results fall within 10-15% of the measured data, although it underestimates the results. In the beginning of the test, the analytical results divert mostly from the experimental results, but with time the measured and predicted values started to conform more with each other. This difference in predicted results with measured values could happen for several reasons. One of the main factors is the degree of variability of the bubble size in real measurement compared to the analytical model. In the analytical model, all bubbles are considered spherical in shape which is not true in real phenomena. Bubble radius also changes due to coalesce and breakage. Moreover, the effect of turbulent mixing inside the pipe after the injector is not included in the model. The velocity of the bubble and water is assumed to be the same with having a mixture velocity which might be another source of error. This average velocity is added in the model without considering the time-dependent velocity change due to the turbulent flow inside the pipe.
Based on , the standard aeration efficiency (SAE) and the standard oxygen transfer rate (SOTR) are determined. All values of , SOTR and SAE are presented in Table 3 for both the experiment and analytical results. By comparing the experimental results, the highest SAE was achieved for the third test, while the maximum oxygen transfer rate was observed for test-1. This is also valid for the analytical results. According to Equation (6), a semi log-based graph is plotted with respect to time in Figure 6. The slope of the curve represents the oxygen transfer coefficient ( ) of the aerator module which is an important parameter in determining the efficiency of the aerator. Oxygen transfer coefficient expresses the relationship between the rate of mass transfer and the concentration deficit which is the driving force to dissolve oxygen in the water. In this case, the values obtained through experiment under standard conditions for test-1, test-2, and test-3 are 0.067 min −1 , 0.048 min −1 , and 0.053 min −1 respectively. The highest oxygen transfer coefficient is achieved through test-1. In the analytical model, a similar log deficit plot is generated to determine the . The analytical model are −0.056 min −1 , 0.042 min −1 , and 0.047 min −1 for test-1,-2, and -3 respectively. In all cases, the analytical model is less than the experimental values and the difference in them are between 10-16%.  Table 2: (a) Test-1; (b) Test-2; (c) Test-3.

Parametric Study
Here, the effects of different hydraulic and flow parameters have been tested with the discrete bubble model. This study gives an insight into which parameters are affecting the oxygen transfer coefficient ( ), SOTR and SAE, and will help find an optimum input for future, larger-scale systems for aeration purposes. In all cases, test-1 conditions are used while varying test-1 parameters are used to investigate effects on aeration efficiency in varying these parameters.

Effect of Pipe Length
One of the main purposes of this study was to determine the effect of coiled pipe length in oxygen transfer efficiency. The coiled pipe will start after the injector and will hang over the aeration basin. According to [47], in the helical tube, bubbles dispersed mostly around the center of the coil and a drag reduction was observed on the bubble surface. In the analytical model, the same inputs are used from test-1 with only a change in the pipe length.
In Figure 7a, the DO concentration profile versus time is plotted to show the effect of pipe length. The profile conforms to a general non-steady aeration curve where the rate of change in dissolved oxygen levels decrease with time. This is due to a reduction in saturation deficit level with time. It can be observed that, with an increase in pipe length, the system takes less time to reach the saturation level. Since, with an increase in pipe length, bubbles are taking a much longer retention time to transfer oxygen with water, the oxygen transfer rate is also increased. Moreover, the bubbles remain in the center and thus less likely to collide with the pipe surface. This helps them to maintain a spherical shape for a longer period of time and increases the overall mass transfer rate. This mass transfer rate increase comes with a small increase in the friction head loss which causes less pump power rise. In effect, the overall oxygenation transfer efficiency increases. However, this positive effect of pipe length increase dampens after certain intervals (in this case, from 80 ft. to 100 ft.) (a) (b) Figure 6. Determination of oxygen transfer coefficient for three test conditions defined in Table 2: (a) Test-1; (b) Test-2; (c) Test-3.

Parametric Study
Here, the effects of different hydraulic and flow parameters have been tested with the discrete bubble model. This study gives an insight into which parameters are affecting the oxygen transfer coefficient (K L a), SOTR and SAE, and will help find an optimum input for future, larger-scale systems for aeration purposes. In all cases, test-1 conditions are used while varying test-1 parameters are used to investigate effects on aeration efficiency in varying these parameters.

Effect of Pipe Length
One of the main purposes of this study was to determine the effect of coiled pipe length in oxygen transfer efficiency. The coiled pipe will start after the injector and will hang over the aeration basin. According to [47], in the helical tube, bubbles dispersed mostly around the center of the coil and a drag reduction was observed on the bubble surface. In the analytical model, the same inputs are used from test-1 with only a change in the pipe length.
In Figure 7a, the DO concentration profile versus time is plotted to show the effect of pipe length. The profile conforms to a general non-steady aeration curve where the rate of change in dissolved oxygen levels decrease with time. This is due to a reduction in saturation deficit level with time. It can be observed that, with an increase in pipe length, the system takes less time to reach the saturation level. Since, with an increase in pipe length, bubbles are taking a much longer retention time to transfer oxygen with water, the oxygen transfer rate is also increased. Moreover, the bubbles remain in the center and thus less likely to collide with the pipe surface. This helps them to maintain a spherical shape for a longer period of time and increases the overall mass transfer rate. This mass transfer rate increase comes with a small increase in the friction head loss which causes less pump power rise. In effect, the overall oxygenation transfer efficiency increases. However, this positive effect of pipe length increase dampens after certain intervals (in this case, from 80 ft. to 100 ft.) Water 2020, 12, x FOR PEER REVIEW 13 of 23 (c) Figure 6. Determination of oxygen transfer coefficient for three test conditions defined in Table 2: (a) Test-1; (b) Test-2; (c) Test-3.

Parametric Study
Here, the effects of different hydraulic and flow parameters have been tested with the discrete bubble model. This study gives an insight into which parameters are affecting the oxygen transfer coefficient ( ), SOTR and SAE, and will help find an optimum input for future, larger-scale systems for aeration purposes. In all cases, test-1 conditions are used while varying test-1 parameters are used to investigate effects on aeration efficiency in varying these parameters.

Effect of Pipe Length
One of the main purposes of this study was to determine the effect of coiled pipe length in oxygen transfer efficiency. The coiled pipe will start after the injector and will hang over the aeration basin. According to [47], in the helical tube, bubbles dispersed mostly around the center of the coil and a drag reduction was observed on the bubble surface. In the analytical model, the same inputs are used from test-1 with only a change in the pipe length.
In Figure 7a, the DO concentration profile versus time is plotted to show the effect of pipe length. The profile conforms to a general non-steady aeration curve where the rate of change in dissolved oxygen levels decrease with time. This is due to a reduction in saturation deficit level with time. It can be observed that, with an increase in pipe length, the system takes less time to reach the saturation level. Since, with an increase in pipe length, bubbles are taking a much longer retention time to transfer oxygen with water, the oxygen transfer rate is also increased. Moreover, the bubbles remain in the center and thus less likely to collide with the pipe surface. This helps them to maintain a spherical shape for a longer period of time and increases the overall mass transfer rate. This mass transfer rate increase comes with a small increase in the friction head loss which causes less pump power rise. In effect, the overall oxygenation transfer efficiency increases. However, this positive effect of pipe length increase dampens after certain intervals (in this case, from 80 ft. to 100 ft.) , SOTR, and SAE are plotted in the graph from Figure 7b-d with respect to pipe length. In every case, the linear trend line is obtained which passes through most of the values. If the pipe length is doubled, the SOTR is also increased by a factor of 1.4. For a 100 ft. section of pipe, the oxygenation efficiency is 2.6 times better than the 20 ft. pipe. Values of , SOTR, and SAE are listed in Table 4.

Effect of Pipe Diameter
In the current research, only 1-inch diameter pipe is being used for all experiments. However, one-inch pipe usually has a higher frictional loss and thus affects the SAE in a negative way. In Figure  8, effects of diameter change in dissolved oxygen rate and in efficiency are shown. K L a, SOTR, and SAE are plotted in the graph from Figure 7b-d with respect to pipe length. In every case, the linear trend line is obtained which passes through most of the values. If the pipe length is doubled, the SOTR is also increased by a factor of 1.4. For a 100 ft. section of pipe, the oxygenation efficiency is 2.6 times better than the 20 ft. pipe. Values of K L a, SOTR, and SAE are listed in Table 4.

Effect of Pipe Diameter
In the current research, only 1-inch diameter pipe is being used for all experiments. However, one-inch pipe usually has a higher frictional loss and thus affects the SAE in a negative way. In Figure 8 From Figure 8, it is evident that the dissolved oxygen level rises with time in the aeration tank. The DO level rises rapidly at first, then this rate of increase decreases as it gets close to saturation. It follows that increasing the diameter helps the system reach the saturation more quickly. This can be attributed to the fact that increasing the pipe diameter decreases the mixture velocity. In real phenomena, there exists a slip velocity due to the density difference between two phase flow. The gas bubbles do not exactly travel with the liquid [48]. Additionally, in smaller diameter pipe more bubbles tend to travel towards the pipe wall. Oxygen transfers to the water through diffusion. To overcome the resistance in gas-water interface, bubbles need sufficient time. Since the water velocity decreases with the diameter increase, the gas bubbles get more time for mass transfer to overcome the diffusion resistance. This positive effect can also be observable while determining the oxygenation coefficient as it increases with the diameter increase. In Figure 8b, c, , and are plotted against the diameter and are fitted with a trend line. This trend line equation is a simple mathematical model of the relationship between the diameter with oxygenation coefficient and SOTR. SAE also attains higher values as the system is able to transfer oxygen faster without changing the operating condition of the pump.

Effect of Gas Flow Rate
This section will discuss the effect of increasing the suction gas flow rate while keeping the motive flow rate and the line pressure constant. In this case, the test-1 condition is used only in changing the gas suction rate.
It can be observed from Figure 9 that in all cases of simulation, gas suction flow rates increase significantly the impact of dissolved oxygen level. With the increase in gas flow rate, the dissolved oxygen attains the saturation concentration more rapidly. Gas flow rate basically increases the mass quality in Equation (24). This positive impact can be explained in two ways: firstly, an increase in gas flow rate generally creates more turbulence in the liquid interface. At this point, mass transfer is From Figure 8, it is evident that the dissolved oxygen level rises with time in the aeration tank. The DO level rises rapidly at first, then this rate of increase decreases as it gets close to saturation. It follows that increasing the diameter helps the system reach the saturation more quickly. This can be attributed to the fact that increasing the pipe diameter decreases the mixture velocity. In real phenomena, there exists a slip velocity due to the density difference between two phase flow. The gas bubbles do not exactly travel with the liquid [48]. Additionally, in smaller diameter pipe more bubbles tend to travel towards the pipe wall. Oxygen transfers to the water through diffusion. To overcome the resistance in gas-water interface, bubbles need sufficient time. Since the water velocity decreases with the diameter increase, the gas bubbles get more time for mass transfer to overcome the diffusion resistance. This positive effect can also be observable while determining the oxygenation coefficient as it increases with the diameter increase. In Figure 8b,c, K L a, and are plotted against the diameter and are fitted with a trend line. This trend line equation is a simple mathematical model of the relationship between the diameter with oxygenation coefficient and SOTR. SAE also attains higher values as the system is able to transfer oxygen faster without changing the operating condition of the pump.

Effect of Gas Flow Rate
This section will discuss the effect of increasing the suction gas flow rate while keeping the motive flow rate and the line pressure constant. In this case, the test-1 condition is used only in changing the gas suction rate.
It can be observed from Figure 9 that in all cases of simulation, gas suction flow rates increase significantly the impact of dissolved oxygen level. With the increase in gas flow rate, the dissolved oxygen attains the saturation concentration more rapidly. Gas flow rate basically increases the mass quality in Equation (24). This positive impact can be explained in two ways: firstly, an increase in gas flow rate generally creates more turbulence in the liquid interface. At this point, mass transfer is governed by the disturbance and rate of restoration of the liquid film. This is described by the Danckwertz surface rejuvenation theory [49]: where, K L is the liquid film coefficient discussed in Equation (10), D L is the diffusion coefficient of oxygen (m 2 ·h −1 ) and r is the rate of restoration of liquid the film (h −1 ). Since an increase in gas flow rate increases the r value in Equation (31), K L also increases. This positive effect of K L also increases the oxygen transfer coefficient and SOTR.
Water 2020, 12, x FOR PEER REVIEW 16 of 23 governed by the disturbance and rate of restoration of the liquid film. This is described by the Danckwertz surface rejuvenation theory [49]: where, is the liquid film coefficient discussed in Equation (10), is the diffusion coefficient of oxygen (m 2 ·h −1 ) and is the rate of restoration of liquid the film (h −1 ). Since an increase in gas flow rate increases the r value in Equation (31), also increases. This positive effect of also increases the oxygen transfer coefficient and SOTR. Secondly, a rise in the gas flow rate also increases the number of bubbles at any time. This increase in gas bubbles raises the total available interfacial area to the surrounding medium for diffusion [50][51][52]. In short, this creates favorable conditions for and SOTR. In this discussion, SAE also increases since no flow parameters have been changed. However, in real life this might not be true as the increase in pressure differential and motive flow rate is required to increase the suction gas flow rate for existing venturi. Another way that a gas suction increase is possible is if a new and improved venturi aerator was available. This is currently other promising research in the field of venturi aeration.

Effect of Liquid Flow Rate
In this section, the discussion is made on the effect of liquid flow rate in overall performance of the aeration system. This effect is found ignoring other changes associated with liquid flow rate change like pressure head decreases in pump. The rest of the other conditions of this simulation are same as in test-1.
In Figure 10a, it can be seen that a liquid flow rate increase has adverse effects on the dissolve oxygen rise with time. As the liquid flow rate increases, the system takes more time to reach saturation. On the other hand, no observable trend has been seen for the oxygen transfer coefficient and SOTR. Both the and the SOTR plots against the liquid flow rate have a downward slope Secondly, a rise in the gas flow rate also increases the number of bubbles at any time. This increase in gas bubbles raises the total available interfacial area to the surrounding medium for diffusion [50][51][52]. In short, this creates favorable conditions for K L a and SOTR. In this discussion, SAE also increases since no flow parameters have been changed. However, in real life this might not be true as the increase in pressure differential and motive flow rate is required to increase the suction gas flow rate for existing venturi. Another way that a gas suction increase is possible is if a new and improved venturi aerator was available. This is currently other promising research in the field of venturi aeration.

Effect of Liquid Flow Rate
In this section, the discussion is made on the effect of liquid flow rate in overall performance of the aeration system. This effect is found ignoring other changes associated with liquid flow rate change like pressure head decreases in pump. The rest of the other conditions of this simulation are same as in test-1.
In Figure 10a, it can be seen that a liquid flow rate increase has adverse effects on the dissolve oxygen rise with time. As the liquid flow rate increases, the system takes more time to reach saturation. On the other hand, no observable trend has been seen for the oxygen transfer coefficient and SOTR.
Both the K L a and the SOTR plots against the liquid flow rate have a downward slope which shows the performance degradation of the system with the flow rate. The reasons should be opposite those of the gas flow rate discussion. One factor is that the increasing liquid flow rate decreases the average gas mass fraction at any time. Therefore, less interfacial area is available for mass transfer through diffusion. Additionally, a higher liquid flow rate decreases the bubble contact time. All of these parameters also affect the aeration efficiency of the system negatively.
Water 2020, 12, x FOR PEER REVIEW 17 of 23 which shows the performance degradation of the system with the flow rate. The reasons should be opposite those of the gas flow rate discussion. One factor is that the increasing liquid flow rate decreases the average gas mass fraction at any time. Therefore, less interfacial area is available for mass transfer through diffusion. Additionally, a higher liquid flow rate decreases the bubble contact time. All of these parameters also affect the aeration efficiency of the system negatively.

Proposed System
Based on all the previous analysis done in this research, a new system is proposed with the aim to reach an SAE of over 1.5. The proposed system is considered for an assumed real-life problem. Therefore, the size of the system is much larger than the system tested. The problem considered here is that the proposed system needs to aerate an aquaculture (catfish) pond the size of two acres. The average depth of the system is four feet. This gives approximately 26,042 gallons of water to aerate. This analysis is still done based on the assumption of clean water. In the previous section, it was observed that diameter and airflow rate increase have the maximum impact on increasing the efficiency of the aeration system. However, the airflow rate increase comes with higher pressure differentials and motive flow which in turn negatively impacts the SAE. In that case, pipe length increase plays a better role, since frictional head loss is much less through PVC or PEX pipe. A pump size of 5 hp is chosen in this case after calculating the total head developed by the system. Finally, multiple simulations have been done in the analytical model to determine the optimum system to aerate this catfish pond. The input parameters that are found to be useful to get an SAE over 1.5 are given in Table 5. Figure 11a,b represent the dissolved oxygen level and log deficit concentration of oxygen with respect to time in the catfish pond for the proposed setup respectively.

Proposed System
Based on all the previous analysis done in this research, a new system is proposed with the aim to reach an SAE of over 1.5. The proposed system is considered for an assumed real-life problem. Therefore, the size of the system is much larger than the system tested. The problem considered here is that the proposed system needs to aerate an aquaculture (catfish) pond the size of two acres. The average depth of the system is four feet. This gives approximately 26,042 gallons of water to aerate. This analysis is still done based on the assumption of clean water. In the previous section, it was observed that diameter and airflow rate increase have the maximum impact on increasing the efficiency of the aeration system. However, the airflow rate increase comes with higher pressure differentials and motive flow which in turn negatively impacts the SAE. In that case, pipe length increase plays a better role, since frictional head loss is much less through PVC or PEX pipe. A pump size of 5 hp is chosen in this case after calculating the total head developed by the system. Finally, multiple simulations have been done in the analytical model to determine the optimum system to aerate this catfish pond. The input parameters that are found to be useful to get an SAE over 1.5 are given in Table 5. Figure 11a,b represent the dissolved oxygen level and log deficit concentration of oxygen with respect to time in the catfish pond for the proposed setup respectively. The Mazzei Injector has a catalogue of their different venturi injectors with the best possible air suction rates at different pressure differentials and motive flow rates. From there, a four-inch injector was chosen which is capable of producing an airflow rate of 7 scfm at the pressure differentials and water flow rates described in Table 5. According to these input values, the analytical model gives an oxygen transfer coefficient of 0.045 min −1 , an SOTR of 4.38 lbs O2/hr, and an SAE of 1.8. The possible setup of the proposed system: 1. The proposed system will be a closed loop system where a pump will be set up on the ground. The water intake of the pump will be at the bottom of the catfish pond. 2. The water will reach the main discharge line located on the ground. There will be a pressure gauge to measure the pressure of the upcoming fluid before it enters into the injector. 3. A pressure reducing valve will be attached to the proposed system in the discharge side of the injector to maintain the desired pressure differential. Therefore, maximum air suction can be achieved. 4. The coiled pipe will start on the discharge side of injector and will hang in from the ground. The pressurized water will be evacuated from the pond on top surface. 5. Finally, the water from the discharge side will be ejected into the catfish pond with high velocity.
A schematic of the proposed system is shown in Figure 12. The Mazzei Injector has a catalogue of their different venturi injectors with the best possible air suction rates at different pressure differentials and motive flow rates. From there, a four-inch injector was chosen which is capable of producing an airflow rate of 7 scfm at the pressure differentials and water flow rates described in Table 5. According to these input values, the analytical model gives an oxygen transfer coefficient of 0.045 min −1 , an SOTR of 4.38 lbs O 2 /hr, and an SAE of 1.8.
The possible setup of the proposed system: 1.
The proposed system will be a closed loop system where a pump will be set up on the ground. The water intake of the pump will be at the bottom of the catfish pond.

2.
The water will reach the main discharge line located on the ground. There will be a pressure gauge to measure the pressure of the upcoming fluid before it enters into the injector.

3.
A pressure reducing valve will be attached to the proposed system in the discharge side of the injector to maintain the desired pressure differential. Therefore, maximum air suction can be achieved.

4.
The coiled pipe will start on the discharge side of injector and will hang in from the ground. The pressurized water will be evacuated from the pond on top surface.

5.
Finally, the water from the discharge side will be ejected into the catfish pond with high velocity. A schematic of the proposed system is shown in Figure 12.
In Table 6, efficiency ranges for various types of aeration equipment obtained from [53] are listed. However, SAE values in Table 6 can be used only as a reference; exact values are equipment and application specific. In this present study, only one aerator is considered. Since multiple aerators attached in parallel form can cause more air to be injected into water, this type of setup can also be used to increase the efficiency of the overall aeration system [54,55]. Future studies could focus on the effect of multiple parallel aerators on aeration efficiency. In Table 6, efficiency ranges for various types of aeration equipment obtained from [53] are listed. However, SAE values in Table 6 can be used only as a reference; exact values are equipment and application specific. In this present study, only one aerator is considered. Since multiple aerators attached in parallel form can cause more air to be injected into water, this type of setup can also be used to increase the efficiency of the overall aeration system [54,55]. Future studies could focus on the effect of multiple parallel aerators on aeration efficiency.

Conclusions
In this study, a novel closed loop aeration system is proposed, consisting of a venturi air injector and a coiled pipe coupled with a pump. Several conclusions have been made from this study which are summarized below: 1. The experimental study of the current closed loop venturi aspirated aeration system achieved an SOTR in the range of 0.050-0.068 lbs O2/h and an SAE in the range of 0.42-0.63 lbs O2/hp-h. With the analytical study, SOTR is found to be similar for this system and is 0.047-0.059 lbs O2/h and SAE in the range of 0.37-0.55 lbs O2/hp-h. In all cases, the difference between experimental and analytical results are within 15%. The analytical model underestimates the experimental Figure 12. A schematic of the proposed system with one aerator.

Conclusions
In this study, a novel closed loop aeration system is proposed, consisting of a venturi air injector and a coiled pipe coupled with a pump. Several conclusions have been made from this study which are summarized below:

1.
The experimental study of the current closed loop venturi aspirated aeration system achieved an SOTR in the range of 0.050-0.068 lbs O 2 /h and an SAE in the range of 0.42-0.63 lbs O 2 /hp-h. With the analytical study, SOTR is found to be similar for this system and is 0.047-0.059 lbs O 2 /h and SAE in the range of 0.37-0.55 lbs O 2 /hp-h. In all cases, the difference between experimental and analytical results are within 15%. The analytical model underestimates the experimental values, likely because it cannot account for turbulent mixing, bubble breakup, and bubble coalescing.

2.
The hydraulic parameters considered herein are the pipe diameter and pipe length. Both significantly impact the system's aeration efficiency. Increasing the pipe length by a factor of five increases the SAE by 2.5 times. For pipe diameter increase, the effect is more evident, since increasing the diameter by 4 times increases the SAE by 4.1 times. The linear trend line is obtained with a simple mathematical relationship between SAE against the pipe length and diameter.