Стаховська А. І. Особливості клітинно-автоматної моделі потоку рідини FHP-M1 // Міжнародний науковий журнал "Інтернаука". — 2017. — №8.
Стаховська Анастасія Ігорівна
Студентка групи БС-32
Факультету Біомедичної Інженерії
Національного Технічного Університету України
«Київський Політехнічний Інститут імені Ігоря Сікорського»,
ОСОБЛИВОСТІ КЛІТИННО-АВТОМАТНОЇ МОДЕЛІ ПОТОКУ РІДИНИ FHP-M1
Клітинний автомат (КА) був запропонований фон-Нейманом у середині минулого століття[1]. На основі цієї моделі підтверджувалась думка про те, що людина може створити пристрій, якому присутні його властивості, перед усім, здатність створювати собі подібного. Тому КА довгий час розглядався як модель самовідтворення і пояснювався як спрощена модель деякої біологічної структури, яка складається із множини клітин. Кожній клітини ставиться у відповідність кінцевій автомат, який називається елементарним автоматом.
Він може знаходиться в одному зі станів: 0 або 1, та змінювати цей стан в залежності від стану клітин деякого особистого оточення, яке називається сусідством. Алгоритм обчислення наступного стану в залежності від стану сусідів у всіх клітин однаковий. Всі клітини виконують перехід у новий стан одночасно, тобто синхронно та паралельно. При цьому спостерігається зміна глобальної чорно - білої картини розподілу станів у просторі КА. Така картина називається конфігурацією КА.
Ітеративна зміна конфігурації при переходах всіх елементарних станів у нові стани називається еволюцією КА. Еволюціонуючи, КА моделює просторову динаміку, яка може мати завершення, повторюватися періодично або змінюватися хаотично. Дослідження цієї "моделі світу" показали, що, незважаючи на простоту кожної клітини, їх кооперативна робота моделює дуже складні і різноманітні процеси, які іноді неможливо (або, принаймні, невідомо як) описати іншим способом.
В середині 80-х років минулого століття стався сплеск інтересу до КА, у зв'язку з побудовою нових моделей обчислень. Були запропоновані КА, еволюція яких моделює процеси дифузії, поділу фаз, реакційно-дифузійні процеси, знамениту реакцію Белоусова-Жаботинського, руху солітонів, освіта дисипативних структур. Перехід від булевих просторів до звичних безперервним просторовим функціям (макровелічінам) проводиться шляхом осереднення значень станів клітин за деякою заданою околиці.
Справжнім потрясінням основ в моделюванні просторової динаміки стала поява клітинно-автоматної газової динаміки, названої "Gas-Lattice", що в перекладі на російську звучить як "гратковий газ". Дуже важливим є той факт, що було суворо математично доведено відповідність Gas-Lattice моделей рівняння Нав'є-Стокса [2].
Передбачається, що цей новий підхід до моделювання природних явищ стане доповненням до традиційних моделей математичної фізики, заснованим на диференціальних рівняннях в приватних похідних і добре розвиненим чисельних методів їх вирішення на комп'ютерах, обчислювальна потужність і архітектура яких змінюється швидше, ніж створюється необхідне математичне забезпечення.
Під клітинним автоматом КА моделі FHP-I будемо розуміти трійку елементів (W, A, N), де
(1)
Де,W – це множина клітин, задана їх координатами в деякому дискретному просторі.
Кожній клітині поставленний у відповідність кінцевий автомат А, який називається елементарним автоматом[3].
Внутрішнім станом автомату є булевий ветор. Кожній клітині
Зіставлені деякі координати x(w) та y(w) на декартовій площині. Тому, між будь-якими двома клітинами та можливо підрахувати відстань d(). Для кожної клітини визначена певна впорядкована множина:
(2)
Елементи якого перебувають у відношенні сусідства з кліткою w і називаються її сусідніми клітинами, або сусідами. Константа b визначає кількість нетотожні сусідів кожної клітини. Кожна клітина є сусідом сама собі.
Вхідний стан елементарного автомату A в клітці поставлено у відповідність внутрішнім станам сусідів цієї клітини. Таким чином, структура безлічі клітин W клітинного автомата представляється графом, в якому вершинами є клітини, а ребра відповідають відношенню сусідства. Цей граф має регулярну структуру і ступеня вершин рівні b.
Стан клітини представлено вектором s(w) з булевими компонентами. Безліч станів s(w) всіх клітин в один і той же момент часу t називається глобальним станом клітинного автомата КА На рисунку 1.2.1 зображена клітина w, вектори швидкостів ній частинок і її сусіди.
Рисунок 1. Вектори швидкостей частин. Нумерація сусідів
Таким чином, кількість сусідів кожної клітини w моделі FHP-Ι рівно семи, одним із сусідів є сама клітина w, тобто .
Стан клітини w в кожен момент дискретного часу t однозначно визначається набором частинок, що в ній знаходяться. Вектор швидкості кожної з них або спрямований в бік однієї з сусідніх клітин (при i=1,…,6), або дорівнює нулю (при i=0).
Таким чином, маса частинок в клітці w дорівнює:
(3)
Де: b = 6 - кількість можливих напрямків вектора швидкості; - i-й компонент вектора станів s.
Фізична інтерпретація значень компонентів вектора наступна: визначає частку з одиничною масою в клітці w, вектор швидкості якої спрямований в бік сусіда.
На рисунку 1.2.2 зображений вектор стану клітини.
Рисунок 2. Вектор стану клітини
Модельний імпульс p в клітині є сума всіх імпульсів, спрямованих до сусідів де i=1,..,b, а b=6:
(4)
Використовуючи (3) та рисунок 1.2.1, достатньо просто підрахувати проекції та імпульсу p на декартовій площині осі Ox та Oy:
(6)
Де – частица у клітині w, з вектором швидкості, направленим до сусіда
Список літератури: