A CFD Validation Effect of YP/PV from Laboratory-Formulated SBMDIF for Productive Transport Load to the Surface

: Several technical factors contribute to the ﬂow of cuttings from the wellbore to the surface of the well, some of which are fundamentally due to the speed and inclination of the drill pipe at different positions (concentric and eccentric), the efﬁcacy of the drilling mud considers plastic viscosity (PV) and yield point (YP), the weight of the cuttings, and the deviation of the well. Moreover, these overlaying cutting beds breed destruction in the drilling operation, some of which cause stuck pipes, reducing the rate of rotation and penetration. This current study, while it addresses the apropos of artiﬁcial intelligence (AI) with symmetry, employs a three-dimensional computational ﬂuid dynamic (CFD) simulation model to validate an effective synthetic-based mud-drilling and to investigate the potency of the muds’ ﬂow behaviours for transporting cuttings. Furthermore, the study examines the ratio effects of YP/PV to attain the safe transport of cuttings based on the turbulence of solid-particle suspension from the drilling ﬂuid and the cuttings, and its velocity– pressure inﬂuence in a vertical well under a concentric and eccentric position of the drilling pipe. The resulting CFD analysis explains that the YP/PV of SBM and OBM, which generated the required capacity to suspend the cuttings to the surface, are symmetric to the experimental results and hence, the position of the drill pipe at the concentric position in vertical wells required a lower rotational speed. A computational study of the synthetic-based mud and its potency of not damaging the wellbore under an eccentric drill pipe position can be further examined.


Introduction
An effective hole cleaning [1] comes with a high expectation of formulating a good drilling mud with accurate rheological properties to remove cuttings from the wellbore. Cuttings [2] from boreholes are intended to rise to the surface while drilling operations are ongoing. This process of lifting cuttings is crucial for drilling operators to avoid some drilling challenges faced, such as lost circulation, stuck pipes, poor rate of penetration, and drilling bit wear. A machine learning approach had been adopted to predict stuck pipe incidents in the drilling operation to decrease non-productive time [3].
Some conventional wells are challenged with this process as well as horizontal wells at an inclined position [4,5], posing a threat to reaching the desired depth of drilling. A higher rate of penetration in a wellbore means that the cuttings from the wellbore are lifted mainly because the mud under operation is assisting the drilling bit to cut fresh rock with little exertion of force; this further explains why it is necessary to keep the mud design in an accurate state, suitable for the wellbore and the effective removal of cuttings.

•
To validate the effective synthetic-based drilling mud formulated from the drilling fluid laboratory at Nazarbayev University.

•
To validate the distribution of synthetic-based mud particles and cuttings in the wellbore.

•
To investigate how rotating the drill pipe affects drilling mud transport capacity.

•
To examine the effects of fluid hydraulics on transport cuttings based on the turbulence of solid-particle suspension. Table 1. Effect of YP/PV on safe transport.

Author Yield Point (YP) Plastic Viscosity (PV) Transport Index (TI) Average TI Mud Type
Wayo's Experiment [6] 136 The novel results of this study will lead to a deeper understanding of the significance of synthetic-based and oil-based mud, perfect for the safe transport of cuttings from the wellbore to the surface. For effective fluid circulation, operators will be able to decide where on the annulus to modify pressure and velocity for effective drilling optimization. However, the arduous challenges combating the oil and gas industry keenly access AI symmetry techniques to simulate, model, and better still, automate the solutions for drilling optimization.

Empirical Data
The crucial response of CFD modelling will be a validation for the experimental data collated from the Nazarbayev University drilling laboratory. This is accompanied by a fairer comparison of data from previous research conducted by other authors [38]. Table 1 showcases rheological data considering only two properties that are a determinant of the safe transporting of cuttings [39] from the wellbore to the surface. An API standard [40] for drilling rheology explains that the ratio of yield point (YP) and plastic viscosity (PV) can be a determinant of an effective drilling fluid capable of lifting cuttings from the well. The ratio of YP and PV gives the transport index (TI), and for a safe index, the ratio is stipulated to be around 0.75-1.5 [6,35]. Earlier studies formulated several drilling fluids, but Table 1 gives structured data to reduce data noise. However, the data emanated from two different types of mud design, synthetic-based and oil-based mud. Each of these contains barite and bentonite as solid [41] particles and oil and water as liquid fluids, as demonstrated in Figure 1.

