# Analytical Solution of Slow Squeeze Flow of Slightly Viscoelastic Fluid Film between Two Circular Disks Using Recursive Approach

Authors: Muhammad Memon; Asif Ali Shaikh; Abdul Majeed Siddiqui; Laveet Kumar
Journal: Mathematical Problems in Engineering
 (2022)

---

## Abstract

This research investigates an analytical solution for the slow squeeze flow of the slightly viscoelastic fluid film between two circular disks in which the upper disk approaches the lower disk with a constant velocity, and the lower disk is kept stationary. The determination of the study is to identify the behavior of the differential type fluid on the steady squeezing flow using Langlois recursive approach. The governing equations for the axisymmetric flow are expressed in cylindrical coordinates and yield the nonlinear system of partial differential equations. The analytical solution of the resulting equations with nonhomogeneous boundary conditions is obtained by the Langlois recursive approach. Flow variables such as stream function, velocity profiles, pressure distribution, shear, normal stresses, and normal force acting on the disk are determined. These flow variables are nondimensionalized by using suitable dimensional quantities. The influence of slightly viscoelastic parameterβ, radial distance r, and aspect ratio on velocity components, pressure distribution, and normal squeeze force is examined mathematically and portrayed graphically. The results illustrate that the axial and radial velocities increase at the higher values of the slightly viscoelastic parameter β, which confirms the shear thickening behavior. The obtained solutions of the flow variables satisfied the existing solutions on squeeze flow of viscous fluid upon vanishing the slightly viscoelastic parameters. This solution could elucidate the classical lubrication problems, particularly in load and thrust bearing characteristics of the human body joints, the compression molding process of materials, etc.


---

## Body

## 1. Introduction

The squeezed flow has many real-life applications, such as dampers, motor bearing, lubrication, chewing between teeth or gums, and the compression molding processes of metals and polymers (filled or unfilled). The valves and arthritic joints are interesting real-world applications in biology and biotechnology of essential squeezed flows [1]. The entire synovial fluid in the human knee joint is not pumped immediately from the two sides of the joint during the heel-strike processes and toe-off period. In the presence of liquid viscous resistance, the contact surfaces require a specific period. During this interval, the pressure is generated, and the squeezed fluid film supports the force [2]. Several flows are often found in traditional lubrication products, such as cams, engine connecting rod bearings, and gears. In such cases, viscoelastic additives are capable of increasing the load-bearing strength of lubricants [3].The study of squeezing flow can be traced down back to the 19th century. The work by Stefan [4], Reynolds [5], and Scott [6] can be considered as pioneering work. Earlier experimental works on squeezing flow between two circular disks were reported in [7–9]. In the beginning, creeping flow, limiting to Newtonian fluids, was studied using nonclassical lubrication approximation. Later, inertia was included in the Newtonian creeping flow by Jackson [10], and Kuzma [11] identified the error in the inertia term of the work done by Jackson [10]. Jones and Wilson [12] found the squeezing force of liquids with inertia for large Reynolds values. Numerical solutions for squeezing flow between parallel disks have been conducted by Hamza and Macdonald [13].Further experimental results for slow squeezing are highlighted in power-law fluids, while rheological models are applied for fast squeezing flow, which describes the overshoot phenomena [14]. The applied forces and spaces between the disks in the squeezing flow of the power-law fluid were found by Lieder [14, 15] in the solution of shear-thinning polymers. From these results, it has been observed that force and distance parameters deviate from the Scott equation. Later on, McClelland and Finlayson [16] extended Oliver’s work and proposed a model that agrees with experimental data at fast squeezing rates by combining the effects of normal stress. The experimental results on a high viscosity of low-density polyethylene in constant load tests were compared with Criminale–Ericksen–Filbey fluid and found a close relationship between the results [3]. However, it has been analyzed that results produced in the case of slow squeezing are better as compared to fast squeezing. Kramer [17] repeated Tanner’s work and used only the lodge rubber-like fluid Maxwell model. Phan–Thien and Tanner [18] extended the work of [9, 19] by including the inertia effect but neglecting the body forces and edge effects of the motion of the nonlinear Maxwell model. Phan–Thien et al. [20] found the solution of creeping squeeze film flow of inelastic fluids such as the Carreau model and modified Phan–Thien–Tanner (MPTT) model by numerical method considering stress overshoot phenomena. Yousefi et al. [2] assumed the constitutive equation of the MPTT model by considering a synovial fluid and showed a valuable contribution to knee joint lubrication problems. Lee et al. [9] obtained a finite element solution to squeeze a convicted Maxwell fluid under the constant force on the upper disk with the inertia of the fluid. Muravleva has studied the squeeze flow of the Bingham fluid [21–23] in the plane and axisymmetric geometries by the numerical technique, that is, augmented Lagrangian method and asymptotic solution. After that, Singeetham and Puttana [24] extended the work of Muravleva and found the analytical solution by matching the asymptotic expansion technique of plane squeeze flow of Herschel–Bulkley and Casson fluid between two disks. The solution is divided into three regions, that is, shear stress dominant, pseudo-plug, and central pseudo-plug plastic regions due to the yielded stress. The behavior of squeezing flow of MHD Casson fluid with slip and no-slip conditions is investigated in [25–27]. The behavior of Newtonian [28–30] and non-Newtonian fluid [31–34] on different geometries is examined numerically and analytically. The theoretical and experimental researches of squeezing flows have been studied by many researchers [35, 36].This research study explores the squeeze flow of the differential type fluids due to their significant features of rod climbing, shear thinning, shear thickening effects, and normal stress. The constitutive equation of a special class of third-order fluid named the slightly viscoelastic fluid has applications in journal bearing and slide bearing [37, 38]. Governing equations of the squeeze flow of slightly viscoelastic fluid film are the nonlinear system of partial differential equations in the axisymmetric form with nonhomogeneous boundary conditions. Such types of differential equations are complex to solve numerically. Therefore, in this research investigation, the recursive approach of Langlois [39] has been used to linearise these equations, and the analytical solution has been obtained on specified boundary conditions.The aim of this research is to deliver the analytical solution of squeeze flow of steady incompressible viscoelastic fluid between two circular disks using the Langlois recursive approach. This approach is successfully applied by Ullah et al. [40–42] for the creeping flow of slightly viscoelastic fluid and Maxwell fluid through a porous plane slit with uniform reabsorption, no-slip wall, and slip wall. Bhatti et al. [43] also used this technique for the second-order fluid flowing through a porous circular tube. Thus, we have employed the recursive approach to study the slow axisymmetric squeezed flow of slightly viscoelastic fluid film between two disks. The analytical expressions for the velocity field, pressure distribution, shear, and normal stress for squeezing flow have been obtained. The effects of the different physical parameters on the motion of fluid have been illustrated graphically. This study was considered to deliver answers to the following related research questions:(1)
What is the variation in velocity components, and pressure distribution of the squeeze flow due to the radial and axial direction?(2)
At various radial points, what is the impact on velocity profile, pressure distribution, and normal squeeze force due to the rise of slightly viscoelastic parameterβ?This paper is structured as follows. In Section2, the governing equations of motions with associated boundary conditions have been presented. In Section 3, the geometry of the problem and the model of the governing equations in cylindrical coordinates have been described. The construction of the velocity field for first, second, and third-order approximation problems has been explained in Section 4. In Section 5, the computations of the velocity field are given. The pressure distribution for squeezed flow has been obtained in section 6. The shear, normal stresses, and normal force are presented in Sections 7 and 8, respectively. The effects of different physical parameters have been investigated in section 9. Finally, in section 10, the concluding remarks are enlightened.


