Редактирование: Тепловые колебания в одномерном кристалле
Внимание! Вы не авторизовались на сайте. Ваш IP-адрес будет публично видимым, если вы будете вносить любые правки. Если вы войдёте или создадите учётную запись, правки вместо этого будут связаны с вашим именем пользователя, а также у вас появятся другие преимущества.
Правка может быть отменена. Пожалуйста, просмотрите сравнение версий, чтобы убедиться, что это именно те изменения, которые вас интересуют, и нажмите «Записать страницу», чтобы изменения вступили в силу.
Текущая версия | Ваш текст | ||
Строка 24: | Строка 24: | ||
Принимались периодические граничные условия: | Принимались периодические граничные условия: | ||
− | <math> {u_{k+N} =u_k, ~~~ N>>1 } ~~~~~ (2) </math> | + | <math> {u_{k+N} =u_k, ~ ~~ N>>1 } ~~~~~ (2) </math> |
где <math> N</math> – число независимых частиц. | где <math> N</math> – число независимых частиц. | ||
Строка 41: | Строка 41: | ||
<math> K=m\left \langle {\dot u_k^2} \right \rangle ~~~~~ (4) </math> | <math> K=m\left \langle {\dot u_k^2} \right \rangle ~~~~~ (4) </math> | ||
− | Для анализа зависимости кинетической энергии от параметра α были проведены численные эксперименты для одномерного кристалла с нелинейным взаимодействием частиц (1) с начальными условиями (3) и граничными условиями (2). Здесь и далее кинетическая энергия была обезразмерена по начальной заданной кинетической энергии, а время расчётов – по отношению к периоду колебаний частицы вблизи положения равновесия. | + | Для анализа зависимости кинетической энергии от параметра α были проведены численные эксперименты для одномерного кристалла с нелинейным взаимодействием частиц (1) с начальными условиями (3) и граничными условиями (2). Здесь и далее кинетическая энергия была обезразмерена по начальной заданной кинетической энергии, а время расчётов – по отношению к периоду колебаний частицы вблизи положения равновесия . |
На рис. 1 представлены результаты численного моделирования задачи об определении кинетической температуры одномерного кристалла с нелинейным взаимодействием. | На рис. 1 представлены результаты численного моделирования задачи об определении кинетической температуры одномерного кристалла с нелинейным взаимодействием. | ||
− | [[File:Simonov 1.png |300px| | + | [[File:Simonov 1.png |left|300px| а) <math> \alpha=0 </math>]] [[File:Simonov 2.png |300px| б) <math> \alpha=1 </math>]] |
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
==Колебания кинетической температуры== | ==Колебания кинетической температуры== | ||
+ | Данный раздел посвящен исследованию колебаний амплитуды синусоидального профиля кинетической температуры одномерного кристалла и их зависимости от величины введенной нелинейности взаимодействия частиц. | ||
Как показано ранее, увеличение нелинейности взаимодействия частиц приводит к увеличению скорости затухания кинетической энергии в одномерном кристалле. Высокочастотные колебания наблюдают и при задании начального синусоидального распределения температуры по кристаллу. Синусоидальное распределение температуры можно представить как случай неравномерного нагрева материала по его длине. В этом случае инициализируется колебательный процесс, связанный со стремлением кинетической температуры равномерно распределиться по кристаллу. При этом колебания амплитуды профиля заданной температуры со временем будут затухать. | Как показано ранее, увеличение нелинейности взаимодействия частиц приводит к увеличению скорости затухания кинетической энергии в одномерном кристалле. Высокочастотные колебания наблюдают и при задании начального синусоидального распределения температуры по кристаллу. Синусоидальное распределение температуры можно представить как случай неравномерного нагрева материала по его длине. В этом случае инициализируется колебательный процесс, связанный со стремлением кинетической температуры равномерно распределиться по кристаллу. При этом колебания амплитуды профиля заданной температуры со временем будут затухать. | ||
− | Для исследования распределения кинетической температуры в нелинейном одномерном кристалле (1) с периодическими граничными условиями (2), рассмотрены следующие начальные условия: | + | Для исследования распределения кинетической температуры в нелинейном одномерном кристалле (1) с периодическими граничными условиями (2), рассмотрены следующие начальные условия: |
− | |||
− | |||
− | где <math> n </math> – номер частицы, <math> a </math> – расстояние между ними, | + | где <math> n </math> – номер частицы, <math> a </math> – расстояние между ними, – случайная величина с равномерным распределением, а <math> σ(х) </math> - начальные распределения скоростей частиц в цепочке. Величину <math> σ(х) </math> определим так, чтобы получить начальное распределение температуры по синусоидальному закону: |
− | <math> | + | где <math> A </math> и <math> \kappa </math> – константы, <math> \kappa=2\pi / N </math>, <math> σ(х) </math> – значение средней кинетической температуры в одномерном кристалле. |
− | + | Определен способ вычисления колебаний амплитуды синусоидальной температуры. Исследована амплитуда колебаний кинетической температуры через разложение кинетической температуры одномерной цепочки в ряд Фурье. Для этого вычислен первый коэффициент при разложении кинетической температуры по синусу: | |
− | |||
− | + | ==Сравнение зависимостей параметров двух законов затухания от нелинейности== | |
− | |||
− | |||
− | |||
− | |||
− | + | * В двух рассмотренных задачах были определены законы затухания колебаний кинетической энергии и амплитуды кинетической температуры в одномерном нелинейном кристалле. Закон затухания для обеих величин определялся в виде степенных функций. Определены два параметра для каждого из законов затухания и установлены их зависимости от введенной нелинейности. | |
− | |||
− | |||
− | + | * Из представленных результатов можно сделать вывод о том, что зависимость параметра a в законе затухания для колебания кинетической энергии и амплитуды кинетической температуры имеют схожий вид. Результаты, представленные на рис. 6, показывают, что показатели степени (b) в законе затухания в двух рассмотренных задачах имеют схожий вид зависимости от параметра нелинейности, однако значения параметров отличаются друг от друга на 20%. Коэффициент а, для двух рассмотренных величин имеет нелинейную зависимость от параметра α, при этом значения коэффициентов двух законов затухания отличаются на 30%. Таким образом, установлено различие в законах затухания при рассмотрении нелинейного взаимодействия | |
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
==Заключение== | ==Заключение== | ||
− | |||
* В работе исследован одномерный кристалл с нелинейным взаимодействием частиц (кубическая нелинейность в выражении для взаимодействия силы). Для проведения исследования введен безразмерный параметр <math> \alpha </math>, характеризующий отношение влияния линейного и нелинейного взаимодействия между частицами одномерного кристалла. Рассмотрены две постановки начальных условий: равномерное задание температуры и синусоидальный профиль температуры в одномерном кристалле. | * В работе исследован одномерный кристалл с нелинейным взаимодействием частиц (кубическая нелинейность в выражении для взаимодействия силы). Для проведения исследования введен безразмерный параметр <math> \alpha </math>, характеризующий отношение влияния линейного и нелинейного взаимодействия между частицами одномерного кристалла. Рассмотрены две постановки начальных условий: равномерное задание температуры и синусоидальный профиль температуры в одномерном кристалле. | ||
* Для первой постановки задачи исследована кинетическая энергия кристалла и ее зависимость от параметра нелинейности <math> \alpha </math>. Получено, что зависимость кинетической энергии одномерного кристалла от времени носит колебательный характер и с увеличением нелинейности скорость затухания колебаний увеличивается. Предложена аппроксимация закона затухания кинетической энергии кристалла для различного значения параметра нелинейности α. Введённая аппроксимация дает хорошее соответствие для малой нелинейности, при большой нелинейности наблюдается расхождение графика колебаний энергии и закона затухания. | * Для первой постановки задачи исследована кинетическая энергия кристалла и ее зависимость от параметра нелинейности <math> \alpha </math>. Получено, что зависимость кинетической энергии одномерного кристалла от времени носит колебательный характер и с увеличением нелинейности скорость затухания колебаний увеличивается. Предложена аппроксимация закона затухания кинетической энергии кристалла для различного значения параметра нелинейности α. Введённая аппроксимация дает хорошее соответствие для малой нелинейности, при большой нелинейности наблюдается расхождение графика колебаний энергии и закона затухания. | ||
Строка 136: | Строка 81: | ||
# Кривцов А.М., Ле-Захаров А.А. Исследование процесса теплопроводности кристаллах с дефектами методами молекулярной динамики// ДАН. 2008. Т. 420. № 1. С. 45 – 49. | # Кривцов А.М., Ле-Захаров А.А. Исследование процесса теплопроводности кристаллах с дефектами методами молекулярной динамики// ДАН. 2008. Т. 420. № 1. С. 45 – 49. | ||
# E. Brown, L. Hao, J.C. Gallop, J.C. Macfarlane. Ballistic thermal and electrical conductance measurements on individual multiwall carbon nanotubes. Appl. Phys. Lett. (2005) 87, 023107. | # E. Brown, L. Hao, J.C. Gallop, J.C. Macfarlane. Ballistic thermal and electrical conductance measurements on individual multiwall carbon nanotubes. Appl. Phys. Lett. (2005) 87, 023107. | ||
− | # Zhaohui Wang, Jeffrey A. Carter, Alexei Lagutchev, Yee Kan Koh, Nak-Hyun Seong, David G. Cahill, Dana D. Dlott. Ultrafast Flash Thermal Conductance of Molecular Chains. Science (2007) Vol. 317 no. 5839 pp. 787-790. | + | # Zhaohui Wang, Jeffrey A. Carter, Alexei Lagutchev, Yee Kan Koh, Nak-Hyun Seong, David G. Cahill, Dana D. Dlott. Ultrafast Flash Thermal Conductance of |
+ | Molecular Chains. Science (2007) Vol. 317 no. 5839 pp. 787-790. | ||
# C. W. Chang, D. Okawa, H. Garcia, A. Majumdar, and A. Zettl. Breakdown of Fourier’s Law in Nanotube Thermal Conductors. Phys. Rev. Lett. (2008) 101, 075903. | # C. W. Chang, D. Okawa, H. Garcia, A. Majumdar, and A. Zettl. Breakdown of Fourier’s Law in Nanotube Thermal Conductors. Phys. Rev. Lett. (2008) 101, 075903. | ||
# Giardina C., Livi R., Politi A., Vassalli M. // Phys. Rev. Lett. 2000. V. 84. № 10. P. 2144–2147. | # Giardina C., Livi R., Politi A., Vassalli M. // Phys. Rev. Lett. 2000. V. 84. № 10. P. 2144–2147. | ||
Строка 143: | Строка 89: | ||
# Кривцов А.М. Колебания энергий в одномерном кристалле ДАН. 2014, том 458, № 3, с. 279–281. | # Кривцов А.М. Колебания энергий в одномерном кристалле ДАН. 2014, том 458, № 3, с. 279–281. | ||
# Andrei A. Gusev, Sergey A. Lurie, Wave-relaxation duality of heat propagation in Fermi–Pasta–Ulam chains // Mod. Phys. Lett. B August 2012, Vol. 26, No. 22 | # Andrei A. Gusev, Sergey A. Lurie, Wave-relaxation duality of heat propagation in Fermi–Pasta–Ulam chains // Mod. Phys. Lett. B August 2012, Vol. 26, No. 22 | ||
− | |||
− |