# Mathematical Modeling of Flattening Process on Rough Surfaces in Thermal Spray

H.Fukanuma

Plasma Giken Co.,Ltd.,Saitama,Japan

## Abstract

Thermal spray layers are formed on rough surfaces;however,the flattening process on rough surfaces has not yet been clarified.A mathematical flattening model which takes into account the roughness of the substrate or previously coated layers is proposed in this paper.As a result of surface roughness,the flattening degree and the flattening time decrease with increasing surface roughness in this model.In addition,the characterization of surface roughness is introduced for the flattening model.Several calculated cases of the flattening model are shown.

THE FLATTENING PROCESS of impinged molten particles onto substrates affects the mechanical and physical properties of thermal sprayed coatings.The flattening process is one of the most important processes in thermal spraying,as it helps one to understand the physics of thermal spray and the characteristics of the sprayed layers.Many flattening models of molten or liquid particle impingement on flat surfaces have been proposed [1-6].These mathematical models were concerned with the flattening processes on smooth surfaces.In addition to these models,many experimental investigations of splat formation on both smooth and rough surfaces have been reported[7-16].

Flattening processes on rough surfaces have been investigated recently.Moreau et al [14] have shown that the flattening degree and the flattening time of molybdenum splats on both roughened glass and molybdenum surfaces were less than on smooth surfaces.Leger et al [15] have shown that the flattening degree of zirconia splats on roughened steel surfaces decreases with increasing roughness in a relationship between the flattening degree and the Reynolds number.Bianchi et al [16] have shown that splat morphology on rough surfaces is different from that on smooth surfaces.

How does the surface roughness influence the flattening process in thermal spray coatings? First of all,we must consider the geometry of the surface roughness.Let’s consider an idealized roughened surface as shown in Fig.1-a.The cones which have height a and bottom radius b are aligned with distance l between the bottom centers of adjacent cones on the smooth,flat surface.

We can assume that an impinging molten particle begins to spread in the radial direction as a thin disk of thickness h at the impinging point of the particle onto the substrate,and that the thickness h is constant during flattening as shown in Fig.2.

We assume that a molten particle impacts onto the idealized rough surface,and we classify the roughness into three cases:(1)the cone’s height a is much larger than the disk thicknessh,(2)the cone’s height is much smaller than the disk thickness,and(3)the cone’s height is approximately equal to the disk thickness.

In the spreading process in case(1),it can be considered that the fluid in the impinging particle,mostly flows along the valleys because the fluid has a velocity component perpendicular to the surface.A part of the fluid encountering the mountains climbs the slopes,then the liquid probably flies off the mountain,and a splash occurs if the radial velocity is large enough compared to the fluid surface tension.On the other hand,the fluid probably travels around the mountains and goes down the valleys if the radial velocity is not large enough for the fluid to climb over the mountains.The perpendicular velocity is largest at the disk center and is smaller farther from the center.If the impinging velocity is large and roughness is much large compared to the splat thickness,splashing possibly begins at impact and some of the fluid splashes away.This is because the radial velocity at impact is 2-5 times the impinging velocity,as shown by Heymann [2] and will be shown in this article,so the radial flow can climb and fly off the mountains in the short time after impact.Thereafter,because the following radial velocity rapidly decreases,and the velocity near the top of the particle becomes quite small when it comes to the surface,the small splat remains.

The mountains’ geometric shapes influence the splash

phenomena.For example,the pyramid mountain as shown in Fig.1-b causes much more splash than the cone mountain,because the flat slope area can carry more fluid than the cone.If the mountains’ arrangement is symmetrical to the axis perpendicular to the surface through the particle center,the splat shape is also symmetrical,but the rim line is both convex and concave. If the particle impinges directly onto the top of the high large mountain,splashing occurs easily,because radial flow is caused at the upper part of the particle and the radial flow direction becomes upward as shown in Fig.3.The actual surface roughness of a blasted substrate or previously sprayed layer is much more complicated than these ideal models,and thus the actual splats will be irregular.The splashing and the splat morphology strongly depend on the topography of the rough surface.When we solve the flattening process in this case,we must consider the interaction of the flowing fluid and every mountain.It is almost impossible to construct the equations to solve the problem,and to find the conditions which cause splashing.

