Кинетические закономерности потери устойчивости и развития дендритных форм при росте кристалла из раствора тема диссертации и автореферата по ВАК РФ 01.04.14, кандидат физико-математических наук Мартюшев, Леонид Михайлович

• Мартюшев, Леонид Михайлович
• кандидат физико-математических наук
• 1998

Оглавление диссертации кандидат физико-математических наук Мартюшев, Леонид Михайлович

СОДЕРЖАНИЕ

ПЕРЕЧЕНЬ УСЛОВНЫХ ОБОЗНАЧЕНИЙ, СИМВОЛОВ,

ЕДИНИЦ И ТЕРМИНОВ

ВВЕДЕНИЕ

ГЛАВА 1. ЗАКОНОМЕРНОСТИ И МЕТОДЫ ИЗУЧЕНИЯ ПОТЕРИ УСТОЙЧИВОСТИ И ДЕНДРИТНОГО РОСТА

1.1. Механизмы потери устойчивости

1.2. Аналитические подходы к проблеме потери устойчивости растущего кристалла

1.3. Методы компьютерного моделирования потери устойчивости и развития дендритной структуры

ГЛАВА 2. ПРОИЗВОДСТВО ЭНТРОПИИ ПРИ УСТОЙЧИВОМ И НЕУСТОЙЧИВОМ РОСТЕ КРИСТАЛЛА ИЗ РАСТВОРА

2.1. Вывод выражения для производства энтропии

2.2. Производство энтропии при росте шара из раствора

2.3. Производство энтропии при росте бесконечного цилиндра

из раствора

2.4. Производство энтропии при росте бесконечной плоскости

из раствора

2.5. Общий принцип развития кристаллизации при

неравновесном росте

Выводы к главе 2

ГЛАВА 3. КОМПЬЮТЕРНЫЙ АНАЛИЗ ОБРАЗОВАНИЯ ДЕНДРИТОВ ПРИ РОСТЕ КРИСТАЛЛА ИЗ РАСТВОРА

3.1. Построение компьютерной модели

3.2. Влияние кинетического коэффициента кристаллизации

на потерю устойчивости и развитие дендритных форм

3.3. Кинетические особенности растущих из раствора кристаллов. Поведение скорости роста и производства энтропии со временем

3.4. Экспериментальная проверка временных осцилляций производства энтропии при росте дендрита

3.4.1. Описание экспериментальной установки

3.4.2. Результаты эксперимента и обсуждение

Выводы к главе 3

ГЛАВА 4. КОМПЬЮТЕРНЫЙ АНАЛИЗ ОБРАЗОВАНИЯ ДЕНДРИТОВ ПРИ РОСТЕ КРИСТАЛЛА ИЗ РАСТВОРА В УСЛОВИЯХ ФАЗОВОГО РАССЛАЕНИЯ

4.1. Модификация компьютерной модели

4.2. Влияние относительного содержания фазово-расслаива-ющейся примеси на морфологию растущего кристалла

4.2.1. Квазиустойчивый гранный рост в фазово-расслаивающейся среде

4.2.2. Рост скелетного кристалла в фазово-расслаивающейся среде

4.2.3. Рост дендрита в фазово-расслаивающейся среде

4.3. Выводы по действию фазово-расслаивающейся примеси на морфологию растущего кристалла

4.4. Автомодельное решение в пределе кинетического режима роста кристалла с фазово- расслаивающейся примесью

ЗАКЛЮЧЕНИЕ

СПИСОК ИСПОЛЬЗОВАННЫХ ИСТОЧНИКОВ

ПЕРЕЧЕНЬ УСЛОВНЫХ ОБОЗНАЧЕНИЙ СИМВОЛОВ, ЕДИНИЦ

И ТЕРМИНОВ

1г площадь поверхности кристалла; 5 - энтропия;

N - число частиц в твердой фазе; Р - производство энтропии; / - энтальпия;

1лт - химический потенциал соли в растворе; /и - химический потенциал соли в кристалле; У1т - гармоника сферической функции; Кк - функция Ханкеля;

Я* - радиус критического зародыша образования кристаллической фазы; Яс - критический радиус потери устойчивости;

Кь- радиус потери устойчивости, найденный термодинамическим методом;

у - свободная поверхностная энергия; Г - коэффициент поверхностного натяжения; Т - абсолютная температура; В - коэффициент диффузии;

сф - первоначальная концентрация соли в растворе;

С - постоянная концентрация растворенного вещества в кристалле;

с - зависящая от координат концентрация соли в растворе;

с0 - равновесная концентрация соли у плоской поверхности раздела;

с5 - равновесная концентрация соли у поверхности общего типа;

сК - концентрация на поверхности неискаженной сферы;

¿о - капиллярная длина;

Ъ - кинетический коэффициент кристаллизации; Ь- размер квадратной ячейки моделирования; х- временной шаг моделирования; М5 приближение - приближение Маллинза-Секерки.

Введение диссертации (часть автореферата) на тему «Кинетические закономерности потери устойчивости и развития дендритных форм при росте кристалла из раствора»

ВВЕДЕНИЕ

