Next Article in Journal
Model Membrane Systems Used to Study Plasma Membrane Lipid Asymmetry
Next Article in Special Issue
Control Problem Related to 2D Stokes Equations with Variable Density and Viscosity
Previous Article in Journal
Group Theory: Mathematical Expression of Symmetry in Physics
Previous Article in Special Issue
A 3D Non-Stationary Micropolar Fluids Equations with Navier Slip Boundary Conditions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Exact Solutions to the Navier–Stokes Equations with Couple Stresses

by
Evgenii S. Baranovskii
1,*,
Natalya V. Burmasheva
2,3 and
Evgenii Yu. Prosviryakov
2,4
1
Department of Applied Mathematics, Informatics and Mechanics, Voronezh State University, 394018 Voronezh, Russia
2
Sector of Nonlinear Vortex Hydrodynamics, Institute of Engineering Science, Ural Branch of the Russian Academy of Sciences, 620049 Ekaterinburg, Russia
3
Ural Institute of Humanities, Ural Federal University, 620002 Ekaterinburg, Russia
4
Institute of Fundamental Education, Ural Federal University, 620002 Ekaterinburg, Russia
*
Author to whom correspondence should be addressed.
Symmetry 2021, 13(8), 1355; https://doi.org/10.3390/sym13081355
Submission received: 21 June 2021 / Revised: 16 July 2021 / Accepted: 19 July 2021 / Published: 26 July 2021
(This article belongs to the Special Issue Mathematical Fluid Dynamics and Symmetry)

Abstract

:
This article discusses the possibility of using the Lin–Sidorov–Aristov class of exact solutions and its modifications to describe the flows of a fluid with microstructure (with couple stresses). The presence of couple shear stresses is a consequence of taking into account the rotational degrees of freedom for an elementary volume of a micropolar liquid. Thus, the Cauchy stress tensor is not symmetric. The article presents exact solutions for describing unidirectional (layered), shear and three-dimensional flows of a micropolar viscous incompressible fluid. New statements of boundary value problems are formulated to describe generalized classical Couette, Stokes and Poiseuille flows. These flows are created by non-uniform shear stresses and velocities. A study of isobaric shear flows of a micropolar viscous incompressible fluid is presented. Isobaric shear flows are described by an overdetermined system of nonlinear partial differential equations (system of Navier–Stokes equations and incompressibility equation). A condition for the solvability of the overdetermined system of equations is provided. A class of nontrivial solutions of an overdetermined system of partial differential equations for describing isobaric fluid flows is constructed. The exact solutions announced in this article are described by polynomials with respect to two coordinates. The coefficients of the polynomials depend on the third coordinate and time.

1. Introduction

The overwhelming majority of studies dealing with fluid flows are based on the application of the conventional Navier–Stokes equations supplemented by the incompressibility condition [1,2]. The Navier–Stokes equations are derived from the postulates (hypotheses) of the Newtonian mechanics of continua, each particle of which is viewed as a material point. When the representative volume of a continuum is substituted by a material point, it is considered by default to have three degrees of freedom (translational degrees of freedom). The application of this approach imposes restrictions on studying changes in fluid viscosity, friction coefficients, and other surface effects [3,4,5,6,7,8]. The difference between the experimentally and theoretically obtained results reported in the pioneering study [8] stems from ignoring the rotational (orientational) degrees of freedom of the representative volume of a continuum. By taking into account additional degrees of freedom of the elementary volume of a deformable medium (continuum), one finds that the Cauchy stresses do not balance each other. The stress tensor becomes asymmetric in this case since additional stresses occur due to taking into account the deformation properties of the vortex velocities of elementary fluid volumes [9,10,11]. These media are currently termed micropolar [9,10,11,12,13,14]. As applied to elastic bodies, media with additional tangential stresses were first described in [15]. It can be stated that micropolar fluids began to be studied as late as in the mid-1960s [3,4].
In [4], not only were the Navier–Stokes equations derived to describe fluids with a representative volume having six degrees of freedom but also the first exact solutions were constructed and studied. Since [4] was published, steady-state and non-steady-state flows of micropolar viscous incompressible fluids have been studied in an exact formulation for unidirectional flows, e.g., in [4,8,16,17]. The exact solutions for Newtonian fluids were extended to micropolar media for Couette flows, the first and second Stokes problems, the Poiseuille flow, and their combinations and modifications.
There still appears to be no publications discussing classes of exact solutions to the Navier–Stokes equations for shearing and three-dimensional micropolar viscous incompressible fluids. In this paper, the Lin–Sidorov–Aristov family found in [18,19,20] is chosen as a basis for constructing exact solutions. This class prescribes functional variable separation for the velocity field described by linear forms with respect to two (horizontal or longitudinal) coordinates x and y. The coefficients of these forms depend on the third coordinate (vertical or transverse) z and time t. In the Lin–Sidorov–Aristov class, pressure is a quadratic form similar to the velocity field. Various extensions and modifications of the Lin–Sidorov–Aristov family can be found in [21,22,23,24,25]. It was shown in [21,22] that the exact solution for the velocity field can nonlinearly depend on two coordinates.
In this paper, in view of the relevance of the research and on account of the insufficient completeness of the exact integration of micropolar fluid motion equations, we construct new exact solutions to the generalized Navier–Stokes equations for unidirectional, shearing, and three-dimensional flows.

