Spotlight #71: Симуляції з фізики, хімії та інженерії

Кожна симуляція на цій платформі починається з конкретного питання: як насправді виглядає це рівняння, коли ви його запускаєте? Нещодавно доповнені симуляції, представлені в цьому огляді, поділяють цю якість. Вони охоплюють фізичні явища, що трапляються в бакалаврських та магістерських програмах з фізики, хімії та інженерії, але рідко демонструються інтерактивно — зазвичай тому, що обчислювальна вартість історично була заборонною для браузерного контексту. WebGL та оптимізовані чисельні методи це змінили. Кожна з шести описаних тут симуляцій працює з інтерактивною швидкістю на будь-якому сучасному пристрої, без жодних плагінів.

I. Дифракція Брегга

Рентгенівська кристалографія — техніка, яка розкрила структуру ДНК, форми активних центрів білків та розташування атомів у кожному матеріалі, властивості якого ми розуміємо глибоко. Ця техніка спирається на формулу, настільки елегантну, що її вивів Вільям Лоуренс Брегг у віці 22 років:

nλ = 2d sinθ

де:
  n      : порядок відбиття (додатне ціле число)
  λ      : довжина хвилі рентгенівського випромінювання
  d      : відстань між площинами кристала (міжплощинна відстань)
  θ      : кут ковзного падіння

Для ГЦК:  відбиття дозволені, коли h, k, l усі непарні або всі парні
Для ОЦК:  відбиття дозволені, коли h + k + l парне
Проста кубічна: усі hkl дозволені

Симуляція обчислює порошковий дифракційний патерн для гранецентрованої кубічної (ГЦК), об'ємноцентрованої кубічної (ОЦК) та простої кубічної кристалічних структур, підсумовуючи внески всіх дозволених відбиттів Брегга до заданого користувачем максимального 2θ. Інтенсивності піків враховують структурний фактор, фактор Лоренца-поляризації та гаусове розширення для моделювання скінченного розміру кристалітів (рівняння Шеррера: β = Kλ / L cosθ). Результат відображається як стовпчикова діаграма, накладена на симуляцію кілець Дебая-Шеррера.

Найповчальніша риса — перемикач кристалічної структури. Перемкніть з ГЦК на ОЦК, і систематичні погашення одразу змінюються: ГЦК пригнічує відбиття зі змішаною парністю індексів Міллера (тож відбиття (1,0,0) і (1,1,0) зникають), тоді як ОЦК пригнічує відбиття з непарним h+k+l. Це робить зв'язок між геометрією елементарної комірки й дифракційним патерном безпосередньо відчутним, а не суто алгебраїчним. Повзунки керують довжиною хвилі λ та параметром решітки a, дозволяючи користувачам відтворювати опубліковані дифракційні патерни для реальних матеріалів, підбираючи положення піків.

II. Довжина Дебая

У плазмі вільні заряди перерозподіляються, щоб екранувати будь-яке електростатичне збурення на характерній відстані, що зветься довжиною Дебая. Формула має вигляд:

λD = √(ε0 kB T / n q²)

Екранований потенціал навколо пробного заряду:
  φ(r) = (q / 4πε0 r) exp(−r / λD)

Це потенціал Юкави: закон Кулона, помножений на
експоненційний спад на масштабі λD.

Симуляція розміщує пробний заряд у центрі 2D сітки плазми й рендерить результуючий електростатичний потенціал як теплову карту в псевдокольорах, зі шкалою кольору, що охоплює від сильного центрального поля до екранованої периферійної області. Повзунки керують електронною температурою T та густиною числа частинок n. Коли ви їх пересуваєте, теплова карта перемальовується в реальному часі, показуючи, як λD масштабується як √(T/n): підвищення температури збільшує довжину Дебая (гарячіші електрони можуть відлітати далі від зони екранування), тоді як підвищення густини її зменшує (більше електронів доступно, щоб екранувати збурення на коротшій відстані).

Супровідний радіальний графік у логарифмічному масштабі підтверджує експоненційний спад і показує перехід від неекранованого закону Кулона 1/r на коротких відстанях до експоненційно пригніченого хвоста Юкави на великих r. Саме це подвійне відображення — теплова карта й радіальний профіль одночасно — найчастіше запитувана функція серед викладачів фізики плазми, які використовують симуляцію.

III. Цикл Карно

Цикл Карно — теоретична верхня межа ефективності теплового двигуна. Жоден двигун, що працює між гарячим резервуаром при температурі TH і холодним резервуаром при TC, не може перевищити ефективність Карно:

ηКарно = 1 − TC / TH

Чотири етапи циклу на P-V діаграмі:
  1. Ізотермічне розширення при TH:  поглинається QH, газ виконує роботу
  2. Адіабатичне розширення:           Q = 0, температура падає до TC
  3. Ізотермічне стиснення при TC: віддається QC, робота виконується над газом
  4. Адіабатичне стиснення:         Q = 0, температура зростає до TH

Чиста робота = площа, обмежена циклом на P-V діаграмі
Ефективність η = W_чист / QH

Симуляція анімує робочий газ через усі чотири стадії на живій P-V діаграмі, з точкою стану, що прокреслює замкнений цикл під час роботи анімації. Повзунки керують TH і TC, а супровідна T-S діаграма показує цикл у просторі ентропії (де цикл Карно є прямокутником, що робить геометричну інтерпретацію ефективності одразу зрозумілою). Живий підрахунок показує поточні значення QH, QC, Wчист і η, оновлюючись з кожним циклом. Панель ентропії — та риса, яка робить цю симуляцію по-справжньому кориснішою за статичну діаграму з підручника: студенти можуть побачити, чому оборотність термодинамічно особлива.

IV. Реакція Бєлоусова-Жаботинського

Реакція Бєлоусова-Жаботинського (БЖ) — родина осцилюючих хімічних реакцій, що спонтанно генерують просторові патерни: концентричні кільця, обертові спіралі й біжучі хвилі в тому, що виглядає як добре перемішаний розчин. Механізм включає окисно-відновну реакцію, що осцилює між окисним і відновним станами, поєднану з автокаталітичним зворотним зв'язком. Модель Орегонатора зводить механізм до трьох речовин:

ОДР Орегонатора (спрощено):
  ∂u/∂t = Du∇²u + (1/ε)(u − u² − f v (u − q)/(u + q))
  ∂v/∂t = Dv∇²v + u − v

де:
  u : концентрація HBrO2 (автокаталізатор)
  v : концентрація Ce(IV) (окиснений каталізатор)
  ε : відношення масштабів часу (~0,04)
  f : стехіометричний параметр (~1,4)
  q : відношення констант швидкості (~0,002)

Симуляція інтегрує рівняння Орегонатора на 2D сітці за допомогою явної схеми скінченних різниць. Починаючи з випадкових початкових умов, система спонтанно розвиває спіральні хвильові патерни після короткого перехідного процесу — точно так, як спостерігається в реакції БЖ у чашці Петрі. Кольорове поле відображає концентрацію HBrO₂ у класичну синьо-червону осциляцію реакції БЖ, каталізованої церієм.

Найяскравіша риса — взаємодія спіральних хвиль. Дві спіралі, що обертаються в протилежних напрямках і наближаються одна до одної, анігілюють, тоді як спіраль, що стикається з межею, відбивається. Ці правила зіткнення суто емерджентні з локальної реакційно-дифузійної динаміки, без явного правила, що визначає, що відбувається на межах або коли хвилі зустрічаються. Ця емерджентність глобального порядку з суто локальних взаємодій — властивість, яка робить реакцію БЖ канонічним прикладом самоорганізації в хімії, далекій від рівноваги.

V. Принцип Бернуллі

Рівняння Бернуллі — одне з найуживаніших і найбільш неправильно зрозумілих рівнянь у вступній фізиці. Це не окремий закон аеродинаміки; це твердження про збереження для нев'язкого, нестисливого, стаціонарного потоку вздовж лінії течії:

P + ½ρv² + ρgh = константа вздовж лінії течії

Для горизонтального потоку (h = const):
  P1 + ½ρv1² = P2 + ½ρv2²

Рівняння нерозривності (нестисливий потік):
  A1 v1 = A2 v2

Тобто вужчий переріз означає вищу швидкість, що означає нижчий тиск.

Симуляція моделює потік крізь канал, поперечний профіль якого може малювати сам користувач. Поле швидкості розв'язується методом функції течії на 2D сітці, з лініями течії, відображеними як анімовані трасувальні шляхи. Поле тиску обчислюється зі швидкості через рівняння Бернуллі й показується як кольорове накладення, тож студенти одразу бачать падіння тиску в звуженій області. Режим витратоміра Вентурі додає індикатори відбору тиску в двох точках і показує перепад тиску та обчислену витрату — точно так, як працював би фізичний витратомір Вентурі.

Симуляція також демонструє обмеження Бернуллі: пересуньте повзунок числа Рейнольдса достатньо високо, і потік відривається від стінок каналу, створюючи зони рециркуляції, де рівняння Бернуллі перестає діяти. Ця чесна демонстрація припущень рівняння — стаціонарність, нев'язкість, нестисливість — робить симуляцію кориснішою для навчання критичного мислення, ніж версія, яка просто підтверджує формулу.

VI. Броунівський рух

Броунівський рух був однією з чотирьох статей Ейнштейна 1905 року, і саме він надав перший кількісний доказ атомної гіпотези. Пилок, що спостерігається під мікроскопом, прокреслює хаотичний, на вигляд випадковий шлях, тому що його штовхають окремі молекули води. Ейнштейн показав, що середньоквадратичне зміщення зростає лінійно з часом:

⟨r²(t)⟩ = 2dDt

де:
  d : просторова розмірність (2 для 2D руху)
  D : коефіцієнт дифузії = kBT / (6πηR)
  η : в'язкість рідини
  R : радіус частинки
  kB : стала Больцмана
  T : температура

Шлях — вінерівський процес: прирости Δr ~ N(0, 2DΔt)

Симуляція відображає кілька броунівських частинок одночасно, прокреслюючи їхні шляхи як згасаючі сліди на 2D полотні. Живий графік середньоквадратичного зміщення відносно часу зростає лінійно, підтверджуючи співвідношення Ейнштейна. Коефіцієнт дифузії визначається підбором нахилу до кривої СКЗ і порівнюється з теоретичним значенням Стокса-Ейнштейна для обраного розміру частинки й температури. Повзунки керують температурою, в'язкістю й радіусом частинки, а коефіцієнт дифузії безперервно оновлюється, роблячи залежності 1/R і T безпосередньо спостережуваними.

Опційний режим окремої частинки наближує одну траєкторію й показує розподіл зміщень на обраних користувачем інтервалах часу, підтверджуючи, що зміщення гаусові з дисперсією, що лінійно зростає — визначальна статистична ознака вінерівського процесу та математична основа всього сучасного стохастичного числення.

Спільні риси

Розглядаючи ці шість симуляцій разом, варто явно виділити кілька принципів дизайну, що повторюються. По-перше, кожна симуляція показує щонайменше два представлення одночасно: просторову візуалізацію та кількісний графік. Цей підхід із подвійним представленням означає, що користувачі можуть переходити між інтуїцією (патерном) і вимірюванням (кривою), не перемикаючись між різними інструментами чи сторінками.

По-друге, кожна симуляція має щонайменше один повзунок чи інтерактивний параметр, який змінює явище якісно, а не лише кількісно. Перемикання між ГЦК і ОЦК у симуляції Брегга не просто зсуває піки — воно змінює, які піки взагалі існують. Підвищення температури в симуляції броунівського руху не просто змушує частинки рухатися швидше — воно змінює коефіцієнт дифузії конкретним, вимірюваним чином. Саме в цих якісних переходах відбувається найглибше навчання.