 Open Access
 Total Downloads : 13
 Authors : Sharanbasappa, E. Sujith Prasad
 Paper ID : IJERTCONV3IS19142
 Volume & Issue : ICESMART – 2015 (Volume 3 – Issue 19)
 Published (First Online): 24042018
 ISSN (Online) : 22780181
 Publisher Name : IJERT
 License: This work is licensed under a Creative Commons Attribution 4.0 International License
Performance Analysis of NonNewtonian Thermoelastohydrodynamic Journal Bearings
Sharanbasappa

ech Machine Design T.John Institute of Technology Bengaluru, India
E. Sujith Prasad
Professor,Dept of Mechanical engg T.John Institute of Technology Bengaluru, India
Abstract Todays compact design of rotating systems requires their supporting journal bearings to operate under more severe operating conditions in which fluidfilm thickness becomes thin enough. The strong dependency of lubricant viscosity on fluid film temperature necessitates a better knowledge of fluidfilm temperature rise and its variation across and along the fluid film for the determination of realistic lubricant viscosity. In such case, thermal consideration becomes most important and realistic condition in the analysis and design of journal bearings operating under full lubrication condition. In addition to this, the structural distortion due to thermoelastic deformation as a result of hydrodynamic fluidfilm pressure.
Further, high molecular weight polymers such as polymethacrylate or polyisobutylene are added to mineral lubricants in order to conserve optimum pro perties under different operating conditions. These polymer added oil no longer behaves as Newtonian fluid and there exists nonlinear relationships between shear stress and shear strain rate. Therefore, a study that predicts the influence of nonNewtonian behavior of the lubricant on journal bearing performance is also more essential. Due to the present day trend towards higher speed as well as the use of unconventional low viscosity lubricants such as water, liquid metal and synthetic oils etc., it is possible to arrive at such a situation where flow is laminar but fluid inertia forces cannot be neglected. Thus, a nonNewtonian thermoelastohydrodynamic (TEHD) analysis of journal bearing which couples the combined influence of surface roughness, thermal, elastic distortion of bearing, nonNewtonian behavior of lubricant and fluid inertia effects is more realistic and appropriate. The present work predicts the influence of surface roughness and fluidinertia on nonNewtonian thermoelastohydrodynamic (TEHD) performance of hydrodynamic journal bearing systems under more realistic operating condition of bearing by considering bearing flexibility, thermal and nonNewtonian behavior of lubricant.

INTRODUCTION

