Next Article in Journal
Characterization of the Internal Stress Evolution of an EB-PVD Thermal Barrier Coating during a Long-Term Thermal Cycling
Previous Article in Journal
Synthesis and Study of the Optical Properties of a Conjugated Polymer with Configurational Isomerism for Optoelectronics
Previous Article in Special Issue
Influence of Latent Heat of Fusion on the Melt Pool Shape and Size in the Direct Laser Deposition Process
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Fast Prediction Model for Liquid Metal Transfer Modes during the Wire Arc Additive Manufacturing Process

1
Institute of Advanced Structure Technology, Beijing Institute of Technology, Beijing 100081, China
2
School of Mechanical Engineering, Beijing Institute of Technology, Beijing 100081, China
*
Authors to whom correspondence should be addressed.
Materials 2023, 16(7), 2911; https://doi.org/10.3390/ma16072911
Submission received: 4 March 2023 / Revised: 28 March 2023 / Accepted: 3 April 2023 / Published: 6 April 2023

Abstract

:
The liquid metal transfer mode in wire arc additive manufacturing (WAAM), plays an important role in determining the build quality. In this study, a fast prediction model based on the Young–Laplace equation, momentum equation, and energy conservation, is proposed, to identify the metal transfer modes, including droplet, liquid bridge, and wire stubbing, for a given combination of process parameters. To close the proposed model, high-fidelity numerical simulations are applied, to obtain the necessary inputs required by the former. The proposed model’s accuracy and effectiveness are validated by using experimental data and high-fidelity simulation results. It is proved that the model can effectively predict the transition from liquid bridge, to droplet and wire stubbing modes. In addition, its errors in dripping frequency and liquid bridge height range from 6% to 18%. Moreover, the process parameter windows about transitions of liquid transfer modes have been established based on the model, considering wire feed speed, travel speed, heat source power, and material parameters. The proposed model is expected to serve as a powerful tool for the guidance of process parameter optimization, to achieve high-quality builds.

1. Introduction

Wire arc additive manufacturing (WAAM) technology has gained increasing attention thanks to its high deposition rate, low equipment and material consumption, and large forming size [1,2]. However, it still suffers from difficulties in controlling morphology, including hump defects, surface roughness, and deviations in shape and size, and difficulties in controlling mechanical properties, such as coarse microstructures caused by high heat input [3,4]. The metal transfer behavior is one of the most important issues affecting the morphology and mechanical properties of parts formed by WAAM. It has been demonstrated that unstable metal transfer can lead to uneven surfaces and internal porous defects in the solidified parts and, therefore, significantly affect the forming quality and mechanical properties [5,6].
In general, three types of metal transfer modes are observed in experiments: the droplet mode, the liquid bridge mode, and the wire stubbing mode [7,8]. Among them, the liquid bridge mode is most favorable [8,9]. At this point, the liquid metal enters the melt pool steadily and forms a liquid bridge connecting the melt pool and the wire. As a result, it leads to a smooth melt pool surface and, therefore, a better forming quality. However, achieving such a stable liquid bridge is challenging in practice. For example, when the distance between the wire and the substrate is large enough, a droplet pattern develops, causing surface ripples. Alternatively, if the distance is too small, a stubbing mode tends to develop, where unmelted wire causes severe disturbance to the melt pool, reducing the forming quality and resulting in metallurgical defects. Other factors that affect liquid metal transfer patterns include the wire feed speed, travel speed, and input power, among others. For controlling the metal transfer pattern, these process parameters must be carefully selected, based on an accurate model, considering multiple physical fields and violent multi-phase flow.
Much work has been devoted to the metal transfer modes in wire feed additive manufacturing processes. Abioye et al. [8] analyzed the influence of heat input, travel speed, and wire feed speed on the deposition cladding layer, and established the process parameter window of energy per unit length of track and wire deposition volume per unit length of track, to predict the liquid bridge transfer and droplet transfer modes. Luo et al. [10] used the arc information in the manufacturing process to identify the droplet transfer mode, and established a quadratic polynomial relationship between droplet transfer frequency and arc power. Scotti et al. [11] addressed the importance of the distance between the wire and substrate in controlling different metal transfer behaviors. Wu et al. [12] reported that when the distance between the electrode and the substrate increases, the liquid bridge transfer mode changes to the droplet transfer mode. Zhao et al. [13] pointed out that the liquid bridge mode had good comprehensive performance, while the wire stubbing mode was unstable and could cause unmelted defects.
To better understand the physical mechanism of heat transfer and fluid flow behavior during metal transfer, many numerical investigations have been conducted. Tang et al. [14] used the level-set method to identify the free surface, and investigated the periodic impact of droplets on the melt pool. Hu et al. [15] investigated the liquid bridge transfer mode, considering the influence of the wire feed. They proposed a dimensionless parameter called “slenderness number”, to roughly estimate whether the liquid bridge can maintain a stable configuration. Bishal et al. [16] used a volume of fluid (VOF) model to study the metal transfer behaviors and presented a corrected formula of droplet detachment in the WAAM process. Chen et al. [17] numerically studied metal transfer modes at different wire feed speeds, and found that the periodic flow pattern of the melt pool, caused by metal transfer impact, leads to ripple and even hump defects, resulting in uneven forming quality. Despite the effectiveness in analyzing metal transfer behaviors, the numerical models above are limited by high computational cost, especially when used in the optimization process. Therefore, a proper model to rapidly predict metal transfer mode, is urgently needed.
This work proposes a prediction model to identify the metal transport mode (droplet, liquid bridge, or wire stubbing) in the WAAM problem. Based on the balance of surface tension and pressure difference, the model can provide the morphology of the liquid bridge and the critical condition for turning into droplet mode. Derived based on energy conservation, the model can also predict whether the liquid bridge mode will turn into a wire stubbing mode. The prediction model requires only a few high-fidelity numerical simulations to calibrate the input parameters, and can accurately predict similar problems for a given combination of process parameters. Validation cases suggest that the prediction model is computationally efficient compared to numerical simulations, with sufficient accuracy, which is very beneficial for optimizing process parameters. Finally, a series of calculations are performed based on the proposed model. A window of process parameters for the WAAM problem is presented, considering wire feed speed, travel speed, heat source power, and material parameters. In contrast to high-fidelity numerical models, the presented approach offers significantly higher efficiency, since it does not need to solve the entire deposition process. Compared to experiments, it can predict liquid metal transfer modes with acceptable accuracy, in an economical and efficient manner. Therefore, the proposed model can serve as a powerful tool for optimizing process parameters when a large number of liquid transfer mode predictions are necessary.

2. Prediction Model for Liquid Metal Transfer Modes

2.1. Brief Introduction to Metal Transfer Modes