Теоретическим и экспериментальным вопросам образования сложных пространственно-временных структур, возникающих при неравновесных процессах, уделяется в настоящее время очень большое внимание. В литературе исследуется много различных систем (гидродинамических, химических, биологических и других), у которых при некоторых критических значениях термодинамических потоков и сил, происходит скачкообразное изменение свойств и симметрии (так называемые неравновесные фазовые переходы) [1-3]. Одной из интересных и практически важных для изучения систем является неравновесный рост кристаллов из раствора/расплава, в результате которого при определенных размерах кристалла и степени пересыщения/переохлаждения возникают скелетные и дендритные (древовидные) формы [4]. Данные неравновесные переходы являются типичными представителями неравновесных фазовых переходов, изучаемых синергетикой (неустойчивость Бенара, возникновение турбулентности и другие) [1,2] и характеризуются разнообразием структур, возникающих в результате развития неустойчивости при росте кристалла из различных растворов/расплавов. Особенностью этого класса неравновесных фазовых переходов может служить то, что при росте кристалла возникшая пространственно-временная структура как бы застывает ("фиксируется") в пространстве. Другой их характерной особенностью является наличие двух одновременных переходов: обычного термодинамического перехода (раствор - кристалл) и неравновесного фазового перехода.

Важность изучения данных явлений определяется двумя причинами. Первая- это практическое значение. Известно, что дендритно затвердевает

большинство расплавов при охлаждении, что существенно и не всегда положительно влияет на физические и механические свойства слитка [5,6]. Изучение морфологии (внешней формы) образовавшейся структуры также может много рассказать о тех неравновесных условиях при которых происходил рост, что практически важно для таких разделов как геофизика [7], астрофизика [8], медицина [9,10], метеорология [11]. Вторая причина -теоретического характера. Образование дендрита и развитие его ветвей является достаточно простым (в сравнении, например, с турбулентностью) и хорошо экспериментально изученным процессом и поэтому может служить и часто служит основой для развития и разработки теоретических подходов, математических методов, выявления общих закономерностей, которые затем могут быть использованы при изучении других самоорганизующихся в неравновесных условиях систем.

Несмотря на то, что дендритный рост интенсивно изучается уже, по крайней мере, последние двадцать лет (см. например обзоры [4,12-15]), многие вопросы либо не рассмотрены до конца, либо не рассматривались вообще.

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

фундаментальные свойства явления не столь сразу заметные в более "подробных и искушенных" методах.

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

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

Здесь мы отметили лишь некоторые наиболее существенные с нашей точки зрения вопросы, разрешение которых позволит более глубоко понять особенности неравновесной кристаллизации.

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

Научная новизна:

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

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

- С помощью т-модели проанализированы потеря устойчивости и развитие дендритных форм при росте из раствора при изменении кинетического коэффициента кристаллизации. Найдены зависимость скорости роста кристалла и производства энтропии от времени.

- Экспериментально, с помощью поляризационно- интерференционной микроскопии подтверждены полученные с помощью т-модели результаты о наличии периодических осцилляций массы растущего дендрита со временем.

- С помощью т-модели впервые проанализирован квазиустойчивый гранный рост, скелетный рост, и рост дендрита в растворе с фазово-расслаивающейся примесью.

На защиту выносятся:

- полученное выражение для производства энтропии;

- расчеты по изменению производства энтропии при росте и потере устойчивости шара, бесконечного цилиндра и бесконечной плоскости из пересыщенного раствора;

- предложенный принцип отбора режима течения процесса при неравновесном росте кристалла из раствора;

- Результаты расчетов и опытов по обнаружению периодических осцилляций скорости роста первичной ветви и массы растущего дендрита, частота которых связана с частотой появления вторичных ветвей, а также предложенный механизм возникновения данных осцилляций;

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

Заключение диссертации по теме «Теплофизика и теоретическая теплотехника», Мартюшев, Леонид Михайлович

4.3. Выводы по действию фазово-расслаивающейся примеси на морфологию растущего кристалл

Обсудим теперь общую схему кинетического действия фазово-расслаивающейся примеси, используя результаты, полученные в разделе 4.2.

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

4.4. Автомодельное решение в пределе кинетического режима роста кристалла с фазово- расслаивающейся примесью

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

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

Данная постановка интересна по следующей причине. Пусть мы хотим изучить морфологию растущего кристалла, для этого необходимо разбить весь объем изучаемой области на маленькие ячейки, размер которых выбирается из соображений необходимой степени подробности в описании. Внутри каждой такой ячейки характеристики кристалла будут однородными. Далее в рамках приведенной выше модели, мы рассчитываем кристалл. В результате мы получаем некоторый pattern/ кристалл с масштабом неоднородности в выбранную нами ячейку. Используемые в расчетах правила при выбранных ограничениях не зависят от выбираемого масштаба. Действительно, если рост происходит при кинетическом режиме, то концентрация у поверхности кристалла во всех точках одинакова и равна концентрации в растворе [19]. Динамика роста однозначно определяется особенностями роста конкретного кристалла (граничными законами / уравнениями), так как мы остаемся в рамках макроскопического описания, то в широком диапазоне масштаба они одни и те же. На диффузионное поле примеси также не сказывается масштаб описания, так как ее коэффициент диффузии равен нулю и распределение ее в объеме одно и тоже и равно начальному, а оттеснение примеси происходит по одним и тем же законам благодаря масштабной инвариантности динамики роста кристаллической поверхности. Допустим, что теперь нам необходимо рассчитать рост того же кристалла с меньшим, либо большим характерным размером (более ранние или более поздний стадии роста, соответственно). Снова разбиваем уже новый характерный размер на ячейки, причем число ячеек оставляем равным предыдущему (т.е. мы хотим исследовать данный характерный размер кристалла с той же степенью подробности, что и в предыдущем случае). Но так как правила расчета кристалла при этом не меняются , то в результате расчета мы, очевидно, получим структуру, абсолютно подобную предыдущей. Изменяется лишь масштаб неоднородностей, а морфология кристалла в целом остается без изменений.