CFD-Eulerian-Eulerian Model
CFD can solve complex radical problems, and this study finds Matlab-CFD an insightful choice for the demonstration of solid-particle dispersion tracking, particle-particle collision, and multiphase flow velocity-pressure monitoring. In modelling this current study, Eulerian-Eulerian (EE) model analysis was considered for the fluid flow of various particle size distributions in the wellbore. Since this chosen model tends to handle particle-to-particle interactions, turbulence, and solid volume fractions all in the same domain, it is regarded as the best, however, Figure 2 provides the algorithm for computation.

CFD-Eulerian-Eulerian Model
CFD can solve complex radical problems, and this study finds Matlab-CFD an insightful choice for the demonstration of solid-particle dispersion tracking, particle-particle collision, and multiphase flow velocity-pressure monitoring. In modelling this current study, Eulerian-Eulerian (EE) model analysis was considered for the fluid flow of various particle size distributions in the wellbore. Since this chosen model tends to handle particleto-particle interactions, turbulence, and solid volume fractions all in the same domain, it is regarded as the best, however, Figure 2 provides the algorithm for computation.
CFD can solve complex radical problems, and this study finds Matlab-CFD an insightful choice for the demonstration of solid-particle dispersion tracking, particle-particle collision, and multiphase flow velocity-pressure monitoring. In modelling this current study, Eulerian-Eulerian (EE) model analysis was considered for the fluid flow of various particle size distributions in the wellbore. Since this chosen model tends to handle particle-to-particle interactions, turbulence, and solid volume fractions all in the same domain, it is regarded as the best, however, Figure 2 provides the algorithm for computation.  The trajectory of the concentrated particles under this model reduces the operational cost for predicting or forecasting solutions to flow impedance in the wellbore or the essence of pressure drops during circulation. Moreso, the limitation of this current study under the EE model is based on the assumption that the particles are all considered spherical. For this reason, emphasis is made on transport equations such as continuity of flow and momentum balance. Particle-to-particle interactions [43,44] and their rigorous movement (turbidity) make a sense of collisions with their bodies which are a result of momentum and pressure, and this, however, is the parameter for EE modelling. Notwithstanding, the model for this study prioritizes the fluid phase under solid-liquid exchanges.

Drag Adaptation
The resistance to particle-free flow is attributed to drag. Previous studies in empirical analysis and on fields have tested the drag force as a means of capitulating the free flow of solid fluids in its domain. Others are of the view that the particles' maximum velocity trajectory is influenced by their shape and size [20]. This current study limits its predictive novelty to spherical particle analysis. The coefficient of drag is a determinate of particle sphericity and a function of the Reynolds number; the lower the number, the laminar the flow and the higher the number showcases the turbidity of the flow.
Particles under the influence of drag force in a microscopic view are easy to analyse and are not in large quantity. It is hard to determine the influence of drag by a large volume of particles and in different sizes. The influence of particle sizes and shapes from the previous empirical results is considered as a whole. However, the moderated drag models published from Epelle-Gerogiorgis (EG), Refs [45][46][47] were considered for adaptation.

Model Assumptions
Satiating the current study's validity, the CFD-Eulerian-Eulerian model is conducted based on an informed assumption below that:

•
The flow of solid-liquid particles under surveillance is in a continuous phase. The walls of the drill pipe are smooth.

Governing Equations
The volume fraction of the solid-liquid assumed uniform flow in the annuli of the wellbore in question is expressed by the continuity equation.
The summation of all momentum acting on the solid-liquid phases describes the force, mass, and velocity of the solid-fluid movement in the wellbore given by: Equations (1) and (2) are models considered from Epelle-Gerogiorgis (EG). The flow continuity and momentum parameters [48][49][50][51] involve the volume fraction solid phase a s , solid phase density ρ s , liquid phase density ρ l , liquid velocity  (2), turbidity force is made relevant because of the high-pressure injection of drilling fluids and the rotational speed of the drilling pipe in the wellbore.
For the efficient review of the solid-liquid exchange coefficient K sl , as adapted from [2,20,31,52], previous studies explain that when the volume fraction of the liquid phase a l > 0.8, then K sl is transformed into: where C D is the drag coefficient, where Re s is the Reynolds number of the solid-particle phase, and, in a case where a l ≤ 0.8, then, The drill pipe is rotating and drilling fluids [53] are released from the drill bit's nozzles to improve the cutting of the formation bed. Therefore, the current study does not anticipate that the fluid will flow in a uniform laminar flow. The pressure and velocity scale under this circumstance would create a turbulence model, where Reynolds numbers are expected to be high. The rate of dissipation (k), and the kinetic energy (ε) for this model were considered dependent on the factors evolving around the transport equation in a single-phase flow: The turbidity of the particles in Figures 3 and 4 is the focus of the investigation, verifying how the physical properties of these formulated muds (synthetic and oil) in the laboratory can be used to validate a CFD simulation to enhance predictive analysis. The turbulence viscosity of the particles under investigation is intermittently linked to the rate of dissipation and the kinetic energy, as adopted in Equation (9):    The modelling of computational fluid flow does not only consider the above equations adopted from previous studies, the iterations of the pressure-velocity profiles of the mud capabilities were also considered using the Eulerian-Eulerian formulas active in the MATLAB FEATool Multiphysics.