Todays compact design of rotating systems requires their supporting journal bearings to operate under more severe operating conditions in which fluidfilm thickness becomes thin enough. As a result, some asperities on the mating surfaces begin to interfere inevitably, hence temperature increases due to asperity interaction and fluid shear action. This temperature rise in the lubricant results into the variation of viscosity along and across the fluidfilm. The strong dependency of lubricant viscosity on fluidfilm temperature necessitates a better knowledge of fluidfilm temperature rise and its variation across and along the fluid
film for the determination of realistic lubricant viscosity. In such case, thermal consideration becomes most important and realistic condition in the analysis and design of journal bearings. In addition to this, the structural distortion due to thermoelastic deformation as a result of hydrodynamic fluid film pressure may significantly cause the change of the fluid film thickness and hence the bearing performance.
From last few decades, the inclusion of surface roughness effects into the analysis of bearing lubrication problems has been conveniently carried out using statistical method. The following subsection details the inclusion of surface roughness into the analysis of lubrication problems through statistical analysis.
Due to unprecedented technological advances during last few decades, the operating conditions of rotating machines are becoming more stringent and exact. This has necessitated many new developments in the area of hydrodynamic journal bearings, since large number of machine parts use hydrodynamic journal bearing and operates continuously at high speed. As a consequence, a considerable amount of research activities have been devoted towards the study of hydrodynamic journal bearing systems and analysis has been carried out by incorporating many physical effects to get much more realistic performance data.
Furthermore, most of the available studies which deals with surface roughness effects in the analysis of hydrodynamic journal bearings have been carried out with the general assumption of negligible fluid inertia forces. However, due to the present day trend towards higher speed as well as the use of unconventional low viscosity lubricants such as water, liquid metal and synthetic oils etc., the fluid inertia becomes important for some range of moderately large Reynolds number. When the reduced
STATISTICAL ANALYSIS OF SURFACE ROUGHNESS
Surface roughness in the form of fine irregularities exists on all finished surfaces of tribological elements and it is inherent to the characteristics of surface finishing process such as turning, grinding, polishing, sand blasting, etc. Fortunately, for many rough surfaces of engineering importance, the surface height, curvature and peak height distributions are random in nature and are represented by a Gaussian probability distribution
SURFACES WITH DIRECTIONAL PATTERNS
Many engineering surfaces exhibits roughness with directional patterns those resulting from different manufacturing process. The directional patterns (roughness orientations) of the rough surfaces may be isotropic or anisotropic depending upon the nature of the surface finishing process. The anisotropic roughness orientations are generally classified as transverse and longitudinal roughness patterns depending upon the principle directions of the surface textures with the sliding direction of opposing surfaces.
To study the effects of roughness orientations on lubrication problems, Peklenik used a parameter known as
surface pattern parameter, , which is also popularly known
as Peklenik number. This surface pattern parameter or Peklenik number is defined as
LITERATURE SURVEY
The literature review presented in this thesis mainly concerns with theoretical investigations of mixed lubrication performance of hydrodynamic journal bearing including the influence of surface roughness coupled with the effects of bearing shell flexibility, nonNewtonian behavior of lubricant and variation in viscosity of the lubricant due to temperature rise in fluid film. Over the last few decades numerous studies on hydrodynamic journal bearings have been carried out and reported in literature. Since the available literature in the area of these journal bearings is quite vast and abundant it is rather difficulty to present all these information. Therefore, few important studies which are most relevant to the present work are reviewed in the following section.
Among the many studies related to the non Newtonian behavior of lubricant on performance of journal
0.5x
0.5 y
(1)
bearing, Dien and Elrod [5] derived a generalized Reynolds equation for the nonNewtonian fluids using perturbation technique. They selected the power law model as the
where
0.5x
and
0.5 y
are the autocorrelation lengths
constitutive relation for the nonNewtonian fluids and carried out a study on slider and journal bearings. Using this
of x – profile (i.e. along the sliding direction of opposing
surface) and y profile (i.e.in the direction perpendicular to
the sliding direction of opposing surface) at which their value reduces to 50 percent of its original value.
y
Sliding

Transversely oriented
y
Sliding

Isotropically oriented
y
generalized Reynolds equation, Wada and Hayashi [7,8] derived a modified Reynolds equation for nonNewtonian lubricants obeying cubic shear stress law model [7]. It has been clarified experimentally that the flow characteristics of pseudoplastic fluids are approximately expressed by a cubic shear stress model containing first and third powers of shear stress [8]. The experimentally measured pressure distribution and load carrying capacity of a journal bearing were found to match with theoretical results obtained from cubic shear stress model. Tayal solved the NavierStokes equations by finite element method. They used power law model to study the static performance characteristics of journal bearing. Their results demonstrated that the static performance characteristics of bearing are significantly affected by non Newtonian behavior of lubricant.
The literature review presented in this subsection clearly indicates that the values of critical mass and threshold speed obtained from linear analysis do not set correct stability margins. On the other hand the nonlinear analysis of transient response generally predicts a higher value of stability parameters than that obtained from linear analysis. This in fact depends on the consideration of various factors such as bearing/rotor flexibility effect, nonNewtonian effect, loading conditions etc. in the analysis.
METHODOLOGY
AVERAGE FLUIDFILM THICKNESS
The coordinate system and geometry of the journal bearing system along with rough surface profile is shown in
Fig. 2.1. The local fluidfilm thickness hL
at any point

