Acta Polytechnica https://doi.org/10.14311/AP.2023.63.0439 Acta Polytechnica 63(6):439–450, 2023 © 2023 The Author(s). Licensed under a CC-BY 4.0 licence Published by the Czech Technical University in Prague SEMI-ANALYTICAL APPROACH-BASED STUDIES OF THE SQUEEZE FILM LUBRICATION BETWEEN ROUGH POROUS ANNULAR DISCS: RABINOWITSCH FLUID MODEL Amit Kumar Rahula,∗, Manoj Kumar Singha, Ravi Tiwarib, Sourabh Paulb, Pentyala Srinivasa Raoc, Rohahn Biswasd a Vellore Institute of Technology, School of Advanced Sciences (SAS), Division of Mathematics, Chennai 600127, Tamil Nadu, India b Vellore Institute of Technology, School of Electronics Engineering (SENSE), Chennai 600127, Tamil Nadu, India c Indian Institute of Technology, Department of Mathematics & Computing, Dhanbad, 826004, India d Vellore Institute of Technology, School of Electrical Engineering (SELECT), Chennai 600127, Tamil Nadu, India ∗ corresponding author: akrahulism@gmail.com Abstract. In recent years, there has been much interest in the effects of porosity and surface roughness (SR) or geometric irregularities between two moving plates under hydrodynamic lubrication. Porous bearings are used extensively in wide range of equipment, including computers, office equipment, home appliances, electric motors, and vehicles. In light of the importance of the aforementioned applications, we explored how SR and porous materials affect annular discs under the condition of a squeeze film. A five-point Gauss quadrature integral formula has been used to examine the characteristics of annular discs and a small perturbation method has been used to discretise the governing Rabinowitsch fluid flow (RFF) equations. The impact of nonlinear parameters on the behaviour of porosity and SR have been visualised in terms of film pressure (FP), load carrying capacity (LCC), and squeeze response time (SRT) of annular discs. Under the conditions of pseudoplastic and dilatant fluids, the effects of SR and porous materials between annular discs have been estimated in the form of the film pressure, LCC, and SRT and are presented in this manuscript as tables and graphs. According to the findings, the performance of an annular disc is significantly affected by porous material and radial roughness patterns. In addition, when RFF is carried through a rough surface and porous media, the performance is found to improve for dilatant fluids but suffer for pseudoplastic fluids. Keywords: Annular discs, squeeze film, Rabinowitsch fluid model, surface roughness, porous wall. 1. Introduction Squeeze film (SF) is used in a variety of applications, such as turbomachinery, disc clutches, and viscous- lock systems. Recently, there has been a significant increase in interest in this area. The performance of annular plates under SF has been studied under New- tonian fluids (NF) and non-Newtonian fluids (NNF) by Allen and McKillop [1], and Naduvinamani et al. [2]. It has been discovered that squeeze film (SF) increases the efficiency and reliability of bearings while also pro- longing their lifespan. The NF is a broad term for the linear relationship between shear strain and the rate of change. Experimental research has shown that novel fluids can be constructed from an NF to include some additional components that behave as an NNF. They consist of high-molecular-weight polymers, viscosity index improvers, polyisobutylene, etc. (Spike [3]). The non-Newtonian fluids (NNF) have a non-linear relationship with both the rate of shear strain and the shear stress. Different NNF models, such as the power- law fluid model, Rabinowitsch fluid model (RFM), and couple-stress fluid model, have been tested by a number of tribologists to investigate the performance of bearings. Lin and Hung [4] have investigated the performance of circular plates lubricated with a couple stress fluid model and a power-law fluid model by Wang et al. [5]. Several researchers have studied the demonstration of the RFM for different types of bearings. Wada and Hayashi [6, 7] were the first to study theoretically and experimentally the mechanism of journal bearing in the presence of pseudoplastic fluids. Furthermore, Lin [8] and Siddangouda et al. [9] investigated the performance of parallel annular discs and static characteristics of an inclined plane slider bearing under the condition of RFM. These studies suggest that under RFM the bearing performance is influenced by the characteristics of the dilatant fluids, whereas the pseudo-plastic behaviour results in a lower concentration of bearing stability compared to NF. For five decades, permeable materials have been widely used in industry to improve bearing perfor- mance (Bujurke et al. [10], Morgan and Cameron [11]). With the use of Darcy’s model, Morgan and Cameron [10] were the first researchers who inves- 439 https://doi.org/10.14311/AP.2023.63.0439 https://creativecommons.org/licenses/by/4.0/ https://www.cvut.cz/en A. K. Rahul, M. K. Singh, R. Tiwari et al. Acta Polytechnica tigated the porous effect on bearing surfaces. Patel et al. [12] have considered double-layer porous surface with a curved squeeze film. Bhat and Deheri [13] con- ducted a study on squeeze film behaviour in porous annular discs that are lubricated with magnetic fluid. Recently, Walicka et al. [14], Rao and Rahul [15– 17], and Rahul et al. [18] have applied the RFM on a curvilinear and SF bearing to study the characteristics of the dilatant fluid between the porous media. These researchers were concerned with how the porous wall affected the performance of the bearings. It was con- cluded that the Rabinowitsch fluid property, which exists in porous media, has a substantial impact on the bearings’ characteristics and lengthens their lifespan. The goal of many researchers is to show how SR affects the thin film lubrication bearing. The bear- ing performance, including pressure, load, friction, drags, and other factors, is influenced by this con- dition’s surface roughness. Various methodologies have been introduced to examine the impact of SR on the performance of bearings. Christensen [19, 20] first established the stochastic theory of thin film lu- brication of rough bearing surfaces. The combined effects of SR and porosity with MHD between annular discs have been taken into consideration by Baksh and Naganagowda [21]. By analysing the permeabil- ity influence on the rotor linear stability, D’Agostino et al. [22] examined the unstable oil film forces in porous bearings. Vashi et al. [23], have considered the Neuringer-Rosenzweig model to investigate the per- formance of circular stepped plates in the existence of couple stress, porosity, and SR Siddangouda et al. [24] investigated the combined impact of SR and viscosity change brought on by additives on long journal bear- ing. According to Munshi et al. [25], a rough porous sine film slider bearing with ferrofluid lubrication was subjected to a sinusoidal magnetic field. Shimpi and Deheri [26] investigated the deforma- tion effect of a magnetic fluid-based squeeze film in rough rotating curved porous annular plates. Patel et al. [27] investigated the lubrication of a rough porous hyperbolic slider bearing with slip velocity using fer- rofluid. Patel and Deheri [28] studied the joint effect of slip velocity and roughness on the ferrofluid lubri- cation of a curved rough annular squeeze film using the Jenkins model. Vashi et al. [29] studied the longi- tudinally rough porous circular stepping plates based on the Neuringer-Roseinweig model in the presence of couple stress. The influence of a porous structure, slip velocity, and Rosensweig’s viscosity on the ferrofluid- based squeeze film on porous curved annular plates was explored by Patel et al. [30]. Pressure genera- tions in a rough conical bearing using non-Newtonian Rabinowitsch fluid with variable viscosity have been investigated by Rao and Rahul [31]. A hydrostatic con- ical bearing that is externally pressured and lubricated with RFM has been studied by Walicka et al. [32] for how wall porosity and surface roughness affect steady performance. Recently, Rahul and Rao [33, 34] inves- tigated the behaviour of annular discs and circular stepping plates in the presence of a viscosity variable of Rabinowitsch fluid with a rough and porous surface. They have discussed that a one-dimensional azimuthal (Radial) roughness pattern on the rough porous circular plate increases (decreases) the load-carrying capacity and the squeeze film time as compared to the corre- sponding smooth case. Moreover, in circular stepped plates with the presence of the porous wall, the load capacity, and squeeze time decreases compared to non-porous case. 1.1. Research gap The above researchers concluded that transverse roughness patterns on the bearing surface increase FP and LCC. Despite the much anticipated future of squeeze film bearings, not enough research has been done. The majority of the work has been per- formed without taking into account the impacts of porous discs, the Rabinowsitch fluid model, and sur- face roughness, according to the literature review. The following issues weren’t covered in the papers that were previously published. a) Analysis of the characteristics of a squeeze-film of annular discs using the RF model to determine the effectiveness of the porosity and surface roughness. b) Investigation of the effect ω, H∗ 0 , R∗ 0, c∗ on the trend P ∗ for a constant Ω. Determining the in- fluence of ω, H∗ 0 , R∗ 0, c∗ on the trend W ∗ for a constant h∗ and Ω. c) The impact of ω, H∗ 0 , R∗ 0, c∗ on t∗ for a constant h∗ and Ω. Evaluation of the disc’s performance in various values of ω and Ω. We highlighted the below interesting points, which have not been addressed in the study of Lin [8]. a) It has been considered to apply the idea of a stochas- tic process to the problem of SR in annular discs. In conjunction with abrasive bearing surfaces, two distinct models of hydrodynamic lubrication are created. The first of these models is associated with a 1-D, azimuthal roughness, and the second model applies to a 1-D radial roughness. b) It has been proposed that disc coating has a two- part geometry. The large-scale portion of the film geometry, including any long wavelength distur- bances, is measured in the first component, which will be referred to as the thickness of the nominal film. The second component of the film geometry is the part caused by SR, which is calculated from the nominal level and is thought to change at random. c) By using Christensens stochastic theory, an aver- aged non-linear modified RE equation has been derived. d) Probability density function, expectation, average, and variance are employed to reveal the impacts of c∗ for P ∗, W ∗ and t∗. 440 vol. 63 no. 6/2023 Squeeze film lubrication between rough porous annular discs (a). 1Detail  2Radial roughness Azimuthal roughnessDetail (b). Figure 1. An illustrative picture showing the physical characteristics of a squeeze film design (A), with details of surface roughness of porous rough annular disks (B). In the world of engineering and material science, the investigation of porous surfaces with particular rough- ness effects on annular discs is a fresh and fascinating topic of study. In annular discs, the combination of the notions of surface porosity and roughness is new and has significant implications for numerous appli- cations. Here, we describe the novelty and possible objectives of this research work in this area. 1.2. Novelty Exploring the interactions between surface porosity and roughness within annular discs is exclusive. In- vestigators can produce materials and structures with different characteristics and capacities by combining these two features, which is not possible by researching them separately. Incorporating porous and rough sur- faces into circular discs can increase the effectiveness of heat transfer. Applications like heat exchangers, where an improved thermal performance can result in energy savings, are particularly interesting in this area. It is novel to understand how the interaction of porosity and roughness affects fluid flow within annular discs. The knowledge gained from this may be useful for fluidic devices, pumps, and turbines. To create specific porous and rough structures in annular discs for specialised purposes, such as biomedical de- vices, aircraft components, or filtration systems, novel surface modification techniques can be created. 1.3. Objective One of the primary objectives is to optimise the de- sign of porous and rough annular discs for specific applications. Researchers can aim to maximise heat transfer efficiency, minimise pressure drop, or enhance filtration performance, depending on the application requirements. Simulations can be done by developing an analytical model to predict how different combi- nations of porosity and roughness affect fluid flow patterns and structural integrity within annular discs. Overall, the objectives of the research in this area revolve around optimising the design, understanding the behaviour, and exploring the applications of an- nular discs with integrated porous and rough surfaces to meet the specific needs of various industries and technologies. 2. Problem formulation The squeeze film, produced by two porous annular discs with inner and outer radii r̂1 and r̂2, moving toward one another at a squeeze velocity (−∂ĥ/∂t), is shown in Figure 1. Porous media with thickness Ĥ0 and capillary radii R̂0 are present in the upper disc. 441 A. K. Rahul, M. K. Singh, R. Tiwari et al. Acta Polytechnica 3. Mathematical formulations and solutions The RFM is also known as the non-linear cubic power that correlates the shear stress τr̂ẑ and strain rate (∂û/∂ẑ). It is described mathematically as τr̂ẑ = µ̂ ∂û ∂ẑ − κτ3 r̂ẑ, (1) where µ̂ is the dynamic viscosity of the NNF. The non-linear component κ, which defines the feature of RFM, can be used to discriminate between three different types of fluids. • The fluids are referred to as pseudo-plastic if κ > 0, • The fluids are referred to as Newtonian if κ = 0, • The fluids are referred to as dilatant if κ < 0. The following continuity and momentum equations have been used to govern the governing equations for the RFM (Lin [8]; Rahul et al. [14, 15]). 1 r̂ ∂(r̂û) ∂r̂ + ∂ŵ ∂z = 0, (2) ∂p ∂r̂ − ∂τr̂ẑ ∂ẑ = 0, (3) ∂p ∂ẑ = 0. (4) Rough porous annular disc boundary conditions are described by Equations 5 and 6. a) At the upper surface: ẑ = 0 : û(r̂, 0, t̂) = 0, ŵ(r̂, 0, t̂) = ŵprs (5) b) At the lower surface: ẑ = ĥ : û(r̂, ĥ, t̂) = 0, ŵ(r̂, ĥ, t̂) = −∂ĥ ∂t̂ , (6) where û(r̂, ĥ, t̂) and ŵ(r̂, ĥ, t̂) are the components of velocity in r̂ and ẑ directions, respectively. ŵprs is the velocity of a through-flow on the upper boundary of the porous layer. The velocity component û is determined by sub- stituting Equation 1 by the above Equation 3 and integrating while taking boundary conditions Equa- tions 5 and 6 into consideration. û(r̂, ẑ) = 1 2µ̂ [ G(0)(ĥ, ẑ)∂p ∂r̂ + κ ( ∂p ∂r̂ )3 G(1)(ĥ, ẑ) ] , (7) By substituting Equation 7 into Equation 2 and integrating with respect to ẑ, we obtained 1 r̂ ∂ ∂r̂ { r̂ ẑ=ĥ∫ ẑ=0 1 2µ̂ [ G(0)(ĥ, ẑ)∂p ∂r̂ + κ ( ∂p ∂r̂ )3 G(1)(ĥ, ẑ) ] dẑ } = − ẑ=ĥ∫ ẑ=0 ∂ŵ ∂ẑ dẑ. (8) The modified NLRE for the rough porous annular discs is given by applying the boundary conditions (5) and (6) into Equation 8 as follows: 1 r̂ ∂ ∂r̂ [ r̂ĥ3 {( ∂p ∂r̂ ) + 3 20κĥ 2 ( ∂p ∂r̂ )3 }] = 12µ [ ∂ĥ ∂t̂ − ŵprs ] . (9) 3.1. Porous wall Let the non-Newtonian RF flow within the porous layer follow the modified Darcy’s law and the layer be isotropic and homogeneously distributed. The capillary network that makes up these homogeneous, isotropic porous layers has an average radius of R̂0 and porosity ψ̂. The axial and radial components of the velocity through the porous wall transform into ûp = ψ̂ µ̂ ( −∂p ∂r̂ ) + ψ̂ µ̂ κR̂2 0 6 ( −∂p ∂r̂ )3 , (10) ŵp = ψ̂ µ̂ ( −∂p ∂ẑ ) + ψ̂ µ̂ κR̂2 0 6 ( −∂p ∂ẑ )3 , (11) where ϕ̂ is the coefficient of porosity, ûp, ŵp and ψ̂ = ϕ̂R̂2 0/8 stand for the velocity components and permeability of porous layer, respectively. At the porous-fluids functionality, the cross-velocity compo- nent ŵprs is always continuous and equals ŵp. i.e. (ŵp)at ẑ=0 = (ŵprs)at ẑ=0. (12) Substituting Equation 12 into Equation 9, the result obtained is known as the modified NLRE which is expressed as: 1 r̂ ∂ ∂r̂ [ r̂ĥ3 {( ∂p ∂r̂ ) + 3 20κĥ 2 ( ∂p ∂r̂ )3 }] = 12µ̂ [ dĥ dt̂ − ψ̂ µ̂ ( −∂p ∂ẑ ) + κR̂2 0 6 ( −∂p ∂ẑ )3 ] ẑ=0 . (13) The modified version of Darcy’s law satisfied Equa- tion 2 for the porous layer (Walicka et al. [13], Rao and Rahul [14, 15]) 1 r̂ ∂(r̂ûp) ∂r̂ + ∂ŵp ∂ẑ = 0. (14) By substituting Equation 10 and Equation 11 into Equation 14, we obtain ∂ ∂ẑ [( −∂p ∂ẑ ) + κR̂2 0 6 ( −∂p ∂ẑ )3 ] = −1 r̂ ∂ ∂r̂ [ r̂ {( −∂p ∂r̂ ) + κR̂2 0 6 ( −∂p ∂r̂ )3 }] , (15) 442 vol. 63 no. 6/2023 Squeeze film lubrication between rough porous annular discs Integrating Equation 15 with respect to ẑ over the porous layer region [−Ĥ0, 0], [( −∂p ∂ẑ ) + κR̂2 0 6 ( −∂p ∂ẑ )3 ] ẑ=0 = −1 r̂ ∂ ∂r̂ r̂ 0∫ −Ĥ0 [( −∂p ∂r̂ ) + κR̂2 0 6 ( −∂p ∂r̂ )3 ] dẑ, (16) since [( −∂p ∂ẑ ) + κR̂2 0 6 ( −∂p ∂ẑ )3 ] ẑ=−Ĥ0 = 0, (17) For small Ĥ0, we obtain [( −∂p ∂ẑ ) + κR̂2 0 6 ( −∂p ∂ẑ )3 ] ẑ=0 ≈ −Ĥ0 r̂ ∂ ∂r̂ [ r̂ {( −∂p ∂r̂ ) + κR̂2 0 6 ( −∂p ∂r̂ )3 }] ẑ=0 . (18) Using the Morgan-Cameron approximation [11], the Equation 18 is used in Equation 13, the modified NLRE, which is defined as 1 r̂ ∂ ∂r̂ [ r̂ { f1(ĥ, ϕ̂, R̂0, Ĥ0)∂p ∂r̂ + 3 20κf2(ĥ, ϕ̂, R̂0, Ĥ0) ( ∂p ∂r̂ )3}] = 12µ̂ ( dĥ dt̂ ) , (19) where f1(ĥ, ϕ̂, R̂0, Ĥ0) = ĥ3 − 3 2 ϕ̂R̂ 2 0Ĥ0, f2(ĥ, ϕ̂, R̂0, Ĥ0) = ĥ5 + 5 3 ϕ̂R̂ 4 0Ĥ0. 3.2. Surface roughness In the case of surface roughness, the concept of film thickness is divided into two parts. ĥ = ĥ(r̂) + ĥs(r̂, θ, ξ), (20) where ĥ(r̂) represents the apparent smooth part of the film geometry, ĥs = δ1 + δ2 provides the arbi- trary region resulting from SR irregularities measured from the apparent level, and ξ represents the irregular variable that represents the apparent positive definite roughness arrangement. The Gaussian distribution is used to calculate the roughness profile heights that are valid up to three standard deviations or more for a range of lubricated rough surfaces. A likelihood rough distribution function is defined by Christensen [16]: f(ĥs) = { 35 32c7 (c2 − ĥ2 s)3, −c ≤ ĥs ≤ +c 0, elsewhere (21) where c denotes the maximum deviation from the mean film thickness, i.e., c = ±3σ, where σ is the standard deviation. Using expected values in Equa- tion 21, we obtain the following equation of the aver- aged NLRE: 1 r̂ ∂ ∂r̂ [ r̂ ( E { f1(ĥ, ϕ̂, R̂0, Ĥ0) } ∂E(p) ∂r̂ + 3 20κE { f2(ĥ, ϕ̂, R̂0, Ĥ0) }( ∂E(p) ∂r̂ )3 ] = 12µ̂dE(ĥ) dt̂ . (22) The expectation operator E(·) is expressed as fol- lows: E(·) = +c∫ −c (·)f(hs)dhs. (23) By using stochastic theory, Christensen [16] has proposed two different types of 1-D roughness pat- terns: azimuthal and radial patterns. The patterns of roughness are confined edges and valleys flowing in the same direction. The film thickness region for the one-dimensional radial and azimuthal roughness is represented as follows: ĥ = { ĥ(r̂) + ĥs(θ, ξ), radial roughness ĥ(r̂) + ĥs(r̂, ξ), azimuthal roughness. (24) For these roughness patterns, the modified stochas- tic NLRE is given by 1 r̂ ∂ ∂r̂ [ r̂ { G1(ĥ, ϕ̂, R̂0, Ĥ0)∂E(p) ∂r̂ + 3 20κG2(ĥ, ϕ̂, R̂0, Ĥ0) ( ∂E(p) ∂r̂ )3}] = 12µ̂dE(ĥ) dt̂ . (25) where G1(ĥ, ϕ̂, R̂0, Ĥ0, c) ={ E{f1(ĥ, ϕ̂, R̂0, Ĥ0)}, radial roughness[ E { 1/f1(ĥ, ϕ̂, R̂0, Ĥ0) }]−1 , azimuthal roughness. G2(ĥ, ϕ̂, R̂0, Ĥ0, c) ={ E{f2(ĥ, ϕ̂, R̂0, Ĥ0)}, radial roughness[ E { 1/f2(ĥ, ϕ̂, R̂0, Ĥ0) }]−1 , azimuthal roughness. (26) 443 A. K. Rahul, M. K. Singh, R. Tiwari et al. Acta Polytechnica Currently, the problem is minimising the construct- ing methods for estimating the left side of Equation 25 according to a specified roughness pattern. The mea- surement of mean FP involves the evaluation of the standard estimate of distinct film thicknesses. The probability density function is incorporated by Equa- tion 21. The accompanying expected estimations of film thickness is discussed in the work by Walicka et al. [23]. E(ĥ) = ĥ, E(ĥ2) = ĥ2 ( 1 + ∆2 9 ) , E(ĥ3) = ĥ3 ( 1 + ∆2 3 ) , E(ĥ4) = ĥ4 ( 1 + 2∆2 3 + ∆4 33 ) , E(ĥ5) = 5ĥ5 ( 1 + 2∆2 5 + ∆4 33 ) , E(ĥ6) = 5ĥ6 ( 1 5 + ∆2 93 + ∆4 11 + ∆6 429 ) , E(ĥ−1) = 1 ĥ [ 35 32 1 ∆7 { (∆2 − 1) ln (1 + ∆ 1 − ∆ ) − 2∆ 15 ( 15 − 40∆2 + 33∆4 )}] , E(ĥ−2) = 1 ĥ2 [ 35 32 1 ∆7 { 6(∆2 − 1) ln (1 + ∆ 1 − ∆ ) − 4∆ 5 ( 15 − 25∆2 + 8∆4 )}] , E(ĥ−3) = 1 ĥ3 [ 35 32 1 ∆7 { 3(5 − ∆2)(∆2 − 1) ln (1 + ∆ 1 − ∆ ) +2∆ ( 15 − 13∆2 )}] , E(ĥ−4) = 1 ĥ4 [ 35 32 1 ∆7 { 4(5 − 3∆2) ln (1 + ∆ 1 − ∆ ) − 8∆ 3 ( 15 − 4∆2 )}] , E(ĥ−5) = 1 ĥ5 [ 35 32 1 ∆7 { 3(∆2 − 5) ln (1 + ∆ 1 − ∆ ) + 2∆ 1 + ∆2 ( 15 − 13∆2 )}] , E(ĥ−6) = 1 ĥ6 [ 35 32 1 ∆7 { 3 ln (1 + ∆ 1 − ∆ ) + 2∆ 5(1 − ∆2) ( 15 − 25∆2 + 8∆4 )}] . (27) where ∆ = c ĥ . (28) The predicted values of film thickness feature E ( ĥ−i ) , i = 1, 2, . . . , 6, are not suitable for numerical purposes for small ∆ values, as there are variations when small amounts are involved, resulting in a loss in large digits. Hence, expansions of Taylors powers ∆ are proposed. The expansion of these film thickness functions in series E(ĥ−1) = 1 ĥ [ 1 + ∞∑ n=1 105∆2n (2n + 1)(2n + 3)(2n + 5)(2n + 7) ] , E(ĥ−2) = 1 ĥ2 [ 1 + ∞∑ n=1 105∆2n (2n + 3)(2n + 5)(2n + 7) ] , E(ĥ−3) = 1 ĥ3 [ 1 + ∞∑ n=1 105(n + 1)∆2n (2n + 3)(2n + 5)(2n + 7) ] , E(ĥ−4) = 1 ĥ4 [ 1 + ∞∑ n=1 35(n + 1)∆2n (2n + 5)(2n + 7) ] , E(ĥ−5) = 1 ĥ5 [ 1 + ∞∑ n=1 35(n + 1)(n + 2)∆2n 2(2n + 5)(2n + 7) ] , E(ĥ−6) = 1 ĥ6 [ 1 + ∞∑ n=1 7(n + 1)(n + 2)∆2n 2(2n + 7) ] . (29) The pressure boundary conditions for the NLRE are given by Equation 30p̂ = 0 at r̂ = r̂1 r̂2 , p̂ = 0 at r̂ = 1. (30) Introducing the dimensionless variables and param- eters: r∗ = r̂ r̂2 , h∗ = ĥ ĥ0 , H∗ = Ĥ ĥ0 , R∗ = R̂ ĥ0 , Ω = r̂1 r̂2 , c∗ = c ĥ0 , P ∗ = pĥ0 3 µ̂0r̂2 2(−dĥ/dt̂) , ω = κµ̂2 0r̂ 2 2(−dĥ/dt̂) ĥ4 0 . (31) Under the above mentioned dimensionless variables and parameters, the NLRE Equation 25 is expressed as: ∂ ∂r∗ [ r∗G∗ 1(h∗, R∗ 0, H ∗ 0 , c ∗) ( ∂P ∗ ∂r∗ ) + 3 20ωG ∗ 2(h∗, R∗ 0, H ∗ 0 , c ∗) ( ∂P ∗ ∂r∗ )3 ] = −12r∗, (32) where G1(h∗, R∗, H∗, c∗) ={ E{f1(h∗, R∗ 0, H ∗ 0 )}, radial roughness [E {1/f∗ 1 (h∗, R∗ 0, H ∗ 0 )}]−1 , azimuthal roughness (33) 444 vol. 63 no. 6/2023 Squeeze film lubrication between rough porous annular discs G2(h∗, R∗, H∗, c∗) ={ E{f2(h∗, R∗ 0, H ∗ 0 )}, radial roughness [E {1/f∗ 2 (h∗, R∗ 0, H ∗ 0 )}]−1 , azimuthal roughness (34) where ω represents the non-dimensional non-linear factor that specifies the RF behaviour. If, ω = 0 the fluids will behave like NF; ω > 0, the fluids will decompose as pseudoplastic fluids and ω < 0, the fluids will be dilatant fluids. The boundary condition in terms of dimensionless form is: P ∗ = 0 at r∗ = Ω, P ∗ = 0 at r∗ = 1. } (35) The value of the pseudoplastic coefficient κ depends on the type and quantity of additives which can be determined experimentally [7]. Thus, the values of ω can be calculated with the appropriate value of κ. However, for the validity of the present analysis, the value of ω is restricted to |ω| < 0.01. It is discov- ered that the dimensionless non-Newtonian averaged modified NLRE is nonlinear in terms of P ∗. A small perturbation method based on Lin [8], Rahul and Rao [14, 15, 24] is used to get approximate analytical solutions. Therefore, the classical perturbation method is used to solve it. The perturbation series for P can be expressed in the form: P ∗ = P ∗ 0 + ωP ∗ 1 + ω2P ∗ 2 + . . . (36) For ω ≪ 1, it is sufficient, for the analysis, to consider the first order term in ω as follows: For small values of the nonlinear parameter ω, the FP is perturbed: P ∗ = P ∗ 0 + ωP ∗ 1 , (37) For the higher values of ω, second and higher-order terms can be considered to increase the accuracy of the results. However, for the higher values of ω, it is more appropriate to adopt a numerical solution procedure such as the finite element method to solve the Reynolds equation. To obtain the FP P ∗ 0 and P ∗ 1 , we substitute Equa- tion 37 into the modified NLRE Equation 32 ∂ ∂r∗ [ r∗G∗ 1(h∗, R∗ 0, H ∗ 0 , c ∗) ( ∂P ∗ ∂r∗ )] + 12r∗ = 0, (38) ∂ ∂r∗ [ r∗G∗ 1(h∗, R∗ 0, H ∗ 0 , c ∗) ( ∂P ∗ ∂r∗ )] = − 3 20G ∗ 2(h∗, R∗ 0, H ∗ 0 , c ∗) ∂ ∂r∗ [ r ( ∂P ∗ 0 ∂r∗ )3 ] . (39) After solving Equation 38 and Equation 39, the perturbed dimensionless film pressures P ∗ 0 and P ∗ 1 are: P ∗ 0 = 3 G∗ 1(h∗, R∗ 0, H ∗ 0 , c ∗) [ 1 − r∗2 + ( Ω2 − 1 log Ω ) log r∗ ] , (40) P ∗ 1 = − A4 G∗ 1(h∗, R∗ 0, H ∗ 0 , c ∗) log r∗ − 3 20 G∗ 2(h∗, R∗ 0, H ∗ 0 , c ∗) G∗ 1(h∗, R∗ 0, H ∗ 0 , c ∗) {[ 54(1 − r∗4) − 54A4(1 − r∗2) − 18A2 4 log r∗] +1 2A 3 4 ( 1 r∗2 − 1 )} , (41) where A4 = −3(1 − Ω2) log Ω A2 = 3 20 G∗ 2(h∗, R∗ 0, H ∗ 0 , c ∗) G∗4 1 (h∗, R∗ 0, H ∗ 0 , c ∗) 1 log Ω { 54(1 − r∗4) − 54A4(1 − Ω2) −18A2 4 logS + 1 2A 3 4 ( 1 Ω2 − 1 )} . (42) The FP is integrated into the domain [r̂1, r̂2] to get the LCC of annular discs. W = 2π r̂2∫ r̂1 pr̂dr̂. (43) The dimensionless form of Equation 43 is given below: W ∗ = Ŵ ĥ3 0 µ̂0r̂4 2(−dĥ/dt) = 2π 1∫ Ω P ∗r∗dr∗. (44) Introducing the dimensionless SRT as: t∗ = Wĥ2 0 µ̂0r̂4 2 t̂, (45) Substituting Equation 45 into Equation 44 yields the ordinary differential equation that incorporates the film height, which varies with SRT. dh∗ dt∗ = − 1 W ∗ . (46) 445 A. K. Rahul, M. K. Singh, R. Tiwari et al. Acta Polytechnica 0.62 0.64 0.66 0.68 0.70 0.72 13.0 13.5 14.0 14.5 15.0 r* P * Ω = 0.0 HNewtonianL Ω = 0.0002 Ω = -0.0002H0 * = R0 * = c* = 0.0 HRedL Lin@8D H0 * = R0 * = 0.2 HGreen, AzimuthalL H0 * = R0 * = 0.2 HBlue, RadialL Red, Green, Blue 0.4 0.5 0.6 0.7 0.8 0.9 1.0 0 5 10 15 20 25 r* P * Figure 2. Variations in FP (P ∗) with a coordinate axis (r∗) for different values of ω, H∗ 0 and R∗ 0 with Ω = 0.4 and c∗ = 0.1 and h∗ = 0.3 . The non-dimensional film thickness t∗ = 0 is ini- tially under the condition of h∗ = 1. Integrate this equation while using the initial condition: t∗ = 1∫ h∗ 1 W ∗ dh ∗. (47) With the use of the Gaussian quadrature inte- gration algorithm, we were able to determine the properties of the discs by applying the boundary conditions and including the variable. 4. Results and Discussion The characteristics of a squeeze film between porous rough annular discs under RFM are investigated and the results of this study are represented in the form of graphs and tables. In the beginning of this study, Darcy’s law and Morgan-Cameron approximation are chosen for NNF in the porous matrix condition with homogeneous capillary tubes and the Christensen stochastic theory for surface roughness (SR) on an- nular discs. For this given special case, the value of ω → 0. ∂ ∂r∗ [ r∗ { G∗ 1(h∗, R∗ 0, H ∗ 0 , c ∗) ( ∂P ∗ ∂r∗ )}] = −12r∗. (48) As the porous, roughness, and viscosity variation parameters tend to zero (H∗ 0 → 0, R∗ 0 → 0, c∗ → 0), the results are obtained, as was also in the case of Lin [8]. Also, one can obtain the corresponding cases from the specific values of the parameters. • H∗ 0 = R∗ 0 = 0: the non-porous case; • c∗ = 0: smooth surface; • ω = 0: the NF case; and Some representative values being used for RF, poros- ity, and SR can be observed in many articles such as Lin [8], Walicka et al. [32], and Siddangouda et al. [9]. Ω = 0.0 HNewtonianL Ω = 0.0002 Ω = -0.0002 c* = 0.2 HGreen, AzimuthalL c* = 0.2 HBlue, RadialL c* = H0 * = R0 * = 0 HRedL Lin@8D Green, Red, Blue 0.4 0.5 0.6 0.7 0.8 0.9 1.0 0 5 10 15 20 25 r* P * Figure 3. Variations in FP (P ∗) with a coordinate axis (r∗) for different values of ω and c∗ with Ω = 0.4, h∗ = 0.3 and H∗ 0 = 0.1 and R∗ 0 = 0.1. Radial roughness H0 * = R0 * = 0.0 HBlueL HNon - PorousL H0 * = R0 * = 0.15 HMagentaL H0 * = R0 * = 0.0 HRedL HNon - PorousL H0 * = R0 * = 0.15 HGreenL Azimuthal roughness Ω = 0.0 HNewtonianL Ω = 0.0002 Ω = -0.0002 Red, Blue, Green, Magenta W = 0.4 0.20 0.25 0.30 0.35 0.40 0.45 0.50 0 50 100 150 h* W * Figure 4. Variations in LCC (W ∗) with h∗ for dif- ferent values of ω, H∗ 0 and R∗ 0 with c∗ = 0.1 and Ω = 0.4. In addition, the following set of parameters is used when performing numerical calculations: RF non- linear parameter: ω ∈ [−0.0002, 0.0002], roughness parameter: c∗ = 0, 0.1, 0.2, 0.3, 0.4 relative porous parameter: H∗ 0 = 0.0, 0.1, 0.2 and R∗ = 0.0, 0.1, 0.2. In this section, the solid lines demonstrate the dila- tant fluids, dot-dashed lines represent the pseudoplas- tic fluids, and the solid lines with a solid circle marker represent the Newtonian fluids. Figure 2 describes how the FP distribution over the porous hydrodynamic annular discs varies along with different r∗ values of the porous and RF parameters under fixed Ω = 0.4, h∗ = 0.3 and c∗ = 0.1. It is observed that FP de- creases with respect to increasing porosity, which is a similar trend to that observed by Lin [8]. In addition, the effect of RF on porous wall increases the FP for dilatant fluids whereas pseudoplastic fluids behaviour shows a reverse trend as compared to NF. Figure 3 shows the effect of squeeze FP (P ∗) along with the coordinate axis (r∗) for different values of ω and c∗ with Ω = 0.4, h∗ = 0.3, H∗ 0 = 0.1 and R∗ = 0.1. It is shown that FP decreases with c∗ for both types of roughness and finds a similar trend for RF with the smooth surface case of Lin [8]. Figures 4 and 5 illustrate the variation of LCC 446 vol. 63 no. 6/2023 Squeeze film lubrication between rough porous annular discs c* = 0.0 HBlueL Smooth c* = 0.2 HGreenL Azmuthal roughness c* = 0.2 HMagentaL Radial roughness Ω = -0.0002 Ω = 0.0 HNewtonianL Ω = 0.0002 Green, Blue, Magenta W = 0.4 0.20 0.25 0.30 0.35 0.40 0.45 0.50 0 50 100 150 h* W * Figure 5. Variations in LCC (W ∗) with h∗ for dif- ferent values of ω and c∗ with H∗ 0 = R∗ 0 = 0.1 and Ω = 0.4. Azimuthal roughness H0 * = R0 * = 0.0 HRedL HNon - PorousL H0 * = R0 * = 0.15 HGreenL Radial roughness H0 * = R0 * = 0.0 HBlueL HNon - PorousL H0 * = R0 * = 0.15 HMagentaL Ω = 0.0 HNewtonianL Ω = 0.0002 Ω = -0.0002 Red, Blue, Green, Magenta W = 0.4 0.20 0.25 0.30 0.35 0.40 0.45 0.50 0 20 40 60 80 100 120 140 h* t* Figure 6. Variations in SRT (t∗) with h∗ for different values of ω, H∗ 0 R∗ 0 with c∗ = 0.1 and Ω = 0.4. (W ∗) as a function of film thickness (h∗), influenced by the presence of porosity (H∗ 0 and R∗ 0) and surface roughness (c∗) between annular discs. These varia- tions are observed using the Rabinowitsch fluid model parameter (ω), while maintaining a constant radius ratio Ω = 0.4. From these Figures 4 and 5, it is con- cluded that the effect of roughness plays a vital role in the LCC of the annular discs. Compared to the cases of the non-porous and smooth surface of discs, it is concluded that the impact of porosity and rough- ness reduces the LCC with increased film thickness. LCC also considerably decreases with the increasing porous parameter. Mostly in the case of roughness patterns, LCC is much more prominent as compared to a smooth surface for the azimuthal roughness. For lower and higher values of radius ratio, moreover, the effect of the dilatant fluids provides a larger LCC influence on the pseudoplastic fluids and is expected to yield inverse outcomes as compared to NF. Fig- ures 6 demonstrates the the disparity of dimensionless SRT (t∗) concerning dimensionless film thickness (h∗) for distinct values of the non-linear parameter (ω), porous parameters (H∗ 0 , and R∗ 0) under the radius ratio Ω = 0.4. The impact of the porous media on discs, SRT is hugely prominent for the non-porous (H∗ 0 = 0.0 , R∗ 0 = 0.0) as compared to the porous c* = 0.2 HGreenL Azimuthal roughness c* = 0.0 HBlueL Smooth c* = 0.2 HMagentaL Radial roughness Ω = -0.0002 Ω = 0.0 HNewtonianL Ω = 0.0002 Green, Blue, Magenta W = 0.4 0.20 0.25 0.30 0.35 0.40 0.45 0.50 0 20 40 60 80 100 120 140 h* t* Figure 7. Variations in SRT (t∗) with h∗ for different values of ω, and c∗ with H∗ 0 = R∗ 0 = 0.1 and Ω = 0.4. case of discs. Influence of RF in the porous medium between annular discs, the SRT is lower for pseudo- plastic and higher for dilatant fluids as compared to NF. Figures 7 shows the behaviour of SRT t∗ against h∗ with ω and c∗ under the radius ratio Ω = 0.4. It is found that t∗ increases (decreases) for the azimuthal (Radial) roughness patterns. The SR accounted for more SRT with the lower radius ratio with respect to higher values of the radius ratio. Tables 1 and 2 illustrate the numerical comparison of LCC with the results of Lin [8] (with H∗ 0 = R∗ 0 = 0.0, h∗ = 0.2) and the present analysis for different values of ω, c∗ and Ω. It can be concluded that increas- ing the roughness and porous parameters leads to a decrease in LCC. Moreover, dilatant fluids (−0.0002) showed a higher LCC as compared to pseudoplastic fluids (+0.0002). 5. Conclusions Based on the RFM, the performance of annular discs under the condition of SR and the porous wall is studied in the present work proposing two different theories. • The modified Darcy’s law and Morgan-Cameron approximation for porous wall; • Christensen’s stochastic theory for SR. The non-linear Reynolds equation (NLRE) is con- structed using the above theories. To solve the NLRE, we used the techniques listed below: • Small perturbation techniques; • Five-point Gauss quadrature integral formula. The annular discs are studied in two physical situ- ations with surfaces: one with porous walls of equal thickness and radius of the capillarity tube, another with SR. FP, LCC, and SRT are calculated for dif- ferent parameters of porous and SR. In the case of non-porous and smooth surfaces, the FP distribution, LCC, and SF are identical to the literature Lin [8]. The conclusions are as follows: 447 A. K. Rahul, M. K. Singh, R. Tiwari et al. Acta Polytechnica Ω = 0.2 Ω = 0.4 ω Lin [8] Present analysis Lin [8] Present analysis c∗ = 0.0 c∗ = 0.0 c∗ = 0.4 c∗ = 0.4 c∗ = 0.0 c∗ = 0.0 c∗ = 0.4 c∗ = 0.4 (Radial) (Azimuthal) (Radial) (Azimuthal) (Radial) (Azimuthal) (Radial) (Azimuthal) −0.0002 329.088 329.088 69.7214 459.684 138.991 138.991 32.6894 197.561 0 250.804 250.804 66.5346 364.434 120.365 120.365 31.9312 174.899 +0.0002 172.519 172.519 63.3478 269.185 101.74 101.74 31.173 152.237 Table 1. Numerical comparison of LCC (W ∗) of Lin [8] (with H∗ 0 = R∗ 0 = 0.0, h∗ = 0.2) and present analysis (with H∗ 0 = R∗ 0 = 0.1, h∗ = 0.2) for different values of ω, c∗ and Ω. Ω = 0.2 Present analysis Ω = 0.4 Present analysis ω Lin[8] H∗ 0 = R∗ 0 = 0.2 Lin[8] H∗ 0 = R∗ 0 = 0.2 H∗ 0 = 0.0 H∗ 0 = 0.0 H∗ 0 = 0.2 H∗ 0 = 0.2 H∗ 0 = 0.0 H∗ 0 = 0.0 H∗ 0 = 0.2 H∗ 0 = 0.2 R∗ 0 = 0.0 R∗ 0 = 0.0 R∗ 0 = 0.2 R∗ 0 = 0.2 R∗ 0 = 0.0 R∗ 0 = 0.0 R∗ 0 = 0.2 R∗ 0 = 0.2 (Radial) (Azimuthal) (Radial) (Azimuthal) (Radial) (Azimuthal) (Radial) (Azimuthal) −0.0002 329.088 329.088 66.4727 72.9544 138.991 138.991 31.5269 34.5767 0 250.804 250.804 64.9249 71.155 120.365 120.365 31.1587 34.1486 +0.0002 172.519 172.519 63.377 69.3557 101.74 101.74 30.7904 33.7205 Table 2. Numerical comparison of LCC (W ∗) of Lin [8] (with h∗ = 0.2, c∗ = 0.0) and present analysis (with h∗ = 0.2, c∗ = 0.1) for different values of ω, H∗ 0 , R∗ 0, and Ω. a) Annular discs with the presence of a porous wall: The FP, LCC, and SRT decrease as compared to non-porous discs, as in the case of Lin [8]. b) Annular discs with the presence of SR: The FP, LCC, and SRT decrease as compared with smooth discs, as in the case of Lin [8]. All these characteris- tics show a lower value for radial surface roughness (SR). c) Annular discs with RFM: The performance of an- nular discs is better for dilatant fluids and lower for pseudoplastic fluids as compared to NF. d) A radial roughness pattern on the porous surface of an annular disc can promote better fluid mixing and increase heat transfer efficiency. Radial rough- ness, on the other hand, can cause higher pressure drop and increased flow resistance, which can be a disadvantage in certain applications where energy efficiency and low-pressure losses are critical. e) An annular disc’s circumferential azimuthal rough- ness pattern can influence the flow of fluid in one direction, the amount of turbulence present inside the disc, and the flow velocity in various annular disc sections. 6. Future research In conclusion, the study of porous roughness annular discs is an important are of research with several ap- plications in fluid dynamics, materials science, and engineering. Several promising directions stand out as having high potential for future research in this area, including advanced materials and fabrication techniques, optimisation of porous structure, numeri- cal modelling, and simulation, applications in energy and environmental engineering, surface modification techniques, data-driven approaches, etc. We can learn much more about porous roughness annular discs and unlock their potential to solve challenging problems across a range of industries by incorporating these research directions. Researchers can help make future engineering solutions more effective and sustainable by improving their design, performance, and applica- tions. List of symbols Acronyms FP Film pressure NF Newtonian fluids NLRE Non-linear Reynolds equation NNF Non-Newtonian fluids RFM Rabinowitsch fluid model SF Squeeze film SR Surface roughness SRT Squeeze response time Symbols r̂1, r̂2 The outer/inner radius of annular discs [m] Ω The ratio of inner and outer radius, Ω = r̂1/r̂2 [–] ĥ0, ĥ1 Inlet/outlet film thickness [m] h Nominal film height [m] ĥ, h∗ Film thickness (ĥ = h+ hs), h∗ = ĥ/ĥ0 [m] hs Film height deviation from the nominal level [m] E(·) Expected value [–] c, c∗ The greatest asperity variation from the nominal height, c∗ = c/ĥ0 [–] dĥ/dt̂ Squeeze velocity [m s−1] 448 vol. 63 no. 6/2023 Squeeze film lubrication between rough porous annular discs r̂, ẑ, r∗ Coordinates of discs, r∗ = r̂/r̂2 [–] ŵprs Through-flow velocity on the upper boundary of the porous layer [m s−1] û, ŵ Velocity components in r̂ and ẑ directions [m s−1] ûp, ŵp Axial and radial velocity component of the porous region. [m s−1] Ĥ0, H ∗ 0 Porous pad thickness, H∗ 0 = ϕ̂Ĥ0/ĥ0 [m] R̂0, R ∗ 0 Porous pad thickness, H∗ 0 = ϕ̂R̂0/ĥ0 [m] p Film pressure in the porous region [Pa] p, P ∗ Squeeze film pressure, P ∗ = pĥ3 0/µ̂0r̂2 2(−∂ĥ/∂t̂) [Pa] W,W ∗ Load-carrying capacity, W ∗ = W ĥ3 0/µ̂0r̂4 2(−∂ĥ/∂t̂) [N] W,W ∗ Squeeze response time, t∗ = Wĥ2 0t̂/µ̂0r̂ 4 2 [s] µ̂ Dynamic viscosity of the NNF [Po] τr̂ẑ Element of the stress tensor [Pa] ψ̂ Permeability of the porous region [m2] µ̂ Coefficient of porosity [pu] ϕ̂ Dynamic viscosity of the NNF [Po] κ Nonlinear factor accounting for the Rabinowitsch fluid model [–] ω Dimensionless non-linear factor, ω = κµ̂2 0r̂2 2(−∂ĥ/∂t̂)/ĥ4 0 [–] θ Circumferential co-ordinate [°] ξ Random variable for surface roughness [–] References [1] C. Allen, A. McKillop. An investigation of the squeeze film between rotating annuli. ASME Journal of Lubrication Technology 92(3):435–441, 1970. https://doi.org/10.1115/1.3451435 [2] N. Naduvinamani, A. Siddangouda, A. Kadadi, S. Biradar. Effect of pressure dependent viscosity on squeeze film characteristics of micropolar fluid in convex curved plates. Tribology-Materials, Surfaces & Interfaces 9(3):154–158, 2015. https://doi.org/10.1080/17515831.2015.1107239 [3] H. Spikes. The behaviour of lubricants in contacts: current understanding and future possibilities. Proceedings of the Institution of Mechanical Engineers, Part J: Journal of Engineering Tribology 208(1):3–15, 1994. https: //doi.org/10.1243/PIME_PROC_1994_208_345_02 [4] J. Lin, C. Hung. Combined effects of non-Newtonian rheology and rotational inertia on the squeeze film characteristics of parallel circular discs. Proceedings of the Institution of Mechanical Engineers, Part J: Journal of Engineering Tribology 222(4):629–636, 2008. https://doi.org/10.1243/13506501JET39 [5] Y. Wang, J. H. Wu, L. Xu. Influence of power-law fluid on transient performance of liquid film seal based on the time-dependent non-Newtonian dynamic reynolds equation. Tribology International 159:106984, 2021. https://doi.org/10.1016/j.triboint.2021.106984 [6] S. Wada, H. Hayashi. Hydrodynamic lubrication of journal bearings by pseudo-plastic lubricants: part 1, theoretical studies. Bulletin of JSME 14(69):268–278, 1971. https://doi.org/10.1299/jsme1958.14.268 [7] S. Wada, H. Hayashi. Hydrodynamic lubrication of journal bearings by pseudo-plastic lubricants: part 2, experimental studies. Bulletin of JSME 14(69):279–286, 1971. https://doi.org/10.1299/jsme1958.14.279 [8] J.-R. Lin. non-Newtonian squeeze film characteristics between parallel annular disks: Rabinowitsch fluid model. Tribology international 52:190–194, 2012. https://doi.org/10.1016/j.triboint.2012.02.017 [9] A. Siddangouda, N. Naduvinamani, S. Siddapur. Effect of surface roughness on the static characteristics of inclined plane slider bearing: Rabinowitsch fluid model. Tribology-Materials, Surfaces & Interfaces 11(3):125–135, 2017. https://doi.org/10.1080/17515831.2017.1347745 [10] N. Bujurke, M. Jagadeeswar, P. Hiremath. Analysis of normal stress effects in a squeeze film porous bearing. Wear 116(2):237–248, 1987. https://doi.org/10.1016/0043-1648(87)90236-5 [11] V. Morgan, A. Cameron. Mechanism of lubrication in porous metal bearings. In International Proceeding Conference Lubrication and Wear, Institute of Mechanical Engineering, vol. 89, pp. 151–157. Inst. Mech. Eng London, 1957. [12] N. C. Patel, J. R. Patel, G. M. Deheri. A comparative study of ferrofluid lubrication on double-layer porous squeeze curved annular plates with slip velocity. Acta Polytechnica 62(4):488–497, 2022. https://doi.org/10.14311/AP.2022.62.0488 [13] M. Bhat, G. Deberi. Squeeze film behaviour in porous annular discs lubricated with magnetic fluid. Wear 151(1):123–128, 1991. https://doi.org/10.1016/0043-1648(91)90352-U [14] A. Walicka, E. Walicki, P. Jurczak, J. Falicki. Curvilinear squeeze film bearing with rough surfaces lubricated by a Rabinowitsch–Rotem–Shinnar fluid. Applied Mathematical Modelling 40(17-18):7916–7927, 2016. https://doi.org/10.1016/j.apm.2016.03.048 [15] P. Rao, A. Rahul, S. Agarwal. Effect of non-Newtonian lubrication of squeeze film conical bearing with the porous wall operating with Rabinowitsch fluid model. Proceedings of the Institution of Mechanical Engineers, Part J: Journal of Engineering Tribology 232(10):1293–1303, 2018. https://doi.org/10.1177/1350650117749735 [16] P. Rao, A. Rahul. Combined effect of viscosity variation and non-Newtonian Rabinowitsch fluid in wide parallel rectangular-porous plate with squeeze-film characteristics. Meccanica 54(15):2399–2409, 2019. https://doi.org/10.1007/s11012-019-01092-2 [17] P. S. Rao, A. Kumar Rahul. Effect of viscosity variation on non-Newtonian lubrication of squeeze film conical bearing having porous wall operating with Rabinowitsch fluid model. Proceedings of the Institution of Mechanical Engineers, Part C: Journal of Mechanical Engineering Science 233(7):2538–2551, 2019. https://doi.org/10.1177/0954406218790 [18] A. K. Rahul, M. K. Singh, S. Saha. Squeeze film lubrication analysis and optimization of porous annular disk with viscosity variation of non-Newtonian fluid: Rabinowitsch fluid model. Tribology International 179:108060, 2023. https://doi.org/10.1016/j.triboint.2022.108060 449 https://doi.org/10.1115/1.3451435 https://doi.org/10.1080/17515831.2015.1107239 https://doi.org/10.1243/PIME_PROC_1994_208_345_02 https://doi.org/10.1243/PIME_PROC_1994_208_345_02 https://doi.org/10.1243/13506501JET39 https://doi.org/10.1016/j.triboint.2021.106984 https://doi.org/10.1299/jsme1958.14.268 https://doi.org/10.1299/jsme1958.14.279 https://doi.org/10.1016/j.triboint.2012.02.017 https://doi.org/10.1080/17515831.2017.1347745 https://doi.org/10.1016/0043-1648(87)90236-5 https://doi.org/10.14311/AP.2022.62.0488 https://doi.org/10.1016/0043-1648(91)90352-U https://doi.org/10.1016/j.apm.2016.03.048 https://doi.org/10.1177/1350650117749735 https://doi.org/10.1007/s11012-019-01092-2 https://doi.org/10.1177/0954406218790 https://doi.org/10.1016/j.triboint.2022.108060 A. K. Rahul, M. K. Singh, R. Tiwari et al. Acta Polytechnica [19] H. Christensen. Stochastic models for hydrodynamic lubrication of rough surfaces. Proceedings of the institution of mechanical engineers 184(1):1013–1026, 1969. https: //doi.org/10.1243/PIME_PROC_1969_184_074_02 [20] H. Christensen. Some aspects of the functional influence of surface roughness in lubrication. Wear 17(2):149–162, 1971. https://doi.org/10.1016/0043-1648(71)90025-1 [21] S. Alla Baksh, H. Bannihalli Naganagowda. Study of surface roughness with MHD and couple stress fluid on porous curved annular plates 62(6):574–588, 2022. https://doi.org/10.14311/AP.2022.62.0574 [22] V. D’Agostino, A. Ruggiero, A. Senatore. Unsteady oil film forces in porous bearings: analysis of permeability effect on the rotor linear stability. Meccanica 44(2):207–214, 2009. https://doi.org/10.1007/s11012-008-9180-0 [23] Y. D. Vashi, R. M. Patel, G. B. Deheri. Neuringer-Roseinweig model based longitudinally rough porous circular stepped plates in the existence of couple stress. Acta Polytechnica 60(3):259–267, 2020. https://doi.org/10.14311/AP.2020.60.0259 [24] A. Siddangouda, T. Biradar, N. Naduvinamani. Combined effects of surface roughness and viscosity variation due to additives on long journal bearing. Tribology-Materials, Surfaces & Interfaces 7(1):21–35, 2013. https://doi.org/10.1179/1751584X13Y.0000000024 [25] M. M. Munshi, A. R. Patel, G. B. Deheri. A study of ferrofluid lubrication based rough sine film slider bearing with assorted porous structure. Acta Polytechnica 59(2):144–152, 2019. https://doi.org/10.14311/AP.2019.59.0144 [26] M. Shimpi, G. Deheri. Magnetic fluid based squeeze film in rough rotating curved porous annular plates: Deformation effect. International Journal of Mathematical and Computational Sciences 7(8):1370– 1380, 2013. https://doi.org/10.1155/2012/148281 [27] S. Patel, G. Deheri, J. Patel. Ferrofluid lubrication of a rough porous hyperbolic slider bearing with slip velocity. Tribology in Industry 36(3):259–268, 2014. [28] J. Patel, G. Deheri. Combined effect of slip velocity and roughness on the Jenkins model based ferrofluid lubrication of a curved rough annular squeeze film. Journal of Applied Fluid Mechanics 9(2):855–865, 2016. https: //doi.org/10.18869/acadpub.jafm.68.225.24447 [29] Y. D. Vashi, R. M. Patel, G. B. Deheri. Neuringer-Roseinweig model based longitudinally rough porous circular stepped plates in the existence of couple stress. Acta Polytechnica 60(3):259–267, 2020. https://doi.org/10.14311/AP.2020.60.0259 [30] N. C. Patel, J. R. Patel, G. Deheri. An effect of a porous structure, slip velocity and Rosensweig’s viscosity on the ferrofluid based squeeze film in porous curved annular plates. Journal of Nanofluids 12(2):498– 505, 2023. https://doi.org/10.1166/jon.2023.1906 [31] P. S. Rao, A. K. Rahul. Pressure generation in rough conical bearing using non-Newtonian Rabinowitsch fluid with variable viscosity. Industrial Lubrication and Tribology 71(3):357–365, 2019. https://doi.org/10.1108/ILT-01-2018-0035 [32] A. Walicka, E. Walicki, P. Jurczak, J. Falicki. Influence of wall porosity and surfaces roughness on the steady performance of an externally pressurized hydrostatic conical bearing lubricated by a Rabinowitsch fluid. International Journal of Applied Mechanics and Engineering 22(3):717–737, 2017. https://doi.org/10.1515/ijame-2017-0045 [33] A. K. Rahul, P. S. Rao. Rabinowitsch fluid flow with viscosity variation: application of porous rough circular stepped plates. Tribology International 154:106635, 2021. https://doi.org/10.1016/j.triboint.2020.106635 [34] A. K. Rahul, M. K. Singh, S. Paul, et al. Performance analysis of annular disks with non-Newtonian Rabinowitsch fluid model: Influence of squeeze film pressure, surface roughness, porosity and viscosity variation. International Journal of Modern Physics B p. 2450268, 2023. [Online Ready]. https://doi.org/10.1142/S0217979224502680 450 https://doi.org/10.1243/PIME_PROC_1969_184_074_02 https://doi.org/10.1243/PIME_PROC_1969_184_074_02 https://doi.org/10.1016/0043-1648(71)90025-1 https://doi.org/10.14311/AP.2022.62.0574 https://doi.org/10.1007/s11012-008-9180-0 https://doi.org/10.14311/AP.2020.60.0259 https://doi.org/10.1179/1751584X13Y.0000000024 https://doi.org/10.14311/AP.2019.59.0144 https://doi.org/10.1155/2012/148281 https://doi.org/10.18869/acadpub.jafm.68.225.24447 https://doi.org/10.18869/acadpub.jafm.68.225.24447 https://doi.org/10.14311/AP.2020.60.0259 https://doi.org/10.1166/jon.2023.1906 https://doi.org/10.1108/ILT-01-2018-0035 https://doi.org/10.1515/ijame-2017-0045 https://doi.org/10.1016/j.triboint.2020.106635 https://doi.org/10.1142/S0217979224502680 Acta Polytechnica 63(6):439–450, 2023 1 Introduction 1.1 Research gap 1.2 Novelty 1.3 Objective 2 Problem formulation 3 Mathematical formulations and solutions 3.1 Porous wall 3.2 Surface roughness 4 Results and Discussion 5 Conclusions 6 Future research List of symbols References