Next Article in Journal
DFT Study of Adsorption Behavior of Nitro Species on Carbon-Doped Boron Nitride Nanoribbons for Toxic Gas Sensing
Next Article in Special Issue
Instability of Liquid Film with Odd Viscosity over a Non-Uniformly Heated and Corrugated Substrate
Previous Article in Journal
Bifunctional Bicarbazole-Benzophenone-Based Twisted Donor–Acceptor–Donor Derivatives for Deep-Blue and Green OLEDs
Previous Article in Special Issue
Influence of Stress Jump Condition at the Interface Region of a Two-Layer Nanofluid Flow in a Microchannel with EDL Effects
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Exact Solutions for Non-Isothermal Flows of Second Grade Fluid between Parallel Plates

by
Evgenii S. Baranovskii
Department of Applied Mathematics, Informatics and Mechanics, Voronezh State University, 394018 Voronezh, Russia
Nanomaterials 2023, 13(8), 1409; https://doi.org/10.3390/nano13081409
Submission received: 20 March 2023 / Revised: 12 April 2023 / Accepted: 17 April 2023 / Published: 19 April 2023
(This article belongs to the Special Issue Advances of Nanoscale Fluid Mechanics)

Abstract

:
In this paper, we obtain new exact solutions for the unidirectional non-isothermal flow of a second grade fluid in a plane channel with impermeable solid walls, taking into account the fluid energy dissipation (mechanical-to-thermal energy conversion) in the heat transfer equation. It is assumed that the flow is time-independent and driven by the pressure gradient. On the channel walls, various boundary conditions are stated. Namely, we consider the no-slip conditions, the threshold slip conditions, which include Navier’s slip condition (free slip) as a limit case, as well as mixed boundary conditions, assuming that the upper and lower walls of the channel differ in their physical properties. The dependence of solutions on the boundary conditions is discussed in some detail. Moreover, we establish explicit relationships for the model parameters that guarantee the slip (or no-slip) regime on the boundaries.

1. Introduction

Many real fluids and fluid-like materials used in nanotechnologies belong to the class of fluids of complexity N (see [1,2]). For these fluids, the Cauchy stress tensor T is given by the relation
T = p I + F ( A 1 , , A N ) ,
where
  • p is the pressure;
  • I is the identity tensor;
  • F is a frame indifferent response function;
  • A 1 , , A N are the first N Rivlin–Ericksen tensors:
    A 1 = def v + ( v ) , A j = def d d t A j 1 + A j 1 v + ( v ) A j 1 , j = 2 , , N ;
  • v is the velocity field;
  • v denotes the velocity gradient;
  • ( v ) denotes the transpose of the velocity gradient;
  • the differential operator d / d t is the material time derivative,
    d d t A j 1 = def t A j 1 + ( v · ) A j 1 .
If F is a polynomial of degree N, then the corresponding fluid is called a fluid of grade N.
An incompressible Newtonian fluid
T = p I + μ A 1 , μ > 0 ,
is a fluid of grade 1. Fluids with shear-dependent viscosity, for which the constitutive equation is given by the equality
T = p I + μ ( A 1 ) A 1 ,
belong to the class of fluids of complexity 1.
In the present paper, we deal with the second grade fluids:
T = p I + μ A 1 + α 1 A 2 + α 2 A 1 2 ,
where μ is the viscosity coefficient, μ > 0 , while α 1 and α 2 are the normal stress moduli.
Note that if the equality α 1 + α 2 = 0 holds, then one can rewrite (1) as follows:
T = p I + μ A 1 + α 1 d d t A 1 + α 1 ( A 1 W W A 1 ) ,
where W is the vorticity tensor defined by
W = def 1 2 ( v ( v ) ) .
The nonlinear constitutive relations (1) and (2), as well as their various modifications, are often used in the dynamics modeling of nanoscale fluids (see, for example [3,4,5,6]).
Many dilute polymer solutions belong to the class of nanofluids that obey (1). It is well known that the addition of a small amount of polymer nanoparticles to water almost does not change the physical characteristics of the solution, such as the density and the viscosity, but the fluid gains some relaxation properties. An important consequence is that the friction drag drastically decreases for both internal problems (flow in pipes) and external problems (flow past bodies). This effect was discovered by Toms [7] and stimulated in a series of experimental and theoretical studies of the dynamics of aqueous solutions of polymers (see [8,9,10,11,12,13,14,15] and the references therein).
A model for the motion of polymer solutions, considering their relaxation properties, was proposed by Voitkunskii, Amfilokhiev, and Pavlovskii [16]. Using ideas of the hereditary theory of viscoelasticity [2,17], these authors introduced the Maxwell-type relationship between the Cauchy stress tensor T and the deformation rate tensor D :
T = p I + 2 μ D ( v ( · , t ) ) + 2 χ τ t exp t s τ d d s D ( v ( · , s ) ) d s ,
where
  • D ( v ) = def ( v + ( v ) ) / 2 ;
  • τ is the shear stress relaxation time, τ > 0 ;
  • χ is the relaxation viscosity coefficient, χ > 0 .
