Авторефераты по всем темам >>
Авторефераты по техническим специальностям
На правах рукописи
Александров Юрий Владимирович
МОДЕЛИРОВАНИЕ УПРУГИХ СВОЙСТВ ДВУМЕРНЫХ ЭЛЕКТРИЧЕСКИ СТАБИЛИЗИРОВАННЫХ КОЛЛОИДНЫХ КРИСТАЛЛОВ
Специальность: 05.13.18 - Математическое моделирование, численные методы и комплексы программ
Автореферат диссертации на соискание ученой степени кандидата физико-математических наук
Ульяновск - 2012
Работа выполнена на кафедре Системы автоматизированного проектирования Ульяновского государственного технического университета.
Научный руководитель - кандидат физико-математических наук, доцент Дышловенко Павел Евгеньевич
Официальные оппоненты:
Самохвалов Михаил Константинович, д.ф-м.н., профессор, зав. кафедрой Проектирование и технология электронных средств Ульяновского государственного технического университета;
Нагорнов Юрий Сергеевич, к.ф-м.н., доцент, доцент кафедры высшей математики Ульяновского государственного педагогического университета им. И.Н. Ульянова.
Ведущая организация - Ульяновский филиал Института радиотехники и электроники им. В.А. Котельникова Российской академии наук.
Защита состоится 23 мая 2012 г. в 10 ч. 00 мин. на заседании диссертационного совета Д 212.277.02 при Ульяновском государственном техническом университете по адресу: 432027, г. Ульяновск, ул. Северный Венец, 32, ауд. 211 (Главный корпус).
С диссертацией можно ознакомиться в библиотеке Ульяновского государственного технического университета.
Автореферат разослан л_______________2012 г.
Ученый секретарь диссертационного совета д.т.н., профессор В. Р. Крашенинников
ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ
Актуальность исследования Электрически стабилизированные коллоидные системы - это простейшие представители обширного класса суспензий заряженных объектов. Примерами таких систем чрезвычайно разнообразны и имеют многочисленные технологические применения [1]. Эти системы также играют огромную роль в молекулярной биологии, поскольку практически все протеины в каждой живой клетке, также как и сама молекула ДНК, являются заряженными макромолекулами, растворенными в соленой воде.
За последние десятилетия предложено несколько способов описания таких систем, отличающихся разной степенью детальности. Теория на основе нелинейного дифференциального уравнения Пуассона-Больцмана занимает в этой связи центральное место, поскольку, с одной стороны, с хорошей точностью описывает многие особенности электрического взаимодействия коллоидных систем, а с другой стороны, является базой, относительно которой проверяются все остальные теории. Значение этой теории особенно возросло в связи с появлением технологических возможностей получения коллоидных систем с частицами все меньшего размера, соизмеримого с длиной Дебая. Для таких систем становятся существенными нелинейные эффекты, не описываемые линеаризованными теориями.
Несмотря на то, что теоретические основы описания электрически стабилизированных коллоидных систем на основе уравнения ПуассонаБольцмана хорошо разработаны [2, 3], применение этой теории ограничивается в основном простейшими системами, что связано со сложностью численного решения нелинейного дифференциального уравнения. В связи с этим особое значение приобретает разработка точных и универсальных методов математического моделирования физически интересных систем, то есть систем с разнообразием электрических свойств, сложной геометрией и большим числом частиц.
Важным примером таких систем являются электрически стабилизированные коллоидные кристаллы, то есть системы, в которых частицы пространственно упорядочены. Как и обычные кристаллы, они обладают определенными свойствами упругости, но, в отличие от них, являются системами с начальным напряжением. Математическое моделирование упругих свойств электрически стабилизированных коллоидных кристаллов дополняет натурный эксперимент и позволяет получить сведения о силовых постоянных и модулях упругости этих сред.
Сведения об упругих свойствах важны для технологических применений.
Кроме того, моделирование позволяет получить информацию о характере взаимодействия в электрически стабилизированных системах, об их акустических свойствах и о фазовых переходах в них, а также о свойствах эффективных взаимодействий частиц. Все вышесказанное обосновывает актуальность темы диссертационного исследования.
Цель и задачи исследования Целью работы является исследование силовых и упругих постоянных двумерных электрически-стабилизированных коллоидных кристаллов средствами математического моделирования.
Для достижения поставленной цели решаются следующие задачи:
1. Разработка методики компьютерного эксперимента по определению силовых и упругих постоянных двумерных электрическистабилизированных коллоидных кристаллов.
2. Создание программного комплекса, реализующего методику определения силовых и упругих постоянных.
3. Проведение экспериментов по определению зависимостей силовых и упругих постоянных двумерных электрически-стабилизированных коллоидных кристаллов от параметров моделей.
4. Анализ результатов экспериментов, формулировка выводов по их итогам.
Методы исследования При решении поставленных задач применялись методы теории дифференциальных уравнений, математического моделирования, вычислительной математики, теории упругости, а также средства программирования на языках высокого уровня и библиотеки программ.
Научная новизна положений, выносимых на защиту 1. В рамках классической теории на основе нелинейного дифференциального уравнения Пуассона-Больцмана построены четыре модели двумерных электрически стабилизированных коллоидных структур, отличающихся: 1) наличием пространственной периодичностью, 2) включением большого число взаимодействующих частиц.
2. Разработан новый алгоритм определения силовых и упругих постоянных коллоидных кристаллов в рамках построенных моделей, в котором учтены: 1) свойства симметрии моделей, 2) вклады ближайших соседей высоких порядков.
3. Разработан новый программный комплекс, позволяющий проводить вычислительные эксперименты по определению силовых и упругих постоянных двумерных электрически стабилизированных коллоидных кристаллов в рамках теории на основе уравнения Пуассона-Больцмана при различных межчастичных расстояниях, размерах и электрических параметрах макроионов.
4. Обнаружены существенные отклонения от соотношений Коши для упругих постоянных коллоидных кристаллов рассматриваемых типов, что свидетельствует о многочастичном характере эффективного взаимодействия макроионов в таких системах.
5. Впервые установлено, что упругие свойства двумерных электрически стабилизированных коллоидных кристаллов практически полностью исчерпываются вкладами ближайших соседей 1-го и 2-го порядков.
Вклады соседей более высокого порядка пренебрежимо малы.
Достоверность результатов, представленных в диссертации, обеспечивается корректным применением классической теории на основе уравнения Пуассона-Больцмана, применением апробированного математического аппарата и численных методов, проверкой предельных и специальных случаев, сопоставлением результатов с литературными данными.
Практическая ценность результатов работы Разработанные алгоритмы и программы могут быть использованы при разработке технологических процессов создания фотонных и фононных кристаллов, а также новых перспективных материалов на основе глобулярных структур.
Реализация результатов работы Результаты диссертационной работы использованы при выполнении гранта РФФИ №09-01-97012 Математическое моделирование упругих и решеточных свойств наноразмерных коллоидных кристаллов.
Апробация работы Основные результаты работы докладывались на 42-й, 43-й, 44-й и 45-й научно-технических конференциях УГТУ Вузовская наука в современных условиях (Ульяновск, УГТУ, 28 января - 4 февраля 2008 г., 26-31 января 2009 г., 1-7 февраля 2010 г., 24-29 января 2011 г.), Всероссийской конференций Проведение научных исследований в области обработки, хранения, передачи и защиты информации (Ульяновск, 1-5 декабря 2009 г.), 7-й Всероссийской научной конференции с международным участием Математическое моделирование и краевые задачи (Самара, 3-6 июня 2010 г.), 3-й Российской научно-технической конференции аспирантов, студентов и молодых ученых ИВТ-2011 (Ульяновск, 24-25 мая 2011 г.), 7-й Всероссийской научнопрактической конференции (с участием стран СНГ) Современные проблемы создания и эксплуатации радиотехнических систем (Ульяновск, 22-сентября 2011 г.) Получены два свидетельства о государственной регистрации программы для ЭВМ Двумерный коллоидный кристалл с квадратной решеткой в модели уравнения Пуассона-Больцмана (quad5), № 2011616981, М.:
РОСПАТЕНТ, 08.09.2011 и Численное решение уравнения ПуассонаБольцмана для двумерных коллоидных кристаллов с квадратной решеткой, № 2011617862, М.: РОСПАТЕНТ, 07.10.2011.
Публикации По теме диссертации опубликовано 13 печатных работ, одна из которых в издании из перечня ВАК, два зарегистрированных программных продукта, и десять в других изданиях, включая тематические сборники и материалы международных и всероссийских научно-технических конференций.
Структура и объем диссертации Диссертация состоит из введения, четырех глав, заключения, списка литературы, включающего 126 наименований, и 3 приложений. Общий объем диссертации составляет 138 страниц и содержит 16 таблиц и 50 рисунков.
ОСНОВНОЕ СОДЕРЖАНИЕ РАБОТЫ
Во введении обоснована актуальность темы исследования, сформулирована цель работы и основные решаемые задачи, представлены научные положения, выносимые на защиту и отмечена их новизна, охарактеризована достоверность результатов и их практическая значимость.
В первой главе на основе литературных источников дана характеристика электрически стабилизированных коллоидных систем и способов описания электростатических взаимодействий в таких системах.
Проанализирован подход к описанию электрических свойств коллоидов на основе нелинейного дифференциального уравнения Пуассона-Больцмана (ПБ), его обоснованность, основные приближения и связь с другими подходами, в том числе линеаризованными теориями. Отмечается, что точные численные решения были получены только для простых систем с небольшим числом частиц и простой геометрией; часто высокая точность решения достигалась за счет использования специфических свойств системы и ее симметрии.
Также проведен обзор литературы по упругим свойствам электрически стабилизированных коллоидных кристаллов.
Ввиду того, что коллоидные кристаллы являются системами с начальным напряжением, описание их упругих свойств имеет ряд особенностей по сравнению с обычными кристаллами. Поэтому отдельный параграф посвящен справке по теории упругости сред с начальным напряжением и уточнению определений и обозначений.
Во второй главе дается описание четырех математических моделей исследуемых двумерных коллоидных кристаллов, разрабатывается методика вычислительного эксперимента по определению силовых постоянных этих кристаллов, дается описание структуры матриц силовых постоянных, проводится анализ свойств симметрии моделей для сокращения требуемого объема вычислений.
Модель коллоидного кристалла представляет собой систему частиц, находящуюся в жидком электролите. Частицы являются абсолютно твердыми диэлектриками и могут рассматриваться как бесконечно длинные круглые стержни (в трех измерениях) либо диски (в плоской двумерной задаче). Электролит является бинарным симметричным одновалентным с валентностям ионов,. В положении равновесия частицы находятся в узлах своей решетки Бравэ, при этом используются два типа двумерных решеток: гексагональная и квадратная. Соответствующие модели показаны на рисунке 1.
Рисунок 1. Модели двумерных коллоидных кристаллов с гексагональной и квадратной кристаллическими решетками. R - радиус частицы, d - параметр решетки Область определения задачи включает в себя ячейки Вигнера-Зейтца частиц, являющихся ближайшими соседями произвольно выбранной центральной частицы. В модели с гексагональной решеткой учитываются соседи до третьего порядка, а с квадратной - до пятого порядка включительно.
Для электрического потенциала в области электролита справедливо уравнение ПБ [3], безразмерная форма которого для рассматриваемых моделей имеет вид. (1) При этом единицей измерения длины является длина Дебая, а электрический потенциал измеряется в единицах, где - равновесная концентрация -й компоненты в объеме электролита, - элементарный заряд, - относительная диэлектрическая проницаемость электролита, - электрическая постоянная, - валентность -й компоненты электролита, - постоянная Больцмана, - абсолютная температура.
Рассматриваются две модели поведения заряда на поверхности частиц:
модель постоянного потенциала (ПП) и модель постоянной плотности заряда (ПЗ). Выбор модели определяет граничные условия на поверхности частиц.
Для модели ПП это граничное условие имеет вид, (2) где - заданный потенциал на границе. Для модели ПЗ рассматривается приближение большой, по сравнению с материалом частиц, диэлектрической проницаемости электролита, что типично для водных растворов. В этом случае граничное условие в безразмерной форме имеет вид, (3) где - нормальная компонента электрического поля, а - постоянная плотность заряда на поверхности частицы. На внешней границе области определения действуют однородные граничные условия Неймана вида. (4) Таким образом, в обеих моделях, как ПП, так и ПЗ, задача вычисления электрического потенциала в системе сводится к краевой задаче для уравнения (1) исключительно для области электролита и не включает области внутри частиц.
Характеристика математических моделей двумерных коллоидных кристаллов представлена в таблице 1.
Таблица Характеристика математических моделей коллоидных кристаллов гексагональная квадратная Название характеристики ПП ПЗ ПП ПЗ Ближайшие соседи До 3-го порядка До 5-го порядка 19 ячеек Вигнера-Зейтца 25 ячеек Вигнера-Зейтца Область определения (гексагоны) за вычетом области (квадраты) за вычетом частиц области частиц Уравнение ПБ для 1:1 электролита Дифференциальное уравнение Граничные условия на частицах Граничные условия на внешних границах Параметры модели Силовые постоянные определяются как коэффициенты при квадратичном члене в разложении энергии кристалла по малым смещениям частиц из их положений равновесия:
, (5) где есть -компонента смещения частицы из ее положения равновесия, задаваемого вектором, и - векторы решетки Бравэ,,. Методика вычислительного эксперимента по определению силовых постоянных основана на формуле, (6) вытекающей из (5), где есть -компонента силы, действующей на частицу. Существенным преимуществом подхода на основы формулы (6) по сравнению с определение (5) является отсутствие необходимости вычислять энергию системы для множественных смещений различных пар частиц, что значительно снижает объем вычислений, требования к вычислительным ресурсам и к точности решения, а также избавляет от необходимости двойного численного дифференцирования. Вычисление силы в безразмерной форме осуществляется по формуле F, N e E E E2 ch 1I n dl, x, y, (7) где выражение в квадратных скобках является тензором напряжения, связанным с уравнением ПБ. Здесь - произвольный контур, охватывающий частицу, I - единичный тензор, n - вектор внешней единичной нормали к элементу dl контура интегрирования, e - соответствующий единичный базисный вектор декартовой системы координат. В качестве контура интегрирования выбирается граница соответствующей ячейки Вигнера-Зейтца. Процедура численного интегрирования предпочтительна в вычислительном отношении по сравнению с дифференцирование. Дополнительным преимуществом является необходимость интегрировать только по периметру области определения.
Основная идея методики определения силовых постоянных состоит во внесении в равновесную систему возмущения в виде смещения единственной частицы и определении сил, действующих на все частицы системы.
Смещению подвергается только центральная частица. Силовые постоянные согласно (6) получаются затем однократным численным дифференцированием зависимостей сил от смещений.
Систематический учет симметрии моделей позволяет дополнительно снизить объем вычислений. В ходе исследования было показано, что для полного определения силовых постоянных всех рассматриваемых моделей достаточно смещений вдоль только одной из осей, причем только в одном из двух направлений. В работе осуществляется смещение в положительном x направлении оси. Детальный анализ свойств симметрии и структура матриц силовых постоянных для четырех исследуемых моделей приведены в диссертации.
Поскольку при смещении центральной частицы в положительном x направлении оси сохраняется зеркальная симметрия относительно этой оси, при численном решении использовалась только половина (верхняя) исходной области определения, что также снижает объем вычислений и требования к вычислительным ресурсам.
В третьей главе описывается программный комплекс, реализующий алгоритм нахождения силовых и упругих постоянных коллоидных кристаллов, описанный в главе 2, процедура обработки первичных данных и приводятся результаты математического моделирования для силовых постоянных всех исследуемых моделей коллоидных кристаллов.
Алгоритм нахождения силовых и упругих постоянных коллоидных кристаллов реализован в виде комплекса программ. Блок-схема комплекса программ показана на рисунке 2.
В состав комплекса входит программа-драйвер, реализующая методику определения силовых и упругих постоянных, а также программы генерации исходной геометрии, адаптивного перестроения сеток, численного решения краевой задачи для дифференциального уравнения ПБ, обработки данных первичных протоколов и определения силовых постоянных, вычисления упругих постоянных и другие программы постпроцессорной обработки.
Программы комплекса написаны на языках Fortran 90, COMSOL Script, MATLAB.
Уравнение ПБ решалось методом конечных элементов с использованием нерегулярных треугольных сеток. Программа реализована в двух вариантах: с адаптивным перестроением сеток и с использованием градиентных сеток без адаптации. Использовались лагранжевы элементы второго порядка.
Входными данными являются тип решетки и параметры модели, указанные в таблице 1: радиус частиц R, параметр решетки d и один из двух 0 электрических параметров, или, для модели ПП и ПЗ соответственно.
Для любого из параметров может быть задан желаемый диапазон изменения значений.
После ввода исходных данных осуществляется развертка по смещениям центральной частицы. Смещения совершаются в положительном направлении x оси. Всего проводится 10 одинаковых шагов смещения так, что максимальное смещение составляло 10% от расстояния d 2R (край-край) между ближайшими частицами. На каждом шаге смещения решается краевая задача для дифференциального уравнения ПБ с надлежащими граничными условиями, соответствующими заданной модели. По получении решения дифференциального уравнения вычисляются силы, действующие на частицы системы. Кроме того, вычисляется давление в равновесной конфигурации.
Полученные данные для всех смещений выводятся в файл протокола первичных данных, образец которого находится в приложении к диссертации.
начало ввод физических параметров решетки Наложение граничных условий квадратная гексагональная выбор Построение сетки конечных эл-тов геометрической модели генератор геометрии генератор геометрии квадратной решетки гексагональной решетки Решение уравнения ПБ постоянного постоянной потенциала плотности заряда выбор Нахождение х- и y- компонент сил электрической модели ввод ввод нет Число шагов сдвига достигло максимального да Импорт геометрической модели да Наложение граничных условий Формирование протокола значений компонент сил и давления Построение сетки конечных эл-тов Решение уравнения ПБ Нахождение силовых постоянных Нахождение давления в системе Нахождение упругих постоянных Нахождение х- и y- компонент сил Формирование протокола силовых и упругих постоянных Сдвиг центральной частицы на шаг в положительном направлении оси х конец Рисунок 2. Блок-схема комплекса программ для определения силовых и упругих постоянных двумерных электрически стабилизированных коллоидных кристаллов На рисунке 3 показан пример типичных зависимостей компонент сил, действующих на частицы системы, от смещения центральной частицы.
Значения компонент сил для положительных берутся непосредственно из протоколов первичных данных. Данные для отрицательных смещений формируются из данных для положительных на основании соображений симметрии, изложенных в главе 2. Как показано в диссертации, силовые постоянные являются производными представленных зависимостей, вычисленными в нуле. Кривые на графиках аппроксимировались кубическими полиномами методом наименьших квадратов. Силовые постоянные получаются тогда как коэффициенты при линейном члене в полиномиальной аппроксимации, при этом определяются как сами постоянные, так и их ошибки.
Выходными данными являются значения силовых и упругих постоянных кристалла. Упругие постоянные вычисляются по данным о силовых постоянных согласно формулам, представленным в главе 4. Образец итогового протокола находится в приложении к диссертации.
Рисунок 3. Зависимости компонент сил от смещений центральной частицы для кристалла с квадратной решеткой в модели ПП, Решение краевой задачи в ходе вычислительного эксперимента проводилось с использованием градиентных сеток. Для оценки точности полученного решения проводилось сравнение величины силы, действующей на центральную частицу, при ее максимальном смещении из положения равновесия, с точным решением, в качестве которого выбиралось аналогичное решение, полученное с адаптивным перестроением сеток.
В обоих случаях компоненты силы вычислялись как по границе ячейки Вигнера-Зейтца,, так и по поверхности частицы, и и. Вычисления проводились для тестовой модели кристалла с гексагональной решеткой в модели ПЗ с включением ближайших соседей только 1-го порядка. Адаптивные сетки конечных элементов содержали порядка 500 000 степеней свободы, что соответствовало предельным возможностям имеющейся вычислительной техники. Градиентные сетки насчитывали только около 200 000 степеней свободы.
Анализ данных для тестовой модели показал, что относительная ошибка решения не превышает 0.01% для всех d. Это свидетельствует о том, что градиентные сетки дают для величины силы результаты, по точности сопоставимые с аналогичными результатами для адаптивных сеток, при условии, что интегрирование осуществляется по границе ячейки Вигнера-Зейтца. Это обстоятельство позволяет использовать градиентные сетки в задачах с большим числом частиц, когда адаптивное решение для достижения той же точности требует числа степеней свободы, которое превышает имеющиеся вычислительные возможности. Типичное значение относительного отклонения величин и,, также составляло несколько сотых долей процента. Это позволило численно убедиться в независимости интеграла для силы от контура интегрирования, что, в свою очередь, свидетельствовало о высоком качестве численного решения. Последнее обстоятельство постоянно контролировалось в ходе моделирования.
Рисунок 4. Силовые постоянные центральной частицы, и соседей 1-го A0 A1 Bпорядка, A2 и B2 соседей 2-го порядка для кристалла с квадратной решеткой в модели ПП Были проведены систематические исследования силовых постоянных четырех моделей коллоидных кристаллов в диапазоне значений параметров решетки. В качестве примера на рисунке 4 представлены результаты для кристалла с квадратной решеткой в модели ПП.
В обозначении констант нижний индекс указывает порядок ближайших соседей центральной частицы, которая имеет номер 0. Полные численные и графические данные по силовым постоянным для всех исследуемых моделей приведены в диссертации.
В четвертой главе на основе полученной в ходе вычислительного эксперимента информации о силовых постоянных модельных кристаллов исследуются такие вопросы, как упругие постоянные кристаллов, наличие многочастичных взаимодействий в таких системах и оценка вкладов ближайших соседей различных порядков в их упругие свойства.
Упругие постоянные 1-го порядка - это напряжение в исходной конфигурации при отсутствии деформации. Для кристаллов как с гексагональной, так и с квадратной решеткой напряжение в исходной конфигурации изотропно и определяется осмотическим давлением.
Давление вычисляется с помощью тензора напряжений с помощью выражения, аналогичного (7) [3].
а) б) Рисунок 5. Упругие постоянные 1-го (давление) и 2-го (модули упругости) порядка для кристалла с гексагональной решеткой в модели ПП Упругие постоянные 2-го порядка (модули упругости) могут быть выражены через силовые постоянные и упругие постоянные 1-го порядка [4].
Для обоих типов решеток исследуемых моделей имеется только три нетривиальные упругие постоянные:. Полные выражения, учитывающие вклады соседей различных порядков в выражения для модулей упругости гексагональной и квадратной решеток приведены в диссертации.
В качестве примера на рисунке 5 показаны зависимости упругих постоянных кристалла с гексагональной решеткой в модели ПП. Аналогичные результаты для остальных трех моделей представлены в диссертации.
Информация об упругих постоянных кристалла позволяет сделать некоторые выводы о том, является ли эффективное взаимодействие в системе парным. Если частицы, составляющие кристалл, взаимодействуют посредством центрального парного потенциала, то на упругие постоянные кристалла налагаются дополнительные условия, известные как соотношения Коши [5]. Для всех моделей, рассматриваемых в данной работе, как с гексагональной, так и с квадратной решеткой, имеется одно соотношение Коши вида (8) Для оценки степени присутствия многочастичных эффективных взаимодействий в исследуемых системах вводится параметр. (9) Параметр для систем с парным взаимодействием и отклоняется от 1 в случае присутствия многочастичных эффективных взаимодействий.
Зависимость от параметра решетки для кристалла с гексагональной решеткой в модели ПП и ПЗ показаны на рисунке 6. Аналогичные результаты для квадратной решетки представлены в диссертации.
а) б) Рис. 6. Проверка выполнимости соотношения Коши для кристалла с гексагональной решеткой в модели а) ПП и б) ПЗ Из представленных данных следует, что имеют место отклонения от соотношения Коши для упругих постоянных всех исследуемых моделей.
Преимущественно реализуется соотношение. Однако в случае кристалла с гексагональной решеткой в модели ПЗ при приближении параметра решетки к области контакта наблюдается обратное соотношение. Наибольшие отклонения наблюдаются при малых и средних значениях параметра решетки. При увеличении параметра решетки парный характер эффективного взаимодействия усиливается. Таким образом, эффективное парное взаимодействие в коллоидных кристаллах, вообще говоря, не исчерпывается суммой парных взаимодействий независимо от конкретного вида парного потенциала. Наряду с парными должны учитываться и многочастичные вклады, роль которых, особенно при высоких плотностях, может быть значительна.
Полученные данные о силовых постоянных позволяют провести оценку относительного вклада соседей различных порядков в общее силовое взаимодействие в кристалле. Оценка основана на одном из свойств силовых постоянных, которое может быть записано в виде M Ca C (10) a MОно справедливо для любой фиксированной пары индексов.
При суммировании лишь по ограниченному числу ближайших соседей соотношение будет выполняться лишь приближенно. Были предложены M M S1 параметры S1 и S2 следующего вида: и S2 C Cx x x x C Cx, x x x MMгде - вектор решетки Бравэ для ближайших соседей только первого Mпорядка, а M2 - для ближайших соседей первого и второго порядка вместе.
Близость этих параметров к единице является мерой влияния ближайших соседей указанных порядков.
Зависимость параметров S1 и S2 от параметра решетки d для кристалла с гексагональной решеткой в модели ПП и ПЗ показаны на рисунке 7.
Аналогичные результаты для квадратной решетки представлены в диссертации. На графиках всех моделей наблюдаются минимумы значений, приходящиеся на область средних значений параметра решетки. Наибольшее отклонение параметра S1 от 1 имеет место в кристалле с квадратной решеткой в модели ПЗ. В то же время максимальное отклонение от параметра S2 для всех моделей не превышает 1%.
а) б) Рисунок 7. Параметры S1 и S2 от параметра решетки для для кристалла с гексагональной решеткой в модели а) ПП и б) ПЗ Из полученных данных следует, что силовое взаимодействие в двумерных коллоидных кристаллах всех рассматриваемых типов практически полностью исчерпывается взаимодействием ближайших соседей первого и второго порядков; вклады соседей более высоких порядков пренебрежимо малы для любых плотностей вплоть до полного контакта макроионов. Более детальные заключения для каждой из моделей представлены в диссертации.
В заключении сформулированы основные результаты диссертационной работы.
1. В рамках теории на основе нелинейного дифференциального уравнения Пуассона-Больцмана построены четыре математические модели двумерных коллоидных кристаллов с гексагональной и квадратной кристаллической решеткой с постоянным потенциалом и постоянной плотностью заряда частиц, в которых учтены взаимодействия ближайших соседей высоких порядков.
2. Для предложенных моделей коллоидных кристаллов разработан алгоритм определения силовых и упругих постоянных, основанный на внесении возмущения в равновесное состояние кристалла путем смещения одной из частиц.
3. Показано, что учет симметрии моделей дает возможность значительно, на порядок, снизить объем вычислений при математическом моделировании, а также снизить требования к оперативной памяти.
4. Разработан программный комплекс, реализующий алгоритм определения силовых и упругих постоянных двумерных коллоидных кристаллов в рамках теории на основе уравнения Пуассона-Больцмана.
5. Проведено математическое моделирование упругих свойств двумерных коллоидных кристаллов в рамках предложенных моделей.
В результате получены высокоточные данные о силовых и упругих постоянных кристалла в диапазоне параметров решетки от (почти) контакта частиц до расстояний, на которых взаимодействие ничтожно мало.
6. Обнаружены существенные отклонения от соотношений Коши для упругих постоянных коллоидных кристаллов рассматриваемых типов, что свидетельствует о многочастичном характере эффективного взаимодействия макроионов в таких системах.
7. Получено выражение для оценки относительных вкладов ближайших соседей различных порядков в упругие свойства двумерных коллоидных кристаллов. Показано, что в рассматриваемых моделях упругие свойства практически полностью исчерпываются вкладами ближайших соседей 1-го и 2-го порядков. Вклады соседей более высокого порядка пренебрежимо малы.
ЦИТИРУЕМАЯ ЛИТЕРАТУРА 1. Ролдугин, В. И. Физикохимия поверхности / В. И. Ролдугин. - Долгопрудный: Издательский дом Интеллект, 2008. - 568 с.
2. Дерягин, Б. В. Поверхностные силы / Б. В. Дерягин, Н. В. Чураев, В. М. Муллер. - М. : Наука, 1987. - 398 с.
3. Belloni, L. Colloidal interaction / L. Belloni // J. Phys.: Condens. Matter, 2000. - Vol. 12. - Pp. R549-R587.
4. Wallace D. C. Lattice Dynamics and Elasticity of Stressed Crystals / D. C. Wallace // Rev. Mod. Phys. - 1965. - V. 37. - Pр. 57-67.
5. Barron, T. H. K. Second-order elastic constants of a solid under stress / T. H. K. Barron, M. L. Klein // Proc. Phys. Soc. - 1965. - V.85. - Pp. 523-532.
ОСНОВНЫЕ ПУБЛИКАЦИИ ПО ТЕМЕ ДИССЕРТАЦИИ Научные статьи, опубликованные в изданиях из списка ВАК:
1. Александров, Ю. В. Силовые постоянные двумерного коллоидного кристалла с квадратной статической решеткой / Ю. В. Александров, П. Е. Дышловенко // Известия вузов. Электроника. - 2011. - № 1(87). - С. 9-16.
Свидетельства о государственной регистрации программы для ЭВМ:
2. Дышловенко, П. Е. Двумерный коллоидный кристалл с квадратной решеткой в модели уравнения Пуассона-Больцмана (quad5) / П. Е. Дышловенко, Ю. В. Александров // Свидетельство о государственной регистрации программы для ЭВМ № 2011616981, М.: РОСПАТЕНТ, 08.09.2011.
3. Дышловенко, П. Е. Численное решение уравнения ПуассонаБольцмана для двумерных коллоидных кристаллов с квадратной решеткой / П. Е. Дышловенко, Ю. В. Александров // Свидетельство о государственной регистрации программы для ЭВМ № 2011617862, М.: РОСПАТЕНТ, 07.10.2011.
Публикации в других изданиях:
4. Александров, Ю. В. Наноразмерный коллоидный кристалл с постоянным потенциалом макроионов / Ю. В. Александров, П. Е. Дышловенко // Вузовская наука в современных условиях: тезисы докладов 42-й научно-технической конференции УГТУ. - Ульяновск :
УГТУ, 2008. - С. 135.
5. Александров, Ю. В. Модель электрически стабилизированного двумерного коллоидного кристалла / Ю. В. Александров // Радиоэлектронная техника: межвузовский сборник научных трудов. - Ульяновск : УГТУ, 2008.
- С. 130-132.
6. Александров, Ю. В. Роль граничных условий в задаче определения силовых постоянных двумерного коллоидного кристалла / Ю. В. Александров, П. Е. Дышловенко // Вузовская наука в современных условиях: тезисы докладов 43-й научно-технической конференции УГТУ. - Ульяновск : УГТУ, 2009. - С. 27.
7. Александров, Ю. В. Роль граничных условий при моделировании силовых постоянных электрически стабилизированных коллоидных кристаллов / Ю. В. Александров, П. Е. Дышловенко // Сборник научных трудов всероссийской конференции Проведение научных исследований в области обработки, хранения, передачи и защиты информации, т.3. - Ульяновск : УГТУ, 2009. - С. 214-218.
8. Александров, Ю. В. Эффекты конечности области определения задачи при численном моделировании силовых постоянных коллоидного кристалла / Ю. В. Александров, П. Е. Дышловенко // Радиоэлектронная техника: межвузовский сборник научных трудов. - Ульяновск : УГТУ, 2009.
- С. 147-150.
9. Александров, Ю. В. Силовые постоянные коллоидного наноразмерного кристалла с постоянным потенциалом частиц / Ю. В. Александров, П. Е. Дышловенко // Вузовская наука в современных условиях: тезисы докладов 44-й научно-технической конференции УГТУ. - Ульяновск : УГТУ, 2010. - С. 164.
10. Александров, Ю. В. Моделирование упругих и решеточных свойств электрически стабилизированных коллоидных кристаллов со статической квадратной решеткой / Ю. В. Александров, П. Е. Дышловенко, А. Ф. Низаметдинов // Математическое моделирование и краевые задачи:
М33. Труды седьмой Всероссийской научной конференции с международным участием. Ч. 1: Математические модели механики, прочности и надежности элементов конструкций. - Самара : СамГТУ, 2010. - С.12-15.
11. Александров, Ю. В. Силовые и упругие постоянные двумерных коллоидных кристаллов / Ю. В. Александров, П. Е. Дышловенко // Вузовская наука в современных условиях: тезисы докладов 45-й научно-технической конференции УГТУ. - Ульяновск : УГТУ, 2011. - С. 185.
12. Александров, Ю. В. Упругие свойства коллоидного кристалла с квадратной решеткой в модели постоянного потенциала / Ю. В. Александров // Сборник научных трудов 3-й Российской научно-технической конференции аспирантов, студентов и молодых ученых ИВТ-2011. - Ульяновск : УГТУ, 2011. - С. 18-22.
13. Александров, Ю. В. Моделирование силовых постоянных двумерного коллоидного кристалла с гексагональной статической решеткой / Ю. В. Александров, П. Е. Дышловенко, А. Ф. Низаметдинов, Д. В. Чернятьев // Современные проблемы создания и эксплуатации радиотехнических систем: труды седьмой всероссийской научно-практической конференции (с участием стран СНГ. - Ульяновск : УГТУ, 2011. - С.177-180.
Александров Юрий Владимирович Моделирование упругих свойств двумерных электрически стабилизированных коллоидных кристаллов Автореферат Подписано в печать 17.04.2012. Формат 6084/16.
Усл. печ. л. 1,16. Тираж 100 экз. Заказ 424.
Типография УГТУ, 432027, г. Ульяновск, Северный Венец, 32.
Авторефераты по всем темам >>
Авторефераты по техническим специальностям