In WAAM, the liquid metal transfer mode plays an important role in the melt track’s morphology and the component’s fabrication quality. The three types of transfer mode, droplet, liquid bridge, and wire stubbing, are illustrated in Figure 1, and briefly introduced as follows. In the droplet transfer mode, the liquid metal periodically enters the melt pool as droplets. The liquid bridge mode denotes that the wire and substrate are bridged by continuous liquid metal. In the wire stubbing mode [18,19,20], the solid wire fails to completely melt before entering the melt pool, leading to severe disturbance to the melt pool. Generally, the liquid bridge mode is most favorable, while the other two modes, especially the wire stubbing mode, need to be avoided. The following section will analyze the liquid bridge mode and determine the critical conditions for its transformation into the droplet and the wire stubbing modes.

2.2. Profile of the Liquid Bridge

In the liquid bridge mode, the profile of the liquid metal is shown in Figure 2. The shape of the liquid bridge is mainly dominated by surface tension, gravity, and pressure. To describe the continuous liquid bridge, let us first assume that the input power of the arc is sufficient to melt the feed wire into liquid completely. Meanwhile, the shape of the liquid bridge is considered to be axisymmetric. As the Marangoni force acts mainly in the tangential direction, its effect on the liquid bridge profile is exerted indirectly, by influencing the flow velocity, which is ignored here. The pinch effect caused by the magnetic field is not considered. As a result, the surface tension in the normal direction mainly determines the profile. The profile can thus be described by the Young–Laplace equation [21]:
σ κ = Δ p
where κ is the curvature, σ is the surface tension coefficient, and Δ p is the pressure difference across the gas–liquid surface. For a body with cylindrical symmetry, the curvature reads:
κ = 1 R 1 + R 2 1 / 2 R 1 + R 2 3 / 2
where R is the radius and superscripts and denote first- and second-order derivatives, respectively.
To simplify the modeling, the fluid momentum conservation is written as the following formulation, considering the liquid to be stable:
p = ρ g
where ρ is the density and g is the gravitational vector. Assuming the pressure at the top of the liquid radius is atmospheric pressure p atm , Equation (3) becomes:
p = p 0 + ρ g ( H z )
where p 0 is atmospheric pressure, z is the height coordinate with z = 0 at the substrate surface, g is the gravitational acceleration in the z direction, and H is the height of the liquid bridge.
Substituting Equation (4) and the curvature Equation (2) into Equation (1), yields the description of the liquid bridge profile as:
R 1 + R 2 3 / 2 1 R 1 + R 2 1 / 2 + ρ g ( H z ) σ = 0
The boundary conditions associated with Equation (5) are as follows:
R z = 0 = R 0 d R d z z = 0 = cot ϕ 0 R z = H = R w
where R 0 is the bottom radius of the liquid bridge, ϕ 0 is the angle between the liquid bridge profile and substrate, and R w is the projection radius of the wire in the z = 0 plane. For the studied case, R w reads:
R w = R c / sin θ
where R c is the wire radius, and θ is the angle between the feed wire and the substrate.
Since the solution to Equation (5) cannot be obtained analytically, we proposed Algorithm 1, an iterative approach, to solve it. In the proposed method, R 0 and ϕ 0 are chosen to be input parameters, to obtain the unknown liquid bridge height, H. To calibrate the values of R 0 and ϕ 0 , a high-fidelity numerical model is applied. Since the surface tension coefficient changes with temperature, σ is also calibrated using the numerical model. By combining bisection and quadrature methods, as listed in Algorithm 1, Equation (5) can be solved with these input parameters. It is worth noting, that the proposed method converges quickly when the appropriate input parameters are selected.
Algorithm 1: Algorithm to the liquid bridge profile.
Materials 16 02911 i001

2.3. The Transition from Liquid Bridge to Droplet Mode

In the proposed model, above, a liquid bridge is assumed to form between the wire and the substrate. However, as the distance between the wire and the substrate increases, the liquid bridge mode will transition to the droplet mode. Algorithm 1 is used to determine the height of the liquid bridge. The critical condition for the transition to the droplet mode can be expressed as follows.
H 0 H
where H 0 is the initial distance between the wire tip and the substrate. If Equation (8) is not met, the liquid bridge mode may transition to the droplet mode.
In the droplet mode, the liquid profile is shown in Figure 3, and it is also governed by Equation (5). However, as the liquid becomes a pendant droplet, the boundary conditions change to:
R z = 0 = 0 R z = H = R w
The bottom tip of the droplet is parabolic, and the two principal curvature radii are equal. So near the tip, one has:
z = 1 2 R 2 r 0
where r 0 is the tip radius of the droplet.
The profile of the pendant droplet is obtained through a scheme similar to Algorithm 1. To avoid the singularity, Equation (10) is used as the initial value of numerical integration. Once the droplet profile is predicted, the maximum volume of droplet growth can be achieved as:
V d = π 0 H R 2 ( z ) d z
with V d being the maximum volume of each droplet.
Assuming the arc power is sufficient to melt the wire completely, we can calculate the wire melting volume rate V ˙ w as:
V ˙ w = π R c 2 v w
with v w being the wire feed velocity. Therefore, the frequency of the liquid drop is:
f = V ˙ w V d
In turn, as the droplet frequency increases, the droplet mode gradually shifts to the liquid bridge mode, resulting in improved quality of formation.

2.4. Transition from Liquid Bridge to Wire Stubbing Mode

In the liquid morphology model presented above, it is assumed that the arc power is sufficient to fully melt the wire. However, in practical applications, this is often not the case. If the arc power is too low or the wire feed speed is too fast, it can result in the wire not melting completely and entering the melt pool, which is referred to as wire stubbing mode. In this mode, the tip of the wire collides with the solid part at the bottom of the melt pool, causing significant disruption to the melt pool flow field. As a consequence, irregularities in the solidified track’s morphology can occur. To predict the wire stubbing mode, the following formulations are derived.
For the WAAM process, the wire tip and the heat source affected area are illustrated in Figure 4. Despite the complexity of the arc heat source, it can be approximated into a Gaussian surface heat source, based on experimental and numerical investigations [22,23,24]. In this study, the heat source is defined as:
q input = η P l π r 2 exp 2 x 2 + y 2 r 2
where η is the absorption coefficient, and r is the equivalent radius of the heat source, which is a function of current [25,26,27,28]:
r = 0.533 I 0.2941
with I being the current. The heat source power P l is determined as:
P l = I U
with U being the voltage.
The total power absorbed by the wire can be calculated by integrating the affected area:
E absorb = q input d S = η P l π r 2 exp 2 x 2 + y 2 r 2 d S
Assuming the temperature of the preheated wire to be T p , if the wire is fully melted, the temperature at the wire cross-section needs to reach the solidus temperature, T s . The power needed for melting satisfies:
E melt = ρ π R c 2 v w c p T s T p
where c p is the specific heat capacity.
According to the conservation of energy, the condition for the complete melting of the wire is:
E melt E absorb
Otherwise, the liquid bridge mode will change to the wire stubbing mode.

