"Численные методы интегрирования уравнений движения для одномерной линейной цепочки и частицы в потенциальной яме Леннарда-Джонса"

Материал из Department of Theoretical and Applied Mechanics
Версия от 22:19, 23 января 2020; 93.185.19.71 (обсуждение) (Новая страница: «'''''Курсовой проект по Механике дискретных сред''''' '''Исполни…»)

(разн.) ← Предыдущая | Текущая версия (разн.) | Следующая → (разн.)
Перейти к: навигация, поиск

Курсовой проект по Механике дискретных сред

Исполнитель: Троцкая Валерия

Группа: 3630103/60101

Семестр: осенний семестр 2019-2020 учебного года

Постановка задачи:

  1. Сравнить различные методы интегрирования уравнений движения (Верле, Рунге-Кутта 4 порядка). Реализовать фиксированные, свободные и периодические условия</div>
  2. Численно определить скорость диссоциации частицы в потенциальной яме Леннарда-Джонса

Теоретическая сводка:

  1. Одномерная линейная цепочка

Рассмотрим модель колебаний одинаковых атомов массой m, находящихся в одномерной цепочке. Пусть в этой цепочке находится N атомов, связанных между собой квазиупругой силой с коэффициентом упругости k.

Изображение1.png

Файл:Изображение2.png
Если учитывать взаимодействие только соседних атомов, уравнение движения можно записать в следующем виде:
Файл:Изображение3.png


Для решения уравнения движения воспользуемся численными методами интегрирования:


а. Метод Верле

б. Метод Рунге-Кутта 4 порядка

Файл:Изображение4.png
Файл:Изображение5.png
Файл:Изображение6.pngФайл:Изображение7.pngФайл:Изображение8.png
Файл:Изображение9.png
Файл:Изображение10.png
Файл:Изображение11.png
Файл:Изображение12.png
Файл:Изображение13.png
Файл:Изображение14.png
Файл:Изображение15.png


Для каждого из методов реализуются 3 вида граничных условий:# Фиксированные граничные условия


Файл:Изображение16.png
Файл:Изображение17.png
Файл:Изображение18.png
Файл:Изображение19.png


  1. Свободные граничные условия


Файл:Изображение20.png
Файл:Изображение21.png


  1. Периодические граничные условия


Файл:Изображение22.png
Файл:Изображение23.png


  1. Частица в потенциальной яме Леннарда-Джонса


Уравнение движения частицы в потенциальной яме Леннарда-Джонса:

Файл:Изображение24.png
Файл:Изображение25.png

Файл:Изображение26.png

Скоростью диссоциации будем называть скорость, которую необходимо сообщить частице, чтобы она улетела на бесконечность.

Файл:Изображение27.png

Решение:

Вывод:# Были реализованы различные методы интегрирования уравнения движения одномерной линейной цепочки. Заметим, что метод Верле является симплектическим и сохраняет энергию, в то время как метод Рунге-Кутта энергию не сохраняет.

  1. Была численно найдена скорость диссоциации частицы в потенциальной яме Леннарда-Джонса, которая с заданной точностью совпала с теоретическим значением.