Friction Stir Welding (FSW) Solid-State Temperature Regulation: Analytical Heat Source vs. ALE Simulation Methodology

1. Technical Definition and Core Principles

Friction Stir Welding (FSW) is a solid-state joining process in which a non-consumable rotating tool—comprising a shoulder and a pin—is plunged into the joint interface between two workpieces. Heat is generated exclusively through frictional contact between the tool and the base material, combined with plastic deformation energy. Unlike fusion welding processes, FSW operates entirely below the melting point of the base material, typically maintaining peak temperatures in the range of 0.8 to 0.9 times the absolute melting temperature (T_melt) of the material being joined.

The solid-state temperature regulation mechanism is the critical factor governing FSW process quality. It encompasses the dynamic thermal field distribution, the rate of heat input, heat dissipation pathways, and the resulting microstructural evolution in the stir zone, thermally affected zone (TAZ), and heat-affected zone (HAZ). Proper temperature control ensures adequate material plasticity for defect-free bonding while avoiding thermal damage, excessive grain growth, or residual stress accumulation.

The analytical heat source model represents a simplified, mathematically tractable approach to describing the thermal distribution in FSW. This model typically employs a double-ellipsoidal or Gaussian-type heat source function that approximates the frictional and deformation heat generation zones. The analytical solution provides rapid parametric evaluation of temperature fields but sacrifices spatial fidelity in regions of complex tool-material interaction.

The Arbitrary Lagrangian-Eulerian (ALE) method, by contrast, is a hybrid numerical formulation that combines the advantages of Lagrangian mesh (material tracking) and Eulerian mesh (fluid-like flow resolution). In FSW simulation, ALE handles the severe plastic deformation and material flow around the rotating pin without mesh distortion, enabling accurate prediction of transient temperature fields, flow patterns, and stress states throughout the welding cycle.

2. Comparative Analysis: Analytical Heat Source vs. ALE Method

2.1 Methodological Framework

The analytical heat source approach models the FSW thermal input as a moving volumetric or surface heat source traveling at the welding speed. The heat generation rate is derived from the friction coefficient, tool geometry, rotational speed, and material properties. The governing heat equation is solved either semi-analytically (Rosenthal-type solutions) or through finite element discretization with prescribed boundary conditions.

The ALE method solves the coupled thermo-mechanical problem by dividing the computational domain into a Lagrangian region (where material moves with the mesh, representing the bulk workpiece) and an Eulerian region (where the mesh is stationary while material flows through it, representing the highly deformed stir zone). The interface between these regions is managed through mesh smoothing, remapping, or rezoning algorithms to prevent element distortion.

2.2 Performance Comparison

Parameter	Analytical Heat Source Model	ALE Method
Computational Cost	Low (seconds to minutes)	High (hours to days)
Thermal Field Accuracy	Adequate for far-field; limited near pin	High fidelity throughout domain
Material Flow Prediction	Not available	Full 3D flow field resolution
Stress/Strain Coupling	Limited or absent	Full thermo-mechanical coupling
Tool Geometry Fidelity	Simplified (idealized shapes)	Full 3D tool geometry
Phase Transformation Modeling	Not feasible	Feasible with constitutive extensions
Process Window Optimization	Excellent for rapid screening	Refined validation of candidate windows
Residual Stress Prediction	Approximate	Accurate with proper boundary conditions

2.3 Temperature Distribution Characteristics

Both methods confirm that the peak temperature in FSW occurs at the interface between the pin and the workpiece, typically on the trailing side of the pin due to the asymmetry in material flow. The shoulder contributes to heat input over a broader area, while the pin concentrates heat generation in a smaller volume. The temperature gradient from the stir zone to the un-deformed base material is typically steep, with a thermal gradient of 100–500 °C/mm in the immediate vicinity of the tool.

The analytical model tends to overestimate peak temperatures by 50–150 °C compared to ALE results when calibrated to the same friction coefficient, primarily because it cannot account for the convective heat transport caused by material flow (advection of hot material away from the pin). The ALE method captures this convective cooling effect, resulting in more realistic peak temperature predictions and more accurate representation of the temperature field asymmetry between the leading and trailing sides.

3. Solid-State Temperature Regulation Mechanism

3.1 Heat Generation Sources

In FSW, heat is generated through two primary mechanisms:

• Interfacial Friction Heat: Generated at the contact surfaces between the tool shoulder/pin and the workpiece material. The frictional heat flux is proportional to the friction coefficient (μ), the normal contact pressure (p), and the sliding velocity (v): q = μ · p · v. The friction coefficient is temperature-dependent, typically ranging from 0.2 to 0.4 for aluminum alloys and 0.3 to 0.5 for steels.
• Plastic Deformation Heat: Generated in the bulk material undergoing severe plastic strain around the pin. The deformation heat rate is given by: Q_def = σ · ε̇, where σ is the flow stress and ε̇ is the strain rate. This component is particularly significant in the stir zone where strain rates can exceed 10–100 s^-1.

3.2 Heat Dissipation Pathways

Heat dissipation in FSW occurs through:

• Conduction into the workpiece: The dominant heat sink, governed by the thermal diffusivity of the base material. Aluminum alloys (α ≈ 8.4 × 10^-5 m²/s) dissipate heat more readily than steels (α ≈ 1.2 × 10^-5 m²/s), requiring higher heat inputs for equivalent plasticity.
• Convective transport by material flow: Hot material in the stir zone is advected away from the pin, carrying thermal energy to the trailing edge and base material. This is the mechanism most accurately captured by ALE simulation.
• Heat loss to the tool: A portion of frictional heat conducts into the tool material (typically hardened steel or tungsten carbide). This represents a parasitic heat loss that reduces the effective heat available for plasticization of the workpiece.
• Radiation and convection to ambient: Typically negligible at FSW temperatures but becomes significant for high-temperature materials (titanium, superalloys).

3.3 Temperature Regulation Parameters

The solid-state temperature is regulated by the following process parameters, which define the process window for defect-free FSW:

Process Parameter	Typical Range (Aluminum 6061)	Effect on Peak Temperature	Effect on Stir Zone Quality
Rotational Speed (ω)	500–2000 rpm	Strong positive correlation	Higher ω → finer grains but risk of overheating
Welding Speed (vw)	100–500 mm/min	Strong negative correlation	Lower vw → more heat per unit length
Tool Tilt Angle (θ)	1.5°–4°	Moderate effect via shoulder contact area	Optimizes plug formation and flow balance
Plunge Depth (d)	0.1–0.5 mm (pin protrusion)	Positive correlation	Deeper plunge → more material deformation
Tool Shoulder Diameter	8–15 mm (for 3–6 mm thickness)	Positive correlation via friction area	Larger shoulder → broader TAZ
Pin Diameter/Profile	1.5–4 mm; cylindrical, tapered, concave, threaded	Complex interaction	Profile dictates material flow pattern and mixing

3.4 The Role of the Ratio ω/v_w

The ratio of rotational speed to welding speed (ω/v_w) is a master parameter governing the thermal balance in FSW. A higher ratio means more heat is generated per unit length of weld, leading to higher peak temperatures and broader heat-affected zones. Conversely, a lower ratio reduces thermal input, which can lead to insufficient plasticization and defects such as lack of fusion (voids at the trailing side) or incomplete bonding at the root.

For aluminum alloys, the optimal ω/v_w ratio typically falls in the range of 3–8 (in units of s^-1·mm^-1). Below this range, cold defects predominate; above this range, hot defects (tunnel voids, flash, excessive thinning) become likely. The ALE method is particularly valuable for determining the precise boundaries of this window for specific material-tool combinations.

4. Technical Purpose and Value to Cladding Technology Shanxi Co., Ltd.

4.1 Process Development and Optimization

While the company's primary technology routes—TIG/MIG weld overlay, hydraulic explosive bonding, and explosion welding—are established, the understanding of solid-state temperature regulation through FSW simulation provides several strategic advantages:

• Cross-process thermal understanding: The analytical and ALE simulation methodologies developed for FSW are directly transferable to modeling heat-affected zone (HAZ) behavior in TIG/MIG weld overlay processes. The analytical heat source model (double-ellipsoidal) is the industry standard for weld thermal modeling and is directly applicable to predicting dilution, microstructural evolution, and residual stress in overlay welds.
• Process window definition: For explosion welding and explosive bonding, the understanding of thermal-mechanical coupling from ALE simulation informs the design of flyer plate velocities, stand-off distances, and collision angles that optimize bonding interfaces while controlling thermal effects.
• Material compatibility assessment: Simulation-based temperature prediction enables evaluation of whether dissimilar material combinations (e.g., carbon steel to stainless steel, aluminum to copper) can be joined or clad without exceeding critical temperature thresholds that would cause intermetallic compound formation, cracking, or degradation.

4.2 Qualification and Certification Support

Thermal modeling results serve as essential supporting documentation for welding procedure qualification under standards such as:

• ASME Section IX, Part QW-401 through QW-407: Thermal modeling can demonstrate equivalence of welding parameters within qualified ranges.
• ISO 15614-1 / ISO 15614-16: Procedure qualification requires demonstration of mechanical properties; thermal analysis provides predictive validation.
• GB/T 19866 series: Chinese national standards for welder and procedure qualification reference thermal input limits.
• NB/T 20341: Nuclear industry welding procedure standards require thermal analysis for critical components.

4.3 Customer Value Delivery

Thermal simulation capability positions the company as a technically sophisticated partner capable of:

• Providing predictive quality assurance before physical trials, reducing development time and cost
• Optimizing overlay parameters for minimum dilution while maintaining bond strength
• Designing multi-pass overlay sequences that control cumulative thermal exposure
• Resolving customer-specific technical challenges through customized simulation studies
• Supporting design-for-manufacturability reviews with quantified thermal impact assessments

5. Implementation Framework for Thermal Simulation in Cladding Operations

5.1 Analytical Heat Source Model for Weld Overlay

The double-ellipsoidal heat source model, originally developed by Goldak, Akhlaghi, and Butler, is the standard approach for TIG and MIG weld thermal simulation. The heat flux distribution is defined as:

q(x, y) = (6√3 · Q) / (a · b · c · 2π√π) × {exp[−3((x−vx)t/a)² − 3y²/b²] for x ≥ 0 (front); exp[−3((x−vx)t/c)² − 3y²/b²] for x < 0 (rear)}

Where Q is the effective heat input, a, b, c are the ellipsoidal semi-axes, and v_x is the welding speed. This model is directly applicable to TIG overlay welding simulations for predicting:

• Heat-affected zone width and depth
• Dilution rate as a function of travel speed and current
• Peak temperature distribution across the overlay/base material interface
• Cooling rates (t_8/5) governing microstructural transformation

5.2 ALE Method Application to Cladding Processes

The ALE formulation is particularly valuable for modeling the following cladding-related phenomena:

• Explosive welding interface dynamics: The ALE method captures the high-velocity collision between flyer plate and base plate, the jet formation, oxide disruption, and subsequent bonding. Temperature and pressure fields at the collision interface determine bonding quality.
• Multi-pass overlay weld interaction: Each subsequent overlay pass re-heats the previous pass, creating complex thermal histories. ALE simulation tracks the cumulative thermal-mechanical state and predicts residual stress build-up, distortion, and potential cracking.
• Hydraulic explosive bonding process control: The controlled energy release in hydraulic explosive bonding generates specific pressure and temperature profiles at the bonding interface. ALE simulation validates that these profiles fall within the bonding window for the material combination.

5.3 Simulation Workflow

1. Geometry Setup: Create 3D models of workpiece, tool (for FSW) or weld bead (for overlay), with appropriate boundary conditions representing clamping, backing, and cooling.
2. Material Property Definition: Input temperature-dependent thermal conductivity, specific heat, density, flow stress curves, and phase transformation data. For dissimilar material systems, ensure interface properties are correctly assigned.
3. Heat Source Calibration: For analytical models, calibrate heat input efficiency (typically 60–80% for TIG, 70–85% for MIG) against thermocouple measurements or infrared thermography.
4. Mesh Generation: For ALE simulations, create fine mesh in the interaction zone (stir zone, weld pool, collision interface) with element sizes of 0.1–0.5 mm. Apply mesh smoothing/rezoning algorithms to maintain element quality.
5. Solution and Convergence: Solve the coupled thermo-mechanical problem with appropriate time stepping (0.01–0.1 ms for FSW/explosion; 0.1–1.0 s for overlay). Monitor convergence of temperature, stress, and flow fields.
6. Post-Processing and Validation: Extract temperature histories, residual stress distributions, and deformation patterns. Validate against experimental measurements (thermocouples, DIC, XRD residual stress, metallography).

6. Applicable Standards and Acceptance Criteria

6.1 Process Qualification Standards

Standard	Scope	Relevance to Thermal Simulation
ASME Section IX	Welding and Brazing Qualification	Thermal input limits for procedure qualification; simulation supports parameter range justification
ISO 15614-1	Specification for Approval of Welding Procedures	Thermal modeling supports essential variable identification and equivalence demonstration
ISO 15614-16	Welding Procedure Approval – Friction Stir Welding	Directly applicable; requires demonstration of mechanical properties within process window
GB/T 19866.1-2017	Specification for Approval of Welding Procedures – Part 1: General Requirements	Chinese national standard; thermal analysis supports procedure specification development
NB/T 20341-2011	Nuclear Power Plant Welding Procedure Qualification	Requires detailed thermal analysis for critical nuclear components
ASTM E1019	Standard Practice for Determining Dilution in Weld Overlay	Thermal simulation predicts dilution; results validated by this method
API 579-1/ASME FFS-1	Fitting-Up on Damaged Components	Thermal modeling supports repair qualification for in-service cladding repairs

6.2 Material and Performance Standards

• ASTM A240: Stainless steel clad plate specifications – thermal simulation ensures HAZ properties meet minimum requirements
• ASTM A516 / A517: Carbon steel clad plate – predicts base material property retention
• ASTM B497: Nickel-base alloy clad plate – critical for thermal limits to prevent sensitization
• NACE MR0175 / ISO 15156: Materials for H₂S environments – thermal simulation ensures clad material properties remain within NACE-compliant ranges
• GB/T 20388-2006: Clad steel plate specifications – Chinese standard for clad plate delivery
• GB/T 24511-2017: Explosion welding of metals – process specification and acceptance criteria

6.3 Simulation Validation Criteria

For simulation results to be accepted as qualification-supporting evidence, the following validation criteria should be met:

• Peak temperature prediction within ±50 °C of experimental measurement
• Cooling rate (t_8/5) prediction within ±20% of thermocouple-derived values
• HAZ width prediction within ±15% of metallographic measurement
• Residual stress prediction within ±50 MPa of XRD measurement (for critical applications)
• Dilution prediction within ±3% (absolute) of spectrographic analysis

7. Common Risks and Controls

7.1 Simulation-Specific Risks

Risk	Description	Mitigation Control
Material property uncertainty	Temperature-dependent properties are estimated or extrapolated	Use experimentally measured properties; perform sensitivity analysis; validate at multiple temperatures
Boundary condition idealization	Actual cooling conditions (air, water, backing) are simplified	Use calibrated convection coefficients; include back-bar cooling models; validate against IR thermography
Friction coefficient variability	Oxide films, surface roughness, and lubrication affect μ	Use temperature-dependent μ curves from literature; perform parametric studies
Mesh sensitivity (ALE)	Results may depend on element size and smoothing parameters	Perform mesh convergence studies; use adaptive meshing; validate at multiple resolutions
Phase transformation neglect	Latent heat and volume change during phase changes are omitted	Include phase transformation models (Kinetics-based); validate against dilatometry
Scale-up errors	Simulation validated at lab scale may not predict production scale accurately	Include production-specific boundary conditions; account for workpiece temperature pre-heating

7.2 Process Risks Related to Temperature Control

• Overheating in TIG/MIG overlay: Excessive thermal input causes grain coarsening, sensitization (in austenitic stainless steels), and reduced corrosion resistance. Control: Use thermal simulation to define maximum allowable heat input (J/mm) for each material combination.
• Insufficient bonding temperature in explosive welding: Below the minimum collision velocity/temperature threshold, oxide films are not disrupted, resulting in unbonded or partially bonded interfaces. Control: ALE simulation of collision dynamics validates that interface temperatures exceed the bonding threshold (typically 0.5–0.7 T_melt).
• Intermetallic compound formation: In dissimilar metal cladding (e.g., Al/steel, Cu/steel), excessive interfacial temperatures promote brittle intermetallic layers. Control: Analytical thermal modeling identifies maximum permissible dwell time and peak temperature at the interface.
• Residual stress-induced cracking: High thermal gradients generate tensile residual stresses that may exceed material yield strength. Control: ALE simulation predicts stress distributions and identifies critical locations for stress-relief treatment.

8. Application Scenarios Across Company Technology Routes

8.1 TIG/MIG Weld Overlay

The analytical heat source model is directly applicable to TIG and MIG overlay welding process optimization. Key applications include:

• Dilution control: Predicting dilution rates for multi-pass overlay sequences on carbon steel substrates with stainless steel or nickel-alloy overlay layers. The thermal model identifies the minimum number of passes and optimal travel speeds to achieve dilution below 5% (per ASTM E1019).
• HAZ property prediction: For nuclear-grade clad plates (per NB/T 20341), thermal simulation predicts whether HAZ hardness and toughness will remain within acceptable limits after multi-pass overlay.
• Preheat and interpass temperature optimization: Simulation determines the required preheat temperature and maximum interpass temperature to prevent cold cracking in high-hardness base materials while avoiding excessive grain growth in the overlay.
• Distortion prediction: For large-diameter pipe overlay (per ASME B31.3 or API 5L), thermal-mechanical simulation predicts angular and longitudinal distortion, enabling fixture design to minimize post-weld machining.

8.2 Hydraulic Explosive Bonding

The ALE method is essential for modeling the dynamics of hydraulic explosive bonding, where controlled energy release generates the collision conditions required for solid-state bonding:

• Collision velocity optimization: ALE simulation determines the minimum flyer plate velocity (typically 300–700 m/s) required to achieve bonding for specific material combinations. The simulation captures the shock wave propagation, material jet formation, and interface temperature at collision.
• Stand-off distance calibration: The gap between flyer and base plate (typically 0.5–2.0 mm) directly affects collision velocity and angle. ALE simulation identifies the optimal gap for maximum bonding area and minimum oxide inclusion.
• Energy control validation: Hydraulic explosive bonding uses controlled hydraulic pressure to initiate and contain the explosive reaction. Simulation validates that the pressure profile maintains collision conditions within the bonding window throughout the bonding zone.
• Multi-layer bonding sequences: For multi-layer clad plates (e.g., CS/SS/Ni-alloy), ALE simulation models each bonding step sequentially, ensuring that subsequent layers do not degrade previously bonded interfaces through thermal or mechanical loading.

8.3 Explosion Welding

Traditional explosion welding relies on detonation-driven collision. Thermal-mechanical simulation supports the following aspects:

• Collision geometry design: The angle of collision (typically 5°–20°) determines the wave pattern of the bonding interface (wavy, spiral, or laminar). ALE simulation predicts the collision angle as a function of flyer velocity, base plate velocity, and stand-off distance.
• Bonding window determination: For each material pair, there exists an upper and lower boundary of collision velocity above which the material melts (excessive temperature) and below which bonding does not occur. Simulation maps these boundaries precisely.
• Residual stress analysis: The explosive collision generates complex residual stress states in both plates. ALE simulation predicts these stresses, informing post-weld stress relief requirements and dimensional stability predictions.
• Spall and fracture prediction: At excessive collision velocities, material spall or fracture occurs. Simulation identifies the maximum allowable collision conditions for defect-free bonding.

9. Integration with Quality Management and Certification

9.1 WPS Development Support

Thermal simulation results directly inform the Welding Procedure Specification (WPS) development process:

1. Define essential variables and their qualified ranges based on simulation-predicted thermal effects
2. Establish preheat and interpass temperature limits from HAZ property predictions
3. Determine post-weld heat treatment (PWHT) parameters from residual stress analysis
4. Specify NDE requirements based on predicted defect susceptibility at critical locations

9.2 Certification Documentation

For customer qualification packages (particularly in nuclear, aerospace, and oil/gas sectors), simulation reports serve as:

• Supporting evidence for procedure qualification when physical testing is impractical or uneconomical
• Justification for parameter changes within qualified ranges (per ASME IX QW-401.11)
• Demonstration of design intent and process control capability
• Documentation of thermal impact on clad material properties for API 923 or ASME Section VIII compliance

9.3 Digital Twin and Process Monitoring

Advanced implementation includes real-time thermal monitoring where embedded thermocouples or infrared sensors feed measured temperatures into a digital twin model. Deviations from the simulated thermal profile trigger process adjustments or stop conditions, ensuring that the actual thermal history remains within the qualified envelope.

10. Conclusions and Recommendations

The study of FSW solid-state temperature regulation through analytical heat source and ALE simulation methods provides Cladding Technology Shanxi Co., Ltd. with a rigorous analytical foundation that transcends its primary application in friction stir welding. The thermal modeling methodologies, validation frameworks, and process optimization approaches developed through this research are directly transferable to the company's TIG/MIG weld overlay, hydraulic explosive bonding, and explosion welding operations.

Key recommendations for leveraging this capability:

1. Establish a thermal simulation capability center within the company, equipped with validated material property databases and qualified simulation personnel.
2. Develop standard simulation protocols for each technology route, ensuring consistent methodology and traceable validation against experimental data.
3. Integrate simulation into the WPS development workflow as a mandatory step before physical procedure trials, reducing qualification costs and time-to-market.
4. Build a validated simulation database of thermal histories for common material combinations, enabling rapid quotation and technical support for new customer inquiries.
5. Pursue third-party validation of simulation capabilities through participation in benchmark exercises (e.g., ESWIS, NIST weld data) and publication of peer-reviewed validation studies.
6. Extend simulation to multi-scale modeling, coupling macro-scale thermal analysis with micro-scale phase transformation and precipitation models to predict long-term property stability of clad components in service.

By mastering solid-state temperature regulation through both analytical and numerical simulation methods, the company positions itself at the forefront of technologically sophisticated cladding and overlay manufacturing, delivering higher quality, greater reliability, and more predictable performance for demanding industrial applications across nuclear, energy, chemical processing, and transportation sectors.