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 :
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/m2, D - electric induction vector As/m, ρ - fluid density kg/m3, κ - thermal conductivity coefficient W/mK, v – fluid velocity vector in m/s, φF – dissipation of energy in W/m3, S - stress tensor in the fluid in Pa, u - displacement vector of the solid body in m, t - time in s, cv - 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 2Td = A1 i.e. constitutive equations are assumed for the lubricant in the following form :
We assume: n – dimensionless flow index (0.5, 1.2), m = m(n) fluid consistency coefficient in Pasn.
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 ϵij take the following form:
for i,j = 1,2,3, where: δij – Kronecker unit tensor component (δij = 1 for i = j and δij = 0 for i≠j), K = Λ + 2G/3, G = EY/[2(1+ ϑ)], Λ = EY ϑ/(1+ϑ)(1–2ϑ), whereas EY - Young‘s Modulus in Pa, αT – thermal coefficient of linear expansionin K–1, ϑ – Poisson ratio.
And by virtue of known linear geometrical relations, the strain tensor components Θij(vi) in the oil, and strain components ϵij(ui) in the sleeve are mutually connected with the oil velocity components vi and displacement vector of solid body components ui, respectively:
2 Electro-magneto-hydrodynamic equations
We assume: unsymmetrical, incompressible (for invariant I1 = 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 ϵT between the two surfaces is signifucantly smaller than other dimensions of the considered surfaces. Taking into account the layer boundary simplification 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: h1 = h1(α3), h2 = 1, h3 = h3(α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π, –bm ≤ α3 bs, 0 ≤ α2 ≤ ϵT and 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 vi (α1,α2,α3) and three layer displacement components ui (α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, v2 = 0 for α2 = 0,
and for α2 = ϵT, v3 = 0 for α2 = 0, v3 = 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 T0 have constant value fc on the journal surface and variable unknown values fp (α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/m2K, qc - heat flux density in oil on journal surface in W/m2.
The elastic layer of bearing alloy is laying on the rigid ring in place α2 = c3, and therefore the contact surface of these bodies is not deformed by the pressure. The contact surface α2 = c3 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: T0 - average ambient temperature, qp - heat flux density across the elastic layer located on the sleeve surface, f2(α1, α3) - temperature increases above the temperature T0 on the external surface of the elastic layer for α2 = c2, additionally e2(α1,α3) temperature increases in the excess of temperature T0 on the bottom surface of the elastic layer in the place of α2 = c3.
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, δE of magnetic, electric influence on the oil viscosity and coefficient δT of temperature influence on the oil viscosity Mentioned phenomena are confirmed by virtue of authors own measurements.
5 Sketch of semi-analytical solutions
After double integration of two equations (15) in α2 direction we apply the low heat flux continuity (qp = qc) 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:
where ϵs = c3 – c2, 0 < α1 < 2π, –bm < α3 < +bm. Solution (25) valid if and only if EY - Young modulus and solid body coefficients αT, κ*, ϑ - Poisson ratio, are independent of α2, The hydrodynamic pressure p, the load carrying capacity and temperature T distributions with oil velocity components v1, v2, v3 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 with qc = 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(v1, v3) 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:
For k = 0 the abovementioned system of equations (28)-(31) determines functions: for Newtonian oil taking into account the boundary conditions in steady motion (18), (19). For 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π, –bm < α3 < +bm, whereas we denote:
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):
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) = fcδk0 for α2 = 0, T(k) = fpδk0 for (α2 = ϵT, and after the term ordering, we finally obtain:
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 fp (α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: for i = 1.2.3. In this case we obtain oil velocity components pressure p(0), temperature T(0). Assuming cylindrical coordinates: α1 = ϕ, α2 = r, α3 = z, h1 = R, h2 = h3 = 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(δBB + δEE).
For δ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.
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. Google Scholar
Schweitzer G., Traxler A., Bleuer H., Magnetlager, Springer Verlag, Berlin, Heidelberg, 1993 Google Scholar
About the article
Published Online: 2018-05-30
Citation Information: Open Physics, Volume 16, Issue 1, Pages 285–291, ISSN (Online) 2391-5471, DOI: https://doi.org/10.1515/phys-2018-0040.
© 2018 K. Wierzcholski and A. Miszczak, published by De Gruyter. This work is licensed under the Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 License. BY-NC-ND 4.0