2.5. High-Fidelity Numerical Model for Model Calibration

In the above model, three input parameters must be calibrated using numerical simulations, including the bottom radius of the liquid bridge R 0 , the angle ϕ 0 , and the surface tension coefficient σ , as a function of temperature. The high-fidelity numerical model is set up as follows.
The metal is assumed to be an incompressible mixture in the numerical simulations. The mass conservation equation reads:
· u = 0
where u is the fluid velocity. The momentum conservation equation reads:
ρ u t + ρ u · u = · ( μ u ) + ρ g p ρ g β T T L + K 0 1 f l 2 f l 3 + B u + b
where μ is the dynamic viscosity, β is the fluid thermal expansion coefficient, and T L is the liquidus temperature. K 0 is the Carman–Kozeny coefficient of the mushy zone flowing through a porous media [29], f l is the liquid volume fraction, B is a small constant to avoid the singularity, and b is the body force. The momentum equation considers the buoyancy force under Boussinesq’s assumption and Darcy’s damping effect.
At the gas–liquid interface, the surface tension and Marangoni forces are defined as follows:
f surface = σ n κ + d σ d T [ T n ( n · T ) ]
where n is the normal vector of the surface, and d σ d T is the of surface tension gradient.
The arc pressure at the gas–liquid interface is assumed to follow a Gaussian distribution, as follows [27]:
p arc = μ m I 2 8 π σ a 2 exp x 2 + y 2 2 σ a 2
where μ m is the magnetic permeability, and σ a is the Gaussian pressure parameter.
The recoil pressure caused by metal evaporation is as follows:
p recoil = 0.54 p 0 exp L v M T T v R g T T v
where L v is the latent heat of vaporization, M is the molar mass, and R g is the universal gas constant. The surface tension, arc pressure, and recoil pressure are converted into body force b , by using the volume fraction.
The free surface evolution is solved by the volume of fluid (VOF) method, as follows:
f m t + ( u · ) f m = 0
where f m is the volume fraction of the metal.
The energy conservation equation reads:
t ( ρ h ) + u · ( ρ h ) = · ( k T ) t ( ρ Δ h ) u · ( ρ Δ h ) + Q
where h is the enthalpy, k is the heat conductivity, Δ h is the latent enthalpy change of the alloy, and Q includes the volume heat source and losses. The arc heat source is given in Equation (14). The heat losses are defined as follows.
q loss = ε σ b T 4 T a 4 + h c T T a + 0.82 L v M 2 π M R g T p 0 exp L v M T T v R g T T v
The right-hand side of Equation (27) includes the radiation, convection, and vaporization heat losses. In Equation (27), ε is the emissivity, σ b is the Stefan–Boltzmann constant, T a is the ambient temperature, and h c is the convection coefficient.
The finite volume method (FVM) is used to solve the model. The momentum equation is solved using an operator-splitting scheme [30]. A second-order upwind scheme is used for the convection term [31], and the first-order Euler method is applied for time integration [32]. The validation of the high-fidelity numerical simulation model is detailed in Appendix A (see Figure A1). The numerical results are in good agreement with the experimental data. Once the numerical results are obtained, they can calibrate the bottom radius R 0 , the angle ϕ 0 , and the surface tension coefficient σ .

2.6. Procedures of the Prediction Model

For the proposed fast prediction model, the procedures to predict the liquid metal transfer mode are as follows.
  • Calibrate the bottom radius R 0 , the angle ϕ 0 , and the surface tension coefficient σ , by high-fidelity numerical simulations.
  • Calculate the liquid bridge profile and height H, using Equation (5).
  • If Equation (8) is unsatisfied, the initial distance between the wire and the substrate exceeds the liquid bridge height limit, and the liquid bridge will develop into the droplet transfer mode. Calculate the dripping frequency of droplets using Equation (13).
  • If the energy conservation condition in Equation (19) is not satisfied, the liquid bridge will develop into the wire stubbing mode.

3. Results and Discussion

3.1. Model Calibration by High-Fidelity Numerical Simulations

Before using the proposed model to predict the liquid metal transfer mode, the numerical examples shown in Figure 5 are conducted, to calibrate the prediction model.
In this case, x, y, and z are the travel, the transverse, and the height directions, respectively. The size of the calculation domain is 3.6 cm × 1.2 cm × 1.4 cm. The material used for both substrate and wire is mild steel, with its material parameters [33,34] listed in Table 1 and Table 2. The thermal conductivity and specific heat capacity are temperature-dependent and are linearly interpolated from the values listed in Table 2, while other parameters remain constant.
Three cases are conducted, with a constant travel speed of 1 cm/s, and wire feed speeds of 2 cm/s, 3 cm/s, and 4 cm/s, respectively. The selected wire feed speeds are in the reasonable range of experimental conditions for WAAM. The arc heat source power is 1440 W and the wire diameter is 1.2 mm. The cell size is 0.01 cm, leading to 6,048,000 cells. The temperature distribution, with respect to the z coordinate on the liquid bridge, is fitted according to an average of the results of the three examples. The second-order Gaussian function is used to fit the curve, and the coefficient of determination (denoting the goodness of fit), is 0.99. The fitting results are as follows.
T ( z ) = a 1 e z b 1 / c 1 2 + a 2 e z b 2 / c 2 2
where a 1 = 546 , b 1 = 0.0837 , c 1 = 0.01505 , a 2 = 1882 , b 2 = 0.1095 , and c 2 = 0.8174 . Since a suitable dataset for temperature-dependent surface tension coefficients of the liquid metal is not currently available, a linear fitting is used to obtain the surface tension coefficient:
σ = σ 0 + d σ d T T ( z ) T L
where σ 0 is the surface tension coefficient at the liquidus temperature T L , and d σ d T is given as a constant [34] in Table 1. By substituting Equation (28) into Equation (29), the surface tension used in Equation (5) is calibrated.
Additional cases are required to calibrate the bottom radius R 0 and angle ϕ 0 as input parameters in the prediction model. This is because, unlike the surface tension coefficient, which is not very sensitive to the process parameters, R 0 depends on the wire feed speed and travel speed. Therefore, nine cases were carried out to calibrate the two parameters. The setups for these nine cases are listed in Table 3, while others are the same as the above case.
The R 0 in Table 3 was obtained from each numerical case by taking the average over a period of time. Ten different time instances were selected, and the statistical results were regressed. Here, parameter ξ is introduced, which is defined as the ratio of the wire feed speed to travel speed. A regression relationship between ξ and R 0 can be established:
R ( ξ ) = a ξ 2 + b ξ + c
where the parameters a, b, and c are 2.3 × 10 3 , 2.5 × 10 2 , and 5.5 × 10 2 , respectively. Equation (30) is then used to calibrate the liquid bridge bottom radius in the liquid bridge profile, Equation (5).
The same approach is used for the calibration of angle ϕ 0 , by averaging the numerical results at various time instances. However, from the results of the high-fidelity numerical model, it is found that ϕ 0 is not sensitive to the values taken for the wire feed speed and travel speed. In the given range of parameters, the average of ϕ 0 is around 10.1 degrees, with no significant deviation, and thus is chosen as the input of Equation (5). Considering that the predictions in the two cases are insufficient to illustrate the model’s prediction capability, more predictions are shown in Section 3.4, in the process parameter window investigation.