## 2. The Governing Equations of Motion

The basic equations governing the motion of an isotropic incompressible fluid [42] are given as(1)divV˜=0,(2)divτ˜+ρf˜=ρDV˜Dt,where V˜ denotes the velocity vector, ρ represents the constant fluid density, f˜ represents the body force per unit mass, and the material time derivative is denoted by D/Dt. The constitutive equation of third-order fluid for Cauchy stress by Truesdell and Noll [44] is given as(3)τ˜=−pI˜+μA˜1+α1A˜2+α2A˜12+β1A˜3+β2A˜1A˜2+A˜2A˜1+β3A˜12A˜1,where p represents the indeterminate part of the stress due to constraint of incompressibility, Ι is identity tensor, μ represents the coefficient of viscosity, the material constants are represented by α1,α2,β1,β2 and β3, A˜12=traceA˜1A˜1T and A˜1,A˜3 and A˜3 are the Rivlin-Ericksen tensors and defined as follows:(4)A˜1=gradV˜+gradV˜T,(5)A˜n+1=DDtA˜n+A˜ngradV˜+gradV˜TA˜nforn≥1.In addition, thermodynamics analysis imposes conditions [45] if all the motions of the fluid meet the Clausius‐Duhem inequality, and assuming the Helmholtz free energy is minimum while the fluid is at rest, then the material coefficients met the following restrictions:(6)μ≥0,α1≥0,α1+α2≤24μβ3,β1=β2=0,β3≥0.By making use of (6), the reduced form of (3) yields as follows:(7)τ˜=−pI˜+μA˜1+α1A˜2+α2A˜12+β3A˜12A˜1.The authors [45] presented (7) and studied the thermal effects in the variable viscosity of journal bearings. By considering the material constants equal to zero in (7), which becomes as follows:(8)τ˜=−pI˜+μA˜1+β3A˜12A˜1.Authors [37, 45] claimed that (8) represents the constitutive relation for the slightly viscoelastic fluid and can be considered the special class of differential type fluids or a special type of power-law model. Furthermore, in (8), the slightly viscoelastic parameter β3>0 and β3<0 represent the shear thickening and shear thinning behavior of the fluid, respectively. However, for β=0, (8) reduces to the Newtonian fluid.Substituting (8) in (2), yielding the result in vector form as(9)divV˜=0,(10)ρDV˜Dt=−gradp+μ∇2V˜+βA˜1gradA˜12+βA˜12∇2V˜+ρf˜,where ∇2 is the Laplacian operator and β3=β.


## 3. Problem Statement

The slow squeeze flow of an incompressible slightly viscoelastic fluid in axisymmetric form has been considered in Figure1. The lower disk is kept fixed while the above disk moves with constant velocity v(t) under the force F, which approaches the lower disk. H(t) represents the fluid film thickness, which decreases with time and fluid flows outside the disks radially. The velocity components for the axisymmetric flows are assumed as follows:(11)V˜=ur,θ,z,vr,θ,z,wr,θ,z=ur,z,0,wr,z;v=0,∂∂θ.=0,where ur,z and wr,z are the radial and axial velocity components.Figure 1
Geometry of squeeze flow [21].In order to write the components form of equations (9) and (10), substituting (11) in the first Rivlin tensor A˜1, we get(12)A˜1=gradV˜+gradV˜T=2∂u∂r0∂u∂z+∂w∂r02ur0∂u∂z+∂w∂r02∂w∂z.Using (12) in the definition of A˜12=traceA˜1A˜1T, we obtain(13)A˜12=traceA˜1A˜1T=4∂u∂r2+ur2+∂w∂z2+2∂w∂r+∂u∂z2=Msay.Putting (11)–(13) in terms of (10), we get the following expressions:(14)∇2V˜=∇2u−ur20∇2w,(15)A˜1gradA˜12=∂M∂z∂w∂r+∂u∂z+2∂M∂r∂u∂r0∂M∂r∂w∂r+∂u∂z+2∂M∂z∂w∂z,(16)A˜12∇2V˜=M∇2u−ur20M∇2w,where M=A˜12.(17)gradp=∂p∂r0∂p∂z.Neglecting the inertial term and body force and putting (11)–(17) in (9) and (10), we get the following nonlinear system of partial differential equations:(18)∂u∂r+ur+∂w∂z=0,(19)∂p∂r=μ+βM∇2u−1r2u+β∂M∂z∂w∂r+∂u∂z+2∂M∂r∂u∂r.(20)∂p∂z=μ+βM∇2w+β∂M∂r∂w∂r+∂u∂z+2∂M∂z∂w∂z.Similarly, by putting (12) and (13) into (8), four nonzero components form for the stress tensor are given as(21)τ˜rr=−p+2μ+βM∂u∂r;τ˜zr=μ+βM∂w∂r+∂u∂z,(22)τθθ=−p+2μ+βMur;τ˜zz=−p+2μ+βM∂w∂z,where M=4∂u/∂r2+u/r2+∂w/∂z2+2∂w/∂r+∂u/∂z2.We plan to solve the above set of nonlinear partial differential equations subject to the following nonhomogeneous boundary conditions:(23)w=−vt,u=0atz=Ht;u=0,w=0atz=0,where vt=−dH/dt.The first two conditions in (23) are due to the adherence of the top disk, and the location Ht of the top disk is unknown. In the other two conditions, there is no slip at the bottom disk. Velocity may be constant, or it may vary with time.


## 4. Construction of an Analytical Approximation of Solution by Langlois Recursive Approach

