Numerical Investigation of a Liquid-Gas Ejector Used for Shipping Ballast Water Treatment

Song X, Cao MS, Shin W, Cao WP, Kang S, Park Y. Numerical Investigation of a Liquid-Gas Ejector Used for Shipping Ballast Water Treatment. Mathematical Problems in Engineering 2014. Copyright: Copyright © 2014 Xueguan Song et al. This an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. DOI link to article: http://dx.doi.org/10.1155/2014/259593 Date deposited: 07-11-2014


Introduction
Ballast water is a huge amount of water (sometimes millions of gallons on a large ship) stored in the long-distance passenger and freight ships to provide momentum for ship stability.When the ship reaches its destination, ballast water needs to be discharged and renewed for the next journey.However, the discharges ballast water usually contains a wide range of exotic species such as animals, plants, bacteria, and viruses, which may pose a threat to the local marine environment [1].These foreign materials often include invasive exotic species that can cause extensive ecological damage and economic losses to aquatic ecosystems.The worldwide concern over the potential harm of ecosystems from shipping ballast water has resulted in the development of a range of treatment systems to minimize these environmental risks.Currently, ballast water treatment solutions can be classified into four types: mechanical treatment, physical treatment, chemical treatment, and a combination of them [1][2][3].Mechanical treatment uses fine filters to remove small immature young eggs of invasive species or uses high pressure jet to kill them.Physical treatment uses heat, ultraviolet, sonic, and other radiations to kill organisms.Chemical treatment kills existing organisms with chlorine bleach and other toxic chemicals.The combination of them has two or more treatment methods to clear the ballast water.Each of them has its advantages and disadvantages.
In the conventional ozone treatment systems, ozone gas is directly bubbled into the water which decomposes and reacts with chemicals.It is very effective to kill microorganisms, but not so at killing larger organisms.In order to kill both large organisms and microorganisms, a new ozone treatment system needs to be developed.Figure 1  of the proposed ballast water treatment system in this work.This system still employs ozone to kill the microorganisms.However, ozone is no longer directly pumped into the ballast water but is absorbed and ejected into ballast water by an ejector (see Figure 2).Specifically, a small amount of ballast water from the main pipe is pumped as the primary fluid through the ejector to produce a lower pressure in the suction chamber of the ejector, which absorbs ozone as the secondary fluid into the ejector and then ejects it into the downstream of the main pipe to clean the ballast water.Microorganisms are killed by an appropriate mixed fluid downstream, whilst large organisms are killed by the high-velocity impact of the ejected flow.This design has advantages of high efficiency in producing intense mixing flow and a high interfacial area for generating small bubbles to rapidly treat a large volume of ballast water.However, the main disadvantage is that the efficiency of this system is strongly dependent on the performance of the ejector.If inappropriately designed or equipped with a high back pressure, the ejector will not absorb enough ozone to treat the discharged ballast water, resulting in microorganism survivals and their escape to the local water.
An ejector can be single-phase (e.g., liquid-liquid and gas-gas ejector) or multiphase (e.g., liquid-gas and liquidsteam ejectors) [4][5][6][7][8][9][10][11][12][13].For example, detailed experiments and computational fluid dynamics (CFD) numerical studies were undertaken to understand the hydrodynamic characteristics of the ejector geometry [6].These CFD models provide an insight into the ejector performance under different operating conditions.Test results show that there is an optimum ratio of the nozzle area to the throat area (area ratio), at which the liquid entrainment rate is the highest.Another study investigated the effects of the ejector geometry and the operating conditions on the hydraulic characteristics in a rectangular bubble column with a horizontal flow ejector [7].It is found that the gas suction rate increases with the liquid circulating rate in the multiphase CFD simulation with the mixture model and the experiments.However, the gas suction rate decreases as the liquid level increases in the column and nozzle diameter.However, there is no research reporting the use of a liquid-gas ejector in ballast water treatment system, relevant numerical, or experimental research.The aim of this work is to introduce a new type of the ballast water treatment system based on a liquid-gas ejector and then utilizes both experimental and computational fluid

