Next Article in Journal
Bio-Convection Effects of MHD Williamson Fluid Flow over a Symmetrically Stretching Sheet: Machine Learning
Next Article in Special Issue
A Comparative Study of the Mode-Decomposed Characteristics of the Asymmetricity of a Vortex Rope with Flow Rate Variation
Previous Article in Journal
Modelling Symmetric Ion-Acoustic Wave Structures for the BBMPB Equation in Fluid Ions Using Hirota’s Bilinear Technique
Previous Article in Special Issue
Analysis of Heat Transfer Behavior of Porous Wavy Fin with Radiation and Convection by Using a Machine Learning Technique
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Permanent Solutions for MHD Motions of Generalized Burgers’ Fluids Adjacent to an Unbounded Plate Subjected to Oscillatory Shear Stresses

by
Constantin Fetecau
1,*,
Shehraz Akhtar
2 and
Costică Moroşanu
3
1
Section of Mathematics, Academy of Romanian Scientists, 050094 Bucharest, Romania
2
Department of Mathematics, The Islamia University of Bahawalpur, Rahim Yar Khan Campus, Rahim Yar Khan 64200, Pakistan
3
Department of Mathematics, “Alexandru Ioan Cuza” University, 700506 Iasi, Romania
*
Author to whom correspondence should be addressed.
Symmetry 2023, 15(9), 1683; https://doi.org/10.3390/sym15091683
Submission received: 16 August 2023 / Revised: 30 August 2023 / Accepted: 31 August 2023 / Published: 1 September 2023
(This article belongs to the Special Issue Symmetrical Mathematical Computation in Fluid Dynamics)

Abstract

:
Closed-form expressions have been obtained to characterize the non-dimensional velocity and corresponding non-trivial shear stress in the context of two magnetohydrodynamic (MHD) motions exhibited by incompressible generalized Burgers’ fluids. These motions occur over an infinite plate, which subjects the fluid to oscillatory shear stresses. The obtained solutions represent the first exact analytical solutions for MHD motions of such fluids under the condition of shear stress prescribed along the boundary. The establishment of these solutions relies upon the utilization of a perfect symmetry existing between the governing equations of fluid velocity and shear stress. To validate the results, a comprehensive analysis has been undertaken using two distinct methods. This validation process is further substantiated through graphical representation, demonstrating the congruence between the obtained solutions. Additionally, the convergence of the initial solutions, obtained through numerical techniques, towards their corresponding permanent counterparts has been visually established. This graphical depiction not only substantiates the accuracy of the solutions but also provides insights into the temporal evolution of the system toward its permanent state. An insight to characterize the non-dimensional shear stresses in the context of two values of the magnetic parameter is to identify that the permanent state is reached at an earlier time and the absolute magnitude of fluid velocity is reduced in the presence of an applied magnetic field.

1. Introduction

The incompressible generalized Burgers’ fluids (IGBFs), whose constitutive equations are given by the next relations [1]
T = p ^ I + S , 1 + λ 1 δ δ t + λ 2 δ 2 δ t 2 S = μ 1 + λ 3 δ δ t + λ 4 δ 2 δ t 2 A ,
represent the larger class of rate-type fluids. They contain, as special cases, the incompressible Burgers’, Oldroyd-B, Maxwell, and Newtonian fluids for λ 4   = 0 ,   λ 4 = λ 2   = 0 ,     λ 4 = λ 3 = λ 2   = 0 or λ 4   = λ 3 = λ 2 = λ 1 = 0 , respectively. In the case of the motions to be investigated here, the governing equations for incompressible second-grade fluids can also be obtained as limiting cases of the present equations. In the constitutive Equation (1), T is the Cauchy stress tensor, S is the extra stress tensor, A is the first Rivlin–Ericksen tensor, p ^ I represents the indeterminate spherical stress, μ is the dynamic viscosity of the fluid, λ 1 ,     λ 2 ,     λ 3 and λ 4 are dimensional material constants while δ / δ t represents the upper-convected time derivative.
First, the exact solutions for motions of IGBFs seem to be those of Fetecau et al. [1] in rectangular domains. In the meantime, other exact solutions for isothermal motions of the same fluids have been determined by Tong and Shan [2], Zheng et al. [3], Tong [4], Jamil [5], Khan et al. [6], and Fetecau et al. [7]. The MHD motion of fluids finds diverse uses in hydrology, horticulture, astrological exploration, aerodynamics, and the design of engineering structures. The motion of electrically conducting fluids in the presence of a magnetic field is accompanied by important effects with applications in physics, chemistry, and engineering. Exact solutions for MHD motions of IGBFs through a rectangular duct or over an infinite flat plate have been established by Sultan et al. [8], Khan et al. [9], Abro et al. [10], Alqahtani and Khan [11], and Hussain et al. [12]. However, it Is noteworthy to emphasize that the aforementioned papers address the examination of fluid motions in which velocity profiles are prescribed along the boundary. It is pertinent to acknowledge that numerous real-world practical scenarios entail the specification of shear stress along the boundaries of the flow domain [13,14].
In a seminal work, Renardy [14] demonstrated the necessity of prescribing boundary conditions for stresses at the inflow boundary in order to establish a rigorously well-posed boundary value problem for the analysis of Maxwell fluid flow. Renardy further extended this research to encompass the realm of viscoelastic fluid dynamics [15]. His investigation elucidated that the Jeffrey model maintains well-posedness within a bounded channel configuration, contingent upon the provision of the requisite components of the extra-stress tensor along the entry boundary. Moreover, the conventional “no-slip” boundary condition, while suitable for many scenarios, becomes inadequate when dealing with the motion of polymeric liquids that inherently possess the capacity to undergo boundary-sliding phenomena. Despite this, it is noteworthy that, to the extent of our current understanding, comprehensive analytical solutions encompassing MHD motions of IGBFs remain conspicuously absent within the existing body of literature. Conversely, in the framework of classical Newtonian mechanics, forces emerge as causal agents, engendering subsequent kinematic responses [16]. This contextual back-drop accentuates the significance of prescribing the shear stress distribution along the boundary, an act akin to specifying the exertion of shear forces requisite for inducing motion.
The fundamental objective of this study is to establish a pioneering framework encompassing exact, permanent solutions for MHD motions of IGBFs situated above an unbounded plate that applies oscillatory shear stress to the fluid medium. This is achieved by deft, employing a deep symmetry concerning the equations that govern fluid speed and the force of shear. In the interest of validation, the solutions are presented in dual manifestations, each form meticulously verified for equivalence through graphical analyses. Notably, these analytical solutions possess the remarkable versatility to seamlessly transition to the corresponding solutions to the fluids Burgers’, Oldroyd-B, Maxwell, second grade, and viscous underpinning the same motions. The outcomes gleaned facilitate the determination of the need for time to establish the permanent state and to underline the influence of the magnetic field on velocity. A discernible trend emerges, illustrating that the absolute magnitude of fluid velocity diminishes as the permanent state is reached at an accelerated pace with escalating values of the magnetic parameter, M. As a result, the fluid moves slower in the presence of a magnetic field.

