Zhuojin Pipe Fitting Co., LtdMathematical Model of Penetration Weld Pool in TIG Welding
Literature Overview
The paper by Cao Zhenning, Wu Chuansong, and Wu Lin, published in Transactions of the China Welding Institution (1996, Vol. 17, No. 1, pp. 62–70), presents a comprehensive mathematical model for the weld pool under penetration conditions in TIG welding. Funded by the National Natural Science Foundation of China, this work derives the deformation equations for both the upper and lower surfaces of a through-penetration weld pool using fluid mechanics theory and variational principles. The model employs non-orthogonal body-fitted curvilinear coordinate systems to handle the complex curved boundaries of the weld pool surfaces, and it couples the flow field and thermal field calculations with surface deformation analysis.
Theoretical Framework
The mathematical model is built upon several fundamental theoretical foundations:
- Fluid mechanics theory: The weld pool is treated as a viscous incompressible fluid governed by the Navier-Stokes equations, with appropriate boundary conditions at the free surfaces.
- Variational principles: The surface deformation equations are derived using variational methods, which provide a rigorous mathematical framework for determining the equilibrium shape of the weld pool surfaces under the action of surface tension, buoyancy, electromagnetic forces, and viscous stresses.
- Non-orthogonal curvilinear coordinates: A key innovation is the adoption of body-fitted non-orthogonal curvilinear coordinate systems that conform to the deformed weld pool surfaces. This approach allows accurate numerical treatment of the complex boundary geometry that arises when the weld pool penetrates through the plate.
| Model Component | Mathematical Basis | Key Feature |
|---|---|---|
| Flow field | Navier-Stokes equations | Viscous fluid dynamics |
| Thermal field | Heat conduction equation | Coupled with flow |
| Upper surface deformation | Variational principle | Surface tension, pressure balance |
| Lower surface deformation | Variational principle | Capillary force, viscous stress |
| Coordinate system | Non-orthogonal curvilinear | Body-fitted to surfaces |
Coupled Analysis Approach
The distinguishing feature of this model is the simultaneous computation of:
- The three-dimensional flow field within the weld pool
- The three-dimensional thermal field within the weld pool
- The deformation of the upper weld pool surface
- The deformation of the lower weld pool surface (penetration side)
This coupled approach is essential because the surface deformation directly affects the boundary conditions for the flow and thermal fields, while the flow and thermal fields determine the stresses acting on the surfaces. An iterative solution procedure is required to achieve self-consistency between the surface shapes and the internal field solutions.
The forces acting on the weld pool surfaces include:
- Surface tension: Acts as a restoring force that tends to minimize surface area
- Buoyancy force: Driven by density differences due to temperature gradients
- Electromagnetic force: Lorentz force from the interaction of current density and magnetic field
- Atmospheric pressure: Acts on the exposed upper surface
- Viscous stress: Internal friction at the free surface
Numerical Implementation and Validation
The numerical solution employs the finite volume method on a non-orthogonal grid that conforms to the weld pool boundaries. The non-orthogonal coordinate transformation introduces metric terms into the governing equations, which must be carefully discretized to maintain numerical accuracy. The grid is adapted iteratively as the surface shapes evolve during the solution process.
The model was validated against welding process experiments, and the calculated weld bead profiles showed good agreement with experimental measurements. This validation confirms that the mathematical formulation and numerical implementation are physically sound and capable of predicting practical welding outcomes.
Engineering Significance
This type of mathematical modeling has several important applications in welding engineering:
- Process optimization: By systematically varying welding parameters in the model, optimal conditions for achieving desired weld geometry can be identified before physical experimentation.
- Quality prediction: The model can predict weld pool dimensions, penetration depth, and surface profile for given parameters, enabling quality assessment without destructive testing.
- Troubleshooting: When welding defects occur in practice, the model can help identify whether the root cause lies in parameter selection, material properties, or geometric configuration.
- Equipment design: Understanding the weld pool dynamics aids in designing torch holders, backing fixtures, and positioning systems that accommodate the expected weld pool behavior.
| Application Area | Benefit | Limitation |
|---|---|---|
| Parameter optimization | Reduced trial-and-error | Requires accurate boundary conditions |
| Defect prediction | Proactive quality control | May not capture all defect mechanisms |
| Equipment design | Better fixture design | Simplified boundary conditions |
| Process development | Faster qualification | Validation still required |
Study Insights and Reflections
This 1996 paper represents a significant milestone in computational welding science, particularly for the Chinese welding research community. The successful implementation of non-orthogonal body-fitted coordinates for the complex weld pool geometry demonstrates sophisticated numerical methods expertise. The coupling of surface deformation with internal field calculations represents a level of physical fidelity that was uncommon in welding models at that time.
The model's strength lies in its rigorous mathematical foundation and its ability to capture the essential physics of weld pool behavior under penetration conditions. However, as with all computational models, its predictive accuracy depends on the quality of input data including surface tension coefficients, electrical conductivity, and thermophysical properties as functions of temperature. The model also assumes axisymmetric conditions, which limits its applicability to non-axisymmetric geometries or complex joint configurations.
The work by this research group at Shandong University of Technology and Harbin Institute of Technology laid important groundwork for subsequent advances in computational welding science in China, including more sophisticated models incorporating multiphase flow, solidification microstructure, and residual stress analysis.
Concluding Remarks
This paper presents a rigorous and physically grounded mathematical model for penetration weld pools in TIG welding that successfully couples flow field, thermal field, and surface deformation calculations using non-orthogonal body-fitted coordinates. The validated agreement between predicted and experimental weld profiles demonstrates the model's practical utility for welding process optimization and quality prediction, establishing a computational methodology that continues to influence modern welding simulation approaches.