Scientific journal
Fundamental research
ISSN 1812-7339
"Перечень" ВАК
ИФ РИНЦ = 1,674

STATISTICAL MODEL OF EDUCATION CORROSION SPOTS ON THE METAL

Platonova E.S. 1 Buchinskas V. 2 Yurov V.M. 3
1 Karaganda State Technical University
2 Vilnyussky Technical University named after Gediminas
3 Karaganda State University named after E.A. Buketov
The paper deals with the statistical model of the formation of corrosion spots on the metal when the oxygen molecule is near the reactive surface of the metal atom with a certain probability. The case, when this probability is defined as the ratio of the number of particles in the reactive zone atoms to the total number of particles in a selected circumferentially around this atom certain radius. Shows the logarithmic dependence of the area of corrosion spots on the metal surface imperfection. The case where the probability of hitting the oxygen molecule near the reactive metal atom is the ratio of the binding energy of the oxygen molecule to the metal atom to the total energy of formation of corrosion spots. Shows the logarithmic dependence of the area of corrosion spots on the chemical equilibrium constant of the oxidation process. Obtained in formulas contain experimentally determined values.
corrosion
corrosion spot
metal
oxidation
defect
1. Azarenkov N.A., Litovchenko S.V., Nekljudov I.M., Stoev P.I. Korrozija i zashhita metallov. Chast 1. Himicheskaja korrozija metallov. Harkov: HNU, 2007. 187 р.
2. Vatanabje S., Ikeda N. Stohasticheskie differencialnye uravnenija i diffuzionnye processy. M.: Nauka, 1986. 448 р.
3. Dynkin E.B., Jushkevich A.A. Teoremy i zadachi o processah Markova. M.: Nauka, 1967. 426 р.
4. Ito K., Makkin G. Diffuzionnye processy i ih traektorii. M.: Mir, 1968. 480 р.
5. Kidjarov B.I. Kinetika obrazovanija kristallov iz zhidkoj fazy. Novosibirsk: Nauka, 1979. 134 р.
6. Skripov V.P. Metastabilnaja zhidkost. M.: Nauka, 1972. 312 р.
7. Stratonovich R.L. Nelinejnaja neravnovesnaja termodinamika. M.: Nauka, 1985. 480 р.
8. Baroux B. The kinetics of pit generation on stainless steel // Corrosion science. 1988. Vol. 28, no. 10. рр. 969–986.
9. Kondo J. Prediction of fatigue crack initiation life based on pit growth // Corrosion science. l989. Vol. 45, no. 1 рр. 7–11.
10. Provan J.W., Rodriguez E.S. Development of a Markov description of pitting corrosion. Corrosion (USA). 1989. Vol. 45, no. 3. рр. 178 192.
11. Strutt I.E., Nicholls and Barbier B. The prediction of corrosion by statistical analysis of corrosion profiles // Corrosion science. 1985. Vol. 25, no. 5. рр. 305–3l6.

Вопросам коррозии металлов посвящено огромное количество работ, из которых отметим лишь работы [1, 9–11], в которых приводится обширная библиография. Несмотря на это работы в области теории коррозии продолжают расти с ростом различных типов конструкционных металлических материалов, применяемых в различных областях промышленного производства.

В самом общем случае коррозию металла можно представить как зарождение и рост новой фазы (окисленного металла). Критический зародыш новой фазы образуется последовательно в серии случайных актов присоединения и отрыва атомов (молекул) друг от друга. Поэтому зародышеобразование – случайный процесс во времени и пространстве. Это предопределяет вероятностный характер параметров, которые описывают кинетику образования зародышей в процессе коррозии или роста кристалла [11].

Марковские процессы и кинетика образования зародышей новой фазы

Количественное описание случайного процесса дается его функцией распределения, удовлетворяющей кинетическому уравнению. В общем случае кинетическое уравнение является сложным интегро-дифференциальным уравнением, решить которое невозможно. Однако, если рассматривать случайный процесс как марковский, то кинетическое уравнение переходит в дифференциальное, которое имеет более простой вид [10].

В частном случае пуассоновского процесса гибели и размножения с конечным числом состояний получается система дифференциальных уравнений [5–6]:

platon01.wmf (1)

Здесь λ0 – вероятность перехода системы из состояния E0 в E1 и т.д.; μ1 – вероятность перехода из состояния E1 в E0 и т.д. Вероятность перехода из En в En–1 полагается равной нулю (μn = 0), т.е. состояние En для такой системы – поглощающее; pi(t) – вероятность нахождения системы в состоянии Ei. Система уравнений (1) является системой уравнений Колмогорова [3]. Общее решение настолько громоздко, что не позволяет провести анализ даже численными методами.

Приведенный выше пример показывает, что большинство исследователей идет по пути решения уравнений диффузионного типа, математическая теория которых разработана достаточно полно и которые описывают процесс случайных (броуновских) блужданий, а также теории надежности, порядковых статистик, массового обслуживания и ряд других.

К числу основных кинетических уравнений, вытекающих из дифференциальных уравнений Колмогорова, для случайных марковских процессов относятся уравнения Смолуховского – Чепмена и Фоккера – Планка. В случае многомерного вектора состояния platon02.wmf, уравнение Фоккера – Планка имеет вид [9]

platon03.wmf (2)

Получить решение уравнения (2) в явном виде удается лишь в частных случаях, например, если величина K линейна по переменным q, а величина Qjk не зависит от q. Уравнение (2) имеет важную роль в нестационарной теории замедления нейтронов, однако его линеаризация приводит к тому, что его решение сводится к отысканию соответствующей функции Грина при различных граничных условиях.

Для решения стохастических дифференциальных уравнений диффузионного типа, к которым также относятся уравнения (1) и (2), развито несколько подходов, из которых необходимо отметить подходы Ито [4] и Стратоновича [7]. Решая записанные для процесса Маркова y(t) стохастические уравнения типа Ланжевена, можно находить соответствующие ему нестационарные или стационарные моменты и корреляторы, одновременные или многовременные стохастические уравнения Ито [7]:

platon04.wmf (3)

где ξ(s)(t) – дельта-коррелированные по времени случайные функции, имеют недостатком то, что с входящими в него выражениями типа σ(y)ξ(t) при различных выкладках нельзя обращаться по обычным правилам, пригодным для гладких функций. Это обстоятельство неудобно с практической точки зрения.

Из сказанного выше отметим следующее:

– формальный вероятностный характер описания процесса новой фазы как случайного марковского стационарного или нестационарного процесса Пуассона требует более строгого обоснования, поскольку ряд предположений (отсутствие последействия, ординарность и др.) не вытекает из физической картины явления коррозии;

– так как экспериментально наблюдаемыми являются (как и в квантовой теории) не сами функции распределения, а их моменты и корреляторы, то целесообразно проведение анализа статистических закономерностей начальной стадии коррозии

Статистическая модель образования коррозионного пятна

Наиболее распространенным случаем коррозии металла является его взаимодействие с молекулами кислорода. Очевидно, что это взаимодействие начинается со «слабых» мест поверхности металла или с его дефектами. Мы изложим нашу модель образования коррозионного пятна с позиции статистической физики.

Рассмотрим поверхность металла с числом дефектов m. Пусть расстояние между дефектами одинаково и равно R. Опишем вокруг каждого дефекта 0 окружность радиусом R. Пусть плотность числа дефектов в этой окружности равна n0, тогда вероятность W0(r) того, что ближайшая частица кислорода попадет на расстояние r от частицы 0, нетрудно получить на основе статистической физики и она равна

platon05.wmf (3)

Вероятность нахождения N0 частиц кислорода в зоне дефекта 0 радиусом r равна, очевидно,

platon06.wmf (4)

Вероятность (4) определим, с другой стороны, как отношение числа частиц N0 в зоне дефекта к общему числу частиц в выделенной окружности – Q0:

platon07.wmf (5)

Для системы из m дефектов имеем

platon08.wmf (6)

Для всего металла с числом дефектов 0, 1, 2, …, m имеем

platon09.wmf (7)

Cистема уравнений (6) и (7) представляет собой систему трансцендентных уравнений, решить которую можно только приближенными или численными методами.

В связи с этим можно сделать численную оценку, основываясь на реальной ситуации и уравнении (1) системы (6):

platon10.wmf

Соответствующая оценка дает, что первый член левой части уравнения и первые два члена правой части – пренебрежимо малы. В результате получим

platon11.wmf (8)

Учитывая, что πr2 = S – площади коррозионного пятна и n0N0 = const, из (8) имеем

platon12.wmf (9)

Последнее выражение показывает логарифмическую зависимость площади коррозионного пятна от «дефектности» поверхности металла.

Учет энергии коррозионного процесса

Вероятность (4) можно определить, с другой стороны, и как отношение энергии связи молекулы кислорода Е0 с атомом металла к общей энергии образования коррозионного пятна Еобщ. Учитывая, что [1]

platon13.wmf (10)

где Kp – константа химического равновесия, проводя, аналогичные приведенным выше вычисления, для площади коррозионного пятна, получим

platon14.wmf (11)

Заключение

Несмотря на простоту формул (9) и (11) они могут быть полезными для исследования процессов коррозии любых конструкционных материалов, поскольку включают в себя легко определяемые в эксперименте параметры.

Работа выполнена по программе МОН РК 055 «Научная и/или научно-техническая деятельность», подпрограмма 101 «Грантовое финансирование научных исследований».

Рецензенты:

Портнов В.С., д.т.н., профессор, начальник УПО, Карагандинский государственный технический университет, г. Караганда;

Жакатаев Т.А., д.т.н., старший преподаватель кафедры ММ и Н, Карагандинский государственный технический университет, г. Караганда.

Работа поступила в редакцию 10.04.2015.