Experimental Setup
An experimental facility for the real-world application is shown in Figure 3.In the experiments, air is used in place of ozone as the secondary fluid; water is employed in place of ballast water as the primary fluid, and a check valve is installed on suction port to prevent back-flow, for demonstration purposes.The whole system is composed of the main pipe, the pump, and the ejector.The ejector is horizontally equipped with the main pipe, which maintains a constant pressure and allows water to circulate.The pressure of the fluid circulating in the main pipe is monitored to ensure a constant back pressure.A small amount of fluid circulating in the main pipe enters the ejector nozzle by a pump; the flow rate and the pressure of the fluid coming out of pump can be measured using a water flow meter and a pressure gauge.As the primary fluid is pumped through the ejector nozzle, the pressure of the secondary fluid entering through suction chamber ( suction ) can be measured using the pressure gauge.Another pressure gauge measures the pressure of the waterair mixture fluid at the elbow..The principle of operation of this ejector is straightforward.It begins with a high velocity liquid from the pump and accelerates the liquid at the throat where it produces a low pressure zone that draws in and entrains a suction fluid.After passing through the throat, a mixed fluid expands and the velocity is reduced which results in recompressing the mixed fluids by converting the kinetics energy back into the pressure energy.Therefore, the liquidgas multiphase flow is taken into account in this study.The governing equations for the two-phase flow are then given as follows.

Governing Equation
The continuity equation for each phase is The momentum conservation equation for each phase is where   ,   , and   are the void fraction, density, and viscosity of the phase , respectively.  represents the sum of interfacial forces including the drag force   , lift force   , wall lubrication force  WL , turbulent dispersion force  TD , and virtual mass force  VM .The virtual mass force has been neglected in this work, since its influence is insignificant as compared with the other drag and nondrag forces for the steady-state multiphase flows [14].  denotes the phase generation rate between different phases.By assuming the continuity and momentum equations are identical for all phases, the mixture momentum equation can be defined as where the mixture density   is defined as In this work, the number of phases  = 2. Since the sum of all phases must occupy the whole domain volume, the following constraint must be satisfied: Similarly, the mass sources between phases must conserve: (6)

Turbulent Closures.
The Reynolds number inside of the ejector is much greater than 10 5 , which shows that the flow through the ejector is turbulent.To deal with the turbulence flow, the Reynolds-averaged Navier-Stokes equations (RANS) are utilized.In this case, simulation based on RANS equations substantially reduces the computational effort in comparison with the direct numerical simulation (DNS), but the averaging procedure introduces additional unknown terms containing products of the fluctuating components.This acts like additional stresses in the fluid, so that a turbulence model is needed to close the set of equations.
In this work four two-equation turbulence models named standard  − , RNG  − ,  − , and shear stress transport (SST) model are adopted to fulfill the aim.It should be noticed that the SST model is based on the combination of standard − and − turbulence models in the formulation of Wilcox.The hybrid procedure consists of the -equation and a special form of the -equation, which enables changing the value of a blend factor  1 switching between a -equation ( 1 = 1) and a -equation ( 1 = 0) [15].

Computational Domain and Boundary Conditions.
Figure 4 shows the grid model in ANSYS CFX 14.5.For accurate simulations, the CFD model is created at a 1 : 1 scale with the exact geometry as shown in Figure 3.The continuous phase from the nozzle section (inlet 1) utilizes water at 25 ∘ C, and the dispersed phase from the suction section utilizes air at 25 ∘ C. Since the volumetric flow rate of water is set and known as constant in the experiments, a fixed flow rate of 7.009 kg/s is used as a boundary condition at inlet 1.
The inlet condition with atmospheric pressure is applied at inlet 2, preventing air and/or water from flowing out of the ejector which coincides with the physical experiment because a check valve is installed inside of the suction pipe.Inlet 3 is defined as an inlet boundary condition with a constant flow rate.The main outlet is set as the opening condition with various back pressures.In addition to the inlet and outlet conditions, all other surfaces are defined as walls under no-slip boundary conditions.The simulation is initialized with full air; that is,  air = 1.With regard to the solver algorithm, high resolution is used for the advection scheme and turbulence numerics.The convergence criteria for continuity, momentum, and transport equations are set to 10 −4 , as multiphase simulation is numerically difficult to converge.

