Ultrasonic Propagation Simulation and Defect Echo Prediction in Anisotropic Weld Overlay Structures

1. Definition and Fundamental Principles

1.1 Anisotropy in Weld Overlay Structures

In bimetallic cladding and weld overlay manufacturing, the repeated deposition of dissimilar alloy layers through TIG/MIG welding, hydraulic explosive bonding, or explosion welding creates complex microstructural gradients. These gradients result in anisotropic material properties — meaning acoustic wave velocity, attenuation, and reflection characteristics vary significantly depending on the direction of propagation relative to the weld pass orientation, grain elongation, and layer interfaces. Unlike isotropic homogeneous materials where ultrasonic wave behavior is direction-independent, anisotropic weld overlay structures exhibit:

• Velocity anisotropy: Longitudinal and shear wave velocities differ by up to 5–15% depending on the propagation angle relative to the weld bead direction and layer boundary planes.
• Mode conversion: At inclined interfaces between dissimilar overlay layers, longitudinal waves convert partially into shear waves and vice versa, generating complex echo patterns.
• Beam skew and refraction: The ultrasonic beam path deviates from the normal propagation direction, causing apparent defect position errors if not corrected.
• Attenuation anisotropy: Acoustic energy loss varies with propagation direction due to grain boundary scattering and interlayer reflections.

1.2 Simulation Methodology

The simulation framework for ultrasonic propagation in anisotropic weld overlay structures typically employs:

• Finite Element Method (FEM): Full-wave simulation using software such as COMSOL Multiphysics or ABAQUS, incorporating Christoffel equations for elastic wave propagation in anisotropic media. The stiffness tensor (C_ijkl) is defined based on measured or calculated elastic constants for each overlay layer.
• Ray Tracing Methods: Geometric acoustics approaches (e.g., using CIVA, ONYX, or Zemax) that compute ray paths through layered anisotropic media, accounting for Snell's law at each interface with direction-dependent velocities.
• Transfer Matrix Method: Efficient computation of through-transmission and reflection coefficients across multilayer stacks with anisotropic properties.
• Hybrid Approaches: Combining ray tracing for near-field geometry with FEM for complex scattering from defects and interface irregularities.

The Christoffel equation governs wave propagation in anisotropic media:

C_ijkl · n_j · n_l · u_i = ρ · v² · u_i

where C_ijkl is the fourth-order elastic stiffness tensor, n_j is the propagation direction unit vector, u_i is the displacement amplitude, ρ is density, and v is the phase velocity. For each overlay layer, the stiffness tensor is oriented according to the local material texture (grain orientation).

2. Category and Business Positioning

2.1 Technical Classification

This capability belongs to the Non-Destructive Testing (NDT) and Quality Assurance domain within Cladding Technology Shanxi Co., Ltd.'s technical portfolio. It specifically addresses the critical gap in ultrasonic inspection reliability for multilayer weld overlay structures where conventional NDT techniques suffer from high false-indication rates and unreliable defect sizing due to anisotropic scattering effects.

2.2 Business Positioning

• Quality Assurance Differentiator: Provides scientific justification for NDT acceptance/rejection decisions, reducing customer disputes and warranty claims.
• WPS Qualification Support: Enables demonstration of adequate inspection coverage and defect detection capability during Welding Procedure Specification (WPS) qualification testing per ASME Section IX.
• Customer Confidence Building: Offers predictive simulation reports alongside physical NDT results, demonstrating engineering rigor and traceability.
• Technology Transfer Value: Positions the company as a knowledge leader in advanced NDT methodology for clad products, supporting premium pricing and long-term partnerships.

3. Technical Purpose and Value

3.1 Primary Objectives

1. Defect Echo Prediction: Predict the amplitude, time-of-flight, and waveform shape of ultrasonic echoes from known defect types (porosity, lack of fusion, cracks, interlayer delamination) at specific locations within anisotropic overlay structures.
2. Inspection Planning Optimization: Determine optimal probe selection, beam angle, frequency, scan pattern, and evaluation reference level (ERL) for reliable defect detection in specific clad configurations.
3. False Indication Identification: Distinguish between true defect signals and artifacts generated by layer interface echoes, grain scattering, and geometric beam focusing.
4. Acceptance Criteria Calibration: Establish scientifically defensible amplitude thresholds and sizing methodologies that account for anisotropic attenuation and beam behavior.

3.2 Quantifiable Value

Value Dimension	Impact Metric	Typical Improvement
False positive rate reduction	Percentage decrease in non-relevant indications	40–65% reduction
Inspection coverage confidence	Probability of defect detection (POD)	Increase from 0.65 to 0.90+
Re-inspection cost savings	Reduction in repeat inspections	25–40% cost reduction
Customer acceptance rate	First-pass acceptance of NDT reports	15–25% improvement
WPS qualification cycle time	Time to demonstrate NDT adequacy	30–50% reduction

4. Key Process and Implementation Points

4.1 Material Characterization for Simulation Input

Accurate simulation requires precise knowledge of the anisotropic elastic properties of each overlay layer. The following parameters must be characterized:

Parameter	Measurement Method	Typical Values (309L/316L Overlay on C-Steel)
Longitudinal wave velocity (along weld)	Pulse-echo ultrasonic measurement	5,800–6,200 m/s
Longitudinal wave velocity (perpendicular to weld)	Pulse-echo ultrasonic measurement	5,600–5,950 m/s
Shear wave velocity (along weld)	Shear wave transducer measurement	3,200–3,450 m/s
Shear wave velocity (perpendicular to weld)	Shear wave transducer measurement	3,100–3,350 m/s
Density	Archimedes method / literature	7,900–8,050 kg/m³
Acoustic attenuation coefficient	Through-transmission method	0.5–2.5 dB/mm (frequency-dependent)
Grain size (equivalent)	ET measurement / metallography	ASTM 4–8 (0.2–0.1 mm)
Layer thickness and interface roughness	UT thickness gauge / macrograph	0.5–3.0 mm per pass

4.2 Simulation Workflow

1. Geometry Definition: Model the clad structure geometry including base material, each overlay layer, dilution zone, and any known or assumed defect configurations (flat-bottom holes, side-drilled holes, cracks).
2. Material Property Assignment: Assign the anisotropic stiffness tensor to each layer based on measured wave velocities and crystallographic texture data.
3. Probe Modeling: Define transducer parameters including center frequency, bandwidth, aperture, focal length, wedge material, and angle. Include beam divergence and near-field effects.
4. Boundary Condition Setup: Define excitation source, coupling medium (water/air), and radiation conditions at model boundaries.
5. Solver Configuration: Select appropriate time-step (for FEM) or ray density (for ray tracing) to ensure convergence and computational efficiency.
6. Simulation Execution: Run the simulation for the defined scan path and record time-amplitude-depth (TAD) curves and A-scan/B-scan/C-scan outputs.
7. Validation: Compare simulation predictions against physical NDT measurements on qualification test specimens (V-blocks, reference blocks with artificial defects).
8. Iterative Refinement: Adjust material properties and model assumptions until simulation predictions match experimental data within acceptable tolerances (typically ±10% in amplitude, ±5% in TOF).

4.3 Defect Echo Prediction Parameters

Defect Type	Typical Size (Reference)	Predicted Echo Amplitude (as % of Reference)	Key Distinguishing Features
Porosity (spherical)	Ø 2 mm	40–70%	Single peak, narrow pulse width
Lack of fusion (planar)	5 × 0.5 mm	70–100%	Broad pulse, possible multiple echoes
Crack (interlayer)	10 × 0.1 mm	85–120%	Sharp peak, high frequency content
Slag inclusion	Ø 1.5 mm	25–55%	Multiple scattered echoes
Interlayer delamination	20 × 0.2 mm	90–130%	Strong reflection, mode-converted echoes

4.4 Validation Against Physical Testing

Simulation predictions must be validated against physical ultrasonic testing of qualification specimens containing artificial reference defects. Acceptance criteria for validation include:

• Time-of-flight prediction accuracy: within ±5% of measured TOF
• Amplitude prediction accuracy: within ±10 dB of measured amplitude
• Defect position prediction accuracy: within ±2 mm lateral and ±1 mm depth
• Waveform shape correlation: cross-correlation coefficient > 0.85

5. Applicable Standards and Acceptance Criteria

5.1 Governing Standards

Standard	Scope	Relevance to Simulation Work
GB/T 11345-2013	Ultrasonic testing of welds — General rules	Defines UT inspection procedures, equipment requirements, and acceptance methodology for weld overlay joints
GB/T 12606-2010	Ultrasonic testing of welds — Procedures	Specifies scan patterns, reference levels, and defect evaluation criteria applicable to overlay welds
NB/T 47013.3-2015	Non-destructive testing of pressure vessels — Ultrasonic testing	Pressure vessel specific UT requirements including acceptance criteria for weld overlay on pressure equipment
ASME BPV Section V, Article 4	Ultrasonic Examination	Defines UT methods, equipment calibration, and acceptance for pressure vessel weld overlay
ASME BPV Section VIII, Div. 2	Rules for Construction — Alternative Rules	Requires demonstration of NDT capability including POD/POD analysis for clad pressure vessels
ASME Section IX, QW-452	Welding procedure qualification — NDT	Requires proof of adequate NDT coverage during WPS qualification
ASTM E164/E164M	Standard practice for calibrating UT equipment	Reference block calibration procedures used to validate simulation predictions
ASTM E2309/E2309M	Standard practice for UT of welds using phased array	Phased array UT methodology for complex geometries including overlay welds
ISO 17640:2017	UT of welds — General rules	International framework for UT inspection planning and execution on weld overlay
ISO 13588:2016	UT of welds — Procedures	Detailed scan procedures and acceptance criteria for weld inspection
NACE SP0775	UT inspection of weld overlay cladding	Corrosion protection specific NDT requirements for overlay cladding
API 579-1/ASME FFS-1	Fitness-for-service assessment	Requires UT data for defect characterization in service evaluation of clad components

5.2 Acceptance Criteria for Simulation Results

• Defect detection capability: Simulation must demonstrate detection of all reference defects at or above the relevant standard's minimum detectable size (typically 0.5 mm planar or 1.5 mm volumetric for most standards).
• Signal-to-noise ratio: Predicted defect echo amplitude must exceed the evaluation reference level by at least 6 dB for reliable detection.
• False indication rate: Simulation-predicted interface echoes and scattering artifacts must be distinguishable from defect echoes with documented criteria.
• Positioning accuracy: Predicted defect location must be within the standard's acceptable tolerance for sizing (typically ±2 mm for planar defects).

6. Common Risks and Controls

6.1 Technical Risks

Risk	Description	Mitigation Control
Over-simplification of anisotropy	Using isotropic assumptions for inherently anisotropic overlay layers leads to significant prediction errors	Mandatory anisotropic characterization of each layer; minimum 3-directional velocity measurements; texture analysis where feasible
Interface modeling error	Incorrect representation of layer boundaries (perfect bonding vs. partial delamination) distorts wave propagation predictions	Include realistic interface impedance models; validate against known interface reflection data; include roughness parameters
Grain scattering underestimation	Coarse-grained weld overlay deposits cause high scattering that may not be captured by smooth-boundary models	Incorporate stochastic scattering models; include grain size and shape distribution in attenuation calculations
Probe-to-surface coupling variability	Surface roughness of overlay welds (typically Ra 25–63 μm) affects coupling efficiency unpredictably	Model multiple surface roughness scenarios; use worst-case coupling assumptions for acceptance criteria
Computational convergence failure	Complex 3D anisotropic models may fail to converge or produce numerical artifacts	Use adaptive mesh refinement; validate 2D cross-sections before full 3D modeling; use multiple solver configurations
Extrapolation beyond validated range	Applying simulation results to geometries or materials outside the validated parameter space	Document validation boundaries explicitly; require re-validation for new geometries exceeding ±20% parameter variation

6.2 Quality Risks

• Uncalibrated simulation software: Ensure simulation tools are validated against known analytical solutions and physical test data before production use. Maintain calibration records per ISO 17025 or equivalent.
• Personnel competency: Simulation operators must demonstrate competency in both ultrasonic physics and numerical modeling. Maintain training records per NB/T 47014 or ASME Section IX qualification requirements.
• Document control: All simulation models, input parameters, and output reports must be version-controlled and traceable to specific production orders or WPS qualifications.

7. Application Across Company Technology Routes

7.1 TIG/MIG Weld Overlay

In TIG and MIG weld overlay processes, multiple thin passes (typically 0.5–2.0 mm per pass) are deposited to build up the cladding layer. This creates highly anisotropic structures with:

• Strong directional grain elongation along the travel direction, creating pronounced velocity anisotropy (up to 8–12% difference between longitudinal and transverse directions).
• Multiple interpass interfaces that generate reflection and mode-conversion echoes, complicating defect detection.
• Progressive dilution gradients from base material to pure overlay alloy, creating continuous property transitions rather than sharp interfaces.

Simulation Application: For TIG/MIG overlay, the simulation model typically includes 5–15 discrete layers representing individual weld passes. Each layer is assigned anisotropic properties based on the specific alloy composition at that dilution level. The simulation helps determine optimal phased array scan configurations (e.g., 64-element linear array, 2.5 MHz center frequency, 45°–70° angle range) that provide adequate coverage of the dilution zone and interpass boundaries where lack-of-fusion defects are most likely.

Key Value: Enables WPS qualification by demonstrating that the proposed UT procedure can detect reference defects (typically 10 mm × 0.2 mm side-drilled hole or equivalent) at all critical locations within the overlay structure, satisfying ASME Section IX QW-452 and NB/T 47013.3 requirements.

7.2 Hydraulic Explosive Bonding

Hydraulic explosive bonding (also known as hydraulic explosive welding or water-jet explosive bonding) creates clad structures through a controlled explosive reaction in a water-filled chamber, achieving metallurgical bonding between base and cladding materials without melting. The resulting structure exhibits:

• Wave-pattern interface morphology with characteristic amplitude and wavelength that creates complex acoustic impedance variations.
• Severe plastic deformation zones near the interface with highly distorted grain structures causing intense scattering.
• Residual stress gradients from the explosive forming process that affect acoustic wave propagation through elastic modulus variations.

Simulation Application: The simulation must model the undulating interface geometry and the highly deformed near-interface microstructure. This is critical for distinguishing between:

1. True bonding defects (unbonded areas, voids at the interface)
2. Interface geometry echoes from the wave pattern
3. Scattering from the deformed grain structure

Simulation results guide the selection of inspection parameters — typically requiring lower frequency (1–2 MHz) with longer wavelength to penetrate the scattering zone while maintaining sufficient resolution for defect detection. The simulation predicts the amplitude reduction due to interface scattering and helps set appropriate evaluation reference levels that avoid excessive false indications from the bonding interface itself.

Key Value: Provides the scientific basis for acceptance criteria in hydraulic explosive bonding NDT, particularly for distinguishing acceptable interface characteristics from true bonding failures. Supports compliance with ASTM F2541 and relevant GB standards for explosion-welded clad products.

7.3 Explosion Welding (Traditional Contact Method)

Traditional explosion welding involves direct contact of cladding and base plates at high velocity (typically 2,000–4,000 m/s) during detonation. The resulting clad structure has distinct characteristics compared to hydraulic explosive bonding:

• More pronounced wave-pattern interface with larger amplitude and wavelength (typically 0.5–5 mm amplitude, 5–20 mm wavelength).
• Higher residual stress levels and more extensive deformation zones.
• Potential for micro-cracking in the deformation zone, particularly at wave crests.
• Thicker deformation zones (up to 0.5–1.0 mm on each side of the interface).

Simulation Application: For explosion-welded clad products, the simulation must account for the larger geometric complexity of the interface. Key applications include:

1. Predicting through-transmission signals to determine minimum detectable unbonded area size
2. Modeling the interaction of UT beams with wave-pattern interfaces to establish acceptable echo patterns vs. defect indications
3. Simulating the effect of the deformation zone's anisotropic grain structure on shear wave propagation for angle-beam inspection
4. Optimizing inspection approach angles to minimize interface interference while maintaining defect coverage

Key Value: Enables reliable qualification of explosion-welded clad plates and pipes per ASTM A240/A270/A516 clad product specifications and ASME SA-240/SA-270 clad plate requirements. The simulation supports the demonstration that NDT procedures can reliably detect minimum-size bonding defects (typically 3 mm diameter unbonded area per common industry acceptance criteria) despite the complex acoustic environment.

8. Contribution to Qualification Building, Product Delivery, and Customer Value

8.1 Qualification Building

• WPS Qualification Support: Simulation results provide objective evidence that the proposed NDT procedure achieves required detection capability, reducing the number of physical qualification coupons needed and accelerating WPS approval timelines.
• Personnel Qualification: Simulation exercises serve as training tools for NDT Level III personnel, enabling competency assessment in complex overlay inspection scenarios without requiring physical test specimens for every scenario.
• Equipment Qualification: Simulation validates that specific UT equipment (frequency, aperture, wedge configuration) meets the detection requirements for specific clad product configurations before equipment procurement or deployment.
• Standard Compliance: Provides documented evidence of NDT capability for ASME Section VIII Div. 2 qualification requirements, which mandate demonstration of adequate inspection coverage through analytical methods.

8.2 Product Delivery Enhancement

• Faster Inspection Turnaround: Simulation-guided inspection parameters reduce the need for iterative physical testing and re-inspection, accelerating product delivery schedules by 20–35%.
• Reduced Scrap Rates: Accurate defect prediction reduces false-positive rejection of acceptable products, directly improving manufacturing yield and reducing material waste.
• Confidence in First-Pass Quality: Simulation-validated inspection procedures increase confidence that delivered products meet specification requirements, reducing post-delivery quality disputes.
• Traceability Documentation: Simulation reports provide a complete technical record of the NDT methodology, supporting traceability requirements for critical applications (nuclear, pressure vessels, aerospace).

8.3 Customer Value

• Engineering Confidence: Customers receive simulation-backed NDT reports that demonstrate scientific rigor beyond simple pass/fail ultrasonic testing, building long-term trust and repeat business.
• Regulatory Compliance Support: Simulation documentation supports customer regulatory submissions to NORSOK, NQA-1, ASME, or other regulatory bodies that require demonstrated NDT capability.
• Cost Optimization: By reducing false indications and re-inspection requirements, simulation-guided NDT reduces total project costs while maintaining or improving quality assurance levels.
• Technology Leadership: Offering simulation-based NDT services positions Cladding Technology Shanxi Co., Ltd. as a technology leader in the clad products industry, supporting premium pricing and differentiation from competitors offering only conventional NDT services.
• Customized Inspection Solutions: Simulation enables tailored NDT procedures for each unique product configuration, rather than applying generic inspection procedures that may be suboptimal for specific clad geometries and materials.

9. Implementation Roadmap

9.1 Phase 1: Foundation (Months 1–3)

1. Acquire and validate simulation software (COMSOL, CIVA, or equivalent)
2. Establish anisotropic property measurement protocols for all standard overlay alloys used
3. Develop and validate baseline simulation models for common clad configurations (309L/316L on carbon steel, Hastelloy C-276 on stainless steel)
4. Train 2–3 personnel to Level II competency in simulation methodology

9.2 Phase 2: Integration (Months 4–6)

1. Integrate simulation into WPS qualification workflow for all three technology routes
2. Develop standardized simulation report templates aligned with relevant standards
3. Validate simulation predictions against physical NDT data from production orders
4. Establish simulation-as-a-service capability for customer presentations

9.3 Phase 3: Optimization (Months 7–12)

1. Expand simulation database to cover all alloy combinations and geometries in the product portfolio
2. Develop automated workflow linking CAD geometry directly to simulation models
3. Implement machine learning enhancement for rapid prediction based on database patterns
4. Publish technical papers and present at industry conferences to build thought leadership

10. Conclusion

Ultrasonic propagation simulation and defect echo prediction in anisotropic weld overlay structures represents a critical technical capability that bridges the gap between theoretical NDT methodology and practical inspection reliability for complex clad products. By providing scientifically rigorous prediction of ultrasonic behavior in anisotropic overlay structures, this capability directly enhances qualification efficiency, product quality assurance, and customer confidence across all three of Cladding Technology Shanxi Co., Ltd.'s manufacturing technology routes — TIG/MIG weld overlay, hydraulic explosive bonding, and explosion welding.

The investment in this simulation capability yields measurable returns through reduced false indication rates, accelerated WPS qualification cycles, improved product acceptance rates, and enhanced market positioning as a technology-driven clad products manufacturer. As industry standards increasingly require demonstrated NDT capability through analytical methods (particularly ASME Section VIII Div. 2), this capability transitions from a competitive advantage to a qualification necessity for access to premium markets in nuclear, petrochemical, and aerospace applications.