CLADDING TECHNOLOGY SHANXI CO., LTDNumerical Simulation of Molten Pool Surface Deformation in TIG Welding
Literature Overview
This paper, authored by Zhang Jianbao and Wang Hong from the School of Materials Science and Engineering at Taiyuan University of Science and Technology, was published in Hot Working Technology in 2016. The work presents a numerical simulation study of the molten pool surface deformation during TIG welding, focusing on the complex interplay between electromagnetic forces, surface tension gradients, and fluid flow dynamics that determine the pool shape and weld bead geometry.
Core Technical Content
The numerical simulation employs computational fluid dynamics (CFD) coupled with electromagnetic field calculations to model the molten pool during TIG welding. The governing equations include the Navier-Stokes equations for fluid flow, the energy equation for heat transfer, and the electromagnetic force equations for Lorentz force and surface tension force calculations.
Governing Equations and Boundary Conditions
| Equation Type | Governing Equation | Boundary Condition |
|---|---|---|
| Mass conservation | ∂ρ/∂t + ∇·(ρv) = 0 | No-flow at free surface |
| Momentum | ρ(∂v/∂t + v·∇v) = -∇p + μ∇²v + F_L + F_ST | Surface tension at free surface |
| Energy | ρc_p(∂T/∂t + v·∇T) = ∇·(k∇T) + Q_arc | Adiabatic at free surface |
| Electromagnetic | F_L = J × B | Current distribution from arc model |
The surface tension force is calculated using the Brinkman penalty method, which models the free surface as a thin layer of viscous material. The surface tension coefficient is temperature-dependent, following the empirical relationship γ = γ_0(1 - T/T_m), where γ_0 is the surface tension at room temperature and T_m is the melting temperature.
Simulation Parameters
| Parameter | Value | Source |
|---|---|---|
| Welding current | 150-250 A | Typical TIG range |
| Travel speed | 3-8 mm/s | Production range |
| Arc radius | 1.5-2.5 mm | Arc model |
| Mesh size | 0.1-0.5 mm | Near pool region |
| Time step | 0.001-0.01 s | Stability requirement |
| Material | 304 stainless steel | Common test material |
Interpretation of Technical Points
The simulation results reveal that the molten pool surface deformation is primarily governed by the competition between electromagnetic depression force and surface tension restoration force. At higher currents, the electromagnetic force dominates, creating a deeper pool depression. However, the surface tension gradient, driven by temperature variations across the pool surface, generates Marangoni convection that significantly influences the pool shape and weld bead profile.
The Marangoni effect is particularly important because the surface tension of most metals decreases with increasing temperature. This creates a surface tension gradient that drives fluid flow from the hot center of the pool toward the cooler edges. This outward flow at the surface is balanced by inward flow at the bottom of the pool, creating a circulation pattern that affects heat distribution and penetration depth.
Effect of Process Parameters on Pool Shape
| Parameter Variation | Pool Depth Change | Pool Width Change | Bead Profile |
|---|---|---|---|
| Increase current | Increases | Increases | Deeper, wider |
| Increase travel speed | Decreases | Decreases | Shallower, narrower |
| Add surfactant | Decreases | Increases | Shallower, wider |
| Increase torch angle | Asymmetric | Asymmetric | Asymmetric profile |
| Increase gap | Increases | Unchanged | Deeper with root reinforcement |
The numerical simulation provides insights that are difficult to obtain experimentally. The internal flow patterns, temperature distributions, and force balances within the pool are invisible to conventional measurement techniques. This information is critical for understanding weld formation mechanisms and predicting weld quality.
Engineering Practice Integration
In cladding and overlay welding operations, the molten pool shape directly determines the dilution level between the overlay material and the base metal. For nickel-based alloy cladding on carbon steel, excessive dilution reduces the corrosion resistance of the overlay layer, while insufficient dilution compromises the metallurgical bond strength. The numerical simulation can predict dilution levels under various process conditions, enabling optimization before actual welding trials.
The simulation approach is particularly valuable for developing welding procedures for critical applications where experimental trials are expensive or time-consuming. For example, when developing a cladding procedure for a large hydrogenation reactor, the cost of each trial weld can be substantial due to the material costs and fabrication time. Numerical simulation can narrow the parameter search space, reducing the number of experimental trials required for procedure qualification.
Application to Cladding Process Development
| Cladding Application | Critical Pool Parameter | Simulation Insight |
|---|---|---|
| Inconel 625 on carbon steel | Dilution ratio | Pool depth determines dilution level |
| 316L on 304 SS | Heat input | Pool temperature affects grain growth |
| Copper on steel | Bond strength | Pool shape affects interface metallurgy |
| Multi-layer cladding | Interpass temperature | Pool cooling rate affects phase transformation |
The simulation can also be used to predict the effects of process disturbances on weld quality. For example, variations in travel speed due to operator inconsistency or mechanical issues can be simulated to determine their impact on pool shape and weld bead geometry. This information can be used to set tolerance limits for process parameter variations in production welding.
Study Insights and Reflections
The numerical simulation of TIG welding molten pools has advanced significantly over the past two decades, but challenges remain in accurately modeling the arc-plasma interaction and the complex physics at the pool surface. The authors' work contributes to the ongoing refinement of these models, particularly in the treatment of surface tension effects and electromagnetic force calculations.
For practical engineering applications, the key value of numerical simulation lies in its ability to provide design guidance and process optimization insights. However, simulation results must always be validated against experimental data before being relied upon for production decisions. The ideal approach combines simulation-based optimization with experimental verification, using simulation to narrow the parameter space and experiments to confirm the optimal settings.
The methodology presented in this paper can be extended to model more complex welding scenarios, including multi-pass welding with residual stress effects, welding of dissimilar materials with different thermal properties, and cladding operations with significant dilution concerns. As computational resources continue to improve, real-time simulation coupled with process monitoring could enable adaptive welding control systems that adjust parameters in response to changing conditions.