# Three-Dimensional Fluctuating Flow of a Second-Grade Fluid along an Infinite Horizontal Plate with Periodic Suction

Authors: Mehwish Zafar; Muhammad Afzal Rana; Atifa Latif
Journal: Mathematical Problems in Engineering
 (2022)

---

## Abstract

Laminar flow control plays a vital role in reducing drag resulting enhancement in the vehicle power by a significant amount. Theoretical and experimental studies have proven that the transition from laminar to turbulent flow (causing the drag coefficient to enhance) may be delayed/prevented by the suction of the fluid from the boundary layer to the wall. The purpose of this work is to study the effects of periodic suction on the unsteady 3-dimensional fluctuating flow of a second grade incompressible fluid flowing laminarly over a horizontal porous infinite plate. The plate is subjected to a sinusoidal transverse suction velocity while the free stream velocity oscillates in time about a constant mean. The flow turns out to be 3-dimensional because of transverse direction of the periodic suction velocity. By the series expansion method, analytic expressions for transient velocity, skin friction components, and pressure are attained. Impact of nondimensional parameters evolving in the mathematical model of the problem on these physical quantities are visualized graphically and discussed analytically. The velocity componentu is found to be rising with a growth in suction parameter α. The pressure is noted to be growing with a growth in Reynolds number Re. Further, due to suction at the plate and transient free stream velocity, the pressure increases near the plate and then reaches to a steady value far-off the plate. The skin friction along the main flow is noticed to be decreasing with the rise in α for 0.5≤Re≤1.4, however, an exponential rise is observed for Re>1.4. The skin friction along the cross flow is noted to be enhancing for the rising values of suction parameter and elastic parameter. The present results have excellent agreement with previous published results in the limiting sense. This study is beneficial in designing and manufacturing laminar flow control system to enhance the vehicle power requirement.


---

## Body

## 1. Introduction