3.2. Liquid Bridge Profile and Critical Height Validation

Two cases are presented to validate the profile predicted by the proposed model. In the first case, the wire feed speed is 2 cm/s, and the calibrated bottom radius R 0 is 0.12 cm. In the second case, the wire feed speed is 1.34 cm/s, and the calibrated bottom radius R 0 is 0.10 cm. Other setups are the same as the numerical cases used in Section 3.1.
The high-fidelity numerical and proposed prediction models’ results are presented in Figure 6. It is demonstrated that the predicted profiles are similar to the results obtained by the high-fidelity numerical model. The profile of the neck region of the liquid bridge deviates slightly compared to the numerical results. This deviation may be because the proposed prediction model neglects some surface forces acting on the fluid. However, the overall profile and bridge height results agree well with the numerical results. For the liquid bridge heights, the height obtained by the prediction model is 0.10 cm in the first case, while the height calculated by the high-fidelity numerical simulation is 0.09 cm and the relative error is 11.1%. In the second case, the height of the liquid bridge calculated by the proposed model is 0.085 cm, while the height achieved by the numerical model is 0.08 cm and the relative error is 6.3%.
For the predicted liquid bridge profiles, the critical condition for the occurrence of the droplet mode can be determined. According to Equation (8), if the distance between the substrate and wire tip, H 0 , exceeds the liquid bridge height H, the liquid bridge changes to the droplet mode. Two cases are presented here, to verify the proposed model’s capability in predicting the liquid bridge’s critical height. In the first case, the initial distance between the substrate and wire tip, H 0 , is 0.12 cm. In the second case, H 0 is 0.14 cm. The wire feed speed is 4 cm/s in both cases. The numerical results are shown in Figure 7.
The numerical results suggest that the first case is droplet mode, and the second is liquid bridge mode. From the proposed model, the critical height H, is 0.12 cm, and the two cases should be droplet and liquid bridge modes, respectively, according to Equation (8). A stable liquid bridge cannot be formed when the distance between the wire and the substrate exceeds the predicted maximum liquid bridge height. It is demonstrated that the predicted results agree well with the high-fidelity numerical simulations.

3.3. Droplet Profile and Frequency Validation

The proposed model can be used to predict the droplet profile and dripping frequency for the droplet mode. Therefore, two cases are presented here to verify the droplet profile and dripping frequency.
The first case is the dripping of water droplets. The predicted results have been compared to the experimental profiles by Zhang et al. [35]. In the water droplet experiment, the water density is 1000 kg/m 3 , the surface tension coefficient is 0.0728 N/m, the outer radius of the tube is 0.16 cm, and the tube flow rate is 1 mL/min. The profiles of droplets at different time instances are shown in Figure 8, where the predicted profiles obtained by Equation (5) are in good agreement with the experimental data [35]. The comparison shown in Figure 8, demonstrates that the proposed model can predict the droplet profiles accurately at different time instances.
As detailed in Appendix B, we conducted WAAM experiments to investigate the dripping frequency. The experiments were conducted using the process parameters listed in Table A2, with a duration of 1 min. The results obtained are plotted in Figure 9, which reveals that all three experimental cases exhibited droplet modes. In this mode, the liquid metal droplet grows gradually as the wire is fed. However, when the surface tension is insufficient to maintain the droplet shape, it separates from the wire, causing the periodic impact of these droplets to yield ripples on the solidified specimen. Since the droplets cannot be observed directly during the experiment, their number and dripping frequency are calculated indirectly, based on the spacing of the ripples on the solidified specimen. A comparison between the experimental and theoretical results is presented in Table 4. The experimental findings indicate that both the number of droplets and dripping frequency increase with the wire feed speed, which is consistent with Equation (13). The predictions are in good agreement with the experimental trend, although the number of droplets and dripping frequency from the experiment surpasses the predicted values. The relative error of the dripping frequency ranges from 10.2% to 18.9%, which may be attributed to unstable conditions during the experiment, such as arc pressure oscillation and equipment vibration. Additionally, the neglected electromagnetic forces and pinch effects in the prediction model could have accelerated the dripping frequency of droplets. Despite these discrepancies, the proposed model still fulfills the process parameter control requirements to a certain extent.

3.4. Process Parameter Window Considering Liquid Bridge–Droplet Transition

In the WAAM process, the liquid bridge mode is preferred, but identifying the critical conditions for the transition between the liquid bridge and droplet modes through experiments can be expensive. High-fidelity numerical simulations are also computationally intensive. For example, the computational time of a typical numerical case established in this study is approximately 24 h, using a workstation with 40 cores CPU (Intel Xeon 6248). Not to mention, that at least dozens of such examples are required to obtain a process parameter window. Therefore, the proposed fast prediction model for obtaining a process parameter window is crucial.
With the proposed prediction model, it is possible to determine the critical condition for the droplet mode, as shown in Equation (8). This model allows for the presentation of a process parameter window that considers the liquid bridge–droplet transition, with a focus on two key process parameters: the initial distance between the substrate and wire tip ( H 0 ) and the wire feed speed to travel speed ratio ( ξ ). The process parameter window generated by the proposed model is illustrated in Figure 10. The process parameters can be clearly divided into two zones in the figure: the upper area represents the droplet mode zone, and the lower area represents the liquid bridge mode zone. The black line is obtained using Equation (5), while the shadowed error bar zone is obtained by selecting different calibrated R 0 values, taking into account the standard deviation in Table 3.
A total of eleven high-fidelity numerical cases were conducted, to verify the accuracy of the proposed prediction model, with their setups summarized in Table 5. Figure 10 presents a comparison of the numerical simulation results and the predictions of the proposed model, where the triangles and circles represent the droplet and liquid bridge modes, respectively. From the comparison, one can find that the predictions from the proposed model are in good agreement with their counterparts from the numerical simulations. However, the proposed model takes only a few minutes to complete the prediction, while high-fidelity numerical simulations often require a couple of days. The proposed model is also compared with experiments conducted by Wang et al. [36]. Fourteen experimental results, using wires with a diameter of 2.4 mm, were chosen for comparison. In these experiments, the distance H 0 , between the substrate and wire, ranged from 0–3 mm, and the wire feed speeds were 0.8 m/min and 1.1 m/min. Details of the setups can be found in Figure 17 in [36]. Figure 11 presents a comparison between the experimental results and the predictions of the proposed model, with triangles, circles, and stars representing the droplet, liquid bridge, and intermediate states observed from experiments, respectively. The predictions are in good agreement with experiments. It is demonstrated that the proposed model can efficiently predict the critical conditions of the liquid bridge–droplet transition and provide guidance for practical engineering applications.