Longitudinally oriented
Sliding
inside the clearance space between rough journal and bearing surfaces is sum of the nominal fluidfilm thickness ( h ) and the roughness amplitudes of journal and bearing surfaces ( J , b ) measured from their mean levels as shown in Fig.
2 may be expressed as
Fig. 1.3 – Contact areas of different roughness orientations on
unwrapped bearing surface
hL h J b h
(2)
AVERAGEREYNOLDS EQUATION
For the laminar flow of an incompressible, non Newtonian lubricant, the modified average Reynolds
where J b is the combined roughness
amplitude of two surfaces. The nominal fluidfilm thickness ( h ) is defined as the distance between the mean levels of two
equation in terms of flow factors, average fluidfilm thickness and inertia term can be expressed as (see Appendix1)
f s T e
f s T e
surfaces as shown in Fig.2.
h 3 p
h 3 p h h R*
x f
T – h 2Gx
h 2Gy
The average fluidfilm thickness hT , which is equal to the expected or mean value of local fluidfilm thickness ( hL ),
12
12 2 2 t 12 T
Where
Where
is the local average (expected) viscosity
is the local average (expected) viscosity
T
hT E[h ]
(h ) ( ) d
corresponds to the local average temperature It may be noted that for the flow of lubricant through the clearance space between two smooth surfaces, the pressure flow factors
(3)
where ( ) is the probability density function of
x y 1 and shear flow factor s 0 . Then the
combined roughness .
As most of the engineering surfaces follow Gaussian height distribution, the present work assumes Gaussian height distribution of the surfaces. The probability density function of this Gaussian distribution is expressed as
modified Reynolds equation (Eq. (5)) reduces to the generalized Reynolds equation governing the flow of lubricant in the clearance space of a smooth journal bearing.
The parameters Gx and Gy appearing in last term of Eq. (5) are known as inertia functions in x and y directions.
( )
1
2
e 2 2 2
(4)
FINITE ELEMENT FORMULATION
Lubricant flow field is discretized by using fournoded quadrilateral isoparametric elements. The fluidfilm pressure variation is assumed to vary linearly over an element. Then,
where is the combined rms or standard deviation of
roughness and is expressed as ( 2 2 )1 2 .
the unknown pressure is expressed in nondimensional form as
J b
l
l
ne
p p j N j j1
n
n
l
l
Where N j is elemental shape function and e is
number of nodes per element of twodimensional flowfield solution domain.
Using Galerkins orthogonality conditions for the modified average Reynolds equation (Eq. (5)) and following the usual assembly procedure for all the elements, the global system equation for the entire lubricant flow field is expressed as
F
p QRH XJ RXJ ZJ RZJ G
Fig. Bearing geometry and surface profile
where,
F
p
Q
= Assembled fluidity matrix,
= Nodal pressure vector,
= Nodal flow vector,
RH
J
J
J
J
hydrodynamic terms,
= Column vectors due to

At the trailing edge of the positive region traditional Reynolds boundary conditions are employed,
p = p =0.0
RX
,RZ
= Global right hand side vectors due to journal center velocities

MEAN VELOCITIES DUE TO PRESSURE INDUCED
FLOW
G
= Global right hand side vector due to inertia functions
The mean velocity components due to pressure induced flow of lubricant in x and y directions are expressed (Appendix1) as
For an
eth element, the elements of the above matrices
2 p *
are defined as follows:
Um
h
12 x
f Re hT G
12 x
e h 3
Ni N j
h 3 Ni N j 2 p *
Fij
x 12
y 12
J d d
V h
f Re hT G
(2.6b)
m
12
y
12 y
Re h

s Ni
J d d
Hi 2 T
The parameters Gx
and Gy appearing in Eqs. (2.5)
(2.6c)
1 h
and (2.11) are called inertia functions in x and y directions. These inertia functions are expressed in nondimensional
e
e
RX e
1 erf Ni cos d d
form (Appendix1) as
Ji 2
2
2
s s s
s s s
A Gx hT 12U m 1 s U m
6 hT
m Um 1 2
(2.6d)
5
U 5
3 1
R e
1 erf h N
sin d d
hT 2 1 s U
s 6V h
Um hT 12U
(1 s ) V m
ZJi
2 1
2 i
3 m 5 m T 2 5 m
Ae
(2.6e)
e R 2 N
N
(2.12a)
Gi e hT Gx i Gy i J d d
12
12 hT
Um
hT 12
s V m
V m 12
m
m
s hT
(2.6f)
Gy V m
5 2
U 1
2 5
2 5 Um 1
h 3 p
p
R*
Qie h 3
x
y
hT s e h 2 Gx Gy Nide
h
12 V
e
12a
12
T V m s
V mhT m
2 5
l
l
where i , j =1,2 ne (number of nodes per element).
BOUNDARY CONDITIONS
The boundary conditions relevant to the lubricant flow field (Eq. 2.6a) are as described below:

At external boundary of the bearing, the pressure is set to the atmospheric pressure,
p  1.0 = 0.0
1.1.1 FLUIDFILM VELOCITY COMPONENTS
The flow of the lubricant between two rough surfaces of bearing and journal can be modeled by an equivalent flow model as shown in Fig. 2.4. The equivalent model is defined as two smooth surfaces separated by a clearance equal to the average fluidfilm thickness ( hT ). Then, the mean or expected velocity components can be obtained by modifying the Poiseuille and Couette terms in the expression of local velocity components using Patir and Chengs flow factors. They can be expressed in nondimensional form as
T
T
h 2 p 2

At the supply hole,
p ps 1
u x 12 z
z z h
sz

Periodicity is enforced where the bearing wraps around
h 2 p 2
The viscosity of nonNewtonian lubricant is
v y 12 z

z
(2.13b)
described by an apparent viscosity ( a ) and is defined as a function of shear strain rate () .
where
z hT
in the equivalent flow model.
a
The fluidfilm velocity component across the fluid film is obtained from the continuity equation and it may be expressed in nondimensional form as
1.1.3 TEMPERATUREVISCOSITY RELATION
w~ w zu hT

zv hT
z z
hT udz hT vdz
The viscosity is assumed to be dependent on temperature and is defined by the exponential law [26],
(2.13c)
0
0
r expv (Tf
Tr )
where r
is the viscosity of the lubricant at reference
1.1.2 NONNEWTONIAN LUBRICANT MODEL
temperature
Tr . The above relation may be expressed in
nondimensional form as
The available studies indicate that the power law model accurately approximates the physical behavior of many polymerthickened lubricating oils
BC 1
Entry Hole
Lubricant
BC 4 BC 5
[5]. Therefore, the power law model is used in thepresent study to express the nonlinear relationship R2
R1
between shear stress and shear strain rate of non
Fluid Journal Interface
Fluid Bearing Interface
Newtonian lubricants. The constitutive relation for a simple nonNewtonian lubricant is expressed in non
dimensional form as
BC 6
BC 7
m()n
BC 2
BC 3
Journal
BC 5
(2.14a)
where n and m are known as the power law index
Fig. 3 Thermal boundaries of bigend bearing system
and consistency index respectively. The shear strain
Tf
273.12 / Tr
rate ( ) is made independent of direction by
exp v 1
considering it as a function of the second strain
r
1 273.12 / Tr
invariant of shear strain rate, nondimensional form as
I 2 and is expressed in
where
(2.20a)
v is the nondimensional temperatureviscosity
2
2 1 2
coefficient determined from the known values of viscosity at
1
u
1 v
specified temperature for a particular lubricant and is
h
z
h
z
expressed as
T
T
v v
Tr
(1 273.14 / Tr )
Differentiation of Eqs. (2.13a) and (2.13b) and substitution into the above equation yields
For the nonNewtonian lubricant, the temperature viscosity relation may be expressed as
2 2 12
Tf
273.12 / Tr
p h 2 1
p h 2
1
a exp v 1
z
1 s
z
1 273.12 / Tr
x h
2 h
h y h
2
T T
T
T
where a
(2.20b)
is the apparent viscosity of nonNewtonian
lubricant. It may be noted that, for Newtonian lubricant
a 1
and hence for Newtonianthermal analysis Eq.
(2.20b) reduces to Eq. (2.20a). On the other hand for non
Newtonian isothermal analysis Tf 1 and Eq. (2.20b) reduces to a , the apparent viscosity of lubricant.
PERFORMANCE CHARACTERISTICS
After establishing the steadystate matched solutions for the pressure, thermal and elastic deformation fields, the static and dynamic performance characteristics of journal bearing system are computed. In the present work the performance characteristics of hydrodynamic journal bearing is studied under both fully lubricated lubrication conditions of the bearing. For the fully lubricated condition of the bearing, static performance characteristics are computed and presented.
Static Performance Characteristics
The static performance characteristic includes the load carrying capacity of fluidfilm pressure, bearing flow, frictional torque at journal surface due to shear viscous of lubricant, midfilm temperature, etc. and dynamic performance characteristics includes fluidfilm stiffness and damping coefficients and critical journal mass.
SOLUTION PROCEDURE
A computer program based on the analysis is developed to obtain the performance characteristics of hydrodynamic journal bearings. In case of a nonNewtonian TEHD analysis of a journal bearing, the fluidfilm thickness, pressure, temperature, lubricant viscosity and bearing deformation are interdependent variables. Therefore, it requires the simultaneous solution of all the governing equations to establish the matched solution of modified average Reynolds, elasticity, energy and heat conduction equations.
RESULTSAND DISCUSSIONS
The mathematical models and solution algorithms presented in the previous chapters have been used to compute the static and dynamic performance characteristics of a journal bearing. The influence of surface roughness under the realistic operating conditions of bearing is presented and discussed by considering the bearing deformation, crossfilm viscosity variation due to nonNewtonian behavior of the lubricant and rise in fluidfilm temperature and fluid inertia. The variations in the bearing performance characteristics are presented for various values of surface roughness parameter,
(i.e. for the inverse rms values of combined roughness heights) and roughness orientations such as

Transversely oriented roughness pattern ( = 1/6)

Isotropically oriented roughness patterns (=1)

Longitudinally oriented roughness patterns ( = 6)
The static performance characteristics of journal bearing are computed for the geometric and operating parameters of the bearing and the results are computed for these nondimensional values unless otherwise mentioned in the text or figures.
Further, since the influence of surface roughness on lubrication process of bearing under the fully lubricated condition is more clearly understandable, the results showing the influence of surface roughness, bearing flexibility, thermal, nonNewtonian behavior of lubricant and fluid inertia on fully lubricated bearing performances are presented and discussed in the first part. The computed results of performance characteristics are presented and discussed in the following sections.
PERFORMANCE CHARACTERISTICS UNDER FULL
LUBRICATION
The results of static performanc characteristics of hydrodynamic journal bearing as computed from IHD,THD and TEHD analysis of bearing operating under fully lubricated condition of bearing are presented in Fig 4 These results are discussed in the following sections. As seen from Fig. 4.12, the transverse, isotropic and longitudinal roughness patterns enhance the load carrying capacity of twosided type rough bearing as compared to the corresponding smooth bearing for both Newtonian and nonNewtonian lubricant cases and this trend is also observed for both fluid inertialess and inertia solutions. Irrespective of the analysis cases considered, the longitudinal roughness, which restrict the dominant pressure induced flow of lubricant in axial direction, is observed to provide maximum enhancement in the load carrying capacity. Though the transverse roughness pattern with longer asperities in axial direction enhance the pressure induced axial flow, it is observed to provide a marginal enhancement in the load carrying capacity of bearing by restricting the pressure induced circumferential flow of lubricant at converging section of the bearing.
1, 0.5,Cd 0, 6,Vrj 0.5
Line
1, 0.5,Cd 0, 6,Vrj 0.5
Line
62
52
Smooth
Line
Smooth
Line
Colour
Colour
1 / 6 e
1 / 6 e
R 0
R 0
*
*
1
1
Type
Type
R 0.6
R 0.6
*
*
42
e
e
6
6
Fo
32
22
12
a
a
2
1 3 5 7 9
13
1, 0.5,Cd 0, 6,Vrj 0.5
1, 0.5,Cd 0, 6,Vrj 0.5
11
9
Fo
LiSmnoeoth
LiSmnoeoth
Colo1u/ 6r
Colo1u/ 6r
Lin* e
Ty
Lin* e
Ty
7
1
6
1
6
5
3
1
b 1 3 5 7 9
13
1, 0.5,Cd 0.01, 6,Vrj 0.5
n=
n=
n=
n=
11 13 15 17
n
n=
n
n=
pe0.2
pe0.2
Re 0
Re 0
R
R
*
*
e
e
11 13 15 17
became considerable when this viscosity variation is considered.

