Abstract

Fibre suspension has garnered considerable attention in turbulent flows that are used in many industries. Solid particles, such as dust particles, notably affect the turbulent flow field in a rotational frame. In assessing their impacts, the dusty turbulent flow for fibre suspensions needs to be studied in a frame of rotation that can be substantially applied in many industries. This study, therefore, aims to build a theoretical model for the energy motion of dusty turbulent flow of fibre suspensions in a rotational frame. The turbulence momentum equation was considered to formulate the model in presence of dusty fluid rotating flow of fibre suspensions. The newly derived equation was derived in second-order correlation tensors Fi,j, Wi,j, Gi,j, Si,j, Xi,j, Yi,j, Qi,j, and Ri,j at any two points in the flow domain, where the tensors were expressed as space, time, and distance functions. The developed model is a considerable improvement because it takes into account all of the potential influential parameters that could affect the motion of turbulent energy, such as dust particles, suspending particles (fibres), and rotating frame. However, the impact of these parameters on turbulence energy motion must be evaluated in order to assess the performance of turbulence systems utilized in a variety of industries, such as paper manufacturing. The present theoretical development will contribute to open up experimental and numerical research opportunities for the advancement of the industry, science, and technology.

1. Introduction

The flows of turbulent are often found in moving objects, for example, various moving vehicles, which affect flow resistance to the bodies of the objects. Turbulent flow motions define the fluid flow behaviours based on mass conservation, energy, and momentum laws [1]. In the second-order pressure-velocity correlation, Hinze and Uberoi [2] derived the turbulent motion at any two points in the flow domain. The correlations between velocity and pressure fluctuations were taken into account in his study. However, the effect of any potential turbulent motion parameters was not considered. The turbulent flows are greatly affected by the introduction of fibres into the flow. In several parts of the industry, fibre injections into the turbulent flow are used, for example, in environmental engineering, paper manufacturing, chemical engineering, the manufacture of composite materials, and the textile industry. The property of the fibre suspension has a notable effect on product quality. The complexity in turbulent fibre motion arises mainly from the translational and revolutionary dynamics of the particles, which are driven by forces and torques, depending on the form and the orientation of the particles [3]. Over the last thirty years, numerous studies have been undertaken on turbulent flows of fibre suspensions. A preliminary study on turbulent motion of fibre (ice crystals) suspension was conducted by Cho et al. [4]. Atmospheric turbulence associated with their effect on fibre orientation was investigated in their study. The mean fibre orientation was not significantly affected by the atmospheric turbulence. Fluid motion equations were solved with a spectral approach by Lin et al. [5], and the trajectories of fibre were determined on slim body theory. Although a key parameter to define the fibre orientations on a broad scale was found to be the Stokes number, its orientation effect on the fibre is not significant and the ratio of fibre direction has a minimum influence on the orientation distribution of fibre.

Particle motions in the flows of turbulent are a fundamental issue for applications [6], for example, atmospheric clouds, and spray combustion engines. In an accurate understanding between dust particle collisions and growth [79], humid terrestrial clouds, cool stars [10], brown dwarves, and cloud production of planets, the interactions between turbulent and the particles have been considered. The evolution of the particles’ size is dependent on the particle collision rate which can be significantly increased by turbulent motions [11]. In the presence of small-sized particles, Bhatti et al. [12] studied the electro-osmotic flow of non-Newtonian fluid moving in a sinusoidal fashion in a Darcy–Brinkman–Forchheimer medium. The Jeffrey fluid model was used to investigate non-Newtonian effects. The Debye length approximation, Poisson–Boltzmann equation, and ionic Nernst–Planck equation were employed to construct the mathematical model. The effect of parameters on both fluid and particle phases was also discussed. The present effects have been proven to be advantageous in the design and manufacture of microfluidic devices for determining the transport mechanism.

A precise and exact calculation of the collision rate includes an accurate understanding of the collision rate impacts caused by turbulence and preferential or clustering particle concentrations. Several studies [1321] have been carried out on the motion of the dusty turbulent flow. The characteristics of the dust particles [2225] mainly rely on the particle size of the turbulence scale. The particles compete over fluid impacts and collisions on turbulent flow into the neighbourhoods of the particle. In a study conducted by [26], entropy generation and irreversibility processes that occur as a result of partial slip on magnetic dusty fluid caused by a peristaltic wave via a porous conduit were studied. The results demonstrated that the presence of dust particles in the fluid causes the flow to decelerate. In the presence of a strong magnetic field, Mahanthesh [27] investigated the effect of Hall current on the two-phase boundary layer flow of an electrically conducting dusty fluid over a permeable stretched sheet. Using appropriate similarity transformations, the governing equations were reduced from a set of partial differential to ordinary differential equations. The author reported that, in controlling the friction factor on the sheet, the dust particles’ mass concentration can be employed as a control parameter. It was also revealed that suction and injection have opposing effects on the formation of the momentum boundary layer. In the presence of dust particles, the effect of Hall current on the time-dependent flow of a nanofluid was also investigated by Gireesha et al. [28]. The boundary layer approximation notion was used for modeling the governing equations for nanofluid as well as dust phases. Increases in heat transfer rate were found to be proportional to Hall current and unsteadiness.

