Reducing Friction with a Liquid Film on the Body Surface

A flow of a thin layer of liquid is simulated on a flat surface of a body located in a stream of air. Liquid film on the surface of the body reduces frictional resistance and can be used as a boundary layer control element. The paper presents a mathematical model of the film flow on a half-plane, located at an angle to the horizon. The fluid flow is determined by the force of gravity and friction from the external air current. A model of an incompressible viscous fluid is used in the boundary-layer approximation. The terms of the motion equation are averaged over the film thickness according to the Leibniz rule. In the cross section of the film, a quadratic law is adopted for the distribution of the longitudinal velocity, taking into account friction on the film surface. An analytical solution of the problem is obtained in the form of series in powers of the small parameter for determining the film thickness and the average longitudinal velocity along the length of the plate. It is shown that the friction decreases with flow around a half-plane with a film of liquid on the surface.


Introduction
Problems of liquid film flow are quite often encountered in technical applications and are discussed by many authors.Paper [1] presents the results of investigation of films for a gas-liquid two-phase flow in short horizontal pipes.The obtained results are compared with the experimental data.In paper [2], the dynamic characteristics of the film flow with a counter-flow of air are given.The calculated data is experimentally confirmed.
A three-dimensional dynamic model of film flow on a cylindrical surface is considered in paper [3].The film is affected by force of gravity and surface tension, as well as the strength of resistance from the counter flow of gas.In paper [4], the solution of the conjugate problem for the air flow over a flat plate with a thin liquid film on the surface is presented.The film is formed by extruding liquid from a porous wall.In paper [5], the film is formed from the external medium.A mathematical model is obtained for the zero approximation.The results of theoretical and numerical studies of liquid films are presented in papers [6,7], and experimental studies have been performed in [8][9][10].
A model of a thin film flow along an inclined plane is presented in paper [11].The terms of Navier-Stokes equations and the energy equations are averaged over the depth of the liquid.A mathematical model of thin film evaporation was proposed in paper [12].In paper [13], experimental results and numerical modeling of the night flow on a cooled wall are presented.The measurements were carried out on a flat plate.The problem of condensate film stability on a vertical surface is analytically studied in paper [14].An analytical model for the laminar flow of a condensate film on the outer surface of an isothermal vertical pipe is presented in [15].Despite the large number of works on the flow of liquid films, there is still a considerable interest in this phenomenon.Liquid film on the surface of the body reduces friction resistance and can be used as a control element for the boundary layer.In this study, the flow of a thin layer of liquid on a flat surface of the body in the airflow is simulated for its ability to reduce frictional resistance.

Materials and Methods
As an element of the surface, we choose a plane and consider the steady-state flow of a liquid film under the action of gravity and the incoming air flow (Figure 1).The film is formed by squeezing the liquid through the porous surface of the plate, which is located at α angle to the horizon.The liquid flows evenly into the film at V k speed.surface is analytically studied in paper [14].An analytical model for the laminar flow of a condensate film on the outer surface of an isothermal vertical pipe is presented in [15].Despite the large number of works on the flow of liquid films, there is still a considerable interest in this phenomenon.Liquid film on the surface of the body reduces friction resistance and can be used as a control element for the boundary layer.In this study, the flow of a thin layer of liquid on a flat surface of the body in the airflow is simulated for its ability to reduce frictional resistance.

Materials and Methods
As an element of the surface, we choose a plane and consider the steady-state flow of a liquid film under the action of gravity and the incoming air flow (Figure 1).The film is formed by squeezing the liquid through the porous surface of the plate, which is located at α angle to the horizon.The liquid flows evenly into the film at k V speed.In the boundary-layer approximation, the mathematical formulation of the problem has the form The boundary conditions of the problem will have the following form where μ kinematic and dynamic viscosity, g-acceleration of gravity, − τ tangential stress at the boundary between the film-the incoming air flow.

Methods of Solution
Combining the equation of motion (8) with the continuity equation, obtain the conservative form In the boundary-layer approximation, the mathematical formulation of the problem has the form The boundary conditions of the problem will have the following form where p-pressure, ρ-density, ν ,µ-kinematic and dynamic viscosity, g-acceleration of gravity, τ-tangential stress at the boundary between the film-the incoming air flow.

Methods of Solution
Combining the equation of motion (8) with the continuity equation, obtain the conservative form We carry out the averaging of the terms entering into Equation ( 5) over the thickness of the film.
Lubricants 2018, 6, 25 3 of 7 For this purpose, we integrate (10) with respect to y from 0 to δ, using the Leibniz rule (derived from an integral with a variable upper limit) We use the definition of the mean longitudinal flow velocity The integral of the first term of Equation ( 5) is written in the form (the complete derivative can be replaced by a partial one) Since the thickness of the film is sufficiently small, the viscous forces will play a decisive role in the formation of the flow.In this case, we can take for the longitudinal velocity the quadratic law of variation of the longitudinal velocity in the cross section of the film.We use the expression for the velocity of a film of constant thickness that flows along a flat vertical wall, taking into account friction on the surface Average flow velocity Taking into account (10), the longitudinal velocity (9) can be rewritten For the problem under consideration, the average flow velocity < u >, the friction on the surface of the film τ, and the thickness of the film δ are functions of the coordinate x.
Using (11) we write the following expressions Using the obtained expressions ( 11)-( 15), we rewrite the equation of motion in the following form We write down the law of conservation of mass, for which we integrate the continuity Equation ( 2) with y from 0 to δ Taking into account expressions ( 12) and (17), Equation ( 16) can be rewritten If we exclude a neighborhood of a point x = 0 , then v(δ) << v k , in the first approximation, the continuity equation takes the form We introduce dimensionless variables and, taking into account (20), Equations ( 18) and ( 19) take the form initial conditions of the task Calculation of friction is feasible according to the approximate model.The film thickness is sufficiently small and it can be neglected, then the friction value is determined from the Blasius problem For the speed of the oncoming stream U ∞ = 5 m/s, v k = 8 × 10 −4 m/s and the length of the wall l = 0.5 m.Reynolds number Re = 1.68 × 10 5 , which corresponds to laminar flow.If the injection of the mass is carried out by slightly receding from the edge of the plate, then the distribution of friction on the surface of the film can be considered a constant value and take its average value τ = 1 0.5 0.5 0 τ w (x)dx.
Moving to dimensionless quantities, we obtain T = 0.67 = const.
Note that for the initial data (liquid-water at temperature t = 20 Consequently, the solution of Equations ( 21) and ( 22) can be represented in the form of series in powers of the small parameter The coefficients α m and β m are determined by substituting the expansions (24) into ( 21) and ( 22) and equating the coefficients to zero for identical powers of ξ.

Results and Discussions
For α = π/6 we get for α = π/2 the solution takes the form Graphs (Figures 2-4) present comparative flow characteristics, where 1 998 kg m ρ = dimensionless coordinate is a small quantity 0 1.2 10 Consequently, the solution of Equations ( 21) and ( 22) can be represented in the form of series in powers of the small parameter The coefficients m α and m β are determined by substituting the expansions (24) into (21) and ( 22) and equating the coefficients to zero for identical powers of ξ .

kg m ρ =
dimensionless coordinate is a small quantity 0 1.2 10 Consequently, the solution of Equations ( 21) and ( 22) can be represented in the form of series in powers of the small parameter The coefficients m α and m β are determined by substituting the expansions (24) into (21) and ( 22) and equating the coefficients to zero for identical powers of ξ .Then we consider air stream flow by a flat plate (the Blasius problem), but with a movable wall , the results are between 1 and 3.The friction reduction, calculated from the average value, is 14.6%.Then we consider air stream flow by a flat plate (the Blasius problem), but with a movable wall

Results and Discussions
Such a formulation allows a self-similar solution for u 0 = const (for u 0 = f (x) the problem it is unjustifiably more complicated).
The results are shown in the graph (Figure 5) in the form of friction distribution on the plate surface, where curve 1, as before, denotes the slope angle α = π/2, and distribution 3 is the solution of the classical Blasius problem.Then we consider air stream flow by a flat plate (the Blasius problem), but with a movable wall , the results are between 1 and 3.The friction reduction, calculated from the average value, is 14.6%.For the slope of the plate α = π/6, the results are between 1 and 3.The friction reduction, calculated from the average value, is 14.6%.

Figure 1 .
Figure 1.Flow pattern of the film: − v u, components of the velocity vector,

Figure 1 .
Figure 1.Flow pattern of the film: u, v-components of the velocity vector, V k -rate of fluid flow, V ∞ -speed of the incoming flow, δ-film thickness.

Figure 2 .
Figure 2. The average longitudinal velocity.Figure 2. The average longitudinal velocity.

Figure 2 .
Figure 2. The average longitudinal velocity.Figure 2. The average longitudinal velocity.

Figure 4 .
Figure 4.The fluid velocity at film surface.
more complicated).Let us calculate the average velocities on the surface of the film (for ) and perform the numerical solution of problem (27).The results are shown in the graph (Figure5) in the form of friction distribution on the plate surface, where curve 1, as before, denotes the slope angle 2 / π α = , and distribution 3 is the solution of the classical Blasius problem.For the slope of the plate 6 / π α =

Figure 5 .
Figure 5. Tangential stress on the inter-phase surface.

Figure 4 .
Figure 4.The fluid velocity at film surface.

Figure 4 .
Figure 4.The fluid velocity at film surface.
more complicated).Let us calculate the average velocities on the surface of the film (for ) and perform the numerical solution of problem (27).The results are shown in the graph (Figure5) in the form of friction distribution on the plate surface, where curve 1, as before, denotes the slope angle

Figure 5 .
Figure 5. Tangential stress on the inter-phase surface.Figure 5. Tangential stress on the inter-phase surface.

Figure 5 .
Figure 5. Tangential stress on the inter-phase surface.Figure 5. Tangential stress on the inter-phase surface.