2. The Navier–Stokes Equations with Couple Stresses

The incompressible Navier–Stokes equations with tangential couple stresses are written as follows [3,4]:
d V d t = P + ν Δ V μ Δ 2 V , · V = 0 ,
where V ( t , x , y , z ) = ( V x , V y , V z ) is the velocity vector; P = p / ρ 0 is the pressure normalized to the fluid density ρ 0 ; ν is kinematic viscosity ( ν > 0 ) ; μ is the couple stress viscosity parameter ( μ > 0 ) ; the symbols ∇, Δ , and Δ 2 denote, respectively, the Hamilton operator, the Laplace operator, and the biharmonic operator:
= x , y , z , Δ = 2 x 2 + 2 y 2 + 2 z 2 , Δ 2 = 4 x 4 + 4 y 4 + 4 z 4 + 2 4 x 2 y 2 + 2 4 x 2 z 2 + 2 4 y 2 z 2 ;
where d d t = t + ( V · ) is the Lagrangian derivative; t is the local derivative, ( V · ) is the convective derivative [1]. The contribution of the couple stresses μ Δ 2 V is due to the possibility of fluid particles to have rotational interaction (the representative fluid volume has rotational degrees of freedom) [3,4,13,16,17,26,27,28,29]. In the limit case μ = 0 , system (1) passes to the standard Navier–Stokes equations, describing flows of a Newtonian fluid.
Note that model (1) can be considered a suitable μ -approximation of the Navier–Stokes equations [30]. Ladyzhenskaya [31] established the global unique solvability of an initial boundary value problem for (1) with the following boundary conditions:
V | Γ = 0 , Δ V × n | Γ = 0 ,
where Γ is the boundary of the flow domain and n denotes the exterior unit normal to the surface Γ .
We write system (1) in coordinate form as follows:
V x t + V x V x x + V y V x y + V z V x z = P x + ν 2 V x x 2 + 2 V x y 2 + 2 V x z 2 μ 4 V x x 4 + 4 V x y 4 + 4 V x z 4 + 2 4 V x x 2 y 2 + 2 4 V x x 2 z 2 + 2 4 V x y 2 z 2 ,
V y t + V x V y x + V y V y y + V z V y z = P y + ν 2 V y x 2 + 2 V y y 2 + 2 V y z 2 μ 4 V y x 4 + 4 V y y 4 + 4 V y z 4 + 2 4 V y x 2 y 2 + 2 4 V y x 2 z 2 + 2 4 V y y 2 z 2 ,
V z t + V x V z x + V y V z y + V z V z z = P z + ν 2 V z x 2 + 2 V z y 2 + 2 V z z 2 μ 4 V z x 4 + 4 V z y 4 + 4 V z z 4 + 2 4 V z x 2 y 2 + 2 4 V z x 2 z 2 + 2 4 V z y 2 z 2 ,
V x x + V y y + V z z = 0 .
System (2)–(5) consists of four scalar equations for the determination of pressure P and three velocity vector projections V x , V y , and V z . The integration of this system is complicated due to quadratic nonlinearity and the presence of the fourth derivatives of the velocity vector projections V x , V y , V z involved in the generalized Navier–Stokes equations (2)–(4).

3. Unidirectional Flows