In case(2),flattening processes do not cause splashing because the small mountains cannot produce the velocity or

the force in the vertical directions that is necessary to cause the fluid to take off from the surface.This is similar to a deep river flowing on a slightly rough river bed.In this case,it is satisfactory to consider the flattening process to be like the one on flat surfaces.

In case(3),the spreading phenomena are explained by case(1)or(2)depending on the impinging and surface conditions.If the radial velocity is not large,splashing probably does not take place.If splashing should occur,it is weak.If radial velocity is large,or the mountain slopes’angle makes the flow velocity in parallel to the slopes large enough for the flow to fly off the surface,then splashing occurs.

In addition to surface roughness,the wetability between a splat and the substrate or a previously sprayed layer greatly affects the flattening process,in particular the splashing process.Even if the surface is smooth,splashing can occur as a splat does not wet the substrate at all,or very slightly.The influence of wetabilty is not considered in this paper.

## Modeling

This study refers to the modeling on a rough surface in case(2).It is difficult to take the roughness into considerations in case(1),because of the lack of understanding of the characterization and the topography of such rough surfaces.Also,we need to clarify the particle impact process onto rough surfaces at the moment of collision.It is also difficult to treat case(3);however,in a case without splashing,we might treat that as an extension of case(2).

It is nearly impossible to consider the interaction between the expanding splat and every mountain in the flattening process on a rough surfaces.A better way is to take the average of the fluid flow field on a small area of rough surface in order to characterize the flow field in the valleys.The flow field of the lower fluid as shown in Fig.4,which flows in the valleys is very different from the flow field of the upper fluid,which flows over the tops of the mountains.The lower fluid flows much less than the upper fluid because the mountains impede the lower part of the fluid flow.If the lower flow has negligible affect on the upper flow,we can classify the fluid of

the disk into two parts:the fluid that flows in the valleys and the fluid that flows over the mountains.We assume that the flow field in the valleys between γ and γ+Δγ,although complicated,can be represented by the flow field in the fluid area Δγ×σ,where σ is obtained by taking the average of the fluid depth in the valleys.The entire fluid in the valleys can be expressed as a thin disk by adding all elements in Δγ from the center to outer rim of the splat.Thus,the flattening model can be represented by two disk layers on a smooth surface instead of the complicated splat on the roughened surface.

How do we characterize the surface roughness which affects the flattening process? We can introduce two parameters which influence the splat spreading process on the rough surface:one parameter is the volume of the valley space contained in a unit area on the surface,the other is a parameter which gives an expression of the fluid’s”flowability”in the valley.

The volume of the valleys is a function of the mountains’ height,shape,and arrangement.In the case of the cones surface in Fig.1-a,because the valley volume ψ is a function of a cone height a,radius of cone bottom circle b and distance l between the bottom centers of adjacent cones,the volume can be obtained on unit surface area as

When b = l,for example,when the bottom circle rims contact each other,we obtain ψ = (1-0.302)a.We can similarly obtain the volume in an idealized case in which i depends on surface topography.Then,Eq.1 is expressed as

where i is a parameter which is smaller than 1,since b≤1.The volume ψ is expressed as the product of mountain height and a coefficient as shown Eq.2.Substituting (1-i)=j into Eq.2,the expression becomes

Valley volume ψ can be expressed by two parameters j and a.In general,j is smaller than 1,since 0 ≤ i ≤ 1.

On the other hand,determining the parameter which expresses fluid mobility depends on the shape of the mountains and their arrangement.We consider that the flow field in the splat is represented in cylindrical coordinates (0,γ,θ,z),and if the flow field in a splat on a smooth,flat surface is U ( Uγ, Uθ Uz ),we assume that the flow field in the valleys on a rough surface is k × Uγ,where k is a coefficient smaller than 1.The flow on a smooth surface is expressed whenk = 1.We have now converted the problem of the flattening process on rough surfaces into the flattening process on smooth,flat surfaces by introducing parameters which characterize the rough surface.The parameters are a,j,and ψ instead of a and ψ.