2. Setting the Problem and Governing Equations

Let us assume that an electrical conducting IGBF is at rest over an infinite horizontal flat plate. At the moment t = 0 + , the plate begins to apply an oscillatory shear stress S cos ( ω t ) or S sin ( ω t ) to the fluid while an applied magnetic field of strength B acts vertical to the plate. Here, the two constants S and ω are the amplitude and the frequency of the oscillations, respectively. The fluid, whose magnetic Reynolds number is assumed to be small enough, is finitely conducting. Henceforth, it is reasonable to disregard the induced magnetic field and the associated Joule heating stemming from the external magnetic field. It is a well-established fact that fluids exhibiting metallic properties and ionized liquids, in particular, tend to possess a magnetic Reynolds number of diminished magnitude, as outlined in reference [17]. Moreover, in our analysis, we posit the absence of any superfluous electric charge distribution, and we deliberate the omission of Hall effects owing to the moderate levels characterizing the magnetic parameter.
Owing to the shear, the fluid is gradually moved. Since the plate is limitless, the velocity vector u corresponding to such fluid motions is of the form [1,9]
u = u ( x , t ) = u ( x , t ) k ,
where k is the unit vector along the z-direction of a convenient Cartesian coordinate system x, y, and z whose x-axis is vertical to the plate. We further posit that the extra-stress tensor S (as well as the velocity vector) is exclusively dependent on variables x and t. The condition of incompressibility is satisfied, and the momentum balance, under the absence of a pressure gradient in the z-direction, is succinctly represented by the ensuing partial differential equation [9,10]
ρ u ( x , t ) t = τ ( x , t ) x σ B 2 u ( x , t ) ;     x > 0 ,     t > 0 .
In the last relation ρ is the fluid density, τ ( x , t ) is the non-null shear stress, while σ is the electrical conductivity. Introducing the velocity vector.
u ( x , t ) from Equation (2) in (1), one obtains the next relation
1 + λ 1 t + λ 2 2 t 2 τ ( x , t ) = μ 1 + λ 3 t + λ 4 2 t 2 u ( x , t ) x ;     x > 0 ,     t > 0 ,
between the dimensional fluid velocity u ( x , t ) and the shear stress τ ( x , t ) . In the following, since we have to solve motion problems in which the shear stress is prescribed on the plate, the next initial and boundary conditions
τ ( x , 0 ) = τ ( x , t ) t t = 0 = 2 τ ( x , t ) t 2 t = 0 = 0 ;     x 0 ,
τ ( 0 , t ) = S cos ( ω t )     or     τ ( 0 , t ) = S sin ( ω t ) ,     lim x τ ( x , t ) = 0 ;     t > 0 ,
will be used. The last condition from the relations (6) assures us that there exists no shear in the free stream. We also assume that the fluid is quiescent at infinity. Consequently,
lim x u ( x , t ) = 0 ;     t > 0 .
The non-dimensional forms of the relations (3), (4), and (6), namely
u ( x , t ) t = τ ( x , t ) x M u ( x , t ) ;     x > 0 ,     t > 0 ,
1 + α t + β 2 t 2 τ ( x , t ) = 1 + γ t + δ 2 t 2 u ( x , t ) x ;     x > 0 ,     t > 0 ,
τ ( 0 , t ) = cos ( ω t )     or     τ ( 0 , t ) = sin ( ω t ) ,     lim x τ ( x , t ) = 0 ;     t > 0 ,
were obtained using the following dimensionless variables, functions, and parameter
x = x S μ ν ,     t = S μ t ,     u = u ρ S ,     τ = τ S ,     ω = μ S ω
and eliminating the star notation. In above relations ν = μ / ρ is the kinematic viscosity of the fluid while the magnetic parameter M and the non-dimensional constants α ,   β ,   γ and δ are defined by the next relations
M = ν S σ B 2 ,     α = S μ λ 1 ,     β = S μ 2 λ 2 ,     γ = S μ λ 3 ,     δ = S μ 2 λ 4 .
In the following, opposite to the usual line from the literature, we eliminate the fluid velocity u ( x , t ) between Equations (8) and (9) and obtain the following partial differential equation
1 + α t + β 2 t 2 τ ( x , t ) t = 1 + γ t + δ 2 t 2 2 τ ( x , t ) x 2 M 1 + α t + β 2 t 2 τ ( x , t ) ;     x > 0 ,       t > 0 ,
for the dimensionless shear stress τ ( x , t ) . It is worth pointing out the fact that by eliminating the shear stress τ ( x , t ) between the same Equations (8) and (9), one obtains for the dimensionless velocity field u ( x , t ) the governing equation
1 + α t + β 2 t 2 u ( x , t ) t = 1 + γ t + δ 2 t 2 2 u ( x , t ) x 2 M 1 + α t + β 2 t 2 u ( x , t ) ;     x > 0 ,       t > 0 ,
which is identical to the form with Equation (13) for the shear stress.
The two motion problems, as well as Stokes’ problems for the same fluids, become permanent or steady in time. Let us denote by u c   ( x , t ) ,     τ c   ( x , t ) and u s   ( x , t ) ,     τ s   ( x , t ) the dimensionless starting solutions corresponding to cosine or sine oscillations of shear stress on the boundary. They can be written as sums of permanent (steady state) and transient components, namely
u c ( x , t ) = u c p ( x , t ) + u c t ( x , t ) ,     τ c ( x , t ) = τ c p ( x , t ) + τ c t ( x , t ) ;     x > 0 ,     t > 0 ,
u s ( x , t ) = u s p ( x , t ) + u s t ( x , t ) ,     τ s ( x , t ) = τ s p ( x , t ) + τ s t ( x , t ) ;     x > 0 ,     t > 0 .
The starting solutions describe the fluid motion some time after its initiation. After that time, when the numerical values of the transient components u c t ( x , t ) ,     τ c t ( x , t ) or u s t ( x , t ) ,     τ s t ( x , t ) are small enough and can be neglected, the fluid behavior can be characterized by the permanent solutions u c p ( x , t ) ,     τ c p ( x , t ) or u s p ( x , t ) ,     τ s p ( x , t ) . This juncture marks a pivotal phase for the attainment of the state of permanence or equilibrium. Within the realm of practical application, this temporal juncture holds profound significance for researchers engaged in experimental endeavors seeking to ascertain the precise moment of transition of motion towards a state of equilibrium. To determine this critical temporal threshold for a given motion, the enduring solutions prove sufficient. Thus, in the subsequent section, we proffer analytic expressions exclusively pertaining to these enduring solutions. To ensure their veracity, these solutions are presented in two congruent manifestations. The requisite temporal interval for achieving the state of permanence can be derived by juxtaposing them against the initial solutions. To be more precise, this temporal threshold signifies the point at which the visual representations of the initial solutions harmonize with the corresponding constituents of enduring permanence, as seen in the diagrams.