The coupled system of nonlinear partial differential equations (18)–(22) subject to nonhomogeneous boundary conditions (23) is not easy to be solved exactly. This nonlinear system is solved by using the recursive approach suggested by Langlois [39], which linearises the governing equations of motion given in equations (18)–(22). To obtain the approximate analytical solution of equations (18)–(22) subject to boundary conditions (23), the following equations for the velocity field, pressure, and stresses are sought in the form of perturbation of the state of the rest.(24)ur,z=∑i=13εiuir,z,(25)wr,z=∑i=13εiwir,z,(26)pr,z=constant+∑i=13εipir,z,(27)τ˜r,z=∑i=13εiτ˜ir,z,where ε is a small dimensionless number, these assumptions lead to the systems of linear partial differential equations corresponding to the boundary conditions for every set ui,wi,pi,τi,i=1,2,3, so that ur,z,wr,z,pr,z,τr,z as given by equations (24)–(27) provide the solution to equations (18)–(22). The first-order equations corresponding to u1,w1,p1,τ1 are describing the solution of the governing equation of the Newtonian fluid subject to nonhomogeneous boundary conditions. The second-order equations corresponding to u2,w2,p2,τ2 are similar except the nonhomogeneous part containing the terms of u1,w1 with homogeneous boundary conditions. The third-order equations corresponding to u3,w3,p3,τ3 are identical except the nonhomogeneous part containing the terms of u1,w1,u2,w2 with homogeneous boundary conditions; hence, third-order equations give the contribution of the slightly viscoelastic term.Substituting equations (24)–(27) into equations (19)–(21) and collecting the coefficients of equal powers of ε, we get the following first, second, and third-order boundary value problems, and the aim is to solve these problems Oε,Oε2, and Oε3.


### 4.1. First-Order Problem with Associated Nonhomogeneous Boundary Conditions

On equating the terms ofε from equations (18)–(22), the following linear system of partial differential equations subject to nonhomogeneous boundary conditions is obtained.(28a)∂u1∂r+u1r+∂w1∂z=0,(28b)∂p1∂r=μ∇2u1−1r2u1,(28c)∂p1∂z=μ∇2w1,(29)τ˜1rr=−p1+2μ∂u1∂r;τ˜1zr=μ∂w1∂r+∂u1∂z;τ˜1θθ=−p1+2μu1r;τ˜1zz=−p1+2μ∂w1∂z,subject to boundary conditions(30)w1=−Vt,u1=0atz=Ht;u1=0,w1=0atz=0.


### 4.2. Second-Order Problem with Associated Homogeneous Boundary Conditions

On equating the terms ofε2 from equations (18)–(22), the following linear system of partial differential equations subject to nonhomogeneous boundary conditions is obtained.(31)∂u2∂r+u2r+∂w2∂z=0;∂p2∂r=μ∇2u2−1r2u2;∂p2∂z=μ∇2w2,(32)τ˜2rr=−p2+2μ∂u2∂r;τ˜2zr=μ∂w2∂r+∂u2∂z;τ˜2θθ=−p2+2μu2r;τ˜2zz=−p2+2μ∂w2∂z,subject to homogeneous boundary conditions(33)w2=0,u2=0atz=Ht;u2=0,w2=0atz=0.


### 4.3. Third-Order Problem with Associated Homogeneous Boundary Conditions

On equating the terms ofε3 from equations (18)–(22), the following linear system of partial differential equations subject to nonhomogeneous boundary conditions is obtained.(34a)∂u3∂r+u3r+∂w3∂z=0,(34b)∂p3∂r=μ∇2u3−1r2u3+βM2∇2u1−1r2u1+β∂M2∂z∂w1∂r+∂u1∂z+2∂M2∂r∂u1∂r,(34c)∂p3∂z=μ∇2w3+βM2∇2w1+β∂M2∂r∂w1∂r+∂u1∂z+2∂M2∂z∂w1∂z,(35a)τ˜3rr=−p3+2μ∂u3∂r+2βM2∂u1∂r;τzr3=μ∂w3∂r+∂u3∂z+βM2∂w1∂r+∂u1∂z,(35b)τ˜3θθ=−p3+2μu3r+βM2u1r;τ˜3zz=−p3+2μ∂w3∂z++βM2∂w1∂z,where M2=4∂u1/∂r2+u1/r2+∂w1/∂z2+2∂w1/∂r+∂u1/∂z2,subject to homogeneous boundary conditions(36)w3=0,u3=0atz=Ht;u3=0,w3=0atz=0.


## 4.1. First-Order Problem with Associated Nonhomogeneous Boundary Conditions

On equating the terms ofε from equations (18)–(22), the following linear system of partial differential equations subject to nonhomogeneous boundary conditions is obtained.(28a)∂u1∂r+u1r+∂w1∂z=0,(28b)∂p1∂r=μ∇2u1−1r2u1,(28c)∂p1∂z=μ∇2w1,(29)τ˜1rr=−p1+2μ∂u1∂r;τ˜1zr=μ∂w1∂r+∂u1∂z;τ˜1θθ=−p1+2μu1r;τ˜1zz=−p1+2μ∂w1∂z,subject to boundary conditions(30)w1=−Vt,u1=0atz=Ht;u1=0,w1=0atz=0.


## 4.2. Second-Order Problem with Associated Homogeneous Boundary Conditions

On equating the terms ofε2 from equations (18)–(22), the following linear system of partial differential equations subject to nonhomogeneous boundary conditions is obtained.(31)∂u2∂r+u2r+∂w2∂z=0;∂p2∂r=μ∇2u2−1r2u2;∂p2∂z=μ∇2w2,(32)τ˜2rr=−p2+2μ∂u2∂r;τ˜2zr=μ∂w2∂r+∂u2∂z;τ˜2θθ=−p2+2μu2r;τ˜2zz=−p2+2μ∂w2∂z,subject to homogeneous boundary conditions(33)w2=0,u2=0atz=Ht;u2=0,w2=0atz=0.


## 4.3. Third-Order Problem with Associated Homogeneous Boundary Conditions

