Редактирование: Моделирование распространения тепла в треугольной кристаллической решетке
Внимание! Вы не авторизовались на сайте. Ваш IP-адрес будет публично видимым, если вы будете вносить любые правки. Если вы войдёте или создадите учётную запись, правки вместо этого будут связаны с вашим именем пользователя, а также у вас появятся другие преимущества.
Правка может быть отменена. Пожалуйста, просмотрите сравнение версий, чтобы убедиться, что это именно те изменения, которые вас интересуют, и нажмите «Записать страницу», чтобы изменения вступили в силу.
Текущая версия | Ваш текст | ||
Строка 4: | Строка 4: | ||
=== Кристаллическая решетка === | === Кристаллическая решетка === | ||
Рассматривается плоская кристаллическая скалярная решетка. Атомы этой решетки расположены в вершинах равных правильных треугольников и имеют одинаковую массу m. Каждый атом имеет одну степень свободы - он может двигаться вдоль прямой, перпендикулярной плоскости решетки. Схематически решетка представлена на рисунке 1. | Рассматривается плоская кристаллическая скалярная решетка. Атомы этой решетки расположены в вершинах равных правильных треугольников и имеют одинаковую массу m. Каждый атом имеет одну степень свободы - он может двигаться вдоль прямой, перпендикулярной плоскости решетки. Схематически решетка представлена на рисунке 1. | ||
− | [[File:Lattice | + | [[File:Lattice 1|thumb|left|Внешний вид кристаллической решетки]] |
− | Процесс распространения тепловой энергии в кристалле проходит на молекулярном уровне, как следствие материалы с правильной регулярной структурой могут обладать особыми свойствами теплопроводности, которые могут зависеть от направления. Для исследования свойств теплопроводности в каком-то одном из направлений необходимо смоделировать бесконечную плоскую решетку, задать температуру вдоль выделенного направления (при этом во всех точках, лежащих на одной прямой, перпендикулярной выделенному направлению температура будет одинаковой) и наблюдать за скоростью, с которой тепловая энергия будет распределяться в кристалле. | + | Процесс распространения тепловой энергии в кристалле проходит на молекулярном уровне, как следствие материалы с правильной регулярной структурой могут обладать особыми свойствами теплопроводности, которые могут зависеть от направления. Для исследования свойств теплопроводности в каком-то одном из направлений необходимо смоделировать бесконечную плоскую решетку, задать температуру вдоль выделенного направления (при этом во всех точках, лежащих на одной прямой, перпендикулярной выделенному направлению температура будет одинаковой) и наблюдать за скоростью, с которой тепловая энергия будет распределяться в кристалле. Чем быстрее температура выравнивается, тем больше теплопроводность кристалла в данном направлении. |
=== Периодичные начальные условия === | === Периодичные начальные условия === | ||
К сожалению, расчет бесконечной решетки (и как следствие бесконечного числа частиц) потребовал бы бесконечных вычислительных мощностей или бесконечного времени на расчет. Если ограничиться достаточно большой решеткой, но конечных размеров, это приведет к ухудшению точности в связи с влиянием эффектов на границе и потребует значительных ресурсов при вычислениях. Из этой ситуации позволяет выйти периодичность кристаллической решетки. Если задать начальное одномерное распределение температуры, которое можно задать также периодичным. Таким образом мы получаем, что всю бесконечную плоскость можно разбить на конечные участки, на каждом из которых будут одинаковые начальные условия и одинаковая геометрия решетки. Это позволяет заменить задачу про бесконечную плоскость на задачу о конечном участке с периодичными граничными условиями. В качестве периодичной функции температуры вдоль выделенного направления использован синус. Распределение температуры в пространстве задавалось законом: | К сожалению, расчет бесконечной решетки (и как следствие бесконечного числа частиц) потребовал бы бесконечных вычислительных мощностей или бесконечного времени на расчет. Если ограничиться достаточно большой решеткой, но конечных размеров, это приведет к ухудшению точности в связи с влиянием эффектов на границе и потребует значительных ресурсов при вычислениях. Из этой ситуации позволяет выйти периодичность кристаллической решетки. Если задать начальное одномерное распределение температуры, которое можно задать также периодичным. Таким образом мы получаем, что всю бесконечную плоскость можно разбить на конечные участки, на каждом из которых будут одинаковые начальные условия и одинаковая геометрия решетки. Это позволяет заменить задачу про бесконечную плоскость на задачу о конечном участке с периодичными граничными условиями. В качестве периодичной функции температуры вдоль выделенного направления использован синус. Распределение температуры в пространстве задавалось законом: | ||
::<math> | ::<math> | ||
− | T = T_0 ( | + | T = T_0 (sin (\omega x)+1.05) |
</math> | </math> | ||
в первом случае и | в первом случае и | ||
::<math> | ::<math> | ||
− | T = T_0 ( | + | T = T_0 (sin (\omega y)+1.05), |
</math> | </math> | ||
где x,y - декартовы координаты в плоскости решетки, слагаемое 1.05 вводится для строгой положительности температуры во всех точках | где x,y - декартовы координаты в плоскости решетки, слагаемое 1.05 вводится для строгой положительности температуры во всех точках | ||
Строка 20: | Строка 20: | ||
=== Периодические граничные условия === | === Периодические граничные условия === | ||
Бесконечную плоскость треугольной кристаллической решетки можно разбивать на фигуры различных форм и размеров. Но наиболее простыми и оптимальными вариантами являются шестиугольники и прямоугольники (конечно, две из сторон такого прямоугольника не будут прямыми из-за геометрии решетки). Поскольку начальные условия задачи имеют одно выделенное направление, то единственным подходящим вариантов ячейки периодичности решетки является прямоугольник. | Бесконечную плоскость треугольной кристаллической решетки можно разбивать на фигуры различных форм и размеров. Но наиболее простыми и оптимальными вариантами являются шестиугольники и прямоугольники (конечно, две из сторон такого прямоугольника не будут прямыми из-за геометрии решетки). Поскольку начальные условия задачи имеют одно выделенное направление, то единственным подходящим вариантов ячейки периодичности решетки является прямоугольник. | ||
− | [[File:lattice.png|200px|thumb| | + | [[File:lattice.png|200px|thumb|left|Внешний вид ячейки периодичности]] |
Синим контуром обведены частицы, которые являются уникальными (не повторяющимися), а черным контуром – частицы, данные о которых хранятся и обрабатываются в программе. Размеры решетки, используемой для расчетов много больше представленных на рисунке и составляют порядка сотни частиц. | Синим контуром обведены частицы, которые являются уникальными (не повторяющимися), а черным контуром – частицы, данные о которых хранятся и обрабатываются в программе. Размеры решетки, используемой для расчетов много больше представленных на рисунке и составляют порядка сотни частиц. | ||
− | == Описание алгоритма == | + | === Описание алгоритма === |
В данной программе решается задача для двух направлений – вдоль одного из направлений связей и перпендикулярно ему. Основные этапы расчета содержат: | В данной программе решается задача для двух направлений – вдоль одного из направлений связей и перпендикулярно ему. Основные этапы расчета содержат: | ||
− | Создание начальной конфигурации решетки, массива температур. Для каждого ряда (или столбца) высчитывается его температура по формуле | + | # Создание начальной конфигурации решетки, массива температур. Для каждого ряда (или столбца) высчитывается его температура по формуле |
::<math> | ::<math> | ||
− | T_0 | + | T = T_0 sin (\omega x). |
</math> | </math> | ||
− | Далее всем точкам этого ряда (столбца) задаются случайные скорости, такие, что значение квадрата их модуля не превышает температуры данного ряда (столбца) | + | Далее всем точкам этого ряда (столбца) задаются случайные скорости, такие, что значение квадрата их модуля не превышает температуры данного ряда (столбца) |
::<math> | ::<math> | ||
− | V_0 = random(-\sqrt{ | + | V_0 = random(-K \sqrt{T(x)}; K \sqrt{T(x)}} |
− | </math>. Положение всех точек в начальный момент времени соответствует равновесному положению. | + | </math>. Положение всех точек в начальный момент времени соответствует равновесному положению. Такие начальные условия не являются реалистичными, поскольку вся энергия системы является кинетической. Но, как показывает эксперимент, через некоторое время колебания потенциальной и кинетической энергии затухают и их значения во все последующие моменты времени практически равны (что и наблюдается в реальных системах). Поэтому такой подход допустим с оговоркой, что при временах меньше 2T0 (характерное время колебания системы) результаты нельзя использовать напрямую для прогнозирования поведения реальных систем. |
− | + | # Расчет одного шага по времени. На этом этапе проводилось интегрирование системы уравнений движения методом Эйлера, который имеет ошибку порядка O(h). Для удобства решения задачи, помимо «уникальных» частиц в массиве содержатся также их копии, которые образуют внешний слой рассматриваемого участка. Интегрирование производилось только для «уникальных» частиц. Остальные частицы (внешний слой) получались копированием «уникальных». Для этого использовались формулы | |
− | Расчет одного шага по времени. На этом этапе проводилось интегрирование системы уравнений движения методом Эйлера, который имеет ошибку порядка O(h). Для удобства решения задачи, помимо «уникальных» частиц в массиве содержатся также их копии, которые образуют внешний слой рассматриваемого участка. Интегрирование производилось только для «уникальных» частиц. Остальные частицы (внешний слой) получались копированием «уникальных». Для этого использовались формулы | ||
::<math> | ::<math> | ||
\ddot{u} = \frac{c}{m}(u_{i+1,j+1}+u_{i+1,j-1}+u_{i+1,j}+u_{i,j+1}+ u_{i,j-1}+u_{i-1,j} - 6 u_{i,j});\text{ } j =2n | \ddot{u} = \frac{c}{m}(u_{i+1,j+1}+u_{i+1,j-1}+u_{i+1,j}+u_{i,j+1}+ u_{i,j-1}+u_{i-1,j} - 6 u_{i,j});\text{ } j =2n | ||
Строка 43: | Строка 42: | ||
</math> | </math> | ||
# Осреднение по множеству реализаций. Поскольку температура – величина статистическая, то неправомерно считать ее как квадрат скорости частицы на некий коэффициент. Поэтому для достижения более высокой точности в программе создается несколько наборов начальных скоростей для частиц, каждый из которых отвечает начальному распределению температуры. При подсчете температуры используется осреднение значений по всем реализациям. | # Осреднение по множеству реализаций. Поскольку температура – величина статистическая, то неправомерно считать ее как квадрат скорости частицы на некий коэффициент. Поэтому для достижения более высокой точности в программе создается несколько наборов начальных скоростей для частиц, каждый из которых отвечает начальному распределению температуры. При подсчете температуры используется осреднение значений по всем реализациям. | ||
− | # Вычисление амплитуды. Амплитуда высчитывается как скалярное произведение температуры в данный момент времени и синуса той же частоты, что и в начальных условиях (интегрирование по периоду). | + | # Вычисление амплитуды. Этот этап необходим непосредственно для определения теплопроводности в данном направлении. Амплитуда высчитывается как скалярное произведение температуры в данный момент времени и синуса той же частоты, что и в начальных условиях (интегрирование по периоду). |
::<math> | ::<math> | ||
− | A = \int_0^l T(x) | + | A = \int_0^l T(x) sin(\omega x) dx |
</math> | </math> | ||
− | == Результаты == | + | === Результаты === |
− | |||
− | |||
− | |||
− | |||
− | |||
В результате численного эксперимента были получены следующий графики: | В результате численного эксперимента были получены следующий графики: | ||
+ | # График зависимости амплитуды от времени при T_max = 10 T0, dt = 0.001 T0, ось вдоль связей, 200 точек. | ||
− | + | [[File:Oy_10T0.jpg|thumb|амплитуда синуса по y]] | |
− | |||
− | [[File: | ||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | [[File: | + | # График зависимости амплитуды от времени при T_max = 10 T0, dt = 0.001 T0, ось перпендикулярна одному из направлений связей, 200 точек. |
+ | [[File:Ox_10T0.jpg|thumb|left|График зависимости амплитуды синуса температур от времени (вдоль оси x)]] | ||
+ | В результате анализа графиков было получено, что тепло распространяется в двух исследуемых направлениях одинаково, при том, что геометрия решетки не является симметричной относительно поворота на 90 градусов. Возможно, из этого эксперимента следует, что теплопроводность треугольной решетки вообще во всех направлениях одинакова.(но это не точно) | ||
− | == Ссылки == | + | === Ссылки === |
#[[Курсовые_работы_по_ВМДС:_2018-2019]] | #[[Курсовые_работы_по_ВМДС:_2018-2019]] | ||
#Автор работы - [[Давыдова Алена]] | #Автор работы - [[Давыдова Алена]] | ||
#Руководитель - [[Виталий Андреевич Кузькин]] | #Руководитель - [[Виталий Андреевич Кузькин]] | ||
− | |||
− | == Список использованной литературы == | + | === Список использованной литературы === |
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− |