Розглянемо конвективну дифузію тих розчинних речовин, які нейтральні до наявних у ґрунті породам, тобто надалі сорбцією й іншими видами поглинання забруднюючі підземні води компонентів знехтуємо й будемо виходити з наступної системи рівнянь фільтрації й конвективної дифузії (гідравлічної дисперсії):
При конвективній дифузії речовин, що забруднюють підземні води, на вході АВ фільтраційного потоку можна прийняти одне з наступних граничних умов:
а) задана концентрація розчиненого у водоймі (ріці) речовини
б) задана умова Данквертса, що враховує як конвективний, так і дифузійний механізми відводу речовини на водопроникній ділянці границі області фільтрації
де n- нормаль до границі; vn - нормальна складова швидкості фільтрації.
На водонепроникних ділянках границі області фільтрації й на криві депресії виконується умова
На ділянці виходу фільтраційного потоку (CD) можна приймати одне з наступних граничних умов:
а)задана концентрація речовини, що дифундує, або задана умова Данквертса (такі умови приймаються, якщо не спостерігається інтенсивного відводу вод на виході фільтраційного потоку)
б) задана умова, що враховує тільки конвективний перенос через границю (у випадку інтенсивного відводу вод на виході фільтраційного потоку)
При конвективній дифузії солей і гіпсів, що залягають на певній глибині T фільтраційного потоку, на границі із сіллю або гіпсом звичайно приймається умова
де
Початкові умови засолення підземних вод, мають вигляд
де c0 - задана концентрація речовини, що дифундує, в області фільтрації в момент часу до настання процесу забруднення (засолення) або промивання підземного середовища.
Трудність, що виникає при рішенні стаціонарних і нестаціонарних крайових завдань, що описують двовимірні процеси, зв'язана не тільки з видом рівнянь у частинних похідних і з видом граничних умов, а головним чином з видом (геометрією) області, у якій відшукується рішення. У зв'язку із цим у рівняннях конвективної дифузії й у наведених вище граничних умовах доцільно перейти до нових незалежних змінних координат комплексного потенціалу ω, що, як відомо, має вигляд багатокутника зі сторонами, паралельними прямокутній системі координат.
Нехай відома характеристична функція течії
яку можна знайти, наприклад, методом конформних відображень. Тоді, зробивши в рівнянні конвективної дифузії (1.49) заміну змінних
Взявши середню величину
Перший тип крайових завдань виникає при фільтрації забруднених вод у відкриті водойми (водоймища), коли в останні підтримується задана концентрація речовин. Ці задачі формулюються в такий спосіб: потрібно знайти рішення
задовольняючій або граничній умовам виду (перша задача)
або умовам, що враховують механізм дифузійного відводу речовини від границі на вході фільтраційної течії (друга задача):
і початковій умові
Безпосередньою перевіркою легко переконатися, що рішенням двовимірних крайових завдань (1.58), (1.59), (1.61) і (1.58), (1.60), (1.61) будуть функції
Підставляючи це рішення у вигляді суми рішень стаціонарного й нестаціонарного завдань і застосовуючи метод поділу змінних, одержимо рішення нестаціонарних завдань конвективної дифузії, які після розподілу на c1 і введення безрозмірних величин
де власні значення
Коефіцієнти
Другий тип крайових задач конвективної дифузії підземної води, речовин що забруднять, характеризується гарничною умовою, що приймається на виході фільтраційного потоку, коли спостерігається інтенсивний відвід із дренажного каналу CD. У цьому випадку рішенням стаціонарних задач буде стала, значення якої залежить від крайової умови на вході фільтраційного потоку.
Тому перейдемо до розгляду нестаціонарних завдань. Осереднюючи швидкість фільтрації по просторовим змінним, приходимо до наступних двох крайових завдань: Потрібно знайти рішення
задовольняючим крайовим умовам:
а у випадку обліку механізму дифузійного відводу речовини на вході фільтраційного потоку (друга крайова задача) потрібно знайти вирішення
задовольняючим крайовим умовам:
Застосування методу Фур'є до крайової задачі(1.72)-(1.73) дає вирішення
де
Коефіцієнти
Рішення крайової задачі (1.74)-(1.75) одержуємо в наступному виді:
де коефіцієнти
а власні значення λn визначаються з рівняння
λn =