CFD Model Implementation
The Eulerian-Eulerian model proposed in this study was implemented to define the efficacy of mud, demonstrating its tendency to carry cuttings in the wellbore to the surface. The finite element method was adapted while conducting the simulation using MATLAB FEATool Multiphysics. This pattern of computation focuses on the continuity of flow, momentum, and transport equations. Table 2, based on previous experimental analysis, gives explicit parameters needed for the simulation, however, the basis of the two-dimensional flow regime was factored on several factors including mud-nozzle inlet velocity, the velocity of the cuttings in the wellbore and the volume of force acting in the x-y direction.

Geometry, Grid & Boundary
A cross-sectioned wellbore coupled with drilling pipe and casing was used for this analysis; liquid-solid interactions between the outer walls of the drill pipe and the inner walls of the casing were the interest of the investigation. The trajectory of the particle's collision with the walls of the casing and drill pipe were set to the limit boundary. The space in between the casing and the drill pipe is termed the annulus. Additionally, the domain for this simulation was set up in both concentric and eccentric locations using Equation (10), where E is the distance between the centres of the inner (R i ) and outer (R o ) tubes of the annulus, as demonstrated in Equation (10) and Figure 5. The grid for the analysis maintained the shape and size set by the geometry to support running the set of equations governing this analysis. The grid was set at 0.03 to discretize the domain efficiently. Also, Figure 6a,b exposes the true cross-sectional view and boundaries of the wellbore, and in further assumptions stated above, the right wall (boundary 2) was considered the outer wall of the drill pipe which was assumed smooth. Boundary 4, being the left wall was named the inner wall of the casing, while boundary 3 was classified as the annulus of the drilling column; it is termed to be the inlet/outlet boundary. No stress or neutral boundary was set to have boundary 1; this was because the focus of the investigation was limited to finding the behaviour of liquid-solid interactions between the casing and drill pipe, so long as an established assumption of cuttings was already in the loop of investigation.
Symmetry 2022, 14, x FOR PEER REVIEW 9 Boundary 4, being the left wall was named the inner wall of the casing, while bound was classified as the annulus of the drilling column; it is termed to be the inlet/o boundary. No stress or neutral boundary was set to have boundary 1; this was bec the focus of the investigation was limited to finding the behaviour of liquid-solid int tions between the casing and drill pipe, so long as an established assumption of cut was already in the loop of investigation.

Results and Discussion
The computational liquid-solid modelling considered the schematic flowchart in the above Figure 2 as the basis to determine a two-dimensional fluid particle distribution, dispersion, and turbidity. These cuttings emanating from the formation beds were simulated using the pressure (80 psi) from the drilling fluids under review to force out, or clean the wellbore, thereby maintaining constant hydrostatic pressure to carry these cuttings to the surface.
Simulated fluids, investigating their rheological properties [54], were simply classified under non-Newtonian fluids. The potentiality of these fluid flows (synthetic and oil-based muds) was tested through a cross-sectioned vertical eccentric annulus. The CFD model development examined the turbulence of the particles holding the assumptions made for its success.
CFD analysis was tuned to hold other parameters constant, as demonstrated in Table 2, while other parameters had different values due to the disparities in fluid rheology [55]. Moreso, the experimental results from Table 1 (results emanating from the laboratory drilling fluids) were keenly used to validate the CFD 2D simulation for this study.
The transport index from the experimental analysis satisfies the API standards [40] for drilling fluids capacity. This, in its original form, had the cutting transport ratio adapted from Rooki's model [29], used for the CFD simulations as demonstrated in Equation (11), where R T is the cuttings transport ratio, V t is the velocity for cuttings transportation, and V a is the velocity of the fluid in the annulus.
The fluid dynamic analysis and the experimental results were duly compared using the transport index or ratio as the basis of analysis, as demonstrated in Figure 7.
for drilling fluids capacity. This, in its original form, had the cutting transport ratio adapted from Rooki's model [29], used for the CFD simulations as demonstrated in Equation (11), where is the cuttings transport ratio, is the velocity for cuttings transportation, and is the velocity of the fluid in the annulus. The fluid dynamic analysis and the experimental results were duly compared using the transport index or ratio as the basis of analysis, as demonstrated in Figure 7.  The non-Newtonian fluids experimented with by Wayo, Okon and Murtaza were analysed by calculating each cutting transport ratio considering yield point and plastic viscosity. The CFD analysis was based on Equation (11); the ratio of transport cutting velocity and the velocity of the fluid in the annulus was defined. Moreover, each model analysis paid attention to the safe transport of cuttings ranging from 0.75 to 1.5. This explains that the higher the safe cutting range, the greater the drag force with a greater proportion of solid cuttings transported through the annulus. The experimental data and the CFD simulations are both agreeable.

Effects of Particle Diameter and Density
The size of the particle and density from the synthetic-based mud had a substantial effect on the cuttings [56] carriage and an impactful increase in the development of turbulent flow. However, fluid motion often has traces of particles in its direction, and the collisions of these particles in an opposite direction. The response time of each collision creates a rise in turbulences and alters the flow patterns. The highest appreciable particle diameter under investigation was 1 mm for both synthetic and oil-based muds. Also, the density of fluids in their entirety, whose motion was in a turbulent flow, was tolerable at over 14,000 kg/m 3 .
The velocities acting on the same particle diameter from different fluids, as demonstrated in grey in Figure 8, spell out the increasing flow of the said particles rushing out to fill the void spaces in the annulus. Moreso, particles from cuttings are noted to have high densities and larger particle diameters (7 mm), and the velocities of the cutting particles (60 m/s) expound the disposition of typical fluid turbulence. ticles (60 m/s) expound the disposition of typical fluid turbulence.
The deposition of smaller particles in high velocities that turn to fill up the void spaces act swiftly to create higher turbidity. Spreading, or the transfer of particles to previously particle-free areas was demonstrated by the 1 mm diameter particles, which further demonstrates the carrier fluid's ability to act continually. Additionally, the velocity profiles demonstrate that when compared to smaller particles, larger particles have a higher propensity to travel along the drill pipe's rotating path. As shown in Figures 8a,b, the combination of particles and cuttings from the corresponding synthetic-based and oil-based muds predict the movement of the particles in various directions. Figure 8 generally shows the movement of concentrated particles from the right to the mid-left (towards the formation walls) during the drilling operations. The deposition of smaller particles in high velocities that turn to fill up the void spaces act swiftly to create higher turbidity. Spreading, or the transfer of particles to previously particle-free areas was demonstrated by the 1 mm diameter particles, which further demonstrates the carrier fluid's ability to act continually. Additionally, the velocity profiles demonstrate that when compared to smaller particles, larger particles have a higher propensity to travel along the drill pipe's rotating path.
As shown in Figure 8a,b, the combination of particles and cuttings from the corresponding synthetic-based and oil-based muds predict the movement of the particles in various directions. Figure 8 generally shows the movement of concentrated particles from the right to the mid-left (towards the formation walls) during the drilling operations.
Most appreciable, the smaller particles (1 mm) were compelled to travel toward the formation's walls by the influence of the drilling pipe rotation.

Effects of Pipe Rotation and Velocity Profiles on SBM/OBM
By further simulating its effectiveness, fluid rheology was evaluated to define its potency, by keeping the pipe rotation constant and the velocity of the cuttings and mud also persistent. Figure 9a,b depict a CFD simulation of both synthetic and oil-based muds. It could be identified that the velocity profiles of these mud particles collide; concentrating on the red arrows, particle-particle collision is imminent. However, the explanation for setting boundaries in Figure 6b applies to Figure 9.
At lower velocities, there are higher tendencies of increased frictional resistance, because the depositions of the cuttings and mud particles on the floors of the formation bed at the lower annulus will adversely reduce the rotation of the drill pipe. Nonetheless, Figure 9 is notably displaying mud particles and cuttings above the annulus; this explains that the velocity of these particles under surveillance is, at its best, not increasing frictional resistance.

Effects of Pipe Rotation and Velocity Profiles on SBM/OBM
By further simulating its effectiveness, fluid rheology was evaluated to define its potency, by keeping the pipe rotation constant and the velocity of the cuttings and mud also persistent. Figure 9a,b depict a CFD simulation of both synthetic and oil-based muds. It could be identified that the velocity profiles of these mud particles collide; concentrating on the red arrows, particle-particle collision is imminent. However, the explanation for setting boundaries in Figure 6b applies to Figure 9. At lower velocities, there are higher tendencies of increased frictional resistance, because the depositions of the cuttings and mud particles on the floors of the formation bed at the lower annulus will adversely reduce the rotation of the drill pipe. Nonetheless, Figure 9 is notably displaying mud particles and cuttings above the annulus; this explains that the velocity of these particles under surveillance is, at its best, not increasing frictional resistance.
Mud particles comprehend the fact that the rotation of the drill pipe rapidly moves these particles toward the formation walls; though we assumed the pipe is smooth and has no slip. The slip velocity for these profiles was considered not to have any adverse effect on the simulation. The EG model elucidates the relevance of slip velocity fluctuations being the results of the speed of the drill pipe rotation, as illustrated in Figure 10. Hence, a lower slip velocity is a ramification of an increased drill pipe rotation. a b c d Figure 9. Velocity Profiles of (a) synthetic-based mud, (b) oil-based mud.
Mud particles comprehend the fact that the rotation of the drill pipe rapidly moves these particles toward the formation walls; though we assumed the pipe is smooth and has no slip. The slip velocity for these profiles was considered not to have any adverse effect on the simulation. The EG model elucidates the relevance of slip velocity fluctuations being the results of the speed of the drill pipe rotation, as illustrated in Figure 10. Hence, a lower slip velocity is a ramification of an increased drill pipe rotation. At lower velocities, there are higher tendencies of increased frictional resistance, because the depositions of the cuttings and mud particles on the floors of the formation bed at the lower annulus will adversely reduce the rotation of the drill pipe. Nonetheless, Figure 9 is notably displaying mud particles and cuttings above the annulus; this explains that the velocity of these particles under surveillance is, at its best, not increasing frictional resistance.
Mud particles comprehend the fact that the rotation of the drill pipe rapidly moves these particles toward the formation walls; though we assumed the pipe is smooth and has no slip. The slip velocity for these profiles was considered not to have any adverse effect on the simulation. The EG model elucidates the relevance of slip velocity fluctuations being the results of the speed of the drill pipe rotation, as illustrated in Figure 10. Hence, a lower slip velocity is a ramification of an increased drill pipe rotation. a b c d Figure 10. (a-d) Slip velocity profiling at different drill pipe rotations, adapted with permission from Ref. [42]. 2018, Epelle, E.I.
Although it has been established that the velocity profile of the cuttings and the mud particles are not in any way increasing the frictional resistance, closing up the space between the walls of the formation and the assumed smooth drill pipe walls would create an eccentric position. Simulation from authors [42,57,58] and this study predict that the particle-particle collisions will intensify; cuttings from the bed would however rise because smaller particles will tend to occupy void spaces causing drag and lift. An eccentric position will adversely create a stuck pipe [59] if the drill pipe rotation for this simulation remains constant; this is because proven fluid viscosity varies, and this is aligned with previous studies delineating that they are non-Newtonian; their shear rates openly lead to cuttings' carriage. For the avoidance of thought and to optimise the fluids carriage rheology, the drill pipe rotation should be maintained constant only at the concentric position, or increased drill pipe rotation at an eccentric position.
The higher the shear region, the lower the viscosity and the lesser the synthetic and oil-based mud transport capacity; this is for the case of an eccentric position. Inversely, the lower the shear region, the higher the viscosity, then the greater the propensity of mud being able to carry or transport cuttings. This is also for the case of concentric position.
The positions of the cuttings in the annulus have no constraints on the fluid transport capacity under a concentric position. This is because the shearing rate at the most extremity is termed pseudo-uniform. This cannot be said the same for eccentric positions because the flow cutting and the dispersion of cuttings are seemly, not uniform.
The ideal situation for optimizing the fluid rheology transport capacity is to keep drill pipe rotation at a favourable speed in a concentric position. Field practice would have the combination of the two positions at their best. For drilling optimization, as purported by previous studies [60,61], the efficiency of transporting cuttings to the surface of the wells is to increase the propitious speed of the drill pipe in an eccentric position, as illustrated in Figure 11. position will adversely create a stuck pipe [59] if the drill pipe rotation for this simulation remains constant; this is because proven fluid viscosity varies, and this is aligned with previous studies delineating that they are non-Newtonian; their shear rates openly lead to cuttings' carriage. For the avoidance of thought and to optimise the fluids carriage rheology, the drill pipe rotation should be maintained constant only at the concentric position, or increased drill pipe rotation at an eccentric position.
The higher the shear region, the lower the viscosity and the lesser the synthetic and oil-based mud transport capacity; this is for the case of an eccentric position. Inversely, the lower the shear region, the higher the viscosity, then the greater the propensity of mud being able to carry or transport cuttings. This is also for the case of concentric position.
The positions of the cuttings in the annulus have no constraints on the fluid transport capacity under a concentric position. This is because the shearing rate at the most extremity is termed pseudo-uniform. This cannot be said the same for eccentric positions because the flow cutting and the dispersion of cuttings are seemly, not uniform.
The ideal situation for optimizing the fluid rheology transport capacity is to keep drill pipe rotation at a favourable speed in a concentric position. Field practice would have the combination of the two positions at their best. For drilling optimization, as purported by previous studies [60,61], the efficiency of transporting cuttings to the surface of the wells is to increase the propitious speed of the drill pipe in an eccentric position, as illustrated in Figure 11.

Effects of Fluid Rheology
The properties of drilling fluid rheology, i.e., yield point and plastic viscosity, at any given point in time must have cuttings transport or carriage capabilities otherwise it tends to lose its significance as a drilling fluid. These listed properties cannot be said the same for water which is less viscous. Fluid rheology under surveillance for any given position Figure 11. Streamlines at an eccentric position with RPM ≥ 80 adapted with permission from Ref. [42]. 2018, Epelle, E.I.

Effects of Fluid Rheology
The properties of drilling fluid rheology, i.e., yield point and plastic viscosity, at any given point in time must have cuttings transport or carriage capabilities otherwise it tends to lose its significance as a drilling fluid. These listed properties cannot be said the same for water which is less viscous. Fluid rheology under surveillance for any given position (concentric and eccentric) demonstrated the capacity to shear and transport cuttings to the surface. Maintaining the drill pipe rotation at 80 RPM, particles in the annulus were seen to have a pseudo-uniform movement in a concentric position rather than in an eccentric position. Further, the vertical well under simulation played a major role in reducing friction resistance, since the walls of the formation are not inclined to challenge the behaviour of the particles.
Both the synthetic and oil-based muds exhibited greater proportionate means of transporting cuttings from the well to the surface.

Conclusions
The novelty of this paper presents the opportunity to simulate experimental results emanating from our previous studies to support the efficiency of the drilling fluids transport index considering some fluid rheology (yield point and plastic viscosity). While defining terms under the Eulerian-Eulerian Model, the CFD analysis emphasized liquid-solid particle dispersion and collisions, their velocity profiles, and their effects on limited annulus space. The simulation revealed that:

•
The CFD simulation validates the experimental results; the symmetrical results purport that the transport index or transport capacity for these two different drilling fluids (synthetic-based and oil-based muds) under surveillance were agreeably efficient and suitable for drilling operations.

•
The smaller particles swiftly occupied the void spaces aiding the lift of cuttings to the surface. Particle diameters of 1 mm were spotted as suitable particle sizes for the enhancement of cuttings carriage. Also, the distribution of cuttings at any point of the annulus was only seen as efficient at the concentric positions. • An 80 RPM pipe rotation maintained throughout the simulation favoured particles in the concentric annulus; movement of these particles was seen as uniformly distributed.

•
The turbulence of the solid particles was caused by the velocity profiles of the smaller particles in question.

Acknowledgments:
We are thankful to Nazarbayev University for the opportunity given to us to further share our work under the competitive grant in the School of Mining and Geosciences. Notwithstanding, we passionately acknowledge all authors cited in this piece of work and we are most grateful for their extensive research work which supports knowledge transfer.

Conflicts of Interest:
We hereby state that there are no known competing financial interests or personal ties that could have influenced the research presented in this study.