Studies are in progress to develop a laminar control system for the purpose of tumbling vehicle power requirements [1]. It has been proved experimentally and theoretically that the laminarization of the boundary layer over a profile lessens the drag resulting enhancement in the vehicle power requirement by a very substantial amount. For the sake of artificially controlling the boundary layer, many methods have been developed since World War II [1]. Theoretical and experimental studies have proven that the transition from laminar to turbulent flow (causing the drag coefficient to enhance) may be delayed/prevented by the suction of the fluid from the boundary layer to the wall. Geometry and configuration of the slits through which the suction takes place is very important and has been deliberated broadly by various scholars. Holes of numerous sizes and slots have received extensive consideration. There is an impervious segment between the suction openings for all configurations, and it is hypothetically apposite to imagine such an area with openings as a region of changing continuous suction. The transverse sinusoidal suction velocity approach is one of the possible suction distributions. This work is oriented with the behavior of such a sinusoidally changing suction distribution over a surface and its impact on the flow.The oscillatory flows are ever momentous from the technological point of view and hence possess various real-world applications. Oscillatory unsteady flows, therefore, play a vital role in aerospace technology, chemical engineering, and turbo machinery. Lighthill [2] initiated the study of such flows in 1954 and presented the free stream oscillations effects on an incompressible viscous fluid flow over an infinite plane. In 1955, Stuart [3] extended further this work for an oscillatory 2-dimensional laminar flow over an infinite porous plane with constant suction. It was observed that the flow in the neighbor of the plate is reversed constructing the boundary layer to detach from it. Further, in 1966, Messiha [4] analyzed oscillatory two-dimensional flow when suction depending on time is subjected to the plane. Soundalgekar [5] investigated the effects of free convection on the flow past an infinite vertical oscillatory plate with the wall temperature in 1979. Raptis [6], in 1985, explored the unsteady flow problem through a porous medium circumscribed by a porous infinite plane subjected to a temperature (variable) but on a constant suction. Singh and Rana [7] performed study of the heat transfer and a viscous fluid flow through a porous medium circumscribed by a plane at constant temperature and presented effects of periodic suction velocity on the heat transfer in 1992. Singh [8] analyzed an incompressible viscous fluid flow through a porous medium circumscribed by a porous infinite plane in the presence of periodic suction with free stream velocity oscillating with respect to time about a nonzero constant mean in 1992. Considerable prevention of separation in the porous medium was noted. Later on, in 1993, Helmy [9] explored the electrically conducting fluid flow across an infinite plate wall with variable suction. In 2006, Guria and Jana [10] analyzed Couette unsteady flow of an incompressible viscous fluid when the upper plate (moving uniformly) is subjected to injection (uniform) while the lower plate (stationary) is subjected to a periodic suction. Decrease in flow velocity (main) with an increase in frequency parameter was noted. However, a rise in cross flow velocity magnitude with a grow in frequency parameter was presented. In 2010, Ahmed [11] studied unsteady three-dimensional flow of an incompressible viscous fluid past a porous vertical infinite plate with heat transfer and viscous dissipation. Workers [12–14] scrutinized transient 3-dimensional flow of viscous fluid under different physical conditions. Gersten and Gross [15] examined heat flow along a plane wall with periodic suction under the asymptotic flow conditions for downstream and determined the wall shear stress components. The researchers [16–18] analyzed technically the fluctuating flow of some Newtonian fluids through porous mediums with different rotating geometries, and they established adequate results. Mabood et al. [19] examined the effects of thermal radiation and melting heat transfer in a stagnation point flow towards a shrinking/stretching surface. Numerical solutions based on the Runge–Kutta fourth-fifth order method for the resulting nonlinear problems were obtained.All above mentioned studies have been performed for Newtonian fluids. However, majority of fluids used in industries are non-Newtonian ones. Therefore, without the study of non-Newtonian fluids, the problems associated with industries issues remain unaddressed. Shoaib et al. [20] dissected steady 3-D flow of a second grade fluid over an infinite horizontal plane subjected to a sinusoidal transverse suction velocity. They found the existence of backflow. Rana and Latif [21] presented 3-dimensional free convective second-grade fluid flow with periodic permeability and heat transfer through a porous medium. They noted a rise in pressure because of fluid thickening. Bhatti et al. [22] analyzed electro-osmotic flow of Jeffery fluid in the presence of small particles moving in the sinusoidal form in a Darcy–Brinkman–Forchheimer medium and transversely applied magnetic field. Riaz et al. [23] presented the exact solutions of non-Newtonian multiphase fluid flow through peristaltic pumping characteristics in an annulus having complaint walls and applied magnetic field. It was concluded that applied magnetic field decreases the velocity of both the fluid and the particles flow. Gireesha et al. [24] analyzed the effect of nanoparticles on a steady MHD flow and heat and mass transfer of Eyring–Powell fluid in two lateral directions over a convectively heated stretching sheet. Similarity transformations were employed to reduce the governing partial differential equations into a set of nonlinear ordinary differential equations which were solved numerically using the fourth-fifth order Runge–Kutta–Fehlberg method with the shooting technique. Mahanthesh et al. [25] scrutinized numerically the combined effects of nonlinear thermal convection and radiation in 3D boundary layer flow of Maxwell nanofluid over a stretching sheet. The solutions were computed via the homotopy procedure. It was noted that the nanoparticle volume fraction and temperature profiles are stronger for the case of solar radiation in comparison with the problem without radiation. Shafiq et al. [26] examined the second grade bioconvective nanofluid flow with the buoyancy effect and chemical reaction. The Brownian motion and thermophoretic mechanisms along with Newtonian heating were also considered. Mass flux was found to be enhancing function of both Brownian diffusion parameter and Lewis number. Mabood et al. [27] studied rheological aspects on the micropolar fluid model by simulating the reactive flow from a continuously moving flat plate. The model was developed employing Brownian motion, thermophoretic diffusion, and chemically reactive species. It was found that the presence of thermophoresis and Brownian motion is more effective to improve the heat transportation phenomenon. Shamshuddin et al. [28] scrutinized the magnetized flow of Casson nanofluid past a Riga surface with nonlinear radiative, uneven heat sink/source, thermophoretic movement, and chemical reaction. The analysis revealed that an increment in the Casson term causes a rise in the temperature profile for CuO and MgO nanofluid and dominant behavior is noted in case of CuO nanofluid on comparing with MgO nanofluid.Above mentioned studies were carried out for steady flows of non-Newtonian fluids; however, oscillatory unsteady flows play a vital role in aerospace technology, turbo machinery, and chemical engineering. The objective of this study is to analyze periodic suction velocity effect on the 3-D flow of a second-grade incompressible fluid flowing laminarly over a horizontal infinite plane with a free stream oscillatory velocity. To the best of authors’ information, no such investigation has been carried out. Basically, laminar flow control performs a vital role in reducing drag resulting enhancement in the vehicle power by a significant amount; therefore, the present study is relevant to various engineering sciences applications, particularly, relevant to aeronautical engineering.


