I wanted to try and take my current knowledge of statistical mechanics (first quarter undergraduate course completed, beginning researcher in far from equilibrium statistical mechanics, basic knowledge of mean field approximation techniques) and apply it to create a rough estimate of corrosion. I want to try and keep this model simple enough that it might be usable in a demonstration to a scientifically/mathematically literate background (not assuming statistical mechanics experience) to showcase some of the interesting ways of applying stat mech to a phenomena.
The tricky part about a problem like this is trying to find the right approximations, so I wanted to try proposing several ideas I had on solving the problem, similar problems I've encountered that can serve as a good inspiration, and difficulties one might encounter. I'll list a few general questions I had about my problem and then follow with my thought process after at the very bottom. Essentially the primary questions I am going to ask are hopefully well encapsulated and (probably) have distinct methodologies to arrive at the conclusion.
What types of values can we determine from statistical methods that can then be used to determine the kinetics of the reaction?
What types of basic model would be a good starting point to work on these calculations?(idea example: would we expect "state condensation" around certain energy levels (something like a bose-einstein model))
More context for what I'm trying to do follows:
Predict (with some accuracy) corrosion
Use variables of interest to be able to then create an approximation of reaction kinetics.
Constants are simple physical constants and basic information available in tables (heats of formations/free energy)
Pure Solid in a well mixed corroding fluid
Simple shapes (spherical or square)
We can represent energy levels from our possible chemical species
Surface rust; weathering operates in a simple, predictable manner, (if this is too complex I'd be willing to make the assumption that it is a weathering resistant material: http://en.wikipedia.org/wiki/Weathering_steel)
The energy level of each possible type of rust is known.
The environment is time-invariant
Difficulties (posed here as research questions as oppose to direct questions meant for SE community, can still be useful look at when trying to answer question 2):
How does geometry effect the calculation of our variables of interest? (example surface rust effecting lower layers? Let's ditch certain types of rusting that are hard to predict like pitting and exfoliating (unless we need to keep them).
Is it justifiable to just use energy levels corresponding to the energy of each chemical species we'd encounter or are there other energy levels to consider?
Would we need to perturb these energy levels as the system evolves?
Idea: Would statistical mechanical models of batteries work well for modelling this system? How would the math of those models need to be adjusted to account for differences in the two models? Or are they completely irreconcilable with each other?