Numerical Simulation of the Thin Film Coating Flow in Two-Dimension

Thin film coating is a process of making liquid film cover and deposit base body surface by the way of dipping, spraying, sliding or spin coating, which is a kind of modern surface engineering. It plays an important role in the actual processing, such as improving the surface properties, fine processing, and new surface properties. Analysis of the influence of substrating morphology and fluid flow properties itself on coating fluid motion has an important significance to optimize the thin film coating and improve the quality of the final film. The influence from uneven substrate surface’s geometry configuration on internal motion of the flow field in slip-coating is analyzed by using the FLUENT software as a calculation platform. A two-dimension model of slip coating under isosceles triangle and isosceles trapezoid substrate was established, and thin film coating fluid motions under different configuration parameters were simulated. It is pointed out that the key factor determining the turbulence generation and evolution is the parameter of substrating surface nature. The effects of the change of Reynolds number on turbulent appearance and action area are studied. The velocity contours of fluid field on different substrate surfaces are shown, and the impact of substrate geometry on the backwater region is analyzed.


Introduction
Thin film coating is a kind of modern surface engineering, which generally refers to the process of covering substrate surface with a layer of liquid film.At present, coating technology including electroplating, painting, thermal spraying and va-Z.Y. Liang, H. Zhou DOI: 10.4236/ojfd.2017.73021331 Open Journal of Genetics por deposition, compared with the heat treatment, bead welding and other surface engineering, has several advantages such as less constraint conditions, large space of selection technology and material type and so on, and is used in practical engineering more and more widely.It can play out in three ways such as fine processing, optimization of surface properties, and making new surface properties [1].Specifically, thin film coating can be used to improve the coated substrate surface rust-proof property, anti-corrosion property, abrasions resistance, heat resistance and other characteristics, and has been applied in many fields.In the machinery industry, bridge engineering and electric power industry, thin film coating technique is used to improve micro damage phenomenon of closely matched components [2]; In the field of Aeronautics, film coating technology is often used in aviation engine parts to improve engine efficiency to prolong its service life [3]; In the field of ship protection, application of coating technology can greatly improve the imputrescibility and rust preventing characteristics of shipboard equipment [4].
With the development of material science, liquid thin film begins to take a more and more important role, and consequential quality requirements of the coated thin film are more and more high.In modern industrial applications, thin film coating process needs to meet some special requirements such as: the shape of thin film flow surface is complicated, the fluid coating needs to be carried out in a great disturbance; problem of the actual flow of the film is nonlinear; rheological properties of thin film flow cannot be changed arbitrarily; it should meet the demand of high speed coating industrial production.According to the form of thin film defect, it could be divided into two types, discontinuity and continuity [5].
In the actual production, due to factors of mechanical equipment, operation, production process and fluid properties, coating film maybe have some defects which cannot be completely eliminated such as folds, ripple and bubble [6] [7] [8].Coating defects will certainly affect the coating quality, thereby affecting the properties of the thin film.As the requirements of coating quality and aesthetic appearance of products continuously increase, how to control or reduce the defects of film coating appearance has become an urgent problem to solve, and this needs to do some numerical simulation in the coating fluid mechanics research.
In this paper, research purposes are to analyze the influence of substrating morphology and fluid flow properties itself on coating fluid motion to optimize the thin film coating and improve the quality of the final film.

Coating Fluid Simulation Theory
In fluid dynamics foreign scholars have carried out much research work on the film flow.However, in this field, most previous studies about the flow of the film focused on the flat surface problem.In recent years, research of the liquid film flow characteristics on surface of basal with specific geometry morphology begins.The reason for this study is the surface of the thin film coating substrate is never perfectly flat but quite complex in the actual production.

Couette Flow Model
The film flow on smooth surface can be approximated as a shear flow between two infinite flat plates with rigid surface and specific geometry morphology, which is called the Couette flow.This is a basic research model of the thin film coating.
Two-dimension Couette model is shown in Figure 1.The lower flat plate is fixed, the upper mobile flat plate moves at a constant velocity v along the x direction.The velocity distribution between the two flat plates is as follows: The boundary conditions are as follows: where u is the distribution velocity, μ is fluid viscosity, p is the pressure, h is the distance between the parallel plates, x, y represents the spatial coordinates respectively, and v is a constant velocity along the x direction.
The velocity distribution can be generated by substituting Formula (1-2) into Formula (1-1) as follows: ( ) The velocity distribution of dimensionless form is expressed as In the above formula, D is the dimensionless pressure gradient In the formula, u is the velocity vector, ρ is the density, t is the time, T is the temperature, k is heat conduction coefficient, p c is specific heat capacity, and T S is viscous dissipation.
In Figure 2 Usually, Formula (1-5) is called the momentum conservation equation, Formula (1-6) is called the mass conservation equation, and Formula (1-7) is called the  phenomenon of thin film coating fluid occurs, accompanied by the vortices.Increasing unevenness r will make the eddy current phenomenon more and more obvious, and effect the location of vortex generation.Here defined unflatness of the basal plate with triangular groove is

Influence of Roughness of Trapezoid Substrate on the Flow Field
Thin increasing the roughness of R will expand the scope of the vortex.Here the definition of basal plate roughness with isosceles trapezoid is

Analysis of the Influence of Reynolds Number on the Flow Field of the Thin Film Coating
In theory, film coating fluid should maintain in the laminar or nearly laminar flow state, which means that the flow must be carried out as far as possible at a low Reynolds number range [9].Low Reynolds number means that the viscous force dominates flow process, in other words, that the effect of inertia force can be ignored to some extent [10].As can be seen from the Figure 4 and Figure 5, As can be seen from Figure 5, at a low Reynolds range of 1 -100, when Reynolds number is small, near the basal plate position the volute is small, and flow state is relatively stable.With the increase of Reynolds number, the eddy is more and more obvious to see.But when the Reynolds number increases to a certain extent, volute changes will not be obvious.In other words, compared with the unevenness, the Reynolds number for the eddy current generation and development is not the dominant factor.

Analysis of Thin Film Coating Fluid Velocity
Velocity contours of film coating fluid on the triangle basal plate and trapezoid substrate at the same low Reynolds number are respectively shown as Figure 6 and Figure 7.Here selected to study is the fluid field on basal plate with high roughness that is greater than 0.5, since there is no need for the research on substrate with a low unevenness (such as r < 0.

Conclusion
We proposed a two-dimension model of slip coating under isosceles triangle and isosceles trapezoid substrate, and thin film coating fluid motions under different configuration parameters were simulated.It is pointed out that the key factor determining the turbulence generation and evolution is the parameter of substratum surface nature, with the increase of basal plate roughness, flow separation phenomenon of thin film coating fluid occurs, accompanied by the vortices.
Increasing unevenness r will make the eddy current phenomenon more and more obvious, and effect the location of vortex generation.The effects of the change of Reynolds number on turbulent appearance and action area are studied, compared with the unevenness, the Reynolds number for the eddy current generation and development is not the dominant factor.The velocity contours of fluid field on different substrate surfaces are shown, and the impact of substrate geometry on the backwater region is analyzed.The center dead zones exists can be taken as a reference to judge whether the vortex exists.
: (a) the physical model with isosceles triangle; (b) the computational model with isosceles triangle; (c) the physical model with isosceles trapezoid groove; (d) the computational model with isosceles trapezoid groove.

Figure 2 .
Figure 2. Physical model and computational model of the backplane silhouette is a cyclical fluctuations curve.
3) by the reason that the velocity contours of fluid field on substrate at low trapezoid substrate are similar because its state is laminar at low Reynolds number range, while the flow state on the basal plate with larger roughness is turbulent accompanied by segregation phenomenon appearing.As can be seen from the diagram, the fluid velocity close to the bottom is zero, and the overall velocity distribution is symmetrical (in the range of low Reynolds).Whether triangular or trapezoidal substrate, there is an oval center dead zones (backwater region), where the fluid is in static state, and velocity of the surrounding fluid is not equal to zero.Moreover, due to the influence of unevenness, dead zone center arises in a slightly different position, along with the increase of roughness, location has been moved up slightly.There is no center backwater region in the flow field of laminar state, so the center dead zones can be taken as a reference to judge whether the vortex exists or not.