# A multigrid method for EHL line contact problem with Grease as lubricant

## B. Vishwanath Awati1, Shankar Naik2

1, 2Department of Mathematics, Rani Channamma University, Belagavi, 591156, India

1Corresponding author

Mathematical Models in Engineering, Vol. 3, Issue 2, 2017, p. 126-134. https://doi.org/10.21595/mme.2017.19498
Received 4 October 2017; accepted 4 December 2017; published 31 December 2017

Copyright © 2017 JVE International Ltd. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Views 119
Abstract.

The paper presents an isothermal Elastohydrodynamic lubrication (EHL) of line contact problem with grease as a lubricant, described by Herschel-Bulkley model. The simulation of the governing problem is analyzed using Multigrid method with full approximation scheme (FAS) for varying load, speed and rheological index $n$ Herschel-Bulkley model. The EHL problem considered comprises of Reynolds equation and film thickness understudy for different parameters of interest. The minimum film thickness and pressure spike decrease with a decrease in rheological index $n$. The present paper traces out the comparison of minimum film thickness with the standard results of Dowson and Higginson. The results obtained are comparable and are presented in terms of graphs and tables.

Keywords: elastohydrodynamic lubrication, line contact, grease lubricant, multigrid method, FAS.

#### 1. Introduction

EHL is the study of hydrodynamic lubrication in rolling/sliding contacts when the elastic deformation of the solids is significant. EHL is usually found in many machines parts, for example, gear teeth contacts, cams and rolling element bearings during variations in load etc. EHL problems constitute the following characteristics such as high pressure, deformation of the contacted surfaces and shear rates. The mathematical model of EHL problem depends on a proper fluid rheology model for the lubricant. Many researchers have utilized deliberate methodologies for the numerical solutions of EHL line and point contact problems [1-4]. As most of the roller bearings are lubricated by grease, therefore the EHL equations are solved using grease as a lubricant in order to determine its properties between contacting surfaces. A complete review of EHL theory has been discussed by Lugt et al. .

Grease is commonly used as a lubricant in the application of bearing surfaces, it is a choice of lubricant for 80-90 % of bearings. The essential use of grease in bearing components is to ensure long life and low erosion. The main advantage of using grease inherits sealing action; low friction (provided no leakage). The disadvantage of grease is that it has short period compared to other lubricants. Kauzlarich and Greenwood  determined experimentally the flow curves for grease which are found to be well correlated with the Herschel-Bulkley model. Wen and Ying  examined a hypothetical and test investigation of EHL with thixotropy of grease. Cheng  analyzed the influence of pressure and film thickness dispersion in grease EHL, based on its theoretical and numerical analysis in line contact problem using Herschel-Bulkley model for the rheological characteristics of the grease. The effects of thickener structure on film thickness and traction force of lubricating greases in EHL point contact were discussed by Kaneta et al.  using optical interferometry technique. Cann and Lubrecht  presented the bearing execution limits with grease lubrication, the interaction of bearing design, working conditions and grease properties. Many aspects of grease lubrication have been discussed by Lugt . Cen et al.  examined the film thickness of oil greased up contacts at low speeds. The actual operating conditions with an assortment of greases for the accurate prediction of film thickness is of immense importance to ensure controlled separation of contacting surfaces, leading to greater life extent of the mechanical components. In this article we consider an isothermal EHL problem with grease as a lubricant for its film formation, adapting multigrid FAS for its solution.

#### 2. Mathematical formulation

The Herschel-Bulkley model is one of the most important models describing the non-Newtonian rheology of the grease as a lubricant. The model is generally written as $\tau ={\tau }_{0}+k{\stackrel{˙}{\gamma }}^{n}$, where $\tau$ is the shear stress, ${\tau }_{0}$ the yield stress, $\stackrel{˙}{\gamma }$ the shear rate, $n$ the rheological index and $k$is the plastic viscosity. The rheological characteristics of grease depend on the index $n\le$ 1, ${\tau }_{0}$ and $k$ are functions of pressure and temperature. If $\tau >{\tau }_{0}$ the Herschel Bulkley model behaves as a fluid, otherwise it behaves as a solid. If $n=$ 1 and ${\tau }_{0}=$ 0, this model reduces to the Newtonian fluid. For an incompressible, steady and isothermal flow, the modified Reynolds equation governing grease EHL in its dimensionless form is :

(1)
$\frac{d}{dX}\left(\xi {\left(\frac{dP}{dX}\right)}^{\frac{1}{n}}\right)-\frac{dH}{dX}=0,$

where $\xi =\lambda {H}^{2+1/n}/{\eta }^{1/n}$ and $\lambda ={p}_{H}^{1/n}{b}^{2+1/n}/2U\left(2+1/n\right){R}^{\left(1+1/n\right)}{2}^{1/n}{\eta }_{0}^{1/n}$. The non-dimensional parameters for magnitude analysis of Reynolds equation are:

(2)

The film thickness at any point on the contacted surface is:

(3)
$H\left(X\right)={H}_{0}+\frac{{X}^{2}}{2}-\frac{1}{\pi }\underset{{X}_{in}}{\overset{{X}_{out}}{\int }}\mathrm{l}\mathrm{n}\left|X-{X}^{\text{'}}\right|P\left({X}^{\text{'}}\right)d{X}^{\text{'}},$

where ${H}_{0}$ is the central film thickness, ${X}^{2}/2$ defines undeformed contact shape and the integral term represents the elastic deformation of the contact. The non-dimensional viscosity pressure relation given by Roelands :

(4)
$\stackrel{-}{\eta }=\mathrm{e}\mathrm{x}\mathrm{p}\left\{\left[\mathrm{l}\mathrm{n}\left({\eta }_{0}\right)+9.67\right]\left[-1+{\left(1+5.1×1{0}^{-9}P{p}_{H}\right)}^{z}\right]\right\}.$

The equation for applied load in dimensionless form is:

(5)
$W=\underset{{X}_{in}}{\overset{{X}_{out}}{\int }}Pdx=\frac{\pi }{2}.$

The corresponding boundary conditions are

1) Inlet boundary condition at: $P=$ 0 and.

2) Outlet boundary condition at: $X={X}_{out}$, $P=dP/dX=$ 0.

#### 3. Discretization of Reynolds and film thickness equations

The Eq. (1) is discretized using second order finite difference with grid size of $N=$ 32, 64, 128 and the domain of interest is [${X}_{in}$, ${X}_{out}$] = [–2, 1.5]. The discretized form of Reynolds equation is:

(6)
$\frac{{\xi }_{i+1/2}{\left({P}_{i+1}-{P}_{i}\right)}^{1/n}-{\xi }_{i-1/2}{\left({P}_{i}-{P}_{i-1}\right)}^{1/n}}{\delta {X}^{1+1/n}}-\frac{{H}_{i}-{H}_{i-1}}{\delta X}=0,$

where $\delta X={X}_{i}-{X}_{i-1}$. The discretized form of film thickness Eq. (3) is:

(7)
${H}_{i}={H}_{0}+\frac{{X}_{i}^{2}}{2}-\frac{1}{\pi }\sum _{j=1}^{n}{K}_{ij}{P}_{j},$

where:

${K}_{ij}=-\left(i-j+\frac{1}{2}\right)\delta X\left[\mathrm{l}\mathrm{n}\left(\left|i-j+\frac{1}{2}\right|\delta X\right)-1\right]+\left(i-j+\frac{1}{2}\right)\delta X\left[\mathrm{l}\mathrm{n}\left(\left|i-j-\frac{1}{2}\right|\delta X\right)-1\right],$

for $i=$ 0, 1, 2, ..., $n$, $j=$ 0, 1, 2, ..., $n$. The force balance equation in discrete form is:

(8)
$\delta X\sum _{i=1}^{n}\frac{\left({P}_{i}+{P}_{i+1}\right)}{2}=\frac{\pi }{2}.$

The boundary condition in the inlet $P\left({X}_{in}\right)=$ 0 and the boundary in the outlet can be obtained by setting negative pressure equal to zero, the discrete form is:

