# Mathematical modeling of splat formation at off-normal angles in thermal spray

H.Fukanuma,Toda City/J,and C.-J. Li,xian/TJ

## Mathematische Modellierung der Spritzerausbildung beim thermischen Spritzen abweichend von der Bauteilnormalen

In dieser Studie wird ein Modell für die Spritzerbildung bei anormalen Spritzwinkeln beim thermischen spritzen vorgeschlagen.Es zeigt sich,daß die Spritzer elliptisch werden,wenn die Spritzwinkel anormal sind,und daß sich die Ellipse dehnt,je mehr der Winkel von der Lotrechten adweicht.Mit dem Modell wurden Versuchsergebnisse der Spritzerdehnung bei anormalen Aufprallwinkeln verglichen.

## 1. Introduction

This paper aims to clarify the molten particle deformation in off-normal angle impingement under certain physical and mathematical assumptions.The splat formation at off-normal angles in thermal spray is not well understood,although there have been some qualitative investigations of splat morphology and the characteristics of sprayed layers in off-normal angle spray.Montavan et al reported that the elongation ratio increased from approximately 1 to 2 when the spray angle changed from 90° to 30°.However,the average equivalent splat diameter was not significantly influenced by the changes of spray angle[1].In another study,Montavan and Coddet showed that a lower spray angle resulted in a greater slope of the splat thickness[2].Smith investigated the spray angle dependence of deposition efficiency,surface roughness and porosity for several spray materials[3].That particular study is practical and useful,however it did not include information about how the splat morphology depended on spray angle.On the other hand,Sobolev reported a mathematical flattening model taking into account impact angles[4].However,the model is treated under the assumption that splat circularity is held between 45° and 90°.The theory does not include an analysis of splat shapes whose circularity decreases when the spray angle changes from the perpendicular.

The relationship between the splat elongation ratio and spray angles is found in this paper under the assumptions that(i) the splats are elliptic at off-normal spray angles,(ii)a molten particle which impinges onto a substrate at declined angle moves away from the impact point along the substrate and(iii)the intersection of the substrate surface and a vertical line through the particle center is the second focal point of the ellipse,if the impact point is the first focal point.

## 2. Modeling

A molten particle whose diameter is do impinges onto a flat smooth substrate at velocity voand impact angle φ.Then,a radial flow begins to spread from the particle impingement point.The splat thickness h is assumed a constant during the flattening.When the impinging angle is φ,the particle velocity components perpendicular and parallel to the substrate at impact are vo sinφ and vo cosφ,respectively.The particle moves towards the substrate and concurrently moves along the substrate.The splat is a circle at the very moment of impact,then the splat becomes elliptic.The oval shape is growing elongated and larger with elapsed time t as shown in Fig.1.The impact point of the particle is one of the two focal points in the ellipse.The other focal point moves away from this point of initial impact and causes the formation of an oval disk,while,simultaneously,supplying molten liquid into the disk.

We consider two sets of cylindrical coordinates O(r,θ,z) and O'(r’,θ’,z).The origin O of the first set is placed at the impact point,with r and θ on the substrate surface and z extending perpendicularly upward.The origin O’ of the other set is the second focal point that is moving with the particle as shown in Fig.2.If we are on the moving coordinates O'(r’,θ’,z) and observe the particle flattening,the particle does not apparently move laterally and the liquid from the particle flows out radially into the spreading disk.Thus the r’ component ur of the flow fieldU'(ur,uθ,uz) in the splat is assumed as follows;

where C1,C2 and α are constants that are determined by initial and boundary conditions,andR0=d0/2 [5,6].The θ’ and z components in r’ can be negligible because their velocity gradient is negligibly small since the flattening disk is too thin in comparison with its radius.When r’ in the moving coordinates is expressed in the original coordinates r and θ,the next equation holds;

The term rp is the distance between the two origins of the two sets of coordinates and also between the first and the second focal points in the ellipse.We assume that rp is

where “n” is a constant and 0 < n ≤ 1,and "k" is the focal length of the ellipse.The second focal point moves at velocity vo cosφ at impact and the velocity decreases with time because the term exp(-nαt) decreases rapidly with increasing time. If the droplet supplies all its liquid into the spreading oval disk and the liquid flows out through the area around the moving focal point whose radius and thickness are R0 and hrespectively,then,the liquid volume of the oval splat from time t = 0 to t equals the volume that flows out through the boundary r'=R0 for the same time range.Then;

where “R” is the elliptic radius from the initial focal point that is the impact point.The elliptic curve is expressed as,

where g is a constant and is expressed as,

where am is the long radius of the splat that stops its spreading just when it loses its kinetic energy and km is the focal length in the final oval splat,and “e” is eccentricity.e is defined as,

where “a” and “k” are the elliptic long radius and focal length of the splat at timet,respectively.

The term “α” is determined by solving equation(5),and C1 and C2 are determined by the procedure shown in a previous study[5,6].α is a factor that shows how the particle is decelerated when the droplet moves down onto the substrate.On the other hand,exp(-nαt) in equation(4)is also the deceleration factor of lateral motion.The term nα is used instead of justα since the lateral speed of the droplet on the flattening disk is decelerated less than that in the perpendicular direction,since sliding motion has a smaller resistance than that of vertical motion.A lower value of “n” indicates lower resistance and indicates a more elongated ellipse.

The area of the elliptic shaped splat is determined by solving the equation of energy dissipation by friction in the splat that is equal to the kinetic energy that the droplet holds at impingement.When the kinetic energy of the impact particle is Eλ then,

where “μ” is the viscous coefficient of the droplet liquid and “V” is the volume of the flattening disk at time t.

## 3. Experimental Procedure and Results

### 3.1 Experimental Procedure

Copper and nickel powders were sprayed in Xian Jiaotong University and Plasma Giken,respectively.

Copper powder whose particle sizes ranged from 74 to 88μm was sprayed onto a stainless steel substrate with a plasma torch.The substrates were preheated to 250℃ on the back side.Shielding plates were fixed in front of the substrates to obtain splats with minimal errors of angle.Spray angle φ was changed from 60°to 90°in 5°steps.Argon and hydrogen gas were used as the arc gas and input current and voltage were 400A and 60V,respectively.The powder was fed into the plasma jet vertically upward.

Nickel powder ranging from -100 to +325 mesh was sprayed onto a stainless steel substrate by a plasma torch using argon and helium as the arc gas.Arc current and voltage were 800A and 35V,respectively.The substrates were preheated to 350℃ on the back side.The spray angle was changed from 15°to 90°in 15°steps.

### 3.2 Experimental Results

Sprayed copper and nickel splats were collected and photographed.The typical morphology of copper and nickel splats are shown in Fig.3 and Fig.4,respectively.Splat diameters were measured choosing ten samples of copper and one splat of nickel at each angle.The ratio of the long diameter to its orthogonal short diameter was calculated for each.

The Cu and Ni splats at spray angles more than 60°are almost circles.The Cu and Ni splats at spray angles greater than 45°exhibit two different parts that are the central and outer areas.The elongation ratio of the central area is smaller than that of the entire splat.In the photographs,the splat surfaces of the central sections are smooth and the peripheral areas are roughened.

The Ni splats at 30°and 15°look splashed because the peripheral sections were not observed.Thus,the real elongation ratios could be larger than those of the photographs.

## 4. Comparison of The Model to The Experimental Results

From the model,the elongation ratio ER of the splats is defined by the following equation.

ER is determined by substituting km into equation(10)which is solved from equation(4)and(5)as a function of radius am of the splat and impinging angle φ.Then;

The results of the experiment and the model are shown in Fig.5.

The graph shows that the experimental data agree quite well with the theoretical curve.Although the data near the normal angle are larger than the predicted values,this arrives from measurement error because the splats are not ideal ellipses and the peripheral curves are not smooth.As well,the elongation ratio is close to 1 due to the spray angle being nearly perpendicular.

## 5. Conclusions

The results of the experiment show that the elongation ratio increases with an increasing inclined spray angle,and that the model can predict the relation between the elongation ratio and the spray angle.

We need experiments on other spray materials and under more varied condition to clarify the processes of the off-normal spray angles and to confirm the model adaptability.

## 6. Literature

1. G.Montavan,C.Coddet S.Sampath H.Herman,C.C.Berndt:Vacuum Plasma Spray Forming of Astroloy:An Investigation of Processing Parameters,Proceedings of the 7th Thermal Spray Conference,Boston,20-24 June,1994,pp.469-475.

2. G.Montavan and C.Coddet:3-D Profilometries of Vacuum Plasma Sprayed Nickel-Based Alloy Splats Using Scanning Mechanical Microscopy,Proceedings of the 8th National Thermal Spray Conference,Houston,11-15 September 1995,pp.285-289.

3. M.F.Smith,R.A.Neiser and R.C.Dykhuizen:An Investigation of the Effects of Droplet Impact Angle In Thermal Spray Deposition,Proceedings of the 7th National Thermal Spray Conference,Boston,20-24 June 1994,pp.603-608.

4. V.V.Sobolev and J.M.Guilemany:Droplet Flattening Thermal Spraying at Off-Normal Angles,Proceedings of the 15th International Thermal conference,Nice,France,25-29 May 1998,pp.497-502.

5. H.Fukanuma:A Porosity Formation and Flattening Model of an Impinging Molten Particle in Thermal Spray Coatings,Journal of Thermal Spray Technology,Vol.3(1)March(1994),pp.33-44.

6. H.Fukanuma and A.Ohmori:Behavior of Molten Droplets Impinging on Flat Surfaces,Proceedings of the 7th National Thermal Spray Conference,Boston,20-24 June 1994,pp.563-568.