In order to begin constructing exact solutions of system (2)–(5), we model unidirectional flows, i.e., we consider fluid flows with the velocity defined as follows:
V = ( V x ( x , y , z , t ) , 0 , 0 ) .
The velocity field (6) satisfies the incompressibility condition (5) if the following condition
V x x = 0
holds. Therefore V x depends only on two spatial coordinates and time, i.e.,
V x = V x ( y , z , t ) .
By substituting Equation (6) into (3) and (4), we arrive at the conclusion that the pressure field is determined by only one spatial coordinate and time:
P = P ( x , t ) .
In view of (7), Equation (2) also becomes simplified:
V x t = P x + ν 2 V x y 2 + 2 V x z 2 μ 4 V x y 4 + 4 V x z 4 + 2 4 V x y 2 z 2 .
If the function P is independent of the coordinate x but determined only by the time t (i.e., P = P ( t ) ), then the first term in the right-hand side of (8) vanishes:
V x t = ν 2 V x y 2 + 2 V x z 2 μ 4 V x y 4 + 4 V x z 4 + 2 4 V x y 2 z 2 .
Herewith, formally, the flow under study is not isobaric since the pressure is time dependent: P = P ( t ) .
Note that, due to the structure of Equation (9), steady-state flows V = ( V x ( y , z ) , 0 , 0 ) (as distinct from the classical Couette flow) depend on viscosity and velocity V x , in general terms, is not a harmonic function,
2 V x y 2 + 2 V x z 2 = μ ν 4 V x y 4 + 4 V x z 4 + 2 4 V x y 2 z 2 0 .
In what follows, we shall present new exact solutions of Equation (9), which also satisfy the generalized Navier–Stokes equations (1).
Let us now construct an extension of the classical Couette solution found in [32] to the case of couple stresses taken into account, i.e., consider a solution of the following form
V x = U ( z , t ) .
Note that the models for describing Couette flow (and some other flows) with couple stresses taken into account were discussed, e.g., in [4].
The substitution of Equation (10) into (9) produces the following heat-conduction-type equation
U t = ν 2 U z 2 μ 4 U z 4 .
If herewith the flow represented by Equation (10) is steady-state ( V x = U ( z ) ), the solution for velocity U is a combination of linear and exponential functions:
U = C 1 z + C 2 + C 3 e z ν / μ + C 4 e z ν / μ .
Moreover, the exponential terms in Equation (12) appear only when couple stresses are taken into account. If we solve Equation (11) for μ = 0 , the solution acquires the form of the classical Couette solution, where the velocity profile is defined by the linear function
U = C 1 z + C 2 .
Here, C i denote integration constants.
Consider another particular exact solution of Equation (9)
V x = y u 1 ( z , t ) .
By substituting (13) into Equation (9), we obtain
y u 1 t = y ν 2 u 1 z 2 μ 4 u 1 z 4 .
That is, the function u 1 satisfies Equation (11) and, hence, its solution can be represented as Equation (12). It can be easily shown that in view of the linearity of Equation (9), the linear combination of the solutions represented by Equations (10) and (13)
V x = U ( z , t ) + y u 1 ( z , t )
also satisfies Equation (9).
We now add a nonlinear (in the y coordinate) term to (14), i.e., consider a solution of the following form
V x = U ( z , t ) + y u 1 ( z , t ) + y 2 2 u 2 ( z , t ) .
Substitute the last sum into (9). Some simple transformations yield
U t + y u 1 t + y 2 2 u 2 t = ν u 2 + 2 U z 2 + y 2 u 1 z 2 + y 2 2 2 u 2 z 2 μ 4 U z 4 + y 4 u 1 z 4 + y 2 2 4 u 2 z 4 .
Applying the method of undetermined coefficients to this equation, we obtain the following system of equations:
U t = ν u 2 + 2 U z 2 μ 4 U z 4 , u 1 t = ν 2 u 1 z 2 μ 4 u 1 z 4 , u 2 t = ν 2 u 2 z 2 μ 4 u 2 z 4 .
The latter two equations in this system are equations of the form of (11) and the equation for determining the homogeneous term U in (15) is not isolated any longer. Thus, the consideration of nonlinearity in the variable results in the fact that the solution represented by (15) stops being a superposition of the previously presented solutions. As the power at the y coordinate increases, this tendency holds true. In other words, a unidirectional flow is described by the following exact solution
V x = U ( z , t ) + k = 1 n y k k ! u k ( z , t ) .
This exact solution, as all those presented in what follows, need to be studied for stability. However, studying for stability questions is not the aim of this paper.
Note that unidirectional flows of non-Newtonian fluids have been studied in many works (see, e.g., [33,34,35,36] and the references therein).

4. Exact Solutions for Three-Dimensional Flows

