Исмаилов Б. Р. Техногенді сипаттағы апаттарды модельдеу және бағалау пәнінен дәрістер жинағы


Турбуленттік диффузияның жартылай эмпирикалық теңдеуі анықтамасы



бет30/41
Дата25.04.2022
өлшемі1,37 Mb.
#140838
1   ...   26   27   28   29   30   31   32   33   ...   41
Байланысты:
лекция

2.Турбуленттік диффузияның жартылай эмпирикалық теңдеуі анықтамасы

Онда қандай да бір траекторияны бойлай қоспаның орнын ауыстыру процесі келесідей дифференциялдық теңдеумен сипатталады.


(1.1.1)
Осы теңдеуге салмақты сақтау заңының математикалық өрнегі болып саналатын беріктік теңдеуін қосу керек,
, (1.1.2)
мұндағы V=x, Vy, Vz> - қоспаның бөлшектерінің шапшаң жылдамдығының векторы.
Нәтижеде келесі формуланы аламыз
(1.1.3)
V векторын келесідей түрде көрсетуге болады
(1.1.4)
мұндағы - жеткілікті үлкен уақыт аралығы бойынша орташаланған қоспа бөлшегінің жылдамдығының векторы;
- V қоспа бөлшегінің жылдамдық векторының флуктуациялық құрастырушысының векторы.
Онда (1.1.3) формуланы келесідей етіп қайта жазуға болады:
(1.1.5)
(1.1.6)
(1.1.7)
(1.1.7) теңдеуіне келесідей бастапқы және шектік шарттарды қосу керек болады:
(1.1.8)
кезінде (1.1.9)
(1.1.7)-(1.1.9) моделі қоспаның таралу процесін жеткілікті түрде сипаттайды, оның практика жүзінде қолданылуы көптеген қосымша факторларға байланысты болған V векторы компонентін нақты анықтаудың мүмкінсіздігінің есебінен күрделілік тудырады. Алайда, турбуленттік атмосферада қоспаның шашырау кезінде алар орны бар турбуленттік диффузия процесіне келесідей қатынастар сай келеді:
(1.1.10)
мұндағы Кх, Ку, Кz – сәйкесінше ОХ, ОҮ, OZ координаттық осьтеріне сәйкес турбуленттік диффузияның коэффиценттері.
Осы коэффиценттер эксперименталдық мәліметтерден шыға келе анықталады немесе қоспаның таралуының кері есебінің шешімінің негізінде анықталады.
Сәйкесінше (1.1.7) теңдеу келесідей түрге ие болады:
(1.1.11)
Жазба ыңғайлы болуы үшін Vx, Vy, Vz үстіндегі сызықшаларды алып тастаймыз:
(1.1.12)
Ары қарай жел жылдамдығының бағыты ОХ осіне сай келеді деп болжанады, ал қоспа бөлшектерінің тұну жылдамдығының векторы төмен қарай вертикал бағытталған деп болжанады:
(1.1.13)
Осыдан басқа, егер қоспа ыдырайтын болса немесе басқа заттармен химиялық реакцияға түсетін болса, мұнда α – қоспаның химиялық белсенділігін сипаттайтын коэффицент ( мұндағы Те – уақыт аралығы, осы аралықта қоспа е есе ыдырайды), атмосферада f (t, x, y, z) функцияларымен сипатталатын қоспа көздері болса, онда (1.1.13) теңдеуі келесі түрге түрленеді:
(1.1.14)
(1.1.14) теңдеу турбуленттік диффузияның жартылай эмпирикалық теңдеуі деп аталады және атмосферада белсенді қоспаның таралу мәсеселсін шешу кезінде қолданбалы маңызды мәнге ие болады.
(1.1.14) теңдеу шешімін қанағаттандыруы керек болатын шектік және бастапқы шартты беру кезінде ұйғарымдардан шыға келе қоспаның таралуының t уақыты аралығына тиісті екенін көруге болады, ал осы қоспа таралатын Ω аймағы келесі түрдегі жартылай кеңістік болып саналады

мұндағы z0 = const>0 – төсеніштік беттің кедір бұдырлығының деңгейі (көбінесе кедір бұдырлық деңгейін елемейді және z0=0 деп ұйғарады) [10].
Осы болжамнан шыға келе бастапқы және шектік шартты келесідей түрде беруге болады:
(1.1.15)
(1.1.16)
(1.1.17)
мұндағы
- Ω фондық концентрацияны сипаттайтын қандай да бір функция;
vs – төсеніштік бетке қоспаның құрғақ тұнуының жылдамдығы.
(1.1.17) шектік шарты мынаны білдіреді, яғни турбуленттік және графитациялық Vzq құрауыштардан тұратын жер бетіндегі Р қоспа ағыны vsq тең.
Практикалық мәселелерді қарастыру кезіндегі (1.1.17) шектік шартты қолданудың күрделілігі мынадай, яғни құрғақ тұну vs жылдамдығының мәніне қатысты мәліметтер аз. Сондықтан көбінесе қарапайым (1.1.18) және (1.1.19) шарттарын қарастыру арқылы шектеледі.
Төсеніштік бетпен қоспаны толық жұту жағдайы үшін келесі шектік шарт орынды
z=z0 кезінде q=0(1.1.18)
Бұл жағдай сулы беттің үстінде қоспаның таралуы кезінде ерекше маңызға ие болады.
Төсеніштік бетте қоспаның толық таралуын көрсету үшін басқа шектік шарт қолданылады
z=z0 кезінде (1.1.19)
Бұл жағдай топырақ бетінде қоспаның таралуы кезінде ерекше маңызға ие болады.
(1.1.14) теңдеудегі Vx және Vz коэффициенттері әр түрлі беріледі. Жел жылдамдығын сипаттайтын Vx коэффициентін көбінесе келесідей түрде береді:
(1.1.20)
мұндағы
z1 – өлшемдік биіктік (көбінесе z1=1 деп алады);
Vx1 – z1 биіктігіндегі жел жылдамдығы.


Достарыңызбен бөлісу:
1   ...   26   27   28   29   30   31   32   33   ...   41




©engime.org 2024
әкімшілігінің қараңыз

    Басты бет