2.2 Lubrication Approximation For Flow in A Thin Layer
2.2 Lubrication Approximation For Flow in A Thin Layer
2.2 Lubrication Approximation For Flow in A Thin Layer
2-2lubri.tex
H/L ≪ 1 (2.2.1)
H ′
t = T t′ , x = Lx′ , y = Hy ′ , u = U u′ , v = U v, p = P p′ , (2.2.2)
L
The normalized continuity equation is
∂u′ ∂v ′
+ =0 (2.2.3)
∂x′ ∂y ′
Both terms are equally important, reflecting the holyness of the law of mass concentration.
The longitudinal momentum equation is normalized to
! !
U ∂u′ U 2 ′ ∂u
′
′ ∂u
′
P ∂p′ νU H 2 ∂ 2 u′ ∂ 2 v ′
+ u + v = g sin θ − + + (2.2.4)
T ∂t′ L ∂x′ ∂y ′ ρL ∂x′ H 2 L2 ∂x′2 ∂y ′2
νU
Dividing by H2
, we get
! !
H 2 ∂u U H H ∂u′ ∂u′ g sin θH 2 P H 2 ∂p′ H 2 ∂ 2 u′ ∂ 2 u′
+ u′ ′ + v′ ′ = − + + ′2 (2.2.5)
νT ∂t ν L ∂x ∂y νU ρLνU ∂x′ L2 ∂x′2 ∂y
UH
= O(1) (2.2.6)
ν
2
and
H2
≪1 (2.2.7)
νT
Omitting terms of the order H/L and smaller, the above equation can be approximated to
the leading order by
g sin θH 2 P H 2 ∂p′ ∂ 2 u′
0= − + (2.2.8)
νU ρLνU ∂x′ ∂y ′2
or in dimensional form,
1 ∂p ∂ 2u
0 = g sin θ − +ν 2 (2.2.9)
ρ ∂x ∂y
All inertia terms are inconsequential; the most important balance is among gravity, the
pressure gredient and the dominant viscous stress. This balance also implies a pressure
scale,
ρLνU
P = (2.2.10)
H2
From the transverse momentum equation,
" !# !
H U ∂v ′ U 2 ′ ∂v
′
′ ∂v
′
P ∂p′ H νU ∂ 2v′ ∂ 2v′
+ u + v = −g cos θ − + + (2.2.11)
L T ∂t′ L ∂x′ ∂y ′ ρH ∂y ′ L H2 ∂x′2 ∂y ′2
or
( !)
H2 H 2 ∂v ′ U H H ′ ∂v
′
′ ∂v
′
+ u + v =
L2 νT ∂t′ ν L ∂x′ ∂y ′
!
g sin θH 2 H P H 2 ∂p′ H 2 H 2 ∂ 2 v ′ ∂ 2 v ′
− − + 2 + (2.2.12)
νU L tan θ ρLνU ∂y ′ L L2 ∂x′2 ∂y ′2
H gH cos θ
U= (2.2.16)
L ν
3
The distinguishing feature of negligible inertia is shared by the slow flow through thin
gaps of bearings in the theory of lubrication. Hence (2.2.9) and (2.2.14) can be called the
lubrication approximation.
We leave it as an exercise to show by similar normalization, that the dynamic boundary
conditions on y = h can be approximated to the leading order by
∂u
=0 (2.2.17)
∂y
for the tangential stress, and
p=0 (2.2.18)
for the normal stress. It follows by integrating (2.2.14)) that
which can be inserted in the integrated mass conservation equation (2.1.12) to give
" !#
∂h ∂Q ∂h ρg ∂ ∂h
+ = + h3 sin θ − cos θ =0 (2.2.22)
∂t ∂x ∂t 3µ ∂x ∂x
This is a nonlinear diffusion equation governing the evolution of the fluid depth.
In the special limit of a uniform flow, ∂/∂x ≡ 0. The velocity profile is then
!
ρgh2 y y2
u= sin θ − 2 (2.2.23)
µ h 2h
ρgh3
Q= sin θ (2.2.24)
3µ