1. Introduction
In recent years, additive manufacturing (AM) technologies for the production of metal parts have drawn an enormous surge of industrial interest; however, process reliability and component quality are not enough for mass production [
1]. In opposition to systems able to selectively melt a powder bed, high build rates and larger part volumes can be obtained by exploiting direct energy deposition (DED) technologies [
2,
3]. In this technology, previously trademarked as laser engineered net shaping (LENS), laser metal deposition (LMD), or direct metal deposition (DMD), a heat source generates a melt pool on a metallic surface into which feedstock is deposited with the consequent building of the parts layer by layer. At the moment, several solutions are available, using wire or powder as feedstock; in former versions, the processes could be considered a type of welding technology, whereas the latter ones are very similar to laser cladding technologies.
In view of the introduction of AM machines to the industrial world, it is necessary to study and develop monitoring and control systems capable of guaranteeing series production capability and reproducibility of these innovative processes. For this reason, recent AM technology reviews have repeatedly called for real-time, closed loop process controls and sensors to ensure quality, consistency, and reproducibility across AM machines [
4].
DED systems employing powder feeding are highly sensitive to working conditions; any variation of process parameters during the deposition and melting of the metal powders can influence the quality of the component, compromising the process efficiency [
5]. In particular, the deposition quality can be affected by intrinsic parameters (related to the properties of the substrate and the metallic powder such as geometry, thermal diffusivity, absorptivity, thermal conductivity, and heat capacity) and extrinsic parameters (related to the laser, the powder feeder, and the positioning system). In particular, extrinsic parameters are strictly related to the specific print head of the employed AM machine, in which powder injection nozzles and the laser beam are assembled according different schemes [
6]. Indeed, common DED or cladding powder injection systems can present:
an off-axis configuration, in which a single powder flow laterally passes through the laser beam;
a continuous coaxial configuration, in which a conical powder flow surrounds and interacts with the laser beam;
a discontinuous coaxial configuration, with several powder flows from different injection nozzles (i.e., multiple nozzle deposition) distributed around the laser beam (up to eight in a robotized laser-based direct metal deposition system recently developed [
7]).
The off-axis powder injection implies a strong relation between the deposition rate and the scan direction [
8]; in the continuous coaxial version, as opposed to the discontinuous axial one, the tilting of the powder injection nozzle is restricted, limiting the potentialities of DED [
6]. Multiple nozzle deposition is suitable for parts with high geometrical complexity, even if it is more difficult to ensure flow uniformity and direct the powder to a specific region of interest [
9]. In the case of both continuous coaxial configuration and discontinuous coaxial configuration, the spot of the powder flow in correspondence with the deposition plane usually ranges between 5 mm and 8 mm [
6,
10,
11,
12,
13].
As already stated, the process efficiency can be influenced by the design of the powder injection nozzle, since it is capable of modifying interactions among the metallic particles, laser beam, and melt pool [
14]. Further, the control of extrinsic parameters can be useful for adjusting deposition errors during the process. For this reason, in order to improve the process reliability and quality, closed loop control systems should be introduced in DED machines [
15]. Indeed, in closed loop control, it is possible to exploit different devices to monitor the deposition and melting processes, by using the recorded data as feedback to setup the optimal process parameters.
In this paper, the first part of a larger research project is presented, in view to design a tailored deposition closed loop control to improve the performances of DED systems, guaranteeing their maximum efficiency during the process. In particular, this preliminary study is focused on the performances of a double chamber nozzle solution, controlled through the evaluation of the shape of the powder flow at the exit of the nozzle outlet and the spread of the powder particles on the deposition plane. During experimental tests, the powder flow is monitored with a high-speed camera, in order to confirm its simulating behavior obtained through a Computational Fluid Dynamics CFD model [
16,
17]. The experimental tests demonstrate that, in the analyzed configurations, the shielding gas seems to have a significant influence on the powder distribution and powder flow geometry.
2. Materials and Methods
A new deposition nozzle solution for DED is designed at the Department of Innovative Technologies (DTI) of University of Applied Sciences and Arts of Southern Switzerland (SUPSI), Manno-CH (see
Figure 1a). It consists of a hybrid solution of continuous and discontinuous coaxial systems, with a double chamber coaxial nozzle (see
Figure 1b), where the central powder-gas flow is shielded by an external annular flux of inert gas in order to prevent the oxidation of the melt pool, even in the case of no-hermetic deposition chamber, and to improve the catchment efficiency during the deposition, limiting the spread of the powder flow coming out from the nozzle outlet. Two deposition nozzles with an inclination of 30° are placed in a testing chamber employed to characterize the powder flow coming out from the nozzle outlets (see
Figure 2a). The deposition nozzles are connected with a flexible feeding system so that the control of shielding, carrier, and powder mass flow can be independent and precise (see
Figure 2b).
In order to investigate the performances of the new nozzle solution, an experimental campaign is carried out to analyse the powder spread at 15 mm from the nozzle outlet (i.e., the location of the deposition plane). The chosen carrier and shielding inert gas is argon. The metal powder employed during this experimental investigation is a gas atomized Ti-6Al-4V powder (EOS GmbH-Electro Optical Systems) with a particle size distribution ranging between 45 μm and 105 μm (d50: 70 ± 5 μm, as mentioned by the powder supplier).
As shown in
Figure 3, in order to evaluate the performance of the nozzle solution, three different process conditions are considered:
in the first one, only one nozzle is active without the presence of the substrate (see
Figure 3a);
in the second one, two nozzles are active without the presence of the substrate (see
Figure 3b);
in the third one, two nozzles and a flat substrate placed at 15 mm from the centre of the nozzle outlet are present. In this case, a flat substrate is included in order to analyse the effectiveness of the provided nozzle solution in limiting the spread of the powder particles at the first layer deposition (see
Figure 3c).
In particular, experimental tests are carried out with various carrier, shielding, and powder mass flow rates applying a full factorial 2 × 3 × 2 design, according to the Design of Experiments (DOE) approach, as listed in
Table 1. Each value range is fixed starting from process parameters commonly adopted in DED [
9]. The experimental results are investigated, taking into account a confidence interval of 95% and running three repetitions for each parameter setting.
A high-speed acquisition camera records the variation of the powder flow for every combination of process parameters, and the shape and width variation of the powder flow is detected and measured through image analyses performed with the open software ImageJ [
18] (see
Figure 4).
To correctly analyse the powder flow profile, reducing the presence of undesired floating particles which can deteriorate and worsen the powder flow characterization, 10 images are extracted from each experimental video and subsequently averaged (see
Figure 4a). From the resulting image, the corresponding background is subtracted (see
Figure 4b) and a “binary-mask” filter is applied to highlight the contour and the shape of the powder flow (see
Figure 4c). To characterize the width variation of the powder flow along its extension, 10 equally spaced evaluation planes are traced below the nozzle outlet (i.e., the red lines in
Figure 4c). The 95% of the portion of plane crossing the powder flow is taken into account as the local width of the powder flow, cutting the total powder flow intensity in correspondence to the analysed evaluation plane of 2.5% at its beginning and end. This is required due to the irregular shape of the powder flow along its edges.
The binary-mask filter converts an image in black and white thanks to the preliminary definition of an intensity threshold. This tool is very useful in image processing, nevertheless, in this case, the choice of an intensity threshold is critical since it directly affects the size and extension of the powder flow (see
Figure 5a,b). To objectively determine this value, an iterative method is designed, forcing an equal total intensity of the powder flow through the last five evaluation planes (see
Figure 6). The purpose of this constraint is to suppose an equal quantity of powder particles passing through the evaluation plane for an equal value of total intensity in correspondence with each plane. The equation employed in this case is listed below:
where
Ik is the corresponding total intensity of the powder flow for each evaluated plane. For each experimental video, the intensity thresholds chosen for the binary-mask filter are those satisfying Equation (1) (see
Figure 6).
3. CFD Modeling and Theoretical Assumptions
The deposition simulation is approached by using CFD analysis, as generally occurs for DED processes, modelling particle powders as a discrete phase dispersed into a continuous phase (i.e., the inert gas). This approach is widely employed in the scientific literature to simulate fluid-dynamics problems when a low volume fraction of powder particles is dispersed and dragged by a gas. Carrier and shielding gases are computed as a continuous phase applying the standard κ-ε turbulent flow model available in ANSYS® FLUENT (v15.0) and based on the Navier-Stokes equations. To describe a turbulent flow, the time-averaging governing equations are:
conservation of mass:
where
ρ is the argon density (1.623 kg/m
3),
ui is the gas velocity, and
xi is the gas position.
Conservation of momentum:
where
p is the pressure,
g is the gravitational acceleration,
μ is the molecular viscosity (2.125 × 10
−5 kg/m s), and
μt is the turbulent viscosity.
The most commonly used model to handle this situation is the standard
k-
ε model, in which
k and
ε represent the turbulent kinetic energy and the dissipation of kinetic energy respectively. The conservation of the kinetic energy of turbulence is given by:
whereas the conservation of the dissipation of kinetic energy of turbulence is defined as:
where
C1ε = 1.44,
C2ε = 1.92,
k = 1.0, and
ε = 1.3 are empirical constants;
Prt is the turbulent Prandtl number;
Gk is the generation of turbulence kinetic energy due to the mean velocity gradients; and
Gb is the generator of turbulence kinetic energy due to buoyancy.
The equations provided above govern the continuous phase constituted by the argon carrier gas. To complete the theoretical basis of this CFD analysis, the equations governing the secondary phase have to be provided. The discrete phase representing the powder particles dispersed into the continuous phase is computed by ANSYS
® FLUENT, which integrates the differential equation of a particle’s force balance in a Lagrange coordinate system. The balance of the forces is given by:
where
up is the particle velocity,
u is the fluid phase velocity,
ρ is the fluid density,
ρp is the density of the particles,
g is the gravitational acceleration, and
Fi is an additional acceleration (force/unit particle mass) term. The
FD coefficient is the drag force per powder mass unit and it can be calculated as:
In the equation defining
FD,
μ is the molecular viscosity of the fluid,
dp is the particle diameter,
Re is the relative Reynolds number, and
CD is the drag coefficient, defined as:
where
a1,
a2, and
a3 are empirical constants. The second term on the right of Equation (8) consists of the gravity and buoyancy forces per unit particle mass. Therefore, the particle velocity can be acquired as:
The particle trajectory can be obtained using Equations (8) and (12).
The domain taken into account during the numerical simulation is composed by a control volume large enough to allow the complete development of the powder flow at the exit of the deposition nozzle. In the case of two active nozzles, a symmetry condition is assumed. The theoretical assumptions taken into account for the numerical analysis are:
only the forces of drag, inertia, and gravity are included in the analysis;
collisions among particles are not considered;
the grain size distribution is considered uniform with an average diameter of 70 μm, since it corresponds to the mass median diameter of the powder particle distribution;
the gas-powder flow is assumed to be a steady state flow;
the powder particles are assumed to be spherical in shape;
the substrate “traps” the powder particles reaching the surface.
4. Results
Table 2 summarizes the results concerning the 95% of the portion of the total powder flow width at 15 mm below the nozzle outlet for the two process conditions where the presence of the substrate is not taken into account (see
Figure 3a,b). In the case of only one active nozzle (
Figure 3a), the effects of the shielding gas (
p-value = 0.000) and of the powder feed rate (
p-value = 0.001) are significant. Increasing the powder feed rate from 0.1 to 0.14 g/s, the 95% of the total spread of the powder flow at 15 mm decreases from 4.58 to 4.43 mm in average, obtaining an average absolute variation of −9.8%. This behaviour is significant and demonstrates the effectiveness of the suggested nozzle solution in the powder particle deposition, since it is capable of ensuring a high particle powder concentration focused at the centre of the flow, thus decreasing the width of the powder flow (i.e., the powder flow width at 15 mm from the nozzle outlet).
The influence of the shielding gas on the powder flow profile is consistent for every evaluated combination of process parameters and it is well represented in
Figure 7. The external flow of inert gas compacts the powder particles at the exit of the nozzle and contrasts the effect of the gravity force that tends to pull down the powder particles, deflecting them from the nominal trajectory imposed by the nozzle geometry. Upon increasing the shielding mass flow rate from 0 to 2.92e
−0.4 kg/s, the spread of the powder flow at 15 mm from the nozzle outlet decreases to about −14.6% on average.
In the case of two active nozzles without the presence of the substrate, analysis of variance (ANOVA) demonstrates that the powder feed rate has a negligible effect (
p-value = 0.064) on the spread variation of the powder flow compared to the shielding (
p-value = 0.001) and carrier gas (
p-value = 0.000). Upon increasing the carrier gas from 9.02e
−0.5 to 1.08e
−0.4 kg/s, the spread of the powder flow increases (see
Table 2). This is mainly due to the increase in the carrier inertia that limits the effect of the shielding gas. Moreover, the influence of the shielding gas is not linear, as happens in the process condition with only one active nozzle. In this case, the experimental analysis shows a significant reduction in the powder particles diffusions for medium values of the shielding gas, followed by a critical enlargement at the highest values. For medium values of the shielding gas, in fact, the average 95% of the powder flow width at 15 mm decreases down to −11.8%, whereas for higher mass flow rates the width only decreases down to −5.6%. The reason for this behaviour can be related to the location at which the two powder flows meet each other. Indeed, upon increasing the shielding mass flow rate from 3.9e
−0.4 kg/s to 5.84e
−0.4 kg/s, the zone where the two powder flows meet each other ends up being lower compared to the previous cases, moving away from the nozzle outlet. This phenomenon affects the measurements since for high values of the shielding mass flow, the deposition plane is no longer located at 15 mm from the nozzle outlet.
Concerning the case of two active nozzles with the presence of the substrate at 15 mm from the nozzle outlet, the acquired images are strongly affected by a high concentration of bouncing powder particles that degrade the image analysis, making it impossible to define a reliable image intensity threshold and preventing the application of the image-based method previously discussed. For this reason, only a qualitative analysis is possible for comparing the different spectra of the powder particles distribution in correspondence to the surface of the substrate.
Figure 8 illustrates the effect of the shielding gas for constant values of the carrier mass flow rate and the powder feed rate. In the case of no shielding (see
Figure 8a), the area with the larger particle mass concentration (i.e., the red one) is wider compared with those represented in
Figure 8b,c, in which the shielding gas is 3.9e
−0.4 kg/s and 5.84e
−0.4 kg/s, respectively. In particular, when no shielding is applied, the top edges of the area with the maximum powder concentration are indented, indicating a strong rebound of the powder particles on the metal surface of the substrate. On the contrary,
Figure 8b,c show a red area with more defined edges and then a limited rebound of the powder particles. The zone with the higher particle concentration seems to be qualitatively smaller for higher values of the shielding mass flow rate (see
Figure 8c). Nevertheless, to have a feedback and a quantitative analysis of the influence of the shielding, carrier, and powder mass flow on the deposition efficiency of the process under this process conditions, further experimental tests are required.
5. Discussion
The CFD model is verified, estimating the spread of the powder flow at 15 mm from the nozzle outlet and taking into account the particle mass concentration (see
Figure 9a). To correctly compare the numerical outputs with the experimental results, the nominal particle mass concentration is adequately filtered for each combination of process parameters, since:
the experimental analysis is based on 2D images (as those reported in
Figure 4);
the employment of a “binary-mask” filter requires the application of an intensity threshold to highlight the shape of the powder flow;
the particle mass concentration is provided by the numerical software with no filtering.
The filtering threshold set to cut off the numerical outputs is equivalent in terms of percentages to the intensity threshold employed during the image analysis. After filtering, 10 lines are transversely traced to discretize the particle mass concentration and detect the spread of the powder flow at 15 mm from the nozzle outlet (as shown in
Figure 9b). The numerical powder flow width to compare with the experimental one is estimated by computing the 95% of the total particle mass concentration resulting from the sum of the 10 tracked lines, cutting off the values lower than the 5% of the maximum powder concentration.
The numerical results do not correspond to the experimental ones; from data reported in
Table 3 it is possible to note that the estimated powder flow width is quite different from the results obtained by the experimental characterization.
For both one active nozzle and two active nozzles systems, the employed CFD model is not capable of correctly estimating the powder behaviour. From an analysis of the CFD results, for each combination of process parameters, the shielding gas seems to have a very low influence on the powder flow geometry (see
Figure 10). Indeed, both in the case with no shielding (see
Figure 10a) and in the case with shielding (see
Figure 10b,c), the graphical CFD results show the presence of secondary isolated powder flows that diverge from the central main powder flow. These secondary flows seem to not be affected by the presence of the shielding that should deflect them or limit their extension.
The motivations for such different results could be attributed to two different reasons:
The employment of a correct value to filter all the data provided by the CFD analysis is a critical issue not easy to deduce, but it is required to compare the numerical results with the experimental ones due to the application of a binary-mask filter during the experimental image analysis. Nevertheless, the numerical method employed in this analysis seems to not fit the behaviour of the powder flow recorded during the experimental investigation, probably because both the influence of the air at the exit of the deposition nozzle is not taken into account and the powder grain size distribution is not correctly estimable by its mass median diameter (i.e., 70 μm). To improve the CFD simulation and to obtain a powder flow behaviour more in compliance with the real one, a more complex model has to be taken into account, such as the Eulerian multiphase model. The Eulerian multiphase is a model implemented in ANSYS® FLUENT, which allows for the modeling of multiple separate, yet interacting phases, where Eulerian treatment is used for each phase, in contrast to the Eulerian-Lagrangian treatment that is used only for the discrete phase model. The application of this numerical model together with a more accurate estimation of the distribution of the powder grain size (i.e., Rosin-Rammler particle size distribution) could improve the experimental fitting of the powder behaviour at the expense of more computational time and memory.