Mathematical Model of MIG Welding Molten Pool Surface Shape and Droplet Enthalpy Distribution

Literature Overview

This paper by Sun Junsheng and Wu Chuansong, published in Computational Physics (2001, Vol. 18, Issue 6, pp. 544-548), represents a foundational contribution to the computational modeling of MIG welding physics. The study develops coupled mathematical models for the molten pool surface deformation and the distribution of droplet enthalpy within the weld pool. Authored by researchers from the State Key Laboratory of Modern Welding at Harbin Institute of Technology and Shandong University, this work reflects the depth of Chinese welding research institutions and provides theoretical underpinnings that remain relevant to modern process modeling.

Model Development Framework

The mathematical model addresses two coupled physical phenomena:

The model formulation considers the following physical interactions simultaneously:

• Arc electromagnetic force acting on the molten pool surface
• Surface tension gradients (Marangoni effect) driven by temperature-dependent surface tension
• Buoyancy forces from thermal density variations
• Arc plasma pressure and convective cooling on the pool surface
• Electromagnetic stirring within the pool volume
• Heat input from transferred droplets carrying enthalpy

Numerical Simulation Results

The numerical simulations reveal several important relationships between welding parameters and weld pool behavior:

Effect of welding current on pool surface shape:

• Increasing current increases the electromagnetic force on the pool surface, causing deeper depression at the arc contact point
• Higher currents also increase Marangoni-driven surface flow, spreading the pool laterally
• The combined effect produces a wider, deeper pool with a more pronounced surface depression

Effect of travel speed on enthalpy distribution:

• Higher travel speeds concentrate the enthalpy distribution in the rear portion of the pool
• This asymmetry affects solidification direction and grain structure orientation
• At very high travel speeds, insufficient enthalpy delivery to the leading edge can cause lack of fusion

Effect of electrode diameter on pool geometry:

• Larger electrode diameters produce wider, shallower pools due to increased arc spreading
• Smaller electrode diameters concentrate heat input, producing narrower, deeper pools
• This relationship directly influences weld bead geometry predictions

Engineering Practice Relevance

For welding engineers working on pipe and fitting production, the model provides several practical insights:

1. Weld geometry prediction: Understanding the relationship between arc forces and pool surface shape enables prediction of weld bead profiles without extensive trial welding. This is particularly valuable for establishing welding procedure specifications for new materials or geometries.
2. Porosity prevention: The enthalpy distribution model helps identify conditions where insufficient heat input to the pool front can cause incomplete melting and subsequent porosity formation. This is critical for thin-wall pipe welding where lack of fusion is a common defect.
3. Distortion prediction: The asymmetric enthalpy distribution at high travel speeds correlates with directional thermal distortion. Engineers can use this understanding to design welding sequences that minimize cumulative distortion in pipe assemblies.
4. Parameter optimization: The coupled model allows systematic exploration of parameter combinations to achieve target weld geometries, reducing the need for extensive trial-and-error qualification testing.

Model Limitations and Development Needs

While the model represents significant progress, several limitations should be noted for practical application:

• The model assumes steady-state conditions, which is not applicable to start and end conditions, weld interruptions, or multi-pass welding sequences.
• The material property database must be temperature-dependent and specific to the alloy being welded; generic properties will produce inaccurate results.
• The model does not account for gas shielding dynamics, which can significantly affect pool surface conditions, especially in outdoor or high-airflow environments.
• Solidification modeling is not included, so microstructural predictions require separate analysis.

For modern applications, the model framework should be extended to include:

• Transient thermal analysis for multi-pass welding
• Coupled thermo-mechanical analysis for residual stress prediction
• Phase transformation modeling for microstructural prediction
• Fluid-structure interaction for thin-wall pipe welding where pool oscillation and sagging are concerns

Key Questions and Reflections

The model's treatment of droplet enthalpy as a distributed heat source within the pool is an elegant approach that captures the physical reality of droplet impact and mixing. However, in practice, the droplet transfer mode (short-circuit, globular, spray) varies significantly with welding parameters and has a profound effect on pool dynamics. A complete model should incorporate droplet transfer mode-dependent boundary conditions rather than treating enthalpy distribution as a continuous function.

The Marangoni effect, driven by surface tension gradients, is a dominant mechanism in pool surface flow. In steel welding, the surface tension decreases with temperature (due to sulfur and oxygen surfactants), driving outward flow and producing wide, shallow welds. In aluminum welding, the surface tension increases with temperature (due to surfactant depletion), driving inward flow and producing narrow, deep welds. The model should explicitly account for this material-dependent behavior, as it fundamentally changes the pool dynamics and weld geometry predictions.

Study Insights and Implications

This paper represents an important step in the computational modeling of MIG welding, providing a rigorous physical framework for understanding the coupled phenomena of pool deformation and heat distribution. For practicing welding engineers, the value lies not in the specific numerical results but in the conceptual framework: welding is a multi-physics problem where electromagnetic, thermal, fluid dynamic, and surface tension effects interact in complex ways. Understanding these interactions enables better process design, defect prevention, and quality control. The model also highlights the importance of numerical simulation as a complementary tool to experimental welding, reducing development costs and accelerating process qualification. For the pipe and fitting industry, where welding quality directly affects structural integrity and safety, such computational tools are indispensable for ensuring consistent production quality across diverse geometries and materials.