Numerical examination of concentration-dependent wastewater sludge ejected into a drinking water source

One of the significant water-related health challenges globally is due to pollutant fate. Contaminants endanger the lives of humans, animals, and even plants. The present mathematical analysis explains reactive wastewater sludge ejected into a drinking water source from wastewater treatment plants. The assumption that wastewater sludge follows a power-law constitutive relation leads to nonlinear momentum and concentration equations. The contaminants are assumed to follow a nonlinear irreversible first-order sorption model. The numerical solution of the coupled problem is solved using the Bivariate Spectral Local Linearization Method and validated with the spectral Chebyshev weighted residual method. Profiles are presented for dimensionless flow velocity and concentration. Comprehensive explanations for the obtained results are provided with relevant applications.

obtained in [5][6][7][8][9] will break down since the shear stress-strain relation cannot be linear in the non-Newtonian fluid.As a result of this behavior, Santos et al. 11 presented a power law rheological model to model sludge behavior from a water treatment plant.Csizmadia and Till 12 , in their study of a comparative analysis of the power law by Bingham and Hersschel-Bulkley, presented non-Newtonian models for activated sludge.A similar study by Haldebwang et al. 13 conducted a curve fitting for three different rheological models for both laminar and turbulent cases of sludge in a pipe flow.More exciting results on the rheological properties of wastewater sludge can be found in works reported by Eshtiaghi et al. 14 and Wei et al. 15 and some references cited therein.
Motivated by the studies on the non-Newtonian flow documented in [11][12][13][14] , reaction kinetics, and the channel porosity effect 10 , the focus here is to investigate the combined impacts on the work done in 5 .The outcome is expected to be useful in treating wastewater sludge from most treatment plants and other related industrial and manufacturing processes where sludge is constantly ejected into drinking water sources.Understanding these flow properties will provide more information to ensure the United Nations' drive to meet goal #6-the provision of clean water.A thorough literature survey shows that the study is yet to be reported.Therefore, this study on the non-Newtonian sludge flow using the power law constitutive model is innovative.In the next section, the mathematical formulations for the unsteady fully developed flow with constant vertical penetration are presented for both flow momentum and concentration of pollutants.This way, a set of coupled nonlinear partial differential equations will be presented in Sect."Mathematical analysis" of this paper with detailed mathematical analysis, while Sect."Bivariate spectral local linearization method of solution" is devoted to constructing a bivariate spectral local quasi-linearization method for the approximation to the initial-boundary value system of partial differential equations.

Mathematical analysis
The transient flow problem earlier studied by Chinyoka and Makinde 5 for nonlinear contaminant ejection is revisited.The buoyancy-induced flow along a vertical channel whose x-axis is in flow direction and the y-axis is taken in a perpendicular direction.The flow assumptions are now as follows: (i) the fluid is dilatant, reactive, and incompressible, (ii) fully developed condition implies that ∂ ∂x = 0 , (iii) the leaky channel walls allow constant fluid penetration by suction/injection, (iv) at a time t = 0 with µ = µ 0 and D = D 0 , (v) at t > 0, the fluid properties become,µ = µ(C), D = D(C), S = Q(C) (vi) the nonlinear buoyancy dependence on the concentration of the contaminant is neglected.(vii) slippage is negligible.
Where µ, µ 0 represents the dynamic and constant component of water viscosity C, D, represents the dissolved pollutant concentration and diffusion coefficient, D 0 , S are the referenced diffusion and nonlinear source term components.Due to possibility of leakages of toxic chemical like flocculants from the water treatment plant the model of Chinyoka and Makinde 5 is modified to include the effect of suction and injection.In this regard, the unsteady coupled governing equations the flow of waste water can be written as: Additional terms in (1) arise due to departure from the Newtonian to non-Newtonian flow behaviour of sludge when m = 1 .Since sludges are mostly emulsions due to particle aggregation, the irreversible coalescence process is represented by the first-order reaction rate ∂C SP ∂t .In Eqs.(1) above, P, ρ, u, v 0 , t, x, y, Ŵ, m represent the pressure, sludge density, velocity, constant injection/suction velocity, time, axial, vertical coordinates, pores, and power-law index.g, β are the gravitational acceleration, and coefficient of expansion.The non-symmetrical nature due to wall leakages suggests the initial, no-slip, and non-moving wall boundary conditions.M, a, C w , C 0 are initial condition control parameters, channel half width, and referenced wall concentrations.Following Serato 16 , the solid reactive contaminant C SP in the sludge follows the irreversible first-order sorption model of the form The dependence of dynamic viscosity, diffusion, and nonlinear source on pollutant concentration are defined as (kindly refer to Chinyoka and Makinde 5 ): (1) 2,3) denote the constant component of viscosity and variation coefficients.First-order rate constant, reactive solid contaminant, and external source coefficients respectively.
with the following variables and parameters, For dilatant fluids with positive values of the index m, we get, with initial and boundary conditions, In ( 5), (u, φ, m) , dimensionless velocity and pollutant concentration, power-law index, (R, K, α, Gc) represents suction Reynolds number, constant pressure gradient, viscosity variation parameter and solutal Grashof number, (Sc, γ , ), Schmidt number, concentration dependent parameter, nonlinear source parameter, (β, n, k 1 ) pollutant variation parameter, order of chemical reaction and coefficient of chemical kinetics.