Таким образом, для большого диапазона пространственного масштаба, начиная с достаточно малого (нескольких радиусов критического зародыша) до относительно больших размеров устойчивых многогранников, растущих из раствора, в данной постановке задачи структура оказывается самоподобна. Другими словами развитие кристаллической структуры во времени при рассматриваемых условиях роста является автомодельным [108]. Следовательно, изучив один из масштабов структуры, мы все знаем и о морфологии в других масштабах (при другом увеличении).

Автомодельные режимы при росте кристалла из раствора при определенных условиях были известны и ранее. Классической задачей об автомодельном росте является задача о росте сферического кристалла в диффузионно-лимитируемых условиях (задача Бспуеп-ЮгкаЫу) [109], решенная аналитически. Также автомодельность возникает в моделях типа БЬА [3,74]. В настоящей работе автомодельность возникает при кинетическом режиме кристаллизации в присутствии примеси с нулевым коэффициентом диффузии. В отличии от тривиальной автомодельности при кинетическом росте без примеси, данный случай морфологически будет намного интереснее. Еще раз отметим, что как любое автомодельное решение оно не является лишь частным и изящным решением нелинейной задачи при определенных условиях, а описывает промежуточное асимптотическое поведение широкого класса других более реальных задач кристаллизации. Нахождение и анализ данного автомодельного решения произведем методом компьютерного моделирования.

Для моделирования условий роста, описанного выше, используем модель, построенную на принципах клеточных автоматов. В модели используются следующие правила:

1. Используется квадратная двухмерная решетка, ЫхЫ клеток. Выбор квадратной решетки предполагает возможность моделирования роста кристалла кубической симметрии. Размер клетки в некотором смысле л произволен, сверху он ограничен 1/Ы частью критического размера устойчивости образующегося кристалла, а снизу критическим зародышем [19]. Выбор квадратной и двухмерной решетки сделан в настоящей работе для определенности, по аналогии модель можно построить для трехмерной решетки и другой симметрии. Тестирование программы показало, что выбирать N>100 не имеет смысла, так как "рисунок" кристалла повторяется.

2. Предполагается, что в центре решетки начинает расти кристалл. Так как рассматривается предельный случай - кинетический режим роста [19], концентрации у вершин и в середине граней кристалла одинаковы и равны концентрации в объеме, то для моделирования такого роста можно выбрать предельно простое правило: кристалл за каждый временной шаг последовательно увеличивается, что выражается в закрашивании черным всех ближайших к кристаллу клеток за один цикл (см. рис.4.6), если это не запрещено свойствами второго компонента (см. ниже). Характер рассматриваемого роста диктует необходимость пересматривать за каждый временной шаг соседей всех закрашенных клеток.

3. Правила, моделирующие действие фазово- расслаивающейся примеси, полностью аналогичны изложенным в разделе 4.1. Поэтому здесь не приводятся.

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

Построенная выше модель для анализа кинетического режима роста кристалла кубической симметрии при влиянии фазово-расслаивающейся примеси имеет лишь один управляющий параметр -относительное содержание примеси Сбнач/Сбтв. При Сбнач/Сбтв->0 кристалл является идеально кубическим.

Область раствора вблизи кристаллического зародыша

Л л * л «ч ч > V а я У

V

Кристаллический зародыш (закрашен черным). Каждая клетка имеет четыре ближайших соседа. Стрелки показывают порядок пересчета клеток в течении одной итерации. Клетка закрашивается черным, если это не запрещено свойствами второго компонента и если ее ближайший сосед является черным.

Рис.4.6

Проанализируем поведение системы при постепенном увеличении содержания примеси. При Сбнач/С6тв<0.05 кристалл, оставаясь в целом однородным, в процессе роста начинает захватывать маленькие области фазово-расслоившейся примеси (рис.4.7(а). Эти включения, благодаря вытеснению, скапливаются в центральной части граней растущего кристалл. С увеличением концентрации примесь начинает осаждаться все на более ранних стадиях роста и распределение ее приобретает ярко

Морфология кристалла при изменении относительной концентрации примеси

СШСт ■ °-05 (а); 0.10 (б); 0.22 (в); 0.34 (г); 0.36 (д); 0.41 (е|; 0.55 (ж). Черным цветом показано распределение соли, серым распределение фазово-расслаившейся примеси.

Рис.4.7. выраженный слоевой характер (рис.4.7(б)). При дальнейшем величении Сбнач/Сбтв расстояние между слоями, образованными фазово-расслоившейся примесью, уменьшается, а ее плотность в этих слоях увеличивается (рис.4.7(в)). Одновременно однородная область без примеси становится все более тонкой. Данная область (остов переменной толщины, образованный двумя взаимно перпендикулярными прямыми, проходящими через центр и вершины кристалла) делит растущий кристалл на четыре четверти. Благодаря симметрии при рассмотрении внутренней структуры кристалла достаточно рассмотреть распределение кристаллизующегося и примесного компонента лишь в одной его четверти (далее в работе будет идти речь лишь об одной четверти). Очевидно, что физические свойства всего кристалла должны зависеть не только от количества примеси, но и от направлений ее основного распределения. В частности, одной из основных характеристик, определяющей физические свойства, является направление <110>, проходящее через центр кристалла и делящее рассматриваемую четверть кристалла пополам. По данной прямой осуществляется связь между двумя зеркально симметричными 1/8 частями кристалла. Введем характеристику связности этих двух областей £ как число клеток кристалла, расположенных на этой прямой. Параметр связности, очевидно, в данной модели может меняться от 35 (л/532 + 502/2), что соответствует росту в пределах рассматриваемой сетки 100x100 в отсутствие примеси до 1 ( примесь подавляет рост на самой ранней стадии). Для уже рассмотренных случаев рис.4.7(а-в) имея относительно большие значения, постепенно уменьшается.

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

При незначительном изменении Сбнач/Сбтв до 0.36 внутренняя структура кристалла резко изменяется (рис. 4.7(д)). Фазово-расслоившаяся примесь преимущественно начинает располагаться под углом 90 градусов к граням, соответственно, свойства кристалла приобретают ярко выраженную анизотропность. Параметр связности уменьшается до двух. Налицо ярко выраженная морфологическая перестройка. Последующее увеличение относительной концентрации снова приводит к однородному распределению примеси и изотропным свойствам кристалла ( по рассматриваемой четверти). Параметр связности, резко увеличившись до 9 (при Сбнач/Сбтв=0.38), начинает снова медленно уменьшаться. На рис.4.7(е) представлена типичная структура для этой области изменения управляющего параметра (£=7). В интервале изменения Сбнач/Сбтв от 0.41 до 0.55 характер распределения примеси остается прежним, однако внутренняя структура приобретает все более регулярный вид рис.4.7(ж). Параметр связности уменьшается до одного, а пористость кристалла увеличивается. Граница роста превращается из относительно гладкой в фракталоподобную. При дальнейшем увеличении относительной концентрации примеси (Сбнач/Сбтв>0.56) происходит образование недоразвитых кристаллических структур.

Итак, в настоящем разделе рассмотрен кинетический режим двухмерного роста кристалла кубической симметрии при влиянии фазово-расслаивающейся примеси с коэффициентом диффузии, равным нулю. В таком приближении рост кристалла является автомодельным. Показаны последовательные этапы модификации структуры при увеличении относительной концентрации примеси. Получено, что в результате таких перестроек характер распределения примеси меняется от слоевого до практически однородного. Также обнаружено, что в результате структурных перестроек возможно возникновение анизотропности в распределении примеси в рассматриваемой кристаллической четверти (направление <110>). Полученные автомодельные решения являются асимптотическими для целых классов задач, например задач о кинетическом режиме роста в присутствии примеси с относительно малым коэффициентом диффузии и взаимодействующей с поверхностью ( адсорбция, хемосорбция и т.п.). Автомодельные решения (рис.4.7) можно использовать как эталоны при приближенных расчетах и компьютерном моделировании более сложных задач, они хорошо передают геометрический аспект действия примеси.

ЗАКЛЮЧЕНИЕ

Основными результатами работы являются.

1. Методом нелокальной неравновесной термодинамики получено выражение для производства энтропии при изотермо-изобарической кристаллизации и проведено его критическое сравнение с имеющимися.На основе обобщения имеющихся экстремальных принципов неравновесного роста, предложена гипотеза о направлении кристаллизации при неравновесном росте - принцип максимально возможного значения производства энтропии при росте.

2. С помощью полученного выражения для производства энтропии на примере трех классических задач о потере устойчивости при кристаллизации (рост шара, бесконечного цилиндра и бесконечной плоскости) рассмотрены следствия, к которым приводит предложенная гипотеза. Получено, что в случае роста шара и цилиндра разница между производством энтропии в возмущенном и невозмущенном случае равна нулю в точке отличной от точки потери устойчивости, предсказываемой классической теорией Маллинза -Секерки. В случае роста плоскости обе эти точки совпали. Также получено, что в точке перехода производство энтропии для шара и цилиндра скачкообразно увеличивается.

3. В работе построена компьютерная т-модель, сочетающая в себе достоинства моделей, построенных на принципах клеточных автоматов и численного моделирования. С помощью построенной модели исследована потеря устойчивости и дендритный рост при изменении кинетического коэффициента кристаллизации.

4. С помощью компьютерного моделирования показано, что скорость роста первичной ветви дендрита, растущего из раствора, испытывает малоамплитудные осцилляции, которые связаны с появлением вторичных ветвей. Показано также, что зависимость производства энтропии от времени имеет осциллирующий характер при росте дендрита.

5. Экспериментально на примере кристаллизации М^СЛ из водного раствора с использованием поляризационно-интерференционной микроскопии исследована зависимость массы дендритов (и производства энтропии) от времени роста. Обнаружены осцилляции массы, имеющие неслучайный, периодический характер. Появление данных осцилляций связывается со взаимодействием диффузионных полей существующих и вновь возникающих вторичных ветвей.

6. С помощью компьютерного моделирования исследовано влияние фазово-расслаивающейся примеси на рост кристалла. Обнаружено, что при квазиустойчивом гранном росте действие такой примеси приводит к постепенной модификации структуры, при скелетном росте возможны как постепенная модификация структуры, так и скачкообразные изменения- неравновесный фазовый переход. При росте дендрита в фазово-расслаивающейся среде увеличение концентрации примеси может привести к циклической последовательности неравновесных фазовых переходов- возвратных фазовых переходов.

7. Рассмотрены условия, при которых рост кристалла с примесью можно рассматривать как автомодельный. С помощью компьютерной модели произведен анализ данного случая.

По результатам проведенного исследования опубликовано 10 работ [110-119].

Список литературы диссертационного исследования кандидат физико-математических наук Мартюшев, Леонид Михайлович, 1998 год

СПИСОК ИСПОЛЬЗОВАННЫХ ИСТОЧНИКОВ.

1. Хакен Г. Синергетика. М.:Мир,1980. 404 с.

2. Николис Г., Пригожин И. Самоорганизация в неравновесных системах. От диссипативных структур к упорядоченности через флуктуации. М.: Мир, 1979. 308 с.

3. Федер Е. Фракталы. М.: Мир, 1991. 260 с.

4. Langer J.S. Instabilities and pattern formation in crystal growth // Rev.Mod.Phys. 1980. v.52. N 1. P. 1-28.

5. Чалмерс Б. Теория затвердевания. M.: Металлугрия, 1968. 288 с.

6. Флеминге М. Процессы затвердевания. М.: Мир, 1977. 423 с.

7. Jamtveit В., Andersen Т. Morphological instabilities during rapid growth of metamorphic garnets // Physics and Chemistry of Minerals. 1992. v. 19,

P. 176-184.

8. Додд P.T. Метеориты. Петрология и геохимия. M.: Мир, 1986. 384 с.

9. Sole A. Stagoskopie. Wien: Franz Deuticke, 1960. 92 с Ю.Кристаллографические методы исследования в медицине. Сб. науч. тр.

1-й Всероссийской научно-практической конференции. Москва. 1997.

11.Mason B.J. Snow crystals, natural and man made // Contemporary Physics. 1992. v. 33. N4. P. 227-243.

12. Kessler D., Koplik J., Levine H. Pattern selection in fingered growth phenomena // Adv. Phys. 1988. v.37. N 37. P.255-259.

13. Brener E.A., Melnikov V.I. Pattern selection in two-dimensional dendritic growth // Adv. Phys. 1991. v. 40. N1. P. S3-97.

14. Langer J.S. Issues and opportunities in materials research // Physics Today. 1992. V.45. N 10. P.24-34.

15. Ben-Jacob E., Garik P. The formation of pattern in non-equilibrium growth // Nature. 1990. V.343. P.523-530.

16. Тоффоли Т., Маргулис H. Машины клеточных автоматов. M.: Мир, 1991.280 с.

17. Джексон К. Основные представления о росте кристалла// Проблемы роста кристаллов. Под ред. Н.Н.Шефталя и Е.И.Гиваргизова М.: Мир, 1968. С. 13-26.

18. Мюллер-Крумбхаар X. Моделирование роста кристаллов методом Монте-Карло // Методы Монте-Карло в статистической физике. М.: Мир, 1982. С. 267-325.

19. Современная кристаллография т.З. Образование кристаллов/ Чернов А.А., Гиваргизов Е.И., Багдасаров Х.С. и др. М.:Наука, 1980. 407 с.

19. Чернов А.А. Теория устойчивости гранных форм роста // Кристаллография. 1971. т. 16. вып.4. С. 842-863.

21. Nanev С., Iwanov D. Alteration of the growth form of zinc single crystals as a result of the diffusion non-homogeneity of the supersaturation // J. Crystal. Growth. 1968. v.3.N.4. P.530-534.

22. Nenov D., Stoyanova V. On the formation of ice dendrites from the vapour phase// J. Crystal. Growth. 1977. vr 41. N1. P.73-76.

23. Саратовкин Д.Д. Дендритная кристаллизация. М.: Металлургиздат, 1957.127 с.

24. Бакли Г. Рост кристаллов. М.: ИЛ, 1954. 403 с.

25. Измайлова В.Н., Ребиндер П.А. Структурообразование в белковых системах. М.: Наука, 1974. 268 с.

26. Браун Г.Г., Уолкен Д.Д. Жидкие кристаллы и биологические структуры. М.: Мир, 1982. 198 с.

27. Минц Р.И., Скопинов С.А., Кадушников P.M. и др. Диффузно-ограниченная агрегация в тонких пленках водно-солевых растворов белка на твердотельной подложке. Эксперимент и компьютерная модель // Ж. Физической химии. 1992. т.66. N2. С. 352-355.

28. Минц Р.И., Берг Д.Б. Экспериментальное и компьютерное исследование агрегации молекул в системе кристалл-жидкий кристалл// Ж. Физической химии. 1995. т.69. N1. С. 48-51.

29. Mintz R.I.,.Skopinov S.A, Yakovleva S.V. et al. Formation of oriented films by mucus glycoproteins on the faces of NaCl crystal// Studia biophysics. 1989. v.133. N3. P.221-225.

30. Yasui M., Matsushita M. Morphological changes in dendritic crystal growth of ammonium chloride on agar plate // Journal of Physical Society of Japan. 1992. v.61.N7. P.2327-2332.

31. Минц Р.И., Скопинов C.A., Яковлева C.B. Фракталы в лиотропных системах // Письма в ЖТФ. 1988. т.14. вып.23. С.2204-2207.

32. Измайлова В.Н., Ямпольская Г.П., Сумм Б.Д. Поверхностные явления в белковых системах. М.: Химия, 1988. 238 с.

33. Шинода К. Коллоидные поверхностно-активные вещества. Физико-химические свойства. М.: Мир, 1966. 319 с.

34. Лодиз Р., Паркер Р. Рост монокристаллов. М.: Мир, 1974. 540 с.

35. Гафийчук В.В. Динамика формирования поверхностных структур в системах со свободной границей. Киев: Наукова Думка, 1990. 216 с.

36. Хорвай Г. Модифицированная задача Стефана // Инженерно-физический журнал. 1965. т. 8. N 6. С. 779-800.

37. Chandrasekhar S. Hydrodynamic and hydromagnetic stability, Oxford: Clarendon Press, 1961. 320 p.

38. Гершуни Г.З., Жуковицкий E. M. Конвективная устойчивость несжимаемых жидкостей. М.: Наука, 1972. 392 с.

39. Mullins W.W., Sekerka R.F. Morphological stability of a particle when growth is controlled by diffusion or heat flow // J. Appl.Phys. 1963. v.34, P.323-340.

40. Маллинз В., Секерка P. Морфологическая устойчивость частицы, растущей за счет диффузии или теплоотвода// Проблемы роста кристаллов. М.: Мир, 1968. с. 88-101.

41. Coriell S.R., Parker R.L. Stability of the shape of a solid cylinder growthing in a diffusion field // J. Appl.Phys. 1965. v.36, N 2. P.632-637.

42. Chou H., Cummins H. Z. Evolution of the dendritic instability in solidifying succinonitrile // Phys. Rev. Lett. 1988. v. 61. N 2. P.173-176.

43. Ben-Jacob E., Garik P., Mueller T., Grier D. Characterization of morphology transitions in diffusion-controlled systems// Phys. Rev. A. 1988. v. 36. N 3. P. 1370-1380.

44. Shochet O., Kassner K., Ben-Jacob E. et al. Morphology transitions during non-equilibrium growth. II. Morphology diagram and characterization of the transition // Physica A. 1992. v. 187, P.87-111.

45. Hutter J.L.,.Bechhoefer J., Many modes of rapid solidification in a liquid crystal // Physica A. 1997. v. 239. P. 103-110.

46. Hutter J.L., Bechhoefer J. Three classes of morphology transitions in the solidification of a liquid crystal // Phys.Rev. Lett. 1997. v. 79. N 20. P. 4022-4025.

47. Mu Wang, Nai-ben Ming. Alternating morphology transitions in electrochemical deposition// Phys.Rev. Lett. 1993. V. 71. N 1. P. 113-116.

48. Овруцкий A.M. Расчеты профиля граней послойно растущих кристаллов// Кристаллография. 1978. т. 23. вып. 5, С. 925-930.

49. Овруцкий A.M. Моделирование свободного роста кристалла на стадии развития дендритной формы // ЖЭТФ. 1991 т. 99. вып. 1. С.250-259.

бО.Уманцев А.Р., Виноградов В.В., Борисов В.Т. Математическое моделирование роста дендритов в переохлажденном расплаве// Кристаллография. 1985. т. 30. вып. 3. С.455-460.

51. Sethian J.A., Strain J. Crystal growth and dendritic solidification// J. Computational Physics. 1992 v. 98. P.231-253.

52. Hunt J.D. A numerical analysis of time dependent isolated dendritic growth for conditions near the steady state// Acta metall. mater. 1990. v. 38. N3. P. 411-418.

53. Hunt J.D. A numerical analysis of dendritic and cellular growth of a pure material investigating the transition from "array" and "isolated" growth// Acta metall. mater. 1991. v. 39, N9, P.2117-2133,

54.Уманцев A.P., Виноградов B.B., Борисов B.T. Моделирование эволюции дендритной структуры// Кристаллография. 1986. т. 31. вып.5. С. 1002-1008.

55. Тарабаев JI. П., Машихин А.Ю., Есин В.О. Дендритный рост в переохлажденном расплаве// Расплавы. 1991. N 2. С. 89-100.

56. Галенко П.К., Журавлев В.А. Физика дендритов. М.: Софт-Москва, 1993. 181 с.

57. Galenko Р.К., Krivilyov M.D., Buzilov S.V. Bifurcations in a sidebranch surface of a free-growing dendrite// Phys.Rev.E. 1997. V.55. N1. P.611-619.

58. Проблемы кристаллизации сплавов и компьютерное моделирование. Тезисы Всесоюзной научно-технической конференции: Проблемы кристаллизации сплавов и компьютерное моделирование.ТЗжевск: Из-во Ижевский университет, 1990.154 с.

59. Кристаллизация и компьютерные модели. Труды пятой научно-технической конференции. Ижевск, 6-7 октября 1992 г. Ижевск: Из-во Ижевский университет, 1994. 171 с.

60. Shu-Zu Lu, Hunt J.D., Gilgien P., Kurz W. Cellular and dendritic growth in rapidly solidified Al-Fe and Al-Cu alloys// Acta metall. mater. 1994. v. 42. N5. P.l653-1660.

61. Saito Y., Ueta T. Monte Carlo studies of equilibrium and growth shapes of a crystal// Phys. Rev. A. 1989. v. 40. N6. P.3408-3419.

62. Saito Y., Ueta T. Lattice-gas madel simulation of crystal shapes// J. Crystal. Growth. 1990. v. 99. P. 171-174.

63. Xiao R.F., Alexander J. I. D., Rosenberger F. Morphological evolution of growing crystals: A Monte Carlo simalation// Phys.Rev.A. 1988. v.38. N5. P.2447-2456.

64. Xiao R.F., Alexander J. I. D., Rosenberger F. Growth morphology with anisotropic surface kinetics// J. Crastal. Growth. 1990. v.100. P.313-329.

65. Xiao R.F., Alexander J. I. D., Rosenberger F. Growth morphologies of crystal surface// Phys. Rev. A. 1991. v. 43 .N6. P.2977-2992.

66. Witten T.A., Sander L.M. Diffusion-limited aggregation, a kinetic critical phenomenon// Phys. Rev Lett. 1981. v.47, N 19, P. 1400-1403.

67. Witten T.A., Sander L.M. Diffusion-limited aggregation// Phys. Rev. B. 1983. v. 27. N9. P.5686-5697.

68. Вичек Т. Формирование структур отвердевания в моделях агрегации // Фракталы в физике, под ред. Л.Пьетронеро и Э.Тозати. М.: Мир, 1988. С. 345-349.

69. Meakin P. Universality, nonuiversality, and the effects of anisotropy on diffusion-limited aggregation//Phys.Rev.A. 1986, v.33. N5. P. 3371-3382.

70. Микин П. Некоторые последние достижения в моделировании ограниченной диффузии агрегации и родственных процессов// Фракталы в физике, под ред. Л.Пьетронеро и Э.Тозати. М.:Мир, 1988. С. 283-295.

71. Uwaha М., Saito Y. Fractal-to-compact transition and velocity selection in aggregation from lattice gas// J. Phys. Soc. Japan. 1988. v.57. N 10. P.3285-3288.

72. Uwaha M., Saito Y. Aggregation growth in a gas of finite density: velocity selection via fractal dimension of diffusion-limited aggregation// Phys. Rev. A. 1989. v.40.N8. P.4716-4723.

73. Lebovka N.I., Ivanenko Ya. V., Vygornitskii N.V. Deterministic Eden model of charged-particles aggregation// Europhys. Lett. 1998. v.41, N 1, P. 19-24.

74. Смирнов Б.М. Физика фрактальных кластеров. М.:Наука, 1991. 136 с.

75. Brown S.G.R., Williams Т., Spittle J.A. A cellular automaton model of the steady-state "free" growth of a non-isotermal dendrite// Acta metall. mater. 1994. v. 42. N 8. P.2893-2898.

76. Spittle J.A., Brown S.G.R. A cellular automaton model of steady-state columnar-dendritic growth in binar alloys// J. Materials Science. 1995. v.30. P.3989-3994.

77. Gandin Ch.-A., Rappaz M. A coupled finite element-cellular automaton model for the prediction of dendritic grain structures in solidification processes// Acta metall. mater. 1994. v.42. N 7. P.2233-2246.

78. Gandin Ch.-A., Schaefer R.J., Rappaz M. Analytical and numerical predictions of dendritic grain envelopes// Acta mater. 1996. v.44. N 8. P. 3339-3347.

79. Берг Д.Б., Кадушников P.M., Скопинов C.A. и др. Компьютерная модель процессов диффузно-ограниченной агрегации в лиотропных системах// Изв. АН СССР. Серия физическая. 1991. Т. 55. N9. с. 18531856.

80. McCarthy J.F. Lattice gas cellular automata method for flow in the interdendritic region// Acta metall. mater. 1994. v.42. N 5. P.1573-1581.

81. Kadanof L.P. On two levels// Physics Today. September 1986. P.7-9.

82. Theory and Application of Cellular Automata/ eds. S. Wolfram. Singapore: World Scientific, 1986. 386 p.

83. Лоскутов А.Ю., Михайлов A.C. Введение в синергетику. М.: Наука, 1990.270 с.

84. RothmanD.H., Zaleski S. Lattice-gas models of phase separation: interfaces, phase transitions, and multiphase flow// Rev. Mod. Phys. 1994. v.66. N 4. P.1417-1480.

85. Беркович С .Я. Клеточные автоматы как модель реальности: поиски новых представлений физических и информационных процессов. М.: Изд.-во МГУ, 1993. 112 с.

86. Toffoli Т. Cellular automata as an alternative to ( rather than an approximation of) differential equation in modeling phycics// Physica D 1984. v.lO.P.l 17-127.

87. Трейвус Е.Б. Кинетика роста и растворения кристаллов. JL: Изд-во Ленин, ун-та, 1979. 248 с.

88. Трейвус Е.Б. Введение в термодинамику кристаллогенезиса: Учеб. пособие. Л.: Лен. ун-т, 1990. 150 с.

89. Гленсдорф П., Пригожин И. Термодинамическая теория структуры, устойчивости и флуктуаций. М. Мир, 1973. 280 с.

90.Мелких А. В., Селезнев В. Д. К вопросу об изменении производства энтропии при кинетических фазовых переходах первого рода// Метастабильные состояния и фазовые переходы. Екатеринбург, 1997. С.188-195.

91.Климонтович Ю.Л. Турбулентное движение и структура хаоса. М.: Наука, 1990. 235 с.

92. Климонтович Ю.Л. Статистическая теория открытых систем. М.: ТОО "Янус" 1995. 624 с.

93. Эбелинг В. Образование структур при необратимых процессах. М.: Мир, 1979. 279 с.

94. Кудрявцев И.К. Химические нестабильности. М.: Изд-во МГУ, 1987. 252 с.

95. Weinberg F.,Chalmers В. Further observations on dendritic growth in metals// Canad. J. Phys. 1952.v.30. N 5. P.488-502.

96. Raz E., Lipson S.G., Polturak E. Dendritic growth of ammonium chloride crystals: Measurements of the concentration field and a proposed nucleation model for growth// Phys. Rev. A. 1989. v. 40. N2. P.1088-1095.

97.Huang S.-C., Glicksman M.E. Fundamentals of dendritic solidification -I Steady-state tip growth// Acta Metallurgica. 1981. v.29. P.701-715.

98.Huang S.-C., Glicksman M.E. Fundamentals of dendritic solidification -II Development of sidebranch structure// Acta Metallurgica 1981. v.29. P.717-734.

99.Chan S.K., Reimer H.H., Kahlweit M. On the stationary growth shape of NH4CI dendrites// J. Cryst. Growth. 1976. v. 32. P. 303-315.

100. Dougherty A., Gollub J.P. Sready-state growth of NH4CI from solution// Phys. Rev. A. 1988. v.38. N 6. P.3043-3053.

101. Jun-Ming L., Zhi-Guo L., Zhuang-Chun W. In-situ observations of dendritic growth of ammonium chloride crystals from an aqueous solution system// Scripta Metallurgica et Materialia. 1995. v. 32. N.3. P. 445-450.

102. Аксельрод Е.Г., Беспалов B.B., Добрин B.A и др. Акустическая эмиссия при фазовом переходе первого рода в жидком кристалле // Доклады РАН. 1995. т. 345. N 3. С. 320-323.

103. Дельмон Б. Кинетика гетерогенных реакций. М.: Мир, 1972. 554 с.

104. Дюк В. Обработка данных на ПК в примерах. Спб: Питер, 1997. 240с.

105. Тюрин Ю.Н., Макаров А.А. Статистический анализ данных на компьютере.М.: ИНФРА-М, 1998. 528 с.

106. Cladis Р.Е. New liquid-crystal phase diagram// Phys. Rev. Lett. 1975. v.35.N 1. P.48-51.

107. Cladis P.E. The reentrant nematic, enhanced smectic A phases and molecular comosition// Mol. Cryst. Liq. Cryst. 1981. v. 67, p. 177-192.

108.Баренблатт Г.И. Подобие, автомодельность, промежуточная асимптотика. JI.: Гидрометеоиздат. 1982, 230 с.

109.Казенин Д.А., Макеев А.А., Марков А.В. Автомодельный рост сферического кристалла в бинарном растворе// Теоретические основы химической технологии 1987. т.21. N 6. С. 763-769.

110. Мартюшев Л.М., Селезнев В.Д., Скопинов С.А. Компьютерное моделирование кристаллизации соли на подложке с помощью метода диффузионных потоков// Письма в ЖТФ. 1996. Т. 22. Вып. 4. С. 28-33.

111. Мартюшев Л.М., Селезнев В.Д., Скопинов С.А. Изучение роста скелетного кристалла в двумерной среде с фазовым расслоением с помощью метода диффузионных потоков // Письма в ЖТФ. 1996. Т. 22. Вып. 16. С. 12-17.

112. Martiouchev L.M., Skopinov S.A. Simulation of fractal growth of salts from biofluids on a substrate// Medical & Biological engineering & Computing. 1997. vol. 35. Supplement Part 1. P. 484 ( Abstracts of the World Congress on Medical Physics and Biomedical Engineering. 1997, Nice, France).

113.Бердышева O.B., Мартюшев Л.М. Изучение механизма кристаллизации соли при тезиграфических исследованиях в медицине// Безопасность Биосферы. Первый всероссийский научный молодежный симпозиум. Екатеринбург. 18-20 декабря 1997. Сб. тезисов. С. 202.

1 Н.Кузнецова И.Е., Мартюшев Л.М. Анализ устойчивости экосистемы с использованием термодинамического подхода. Эволюция шарообразной колонии// Безопасность Биосферы. Первый всероссийский научный молодежный симпозиум. Екатеринбург. 1820 декабря 1997. Сб. тезисов. С. 215.

115. Мартюшев Л.М., Селезнев В.Д., Скопинов С.А. Кинетические возвратные фазовые переходы при дендритном росте кристаллов в двумерной среде , с фазовым расслоением // Письма в ЖТФ. 1997. Т. 23. Вып. 13. С. 1-6.

116. Мартюшев Л.М., Скопинов С.А. Изучение роста кристалла в двумерной среде с фазовым расслоением методом компьютерного моделирования// сб. науч. тр. Метастабильные состояния и фазовые переходы. Екатеринбург: Институт Теплофизики УрО РАН. 1997. С. 105-111

117. Мартюшев Л.М., Скопинов С.А. Компьютерное моделирование

кристаллизации соли из биожидкостей// Сб. науч. тр. 1-й Всероссийской научно-практической конференции:

Кристаллографические методы исследования в медицине. Москва. 1997. С. 27-33.

118.Мартюшев JT.M., Селезнев В.Д., Скопинов С.А. Компьютерное моделирование потери устойчивости и развития дендритных форм методом диффузионных потоков// Кристаллография. 1997. Том 42. Вып.5. С. 735-741.

119. Martiouchev L.M., Seleznev V.D., Skopinov S.A. The computer simulation of nonequilibrium growth of crystals in the two-dimensional medium with a phase separating impurity // Journal of Statistical Physics 1998 v.90 (5/6). P.1413-1427.