In Fig.5,we assume the energy dissipation in region I is negligible,because the velocity gradient is most likely small.If the velocity gradient is large in that part,the spherical part on the disk will change shape.The splat formation pictures in the work by Engel [1] and in the previous work by Ohmori and myself [6],however,show that the shape of the spherical part keeps its sphericity quite well.Thus,we consider the flow field in just the disk(region II and III).We can assume the radial flow field Uγ in the disk is:

where Uγ is the radial flow component,C1,C2 and α are constants,t is the time after the impact,h is the disk thickness,σ is the average thickness of the fluid in the valleys,and k is the parameter which indicates fluid mobility in the valleys.Equation 4 and 5 are γ-component of the fluid field in lower and upper part of region II,respectively,and Eq.6 and 7 are γ-component of the fluid field in lower and upper part of region III,respectively.We assume R0 = d0 / 2,where d0 is

the particle diameter at impact.

The θ-component flow Uθ equals 0,because the flow field should be symmetrical around theZ-axis.The internal friction in the perpendicular direction is negligible,because the disk is so thin that the energy dissipated by the velocity gradient in the perpendicular direction is negligible compared to that in the radial direction.

The assumed flow field show that the radial flow velocity will increase proportionally to the square of the height from the fluid bottom,and will increase proportionally to the radius from the center to a certain radius R0,and will decrease inversely to the radius from R0 to the outer rim.We call these region II and III,respectively,as shown in Fig.5.The formulas are chosen so that the velocity gradient is equal to zero on the top surfaces of region II and III.The flow field will depend on time with an exponential exp(-α1),after taking into consideration the fact that the spreading rate rapidly decreases with the passage of time [6].

The disk volume of radius R(t) and thickness h at time t after impact is expressed as

where R0 is the disk radius at time t0.The integration on the right hand side represents the fluid volume which flows out from region II into region III through the boundary at R0 during the time between t0 and t.Substituting the sum of Eq.6 and 7 into Eq.8,the expression becomes

We obtain the disk radius as a function of time by solving Eq.9 for R ,and finding coefficientsC2,α and the disk thickness h.Coefficients C1,C2 and &aplha; are found by solving the equations of the boundary conditions.The boundary conditions are:

(i)when a thin disk of thickness h begins to spread at the impact point in the radial direction parallel to the surface, and the particle keeps moving perpendicular to the surface at the impact velocity ν0 for a time after impingement,the fluid volume ω in Fig.6 is expressed as

where γ is the circle radius given by cutting the particle by planes parallel to the substrate surface,and H is the particle’s moving distance during time t(i.e. H = ν0t in the short time after impact).Since the fluid volume ω converts to the thin disk,namely ω = πhR2,the next equation holds.

Differentiating Eq.11 by t,and since ( H + h ) /d0 << 1,the radial expanding rate of the disk can be expressed as

where d0 is the particle diameter at the time of collision.

Assuming t0 ≅ 0,we obtain

Differentiating Eq.9 by t,the expression becomes

When t = t0,Eq.14 becomes

We obtain the next equation from Eq.13 and 15.

(ii)Because the fluid volume which flows out through the boundary between region II and III is the same as the fluid volume which flows into region III,the next equation holds.

Solving Eq.17,We obtain

(iii)At t = ∞,the radius becomes R = Rm.Substituting t = ∞,t0 = 0 and R = Rm into Eq.9,we obtain

From Eq.16,18 and 19,we obtain

Substituting Eq.20,21 and 22 into Eq.9,and then substituting t0 = 0 into the equation and directing our attention to the relationships R = D / 2, R0 = d0 / 2 , Rm = Dm / 2 and h = d03 / 6Rm2,we obtain an equation which expresses the splat diameter expanding as a function of time as follows,

If d02 / Dm2 << 1,Eq.23 is simplified as

The unknown is still Dm in Eq.23 and 24.If we find Dm the problem has been solved.The parameter k,j and α do not appear in Eq.23 and 24.They must be hidden in Dm.

Let’s find Dm.The next equation holds because of the conservation of energy.

where μ and γ are the viscosity coefficient and the surface tension of the fluid,respectively.In Eq.25,the left hand side is the particle kinetic energy at impact,the integration is the dissipate energy of friction by shear stress in the fluid,and the second term on the right hand side is the increase in surface energy due to the increase in splat surface area during the flattening process.Substituting Eq.4,5,6 and 7 into the integration term in Eq.25,and since the volume element dV is equal to γdrdθdz in cylindrical coordinates,the expression becomes