3. Closed form Expressions for the Dimensionless Permanent Solutions

As we previously mentioned, different exact expressions will be provided for the dimensionless permanent solutions, and their equivalence will be graphically proved.

3.1. Exact Expressions for the Shear Stresses τ c p ( x , t ) , τ s p ( x , t )

The dimensionless permanent shear stresses τ c p ( x , t ) and τ s p ( x , t ) corresponding to the two motions in discussion have to satisfy the governing Equation (13) and the boundary conditions (10). Lengthy but straightforward computations show that these shear stresses can be presented in simple forms
τ c p ( x , t ) = e m x cos ( ω t n x ) ;       x > 0 ,     t R ,
τ s p ( x , t ) = e m x sin ( ω t n x ) ;       x > 0 ,     t R .
In the above relations, the non-dimensional constants m and n have the expressions
m = ω 2 a ω + ( a ω ) 2 + b 2 ( γ ω ) 2 + ( 1 δ ω 2 ) 2 ,       n = ω 2 a ω + ( a ω ) 2 + b 2 ( γ ω ) 2 + ( 1 δ ω 2 ) 2 ,
a = γ ω 2 ( 1 β ω 2 + α M ) + ( 1 δ ω 2 ) [ α ω 2 + ( 1 β ω 2 ) M ] ω 2 ,
b = ( 1 δ ω 2 ) ( 1 β ω 2 + α M ) γ [ α ω 2 + ( 1 β ω 2 ) M ] .
The first exact solution under the form (17) has been provided by Rajagopal [18] for the dimensional velocity of Stokes’ second problem corresponding to the second-grade fluids.
In order to determine equivalent expressions for these permanent shear stresses, we use the dimensionless complex shear stress
τ p ( x , t ) = τ c p ( x , t ) + i τ s p ( x , t ) ;       x > 0 ,       t R ,
where i is the complex unit. This complex shear stress has to satisfy the partial differential Equation (13) and the boundary conditions
τ p ( 0 , t ) = e i ω t ,       lim x τ p ( x , t ) = 0 ;       t R .
Because the governing Equation (13) is linear and the form of boundary conditions (23) is considered, we are in search of a solution with the following structure
τ p ( x , t ) = T ( x ) e i ω t ;       x > 0 ,       t R ,
where T ( ) is an unknown function. Direct computations show that
τ p ( x , t ) = e i ω t p x ;       x > 0 ,       t R ,
while
τ c p ( x , t ) = Re e i ω t p x   ;       x > 0 ,       t R ,
τ s p ( x , t ) = Im e i ω t p x   ;       x > 0 ,       t R .
In the above relations, the dimensionless constant p is defined by the equality
p = ( 1 β ω 2 + i ω α ) ( i ω + M ) 1 δ ω 2 + i ω γ
and simple computations clearly show that τ c p ( x , t ) and τ s p ( x , t ) given by Equations (26) and (27) satisfy the governing Equation (13) and the boundary conditions (10).
The equivalence of their expressions given by Equations (17), (26) and (18), (27), respectively, is proved by means of Figure 1.

3.2. Exact Expressions for the Velocity Fields u c p ( x , t ) , u s p ( x , t )

Once the permanent shear stresses τ c p ( x , t ) and τ s p ( x , t ) are known, we can determine the corresponding velocity fields u c p ( x , t ) and u s p ( x , t ) . They have to satisfy the governing Equations (8) and (9) and the limiting condition (7). Bearing in mind the simple forms of τ c p ( x , t ) and τ s p ( x , t ) from Equations (17) and (18), we are looking for u c p ( x , t ) an expression of the form
u c p ( x , t ) = c cos ( ω t n y ) + d sin ( ω t n y )   e m x ;       x > 0 ,       t R .
Lengthy but straightforward computations show that u c p ( x , t ) from Equality (29) satisfies the governing Equation (9) if and only if
c = ( 1 β ω 2 ) [ ω γ   n ( 1 δ ω 2 ) m ] ω α [ ω γ   m + ( 1 δ ω 2 ) n ] ( m 2 + n 2 ) [ ( ω γ ) 2 + ( 1 δ ω 2 ) 2 ]
d = ω α [ ω γ   n ( 1 δ ω 2 ) m ] + ( 1 β ω 2 ) [ ω γ   m + ( 1 δ ω 2 ) n ( m 2 + n 2 ) [ ( ω γ ) 2 + ( 1 δ ω 2 ) 2 ]
and
u c p ( x , t ) = c 2 + d 2 e m x sin ( ω t n x + φ ) ;       x > 0 ,       t R ,
where φ = arctg ( c / d ) . Similar computations also show that
u s p ( x , t ) = c 2 + d 2 e m x cos ( ω t n x + φ ) ;       x > 0 ,       t R .
Equivalent expressions for u c p ( x , t ) and u s p ( x , t ) , namely
u c p ( x , t ) = Re p i ω + M e i ω t p x ;       x > 0 ,       t R ,
u s p ( x , t ) = Im p i ω + M e i ω t p x ;       x > 0 ,       t R ,
have been obtained by introducing τ c p ( x , t ) and τ s p ( x , t ) from the relations (26) and (27) in Equation (8). Simple computations show that u c p ( x , t ) and u s p ( x , t ) given by Equations (34) and (35) satisfy all governing Equations (8), (9) and (14). The equivalence of the expressions of u c p ( x , t ) and u s p ( x , t ) given by Equations (32), (34) and (33), (35), respectively, has been proved by means of Figure 2. Furthermore, it’s readily noticeable that the dimensionless permanent solutions, which correspond to identical movements of IGBFs when magnetic influences are absent, can be directly derived by setting M = 0 in the overall solutions. These solutions, along with the previously mentioned ones, are new contributions within the existing body of literature.

4. Limiting Cases

As was already mentioned at the beginning of Section 2, similar solutions for incompressible Newtonian, second grade, Maxwell, Oldroyd-B, and Burgers’ fluids performing the same motions can be immediately obtained as limiting cases of the present solutions. In order to avoid repetition, we shall show this thing for the second-grade fluids because their constitutive equations cannot be obtained from those of IGBFs. However, for the present motions, the governing equations corresponding to these fluids can be obtained as limiting cases of present equations.

4.1. Case α = β = δ = 0 (Permanent Solutions for Second-Grade Fluids)

Making α = β = δ = 0 in Equations (17), (18), (26) and (27), one finds the following equivalent expressions,
τ S G c p ( x , t ) = e m 1 x cos ( ω t n 1 x ) ;       x > 0 ,       t R ,
  τ S G s p ( x , t ) = e m 1 x sin ( ω t n 1 x ) ;       x > 0 ,       t R
Respectively,
τ S G c p ( x , t ) = Re e i ω t p 1 x   ;       x > 0 ,       t R ,
τ S G s p ( x , t ) = Im e i ω t p 1 x   ;       x > 0 ,       t R ,
for the dimensionless permanent shear stresses τ S G c p ( x , t ) and τ S G s p ( x , t ) . They correspond to unsteady MHD motions of the incompressible second-grade fluids over an infinite plate that applies the oscillatory shear stresses S cos ( ω t ) or S sin ( ω t ) to the fluid. The expressions of these solutions are identical to those obtained by Fetecau and Morosanu [19] in which K = 0 (i.e., in the absence of porous effects). In the above relations, the new constants m 1 ,     n 1 , and p 1 are given by the relations
m 1 = ω 2 a 1 ω + ( a 1 ω ) 2 + b 1 2 1 + ( γ ω ) 2 ,       n 1 = ω 2 a 1 ω + ( a 1 ω ) 2 + b 1 2 1 + ( γ ω ) 2 ,       p 1 = i ω + M 1 + i ω γ ,
in which
a 1 = γ + M ω 2       and       b 1 = 1 γ M .
Taking now α = β = δ = 0 in the Equations (32)–(35), one obtains the corresponding permanent velocity fields
u S G c p ( x , t ) = c 1 2 + d 1 2 e m 1 x sin ( ω t n 1 x + φ 1 ) ;       x > 0 ,       t R ,
u S G s p ( x , t ) = c 1 2 + d 1 2 e m 1 x cos ( ω t n 1 x + φ 1 ) ;       x > 0 ,       t R ,
u S G c p ( x , t ) = Re p 1 i ω + M e i ω t p 1 x ;       x > 0 ,       t R ,
u S G s p ( x , t ) = Im p 1 i ω + M e i ω t p 1 x ;       x > 0 ,       t R ,
in which
c 1 = ω γ   n 1 m 1 ( m 1 2 + n 1 2 ) [ ( ω γ ) 2 + 1 ] ,       d 1 = ω γ   m 1 + n 1 ( m 1 2 + n 1 2 ) [ ( ω γ ) 2 + 1 ] ,       φ 1 = arctg ω γ   n 1 m 1 ω γ   m 1 + n 1 .
Simple computations show that the expressions of u S G c p ( x , t ) ,   u S G s p ( x , t ) from Equations (43) and (44) are identical to those obtained by Fetecau and Morosanu [19] in their Equation (61). In exchange, the expressions of these entities from Equations (41) and (42) have different forms from those from the above-mentioned reference. However, they are equivalent by transitivity. Baranovskii [20,21] has recently obtained interesting results concerning the movement of incompressible second-grade fluids within rectangular regions.

4.2. Case α = β = γ = δ = 0 (Permanent Solutions for Newtonian Fluids)

Taking α = β = γ = δ = 0 in Equations (17), (18), (32) and (33) or γ = 0 in Equations (36), (37), (41) and (42), the similar solutions
τ N c p ( x , t ) = e q x cos ( ω t r x ) ;       x > 0 ,       t R ,
τ N s p ( x , t ) = e q x sin ( ω t r x ) ;       x > 0 ,       t R ,
u N c p ( x , t ) = 1 M 2 + ω 2 4 e q x sin ( ω t r x + ψ ) ;       x > 0 ,       t R ,
u N s p ( x , t ) = 1 M 2 + ω 2 4 e q x cos ( ω t r x + ψ ) ;       x > 0 ,       t R ,
corresponding to the identical movements of incompressible Newtonian fluids are acquired. Here, the constants q, r and ψ are given by the relations
q = M + M 2 + ω 2 2 ,       r = M + M 2 + ω 2 2 ,     ψ = arctg M + M 2 + ω 2 ω .
The equivalent expressions of these entities have also simple forms, namely
τ N c p ( x , t ) = Re { e i ω t x i ω + M } ;       x > 0 ,       t R ,
τ N s p ( x , t ) = Im { e i ω t x i ω + M } ;       x > 0 ,       t R ,
u N c p ( x , t ) = Re 1 i ω + M e i ω t x i ω + M ;       x > 0 ,       t R ,
u N s p ( x , t ) = Im 1 i ω + M e i ω t x i ω + M ;       x > 0 ,       t R .
Finally, it is interesting to observe that in the absence of magnetic effects, the expressions of τ N c p ( x , t ) and τ N s p ( x , t ) from the Equations (46) and (47), respectively, are identical to the dimensionless forms of the velocity fields obtained by Erdogan [22] in the relations (12) and (17). This is possible since the shear stress in the present case and the fluid velocity in Erdogan’s paper satisfy the same governing equations and boundary conditions.

4.3. The Case ω = 0 (the Plate Applies a Constant Shear Stress S to the Fluid)

Making ω = 0 in the relations (26) and (34) or (51) and (53) one finds the expressions of dimensionless permanent shear stress and velocity fields
τ C p ( x ) = e x M ;       x > 0 ,       t R ,
  u C p ( x ) = 1 M e x M ;       x > 0 ,       t R ,
corresponding to the unsteady motion of IGBFs over an infinite flat plate that applies a constant shear stress S to the fluid. As it was to be expected, these expressions are the same both for Newtonian and non-Newtonian fluids. This is not a surprise because, as it results from the literature [23], the non-Newtonian effects disappear in time.

5. Some Numerical Results and Conclusions

The present study offers closed-form analytical expressions for dimensionless steady-state velocities and non-null shear stress distributions that pertain to isothermal MHD unsteady flows of IGBFs above an infinite flat plate subject to oscillatory shear stress conditions. The obtained solutions also encompass scenarios where the motion arises from a constant shear stress applied at the boundary, thus serving as a limiting case. Collectively, these solutions stand as the inaugural exact derivations for MHD flow behaviors exhibited by such fluids, particularly in scenarios where the shear stress is prescribed at the boundary. Moreover, these formulations can readily be tailored to yield solutions for other fluid models, such as Burgers’, Oldroyd-B, Maxwell, second grade, and Newtonian, all engaging in the same underlying motions. Notably, these derivations recover well-established solutions pertaining to second-grade fluids as limiting instances. Furthermore, akin solutions describing the identical motions of IGBFs sans the influence of magnetic effects can similarly be derived as limiting cases within the broader framework of these general solutions. These findings, noteworthy for their novelty, contribute novel insights to the existing body of scientific literature.
To validate the derived outcomes, we presented all solutions in two distinct formats, with their equivalence being visually demonstrated. In addition, as was to be expected, Figure 3 and Figure 4 clearly show that the starting solutions τ c ( x , t ) and τ s ( x , t ) (numerical solutions) converge to the corresponding permanent solutions τ c p ( x , t ) and τ s p ( x , t ) , respectively, for increasing values of the time t.
From these graphical illustrations, which have been generated for two different values of the magnetic parameter M, it can be inferred that the time required to attain the permanent state diminishes as the magnetic parameter M increases. Consequently, the presence of a magnetic field leads to an earlier achievement of the permanent state for these unsteady motions of IGBFs. In addition, the required time to touch the permanent state for motions due to sine oscillations is shorter than that for cosine oscillations of shear stress on the boundary.
Now, in order to bring to light some characteristics of the fluid motion, Figure 5 and Figure 6 have been included here. Figure 5, which presents the time variations of the permanent velocities u c p ( x , t ) and u s p ( x , t ) at the middle of the channel for increasing values of the magnetic parameter M, shows that the oscillations’ amplitude decreases in the presence of a magnetic field. Moreover, the oscillatory characteristic feature of the two movements, along with the phase discrepancy between them, can be readily noticed.
Figure 6 shows the influence of the same parameter M on the dimensionless permanent solutions τ C p ( x ) and u C p ( x ) corresponding to the unsteady MHD motion of IGBFs induced by the flat plate that applies a constant shear stress S to the fluid. It results in both the fluid velocity absolute value and the shear stress decline for increasing values of the magnetic parameter M. Consequently, the fluid moves slower in the presence of a magnetic field.
Finally, for comparison, the spatial profiles of the dimensionless starting shear stresses τ c ( x , t ) and τ s ( x , t ) (numerical solutions) are depicted adjoining in Figure 7 for α = 0.8 ,   β = 0.5 ,     γ = 0.7 ,     δ = 0.4 ,     ω = π / 6 and M = 0.9 . As before, the oscillatory behavior and the phase difference between the two motions are easily observed and the initial and boundary conditions are clearly satisfied. The three-dimensional pattern of the initial shear stresses for both cases is also depicted using two-dimensional contour plots [24] in Figure 8. This visualization is based on the same set of physical parameter values. Different colors are used to indicate their trajectories. The red and purple colors are used to bring to light the solutions’ maximum and minimum values, respectively. The oscillatory behavior of the solutions is indicated by the alternation of two distinct sets of almost closed trajectories throughout the time t. The increasing values of x lead to a decrease in the oscillation amplitude.
The main results obtained through this study are:
(1)
Concise analytical expressions have been provided for the dimensionless permanent solutions associated with unsteady MHD motions of IGBFs over an unbounded flat plate that applies oscillatory or constant shear stresses upon the fluid;
(2)
These expressions can be promptly tailored to yield comparable solutions for incompressible Burgers’, Oldroyd-B, Maxwell, second grade, and Newtonian fluids performing the same motions, and their correctness has been graphically proved;
(3)
The acquired outcomes have been employed in investigating the magnetic effects on both the steady state and fluid velocity. It was found that the permanent state is more quickly obtained, and the fluid velocity is diminished in the presence of a magnetic field;
(4)
It is pertinent to highlight that the governing Equation (14), which characterizes shear stress, exhibiting an analogous structure to Equation (13) delineating velocity, assumes pivotal significance in obtaining new exact solutions for MHD motions of IGBFs.

Author Contributions

Conceptualization, C.F. and C.M.; methodology, C.F. and C.M.; software, C.F. and S.A.; validation, C.F., C.M. and S.A.; writing—review and editing, C.F., C.M. and S.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Not applicable.

Acknowledgments

The authors wish to extend their sincere appreciation to the reviewers for their meticulous evaluation, valuable insights, and constructive recommendations pertaining to the initial version of the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

TCauchy stress tensor
A First Rivlin–Ericksen tensor
IIdentity tensor
p ^ Hydrostatic pressure
BMagnitude of the applied magnetic field
x , y , z Cartesian coordinates
u ( x , t ) Fluid velocity
MMagnetic parameter
λ 1 ,   λ 2 ,   λ 3 ,   λ 4 Dimensional material constants
α , β ,   γ ,   δ Non-dimensional material constants
μ Dynamic viscosity
ρ Fluid density
ν Kinematic viscosity
τ ( x , t ) Shear stress
ω Frequency of oscillations
σ Electrical conductivity

References

  1. Fetecau, C.; Hayat, T.; Corina, F. Seady-state solutions for some simple flows of generalized Burgers fluids. Int. J. Non-Linear Mech. 2006, 41, 880–887. [Google Scholar] [CrossRef]
  2. Tong, D.; Shan, L. Exact solutions for generalized Burgers’ fluid in an annular pipe. Meccanica 2009, 44, 427–431. [Google Scholar] [CrossRef]
  3. Zheng, L.C.; Zhao, F.F.; Zhang, X.X. An exact solution for an unsteady flow of a generalized Burgers’ fluid induced by an accelerating plate. Int. J. Nonlinear Sci. Numer. Simul. 2010, 11, 457–464. [Google Scholar] [CrossRef]
  4. Tong, D. Starting solutions for oscillating motions of a generalized Burgers’ fluid in cylindrical domains. Acta Mech. 2010, 214, 395–407. [Google Scholar] [CrossRef]
  5. Jamil, M. First problem of Stokes for generalized Burgers’ fluids. Int. Sch. Res. Not. 2012, 2012, 831063. [Google Scholar] [CrossRef]
  6. Khan, I.; Hussanan, A.; Salleh, M.Z.; Tahar, R.M. Exact solutions of accelerated flows for a generalized Burgers’ fluid, I: The case. In Proceedings of the 4th International Conference on Computer Science and Computational Mathematics (ICCSCM 2015), Langkawi, Malaysia, 7–8 May 2015; pp. 47–52. [Google Scholar]
  7. Fetecau, C.; Corina, F.; Akhtar, S. Permanent solutions for some axial motions of generalized Burgers fluids in cylindrical domains. Ann. Acad. Rom. Sci. Ser. Math. Appl. 2015, 7, 271–284. [Google Scholar]
  8. Sultan, Q.; Nazar, M.; Ali, U.; Ahmad, I. On the flow of generalized Burgers’ fluid induced by sawtooth pulses. J. Appl. Fluid Mech. 2015, 8, 243–254. [Google Scholar] [CrossRef]
  9. Khan, M.; Malik, R.; Anjum, A. Exact solutions of MHD second Stokes’ flow of generalized Burgers fluid. Appl. Math. Mech. 2015, 36, 211–224. [Google Scholar] [CrossRef]
  10. Abro, K.A.; Hussain, M.; Baig, M.M. Analytical solution of magnetohydrodynamics generalized Burgers’ fluid embedded with porosity. Int. J. Adv. Appl. Sci. 2017, 4, 80–89. [Google Scholar] [CrossRef]
  11. Alqahtani, A.M.; Khan, I. Time-dependent MHD flow of non-Newtonian generalized Butgers’ fluid (GBF) over a suddenly moved plate with generalized Darcy’s law. Front. Phys. 2020, 7, 214. [Google Scholar] [CrossRef]
  12. Hussain, M.; Qayyum, M.; Sidra, A. Modeling and analysis of MHD oscillatory flows of generalized Burgers’ fluid in a porous medium using Fourier transform. J. Math. 2022, 2022, 2373084. [Google Scholar] [CrossRef]
  13. Renardy, M. Inflow boundary condition for steady flow of viscoelastic fluids with differential constitutive laws. Rocky Mt. J. Math. 1988, 18, 445–453. Available online: http://www.jstor.org/stable/44237133 (accessed on 30 August 2023). [CrossRef]
  14. Renardy, M. An alternative approach to inflow boundary conditions for Maxwell fluids in three space dimensions. J. Non-Newton. Fluid Mech. 1990, 36, 419–425. [Google Scholar] [CrossRef]
  15. Renardy, M. Recent advances in the mathematical theory of steady flow of viscoelastic fluids. J. Non-Newton. Fluid Mech. 1988, 1, 11–24. [Google Scholar] [CrossRef]
  16. Rajagopal, K.R. A new development and interpretation of the Navier-Stokes fluid which reveals why the “Stokes Assumption” is inapt. Int. J. Non-Linear Mech. 2013, 50, 141–151. [Google Scholar] [CrossRef]
  17. Cramer, K.R.; Pai, S.I. Magnetofluid Dynamics for Engineers and Applied Physicists; McGraw-Hill: New York, NY, USA, 1973. [Google Scholar]
  18. Rajagopal, K.R. A note on unsteady unidirectional flows of a non-Newtonian fluid. Int. J. Non-Linear Mech. 1982, 17, 369–373. [Google Scholar] [CrossRef]
  19. Fetecau, C.; Morosanu, C. Influence of magnetic field and porous medium on the steady state and flow resistance of second grade fluids on an infinite plate. Symmetry 2023, 15, 1269. [Google Scholar] [CrossRef]
  20. Baranovskii, E.S.; Artemov, M.A. Steady flows of second grade fluids in a channel. Vestn. St. Petersburg Univ. Appl. Math. Comput. Sci. Control Process 2017, 13, 342–353. (In Russian) [Google Scholar] [CrossRef]
  21. 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]
  22. Erdogan, M.E. A note on an unsteady flow of a viscous fluid due to an oscillating plane wall. Int. J. Non-Linear Mech. 2000, 35, 1–6. [Google Scholar] [CrossRef]
  23. Joseph, D.D. Fluid Dynamics of Viscoelastic Liquids; Springer: New York, NY, USA, 1990. [Google Scholar]
  24. Fullard, L.A.; Wake, G.C. An analytical series solution to the steady laminar flow of a Newtonian fluid in a partially filled pipe, including the velocity distribution and the dip phenomenon. IMA J. Appl. Math. 2015, 80, 1890–1901. [Google Scholar] [CrossRef]
