Syvum Home Page

Home > Print Preview

Heat Transfer Problem :
Maximum temperature in lubricant by viscous heating


Problem.

An oil (of viscosity symbol : mu and thermal conductivity k) acts as a lubricant between two coaxial cylinders. The inner cylinder is stationary and the outer cylinder of radius R rotates at an angular velocity symbol : Omega. The clearance between the cylinders is b, which is small compared to the radii of the cylinders; so, curvature effects can be neglected and the cylindrical system can be approximated by a plane narrow slit (to be solved in Cartesian coordinates) as shown in the figure. Derive an expression for the maximum temperature in the lubricant if both cylinders are at temperature T0. Neglect the temperature dependence of symbol : mu and k, but explicitly take into account the heat generated by viscous dissipation.

figure : coaxial cylinders approximated by plane narrow slit

Figure. Temperature profile for viscous heat generation. The rectangular section in the flow between two coaxial cylinders can be approximated by the plane narrow slit on neglecting the curvature of the bounding surfaces.


Solution.

The clearance between the cylinders is b, which is small compared to the radii of the cylinders; so, curvature effects are neglected and the cylindrical system is approximated by a plane narrow slit. From a shell momentum balance in Cartesian coordinates over a thin rectangular slab of thickness symbol : Deltax, the momentum flux (shear stress) distribution is obtained as

equation : d tau_{xz}/dx = Delta P / L = 0 implies tau_{xz} = K_1 (1)

The flow is solely due to the movement of the top surface and there is no pressure gradient imposed on the system.

On substituting Newton's law of viscosity (equation : tau_{xz} = -mu dv_z/dx), the expression for the velocity profile is

equation : dv_z/dx = -K_1/mu implies v_z = -K_1/mu z + K_2 (2)

The integration constants K1 and K2 are determined using the boundary conditions:

equation : BC 1: x = 0, v_z = 0 implies K_2 = 0 (3)

equation : BC 2: x = b, v_z = V implies V = - K_1/mu b (4)

On substituting the integration constants, the velocity profile is given by

equation : v_z / V = x / b (5)

During viscous flow, fluid layers rub against adjacent layers of fluid and the internal friction produces heat. In other words, mechanical energy is degraded into thermal energy irreversibly. The heat generation by viscous dissipation per unit volume depends on the local velocity gradient and is given by

equation : S = -tau_{xz} dv_z/dx = mu (dv_z/dx)^2 = mu (V/b)^2 (6)

For the system under consideration, S is found to be constant.

Now,

Rate of Heat           In - Out + Generation = Accumulation

At steady-state, the accumulation term is zero. From a thermal energy balance over a thin rectangular slab of thickness symbol : Deltax in the fluid, we get

equation : (A q_x)|_x - (A q_x)|_x+Delta_x + S (A Delta_x) = 0 (7)

where A is the cross-sectional area of the narrow slit and qx is the heat flux for conduction in the x-direction.

Dividing by A symbol : Deltax and taking the limit as symbol : Deltax tends to zero,

equation : lim Delta-x --> O [(q_x)|_x+Delta-x - (q_x)|_x]/Delta-x = S (8)

equation : dq_x/dx = S (9)

On integrating, the heat flux is given by

equation : q_x = S x + C_1 (10)

On substituting Fourier's law (equation : q_x = -k dT/dx), we get

equation : dT/dr = -S/k x - C_1/k (11)

On integrating, the temperature profile is given by

equation : T = -S/(2k) x^2 - C_1/k x + C_2 (12)

The integration constants C1 and C2 are determined using the boundary conditions:

equation : BC 1: x = 0, T = T_0 implies C_2 = T_0 (13)

equation : BC 2: x = b, T = T_0 implies -C_1/k = Sb/(2k) (14)

On substituting the integration constants, the temperature profile is

equation : T = T_0 + Sb^2/(2k) [x/b - (x/b)^2] (15)

The maximum temperature occurs when

equation : dT/dx = Sb^2/(2k) [1/b - 2x/b^2] = 0 implies x = b/2 (16)

Thus, the maximum temperature occurs at the mid-plane of the slit and is obtained as

equation : T_max = T_0 + Sb^2/(8k) = T_0 + mu V^2/(8k) = T_0 + mu (Omega R)^2 / (8k) (17)


Previous
-
Next
-

Contact Info © 1999-2024 Syvum Technologies Inc. Privacy Policy Disclaimer and Copyright