Я не розумію різної поведінки рівняння адвекції-дифузії, коли я застосовую різні граничні умови. Моя мотивація - це моделювання реальної фізичної величини (щільності частинок) при дифузії та адвекції. Щільність частинок повинна зберігатися у внутрішніх приміщеннях, якщо вона не витікає з країв. За цією логікою, якщо я виконую умови кордону Неймана, кінці системи, такі як (зліва та праворуч), систему слід "закрити", тобто якщо потік на кордоні дорівнює нулю, то жодна частинка не може вийти.
Для всіх моделей, наведених нижче, я застосував дискретизацію Кранка-Ніколсона до рівняння адвекційно-дифузійної форми, і всі моделювання мають граничних умов. Однак для першого та останнього рядків матриці (рядки граничної умови) я дозволяю змінювати незалежно від внутрішнього значення. Це дає змогу кінцеві точки бути повністю неявними.
Нижче я обговорюю 4 різні конфігурації, лише одна з них - це те, що я очікував. Наприкінці я обговорюю свою реалізацію.
Дифузія лише обмежує
Тут умови відхилення вимикаються, встановлюючи швидкість на нуль.
Лише дифузія, з = 0,5 (Crank-Niscolson) у всіх точках
Кількість не зберігається, як видно із зменшення області пульсу.
Дифузія лише з = 0,5 (Crank-Niscolson) у внутрішніх точках, і = 1 (повністю неявна) на межахβ
Використовуючи повністю неявне рівняння меж, я досягаю того, чого очікую: жодні частинки не виходять . Ви можете бачити це, зберігаючи ділянку як дифузну частинку. Чому вибір в граничних точках повинен впливати на фізику ситуації? Це помилка чи очікується?
Дифузія та адвекція
Коли додається термін адвекції, значення на кордонах, здається, не впливає на рішення. Однак для всіх випадків, коли межі здаються "відкритими", тобто частинки можуть вийти з меж. Чому це так?
Адвекція та дифузія з = 0,5 (Crank-Niscolson) у всіх точках
Advection and Diffusion з = 0,5 (Crank-Niscolson) у внутрішніх точках, і = 1 (повна неявна) на межахβ
Реалізація рівняння адвекції-дифузії
Починаючи з рівняння адвекційно-дифузійної,
Написання за допомогою Кранка-Нікольсона дає,
Зверніть увагу, що = 0,5 для Crank-Nicolson, = 1 для повністю неявного, і, = 0 для повністю явного.β β
Для спрощення позначень давайте зробимо підстановку,
і перемістіть відоме значення часової похідної праворуч,
Факторинг термінів дає,
яку ми можемо записати у матричній формі як де,
Застосування граничних умов Неймана
NB знову працює над виведенням, я думаю, що я помітив помилку. Я припускав повністю неявну схему ( = 1) при написанні кінцевої різниці граничної умови. Якщо припустити схему Кранка-Нісколсона, тут складність стає надто великою, і я не міг вирішити отримані рівняння для усунення вузлів, які знаходяться за межами домену. Однак, можливо, можливо, є два рівняння з двома невідомими, але я не міг ним керувати. Це, ймовірно, пояснює різницю між першим та другим сюжетами вище. Я думаю, ми можемо зробити висновок, що справедливі лише графіки з = 0,5 в граничних точках.
Якщо припустимо, що потік ліворуч відомий (припускаючи повністю неявну форму),
Якщо писати це як різницю по центру, це дає
Тому
Зверніть увагу, що це вводить вузол який знаходиться поза межами проблеми. Цей вузол можна усунути, використовуючи друге рівняння. Ми можемо записати вузол як
Підстановка значення знайдене з граничної умови, дає такий результат для рядка = 1,
Виконуючи ту саму процедуру для виходу кінцевого рядка (при = ),
Нарешті, робити граничні рядки неявними (встановлення = 1) дає,
Тому за граничних умов Неймана ми можемо записати матричне рівняння ,
де,
Моє теперішнє розуміння
Я думаю, що різниця між першим та другим сюжетами пояснюється зазначаючи помилку, викладену вище.
Щодо збереження фізичної величини. Я вважаю, що причина полягає в тому, що, як зазначалося тут , рівняння адвекції у формі, яку я записав, не дозволяє розповсюджуватись у зворотному напрямку, тому хвиля просто проходить через навіть граничні умови нульового потоку . Моя початкова інтуїція щодо збереження застосовується лише тоді, коли термін адвекції дорівнює нулю (це рішення на ділянці №2, де збережена площа).Навіть при граничних умовах Неймана з нульовим потоком маса все ще може залишити систему, це тому, що правильні граничні умови в цьому випадку є граничними умовами Робіна, в яких загальний потік. задається . Крім того, умова Нойнмана вказує, що маса не може покинути область шляхом дифузії , вона нічого не говорить про адвекцію. По суті, ми чуємо, що це закриті граничні умови дифузії та відкриті граничні умови до адвекції. Для отримання додаткової інформації дивіться відповідь тут: Реалізація градієнта нульового граничного стану в адвекційно-дифузійному рівнянні.
Чи погодились би ви?