(9)

The convergence criteria for pressure relation is:

(10)
$\frac{\sum \left|{P}_{i}^{k+1}-{P}_{i}^{k}\right|}{\sum {P}_{i}^{k+1}}\le \epsilon ,$

where $\epsilon$ is the error.

#### 4. Multigrid method

Multigrid method plays an important role in solving nonlinear mathematical models. Many researchers have developed the generalized ideas in various disciplines underlying multigrid method for the solution of EHL problems. Lubrechet et al.  discussed the multigrid method to solve the EHL problem for both line and point contacts which is a significant breakthrough. The efficiency and robustness are improved by eliminating, both high and low-frequency error components by a series of grid with different mesh levels. Venner [16, 17] and Nurgat et al.  have introduced new relaxation methods for which the technique was greatly improved. The concept of Multigrid method has been discussed in detail by Briggs et al. . For the solution of a nonlinear system of equations, we use multigrid method with FAS scheme. FAS involve an inner and an outer iteration; the outer iteration is the FAS correction scheme, while the inner iteration is usually a standard relaxation method such as Gauss-Seidel iteration scheme.

#### 4.1. Algorithm (FAS)

The system of nonlinear equations can be written as:

(11)
${L}^{h}\left({u}^{h}\right)={f}^{h},$

where $L$ is a nonlinear operator, $u$ is the exact solution, $f$ is a right-hand side function and $h$ is the mesh size of the uniform finer grid ${\mathrm{\Omega }}^{h}$. Let $v$ be an approximation to the exact solution $u$ then:

For simplicity, we consider two grid correction schemes, as follows.

• Applying few number of Gauss-Seidel iteration method on fine grid ${\mathrm{\Omega }}^{h}$ with initial ${v}^{h}$.

• Computation of the restricted residual and restricted current approximation obtained from step 1 (finer to coarser):

• Solve the problem on coarser grid ${\mathrm{\Omega }}^{2h}$:

${L}^{2h}\left({u}^{2h}\right)={L}^{2h}\left({v}^{2h}\right)+{r}^{2h}.$

• Compute the error approximation of the coarser grid:

${e}^{2h}={u}^{2h}-{v}^{2h}.$

• Interpolation of error approximation and correction (coarser to finer):

${v}^{h}←{v}^{h}+{I}_{2h}^{h}{e}^{2h}.$

• Applying few number of Gauss-Seidel iteration method on finer grid ${\mathrm{\Omega }}^{h}$

The schematic representation of restricted (interpolated) operator from finer grid to coarser grid (coarser grid to finer grid) is shown in the Fig. 1(a, b). The usual multigrid V cycle is shown in the Fig. 2.

Fig. 1. a) Restriction, b) interpolation a) b)

Fig. 2. Multigrid V-cycle #### 4.2. Method of solution

The process of multigrid method includes the pressure correction and load balance by adjusting the rigid film thickness, which is carried out on the same grid. For an instant, the iteration procedure consists of the pressure iteration and the adjustment for ${H}_{0}$. The pressure must be relaxed on each grid level. If we denote ${\left(LP\right)}_{i}$ by ${L}_{i}\left(P\right)$ and ignoring superscripts which represent the grid level, the algebraic equation for pressure on each level can be expressed as: ${L}_{i}\left(P\right)={F}_{i}$ ($i=$ 1, 2, ..., $n$-1). For an under-relaxation factor ${c}_{1}$, the process for the pressure iteration can be described as :

(12)
${\stackrel{~}{P}}_{i}={\stackrel{-}{P}}_{i}+{c}_{1}{\delta }_{i},$

where:

$\left\{\begin{array}{l}{\delta }_{i}={\left(\frac{\partial {L}_{i}}{\partial {P}_{i}}\right)}^{-1}{r}_{i},\\ {r}_{i}={F}_{i}-\frac{{\xi }_{i+1/2}{\left({P}_{i+1}-{P}_{i}\right)}^{\frac{1}{n}}-{\xi }_{i-1/2}{\left({P}_{i}-{P}_{i-1}\right)}^{\frac{1}{n}}}{\delta {X}^{1+\frac{1}{n}}}+\frac{{H}_{i}-{H}_{i-1}}{\delta X},\\ \frac{\partial {L}_{i}}{\partial {P}_{i}}=\frac{1}{n}\left[\frac{{\xi }_{i+1/2}{\left({P}_{i+1}-{P}_{i}\right)}^{1/n}-{\xi }_{i-1/2}{\left({P}_{i}-{P}_{i-1}\right)}^{1/n}}{\delta {X}^{1+1/n}}\right]-\left[\frac{\frac{1}{\pi }\left(\stackrel{-}{\rho }{K}_{ij}-{\stackrel{-}{\rho }}_{i-1}{K}_{i-1j}\right)}{\left(\delta X\right)}\right],\end{array}\right\$

and ${\stackrel{-}{P}}_{i}$, ${\stackrel{-}{P}}_{i+1}$ are initial values of pressure in an iterative process and ${\stackrel{~}{P}}_{i-1}$ is the new value of pressure obtained in the iteration. The load balance condition can be achieved by modifying the rigid film thickness ${\stackrel{-}{H}}_{0}$ as follows:

(13)
${\stackrel{-}{H}}_{0}={\stackrel{~}{H}}_{0}+{c}_{1}\left[{G}^{\mathrm{\Delta }}-\frac{\mathrm{\Delta }}{\pi }\sum _{j=1}^{N-1}\left({P}_{j}+{P}_{j+1}\right)\right],$

where ${c}_{1}$ is the relaxation factor and ${G}^{\mathrm{\Delta }}$ is the non-dimensional load on the coarsest grid. The following procedure is employed to determine the pressure cavitation point is given in Zargari et al. . Let ${X}_{c}$ be a point in (${X}_{in}$, ${X}_{out}$), if pressure gradient ${dP/dX|}_{{X}_{c}}=$ 0, then ${X}_{c}$ is the required pressure cavitation point. Otherwise, shift the point ${X}_{c}$ on next interval to the left or right depending on ${dP/dX|}_{{X}_{c}}>$ 0 or ${dP/dX|}_{{X}_{c}}<$ 0 respectively. Next, approximate another cavitation point ${{X}_{c}}^{\text{'}}$ and solve EHL problem. If pressure gradient ${dP/dX|}_{{X\text{'}}_{c}}$ is zero, we get converged solution, otherwise verify the condition:

(14)
${\frac{dP}{dX}|}_{{X}_{c}}×{\frac{dP}{dX}|}_{{X\text{'}}_{c}}<0.$

If the Eq. (14) is not satisfied, repeat above process and find two approximate values ${X}_{c}^{\text{'}}$ and ${X}_{c}^{\text{'}\text{'}}$ for which Eq. (14) is satisfied. A new cavitation point ${X}_{new}$ is obtained by interpolating between the two points ${X}_{c}^{\text{'}}$ and ${X}_{c}^{\text{'}\text{'}}$ which ensures that the pressure gradient nearly equal to zero at ${X}_{new}$ and results in the required smooth solution of EHL problem.

Fig. 3. Pressure and film thickness plot for $W=$ 4E-05, $U=$ 5E-11 a) b)

Fig. 4. Pressure and film thickness plot for $W=$ 3E-05, $U=$ 2.04E-11 a) b)

Fig. 5. Pressure and film thickness plot for $W=$ 4E-05, $U=$ 5E-11 and $n=$ 0.9 a) b)

Fig. 6. Pressure and film thickness plot for $W=$ 3E-05, $U=$ 2.04E-11 and $n=$ 0.9 a) b)

Fig. 7. a) Convergence of pressure for $W=$ 4E-05, $U=$ 5E-11 and $n=$ 0.9, b) convergence of pressure for $W=$ 3E-05, $U=$ 2.04E-11 and $n=$ 0.9 a) b)

Fig. 8. a) Convergence of pressure for $W=$ 4E-05, $U=$ 5E-11 and $n=$ 1.0, b) convergence of pressure for $W=$ 3E-05, $U=$ 2.04E-11 and $n=$ 1.0 a) b)