Mesh and Bubble Independence.
The mesh size is a virtual issue for CFD simulation because it not only significantly influences the computational accuracy but also directly determines the computational effort.Therefore a mesh independence study has to be carried out to find out the desirable mesh size, which can represent a good compromise between computational effort and accuracy.In this work, a total of five types of mesh sizes are used, namely, coarse, medium, fine, very fine, and extremely fine.Table 1 presents a comparison of these mesh sizes with different nodes and elements.Their effect on the elbow pressure and the mass flow of suction air is illustrated in Figure 5.It can be seen from this figure that the mesh size has little effect on the CFD simulation when it increases to fine mesh size and above.Therefore, the very fine mesh size is determined as    the most appropriate size striking a balance between the computational accuracy and the computational effort.As a large amount of discrete air bubbles are involved in this simulation, a bubble size independence study is carried out so as to identify an appropriate bubble size.Figure 6 shows the mass flow of suction air and the elbow pressure with the change of bubble sizes.It is easily concluded that bubble size has negligible effect on the results.Therefore, a classical bubble diameter of 1.0 mm is defined in this work, which is also in agreement with a previous study [9].

Results
The effect of the back pressure on the suction pressure and elbow pressure is investigated in this work, and four turbulence models are also compared with the experiments to find the most appropriate one for the targeted application.
As shown in Figures 7 and 8, the deviation between the CFD results and experimental results is all below 8% for the prediction of the elbow pressure and below 10% for the prediction of the suction pressure.The deviations could be attributed to the following modeling issues: (i) the application of impressible air at 25 ∘ C instead of the real compressible air and the exclusion of heat transfer between the different phases and (ii) the inadequacy for the mixture model and interfacial force model produced to handle the large amount of bubbles.On the whole, the RNG  −  turbulence model is in the best agreement with experiments.Its mean that error is 3.25% for elbow pressure prediction and 3.36% for suction pressure prediction.Next, the SST and  −  models show less accurate results than the RNG  −  model.The standard  −  turbulence model failed to converge for any meaningful results and is thus inappropriate for the CFD analysis of the liquid-gas ejector.Therefore, the RNG  −  model is  chosen for the multiphase ejector simulation and is used for the subsequent analysis for understanding the characteristics of the ejector.In addition, it should be noticed that the residue for each equation falls down below 10 −5 when the back pressure is higher than 2.86 × 10 5 Pa (even though the default value is 10 −4 ).This is due to the fact that the two-phase liquidgas problem becomes a single-phase liquid problem under a high back-pressure condition, where a wall has been placed at portion(s) of the inlet 2 to prevent fluid from flowing out of the domain.
For a pressure comparison, as the back pressure increases from 1.15 × 10 5 Pa to 2.86 × 10 5 Pa, the pressure ahead of the elbow increases almost linearly from 1.9 × 10 5 Pa to 3.8 × 10 5 Pa.The suction pressure shows a similar trend, increasing from −7049 Pa to 92025 Pa.This suggests that the flow rate of absorbed air decreases from 48.0 m 3 /hr to 0.0 m 3 /hr as shown in Figure 9.That is, no air is absorbed into the ejector when the back pressure is greater than 2.8 × 10 5 Pa.In essence, the kinetic energy in the high velocity water jet is mainly transferred into the potential energy (high back pressure), leaving no extra energy to absorb the air at the ejector.
Figure 10 shows the contour of the air volume fraction at three different back pressures.It can be clearly found that as water leaves the nozzle tip, it forms a strong jet that passes through the suction chamber and mixing chamber which produce a lower pressure region to absorb the air from the suction port.A clear phase interface between water and air in the suction chamber means they are not well mixed.Air flows in the annular region between the jet of water and the ejector wall (i.e., coaxial-flow due to the high velocity water jet).As water and air flow in the same path and annularly    downstream, their mixing becomes evident and intensive, which results in a large amount of air bubbles dispersed in the water.The jet intensity decreases dramatically as the back pressure increases because it needs to be transferred to pressure energy to withstand the effect of the back pressure.In turn, the decrease in the jet momentum leads to a pressure drop between the suction chamber and the ambient air.Therefore the flow rate of absorbed air drops sharply.This ejector-based system will be used under a high back pressure condition, so its performance under high back pressure condition is of great concern.Optimization of ejector in the previous work is mainly focusing on the suction chamber and mixing chamber; this usually costs a lot of effort and time as there are too many dimensions.In this work, a quick method which only alters the nozzle shape is proposed to efficiently achieve the goal of optimization.Figure 11 shows the initial nozzle and the newly designed nozzle.The new one has a straight shape and its diameter is the same as the throat diameter of the initial curved shape.The CFD simulation of the new design has been conducted with the same mesh size, bubble size, and boundary conditions as depicted in Section 3.3.
The CFD results are shown in Figure 12.It can be found that at low back pressure condition (  ≤ 1.7 × 10 5 Pa), the suction performance is a little worse than the initial nozzle.However, when they are working under higher back pressures, the advantage of the new straight nozzle is becoming more and more evident.Specifically, at the back pressure of 2 × 10 5 Pa, the absorbed air by the initial design and new design are 20.95Kg/hr and 23.44 Kg/hr, respectively, with an improvement of 11.89%.At the back pressure of 2.5 × 10 5 Pa, the absorbed air by the initial design and new design are 6.915 Kg/hr and 12.99 Kg/hr, respectively, which means that an improvement of 87.85% has been achieved.The initial one will loss the capability of absorbing air at back pressure of about 2.8 × 10 5 Pa, but the new design still can absorb approximately 8 Kg/hr air from the suction port at the same back pressure.The new design will stop absorbing air at the back pressure of 3.13 × 10 5 Pa eventually.This improvement can be attributed to two reasons: (1) the straight nozzle generates a less pressure loss than the curved shape, so a higher motivate energy is maintained through the nozzle; (2) a portion of the drive water diverge to the suction chamber from the nozzle due to the curved shape and block the air flow to mixing chamber, which does not happen in the straight nozzle design.