In order to study the properties of three-dimensional flows of viscous fluids, one needs to have a reserve of exact solutions to the Navier–Stokes equations. The exact solutions inheriting nonlinear properties of the motion equations help to understand better the equations structure in further analytical and numerical integrations [23,24,25,37,38,39,40,41,42,43,44].
As a preliminary, the exact solution of the system of nonlinear partial Equations (2)–(5) is sought within the Lin–Sidorov–Aristov family [18,19,20]:
V x ( x , y , z , t ) = U ( z , t ) + x u 1 ( z , t ) + y u 2 ( z , t ) , V y ( x , y , z , t ) = V ( z , t ) + x v 1 ( z , t ) + y v 2 ( z , t ) , V z ( z , t ) w , w = w ( z , t ) .
Note that the expressions in Equation (16), they can be treated as the Taylor series expansion of the velocity vector components that is confined to linear terms [37].
We substitute the class (16) into (2) and immediately ignore the terms containing the second and fourth derivatives with respect to the x and y variables since the expressions in Equation (16) linearly depend on these coordinates; this gives
( U + x u 1 + y u 2 ) t + ( U + x u 1 + y u 2 ) ( U + x u 1 + y u 2 ) x + + ( V + x v 1 + y v 2 ) ( U + x u 1 + y u 2 ) y + w ( U + x u 1 + y u 2 ) z = = P x + ν 2 ( U + x u 1 + y u 2 ) z 2 μ 4 ( U + x u 1 + y u 2 ) z 4 .
We now make the following simple transformations in this equation by computing the partial derivatives and combining similar terms at the identical powers of the independent variables x and y:
U t + U u 1 + V u 2 + w U z + x u 1 t + u 1 2 + v 1 u 2 + w u 1 z + + y u 2 t + u 1 u 2 + v 2 u 2 + w u 2 z = P x + ν 2 U z 2 μ 4 U z 4 + + x ν 2 u 1 z 2 μ 4 u 1 z 4 + y ν 2 u 2 z 2 μ 4 u 2 z 4 .
The structure of the obtained expression suggests that the pressure should be viewed as a quadratic form of the x and y variables. In other words, we have
P ( x , y , z , t ) = P 0 ( z , t ) + x P 1 ( z , t ) + y P 2 ( z , t ) + + x 2 2 P 11 ( z , t ) + x y P 12 ( z , t ) + y 2 2 P 22 ( z , t ) .
Taking into account (18), when the method of undetermined coefficients is applied, Equation (17) splits into the following three equations:
U t + U u 1 + V u 2 + w U z = P 1 + ν 2 U z 2 μ 4 U z 4 , u 1 t + u 1 2 + v 1 u 2 + w u 1 z = P 11 + ν 2 u 1 z 2 μ 4 u 1 z 4 , u 2 t + u 1 u 2 + v 2 u 2 + w u 2 z = P 12 + ν 2 u 2 z 2 μ 4 u 2 z 4 .
Similar operations with the substitution of Equations (16) and (18) into the other two equations of system (2)–(4) yield, respectively, the following two systems:
V t + U v 1 + V v 2 + w V z = P 2 + ν 2 V z 2 μ 4 V z 4 , v 1 t + u 1 v 1 + v 1 v 2 + w v 1 z = P 12 + ν 2 v 1 z 2 μ 4 v 1 z 4 , v 2 t + u 2 v 1 + v 2 2 + w v 2 z = P 22 + ν 2 v 2 z 2 μ 4 v 2 z 4 ;
w t + w w z = P 0 z + ν 2 w z 2 μ 4 w z 4 , P 1 z = 0 , P 2 z = 0 , P 11 z = 0 , P 12 z = 0 , P 22 z = 0 .
Note that the last equations in system (21) indicate that the overwhelming majority of the coefficients in Equation (18) are independent of the vertical coordinate; consequently, we gain a possibility of correcting the structure of the pressure field:
P ( x , y , z , t ) = P 0 ( z , t ) + x P 1 ( t ) + y P 2 ( t ) + x 2 2 P 11 ( t ) + x y P 12 ( t ) + y 2 2 P 22 ( t ) .
Herewith, the homogeneous component of this field (background pressure P 0 ) is determined by integrating the first equation of system (21) up to the additive function of time,
P 0 = w t d z w 2 2 + ν w z μ 3 w z 3 + c 0 ( t ) .
In addition to the equations contained in systems (19)–(21), the coefficients in Equation (16) must satisfy the following equation
u 1 + v 2 + w z = 0 .
Equation (23) is obtained from the substitution of class (16) into the incompressibility condition (5).
By analogy with the derivation of systems (19)–(21), other solutions with the velocity field arbitrarily dependent on the horizontal coordinates can be constructed [21,22]. For example, the following exact solution of system (2)–(4) is valid:
V x = k = 0 n U k , V y = k = 0 n V k , V z = k = 0 n 1 W k , A = k = 0 n 2 n A k , B = k = 0 n 2 n B k , C = k = 0 n 2 n + 1 C k , P = k = 0 n 2 n + 1 P k .
Here, the forms of U k , V k , W k , P k are determined by the following expressions:
U k = 1 k ! i = 0 k C k i U i ( k i ) ( z , t ) x i y k i , V k = 1 k ! i = 0 k C k i V i ( k i ) ( z , t ) x i y k i , W k = 1 k ! i = 0 k C k i W i ( k i ) ( z , t ) x i y k i , A k = 1 k ! i = 0 k C k i A i ( k i ) ( z , t ) x i y k i , B k = 1 k ! i = 0 k C k i B i ( k i ) ( z , t ) x i y k i , C k = 1 k ! i = 0 k C k i C i ( k i ) ( z , t ) x i y k i , P k = 1 k ! i = 0 k C k i P i ( k i ) ( z , t ) x i y k i ,
where C k i = k ! ( k i ) ! i ! is a binomial coefficient.

5. Exact Solutions for Shearing Flows