Figure 1. Diagrams of shear stresses τ c p ( x , t ) and τ s p ( x , t ) given by Equations (17), (26) and (18), (27), respectively, for α = 0.8 ,     β = 0.5 ,     γ = 0.7 ,   δ = 0.4 ,     ω = π / 6 ,     M = 0.9 and t = 10 .
Figure 1. Diagrams of shear stresses τ c p ( x , t ) and τ s p ( x , t ) given by Equations (17), (26) and (18), (27), respectively, for α = 0.8 ,     β = 0.5 ,     γ = 0.7 ,   δ = 0.4 ,     ω = π / 6 ,     M = 0.9 and t = 10 .
Symmetry 15 01683 g001
Figure 2. Diagrams of the velocities u c p ( x , t ) and u s p ( x , t ) given by Equations (32), (34) and (33), (35), respectively, for α = 0.8 ,     β = 0.5 ,     γ = 0.7 ,   δ = 0.4 ,     ω = π / 6 ,     M = 0.9 and t = 10 .
Figure 2. Diagrams of the velocities u c p ( x , t ) and u s p ( x , t ) given by Equations (32), (34) and (33), (35), respectively, for α = 0.8 ,     β = 0.5 ,     γ = 0.7 ,   δ = 0.4 ,     ω = π / 6 ,     M = 0.9 and t = 10 .
Symmetry 15 01683 g002
Figure 3. Convergence of the starting shear stress τ c ( x , t ) (numerical solution) to its permanent component τ c p ( x , t ) for α = 0.8 ,     β = 0.5 ,     γ = 0.7 ,   δ = 0.4 ,     ω = π / 6 , two values of M and increasing values of the time t.
Figure 3. Convergence of the starting shear stress τ c ( x , t ) (numerical solution) to its permanent component τ c p ( x , t ) for α = 0.8 ,     β = 0.5 ,     γ = 0.7 ,   δ = 0.4 ,     ω = π / 6 , two values of M and increasing values of the time t.
Symmetry 15 01683 g003
Figure 4. Convergence of the starting shear stress τ s ( x , t ) (numerical solution) to its permanent component τ s p ( x , t ) for α = 0.8 ,     β = 0.5 ,     γ = 0.7 ,   δ = 0.4 ,     ω = π / 6 , two values of M and increasing values of the time t.
Figure 4. Convergence of the starting shear stress τ s ( x , t ) (numerical solution) to its permanent component τ s p ( x , t ) for α = 0.8 ,     β = 0.5 ,     γ = 0.7 ,   δ = 0.4 ,     ω = π / 6 , two values of M and increasing values of the time t.
Symmetry 15 01683 g004
Figure 5. The time variations of the midplane permanent velocities u c p ( 0.5 , t ) and u s p ( 0.5 , t ) for α = 0.8 ,     β = 0.5 ,     γ = 0.7 ,     δ = 0.4 ,     ω = π / 6 ,     x = 0.5 and three values of M.
Figure 5. The time variations of the midplane permanent velocities u c p ( 0.5 , t ) and u s p ( 0.5 , t ) for α = 0.8 ,     β = 0.5 ,     γ = 0.7 ,     δ = 0.4 ,     ω = π / 6 ,     x = 0.5 and three values of M.
Symmetry 15 01683 g005
Figure 6. Configurations of the permanent shear stress and velocity fields τ C p ( x ) and u C p ( x ) given by Equations (55) and (56), respectively, for three values of M.
Figure 6. Configurations of the permanent shear stress and velocity fields τ C p ( x ) and u C p ( x ) given by Equations (55) and (56), respectively, for three values of M.
Symmetry 15 01683 g006
Figure 7. Spatial profiles of the starting shear stresses τ c ( x , t ) and τ s ( x , t ) (numerical solutions) for α = 0.8 ,     β = 0.5 ,     γ = 0.7 ,   δ = 0.4 ,     ω = π / 6 , and M = 0.9 .
Figure 7. Spatial profiles of the starting shear stresses τ c ( x , t ) and τ s ( x , t ) (numerical solutions) for α = 0.8 ,     β = 0.5 ,     γ = 0.7 ,   δ = 0.4 ,     ω = π / 6 , and M = 0.9 .
Symmetry 15 01683 g007
Figure 8. Contour profiles of the starting shear stresses τ c ( x , t ) and τ s ( x , t ) (numerical solutions) for α = 0.8 ,     β = 0.5 ,     γ = 0.7 ,   δ = 0.4 ,     ω = π / 6 and M = 0.9 .
Figure 8. Contour profiles of the starting shear stresses τ c ( x , t ) and τ s ( x , t ) (numerical solutions) for α = 0.8 ,     β = 0.5 ,     γ = 0.7 ,   δ = 0.4 ,     ω = π / 6 and M = 0.9 .
Symmetry 15 01683 g008
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

