Corrosion Modleing in Comsol
Corrosion Modleing in Comsol
Corrosion Modleing in Comsol
Motivation
Why Magnesium in Automobiles?
Light weight increases fuel economy High strength to weight ratio makes it a good engineering material Good castability enables high production rate
Challenges
Poor corrosion resistance
Mg is not very suitable for automotive applications which are exposed to corrosion media
High galvanic corrosion
Mg is very anodic to other metals such as Fe, Al poor durability in dissimilar material
joints and applications
Research Focus
Al
2e-
Mg
Cathodic
Objective: To optimize the processing route to obtain the desirable microstructure for the better corrosion resistance
Objective: To develop a numerical tool to predict galvanic corrosion rate, which can provide design specifications
AM50 DC
AM50 PM
AM50 SC
Microstructure effect on corrosion b phase fraction b phase Distribution Al content in the a phase India Science Lab
Comsol Conference, Bangalore, October 29 - 30, 2010
Modeling approach
Model Schematic n 0
Electrolyte
Assumptions
Well mixed, incompressible and electroneutral electrolyte solution Uniform Al content throughout alpha phase No interfacial resistance due to corrosion product is considered here The dissolution reaction is considered at the anode surface
n 0
2 0
b phase (Cathode) a phase (Anode)
n 0
j k n f ( )
Corrosion rate in m.sec-1 mole.Amp-1.sec-1 Equivalent weight in kg.mole-1 Density in kg.m-3 anodic current density in Amp.m-2
The model is capable of explicitly tracking the corroding phase. The model predictions are based on non-linear polarization data.
Electrolyte
Cathode
Anode
Deformed mesh
x
Comsol Conference, Bangalore, October 29 - 30, 2010
50% b phase
90% b phase
as b phase fraction
Schematic
Vin RE
Potentiostat WE CE RE
Sample specimen
j
b a b
E A
Current density, Amp/m2 Conductivity, S/m Potential difference at the two extremes of amplitude of probe vibration, V Amplitude of probe vibration, m
E
SVET experiments are performed on a galvanic couple where individual components are in direct physical and electrical contact, which eliminates IR drop.
50% b phase
90% b phase
Model and SVET experiments : Corrosion rate as b phase fraction The corrosion rate AM50 dc was the lowest despite the highest b phase fraction This behavior motivates us to investigate the effect of b phase distribution
India Science Lab
Comsol Conference, Bangalore, October 29 - 30, 2010
A single boundary is used to specify the anode surface and the cathode surface using the level set function The above formulation can capture the varying ratios of the a phase and the b phase along the depth of the alloy
Acknowledgement: Shashank Tiwari for microstructure
Electrolyte
b phase
Dissolved a phase
a phase preferentially dissolves in electrolyte solution until continuous network of b phase is exposed to electrolyte solution after which corrosion is halted.
India Science Lab
Comsol Conference, Bangalore, October 29 - 30, 2010
a phase dissolution
b phase
India Science Lab
Dissolved a phase
Average anodic current density as surface b phase fraction The corrosion activity is however halted after the a phase present on the surface is dissolved into the electrolyte solution. The total charge passed per square meter area = 1.3844105 C/m2 The total mass of material removed per square meter area = 17.22 g/m2
Discrete b phase
Model prediction Schematic of microstructure Electrolyte
b phase
Dissolved a phase
a phase continuously dissolves in the electrolyte solution as discontinuous b phase network assists corrosion. As discrete b phase is not well supported, b phase eventually spatters off.
References: Ambat, et al., Corrosion Science, 42, 1433-1455 (2000); Song et al., Corrosion Science, 41, 249-273 (1999).
Discrete b phase
Surface b phase fraction evolution
a b a b
Neck formation and neck thinning leading to b phase spattering is captured using the model. The corrosion activity lasts for a lot longer for discrete b phase than in continuous b phase network. The model captures the scenario leading to b phase spattering, but not the actual spattering as the model cannot capture topological changes.
The total charge passed per square meter area = 2.5060 10 C/m2
5
The total mass of material removed per square meter area = 31.17 g/m2
Objective: To optimize the processing route to obtain the desirable microstructure for the better corrosion resistance
Objective: To develop a numerical tool to predict galvanic corrosion rate, which can provide design specifications
Modeling approach
Model Schematic n 0
Electrolyte
Assumptions
Well mixed, incompressible and electro-neutral electrolyte solution No interfacial resistance due to corrosion product is considered here Anode surface is assumed to be corroding and cathode surface is assumed to be non-corroding
n 0
2 0
Mild Steel (Cathode) Magnesium (Anode)
n 0
j k n f ( )
The model is capable of explicitly tracking the corroding phase. The model predictions are based on non-linear polarization data.
Mild Steel
2 mm
The corrosion rate predicted using the numerical model is within +9% of that estimated using SVET analysis for AE44 mild steel couple. The corrosion rate predicted using the numerical model is within -29% of that estimated using SVET analysis for AE44 AA6063 couple.
References: K. B. Deshpande, Corrosion Science, 52 (2010) 2819 2826 K. B. Deshpande, Corrosion Science, 52 (2010) 3514 - 3522 Comsol Conference, Bangalore, October 29 - 30, 2010
The corrosion rate predicted using the numerical model is within -20% of that estimated from immersion technique for AE44 mild steel couple. The corrosion rate predicted using the numerical model is within -47% of that estimated from immersion technique for AE44 AA6063 couple.
Acknowledgement: Sampath Vanimisetti (Image digitization)
Summary
Corrosion rate (mm/y) Galvanic Couples SVET experiments 197 42
Immersion experiments
243 52
The corrosion rate predicted using the numerical model is in good agreement with that estimated from the two experimental techniques.
The galvanic corrosion behavior is demonstrated at both macro and micro scale using the numerical model developed using COMSOL Multi-Physics.