## 2. Mathematical Formulation

Consider an unsteady flow of a second grade incompressible fluid along an infinite porous plate lying horizontally on thex∗z∗−plane with x∗−axis along the plate, being the direction of the flow. The y∗−axis which is directed into the fluid is along normal to the plate. All fluid physical quantities are not dependent of x∗ due to the plate infinite length in x∗−direction [11, 15]. The flow, however, ruins 3-D because of the transverse sinusoidal suction velocity [11, 15]:(1)v∗z∗=V01+ϵcosπz∗l,where V0<0 is suction velocity, l is the wave length of the suction distribution, and ∈≪1 is the amplitude of the suction variation.The continuity equation and momentum equation respectively governing the fluid flow are(2)∇.V⟶=0,(3)ρdV⟶dt=∇.τ˜,where(4)V⟶=u∗y∗,z∗,t∗,v∗y∗,z∗,t∗,w∗y∗,z∗,t∗,is velocity field, and the Cauchy stress tensor τ˜ for the second grade fluid model [29–31] is(5)τ˜=−pI˜+μA˜1+α1A˜2+α2A˜12,while Rivlin–Ericksen tensors, A˜1 and A˜2, are(6)A˜1=gradV⟶T+gradV⟶,A˜2=gradV⟶TA˜1+A˜1gradV⟶+dA˜1dt..For the model (5) to meet the assumption that the specific Helmholtz free energy is a minimum in equilibrium and to be well-matched with the thermodynamics in the sense that all motions meet the Clausius–Duhem inequality [31], then(7)α1≥0,μ≥0,α1+α2=0.The equations (2) and (3) in view of Equations (4)–(7) yield(8)∂v∗∂y∗+∂w∗∂z∗=0,(9)ρ∂u∗∂t∗+ρv∗∂u∗∂y∗+w∗∂u∗∂z∗=ρ∂U∗∂t∗+μ∂2u∗∂y∗2+∂2u∗∂z∗2+α1v∗∂3u∗∂y∗3+w∂3u∗∂z∗∂y∗2+v∗∂3u∗∂y∗∂z∗2+w∗∂3u∗∂z∗3,,(10)ρ∂v∗∂t∗+ρv∗∂v∗∂y∗+w∗∂v∗∂z∗=−∂p∗∂y∗+μ∂2v∗∂y∗2+∂2v∗∂z∗2+α1v∗∂3v∗∂y∗3+w∂3v∗∂z∗∂y∗2+v∗∂3v∗∂y∗∂z∗2+w∗∂3v∗∂z∗3+5∂v∗∂y∗∂2v∗∂y2∗+2∂u∗∂2u∗∂y∗∂y2∗+2∂w∗∂2w∗∂y∗∂y2∗+∂v∗∂2w∗∂z∗∂y2∗+∂u∗∂z∗∂2u∗∂z∗∂y∗+∂u∗∂y∗∂2u∗∂z2∗+∂v∗∂y∗∂2v∗∂z2∗+∂v∗∂z∗∂2v∗∂z∗∂y,(11)ρ∂w∗∂t∗+ρv∗∂w∗∂y∗+w∗∂w∗∂z∗=−∂p∗∂z∗+μ∂2w∗∂y∗2+∂2w∗∂z∗2+α1v∗∂3w∗∂y∗3+w∗∂3w∗∂z∗∂y∗2+v∗∂3w∗∂z∗∂y∗2+w∗∂3w∗∂z∗3+∂w∗∂y∗∂2w∗∂y∗2∂z∗+∂u∗∂2u∗∂y∗∂y∗2∂z∗+5∂w∗∂2w∗∂z∗∂z∗2+∂w∗∂2v∗∂y∗∂z∗2+2∂u∗∂2u∗∂z∗∂z2∗+2∂v∗∂2v∗∂z∗∂z2∗+∂u∗∂2u∗∂z∗∂z2∗+∂w∗∂2w∗∂z∗∂y∗2.The appropriate boundary conditions [11] are(12)u∗=0,v∗=V01+∈cosπz∗l,w∗=0aty∗=0,u∗=U∗t∗=U01+ϵein∗t∗,v∗=V0,w∗=0,p∗=p∞∗aty∗⟶∞,,where n∗ is frequency of free stream oscillation and t∗ denotes time.


