Streamlines in the Two-Dimensional Spreading of a Thin Fluid Film: Blowing and Suction Velocity Proportional to the Spatial Gradient of the Height ()
1. Introduction
In a recent paper [1] we investigated the effect of fluid leak-off (suction) and fluid injection (blowing) at the horizontal base on the two-dimensional spreading under the gravity of a thin film of viscous incompressible fluid. We derived and analysed the plots of the streamlines to understand the fluid flow in the thin film on which the evolution of the surface profile depended. The system of equations was closed by making the assumption that the normal component of the fluid velocity at the base, vn, is proportional to the height of the thin film at that point. With this assumption, an analytical solution for the surface profile could be derived and the streamlines were plotted by solving a cubic equation numerically. In this model, the magnitude of the suction and blowing is greatest near the centre line of the thin film because the height is greatest on the centre line and is least near the moving contact lines where the height vanishes.
In this paper, we will assume that vn is proportional to the spatial gradient of the height. An analytical solution for the surface profile can again be derived. We will plot the streamlines and investigate their properties and compare them with the results obtained in [1]. In this model, the magnitude of the suction and blowing is least near the centre line on which the spatial gradient vanishes and is greatest near the moving contact lines on which the magnitude of the spatial gradient is infinite.
Investigations of the evolution of a thin fluid film with a free surface generally concentrate on the surface profile [2] [3] [4]. In this paper, we investigate the properties and evolution of the thin film by deriving and analysing the streamlines of the fluid flow in the thin film.
Mason and Momoniat [5] made the assumption that vn is proportional to the spatial gradient of the height in the axisymmetric spreading under the gravity of a thin liquid drop with suction and blowing at the base. They did not investigate the streamlines in the liquid drop. We will use Lie group analysis of differential equations to reduce the nonlinear diffusion equation to an ordinary differential equation and derive an analytical solution. This is a powerful and systematic method that has been applied successfully in other investigations in thin fluid film theory. It can be used to analyse the evolution of a thin fluid film from a point or line source of fluid and from a fluid film with an initial non-zero half-width as considered here. A literature review of the relevant papers of thin fluid film theory and of the application of Lie group analysis to thin fluid films was given in [1].
An outline of the paper is as follows. In Section 2 a summary is given of the results needed from [1]. In Section 3 the invariant solution is derived from the Lie group analysis performed in [1]. In Section 4 the evolution with time of the fluid variables is analysed. In Section 5 the streamlines are derived and plotted and in Section 6 an approximate analytical solution for the dividing streamline is obtained. Finally, the conclusions are summarised in Section 7.
A nomenclature is presented in Table 1.
2. Summary of Results
The fluid is viscous and incompressible. The spreading of the thin fluid film is two-dimensional as shown in Figure 1. There is suction (fluid leak-off) or blowing (fluid injection) at the base. The moving contact lines are
and the height of the thin fluid film is
. The fluid film is infinite in the y-direction. The variables are dimensionless.
The thin fluid film approximation is imposed. The x- and z-components of the fluid velocity satisfy
(2.1)
(2.2)
Table 1. Nomenclature for significant variables and parameters.
Figure 1. Spreading of a two-dimensional thin fluid film under gravity illustrating the coordinate system, body force and boundary conditions.
The height
satisfies the nonlinear diffusion equation
(2.3)
where
is the normal component of the fluid velocity at the base. The boundary conditions are that at the moving contact lines the height vanishes,
(2.4)
and the fluid flux vanishes,
(2.5)
The balance law for fluid volume is
(2.6)
where
is the total volume of the thin film per unit length in the y-direction:
(2.7)
An invariant solution is derived using Lie group analysis for differential equations [1] [6]. The nonlinear diffusion Equation (2.3) admits the Lie point symmetry
(2.8)
provided
satisfies
(2.9)
Using the Lie point symmetry (2.8) the nonlinear diffusion Equation (2.3) is reduced to an ordinary differential equation. There are two cases, the general case,
and the special case,
.
When
the invariant solution is obtained by solving the ordinary differential equation
(2.10)
where
(2.11)
subject to the boundary conditions
(2.12)
(2.13)
and to the balance law for fluid volume,
(2.14)
The invariant solutions for the fluid variables are
(2.15)
(2.16)
(2.17)
(2.18)
where
(2.19)
The Lie point symmetry which generates the invariant solution is
(2.20)
when
the invariant solution is obtained by solving the ordinary differential equation
(2.21)
subject to the boundary conditions (2.12) and (2.13) and to the balance law for fluid volume
(2.22)
The invariant solution is of the form
(2.23)
(2.24)
(2.25)
(2.26)
where
is given by (2.19) and
is given by (2.11). The Lie point symmetry which generates the invariant solution is
(2.27)
Unlike (2.20), the Lie point symmetry (2.27) is not a scaling symmetry. It can be verified that the system of equations (2.21) to (2.26) for
can be obtained from the Equations (2.10) to (2.20) for
by letting
.
For the two-dimensional flow of an incompressible fluid the streamlines can be obtained from the stream function. The stream function
satisfies
(2.28)
and the streamlines at time t are the curves
(2.29)
where
is given a range of discrete values. In the
-plane the streamlines at timet are obtained by solving numerically for z as a function of x the cubic equation [1]
(2.30)
for a range of discrete values of
where
(2.31)
We will also consider the streamlines in the
-plane where
and
. In the
-plane the streamlines satisfy the cubic equation [1]
(3.32)
where
(3.33)
In order to close the systems of equations and complete the formulation of the problem an assumption needs to be made on
. In [1] it was assumed that
is proportional to
. In this paper it will be assumed that
is proportional to the spatial gradient
.
3. Invariant Solution
We first consider the general case
. Assume that
(3.1)
where
is a constant. The balance law (2.14) becomes
(3.2)
and by integrating by parts we find that
(3.3)
Hence
(3.4)
The differential Equation (2.10) with
given by (3.4) becomes
(3.5)
which can now be integrated subject to the boundary conditions (2.12) and (2.13). We obtain
(3.6)
and therefore
(3.7)
For
there is one real solution for
which is negative and two complex conjugate solutions which we do not consider. From (2.15) to (2.18) we see that the real negative solution gives a positive solution for
and that the other fluid variables depend only on
. Hence
with
.
In order to determine the physical significance of
consider
. From (2.16),
(3.8)
We introduce the new parameter
defined by
(3.9)
The factor
is a normalisation factor to give
for the special case
. The range of
is
. Equation (3.8) becomes
(3.10)
The normal fluid velocity at the base has the same sign as
and
describes blowing (fluid injection) into the thin film,
describes an impermeable base and
describes suction (fluid leak-off) at the base. The other fluid variables are
(3.11)
(3.12)
(3.13)
(3.14)
From (3.10), (3.11) and (3.14)
(3.15)
We see that the physical significance of the solution is that
is proportional to
which like
is an even function of x.
Consider now the special case,
. Assuming (3.1), the balance law (2.22) becomes
(3.16)
and therefore
and
The differential Equation (2.21) becomes
(3.17)
Integrating and imposing the boundary conditions we obtain
(3.18)
Hence from (2.23) to (2.26),
(3.19)
(3.20)
(3.21)
(3.22)
where
is given by (3.13).
4. Analysis of the Solution
We will determine the evolution in time of the fluid variables, especially
and
and investigate the fluid velocity on the centre line,
.
4.1. Time Evolution of the Fluid Variables
Consider first the thin fluid film approximation [1]. From (3.11) and (3.14) for
,
(4.1)
and from (3.19) and (3.21) for
,
(4.2)
For
,
as
where
(4.3)
for
,
exponentially as
while for
,
algebraically as
. For
,
and therefore remains constant and for
,
as
. Since it is assumed that the thin fluid film approximation is satisfied initially it will remain satisfied for all time for
which is the range of
for which the solution is valid physically. The special solution for
with exponential behaviour and the solution in the range
in which the limit is attained in finite time lie outside the physical range. Interestingly the thin film approximation is satisfied for the whole range of blowing,
, which also has the most interesting streamlines. The spreading due to gravity is the dominant mechanism and is stronger than the effect of blowing on the height. The results are summarised in Table 2.
Consider next the half-width of the base,
, for the range
in which the solution applies. For this range,
is given by (3.11). For
,
as
but in such a way that the thin fluid film approximation remains valid. The two-dimensional thin film shrinks to a line as
due to suction. For
,
remains constant as the fluid evolves and for
,
as
. The behaviour of
for the full range,
is presented in Table 2.
Consider next
for the physical range,
. From (3.14)
Table 2. Time evolution of
,
,
and
for
. The finite time
is defined by (4.3). The thin fluid film approximation is satisfied for all time only for
.
(4.4)
and
(4.5)
We see that
for
and that
as
,
even for blowing when
. The behaviour of
as t evolves for the full range,
, is given in Table 2.
The evolution of
is summarised in Table 2. The time evolution of
for a fixed value of
in the range
is given in Table 3. Also given in Table 3 is the time evolution of
. We see from (4.5) and (3.10) that
and
evolve with time in the same way.
Table 3. Time evolution of
for
and of
for
. The finite time
is defined by Equation (4.3). The thin fluid film approximation is satisfied for all time only for
.
4.2. Fluid Velocity on Centre Line
On the centre line,
, we have
by (2.1). From (2.2), since
by (3.10) and (3.20),
(4.6)
But
(4.7)
and therefore
(4.8)
(4.9)
Equation (4.6) becomes
(4.10)
where
(4.11)
For points in the thin film on the centre line,
. Now
at
and
which lies outside the thin film. There is a stagnation point at
,
at which
and
. There is no stagnation point on the centre line for
. For
,
.
At the point of maximum height of the thin film,
(4.12)
Using (3.11) and (3.14) for
we obtain again (4.5). For
and using (3.19) and (3.21) we obtain
(4.13)
The time evolution of
for
is summarised in Table 3.
5. Streamlines
The tangent vector to a streamline is everywhere parallel to the fluid velocity vector instantaneously. The streamlines in the
-plane are obtained by solving numerically the cubic Equation (2.30) for z,
(5.1)
where the coefficients
,
and
are given in (2.21). Now from (3.11), (3.14) and (3.19), (3.21),
(5.2)
and from (3.10), (3.11) and (3.20), (3.21),
(5.3)
Thus from (2.31),
(5.4)
(5.5)
(5.6)
where
and
.
We will mainly consider the streamlines in the
-plane where
,
, and
,
. In this plane the width of the base and the maximum height are constant and the streamline pattern is independent of time and depends only on the parameter
. The evolution of the streamlines with time can be determined using the evolution with time of the half-width of the base,
, and the maximum height
given in Table 2. In the
-plane the cubic equation is given by (2.32) with coefficients
(5.7)
(5.8)
(5.9)
The coefficients
,
and
have the same time dependent factor. The cubic Equation (2.32) for
is
(5.10)
where
(5.11)
Since the left hand side of (5.10) does not depend on time explicitly,
is a constant. The constant
is given a range of values. The streamlines in the
-plane are obtained by solving numerically the cubic Equation (5.10) for
with
. The pattern of streamlines is independent of time because (5.10) does not depend on time explicitly. It depends only on the parameter
.
For
the base is impermeable and the streamlines in the
-plane are plotted in [1]. In Figure 2, leak-off is considered and the streamlines are plotted in the
-plane for strong suction,
,
and for weak suction,
. The streamlines begin on the free surface. They are perpendicular to the base because there is no slip of fluid at the base. The point
where the centre line meets the base is a stagnation point. Since
is proportional to
, near the centre line the suction is weak and spreading due to gravity dominates. As distance from the centre line increases the suction increases and is strongest near the moving contact lines. The streamlines are strongly curved near the centre line because of spreading due to gravity and are weakly curved near the moving contact lines because suction is more important than spreading.
In Figure 3 the streamlines are plotted in the
-plane for fluid injection (blowing) at the base for
, 15, 30 and 50. The magnitude of the fluid injection velocity
increases as
increases and the blowing becomes stronger. Since
is proportional to
,
and the point (0, 0) is again a stagnation point. A dividing streamline extends from this stagnation point and separates the flow into two regions. In the upper region the streamlines begin and end on the free surface. The upper region shows fluid descending due to
Figure 2. Streamlines for leak-off at the base with (a)
and (b)
plotted in the
-plane.
Figure 3. Streamlines for fluid injection at the base with (a)
, (b)
, (c)
and (d)
, plotted in the
-plane.
gravity and spreading. The fluid flow near the centre line is downwards but the flow is forced upwards as it approaches the free surface because of the stronger injection of fluid as the moving contact lines are approached. The lower region describes the injection of fluid into the thin film. The streamlines begin at the base and end on the free surface. The thickness of the fluid layer at the base increases approximately linearly from zero at the stagnation point to the height of the free surface
as
increases because blowing increases as the moving contact lines are approached. The angle which the dividing streamline makes with the base at the stagnation point increases as
increases and therefore as the strength of the blowing increases.
In order to investigate the evolution in time of the streamlines only
and
are required. The evolution of the streamline pattern with time can be determined using
(5.12)
and the results of Table 2 for the evolution of
and
. In Figure 2(a),
which lies in the range of strong suction,
, and
and
as
. The half-width of the base,
, and the maximum height,
, will decrease steadily with time and the thin fluid film will disappear as a line sink as
. In Figure 2(b),
which lies in the range of weak suction,
, and
and
as
. The half-width of the base,
, will increase steadily with time and the maximum height,
, will decrease steadily with time. The thin fluid film will spread over the whole plane
with decreasing thickness and disappear as
. In Figure 3 fluid injection (blowing) is considered and
lies in the range
for which
and
as
. The half-with of the base will therefore increase steadily with time and the maximum height of the thin film will decrease steadily with time. The angle which the dividing streamline makes with the base at the stagnation point will therefore decrease steadily with time.
The evolution of the streamline pattern with time is illustrated in Figure 4 where the streamlines are plotted in the
-plane at time
and
. They are obtained by solving numerically the cubic Equation (5.1) for z as a function of x at
and
where
,
and
are given by (5.4) to (5.6).
In Figures 2-4,
lies in the range
for which the thin fluid film approximation is satisfied for all time. In the next section an approximate analytical solution for the dividing streamline will be derived and its dependence on
and its evolution with time will be determined.
6. Dividing Streamline
We will derive an approximate solution for the dividing streamline in the
-plane. The solution can be expressed in the coordinates
when required
Figure 4. Streamlines for
at (a)
and (b)
, plotted in the
-plane.
later. Since the dividing streamline passes through the stagnation point (0, 0), the constant
in (5.10). The cubic Equation (5.10) becomes
(6.1)
In oder to transform (6.1) to standard form [7],
(6.2)
let
(6.3)
Equation (6.1) becomes
(6.4)
with
(6.5)
(6.6)
On the dividing streamline,
. We perform an expansion in powers of
and derive a solution for blowing,
, to first order in
. No assumption is placed on
or t. We find that this gives a sufficiently accurate result. We see from Figure 3 that as
increases the dividing streamline moves closer to the centre line and the maximum value
on the dividing streamline decreases. The solution should therefore be more accurate as
increases.
Now from (6.5) and (6.6),
(6.7)
(6.8)
The expansion of the cubic equation in powers of
is therefore
(6.9)
and from (6.3),
(6.10)
The discriminant,
, of the cubic Equation (6.9) is [7]
(6.11)
Consider first the solution to zero order in
. Then
,
and
and therefore the cubic equation has three real roots with two roots coincident and one root different [7]. Equation (6.9) reduces to
(6.12)
Thus
,
,
and from (6.10),
,
and
. For zero order in
we have
. The point (0, 3) lies outside the thin film. The solution for the dividing streamline is the stagnation point (0, 0).
For the solution to first order in
,
since
and the cubic equation has three real and distinct roots [7]. We use the trigonometric method of solution described in [1]. The three roots of (6.9) are
(6.13)
where
is the smallest non-negative angle which satisfies
(6.14)
Thus from (6.7) and (6.10),
(6.15)
Consider now
. For
,
and therefore
since
.
Thus for first order in
,
will be small. Using the series expansion for cosine, (6.14) becomes
(6.16)
and therefore since
,
(6.17)
Consider first
given by (6.15) with
. Using again the series expansion for cosine and (6.17),
(6.18)
and therefore
(6.19)
Consider next
given by (6.15) with
. Now
(6.20)
But
is given by (6.18) and from the series expansion for sine and (6.17)
(6.21)
Thus
(6.22)
and therefore
(6.23)
Finally,
is given by (6.15) with
. Now
(6.24)
and using (6.18) and (6.21) we obtain
(6.25)
Thus
(6.26)
The dividing streamline must pass through the stagnation point (0, 0) and satisfy
. The root
is not a dividing streamline because it does not pass through (0, 0). The dividing streamline in the
plane correct to order
is
(6.27)
(6.28)
Let
be the angle that the dividing streamline makes with the positive
-axis at the stagnation point (0, 0). Then
(6.29)
Thus
increases as
increases. In Figure 3 the streamlines are plotted in the
plane for
and 50. From (6.29) the corresponding values of
are
,
,
and
, which agree well with the slope of the dividing streamlines at (0, 0) in Figure 3.
Let the point of intersection, P, of the dividing streamline with the surface of the thin film in the
-plane be
. Then from (3.14),
(6.30)
and therefore from (6.27) for
,
(6.31)
Thus neglecting terms
,
satisfies the quadratic equation
(6.32)
If terms of order
are neglected then the point P is
(6.33)
This approximation for
and
is consistent with (6.29) for
because approximately
(6.34)
As
increases and the strength of the blowing increases,
decreases, verifying that the dividing streamline moves closer to the centre line. The product
occurs in the expansions (6.27) and (6.28) which are valid because
is of order of magnitude unity. The approximation correct to order
is obtained by solving the quadratic Equation (6.32). The positive root of (6.32) can be expressed as
(6.35)
The evolution with time can again be determined by using the transformation (5.12). Let
be the angle which the dividing streamline in the
-plane makes with the positive x-axis at the stagnation point. Then
(6.36)
and using (6.29) with (3.11) and (3.14) at
for
and
we obtain
(6.37)
The asymptotic solution for large values of
is
(6.38)
Thus
decreases as time increases. For large values of
corresponding to strong blowing the gradient of the dividing streamline at the stagnation point decreases with time like
and is independent of
. This makes quantitative the conclusion reached in Section 5 on the evolution of
with time.
Let
be the point of intersection of the dividing streamline with the surface of the thin film in the
-plane. Then
(6.39)
where
and
are given by (3.11) and by (3.14) at
. Neglecting terms of order
and using (6.33),
(6.40)
and therefore for large
,
(6.41)
This approximation for
and
is consistent with (6.38) for
because approximately
(6.42)
Keeping terms of order
we obtain from (6.30),
(6.43)
where
is given by (6.35) and
(6.44)
In the
-plane the dividing streamline is obtained by solving numerically for z the cubic Equation (5.1) with
and
,
and
given by (5.4) to (5.6). From (6.27) and (6.28) the approximate solution is, correct to second order in
,
(6.45)
(6.46)
In Figure 5 the approximate solution, (6.45) and (6.46), for the dividing streamline is plotted in the
plane for
and times ranging from
to
. We see clearly that as t increases the angle
that the dividing streamline makes with the x-axis steadily decreases because spreading is stronger than blowing near the centre line. In Figure 6 the approximate solution, (6.45) and (6.46), for the dividing streamline is compared with the numerical solution of (5.1) with
. The approximate solution is more accurate as
increases. The approximate solution (6.43) for
is also compared with the numerical solution and its accuracy also increases as
increases. Figure 6 also
Figure 5. Approximate solution for the dividing streamline in the
-plane given by (6.45) and (6.46) for
at time
,
,
and
.
Figure 6. Approximate solution for the dividing streamline in the
-plane given by (6.45) and (6.46) (----) compared with the numerical solution (—) at time
and (a)
, (b)
, (c)
and (d)
. The point of intersection, P, of the dividing streamline with the surface is given approximately by (6.43).
illustrates the increase in the base,
, with time and the decrease during blowing of the maximum height although the rate of decrease of
decreases as
increases. From (3.14) for large values of
,
(6.47)
7. Conclusions
Analysis of the streamlines inside the thin fluid film contributed to the understanding of the effect suction and blowing at the base has on the fluid flow inside the thin film. The properties could be analysed both analytically and numerically
because an invariant solution could be derived when
is proportional to
.
The properties of the fluid flow in this study and in [1] where
is proportional to
can be expected to be present in other models for suction and blowing. If more realistic models cannot be solved analytically these properties will be a useful guide in a numerical solution. In both studies, there is a range of suction and blowing for which the thin fluid film approximation is satisfied and outside this range, the solutions are not physically acceptable. In both models a dividing streamline exists, starting at a stagnation point, which separates the flow into two regions, a lower region at the base consisting of rising fluid and an upper region consisting mainly of descending fluid. Both models have fluid variables which tend to zero or infinity in a finite time. The transition solution to an infinite limiting time is an exponential solution.
For the study presented here the thin fluid film approximation was satisfied for
, that is, for moderate to weak suction and the whole range of blowing. The exponential solution, generated by a special case of the Lie point symmetry, and the solutions which tended to zero or infinity in a finite time, lie outside this range. The study complemented the investigation in [1] for which the thin fluid film approximation was satisfied for
, that is, for the whole of suction and for weak to moderate blowing. There was unexpected behaviour in the evolution of
and
caused by the relative importance of spreading due to gravity and suction or blowing. The base half-width
for all blowing and even for weak suction for which
. The maximum height
for all suction and also for all blowing for which
, although
decreases as
increases. Because blowing is weakest at the centre line and strongest at the moving contact lines the dividing streamline generated an approximate V-shaped surface along the length of the two-dimensional thin fluid film, with vertex of each section at a stagnation point. This compares with the dividing streamline in [1] which also passed through a stagnation point and generated an approximately horizontal surface that dipped down at the moving contact lines.
Further research can be undertaken. Surface tension was neglected and there was no slip at the fluid/base interface. The streamlines in a thin fluid film with surface tension or slip could be investigated to give a better understanding of the fluid film flow inside the thin film.
Acknowledgements
DPM thanks the National Research Foundation, Pretoria, South Africa, for financial support under grant number 132189.