Let us now consider an important particular case of isothermal flows of viscous fluids defined by Equations (16)–(23), namely, shearing flows. To perform this, we assume that the vertical velocity V z is zero, i.e.,
V z ( z , t ) = w ( z , t ) = 0 .
Using this assumption, we derive from (19) and (20) the following two systems:
u 1 t + u 1 2 + v 1 u 2 = P 11 + ν 2 u 1 z 2 μ 4 u 1 z 4 , u 2 t + u 1 u 2 + v 2 u 2 = P 12 + ν 2 u 2 z 2 μ 4 u 2 z 4 , v 1 t + u 1 v 1 + v 1 v 2 = P 12 + ν 2 v 1 z 2 μ 4 v 1 z 4 , v 2 t + u 2 v 1 + v 2 2 = P 22 + ν 2 v 2 z 2 μ 4 v 2 z 4
and
U t + U u 1 + V u 2 = P 1 + ν 2 U z 2 μ 4 U z 4 , V t + U v 1 + V v 2 = P 2 + ν 2 V z 2 μ 4 V z 4 .
All the coefficients in Equation (22) for the pressure depend only on the time coordinate t. The scientific literature on hydrodynamics reports that shearing flows are used in oceanology, atmosphere physics, astrophysics, and other sciences studying large-scale flows [23,24,25,45,46]. Large-scale flows are not only fluid motions at the planetary or galactic scale. The class of large-scale fluxes encompasses all flows with geometric anisotropy (thin layer approximation). In this case, if one filters out surface waves and neglects free surface deformation, the hydrodynamics of large-scale flows (fluid motions in thin layers) are describable by shearing flows V ( x , y , z , t ) = ( V x , V y , 0 ) .
Note that, when shearing flows are considered ( V z ( z , t ) = 0 ), the problem of integrating the generalized Navier–Stokes Equations (2)–(5) becomes noticeably more complicated; namely, we obtain an overdetermined system since it follows from Equation (4) that the pressure function must be specified by boundary conditions. This is what motivates us to start from constructing exact solutions for three-dimensional and unidirectional fluid flows. In other words, it has been demonstrated that the intermediate reduction of some components of the velocity field to zero does not always simplify the problems; it may even complicate them significantly. In what follows, we demonstrate that there exists a nontrivial exact solution in the Lin–Sidorov–Aristov class for shearing flows of the form ( V x ( x , y , z , t ) , V y ( x , y , z , t ) , 0 ) by the mathematics developed in [23,24,25,47,48,49].
In the case of shearing flows, Equation (23) allows the number of unknowns in system (25) to be reduced by using the relation between the spatial accelerations:
v 2 = u 1 .
Equation (23) follows from the incompressibility condition (5).
System (25) is transformed in view of (27) as follows:
u 1 t ν 2 u 1 z 2 + μ 4 u 1 z 4 + u 1 2 + v 1 u 2 = P 11 , u 1 t ν 2 u 1 z 2 + μ 4 u 1 z 4 ( u 1 2 + v 1 u 2 ) = P 22 , u 2 t ν 2 u 2 z 2 + μ 4 u 2 z 4 = P 12 , v 1 t ν 2 v 1 z 2 + μ 4 v 1 z 4 = P 12 .
Comparing the first two equations of this system with one another, we can conclude that they will be simultaneous only if the following condition
u 1 2 + v 1 u 2 = P 11 + P 22 2 .
holds.
Let us define a linear differential operator L ν , μ by the following formula
L ν , μ = t ν 2 z 2 + μ 4 z 4 .
In view of the consistency condition (29), system (28) can be represented as follows:
L ν , μ u 1 = P 11 P 22 2 , L ν , μ u 2 = P 12 , L ν , μ v 1 = P 12 .
Similar results for isobaric two-dimensional and quasi-two-dimensional flows were obtained earlier (see, e.g., [22,47,50]), but the possibility of the rotational interaction of fluid particles was not then taken into account. Herewith, system (30) generalizes the results reported in [47,48,49] and it can be reduced to them by passing to the limit as μ 0 in the expression for the operator L ν , μ :
L ν , 0 = t ν 2 z 2 .
Note that the consistency condition (29) for the overdetermined system of Equations (25)–(27) can be obtained on other grounds. Let us write the Navier–Stokes equations (2)–(4) for shearing flows as follows:
V x t + V x V x x + V y V x y = P x + ν 2 V x x 2 + 2 V x y 2 + 2 V x z 2 μ 4 V x x 4 + 4 V x y 4 + 4 V x z 4 + 2 4 V x x 2 y 2 + 2 4 V x x 2 z 2 + 2 4 V x y 2 z 2 ,
V y t + V x V y x + V y V y y = P y + ν 2 V y x 2 + 2 V y y 2 + 2 V y z 2 μ 4 V y x 4 + 4 V y y 4 + 4 V y z 4 + 2 4 V y x 2 y 2 + 2 4 V y x 2 z 2 + 2 4 V y y 2 z 2 .
These equations are closed by the incompressibility condition
V x x + V y y = 0 .
We differentiate Equation (31) with respect to the variable x and Equation (32) with respect to y and add up the obtained expressions. Some simple transformations and the use of the incompressibility condition (33) result in the following relationship (the consistency condition):
V x y V y x V x x V y y = 1 2 2 P x 2 + 2 P y 2 .
The particular case of the condition represented by Equation (34) for isobaric flows was obtained in [22].
Substituting expressions (16) and (18) into Equation (34), we obtain
( U + x u 1 + y u 2 ) y · ( V + x v 1 + y v 2 ) x ( U + x u 1 + y u 2 ) x · ( V + x v 1 + y v 2 ) y = = 1 2 2 x 2 + 2 y 2 P 0 + x P 1 + y P 2 + x 2 2 P 11 + x y P 12 + y 2 2 P 22 .
Elementary transformations for calculating partial derivatives produce the following relation
u 2 v 1 u 1 v 2 = P 11 + P 22 2 .
Finally, by substituting Equation (27) into the latter equation, we arrive at the consistency condition (29).
By analogy with the results reported in [22], the class of exact solutions (16) for system (19)–(21) can be expanded. Note the following velocity field
V x ( i ) = k = 0 n U k ( i ) ( z , t ) y k k ! , V y ( i ) = V ( i ) ( z , t )
and the pressure field
P ( x , y , z , t ) = P 0 ( z , t ) + x P 1 ( t ) + y P 2 ( t ) + x 2 2 P 11 ( t ) + x y P 12 ( t ) + y 2 2 P 22 ( t )
satisfy the reduced Navier–Stokes equation system and the incompressibility condition, i.e., system (25)–(27). If a rotational transformation is performed for the coordinates and velocities according to the following rule:
x x cos φ y sin φ , y x sin φ + y cos φ , V x V x cos φ V y sin φ , V y V x sin φ + V y cos φ ,
then we obtain a family of exact solutions of the form (24) with a nonlinear dependence on two horizontal coordinates (the x and y coordinates).