## 3. Dimensionless Equations

Introducing the nondimensional variables,(13)y=y∗l,z=z∗l,u=u∗U0,v=v∗U0,w=w∗U0,U=U∗U0,t=n∗t∗,n=n∗l2ϑ,p=p∗ρU02,p∞=p∞∗ρU02,,(14)α=V0U0Suction parameter,Re=U0lϑReynolds number,K=α1ρl2Elastic parameter,,then the equations (8)–(11) become(15)∂v∂y+∂w∂z=0,(16)nRe∂u∂t+v∂u∂y+w∂u∂z=nRe∂U∂t+1Re∂2u∂y2+∂2u∂z2+Kv∂3u∂y3+w∂3u∂z∂y2+v∂3u∂y∂z2+w∂3u∂z3,(17)nRe∂v∂t+v∂v∂y+w∂v∂z=−∂p∂y+1Re∂2v∂y2+∂2v∂z2+Kv∂3v∂y3+w∂3v∂z∂y2+v∂3v∂y∂z2+w∂3v∂z3+5∂v∂y∂2v∂y2+2∂u∂y∂2u∂y2+2∂w∂y∂2w∂y2+∂u∂z∂2u∂z∂y+∂v∂z∂2w∂y2+∂u∂y∂2u∂z2+∂v∂y∂2v∂z2+∂v∂z∂2v∂z∂y,(18)nRe∂w∂t+v∂w∂y+w∂w∂z=−∂p∂z+1Re∂2w∂y2+∂2w∂z2+Kv∂3w∂y3+w∂3w∂z∂y2+v∂3w∂y∂z2+w∂3w∂z3+5∂w∂z∂2w∂z2+2∂u∂z∂2u∂z2+∂w∂y∂2w∂y∂z+∂w∂z∂2w∂y2+2∂v∂z∂2v∂z2+∂u∂z∂2u∂y2+∂u∂y∂2u∂y∂z+∂w∂y∂2v∂z2.The corresponding boundary conditions (12) reduce to(19)u=0,v=α1+ϵcosπz,w=0aty=0,u=1+ϵeint,v=α,w=0,p=p∞asy⟶∞..It is wealth mentioning that the set of nonlinear partial differential equations (16)–(18) reduce to momentum equations of [11] for non-Newtonian parameter K⟶0 and momentum equations of [21] under identical physical conditions.


## 4. Solution

