## Abstract

The topic of the presented paper aims to demonstrate a new principle of hydrodynamic lubrication in mechanical, thermal and electro-magnetic fields. Up till now, when dealing with the hydrodynamic theory lubrication, many authors of scientific papers have assumed the constant oil dynamic viscosity value without variations caused by temperature crosswise the film thickness. Simultaneously, due to the numerous AFM measurements, it appears that oil temperature gradients and oil viscosity changes in the bearing gap height directions cannot be omitted. Therefore, in this paper, the problem of the viscosity changes across the lubricant thin layer was resolved as the main novelty in principles of mechanical thermal lubrication. The method of solving the mentioned problem was manifested by a general model of semi-analytical solutions of isothermal electro-magneto-elastohydro-dynamic and non-Newtonian, lubrication problem formulated for two deformable rotational surfaces in curvilinear, co-ordinates.

## 1 Introduction and governing equations

This paper presents a semi analytical method of solution of the asymmetrical, laminar, steady, non-Newtonian lubrication flow problem between two rotational and deformable, curvilinear orthogonal movable surfaces in conjugated elasto-hydro-electro-magnetic fields. The fluid flow between the two above mentioned solid surfaces in the electromagnetic field will be described by the 3 momentum equations of equilibrium in a vector form, a fluid continuity equation, and by equation of a conservation of energy equation in a scalar form. Hence, we obtain the following system [1, 2]:

The above-mentioned system of equations is completed by magneto-thermo-elasticity equations describing problem in stresses for two solid surfaces restricting the thin fluid layer. This system consists of the three partial differential equations (4) in a vector form. Heat conductivity equation in a solid body (5) is added to this set of equations and we obtain a system in the following form [2, 3, 4]:

Moreover, we add the Maxwell and Ohm equations as well for the two surfaces as for the thin boundary liquid layer between two surfaces. Thus we get such equations [5]:

We assume the following notations: *μ _{o}* – magnetic permeability in vacuum H/m, T – fluid temperature in K, T

^{*}– solid body temperature in K,

**B**– magnetic induction vector in T,

**N**- magnetization vector A/m,

**E**– electric intensity vector V/m,