On equating the terms ofε3 from equations (18)–(22), the following linear system of partial differential equations subject to nonhomogeneous boundary conditions is obtained.(34a)∂u3∂r+u3r+∂w3∂z=0,(34b)∂p3∂r=μ∇2u3−1r2u3+βM2∇2u1−1r2u1+β∂M2∂z∂w1∂r+∂u1∂z+2∂M2∂r∂u1∂r,(34c)∂p3∂z=μ∇2w3+βM2∇2w1+β∂M2∂r∂w1∂r+∂u1∂z+2∂M2∂z∂w1∂z,(35a)τ˜3rr=−p3+2μ∂u3∂r+2βM2∂u1∂r;τzr3=μ∂w3∂r+∂u3∂z+βM2∂w1∂r+∂u1∂z,(35b)τ˜3θθ=−p3+2μu3r+βM2u1r;τ˜3zz=−p3+2μ∂w3∂z++βM2∂w1∂z,where M2=4∂u1/∂r2+u1/r2+∂w1/∂z2+2∂w1/∂r+∂u1/∂z2,subject to homogeneous boundary conditions(36)w3=0,u3=0atz=Ht;u3=0,w3=0atz=0.


## 5. Computation of Velocity Field

### 5.1. Velocity Field for First-Order Problem

In this subsection, we compute the velocity field of first-order approximation by rewriting the system of partial differential equations (28a)–(28c) with associated boundary conditions (30) in terms of stream function. The radial and axial velocity components of the axisymmetric flow can be expressed in the form of scalar stream function as follows:(37)u1=−1r∂ψ1∂z,w1=1r∂ψ1∂r.It should be noted that (28a) is identically satisfied by (37). By differentiating (28b) and (28c) partially with respect to z and r, respectively, and by eliminating the pressure, the compatibility equation of first-order approximation is obtained as(38)E4ψ1r,z=0,where E2=∂2/∂r2−1/r∂/∂r+∂2/∂z2 and E4∗=E2E2∗. Furthermore, by using (37) in (30), boundary conditions are obtained in terms of stream function as follows:(39)1r∂ψ1∂r=−V,∂ψ1∂z=0atz=Ht;∂ψ1∂r=0,∂ψ1∂z=0atz=0.In order to obtain the solution for the boundary value problem given in (38) and (39) by using the inverse method [46], assuming the following solution for stream function, a priori:(40)ψ1r,z=r2ϕ1z,where ϕ1z is an unknown function, which needs to be determined. Thus, using Equation (40) into (38) and (39) yields(41)d4ϕ1dz4=0,(42)ϕ1z=−V2,dϕ1dz=0atz=Ht;ϕ1z=0,dϕ1dz=0atz=0.The unknown functionϕ1z can be obtained from (41) by integrating and using associated conditions (42).(43)ϕ1z=VzH3−32zH2.The steam function, redial velocity, and the axial velocity can be obtained by mathematical simplification, using equation (43) into (40) and (37).(44)ψ1r,z=r2VzH3−32zH2.(45)u1r,z=3rVHzH−zH2.(46)w1r,z=−3VzH2−23zH3.It is noticed that the first-order velocity components (45) and (46) are in coherence with those obtained in [4] for the creeping squeeze flow of viscous fluid between two disks.


### 5.2. Velocity Field for Second-Order Problem

We determine the velocity field of second-order approximation from equations (31) and (32) subject to homogeneous boundary conditions (34a) by reducing the system of partial differential equations in terms of stream function. Defining the stream function for second-order approximation is as follows:(47)u2=−1r∂ψ2∂z,w2=1r∂ψ2∂r.It is noted that (31) is identically satisfied by using (47). Therefore, using (47) into (31) and eliminating the pressure, it reduces into the following form:(48)E4ψ2r,z=0.Furthermore, associated boundary conditions (33) are reduced in terms of stream function as follows:(49)1r∂ψ2∂r=0,∂ψ2∂z=0atz=Ht;∂ψ2∂r=0,∂ψ2∂z=0atz=0.The inverse solution of (48) corresponding to conditions (49) for any assumed form ψ2 is zero due to the homogeneous boundary conditions. Therefore, we get ψ2r,z=0.The radial and axial velocity components for second-order approximation are obtained by substituting stream function into (47). We get u2r,z=0 and w2r,z=0.


### 5.3. Velocity Field for Third-Order Problem

In this subsection, we compute the third-order approximations of the velocity field from equations (34a–34c) corresponding to homogeneous boundary conditions (36), and by using the first-order solution in (34a), we get(50)∂u3∂r+u3r+∂w3∂z=0.(51)∂p3∂r−μ∇2u3−1r2u3=−108V3r3H73−12zH+12zH2β−256V3rH534zH4−68zH3+41zH2−7zHβ.(52)∂p3∂z−μ∇2w3=432V3r2H62zH−6zH2+4zH3β−432V3H430zH5−75zH4+60zH3−15zH2β.Defining the stream function for the third-order problem is as follows:(53)u3=−1r∂ψ3∂z,w3=1r∂ψ3∂r.(50) is identically satisfied, and using (53) in equations (51) and (52) by eliminating the pressure, the following equation is obtained:(54)μ1rE4ψ3=216V3rβH67−90zH+228zH2−152zH3+216V3r3βH86−12zH.Associated boundary conditions in terms of stream function are written as follows:(55)1r∂ψ3∂r=0,∂ψ3∂z=0atz=Ht;∂ψ3∂r=0,∂ψ3∂z=0atz=0.Similarly, an inverse solution is sought for stream functionψ3, a priori as follows:(56)ψ3r,z=r2ϕ3z+r4ξ3z,where ϕ3z and ξ3z are unknowns and to be determined. On substituting (56) into (54), it takes the following form:(57)μrϕ3ivz+16ξ3ivz+r3ξ3ivz=216V3r3βH86−12zH+216V3rβH67−90zH+228zH2−152zH3.The following equations (58) and (59) are obtained by comparing the coefficients of r and r3 in (57).(58)μξ3ivz=1296V3H8β1−2zH.(59)μϕ3ivz+16ξ3″z=216V3βH67−90zH+228zH2−152zH3.Using the assumption (56), the corresponding boundary conditions (55) reduce in (60).(60)ξ3H=0,ξ3′H=0,ξ30=0,ξ3′0=0;ϕ3H=0,ϕ3′H=0,ϕ30=0,ϕ3′0=0.The solutions of an ordinary differential (58) and (59) corresponding to boundary conditions (60) are obtained as(61)ξ3z=54V3β5μH4zH2−4zH3+5zH4−2zH5.(62)ϕ3z=27V3β175μH6−37zH2+48zH3+315zH4−826zH5+700zH6−200zH7.The expressions for stream function and velocity components are obtained by using equations (61) and (62) into equations (53) and (56) as follows:(63)ψ3r,z=54V3r4β5μH4zH2−4zH3+5zH4−2zH5+27V3r2β175μH6−37zH2+48zH3+315zH4−826zH5+700zH6−200zH7.(64)u3r,z=54V3β175μH337zH−72zH2−630zH3+2065zH4−2100zH5+700zH6r+54V3β175μH5−70zH+420zH2−700zH3+350zH4r3.(65)w3r,z=54V3β175μH2−37zH2+48zH3+315zH4−826zH5+700zH6−200zH7+54V3β175μH4140zH2−560zH3+700zH4−280zH5r2.It is observed that the third-order approximate solution of the velocity profile and stream function contains the terms of the slightly viscoelastic parameterβ, which is the key feature of the present study.