Using the smallness of the parameter τ , Pavlovskii [18] performed the asymptotic expansion of the integral term from (3) with respect to τ 0 + . Retaining only the first term of this expansion, he obtained
T = p I + 2 μ D + 2 χ d d t D ( v ) .
Clearly, the last relation is a simplified version of (2) with α 1 = χ under the assumption that the product α 1 ( A 1 W W A 1 ) is small compared to the other terms in equality (2) and can be dropped.
The adequacy of the rheological model (4) has been supported by experimental studies. In particular, (4) is considered as a suitable constitutive relationship for low-concentrated aqueous solutions of polyethylenoxide, polyacrylamide, and guar gum [10,11]. The analysis of exact solutions of the corresponding nonlinear motion equations confirms that polymer nanoparticles added to water even in small amounts have a significant influence on the flow pattern [19].
A fluid modeled by (1) is compatible with the thermodynamic laws and stability principles if the following restrictions are imposed on the material constants α 1 and α 2 :
α 1 0 , α 1 + α 2 = 0 ,
for details, see [20]. Moreover, Fosdick and Rajagopal [21] showed that for arbitrary values of the sum α 1 + α 2 , with α 1 < 0 , a fluid totally filling a bounded domain and adhering to the boundary of this domain exhibits an anomalous behavior not expected with real fluids. For a detailed discussion on the physical background, we refer readers to the critical and extensive historical review of second grade (and higher-order) fluid models [22].
Assuming that (5) holds, we introduce the notation α = α 1 = α 2 and rewrite (1) as follows:
T = p I + μ A 1 + α A 2 α A 1 2 .
The aim of the present paper is to obtain exact solutions for non-isothermal steady-state flows of the fluid (6) in a flat infinite channel with impermeable solid walls.
The pointed feature of our work is that different types of boundary conditions on channel walls are used. In addition to the standard no-slip boundary condition v = 0 , we will consider the threshold slip conditions, which include Navier’s slip conditions as a limit case, as well as mixed boundary conditions, which are suitable for the case when the upper and lower walls have different physical properties. Importance of the wall slip effect and its influence on various characteristics of fluid flows, especially in the case of non-Newtonian fluids, are mentioned in many studies (see [23,24,25,26,27] and the references therein). In particular, as noted in [27], the study of wall slip is very important because it can be used to determine the true rheology of complex fluids by correcting data for slipping effects and explaining a mismatch of rheological data that are obtained from rheometers having different geometries.
Another feature of this paper is that we take into account the fluid energy dissipation (mechanical-to-thermal energy conversion) in the heat transfer equation. In many studies (see, e.g., [28,29]), the influence of the Rayleigh dissipation function is neglected because the mathematical analysis of heat and motion equations are considerably simplified due to artificially vanishing the term involving a quadratic function of space derivatives of the velocity field. However, from the physical point of view, it is more interesting not to use this simplifying assumption and keep all nonlinearities in the original equations [30,31,32].
For each boundary value problem under consideration, we construct exact solutions which determine the velocity field, the temperature, and the pressure in the flow region. Since the used boundary conditions allow for various types of “fluid–solid walls” interactions, we establish explicit relationships for model parameters that guarantee the slip/no-slip regime on the channel walls. Note that the obtained results are valid and new for a Newtonian fluid too, which can be considered as the limit case of a second grade fluid as α 0 + .
The present paper is a continuation of [33,34], in which analogous boundary value problems were considered for isothermal flows. It should be mentioned that many exact solutions for steady and time-dependent motions of the second grade fluids have been established by different authors. The first exact solutions for unsteady flows of this class of non-Newtonian fluids seem to be those of Ting [35], both in rectangular and cylindrical domains. In particular, he showed that solutions are unbounded when α 1 < 0 . Ting’s results were extended by Coleman et al. [36], who performed a mathematical analysis (instability, uniqueness, and nonexistence theorems) of initial boundary value problems describing non-steady simple shearing flows of second grade fluids provided that α 1 < 0 . Hron et al. [37] investigated exact solutions for steady-state flows of fluids of complexity 2 in a plane channel and a cylindrical pipe and flows between two rotating concentric cylinders subject to Navier’s slip boundary condition. Exact solutions for the velocity field corresponding to the second problem of Stokes were obtained in [38] by the Laplace transform method. Fetecau et al. [39] analytically studied the magnetohydrodynamic (MHD) flow of second grade fluids with Caputo–Fabrizio time fractional derivatives over a moving infinite flat plate. In [40], it was shown that the governing equations for the fluid velocity and non-trivial shear stress corresponding to some isothermal MHD unidirectional motions of second grade fluids through a porous medium have identical forms. Fetecau and Vieru [41] provided the first exact general solutions for isothermal MHD flows of incompressible second grade fluids between infinite horizontal parallel plates embedded in a porous medium. Note also that there are numerous mathematical studies concerning the existence and uniqueness of solutions to the motion equations of second grade fluids [42,43,44,45,46] as well as optimal control flow problems [47,48,49,50]. The literature on these fluids continues to grow, providing a deeper understanding of the physical processes and support for modern technological advances, in particular for nanotechnologies.

2. Statements of Boundary Value Problems

The non-isothermal steady flow of a fluid with constant density is governed by the following system of equations:
ρ ( v · ) v = div T + ρ g ,
· v = 0 ,
ρ ( v · ) θ κ 2 θ = ω + Φ ,
where
  • ρ is the fluid density, ρ > 0 ;
  • v = ( v 1 , v 2 , v 3 ) is the velocity vector;
  • T is the Cauchy stress tensor;
  • g = ( g 1 , g 2 , g 3 ) is the external force per unit mass;
  • θ is the temperature;
  • κ is the thermal conductivity, κ > 0 ;
  • ω is the heat source intensity;
  • Φ is the Rayleigh function that determines the fluid energy dissipation (mechanical-to-thermal energy conversion) according to the formula
    Φ = 1 2 c p T : A 1 ;
  • the colon symbol: denotes the scalar product of tensors;
  • c p is the heat capacity of the fluid, c p > 0 ;
  • the operators div and ∇ are the divergence and the gradient, respectively, with respect to the space variables x, y, z;
  • 2 = def · = 2 x 2 + 2 y 2 + 2 z 2 .