#### 5. Results and discussion

The paper presents, the computation of minimum film thickness of grease EHL with different load and speed using multigrid method with FAS, as given in section 4.2. The solution has been calculated on the domain –2.5 < $X$ < 1.5 with different grid size $N$. Table 1 and 2 presents the calculated minimum film thickness for different rheological index ($n$), load ($W$), speed ($U$) and compared with the standard result of Dowson and Higginson . The minimum value of film thickness increases as the index $n$ increases. The pressure profiles have slight variation, corresponding to Newtonian case and pressure spike decreases as index $n$ decreases. Figs. 3 and 4 represents pressure and film thickness profile for different load $W=$ 4E-5, 3E-5, speed $U=$ 5E-11, 2.04E-11, index $n=$ 1.0, 0.95, 0.9, 0.85, 0.8 on the grid size of $N=$ 128. Figs. 5 and 6 represent pressure and film thickness profile for different load ($W$), speed ($U$) taking $n=$ 0.9 over the grid size of $N=$ 32, 64, 128.

Table 1. The minimum film thickness for different rheological index and grid

 $n$ $N$ $W$ $U$ $G$ ${H}_{min}$  ${H}_{min}$ (calculated) 1.0 32 3E-05 2.04E-11 5373 0.384185 3.30E-01 0.95 32 2.12E-01 0.9 32 1.30E-01 0.85 32 7.49E-02 0.8 32 4.48E-02 1.0 64 0.384185 3.34E-01 0.95 64 2.22E-01 0.9 64 1.43E-01 0.85 64 9.03E-02 0.8 64 5.14E-02 1.0 128 0.384185 3.24E-01 0.95 128 2.04E-01 0.9 128 1.29E-01 0.85 128 7.57E-02 0.8 128 4.30E-02

Table 2. The minimum film thickness for different rheological index and grid

 $n$ $N$ $W$ $U$ $G$ ${H}_{min}$  ${H}_{min}$ (calculated) 1.0 32 3E-05 2.04E-11 5373 0.384185 2.56E-01 0.95 32 1.65E-01 0.9 32 1.02E-01 0.85 32 5.93E-02 0.8 32 3.51E-02 1.0 64 0.384185 2.61E-01 0.95 64 1.75E-01 0.9 64 1.14E-01 0.85 64 7.28E-02 0.8 64 4.01E-02 1.0 128 0.384185 2.49E-01 0.95 128 1.53E-01 0.9 128 9.58E-02 0.85 128 5.95E-02 0.8 128 3.47E-02

The convergence of pressure with different grid over a number of iteration is plotted in the Fig. 7(a, b) for $W=$ 4E-5, 3E-5, $U=$5E-11, 2.04E-11, $n=$ 0.9 and Fig. 8(a, b) for $W=$ 4E-5, 3E-5, $U=$ 5E-11, 2.04E-11, $n=$1.0. In this method the convergence of pressure profile and film thickness is obtained up to $n=$ 0.8. As the mesh points on the finer grid increases, the converging solution decreases. For decrease in the value of index $n$ and increase in the mesh points of finer grid, the pressure spike shrinks and moves toward the outlet.

#### 6. Conclusions

In the present paper, the Reynolds equation and elastic equation is solved to obtain a minimum film thickness of EHL in line contact with grease as a lubricant and can be concluded as: 1) Minimum film thickness increases as the rheological index $n$ increases. 2) Pressure spike decreases as the rheological index $n$ decreases. 3) The advantage of multigrid FAS method incorporates naturally and overcomes the drawbacks of Newton-Raphson scheme.

1. Dowson D., Higginson G. R. Elastohydrodynamic Lubrication, The Fundamentals of Roller and Gear Lubrication. Pergamon Press, Oxford, 1966. [Search CrossRef]
2. Hamrock B. J., Dowson D. Isothermal elastohydrodynamic lubrication of point contact. 1. Theoretical formulation. Journal of Tribology, Vol. 98, 1976, p. 223-229. [Search CrossRef]
3. Hamrock B. J., Dowson D. Isothermal elastohydrodynamic lubrication of point contact. 2. Ellipticity parameters results. Journal of Tribology, Vol. 98, 1976, p. 375-83. [Search CrossRef]
4. Venner C. H., Lubrecht A. A. Multilevel Methods in Lubrication. Elsevier, Amsterdam, 2000. [Search CrossRef]
5. Lugt P. M., Morales Espejel A review of elasto-hydrodynamic lubrication theory. Tribology Transactions, Vol. 54, Issue 3, 2011, p. 470-496. [Search CrossRef]
6. Kauzlarich J. J., Greenwood J. A. Elastohydrodynamic lubrication with Herschel-Bulkley model grease. ASLE Transactions, Vol. 15, Issue 4, 1972, p. 267-277. [Search CrossRef]
7. Zhu W. S., Ying T. N. A Theoretical and experimental study of EHL lubricated with grease. Journal of Tribology, Vol. 110, Issue 1, 1988, p. 38-43. [Search CrossRef]
8. Cheng J. Elastohydrodynamic grease lubrication theory and numerical solution in line contacts. Tribology Transactions, Vol. 37, Issue 4, 1991, p. 711-718. [Search CrossRef]
9. Kaneta M., Ogata T., Takubo Y., Naka M. Effects of thickener structure on grease elastohydrodynamic lubricant films. Journal of Engineering Tribology, Vol. 214, 2000, p. 327-336. [Search CrossRef]
10. Cann P. M. E., Lubrecht A. A. Bearing performance limits with grease lubrication: the interaction of bearing design, operating conditions and grease properties. Journal of Physics D: Applied Physics., Vol. 40, 2007, p. 5446-5451. [Search CrossRef]
11. Lugt P. A review on grease lubrication in rolling bearings. Tribology Transactions, Vol. 52, 2009, p. 470-480. [Search CrossRef]
12. Cen H., Lugt P. M., Morales Espejel On the film thickness of grease-lubricated contacts at low speeds. Tribology Transactions, Vol. 57, 2014, p. 668-678. [Search CrossRef]
13. Wen S., Ping Huang Principles of Tribology. Wiley, 2011. [Publisher]
14. Roelands C. J. A. Correlational Aspects of the Viscosity-Temperature-Pressure Relationship of Lubricating Oils. Doctoral Dissertation, TU Delft, Delft University of Technology, 1966. [Search CrossRef]
15. Lubrecht A. A., Ten Napel W. E., Bosma R. The influence of longitudinal and transverse roughness on the elastohydrodynamic lubrication of circular contacts. Transactions of the ASME, Journal of Tribology, Vol. 110, 1988, p. 421-426. [Publisher]
16. Venner C. H. Multilevel Solution of the EHL Line and Point Contact Problems. Doctoral Dissertation, University of Twente, 1991. [Search CrossRef]
17. Venner C. H. Higher-order multilevel solvers for the EHL line and point contact problems. ASME Journal of Tribology, Vol. 116, Issue 4, 1994, p. 741-750. [Search CrossRef]
18. Nurgat E., Berzins M., Laurence S. Solving EHL Problems Using Iterative, Multigrid and Homotopy Methods. Journal of Tribology, Vol. 121, 1999, p. 28-33. [Search CrossRef]
19. Briggs W. L., Van E. H., Mccormick S. F. A Multigrid Tutorial. Second edition, Society for Industrial and Applied Mathematics, Philadelphia, PA, 2000. [Publisher]
20. Awati V. B., Shankar N., Mahesh Kumar N. Multigrid method for the solution of EHL line contact with bio-based oils as lubricants. Applied Mathematics and Nonlinear Sciences, Vol. 2, 2016, p. 359-368. [Search CrossRef]
21. Zargari E. A., Jimack P. K., Walkley M. A. An investigation of the film thickness calculation for elastohydrodynamic lubrication problems. International Journal for Numerical Methods in Fluids, 2007, p. 1-6. [Search CrossRef]

#### Cited By

 Australian Journal of Mechanical Engineering Vishwanath B. Awati, Shankar Naik 2020