The influence of fluidinertia is to increase the bearing load carrying capacity and lubricant side leakage while the influence of nonNewtonian behavior of lubricant is to reduce these parameters.

The transverse, isotropic and longitudinal roughness patterns enhance the load carrying capacity of two sided type rough bearing as compared to the corresponding smooth bearing for both Newtonian and nonNewtonian lubricant cases and this trend is also observed for both fluid inertialess and inertia solutions.

The longitudinal roughness, which restrict the dominant pressure induced flow of lubricant in axial direction, is observed to provide maximum enhancement in the load carrying capacity.

the transverse roughness pattern with longer asperities in axial direction enhance the pressure induced axial flow, it is observed to provide a marginal enhancement in the load carrying capacity of bearing by restricting the pressure induced circumferential flow of lubricant at converging section of the bearing.


11
9
o
7
5
o
7
5
F
Smooth
Smooth
Line Colour 1 / 6
1
6
3
Line
e
e
TyRp* e0
e
e
R* 0.2
n n=0.
REFERENCES

Nayak P. R., Random Process Model of Rough Surfaces, Trans, ASME, J. Lubr. Technol., Vol. 93, 1971, pp.398407.

Thomas T. R., Rough surfaces, Longman, London, UK,1982

Whitehouse D. J., Handbook of Surface Metrology, Institute of Physics Publishing, Bristal, UK, 1994.

Peklenik ., New Developments in Surface Characterization and Measurements by Means of Random Process Analysis, Proc.
Fig.c1 4.12 – Variation of lFooad carrying capacity
IMechE., Part 3K, Vol.182, 196671968,
1 3 5 7 9 11 13 15 17 [5] Dien I. K., and Elrod H. G., A Generalised SteadyState Reynolds
( ) of rough bearing with nonNewtonian
Equation for NonNewtonian Fluids, with application to Journal Bearings, ASME, J. Lub. Tech., Vol. 105, 1983, pp. 385390.
CONCLUSIONS
The present study investigate the influence of surface roughness on nonNewtonian thermoelastohydrodynamic (TEHD) performance of hydrodynamic journal bearing including fluidinertia effects. The modified form average Reynolds equation is derived in terms of Patir and Chengs flow factors and inertia functions to include the surface roughness and fluidinertia. The mean pressure induced velocity components are also modified to include surface roughness in fluidinertia analysis. The expressions for pressure derivatives which are required to compute fluidfilm dynamic coefficients are also developed.
1. In general influence of fluidinertia on performance characteristics of bearing is insignificant when viscosity variation due to rise in fluidfilm temperature and non Newtonian behavior of lubricant is not considered. It

Buckholz R. H., and Lin J. F., The Effect of journal Bearing Misalignment on Load and Cavitation for NonNewtonian Lubricants, ASME, J. Tribol., Vol.108, 1986, pp. 645654.

Wada S., and Hayashi H., Hydrodynamic Lubrication of Journal Bearings by Pseudoplastic Lubricants; Part 1: Theoretical Studies, Bull JSME, Vol. 14(39), 1971, pp. 268278.

Wada S., and Hayashi H., Hydrodynamic Lubrication of Journal Bearings by Pseudoplastic Lubricants; Part 2: Experimental Study, Bull JSME, Vol. 14(39), 1971, pp. 279286.

Tayal S. P., Sinhasan R., and Singh D. V., Analysis of Hydrodynamic Journal Bearing with NonNewtonian Power Law Lubricants by the Finite Element Method, WEAR, Vol. 71, 1981, pp.1527.

Paranjpe R. S., Analysis of NonNewtonian Effects in Dynamically Loaded Finite Journal Bearings Including Mass Conserving Cavitation, ASME, J. Tribol., Vol.114, 1992, pp. 736746.