Fetecau, C.; Akhtar, S.; Moroşanu, C. Permanent Solutions for MHD Motions of Generalized Burgers’ Fluids Adjacent to an Unbounded Plate Subjected to Oscillatory Shear Stresses. Symmetry 2023, 15, 1683. https://doi.org/10.3390/sym15091683

AMA Style

Fetecau C, Akhtar S, Moroşanu C. Permanent Solutions for MHD Motions of Generalized Burgers’ Fluids Adjacent to an Unbounded Plate Subjected to Oscillatory Shear Stresses. Symmetry. 2023; 15(9):1683. https://doi.org/10.3390/sym15091683

Chicago/Turabian Style

Fetecau, Constantin, Shehraz Akhtar, and Costică Moroşanu. 2023. "Permanent Solutions for MHD Motions of Generalized Burgers’ Fluids Adjacent to an Unbounded Plate Subjected to Oscillatory Shear Stresses" Symmetry 15, no. 9: 1683. https://doi.org/10.3390/sym15091683

APA Style

Fetecau, C., Akhtar, S., & Moroşanu, C. (2023). Permanent Solutions for MHD Motions of Generalized Burgers’ Fluids Adjacent to an Unbounded Plate Subjected to Oscillatory Shear Stresses. Symmetry, 15(9), 1683. https://doi.org/10.3390/sym15091683

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