Abstract
This paper presents a novel analytical solution for journal-bearing viscoelastic lubrication using the perturbation method. The nonlinear Giesekus model was used for the constitutive equations to study the effects of fluid elasticity, shear-thinning viscometric functions, and strain-hardening elongational viscosity of viscoelastic lubrication. The investigation focuses on the impact of characteristic parameters such as mobility factor, eccentricity ratio, and Weissenberg number on the fluid film pressure distribution, load capacity, and shear stress. Although distinguishing between the normal stress differences and extensional viscosity in mixed viscoelastic flows is complicated, we investigated the role and contribution of these two factors. By increasing the elasticity of the fluid, the portion of both mentioned parameters increases consequently. Furthermore, analyses and comparisons show the contributions of the first normal stress and elongational viscosity to the load capacity of the bearing through the stress ratio and flow type parameter for the first time. The research findings indicate that fluid elasticity enhances the load capacity of the bearing compared to a Newtonian lubricant with the same effective viscosity. Moreover, the bearing load capacity is divided into two regions. In the linear region, the mobility factor and Weissenberg numbers have minimal effects leading to a linear increase in the load distribution, and in the exponential region, the load capacity changes are considerable. This research provides valuable insights into the behavior of viscoelastic lubrication in journal-bearing systems.
Similar content being viewed by others
Introduction
Despite decades of research, finding an exact analytical solution for the Reynolds equation governing journal-bearing lubrication remains a challenging problem. The primary objectives of lubrication include reducing friction, preventing heat generation, and minimizing wear on rotating components. To achieve these objectives, a thin film of lubricant separates the shaft from the bearing. Tower1,2,3 conducted experiments on lubrication with meticulous precision inspired by an unexpected accident or mistake. This ultimately led to the development of hydrodynamic lubrication, the foundation of modern lubrication theory. He also calculated the pressure distribution of the fluid film and reported the results.
Reynolds4 used the brief forms of NavierâStokes and continuity equations to obtain a second-order differential equation for the pressure distribution of lubrication film. Sommerfeld5 was the first to introduce a series of dimensionless numbers relevant to bearings, among which the Sommerfeld number stands out as particularly important. Solving the lubrication equation is a complex problem that has captured the interest of numerous researchers. Dowson6 studied the generalized Reynolds equation for the lubricant layer. Raimondi and Boyd7,8,9 used an iterative method to solve the Reynolds equation numerically. And in 1958, they presented three papers on the numerical solution of the Reynolds equation for journal bearings. A large volume of data was available for designers for the first time.
Using the finite difference method, Raimondi10 investigated the problem of hydrodynamic lubrication of journal bearings for compressible lubricants. The results of Raimondiâs numerical solution were very accurate and had an acceptable agreement with the experimental results. Reddi11 solved the problem of lubrication for the incompressible lubricant using the finite element numerical method. He described the advantages of finite element methods over other numerical methods. Reddi and Chu12 then generalized the solutions for the compressible lubricant using the finite element method.
Malik13 used the finite element method to investigate gas-bearing lubrication. He applied a slip condition to the surfaces and examined the different parameters of the bearing. Dinariev14 extended the Reynolds theory for the hydrodynamic lubrication problem in the presence of a viscoelastic lubricant. Based on his qualitative and numerical analysis, the relaxation time effect expanded the pressure distribution and reduced the bearing capacity of a lubricating layer. Urreta et al.15 summarize the work developing hydrodynamic lubricated journal bearings with magnetic fluids. Sfyris and Chasalevris16 obtained an exact solution for the Reynolds equation while assuming a Newtonian fluid and considering time dependency. They implemented this solution for a finite journal bearing. They used their exact solution results and examined the characteristics of journal bearing while assuming a Newtonian fluid17. Rao et al.18 analytically solved the Reynolds equation for sliding and journal bearings while considering slip boundary conditions. Vignolo et al.19 solved the Reynolds equation for a finite-length bearing. They used the perturbation method and regarded the term \(\upvarepsilon ={\left(L/D\right)}^{2}\) as the perturbation parameter. Gong et al.20 continued Gustavoâs solution and solved the unsteady Reynolds equation.
Another common subject in bearing studies is misalignment. Jang and Khonsari21 studied the misaligned journal bearing and investigated the effect of the deviation on the static and dynamic characteristics of the bearing. Most studies dealing with bearing lubrication start with the well-known Reynolds equation. This equation is known as the classic lubrication theory of fluid film and is based on the Newtonian fluid. However, in many applications of lubricants in different industries, the approximation of the Newtonian fluid model cannot be used for the lubricants.
Based on the literature, using a few long-chain polymeric additives could improve the lubrication with much higher efficiency. In the mid-1950s, converting Newtonian lubricants to non-Newtonian lubricants by adding mineral oils became an accepted method22. A decade of studies has determined the effects of lubricants on polymer structures. In most cases, polymeric lubricants perform better than Newtonian ones23. One of the early investigations on finite-width journal bearing using non-Newtonian fluid was done by Tayal et al.24. They used the Prandtl model to study the shear thinning and shear thickening of lubricant in cylindrical coordinates. Horowitz and Steidler25 examined the effect of non-Newtonian oil on finite bearing characteristics by employing a logarithmic function for viscosity alternations. Tanner26 focused on a power-law fluid for a short journal bearing.
Wada and Hayashi27,28 investigated the necessary parameters of a journal bearing for a pseudo-plastic fluid assuming a finite width. Swamy29 computed the load capacity of non-Newtonian lubricants in limited-length journal bearings and showed that using non-Newtonian lubricants increases the load capacity. They examined the damping properties of non-Newtonian lubricants in journal bearings besides the effect of non-Newtonian lubricants on fluid instability. They reported that non-Newtonian lubricants are more stable than Newtonian ones30. Raghunandana et al.31 employed the non-Newtonian model developed by Dien and Elrod to find the stability margin for different support parameters of non-Newtonian lubricants.
Das et al.32,33 examined the performance of a hydrodynamic journal bearing under a micropolar lubricant. Abdel-Rahman34 studied the flow of non-Newtonian lubricant through a conical bearing in an external magnetic field. Elsharkawy35 used a finite difference scheme to solve a modified form of the Reynolds equation to investigate the effect of geometry, pressure distribution, load carrying capacity, side leakage flow, and friction factor. The results demonstrated that additives increase the load-carrying capacity but decrease friction and side leakage coefficient. Tian et al.36 simulated the non-Newtonian fluid between eccentric cylinders with the finite element method. In their study, the flow in the annular gap between the eccentric rotating cylinders was a shear-extensional controllable flow. Sakim et al.37 showed that couple stresses improve the load capacity and reduce friction, while permeability and deformation of a porous elastic journal bearing have reverse trends. Chetti and Zouggar38 presented a numerical study of the effect of elastic deformation on the static characteristics of a circular journal bearing operating with non-Newtonian fluids obeying the power law model. Li et al.39 numerically solved the generalized Reynolds equation, heat conduction, and energy equations to study variations in bearing performance due to misalignment. They also showed that thermal effect and surface roughness considerably influence lubrication performance. Gwynllyw and Phillips40 investigated the effect of relaxation time and gap size considering the PTT and the Oldroyd-B constitutive equation. Gertzos et al.41 simulated a bearing with Bingham fluid as the lubricant using the commercial software Fluent. Lin42 examined the non-Newtonian effect of couple stresses of short journal bearing. The Study showed that decreasing values of the system parameter for constant couple stress parameters can shift sub-critical bifurcation into super-critical bifurcation. Wierzcholski43 used the Rivlin-Ericksen constitutive equation to model and solve the lubrication problem for a micro-bearing. Guemmadi and Ouibrahim44 investigated the behavior of the generalized Maxwell viscoelastic fluid for lubrication of a journal bearing by applying a finite volume method.
Tichy45 used Maxwell's convective equation to study the effect of the Deborah number on the pressure distribution. He solved the problem concerning the perturbation theory, and the Deborah number was considered as perturbation parameter. Huang et al.46 calculated the pressure distribution for sliding and journal bearings, using a second-order fluid model. By varying the thickness of the fluid film, they discovered that the normal stress distribution was strictly dependent on the thickness. When it decreased, the normal stresses for a second-order fluid increased. Akyildiz and Bellout47 investigated the pressure distribution of a slider bearing for a PTT lubricant. They focused on the effect of the Deborah number on the pressure field. Kumar and Sharma48 evaluated a conical hybrid journal bearing with micro-grooves while considering the shear-thinning and piezo-viscous behaviors of the lubricant. Chetti et al.49 presented a theoretical study of the effects of elastic deformation and viscosity variation with pressure on the performance characteristics of a circular journal bearing lubricated with non-Newtonian fluids, considering the Barus law and the power law model. Mokhtari et al.50 looked into the Deborah number's impact on the pressure distribution and other characteristics using a FENE-P lubricant. Ahmed and Biancofiore51 proposed a new modeling technique based on lubrication theory, considering the viscoelastic effects. As a result, they obtained a modified equation for the pressure, i.e., the viscoelastic Reynolds (VR) equation. They also52 extended the VR approach to the non-linear finitely extensible non-linear elastic (FENE) type constitutive relations that account for the finite extension of the polymer chains and shear thinning. Soni53 investigated the numerical solution to predict the dynamic performance of finite bearing considering the combined influence of turbulence regime and non-Newtonian flow. Agrawal and Sharma54 examined the performance of micro-grooved hole-entry hybrid spherical thrust bearing (HSTB), considering the non-Newtonian behavior of the lubricant. Hashemabadi and Mirnajafizadeh55 solved the lubrication problem analytically by considering the SPTT model as representing fluid. They investigated the effects of some variables on velocity and pressure distribution. They also examined the impacts of changing the slope of the slider bearing.
Li et al.56 considered the upper-convected Maxwell model (UCM) to study the lubrication problem. In addition to the Deborah number, they also took interest in eccentricity (ε) effects. They showed that the viscoelastic property of the fluid causes increase of the lubricantâs pressure distribution and positively affects the lubrication procedure. Nessil et al.57 analyzed various aspects of lubrication theory, including heat transfer for a power-law fluid, and discussed the hydrodynamic parameters of the bearing. They concluded that the power-law index (n) plays a significant role in temperature distribution. Further, they examined the load capacity, pressure distribution, temperature distribution, and frictional force. Li58 conducted a study on non-Newtonian lubrication. The PTT constitutive equation represented the lubricant fluid, and they considered the Deborah number as the perturbation parameter. Soni and Vakharia59 studied turbulence effects in addition to the non-Newtonian influence of lubricants in cylindrical coordinates for a finite journal bearing. Ng and Pan's linear turbulence model was applied in the finite element method to analyze the problem. Abbaspur et al.60 analyzed the viscoelastic lubrication of a thrust bearing using the Giesekus constitutive equation. In their work, Cartesian coordinates were employed, while the present investigation necessitated the use of polar coordinates due to the journal bearing geometry. They used perturbation theory and considered the mobility factor as the perturbation parameter. From the geometry perspective, it is worth mentioning that the current study focused on the journal bearing (with the rotational motion of shaft), whereas in the Abbaspur et al.60 the geometry was sliding bearing with reciprocating motion. in the current study for the first time, the role of two parameters for journal bearings: the elongational viscosity and the first normal stress difference compared.
Most models used to investigate the lubrication problem for a non-Newtonian fluid were simple and linear models. They could not predict the behavior of these fluids precisely. For more reliable results, nonlinear and complex models such as the Giesekus, FENE-P61, and Phan-Thein-Tanner (PTT) model62 are better choices. Viscoelasticity is the property of a material to demonstrates both viscous and elastic properties under the same conditions when it undergoes deformation. Broadly speaking, viscoelasticity is divided into two major fields: linear and nonlinear. Linear viscoelasticity is the field of rheology devoted to the study of viscoelastic materials under very small strain or deformation where the displacement gradients are very small, and a linear relationship between stress and rate of strain for linear materials can describe the flow regime. In principle, the strain has to be sufficiently small so that the material's structure remains unperturbed by the flow history63,64. The linear viscoelastic models have several limitations. For example, they cannot describe strain rate dependence of viscosity or normal stress difference phenomena since they are nonlinear effects. Due to the restriction to infinitesimal deformations, the linear models may be more appropriate for describing viscoelastic solids rather than viscoelastic fluids. Despite the limitations of the linear viscoelastic models and despite being of less interest to the study of flow where the material is usually subject to large deformation, they are very important in the study of viscoelasticity for several reasons65,66.
-
1.
They are used to characterize the behavior of viscoelastic materials at small deformations.
-
2.
They serve as a motivation and starting point for developing nonlinear models since the latter are generally extensions to the linear.
-
3.
They are used for analyzing experimental data obtained in small deformation experiments and for interpreting important viscoelastic phenomena, at least qualitatively.
On the other hand, in reality many of fluids are nonlinear, with large deformations, and with nonlinear response in the presence of such deformations. Nonlinear viscoelastic behavior is usually exhibited when the deformation is large and most of the time when the material changes its properties under . It is worth mentioning that they can describe strain rate dependence of viscosity or normal stress phenomena since these are nonlinear effects. For these reasons, nonlinear viscoelastic mathematical models are needed. Among the nonlinear models, the Giesekus model obtained the most precise outcome due to its generality for the lubrication problem. The Giesekus viscoelastic model is based on a molecular configuration and describes non-linear viscoelastic properties. Also, this model describes many parameters such as linear regions (power-law) of viscosity, the first and second normal stress coefficients, the continuous strain-hardening elongational viscosity with a finite asymptote, the complex viscosity, and start-up curves.
This Study presents an analytical solution for hydrodynamic lubrication of a journal bearing with variable film thickness using Giesekus fluid for the first time. The equations were solved using the perturbation method, and the mobility factor was employed as the perturbation parameter. The effects of the mobility factor, the Weissenberg number, eccentricity of pressure distribution, the first and second normal stress differences, and load capacity are discussed. In this paper, the contribution of each normal force via stress ratio (SR) analysis (ratio of normal stress arising from elongational viscosity to the first normal stress difference) distinguish which can be considered the main innovative aspect of the present study. Also the flow type parameter to distinguish the flow type from kinematic point of view is used. It is a term that can be used to describe certain characteristics or properties related to the flow behavior of fluids, particularly in the field of fluid dynamics and rheology. This parameter provides information about how a fluid flows or behaves under different conditions. Depending on the context, various flow type parameters may be relevant. It helps determine flow regimes, predict flow patterns, and design efficient systems.
The main innovative aspects of the present study can be summarized as follows:
-
A new analytical solution for nonlinear viscoelastic lubrication in journal bearings is presented.
-
The role and contribution of the first normal stress difference and the elongational viscosity on the load capacity are investigated via the stress ratio and flow type parameter.
-
The effects of mobility factor, fluid elasticity, and geometrical properties on journal-bearing lubrication via viscoelastic fluids are studied in detail.
Governing equations and constitutive equation
The governing equations for incompressible viscoelastic fluids are the continuity and momentum.
equations, given by67:
where \(\widetilde{{\varvec{V}}}\) is the velocity vector, \(\widetilde{\rho }\) is the density, \(\widetilde{\mathrm{p}}\) is the pressure, and \(\widetilde{{\varvec{\tau}}}\) is the stress tensor. Here, the superscript ââ~ââ is used to characterize the quantities with dimension, and \(\overrightarrow{\nabla }\) is the gradient operator. The Giesekus model was used for the constitutive equation to model the stress field64:
where \(\alpha\) is the mobility factor, \(\widetilde{\lambda }\) is the relaxation time, \({\widetilde{\eta }}_{0}\) is the viscosity at zero shear rate, and \(\widetilde{{\varvec{D}}}\) is the deformation rate defined as:
The superscript âTâ is the transpose operator, and \(\stackrel{\nabla }{\widetilde{{\varvec{\tau}}}}\) is the upper convected derivative of the stress tensor:
The Giesekus model is based on molecular concepts and represents many of the characteristics of viscoelastic fluids well. In this model, the viscoelastic component of the extra stress tensor is depicted with the parameters \({\widetilde{\eta }}_{0}\), \(\widetilde{\lambda }\), and \(\alpha\). This model has gained prominence because it describes the power-law regions for viscosity and normal-stress coefficients and reasonably explains the elongational viscosity and complex viscosity61. It should be noted that the values of \(\widetilde{\lambda }\) for lubricants used in bearings are in the range of \({10}^{-4}-{10}^{-6}\)45,56 and physically acceptable mobility factor (\(\alpha\)) values are \(0<\alpha <0.5\)68.
The geometry of the journal bearing system used in this study is shown in Fig. 1. A journal bearing consists of an approximately cylindrical body around a rotating shaft, used either to support a radial load or simply as a guide for smooth transmission of torque. It involves a stationary sleeve (or bushing) with a complete 360⦠arc or various arrangements of a partial arc or arcs in a housing structure, as shown in Fig. 1a. The inner surface is commonly lined with a soft bearing material such as lead or tin babbitt, bearing bronze, or a plastic. When motion is initiated, the shaft first rolls up the wall of the sleeve in the direction opposite to rotation due to metal-tometal friction between the steel shaft and the bearing bore.With an adequate lubricant supply, a supporting wedge-shaped film of lubricant is almost immediately formed to lift the journal into its steady-state position. Figure 1b shows the components of pressure projected along the line of centers (radial direction R and tangential direction perpendicular to the line of centers T). As the figure depicts, the circular shaft rotates with an angular velocity \(\widetilde{\omega }\). The radius of the journal is \(\widetilde{R}\) which is smaller than the radius of the bearing surrounding it.
The thickness of the lubricant layer varies according to the angular position θ and can be expressed as follows:
where ε is the relative eccentricity, \(\widetilde{\delta }\) is the radial gap size (the difference between the radius of bearing and journal), and e is the eccentricity. As the figure depicts, the space between the bearings and the journal is very tight and varies between \(\widetilde{H}-\widetilde{\delta }\) and \(\widetilde{H}+\widetilde{\delta }\).
As an assumption, the cavitation is not considered in this study, However, it's important to note that there may be concerns about cavitation or other fluid-related issues in some specialized bearing applications, such as high-speed or high-load scenarios. In such cases, engineers and lubricant manufacturers may take steps to address these concerns, which could include selecting lubricants with specific additives or properties to mitigate the risk of cavitation or other issues related to fluid dynamics. Using updated cavitation models for solving the governing equation can be a new study, and we will start this issu.
The first step to analyze the lubrication problem in the journal bearing is investigating the fluid flow in a Cartesian coordinate system. In order to non-dimensionalize the governing and constitutive equations, the following dimensionless parameters are defined:
where Wi is the Weissenberg number, which expresses the ratio of the elastic force to the viscous force. De called the Deborah number and in order to characterize the intrinsic fluidity of a material or how ââfluid-likeââ the material is, this number has been introduced. After deriving the governing and constitutive equations by applying the appropriate assumptions60, the following non-dimensionalized equations can be derived:
The dimensionless form of the fluid film thickness can be written as
Perturbation solutions
The implemented method to linearize the equations is the perturbation method. The mobility factor, \(\alpha\), was employed as the perturbation parameter because its value is always less than one. Therefore, it is a good choice because we can derive solutions that predict high-order nonlinearity without serious convergence problems. For \(\alpha =0\), the solution is simplified to the UCM lubrication problem, so the high-order terms specify the deviations from the quasilinear to the nonlinear response. Hence, the variables in the Eqs. (8) to (12) expanded as follows:
Reynolds equations
Considering the exhibited geometry, the velocity boundary conditions can be written as
The dimensionless Reynolds equations in a polar coordinate system where \(\widetilde{x}=\widetilde{R}\theta\) and \(d\widetilde{x}=\widetilde{R} d\theta\) can be obtained from zero to second order as follows5,69:
Zero-order Reynolds equation:
The first-order Reynolds equation:
The second-order Reynolds equation:
In order to solve Reynolds equations using the change of variables, the following assumptions proposed by Somerfield5,69 are considered:
The pressure boundary conditions for the present problem are written as:
An unwrapped schematic of the film shape is shown in Fig. 2. According to the figure, \({h}_{min}\) = \(\left(1+\varepsilon cos{\theta }_{m}\right),\) and \({\theta }_{m}\) represents the circumferential angle at which \(d\mathrm{p}/d\theta =0\).
Finally, after solving Reynolds equations, by considering Sommerfeld's assumption and pressure boundary conditions, the final form of Eq. (14) can be expressed as follows:
Due to the large size of the solution of pressure distribution, \({p}_{1}\left(\theta \right)\) and \({p}_{2}\left(\theta \right)\) are presented in the appendix in supplementary material.
Load-carrying capacity
Load bearing capacity is the maximum ability of a journal bearing to take loading before failure occurs. Once the pressure distribution is determined, the load capacity can be calculated. The radial and tangential nondimensional load capacity obtained by the integration of nondimensional pressure along and perpendicular to the line of centers are expressed in Eq. (25):
After integration and simplification, the different orders of load-carrying capacity are written as follows:
Important parameters
The first and second normal stress differences are two important rheological properties of viscoelastic liquids. The shear flow of polymeric solutions is not isotropic due to the alignment of polymeric chains along the streamlines, and it causes an imbalance between the normal forces:
where \({N}_{1}\) and \({N}_{2}\) are the first and the second normal stress differences, and \({\widetilde{\Psi }}_{1}\) and \({\widetilde{\Psi }}_{2}\) are their corresponding coefficients. In Eqs. (27) and (28), the subscript 1 on Ï denotes the shear flow direction, 2 denotes the direction of the variation of velocity profile, and 3 denotes the third right-hand direction. The viscosity and the coefficients of the first and the second normal stress differences are usually known as the viscometric functions. The viscometric functions of the Giesekus model in the simple shear flow can be expressed as follows64:
where
where \(\widetilde{\dot{\gamma }}\) is the shear rate, and \(\widetilde{\kappa }\) is the retardation time constant of the model.
The viscoelastic liquids exhibit a strain-hardening resistance against elongational deformations. The elongational viscosity of the Giesekus model in the steady uniaxial extensional deformation can be obtained from the following relationship:
where \(\dot{\varepsilon }\) is the strain rate, and \({\upeta }_{\mathrm{E}}\) is the extensional viscosity. The combined effects of normal stress differences and extensional viscosity in such flows are rather complex. To compare the role of elongational viscosity and the first normal stress difference, we defined a parameter that is known as the stress ratio, SR:
where subscripts E and S denote the contribution of elongational and shear effects in a mixed flow, respectively.
The flow type parameter, ξ, is another useful dimensionless factor to identify the contribution of different types of deformations in a mixed flow:
where \(\Vert {{\varvec{D}}}_{{\varvec{j}}}\Vert\) and \(\Vert {{\varvec{\Omega}}}_{{\varvec{j}}}\Vert\) are the magnitudes of the rate of deformation and vorticity tensors:
Unlike the stress ratio, ξ is a kinematic parameter and may vary within the range [â 1, 1], in which ξâ=ââ 1 characterizes a solid-like rotational flow, ξâ=â1 characterizes a pure extensional flow, and ξâ=â0 characterizes a simple shear flow.
Validation
Upper convective maxwell validation
The work of Liu and Grecov70 was used to validate the present solution. They employed the UCM model for viscoelastic fluid within a journal bearing. The UCM model is a differential generalization of the Maxwell model for the case of large deformations based on the upper-convected time derivative. According to the UCM the shear viscosity and the first normal stress difference are independent of shear rate and hence the model fails to describe the behavior of most viscoelastic fluids. Furthermore, it predicts a steady-state elongational viscosity that becomes infinite at a finite elongation rate, which is obviously far from physical reality. A major limitation of the UCM is that it do not allow for strain dependency and second normal stress difference. To account for strain dependent viscosity and non-zero second normal stress difference among other phenomena, more sophisticated models such as Giesekus which introduce additional parameters should be considered. This equation, Giesekus model, however, has rarely been used because of the theoretical and experimental complications it introduce.
The model can be written as:
To verify the precision of the perturbation method, a comparison was made with the UCM modelâs results. There are resemblances between the behaviors of the Giesekus model (Eq. 3) and the UCM model. For \(\alpha =0,\) the model is simplified to the isotropic UCM model71. In the study by Kai Liu and Dana Grecov, two parameters are available for validation along the cylinder: pressure distribution and shear stress. Given the relevance of these parameters to the context of journal bearings, the pressure distribution have chosen. This selection is based on the crucial role pressure distribution plays in journal bearings, as it directly influences their load-carrying capacity, friction reduction, lubrication efficiency, stability, wear and tear, and overall performance.
Figure 3 illustrates the comparison of the pressure distribution between Giesekus and UCM models. The results of Liu and Grecov70 are consistent with a perturbation method. Another point to note is that as the Weissenberg number decreases (Fig. 3), the pressure distribution decreases, similar to the work of Liu and Grecov70. However, on closer inspection, it is noticeable that when reducing the Weissenberg number, the validation becomes much more accurate. Therefore, the outputs were taken with lower Weissenberg numbers.
Newtonian lubricant validation
The exact solution of the pressure distribution equation for a journal bearing assuming Newtonian lubricant (\({\mathrm{p}}_{\mathrm{N}}\)) is69:
The case of \(\alpha =0\) and \(\widetilde{\lambda }=0\) in Giesekus model (Eq. 3) is the Newtonian lubrication scenario. The results of the perturbation solution, such as the load capacity of the bearing, have excellent agreement with the exact solution when \(\alpha\) and \(\widetilde{\lambda }\) are equal to zero. Therefore, according to the two validations above, the mobility factor can be a good choice as a perturbation parameter for solving the equation using the perturbation method.
Results and discussion
Pressure distribution
Generally, we should consider 0â<âαâ<â0.5 for realistic properties64. After applying the perturbation theory, considering α as the perturbation parameter for the fourth order, and comparing the results with some works40,70, it has been found that the best range in which the desired answers are correct is 0â<âαâ<â0.3. Figure 4 displays the effects of the eccentricity ratio on the pressure distribution. As the eccentricity ratio increases, the fluid-bearing pressure increases. Due to the continuity equation and Bernoulli's principle, the amount of pressure rises significantly to the apex at approximately \(\theta =180\) because the radial gap size, \(\widetilde{\delta }\), reaches its minimum value in this region. After crossing the utmost point, it experiences a decreasing trend until it finally reaches zero. An issue to notice is that the pressure distribution is symmetric and possesses positive amounts in the convergence region (\(0^\circ <\theta <180^\circ\)) of the journal bearing. Then the fluid film flow enters the diverging area \(180^\circ <\theta <360^\circ\) and yields a negative pressure until the value finally reaches zero at \(\theta =360^\circ\).
Figures 5 and 6 emphasize how the mobility factor and the Weissenberg number alter the pressure distribution. With the increase of the mobility factor and the Weissenberg number, the pressure value decreases before \(\theta =180^\circ\), and an upward trend follows. Evidently, increasing the mobility factor decreases the pressure distribution. It is noteworthy that since the mobility factor in the Giesekus model is in the range of \(0<\alpha <1\), the model behaves as a shear-thinning fluid71. Hence, it is expected that when increasing the value of \(\alpha\), the fluid viscosity (\({\eta }_{0}\)) decreases, and the pressure distribution simultaneously decreases. Another point is that the pressure distribution behavior for all of the parameters is almost identical, but there is a slight difference; in particular, it is obvious that the mobility factor is changing.
Shear stress variations
Figures 7, 8 and 9 show shear stress versus angle variations on the inner cylinder for different parameters. Figure 7 illustrates the influence of the mobility factor on shear stress. Initially, a decrease in shear stress is seen as the angle grows larger until it hits its maximum at \(180^\circ <\theta <210^\circ\). At \(\theta \approx 180^\circ\), specifically relating to low mobility factors, a noticeable difference is identifiable in the range of \(180^\circ <\theta <210^\circ\). Beyond this zone, the behavior of \(\tau\) and its values for all the mobility factors are almost the same. Increasing the mobility factor in the Giesekus model leads to an increase in shear stress because it enhances the material's elastic response and its ability to resist deformation, store energy, and generate shear stress as a result.
Figure 8 displays the Weissenberg numberâs effects on shear stress. As it can be seen, the shear stress profile shows roughly symmetric behavior for various Weissenberg numbers. At \(\theta =120^\circ\), there is a gradual decrease for \(\tau\). Next, it falls dramatically until it peaks at \(180^\circ <\theta <210^\circ\). In the end, it follows a symmetric distribution. Furthermore, the shear stress profile changes substantially when the Weissenberg number is further increased. The material can store deformation energy and release it over time. As the material experiences shear deformation, the elastic component of the material resists the change in shape. This resistance to deformation results in the generation of shear stress within the material. In materials with a high Weissenberg number, the elastic response causes a delay in reaching a steady-state deformation. The material "remembers" the deformation, and the shear stress continues to increase until the material fully responds to the applied deformation. The elastic component of the material stores energy during deformation, and this stored energy contributes to an increase in shear stress until the material reaches equilibrium.
Figure 9 illustrates the shear stress variations concerning the eccentricity ratios. The results for \(\varepsilon <0.5\) are interpreted roughly as in Fig. 8. At \(\varepsilon <0.5\), the shear stress variations are dissimilar to other characteristics. It first moves upward, then downward. Additionally, it is intriguing that the maximum and the minimum values occur in a tight range of angle change, and there is a considerable shift in shear stress in the range of \(170^\circ <\theta <190^\circ\). This change causes an asymmetric shear stress distribution at \(\varepsilon <0.5\). As the eccentricity ratio increases, the journal moves away from the central position within the bearing. This means that the lubricating oil film between the journal and the bearing surface becomes thinner on one side (in the current study it happen around θâ=â180) and thicker on the other side due to the offset. When there is a variation in shear rate, it results in increased shear stress within the lubricant. In areas where the lubricant film is thinner, the shear rate is higher, leading to higher shear stresses. In contrast, where the film is thicker, the shear rate is lower, resulting in lower shear stresses. Higher eccentricity ratios can create higher shear stress regions in the bearing, which are essential for load-carrying capacity but must be managed to prevent excessive wear and overheating.
Apparent First normal stress difference
One of the most relevant effects of elasticity on a flow is related to the apparent first normal stress difference. In a journal bearing, the flow type is mixed (elongation and shear), so the total amount of normal stress difference resulting from shear and elongation is considered as the apparent first normal stress difference. Figures 10, 11 and 12 demonstrate the profiles of the apparent first normal stress difference on the inner cylinder for various relaxation times. The apparent first normal stress difference for polymeric fluids is practically always positive. Generally, this behavior can be noticed in all the apparent first normal stress difference figures. For a first approximation, this indicates that polymeric fluids exhibit extra tension in addition to the shear stresses along the streamlines in the â1â direction64.
For all three influential parameters (mobility factor, Weissenberg number, and eccentricity ratio), the apparent first normal stress difference first has a descending trend followed by an increasing trend. At first glance, augmenting each of the three parameters magnifies the amount of apparent first normal stress difference. Moreover, the apparent first normal stress difference profiles are almost flat for low relaxation times. Changes to the stress profiles due to changes in influential parameters are clearly nonlinear and display significant viscoelastic effects.
According to Figs. 10 and 12, the trend of change is similar for different values of the mobility factor and Weissenberg number. When increasing these quantities, the distance between extremum points decreases, and the highest value of the apparent first normal stress difference also occurs at \(\theta \approx 180^\circ\). Figure 12 exhibits that for \(\varepsilon \le 0.4\), the trend is similar, and no notable difference is seen. However, after passing this value, the variations of the apparent first normal stress difference change substantially (especially for \(180^\circ <\theta <210^\circ\)).
Load capacity of bearing
Figures 13 and 14 present the evolution of load capacity on the inner cylinder for various mobility factors and the Weissenberg number in terms of the eccentricity ratio. Changes for both parameters are initially linear (\(\varepsilon \le 0.4\)) and have predictable trends. However, changes in the mobility factor and the Weissenberg number do not impact the bearing capacity distribution in this area. For \(\varepsilon >0.4\), the trend of changes is exponential, and an effect of changes in parameters is observable, especially for the Weissenberg number. The results point out the dominant influence of elastic properties with respect to viscosity effects, as mentioned before.
Consequently, this makes the journal bearing work more efficiently and safely since high load capacity maintains a minimum thickness of lubricant that is necessary between two relatively moving surfaces. This is a beneficial impact of viscoelasticity on journal-bearing performance. An increase in eccentricity enhances the load capacity. This is consistent with the experimental observations of higher load capacity load for higher eccentricity10. According to the Fig. 13, as the mobility factor increases, it implies that the fluid can respond more quickly to changes in shear deformation. In a journal bearing, where there is relative motion between the journal and the bearing surface, the lubricating fluid experiences shear deformation due to this motion. A faster response to shear means that the lubricating fluid can build up a thicker lubricating film more quickly between the journal and the bearing surface. With a thicker lubricating film and reduced friction, the bearing can support higher loads without excessive heat generation or wear. This can lead to an increase in load-carrying capacity. Figure 14 shows the elastic effects can enhance the load carrying capacity in certain situations. When a viscoelastic fluid is subjected to deformation or shear, the polymer chains in the fluid can stretch and store elastic energy. This elastic energy can help support and distribute the load, increasing the load carrying capacity which cannot be seen in Newtonian lubricants.
Stress ratio
As mentioned before, the stress ratio (SR) is defined as the ratio of stress resulting from the elongational viscosity to the first normal stress difference. A diagram of SR along the angular coordinate θ is shown in Figs. 15, 16, 17 and 18 according to the variation of \(\upvarepsilon\), Wi, α, and De. As the figures illustrate, two areas are very important: \(120^\circ <\uptheta <180^\circ\) and \(180^\circ <\uptheta <240^\circ\). Figure 15 depicts how the stress ratio, SR, changes with the alteration of the eccentricity ratios. Decreasing \(\upvarepsilon\) leads the fluid flow to a narrowing gap; therefore, the surface areas in the zone are under remarkable stress.
According to Fig. 16, increasing the Weissenberg number results in the decrease of SR. Consequently, it reaches zero at Wiâ=â1. The impact of mobility factors and Deborah numbers on SR is shown in Figs. 17 and 18. Just under \(\uptheta =240^\circ\). Increasing \(\mathrm{\alpha }\) and particularly De have significant influences on SR. If the Deborah number increase, polymer molecules that are distorted by the flow will not have time to relax during the time scale of the process or experiment. Therefore, according to Fig. 18, at about \(\uptheta =240^\circ\), increasing the Deborah number has resulted in an overwhelming growth in the value of SR. Consequently, in this zone, the contribution of normal stress differences and extensional viscosity is higher than in other areas. According to Eq. (34), there is a direct correlation between SR and extensional viscosity \(({\upeta }_{\mathrm{E}})\). However, overall, for \(\kappa =0\) and \(\mathrm{\alpha }\)â 0, the following asymptotic formula holds64:
As a result, the magnitude of the elongational viscosity, which characterizes a fluid's resistance to elongational deformation, is higher than the first normal stress difference in the present problem. Moreover, because of their long chain molecules, polymers exhibit stiff resistance to any elongational deformation; therefore, the elongational viscosity of a polymer is generally quite high.
Flow type parameter
Figure 19 presents the impacts of the mobility factor, eccentricity ratio, and Weissenberg number on flow type parameters for the Giesekus model at different axial positions. Overall, it is clear that in all contours, the general type of flow is simple shear flow (ξâ=â0). Nevertheless, there are two symmetrical points around \(\uptheta =180^\circ\) at which we have pure extensional flow (ξâ=â1). To be more precise, regarding Fig. 19a, increasing the mobility factor shows no remarkable change in the fluid flow, and approximately, the positions of two points at which we have extensional flow are stable at almost \(\uptheta =130^\circ \mathrm{ and }220^\circ\). On the other hand, as Fig. 19b depicts, increasing the figure of eccentricity ratio leads to the fluid within the journal bearing experiencing a pure extensional flow close to \(\uptheta =180^\circ\). Surprisingly, an increase in the Weissenberg number causes a new location to appear around \(\uptheta =90^\circ\) with another extensional flow area. The flow type parameter, is crucial in determining the lubrication regime and, consequently, the performance and reliability of journal bearings. Proper design and lubricant selection are essential to ensure that the bearing operates within the desired hydrodynamic regime to minimize wear and maximize load-carrying capacity.
Conclusions
The present study looked at the effect of viscoelastic fluid lubrication on the performance of journal bearings by employing the Giesekus model. The major variables of the Giesekus fluid flow in a journal bearing were developed by applying the perturbation theory.
The main conclusions drawn from the present study are:
-
As the eccentricity ratio increases, the pressure distribution of the fluid increases, while with the increase of the mobility factor and the Weissenberg number, the pressure distribution decreases.
-
Mobility factor variations have no significant effects on the bearing pressure distribution.
-
The first normal stress difference for variation of all three effective parameters (the mobility factor, the Weissenberg number, and the eccentricity ratio) initially shows a downward trend and then an upward trend while by increasing all three parameters, the first normal stress difference also increases.
-
The second normal stress difference alternations for the mobility factor, the Weissenberg number, and the eccentricity ratio are similar. For large values of independent variants, the dependent quantity fluctuates more evidently. The extreme amounts are near the narrowest region of a journal bearing.
-
The bearing load capacity can be divided into two parts: a linear region where the mobility factor and Weissenberg numbers have little effect and the bearing load distribution increases linearly, and a second area where load capacity changes are exponential.
-
When the Weissenberg number is zero, it implies that elastic effects are negligible, and the fluid can be treated as a simpler, Newtonian fluid. As a result, you can observe that increasing the Weissenberg number significantly affects bearing load capacity when using a viscoelastic fluid compared to a Newtonian fluid.
-
When a viscoelastic fluid experiences deformation or shear, the polymer chain within the fluid have the capacity to elongate and retain elastic energy. This elastic energy can help support and distribute the load, increasing the load carrying capacity and a feature not observed in conventional Newtonian lubricants.
-
Increasing the mobility factor in the Giesekus model leads to an increase in shear stress than Newtonian lubricants because it enhances the material's elastic response and its ability to resist deformation, store energy, and generate shear stress as a result.
-
Greater eccentricity ratios can lead to the formation of elevated shear stress areas within the bearing. While these regions are vital for supporting loads, it's imperative to effectively control them to avoid excessive wear and overheating.
-
Polymers exhibit a stiff resistance to any elongational deformation due to their long chain molecules; hence, the elongational viscosity of a polymer is naturally high. Therefore, this subject can be monitored particularly where the flow has higher viscoelasticity, such as \(\theta =240^\circ\) and Deâ=â1.
-
Generally, the fluid flow is similar to a simple shear flow except around \(\theta =180^\circ ,\) where we have a pure extensional flow.
Abbreviations
- \(\widetilde{D}\) :
-
Deformation rate
- \(\mathrm{De}\) :
-
Deborah number
- e :
-
Eccentricity
- \(\widetilde{H}\) :
-
The thickness of the lubricant layer
- \({N}_{1}\) :
-
First normal stress difference (Pa)
- \({N}_{2}\) :
-
Second normal stress difference (Pa)
- \(\widetilde{\mathrm{p}}\) :
-
Pressure (Pa)
- \(\widetilde{R}\) :
-
Radius of the journal (m)
- \(\mathrm{SR}\) :
-
Stress ratio
- \(\widetilde{U}\) :
-
Surface sliding speed parameter (m/s)
- \(\widetilde{V}\) :
-
Vector of velocity (m/s)
- \(u,v\) :
-
Velocities along x and y directions, respectively
- \(W\) :
-
Load-carrying capacity
- Wi:
-
Weissenberg number
- \(x\),\(y\) :
-
Dimensionless coordinate
- \(\alpha\) :
-
Mobility factor
- \(\widetilde{\dot{\gamma }}\) :
-
Shear rate (s-1)
- \(\widetilde{\delta }\) :
-
Radial gap size (m)
- É:
-
Relative eccentricity
- \(\widetilde{\dot{\upvarepsilon }}\) :
-
Strain rate (s-1)
- \({\widetilde{\eta }}_{0}\) :
-
Viscosity (Pa.s)
- \({\widetilde{\upeta }}_{\mathrm{E}}\) :
-
Extensional viscosity (Pa.s)
- θ :
-
Angular position
- \(\widetilde{\kappa }\) :
-
Retardation time (s)
- \(\widetilde{\lambda }\) :
-
Relaxation time (s)
- ξ:
-
Flow type parameter
- \(\widetilde{\rho }\) :
-
Density (kg/m3)
- \(\widetilde{\tau }\) :
-
Stress tensor (Pa)
- \(\stackrel{\nabla }{\widetilde{\tau }}\) :
-
Upper-convected derivative (Pa)
- \(\sigma\),\(\tau\) :
-
Dimensionless normal stress components
- \(\gamma\) :
-
Dimensionless shear stress
- \(\widetilde{\Psi }\) :
-
The coefficient of normal stress difference
- \(\widetilde{\omega }\) :
-
Angular velocity
- \({\Omega }_{j}\) :
-
Rate of vorticity tensor
- \(\overrightarrow{\nabla }\) :
-
Gradient operator
References
Tower, B. First report on friction experiments. Proceed. Institut. Mech. Eng. 34, 632â659 (1883).
Tower, B. Second report on friction experiments. Proc. Inst. Mech. Engrs. 36, 38â70 (1985).
Tower, B. Third report on friction experiments 173â205 (Proc. Inst. Mech. Engrs, 1888).
Reynolds, O. IV. On the theory of lubrication and its application to Mr. Beauchamp towerâs experiments, including an experimental determination of the viscosity of olive oil. Philosophical Trans. R. Soc. London 177, 157â234 (1886).
Philipzik, W. Zur hydrodynamischen Theorie der Schmiermittelreibung. ZAMM-J. Appl. Math. Mech. /Zeitschrift für Angewandte Mathematik und Mechanik 36, 51â60 (1956).
Zhang, X., Zhang, Y., Li, J. & Jin, Z. Scientific contributions of Duncan Dowson to Tribological Studies of Artificial Joints. Proceed. Institut. Mech. Eng. J J. Eng. Tribol. 236, 1554â1569 (2022).
Raimondi, A. A. & Boyd, J. A solution for the finite journal bearing and its application to analysis and design: I. A S L E Trans. 1, 159â174. https://doi.org/10.1080/05698195808972328 (1958).
Raimondi, A. & Boyd, J. A solution for the finite journal bearing and its application to analysis and design: II. ASLE Transactions 1, 175â193 (1958).
Raimondi, A. & Boyd, J. A solution for the finite journal bearing and its application to analysis and design: III. ASLE Trans. 1, 194â209 (1958).
Raimondi, A. A numerical solution for the gas lubricated full journal bearing of finite length. ASLE Trans. 4, 131â155 (1961).
Reddi, M. M. Finite-element solution of the incompressible lubrication problem. J. Lubricat. Technol. 91, 524â533 (1969).
Reddi, M. & Chu, T. Finite element solution of the steady-state compressible lubrication problem. J. Lubricat. Technol. 92, 495â502 (1970).
Malik, M. Theoretical considerations of molecular mean free path influenced slip in self-acting gas-lubricated plain journal bearings. Proceed. Institut. Mech. Eng. C J. Mech. Eng. Sci. 198, 25â31 (1984).
Dinariev, O. Y. Plane hydrodynamic problem for viscoelastic lubrication. J. Eng. Phys. Thermophys. 74, 771â776 (2001).
Urreta, H. et al. Hydrodynamic bearing lubricated with magnetic fluids. J. Intell. Mater. Syst. Struct. 21, 1491â1499 (2010).
Sfyris, D. & Chasalevris, A. An exact analytical solution of the Reynolds equation for the finite journal bearing lubrication. Tribol. Int. 55, 46â58 (2012).
Chasalevris, A. & Sfyris, D. Evaluation of the finite journal bearing characteristics, using the exact analytical solution of the Reynolds equation. Tribol. Int. 57, 216â234 (2013).
Rao, T., Rani, A., Nagarajan, T. & Hashim, F. Analysis of slider and journal bearing using partially textured slip surface. Tribol. Int. 56, 121â128 (2012).
Vignolo, G. G., Barilá, D. O. & Quinzani, L. M. Approximate analytical solution to Reynolds equation for finite length journal bearings. Tribol. Int. 44, 1089â1099 (2011).
Gong, R.-Z., Li, D.-Y., Wang, H.-J., Han, L. & Qin, D.-Q. Analytical solution of Reynolds equation under dynamic conditions. Proceed. Institut. Mech. Eng. J J. Eng. Tribol. 230, 416â427 (2016).
Jang, J. Y. & Khonsari, M. M. On the characteristics of misaligned journal bearings. Lubricants 3, 27â53 (2015).
Barnes, H. A., Hutton, J. F. & Walters, K. An introduction to rheology. Vol. 3 (Elsevier, 1989).
Bates, T. W., Williamson, B., Spearot, J. A. & Murphy, C. K. A correlation between engine oil rheology and oil film thickness in engine journal bearings. SAE transactions, 94â112 (1986).
Tayal, S., Sinhasan, R. & Singh, D. Analysis of hydrodynamic journal bearings having non-newtonian lubricants (prandtl model) by a finite element method. J. Mech. Eng. Sci. 23, 63â68 (1981).
Horowitz, H. & Steidler, F. Calculated performance of non-Newtonian lubricants in finite width journal bearings. ASLE TRANSACTIONS 4, 275â281 (1961).
Tanner, R. A short-bearing solution for pressure distribution in a non-Newtonian lubricant. (1964).
Wada, S. & Hayashi, H. Hydrodynamic lubrication of journal bearings by pseudo-plastic lubricants: Part 1, theoretical studies. Bulletin of JSME 14, 268â278 (1971).
Wada, S., Hayashi, H. & Haga, K. Behavior of a Bingham solid in hydrodynamic lubrication: Part 1, general theory. Bulletin JSME 16, 422â431 (1973).
Swamy, S., Prabhu, B. & Rao, B. Calculated load capacity of non-Newtonian lubricants in finite width journal bearings. Wear 31, 277â285 (1975).
Swamy, S., Prabhu, B. & Rao, B. Stiffness and damping characteristics of finite width journal bearings with a non-Newtonian film and their application to instability prediction. Wear 32, 379â390 (1975).
Raghunandana, K., Majumdar, B. C. & Maiti, R. Stability of flexibly supported oil journal bearings using non-Newtonian lubricants linear perturbation analysis. J. Trib. 123, 651â654 (2001).
Das, S., Guha, S. & Chattopadhyay, A. On the steady-state performance of misaligned hydrodynamic journal bearings lubricated with micropolar fluids. Tribol. Int. 35, 201â210 (2002).
Das, S., Guha, S. & Chattopadhyay, A. Linear stability analysis of hydrodynamic journal bearings under micropolar lubrication. Tribol. Int. 38, 500â507 (2005).
Abdel-Rahman, G. M. Flow of a non-Newtonian power law through a conical bearing in an applied magnetic field. Appl. Math. Computat. 159, 237â246 (2004).
Elsharkawy, A. A. Effects of lubricant additives on the performance of hydrodynamically lubricated journal bearings. Tribol. Lett. 18, 63â73 (2005).
Tian, G., Wang, M., Wang, X. & Jin, G. Flow between eccentric cylinders: a shear-extensional controllable flow. Korea-Australia Rheol. J. 28, 139â148 (2016).
Sakim, A., Nabhani, M. & El Khlifi, M. Non-Newtonian effects on porous elastic journal bearings. Tribol. Int. 120, 23â33 (2018).
Chetti, B. & Zouggar, H. Steady-state performance of a circular journal bearing lubricated with a non-Newtonian fluid considering the elastic deformation of the liner. Proceed. Institut. Mech. Eng. J J. Eng. Tribol. 233, 1389â1396 (2019).
Li, B. et al. Thermohydrodynamic lubrication analysis of misaligned journal bearing considering the axial movement of journal. Tribol. Int. 135, 397â407 (2019).
Gwynllyw, D. R. & Phillips, T. N. The influence of Oldroyd-B and PTT lubricants on moving journal bearing systems. J. Non-newtonian Fluid Mech. 150, 196â210 (2008).
Gertzos, K., Nikolakopoulos, P. & Papadopoulos, C. CFD analysis of journal bearing hydrodynamic lubrication by Bingham lubricant. Tribol. Int. 41, 1190â1204 (2008).
Lin, J. R. Weakly nonlinear bifurcation analysis of a short journal bearing lubricated with non-Newtonian couple stress fluids. J. Chinese Institute Eng. 31, 721â727 (2008).
Wierzcholski, K. The viscoelastic lubrication problem of micro-bearing. Tribologia 3, 231â240 (2010).
Guemmadi, M. & Ouibrahim, A. in Key Engineering Materials. 64â69 (Trans Tech Publ).
Tichy, J. Non-Newtonian lubrication with the convected Maxwell model. (1996).
Huang, P., Li, Z.-H., Meng, Y.-G. & Wen, S.-Z. Study on thin film lubrication with second-order fluid. J. Tribol. 124, 547â552. https://doi.org/10.1115/1.1467636 (2002).
Akyildiz, F. T. & Bellout, H. Viscoelastic lubrication with phan-thein-tanner fluid (ptt). J. Trib. 126, 288â291 (2004).
Kumar, A. & Sharma, S. C. Optimal parameters of grooved conical hybrid journal bearing with shear thinning and piezo-viscous lubricant behavior. J. Tribol. 141, 071702 (2019).
Chetti, B., Hemis, M., Tahar, O. & Smara, M. Combined effects of elastic deformation and piezo-viscous dependency on the performance of a journal bearing operating with a non-Newtonian fluid. Proceed. Institut. Mech. Eng. J J. Eng. Tribol. 236, 2457â2467 (2022).
Mokhtari Sharghi, M., Norouzi, M., Akbarzadeh, P. & Abbaspur, A. An investigation on nonlinear viscoelastic lubrication using FENE-P constitutive equation. J. Brazilian Soc. Mech. Sci. Eng. 44, 196 (2022).
Ahmed, H. & Biancofiore, L. A new approach for modeling viscoelastic thin film lubrication. J. Non-Newtonian Fluid Mech. 292, 104524 (2021).
Ahmed, H. & Biancofiore, L. Modeling polymeric lubricants with non-linear stress constitutive relations. Available at SSRN 4450394
Soni, S. Performance of finite bearing under the combined influence of turbulent and non-Newtonian lubrication. Proceed. Institut. Mech. Eng. J J. Eng. Tribol. 235, 2350â2363 (2021).
Agrawal, N. & Sharma, S. C. Micro-grooved hybrid spherical thrust bearing with Non-Newtonian lubricant behaviour. Int. J. Mech. Sci. 240, 107940 (2023).
Hashemabadi, S. & Mirnajafzadeh, S. Analytical solution of simplified Phan-Thien-Tanner fluid between nearly. J Appl Sci 7, 1271â1278 (2007).
Li, X. K., Luo, Y., Qi, Y. & Zhang, R. On non-Newtonian lubrication with the upper convected Maxwell model. Appl. Math. Modell. 35, 2309â2323 (2011).
Nessil, A., Larbi, S., Belhaneche, H. & Malki, M. Journal bearings lubrication aspect analysis using non-Newtonian fluids. Advances in Tribology 2013 (2013).
Li, X. K. Non-newtonian lubrication with the phan-thienâtanner model. J. Eng. Math. 87, 1â17 (2014).
Soni, S. & Vakharia, D. P. Static analysis of finite hydrodynamic journal bearing in turbulent regime with non-Newtonian lubricant. Tribol. Online 10, 246â261 (2015).
Abbaspur, A., Norouzi, M., Akbarzadeh, P. & Vaziri, S. A. Analysis of nonlinear viscoelastic lubrication using Giesekus constitutive equation. Proceed. Institut. Mech. Eng. J J. Eng. Tribol. 235, 1124â1138 (2021).
Bird, R. B. & Wiest, J. M. Constitutive equations for polymeric liquids. Annual Rev. Fluid Mech. 27, 169â193 (1995).
Thien, N. P. & Tanner, R. I. A new constitutive equation derived from network theory. J. Non-Newtonian Fluid Mech. 2, 353â365 (1977).
Carreau, P. J., De Kee, D. C. & Chhabra, R. P. Rheology of polymeric systems: principles and applications. (Carl Hanser Verlag GmbH Co KG, 2021).
Bird, R. B., Armstrong, R. C. & Hassager, O. Dynamics of polymeric liquids. Vol. 1: Fluid mechanics. (1987).
Tuna, N. Finlayson. J. Rheol. 2879, 93 (1984).
Zhou, C. & Kumar, S. Thermal instabilities in melt spinning of viscoelastic fibers. J. Non-newtonian Fluid Mech. 165, 879â891 (2010).
White, F. M. & Majdalani, J. Viscous fluid flow. Vol. 3 (McGraw-Hill New York, 2006).
Yoo, J. & Choi, H. C. On the steady simple shear flows of the one-mode Giesekus fluid. Rheologica acta 28, 13â24 (1989).
Khonsari, M. M. & Booser, E. R. Applied tribology: bearing design and lubrication. (John Wiley & Sons, 2017).
Liu, K. & Grecov, D. Rheological and flow modelling of viscoelastic fluids between eccentric cylinders. Appl. Math. Modell. 35, 1603â1615 (2011).
Cherizol, R., Sain, M. & Tjong, J. Review of non-Newtonian mathematical models for rheological characteristics of viscoelastic composites. Green Sustain. Chem. 5, 6 (2015).
Acknowledgements
This work was supported by Brain Pool Program through the National Research Foundation of Korea (NRF), which is funded by Ministry of Science and ICT (NRF-2022H1D3A2A01091637). This work was also supported by the National Research Foundation of Korea (NRF) grant, which is funded by the Korean government (MSIT) (No. 2020R1A5A8018822, No. 2021R1C1C2009287). The authors also appreciate the financial support from Shahrood University of Technology.
Author information
Authors and Affiliations
Contributions
Ali Abbaspur: Implementation of analytical solution, presenting the results and physical dissuasion, and writing the paper. Mahmood Norouzi: Supervisor, Presenting the idea of the problem, Presenting the method of the solution, Conducting the research, scientific discussion, and editing. Pooria Akbarzadeh: Advisor, Conducting the research, scientific discussion, and editing. Seyyed Amirreza Vaziri: Editing, presenting some parts of the results. Melika Mokhtari Sharghi: Editing, presenting some parts of the results. Kyung Chun Kim: Scientific discussion and interpretation of the results, modification, and English editing. Mirae Kim: Scientific discussion and interpretation of the results, modification, and English editing.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher's note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Abbaspur, A., Norouzi, M., Akbarzadeh, P. et al. An analytical study on nonlinear viscoelastic lubrication in journal bearings. Sci Rep 13, 16836 (2023). https://doi.org/10.1038/s41598-023-43712-8
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-023-43712-8