### 5.4. Stream Function and Velocity Field Correct to Third-Order Approximations

The expressions for stream function correct to third-order approximation obtained by adding equations (44) and (63) as follows:(66)ψr,z=r2VzH3−32zH2+54V3r4β5μH4zH2−4zH3+5zH4−2zH5+27V3r2β175μH6−37zH2+48zH3+315zH4−826zH5+700zH6−200zH7+oε3.In order to obtain the velocity components, correct to third-order by substituting equations (45) and (64) into (24), we get(67)ur,z=3rVHzH−zH2+54V3β175μH337zH−72zH2−630zH3+2065zH4−2100zH5+700zH6r+54V3β175μH5−70zH+420zH2−700zH3+350zH4r3+oε3.Similarly, by substituting equations (46) and (65) into (25), we get(68)wr,z=−3VtzH2−23zH3+54V3β175μH2−37zH2+48zH3+315zH4−826zH5+700zH6−200zH7+54V3β175μH4140zH2−560zH3+700zH4−280zH5r2+oε3a.Results achieved from this proposed approach agree with the results presented in the literature [47, 48] for viscous fluids when β = 0.


## 5.1. Velocity Field for First-Order Problem

In this subsection, we compute the velocity field of first-order approximation by rewriting the system of partial differential equations (28a)–(28c) with associated boundary conditions (30) in terms of stream function. The radial and axial velocity components of the axisymmetric flow can be expressed in the form of scalar stream function as follows:(37)u1=−1r∂ψ1∂z,w1=1r∂ψ1∂r.It should be noted that (28a) is identically satisfied by (37). By differentiating (28b) and (28c) partially with respect to z and r, respectively, and by eliminating the pressure, the compatibility equation of first-order approximation is obtained as(38)E4ψ1r,z=0,where E2=∂2/∂r2−1/r∂/∂r+∂2/∂z2 and E4∗=E2E2∗. Furthermore, by using (37) in (30), boundary conditions are obtained in terms of stream function as follows:(39)1r∂ψ1∂r=−V,∂ψ1∂z=0atz=Ht;∂ψ1∂r=0,∂ψ1∂z=0atz=0.In order to obtain the solution for the boundary value problem given in (38) and (39) by using the inverse method [46], assuming the following solution for stream function, a priori:(40)ψ1r,z=r2ϕ1z,where ϕ1z is an unknown function, which needs to be determined. Thus, using Equation (40) into (38) and (39) yields(41)d4ϕ1dz4=0,(42)ϕ1z=−V2,dϕ1dz=0atz=Ht;ϕ1z=0,dϕ1dz=0atz=0.The unknown functionϕ1z can be obtained from (41) by integrating and using associated conditions (42).(43)ϕ1z=VzH3−32zH2.The steam function, redial velocity, and the axial velocity can be obtained by mathematical simplification, using equation (43) into (40) and (37).(44)ψ1r,z=r2VzH3−32zH2.(45)u1r,z=3rVHzH−zH2.(46)w1r,z=−3VzH2−23zH3.It is noticed that the first-order velocity components (45) and (46) are in coherence with those obtained in [4] for the creeping squeeze flow of viscous fluid between two disks.


## 5.2. Velocity Field for Second-Order Problem

We determine the velocity field of second-order approximation from equations (31) and (32) subject to homogeneous boundary conditions (34a) by reducing the system of partial differential equations in terms of stream function. Defining the stream function for second-order approximation is as follows:(47)u2=−1r∂ψ2∂z,w2=1r∂ψ2∂r.It is noted that (31) is identically satisfied by using (47). Therefore, using (47) into (31) and eliminating the pressure, it reduces into the following form:(48)E4ψ2r,z=0.Furthermore, associated boundary conditions (33) are reduced in terms of stream function as follows:(49)1r∂ψ2∂r=0,∂ψ2∂z=0atz=Ht;∂ψ2∂r=0,∂ψ2∂z=0atz=0.The inverse solution of (48) corresponding to conditions (49) for any assumed form ψ2 is zero due to the homogeneous boundary conditions. Therefore, we get ψ2r,z=0.The radial and axial velocity components for second-order approximation are obtained by substituting stream function into (47). We get u2r,z=0 and w2r,z=0.


## 5.3. Velocity Field for Third-Order Problem

In this subsection, we compute the third-order approximations of the velocity field from equations (34a–34c) corresponding to homogeneous boundary conditions (36), and by using the first-order solution in (34a), we get(50)∂u3∂r+u3r+∂w3∂z=0.(51)∂p3∂r−μ∇2u3−1r2u3=−108V3r3H73−12zH+12zH2β−256V3rH534zH4−68zH3+41zH2−7zHβ.(52)∂p3∂z−μ∇2w3=432V3r2H62zH−6zH2+4zH3β−432V3H430zH5−75zH4+60zH3−15zH2β.Defining the stream function for the third-order problem is as follows:(53)u3=−1r∂ψ3∂z,w3=1r∂ψ3∂r.(50) is identically satisfied, and using (53) in equations (51) and (52) by eliminating the pressure, the following equation is obtained:(54)μ1rE4ψ3=216V3rβH67−90zH+228zH2−152zH3+216V3r3βH86−12zH.Associated boundary conditions in terms of stream function are written as follows:(55)1r∂ψ3∂r=0,∂ψ3∂z=0atz=Ht;∂ψ3∂r=0,∂ψ3∂z=0atz=0.Similarly, an inverse solution is sought for stream functionψ3, a priori as follows:(56)ψ3r,z=r2ϕ3z+r4ξ3z,where ϕ3z and ξ3z are unknowns and to be determined. On substituting (56) into (54), it takes the following form:(57)μrϕ3ivz+16ξ3ivz+r3ξ3ivz=216V3r3βH86−12zH+216V3rβH67−90zH+228zH2−152zH3.The following equations (58) and (59) are obtained by comparing the coefficients of r and r3 in (57).(58)μξ3ivz=1296V3H8β1−2zH.(59)μϕ3ivz+16ξ3″z=216V3βH67−90zH+228zH2−152zH3.Using the assumption (56), the corresponding boundary conditions (55) reduce in (60).(60)ξ3H=0,ξ3′H=0,ξ30=0,ξ3′0=0;ϕ3H=0,ϕ3′H=0,ϕ30=0,ϕ3′0=0.The solutions of an ordinary differential (58) and (59) corresponding to boundary conditions (60) are obtained as(61)ξ3z=54V3β5μH4zH2−4zH3+5zH4−2zH5.(62)ϕ3z=27V3β175μH6−37zH2+48zH3+315zH4−826zH5+700zH6−200zH7.The expressions for stream function and velocity components are obtained by using equations (61) and (62) into equations (53) and (56) as follows:(63)ψ3r,z=54V3r4β5μH4zH2−4zH3+5zH4−2zH5+27V3r2β175μH6−37zH2+48zH3+315zH4−826zH5+700zH6−200zH7.(64)u3r,z=54V3β175μH337zH−72zH2−630zH3+2065zH4−2100zH5+700zH6r+54V3β175μH5−70zH+420zH2−700zH3+350zH4r3.(65)w3r,z=54V3β175μH2−37zH2+48zH3+315zH4−826zH5+700zH6−200zH7+54V3β175μH4140zH2−560zH3+700zH4−280zH5r2.It is observed that the third-order approximate solution of the velocity profile and stream function contains the terms of the slightly viscoelastic parameterβ, which is the key feature of the present study.