Conclusion
Shipping ballast water treatment is a critical process to minimize the impact of ballast water discharge on the marine environment.The paper has presented a novel ballast water treatment system with a liquid-gas two-phase ejector.Numerical simulation and experimental tests are conducted to understand the performance of the ejector.Pressures of the ejector were first examined on an experimental bench to validate CFD models with different turbulence models, followed by extensive CFD analysis.Test results on the liquidgas ejector have clearly shown that RNG  −  turbulence model with the very fine grid size is most effective striking a balance between accuracy and computational costs.CFD analysis also indicates that the initial liquid-gas ejector may lose its absorption capability when the back pressure exceeds 2.8 bars in the ballast water treatment system.With the straight nozzle altered, the new design can efficiently improve the ejector performance at high back pressures and will not lose the absorption capability until 3.13 bar back pressure.In the work, in-depth study on the ejector nozzle and other geometrical parameters will be carried out to find the optimal geometry for the ballast water treatment.The multiphase reaction flow will also be taken into account for understanding the mixing characteristics and the treat efficiency of ballast water.

Figure 1 :
Figure 1: Schematic of the proposed ejected ozone ballast water treatment system.

Figure 2 :
Figure 2: The liquid-gas ejector for the ballast water treatment.

Figure 4 :
Figure 4: CFD model of ejector and pipe.

Figure 5 :
Figure 5: Mesh effect on the results.

Figure 6 :
Figure 6: Bubble effect on the results.

5 CFD
results with RNG k- model CFD results with k- model Fitted curve of RNG k- results Fitted curve of RNG k- results

Figure 7 :
Figure 7: Comparison of  elbow with back pressures.
Back pressure (Pa) Suction pressure (Pa) CFD results with SST model Fitted curve of SST results Ambient pressure (1 atm) Fitted curve of experimental results ×10 5 CFD results with RNG k- model CFD results with k- model Fitted curve of RNG k- results Fitted curve of k- results

Figure 8 :
Figure 8: Comparison of  suction with back pressures.
Flow rate of air (kg/hr) CFD results with SST model Fitted curve of SST results ×10 5 CFD results with RNG k- model CFD results with k- model Fitted curve of RNG k- results Fitted curve of k- results

Figure 9 :
Figure 9: Absorbed air flow rate with back pressures.

Figure 10 :
Figure 10: Volume fraction of air with different back pressures at the middle plane.

Figure 11 :
Figure 11: Initial nozzle and new nozzle.
Flow rate of air (kg/hr) CFD results of initial nozzle Fitted curve of initial nozzle CFD results of new nozzle Fitted curve of new nozzle ×10 5

Figure 12 :
Figure 12: Comparison of air mass flow for different nozzles.

Table 1 :
Five types of mesh sizes and their nodes/elements.