Bivariate spectral local linearization method of solution
In the following section, the procedure for the computable form of the system of Eqs. ( 6), ( 7) will be presented based on bivariate spectral local linearisation approximation theory.Following Bellman and Kalaba 17 , including [18][19][20] we set Such that the quasi-linearized version of (6) becomes where the coefficients used in (9) are defined as It is well known that set of polynomials is dense in the set of continuous functions, therefore, to obtain the solution of the (6) subject to (7), we seek a series solution that is based on Lagrange cardinal polynomial, L p y L q (τ ) , approximation of the form: (5) a 0,r (y, t) ∂ 2 u r+1 ∂y 2 + a 1,r (y, r) At every Chebyshev-Gaus-Lobatto points,(y i , τ j ) , derivatives appearing in ( 9) can be computed easily based on the following Chebyshev differentiation matrix procedure: In (14), D and d j,q = T 2 d j,q represent the Chebyshev differentiation matrices for (N y + 1) × (N y + 1) and (N τ + 1) × (N τ + 1) respectively while U j and j are given by: with superscript T and I denoting the transpose and identity matrices, respectively.Substituting (15) and ( 14) into (9) we get The transformed boundary conditions are: The vector U r+1,N τ and r+1,N τ correspond to the initial condition given in Eq. (7).Matrices (15) are solved iteratively until suitable results are obtained.The accuracy of the computation is validated by using the regular bivariate Chebyshev collocation method and the results are presented as Tables 1 and 2 in the results and discussions section.