6. Conclusions

The paper has presented families of exact solutions to the Navier–Stokes equations for describing viscous fluid flows with rotational interaction of fluid particles taken into account. Models describing unidirectional vertical vortex flows have been separately discussed. The class proposed for such flows is characterized by a polynomial dependence of velocity on one of the horizontal coordinates. The polynomial can be of any order; the coefficients in the polynomial representation are arbitrarily related to the vertical coordinate and time. It has been demonstrated that quadratic and higher-power polynomial solutions cannot be obtained by superposition of lower-power solutions.
The solution families for describing three-dimensional flows have been considered. The Lin–Sidorov–Aristov family of exact solutions is chosen as the basic class and extended to the case of arbitrary-power dependences of the velocity field on the horizontal coordinates. The case of shearing flows has been separately discussed. It has been shown that the system of constitutive relations for spatial accelerations is reducible to a system of operator equations for which the solutions must satisfy a consistency condition of a special form. The consistency condition itself is derived in several ways.
It has also been shown that the proposed classes are a generalization of our previous results. Moreover, in the passage to the limit, which enables one to ignore the possibility of rotational interaction among fluid particles, these classes coincide with previously published results.

Author Contributions

Conceptualization, N.V.B. and E.Y.P.; methodology, N.V.B. and E.Y.P.; writing—original draft preparation, E.S.B., N.V.B. and E.Y.P.; writing—review and editing, E.S.B., N.V.B. and E.Y.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Landau, L.D.; Lifshitz, E.M. Course of Theoretical Physics. Volume 6. Fluid Mechanics; Pergamon Press: Oxford, UK, 1987. [Google Scholar]
  2. Temam, R. Navier–Stokes Equations. Theory and Numerical Analysis; North-Holland: Amsterdam, The Netherlands, 1979. [Google Scholar]
  3. Aero, E.L.; Bulygin, A.N.; Kuvshinskii, E.V. Asymmetric hydromechanics. J. App. Math. Mech. 1965, 29, 333–346. [Google Scholar] [CrossRef]
  4. Stokes, V.K. Couple stresses in fluids. Phys. Fluids. 1966, 9, 1709–1715. [Google Scholar] [CrossRef]
  5. Stokes, V.K. Theories of Fluids with Microstructure. An Introduction; Springer: Berlin, Germany, 1984. [Google Scholar] [CrossRef]
  6. Stokes, V.K. Effects of couple stresses in fluids on hydromagnetic channel flows. Phys. Fluids. 1968, 11, 1131–1133. [Google Scholar] [CrossRef]
  7. Stokes, V.K. On some effects of couple stresses in fluids on heat transfer. J. Heat Transfer. 1969, 91, 182–184. [Google Scholar] [CrossRef]
  8. Devakar, M.; Iyengar, T.K.V. Stokes’ problems for an incompressible couple stress fluid. Nonlinear Anal. Model. Control. 2008, 13, 181–190. [Google Scholar] [CrossRef]
  9. Eringen, A.C. Simple microfluids. Int. J. Eng. Sci. 1964, 2, 205–217. [Google Scholar] [CrossRef]
  10. Eringen, A.C.; Suhubi, E.S. Nonlinear theory of simple micro-elastic solids. Int. J. Eng. Sci. 1964, 2, 189–203. [Google Scholar] [CrossRef]
  11. Eringen, A.C. Theory of micropolar fluids. J. Math. Mech. 1966, 16, 1–18. [Google Scholar] [CrossRef]
  12. Toupin, R.A. Elastic materials with couple-stresses. Arch. Ration. Mech. Anal. 1962, 11, 385. [Google Scholar] [CrossRef] [Green Version]
  13. Eremeev, V.A.; Zubov, L.M. Fundam. Mech. A Viscoelastic Micropolar Fluid; SSC RAS: Rostov, Russia, 2009. (In Russian) [Google Scholar]
  14. Joseph, S.P. Some exact solutions for incompressible couple stress fluid flows. Malaya J. Mat. 2020, 648–652. [Google Scholar] [CrossRef]
  15. Cosserat, E.; Cosserat, F. Théorie des Corps Déformables; A. Hermann et Fils: Paris, France, 1909. [Google Scholar]
  16. Ahmad, F.; Nazeer, M.; Ali, W.; Saleem, A.; Sarwar, H.; Suleman, S.; Abdelmalek, Z. Analytical study on couple stress fluid in an inclined channel. Sci. Iran. 2021. [Google Scholar] [CrossRef]
  17. Srinivas, J.; Ramana Murthy, J.V. Thermal analysis of a flow of immiscible couple stress fluids in a channel. J. Appl. Mech. Tech. Phys. 2016, 57, 997–1005. [Google Scholar] [CrossRef]
  18. Lin, C.C. Note on a class of exact solutions in magneto-hydrodynamics. Arch. Ration. Mech. Anal. 1957, 1, 391–395. [Google Scholar] [CrossRef]
  19. Sidorov, A.F. Two classes of solutions of the fluid and gas mechanics equations and their connection to traveling wave theory. J. Appl. Mech. Tech. Phys. 1989, 30, 197–203. [Google Scholar] [CrossRef]
  20. Aristov, S.N. Eddy Currents in Thin Liquid Layers. Ph.D. Thesis, Institute of Automation and Control Processes, Vladivostok, Russia, 1990. (In Russian). [Google Scholar]
  21. Prosviryakov, E.Y. New class of exact solutions of Navier–Stokes equations with exponential dependence of velocity on two spatial coordinates. Theor. Found. Chem. Eng. 2019, 53, 107–114. [Google Scholar] [CrossRef]
  22. Zubarev, N.M.; Prosviryakov, E.Y. Exact solutions for layered three-dimensional nonstationary isobaric flows of a viscous incompressible fluid. J. Appl. Mech. Tech. Phys. 2019, 60, 1031–1037. [Google Scholar] [CrossRef]
  23. Burmasheva, N.V.; Prosviryakov, E.Y. Exact solution of Navier–Stokes equations describing spatially inhomogeneous flows of a rotating fluid. Tr. Instit. Mat. I Mekh. UrO RAN. 2020, 26, 79–87. (In Russian) [Google Scholar] [CrossRef]
  24. Burmasheva, N.V.; Prosviryakov, E.Y. A class of exact solutions for two-dimensional equations of geophysical hydrodynamics with two Coriolis parameters. Izv. Irkutsk. Gos. Univ. Ser. Mat. 2020, 32, 33–48. (In Russian) [Google Scholar] [CrossRef]
  25. Burmasheva, N.V.; Prosviryakov, E.Y. Isothermal layered flows of a viscous incompressible fluid with spatial acceleration in the case of three Coriolis parameters. DReaM 2020, 3, 29–46. (In Russian) [Google Scholar] [CrossRef]
  26. Öz, Y. Rigorous investigation of the Navier–Stokes momentum equations and correlation tensors. AIP Adv. 2021, 11, 055009. [Google Scholar] [CrossRef]
  27. Münch, I.; Neff, P.; Madeo, A.; Ghiba, I.-D. The modified indeterminate couple stress model: Why Yang et al.’s arguments motivating a symmetric couple stress tensor contain a gap and why the couple stress tensor may be chosen symmetric nevertheless. Z. Angew. Math. Mech. 2017, 97, 1524–1554. [Google Scholar] [CrossRef] [Green Version]
  28. Hadjesfandiari, A.R.; Dargush, G.F. Evolution of Generalized Couple-Stress Continuum Theories: A Critical Analysis. arXiv 2014, arXiv:1501.03112. [Google Scholar]
  29. Reggiani, P.; Sivapalan, M.; Hassanizadeh, S.M.; Gray, W.G. Coupled equations for mass and momentum balance in a stream network: Theoretical derivation and computational experiments. Proc. R. Soc. Lond. A. 2001, 457, 157–189. [Google Scholar] [CrossRef]
  30. Ladyzhenskaya, O.A. On nonstationary Navier–Stokes equations. Vestn. Leningr. Univ. 1958, 19, 9–18. [Google Scholar]
  31. Ladyzhenskaya, O.A. On some gaps in two of my papers on the Navier–Stokes equations and the way of closing them. J. Math. Sci. 2003, 115, 2789–2791. [Google Scholar] [CrossRef]
  32. Couette, M. Études sur le frottement des liquids. Ann. Chim. Phys. 1890, 21, 433–510. [Google Scholar]
  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. (In Russian) [Google Scholar] [CrossRef]
  34. Domnich, A.A.; Baranovskii, E.S.; Artemov, M.A. A nonlinear model of the non-isothermal slip flow between two parallel plates. J. Phys. Conf. Ser. 2020, 1479, 012005. [Google Scholar] [CrossRef]
  35. Fetecau, C.; Vieru, D.; Abbas, T.; Ellahi, R. Analytical solutions of upper convected Maxwell fluid with exponential dependence of viscosity under the influence of pressure. Mathematics 2021, 9, 334. [Google Scholar] [CrossRef]
  36. Fetecau, C.; Vieru, D. Symmetric and non-symmetric flows of Burgers’ fluids through porous media between parallel plates. Symmetry 2021, 13, 1109. [Google Scholar] [CrossRef]
  37. Aristov, S.N.; Knyazev, D.V.; Polyanin, A.D. Exact solutions of the Navier–Stokes equations with the linear dependence of velocity components on two space variables. Theor. Found. Chem. Eng. 2009, 43, 642–662. [Google Scholar] [CrossRef]
  38. Molchanov, A.M. Numerical Methods for Solving the Navier–Stokes Equations. OSF Preprints. 2019. Available online: https://osf.io/zf3j2/ (accessed on 20 June 2021).
  39. Chen, Z.J.; Li, Z.Y.; Xie, W.L.; Wu, X.H. A two-level variational multiscale meshless local Petrov–Galerkin (VMS-MLPG) method for convection-diffusion problems with large Peclet number. Comput. Fluids. 2017, 1, 1–10. [Google Scholar] [CrossRef]
  40. Li, S. Time Advancement of the Navier–Stokes equations: P-adaptive exponential methods. J. Flow Cont. Meas. Visual. 2020, 8, 63–76. [Google Scholar] [CrossRef]
  41. Li, S.-J.; Wang, Z.J.; Ju, L.; Luo, L.-S. Fast time integration of Navier–Stokes equations with an exponential-integrator scheme. In Proceedings of the AIAA Aerospace Sciences Meeting, Kissimmee, FL, USA, 8–12 January 2018. [Google Scholar] [CrossRef]
  42. Li, L.; Yang, Z.; Dong, S. Numerical approximation of incompressible Navier–Stokes equations based on an auxiliary energy variable. J. Comput. Phys. 2019, 388, 1–22. [Google Scholar] [CrossRef] [Green Version]
  43. Chen, H.; Sun, S.; Zhang, T. Energy stability analysis of some fully discrete numerical schemes for incompressible Navier–Stokes equations on staggered grids. J. Sci. Comput. 2018, 75, 427–456. [Google Scholar] [CrossRef]
  44. Dong, S. Multiphase flows of N immiscible incompressible fluids: A reduction-consistent and thermodynamically-consistent formulation and associated algorithm. J. Comput. Phys. 2018, 361, 1–49. [Google Scholar] [CrossRef] [Green Version]
  45. Ingel, L.K.; Kalashnik, M.V. Nontrivial features in the hydrodynamics of seawater and other stratified solutions. Phys. Usp. 2012, 55, 356–381. [Google Scholar] [CrossRef]
  46. Beskin, V.S. Axisymmetric steady flows in astrophysics. Phys. Usp. 2003, 46, 1209–1214. [Google Scholar] [CrossRef]
  47. Aristov, S.N.; Prosviryakov, E.Y. Inhomogeneous Couette flow. Rus. J. Nonlin. Dyn. 2014, 10, 177–182. (In Russian) [Google Scholar] [CrossRef] [Green Version]
  48. Aristov, S.N.; Prosviryakov, E.Y. Stokes waves in vortical fluid. Rus. J. Nonlin. Dyn. 2014, 10, 309–318. (In Russian) [Google Scholar] [CrossRef]
  49. Aristov, S.N.; Prosviryakov, E.Y. Nonuniform convective Couette flow. Fluid Dyn. 2016, 51, 581–587. [Google Scholar] [CrossRef]
  50. Shmyglevskii, Y.D. On isobaric planar flows of a viscous incompressible liquid. USSR Comput. Math. Math. Phys. 1985, 25, 191–193. [Google Scholar] [CrossRef]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Baranovskii, E.S.; Burmasheva, N.V.; Prosviryakov, E.Y. Exact Solutions to the Navier–Stokes Equations with Couple Stresses. Symmetry 2021, 13, 1355. https://doi.org/10.3390/sym13081355

AMA Style

Baranovskii ES, Burmasheva NV, Prosviryakov EY. Exact Solutions to the Navier–Stokes Equations with Couple Stresses. Symmetry. 2021; 13(8):1355. https://doi.org/10.3390/sym13081355

Chicago/Turabian Style

Baranovskii, Evgenii S., Natalya V. Burmasheva, and Evgenii Yu. Prosviryakov. 2021. "Exact Solutions to the Navier–Stokes Equations with Couple Stresses" Symmetry 13, no. 8: 1355. https://doi.org/10.3390/sym13081355

APA Style

Baranovskii, E. S., Burmasheva, N. V., & Prosviryakov, E. Y. (2021). Exact Solutions to the Navier–Stokes Equations with Couple Stresses. Symmetry, 13(8), 1355. https://doi.org/10.3390/sym13081355

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