Analysis of the Possibilities of Energy Recovery from Gravity Flows in Pipelines in a Copper Ore Enrichment Plant

: Considering the need to reduce climate change and energy costs, the possibilities of energy recovery have been sought for many years, even if this should involve signi ﬁ cant complexity of the technological process. A tt empts to recover energy from a ﬂ owing medium containing a large amount of a solid phase are rather rare. This is mainly due to unfavorable erosion and sludge phenomena. The aim of this work is to analyze the possibility of recovering energy from gravity ﬂ ows in the pipelines of a copper ore enrichment plant. The analysis will be carried out for a system operating in a plant in Poland. An energy recovery device will be designed taking into account the ﬂ ow of a two-phase medium. For this purpose, ﬂ ow calculations will be performed. The calculation results obtained for the two-phase ﬂ ow will be veri ﬁ ed based on the experimental results presented in the literature. The simulations will make it possible to determine the power of the designed device.


Introduction
Considering the need to reduce climate change and energy costs, the possibilities of energy recovery have been sought for many years, even if this should involve significant complexity in its technological process.The application of a simple system preheating the mine intake air in arctic and subarctic regions is proposed in [1].The idea is to use waste heat from the exhaust gases of diesel electricity generators.Most often, however, there is a need to cool the mine, which results from the auto-compression of mine air and heat emissions from ore beds and mine machinery.The benefits of reutilizing the water of the bulk air-coolers for mine air pre-conditioning during the wintertime are studied in [2].The application of dynamic empirical models for mine ventilation purposes brings many benefits [3].A gravity heat pipe for heat recovery from return air in a mine is analyzed in [4].The paper presents theoretical guidelines and improvement solutions for practical engineering applications, which are formulated based on the presented numerical calculations.A medium-size excavator is modified and an energy recovery system is introduced in [5].The proposed system was designed to recover the potential energy of the boom, using a hydraulic accumulator as a storage device.The numerical results show that the designed energy recovery system is able to provide about 9% of mechanical energy.In [6], a new method of magnesium extraction is proposed.Compared to the traditional method, the Mg leaching efficiency increased by 36.42%, the washing time was significantly shortened and the energy consumption was reduced.Energy can also be recovered from waste using processes such as combustion, gasification and pyrolysis [7,8].Hydrogen is one of the most important green energy sources, although it is difficult to store.The recent developments in the area of photoreactor development are presented in [9].Cu-based catalysts for H2 production are characterized by various techniques in [10].Greater energy density in portable proton exchange membrane fuel cell applications can be achieved by using methanol instead of hydrogen as a fuel [11].In [12], the authors present a turbine which, when installed, does not reduce the mass flow rate in the treated sewage water distribution line.The turbine power is determined using computational fluid dynamics (CFD).The energy potential of water supply systems is determined in [13].The analyses take account of hydropower recovery and thermal recovery.The proposed method can be applied for hydroelectricity production, as well as for heat and cold recovery.The authors of [14] present the use of a propeller-driven pyroelectric turbine to recover approximately 10% of the energy needed to pump water through the filtration system in an existing building.Energy recovery by means of water turbines at the outlets of wastewater treatment plants is proposed in [15].A comparison between a system consisting of a hybrid Pump-As-Turbine (PAT) and a pilot solar system with a traditional diesel generator in an off-grid farm in southern Spain is shown in [16].The PAT was installed on a bypass to ensure continuous operation of the irrigation network.It only represented a maximum of 8% of the total cost of the environmental protection system.
Much less frequent attempts are made to recover energy from a medium containing a solid phase.The presence of solids causes erosion, which makes operation difficult, and the amount of energy that can be captured from the flowing multiphase mixture is smaller compared to the single-phase flow [17,18].The aim of this work is to analyze the possibility of recovering energy from gravity flows occurring in the pipelines of a copper ore enrichment plant.The analysis will be carried out for a system operating in a plant in Poland.An energy recovery device will be designed taking into account the flow of a two-phase medium.For this purpose, flow calculations will be performed.
The literature offers a great many examples of numerical research focused on determining the power output of water turbines and on the analysis of flows through pipelines, particularly in single-phase simulations, predominantly concerning water flows.However, there are fewer publications addressing two-phase flows with a small amount of the solid phase.With regard to water turbines, the most common investigations revolve around studying the impact of the solid phase on erosion in pipeline components or turbine blades.The authors of [18] present a simulation of the Pelton turbine with the flow of solid particles suspended with the solid phase maximum volume concentration of 2% only.However, it is challenging to find, in the literature, examples of simulations of water turbines within pipelines where the suspension concentration exceeds 50%.The numerical models for such scenarios still need more adequate validation; they are often limited to very basic geometries, such as flows through a pipe with bends.
When it comes to numerical models for multiphase simulations, several approaches can be distinguished.The Lagrangian method analyzes the motion of each solid-phase particle as a point with a certain mass.In this approach, there are no collisions between particles because they are simulated as points.In [19], the researchers show that the results obtained using the Lagrangian approach are worse compared to the Eulerian approach, and the simulation time is longer for the latter.This method can be applied in flows with a volume concentration of up to 10%.An extension of the Lagrangian method is the Discrete Element Method.In this approach, physical particles are simulated, each with its own mass and diameter.It is the most accurate method for modeling the dispersed solid phase in a liquid, but it becomes computationally expensive at a significant number of particles [20,21].On the other hand, Eulerian approaches to multiphase-flow modeling can be distinguished further.If the phases are immiscible and there is a large multiphase surface area, the Volume of Fluid (VOF) method can be applied [22].Various turbulence models are also tested in this work.If the phases are dispersed among each other, the Two Fluid Model (TFM) should be used instead of the VOF method [17,20,[23][24][25].
The TFM approach makes use of a closure model, utilizing the kinetic theory of granular flows with solid particles.This model incorporates parameters such as the solid phase stress, the interphase momentum transfer, the interactions between particles and between particles and the continuous medium.An accurate prediction of multiphase flows with the use of this model relies heavily on the interphase frictional force.It is this force that shapes the velocity field and relative velocities between phases.This article will present the validation of the Two Fluid Model based on experimental data for a straight pipe and a pipe with a bend.Subsequently, for the designed turbine structure, single and multiphase simulations will be performed using the TFM.

Visual Inspection
The purpose of the site visit was to select a location for the energy recovery device and to obtain the necessary data.Information was obtained regarding the operating hours of this part of the installation from the day of its construction to the on-site inspection.On this basis, the average percentage of the installation's operational time was determined.Measurements were made of the maximum height difference and the distance of the pipeline from other elements of the system.Diameters, thicknesses and fragments of the external geometry of the pipeline were also measured to verify the received technical documentation.Laser distance meter KRT706500 with an accuracy of 2 mm and thickness gauge PCE-TG 75 with an accuracy of ±0.5% of the measured value were used to determine the pipeline geometry.The uncertainty analysis should include the uncertainty of the experimenter caused by increased humidity and high dustiness of the building area.
The lower section of the outflow pipeline, which is shown in Figure 1, was chosen as the installation site for the energy recovery device.The location of the proposed energy recovery device is shown in Figure 2. The selected location allows for the use of a piping system that bypasses the energy recovery system.The purpose of the so-called bypass is to ensure the continuity of operation of the installation in the event of a failure or service of the energy recovery system.Additionally, the bypass protects the energy recovery system against damage during its start-up.A decision was made then to design an inexpensive, simple prototype of a turbine.
The assumption is that the rotor wear and tear will proceed quickly due to strong friction, causing considerable erosion of the blades.Figure 2 presents the pipeline with the proposed turbine consisting of three blades, the shaft casing and eight stationary blades upstream and downstream of the rotor, all of which make up the supporting structure of the proposed solution.The turbine was installed in the lower section of the pipeline.Elements of the designed device that are in contact with the medium flowing in the pipeline, with a significant relative flow velocity in the boundary layer, must be covered with coatings of wear-resistant materials to protect against excessive wear resulting from erosion/abrasion.As part of the work carried out, two groups of wear-resistant materials suitable for use in this application were selected: wear-resistant coatings made of polyurethane (PU) and wear-resistant coatings made of silicon carbide (SiC).

Mathematical Model
The continuous and the dispersed phase in a multiphase flow will be described using the Eulerian two-phase model.In this model, space can be occupied only by one phase.
Based on that, the concept of volume fraction i α of phase i is introduced The laws of mass and momentum conservation are written and solved individually for each phase i.The continuity equation for phase i has the following form The general form of the momentum conservation principle in fluid dynamics is described by Equation (3).In the analysis of the single-phase flow, the changes in momentum of a fluid element depend only on gravitational forces, external pressure and internal viscous forces within the fluid.These forces can be divided into body forces per unit volume b f  , exemplified by the gravitational force, and stresses σ, exemplified by the fluid pressure or shear stresses caused by the fluid viscosity: For a fluid, the stress tensor σ is often written as the sum of normal stresses and shear stresses, given by Equation ( 4).
Equation ( 3) can be written for each of the phases involved.
where g is gravitational acceleration [m/s 2 ]; D A is the liquid-solid phase interphase momentum exchange coefficient [kg/(m 3 s)]; vm C is the virtual mass force coefficient [21], The first term on the right side of Equations ( 5) and ( 6) represents surface forces associated with pressure.The second terms represent shear stresses induced by the fluid viscosity.The third terms represent the body forces resulting from gravitational interaction.The subsequent terms describe the interactions between the phases being described.In the fourth terms, there are the so-called interphase drag forces, which influence the slip and relative velocity of each phase.The next forces (terms 5) are associated with acceleration.The inertia of the surrounding fluid affects the acceleration rate of a particle submerged in the flow.In the inviscid flow theory, this effect can be shown as "virtual mass" or "added mass".The final, sixth components relate to the lift forces.In situations where the continuous phase flow field exhibits non-uniform or swirling characteristics, particles within the flow experience a phenomenon known as the lift force.This lift force is perpendicular to the direction of the relative velocity between the particles and the surrounding fluid.

The interphase drag coefficient D
A is determined as follows [20]: where Re s is the particle phase Reynolds number defined as The solid phase stress tensor s τ is expressed as where s ξ is the solid phase bulk viscosity [kg/(m•s)]; s μ is effective granular viscosity [kg/(m•s)]; and I is the identity tensor [-].
Two regimes of the solid phase flow are possible: kinetic and frictional.The kinetic regime occurs below the concentration limit ,max s α , which for spherical particles is 0.624.
In the kinetic regime, collisions of particles dominate in the flow and behave like an ideal gas in the kinetic theory of gases.On the other hand, if the movement is primarily influenced by contact between particles (close to the concentration ,max s α ), empirical formulae borrowed from the field of soil mechanics are employed [26].
Two regimes were included in this paper.Below critical volume concentration crit α with the value of 0.55, the kinetic regime was activated and the solid phase bulk viscosity s ξ in Equation ( 10) was calculated using the following formula [27] 2 0

4
(1 ) 3 where 0 g is the particle radial distribution function [-]; s θ is the granular temperature [m 2 /s 2 ]; and e is the restitution coefficient, with a value of 0.9 [-].
Next, the particle radial distribution 0 g was calculated using Equation ( 12): In the kinetic regime, effective granular viscosity s μ was computed as a sum of col- lisional and kinetic contributions The collisional term was calculated using the following equation [28]: while the kinetic term is expressed as Granular temperature can be computed using the granular energy transport equation: The dissipation of granular energy γ can be expressed as The dissipation of granular energy due to the interphase drag p J was calculated using the following formula: Above the critical value of volume concentration crit α , the frictional regime was ac- tivated and the solid phase bulk viscosity s ξ was set as 0 in Equation ( 10) [26].Then, the solid pressure term in this equation was calculated using the following formula: where F is a constant with a value of 0.05 [kg/m•s 2 ]; k is a constant with a value of 2.0 [-] and h is a constant with a value of 5.0 [-].
Then, in the frictional regime, effective granular viscosity s μ in Equation ( 10) is equal to frictional viscosity defined as where φ is the internal friction angle with the value of 28.5°; p D is the strain rate tensor The turbulence modeling process involved the application of the k ε − realizable turbulence model independently for each phase.This approach enabled a comprehensive analysis of the turbulence behavior within the fluid system under consideration.The Simcenter Star CCM+ version 2210 [29] software facilitated various tasks, including the mesh generation, the simulation performance and the post-simulation analysis.
It is noteworthy to emphasize that all computations were executed under steadystate conditions, leveraging the capabilities of the SIMPLE solver, also known as the Segregated solver within the software framework.By choosing this solver, efficient convergence and stability were ensured throughout the simulation process.Furthermore, the equations were discretized using second-order schemes, enhancing the accuracy of numerical solutions while maintaining computational efficiency.Multiphase flows with the liquid and the solid phase are not an easy task in terms of convergence.
To obtain the simulation solution in a multiphase simulation, the following underrelaxation factors (URFs) were selected: The calculations were performed using a dedicated Prometheus computing cluster located at the Academic Computer Center Cyfronet AGH.A node consisting of 256 cores was used for the calculations.Regarding the simulations of the pipeline with the proposed turbine, for the single-phase simulation, the average computation time and number of iterations were 20 min and 1800 iterations.For the multiphase simulation, considering the slurry, the average computation time was 4.5 h and 140,000 iterations (the effect of the low value of URF coefficients of most equations).

Validation of the Numerical Approach
The flow of the liquid phase with the solid phase in the pipeline under consideration was analyzed numerically.The first part of this section presents experimental verification of the numerical model based on the flow through a straight pipe [30] and through a pipe with a 90° knee [17,24].The numerical model verification started with the flow through a straight pipe [17,30].The pipe length was 100•D, as shown in Figure 3.The diameter D = 0.130 m.Next, discretization was performed on this geometry using a polyhedral mesh.The polyhedral mesh type was selected considering that the combination of the pipeline with the proposed turbine intended for energy recovery had a rather complex structure difficult to represent using a hexahedral mesh.The mesh density was modified using only one parameter, referred to as "Base Size".In every mesh, the number of inflation layers was set to 6, and the coefficient of their growth rate from layer 1 was set to 1.2.Three meshes were generated to assess the impact of the concentration of the cells on quantitative results (here: the pressure drop between the inlet and the outlet).Figure 4   Next, the boundary conditions were assumed.At the inlet, velocity was assumed as 1.0, 2.0, 3.0, 4.0, 5.4 and 6.0 m/s for the liquid phase and the suspension, and the suspension volume fraction was assumed as 0.4.In all simulations conducted in this work, no special velocity profile was assumed at the inlet.Each cell had exactly the same velocity value.Moreover, it was assumed in the context of two-phase flows that the mixture at the inlet was homogeneous.A gravitational vector was assumed in the direction perpendicular to the flow, with a value of 9.81 m 2 /s.The density of the water and sand making up the suspension was adopted as 998.2 and 2650.0 kg/m 3 , respectively.The diameter of the sand grains was set as 2.7 × 10 −4 m [17,30].
After the simulations, the first to be presented was the correlation of the suspension content Cs at a given relative height zn.The relative height was defined as the quotient of the vertical coordinate and the pipe diameter zn = z/D, where 0 means the pipe bottom and 1 is the pipe top.The distribution was measured at length L = 100•D from the inlet.The chart presents the results for three meshes compared to experimental data [30].
The chart in Figure 5 shows good correlation for points located below a relative height zn = 0.77.The error level does not exceed 10%.Analyzing the impact of the computational mesh on the vertical distribution of the suspension, it can be seen that the ME-DIUM and FINE meshes give practically the same results, whereas the COARSE mesh type displays only a slight difference in the distribution compared to the MEDIUM and FINE meshes.
The next chart in Figure 6 presents the pressure gradient defined as Δp/L depending on the inlet velocity v.The pressure drop was calculated as the difference in total pressures between the pipe inlet and outlet.In Figures 5 and 6, it can be seen that all the meshes show good convergence in terms of the order of magnitudes and goodness of fit to experimental data.The COARSE-type mesh differed slightly from the MEDIUM and FINE meshes, which showed very good convergence in terms of values.Based on the charts presented, it was decided that the sizing generated for the MEDIUM-type mesh was sufficient to perform further simulations (flow through the pipe with a knee and through the system with the turbine).
Then, on a few selected sections, the sand volume fraction was presented graphically (cf. Figure 7).The sections were created, respectively, for 0.0, 0.25, 0.5, 0.75 and 1.0 of the pipe relative length x/L, where x was the length coordinate.The above-presented distributions of the suspension content in individual sections illustrate the course of the sedimentation phenomenon.Starting from relative length x/L = 0.5, considerable separation of layers and sand sedimentation can be observed on the pipe's lower part due to gravitation.Figure 7e shows the suspension distribution at the pipe outlet.It can be noticed that in the pipe's upper part, there is no sudden growth in the suspension concentration, as shown by the measurement results presented in Figure 5.
The next stage of the verification was to determine correlations for the flow through a pipe with a 90° bend [24].In this case, the diameter was D = 0.053 m, the bend radius was R = 0.1488 m, and the total length of the pipe with the bend was 7.234 m.
The mesh generated for this pipe was similar to the MEDIUM-type mesh for the straight pipe, as shown in Figure 4b.The mesh is visualized in Figure 8 below.
In this verification, one velocity of the suspension at the inlet was simulated −3.56 m/s.The sand volume content was assumed at the level of 0.0882 [24].Gravitation was taken into account in the simulation.The density of the water and sand making up the suspension was adopted as 998.2 and 2650.0 kg/m 3 , respectively.The diameter of the sand grains was set as 4.5 × 10 −4 m [17,24].
After the simulations, the suspension normalized content was presented: Csn = Cs/Cs,i, where Cs,i is the suspension content at the inlet, totaling 0.0882 [24].Normalized height zn is defined in the same way as for the straight pipe.The measurements were performed for the three sections (C1, C2 and C3) marked in Figure 9.The figures presented above show fairly good correlations.The general trend in the suspension distribution has been preserved.Figure 10a-c for zn lower than 0.4 show much worse goodness of fit compared to values higher than 0.4.Globally, however, the simulations predicted the forming velocity field in a two-phase flow quite well.The next verification was prepared for the pressure drop.In this case, the pressure drop was determined based on measuring points P1-16 presented in Figure 9.The chart presents the normalized pressure drop Δpn = Δp/(ρv 2 /2) as a function of the distance from point P1.Pressure at point P1 was adopted as the reference pressure.The correlation is presented in the figure below.
As presented by the chart in Figure 11, the correlation is very good.Relative errors are less than 10% for all the measuring points.Slight deviations appear for pressure drops measured in the bend area, which is due to local losses and eddies within the bend.The next figure illustrates the suspension distribution in a few selected sections of the pipe with a bend.The following sections were selected: the pipe inlet, the bend inlet, the middle of the bend arc length, the bend outlet, and the section for point P15, based on Figure 9. Figure 12 shows that the heavier fraction, sand in this case, is subject to sedimentation upstream of the bend due to gravity (Figure 12b), then under the influence of the resultant force of gravity and the centrifugal force, it is pressed into the lateral outer wall of the pipe (Figure 12c,d), to finally settle at the pipe bottom wall again due to gravitation, with the separation of layers in the flow (Figure 12e).Summing up this section of the paper, both cases the suspension flow through a straight pipe and the flow through a pipe with a bend show that the defined mathematical model correlates well with the experimental data and can be applied for the suspension flow through a pipeline, where changes occur in the flow direction.

Numerical Simulation of a Pipeline with a Turbine
This section presents simulations of a pipeline with a turbine using the single-and the two-phase approach.The main goal here was to make a quantitative comparison of the pressure drop, efficiency and power for a pipeline with pure water and with a mixture of water with sand.In each approach, the inlet velocity and a boundary condition of the pressure outlet type were assumed.The velocities were selected to establish a point at which the turbine efficiency began to be higher than zero, and then find the maximum point to see how high efficiency would be when this maximum point were crossed.The single-phase simulation was performed without taking account the gravity equation as it has no impact on the solution of the system of equations.The so-called multiple reference frame approach was adopted into the rotational domain, with the rotational speed adopted at 60 rpm (the rotational speed adopted while selecting the angle of incidence of the blades).The accuracy of the proposed method was previously tested using Laser Doppler Anemometry [22].For the multiphase analysis, the inlet volume contents of water and of the solid phase were assumed with their respective values defined as 0.487 and 0.513.The grain average size in the suspension was adopted as 5 × 10 −4 m.The water and the solid phase density were assumed as 998 kg/m 3 and 2660 kg/m 3 , respectively.In addition to sand, the solid phase also included copper ore fractions.Figure 13 presents a visualization of the boundary conditions in the pipeline under analysis.
The geometry was then discretized using a polyhedral mesh (cf. Figure 14).The same settings were applied as for the MEDIUM-type mesh in the straight pipe.Additionally, the mesh was made denser in the area of the blades, as this is where key flow phenomena, such as layer separation, secondary vortices and large velocity gradients occur.These phenomena have an impact on the momentum transfer from the fluid to the rotor, and simulating them correctly is crucial for the analysis accuracy.Using Equation ( 23), a chart was made showing the turbine efficiency based on the quantitative data gathered after the simulations performed for both cases, i.e., the singleand the two-phase approach (cf. Figure 15).The chart indicates that the efficiency for the flow with the suspension was more than 20% lower for the turbine maximum working point compared to the single-phase flow.The optimal output for the designed turbine in the slurry flow was about 0.6 m 3 /s.Moreover, the point of maximum efficiency was shifted toward a higher volume output, and the characteristic curve of the useful operating range was much narrower.Moment M acting on the turbine is expressed by Equation ( 24): ( ) where rf-is position of the cell's face relative to origin of coordinate system describing the MRF,  Next, calculations were made of the pressure drops measured from the inlet to the outlet in the pipeline.The results are presented in the form of the height of the liquid column losses h in Figure 16 according to Equation (25), where the densities of the water and of the suspension were adopted as 998.2 and 1850.6 kg/m 3 , respectively.These drops are necessary to calculate the power that the turbine is able to generate depending on the volume output, taking account of flow-related pressure losses.As indicated by the chart, the pressure drops for the slurry flow are by more than 1.2 m bigger compared to the flow of pure water.For the optimal point, the pressure drop totalled 1.42 m for the slurry flow.Formula (26) was then used to calculate mechanical power P that the turbine is able to generate depending on the output, efficiency and flow-related losses.The power is presented in Figure 17.For the optimal efficiency in the slurry flow, power totaled 36,590 W, whereas for the same output for pure water it is 37,878 W. Despite the lower efficiency for the slurry flow and higher flow losses, power is compensated for owing to the almost twice higher density of the medium.
( ) where h Δ is the difference in the height of the liquid column expressed in meters.The next step was to present a 3D visualization of the suspension concentration in the entire pipeline (cf. Figure 18) for the optimal working point, i.e., for the output of 0.6 m 3 /s.The distribution illustrates how the liquid phase and the solid phase behave depending on the pipeline inclination angle, the bends on the flow path and the rotor (cf. Figure 19).The solid phase separation can be seen on the bends upstream of the turbine.By its rotation, the turbine pushes the suspension into the upper part of the pipe.Due to the gravitational force, sand sedimentation re-occurs in the pipe lower part.A close-up of the turbine area is presented below to show the behavior of the solid phase more accurately.It can be seen how the rotating rotor sweeps away the sand and pushes it into the upper part of the pipe.Figure 22 shows the distribution of the y+ parameter on the surfaces that compose the turbine's construction.The graph is shown for pure water flow through the pipeline for a flow rate of 0.6 m 3 /s.Classically, this value should oscillate around y+ = 1.However, achieving such a low value would involve a very large number of elements in the wall layer, which would increase the size of the computational grid considerably.In addition, the All+ Wall Treatment in the turbulence model in all simulations was used [31].This means that the values presented in Figure 22 were acceptable values for this turbulence model.
Summing up, a simple-structure turbine used in the pipeline under analysis makes it possible recover more than 35 kW of power for the optimal working point.A problem that remains is the solution lifetime.However, considering the costs of the replacement of the blades, which are not expensive to make, the system may bring financial benefits from the recovery of electric energy.

Conclusions
This paper presents an idea of a simple turbine solution to recover mechanical energy from a pipeline with a slurry flow involving a high concentration of sand.The efficiency of the mechanical energy recovery was evaluated using CFD.The numerical model of the turbine was verified against experimental data for both a straight pipe and a pipe with a knee.The analysis of the results, including sand concentration distributions and pressure drops, indicates that the numerical model simulated the turbine behavior well compared to experimental data.This suggests that the model is suitable for predicting the turbine efficiency.
Further research involved simulating both pure water and slurry flows through the turbine.It turned out that despite a 20% lower efficiency of the turbine at the flow with a sand suspension, the powers obtained in the two cases were similar.This was the effect of the slurry density, which was by more than 80% higher and thus effectively compensated for the turbine lower efficiency during the slurry flow.At a constant water output of about 0.6 m 3 /s, for which the turbine efficiency at the slurry flow reached the maximum value, the mechanical power on the turbine shaft was about 35 kW.These results suggest that a simple turbine can be an effective tool for recovering mechanical energy from slurry pipelines.More experiments and simulations could be performed to further understand and optimize the process.

Figure 2 .
Figure 2. Pictorial view: (a) top view and (b) side view, with an opening to show the installed turbine; (c) front view and (d) isometric view of the turbine with basic dimensions.
as 0.5 [-]; L C is the lift force coefficient, set as 0.25 [-]; l τ is the liquid phase stress tensor [kg/m•s 2 ]; s τ is the solid phase stress tensor [kg/m•s 2 ] and the subscript i indicates the phase (i = l liquid or i = s solid).
where l μ is the liquid dynamic viscosity [kg/(m•s)]; s d is the solid particle mean diam- eter [m]; tr α is the transition volume fraction, set as 0.2 [-]; D C is the particle drag co- efficient [-].Next, the particle drag coefficient D C was calculated using the following formula [20]: 0.687 24 (1 0.15 Re ) 0 Re 1000 Re 0.44 Re 1000

[s − 1
] and ,max s u is the maximum solid viscosity, with a value of 1000 [kg/(m•s)].

Figure 3 .
Figure 3. Verification of the numerical model on a straight section of the pipeline.
presents the computational mesh visualization for three different settings of the mesher.The number of the generated elements totaled 0.3 M, 0.7 M and 1.5 M.

Figure 5 .
Figure 5. Distribution of the suspension content Cs at the pipe outlet measured along normalized height zn.

Figure 8 .
Figure 8. Computational mesh for the pipe with a knee.

Figure 9 .
Figure 9. Geometry of the pipe with a knee used for the verification[24].

Figure 10 .
Figure 10.Distribution of the suspension normalized content Cs at the pipe outlet measured along normalized height zn; (a), (b) and (c) present results for sections C1, C2 and C3, respectively.

Figure 11 .
Figure 11.Normalized pressure drop depending on distance from measuring point P1.
where M-moment on the turbine [N•m], ω -angular velocity [rad/s], Q  -volume output [m 3 /s], , T in p -total pressure upstream of the turbine and [Pa] , T out p -total pressure downstream of the turbine [Pa].
[N/m 2 ], and a -vector describing the revolving axis used in MRF; subscript f indicates the cell's face.

Figure 15 .
Figure 15.Chart illustrating the turbine efficiency measured in the pipeline depending on the volume output and presented for the single-and the two-phase simulation.
difference in pressure drops between the pipeline inlet and outlet [Pa], g is gravitational acceleration [m/s 2 ] and ρ is the fluid density [kg/m 3 ].

Figure 16 .
Figure 16.Pressure losses expressed in the form of the height of the liquid column for the singleand the two-phase flow.

Figure 17 .
Figure 17.Changes in the turbine power for the single-and the two-phase flow depending on the volume output.

Figure 18 .
Figure 18.Visualization of the suspension concentration in the pipeline: (a) side view; (b) front view.

Figure 19 .
Figure 19.Three-dimensional visualization of the suspension concentration in the pipeline-closeup of the turbine area.The suspension distribution is shown in the figure below in sections S1 and S2 to illustrate the behavior of the suspension in the pipe section (Figure20).Section S1 shows the separation of layers in the mixture due to the centrifugal force acting in the bends and due to gravitation, whereas section S2 shows the sand ejection toward the pipe top wall due to the impact of the rotor.

Figure 20 .
Figure 20.Distribution of the suspension in sections S1 and S2: (a) visualization of the sections, (b) section S1, and (c) section S2.

Figure 21
Figure21below presents a three-dimensional visualization of the velocity field using streamlines.The velocity is defined by Equation (27) and presents the weighted average for the solid and the liquid phase, where the weights are volume fractions of the phases.

C
is the water content [-], W v  is the velocity vector for the liquid phase [m/s], s C is the sand content and s v  is the velocity vector for the solid phase [m/s].

Figure 21 .
Figure 21.Velocity field visualization illustrated by streamlines in the turbine.

Figure 22 .
Figure 22.Distribution of the y+ parameter on the surface of the turbine.