Конформні перетворення й моделювання масо переносу

3.5.     Конформні перетворення й моделювання масо переносу

Процес міграції розчинних речовин при фільтрації підземних вод, як відомо, описується системою рівнянь:

(3.62)

; (3.63)

(3.64)

 або в скалярній формі

 (3.65)

(3.66)

 (3.67)

де v = {v_x,v_y,v_z} - вектор швидкості фільтрації, м/доб; (x,y,z,t) - потенціал фільтрації; c(x, y, z, t) і N(x, y, z, t) - концентрація дифундуючої речовини відповідно в рідинній і твердій фазах, г/см , D - коефіцієнт конвективної дифузії, м³/доб; σ - активна (або ефективна) пористість середовища, у якому протікає фільтрація розчину; t - година (у добах);  - оператор Гамільтона; γ - стала швидкості масообміну; β - коефіцієнт розподілу речовини між; фазами в умовах рівноваги по лінійній ізтермі Генрі c_p = βN. У багатьох практичних завданнях можна обмежитися дослідженням процесу масопереносу розчинних у фільтраційному потоці речовин тільки на основі рівнянь, що описують конвективний процес, а саме:

(3.68)

(3.69)

причому масообмін визначається такою досить розповсюдженою залежністю:

 (3.70)

де - концентрація граничної насиченості.

Наведені рівняння описують, як правило, міграцію й фізичну трансформацію (сорбцію, десорбцію) консервативних водорозчинних речовин.

Якщо дослідити масоперенос при плоско-вертикальній і плановій усталеній або квазіусталеній фільтрації підземних вод, то для моделювання цього процеса доцільно застосувати конформне перетворення рівнянь масопереноса до криволінійних змінних - координатам точок області комплексного потенціала фільтрації.

У разі плоско-вертикальної (профільної) фільтрації рівняння рухові підземних вод запишуться у вигляді

(3.71)

де χ - коефіцієнт фільтрації, h - напір, який визначається рівністю

(3.72)

причому вісь Oy спрямовано вертикально вниз, p - тиск, ρ - щільність, g - прискорення сили тяжіння.

У разі планової напорної фільтрації відповідні рівняння записуються так:

Φ = -χTh;

(3.73)

а в разі планової безнапорної фільтрації

(3.74)

У рівняннях (3.73), (3.74) через T позначено потужність напірного водоносного шару: q - вектор питомої фільтраційної витрати (м²/доб), a h - напір, який у даному випадку визначається таким рівнянням:

(3.75)

де Z - вертикальна координата точки фільтраційного потоку.

Припущення, що для фільтраційних течій, що розглядаються, можна побудувати область комплексного потенціала ω =  + iψ , де ψ - функція течії, і що відома характеристична функція течії

z = x + iy = F(ω) = F₁( , ψ) + i₂( ,ψ), (3.76)

за допомогою заміни

x = F₁( ,ψ);  y = F₂( ,ψ) (3.77)

перетворимо рівняння конвективної дифузії до нових незалежних змінних  і ψ. У результаті такої заміни рівняння конвективної дифузії в разі плоско-вертикальної фільтрації запишеться у вигляді

(3.78)

у разі планової напорної фільтрації - у такому вигляді:

 (3.79)

а в разі планової безнапорної фільтрації перетворюється до вигляду

(3.80)

Якщо в рівняннях (3.78)-(3.80) D = 0, отримаємо рівняння конвективного масопереносу без враховування дифузійних процесів, що перетворені до нових змінних , ψ або Φ, Ψ або Φ*, Ψ* відповідно для випадків плоско-вертикальної, планової напорної і планової безнапорної фільтрації, а саме:

(3.81)

 (3.82)

(3.83)

Отже, у разі нехтування дифузійними процесами питання про визначення концентрації речовин, що забруднюють підземні води, зводиться до розв'язання відповідної фільтраційної задачі та одного з рівнянь (3.81)-(3.83) із одною додатковою (початковою) умовою, яка задається залежно від фізичної постановки задачі.

Важливою характеристикою при дослідженні процесу забруднення підземних вод є час, протягом якого в даній точці області  (або області z) концентрація розчинної речовини досягає визначеної величини. Крім того, виникає питання про визначення часу, протягом якого концентрація розчинної речовини досягає в даній точці максимального значення. Основні диференціальні рівняння, з яких визначаються ці характеристики, а також фронт просування речовини (домішку) у фільтарційному потоці будуть наведені нижче.

