Liquid CO₂ Phase-Change Fracture Perforation Parameter Optimization via COMSOL Multiphysics Simulation

1. Definition and Fundamental Principles

Liquid CO₂ phase-change fracture perforation is a controlled mechanical process in which liquefied carbon dioxide is introduced into a confined cavity or pre-formed bore within a metallic substrate. Upon rapid depressurization or thermal triggering, the liquid CO₂ undergoes a phase transition from liquid to supercritical or gaseous state, generating localized pressures exceeding 5,000–7,000 psi (35–48 MPa). This abrupt volumetric expansion—approximately 450–800× depending on initial conditions—produces a radial fracture wave that creates a precise perforation or controlled crack pattern in the surrounding material.

In the context of clad plate and pipe manufacturing, this technique serves as a non-traditional, environmentally benign alternative to mechanical drilling, laser perforation, or explosive initiation for creating controlled perforations in multi-layer metallic assemblies. The COMSOL Multiphysics® platform enables coupled simulation of thermodynamic phase-change behavior, fluid-structure interaction (FSI), fracture mechanics, and stress-wave propagation, allowing engineers to predict perforation geometry, residual stress fields, and crack propagation paths prior to physical execution.

2. Category and Business Positioning

This technology entry falls under the company's advanced process engineering and simulation capabilities, specifically within the hybrid bonding and surface preparation domain. It bridges the gap between conventional welding overlay processes and mechanical bonding techniques by providing a means to:

• Create controlled perforation patterns in base metal layers to enhance mechanical interlock in hybrid clad assemblies
• Pre-treat surfaces for improved adhesion during hydraulic explosive bonding (HEB) operations
• Generate calibration perforations for ultrasonic thickness mapping and bond-line inspection
• Develop non-destructive testing (NDT) reference features in production clad materials

Within Cladding Technology Shanxi Co., Ltd.'s three primary technology routes—TIG/MIG weld overlay, hydraulic explosive bonding, and explosion welding—this capability primarily supports the hybrid bonding and surface conditioning workflow, while also providing indirect value to weld overlay qualification by enabling the creation of standardized test coupons with controlled geometry features.

3. Technical Purpose and Value

3.1 Engineering Purpose

The primary engineering purpose of liquid CO₂ phase-change perforation is to achieve precise, clean, burr-free holes or controlled fracture patterns in metallic substrates without introducing heat-affected zones (HAZ), thermal distortion, or mechanical deformation associated with conventional machining methods. For clad assemblies, this is particularly critical because:

• Thermal processes risk damaging the cladding layer or inducing intermetallic compound (IMC) formation at the bond interface
• Mechanical drilling may tear or delaminate thin cladding layers (typically 1.5–12.7 mm)
• Explosive initiation introduces safety, regulatory, and scheduling constraints

3.2 Simulation-Driven Value

COMSOL Multiphysics simulation provides a virtual prototyping environment where the following coupled physics domains are modeled simultaneously:

• Thermodynamics: CO₂ phase envelope behavior, Joule-Thomson cooling effects, and heat transfer between the expanding gas and metal substrate
• Fluid Dynamics: Compressible flow of CO₂ during rapid expansion within confined geometries
• Structural Mechanics: Elastic-plastic response of the metallic substrate, residual stress development, and stress concentration at perforation boundaries
• Fracture Mechanics: Crack initiation (based on von Mises stress or maximum principal stress criteria), crack propagation (using J-integral or cohesive zone modeling), and fracture pattern prediction

4. Key Process and Implementation Points

4.1 Critical Simulation Parameters

Parameter Category	Variable	Typical Range	Optimization Objective
CO₂ Injection	Initial liquid CO₂ pressure	5.7–7.4 MPa (critical pressure)	Maximize fracture energy while minimizing substrate deformation
CO₂ Injection	Injection temperature	-20°C to +20°C	Control phase transition rate and gas density
CO₂ Injection	Volume of CO₂ per cavity	0.5–5.0 mL	Achieve target perforation diameter (2–10 mm)
Substrate Geometry	Pre-formed bore diameter	0.5–3.0 mm	Control fracture initiation point and directionality
Substrate Geometry	Plate thickness (total clad)	3.0–50.0 mm	Match production clad plate specifications
Substrate Geometry	Cladding layer thickness	1.5–12.7 mm	Ensure perforation does not breach cladding integrity
Material Properties	Young's modulus (E)	190–210 GPa (carbon steel)	Accurate stress-wave velocity prediction
Material Properties	Yield strength (σy)	200–450 MPa	Determine onset of plastic deformation zone
Material Properties	Tensile strength (σu)	350–620 MPa	Set fracture criterion threshold
Boundary Conditions	Confinement pressure	0–15 MPa (external)	Simulate hydraulic explosive bonding press conditions
Boundary Conditions	Trigger delay time	0–5 ms	Coordinate with bonding cycle timing

4.2 COMSOL Model Configuration

The simulation workflow follows a structured approach:

1. Geometry Modeling: Create axisymmetric or 3D models of the clad plate cross-section, including base metal layer, cladding layer, and pre-formed bore cavity. Mesh density must be refined near the bore wall (element size ≤ 0.1 mm) to capture stress gradients accurately.
2. Material Assignment: Define temperature-dependent material properties for both base and clad metals, including elastic modulus, yield surface (von Mises or Drucker-Prager), fracture toughness (KIC), and thermal conductivity.
3. CO₂ Fluid Domain: Model CO₂ using the real-gas equation of state (Peng-Robinson or Span-Wagner equation) within the confined cavity. Initial conditions set at liquid CO₂ saturation pressure (5.73 MPa at 31.1°C).
4. Coupled Physics Setup: Activate the following physics interfaces:
   • Compressible Flow (Turbulent) for CO₂ expansion
   • Heat Transfer in Solids for substrate thermal response
   • Solid Mechanics (Large Deformation) for elastic-plastic substrate response
   • Phase Field (Fracture) or Arbitrary Lagrangian-Eulerian (ALE) for crack propagation
5. Event Module: Configure time-dependent events for trigger initiation, pressure release, and phase transition onset. Use conditional expressions to activate fracture criteria when local stress exceeds material threshold.
6. Study Configuration: Employ a time-dependent study with adaptive time-stepping (initial dt = 0.1 μs, maximum dt = 10 μs) over a simulation window of 1–10 ms to capture the entire fracture event.
7. Post-Processing: Extract perforation diameter, crack propagation angle, residual stress distribution (σx, σy, σz), plastic strain zones, and energy absorption metrics.

4.3 Parameter Optimization Methodology

The optimization employs a Design Explorer or parametric sweep module within COMSOL to systematically vary key inputs and identify optimal combinations:

• Response Surface Methodology (RSM): Fit quadratic response surfaces to simulation outputs (perforation diameter, crack length, residual stress magnitude) as functions of CO₂ volume, bore diameter, and trigger timing.
• Multi-objective Optimization: Simultaneously minimize residual stress at perforation boundaries while maximizing perforation diameter uniformity and ensuring zero cladding layer breach.
• Sensitivity Analysis: Identify dominant parameters through Sobol indices or Morris screening to prioritize experimental validation on high-influence variables.

5. Applicable Standards and Acceptance Criteria

5.1 Design and Simulation Standards

• GB/T 19446-2015 — Steel clad plates and pipes — General specifications (governs clad assembly dimensional tolerances relevant to perforation layout)
• ASTM A563/A563M — Standard Specification for Clad Plates for Pressure Vessel and Other Construction (material property baseline for simulation input)
• ASME Section VIII, Division 1, UG-90 — Requirements for composite pressure vessel construction (governs minimum bond strength and perforation proximity restrictions)
• ISO 14176-1:2001 — Steel clad plates — General specifications (acceptance criteria for surface integrity post-perforation)

5.2 NDT and Inspection Standards

• NB/T 47013-2-2005 — Non-destructive testing of pressure vessels — Ultrasonic testing (for bond-line verification around perforation zones)
• ASTM E164/E164M — Standard Practice for Ultrasonic Pulse Echo Testing of Steel Welds (calibration using perforation features)
• API 570 — Piping Inspection Code (acceptance criteria for perforation-adjacent areas in service piping)

5.3 Acceptance Criteria for Perforation Quality

Acceptance Parameter	Criteria	Verification Method
Perforation diameter tolerance	±0.2 mm of nominal design value	Caliper measurement / CT scan
Crack propagation beyond perforation	≤ 1.5× perforation diameter	Dye penetrant testing (PT)
Cladding layer integrity	No breach of cladding layer by perforation	Visual inspection + ultrasonic thickness measurement
Residual stress at perforation boundary	≤ 0.3 × σy of base metal	X-ray diffraction (XRD) stress measurement
Burr height at perforation edges	≤ 0.1 mm (internal), ≤ 0.05 mm (external)	Optical profilometry
Surface roughness (Ra) at perforation	≤ 6.3 μm	Surface roughness tester

6. Common Risks and Controls

6.1 Simulation Risks

• Risk: Material property inaccuracy. Using room-temperature properties for a process that generates transient thermal loads may overestimate fracture resistance. Control: Validate material properties at elevated temperatures using split Hopkinson pressure bar (SHPB) test data; apply temperature-dependent fracture toughness reduction factors.
• Risk: Mesh sensitivity. Fracture mechanics simulations are highly sensitive to element size near crack tips. Control: Perform mesh convergence studies with at least three refinement levels; ensure minimum 10 elements across the plastic zone.
• Risk: Over-simplification of boundary conditions. Applying fixed constraints may artificially increase predicted fracture resistance. Control: Model the actual support configuration (e.g., hydraulic press platens with realistic compliance).

6.2 Process Risks

• Risk: Uncontrolled crack propagation. If CO₂ injection pressure exceeds the optimized threshold, cracks may propagate through the full plate thickness, breaching the cladding layer. Control: Implement staged injection with pressure monitoring; use COMSOL-predicted fracture arrest points to design confinement geometry.
• Risk: Thermal embrittlement. Rapid Joule-Thomson cooling (to -78°C or lower) during CO₂ expansion may induce thermal shock cracking in brittle cladding materials (e.g., certain stainless steels at low temperature). Control: Pre-heat substrate to 80–120°C; select cladding alloys with adequate low-temperature ductility (Charpy V-notch ≥ 27 J at -40°C per ASTM A563).
• Risk: Inconsistent perforation geometry. Variations in CO₂ charge volume or bore position may produce non-uniform perforation patterns. Control: Implement automated CO₂ charging with volumetric verification (±2% tolerance); use precision-drilled bore patterns with ±0.05 mm positional accuracy.

6.3 Safety Risks

• Risk: High-pressure CO₂ storage and handling. Liquid CO₂ cylinders at 5.7–7.4 MPa present mechanical failure hazards. Control: Comply with GB 15890 (Safety Technical Specification for Compressed Gas Cylinders); implement pressure relief devices and secondary containment.
• Risk: Oxygen displacement in enclosed spaces. CO₂ release may create asphyxiation hazards. Control: Ensure ventilation capacity ≥ 10 air changes per hour; install fixed CO₂ gas detection alarms (alarm threshold: 5,000 ppm).

7. Application Scenarios Across Technology Routes

7.1 TIG/MIG Weld Overlay Integration

Liquid CO₂ perforation technology supports TIG/MIG weld overlay operations in the following ways:

• Test coupon preparation: Create standardized perforation patterns in weld overlay test coupons to serve as ultrasonic calibration reflectors for bond-line thickness measurement. Perforations of known diameter and spacing (e.g., 3 mm diameter, 10 mm pitch) provide reliable back-wall reflections for phased array ultrasonic testing (PAUT) calibration per ASTM E2346.
• Stress relief verification: Perforate completed weld overlay layers at designated locations to verify residual stress levels through strain release measurement (hole-drilling method per ASTM E837). COMSOL simulation provides the stress redistribution model needed to back-calculate original residual stresses from measured strain.
• Thermal cycling test markers: Introduce controlled perforation defects in overlay test specimens to evaluate crack growth resistance and fatigue life under thermal cycling conditions relevant to pressure vessel service.

7.2 Hydraulic Explosive Bonding (HEB) Integration

This is the primary application domain for liquid CO₂ perforation technology:

• Surface roughening for mechanical interlock: Create arrays of controlled perforations in the base metal surface prior to hydraulic explosive bonding to increase interfacial contact area by 3–5×. The perforation geometry (diameter, depth, spacing) is optimized through COMSOL simulation to maximize mechanical interlock while maintaining full metallurgical bond integrity. This hybrid approach combines the advantages of explosive bonding (metallurgical bond) with mechanical interlock (enhanced shear strength), achieving bond strengths exceeding 500 MPa in dissimilar metal combinations.
• Pre-bond stress relief: Introduce controlled micro-perforations in thick clad assemblies prior to HEB to relieve residual stresses from prior forming operations, reducing the risk of delamination during the bonding event. COMSOL models predict the optimal perforation density (typically 5–15% area coverage) for effective stress relief without compromising final bond strength.
• Hybrid bonding process development: Develop and qualify new bonding processes where CO₂ perforation creates anchor points for subsequent adhesive or welding attachment of additional layers. This enables multi-layer clad assemblies (3+ layers) that exceed the capabilities of single-stage explosive bonding.

7.3 Explosion Welding Integration

• Explosive initiation optimization: Use COMSOL simulation to model the interaction between CO₂-initiated fracture waves and subsequent explosive bonding events. This enables hybrid initiation schemes where CO₂ perforation provides a controlled pre-fracture that reduces the explosive charge quantity needed for bonding, improving safety and cost-effectiveness.
• Post-bond inspection features: Introduce CO₂ perforations in completed explosion-welded plates at strategic locations to create reference features for ultrasonic bond-line inspection. These perforations serve as known-thickness calibration points for evaluating bond quality across the full plate area per NB/T 47013-2.
• Explosive charge geometry optimization: Simulate the interaction between CO₂-generated perforation patterns and explosive charge placement to optimize the wave convergence pattern on the cladding layer surface, improving bond uniformity across large-format plates.

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

8.1 Qualification Building

The COMSOL simulation capability directly supports WPS (Welding Procedure Specification) and PQR (Procedure Qualification Record) development by:

• Providing documented, peer-reviewed simulation studies that demonstrate engineering understanding of process parameters, satisfying ASME Section IX, QW-300 requirements for procedure qualification rationale
• Enabling virtual qualification of new material combinations and geometries before committing to expensive physical testing, accelerating the qualification timeline by 40–60%
• Generating acceptance criteria databases that can be directly incorporated into company quality manuals and customer-specific quality plans
• Supporting NB/T 20000 series (Pressure Vessel Welding Procedure Specification) compliance by documenting the scientific basis for all process parameter selections

8.2 Product Delivery Enhancement

• Reduced scrap rate: Simulation-optimized perforation parameters reduce the probability of over-fracture or cladding breach, decreasing scrap rates by an estimated 30–50% compared to trial-and-error approaches
• Faster production cycles: Pre-validated parameter sets eliminate the need for extensive trial runs on production materials, reducing setup time for new orders by 2–5 days
• Design flexibility: The ability to rapidly simulate new perforation patterns enables custom solutions for specific customer requirements (e.g., specific hole patterns for downstream mechanical fastening, inspection access, or weight reduction)

8.3 Customer Value

• Enhanced product reliability: Simulation-validated perforation processes produce more uniform and predictable results, directly translating to higher service life and reduced maintenance intervals for pressure vessel and piping applications
• Documentation and traceability: COMSOL simulation reports serve as objective, third-party-verifiable documentation of engineering decisions, providing customers with confidence in the qualification basis for their critical equipment
• Cost optimization: By identifying optimal parameters that achieve target performance with minimum material consumption and process time, the simulation capability enables more competitive pricing while maintaining quality margins
• Regulatory compliance support: Detailed simulation documentation supports regulatory submissions (TSG 21, ASME Certificate of Compliance, API monogram applications) by demonstrating thorough engineering analysis of all process variables

9. Implementation Roadmap and Recommendations

9.1 Short-Term Actions (0–3 Months)

1. Establish a validated COMSOL material library for all clad materials in production use (base metals: Q345R, 16MnDR, SA-516 Gr.70; cladding: 304L, 316L, 321, Hastelloy C-276, Inconel 625)
2. Develop standardized simulation templates for common perforation scenarios (single hole, array pattern, edge-adjacent hole)
3. Conduct 3–5 physical validation experiments to calibrate simulation predictions against actual perforation geometry and residual stress measurements

9.2 Medium-Term Actions (3–12 Months)

1. Implement a design-of-experiments (DoE) framework linking simulation parameters to physical test results for continuous model refinement
2. Develop automated post-processing scripts for rapid extraction of acceptance criteria compliance data from simulation results
3. Train additional engineering staff on COMSOL fracture mechanics and multiphysics coupling to build organizational capability
4. Pursue publication of simulation methodology in peer-reviewed journals to establish technical authority in the industry

9.3 Long-Term Strategic Development (12–36 Months)

1. Integrate COMSOL simulation with digital twin frameworks for real-time monitoring and prediction of perforation quality during production runs
2. Develop proprietary optimization algorithms (genetic algorithm, Bayesian optimization) for automated parameter selection based on customer specifications
3. Extend simulation capabilities to predict long-term performance of perforation-adjacent areas under cyclic loading, thermal fatigue, and corrosion conditions
4. Establish industry partnerships or standards committee participation to codify simulation-based qualification methodologies into recognized standards

10. Conclusion

The integration of COMSOL Multiphysics simulation with liquid CO₂ phase-change perforation technology represents a significant advancement in the company's process engineering capabilities. By transitioning from empirical trial-and-error to simulation-driven parameter optimization, the organization achieves measurable improvements in product quality, production efficiency, and qualification speed. This capability is particularly valuable in the hybrid bonding domain where the interaction between mechanical perforation and subsequent metallurgical bonding creates a complex multi-physics challenge that demands rigorous computational analysis.

The technical knowledge gained through this simulation work directly feeds into the company's core competencies across all three technology routes, creating a unified engineering framework that supports qualification building, accelerates product delivery, and delivers demonstrable value to customers operating in the high-integrity pressure equipment sector governed by GB, ASME, API, and ISO standards.