The motions between a particle of the fluid and the suspended fibre rely on fluid dynamics fundamental of turbulence characteristics with the velocity pressure correlations of fluctuating components. A good number of [15, 2935] studies have been undertaken in rotating systems to analyze the motion of turbulent flows. The mean rotation induces dynamical influences on turbulent flows through the pressure-stream-rate correlation in the transport equations. The rotation produces the auxiliary body force-like Coriolis and centrifugal force in rotating turbulent flows, accompanied by the turbulent structures. Consequently, the momentum transfer mechanism suffers from complexity. However, most research studies did not take into account or address the rotational impacts of the system at different angles, owing to rotational frame complications. Coriolis force plays a critical role in a turbulent flow rotating system [36, 37]. To measure the impact of the Coriolis force, You et al. [30, 35] directed an experimental-based study on the rotating turbulent flow. The velocity for particle images was considered in this study to calculate the flow profile. Because of the Coriolis effect, the flow domains between the leading and trailing sides were found to be different [30]. Coriolis force widened the vortex nearby the front side and intimidated the vortex nearby the trailing end. The force of Coriolis showed major effects in the flow field on the vortices. It was also revealed that not only the secondary flow but also the Coriolis force influences the flow domain which must not be overlooked.

Due to the simplification of laminar flow, it is shown through comprehensive literature surveys that several research studies were undertaken on assessing the behaviour of laminar fibre motion. However, the turbulent fibre motion studies are still inadequate because of the complexities of fibre motion and turbulence. Flow phenomena generated in the flow domain are complex if it includes any type of solid particle such as dust particles into the turbulent fibre flow. The dust particles can range from 1 to 400 μm diameters, with particles more than 100 μm settling down towards the formation source [38]. The overall range of the particle sizes can be classified into three categories: lower than 1 μm, 1–20 μm, and more than 20 μm, which are referred to as ultrafine, fines, and large particles, respectively. However, the smaller particles are comparatively difficult to remove or separate from the airstream because they have a stronger inclination to stay in suspension.

In a frame of rotation, the dusty turbulent flow becomes more complex. However, some mathematical models were developed on the energy motion of turbulent flow in second-order correlation tensors by considering dust particles [13], fibre suspension [39], dust particles with fibre suspensions [40], a rotating frame [32], and dust particles in a frame of rotation [15]. Nevertheless, any of these articles did not consider all of the parameters that can influence turbulence energy motion, such as dust particles, fibre suspension, and rotational frame. The main fact behind this is because of its complexity level in mathematical modeling. As a result, the previous models cannot be considered fully complete and feasible, and they must be improved by taking into account all of these influential parameters in order to make the model feasible. By considering all of these parameters, the present study thus constructs a model for the energy motion of dusty rotating turbulent flow of fibre suspensions. To address the research gaps and improve the previous models, the present study introduced correlation tensors Wi,j, Xi,j, Yi,j, and Ri,j, which represent the correlations between fluid velocities, angular velocities, and fluid velocities, angular velocities produced due to the rotation, and the suspending fluid velocities in the fluid flow domain. The resulting partial differential equation model can be generalized to time-space fractional order, making it valuable to the broader scientific community. Based on a numerical investigation, the DNS (direct numerical simulation) method can be applied to the rotation of dusty rigid fibres with different lengths relative to the neighbouring fluid in the channel flow. Due to its high precision, the spectral element method (SEM) may be used to solve the newly developed equation, and SEM applications are also suggested for future research.

2. Methodology

A turbulent flow domain was considered in a rotational frame to formulate and develop the present model. The fluid flow domain comprised fibre suspensions in the existence of dust particles. The mathematical formulation of the model was directed with the equation of momentum for viscid turbulent incompressible flow. The mathematical consequences were sequentially added to the momentum equation. To find the relationships among the velocities of fluid, fibres, and dust particles in the rotational frame, any two points, C and D, were chosen at a distance r in the turbulent flow field. Point C was regarded as the origin of the system coordinate since no difference exists in representing the relationships of turbulence momentum between the points if any of the two points is assumed as the source (origin). The independent variables were chosen to differentiate the influences between distance and position. The method of averaging was used to express first-, second-, and third-order pressure-velocity correlations between the two points, C and D. These correlations were demonstrated in a second-order partial differential equation by introducing second-order correlation tensors and following averaging method and tensor properties.

3. Model Formulation and Development

