Effluent Discharges from Two Outfalls on a Sloping Beach

A marine outfall is a long pipeline that continuously discharges large amounts of effluent streams into the sea. As the number of marine outfalls along the coastal areas is growing, a far field mathematical model with two point sources on a sloping beach is used to assess the coastal water quality following discharges from two outfalls. Asymptotic approximation will be made to the concentration at the beach to measure how well the effluent plumes are mixed and diluted in the coastal waters. The result found agrees with the engineering practice of installing a two-port diffuser at the end of a single outfall to minimize its potential environment impacts.


Introduction
Most coastal industrial installations and plants, such as municipal sewage treatment plants [1] [2], power generation stations [3], and seawater desalination plants [4]- [8], dispose of their wastewater effluents through long outfall pipes that stretch far into the ocean.For a modern plant, a multiport diffuser would also be installed at the pipe-end to rapidly dilute the effluent stream.Because of relatively shallow coastal waters, it is observed that the elongated effluent plumes are spreading towards the shoreline and may cause a concentration build-up [9]- [11].Due to the uncertainty in sea conditions, a clear understanding of the mixing processes of effluent plumes is not yet known, and the use of mathematical models has been a key strategy for the basis of sound engineering design and for assessing the potential environmental impacts of marine outfall effluent discharges [2] [3] [5]- [8].
When many marine outfalls are discharging effluents to shallow coastal waters, the adverse long-term impacts are strongly inter-dependent, and compounded from neighbouring outfalls.The interactions of two or more effluent plumes are expected as many outfalls often tend to be closely clustered together along the coastal areas.Newly-constructed coastal plants may need to build two outfalls as a contingency plan for a future increase in the plant's production capacity [6].For two outfalls discharging a given integrated total effluent stream, the waste load can be allocated optimally between them to minimise the impact [12] [13].
As coastal industrial plants are built predominantly on the sloping sandy beaches, a mathematical model using a two-dimensional advection diffusion equation with two point sources is presented.The solution is plotted to graphically study the merging of two effluent plumes from two outfalls.While the far field modelling in this paper involves drastic simplifications, key physical mixing and dispersion processes are represented, and thus the analytical solution remains useful in providing a qualitative understanding and in suggesting general behaviour of the marine outfall effluent discharge plumes in coastal environment [9] [11] [12].

Mathematical Analysis
The beach is considered to be straight and the sea wide, and the outfall's effluent plume is assumed to be vertically well-mixed over the water depth.The coastal (drift) current is assumed to be steady with a speed U and remains in the x -direction parallel to the beach at all times.The dispersion mechanisms are represented by eddy diffusivities, and diffusion in the x -direction is neglected, as the effluent plumes in steady currents become very elongated in the x -direction.The variations in the y -direction of U and coefficient of dispersivity D are assumed as the power functions only of water depth h , and for application, we take U to be proportional to 1 2   0 h and D to 3 2 0 h , where 0 h is an arbitrary reference water depth.These scalings are appropriate for a turbulent shallow water flow over a smooth bed [9]-[11] [14].For simplicity, other complexities such as tidal motions, density and temperature are ignored.
As illustrated in Figure 1, we represent the old outfall as a point source at the position ( ) discharging an effluent stream at a constant rate 0 Q , where α is the source length.Similarly, the new outfall as a point source at ( ) ( ) Q , where ε is the outfall's (offshore) and  (along the shore) separation distances.Without loss of generality, we also assume that if these two out- falls are operated by one plant discharging, then a combined effluent total rate 0 1 As the water depth is gradually decreasing towards the beach at 0 y = , on a uniformly sloping beach with slope m , we formulate ( ) h y my = . Following [9] [11] and applying a linear superposition, the two-dimensional advection-diffusion equation for the far field plume concentration ( ) , c x y from the two point sources is given by with boundary condition 0 hD c y ∂ ∂ = at the beach 0 y = , and the effluent concentration is assumed to be ultimately dissolved far into the sea, 0 c → as y → ∞ .δ is the Dirac delta function.In terms of dimension- less quantities * 0 where Q denotes a reference discharge rate which usually adopts the value of the original discharge rate of the old outfall.By setting A. Purnama , the analytical solution of Equation ( 1) is given by where 3 2 I is a modified Bessel function [15] [16].The model parameter λ represents the effluent plume elongation in the x -direction [9]- [11].In coastal waters, larger values of λ are mostly due to a stronger cur- rent 0 U with less dispersivity 0 D .For the quantitative illustration of the solutions, the values of 0.02 m = and 5 2 λ = will be used in all plots.The other parameters are related to the position of the point sources: α the single point source length, and ε and  the point source's separation distances.
As a higher build-up in concentration is more likely found at the shallow water close to the beach [7] [8] [10] [14], the appropriate measure for assessing the impact of marine effluent discharges from sea outfalls would be the concentration values at the beach.In the limit as * 0 y → and replacing 3 2 I in Equation ( 2) by its asymptotic form [16], we obtain the compounded concentration at the beach ( ) ( ) It is easy to see for effluent discharges from a single point source at ( ) and since 0 1 q = (and 1 0 q = ), the concentration at the beach reduces to ( ) .
The concentration at the beach for a point source length 3 α = and 4 α = is plotted in Figure 2. By diffe- rentiating, the concentration has a maximum value of α = , and this maximum value is reduced by more than 50% to 0.95 when the length is extended to 4 α = .This result agrees with the standard practice of building a longer sea outfall in order to mi- nimize its potential environmental impact in the coastal waters [1] [2].

Two Independent Outfalls
Apart from the effluent discharge rates, the compounded impact of the new outfall is governed by the outfall's separation distances ε and  , and in particular, if the value of  is large (e.g.
), the two outfalls are well separated, and thus it is expected that the contribution of the new outfall located at ( 1 0 , . We first consider the case where the two outfalls are operated independently by two coastal plants, i.e., when 0 1 1 q q + ≠ and when both separation distances 0 >  and 0 ε > .The concentration contour plots of Equation (2), the solution for two point sources when 3 α = and 2 ε = , is shown in Figure 3, which illustrates the merging of plumes from the two point sources for two values of the separation distance 3 =  with 0 1 q = and 1 0.5 q = , and 8 =  with 0 0.5 q = and 1 1 q = , where two separate plumes are clearly shown., where the new point source is discharging with a rate 1 1 q = , double that of the old point source.Next, the compounded concentration at the beach for two point sources as given by Equation ( 3) can be simplified to  , 0 , 0 exp Figure 4 shows the concentration at the beach for the two cases illustrated in Figure 3, where two separated plumes are depicted as two distinctive peaks.For comparison, the concentration for the single point source with 0 1 q = is also shown in Figure 4 by the dotted line.It is worth noting that for * 0 x > , the presence of the new outfall does not change the location of the maximum concentration at the beach [11].
Substituting * 2 5 m x λα = in Equation ( 4), the maximum value of concentration at the beach for two point sources is approximated as ( )

( )
f z has a maximum value of ( ) , which occurs at m z ε α = , and therefore Note that for a given value of m z ε α = , the value of α  can be determined by ( )( ) As shown in Figure 5, the largest maximum value of compounded concentration at the beach for two point sources is always greater than that of the concentration at the beach for a single point source.The contribution of the new outfall is much greater for 1 ε α < .For example, if the new point source length is double that of the old point source length, i.e. 1 ε α = (and α ε α ), then from Equation ( 5), the maximum value of the concentration at the beach is [ ] . If both point sources are discharging at A equal rates, i.e. 1 0 1 q q = , then the maximum value of 1m c is about 18% higher than 0 0m q c , the maximum value of the single point source.However, if 1 0 1.5 q q = , the maximum value of 1m c is increased to about 27% higher than 0 0m q c .

Single Outfall with Two Ports
Modifying the rates of discharge can be achieved through a relatively less expensive method than extending the A. Purnama outfall lengths.For the case where two outfalls is operated by one plant, i.e. when 0 1 1 q q + = , the effluent dis- charge allocation can be arranged between these two outfalls [12].Binomial expansion of Equation ( 5) for 1 ε α < , gives 5 2 That is, the maximum value of 1m c is smaller than that of the single point source value 0m c .This result agrees with the previous finding [12] [13], that if two outfalls are discharging a given integrated total effluent stream, the waste load can be allocated optimally between them to minimise the impact.However, economically it is cheaper to build an outfall and install a two-port diffuser at its pipe-end than build one more sea outfall.
To include the case of the modern engineering practice that installs two ports at the end of a marine outfall [5] [6], we also assume that both port separation distances ε α and α  are small, and thus the value of ( ) can be approximated for 0.05 z ≤ by ( ) . For plotting the contours of the effluent plume, since the point sources are close to each other, and in the limit as * 0 x → and replacing 3 2 I in Equation ( 2) by its asymptotic form [16], we obtain The contours of the solution for a single outfall with two ports are reproduced graphically in Figure 6 for 3 α = with 0 1 0.5 q q = = , when 0.05 z = , which is equivalent to the port separation distance 0.05 α =  for 5 2 λ = .The effluent plumes from these two closely located point sources are immediately merged as they are released, and for * 0 x > , and the combined plumes appears to be spreading as one.This supports the concept that a two-port diffuser will rapidly dilute effluent streams.
Finally, the maximum value of compounded concentration at the beach for two closely located point sources can be approximated to Again, the maximum value of 1m c is smaller than that of the single point source value 0m c .From Figure 7, we noted that as the value of 1 q increases and the port separation distance ε α gets longer, the maximum value becomes smaller than 0m c .For example, if 0.1 ε α = and both ports are discharging at equal rates, i.e. 1 0 0.5 q q = = , then the maximum value of 1m c is about 10% less than 0m c .This result agrees with the modern engineering practice that installing a two-port diffuser at the end of a marine outfall will improve the mixing and dilution of effluent discharge plumes in coastal waters [5] [6] [12].

Conclusions
The solutions for an advection diffusion equation with two point sources are applied to study the interaction and merging of effluent discharge plumes from two outfalls on a sloping beach.As a measure for assessing the impact in the coastal environment, the maximum compounded concentration at the beach is formulated.If the two outfalls are independently operated, then the maximum value of the concentration at the beach can be minimized as long as the new outfall length is more than double the old outfall length, and discharging at a rate smaller than the old outfall.
If two outfalls are operated by one plant, then the integrated total effluent load can be shared between them, and it is found that the maximum value of the concentration at the beach is smaller than that of the single outfall.A similar result is also obtained for a single outfall where a two-port diffuser is installed at the end outfall pipe.However, implementation issues related to the control of discharge rates, reliability and cost effectiveness of the marine outfall are not addressed.
The mathematical formulation presented can be extended to the case of multiple outfalls discharging on a coastal area, and in particular, the result for effluent discharges from multiport diffusers has been reported in [11].

.
This maximum value is inversely proportional to the point source length α[9] [10].A value of 1.96 is obtained for 3

and the top part of Figure 3 Figure 3
represents a situation where two outfalls are relatively close to each other at a small distance apart, represents two outfalls at a slightly longer distance apart,

Figure 2 .
Figure 2. The concentration at the beach for a single point source.

Figure 3 .
Figure 3. Merging contour of two point sources plumes when 3 α = and

Figure 4 .
Figure 4. Compounded concentration at the beach for two point sources.

Figure 5 .
Figure 5. Maximum value of 1 0 0 m m c q c .

Figure 7 .
Figure 7. Maximum value of 1 0 m m c c .