**H**– magnetic intensity vector A/m, (– first derivative of magnetization vector with respect to temperature A/mK,

*σ*- electrical conductivity coefficient S/m,

**J**- electric current density in A/m

^{2},

**D**- electric induction vector As/m,

*ρ*- fluid density kg/m

^{3},

*κ*- thermal conductivity coefficient W/mK,

**v**– fluid velocity vector in m/s,

*φ*– dissipation of energy in W/m

_{F}^{3},

**S**- stress tensor in the fluid in Pa,

**u**- displacement vector of the solid body in m, t - time in s, c

_{v}- specific heat in J/kgK. The symbols with an asterisk are related to the solid body.

The relationship between stress tensor **S** and strain rate tensor 2**T**_{d} = **A**_{1}*i*.*e*. constitutive equations are assumed for the lubricant in the following form [3]:

whereas unit tensor **δ**, strain rate tensor **A**_{1} have the following components: *δ _{ij}*,

*Θ*in s

_{ij}^{–1}. For power law of compressible fluid the apparent viscosity

*η*has the form [2, 3, 4]:

_{p}We assume: *n* – dimensionless flow index (0.5, 1.2), *m* = *m*(*n*) fluid consistency coefficient in Pas^{n}.

The Duhamel Neumann relations between the components *τ _{ij}* of the stress tensor S

^{*}of the elastic body on the sleeve and the strain tensor components

*ϵ*take the following form:

_{ij}for *i*,*j* = 1,2,3, where: *δ _{ij}* – Kronecker unit tensor component (

*δ*= 1 for

_{ij}*i*=

*j*and

*δ*= 0 for

_{ij}*i*≠

*j*),

*K*=

*Λ*+ 2

*G*/3,

*G*=

*E*/[2(1+

_{Y}*ϑ*)],

*Λ*=

*E*

_{Y}*ϑ*/(1+

*ϑ*)(1–2

*ϑ*), whereas

*E*- Young‘s Modulus in Pa,

_{Y}*α*– thermal coefficient of linear expansionin

_{T}*K*

^{–1},

*ϑ*– Poisson ratio.

And by virtue of known linear geometrical relations, the strain tensor components *Θ _{ij}*(

*v*) in the oil, and strain components

_{i}*ϵ*(

_{ij}*u*) in the sleeve are mutually connected with the oil velocity components

_{i}*v*and displacement vector of solid body components

_{i}*u*, respectively:

_{i}## 2 Electro-magneto-hydrodynamic equations

We assume: unsymmetrical, incompressible (for invariant **I**_{1} = 0), steady, pseudo-plastic lubrication in electro-magnetic field for power law model without viscous-elastic properties. Oil apparent viscosity *η _{p}* varies in length, width and gap-height directions and depends on pressure, temperature flow shear ratio, and electro-magnetic field [1, 2, 3, 4, 5]. Inertia forces and terms of energy convections are neglected, pressure is constant in gap height direction. The total distance

*ϵ*between the two surfaces is signifucantly smaller than other dimensions of the considered surfaces. Taking into account the layer boundary simplification

_{T}*i*.

*e*. neglecting in (1)-(7) the negligibly small terms of order

*Ψ*= 0.001 (radial clearance) presenting the quotient of characteristic gap height (

_{0}to the radius R of journal, the problem is made in local curvilinear and orthogonal coordinates (

*α*

_{1},

*α*

_{2},

*α*

_{3}) connected with one of the movable surfaces, where

*α*

_{2}denotes the direction of gap height. The parallel and longitudinal intersections of the cooperating deformed, rotational surfaces have curvilinear and non-monotone generating lines. Hence, Lame coefficients have the form: h

_{1}= h

_{1}(

*α*

_{3}), h

_{2}= 1, h

_{3}= h

_{3}(

*α*

_{3}).

We put physical and geometrical dependencies (9), (10), (11) into expanded boundary simplified equations (1)-(7). Hence, we obtain the system of equations of conservation of momentum, continuity, energy, thermo-elasticity, heat transfer, Maxwell simplified results in the following form of curvilinear orthogonal co-ordinates (*α*_{1}, *α*_{2}, *α*_{3}):

where the length, width and gap- height directions, are limited respectively: 0 < *α*_{1} ≤ 2*π*, –*b _{m}* ≤

*α*

_{3}

*b*, 0 ≤

_{s}*α*

_{2}≤

*ϵ*and

_{T}*i*= 1,2,3. The system of equations (12)((17) contains the following unknown dimensional values: pressure

*p*(

*α*

_{1},

*α*

_{3}), temperature in oil

*T*(

*α*

_{1},

*α*

_{2},

*α*

_{3}) and sleeve

*T*

^{*}, three oil velocity components

*v*(

_{i}*α*

_{1},

*α*

_{2},

*α*

_{3}) and three layer displacement components

*u*(

_{i}*α*

_{1},

*α*

_{2},

*α*

_{3}), for

*i*= 1,2,3 in three curvilinear, orthogonal dimensional directions:

*α*

_{1},

*α*

_{2},

*α*

_{3}.

## 3 Boundary conditions

The lubricant flow in the bearing gap is generated by the rotation of a rotational, curvilinear journal. Hence the boundary conditions for lubricant velocity components have the form:

for *α*_{2} = *ϵ _{T}*,

*v*

_{2}= 0 for

*α*

_{2}= 0,

and for *α*_{2} = *ϵ _{T}*,

*v*

_{3}= 0 for

*α*

_{2}= 0,

*v*

_{3}= 0 for

*α*

_{2}=

*ϵ*,

_{T}where *ω* denotes the angular velocity of rotational journal, in the circumferential direction.

Decrements or increments that are above the characteristic environmental temperature *T*_{0} have constant value *f _{c}* on the journal surface and variable unknown values

*f*(

_{p}*α*

_{1},

*α*

_{3}) on the sleeve surface. Heat flux is transferred from the rotational surface of bearing journal into lubricant. Hence the boundary temperature values are as follows:

where *Δ*f – temperature difference across the film thickness, *υ* – heat transfer coefficient in W/m^{2}K, q* _{c}* - heat flux density in oil on journal surface in W/m

^{2}.

The elastic layer of bearing alloy is laying on the rigid ring in place *α*_{2} = *c*_{3}, and therefore the contact surface of these bodies is not deformed by the pressure. The contact surface *α*_{2} = *c*_{3} is yet deformed by variable temperature, which implies the shape and the total volume change of the rigid ring. These deformations are not considered now.

The temperatures and surfaces localization is presented in Figure 1.

The boundary conditions in the elastic layer have the following form:

where *i* = 1,2,3. Moreover: *T*_{0} - average ambient temperature, *q _{p}* - heat flux density across the elastic layer located on the sleeve surface,

*f*

_{2}(

*α*

_{1},

*α*

_{3}) - temperature increases above the temperature

*T*

_{0}on the external surface of the elastic layer for

*α*

_{2}=

*c*

_{2}, additionally

*e*

_{2}(

*α*

_{1},

*α*

_{3}) temperature increases in the excess of temperature

*T*

_{0}on the bottom surface of the elastic layer in the place of

*α*

_{2}=

*c*

_{3}.

## 4 Oil consistency and viscosity variations

Temperature decreases and magnetic induction, electric intensity increases the oil consistency m and dynamic viscosity coefficient *η*. This phenomenon is described in the following form:

with coefficient *δ _{B}*,

*δ*of magnetic, electric influence on the oil viscosity and coefficient

_{E}*δ*of temperature influence on the oil viscosity Mentioned phenomena are confirmed by virtue of authors own measurements.

_{T}## 5 Sketch of semi-analytical solutions

After double integration of two equations (15) in *α*_{2} direction we apply the low heat flux continuity (*q _{p}* =

*q*) in the contact region between the elastic alloy layer and lubricant, and finally we obtain the external displacement of the sleeve surface in the following dimensional form:

_{c}where *ϵ _{s}* =

*c*

_{3}–

*c*

_{2}, 0 <

*α*

_{1}< 2

*π*, –

*b*<

_{m}*α*

_{3}< +

*b*. Solution (25) valid if and only if

_{m}*E*- Young modulus and solid body coefficients

_{Y}*α*,

_{T}*κ*

^{*},

*ϑ*- Poisson ratio, are independent of

*α*

_{2}, The hydrodynamic pressure p, the load carrying capacity and temperature T distributions with oil velocity components v

_{1}, v

_{2}, v

_{3}in machine slide journal bearing are solved by the means of the quasi half numerical Frobenius small parameter method implemented by the finite difference method using Mathcad 15 Professional Program. A system of partial differential equations (12)-(14) defined in thin space (

*α*

_{1},

*α*

_{2},

*α*

_{3},

*t*) between two movable rotational surfaces and for continuous

*m*(

*n*),

*η*==

*m*(

*n*= 1), has an analytical solution in the form of the following infinite uniform convergent power function series in following form:

where
*q _{c}* =

*q*/(1 –

*g*); for

*i*= 1,2,3;

*k*= 0,1,2,…,

*g*=

*g(n)*≡

*(n*-1)/2 where for 0 <

*n*≤ 3/2 the small parameter

*g*is less than +1/4 and greater than –1/2. In order to achieve the linearization of apparent viscosity (16), we put functions (26a) into Eq. (16) and expand function in closed interval [1,

*n*] or [

*n*, 1] for 0 <

*n*≤ 3/2 using Taylor series in the neighborhood of point

*n*= 1 with respect to the small parameter in the following form:

where *η*_{1}(T, B, E) = *η*_{1}(*α*_{1}, *α*_{2}, *α*_{3}) and functions *η _{prk}* ≡

*η*(

_{prk}*v*

_{1},

*v*

_{3}) for

*k*= 1, 2, … – are the dimensionless expansion coefficients of dynamic viscosity dependent on

*α*

_{1},

*α*

_{2},

*α*

_{3}, whereas for

*k*= 0

*η*

_{pr0}= 1. Such functions after involved calculations are determined in analytical form. Putting the infinite series (26), (27) into the system (12)-(14), multiplying the series by Cauchy methods and equating the terms by the same powers

*k*= 0,1,2,… of the small parameter g, and then using the L. Kronecker symbol

*i*.

*e*.

*δ*

_{k0}= 1 for

*k*= 0,

*δ*

_{k0}= 0 for

*k*= 1,2,…, we obtain the linear differential systems of equations for

*k*= 0,1,2,…

*i*= 1,3 in following form:

with

For *k* = 0 the abovementioned system of equations (28)-(31) determines functions:
*k* = 1,2,… and *i* = 1,2,3 the system of equations (28)-(31) determines the first, second,… corrections:
*T*^{(k)}, *p*^{(k)} caused by the non-Newtonian oil. The boundary conditions for *k* = 0 are as follows:

for *α*_{2} = 0,

for *α*_{2} = *ε _{T}*;

and for k = 1,2,… we have:

The total gap height *ϵ _{T}* includes surface displacement (25). Integrating twice (28) and (29) once with respect to the variable

*α*

_{2}, then the oil velocity components for

*k*= 0,1,2,… by virtue of boundary conditions (32a, b) have the following form:

for *i* = 1, 3; 0 < *α*_{21} < *α*_{2} < *ϵ _{T}*, 0 <

*α*

_{1}< 2

*π*, –

*b*<

_{m}*α*

_{3}< +

*b*, whereas we denote:

_{m}It is easy to see, that:

Imposing the boundary condition (32ab) in point *α*_{2} = *ϵ _{T}* on the solution (34) and taking into account the law of integral differentiation with variable limits of integration, we obtain the following modified Reynolds equation, which determines the unknown function

*p*

^{(k)}(

*α*

_{1},

*α*

_{3}):

with

for i = 1,3; k = 1,2,3.

To obtain the temperature functions *T*^{(k)}(*α*_{1}*α*_{2}, *α*_{3}) for k = 0,1,2,…, we put solutions of oil velocity components (33), (34) for i = 1,3; into the temperature equation (30). After double integration with respect to the variable *α*_{2}*i*.*e*. in film thickness direction, using boundary conditions (19ab) in the form: T^{(k)} = *f _{c}*

*δ*

_{k0}for

*α*

_{2}= 0, T

^{(k)}=

*f*

_{p}*δ*

_{k0}for (

*α*

_{2}=

*ϵ*, and after the term ordering, we finally obtain:

_{T}where *s* ≡ *α*_{2}/*ϵ _{T}*, 0 ≤

*s*≤ 1. Non linear differential equation (38) determines temperature functions T

^{(k)}for k = 0,1,2,… Imposing the condition (19c) on the temperature T

^{(k)}determined by the (38), we obtain the unknown temperature

*f*(

_{p}*α*

_{1},

*α*

_{3}) on the sleeve.

To obtain the oil velocity components (33), hydrodynamic pressure (36) and temperature (38) for Newtonian oil, and taking into account oil viscosity variations crosswise the gap height, we put k = 0 and hence:
*p*^{(0)}, temperature *T*^{(0)}. Assuming cylindrical coordinates: *α*_{1} = *ϕ*, *α*_{2} = *r*, *α*_{3} = *z*, *h*_{1} = *R*, *h*_{2} = *h*_{3} = 1, where *R* is the radius of the cylindrical journal, we put the temperature *T*^{(0)}, determined from the reduced Eq. (38), into the reduced pressure Equation (36). Finally, the Eq. (36) determines the pressure *p*(*ϕ*, *z*), 0 ≤ *ϕ* ≤ 2*π*, –*b*≤ *z* ≤ *b*, and tends to the following form:

0≤ *ϕ* ≤ 2*π*, –*b* ≤ *z* ≤ *b*, *p* = *p*(*ϕ*, *z*), *η*_{1BE} = exp(*δ _{B}B* +

*δ*).

_{E}EFor *δ _{T}* = 0, viscosity is constant in gap height direction and (39) tends to the well known classical form of Reynolds equation.

## 6 Conclusions and results

After semi-analytical solutions and initial numerical calculations it appears that the pressure obtained directly for the oil dynamic viscosity variations crosswise the film thickness caused by the temperature gradients in gap height directions are about 5 to 7 percent different in comparison with the pressure values calculated for constant oil viscosity across the film thickness

*i*.*e*. in the case when the temperature is constant in gap height direction and varies only in bearing length and circumference direction.The mutually direct interactions of influences between hydrodynamic pressure on temperature as well as temperature on pressure, are valid only if the temperature, hence oil dynamic viscosity, varies in the gap height direction.

It is worse to notice that constant temperature in the gap height direction does not correspond with 3D temperature field obtained directly from energy equation. The assumptions of constant temperature and viscosity crosswise the bearing gap are in contradiction with the contemporary achievements connected with new devices such as micro-bearing, nano-bearing, magnetic bearings, artificial joints in humanoid robots, micro-motors. Unfortunately, numerous authors in the field of hydrodynamic lubrication, avoid to assume the oil viscosity variations in the gap height direction.

The performed numerical calculations show that oil, presented in the calculation, increases the pressure and load capacity in the slide journal bearing for Newtonian and non-Newtonian oil, in the presence of magnetic induction field.

## References

[1] Khonsari M.M., Wang S.H., On the Maximum Temperature in Double Layered Journal Bearing, Journal of Tribology, 1991, 113, 464-469.10.1115/1.2920647Search in Google Scholar

[2] Pincus O., Wilcock D.J., Thermal Effects in Fluid Film Bearings, Proceeding of the 6th Leeds-Lyon Symposium on Tribology: Thermal effects on tribology, 1980, 3-23.Search in Google Scholar

[3] Massoudi M., Local non-similarity solutions for the flow of a non-Newtnian fluid over wedge, International Journal of Non-Linear Mechanics, 2001, 36, 961-976.10.1016/S0020-7462(00)00061-5Search in Google Scholar

[4] Makinde O.D., On thermal stability of a reactive third-grade fluid in a channel with convective coolingthe walls, Applied Mathematics and Computation, 2009, 213, 170-176.10.1016/j.amc.2009.03.003Search in Google Scholar

[5] Schweitzer G., Traxler A., Bleuer H., Magnetlager, Springer Verlag, Berlin, Heidelberg, 199310.1007/978-3-662-08448-9Search in Google Scholar

**Received:**2017-11-02

**Accepted:**2017-11-12

**Published Online:**2018-05-30

© 2018 K. Wierzcholski and A. Miszczak, published by De Gruyter

This work is licensed under the Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 License.