Збереження фізичної величини при використанні граничних умов Неймана, застосованих до рівняння адвекції-дифузії


25

Я не розумію різної поведінки рівняння адвекції-дифузії, коли я застосовую різні граничні умови. Моя мотивація - це моделювання реальної фізичної величини (щільності частинок) при дифузії та адвекції. Щільність частинок повинна зберігатися у внутрішніх приміщеннях, якщо вона не витікає з країв. За цією логікою, якщо я виконую умови кордону Неймана, кінці системи, такі як (зліва та праворуч), систему слід "закрити", тобто якщо потік на кордоні дорівнює нулю, то жодна частинка не може вийти.ϕx=0

Для всіх моделей, наведених нижче, я застосував дискретизацію Кранка-Ніколсона до рівняння адвекційно-дифузійної форми, і всі моделювання мають граничних умов. Однак для першого та останнього рядків матриці (рядки граничної умови) я дозволяю змінювати незалежно від внутрішнього значення. Це дає змогу кінцеві точки бути повністю неявними.ϕx=0β

Нижче я обговорюю 4 різні конфігурації, лише одна з них - це те, що я очікував. Наприкінці я обговорюю свою реалізацію.

Дифузія лише обмежує

Тут умови відхилення вимикаються, встановлюючи швидкість на нуль.

Лише дифузія, з β = 0,5 (Crank-Niscolson) у всіх точках

Тільки дифузія (межі Неймана з бета = 0,5)

Кількість не зберігається, як видно із зменшення області пульсу.

Дифузія лише з = 0,5 (Crank-Niscolson) у внутрішніх точках, і = 1 (повністю неявна) на межахβββ

Тільки дифузія (межі Неймана з бета = 0,5 для внутрішніх приміщень, бета = 1 повністю неявна) межі

Використовуючи повністю неявне рівняння меж, я досягаю того, чого очікую: жодні частинки не виходять . Ви можете бачити це, зберігаючи ділянку як дифузну частинку. Чому вибір в граничних точках повинен впливати на фізику ситуації? Це помилка чи очікується?β

Дифузія та адвекція

Коли додається термін адвекції, значення на кордонах, здається, не впливає на рішення. Однак для всіх випадків, коли межі здаються "відкритими", тобто частинки можуть вийти з меж. Чому це так?β

Адвекція та дифузія з = 0,5 (Crank-Niscolson) у всіх точкахβ

Адвекція-дифузія (межі Неймана з бета = 0,5)

Advection and Diffusion з = 0,5 (Crank-Niscolson) у внутрішніх точках, і = 1 (повна неявна) на межахβββ

Адвекція та дифузія (межі Неймана з бета = 0,5 для внутрішніх приміщень, бета = 1 повністю неявна) межі

Реалізація рівняння адвекції-дифузії

Починаючи з рівняння адвекційно-дифузійної,

ϕt=D2ϕx2+vϕx

Написання за допомогою Кранка-Нікольсона дає,

ϕjn+1ϕjnΔt=D[1β(Δx)2(ϕj1n2ϕjn+ϕj+1n)+β(Δx)2(ϕj1n+12ϕjn+1+ϕj+1n+1)]+v[1β2Δx(ϕj+1nϕj1n)+β2Δx(ϕj+1n+1ϕj1n+1)]

Зверніть увагу, що = 0,5 для Crank-Nicolson, = 1 для повністю неявного, і, = 0 для повністю явного.β ββββ

Для спрощення позначень давайте зробимо підстановку,

s=DΔt(Δx)2r=vΔt2Δx

і перемістіть відоме значення часової похідної праворуч,ϕjn

ϕjn+1=ϕjn+s(1β)(ϕj1n2ϕjn+ϕj+1n)+sβ(ϕj1n+12ϕjn+1+ϕj+1n+1)+r(1β)(ϕj+1nϕj1n)+rβ(ϕj+1n+1ϕj1n+1)

Факторинг термінів дає,ϕ

β(rs)ϕj1n+1+(1+2sβ)ϕjn+1β(s+r)ϕj+1n+1Aϕn+1=(1β)(sr)ϕj1n+(12s[1β])ϕjn+(1β)(s+r)ϕj+1nMϕn

яку ми можемо записати у матричній формі як де,Aϕn+1=Mϕn

A=(1+2sββ(s+r)0β(rs)1+2sββ(s+r)β(rs)1+2sββ(s+r)0β(rs)1+2sβ)

M=(12s(1β)(1β)(s+r)0(1β)(sr)12s(1β)(1β)(s+r)(1β)(sr)12s(1β)(1β)(s+r)0(1β)(sr)12s(1β))

Застосування граничних умов Неймана

NB знову працює над виведенням, я думаю, що я помітив помилку. Я припускав повністю неявну схему ( = 1) при написанні кінцевої різниці граничної умови. Якщо припустити схему Кранка-Нісколсона, тут складність стає надто великою, і я не міг вирішити отримані рівняння для усунення вузлів, які знаходяться за межами домену. Однак, можливо, можливо, є два рівняння з двома невідомими, але я не міг ним керувати. Це, ймовірно, пояснює різницю між першим та другим сюжетами вище. Я думаю, ми можемо зробити висновок, що справедливі лише графіки з = 0,5 в граничних точках.ββ

Якщо припустимо, що потік ліворуч відомий (припускаючи повністю неявну форму),

ϕ1n+1x=σL

Якщо писати це як різницю по центру, це дає

ϕ1n+1xϕ2n+1ϕ0n+12Δx=σL

Тому ϕ0n+1=ϕ2n+12ΔxσL

Зверніть увагу, що це вводить вузол який знаходиться поза межами проблеми. Цей вузол можна усунути, використовуючи друге рівняння. Ми можемо записати вузол якϕ0n+1j=1

β(rs)ϕ0n+1+(1+2sβ)ϕ1n+1β(s+r)ϕ2n+1=(1β)(sr)ϕj1n+(12s[1β])ϕjn+(1β)(s+r)ϕj+1n

Підстановка значення знайдене з граничної умови, дає такий результат для рядка = 1,ϕ0n+1j

(1+2sβ)ϕ1n+12sβϕ2n+1=(1β)(sr)ϕj1n+(12s[1β])ϕjn+(1β)(s+r)ϕj+1n+2β(rs)ΔxσL

Виконуючи ту саму процедуру для виходу кінцевого рядка (при = ),jJ

2sβϕJ1n+1+(1+2sβ)ϕJn+1=(1β)(sr)ϕJ1n+(12s(1β))ϕJn+2β(s+r)ΔxσR

Нарешті, робити граничні рядки неявними (встановлення = 1) дає,β

(1+2s)ϕ1n+12sϕ2n+1=ϕj1n+1ϕjn+2(rs)ΔxσL

2sϕJ1n+1+(1+2s)ϕJn+1=ϕJn+2(s+r)ΔxσR

Тому за граничних умов Неймана ми можемо записати матричне рівняння ,Aϕn+1=Mϕn+bN

де,

A=(1+2s2s0β(rs)1+2sββ(s+r)β(rs)1+2sββ(s+r)02s1+2s)

M=(100(1β)(sr)12s(1β)(1β)(s+r)(1β)(sr)12s(1β)(1β)(s+r)001)

bN=(2(rs)ΔxσL002(s+r)ΔxσR)T

Моє теперішнє розуміння

  • Я думаю, що різниця між першим та другим сюжетами пояснюється зазначаючи помилку, викладену вище.

  • Щодо збереження фізичної величини. Я вважаю, що причина полягає в тому, що, як зазначалося тут , рівняння адвекції у формі, яку я записав, не дозволяє розповсюджуватись у зворотному напрямку, тому хвиля просто проходить через навіть граничні умови нульового потоку . Моя початкова інтуїція щодо збереження застосовується лише тоді, коли термін адвекції дорівнює нулю (це рішення на ділянці №2, де збережена площа).

  • Навіть при граничних умовах Неймана з нульовим потоком маса все ще може залишити систему, це тому, що правильні граничні умови в цьому випадку є граничними умовами Робіна, в яких загальний потік. задається . Крім того, умова Нойнмана вказує, що маса не може покинути область шляхом дифузії , вона нічого не говорить про адвекцію. По суті, ми чуємо, що це закриті граничні умови дифузії та відкриті граничні умови до адвекції. Для отримання додаткової інформації дивіться відповідь тут: Реалізація градієнта нульового граничного стану в адвекційно-дифузійному рівнянніϕx=0j=Dϕx+vϕ=0.

Чи погодились би ви?


Здається, що граничні умови виконані неправильно. Чи можете ви показати нам, як ви наклали граничні умови?
Девід Кетчесон

Гаразд Я оновив реалізацію, і я думаю, що я помітив помилку щодо застосування = 0,5 лише в граничних рядках. Я оновив своє "поточне розуміння" внизу питання. Чи є у вас коментар? β
boyfarrell

Отже ... як виглядає розбірливість меж у випадку кордонів Робіна? Ви показали це для меж Неймана, але не для меж Робіна.

Відповіді:


15

Я думаю, що одна з ваших проблем полягає в тому, що (як ви зауважили у своїх коментарях) умови Ноймана - це не умови, яких ви шукаєте , в тому сенсі, що вони не передбачають збереження вашої кількості. Щоб знайти правильну умову, перепишіть PDE як

ϕt=x(Dϕx+vϕ)+S(x,t).

Тепер термін, який відображається в дужках, - це загальний потік, і це кількість, яку потрібно поставити до нуля на межах, щоб зберегти . (Я додав заради загальності та для ваших коментарів.) Прикордонні умови, які вам доведеться встановити, - це (припустимо, що ваш космічний домен )Dϕx+vϕ=0ϕS(x,t)(10,10)

Dϕx(10)+vϕ(10)=0

для лівої сторони та

Dϕx(10)+vϕ(10)=0

для правого боку. Це так звана гранична умова Робіна (зауважте, що Вікіпедія прямо говорить, що це ізоляційні умови для рівнянь адвекційно-дифузійної).

Якщо встановити ці граничні умови, ви отримаєте властивості збереження, які шукали. Дійсно, інтегруючись у космічну область, ми маємо

ϕtdx=x(Dϕx+vϕ)dx+S(x,t)dx

Використовуючи інтеграцію за частинами з правого боку, ми маємо

ϕtdx=(Dϕx+vϕ)(10)(Dϕx+vϕ)(10)+S(x,t)dx

Тепер два центральних терміни зникають завдяки граничним умовам. Інтегруючись у часі, ми отримуємо

0Tϕtdxdt=0TS(x,t)dxdt

і якщо нам дозволено перемикати перші два інтеграли ,

ϕ(x,T)dxϕ(x,0)dx=0TS(x,t)dx

Це показує, що домен ізольовані від зовнішності. Зокрема, якщо , отримаємо збереження .S=0ϕ


Тепер я розумію, чому це спрацювало лише тоді, коли = 0; тому що це означатиме збереження за вашим підходом вище. Яким був би наслідком використання цієї граничної умови на вищезгаданому, чи відобразила б хвиля? Я думав, що це буде неможливо, оскільки в рівнянні немає нічого, що могло б дати мені негативну швидкість? v
boyfarrell

Найкращий спосіб дізнатися - це, мабуть, спробувати! Але якщо це поводиться правильно (і IMO це робить), ви повинні побачити певну кількість яка починає накопичуватися в лівій частині домену: рекламне поштовх штовхає у тому напрямку, але межа закрита. Накопичення припиняється, коли дифузія достатньо велика, щоб збалансувати її. Так ні, не повинно бути відбитої хвилі. ϕϕ
Dr_Sam

@DrSam Просто запитання щодо реалізації. Я розумію вашу думку щодо того, як зробити нуль кількості ліворуч. Але коли ти кажеш "праворуч лише граничний термін", що це означає? Я думав, що граничні умови повинні бути або Нейманом, або Діріхле (або сумішшю обох)?
boyfarrell

@boyfarrell У відповіді ліворуч / праворуч посилалося на виведення правильних граничних умов, а не на спосіб їх реалізації (відредаговано для наочності). Умови Робіна - класичні умови, хоча там менш відомі, ніж Діріхлет та Нойман.
Dr_Sam

Отже, що стосується реалізації, чи вважаєте ви, що я повинен встановити граничні умови Робіна для обох кордонів? Також, якщо рівняння має термін реакції (наприклад, чи потрібно гранична умова також включати цей термін?
ϕt=x(Dϕx+vϕ)+S(x,t)
boyfarrell
Використовуючи наш веб-сайт, ви визнаєте, що прочитали та зрозуміли наші Політику щодо файлів cookie та Політику конфіденційності.
Licensed under cc by-sa 3.0 with attribution required.