As the amplitude∈≪1, we take up solutions for the equations (15)–(18) in the vicinity of the plate in the form:(20)Fy,z,t=F0y+∈F1y,z,t+∈2F2y,z,t+⋯,where F stands for any u,v,w,andp.For∈=0, the problem becomes the 2-dimensional steady flow having constant suction at the plate. Plugging in the series (20) into equations (15)–(19), the zeroth order equations in ∈ are as follows and equation (15) yields(21)dv0dy=0⇒v0=constant.The boundary conditionv00=α⇒v0=α.From equation (16), the zeroth order equation is(22)v0du0dy=1Red2u0dy2+Kv0d3u0dy3.In view ofv0=α and the boundary conditions (19), the zeroth order boundary value problem becomes(23)Kαd3u0dy3+1Red2u0dy2−αdu0dy=0,u0=0aty=0,u0=1aty⟶∞..The differential (23) has order higher than the available boundary conditions as a result of elastic parameter K. For unique solution of the differential equation, three boundary conditions are required. To resolve this issue, we expand solution u0 in the power series of K≪1 such that(24)u0=u00+Ku01+OK2.Employing equation (24) in equation (23) and comparing like terms in O(K0) and O(K) on both sides to get the zeroth order and first order boundary value problems is given by(25)d2u00dy2−αRedu00dy=0,u000=0,u00∞=1,and(26)αRed3u00dy3+d2u01dy2−αRedu01dy=0,u010=0,u01∞=0.Now, solving the boundary value problem (25), we get(27)u00y=1−eαRey.In view of solution (27), the boundary value problem (26) yields(28)d2u01dy2−αRedu01dy=αRe4eαRey,u010=0,u01∞=0.Its solution is(29)u01y=αRe3yeαRey.Hence, we get(30)u0y=1−eαRey+KαRe3yeαRey.When∈≠0, the substitution of equation (20) for u,v,w into the equations (15)–(19) and comparison of the coefficients of ∈ yields partial differential equations that is given by(31)∂v1∂y+∂w1∂z=0,(32)nRe∂u1∂t+v1∂u0∂y+α∂u1∂y=1Re∂2u1∂y2+∂2u1∂z2+in2Reeint+Kα∂3u1∂y3+α∂3u1∂y∂z2+v1∂3u1∂y3,(33)nRe∂v1∂t+α∂v1∂y=−∂p1∂y+1Re∂2v1∂y2+∂2v1∂z2+Kα∂3v1∂y3+∂3v1∂y∂z2,(34)nRe∂w1∂t+α∂w1∂y=−∂p1∂z+1Re∂2w1∂y2+∂2w1∂z2+Kα∂3w1∂y3+∂3w1∂y∂z2.Likewise, boundary conditions are(35)u1=0,v1=αcosπz,w1=0aty=0,u1=eint,v1=0,w1=0,p1=0asy⟶∞..The set of differential equations (31)–(34) describe the unsteady three-dimensional flow. The equation (32) governs the main flow while equations (31), (33), and (34) govern the cross flow.


## 5. Cross Flow, Main Flow Solutions, and Pressure

In order to solve the boundary value problems (32)–(35), we consider the solutions for u1,v1,w1,p1 in the form of following complex notations whose real parts have physical consequence:(36)u1y,z,t=u11yeint+u12ycosπz,(37)v1y,z,t=v11yeint+v12ycosπz,(38)w1y,z,t=−zv11′yeint+1πv12′ysinπz,(39)p1y,z,t=p11yeint+p12ycosπz.We note that the velocity components (37) and (38) identically satisfy the continuity equation (31). Moreover, prime denotes differentiation with respect to y. Setting these equations into equations (32)–(35) and solving the resulting partial differential equations to get the solutions for u1,v1,w1, and p1 given by(40)u1y,z,t=1−e−βyeint+α2Re2π−λλπαRe−π2λαRee−λy+π2λαRee−λ−αRey−λπαRee−π−αRey+KFe−λy+A−Ce−λ−αRey+D−B+Eye−π−αRey−2Eye−λ−αReycosπz,(41)v1y,z,t=απ−λπe−λy−λe−πycosπz,(42)w1y,z,t=αλπ−λe−λy−e−πysinπz,(43)p1y,z,t=αλReππ−λKαReλλ2−π2e−λy+α2λπ−λe−πycosπz,where(44)λ=−αRe2+αRe22+π2,β=−αRe2+αRe22+in2,A=αRe4π/π−λ+αRe3παRe−λπ2−αRe−λ2/π−λ2λαReλ−αRe2−αReαRe−λ−π2,B=αRe4π/π−λ+αRe3λαRe−παRe−π2−π2/π−λαReππ−αRe2−αReαRe−π−π2,C=παRe5αRe−2λπ−λ,D=λαRe5αRe−2ππ−λ,E=πλαRe6π−λ,F=−A+B+C−D.In view of equations (30) and (40), we get(45)uy,z,t=1−eαRey+KyαRe3eαRey+∈1−e−βyeint+α2Re2π−λe−λyπ2λαRee−λ−αRey−λπαRee−π−αRey+KFe−λy+A−Ce−λ−αRey+D−B+Eye−π−αRey−2Eye−λ−αReycosπz.


## 6. Shear Stress Components

Shear stress components inx− and z− directions, respectively, are(46)Cfx=∂u/∂yy=0αRe,or(47)Cfx=−1+KαRe2+ϵβeintαRe+ϵ−λπ−π2λ+παRe−λ2λ−λαRe−ππλπ−λ+ϵKαRe−λ−A+B+C−D+A−CαRe−λ+D−BαRe−π−E,and(48)Cfz=μ∂w/∂yy=0α,or(49)Cfz=μλsinπz.


## 7. Results and Discussions

In the present study, a 3-dimensional transient flow of a second grade fluid over a porous plane having periodic suction and time dependent free stream velocity is modeled and analyzed theoretically. Analytic expressions for the velocity field, pressure, and skin friction are computed by the perturbation technique. On these physical quantities, the effects ofRe (Reynolds number), K (elastic parameter), a (suction parameter), t (time parameter), and n (frequency parameter) are visualized graphically. For the sake of the physical inside into the problem, the velocity field, skin friction components, and pressure are conferred by passing on numerical values to numerous nondimensional parameters arising in the mathematical model of the problem.Figure1 displays the effect of Re, n, t, suction parameter a, and elastic parameter K on the transient velocity u. In Figures 1(a)–1(e), one of the parameters varies while all other nondimensional parameters are fixed. In Figure 1(a), Re = 8, 10, 15, whereas other parameters are fixed i. e., t = 0.5, a = −0.1, n = 5, z = 0, K = 0.05, and ε = 0.05. It is detected that with the upsurge of the Re, the velocity component u also upsurges and grasps its maximum value at the boundary for each value of the Re is taken in this regard. Physically it meant that inertial forces are overriding the viscous forces. Moreover, the boundary layer thickness declines with a growth in Re. Similarly, it is perceived that the velocity component u (Figure 1(b)) increases with a rise in n. This upshot is very effective in the vicinity of the plate. Also, this is motivating to note that all velocity profiles tend to unity, away from the plate. Figure 1(c) displays the influence of the transient number on u. It discloses that an upsurge in time parameter results in a reduction in u. This is quite in agreement with the physical fact that unsteady flow approaches to steady flow after large time. It is evident that u (Figure 1(d)) rises with a growth in suction parameter a. In general, a rise in a leads to decelerate the flow velocity. Of course, in this case, an enhancement in u is because of the transient free stream velocity. Figure 1(e) parades the impact of K on u. It is detected that a growth in K leads to a reduction in u, which looks correct physically as an increase in K results in fluid thickening causing reduction in the velocity.Figure 1
Impact of all nondimensional parameters on main flow. (a) Impact ofRe on u. (b) Impact of n on u. (c) Impact of t on u. (d) Impact of α on u. (e) Impact of K on u.
(a)(b)(c)(d)(e)Impact of Reynolds numberRe, frequency parameter n, time parameter t, suction parameter α, and elastic parameter K on u is presented in Figure 2 for z as abscissa. It is observed that a growth in Re (Figure 2(a)) leads to enhance the amplitude of oscillation and consequently enhancement in main flow which is in agreement with the results discussed in the Figure 1(a). The similar behavior is noted for the suction parameter on u (Figure 2(d)). It is seen from Figure 2(b) that enhancement in the values of n only shifts the location of oscillation. Similar impact of t on u is noted (Figure 2(c)). It is examined from Figure 2(e) that larger value of k rises the amplitude of oscillation.Figure 2
Impact of nondimensional parameters on main flow forz as abscissa. (a) Impact of Re on u. (b) Impact of n on u. (c) Impact of t on u. (d) Impact of α on u. (e) Impact of K on u.
(a)(b)(c)(d)(e)Figure3 illustrates the impact of Re and α on v and w. The magnitude of v shrinks with a growth in Re (Figure 3(a)). Instead, magnitude of v increases for growing values of α, and this influence is significant near the plate (Figure 3(c)), and it appears physically true, since it is well known that the suction impact on the flow is maximum adjacent to the plate. Further, the velocity component v becomes stable as y⟶∞. It is detected from Figure 3(b) that for a fixed value of Re, w rises exponentially adjacent to the plate, achieves its maximum value, and then declines rapidly, and ultimately approaches to 0 as y⟶∞. Adjacent to the plate, a parabolic profile is attained. This velocity component growths with a rise in Re adjacent to the plate but possesses opposite trend away from the plate. Physically, it reflects the dominance of inertial forces over viscous forces adjacent to the plate. However, away from the plate, viscous forces dominate the inertial forces due to oscillatory free steam velocity resulting opposite trend of velocity component w and eventually approaching to zero far away from the plate. It is viewed from Figure 3(d) that, for a fixed value of a, w growths exponentially adjacent to the plate, achieves its optimum (maximum value), then shrinkages rapidly, and lastly approaches to zero as y⟶∞. Adjacent to the plate, a parabolic profile is attained. It is also observed that the suction parameter enhances the velocity component near the plate which is against the general conception (suction parameter leads to decelerate velocity). An enhancement, in fact, in the velocity component near the plate is because of the transient free stream velocity.Figure 3
Impact of nondimensional parameters onv and w. (a) Impact of Re on v (b) Impact of Re on w. (c) Impact of α on v. (d) Impact of α on w.
(a)(b)(c)(d)Figure4 exemplifies the impact of Re and α on v and w against z as abscissa. It is perceived from Figure 4(a) that an increase in Re causes a decline in the amplitude of oscillation resulting reduction in the flow velocity in the y− direction. Further, the impact of Re on w (Figure 4(b)) is to upsurge the amplitude of oscillation and consequently enhancement in velocity component w agree with the results discussed previously (Figure 3(b)). From Figure 4(d), similar influence of α on w is noted. The behavior of α on v is shown in Figure 4(c). It is found that increasing α just results in shift of the axis of oscillation.Figure 4
Impact of nondimensional parameters on v and w forz as abscissa. (a) Impact of Re on v. (b) Impact of Re on w. (c) Impact of α on v. (d) Impact of α on w.
(a)(b)(c)(d)Figure5 displays the effect Re,α, and K on the pressure. In Figure 5(a), Reynolds number is varied, that is, Re = 8, 10, 15 and other existing parameters are fixed, i. e., α=−0.1,z = 0, K = 0.05, and ε = 0.05. The pressure grows near the plate with the rise in Re. It shows the dominance of inertial forces over the viscous forces. Of course, the pressure attains a steady value far away from the plate. In contrast, the impact of α is to decline the pressure near the plate (Figure 5(b)). Moreover, pressure converges to a steady value away from vicinity of the plate, where the suction effect vanishes, for increasing values of suction parameter. From Figure 5(c), it is perceived that the pressure decreases adjacent to the plate when K is elevated. It is consistent with the physical result that the boundary layer is thinned by suction causing reduction in pressure near the plate. Further, similar to influence of α on the pressure far away from the plate, pressure again approaches to a steady value away from vicinity of the plate for increasing values of K.Figure 5
Impact of nondimensional parameters on pressure. (a) Impact ofRe on p. (b) Impact of α on p. (c) Impact of K on p.
(a)(b)(c)Figure6 characterizes the impact of Re,α, and K on pressure against z as abscissa. It is perceived from Figure 6(a) that growth in Re results in enhancement in the amplitude of oscillation causing a rise in the pressure. Similar effect of α (Figure 6(b)) and K (Figure 6(c)) on the pressure is found.Figure 6
Impact of nondimensional parameters on pressure forz as abscissa. (a) Impact of Re on p. (b) Impact of α on p. (c) Impact of K on p.
(a)(b)(c)Figure7 typifies Cfx (nondimensional skin friction component) in the main flow direction against Re for numerous values of frequency parameter n, transient number t, a, and K. Figure 7(a) shows the growing influence on Cfx magnitude for increasing values of n. It is also revealed that there is exponential falloff in the magnitude of this skin friction component for 0.5≤Re≤1.5 up to a certain value depending upon the value of n taken in this regard and then a rapid rise in the magnitude of this friction component is observed. Physically it reflects the dominance role of viscous forces over inertial forces for the range 0.5≤Re≤1.5. Of course, it is worthwhile to point out that the role of viscous forces over inertial forces is overturned when Re>1.5. On the contrary, the role of t regarding its impact on Cfx is inverted (Figure 7(b)). Figure 7(c) portrays the effect of α upon Cfx. It is observed that the magnitude of Cfx decreases for 0.5≤Re≤1.4 and then increases exponentially for rising values of α. It shows the dominant role of viscous forces over the inertial forces in the range 0.5≤Re≤1.4 causing decay in Cfx magnitude. In contrast, for Re>1.4, inertial forces overcome the viscous forces resulting rapid grow in Cfx magnitude. In Figure 7(d), similar behavior of K on Cfx is noted for Re>1.1, and increasing values of K has no significant impact for 0.5≤Re≤1.1.Figure 7
Impact of nondimensional parameters on skin friction the alongx-axis. (a) Impact of n on skin friction along the x-axis. (b) Impact of t on skin friction along the x-axis. (c) Impact of α on skin friction along the x-axis. (d) Impact of K on skin friction along the x-axis.
(a)(b)(c)(d)Figure8 is traced for Cfz (nondimensional skin friction component) in the z-direction versus Re for various values of K and α. It is detected that a growth in Re results in growth in Cfz magnitude in both cases. Further, the enhancement in Cfz magnitude by increasing suction parameter (Figure 8(a)) seems naturally correct as the role of suction is to promote the resistance in flow causing growth in the skin friction. It is also detected that increasing elastic parameter results (Figure 8(b)) in escalation of skin friction magnitude. In fact, an increase in K results in the thickening of fluid causing rise in skin friction.Figure 8
Impact of nondimensional parameters on skin friction the along thez-axis. (a) Impact of α on skin friction along the z-axis. (b) Impact of K on skin friction along the z-axis.
(a)(b)