Let us consider the unidirectional fluid motion between horizontal parallel plates y = h and y = h , assuming that the flow is driven by a constant pressure gradient
p x = ξ , ξ = const , ξ > 0 ,
and
g = ( 0 , g , 0 ) , θ = θ ( y ) , ω = 0 ,
where g is the value of acceleration due to gravity. This means that we deal with the plane Poiseuille flow. Figure 1 shows the chosen coordinate system and the flow geometry.
For such flow, we obviously have
v 1 = u , v 2 = 0 , v 3 = 0 ,
where u = u ( y ) is an unknown function. Then the following equalities hold:
· v = 0 , ( v · ) v = 0 , ( v · ) θ = 0 .
In view of relations (12) and (13), system (7)–(9) reduces to
div T + ρ g = 0 ,
κ θ = Φ ,
where the symbol denotes the differentiation with respect to y.
Assuming that the fluid obeys the constitutive relation (6), we rewrite (14) in the form
div ( μ A 1 + α A 2 α A 1 2 ) = p ρ g .
We will use the nonlinear system (15), (16) for handling second grade fluid flows in the channel h y h . Note that the unknowns of this system are u, p, and θ , while all other quantities are assumed to be given.
Of course, in order to obtain physically important solutions, Equations (15) and (16) must be supplemented with appropriate boundary conditions for the velocity field and the temperature. Experimental data and theoretical works point to different possibilities for the behaviour of fluid flows on solid walls. Along with the standard no-slip condition, various slip conditions are widely used (see, e.g., [23], § 5).
In this paper, we will investigate four boundary value problems describing flows of second grade fluids in the plane channel with impermeable solid walls.
Problem 1. 
Find a triplet ( u , p , θ ) that satisfies system (15), (16) supplemented with the no-slip boundary condition
v ( ± h ) = 0
and the Robin boundary condition for the temperature θ
κ θ ( ± h ) = β θ ( ± h ) ,
where β is a positive coefficient that characterizes the heat transfer on the channel walls.
Problem 2. 
Find a triplet ( u , p , θ ) that satisfies system (15), (16) supplemented with the threshold slip conditions on the plates y = ± h :
v · n = 0 ,
( T n ) tan R 3 σ v tan = 0 ,
( T n ) tan R 3 > σ ( T n ) tan = ( σ + k v tan R 3 ) v tan v tan R 3
and boundary condition (18) for the temperature θ.
Here, and in the succeeding discussion, the following notations are used:
  • n is the exterior unit normal vector on the channel walls;
  • v tan denotes the tangential component of v ;
  • k is the slip coefficient, k > 0 ;
  • σ is the threshold value of the tangential stresses, σ 0 .
Equality (19) represents the impermeability condition on the channel walls. Relations (20) and (21) mean that the fluid slips at a point on the boundary if and only if the magnitude of the tangential traction exceeds the slip threshold σ . These conditions are called the threshold slip conditions as well as the Navier–Fujita slip conditions [51].
Problem 3. 
Find a triplet ( u , p , θ ) that satisfies system (15), (16) under the mixed conditions for the velocity field v and the temperature θ:
v · n = 0 on   the   plate   y = h ,
( T n ) tan R 3 σ v tan = 0 on   the   plate   y = h ,
( T n ) tan R 3 > σ ( T n ) tan = ( σ + k v tan R 3 ) v tan v tan R 3 on   the   plate   y = h ,
v = 0 on   the   plate   y = h ,
θ = 0 on   the   plate   y = h ,
θ = 0 on   the   plate   y = h .
Problem 4. 
Find a triplet ( u , p , θ ) that satisfies system (15), (16) under the mixed conditions for the velocity field v and the temperature θ:
v · n = 0 on   the   plates   y = ± h ,
( T n ) tan R 3 σ v tan = 0 on   the   plate   y = h ,
( T n ) tan R 3 > σ ( T n ) tan = ( σ + k 1 v tan R 3 ) v tan v tan R 3 on   the   plate   y = h ,
( T n ) tan = k 2 v tan on   the   plate   y = h ,
θ = 0 on   the   plates   y = ± h ,
where k 1 > 0 and k 2 0 .
Note that condition (30) states that the fluid slips on the solid wall for any non-zero shear stresses. This situation corresponds to the limit case as σ 0 + for the threshold slip conditions. In the literature, equality (30) is referred to as the Navier slip condition, after Navier [52] who first proposed it. The corresponding slip regime is sometimes referred to as the free slip [23]. However, this condition should not be confused with the perfect slip condition  ( T n ) tan = 0 (see [53,54,55]), which is valid only when k 2 = 0 . As noted in [37], Navier’s slip condition can be considered as a homotopy transformation that links the no-slip boundary condition on the one hand with the no-stick boundary condition on the other hand.

3. Analysis and Exact Solution of Problem 1

