168753 (595691), страница 9
Текст из файла (страница 9)
(3.96)
(3.97)
(3.98)
де D - коефіцієнт конвективної дифузії, м3/доб; c,N - концентрації дифундуючої речовини відповідно в рідинній і твердій фазах (кг/м3 ); vx(x, y, t) , vy(x, y, t) - координати швидкості фільтрації v, м/доб; - пористість або активна пористість ґрунту, у якому здійснюється рух вод і конвективна дифузія розчиненої речовини; t - час, доба;
-стала масообміну (швидкості сорбції);
- коефіцієнт розподілу речовини між; рідинною і твердою фазами в умовах рівноваги між; рідинною і твердою фазами за законом лінійної ізотреми Генрі, який виражається рівністю cp =
N , cp - рівновагова концентрація розчину, яка за величиною дорівнює кількості речовини, що поглинається твердою фазою;
- потенціал швидкості фільтрації;
- коефіцієнт фільтрації. м/доба;
- напір, м; p- тиск, Н/м2=кг/м·c2 ;
- щільність, кг/м3 ;
-прискорення сили тяжіння, м/c2 .
Будемо розглядати конвективну дифузію тих розчинних речовин, які нейтральні до порід, що наявні в ґрунті, тобто сорбцією та іншими видами поглинання речовин, що забруднюють підземні води, будемо нехтувати й розглядати систему рівнянь фільтрації та конвективної дифузії (гідравлічної дисперсії):
(3.99)
(3. 100)
рис. 3.3.
Будемо вважати, що розв'язки фільтраційних завдань для кожної конкретної схеми (мал. 3.3 -3.5) відомі, а також; відомі для цих схем відповідні області комплексного потенціалу (3.6), . Знайдемо розв'язки різних крайових задач для рівняння (3.100).
-
Крайові й початкові умови для шуканої функції с(х, у, t) :
При конвективній дифузії речовин, що забруднюють підземні води, на вході АВ фільтраційного потоку (3.3,3.4) можна прийняти одну із наступних крайових умов:
а)задана концентрація розчиненої у водоймі (річці) речовини
(3. 101)
б)задана умова Данкверста, яка враховує як конвективний, так й дифузійний механізми відведення речовини на водопроникненій ділянці межі області фільтрації:
(3. 102)
рис. 3.4.
рис. 3.5.
де n - нормаль до межі; vn - нормальна складова швидкості фільтрації.
На водонепроникних ділянках межі області фільтрації z та на кривих депресії виконується умова (наприклад, ділянка BC на мал. 3.3,а-г)
(3. 103)
На ділянці виходу фільтраційного потоку (на мал. (3.4) (11.18) - водопроникнені ділянки CD) можна прийняти одну із наступних крайових умов:
а) задана концентрація дифундуючої речовини або задана умова Данкверста(такі умови приймаються, коли не спостерігається інтенсивний відвід вод на виході фільтраційного потоку)
(3. 104)
б) задана умова, яка враховує тільки конвективне перенесення через межу (приймається в разі інтенсивного відведення вод на виході фільтраційного потоку; наприклад,
рис. 3.6.
вихід у дренаж CD на мал. 3.3б-г):
(3. 105)
Умови (3.103) і (3.105) начебто рівнозначні. Алі це не так, бо ці умови треба враховувати разом із межовими умовами для рівнянь фільтрації підземних вод, які різні для водопроникнених і водонепроникних ділянок межі області фільтрації (у першому випадку vn = 0, в іншому vn 0).
При конвективній дифузії солей і гіпсів, що залягають на визначеній глибині T фільтраційного потоку, на межі із сіллю або гіпсом зазвичай приймається умова
(3. 106)
де c - концентрація повної насиченості солі або гіпсу.
Початкові умови, що приймаються при розв'язанні задач про забруднення та засолювання підземних вод, мають вигляд
(3. 107)
або
(3. 108)
де c0 - задана концентрація дифундуючої речовини в області фільтрації в момент години до настання процесу. Складність, що виникає при розв'язанні стаціонарних і нестаціонарних крайових завдань, які описують двовимірні процеси, пов'язана не тільки з виглядом рівнянь у частинних похідних й виглядом крайових умов, але головним чином, із виглядом (геометрією) області, у якій відшукується розв'язок.
Тож, у рівняннях конвективної дифузії та наведених вище крайових умовах перейдемо до нових незалежних змінних - координат області комплексного потенціалу , які, як відомо, мають вигляд многокутника зі сторонами, паралельними прямокутній системі координат.
Нехай відома характеристична функція течії
, (3.109)
яку можна знайти, наприклад, методом конформних відображень. Тоді, виконуючи в рівнянні конвективної дифузії (3.100) заміну змінних і
отримаємо рівняння:
(3. 110)
де
(3. 111)
Після заміни змінних у крайових умовах через підстановку (3.109) отримуємо таке: межові умови (3.101), (3.102) набудуть вигляду відповідно
(3. 112)
(3. 113)
межову умову (3.103) перепишемо у вигляді
(3. 114)
межові умови (3.104), (3.105) набудуть вигляду відповідно
(3. 115)
(3. 116)
(3. 117)
(3. 118)
Межові умови (3.101)-(3.106) і відповідні до них умови (3.112)-(3.118) справджуються як для нестаціонарних, так і для стаціонарних крайових задач. Початкові умови (3.107), (3.108) перетворюються на умови
(3. 119)
(3. 120)
Перейдемо до розгляду завдань конвективної дифузії, розв'язання яких для різних схем фільтрації (мал. 3.5в) будемо шукати в областях комплексного потенціалу, які зображуються у вигляді прямокутників (мал. 3.5г-є).
При заданій концентрації розчинної речовини на межі з водоймами виникає така крайова задача: у прямокутнику ABCD знайти розв'язок рівняння (3.110). що задовольняє межові умови
(3. 121)
(3. 122)
(3. 123)
У разі усталеної конвективної дифузії отримуємо таку крайову задачу:
(3. 124)
(3. 125)
розв'язок якої, вочевидь, не залежить від змінної і має вигляд:
(3. 126)
Якщо враховувати механізм дифузійного відводу речовини на вході фільтраційного потоку, то отримаємо крайову задачу
(3. 127)
(3. 128)
розв'язок якої можна записати у вигляді
(3. 129)
Осереднюючи величину , що є в правій частині рівняння (3.110), по області приведеного комплексного потенціала
і замінюючи її деякою середньою величиною
, розглянемо типи двох нестаціонарних крайових завдань.
-
Перший тип крайових задач
Виникає при фільтрації забруднених вод у відкриті водойми, коли у відкритих водоймах підтримується задана концентрація речовин. Ці задачі формулюються таким чином:
треба знайти розв'язок рівняння
(3. 130)
що задовольняє межові умови вигляду (перша задача)
(3. 131)
або умови, що враховують механізм дифузійного відводу речовини від межі на вході фільтраційної течії (друга задача):
(3. 132)
і початкову умову
(3. 133)
Безпосередньою перевіркою легко переконатися, що розв'язком крайової задачі (3.130), (3.131),(3.133) та (3.130),(3.132),(3.133) будуть функції і
, які є розв'язки відповідних одновимірних крайових задач:
(3. 134)
(3. 135)
(3. 136)
(3. 137)
Якщо підставити в ці рівняння розв'язок у вигляді суми розв'язків стаціонарної і нестаціонарної задач і застосувати метод Фур'є, отримаємо розв'язки нестаціонарних задач конвективної дифузії (3.134) -(3.137), які після ділення на c1 і запровадження безрозмірних змінних та
набувають вигляд:
(3. 138)
(3. 139)
де власні значення визначаються рівняннями
(3. 140)
(3. 141)
коефіцієнти обчислюються за формулами
(3. 142)
(3. 143)
a функції визначаються рівнянням
коли
(3.144)
коли
(3.145)
Графіки цих функцій наведено на мал. 3.7, 3.8.
рис. 3.7.
-
Другий тип крайових задач
Другий тип крайових задач конвективної дифузії забруднюючих підземні води речовин характеризується межовою умовою вигляду (3.105), яку приймається на виході фільтраційного потоку, коли має місце інтенсивний відвід вод із дренажного каналу CD. У цьому разі розв'язком стаціонарних задач буде стала, значення якої залежить від межової умови на вході фільтраційного потоку. Тому будемо розглядати нестаціонарні задачі. Осереднюючи, як і раніше, швидкість фільтрації по просторовим змінним, приходимо до таких двох крайових завдань: треба знайти розв'язок рівняння
, (3.146)
що задовольняє межові умови
(3. 147)
а в разі враховування механізму дифузійного відводу речовини на вході фільтраційного потоку (друга крайова задача) треба знайти розв'язок рівняння
(3. 148)
яке задовольняє межові умови
(3. 149)
Застосування методу відокремлювання змінних (метод Фур'є) призводить до розв'язку (у частках від c1):
(3. 150)
рис. 3.8.
де , функція
визначається рівностями (3.144),(3.145), а коефіцієнти
обчислюються за формулою
(3. 151)
Розв'язок крайової задачі (3.148)-(3.149) отримуємо у вигляді:
(3. 152)
де , коефіцієнти
обчислюються за формулою
(3. 153)