Squeezing flows of Newtonian liquid films an analysis including fluid inertia

A theoretical study is made of the flow behavior of thin Newtonian liquid films being squeezed between two flat plates. Solutions to the problem are obtained by using a numerical method, which is found to be stable for all Reynolds numbers, aspect ratios, and grid sizes tested. Particular emphasis is placed on including in the analysis the inertial terms in the Navier-Stokes equations.Comparison of results from the numerical calculation with those from Ishizawa's perturbation solution is made. For the conditions considered here, it is found that the perturbation series is divergent, and that in general one must use a numerical technique to solve this problem.Nomenclature a half of the distance, or gap, between the two plates - a ₀ the value of a at time t=0 - adot da/dt - ä d² a/dt ² - d³ a/dt ³ - a ⁱ components of a contravariant acceleration vector - f unknown function of z ₀ and t defined in (6) - f ⁱ function defined in (9) f ¹=r ₀ g(z ₀, t) f ²= ₀ f ³=f(z ₀, t) - F force applied to the plates - g unknown function of z ₀ and t defined in (6) - g g/z ₀ - h grid dimension in the z ₀ direction (see Fig. 5) - Christoffel symbol - i, j, k, l indices - k grid dimension in the t direction (see Fig. 5) - r radial coordinate direction defined in Fig. 1 - r ₀ radial convected coordinate - R radius of the circular plates - t time - v _r fluid velocity in the r direction - v _z fluid velocity in the z direction - v fluid velocity in the direction - x ⁱ cylindrical coordinate x ¹=r x²= x³=z - z vertical coordinate direction defined in Fig. 1 - z ₀ vertical convected coordinate - tangential coordinate direction - ₀ tangential convected coordinate - viscosity - kinematic viscosity, / - ⁱ convected coordinate ¹=r₀ ²=₀ ³=z₀ - density