Перераспределение энергии по пространственным направлениям в кристаллах — различия между версиями

Версия 22:50, 20 июня 2017

Выпускная квалификационная работа

Выполнил: студент группы 43604/1 Н.Г. Шварёв

Руководитель: кандидат физ.-мат. наук, доцент кафедры Теоретическая механика В.А. Кузькин

Введение

Количественное описание неравновесных тепловых процессов в кристаллах – одна из актуальных проблем современной физики. Это связано со стремительным развитием нанотехнологий. В связи с этим большой интерес представляют процессы, происходящие в твердых телах при переходе к состоянию термодинамического равновесия. Неравновесное состояние может быть вызвано, к примеру, прохождением ударных волн или быстрым лазерным воздействием. Тогда кинетические энергии теплового движения атомов в разных направлениях могут значительно различаться. Это, в свою очередь, показывает, что кинетическая температура может проявлять тензорные свойства. На фронте ударной волны, распространяющейся вдоль одной из осей, например, оси X, выполняются следующее соотношение:

[math]\widehat T_{xx} \gt \widehat T_{yy} [/math],

где [math]\widehat T_{xx}= \frac{{m\lt v_x^2\gt }}{k_B}, \widehat T_{yy}= \frac{{m\lt v_y^2\gt }}{k_B}[/math] - кинетические температуры вдоль соответствующих направлений, [math]k_B[/math] - постоянная Больцмана.

При переходе к равновесному состоянию в кристалле реализуется два процесса:

1)Выравнивание кинетической и потенциальной энергий;

2)Перераспределение кинетической энергии по пространственным направлениям.

Данная работа посвящена численному описанию перераспределения кинетической энергии по пространственным направлениям в негармонических кристаллах с треугольной кристаллической решеткой.

Далее будет рассматриваться обезразмеренное значение температуры: [math] T_{xx} = \frac{\widehat T_{xx}}{\widehat T_{xx}^o - \widehat T_{yy}^o}, T_{yy} = \frac{\widehat T_{yy}}{\widehat T_{xx}^o - \widehat T_{yy}^o} [/math],

где [math] \widehat T_{xx}^o , \widehat T_{yy}^o [/math] - значения в начальный момент времени [math] \widehat T_{xx}, \widehat T_{yy} [/math] соответственно

Из-за того, что температура прямо пропорциональна кинетической энергии, а в дальнейшем будет происходить рассмотрение только обезразмеренного значения температуры, то понятия температуры и кинетической энергии будут равносильны: [math] T_{xx} = E_x, T_{yy} = E_y [/math]

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

Цели и задачи работы

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

• рассмотрение процесса выравнивания температур;

• рассмотрение влияния нелинейности на поведение системы;

• выделение медленного процесса, вызываемого нелинейностью;

• определение формы выделенного медленного процесса.

Модель двумерного кристалла

Рис.1. Пример треугольной кристаллической решетки 10x10, получаемый в результате работы программы

• Рассматривается Треугольная кристаллическая решетка.

• Для взаимодействия между частицами используется Потенциал Леннард-Джонса.

• Задаются следующие начальные условия:

[math]v_x ≤ v_{max}, v_y = 0, u_x = 0, u_y = 0[/math].

В начальный момент времени рассматриваются случайные начальные скорости вдоль одной оси (оси X), ограниченные некоторым варьируемым максимальным значением [math]v_{max}[/math], нулевые начальные скорости вдоль другой оси (оси Y) и нулевые перемещения вдоль обеих осей.

• Используются периодические граничные условия Борна-Кармана.

Выравнивание температуры

• 

  Рис.2. Поведение кинетической температуры при усреднении по 100 реализациям

Как можно заметить, оба графика стремятся к асимптоте T = ¼ . Это связано с тем, что со временем при переходе к стационарному состоянию кинетическая и потенциальная энергия выравниваются, значит, половина кинетической энергии уходит в потенциальную. А при наличии нелинейности разность [math]T_{xx}-T_{yy}[/math] стремится к нулю, следовательно, половина от оставшейся половины уходит на равное распределение по пространственным направлениям.

На начальном интервале в несколько периодов [math]τ_o[/math] происходит перераспределение кинетической и потенциальной энергии, а далее – перераспределение кинетической энергии по пространственным направлениям.

Redistribution_of_energy.png

• 

  Рис.3. График перераспределения кинетической температуры по направлениям