Покладіть Lх = Lу = 30 і rsса1е = 0.8. Чому дорівнює наведена щільність системи? Оцініть початкову щільність крапельки. Оскільки спочатку частинки ближче один до одного, необхідно вибрати менше значення? T, наприклад? T = 0.01. Чому повна енергія позитивна? Пройдіть принаймні 200 кроків за часом і обчисліть середні температуру і тиск за останніми 100 кроків. Запам'ятайте кінцеву конфігурацію, з тим щоб використовувати її в задачі 6.5. Який характер ваших знімків? Поясніть, чому дане рівноважний стан можна трактувати як газ. p align="justify"> ЗАВДАННЯ 2. Розподіл за швидкостями
а. В якості початкової конфігурації використовуйте для цього завдання кінцеву конфігурацію з завдання 6.4г. Наша мета-обчислити для рівноважного стану ймовірність Р (v)? V того, що швидкість частинки укладена між v і v +? V. Виходячи з початкової конфігурації, оцініть максимальну швидкість частинок. Виберіть інтервали шириною? V, при цьому номер інтервалу k дорівнює v /? V, де v - швидкість частки. Доцільно вибрати? V рівним 0.1 * sqr (Т), де Т-температура. Для реєстрації кількості влучень швидкостей часток у відповідний інтервал з номером k використовуйте масив prob (k): Припустимо, що? V = 0.1, і розглянемо швидкості v = 0.3, 0.49, 0.5, 0.51 і 0.9. Яким значенням k вони відповідають? Напишіть коротку програму для отримання масиву рrоb (k). Визначте рrоb (k) після кожного часового кроку і обчисліть середні значення рrоb (k) принаймні по 100 крокам за часом. Нормуйте рrоb (k), поділивши на число часток і кількість кроків за часом. Зауважимо, що спостерігати траєкторії часток немає необхідності. Надрукуйте таблицю, що містить значення k, v і рrоb (k). p align="justify"> б. Побудуйте графік щільності ймовірності Р (v) як функції від v. Який якісний вид Р (v)? Чому одно найбільш ймовірне значення v? Яка приблизно В«ширинаВ» Р (v)? Даний розподіл ймовірності називається розподілом Максвелла - Больцмана. p align="justify"> в. Визначте розподіл ймовірності для кожної компоненти швидкості. Переконайтеся, що ви розрізняєте позитивні і негативні швидкості. Які найбільш вірогідні значення для x-і у-компонент? Чому рівні середні значення? p align="justify"> ЗАВДАННЯ 3. Температурна залежність внутрішньої енергії
а. Однією з характерних рис методу молекулярної динаміки є те, що повна енергія визначається початковими умовами, а температура є величина похідна, обумовлена ​​тільки після, досягнення системою теплового рівноваги. В результаті важко вивчати систему, що знаходиться при конкретній температурі. Звичайно для досягнення необхідної температури Т f в якості початкової умови беруть рівноважну конфігурацію при температурі Т i , по можливості найменш відрізняється від Т f. Визначають В«масштабнийВ» множник f зі співв...
