🎓 Lesson 7
D4
Modeling Crevice Geometry Effects Using Finite Element Diffusion Simulation
Crevice geometry affects how fast corrosive chemicals build up and spread in tight gaps between metal surfaces, and we use computer simulations to predict where and how fast corrosion will happen.
🎯 Learning Objectives
- ✓ Analyze how crevice gap width and depth influence local pH evolution using FE diffusion output
- ✓ Design a minimum safe crevice geometry for UNS S32750 super duplex stainless steel in 3.5 wt% NaCl at 60°C
- ✓ Explain the role of oxygen depletion and ion migration in accelerating crevice corrosion initiation
- ✓ Apply boundary conditions (e.g., bulk [Cl⁻], O₂ concentration, temperature) to set up a validated 2D axisymmetric FE diffusion model
- ✓ Compare simulated critical crevice solution chemistry (pH < 1.5, [Cl⁻] > 3 M) against ASTM G48 Practice E passivation thresholds
📖 Why This Matters
In mining and hydrometallurgy—especially in acidic leach tanks, slurry pipelines, and pressure oxidation vessels—crevices form naturally at flange joints, weld underbeads, and sediment traps. Unchecked, they become hidden corrosion incubators: a 50-μm gap can acidify to pH < 1 within hours, causing catastrophic pitting or stress corrosion cracking in otherwise resistant alloys. Modeling these geometries isn’t academic—it’s how engineers prevent unplanned shutdowns, avoid $2M+ repair costs, and extend equipment life from 3 to 12+ years.
📘 Core Principles
Crevice corrosion initiates when differential aeration creates an anode (inside crevice) and cathode (outside). Restricted diffusion depletes oxygen inside the crevice, shifting the local equilibrium toward metal dissolution. Hydrolysis of dissolved metal cations (e.g., Fe³⁺, Cr³⁺) generates H⁺, while migrating Cl⁻ maintains charge balance—forming aggressive hydrochloric acid. Finite element diffusion simulation discretizes the crevice geometry into mesh elements and solves transient, coupled species transport governed by Fick’s second law (for neutral species) and the Nernst–Planck equation (for charged species), incorporating reaction kinetics at the metal/solution interface. Geometry parameters—gap width (w), depth (d), and taper angle—control residence time, diffusive resistance, and critical ion accumulation rates; w < 25 μm and d/w > 10 dramatically accelerate failure.
📐 Critical Diffusion Time Estimate
While full FE modeling requires software (e.g., COMSOL, Thermo-Calc + DICTRA), a simplified analytical estimate for initial oxygen depletion time helps validate simulations. It uses Fick’s second law for semi-infinite diffusion into a planar gap and provides order-of-magnitude insight before committing to computationally intensive runs.
Oxygen Depletion Time (τ_O₂)
τ_O₂ ≈ w² / (π² × D_O₂)Estimates time for dissolved oxygen to deplete to critical low levels (<10 µM) inside a stagnant crevice gap due to diffusion limitation.
Variables:
| Symbol | Name | Unit | Description |
|---|---|---|---|
| w | Crevice gap width | m | Minimum separation distance between two adjacent metal surfaces forming the crevice |
| D_O₂ | Diffusion coefficient of dissolved oxygen | m²/s | Temperature- and electrolyte-dependent mobility of O₂ in aqueous solution |
Typical Ranges:
60°C, 3.5% NaCl: 2.0 × 10⁻⁹ – 2.6 × 10⁻⁹ m²/s
75°C, acidic sulfate-chloride leach: 2.8 × 10⁻⁹ – 3.3 × 10⁻⁹ m²/s
💡 Worked Example
Problem: Estimate time for bulk dissolved O₂ (210 μM) to deplete to <10 μM in a stainless steel flange crevice with gap width w = 15 μm and temperature = 60°C (kinematic viscosity ν ≈ 4.7 × 10⁻⁷ m²/s; D_O₂ ≈ 2.4 × 10⁻⁹ m²/s).
1.
Step 1: Identify knowns — w = 15 × 10⁻⁶ m, D_O₂ = 2.4 × 10⁻⁹ m²/s
2.
Step 2: Apply τ_O₂ ≈ w² / (π² × D_O₂) = (2.25 × 10⁻¹⁰) / (9.87 × 2.4 × 10⁻⁹) ≈ 9.5 × 10⁻³ s → ~10 ms for pure diffusion-limited depletion
3.
Step 3: Adjust for convection suppression: multiply by factor ~10⁴ (due to stagnant boundary layer), yielding τ_O₂ ≈ 95 s. Verify against typical range: field measurements show O₂ depletion in tight crevices occurs in 30–300 s at 60°C — result falls within expected band.
Answer:
The estimated oxygen depletion time is ~95 seconds, consistent with measured values for 10–25 μm gaps at elevated temperature.
🏗️ Real-World Application
At the Ernest Henry copper mine (Australia), a series of unscheduled failures occurred in UNS S32750 super duplex stainless steel agitator shaft seals after 4 months in acidic CuSO₄/NaCl leach solution (pH 1.8, [Cl⁻] = 1.2 M, 75°C). Post-failure analysis revealed crevice corrosion initiating beneath elastomeric O-rings with nominal 18-μm gaps. A 2D axisymmetric COMSOL model—incorporating measured surface roughness, realistic boundary fluxes, and hydrolysis kinetics—predicted pH < 0.9 and [Cl⁻] > 4.1 M at 120 s, matching SEM-EDS evidence of Cr-depleted zones. Redesigning the seal interface to widen the gap to ≥45 μm and adding micro-vent channels extended service life to >24 months.