The equation of momentum and the continuity for viscid turbulent incompressible flow are represented by [32, 41]with the fluid velocity components of , fluid density , pressure , suspending fluids’ kinematical viscosity , three-dimensional permutation symbol , rotation vector , dissipation of turbulence , position x, and time t.

For the turbulent flow with dust (solid) particles, the following equations are obtained:with dust (solid) particle velocity components , dust particle mass for single-spherical with radius , Stoke’s formula for the drag , frequency dimension, and dust particle density number . It is noted that the term appears in equation (2) due to the solid dust particles existent in the field of turbulent flow [14]. The fluid is presumed to be air and treated as multiphase since this present study comprises fibre suspensions and dust particle interactions into the dusty turbulent flow in the frame of rotation.

In presence of fibre, the equation of turbulent motion (3) giveswith the apparent viscosity of fibre suspension , the turbulent intensity of fibre suspension , tensor strain rate , and fourth- and second-order tensors for fibre orientation and , respectively. The term occurs in equation (5) because of the fibres poured into the turbulent flow domain [42].

For a rotating frame, the turbulence momentum equation is given by [32]where denotes the Coriolis force, where , , and specify unit vector normal to and , angular velocity, and the angle makes between and , respectively.

Therefore, combining equations (5) and (6), a new equation is obtained as follows:

Equation (7) represents the new energy equation of motion for fibre suspended dusty rotating turbulent flow.

Suppose C and D are any two points in the turbulent flow field and assume c and d be the directions along with the points C and D, respectively. Therefore, and can be used as the components of the velocity along with C and D directions. Suppose that is the mean velocity, time independent, and constant throughout the considered domain.

Thus, and .

All the terms have a value, which can be calculated from the equations of at C point and of at D point.

At point C, the turbulence momentum can be written from equation (7) for :

Since, for incompressible fluids, , equation (8) yields

Multiplying equation (9) by ,where is taken as constants for the differential system at C point.

Likewise, the turbulence momentum equations for are obtained at D point:

By using the condition of incompressible fluids, , equation (11) gives

Multiplying equation (12) by , the following is obtained:where are taken as constants for differential systems at D point.

Summing equations (10) and (13),

In order to represent the relationships of turbulence momentum between D and C point, no difference exists if C or D is assumed as the source (origin) of the system coordinate. Point C was regarded as origin here. Suppose the variables are independent and used to differentiate the influences between distance and position, where .

Therefore, the following relations are obtained:

Using these relations and averaging all the terms, equation (15) yields

Equation (16) yields the mean energy motion of dusty turbulent flow for fibre suspensions in a rotational frame, where the motion of fibres of turbulent flow is at average speed with respect to the system coordinate. Due to the constant derivative, the coefficient term was omitted in this equation. The equation consists of the double correlations of velocity , pressure-velocity , and triple-velocity such as , where the terms are located at a certain distance from each other. The correlations of pressure-velocity and form first-order tensor as the pressure is taken as a scalar. Likewise, triple correlations of velocity, and , make third-order tensors. The double correlations are displayed in Figure 1, while the triple correlations are illustrated in Figure 2, respectively, at points C and D, where r is the distance from C to D.

The correlations of the order first, second, and third can be labelled by , and , respectively. Thus, the correlations of pressure-velocity and velocity givewhere denotes the pressure which is not an index such as dummy or , and therefore, the summation should not be used to .

The terms of , , , and form the correlation tensors of second order those can be represented by , , , and , accordingly, while and form the correlation tensor of third order and expressed as and , respectively. Consequently, the followings are obtained:

Using the above correlation terms, equation (16) yieldswhere the terms of correlations used in equation (19) are to the considering points, C and D.

The double pressure-velocity correlations are zero for incompressible and isotropic turbulence, which means

In the field of isotropic turbulence, the invariance condition under the reflection to C is followed by

Thus, the simplification of equation (19) gives

The terms , , , , and form tensors of second order that can be represented by , , , , and , respectively, and expressed as

Therefore, equation (23) yields

This newly derived equation expresses the energy motion of dusty rotating turbulent flow of fibre suspensions in 2nd order tensors.

4. Results and Discussion

The energy motion equation of the dusty turbulent flow of fibre suspension (19) was derived in a rotational frame. This new equation was developed in the 2nd order tensors of pressure-velocity and velocity-velocity correlation at any two points C and D of the fluid flow field. All of these tensors were expressed as coordinates of time, space, and distance. The present model was compared to the relevant existing models to assess its feasibility which is summarized in Table 1.

Dynamic behaviours of the particles in turbulent flow fields are a major concern related to the applications. Fibre suspensions in a fluid flow or gas have been known for some time to reduce the shear stress generated by the motion of the fluid through a solid surface under the conditions of turbulent flow. In absence of dust particles into the flow domain, , and consequently, the newly developed equation (24) becomes

Equation (25) signifies the energy motion for turbulent flow of fibre suspensions in a rotational frame in the tensors of 2nd order. Here, the correlation tensors, , and express the velocity relationships at any two points C and D (Figure 2) in the flow domain. The tensors, and , represent the fluid velocity correlations, and , indicate the correlations between the suspending fibres and fluid velocities, and , describe the velocity correlations between fluid and dust particles while and defines the correlations between fluid velocities, and angular velocities between fluid and fibre formed due to the rotation.

Fibre suspensions in the fluid have been known sometimes to reduce the shear stress generated by the motion of the fluid through a solid surface under the conditions of turbulent flow. If fibres are not available in the fluid flow domain, the apparent viscosity for suspending fluid; equation (25) thus yields the form

Equation (26) expresses the energy motion of turbulent flow in the correlation tensors in a rotating frame.

Owing to the broad applications to both geophysical and engineering fluid mechanics, turbulent flow is commonly researched in the rotational frame, for example, in combustion systems, rotary shafts, rotational heat exchangers, turbine blade passages, centres of the nuclear reactor, rotated heat exchangers, combustion systems, and rotor shafts. In almost all of these complicated systems, a coordinate with nonrotation is usually used to investigate the problems. For these types of rotating systems, the Coriolis effect performs a critical role in the fluid flow field as mentioned earlier. In the nonrotational frame, the Coriolis effects do not produce in the fluid flow field, which means the term due to the Coriolis effect vanishes, and therefore, equation (26) gives

Equation (27) signifies the energy motion for turbulent flow of 2nd-order correlation tensors.

If no energy dissipation generated by turbulence, , therefore, equation (27) yields

Equation (28) describes the turbulence motion in the 2nd-order tensors which was obtained by Hinze and Uberoi [2].

5. Conclusion and Future Directions

A new energy motion equation was derived for dusty rotating turbulent flows of fibre suspensions using the method of averaging. Fibre suspension in rotating turbulent flow was taken into account in presence of solid (dust) particles to develop the present theoretical model in second-order correlation tensors. The turbulent fibre suspension flow was subjected to both average fluid motion and random motion caused by fluctuating fluid velocity components. Comparing the present theoretical model to previous relevant models revealed a good agreement, as seen in Table 1, in comparison to previous models, the present model represents a significant improvement as it considers all of the potential influential parameters that could affect the motion of turbulent energy, for instance, dust particles, suspending particles (fibres), and rotating frame. However, the impact of these parameters on turbulence energy motion needs to be assessed to measure the performance of the turbulence system used in several industries.

The present model should be optimized in future studies to identify the most influential parameters affecting the turbulence systems. However, taking into account the nonlinearity, discrete nature, constrained parameters, and boundary conditions of turbulence systems, it is extremely difficult to design the most suitable model that accurately represents reality. Considering all of these difficulties, the future study can develop a physical model for a turbulent channel flow that depicts the modeled system in its physical form through both experimental and numerical studies. The difference between angular velocities of fluid and fibre at the fibre centre of mass should be calculated by estimating the slip spin, which enables a direct comparison of fibre rotation as well as fibre orientation with the carrier flow’s local parameters. The findings of the future study could provide physical mechanisms underlying fibre rotation in the closed region of the channel walls, emphasizing how particularly fibres spin in relation to the fluid when swept or ejected. Future research in this field can significantly contribute to natural sciences, technologies, and industries.

Nomenclature

:Fluid velocity components (m s−1)
:Fluid density (kg m3)
:Pressure (Pa)
:Kinematical viscosity of suspending fluids (m2 s−1)
:Three-dimensional permutation symbol
:Rotation vector
:Dissipation of turbulence
x:Position (m)
:Dust (solid) particle velocity components (m s−1)
:Dust particle mass (kg)
:Dust particle radius
:Stoke’s formula for the drag
:Frequency dimension
:Dust particle density number
:Apparent viscosity of fibre suspension (kg m−1s−1)
:Turbulent intensity of fibre suspension
:Tensor strain rate
:Second-order tensors for fibre orientation
:Fourth-order tensors for fibre orientation
:Coriolis force
:Unit vector normal to and
:Angular velocity
:Angle makes between and
C, D:Any two points in the turbulent flow field
c, d:Directions along with the points C and D, respectively
, :Components of the velocity along with C and D directions
:Considered independent variables
r:Distance between any two points C and D in the flow filed
t:Time (s)
, :Fluid velocity correlations
, :Correlations between the suspending fibres and fluid velocities
, :Velocity correlations between fluid and dust particles
:Correlations between the velocities of fluid and fibre
:Angular velocities between fluid and fibre formed due to the rotation.

Data Availability

No data were used to support this study.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The author highly expresses his gratitude to Asian University for Women, Chattogram, Bangladesh, for their support to carry out this study.