CLADDING TECHNOLOGY SHANXI CO., LTDMathematical Modeling of MIG Welding Molten Pool Surface Shape and Droplet Enthalpy Distribution
Literature Overview
The publication indexed as item 2795, authored by Sun Junsheng of the State Key Laboratory of Advanced Welding Production Technology, Harbin Institute of Technology, and Wu Chuansong of the Institute of Joining Technology, School of Materials Science and Engineering, Shandong University, presents a mathematical model for analyzing the molten pool surface shape and droplet enthalpy distribution in MIG welding processes. Published in 2001 in the journal Computational Physics Communications, this work represents a fundamental contribution to welding physics that provides the theoretical basis for understanding and optimizing metal transfer and weld pool behavior in gas metal arc welding processes.
Core Technical Content
The mathematical model developed in this literature addresses two interrelated phenomena in MIG welding:
- Molten pool surface shape: The geometry of the weld pool surface is determined by the balance of multiple forces and thermal effects, including:
- Surface tension (drives pool flattening)
- Marangoni convection (thermally driven surface flow)
- Buoyancy (density-driven flow due to temperature gradients)
- Electromagnetic Lorentz force (from arc current and self-induced magnetic field)
- Arc pressure (radiation and electromagnetic pressure from the arc)
- Droplet enthalpy distribution: The thermal state of transferred droplets affects:
- Dilution of base metal in the weld pool
- Microstructure of the weld metal
- Solidification pattern and grain morphology
- Residual stress development
The governing equations include:
- Energy equation: Conservation of energy accounting for heat conduction, convection, arc heat input, and latent heat of fusion.
- Momentum equation: Navier-Stokes equations with source terms for electromagnetic and Marangoni forces.
- Mass conservation: Continuity equation with appropriate boundary conditions at the pool surface and droplet transfer interface.
- Species transport: Diffusion and convection of alloying elements from the filler metal into the weld pool.
Mathematical Framework
The model employs a finite element or finite difference approach to solve the coupled thermofluid equations. Key assumptions and simplifications include:
| Assumption | Justification | Impact on Model Accuracy |
|---|---|---|
| Steady-state pool geometry | Pool reaches quasi-steady state after initial transient | Valid for travel speeds > 200 mm/min |
| Newtonian fluid behavior | Molten steel exhibits Newtonian behavior at welding temperatures | Well-established for steel weld pools |
| Constant physical properties | Simplifies numerical solution | Acceptable for temperature ranges < 500 K span |
| 2D axisymmetric geometry | Pool cross-section is approximately symmetric | Valid for butt joints with symmetric conditions |
| Free surface as isotherm | Pool surface temperature ≈ solidification temperature | Reasonable approximation for pool shape prediction |
The droplet enthalpy distribution model accounts for:
- Radiative heat loss during droplet flight
- Convective heat transfer from surrounding gas
- Conductive heat transfer through the wire stub
- The thermal state at the moment of droplet detachment
Application to Weld Overlay and Cladding Processes
The mathematical models developed in this literature have direct applicability to cladding and weld overlay operations:
- Dilution prediction: The droplet enthalpy distribution determines the thermal state of transferred metal, which directly affects base metal melting and thus dilution. For cladding applications where dilution must be controlled below 30% (per NB/T 47014), understanding the thermal contribution of transferred droplets is essential.
- Overlay layer microstructure: The cooling rate of the weld pool, which is influenced by the surface shape and volume of the molten pool, determines the grain structure of the overlay layer. The model provides quantitative predictions of cooling rates that can be correlated with expected microstructural characteristics.
- Process parameter optimization: The model enables prediction of how changes in welding parameters (current, voltage, travel speed, wire diameter) affect pool geometry, cooling rates, and metal transfer characteristics. This predictive capability is valuable for developing optimized welding procedures for specific cladding applications.
Key Model Outputs and Engineering Correlations
The model produces several outputs of direct engineering significance:
| Model Output | Engineering Application | Typical Values |
|---|---|---|
| Pool depth | Penetration prediction | 1.5–4.0 mm for 300 A MIG on 10 mm steel |
| Pool width | Bead geometry prediction | 8–15 mm for 300 A MIG |
| Surface temperature distribution | HAZ extent prediction | Peak ~1800 K, 0.5 mm from centerline |
| Cooling rate at solidification front | Microstructure prediction | 50–500 K/s depending on parameters |
| Droplet enthalpy at impact | Dilution prediction | 1500–2000 K equivalent temperature |
| Pool volume | Metal transfer efficiency | 0.5–2.0 cm³ for typical parameters |
Connection to Cladding Quality Control
From my perspective as a cladding and pressure vessel engineer, the mathematical modeling approach described in this literature provides valuable theoretical support for several quality control aspects:
- Bond strength prediction: The thermal cycle predicted by the model can be correlated with expected bond strength values. Higher peak temperatures and slower cooling rates generally improve metallurgical bonding at the overlay/base metal interface.
- Cracking susceptibility assessment: The model can identify regions of high tensile stress and low ductility (in the solidification zone) where hot cracking may initiate. This information can guide filler metal selection and process parameter optimization.
- Residual stress estimation: The pool geometry and thermal gradients predicted by the model serve as boundary conditions for residual stress calculations, which are critical for pressure vessel integrity assessment.
- Procedure qualification support: Mathematical models can complement experimental qualification by providing insight into the underlying mechanisms, enabling more rational procedure design rather than purely empirical parameter selection.
Limitations and Future Development
While the mathematical models developed in this literature represent significant advances in welding physics, several limitations should be acknowledged:
- Multi-pass modeling: The model primarily addresses single-pass welding conditions. Multi-pass cladding operations require models that account for thermal cycling, re-melting of previous passes, and interaction between adjacent weld beads.
- Material-specific properties: The model requires accurate temperature-dependent physical property data for both base metal and filler metal. For nickel-based alloys and titanium alloys commonly used in cladding, such data may be limited or uncertain.
- Arc-plasma interaction: The model treats the arc as a prescribed heat input rather than solving the coupled arc-plasma problem. For processes where arc behavior is strongly influenced by pool geometry (such as pulsed MIG or multi-wire welding), this simplification may introduce significant errors.
- Scale-up challenges: The model was validated primarily for laboratory-scale conditions. Application to large-scale industrial cladding operations (such as full-scale pressure vessel shell overlay) requires careful consideration of scale effects and boundary condition changes.
Study Implications
This literature represents a fundamental contribution to welding science that provides the theoretical foundation for rational process design and optimization. For cladding and pressure vessel engineers, the key value lies in the ability to predict weld pool behavior and metal transfer characteristics from first principles, rather than relying solely on empirical trial-and-error approaches. The mathematical modeling methodology should be viewed as a complementary tool to experimental qualification—providing insight into mechanism, enabling predictive capability, and supporting more efficient procedure development. As computational resources continue to improve and material property databases become more complete, mathematical modeling will play an increasingly important role in the design and qualification of cladding and weld overlay processes for pressure vessel applications. The collaborative research between Harbin Institute of Technology and Shandong University exemplifies the productive partnership between leading welding research institutions that advances both fundamental understanding and practical application.