Solving Eq.26 after substituting Eq.20,21 and 22 into Eq.26,the expression becomes

Substituting Eq.27 into Eq.25,the expression becomes

When surface energy is negligible,Eq.28 is simplified as

Approximating the term in brackets on the right hand side in Eq.29,and substituting ρν0d0 / μinto Eq.29,the expression becomes

Substituting h = 2d03 / 3Dm2 into Eq.30,the expression becomes

Since the fluid volume ψ in valleys on unit area equals the average depth δ of the fluid in valleys on unit surface area, namely δ ＝ ψ ＝ j・a ,substituting δ ＝ ψ ＝ j・a into Eq.31,the expression becomes

We have obtained the equation expressing the relationship between Reynolds number and flattening ratio D (t = ∞) = Dm.the equation includes the parameters characterizing surface roughness.If k = 0,that is,the fluid cannot flow in the valleys,Eq.32 is simplified as

Solving Eq.32 or 33 for Dm / d0 and substituting the solution into Eq.23 or 24,we obtain the deformation ratio of the expanding splat diameter to the particle diameter at impact as a function of time.In particular,since Eq.33 is a cubic equation of (Dm / d0)2,the equation can be solved.Furthermore,if Dm2 / d02,the simplified solution is found as follows,

Equation 34 is the simplified formula of deformation ratios on rough surfaces at k = 0.

If (3ja/2d0)3 is negligibly small,Eq.34 becomes

We obtain more simplified solution of the flattening ratio.This equation holds when surface roughness is much small compared to the particle diameter at impact.

The equation by substituting the solution of Eq.32,33,34 or 35 into Eq.23 or 24 expresses the flattening ratios of the expanding splat with passage of time on rough surfaces.

## Numerical results and discussion

Some numerical results of this model show how roughness affects flattening processes,such as,the influence of surface roughness to flattening time and flattening ratios.The case of a molten nickel particle will be shown and the conditions at impact are shown in Table 1.

The flattening ratios of a spreading nickel splat with different conditions of the rough surface with time are shown in Fig.7.It shows that flattening ratios and spreading rates become smaller with increasing surface roughness,which the parameter a represents.The parametera greatly affects the flattening ratio and flattening time compared to parameter k,which represents fluid mobility in valleys on rough surfaces as shown Fig.8 and 9.It is noted that the product of a and j influences the flattening process,not a or j alone,as shown in Eq.32,33,34 or 35.Figure 8 shows the relationship between the mountain height of surface roughness a and the flattening time t0.9 when the splat diameter D(t) reaches 0.9Dm.The time t0.9 is obtained by substituting D = 0.9Dm into Eq.24 followed by solving for t as follows,

The graph shows that the flattening time decreases with increasing roughness.The influence of the parameter k on the flattening time is less than that of the parameter k.The flattening time is 1.14 μsec on a smooth,flat surface when a 50 μm nickel particle impinges at velocity of 100 m/sec;however,it is 0.86 μsec on a rough surface of 2.2 μm and fluid mobility of k = 0as shown in Fig.8.The deviation of flattening time on a rough surface from that on the smooth surface is approximately 25% less and indicate that a greatly affects the flattening time.

Figure 9 shows that the flattening ratio is affected strongly by the parameter of mountain height a on a rough surface,but affected weakly by the parameter of flowability k.Figure 10 shows that the increasing rates of flattening ratios with Reynolds number become smaller with increasing roughness.

The flattening ratio is strictly a function of a / d0 and Re,not directly of a as shown in Eq.32,33,34 and 35.surface roughness is related to the impinging particle diameter or strictly to the splat thickness,because a larger particle expands with a thicker splat compared to a smaller one with other conditions the same.

When we consider applying this model to very rough surfaces,although actual rough surfaces are much more complicated than the idealized ones,it is difficult to determine the parameters j and a.We can probably determine j and a in the following ways:for example,taking the average a of the mountains’height in a certain area on the actual surface by using an instrument which can measure the topography of rough surfaces,such as a laser microscope;the average a can be applied to the model instead a.Also,using the microscope,the valleys’volume can be determined by putting a hypothetical plane at proper height and angle over an actual surface,the valleys’volume ψ can be calculated under the plane.Since a and the valleys’volume ψ are determined,and from j = a × ψ ,j is determined.