Results and discussion
In this section, numerical results of the dimensionless problems (6) subject to the initial and boundary conditions (7) are presented graphically and in tabular form to describe the effects of various flow parameters on the concentration and flow fields.The numerical results are computed using the quasi-linearized approach based on the bivariate Lagrange polynomials as the trial function and validated with the spectral Chebyshev collocation (11)   u(y, t) ≈ U(y, t) (15) U j = [U(y 0 , τ j ), U(y 1 , τ j ), U(y 2 , τ j ), ..., U(y N y , τ j )] T , for j = 0, 1, 2, ..., N τ , j = [�(y 0 , τ j ), �(y 1 , τ j ), �(y 2 , τ j ), ..., �(y N y , τ j )] T , for j = 0, 1, 2, ..., N τ , (16) a 0,r (y, τ j )D 2 U r+1,j + a 1,r (y, τ j )DU r+1,j + a 2,r (y, τ j ) www.nature.com/scientificreports/method.At this point, the quasi-linear bivariate Langrage method converges faster than SCCM, and the CPU time is more minor even with the power law index greater than unity.In this regard, we have taken a near unity value for the comparative study of the two numerical schemes.As shown in Tables 1 and 2, there is a high level of agreement between the two results.Therefore, the numerical solutions are unique to problems ( 6) and ( 7).
In the obtaining Figs. 1, 2, 3, 4, 5, 6 and 7, we have used the following parameter values R = 1, m = 1.1, Gr = 0.1, α = 0.1, Sc = 0.6, γ = 0.1, = 0.5, β = 0.1, = 0.1, n = 2 .Figure 1 represents the resid- ual error plots when the numerical approximations are returned as solutions to the problem ( 6), (7).The residual errors are of order 10 -9 and 10 -10 for the dimensionless momentum and pollutant concentration, respectively.As a result, the numerical results give excellent residuals of R 1,r y, t and R 2,r y, t in ( 9) and the residual errors converges to 10 -10 and 10 -11 respectively as shown in Fig. 2. Figure 3 reveals the steady state behavior of the numerical solutions.Evidently, the solution approached steady state from t = 2 .However, as t = 5 and above, further solutions are superimposed on t = 2 .Moreover, as seen in Fig. 4, the time-independent solution for pollutant concentration is recovered from t = 2 while superimpositions are seen afterward beyond this point.Figure 5 reveals the non-global nature of the solution that is typical of all nonlinear PDEs with exponential nonlinearity.The solution is seen to blow up with the nonlinear source term parameter.The 3-D solutions of the initial boundary value problem ( 6)-( 7) are presented in Figs. 6 and 7. Evidently, the solutions are well-behaved as the initial and boundary conditions are well satisfied.Figure 8 shows the impact of buoyance on the flow velocity.As seen from the plot, the flow velocity peak increases with increasing solutal Grashof number values.This results from increasing buoyancy force over viscous force around the core area of the flow channel.Figure 9 depicts the influence of nonlinear sorption kinetic on the flow velocity.The result shows that the wastewater velocity is at its peak without sorption.Further increase in the nonlinear sorption coefficient decreases the flow velocity due to adsorption at the porous walls.In Fig. 10, pollutant concentration is observed to be at maximum when k1, the nonlinear sorption coefficient is absent from the pollutant concentration.Further increase in the parameter resulted in decreased concentration due to the filtration process at the leaking walls.The influence of the nonlinear source term on the flow velocity is presented in Fig. 11.The result shows that the exponential dependence on pollutant concentration enhances flow velocity.This is physically correct due to the convective buoyancy effect on the flow velocity.However, a reverse trend is seen in the spatial variation of pollutant concentration with nonlinear source parameters, as seen in Fig. 12.Meanwhile, Figs. 13 and 14 represent suction/injection Reynold's number due to wall porosity.These plots show that the flow retains symmetry if the walls are impervious.Otherwise, the symmetry is lost as the pore opens and becomes wider, and more solid particles in the wastewater are introduced into the channel.Therefore, the flow pattern is distorted through vertical penetration.In bed filtration, Reynold's number is the main application being used since polluted water is injected into the fixed bed where solid particles are first removed when the pores are more prominent, then filtered water is made to pass through several layers with varying porosity until clean water is collected at the suction end.Lastly, Figs. 15 and 16 show the effect of mass diffusivity and molecular diffusion on the flow and concentration profile.From the plot, pollutant diffusion within the flow channel decreases as the Schmidt number increases.As a result, the concentration of pollutants rises around the deposited point within the channel.Consequently, mass diffusivity increases, and the fluid behaves like a viscoelastic fluid with increasing velocity.

Conclusion
In this work, a numerical investigation into the unsteady flow of wastewater sludges into drinking water through a porous channel has been carried out.The pollutant is assumed to be chemically reactive and follows the irreversible nonlinear first-order sorption kinetics.The nonlinear governing partial differential equations are formulated based on the fluid dynamics approach and turned into dimensionless problems.The set of nonlinear initial-boundary value problems is solved by the bivariate quasi-linearized model and validated with the spectral Chebyshev collocation method.The significant contributions to knowledge from the present computations are as follows: effect of suction/injection Reynold's number in flow control is more substantial and cannot be neglected.Irreversible chemical kinetics is significant in the understanding and treatment of wastewater.The influence of solid particles is more pronounced in non-Newtonian flow behaviour wastewater than in Newtonian flow.
In our subsequent formulation, the gravity-driven flow with bioremediation and biochemical oxygen demand effects in contaminated river water flows will be considered under the Monod reaction.

Figure 3 .
Figure 3. Steady state solution for u t, y .

Figure 6 .
Figure 6.Solution for u t, y in 3-D.

Figure 8 .
Figure 8. Impact of solutal Grashof number on u t, y .

Figure 9 .
Figure 9. Impact of nonlinear sorption kinetic parameter on u t, y .

Figure 10 .
Figure 10.Impact of nonlinear sorption kinetic parameter on φ t, y .

Figure 11 .
Figure 11.Impact of nonlinear source parameter on u t, y .

Figure 12 .
Figure 12.Impact of nonlinear source parameter on φ t, y .

Figure 13 .
Figure 13.Impact of suction Reynold's number parameter on u t, y .

Figure 15 .
Figure 15.Impact of Schmidt number on u t, y .

Figure 16 .
Figure 16.Impact of Schmidt number on φ t, y .

Table 2 .
Validation of Bivariate spectral local linearization scheme for dissolved concentration when