Streamlines in the Two-Dimensional Spreading of a Thin Fluid Film: Blowing and Suction Velocity Proportional to the Height ()
1. Introduction
Thin fluid films occur in nature and they also play a significant role in many technological processes such as in coating applications and in the spreading of paints. There are important applications of thin fluid films to the lung [1] and the eye [2]. Reviews of the research up to the year 2000 on the spreading of thin fluid films and liquid drops have been given by Oron et al. [3] and by Davis [4].
Thin fluid flows usually occur under the action of gravity [5], surface tension gradients [4] or on a rotating substrate [6] [7]. The moving contact line at the boundary of the thin fluid film has been studied by O’Brien and Schwartz [8] and by Hocking [9]. Davis and Hocking [10] [11] have considered the spreading of a liquid drop on a porous base. The spreading of an axisymmetric liquid drop with fluid injection or suction at the base has been investigated by Mason and Momoniat [12] who later extended their work to include surface tension [13].
A study of the streamlines of the flow in a thin film can illustrate how the thin fluid film spreads. Momoniat et al. [14] investigated the effect of slot injection or suction on the spreading of a thin film where the surface tension gradient and gravity effects are included. They also investigated the behaviour of the streamlines and how they evolve with suction or blowing. They found that a breakup of the streamlines occurs. Mason and Chung [15] considered the spreading of a thin liquid drop with slip at the base but without fluid injection or suction at the base. The streamlines made clear the properties of the flow. It was found that vortices can occur at the base and close to the axis of the liquid drop.
The motivation for this research is to obtain an improvement on our understanding of the fluid flow inside a spreading thin fluid film with suction or blowing at the base by analysing the streamlines of the flow.
In general a nonlinear diffusion equation for the height of the thin film occurs in the mathematical formulation of thin fluid film flows. Exact analytical solutions for these equations have been found in the form of similarity solutions [16] [17]. The application of Lie group analysis of differential equations has been successful in solving problems in thin fluid film theory [18] [19] [20] and this method will be used here to reduce the nonlinear diffusion equation obtained to an ordinary differential equation. The half-width of the base is obtained during this reduction.
In this paper we investigate the effect of suction and blowing at the base on the spreading under gravity of a two-dimensional thin fluid film. The general form of the suction/blowing velocity at the base is determined from the condition that the partial differential equation and boundary conditions admit an invariant solution in a Lie symmetry analysis. We consider the case in which the suction/blowing velocity at the base is proportional to the height of the thin fluid film. The streamlines of the fluid flow inside the thin film are derived and analysed.
A nomenclature is provided in Table 1.
Table 1. Nomenclature for significant variables and parameters.
2. Mathematical Model
Consider the two-dimensional spreading under gravity of a thin film of viscous incompressible fluid on a fixed horizontal plane. At the base there is either suction/leak-off of fluid from the thin film or blowing/injection of fluid into the thin film. The model is illustrated in Figure 1. The derivation of the thin fluid film equations is well documented in the literature [3] [8]. We follow the derivation of Acheson [21].
The thin film is infinite in the y-direction and symmetric in the
-plane. The equation of the free surface is
(2.1)
where
(2.2)
and the width of the thin film is
(2.3)
The flow is two-dimensional in the
-plane with velocity components and pressure
(2.4)
Figure 1. Two-dimensional thin fluid film showing coordinate system, body force and boundary conditions.
where
(2.5)
The x and z components of the Navier-Stokes equation and the conservation of mass equation are [21]
(2.6)
(2.7)
(2.8)
where g is the acceleration due to gravity.
We introduce characteristic quantities and dimensionless variables and write (2.6) to (2.8) in dimensionless form:
characteristic length in the x-direction
,
characteristic length in the z-direction
,
characteristic velocity in the x-direction U,
characteristic velocity in the z-direction
,
characteristic pressure and stress P,
characteristic half-width
.
The characteristic velocity V is determined from the conservation of mass Equation (2.8) which is not approximated so that the two terms balance while U and P are yet to be specified. The Reynolds number is
(2.9)
where
is the kinematic viscosity. Dimensionless variables are defined by
(2.10)
The x-component of the Navier-Stokes Equation (2.6), becomes
(2.11)
For the thin fluid film to spread in the x-direction the magnitude of the pressure gradient in the x-direction should be sufficiently large to balance the viscous term. Thus
(2.12)
The z-component of the Navier-Stokes Equation (2.7), in dimensionless form is
(2.13)
where (2.12) was used for P in all terms except the body force term. For the thin fluid film to spread the body force term must be sufficiently large to balance the pressure gradient term in the z-direction. Thus we obtain the second condition on P,
(2.14)
Equating (2.12) and (2.14) for P gives for the characteristic velocity U,
(2.15)
The thin fluid approximation is [21]
(2.16)
Equations (2.11) and (2.13) reduce to
(2.17)
(2.18)
and (2.8) expressed in dimensionless variables is
(2.19)
Before stating the boundary conditions we derive two preliminary results. From the Navier-Poisson law for an incompressible fluid [21]
(2.20)
where
and
are components of the Cauchy stress tensor. Expressed in dimensionless variables,
(2.21)
and neglecting terms of order
,
(2.22)
Also, a fluid particle on the free surface of the thin film remains on the free surface as the fluid evolves. Thus
(2.23)
where
is the convective time derivative. Hence
(2.24)
and expressed in terms of dimensionless variables,
(2.25)
We now state the boundary conditions and then comment on them briefly.
(2.26)
(2.27)
(2.28)
(2.29)
(2.30)
Boundary condition (2.26) is the no slip condition for a viscous fluid at a solid boundary. Boundary condition (2.27) balances the normal velocity induced in the thin fluid film with the velocity of fluid injection or suction at the base where
The characteristic value for
is the characteristic velocity in the z-direction, V, and
(2.31)
The boundary condition (2.28) is the balance of the normal stress at the surface of the thin film where
is the atmospheric pressure and the boundary condition (2.29) is that the tangential stress vanishes at the surface.
Consider now the partial differential equation for the surface
which is derived from (2.30). In order to obtain
, we integrate the conservation of mass Equation (2.19) with respect to
from
to
and use the boundary condition (2.27). This gives
(2.32)
But using the formula for differentiation under the integral sign [22],
(2.33)
and therefore (2.32) becomes
(2.34)
Substituting (2.34) into (2.30) we obtain
(2.35)
It remains to obtain
. Integrating (2.18) and imposing the boundary condition (2.28) gives
(2.36)
and Equation (2.17) becomes
(2.37)
Integrating (2.37) twice with respect to
and imposing the no slip boundary condition (2.26) and the no tangential stress boundary condition (2.29) we obtain
(2.38)
Finally substituting (2.38) into (2.35) and integrating we find that
(2.39)
Equation (2.39) is a nonlinear diffusion equation for
with a source/sink term
. The velocity
is not specified but is determined from the condition that the problem admits an invariant solution.
Since the height of the thin fluid film vanishes at the moving contact lines,
, it follows that
(2.40)
The volume flux of fluid in the x-direction per unit breadth in the y-direction is
(2.41)
The characteristic volume flux is UH. Expressed in dimensionless variables and using (2.38),
(2.42)
Consider now the balance law for fluid volume. The total volume of the thin fluid film per unit length in the y-direction is
(2.43)
where the characteristic volume is
. Applying the formula for differentiation under the integral sign [22]
(2.44)
Imposing the boundary conditions (2.40) and substituting the nonlinear diffusion Equation (2.39), (2.44) becomes
(2.45)
At the moving contact lines,
, we take as boundary condition that the volume flux of fluid in the x-direction vanishes, that is
(2.46)
The balance law for fluid volume reduces to
(2.47)
The problem is to solve the nonlinear diffusion Equation (2.39) for
subject to the boundary conditions (2.40) and (2.46) and to the requirement that (2.39) admits a Lie point symmetry which generates an invariant solution.
In the remainder of the paper the overhead bars are suppressed to keep the notation simple, it being understood that dimensionless quantities are used.
3. Lie Point Symmetries and Group Invariant Solutions
When deriving the Lie point symmetries of the nonlinear diffusion Equation (2.39) the variables
and all the partial derivatives of h with respect to t and x are regarded as independent variables. A suffix is used to denote the partial derivatives of
. Equation (2.39) can be written in the form
(3.1)
where
(3.2)
The Lie point symmetry generators of (2.39) are of the form [23]
(3.3)
They are derived from the determining equation
(3.4)
where the second prolongation
of X is defined by
(3.5)
where
(3.6)
(3.7)
with summation over the repeated index k from 1 to 2. The total derivative operators with respect to t and x are
(3.8)
(3.9)
The determining Equation (3.4) is separated by powers and products of the partial derivatives of h. It is found that
(3.10)
where
(3.11)
provided the velocity
satisfies the first order linear partial differential equation
(3.12)
The velocity
therefore cannot be arbitrary for the Lie point symmetries of Equation (2.39) to exist.
Now,
is a group invariant solution of the nonlinear diffusion Equation (2.39) provided
(3.13)
that is, provided
satisfies the first order linear partial differential equation
(3.14)
We first consider the general case in which
,
and
and then consider the special case
,
and
.
3.1. General Case
The differential equations of the characteristic curves of (3.14) are
(3.15)
Integrating of the first pair of terms in (3.15) gives
(3.16)
where
is a constant while integration of the first and last pair of terms gives
(3.17)
where
is a constant. The general solution of (3.14) is
(3.18)
where F is an arbitrary function. Hence, since
, the form of the invariant solution is
(3.19)
where
(3.20)
Consider next the partial differential Equation (3.12) for
. We assume in addition that
. The differential equations of the characteristic curves of (3.12) are
(3.21)
Integration of the first pair of terms gives again (3.16) and integrating the first and last pair of terms gives
(3.22)
where
is a constant. The general solution of (3.12) is
(3.23)
where G is an arbitrary function. Hence the form of the invariant solution is
(3.24)
where
is given by (3.20).
Substituting (3.19) and (3.24) into the nonlinear diffusion Equation (2.39) reduces (2.39) to the ordinary differential equation
(3.25)
Equation (3.25) depends only on the ratio
and does not depend on
and
.
The half-width of the base,
, is determined from the boundary conditions in (2.40). Using (3.19) for
, (2.40) becomes
(3.26)
where
(3.27)
Differentiating (3.26) with respect to t gives
(3.28)
and therefore since
is not a constant function
(3.29)
Thus
and
where
and
are constants and
(3.30)
Adding we obtain
(3.31)
and differentiating (3.31) with respect to t, we find that
(3.32)
Thus
and from (3.31),
. Equation (3.30) becomes
(3.33)
and therefore since from (2.10),
,
(3.34)
The boundary condition (3.26) becomes
(3.35)
and from (3.20) with
and (3.33),
(3.36)
The total volume of the thin fluid film per unit length in the y-direction is
(3.37)
By using (3.19) and (3.36) it can be expressed as
(3.38)
where
(3.39)
The balance law for fluid volume is given by (2.47) and using (3.24) for
it becomes
(3.40)
In order to simplify the ordinary differential Equation (3.25) and the boundary conditions (3.35) we make the change of variables
(3.41)
where
(3.42)
Equations (3.25) and (3.35) become
(3.43)
and
(3.44)
The other variables transform to
(3.45)
(3.46)
(3.47)
(3.48)
where
(3.49)
There is one condition still to be imposed. Since the characteristic length in the z-direction is the initial height of the thin film at the central line
(3.50)
Hence from (3.45),
(3.51)
The ratio
occurs in the differential Equation (3.43) and balance law (3.47) and is either given or determined from the balance law. Let
(3.52)
then
(3.53)
The problem expressed in terms of the transformed variable
and the parameter
is to solve the differential Equation (3.43),
(3.54)
subject to the boundary conditions (3.44) and (2.46),
(3.55)
(3.56)
and the balance law
(3.57)
The remaining quantities are given by
(3.58)
(3.59)
(3.60)
(3.61)
where
(3.62)
and
is given by (3.49). The Lie point symmetry (3.10) which generates the invariant solution is
(3.63)
3.2. Special Case
and
When
the Lie point symmetry (3.10) reduces to
(3.64)
Now
is a group invariant solution of (2.39) provided (3.13) is satisfied, that is, provided
(3.65)
The differential equations of the characteristic curves of (3.65) are
(3.66)
The first pair of terms gives
(3.67)
while the pair consisting of the first and last terms give
(3.68)
where
and
are constants. The general solution of (3.65) is
(3.69)
where F is an arbitrary function. Hence since
,
(3.70)
where
(3.71)
When
, equation (3.12) reduces to
(3.72)
which has the same form as (3.65). Hence
(3.73)
where
is an arbitrary function of
.
Substituting (3.70) and (3.73) into the partial differential Equation (2.39) reduces it to the ordinary differential equation
(3.74)
The half-width of the base,
, is again obtained from the boundary conditions (2.40) which can be written as (3.26) with
(3.75)
It can be verified as for the general case that
and
where
and
are constants with
and that
. Since
we find that
(3.76)
The similarity variable (3.71) becomes
(3.77)
and the boundary conditions (3.26) are
(3.78)
The total volume (3.37), expressed in terms of the invariant solution, is
(3.79)
where
(3.80)
The balance law (2.47) becomes
(3.81)
The remaining condition,
, yields
(3.82)
In order to treat the special case
as the limit
of the general case we transform the special case by making the change of variables (3.41) where now
(3.83)
Expressed in terms of the transformed variables the problem is to solve the differential equation
(3.84)
subject to the boundary conditions (3.55) and (3.56) and to the balance law
(3.85)
and to the scaling condition
(3.86)
The remaining quantities are
(3.87)
(3.88)
(3.89)
(3.90)
where
(3.91)
and
(3.92)
The Lie point symmetry which generates the invariant solution is
(3.93)
Unlike the Lie point symmetry (3.63) which generates the general solution, (3.93) is not a scaling symmetry.
The results (3.84) to (3.93) can be obtained from (3.54) to (3.63) by writing
(3.94)
and using the expansion
(3.95)
for
we obtain
(3.96)
and therefore
(3.97)
4. Normal Velocity at Base Proportional to the Thin Fluid Film Height
An assumption on
needs to be made to close the system of equations. In this section we will investigate the solution for which
is proportional to
.
4.1. Invariant Solution
We will first consider the general case
and then the special case
. When
suppose that
(4.1)
where
is a constant. Then from (3.58) and (3.59)
(4.2)
Thus
is proportional to
and it is greatest at the centre line and vanishes at the moving contact lines,
.
The balance law (3.57) becomes
(4.3)
and therefore
(4.4)
The value
and its relation to the special case
will be considered later.
The differential Equation (3.54) reduces to
(4.5)
The solution of (4.5) which satisfies the boundary conditions (3.55) and (3.56) is
(4.6)
Hence using (3.58) to (3.62),
(4.7)
(4.8)
(4.9)
(4.10)
(4.11)
where
is the Gamma function [22] and
(4.12)
The Lie point symmetry (3.63) which generates the invariant solution becomes
(4.13)
Consider now the special case
and that (4.1) is satisfied. Then from (3.87) and (3.88),
(4.14)
and therefore
is proportional to
. Equation (4.14) agrees with (4.2) in the limit
. From balance law (3.85),
(4.15)
and therefore
. The differential Equation (3.84) reduces to (4.5) and the solution which satisfies the boundary conditions (3.55) and (3.56) is (4.6). Hence from (3.87) to (3.92) we obtain
(4.16)
(4.17)
(4.18)
(4.19)
where
is given by (4.11) and
(4.20)
The Lie point symmetry which generates the invariant solution, (3.93), becomes
(4.21)
In the same way as (4.17) was derived it can be shown that
(4.22)
and it can be verified that (4.16) to (4.20) can be obtained from the general results (4.7) to (4.12) as
.
4.2. Time Evolution of the Invariant Solution
We now investigate the evolution of the thin fluid film with time. The behaviour of the fluid variables depends on the value of the constant
where
.
We first consider the second condition in the thin fluid film approximation (2.16) and determine the range of values of
for which it is satisfied. Return to dimensional variables and define
(4.23)
where
(4.24)
We assume that the thin fluid film approximation is satisfied initially so that
(4.25)
Hence the thin fluid film condition will be satisfied as the fluid evolves if either
remains constant (
) or
as
or as
where
, is some limiting time. If
as
or as
then the thin fluid film approximation will break down.
Return to dimensional variables. From (4.7) and (4.9) for
,
(4.26)
and from (4.16) and (4.18) for
,
(4.27)
For
,
algebraically as
, for
,
exponentially as
while for
,
in the finite time
(4.28)
For
,
remains constant at
while for
,
in the finite time
. The thin fluid film approximation is therefore not satisfied for
and the solution is not valid. For
there is fluid leak-off or suction at the base and the fluid film remains thin as it evolves. For
there is fluid injection or blowing at the base and as long as the blowing is not too strong (
) the thin fluid film condition is satisfied. In the following we will consider only the range
.
Consider next the evolution with time of the half-width of the base,
, given by (4.9) for
and by (4.18) for
. Now
(4.29)
and therefore
for
. The base therefore always expands even for
which describes suction of fluid from the thin film. For
,
algebraically as
,
exponentially as
for
and for
,
in the finite time
defined by (4.28).
The total volume of the thin film per unit length in the y-direction,
, can be analysed in the same way. The evolution of
,
and
is summarised in Table 2.
The remaining two quantities,
and
have additional values of
to analyse. Consider first
which is given by (4.7) for
and by (4.16) for
. For all
,
(4.30)
and
(4.31)
Hence in the vicinity of the moving contact lines,
, the thin fluid film approximation is no longer satisfied. The inclusion of blowing or suction does not eliminate the singularity at
. Also, for
,
Table 2. Time evolution of
,
and
for
. The finite time
is defined by (4.28). The thin fluid film approximation is not satisfied for
.
(4.32)
For
,
and therefore the height of the centre line of the thin fluid film decreases with time even when there is blowing at the base,
. The rate of decrease in the height due to spreading is greater than the rate of increase due to blowing. When
the two rates of change balance and when
the rate of increase due to blowing is greater than the rate of decrease due to spreading. For
,
algebraically as
. For
,
remains constant, while for
,
algebraically as
. For
,
exponentially and for
,
in the finite time
given by (4.28). The time evolution of
is summarised in Table 3.
Finally
which describes suction and blowing at the base, is given by (4.8) for
and by (4.17) for
. Now for
,
(4.33)
For
,
which describes leak-off or suction at the base. For this range of
,
and hence the magnitude of the suction decreases as t increases and
algebraically as
. When
the base is impermeable. When
,
and there is fluid injection or blowing of fluid into the thin film at the base. For
,
Table 3. Time evolution of
for
. The finite time
is defined by (4.28).
so that the blowing becomes weaker as time increases and
algebraically as
. For
,
is constant and the blowing is constant in time. For
,
and the blowing becomes stronger as time increases. For
,
algebraically as
, for
,
exponentially as
and for
,
in the finite time
defined by (4.28). The results for the time evolution of
are summarised in Table 4.
In Figure 2 the evolution of
with time for a range of values of
is presented. Similar figures are obtained for the other fluid variables.
4.3. Stagnation Points on the Centre Line
In this section we investigate the z-component of the fluid velocity on the centre line,
. Since from symmetry
, points on the centre line at which
are stagnation points.
To obtain
we integrate the continuity Equation (2.19) with respect to z from 0 to z and impose the boundary condition (2.27) which gives
Table 4. Time evolution of
for
. The finite time
is defined by (4.28).
Figure 2. The maximum height of the thin fluid,
, plotted against t for a range of values of
.
(4.34)
By using (2.38) for
we obtain
(4.35)
and therefore on the centre line
,
(4.36)
But from (4.7) and (4.9) for
and from (4.16) and (4.18) for
(4.37)
where
(4.38)
and on the centre line,
. Thus
(4.39)
(4.40)
It remains to obtain
. From (4.8) and (4.9) for
and (4.17) and (4.18) for
,
(4.41)
and therefore using (4.38)
(4.42)
Substituting (4.39), (4.40) and (4.42) into (4.36) gives
(4.43)
where
(4.44)
The points on the centre line
at which
satisfy the cubic equation
(4.45)
The roots of the cubic Equation (4.45) which satisfy
are stagnation points of the flow.
We first transform (4.45) to standard form [24]
(4.46)
by making the transformation
(4.47)
The cubic Equation (4.45) becomes
(4.48)
Comparison of (4.48) with (4.46) gives
(4.49)
The roots of (4.48) which satisfy
are stagnation points of the flow. The discriminant
of the cubic Equation (4.48) is
(4.50)
The roots of (4.48) satisfy the following general properties [24].
The thin fluid film approximation breaks down for
.
The real zero for
does not lie in the range
. To show this let
(4.51)
Consider first
. Then
and
. The real zero therefore lies in the range
. Consider next
. Then
and
. The real zero therefore lies in range
. We therefore consider the real roots of (4.48) in the range
for which
. We use the trigonometric solution of the cubic equation when
because the Cardan method gives the solution in complex form even although the roots are real [24].
We look for a solution of the form
(4.52)
Substituting (4.52) into (4.48) gives
(4.53)
Comparing (4.53) with the trigonometric identity
(4.54)
we find that
(4.55)
The condition
is satisfied because
(4.56)
Let
be the smallest non-negative angle which satisfies
(4.57)
and
, that is,
. The general solution for
is
(4.58)
Hence the three real roots of (4.48) for
are, from (4.52) and (4.55),
(4.59)
and since
and
the three real roots of the cubic Equation (4.45) are
(4.60)
4.3.1. Impermeable Base
From (4.57),
and
since
. Thus
(4.61)
which may also be derived directly from (4.45) with
. The root
lies outside the thin fluid film. That (0,0) is a stagnation point follows directly from the boundary conditions.
4.3.2. Weak Blowing
We will derive expansions for
and
for small values of
. Since
when
we expect
to be small and positive when
is small. From (4.57)
(4.62)
and neglecting terms
we obtain form (4.62) the quadratic equation for
.
(4.63)
Hence
(4.64)
From (4.60) with
,
(4.65)
and
(4.66)
Thus
(4.67)
which, although decreasing in magnitude lies outside the thin fluid film.
From (4.60) with
,
(4.68)
But
(4.69)
and substituting (4.66) and (4.69) into (4.68) we obtain
(4.70)
The root
is negative and lies outside the thin fluid film.
From (4.60) with
,
(4.71)
and substituting (4.66) and (4.69) into (4.71) we obtain
(4.72)
The root
lies inside the thin fluid film and the point
is a stagnation point on the centre line
.
In Table 5 the approximate solution (4.72) is compared with the exact solution (4.71) where
satisfies (4.57) for
.
The expansions (4.67), (4.70) and (4.72) can also be derived using perturbation methods [25]. The solution (4.67) is a regular perturbation expansion while (4.70) and (4.72) are singular perturbation expansions.
4.3.3. Moderate Blowing
From (4.57) with
,
and since
it follows that
. Hence from (4.60)
(4.73)
The roots
and
lie outside the thin fluid film. The root
is the limiting value of the root lying inside the thin film. For
, the root
and lies outside the thin film. The point
is the limiting case of a stagnation point on the centre line. The maximum height of the thin film remains constant with time. This is consistent with the results of Table 3.
4.3.4. Strong Blowing
When
, from (4.57),
and since
we have
. Thus from (4.60),
(4.74)
which can be derived directly from (4.45) with
. There are no roots inside the thin fluid film for
.
The stagnation point
on the centre line plays a significant part in the pattern of the streamlines in the thin film which we now consider.
Table 5. The exact and approximate solution for the zero,
, of
on the centre line.
4.4. Streamlines
In order to study further the fluid flow in the thin film we investigate the streamlines of the flow. The tangent vector to a streamline is everywhere parallel to the fluid velocity vector instantaneously. In two-dimensional incompressible flow the streamlines can be obtained from the stream function
which satisfies
(4.75)
Since the stream function is constant along a streamline the streamlines in the
plane at time t can be obtained by plotting
(4.76)
for a range of values of
.
By substituting (2.38) and (4.36) for
and
into (4.75) we obtain
(4.77)
(4.78)
It is readily verified that the compatibility condition
(4.79)
is satisfied. Integrating (4.77) yields
(4.80)
and by substituting (4.80) into (4.78) we find that
(4.81)
Hence
(4.82)
where
. Thus
(4.83)
where
(4.84)
The streamlines in the
plane at time t are given by (4.76). They are obtained by solving numerically for z as a function of x the cubic equation
(4.85)
at a given time t and a range of values of
. The cubic Equation (4.85) for the streamlines is quite general and does not depend on any relation between
and
.
Consider now the invariant solutions (4.7) to (4.9) and (4.16) to (4.18). For both
and
,
(4.86)
(4.87)
(4.88)
The range of the base is
.
The streamlines derived from (4.85) depend on time. To determine the evolution of the streamline pattern with time it would have to be plotted for a range of values of time. Instead we will plot the streamlines in the
plane where
(4.89)
which we will see are independent of time and depend only on the parameter
. Clearly the base and maximum height are independent of time. The results for the evolution in time of
and
derived in Section 4.2 can then be used to determine how the streamline pattern evolves with time.
The cubic Equation (4.85) expressed in terms of
and
is
(4.90)
where using (4.86) to (4.88) and also (4.38)
(4.91)
(4.92)
(4.93)
Equations (4.91) to (4.93) are satisfied for both
and
. The coefficients
and
have the same time dependent factor
. Hence the cubic Equation (4.90) may be written as
(4.94)
where
(4.95)
The streamlines are obtained by giving
a range of values. Since (4.94) does not depend on time explictly the streamline pattern in the
plane depends only on the parameter
and is independent of time.
In Figure 3 the streamlines are plotted in the
plane for an impermeable base with
and for fluid leak-off (suction) at the base with
. For
the spreading is due only to gravity. The streamlines begin and end on the free surface except the streamline
which is the centre line. The point (0,0) is a stagnation point because of the boundary conditions. For fluid leak-off or suction the streamlines begin on the free surface and either end on the free surface or on the base. For stronger suction described by more negative values of
the streamlines are less curved and more streamlines end on the base. The streamlines are always perpendicular to the base because there is no slip at the base. There are no stagnation points in the flow.
In Figure 4 the streamlines are plotted in the
plane for values of
in the range
which describe fluid injection into the thin film at the base. There is a stagnation point on the centre line at
where
(4.96)
where
is given by (4.57). Dividing streamlines emanate from the stagnation point
in the positive and negative
-directions and divide the flow into two regions. In the upper region the streamlines begin and end on the free surface except the streamline along the centre line which ends at the stagnation point. The flow is similar to that for
and is dominated by spreading due to gravity. In the lower region the flow is dominated by fluid injection at the base. This is especially the case near the centre line and for small values of
but as
increases spreading due to gravity causes the streamlines to bend towards the free surface. The maximum thickness of the lower region is
. As the strength of the blowing,
, increases we see that the thickness of the lower region increases. It forms a larger fraction of the thin film in agreement with the
Figure 3. Streamlines for (a) impermeable base with
and (b) leak-off at the base with
, plotted in the
plane.
Figure 4. Streamlines for (a)
, (b)
, (c)
, (d)
, (e)
and (f)
, plotted in the
plane.
numerical vales in Table 5. When
increases to
the lower region completely fills the thin film and there is no longer two regions. When
, the fluid velocity on the centre line is positive except at
where
.
In Figure 5 the streamlines are plotted in the
plane for
,
and for the limit of the thin fluid film approximation,
. The streamline pattern is similar to that for
. The streamlines emanate from the base at right angles because of the no slip boundary condition. Spreading due to gravity causes them to bend towards the free surface. The blowing is stronger
Figure 5. Streamlines for (a)
, (b)
and (c)
, plotted in the
plane.
as
increases and the bending occurs at larger values of
. On the centre line
for
and there is no stagnation point in the flow.
The evolution of the streamline pattern with time can be determined using
(4.97)
and the results in Table 2 and Table 3 for the evolution of
,
and
with time.
In Figure 3,
and
and
as
. The base will steadily increase and the maximum height will steadily decrease for
.
In Figure 4,
and
as
. When
,
as
but for
,
remains constant and the point of maximum height is a stagnation point. The height of the stagnation points above the base on the centre line,
(4.98)
will decrease with time for
and remain constant for
. The maximum thickness of the lower and upper regions will both decrease with time but the ratio of the maximum thicknesses will remain constant as time increases:
Figure 6. Streamlines for
for (a)
and (b)
, plotted in the
-plane.
(4.99)
For
, the maximum height
will remain constant but the base
will expand to
as
.
In Figure 5,
and
and
algebraically as
for
, exponentially for
and in the finite time
defined by (4.28) for
. The thin fluid film ratio
as
or
for
and therefore the base expands to infinity faster than the maximum height but for
,
and
and
tend to infinity at the same rate. For
the thin fluid film approximation is not satisfied for all time.
The evolution of the streamlines with time is illustrated in Figure 6 where the streamlines are plotted in the
-plane for
at times
and
by solving numerically the cubic Equation (4.85) where
,
and
are given by (4.86) to (4.88).
5. Conclusions
We first list the findings of the paper and then comment on the findings.
The thin fluid film equations were derived and it was shown how the characteristic fluid pressure could be determined from the dimensionless partial differential equations and how the characteristic horizontal fluid velocity could be determined from the two expressions obtained for the characteristic fluid pressure.
It was found that the height of the thin fluid film satisfied a nonlinear diffusion equation with a source/sink term consisting of the fluid injection/leak-off velocity,
, at the base.
The nonlinear diffusion equation for the height of the thin fluid film was found to possess a Lie point symmetry provided
which satisfies a first order linear partial differential equation. This leads to the invariant form for
.
The general form of the Lie point symmetry generates an invariant solution with height and half-width which evolve algebraically with time. A special form of the Lie point symmetry was also found which generates a solution for which the height and half-width evolve exponentially with time.
To close the system of equations, we found that one further condition is required. We assumed that
is proportional to the height
.
The time evolution of the invariant solution for the thin fluid film ratio
, the height
, the total volume per unit breadth
and
was determined for both the general and special solution.
The stagnation point on the centre line was found by solving a cubic equation. At this point the upward flow due to fluid injection at the base is balanced by the downward flow due to spreading by gravity.
The fluid flow inside the thin fluid film was investigated by plotting the streamlines which were obtained from the stream function by solving numerically a second cubic equation. Streamlines for both fluid leak-off and fluid injection were plotted. For fluid injection, it was found that, provided the injection velocity is not too strong, there is a dividing streamline which separates the flow into two regions, an upper region with downflow due to gravity and a lower region with upflow due to fluid injection at the base.
We make the following comments on the findings.
When the fluid velocity at the base is proportional to the height of the thin fluid film, a complete analytical solution can be derived and the fluid variables and streamlines can be fully analysed. The thin fluid film approximation is satisfied for all values of time for the whole range of suction and significantly it is also satisfied for blowing at the base which provided the blowing that is not too strong.
There was unexpected behaviour of fluid variables due to the relative importance of spreading by gravity compared with blowing or suction at the base. The base half-width increased with time for all values of
, even for suction when
. The maximum height of the thin film decreased for all time for suction and even for blowing provided the blowing was not too strong (
). Some limits were attained in a finite time and in such a way that the thin fluid film approximation remained satisfied. For sufficiently strong blowing (
) the base half-width and the maximum height of the thin fluid film tended to infinity in time
.
Much of the literature on thin fluid films is concerned with how the surface profile evolves with time. By plotting the streamlines numerically we were able to investigate the fluid flow inside the thin film. The streamline pattern for blowing was much richer than that for suction which did not show any unexpected features. The depth of penetration of the injected fluid at the base and its effect on the fluid flow inside the thin film could be investigated. The streamline pattern showed the existence of a dividing streamline which separated the flow into two regions, a lower region at the base consisting of rising fluid and an upper region consisting of descending fluid. We can expect that this structure will exist for other models of fluid injection. A streamline pattern independent of time was obtained. This was useful because from it and the known time evolution of
and
the streamline pattern at any later time could be determined.
The general solution for
and the special solution for
where both generated by the Lie point symmetry (3.10). The general solution was generated from (3.10) with
and the special solution with
. This unified the two solutions and is the reason why the special solution could be obtained as a limiting case of the general solution. When
, the Lie point symmetry is a scaling symmetry while the special case
is not a scaling symmetry.
The theory developed in the paper up to the end of Section 3 is quite general. It is only in Section 4 that the assumption that
is proportional to
is made. Other models for
can be considered that lead to analytical solutions for the streamlines, for example,
proportional to the spatial gradient of the height which was considered in the axisymmetric spreading of a liquid drop [12].
Acknowledgements
DPM thanks the National Research Foundation, Pretoria, South Africa, for financial support under grant number 132189.