## 5.4. Stream Function and Velocity Field Correct to Third-Order Approximations

The expressions for stream function correct to third-order approximation obtained by adding equations (44) and (63) as follows:(66)ψr,z=r2VzH3−32zH2+54V3r4β5μH4zH2−4zH3+5zH4−2zH5+27V3r2β175μH6−37zH2+48zH3+315zH4−826zH5+700zH6−200zH7+oε3.In order to obtain the velocity components, correct to third-order by substituting equations (45) and (64) into (24), we get(67)ur,z=3rVHzH−zH2+54V3β175μH337zH−72zH2−630zH3+2065zH4−2100zH5+700zH6r+54V3β175μH5−70zH+420zH2−700zH3+350zH4r3+oε3.Similarly, by substituting equations (46) and (65) into (25), we get(68)wr,z=−3VtzH2−23zH3+54V3β175μH2−37zH2+48zH3+315zH4−826zH5+700zH6−200zH7+54V3β175μH4140zH2−560zH3+700zH4−280zH5r2+oε3a.Results achieved from this proposed approach agree with the results presented in the literature [47, 48] for viscous fluids when β = 0.


## 6. Computation of Pressure Distribution

### 6.1. First-Order Approximation of Pressure Distribution

Substituting the first-order approximation of velocity profile (45) and (46) into equations (28b) and (28c), the following equations of pressure distribution are obtained as(69)∂p1∂r=−6μVH3r.(70)∂p1∂z=6μVH22zH−1.Integrating (69) with respect to r, we obtain(71)p1r,z=−3μVH3r2+gzwhere gz is an arbitrary function of z. Differentiating (71) with respect to z and comparing with (70) yield the following(72)g′z=6μVH22zH−1,gz=6μVHzH2−zH+l,where l is the constant of integration. Using equation (72) into (71), the first-order pressure distribution gets as follows:(73)p1r,z=3μVH2zH2−2zH−r2H2+l.


### 6.2. Second-Order Approximation of Pressure Distribution

Subsequently, substituting the second-order approximation of velocity profile into (31), the following equations of pressure distribution are obtained as(74)∂p2∂r=0;∂p2∂z=0.After integrating (74) and simplification, the second-order pressure distribution is obtained p2r,z=m where m is the constant of integration.


### 6.3. Third-Order Approximation of Pressure Distribution

Similarly, substituting the third-order approximation of velocity profile (64) and (65) into equations (34b) and (34c), the following equations of pressure distribution are obtained as(75)∂p3∂r=54V3β175H7−210r3+54V3β175H5r−144+560zH−560zH2.(76)∂p3∂z=54V3β175H4−74+288zH−16660zH2+65240zH3−81200zH4+32480zH5+54V3β175H6r2280−560zH.The solution of the third-order approximation of pressure distribution is obtained from equations (75) and (76) by using the procedure of the first-order pressure distribution as follows:(77)p3r,z=−2835V3β175H7r4+54V3β175H5r2−72+280zH−280zH2+54V3β175H3−74zH+144zH2−166603zH3+16310zH4−16240zH5+162403zH6+n.


### 6.4. Pressure Distribution Correct to Third-Order Approximation

For obtaining the pressure distribution, (78) corrects to third-order approximation, adding equations (73) and (77).(78)pr,z=3μVH2zH2−2zH−r2H2−2835V3β175H7r4+54V3β175H5r2−72+280zH−280zH2+54V3β175H3−74zH+144zH2−166603zH3+16310zH4−16240zH5+162403zH6+K.where K=lε+mε2+nε3. For the evaluation of the constant K, there is required an extra boundary condition at the outer edge of the disk.In this research investigation, it is observed that the edge of the disk of fluid that resides is subjected to atmospheric pressure; therefore, at the free surface, the balance of normal stressτzz is equivalent to atmospheric pressure. There is no loss of generality in taking τzz=0[9].(79)τzz=−p+2μ+βM∂w∂z=0atz=Handr=R.From (68), at z=H, we get ∂w/∂z=0; hence, (79) reduces to pR,H=0. Using this condition in Equation (78) to obtain the constant K=3μVR2/H3+2835V3R4/175H7β+3888V3R2/175H5β.Putting the value ofK in (78), we have pressure distribution correct to third-order approximation:(80)pr,z=3μVH2zHzH−1+R2−r2H2+2835V3β175H7R4−r4+3888V3β175H5R2−r2+54V3βr2175H5280zH−280zH2+54V3β175H3−74zH+144zH2−166603zH3+16310zH4−16240zH5+162403zH6+oε3.


## 6.1. First-Order Approximation of Pressure Distribution

Substituting the first-order approximation of velocity profile (45) and (46) into equations (28b) and (28c), the following equations of pressure distribution are obtained as(69)∂p1∂r=−6μVH3r.(70)∂p1∂z=6μVH22zH−1.Integrating (69) with respect to r, we obtain(71)p1r,z=−3μVH3r2+gzwhere gz is an arbitrary function of z. Differentiating (71) with respect to z and comparing with (70) yield the following(72)g′z=6μVH22zH−1,gz=6μVHzH2−zH+l,where l is the constant of integration. Using equation (72) into (71), the first-order pressure distribution gets as follows:(73)p1r,z=3μVH2zH2−2zH−r2H2+l.


