Математика. Физика
УДК 543.42
Б01: 10.17277/уе81тк.201б.04.рр.ббб-673
ВОЗМОЖНОСТЬ ТЕОРЕТИЧЕСКОГО УЧЕТА ВЛИЯНИЯ КИСЛОТ НА РЕЗУЛЬТАТЫ ОПРЕДЕЛЕНИЯ МАГНИЯ МЕТОДОМ ПЛАМЕННОЙ ЭМИССИОННОЙ ФОТОМЕТРИИ
В. И. Барсуков, О. С. Дмитриев, А. А. Барсуков
Кафедра «Физика», ФГБОУВО «ТГТУ», г. Тамбов, Россия; рку?@ nnn.tstu.ru
Ключевые слова: аэрозоль; вязкость; интенсивность излучения; поверхностное натяжение; растворы; расход; твердые частицы.
Аннотация: Показана возможность теоретического учета мешающего влияния кислот в пламенно-фотометрическом методе анализа по физико-химическим параметрам анализируемых растворов. Выведена формула для определения доли раствора, поступающего в плазму пламени. Отмечено совпадение расчетных и экспериментальных величин.
Обозначения
А, Ь, I, т и г - постоянные; В - постоянная, зависящая от свойств газа и параметров установки; С и СА - концентрации соответственно анализируемого элемента и кислоты; Б0 и Б - средний и текущий диаметры частиц аэрозоля;
dM/Л - интенсивность испарения вещества с поверхности одной частицы; I,10 - интенсивность излучения для кисло-тосодержащих и не содержащих кислоты растворов соответственно; К - коэффициент, зависящий от химической природы кислоты; кр - коэффициент, характеризующий пропускную способность камеры распылителя;
1р - длина капилляра;
Мн, N - число частиц аэрозоля, осевших в камере, и распыляемых в единицу времени соответственно; ^ - число частиц, попадающих в пламя; п - число частиц данного диаметра; пп - число частиц, прошедших камеру; Р(Б) - вероятность появления частиц аэрозоля данного диаметра; АР - скорость газа относительно раствора; 0> - расход раствора;_
Qn - количество раствора, попадающего в пламя в единицу времени; г - удельная теплота парообразования раствора;
гк - радиус капилляра; АТ - разность между температурами пламени и поверхности испарения; и - подъемная сила распылителя; V - расход газа;
Уп - объем частиц данного диаметра; у, у, - практические и теоретические значения определяемой величины соответственно;
а - коэффициент теплообмена; 5 - доля раствора, поступающего в пламя; у - отношение излучения атомов анализируемого элемента в присутствии кислоты к их излучению при отсутствии в растворе кислоты;
П - динамическая вязкость раствора; X - коэффициент теплопроводности газа; ст - среднеквадратическое отклонение; ст0, р - поверхностное натяжение и плотность раствора соответственно;
- вероятность появления частиц данного диаметра;
- критерий разделения.
Введение
Очевидно, что кислоты, содержащиеся в анализируемых растворах, искажают определение элементов по их эмиссионным спектрам в пламенно-фотометрическом методе (ПФМ) [1]. Поэтому в настоящей работе сделана попытка аналитически связать некоторые физико-химические свойства растворов, такие как плотность, вязкость, поверхностное натяжение, кислотность с величиной мешающего действия на определение магния. Рассмотрен вопрос о возможности предварительного определения доли анализируемого раствора, поступающего в пламя.
Исследования проводили на спектрофотометре, конструкция которого описана в [1, 2]. Пламя - воздушно ацетиленовое. Аналитическая линия для магния 285,2 нм. Указаны методики и приборы для определения названных параметров анализируемых растворов. При определении ряда физических величин применили метод наименьших квадратов [3]. Расчеты производились на персональном компьютере.
Теоретический расчет влияния неорганических кислот на определение магния ПФМ
Экспериментальное определение мешающего действия неорганических кислот, таких как азотная, соляная, серная и фосфорная на определение ряда элементов ПФМ достаточно подробно проведено в [1, 2]. Показано, что это, в основном, связано с изменением физико-химических свойств растворов, поступающих в пламя с помощью пневматического распылителя.
В настоящей работе сделана попытка теоретического учета данного влияния. Величина аналитического сигнала, в данном случае интенсивности излучения, зависит от числа частиц, поступающих в пламя, и скорости их испарения. Эту зависимость можно представить в виде
'=4 О <1>
к р 60
где Ы„ = 3 , а интенсивность испарения вещества с поверхности одной частицы определяется как [4]
Ш =а^ дг (2)
С г
При испарении капель малых размеров коэффициент теплообмена определяется равенством а = 2Х / £>0. Таким образом, уравнение (2) примет вид
i А
Ш = 2^, дг, а уравнение (1) I = / к_ Щм или I = /0/Д?).
Ж г и гл2
Количество анализируемого раствора, проходящего через распылитель в единицу времени, может быть представлено либо расчетным методом [1, с. 38] по
држг4
формуле Пуазейля 0 =-, либо определено экспериментально, а средний
8/п
диаметр частиц [1, с. 45] по формуле
* = fj* + 5«'
(
\0,45
iooog у-
V
1,5
п
Расчеты показали, что при содержании магния в анализируемых растворах в количестве 0,05; 0,10 и 0,20 мг/мл степень его ионизации в пламени пропан-бутан-воздух составляет 0,031; 0,021 и 0,015 соответственно, то есть незначительна и ею можно пренебречь. Тогда равенство (1) можно записать в виде
I = I |Ч—ССАк р )
IБ2 СТ0П )
Исходя из этого, относительное изменение интенсивности излучения атомов магния в пламени можно записать
I AQCЬ ксА (3)
У = Т- = 2 т ^ е А . (3)
Т г\2_т
j0 D ст0 n j
0
Тогда для раствора, не содержащего кислоту, интенсивность излучения атомов магния в пламени будет определяться выражением
j = AQ0Cb
в общем случае
AQCb
Значения коэффициентов в уравнении определяли методом наименьших квадратов [3], согласно которому £ (yj - yt )2 = min, а значит и
д(£ (у - y )2) 0,
дх
где yt можно заменить на /0, а yt на I.
Прологарифмировав последнее выражение (4) и взяв частные производные по всем параметрам и приравняв их нулю, получаем систему из четырех уравнений с четырьмя неизвестными A, b, l, m. Решая эту систему уравнений, получаем
значения: A = 5,82-106, b = 0,75, l = 0,51 и m = 0,21. Введя в уравнение (3) величину
т n2 ml
J = J0D ст 5,
AQCb
k с Z
получим новое его выражение J = e р А , прологарифмировав которое, имеем
7
ln J = kрСА. Вторичное логарифмирование этого выражения дает lnln J = = ln k р + Z ln С а. Вновь используя метод наименьших квадратов, получаем систему уравнений:
n ln kp + Z £ ln С а = £ ln ln J,
ln kp £ ln СА + Z £ ln2 СА = £ (lnln J) ln СА.
Решение данной системы дает значение параметра kp для азотной кислоты, равное 0,252; соляной - 0,198; серной - 0,0478 и фосфорной - 0,0524. Значение параметра Z для всех кислот оказался практически постоянной величиной равной 0,335. Очевидно, что параметр kp для всех кислот отрицателен и, по-видимому,
0 1 2 3 4 СА, н 0 1 2 3 4 СА, н
а) б)
0 1 2 3 4 Са, н 0 1 2 3 4 СА, н
б) г)
Рис. 1. Экспериментальное и теоретическое изменения интенсивности излучения магния в присутствии азотной (0,10 мг/мл М§) (а), серной (0,10 мг/мл М£) (б), соляной (0,05 мг/мл М§) (в), фосфорной (0,20 мг/мл М§) (г) кислот в растворе:
о - экспериментальное определение у; А - теоретическое определение у (вычисленное по (5))
связан с активностью кислот в растворах, характер изменения которых схож с изменением свойств самих растворов.
С учетом значений найденных коэффициентов формула (3) примет вид
5,82-10б 6С0,75 -крСА'335 У 10 Д2ст°,21п°,51 '
Из анализа приведенного рис. 1 сделан вывод, что экспериментальные и вычисленные по формуле (5) значения у находятся в хорошем соотношении.
Теоретическое определение доли раствора, поступающего в пламя
Как показали экспериментальные данные [5], при прохождении частиц аэрозоля через камеру распылителя часть их осаждается на стенках распылителя, а часть проходит в горелку. Критерий разделения для потока частиц, набегающих на стенку, с одной стороны, имеет вид
¥ = ЛрШ2, (6)
с другой стороны, он равен
Т = Мн /Мп. (7)
Величину Мн можно представить как
МН = М - Мп. (8)
Подставив выражение (8) в (7), получим Т = (М - Мп)/ Мп, откуда
Мп = М /(1+ Т).
Подсчитаем число частиц, образуемых распылителем в единицу времени. Известно, что распределение частиц по размерам подчиняется закону Гаусса, математическая запись которого имеет вид
( Р - Р)2
1--2—
Р( Р) = ^= е 2°2 . л/2яст
Однако применение этого закона приводит к появлению интеграла, который невозможно вычислить. Поэтому, в первом приближении заменим закон Гаусса равномерным распределением, имеющим следующий вид:
0, -<»<Р < Р0 -3ст; Р(Р) = |с, Р0 -3ст<Р <Р0 + 3ст; 0, Р0 + 3ст < Р
Р0 +3ст
Из теории вероятности известно, что | Р(Р)йР = 0,995 и 1. Величину
Р0 -3ст
ст можно определить из следующих соображений. Распылитель не может образовывать частицы диаметром меньше или равным нулю в силу правила физической реализуемости. Тогда, учитывая, что кривая Гаусса симметрична относительно Р0, можно записать 3ст = Р0 + Ртт, и приняв Ртт = 0, получим
ст = Рр/3. (9)
Р0 +3ст Р0 +3ст
С учетом (9) получим | Р(Р)ёР = ¿Р = 1 или, подставив в пределы Рр-3ст Рр-3ст
Ртах 1
интегрирования значение (9), ¡Р = 1, интегрируя которое, имеем ^ =-•
0 Ртах
С другой стороны, л = Цтп и ^, откуда
ММ
п- = М / Р (10)
"г 14 ' ^тах- 4 '
Ртах
Расход раствора определяется следующим образом д = |п.уг ¿¡Р.
о
Подставив в это равенство значения пг и Уг., и интегрируя его, получаем
, г Ртах
д Г пР3аР = 1.МРтах,
д Ртах 0 6 24 max,
3
или, преобразовав его с учетом того, что Ртах = 2Рр, имеем М = 3д / пРр.
Объем раствора, попадающего в пламя в единицу времени определим по формуле
Вшах у-^3
Q = f ^nndD, (11)
J 6 0
где пп можно представить через критерий разделения как пп = щ /(1 + Т). С учетом формулы (10) оно примет вид
пп = N-1-. (12)
П Дшах(1 + Т)
Подставляя формулу (6) в (12), а (12) в (11), с учетом сделанных выше замечаний, имеем
nN 1 f tdt
Qn = 6DmaX 2 f 1 + ApUt
Этот интеграл является табличным и, вычислив его, получаем
6п = ( 26 У2 - 1п(1 + 4АрШо2)].
(4 АрЩ2 )
2
Обозначив 4А/^ = Б, Т = БриВо получим окончательно
8= Щп = — [- 1п(1 + у)]. (13)
б
Параметр Б определен экспериментально Б = 0,088.
Анализ проведенных исследований показал, что экспериментальные и вычисленные по формуле (13) значения доли раствора, поступающего в пламя, находятся в хорошем соответствии (табл. 1).
Таблица 1
Вычисленные и экспериментально найденные доли кислотосодержащего раствора магния, поступающего в пламя
CS т о н тграция ы, моль £ ?г Плотность раствора, кг/м3 Поверхностное натяжение, мН/м е 2 я s § т 8 2 m х ^CJ Н CS X £ & S й иа , Доля раствора, в пламени, о4 (U о4 о, онт ,ние а и
с & ° о У « S U g э и оа « S И « о CS См о t g и ft О IS « Н rt ^ и & о о ^ <D ¡Т н СЕ W и S § Я и <р л ^ R с g о н и m о к соло ю 5 < 1
1 2 3 4 5 6 7 8 9 10 11
Азотная 0,0 0,5 1,0 2,0 4,0 0,10 0,998 1,019 1.029 1.030 1,128 71,61 69,43 71,05 67,77 66,80 1,033 1,026 1,025 1,024 1,060 0,058 72.1 68.7 69.2 67,5 62.8 7,02 7,04 6,89 7,20 8,32 7,87 6,99 8,21 8,61 7,55 +0,35 -0,06 +1,32 +1,41 -0,77
Продолжение табл. 1
1 2 3 4 5 6 7 8 9 10 11
Соляная 0,5 1,0 2,0 4,0 0,05 1,009 1,016 1,033 1,068 71,84 71,72 71,78 71,13 1,047 1,068 1,102 1,172 0,058 66,8 66,8 65.6 64.7 7.47 7.48 7,56 7,52 8,17 7,59 6,48 7,39 +0,70 +0,11 -1,08 -0,13
Серная 0,5 1,0 2,0 4,0 0,10 1,018 1,037 1,075 1,150 76.57 76,51 76,77 77.58 1,062 1,091 1,145 1,292 0,058 70,4 70,0 66,9 66,0 7,32 7,23 7,22 7,51 8,30 8,39 7,00 8.78 +0,98 +1,16 -0,22 +1,27
Фосфорная 0,5 1,0 2,0 4,0 0,20 1,008 1,018 1,038 1,066 71,75 71,87 71,92 73,22 1,075 1,105 1,171 1,329 0,078 68.3 68,5 68,8 65.4 8,25 8,40 8,17 8,73 8,98 7,98 8,59 9,00 +0,73 -0,42 +0,42 +0,27
Заключение
Таким образом, зная физико-химические параметры анализируемого раствора, можно предварительно определить характер изменения интенсивности излучения при определении магния в кислотосодержащих растворах методом пламенной фотометрии и выбрать условия работы спектрофотометра для получения наибольшей чувствительности.
Список литературы
1. Барсуков, В. И. Пламенно-эмиссионные и атомно-абсорбционные методы анализа и инструментальные способы повышения их чувствительности / В. И. Барсуков. - М. : Машиностроение-1, 2004. - 171 с.
2. Планирование эксперимента при определении микроэлементов в водных растворах методом атомно-абсорбционной спектроскопии / В. И. Барсуков [и др.] // Вестн. Тамб. гос. техн. ун-та. - 2016. - Т. 22, № 1. - С. 114 - 121.
3. Налимов, В. В. Логические основания планирования эксперимента / В. В. Налимов, Т. И. Голиков. - 2-е изд., перераб. и доп. - М. : Металлургия, 1981. - 152 с.
4. Лыков, М. В. Сушка распылением / М. В. Лыков. - М. : Пищепромиздат, 1955. - 124 с.
5. Брицке, М. Э. Анализ металлургических продуктов методом эмиссионной фотометрии пламени / М. Э. Брицке. - М. : Металлургия, 1969. - 224 с.
Possibility of Theoretical Consideration of the Effect of Acids on Results of the Magnesium Determination by Flame Photometric Method
V. I Barsukov, O. S. Dmitriev, A. A. Barsukov
Department "Physics", TSTU, Tambov, Russia; phys@ nnn.tstu.ru
Keywords: consumption; radiation intensity; solid particles; solutions; spray; surface tension; viscosity.
Abstract: The paper shows the possibility of theoretical consideration of the interfering effect of acids in flame photometric method for the analysis of physical-
chemical parameters of the solutions. The formula for determining the proportion of solution entering the flame plasma is obtained. Reasonably good agreement between the calculated and experimental values is observed.
References
1. Barsukov V.I. Plamenno-emissionnye i atomno-absorbtsionnye metody analiza i instrumental'nye sposoby povysheniya ikh chuvstvitel'nosti [Flame emission and atomic absorption methods of analysis and instrumental methods of increasing their sensitivity], Moscow: Mashinostroenie-1, 2004, 171 p. (In Russ.)
2. Barsukov V.I., Grebennikov M.V., Dmitriev O.S., Emel'yanov A.A., Barsukov A.A. [Experimental Design for Identification of Microelements in Aqueous Solutions by Atomic Absorption Spectroscopy], Transactions of Tambov State Technical University, 2016, vol. 22, no. 1, pp. 114-121, doi: 10.17277/vestnik.2016.01.pp.114-121 (In Russ., abstract in Eng.)
3. Nalimov V.V., Golikov T.I. Logicheskie osnovaniya planirovaniya eksperimenta [The logical design of experiment base], Moscow: Metallurgiya, 1981, 152 p. (In Russ.)
4. Lykov M.V. Sushka raspyleniem [Spray drying], Moscow: Pishchepromizdat, 1955, 124 p. (In Russ.)
5. Britske M.E. Analiz metallurgicheskikh produktov metodom emissionnoi fotometrii plameni [Analysis of metallurgical products by emission flame photometry], Moscow: Metallurgiya, 1969, 224 p. (In Russ.)
Möglichkeit der theoretischen Berücksichtigung des Einflusses der Säuren auf die Ergebnisse der Bestimmung des Magnesiums von der Methode der flammenden Emissionsphotometrie
Zusammenfassung: Es ist die Möglichkeit der theoretischen Berücksichtigung des störenden Einflusses der Säuren in der flammenden photometrischen Methode der Analyse nach den physikalisch-chemischen Parametern der analysierenden Lösungen gezeigt. Es ist die Formel für die Bestimmung des Anteiles der Lösung, die ins Plasma der Flamme einfließt, herausgeführt. Es ist das Zusammenfallen der Rechen- und Experimentalgrößen festgestellt.
Possibilité du mesurage théorique de l'influence des acides sur les résultats de la détermination du magnésium par la méthode de la photométrie émissive flambante
Résumé: Est montrée la possibilité théorique du mesurage théorique de l'influence nuisible de l'influence des acides dans la méthode photométrique émissive flambante de l'analyse par les paramètres physico-chimiques des solutions analysées. Est donnée la formule pour déterminer la part de la solution entrant dans le plasma de la flamme. Est notée la coïncidence des valeurs calculées et expérimentales.
Авторы: Барсуков Владимир Иванович - кандидат химических наук, доцент кафедры «Физика»; Дмитриев Олег Сергеевич - доктор технических наук, профессор, заведующий кафедрой «Физика»; Барсуков Александр Андреевич - студент, ФГБОУ ВО «ТГТУ», г. Тамбов, Россия.
Рецензент: Попов Николай Сергеевич - доктор технических наук, профессор кафедры «Природопользование и защита окружающей среды», ФГБОУ ВО «ТГТУ», г. Тамбов, Россия.