Нехай відома концентрація розчинної в фільтраційному потоці речовини як функції координат точок області комплексного потенціала й і години t . Тоді для кожного значення (моменту) часу t можна побудуввати поверхню розподілу концентрації відносно області комплексного потенціала , а отже, й відносно області фільтрації z . Цим самим для кожного моменту часу буде визначено значення концентрації речовини, що розповсюджується в підземних водах, у будь-якій точці області фільтрації або впродовж; лінії, зокрема, впродовж; будь-якої йз ліній чи течії еквіпотенціальних ліній.

Якщо ж припустити, що міграція речовини здійснюється зі сталою концентрацією, то час, протягом якого станеться забруднення визначеної частки області фільтрації, знайдемо таким чином. Нехай відома швидкість фільтрації v(x,y,t) і характеристична функція течії, що отримана у вигляді (3.77). Швидкість розповсюдження розчинної у фільтраційному потоці речовини U(x,y,t) у даному разі дорівнює дійсній швидкості руху підземних вод V(x,y,t), яка зв'язана зі швидкістю фільтрації v(x,y,t) співвідношенням

(3.84)

де через  позначена активна пористість ґрунту (породи). За миттєво протікаючих сорбційних процесах, що визначаються рівністю (3.70), активна пористість замінюється ефективною пористістю середовища, що визначається рівністю

(3.85)

Із (3.84) отримуємо

(3.86)

Після перетворення рівності (3.86) до нових незалежних змінних  та  маємо

(3.87)

Замість рівнянь (3.81)-(3.83) зручно розглядати рівняння

 (3.88)

де  - безрозмірні величини, причому . До рівняння (3.88) легко звести кожне з рівнянь (3.81)-(3.83). Дійсно, якщо в рівнянні (3.88) покласти

(3.89)

то отримаємо рівняння (3.81), якщо в рівнянні (3.88) покласти

(3.90)

то отримаємо рівняння (3.82), а якщо в рівнянні (3.88) покласти

(3.91)

то отримаємо рівняння (3.83).

Якщо розглядається квазістаціонарна фільтрація (фільтраційні характеристики залежатимуть й від години за незмінного розташування ліній потоку), то конвективний масоперенос описується рівнянням

 (3.92)

де  залежна від дослідження переносу при профільній, плановій напорній і плановій безнапорній фільтрації визначається відповідно рівностями

(3.93)

(3.94)

(3.95)

Щоб знайти частинний розв'язок рівняння (3.88), треба задати додаткову умову при = 0 або t = 0 , тобто розлянути задачу Коші. При цьому суттєвою є фізична інтерпретація незалежних координат  it. Тож під часом розв'язання конкретних задач Коші для рівняня (3.88) слід відокремлювати початково-часову та початково-просторову задачі. Перший тип завдань, як правило, виникає під час дослідження процесів очищення або розсолення підземних вод та засолених земель; другий тип завдань (початково-просторові) з'являється зазвичай під час дослідження процесів забруднення або засолення підземних вод та родючих земель.

Похожие работы

Моделювання забруднення ґрунтів

Скачать

35924

3

4

... що концентрація речовини, яка поглинається снігом, пропорційна наземній концентрації. Застосовуючи тепер метод прямих до рівняння (6) для розрахунку забруднення по однорідному профілю, отримаємо рівняння:  (8) 5. Моделювання забруднення ґрунту пестицидами Одним із найбільш важливих інтегральних показників, які відбивають кінцевий результат взаємодії пестицидів, середовища і зовнішніх ...

Вплив на атмосферу підприємств харчової промисловості та шляхи його зниження (на прикладі ВАТ "Жашківський маслозавод")

Скачать

73446

5

6

... параметрів очисного пристрою; проектування та вибір очисного пристрою або фільтра 3.2 Заходи по охороні атмосферного повітря на ВАТ "Жашківський маслозавод" Основними й найбільш дійовими методами боротьби з забрудненням атмосфери на підприємстві є екологічні, діє продумана система заохочувальних і заборонних заходів, які допомагають запобігти забрудненню. Впровадження підприємством певних ...