## 6.2. Second-Order Approximation of Pressure Distribution

Subsequently, substituting the second-order approximation of velocity profile into (31), the following equations of pressure distribution are obtained as(74)∂p2∂r=0;∂p2∂z=0.After integrating (74) and simplification, the second-order pressure distribution is obtained p2r,z=m where m is the constant of integration.


## 6.3. Third-Order Approximation of Pressure Distribution

Similarly, substituting the third-order approximation of velocity profile (64) and (65) into equations (34b) and (34c), the following equations of pressure distribution are obtained as(75)∂p3∂r=54V3β175H7−210r3+54V3β175H5r−144+560zH−560zH2.(76)∂p3∂z=54V3β175H4−74+288zH−16660zH2+65240zH3−81200zH4+32480zH5+54V3β175H6r2280−560zH.The solution of the third-order approximation of pressure distribution is obtained from equations (75) and (76) by using the procedure of the first-order pressure distribution as follows:(77)p3r,z=−2835V3β175H7r4+54V3β175H5r2−72+280zH−280zH2+54V3β175H3−74zH+144zH2−166603zH3+16310zH4−16240zH5+162403zH6+n.


## 6.4. Pressure Distribution Correct to Third-Order Approximation

For obtaining the pressure distribution, (78) corrects to third-order approximation, adding equations (73) and (77).(78)pr,z=3μVH2zH2−2zH−r2H2−2835V3β175H7r4+54V3β175H5r2−72+280zH−280zH2+54V3β175H3−74zH+144zH2−166603zH3+16310zH4−16240zH5+162403zH6+K.where K=lε+mε2+nε3. For the evaluation of the constant K, there is required an extra boundary condition at the outer edge of the disk.In this research investigation, it is observed that the edge of the disk of fluid that resides is subjected to atmospheric pressure; therefore, at the free surface, the balance of normal stressτzz is equivalent to atmospheric pressure. There is no loss of generality in taking τzz=0[9].(79)τzz=−p+2μ+βM∂w∂z=0atz=Handr=R.From (68), at z=H, we get ∂w/∂z=0; hence, (79) reduces to pR,H=0. Using this condition in Equation (78) to obtain the constant K=3μVR2/H3+2835V3R4/175H7β+3888V3R2/175H5β.Putting the value ofK in (78), we have pressure distribution correct to third-order approximation:(80)pr,z=3μVH2zHzH−1+R2−r2H2+2835V3β175H7R4−r4+3888V3β175H5R2−r2+54V3βr2175H5280zH−280zH2+54V3β175H3−74zH+144zH2−166603zH3+16310zH4−16240zH5+162403zH6+oε3.


## 7. Computation of Shear and Normal Stresses

### 7.1. First-Order Shear and Normal Stresses

The expressions for first-order approximations of shear and normal stresses are obtained by substituting first-order velocity approximations (45) and (46) and pressure distribution (73) into (29) as follows:(81)τ˜rr1=3VμH4zH−zH2+r2−R2H2;τ˜zr1=3VμH21−2zHr2,τ˜θθ1=3VμH4zH−zH2+r2−R2H2;τ˜zz1=3VμH2zH2−zH+r2−R2H2.


### 7.2. Second-Order Shear and Normal Stresses

Similarly, the expressions for second-order shear and normal stresses are obtained by substituting second-order velocity approximation and pressure distribution into (32) as follows:(82)τ˜rr2=τ˜zr2=τ˜θθ2=τ˜zz2=0.


### 7.3. Third-Order Shear and Normal Stresses

Similarly, the expressions for third-order shear and normal stresses are obtained by substituting third-order velocity approximation (67) and (68) and pressure distribution (77) into Equation (35a) as follows:(83a)τ˜3rr=9V3β175H54200zH4−8400zH3+4620zH2−140zHr2+432r2−R2−1680zH−zH2.(83b)τzr3=54V3β175H437−144zH+490zH2−1260zH3+1400zH4−560zH5r+9V3β175H3888zH−1728zH2+50960zH3−148680zH4+147840zH5−49280zH6−315r4−R4H4.(83c)τ˜zz3=−9V3β175H5−432r2−R2+1680zH−zH2+1260−zH+5.75zH2−13.33zH3+6.66zH4r2−9V3β175H3444zH−864zH2−23240zH3+71820zH4+72240zH5−24080zH6+315r4−R4H4.(83d)τ˜θθ3=9V3β175H5432r2−R2−1680zH−zH2+210zH−zH2r2+9V3β175H3888zH−1728zH2+38360zH3−110880zH4+110040zH5−36680zH6−315r4−R4H4.


## 7.1. First-Order Shear and Normal Stresses

The expressions for first-order approximations of shear and normal stresses are obtained by substituting first-order velocity approximations (45) and (46) and pressure distribution (73) into (29) as follows:(81)τ˜rr1=3VμH4zH−zH2+r2−R2H2;τ˜zr1=3VμH21−2zHr2,τ˜θθ1=3VμH4zH−zH2+r2−R2H2;τ˜zz1=3VμH2zH2−zH+r2−R2H2.


## 7.2. Second-Order Shear and Normal Stresses

Similarly, the expressions for second-order shear and normal stresses are obtained by substituting second-order velocity approximation and pressure distribution into (32) as follows:(82)τ˜rr2=τ˜zr2=τ˜θθ2=τ˜zz2=0.


## 7.3. Third-Order Shear and Normal Stresses

Similarly, the expressions for third-order shear and normal stresses are obtained by substituting third-order velocity approximation (67) and (68) and pressure distribution (77) into Equation (35a) as follows:(83a)τ˜3rr=9V3β175H54200zH4−8400zH3+4620zH2−140zHr2+432r2−R2−1680zH−zH2.(83b)τzr3=54V3β175H437−144zH+490zH2−1260zH3+1400zH4−560zH5r+9V3β175H3888zH−1728zH2+50960zH3−148680zH4+147840zH5−49280zH6−315r4−R4H4.(83c)τ˜zz3=−9V3β175H5−432r2−R2+1680zH−zH2+1260−zH+5.75zH2−13.33zH3+6.66zH4r2−9V3β175H3444zH−864zH2−23240zH3+71820zH4+72240zH5−24080zH6+315r4−R4H4.(83d)τ˜θθ3=9V3β175H5432r2−R2−1680zH−zH2+210zH−zH2r2+9V3β175H3888zH−1728zH2+38360zH3−110880zH4+110040zH5−36680zH6−315r4−R4H4.