## 8. Remarks

This analysis aims to explore analytic solutions for 3-dimensional fluctuating flow of a second grade fluid with periodic suction. Approximate solutions foru,v,w (velocity field), Cx,Cfz (skin friction components), and pressure based upon perturbation are computed and explored theoretically. As a regular perturbation technique based upon small parameter existing in the governing equations, series solutions obtained by this method always converges. Consequently, analytic solutions for u,v,w computed by this approach in this work also converges. The conclusions of this analysis are as follows:(i)
Velocity componentu grows with Re, frequency parameter n and α, whereas it declines with increasing time parameter and non-Newtonian parameter(ii)
v increases with growing values of Re, and it decreases for growing values of α, and this impact is significant near the plate and ultimately it becomes stable as y⟶∞(iii)
w is noted to be enhancing exponentially adjacent to the plate, attaining its optimum (maximum) value, then decreasing rapidly, and ultimately approaches to 0 as y⟶∞ due to Re and α(iv)
A rise inRe, n, a, and K leads to enhance the oscillation amplitude and consequently enhancement in u whereas enhancement in transient parameter causes an inverse effect(v)
The amplitude of oscillation of cross flow velocityw upsurges due to enhancement in Re and α causing intensification in the flow velocity. In contrast, Re possesses inverse impact on cross flow velocity v(vi)
The influence ofRe is to boost the pressure; however, α and K possess inverse impact on it(vii)
The effect of suction parameterα and non-Newtonian parameter K is to decrease the pressure near the plate(viii)
The growth inRe,α, and K results in enhancement in the amplitude of oscillation causing rise in the pressure(ix)
Cfx rises for growing values of n, and the role of t regarding its impact on Cfx is inverted(x)
It is observed thatCfx decreases for 0.5≤Re≤1.4 and then increases exponentially for rising values of α and K(xi)
The role of viscous forces over the inertial forces is found to be dominant for0.5≤Re≤1.4; in contrast, this role is noted to be inverted for Re>1.4(xii)
Cfz enhances for rising values α,K, and Re(xiii)
A growth inRe causes to decline boundary layer thickness(xiv)
Inertial forces are found to be dominant over viscous forces adjacent to the plate(xv)
Results of [11] are recovered for non-Newtonian parameter K⟶0(xvi)
Results of Ref. [21] are validated under identical physical conditions



---

*Source: 3509721-2022-10-29.xml*