On the other hand,it is very difficult to determine the fluid mobility k in the valleys on both real and idealized rough surfaces.Actually,the model is useful even if the mobility k is treated as k = 0,when the mountain height is small.We need further studies to find practical parameters to characterize real rough surfaces.

The mobility k in valleys seems actually to depend on the flow velocity of fluid,although k is considered as constant in the model.If the velocity is larger,k is expected to become smaller,namely,the mobility become small.

The relationship between surface roughness and splashing has not been understood yet.The mechanism how roughness causes splashing is also not clarified.Impact velocity might affect splashing phenomena or,more strictly,the radial expansion rate at impingement.The radial expansion rate at impact is described in Eq.13.Substituting h = 2d03/ 3Dm2 into Eq.13,then substituting Dm / d0 = 1.06Re1/6 into the obtained equation,the expression becomes

Equation 37 shows the radial expanding rate at impact is 441 m/sec in the case of a molten nickel particle as shown Table 1. In this case,the radial expansion rate is about 4.4 times the impact velocity.Although the splashing conditions are not clarified,the impact velocity might be one of the factors which cause the splashing.A lager impinging velocity has more chance of causing splash if other conditions are the same.We must note the radial flow velocity is much larger than the impact velocity.Larger impact velocity makes the expanding splat thinner,If the splat is thinner,the flattening behavior is more influenced by roughness.Therefore,the parameter jaν0/d0 probably expresses more strictly the tendency to

splashing than impact velocity ν,however,the critical value is not clear.

In conclusion,a formula which describes the flattening process on rough surfaces is obtained in this study by converting the process on a rough surface to that on a smooth,flat surface as a splat consisting of two layers with different flow fields,and introducing parameters characterizing rough surfaces.This model shows the flattening ratio and the flattening time decrease with increasing roughness which have been experimentally confirmed by some studies.

## References

1. O.G.Engel,Journal of Research of the National Bureau of Standards,Vol.54(5),281-298,May(1955).

2. F.J.Heymann,Journal of Applied Physics,Vol.40(13),5113-5122,December(1969).

3. J.Madejski,Int.J.Heat Mass Transfer,Vol.19,1009-1013,(1976).

4. G.Trapaga,J.Szekely,Metallurgical Transactions B,Volume 22B,901-914(1991).

5. H.Fukanuma,Journal of Thermal Spray Technology,Vol.3(1),33-44,March(1994).

6. H.Fukanuma and A.Ohmori,Proceedings of the 7th NTSC,Boston,Massachusetts,563-568,June(1994).

7. J.Yankee and B.J.Plekta,Journal of Thermal Spray Technology,Vol.2(3),271-281,September(1993).

8. Fantassi,M.Vardelle,A.Vardelle,and P.Fauchais,Journal of Thermal Spray Technology,Vol.2(4),379-384,December(1993).

9. Vardelle,A.Vardelle,A.C.Leger,P.Fauchais,and D.Gobin,Journal of Thermal Spray Technology,Vol.4(1),50-57,March(1995).

10. ChangJiu Li,A.Ohmori and Y.Harada,Proceeding of ITSC’95,Kobe,Japan,333-339,May(1995).

11. Fukumoto,S.Katoh and I.Okane,proceeding of ITSC’95,Kobe,Japan,353-358,May(1995).

12. Montavan,S.Sampath,C.C.Berndt,H.Herman and C.Coddet,Proceeding of ITSC’95,Kobe,Japan,365-370,May(1995).

13. Montavan,C.Coddet,Proceeding of 8th NTSC,Houston,Texas,285-289,September(1995)

14. Moreau,P.Gougeon,M.Lamontagne,Journal of Thermal Spray Technology,Vol.4(1),25-33(1995)

15. Leger,M.Vardelle,A.Vardelle,B.Dussoubs,P.Fauchais,Proceedings of the 8th NTSC,Houston,Texas,169-174,September 1995.

16. Bianchi,A.Grimaud,F.Blein,P.Lucchee,and P.Fauchais,Journal of Thermal Spray Technology,Vol.4(1),59-65,March(1995)