## 8. Normal Force on the Upper Disk

The total force applied on the upper disk calculated by integral of the negative ofτzz that is exerted by the slightly viscoelastic fluid film in the positive z-direction at z=H is estimated as(84)F=∫0R2πrpr,Hdr.Substitutingz=H in (80) to obtain pr,H and using in (84), then obtained the following(85)F=3μπVR42H3+2πV3R4175H5β945RH2+972.


## 9. Analysis of Results

The convergence of the solution using the Langlois recursive approach is validated by determining the absolute residual error for different fluid parameters as given in Table1. It can be observed that the absolute residual error of the solution for third-order approximation is approaching to zero that confirmed the convergence of the solution.Table 1
Absolute residual error for Langlois recursive approach solution at various slightly viscoelastic parameters.
Zβ=0.2β=0.4β=0.60.1.01252∗10−84.05095∗10−89.11656∗10−80.16.56195∗10−92.62525∗10−85.90788∗10−80.22.24776∗10−98.9907∗10−92.02283∗100.31.97737∗10−97.91431∗10−91.7818∗10−80.45.54702∗10−92.21957∗10−84.99573∗10−80.58.11939∗10−93.24857∗10−87.31114∗10−80.69.53336∗10−93.81402∗10−88.58308∗10−80.79.76923∗10−93.90807∗10−88.79403∗10−80.88.91384∗10−93.56555∗10−88.0225∗10−80.97.13011∗10−92.85168∗10−86.41547∗10−81.04.63052∗10−91.85154∗10−84.16447∗10−8


## 10. Discussion of Results

In this research study, an approximate analytical solution of the slow squeeze flow of slightly viscoelastic fluid films between two disks is obtained successfully using the recursive approach of Langlois. The following dimensionless variables are introduced to examine the effect of the slightly viscoelastic parameterβ and radial distance r on flow variables such as velocity profile, pressure distribution, and normal force and depicted in Figures 2–7, drawn using the mathematics-based Maple software.(86)r∗=rR,z∗=zH,u∗=HVRu,w∗=wV,p∗=H3μVR2p,H∗=HR,β∗=V2βμHR.Figure 2
Radial velocity at different radial points (r) when aspect ratioH/R=0.4 and β=0.5.Figure 3
Radial velocity at different values ofβ when aspect ratio H/R=0.4 and r=0.5.Figure 4
Axial velocity at different radial points (r) when aspect ratioH/R=0.4 and β=0.5.Figure 5
Axial velocity at different values ofβ when aspect ratio H/R=0.4 and r=0.5.Figure 6
Pressure profiles at different values ofβ when z/H=0.7 and H/R=0.4.Figure 7
Normal force with constant squeezing velocity for different values ofβ.The impact of the radial velocity component at various radial points keeping fixed valuesβ=0.5 and H/R=0.4 is shown in Figure 2. It is observed that the radial velocity increases as the radius increases and velocity is maximum at the axial distance z^=0.5. The effect of the slightly viscoelastic parameter β on the radial velocity with fixed values r=0.5 and H/R=0.4 illustrated in Figure 3 and shows the behavior of the radial velocity is parabolic in the case of viscous fluid [4, 48, 49]. The radial velocity increases and occurs backward flow at the edges with the rise in β. Thus, physically slightly viscoelastic parameter β is inversely proportional to the fluid viscosity. The rise of this parameter β demotes the viscosity, and as a result, the fluid velocity intensified. This profile proved the behavior of the shear thickening fluid. Furthermore, behavior of the radial velocity in the case of viscous fluid, when β  = 0, is taken as a special case which is also illustrated in Figure 3.Figure4 represents the axial velocity profile at different radial values keeping β=0.5 and H=0.4. Axial velocity w rises in contrast to an increase in the value of r. Figure 5 depicts the impact of the slightly viscoelastic parameter β on the axial velocity w. It is noticed that the negative magnitude of axial velocity indicates that the direction of flow is downward; also, reverse flow is occurring in axial velocity due to the rise of β, and it is gradually increased with the maximum magnitude at z=0.75. In Figures 1–4, the shear-thickening behavior of the velocity profile appears at the parameter β.Figure6 demonstrates the effect of β on pressure distribution at H/R=0.4 and z/H=0.7. It shows that the pressure distribution is directly proportional to the thickness of the fluid β and reduces radially approaching to zero at the edge. Physically when the viscoelasticity of the material rises, then pressure distribution increases to squeeze the fluid film. In Figure 7, force swiftly increases at the significant values of the aspect ratio of the film thickness and high values of β.The results obtained from this proposed approach forβ  = 0 agree with the results presented in the literature [47,48]. Furthermore, it also satisfies the physical behavior of the shear thickening fluid.


## 11. Conclusion and Future Recommendations

This research study is focused on the analytical solution of the sluggish squeeze flow of the slightly viscoelastic fluid film between two circular disks where the upper disk is moving with constant velocity and the lower disk is stationary. The recursive approach of Langlois has been applied successfully for the analytical solution of the governing equations of the flow in the axisymmetric form. In order to examine the physical behavior of the flow, analysis is taken up to third-order linearization by this approach. The physical parameters of the flow such as velocity, pressure distribution, normal squeeze force, and stresses are determined as a function of aspect ratio H/R, radial distancer, and viscoelastic parameter β, which are key flow variables of the study. The obtained results by this approach are in good agreement with the results of the squeezing flow of viscous fluid by vanishing slightly viscoelastic parameter which also confirms the physical behavior of shear thickening fluid. The following key points are drawn from the theoretical investigations.(1)
The axial and radial velocities increase with the increment in the radial distance and observed that the maximum velocity occurs at thez = 0.5(2)
The axial and radial velocities increase at the higher values of the slightly viscoelastic parameterβ, and backward flow occurs at the boundaries of the channel due to the slightly viscoelastic term(3)
The pressure distribution elevates with the rise of the slightly viscoelastic parameterβ and crushes at the edges of the channel(4)
The presence of a slightly viscoelastic parameter boosted the normal squeeze force on the upper diskThis study attempted to get an analytical solution for the slow squeeze flow of the slightly viscoelastic fluid film between two circular disks at no-slip conditions by the Langlois recursive approach. As a future recommendation, it is proposed that this method can also be used to get an analytical solution of the nonlinear partial differential equations arising from the squeeze flow of differential type fluids with slip conditions and inertia effects.



---

*Source: 4043909-2022-06-11.xml*