1. Introduction
In hydrodynamics, the description of a fluid’s behavior depends on its flow regime [
1,
2]. At relatively low Reynolds numbers, a laminar flow regime occurs, described by the Navier–Stokes equations, which are based on Newton’s rheological relation for a linearly viscous fluid. Then, with an increase in the flow velocity and at high Reynolds numbers, the rheology of the moving fluid changes: a turbulent flow regime is reached. There is no conventional uniform rheological relation for this flow regime, although work in this direction has been carried out since the times of Reynolds. Currently, modern turbulence models are based on the concept of turbulent viscosity, which has developed since the first theory of turbulence (L. Prandtl’s mixing length theory), up to the present time in the form of modern differential models of turbulence such as the
and
models, as well as other models [
3]. In addition, with an increase in the Reynolds number, the turbulent flow regime may change to another regime when small additives of polymers are added to the fluid flow, leading to the demonstration of the Toms effect. In this case, a “laminar-type” flow regime arises with its own rheological ratio [
4,
5]. Then, with an increase in the Reynolds number, the Virk flow reaches the limit regime [
6] with its rheology.
All the mentioned flow regimes and all types of rheology can be approximately described using one rheological relation, which is a power-law generalization of Newton’s formula for the flow of a viscous fluid, changing the power values in this ratio. Newton’s formula for the shear stress in a longitudinal flow past a flat surface with the velocity
u, depending on the
y transverse coordinate is as follows:
where
is the density,
is the kinematic viscosity, and the formula can be generalized [
7] to the following form
In this expression
is the non-dimensional coefficient, depending on the
nth power value (
). For
and when
, Equation (
2) leads to Newton’s rheological relation and, as a result, to Poiseuille’s formulas for laminar flow in a pipe. For
and when
, this formula leads to the rheological relation for turbulent fluid flow in a pipe and then to the Blasius formula for the resistance coefficient. Equation (
2) can be generalized to the spatial case of flow and the corresponding rheological relation can be presented in tensor form [
8]. As a result, a system of differential equations can be obtained, similar to the system of Navier–Stokes equations, which also allows for the solution of boundary value problems in fluid mechanics [
7], even in relation to the turbulent flow regime.
2. Steady-State Flow in a Circular Cylindrical Pipe at an Arbitrary Value of the n Power
A circular pipe is the most common hydrodynamic object used for tests of theoretical solutions, since a huge amount of experimental material has been accumulated for the flow within this object, which makes it possible to evaluate the quality of the formulas obtained for resistance coefficients and velocity profiles. For the steady-state flow in a straight circular cylindrical pipe with radius
R,
y is the coordinate measured from the wall (
). Let us introduce non-dimensional variables
(non-dimensional coordinate) and
V (non-dimensional velocity) as follows:
where
is the friction velocity, expressed in terms of shear stress
on the wall:
where
is the pipe diameter and
is the longitudinal pressure drop along the pipe of
l length.
Then, considering notation (
3), the formula (
2) for the flow in a pipe in non-dimensional form can be represented as follows:
This non-negative expression can also be associated with the non-negative value of the shear stress
from the motion equation of the continuous medium in stresses, i.e.,
:
Hence follows the differential equation:
where
is the Reynolds number, calculated from the friction velocity. Then:
The boundary condition for this equation is the no-slip condition: , .
The integration of Equation (
7) with this boundary condition leads to the following expression:
from which the non-dimensional velocity profile takes the form:
The non-dimensional velocity averaged over the pipe’s cross section is calculated using the formula [
9,
10]:
Substitution of expression (
8) for the velocity profile into this formula leads to the expression:
where
has the following form for a flow in the pipe:
This parameter is expressed through hyper-geometric functions; its values for different power values are shown in
Table 1.
For other
n values,
can be found in reference materials. It is useful to note that when
,
. The true Reynolds number
, according to the average velocity
, can be calculated [
2] in terms of the
number using friction velocity:
Then, considering Formula (
9), the following can be written:
Hence, if expressing
inversely in terms of
, one obtains
The resistance coefficient
can be represented in terms of the square of the non-dimensional average velocity as
, which, when considering Equation (
9), leads to a formula for this coefficient at an arbitrary
n value:
Formulas (
9)–(
13) make it possible to describe the velocity field and resistance coefficient for a flow in the pipe for any values of the
nth power.
3. Fluid Flow Development in the Pipe with Increasing Reynolds Number
As the Reynolds number increases on the resistance curve, there is a transition in the values of resistance coefficient from one region to another due to a change in rheology.
Figure 1 shows three schemes for the change in the resistance coefficient
for three different cases of flow description implementations corresponding to different values of the
nth power and the
non-dimensional coefficient, respectively. These diagrams indicate: 1—Hagen–Poiseuille curve
,
; 2—Blasius curve,
,
; 3—curve for
,
, 4—curve for
,
, where
m is the roughness parameter, equal to the ratio of the pipe radius
R to the reduced height of the roughness bumps
k (
).
Figure 1a shows the transition from a laminar flow regime to a turbulent one,
Figure 1b shows this transition to the Prandl–Nikuradze curve, which is approximated by curves 2 and 3 with different values of the
nth power; and
Figure 1c shows the transition of curve 2 to curve 4 with a horizontal region with its own roughness value.
The transition of the laminar flow regime with resistance
to turbulent one, when the Blasius formula takes place
is shown in
Figure 1a. From curve 1 corresponding to the
and
, there is a transition to the Blasius curve 2, which corresponds to
and
. The velocity profile (
9) takes the following form for
:
and it becomes a turbulent profile for
:
In dimensional form, taking into account the notation (
3) after the transition to cylindrical coordinates, when
, Formula (
14) takes the form:
and then (
15) takes the form given below:
The Reynolds number (
value) of the transition originating from the laminar regime to the turbulent one can be calculated as the result of the intersection of the Hagen–Poiseuille and Blasius curves:
from which
. This value turns out to be less than 2300, indicated in the references [
1,
2], which is due to the disregard of the transition region from the laminar regime to the turbulent one. If we assume that this region exists, then it would be advisable to accept that
and understand that there is a gap in the resistance curve at this value of the Reynolds number. Thus, for the case of the flow shown in the diagram in
Figure 1a, one obtains:
We shall not consider this case further.
Figure 1b shows a resistance curve, complying with the Prandl–Nikuradze formula:
which is in good agreement with experimental data [
9,
11], and is approximated by two curves corresponding to
and
. For
and
, the Blasius formula
takes place, for
and
the resistance curve, according to the (
13), corresponds to:
The intersection point of these curves corresponds to the Reynolds number, which can be found according to the relation:
from which we obtain
.
The velocity profile, corresponding to the
, complies with expression (
15); for
it takes the following form, according to formula (
8):
Figure 1c shows resistance curves for the flow in rough pipes with different values of the roughness parameter
. For a flow in rough pipes, one can take
, which gives an almost horizontal line in coordinates as follows: (Logarithm of Reynolds number − logarithm of resistance coefficient), i.e., as (
).The rheological power relation (
2) can be presented as follows for
:
Hence, the velocity profile in non-dimensional coordinates, according to (
8), can be presented in the form:
The average non-dimensional velocity value is
taking into account the expression for
, it takes the following form:
For
values of
equal to 11.6, 14.4, and 20.2 can be obtained, which is in good agreement with experimental data [
10]. The connection between Reynolds numbers for average and dynamic velocities for flow in rough pipes is as follows:
The resistance coefficient, according to (
13), is equal to
, and can be rounded to the approximate value in this case:
which is in satisfactory agreement with the experimental data.
The region of a hydraulically smooth pipe in
Figure 1c transitions to the area with a constant value
, depending on the value of the roughness parameter
m. The corresponding Reynolds number
is found as a result of the intersection of the corresponding resistance curves:
which leads to
The Reynolds number
values are presented in
Table 2.
Resistance curves
, when the quadratic law of resistance is satisfied, may have different positions according to each other. As a result, as the flow develops in rough pipes, and the
nth power and
coefficient change:
Thus, with an increase in the values of the Reynolds number, the
nth power changes, which leads to changes in the velocity profiles and resistance coefficients according to formulas (
8) and (
13), respectively.
4. The Flow of Weakly Concentrated Aqueous Solutions of Polymers
In the flow of weakly concentrated aqueous solutions of polymers, in which the Toms effect is observed [
4,
5], the resistance curves have the typical form shown in
Figure 2.
It can be seen in
Figure 2 that the Toms effect is associated with the deviation of the resistance curve from the Prandtl–Nikuradze law with access to an area, which is equidistant to the resistance curve
for the laminar flow regime.The case arising with the resistance curve can be schematically represented in
Figure 3, which shows the “laminar-type” regions of the resistance curves during the flow of polymer solutions, where the Toms effect is observed.
The resistance curve for the flow of polymer solutions, immediately after deviating from the resistance law for a hydraulically smooth pipe, becomes similar to the resistance curve
(laminar regime). Therefore, this resistance law for flows of polymer solutions in the region after deviation from the resistance law for a hydraulically smooth pipe can be represented [
12] in the form:
where
A takes its value depending on the kind of polymer, its concentration, and the radius of the pipe. According to
Figure 3, it is possible to estimate the values of this constant for the flow in pipes of different diameters at a constant volumetric concentration of the polymer solution
(
Table 3).
Then, with an increase in the Reynolds number, there is a deviation from the “laminar-type” resistance curve with a tendency to reach the Virk limit curve [
11]. However, before beginning of the deviation from the resistance law for a hydraulically smooth pipe, the turbulent flow behaves in accordance with the flow laws of ordinary viscous Newtonian fluids (curve 2 in
Figure 2).
With the development of the flow of weakly concentrated aqueous solutions of polymers, the resistance curve consists of several regions, shown schematically in
Figure 4:
The laminar a–1 region (beginning of the flow), which complies with the formula ;
The turbulent 1–2 region with an increase in the flow velocity, which complies with the Blasius formula ;
The “laminar-type” 2–3 region with a further increase in speed, which complies with the formula ;
The limit 3–
b region, which is represented by the Virk curve [
6] for the final stage of the flow, which complies with the formula
.
The sequence of the process of the resistance change during the flow of a fluid with additives of polymers in a pipe is considered in detail as follows:
- 1.
In the region, there is a laminar flow regime. Disregarding the transition region from the laminar regime to the turbulent one, we can approximately set the value as ;
- 2.
In the 1–2 region, where
, a turbulent flow regime develops in a hydraulically smooth pipe. The Reynolds number
of the end of this flow regime is determined as the result of the intersection of the corresponding resistance curves for the Blasius law and the “laminar-type” resistance law:
, and then
- 3.
In the “laminar-type” 2–3 region, where
, there is a deviation from the resistance curve of a hydraulically smooth pipe to the resistance curve for a polymer solution flow with a resistance law
. The value of the “
A” constant determines the type of polymer, its concentration in a solution and the pipe’s diameter [
4,
5]. The
value can be found according to the intersection conditions of the “laminar-type” curve
with the Virk limit curve
, i.e.,
, and then
for different “
A” parameter values; the values of
and
are given in
Table 4;
- 4.
With a further increase in the Reynolds number when , there is a transition to the Virk resistance curve, complying with the law .
Thus, with the development of the flow of weakly concentrated aqueous solutions of polymers, each region has its own values of the nth power, non-dimensional parameters and :
- 1.
In region a–1, when , there are , ,
- 2.
In region 1–2, when , there are , ,
- 3.
In region 2–3, when , there are , ,
- 4.
In region 3–b, when , there are , ,
In accordance with these values, for each region expressions for the resistance coefficients can be written according to (
13), velocity profiles can be written in terms of universal coordinates according to (
3), and Reynolds numbers can be written in terms of friction velocity. These expressions for each region are as follows:
For the flow of polymer solutions in rough pipes with a roughness parameter
m, the Reynolds number
determines the beginning of the transition from the Blasius curve to the region with a constant resistance coefficient, according to (
27). If
, then there will be no transition to the flow regime of the “laminar-type” region and polymer additives will have no effect on the flow (
Figure 4). If it turns out that
, then the roughness does not have time to influence the flow, and the resistance curve during the flow of polymer solutions will contain the same four regions (
29)–(
32) as a hydraulically smooth pipe.
For the resistance coefficients at points 1, 2, and 3 in
Figure 4, there is a transition from one flow regime to another. Moreover, this transition, within the framework of this paper, is shown to be an abrupt one. In fact, there is always a smooth change in the resistance coefficients, and there are always transient flow regimes. The same goes for velocity profiles that transition smoothly from one shape to another. However, these transient processes occur in relatively narrow ranges of Reynolds numbers and in the first approximation they can be disregarded.
Let us change the velocity profiles at points 1, 2 and 3 in
Figure 4 in universal coordinates, which can be considered as follows:
where
y is the coordinate measured from the pipe wall,
is friction velocity,
is kinematic viscosity, and
u is the Reynolds-averaged longitudinal velocity. In this case, the non-dimensional coordinate
, appearing in the expressions for the velocity profiles, will be represented in terms of the Reynolds number for the friction velocity as
. At point 1 at the value of the Reynolds number
by the average velocity, the laminar flow regime (
) changes to a turbulent one (
), and the velocity profile (
14) becomes a turbulent one (
15). So, for
and
, the Reynolds number in terms of friction velocity according to formula (
12) takes the same value, equal to
(for this transitional Reynolds number
). These profiles are shown in
Figure 5.
It can thus be seen how the laminar profile 1 changes dramatically to the turbulent profile 2 when
. The corresponding value of the non-dimensional average velocity according to (
9) for both profiles will be the same:
. Then, with an increase in the Reynolds number, there are power-law profiles 3 and 4 of the turbulent flow regime (
) according to (
15) for
and
, respectively.
When small additives of polymers, for which the Toms effect is observed, are introduced into the flow, point 2 appears on the resistance curve shown in
Figure 4, in which the turbulent velocity profile (
) changes to a “laminar-type” one (
) with the corresponding
A value for this polymer.
Figure 6 at this point shows the change in the turbulent profile (
15) to the profile (
31) when
for
. In both cases for this Reynolds number, the dynamic Reynolds number is the same and equal to
.
It can be seen that profile 1 is abrupt at
and transitions to profile 2. The corresponding value of the dimensionless average velocity according to (
9) for both profiles will be the same, and equal to
. Then, with an increase in the Reynolds number, there are “laminar-type” profiles 3, 4 and 5, according to (
31) for
= 650, 100, and 1625, respectively.
At point 3, the “laminar-type” flow regime (
,
) for
changes to the Virk flow regime (
), and the velocity profile (
31) changes to profile (
32). For both profiles at this Reynolds number, the Reynolds number in terms of the friction velocity according to formula (
12) takes the same value,
. These profiles are shown in
Figure 7.
It can be seen that the “laminar-type” velocity profile 1 abruptly changes to the Virk profile 2. The corresponding value of the non-dimensional average velocity according to (
9) for both profiles will be the same and equal to
. Then, with an increase in the Reynolds number, one obtains Virk profiles 3 and 4, according to (
31) for
and
, respectively.
5. Conclusions
In the present research, a power-law generalization of Newton’s formula for the flow of a viscous fluid is used, in which several variants of rheological relations arise when the value of the exponent changes to be in accordance with the description of various fluid flow regimes, including non-Newtonian ones. This generalization, followed by expressing the corresponding rheological relation in tensor form and substituting it into the equation of motion of a continuous medium under stresses, leads to a system of differential equations of fluid motion, which can be used to describe flows with different rheologies and in different flow regimes (in turbulent ones as well). It should be noted that there are no exact analytical solutions to the problems of simple shear flows of an incompressible fluid in modern papers on the turbulent regime. There are no differential equations for this flow regime in the known papers that allow us to consider the problems of calculating turbulent flows as problems of mathematical physics. By using a system of equations (obtained on the basis of the considered power-law generalization of Newton’s formula), it makes it possible to interpret the problems of calculating turbulent flows as problems of mathematical physics and obtain accurate analytical solutions to problems of simple shear flows. This system can be used in studies of both plane-parallel and three-dimensional flows, which was demonstrated to a certain extent in [
7], which contains a solution for the problem of longitudinal flow around a flat plate and provides estimates of the boundary layer type. In this paper, a system of equations for a 2D fluid flow in the boundary layer of a flat plate is proposed and this system is reduced to one ordinary third-order equation, similarly to the work of Blasius in relation to the laminar boundary layer. The method of the direct reduction of the boundary value problem to the Cauchy problem was used to solve this equation. The results of this solution made it possible to determine the expressions for the resistance coefficients, as well as the thicknesses of the boundary layer, displacement, and momentum loss. These values are comparable with the available experimental data.
The present research shows that the fluid flow development in a pipe and all the considered flow regimes within it, along with their rheology, can be described by means of one rheological relation, which is a power-law two-parameter generalization of Newton’s formula for a viscous fluid. In this case, the change in the flow movement characteristics can be described by changing the values of the nth power and the coefficient in this ratio when the corresponding critical Reynolds numbers are reached. This generalization makes it possible to research the behavior of flows with small concentrations of polymers, in which the Toms effect is observed. For such flows, the turbulent regime changes with an increase in the Reynolds number to another one, leading to “laminar-type” velocity profiles and the corresponding resistance laws. Then, the Virk flow reaches its limit regime with its own rheology with an increase in the Reynolds number. The results of this research allow us to carry out practical calculations of the velocity profiles and resistance coefficients of flows with arbitrary rheology based on a unified approach.
Subsequently, we plan to compare the proposed model with modern approaches such as DNS, and also perform calculations for more complex geometries on the basis of the proposed approach.