Mathematical Modeling of MIG Weld Pool Surface Shape and Droplet Enthalpy Distribution

Literature Overview

This study by Sun Junsheng (Harbin Institute of Technology, State Key Laboratory of Modern Welding Production Technology) and Wu Chuansong (Shandong University, Institute of Joining Technology) was published in 2001 and represents an early foundational work in computational welding mechanics. The research develops a mathematical model for the surface shape of the MIG (Gas Metal Arc Welding) melt pool and the enthalpy distribution within the transferred droplets. The work was supported by the Key Laboratory Visiting Scholar Fund for Higher Education Institutions, reflecting its significance in the Chinese welding research community at the time.

Core Technical Content

The fundamental challenge addressed in this paper is the accurate description of the melt pool geometry and the thermal energy carried by each transferred droplet during GMAW/MIG welding. In conventional welding models, the heat input is often simplified as a Gaussian or double-elliptical surface heat flux, which fails to capture the dynamic nature of the droplet transfer process. This study attempts to bridge that gap by coupling the droplet enthalpy model with the melt pool surface shape prediction.

The mathematical framework considers several key physical phenomena:

• Surface tension effects on the melt pool free surface
• Marangoni convection driven by surface tension gradients
• Radiative and convective heat loss from the pool surface
• The thermal energy content of each individual droplet at the moment of detachment

Key Equations and Assumptions

The droplet enthalpy is calculated by integrating the temperature-dependent specific heat over the droplet volume, accounting for the temperature gradient between the wire tip and the droplet center. The melt pool surface shape is derived from a free-surface boundary condition where the pressure, surface tension, and gravitational forces are in equilibrium.

Interpretation of Technical Points

The enthalpy distribution model is particularly relevant to cladding applications because it determines the effective heat input delivered to the base metal by each droplet. In weld-overlay processes, the dilution ratio between the overlay material and the base metal is one of the most critical quality parameters. For example, in 308L stainless steel cladding on carbon steel (per ASTM A263/A264), the dilution must typically remain below 30% to ensure adequate corrosion resistance of the overlay layer.

The surface shape model provides insight into the pool geometry, which governs the solidification pattern and microstructure of the overlay. A deeper, narrower pool (characteristic of higher current with active gas) tends to produce more dilution, while a shallower, wider pool (inert gas, lower current) reduces dilution but may compromise mechanical bonding.

Connection with Engineering Practice

In bimetal pressure vessel fabrication, particularly for hydrogenation reactors lined with 316L stainless steel or Inconel 625, understanding the heat input distribution is essential for:

1. Predicting the depth of base metal penetration and ensuring adequate overlay thickness
2. Controlling the heat-affected zone (HAZ) to prevent sensitization or martensitic transformation
3. Optimizing multi-pass cladding sequences to achieve uniform overlay composition

Practical Implications for Cladding Process Design

When applying this model to overlay welding, the following considerations arise:

• The droplet enthalpy directly correlates with the local melting rate of the base metal, which determines dilution in each pass
• The melt pool surface shape affects the geometric profile of each overlay pass, influencing the final surface finish and thickness uniformity
• For electroslag welding (ESW) overlay, the analogous principle applies where the slag pool geometry and heat distribution determine the cladding quality

A notable engineering case involves the overlay of Hastelloy C276 on low-alloy steel for sulfuric acid service. The dilution must be kept below 20% to maintain the alloy's resistance to reducing acids. Using the enthalpy model, engineers can predict how changes in wire feed speed, travel speed, and arc current affect the effective heat input per unit length and thus the dilution ratio.

Key Questions and Reflections

The 2001 publication date raises an important question: how well does this model account for the complex physics of modern multi-wire GMAW overlay systems or hot-wire TIG cladding? The droplet enthalpy concept remains valid, but the multi-wire interaction and the significantly different heat input profiles of hot-wire processes would require substantial model extensions.

Another reflection is the practical applicability of the melt pool surface shape prediction. In real cladding operations, the pool is constantly moving, and the dynamic surface shape differs significantly from steady-state predictions. However, for process parameter selection and initial qualification, the model provides a valuable starting point that can be refined through experimental validation.

The study also raises the question of whether the enthalpy distribution should be coupled with the dilution prediction model for more accurate cladding qualification. In practice, engineers often use empirical dilution curves derived from trial welds, but a physics-based approach as presented here could reduce the number of qualification trials required.

Study Insights and Implications

This foundational work reminds us that accurate thermal modeling of the welding process remains the cornerstone of overlay process design. While computational resources and software have advanced dramatically since 2001, the fundamental physics described in this paper—droplet enthalpy, melt pool surface tension, and Marangoni convection—remain unchanged. Modern finite element welding models for cladding still rely on similar heat source formulations, often enhanced with experimental validation data.

For engineers working on bimetal pressure vessel cladding, the key takeaway is that the effective heat input is not simply a function of arc power divided by travel speed. The droplet transfer dynamics add a layer of complexity that must be considered, especially in low-dilution overlay applications where even small changes in heat input can significantly affect the overlay composition and, consequently, the corrosion performance of the finished component.