3.5. Process Parameter Window Considering Liquid Bridge–Wire Stubbing Transition

In the WAAM process, the wire stubbing mode must be avoided as much as possible, because it can lead to poor forming quality. The proposed prediction model can be used to determine the critical conditions for the wire stubbing mode, as shown in Equation (19).
Using the proposed model, we present a process parameter window that considers the transition of the liquid bridge–wire stubbing. The two most critical parameters, i.e., the arc power and the wire feed speed, are considered. The wire feed speed ranges from 0 to 9 cm/s, while the arc power ranges from 0 to 2600 W. The process parameter window depicted by the proposed model is shown in Figure 12. The figure clearly divides the process parameters into two zones. The top left zone corresponds to the liquid bridge mode, while the bottom right zone is for the wire stubbing mode.
Ten high-fidelity numerical cases were conducted, to validate the accuracy of the prediction model. The material parameters and setups used in these cases are listed in Table 1 and Table 6, respectively. Judging whether the numerical results predict the liquid bridge or wire stubbing mode in all ten examples is challenging, as it cannot be directly observed from the high-fidelity numerical simulations. Therefore, the deposition efficiency is defined, and proposed to make such judgments. This parameter represents the ratio of the actual volume increase in the specimen after solidification, to the volume rate of the wire feed, and an ideal deposition efficiency should be 100%. If the wire is not completely melted before it enters the melt pool in the numerical simulations, the deposition efficiency will be less than 100%, indicating the occurrence of wire stubbing mode. In this study, a 90% deposition rate is used as a critical condition to account for any numerical errors. All ten cases are illustrated in Figure 12, with diamonds representing wire stubbing mode and circles denoting liquid bridge mode. The proposed prediction model is shown to agree well with the numerical results, effectively predicting the critical conditions of the liquid bridge–wire stubbing transition; and so could be used to guide engineering practices. It is identified that the wire stubbing mode is less likely to occur with higher arc power and lower wire feed speed.

4. Conclusions

In this study, a fast prediction model for identifying liquid metal transfer modes in wire arc additive manufacturing is proposed. By utilizing the Young–Laplace equation, momentum equation, and energy conservation, the profile and height of the favorable liquid bridge are obtained, which is crucial in determining the critical conditions for transitioning between different modes. Specifically, the critical height of the liquid bridge is used to determine the liquid bridge–droplet transition, while the energy conservation principle is applied to derive the critical condition for the liquid bridge–wire stubbing transition. In addition, the dripping frequency of the liquid drop can be predicted. High-fidelity numerical simulations are used to calibrate the surface tension coefficient, the bottom radius of the liquid bridge, and the angle between the liquid bridge profile and substrate, to close the solution of the proposed model. Numerical cases and experiments were carried out to validate the accuracy of the prediction model, and the results show that the error ranges from 10.2% to 18.9% for the dripping frequency in the droplet mode, and from 6.3% to 11.1% for the liquid bridge height in the liquid bridge mode. Moreover, it is demonstrated that the proposed model can accurately predict the transition from liquid bridge to droplet and wire stubbing modes.
Based on the proposed prediction model, a process parameter window, considering the transition from liquid bridge mode to droplet and wire stubbing modes, is presented. The model is demonstrated to be effective and efficient in optimizing the process parameters. In the WAAM process, the choice of process parameters and the liquid metal transfer modes can have significant impacts on the quality of the final build. Therefore, it is of far-reaching significance to establish an efficient prediction model and corresponding process parameter windows in actual production. The proposed model, therefore, has practical significance in guiding the actual industrial production process and improving the fabrication quality.
In addition, the proposed model in this study has some limitations that should be acknowledged. The model assumes a steady-state liquid transport process, and considers most material properties to be constants. Potential improvements include, using temperature-dependent coefficients and considering the thermocapillary stress, e.g., the semi-analytical model for the apparent slip length of Poiseuille and Couette flows [37]. The influence of electromagnetic forces, the pinch effect, and wire melting rate considering different power sources, are also neglected. Future research could incorporate these factors into the model to enhance its accuracy and applicability.

Author Contributions

Software, Investigation, Writing-original draft, J.O.; Conceptualization, Methodology, Writing-review & editing, M.L.; Conceptualization, Methodology, Writing-review & editing, Supervision, Y.L.; Validation, Investigation, S.P.; Resources, Validation, Investigation, C.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China, under grant numbers 11972086 and 52105287, and the Beijing Institute of Technology Research Fund Program for Young Scholars, under grant number XSQD-202223002.

Data Availability Statement

The main data supporting the work are available within the article. Extra data are available from the corresponding author upon reasonable request.

Conflicts of Interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

Appendix A. Numerical Model Validation

The experiment conducted by Ogino et al. [38] was addressed, using the numerical model for its validation. The material properties used in this numerical model are summarized in Table 1, and the process parameters are set to wire feed speed of 5 cm/s, travel speed of 1 cm/s, and heat source power of 1440 W. The comparison of the deposition layer’s size and shape between the numerical simulation and experimental results is shown in Figure A1, indicating the accuracy of the numerical model. On the transverse cross-section, the average bead width and depth from the experiment are 0.46 cm and 0.26 cm, respectively. The bead width and depth from the numerical simulation are 0.49 cm and 0.24 cm, respectively. The average width and depth errors are found to be 6.5% and 7.6%, respectively, indicating that the morphology and dimensions of the bead from the numerical results agree well with their counterparts from the experimental data.
Figure A1. Transverse cross-sectional view of the molten track, where the left part is from the numerical results and the right part from the experiment [38].
Figure A1. Transverse cross-sectional view of the molten track, where the left part is from the numerical results and the right part from the experiment [38].
Materials 16 02911 g0a1
Table A1. Comparison of melt pool width and depth between experiment data and simulation results.
Table A1. Comparison of melt pool width and depth between experiment data and simulation results.
WidthDepth
Experimental0.46 cm0.26 cm
Numerical simulation0.49 cm0.24 cm
Relative error6.5%7.6%

Appendix B. Experimental Equipment

The WAAM equipment utilized in this study is shown in Figure A2a. The equipment comprises a control system, a computer numerical control (CNC) machine tool, an arc power supply machine, a tungsten inert gas torch installed on the CNC machine tool, two hot-wire resistance power sources, and two wire feeders. As depicted in Figure A2b, the wire feed tubes are situated on both sides of the torch at a 45-degree oblique downward angle, and two wire feeders manage the feed speed. The travel speed of the CNC machine tools is regulated by the motion programs written in the control system.
Figure A2. The WAAM equipment: (a) schematic diagram of WAAM equipment, (b) equipment details, (c) photo of the equipment.
Figure A2. The WAAM equipment: (a) schematic diagram of WAAM equipment, (b) equipment details, (c) photo of the equipment.
Materials 16 02911 g0a2
Regarding the power source, a direct current power source is used, with a current of I = 95 A and a voltage of U = 15.16 ± 0.6 V. The current is obtained by averaging the square wave current over time, and the voltage is measured at the tungsten inert gas torch. The polarity of the power supply is positive for the torch, and zero potential for the substrate. The power of the arc heat source is calculated by multiplying the current by the voltage.
The wire is fed to the bottom of the welding torch, heated by the heat source, melted into liquid metal droplets, and deposited on the substrate. The scanning direction is along the x direction. With the movement of the machine tool, the designed track is deposited to complete the WAAM manufacturing process.
Table A2. Process parameters used in the experiments.
Table A2. Process parameters used in the experiments.
ParameterValue (Unit)
Arc power1440 W
Wire diameter1.2 mm
Wire feed speed0.5–0.83 cm/s
Travel speed0.167 cm/s

References

  1. Hu, Z.; Hua, L.; Qin, X.; Ni, M.; Ji, F.; Wu, M. Molten pool behaviors and forming appearance of robotic GMAW on complex surface with various welding positions. J. Manuf. Process. 2021, 64, 1359–1376. [Google Scholar] [CrossRef]
  2. Jia, C.; Liu, W.; Chen, M.; Guo, M.; Wu, S.; Wu, C. Investigation on arc plasma, droplet, and molten pool behaviours in compulsively constricted WAAM. Addit. Manuf. 2020, 34, 101235. [Google Scholar] [CrossRef]
  3. DebRoy, T.; Wei, H.; Zuback, J.; Mukherjee, T.; Elmer, J.; Milewski, J.; Beese, A.M.; Wilson-Heid, A.d.; De, A.; Zhang, W. Additive manufacturing of metallic components–process, structure and properties. Prog. Mater. Sci. 2018, 92, 112–224. [Google Scholar] [CrossRef]
  4. Geng, H.; Li, J.; Xiong, J.; Lin, X.; Huang, D.; Zhang, F. Formation and improvement of surface waviness for additive manufacturing 5A06 aluminium alloy component with GTAW system. Rapid Prototyp. J. 2018, 24, 342–350. [Google Scholar] [CrossRef]
  5. Aldalur, E.; Suárez, A.; Veiga, F. Metal transfer modes for Wire Arc Additive Manufacturing Al-Mg alloys: Influence of heat input in microstructure and porosity. J. Mater. Process. Technol. 2021, 297, 117271. [Google Scholar] [CrossRef]
  6. Xu, T.; Zhang, M.; Wang, J.; Lu, T.; Ma, S.; Liu, C. Research on high efficiency deposition method of titanium alloy based on double-hot-wire arc additive manufacturing and heat treatment. J. Manuf. Process. 2022, 79, 60–69. [Google Scholar] [CrossRef]
  7. Torkamany, M.J.; Kaplan, A.F.; Ghaini, F.M.; Vänskä, M.; Salminen, A.; Fahlström, K.; Hedegård, J. Wire deposition by a laser-induced boiling front. Opt. Laser Technol. 2015, 69, 104–112. [Google Scholar] [CrossRef]
  8. Abioye, T.; Folkes, J.; Clare, A. A parametric study of Inconel 625 wire laser deposition. J. Mater. Process. Technol. 2013, 213, 2145–2151. [Google Scholar] [CrossRef] [Green Version]
  9. Zhao, J.; Zhang, B.; Li, X.; Li, R. Effects of metal-vapor jet force on the physical behavior of melting wire transfer in electron beam additive manufacturing. J. Mater. Process. Technol. 2015, 220, 243–250. [Google Scholar] [CrossRef]
  10. Luo, Y.; Li, J.; Xu, J.; Zhu, L.; Han, J.; Zhang, C. Influence of pulsed arc on the metal droplet deposited by projected transfer mode in wire-arc additive manufacturing. J. Mater. Process. Technol. 2018, 259, 353–360. [Google Scholar] [CrossRef]
  11. Scotti, A.; Ponomarev, V.; Lucas, W. Interchangeable metal transfer phenomenon in GMA welding: Features, mechanisms, classification. J. Mater. Process. Technol. 2014, 214, 2488–2496. [Google Scholar] [CrossRef] [Green Version]
  12. Wu, B.; Ding, D.; Pan, Z.; Cuiuri, D.; Li, H.; Han, J.; Fei, Z. Effects of heat accumulation on the arc characteristics and metal transfer behavior in Wire Arc Additive Manufacturing of Ti6Al4V. J. Mater. Process. Technol. 2017, 250, 304–312. [Google Scholar] [CrossRef]
  13. Zhao, S.; Yang, L.; Huang, Y.; Xu, S. A novel method to fabricate Ni/WC composite coatings by laser wire deposition: Processing characteristics, microstructural evolution and mechanical properties under different wire transfer modes. Addit. Manuf. 2021, 38, 101738. [Google Scholar] [CrossRef]
  14. Tang, Q.; Pang, S.; Chen, B.; Suo, H.; Zhou, J. A three dimensional transient model for heat transfer and fluid flow of weld pool during electron beam freeform fabrication of Ti-6-Al-4-V alloy. Int. J. Heat Mass Transf. 2014, 78, 203–215. [Google Scholar] [CrossRef]
  15. Hu, R.; Chen, X.; Yang, G.; Gong, S.; Pang, S. Metal transfer in wire feeding-based electron beam 3D printing: Modes, dynamics, and transition criterion. Int. J. Heat Mass Transf. 2018, 126, 877–887. [Google Scholar] [CrossRef]
  16. Silwal, B.; Santangelo, M. Effect of vibration and hot-wire gas tungsten arc (GTA) on the geometric shape. J. Mater. Process. Technol. 2018, 251, 138–145. [Google Scholar] [CrossRef]
  17. Chen, X.; Wang, C.; Ding, J.; Bridgeman, P.; Williams, S. A three-dimensional wire-feeding model for heat and metal transfer, fluid flow, and bead shape in wire plasma arc additive manufacturing. J. Manuf. Process. 2022, 83, 300–312. [Google Scholar] [CrossRef]
  18. Ríos, S.; Colegrove, P.A.; Williams, S.W. Metal transfer modes in plasma Wire+ Arc additive manufacture. J. Mater. Process. Technol. 2019, 264, 45–54. [Google Scholar] [CrossRef]
  19. Jing, C.; Chen, Z.; Liu, B.; Xu, T.; Wang, J.; Lu, T.; Lu, J.; Guo, Y.; Liu, C. Improving mechanical strength and isotropy for wire-arc additive manufactured 304L stainless steels via controlling arc heat input. Mater. Sci. Eng. A 2022, 845, 143223. [Google Scholar] [CrossRef]
  20. Peng, W.; Jiguo, S.; Shiqing, Z.; Gang, W. Control of wire transfer behaviors in hot wire laser welding. Int. J. Adv. Manuf. Technol. 2016, 83, 2091–2100. [Google Scholar] [CrossRef]
  21. Nemchinsky, V.A. Size and shape of the liquid droplet at the molten tip of an arc electrode. J. Phys. D Appl. Phys. 1994, 27, 1433. [Google Scholar] [CrossRef]
  22. Cho, M.H.; Farson, D.F. Understanding bead hump formation in gas metal arc welding using a numerical simulation. Metall. Mater. Trans. B 2007, 38, 305–319. [Google Scholar] [CrossRef]
  23. Abe, T.; Kaneko, J.; Sasahara, H. Thermal sensing and heat input control for thin-walled structure building based on numerical simulation for wire and arc additive manufacturing. Addit. Manuf. 2020, 35, 101357. [Google Scholar] [CrossRef]
  24. Ou, W.; Mukherjee, T.; Knapp, G.; Wei, Y.; DebRoy, T. Fusion zone geometries, cooling rates and solidification parameters during wire arc additive manufacturing. Int. J. Heat Mass Transf. 2018, 127, 1084–1094. [Google Scholar] [CrossRef]
  25. Cho, M.; Farson, D. Simulation study of a hybrid process for the prevention of weld bead hump formation. Weld. J. N. Y. 2007, 86, 253. [Google Scholar]
  26. Pardo, E.; Weckman, D. Prediction of weld pool and reinforcement dimensions of GMA welds using a finite-element model. Metall. Mater. Trans. B 1989, 20, 937–947. [Google Scholar] [CrossRef]
  27. Ou, W.; Knapp, G.; Mukherjee, T.; Wei, Y.; DebRoy, T. An improved heat transfer and fluid flow model of wire-arc additive manufacturing. Int. J. Heat Mass Transf. 2021, 167, 120835. [Google Scholar] [CrossRef]
  28. Ke, W.; Oliveira, J.; Cong, B.; Ao, S.; Qi, Z.; Peng, B.; Zeng, Z. Multi-layer deposition mechanism in ultra high-frequency pulsed wire arc additive manufacturing (WAAM) of NiTi shape memory alloys. Addit. Manuf. 2022, 50, 102513. [Google Scholar] [CrossRef]
  29. Cho, J.H.; Farson, D.F.; Milewski, J.O.; Hollis, K.J. Weld pool flows during initial stages of keyhole formation in laser welding. J. Phys. D Appl. Phys. 2009, 42, 175502. [Google Scholar] [CrossRef]
  30. Issa, R.; Ahmadi-Befrui, B.; Beshay, K.; Gosman, A. Solution of the implicitly discretised reacting flow equations by operator-splitting. J. Comput. Phys. 1991, 93, 388–410. [Google Scholar] [CrossRef]
  31. Van Leer, B. Towards the ultimate conservative difference scheme. IV. A new approach to numerical convection. J. Comput. Phys. 1977, 23, 276–299. [Google Scholar] [CrossRef]
  32. Saad, Y. Iterative Methods for Sparse Linear Systems; SIAM: Philadelphia, PA, USA, 2003. [Google Scholar]
  33. Ushio, M.; Wu, C. Mathematical modeling of three-dimensional heat and fluid flow in a moving gas metal arc weld pool. Metall. Mater. Trans. B 1997, 28, 509–516. [Google Scholar] [CrossRef]
  34. Rao, Z.; Hu, J.; Liao, S.; Tsai, H.L. Modeling of the transport phenomena in GMAW using argon–helium mixtures. Part I–The arc. Int. J. Heat Mass Transf. 2010, 53, 5707–5721. [Google Scholar] [CrossRef]
  35. Zhang, X.; Basaran, O.A. An experimental study of dynamics of drop formation. Phys. Fluids 1995, 7, 1184–1203. [Google Scholar] [CrossRef]
  36. Wang, C.; Suder, W.; Ding, J.; Williams, S. The effect of wire size on high deposition rate wire and plasma arc additive manufacture of Ti-6Al-4V. J. Mater. Process. Technol. 2021, 288, 116842. [Google Scholar] [CrossRef]
  37. Hodes, M.; Kirk, T.; Karamanis, G.; Lam, L.; MacLachlan, S.; Papageorgiou, D. Conformal map and asymptotic solutions for apparent slip lengths in the presence of thermocapillary stress. In Proceedings of the 23rd National Heat and Mass Transfer Conference and 1st International ISHMT-ASTFE Heat and Mass Transfer Conference, Thiruvananthapuram, India, 17–20 December 2015; pp. 17–20. [Google Scholar]
  38. Ogino, Y.; Asai, S.; Hirata, Y. Numerical simulation of WAAM process by a GMAW weld pool model. Weld. World 2018, 62, 393–401. [Google Scholar] [CrossRef]
Figure 1. Metal transfer modes in the WAAM process: (a) droplet transfer mode, (b) liquid bridge mode, and (c) wire stubbing transfer mode.
Figure 1. Metal transfer modes in the WAAM process: (a) droplet transfer mode, (b) liquid bridge mode, and (c) wire stubbing transfer mode.
Materials 16 02911 g001
Figure 2. Schematic diagram of the liquid bridge metal transfer mode.
Figure 2. Schematic diagram of the liquid bridge metal transfer mode.
Materials 16 02911 g002
Figure 3. Schematic diagram of the liquid droplet metal transfer mode.
Figure 3. Schematic diagram of the liquid droplet metal transfer mode.
Materials 16 02911 g003
Figure 4. Schematic diagram of interaction area between the heat source and wire.
Figure 4. Schematic diagram of interaction area between the heat source and wire.
Materials 16 02911 g004
Figure 5. High-fidelity numerical model used for calibration of the prediction model.
Figure 5. High-fidelity numerical model used for calibration of the prediction model.
Materials 16 02911 g005
Figure 6. Comparison of the predicted profile (dashed line) and the numerical results for (a) wire feed speed of 2 cm/s, (b) wire feed speed of 1.34 cm/s.
Figure 6. Comparison of the predicted profile (dashed line) and the numerical results for (a) wire feed speed of 2 cm/s, (b) wire feed speed of 1.34 cm/s.
Materials 16 02911 g006
Figure 7. Verification of critical condition for the transition from liquid bridge to droplet mode: (a) initial distance H 0 = 0.14 cm, (b) initial distance H 0 = 0.12 cm.
Figure 7. Verification of critical condition for the transition from liquid bridge to droplet mode: (a) initial distance H 0 = 0.14 cm, (b) initial distance H 0 = 0.12 cm.
Materials 16 02911 g007
Figure 8. Droplet profile predictions. The pictures in black and white are from experiments [35], and the red dashed lines are predicted results: (a) t = 2092 ms, (b) t = 2990 ms, (c) t = 3075 ms.
Figure 8. Droplet profile predictions. The pictures in black and white are from experiments [35], and the red dashed lines are predicted results: (a) t = 2092 ms, (b) t = 2990 ms, (c) t = 3075 ms.
Materials 16 02911 g008
Figure 9. WAAM experiment for dripping frequency under different wire feed speeds: (a) wire feed speed of 0.5 cm/s, (b) wire feed speed of 0.67 cm/s, (c) wire feed speed of 0.83 cm/s.
Figure 9. WAAM experiment for dripping frequency under different wire feed speeds: (a) wire feed speed of 0.5 cm/s, (b) wire feed speed of 0.67 cm/s, (c) wire feed speed of 0.83 cm/s.
Materials 16 02911 g009
Figure 10. Process parameter window considering liquid bridge–droplet transition.
Figure 10. Process parameter window considering liquid bridge–droplet transition.
Materials 16 02911 g010
Figure 11. Liquid bridge and droplet modes predicted by the proposed model and observed from experiments [36].
Figure 11. Liquid bridge and droplet modes predicted by the proposed model and observed from experiments [36].
Materials 16 02911 g011
Figure 12. Process parameter window considering liquid bridge–wire stubbing transition.
Figure 12. Process parameter window considering liquid bridge–wire stubbing transition.
Materials 16 02911 g012
Table 1. Material properties of the mild steel.
Table 1. Material properties of the mild steel.
Material PropertiesValueUnits
Density, ρ 7200kg m 3
Ambient pressure, p 0 0.1MPa
Solidus temperature, T s 1750K
Liquidus temperature, T L 1800K
Vaporization temperature, T v 3143K
Viscosity, μ 6 × 10 3 Pa · s
Thermal conductivity, k27 (293 K)W m 1 K 1
Specific heat capacity, c p 710 (293 K)J kg 1 K 1
Latent heat of fusion, L m 2.47 × 10 5 J kg 1
Latent heat of vaporization, L v 6.08 × 10 6 J kg 1
Molar mass, M5.2 × 10 3 kg mol 1
Surface tension coefficient, σ 0 1.2N m 1
Surface tension gradient, d σ d T −1 × 10 4 N m 1 K 1
Convection coefficient, h c 10W m 2 K 1
Table 2. Temperature dependent properties of the mild steel.
Table 2. Temperature dependent properties of the mild steel.
Temperature (K)29380013001800
Thermal conductivity
(W m 1 K 1 )27.035.443.752.0
Specific heat capacity
(J kg 1 K 1 )710.0740.2770.1800.0
Table 3. Numerical cases for calibrating the bottom radius.
Table 3. Numerical cases for calibrating the bottom radius.
CaseWire Feed Speed (cm/s)Travel Speed (cm/s) ξ Bottom Radius R 0 (cm)Standard Deviation (cm)
10.500.750.670.0850.0195
20.500.501.000.1100.0241
30.750.501.500.1160.0231
41.801.001.800.1370.0245
52.201.002.200.1430.0160
62.000.752.670.1490.0191
74.001.253.000.1700.0125
82.000.54.000.1770.0137
92.250.54.500.2060.0191
Table 4. Comparison of droplet number and dripping frequency between experimental data and theoretical results.
Table 4. Comparison of droplet number and dripping frequency between experimental data and theoretical results.
Case 1Case 2Case 3
Wire feed speed0.50 cm/s0.67 cm/s0.83 cm/s
Experimental number of droplets71013
Theoretical number of droplets (taking floor)6810
Experimental dripping frequency0.117 Hz0.167 Hz0.217 Hz
Theoretical dripping frequency0.105 Hz0.141 Hz0.176 Hz
Relative error (for frequency)10.2%15.6%18.9%
Table 5. Numerical cases to validate the prediction model for liquid bridge–droplet transition.
Table 5. Numerical cases to validate the prediction model for liquid bridge–droplet transition.
CaseWire Feed Speed (cm/s)Travel Speed (cm/s) H 0 (cm)Power (W)
12.001.250.101440
21.000.500.101440
31.500.500.101440
42.001.250.111440
51.000.750.111440
64.001.500.121440
71.500.500.121440
84.001.000.121440
92.200.700.131440
103.501.000.141440
114.001.000.141440
Table 6. Numerical cases to validate the prediction model for the liquid bridge–wire stubbing mode transition.
Table 6. Numerical cases to validate the prediction model for the liquid bridge–wire stubbing mode transition.
CaseWire Feed Speed (cm/s)Travel Speed (cm/s)Power (W)
12.001.001440
24.001.001440
35.001.001440
46.001.001440
52.501.001100
63.501.001100
74.501.001100
82.001.00800
93.001.00800
104.001.00800
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Ouyang, J.; Li, M.; Lian, Y.; Peng, S.; Liu, C. A Fast Prediction Model for Liquid Metal Transfer Modes during the Wire Arc Additive Manufacturing Process. Materials 2023, 16, 2911. https://doi.org/10.3390/ma16072911

AMA Style

Ouyang J, Li M, Lian Y, Peng S, Liu C. A Fast Prediction Model for Liquid Metal Transfer Modes during the Wire Arc Additive Manufacturing Process. Materials. 2023; 16(7):2911. https://doi.org/10.3390/ma16072911

Chicago/Turabian Style

Ouyang, Jiaqi, Mingjian Li, Yanping Lian, Siyi Peng, and Changmeng Liu. 2023. "A Fast Prediction Model for Liquid Metal Transfer Modes during the Wire Arc Additive Manufacturing Process" Materials 16, no. 7: 2911. https://doi.org/10.3390/ma16072911

APA Style

Ouyang, J., Li, M., Lian, Y., Peng, S., & Liu, C. (2023). A Fast Prediction Model for Liquid Metal Transfer Modes during the Wire Arc Additive Manufacturing Process. Materials, 16(7), 2911. https://doi.org/10.3390/ma16072911

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop