13.02.2013 Views

Mechanics of Fluids

Mechanics of Fluids

Mechanics of Fluids

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

236 Laminar flow between solid boundaries<br />

Similarly there is a net outflow (∂Qz/∂z)δz from the box in the z direction.<br />

If, at the particular position considered, h does not vary with time (i.e. if<br />

both the speed and the load are steady), the total net outflow <strong>of</strong> lubricant<br />

from the box must be zero.<br />

∴ ∂Qx<br />

∂x<br />

∂Qz<br />

δx + δz = 0 (6.59)<br />

∂z<br />

To obtain ∂Qx/∂x we make appropriate substitutions in eqn 6.23:<br />

� �<br />

∂Qx ∂<br />

= δz −<br />

∂x ∂x<br />

∂p h<br />

∂x<br />

3<br />

��<br />

�Rh<br />

+<br />

12µ 2<br />

�<br />

�<br />

�R ∂h 1 ∂ h<br />

= −<br />

2 ∂x 12 ∂x<br />

3<br />

��<br />

∂p<br />

δz (6.60)<br />

µ ∂x<br />

where the negligible changes <strong>of</strong> elevation allow us to drop the asterisk from<br />

p∗ and δz takes the place <strong>of</strong> b. The corresponding equation for the z (axial)<br />

direction, in which there is no relative velocity between journal and bush, is<br />

∂Qz<br />

∂z =−1<br />

�<br />

∂ h<br />

12 ∂z<br />

3<br />

�<br />

∂p<br />

δx (6.61)<br />

µ ∂z<br />

Substituting eqn 6.60 and 6.61 into eqn 6.59, we obtain<br />

�<br />

�R ∂h 1 ∂ h<br />

δxδz −<br />

2 ∂x 12 ∂x<br />

3<br />

�<br />

∂<br />

δxδz −<br />

µ ∂x<br />

1<br />

�<br />

∂ h<br />

12 ∂z<br />

3<br />

�<br />

∂p<br />

δxδz = 0<br />

µ ∂z<br />

If the axes <strong>of</strong> journal bush are parallel, h will be independent <strong>of</strong> z, and<br />

dependent only on x. With µ again assumed constant throughout, the<br />

equation thus becomes<br />

6µ�R dh<br />

�<br />

∂<br />

− h<br />

dx ∂x<br />

3 ∂p<br />

∂x<br />

�<br />

− h 3 ∂2p = 0 (6.62)<br />

∂z2 The position z = 0 may be taken mid-way between the ends <strong>of</strong> the bearing;<br />

for a given value <strong>of</strong> x, the pressure varies from, say, pc at z = 0 to atmospheric<br />

at each end <strong>of</strong> the bearing (where z =±L/2). That is, p/pc varies<br />

symmetrically along the length <strong>of</strong> the bearing from zero to 1 and back to zero.<br />

It is reasonable to assume that the form <strong>of</strong> this variation is the same whatever<br />

the value <strong>of</strong> x; that is, that f = p/pc is a function <strong>of</strong> z only. (We shall show<br />

later that for a short bearing this assumption is justified.) Substituting pc for<br />

p in eqn 6.62 gives<br />

6µ�R dh d<br />

−<br />

dx dx<br />

and division by h 3 pc then yields<br />

�<br />

h<br />

�<br />

f − h<br />

dx<br />

3 d<br />

pc<br />

2f = 0 (6.63)<br />

dz2 3 dpc<br />

a − b 2 f − d2f = 0 (6.64)<br />

dz2

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!