First we calculate the Rivlin–Ericksen tensors A 1 and A 2 :
A 1 = 0 u ( y ) 0 u ( y ) 0 0 0 0 0 , A 2 = 0 0 0 0 2 ( u ( y ) ) 2 0 0 0 0 .
Next, using these equalities and (6), we obtain
A 1 2 = ( u ( y ) ) 2 0 0 0 ( u ( y ) ) 2 0 0 0 0 ,
T = p ( x , y , z ) I + α ( u ( y ) ) 2 μ u ( y ) 0 μ u ( y ) α ( u ( y ) ) 2 0 0 0 0
and rewrite (16) in the form
div α ( u ( y ) ) 2 μ u ( y ) 0 μ u ( y ) α ( u ( y ) ) 2 0 0 0 0 = p ( x , y , z ) ρ g .
The last equation is equivalent to the following system:
μ u ( y ) = p ( x , y , z ) x ,
α [ ( u ( y ) ) 2 ] = p ( x , y , z ) y + ρ g ,
0 = p ( x , y , z ) z .
From (10) and the first equality of (32) it follows that
Φ = μ c p ( u ( y ) ) 2 .
Using this equality, we rewrite (15) as follows:
κ θ ( y ) = μ c p ( u ( y ) ) 2 .
Note that system (34)–(37) can be considered as a starting point for solving all the boundary value problems that are stated in this paper.
Let us construct the exact solution to Problem 1. First, we will find the pressure p. In view of (36), the pressure is independent of z, that is, p = p ( x , y ) . Moreover, taking into account condition (11), we conclude that the pressure should be sought in the form
p ( x , y ) = ξ x + ϕ ( y )
with an unknown function ϕ = ϕ ( y ) .
From (35) it follows that
ϕ ( y ) = α ( u ( y ) ) 2 ρ g y + C ,
where C is a constant. By setting C = ρ g h , we obtain
ϕ ( y ) = α ( u ( y ) ) 2 + ρ g ( h y ) .
From (11) and (34) it follows that
μ u ( y ) = ξ y + C 0 ,
where C 0 is a constant.
In view of the physical meaning of Problem 1, the velocity field is symmetric with respect to the plane y = 0 , that is, the function u = u ( y ) is even. Hence, we have u ( 0 ) = 0 . Setting y = 0 in (39), we obtain C 0 = 0 . Therefore,
u ( y ) = ξ μ y .
Substituting (40) into (38), we arrive at
ϕ ( y ) = α ξ 2 μ 2 y 2 + ρ g ( h y ) .
Next, we shall find u. From (40) it follows that
u ( y ) = ξ 2 μ y 2 + C 1 .
It is clear that the value of the constant C 1 must be chosen such that the no-slip condition (17) on the channel walls is satisfied. Since u is an even function, it suffices to verify that the boundary condition holds on the upper wall. Setting y = h in (41), we find C 1 = ξ h 2 / ( 2 μ ) .
We now know the function u and hence, in order to find the temperature distribution in the channel, we can solve the boundary value problem (18), (37) with respect to θ .
Thus, we have obtained the exact solution to Problem 1:
u ( y ) = ξ 2 μ ( y 2 h 2 ) , p ( x , y ) = ξ x + α ξ 2 μ 2 y 2 + ρ g ( h y ) , θ ( y ) = ξ 2 12 c p μ κ y 4 + ξ 2 h 3 ( h β + 4 κ ) 12 c p μ β κ .

4. Analysis and Exact Solution of Problem 2

Now let us consider the case of the threshold slip boundary conditions (19)–(21).
With the help of arguments similar to those given in the previous section, we can verify the validity of relations (38), (40), and (41) for Problem 2.
Next, directly calculating
T n | y = ± h = p ( x , ± h ) 1 0 0 0 1 0 0 0 1 0 ± 1 0 + α ( u ( ± h ) ) 2 μ u ( ± h ) 0 μ u ( ± h ) α ( u ( ± h ) ) 2 0 0 0 0 0 ± 1 0 = ± μ u ( ± h ) p ( x , ± h ) ± α ( u ( ± h ) ) 2 0 ,
and using (40), we derive
( T n | y = ± h ) tan R 3 = | μ u ( ± h ) | = ξ h .
Taking into account the last relation and (19)–(21), we separately consider two cases: ξ h σ and ξ h > σ .
If the inequality ξ h σ is valid, then in view of equality (20), there is no boundary slip and the velocity field in the channel is determined by formula (42) as in Problem 1.
If the inequality ξ h > σ holds, then the fluid velocity on the channel walls is non-zero. In view of (21), the following equalities hold:
± μ u ( ± h ) = ( σ + k u ( ± h ) ) .
Taking into account (40), we obtain
ξ h = ( σ + k u ( ± h ) ) ,
and hence
u ( ± h ) = ξ h σ k .
Using this equality and (41), we find the function u:
u ( y ) = ξ 2 μ ( y 2 h 2 ) + ξ h σ k .
In both cases, to obtain the temperature distribution in the channel, it is sufficient to solve equation (37) under boundary condition (18) with respect to the function θ .
Combining the solutions constructed for each of the above-mentioned cases, we obtain the general solution to Problem 2, which satisfies all imposed conditions:
u ( y ) = ξ 2 μ ( y 2 h 2 ) + ( 1 H ( σ ξ h ) ) ( ξ h σ ) k , p ( x , y ) = ξ x + α ξ 2 μ 2 y 2 + ρ g ( h y ) , θ ( y ) = ξ 2 12 c p μ κ y 4 + ξ 2 h 3 ( h β + 4 κ ) 12 c p μ β κ ,
where H is the Heaviside step function defined by
H ( s ) = def 0 if s < 0 , 1 if s 0 .

5. Analysis and Exact Solution of Problem 3

For flow models with mixed boundary conditions on the channel walls, one must keep in mind that the velocity field is not symmetric with respect to the plane y = 0 , and hence relation (40) may not hold. Therefore, we turn to relation (34), from which, after integrating with respect to y, we find
μ u ( y ) = ξ 2 y 2 + C 1 y + C 2 ,
where C 1 and C 2 are some constants.
Using boundary condition (24), we derive from (43) that
C 2 = ξ 2 h 2 + C 1 h .
Substituting C 2 into (43), we arrive at the relation
μ u ( y ) = ξ 2 ( y 2 h 2 ) + C 1 ( y + h ) .
Now let us find the value of the constant C 1 based on the threshold slip condition on the upper wall of the channel. Differentiating both sides of identity (43) with respect to y, we obtain
μ u ( y ) = ξ y + C 1 .
Taking into account (33) and (45), we derive
( T n | y = h ) tan = μ u ( h ) 0 0 = ξ h + C 1 0 0 .
Let us consider separately the two cases: the no-slip regime and the slip regime on the wall y = h .
If the fluid adheres to the wall y = h (mathematically, this means that u ( h ) = 0 ), then from (44) it follows that C 1 = 0 . In view of condition (22), this regime is realized if
( T n | y = h ) tan R 3 = ξ h σ .
The slip regime arises if
( T n | y = h ) tan R 3 = | ξ h + C 1 | > σ .
In view of condition (23), the following equality holds:
μ u ( h ) = ( σ + k | u ( h ) | ) s g n ( u ( h ) ) .
Using (44) and (45), we rewrite the last equality as follows:
ξ h + C 1 = σ + 2 k h | C 1 | μ s g n ( C 1 ) .
This implies, in particular, that C 1 > 0 . Therefore, (47) reduces to
ξ h + C 1 = σ + 2 k h C 1 μ ,
where C 1 = μ ( ξ h σ ) / ( μ + 2 h k ) . Combining this with C 1 > 0 , we arrive at the inequality ξ h > σ indicating the slip regime.
After finding the velocity component u, one can derive the temperature θ from system (15), (25), and (26).
Thus, we have obtained the general solution of Problem 3, which is suitable for any admissible values of the model parameters:
u ( y ) = ξ 2 μ ( y 2 h 2 ) + ψ ( y + h ) , p ( x , y ) = ξ x + α μ 2 ( ξ y ψ μ ) 2 + ρ g ( h y ) , θ ( y ) = ξ y ψ μ 4 12 c p μ κ ξ 2 ( ψ μ + h ξ ) 3 3 c p μ κ ξ y + ψ 4 μ 4 + 18 ψ 2 h 2 μ 2 ξ 2 + 8 ψ h 3 μ ξ 3 + 5 h 4 ξ 4 12 c p μ κ ξ 2 ,
where
ψ = def ( 1 H ( σ ξ h ) ) ( ξ h σ ) μ + 2 h k .

6. Analysis and Exact Solution of Problem 4

Obviously, for solving Problem 4 one can use relations (43), (45), and (46). Two cases are possible: either the no-slip condition holds on the plate y = h , or the slip regime is realized on this plate.
First, let us find the solution for the first case. Substitute y = h into (43). Since u ( h ) = 0 , we see that
C 2 = ξ h 2 2 C 1 h ,
and hence
u ( y ) = ξ 2 μ ( y 2 h 2 ) + C 1 μ ( y h ) .
Further, let us choose the value of the constant C 1 such that the Navier slip condition
μ u ( h ) = k 2 u ( h )
is satisfied on the lower wall of the channel. Taking into account (45) and (48), we rewrite (49) as follows
ξ h C 1 = 2 h k 2 C 1 μ ,
where
C 1 = ξ h μ μ + 2 k 2 h .
Let us now determine relations on the model parameters under which the above case is realized. In view of (28), the following inequality
( T n ) tan R 3 σ
holds for y = h . Using (46) and (50), we conclude that (51) is true if
ξ h σ ˜ = def σ 1 μ 2 ( μ + k 2 h ) .
Let us now consider the case when the slip regime is realized on the upper wall of the channel. Then, boundary conditions (27)–(30) reduce to the following system:
μ u ( h ) = ( σ + k 1 u ( h ) ) ,
μ u ( h ) = k 2 u ( h ) .
Using (43) and (45), one can rewrite (53) and (54) as follows:
ξ h + C 1 = σ + k 1 ξ h 2 2 μ + C 1 h μ + C 2 μ , ( ξ h + C 1 ) = k 2 ξ h 2 2 μ C 1 h μ + C 2 μ .
Solving this system, we find the values of the constants C 1 and C 2 :
C 1 = | ( ξ h ( k 1 k 2 ) + σ k 2 ) μ 2 h k 1 k 2 + μ ( k 1 + k 2 ) , C 2 = | 2 h 3 ξ k 1 k 2 + 3 h 2 μ ξ ( k 1 + k 2 ) + 4 h μ 2 ξ 2 h μ σ k 2 2 μ 2 σ 4 h k 1 k 2 + 2 μ ( k 1 + k 2 ) .
Taking into account (29), it is easy to check that the case under consideration is realized if the following inequality holds:
ξ h > σ ˜ .
When the function u is obtained, one can derive the temperature θ by solving (15) with the Dirichlet boundary condition (31).
Finally, summarizing the results of this section, we write the solution of Problem 4 in the explicit form: if the model parameters satisfy relation (52), then
u ( y ) = ξ 2 μ ( y 2 h 2 ) + λ 1 ( y h ) , p ( x , y ) = ξ x + α μ 2 ( ξ y λ 1 μ ) 2 + ρ g ( h y ) , θ ( y ) = ( ξ y λ 1 μ ) 4 12 c p μ κ ξ 2 λ 1 λ 1 2 μ 2 + h 2 ξ 2 3 c p κ ξ y + λ 1 4 μ 4 + 6 λ 1 2 h 2 μ 2 ξ 2 + h 4 ξ 4 12 c p μ κ ξ 2 ,
where
λ 1 = def ξ h μ + 2 k 2 h ;
otherwise, when inequality (55) holds, the exact solution is determined by the following formulas:
u ( y ) = ξ 2 μ ( y 2 h 2 ) + λ 2 y + ξ h 2 ( k 1 + k 2 ) + 2 ξ h μ σ h k 2 μ σ 2 h k 1 k 2 + μ ( k 1 + k 2 ) , p ( x , y ) = ξ x + α μ 2 ( ξ y λ 2 μ ) 2 + ρ g ( h y ) , θ ( y ) = ( ξ y λ 2 μ ) 4 12 c p μ κ ξ 2 λ 2 λ 2 2 μ 2 + h 2 ξ 2 3 c p κ ξ y + λ 2 4 μ 4 + 6 λ 2 2 h 2 μ 2 ξ 2 + h 4 ξ 4 12 c p μ κ ξ 2 ,
where
λ 2 = def ξ h ( k 1 k 2 ) + σ k 2 2 h k 1 k 2 + μ ( k 1 + k 2 ) .

7. Conclusions

In this work, we have studied the non-isothermal steady-state flow of a second grade fluid in the channel h y h with impermeable solid walls, taking into account the fluid energy dissipation (mechanical-to-thermal energy conversion) in the heat transfer equation. It is assumed that the flow is created by a constant pressure gradient p / x = ξ . We have established exact solutions of the nonlinear governing equations for the velocity vector, the pressure, and the temperature under the no-slip boundary conditions and threshold-type slip boundary conditions, which include Navier’s slip condition as a limit case. Moreover, we analytically solved two problems for channel flows with mixed boundary conditions, assuming that the upper and lower walls of the channel differ in their physical properties. The obtained solutions show that the pressure in the channel significantly depends on the normal stress coefficient α , especially in those layers where the change in the flow velocity in the transverse direction to the flow is large. At the same time, the velocity field is independent of α , and therefore coincides with the velocity field that occurs in the case of a Newtonian fluid ( α = 0 ). In the analysis of flows with threshold slip, the key point is the value of ξ h . If ξ h exceeds a given threshold value σ , then the slip regime holds at solid surfaces, otherwise the fluid adheres to the walls of the channel. If it is assumed that on one part of the boundary Navier’s condition is provided, while on the other one the threshold slip condition holds, then, for the slip regime, the associated threshold value σ ˜ is reduced to a certain extent, but not more than twice. An interesting feature of the obtained results is that the temperature distribution is given by a fourth-degree polynomial, and not by a quadratic function. This is due to the fact that when deriving the heat transfer equation, the simplifying assumption that the viscous energy dissipation function is identically equal to zero is not used. The proposed approach leads to a more delicate description of the heat and mass transfer in second grade fluids as well as a deep understanding of the related physical processes. Finally, note that the exact solutions obtained in the present paper can be applied to testing numerical, asymptotic, and approximate analytical methods of solving boundary value problems that describe non-isothermal flows of nanofluids.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The author thanks four anonymous reviewers for their valuable comments, which led to the improvement of the paper.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Rivlin, R.S.; Ericksen, J.L. Stress-deformation relations for isotropic materials. J. Ration. Mech. Anal. 1955, 4, 323–425. [Google Scholar] [CrossRef]
  2. Cioranescu, D.; Girault, V.; Rajagopal, K.R. Mechanics and Mathematics of Fluids of the Differential Type; Springer: Cham, Switzerland, 2016. [Google Scholar] [CrossRef]
  3. Abbas, S.Z.; Waqas, M.; Thaljaoui, A.; Zubair, M.; Riahi, A.; Chu, Y.M.; Khan, W.A. Modeling and analysis of unsteady second-grade nanofluid flow subject to mixed convection and thermal radiation. Soft Comput. 2022, 26, 1033–1042. [Google Scholar] [CrossRef]
  4. Imran, M.; Yasmin, S.; Waqas, H.; Khan, S.A.; Muhammad, T.; Alshammari, N.; Hamadneh, N.N.; Khan, I. Computational analysis of nanoparticle shapes on hybrid nanofluid flow due to flat horizontal plate via solar collector. Nanomaterials 2022, 12, 663. [Google Scholar] [CrossRef] [PubMed]
  5. Shah, F.; Hayat, T.; Momani, S. Non-similar analysis of the Cattaneo–Christov model in MHD second-grade nanofluid flow with Soret and Dufour effects. Alex. Eng. J. 2023, 70, 25–35. [Google Scholar] [CrossRef]
  6. Hosseinzadeh, K.; Mardani, M.R.; Paikar, M.; Hasibi, A.; Tavangar, T.; Nimafar, M.; Ganji, D.D.; Shafii, M.B. Investigation of second grade viscoelastic non-Newtonian nanofluid flow on the curve stretching surface in presence of MHD. Results Eng. 2023, 17, 100838. [Google Scholar] [CrossRef]
  7. Toms, B.A. Some observations on the flow of linear polymer solutions through straight tubes at large Reynolds numbers. Proc. First Int. Congr. Rheol. 1948, 2, 135–141. [Google Scholar]
  8. Barenblatt, G.I.; Kalashnikov, V.N. Effect of high-molecular formations on turbulence in dilute polymer solutions. Fluid Dyn. 1968, 3, 45–48. [Google Scholar] [CrossRef]
  9. Pisolkar, V.G. Effect of drag reducing additives on pressure loss across transitions. Nature 1970, 225, 936–937. [Google Scholar] [CrossRef]
  10. Amfilokhiev, V.B.; Pavlovskii, V.A.; Mazaeva, N.P.; Khodorkovskii, Y.S. Flows of polymer solutions in the presence of convective accelerations. Trudy Leningrad. Korablestr. Inst. 1975, 96, 3–9. [Google Scholar]
  11. Amfilokhiev, V.B.; Pavlovskii, V.A. Experimental data on laminar-turbulent transition for flows of polymer solutions in pipes. Trudy Leningrad. Korablestr. Inst. 1976, 104, 3–5. [Google Scholar]
  12. Fu, Z.; Otsuki, T.; Motozawa, M.; Kurosawa, T.; Yu, B.; Kawaguchi, Y. Experimental investigation of polymer diffusion in the drag-reduced turbulent channel flow of inhomogeneous solution. Int. J. Heat Mass Transf. 2014, 77, 860–873. [Google Scholar] [CrossRef]
  13. Han, W.J.; Dong, Y.Z.; Choi, H.J. Applications of water-soluble polymers in turbulent drag reduction. Processes 2017, 5, 24. [Google Scholar] [CrossRef]
  14. Pukhnachev, V.V.; Frolovskaya, O.A. On the Voitkunskii–Amfilokhiev–Pavlovskii model of motion of aqueous polymer solutions. Proc. Steklov Inst. Math. 2018, 300, 168–181. [Google Scholar] [CrossRef]
  15. Frolovskaya, O.A.; Pukhnachev, V.V. The problem of filling a spherical cavity in an aqueous solution of polymers. Polymers 2022, 14, 4259. [Google Scholar] [CrossRef] [PubMed]
  16. Voitkunskii, Y.I.; Amfilokhiev, V.B.; Pavlovskii, V.A. Equations of motion of a fluid, with its relaxation properties taken into account. Trudy Leningrad. Korablestr. Inst. 1970, 69, 19–26. [Google Scholar]
  17. Lokshin, A.A. Tauberian Theory of Wave Fronts in Linear Hereditary Elasticity; Springer: Singapore, 2020; ISBN 978-981-15-8577-7. [Google Scholar]
  18. Pavlovskii, V.A. On the theoretical description of weak water solutions of polymers. Dokl. Akad. Nauk. 1971, 200, 809–812. [Google Scholar]
  19. Baranovskii, E.S. Flows of a polymer fluid in domain with impermeable boundaries. Comput. Math. Math. Phys. 2014, 54, 1589–1596. [Google Scholar] [CrossRef]
  20. Dunn, J.E.; Fosdick, R.L. Thermodynamics, stability and boundedness of fluids of complexity 2 and fluids of second grade. Arch. Ration. Mech. Anal. 1974, 56, 191–252. [Google Scholar] [CrossRef]
  21. Fosdick, R.L.; Rajagopal, K.R. Anomalous features in the model of “second order fluids”. Arch. Ration. Mech. Anal. 1979, 70, 145–152. [Google Scholar] [CrossRef]
  22. Dunn, J.E.; Rajagopal, K.R. Fluids of differential type: Critical review and thermodynamic analysis. Int. J. Eng. Sci. 1995, 33, 689–729. [Google Scholar] [CrossRef]
  23. Rajagopal, K.R. On some unresolved issues in non-linear fluid dynamics. Russ. Math. Surv. 2003, 58, 319–330. [Google Scholar] [CrossRef]
  24. Xu, H.; Cui, J. Mixed convection flow in a channel with slip in a porous medium saturated with a nanofluid containing both nanoparticles and microorganisms. Int. J. Heat Mass Transf. 2018, 125, 1043–1053. [Google Scholar] [CrossRef]
  25. Wang, G.J.; Hadjiconstantinou, N.G. Universal molecular-kinetic scaling relation for slip of a simple fluid at a solid boundary. Phys. Rev. Fluids 2019, 4, 064201. [Google Scholar] [CrossRef]
  26. Wilms, P.; Wieringa, J.; Blijdenstein, T.; van Malssen, K.; Hinrichs, J. Wall slip of highly concentrated non-Brownian suspensions in pressure driven flows: A geometrical dependency put into a non-Newtonian perspective. J. Non-Newton. Fluid Mech. 2020, 282, 104336. [Google Scholar] [CrossRef]
  27. Ghahramani, N.; Georgiou, G.C.; Mitsoulis, E.; Hatzikiriakos, S.G. J.G. Oldroyd’s early ideas leading to the modern understanding of wall slip. J. Non-Newton. Fluid Mech. 2021, 293, 104566. [Google Scholar] [CrossRef]
  28. Rana, M.A.; Latif, A. Three-dimensional free convective flow of a second-grade fluid through a porous medium with periodic permeability and heat transfer. Bound. Value Probl. 2019, 44. [Google Scholar] [CrossRef]
  29. Sehra; Haq, S.U.; Jan, S.U.; Bilal, R.; Alzahrani, J.H.; Khan, I.; Alzahrani, A. Heat-mass transfer of MHD second grade fluid flow with exponential heating, chemical reaction and porosity by using fractional Caputo-Fabrizio derivatives. Case Stud. Therm. Eng. 2022, 36, 102104. [Google Scholar] [CrossRef]
  30. Baranov, A.V. Nonisothermal dissipative flow of viscous liquid in a porous channel. High Temp. 2017, 55, 414–419. [Google Scholar] [CrossRef]
  31. Goruleva, L.S.; Prosviryakov, E.Y. A new class of exact solutions to the Navier–Stokes equations with allowance for internal heat release. Opt. Spectrosc. 2022, 130, 365–370. [Google Scholar] [CrossRef]
  32. Privalova, V.V.; Prosviryakov, E.Y. A new class of exact solutions of the Oberbeck–Boussinesq equations describing an incompressible fluid. Theor. Found. Chem. Eng. 2022, 56, 331–338. [Google Scholar] [CrossRef]
  33. Baranovskii, E.S.; Artemov, M.A. Steady flows of second-grade fluids in a channel. Vestn. S.-Peterb. Univ. Prikl. Mat. Inf. Protsessy Upr. 2017, 13, 342–353. [Google Scholar] [CrossRef]
  34. Baranovskii, E.S.; Artemov, M.A. Steady flows of second grade fluids subject to stick-slip boundary conditions. In Proceedings of the 23rd International Conference Engineering Mechanics, Svratka, Czech Republic, 15–18 May 2017; pp. 110–113. [Google Scholar]
  35. Ting, T.W. Certain unsteady flows of second grade fluids. Arch. Ration. Mech. Anal. 1963, 14, 1–26. [Google Scholar] [CrossRef]
  36. Coleman, B.D.; Duffin, R.J.; Mizel, V.J. Instability, uniqueness, and nonexistence theorems for the equation ut = uxxuxtx on a strip. Arch. Ration. Mech. Anal. 1965, 19, 100–116. [Google Scholar] [CrossRef]
  37. Hron, J.; Le Roux, C.; Malek, J.; Rajagopal, K.R. Flows of incompressible fluids subject to Navier’s slip on the boundary. Comput. Math. Appl. 2008, 56, 2128–2143. [Google Scholar] [CrossRef]
  38. Nazar, M.; Fetecau, C.; Vieru, D.; Fetecau, C. New exact solutions corresponding to the second problem of Stokes for second grade fluids. Nonlinear Anal. Real World Appl. 2010, 11, 584–591. [Google Scholar] [CrossRef]
  39. Fetecau, C.; Zafar, A.A.; Vieru, D.; Awrejcewicz, J. Hydromagnetic flow over a moving plate of second grade fluids with time fractional derivatives having non-singular kernel. Chaos Solit. Fractals. 2020, 130, 109454. [Google Scholar] [CrossRef]
  40. Fetecau, C.; Vieru, D. On an important remark concerning some MHD motions of second grade fluids through porous media and its applications. Symmetry 2022, 14, 1921. [Google Scholar] [CrossRef]
  41. Fetecau, C.; Vieru, D. General solutions for some MHD motions of second-grade fluids between parallel plates embedded in a porous medium. Symmetry 2023, 15, 183. [Google Scholar] [CrossRef]
  42. Cioranescu, D.; Ouazar, E.H. Existence and uniqueness for fluids of second-grade. Nonlinear Partial. Differ. Equ. Their Appl. 1984, 109, 178–197. [Google Scholar]
  43. Cioranescu, D.; Girault, V. Weak and classical solutions of a family of second grade fluids. Int. J. Non-Linear Mech. 1997, 32, 317–335. [Google Scholar] [CrossRef]
  44. Le Roux, C. Existence and uniqueness of the flow of second-grade fluids with slip boundary conditions. Arch. Ration. Mech. Anal. 1999, 148, 309–356. [Google Scholar] [CrossRef]
  45. Baranovskii, E.S. Existence results for regularized equations of second-grade fluids with wall slip. Electron. J. Qual. Theory Differ. Equ. 2015, 91. [Google Scholar] [CrossRef]
  46. Baranovskii, E.S. Weak solvability of equations modeling steady-state flows of second-grade fluids. Differ. Equ. 2020, 56, 1318–1323. [Google Scholar] [CrossRef]
  47. Baranovskii, E.S. Solvability of the stationary optimal control problem for motion equations of second grade fluids. Sib. Electron. Math. Rep. 2012, 9, 554–560. [Google Scholar]
  48. Chemetov, N.V.; Cipriano, F. Optimal control for two-dimensional stochastic second grade fluids. Stoch. Proc. Their Appl. 2018, 128, 2710–2749. [Google Scholar] [CrossRef]
  49. Baranovskii, E.S. Optimal boundary control of the Boussinesq approximation for polymeric fluids. J. Optim. Theory Appl. 2021, 189, 623–645. [Google Scholar] [CrossRef]
  50. Almeida, A.; Chemetov, N.V.; Cipriano, F. Uniqueness for optimal control problems of two-dimensional second grade fluids. Electron. J. Differ. Equ. 2022, 22. [Google Scholar]
  51. Le Roux, C.; Tani, A. Steady solutions of the Navier–Stokes equations with threshold slip boundary conditions. Math. Meth. Appl. Sci. 2007, 30, 595–624. [Google Scholar] [CrossRef]
  52. Navier, C.L.M.H. Mémoire sur le lois du mouvement des fluides. Mémoires L’académie Sci. L’institut Fr. 1823, 6, 389–416. [Google Scholar]
  53. Solonnikov, V.A.; Ščadilov, V.E. On a boundary value problem for a stationary system of Navier–Stokes equations. Proc. Steklov Inst. Math. 1973, 125, 186–199. [Google Scholar]
  54. Burmasheva, N.V.; Larina, E.A.; Prosviryakov, E.Y. A Couette-type flow with a perfect slip condition on a solid surface. Vestn. Tomsk. Gos. Univ. Mat. Mekh. 2021, 74, 79–94. [Google Scholar] [CrossRef]
  55. Gkormpatsis, S.D.; Housiadas, K.D.; Beris, A.N. Steady sphere translation in weakly viscoelastic UCM/Oldroyd-B fluids with perfect slip on the sphere. Eur. J. Mech. B Fluids 2022, 95, 335–346. [Google Scholar] [CrossRef]
Figure 1. Flow configuration.
Figure 1. Flow configuration.
Nanomaterials 13 01409 g001
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Baranovskii, E.S. Exact Solutions for Non-Isothermal Flows of Second Grade Fluid between Parallel Plates. Nanomaterials 2023, 13, 1409. https://doi.org/10.3390/nano13081409

AMA Style

Baranovskii ES. Exact Solutions for Non-Isothermal Flows of Second Grade Fluid between Parallel Plates. Nanomaterials. 2023; 13(8):1409. https://doi.org/10.3390/nano13081409

Chicago/Turabian Style

Baranovskii, Evgenii S. 2023. "Exact Solutions for Non-Isothermal Flows of Second Grade Fluid between Parallel Plates" Nanomaterials 13, no. 8: 1409. https://doi.org/10.3390/nano13081409

APA Style

Baranovskii, E. S. (2023). Exact Solutions for Non-Isothermal Flows of Second Grade Fluid between Parallel Plates. Nanomaterials, 13(8), 1409. https://doi.org/10.3390/nano13081409

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop