УДК 539.37:537.221
ВЗАИМОСВЯЗЬ СТРУКТУРЫ РАСТУЩЕГО ПОЛИКРИСТАЛЛА ЛЬДА С ПАРАМЕТРАМИ СПЕКТРА ЭЛЕКТРОМАГНИТНОЙ ЭМИССИИ В ДИАПАЗОНЕ 20-104 Гц
© А.А. Шибков, М.А. Желтое, М.А. Татарко, Е.А. Мишалуев
Shibkov A.A, Zheltov M.A., Tatarko M.A., Mishaluev E.A. Relationship between the polycrystal ice growth and the electromagnetic emission parameters in the frequency range 20-104 Hz. Relationship between the polycrystal structure of the growth ice and the electromagnetic emission spectral parameters in the frequency range 20-104 Hz was investigated. We discovered that the of distilled water freezing is accompanied by electromagnetic wave emission (EME) in 20-104 Hz range. EME is the sequence of electric potential discrete pulses induced on a capacitive probe disposed close to freezing sample. We found good correspondence between time dependence of EME pulse amplitude sum and solid phase fraction volume. Besides there existed quantitative correlation between the total amount of pulses, accumulated in certain moment of time and number of grains of the growing ice polycrystal and also between distribution function of pulse amplitudes and ice grains sizes. We assume that the EME pulses correspond the growth of each grain and that the sequence of pulses displays stepped nature of the polycrystal growth.
Систематическое изучение электрических явлений, сопутствующих фазовым превращениям, особенно кристаллизации прозрачных диэлектриков, началось после открытия в 1945 г. Рибейро [1], а затем Воркманом и Рейнольдсом [2] термодиэлектрического эффекта - возникновения значительной, до сотен вольт, неравновесной межфазной разности потенциалов в системе лед-вода. К настоящему времени обнаружен ряд нестационарных электромагнитных явлений (импульсное радиоизлучение в диапазоне частот 0,1-1,0 МГц, криолюминесценция в инфракрасной и видимой области спектра и др.), сопровождающих замерзание водных растворов и ряда неорганических и органических жидкостей [3 - 7]. В [8] нами зарегистрирована низкочастотная (0,1 -105 Гц) электромагнитная эмиссия (ЭМЭ) при кристаллизации дистиллированной воды. Она имела непрерывную (в виде шума) и дискретную (в виде одиночных импульсов) компоненты и была обусловлена флуктуациями фронта кристаллизации. Таким образом, рост кристаллов льда из расплава (воды) сопровождается весьма сложным широкополосным электромагнитным излучением в диапазоне 0-10ы Гц.
В работах, посвященных изучению электрических явлений при кристаллизации, однако, отсутствует структурный аспект. В то же время хорошо известно, что поликристаллический лед, обычно образующийся в экспериментах такого рода, характеризуется широким спектром структурных элементов различного масштаба, зависящим от термодинамических условий охлаждения. Представляется необходимым сопоставление параметров спектра сигналов ЭМЭ со спектром событий эволюции элементов структуры кристаллизующегося льда, определяемого характерными временами зарождения и роста зерен различной геометрии - игл, дендритов, элементов финговой структуры, их динамического взаимодействия, столкновения, трения и т. д. Поэтому цель настоящей работы заключалась в сопоставлении параметров спектра'ЭМЭ в ранее неисследованном диапазоне
частот 20-104 Гц с особенностями структуры растущего из дистиллированной воды поликристалла льда.
1. МЕТОДИКА ЭКСПЕРИМЕНТА
На предметное стекло (1) микроскопа (2), термостатированного при -10° С, шприцем вводилась проба дистиллированной воды (3) объемом 1 мл в форме диска диаметром 20 мм. Потенциал нестационарного электрического поля измерялся емкостным зондом площадью 2 см2 (4) с отверстием для наблюдения в проходящем свете процесса кристаллизации (рис. 1). Канал регистрации электрического сигнала, состоящий из широкополосного предусилителя (5) (Явх = 5109 Ом, Свх = 20 пФ), аналого-цифрового преобразователя с тактовой частотой 25 кГц (6) и компьютера 1ВМ РС 486 (7), позволял непрерывно в ходе кристаллизации измерять и статистически обрабатывать сигналы ЭМЭ с амплитудой выше 3 мкВ, в полосе частот 20-104 Гц.
Рис. 1. Схема экспериментальной установки: 1 - предметное стекло, 2 - микроскоп, 3 - образец, 4 - емкостный зонд в форме плоского кольца, 5 - предусилитель, 6 - АЦП, 7 - компьютер, 8 - стальной экран, 9 - зеркало, 10 - источник света, 11 - криостат.
Подобная экспериментальная методика отличается от использовавшейся в [4] для регистрации радиоизлучения селективными вольтметрами своей широкопо-лосностью и возможностью анализировать реальную форму электрических импульсов. Кроме того, отсеивание постоянной составляющей сигнала позволяет существенно повысить чувствительность к более тонким электрокинетическим явлениям, по сравнению с электрометрическими измерениями межфазных потенциалов [9]. Наконец, бесконтактность метода избавляет от необходимости учитывать искажения сигнала на низких частотах вследствие процессов перезарядки барьерных емкостей типа элекгрод/расплав, элекгрод/лед из-за разной природы носителей в этих проводниках [10].
Визуальные наблюдения показывают, что кристаллизация воды происходит преимущественно за счет гетерогенного зарождения монокристаллических игл, реже дендритов, которые растут от периферии к центру со скоростью ~10'3 м/с и "сталкиваются" в центральной области, образуя характерную игольчатую структуру поликристаллического льда (рис. 2). Кинетическую картину кристаллизации следует дополнить образованием и движением со скоростями -10'2 м/с многочисленных ступенек роста высотой ~(10-50) мкм на поверхности игл. Учитывая размеры этих элементов структуры (длина иглы (1 -6)-10'3 м, расстояние между крупными ступеньками (2-4)-10"1 м), следует ожидать, что частотный спектр событий, связанных с их эволюцией и взаимодействием друг с другом, находится, по крайней мере, в интервале Ю '-Ю3 Гц.
б)
Рис. 2. Типичная структура поликристаллического льда на начальной (а) и конечной (б) стадиях кристаллизации.
Рис. 3. Частотная зависимость электропроводимости ст (1) и диэлектрической проницаемости е (2) при -10 С [11].
Рис. 4. Частотные зависимости максвелловского времени релаксации тм(1), 1§8 (2), коэффициента передачи К (3), построенные поданным работы [11] и коэффициента передачи *■(4) - по данным калибровки канала регистрации ЭМЭ.
Диэлектрические свойства льда обнаруживают удивительную особенность, заключающуюся в том, что он наиболее прозрачен именно в низкочастотной области спектра. Действительно, как известно, электропроводность льда о в интервале V = 10-104 Гц возрастает с ростом частоты V, а диэлектрическая проницаемость е падает [11] (рис. 3), так что частотная зависимость tg5 (который в низкочастотной области стремится к предельному выражению = а/27СУеео [12]) имеет минимум вблизи 102 Гц, а максвелловское время релаксации тм = о/еео стремится к значению —5* 10‘2 с в области очень низких частот, при V ^ 10 Гц (рис. 4). Если характерное время электрокинетического собы-
тая при кристаллизации Д/ < тм, то лед оказывается “прозрачным” в этом частотном диапазоне и временные изменения потенциала зонда ср(Г) без искажения воспроизводят изменения потенциала электрического поля 11(1), генерируемого электрически активной меж-фазной границей; если Д/ > тм, то у-дШдг, а коэффициент ослабления сигнала вследствие его внутреннего дифференцирования К = ф/£У~тм/Д/. Для построения частотной зависимости К(у) двойной электрический слой у межфазной границы лед-вода моделировался вмороженной в лед дипольной антенной в виде плоского конденсатора, питающегося от генератора сигналов Г6-15. На рис. 4 приведена частотная зависимость коэффициента передачи сигнала К* = ф*/(/*, где ф* - потенциал на зонде, I/ - потенциал на передающей антенне. Видно ее количественное соответствие с кривой К(у), построенной обработкой частотных зависимостей о(у) и е(у), взятых из работы [11].
Прозрачность льда в области низких частот позволяет считать, что спектр сигналов ЭМЭ в этом частотном диапазоне непосредственно несет информацию о нестационарных процессах эволюции электрически активного фронта кристаллизации.
2. РЕЗУЛЬТАТЫ И ОБСУЖДЕНИЕ
По форме и амплитуде все сигналы ЭМЭ, сопровождающие кристаллизацию воды, можно разбить на три характерные группы: дискретные импульсы I типа, треугольной формы с временем фронта Г/ » 10 мс и временем спада т « 50 мс и средней амплитудой 40 мкВ образуют самую многочисленную группу - более 90 % всех зарегистрированных импульсов; импульсы И тана, представляющие собой серию из двух разнополярных импульсов по форме близких к импульсам I типа, наблюдались с вероятностью <10 %; наконец, в редких случаях, особенно на поздних стадиях кристаллизации, возникали короткие, длительностью ~1 мс, почти прямоугольные импульсы различной полярности амплитудой -20-30 мкВ и временем фронта 40 мкс - импульсы Ш типа (рис. 5).
С целью выявления соответствия спектра импульсов ЭМЭ с геометрическими характеристиками игольчатой структуры кристаллизующегося льда, последняя фотографировалась через каждые 30 с, и считывалось число импульсов каждого типа, накопленных к моменту времени Г, которое затем сравнивалось с полным числом игл М,(/). Обнаружено количественное соответствие между /V, и числом импульсов I типа А^, а также между функциями распределения размеров зерен Ь и амплитуд импульсов ф (рис. 6, 7). Можно предположить поэтому, что эволюция каждой иглы сопровождается генерированием импульса I типа. Следует отметить, что длительность этих импульсов значительно меньше характерных времен роста зерен (-3-10 с) или движения на ее поверхности крупных ступенек (-0,1-1 с). Поэтому источниками импульсов I типа могут быть наиболее быстрые события в "биографии" этих объектов, например, зарождение и остановка у препятствия.
Рассмотрим теперь возможные механизмы формирования дискретных импульсов ЭМЭ. Существующие теоретические модели термодиэлектрического эффекта, основанные на анализе реакции гидролиза на фронте кристаллизации [13 - 17], являются динамическими и предсказывают зависимость межфазной разности
б)
в)
Рис. 5. Типичные импульсы ЭМЭ, сопровождающие кристаллизацию дистиллированной воды: а - импульс I типа, б - импульс II типа, в - импульс, построенный интегрированием по времени импульса II типа, г - импульс III типа.
Рис. 6. Временные зависимости числа импульсов I типа N1 (1) и числа зерен льдаМ(2) в ходе кристаллизации.
Рис. 7. Функция распределения амплитуд импульсов I типа (1) и размеров зерен (2) в безразмерных координатах: для импульсов х = ф/ф - безразмерная амплитуда, <р = (40 ± 3) мкВ -
средняя амплитуда импульса; для зерен х = Ь1Ь - безразмерная длина зерна, 1^ = (4,4 ± 0,5) мм - средняя длина зерна.
потенциалов и от скорости роста и твердой фазы типа и~и", где п = 2 [17]. Поэтому эффекты нестабильности фронта кристаллизации, обусловленные, например, стартом и остановкой ледяной иглы или крупной ступенькой роста, должны приводить к временным изменениям 0(1) по форме, близкой к прямоугольному импульсу (рис. 5, в). Отсеивание постоянной составляющей приводит к дифференцированию подобного сигнала, обнаруживаемого в виде пары последовательных разнополярных импульсов с возрастающей амплитудой -сигналы II типа (рис. 5, б). Учитывая характерное время между разнополярными составляющими -100 мс, этот вид сигнала обусловлен, по-видимому, зарождением достаточно крупных ступенек роста и их последующей остановкой, например, в результате зависания на краю иглы. Таким образом, возможной причиной возникновения сигнала ЭМЭ I и II типов является скачкообразное изменение межфазной разности потенциалов вследствие скачков скорости фронта кристаллизации, обусловленных зарождением и блокировкой зерен и крупных ступенек роста. Поэтому функция распределения амплитуд импульсов I типа, вызванных преимущественно столкновением отдельных ледяных игл, отражает распределение их межфазных потенциалов и, следовательно, скоростей в моменты контактов и, в конечном счете, однозначно связана с функцией распределения зерен поликристаплического льда по размерам (рис. 7).
Другой причиной генерирования дискретных сигналов ЭМЭ I и II типов может быть разделение зарядов
в области столкновения игл вследствие трения или, как видно из рис. 2, "микроиндентированием" неподвижной иглы или дендрита растущей иглой за счет, например, псевдопьезоэлектрического эффекта в неоднородном поле упругих напряжений в области контакта [18], движения заряженных дислокаций, частичной режек-ции и микроразрушения. Количественная оценка вклада указанных контактных явлений в формирование дискретных импульсов ЭМЭ требует проведения детального экспериментального моделирования соответствующих ситуаций. Но и в этих случаях амплитуды импульсов ЭМЭ будут пропорциональны скоростям игл в моменты их механического контакта, что связывает функции распределения амплитуд ЭМЭ и размеров зерен.
Наконец, наиболее высокочастотные (в диапазоне (1-10)103 Гц) импульсы ЭМЭ - импульсы III типа - вызваны хрупким разрушением некоторых игл на поздних стадиях кристаллизации вследствие роста дилатацион-ных напряжений: именно на этой стадии наблюдались случаи отламывания и всплывания некоторых игл.
Таким образом, в работе проведено сопоставление спектра дискретной ЭМЭ с кинетическими особенностями эволюции преимущественно игольчатой структуры льда. Показано, что основная часть импульсов находится в диапазоне частот -20-200 Гц и обусловлена динамическим взаимодействием зерен, причем их амплитудно-частотный спектр отражает кинетику формирования структуры зерен поликристаплического льда, а другая - в диапазоне частот 103-104 Гц - обусловлена сопутствующими процессами трещинообра-зования.
ЛИТЕРАТУРА
1. Ribeiro J.C. //An. Acad. Brasil S:. 1945. V. 17. №3. P. 2.
2. Workman E.J., Reynolds S.E. II Phys. Rev. 1950. V. 78. P. 254.
3. Качурин Л.Г.. Бекряев В.И. II ДАН СССР. 1960. Т. 130. С. 57.
4. Качурин Л.Г., Колее С.А., Псаломщиков В.Ф. II ДАН СССР. 1982.
Т. 267. С. 347.
5. Берри Б.Л., Грибов В.А. II Вестн. МГУ. Сер. 5. География 1982.
B. 2.С. 15.
6. Татарченко В.А., Умаров Л.М. II Кристаллография. 1980. Т. 25.
C. 1311.
7. Трохаи А.М., Лапшин А.И., Гудзенко О.И. II ДАН СССР. 1984 Т. 275. С. 83.
8. Головин Ю.И., Шибкое А.А., Желтое М.А. и др. // Изв. РАН. Сер. физическая. 1997. Т. 61. С. 913.
9. Качурин Л.Г., Бекряев В.И., Псаломщиков В.Ф. И ДАН СССР. 1967. Т. 174. С. 1922.
10. Петренко В.Ф., Чеснаков В.А. IIФТТ. 1990. Т. 32. С. 2655.
11. Маэно Н. Наука о льде. М.: Мир, 1988.
12. Кузнецов А.А. Биофизика низкочастотных электромагнитных воздействий. М.: МФТИ, 1994.
13. Чернов А.А., Мельникова A.M. II Кристаллография. 1971. Т. 16. С. 477.
14. Чернов А.А., Мельникова A.M. II Кристаллография. 1971. Т. 16. С. 488.
15. Bronshteyn V.A., Chernov А.А. // J. Cryst. Growth. 1991. V. 112. P. 129.
16. Качурин Л.Г. II Электрохимия. 1970. Т. 6. С. 1294.
17. Качурин Л.Г. И Поверхностные явления в жидкостях. Л.: ЛГУ, 1975. С. 137.
18. Евтушенко А.А., Петренко В.Ф. II ФТТ. 1991. Т. 33. С. 1509.
БЛАГОДАРНОСТИ: Работа выполнена при финансовой поддержке Российского Фонда Фундаментальных Исследований (проект № 98-02-17054).
Поступила в редакцию 25 июля 1998 г.