1.Техногенді сипаттағы апаттарды модельдеуде қоспаның бұлт ішінде таралу есебін сандық шешу тәсілі
Белсенді қоспалардың бастапқы және шекаралық шарттары берілген дифференциалдық теңдеулер жиынтығы болып табылады. Бірақ тікелей есептеулерді іске асыру үшін оның аналитикалық немесе сандық шешімін табу қажет. Автор ұсынып отырған бұл есепті сандық шешу тәсілі күрделілігі жағынан ойға қонымды және осы параграфта сөз болады (1 суретті қараңыз). Байқайтынымыздай бұл параграфта келтірілген формулалардың көпшілігі белгіленуіне дейін алдыңғы зерттеме жұмыстарында келтірілгендерге сай келеді. Мұндай қайталану мұқтаждығы себепті және ерекше алып қарағанда жазылып отырған материалдың түсінікті болуын арттыру үшін қажет.
Сурет 1. Белсенді қоспаның бұлт ішінде таралуын сандық есептеудің жалпы сызбасы
Сәйкесінше OX, OY, OZ осьтерін бойлай Δx, Δy, Δz қадамы бар біркелкі (xi, yi, zi), (i, j, k Z) тормен жабамыз және ең бір кіші [t, t+ Δt] уақыттық аралықта осы тордың торабынан шешімін іздейміз. Мәселелерінің шешімі екі өзара байланысқан ішкі мәселелердің шешіміне орайластырылады – яғни бұлтты түзетін сулы будың тиаралу процесін сипаттайтын мәселелердің шешілуіне және бұлт ішінде белсенді қоспаның таралуы прцоесін сипаттайтын мәселелердің шешілуіне орайластырылады(сурет 1).
Ыдырау әдісін қолдана отырып, мәселені 3 ішкі мәселеге бөлшектейміз:
1) ауырлық күші мен жел әсерінен бұлт бөліктерінің ауысу мәселесі:
(1.1)
(1.2)
, (1.3)
кезінде (1.4)
2) бұлт түзетін сулы булардың диффузия мәселесі:
(1.5)
(1.6)
(1.7)
кезінде (1.8)
3) бұлттың сулы буын алу мәселесі:
(1.9)
(1.10)
[4] сай кіші уақыттық [t, t+Δt] аралығында осындай ыдырауға рұқсат етіледі және кезінде шешімдеріне ұқсас шешімге алып келеді.
Будың сулы буларын тасымалдаудың модельдерінде жел жылдамдығы және бұлттың сөну жылдамдығы тұрақты шамалар болып саналады, бұлттың сулы буының қозғалу жылдамдығы t уақыт сәтіндегі нүктесінде t+Δt уақыт сәтіндегі нүктесіндегі сулы будың қозғалыс жылдамдығына тең болады.
(1.11)
Сонда
(1.12)
Осы зерттеудлің 2 тарауы 6 параграфында келтірілген бос атомферада қопслардың бөлшектерін тасымалдауға арналған пікірге ұқсас пікір жүргізе отырып келесіні жазуға болады:
(1.13)
Тор торабтарында мәселелерді нақты шешу үшін келесі қатынастар орындалуы керек
(1.14)
мұндағы Δi, Δj, Δk – бүтін сандар.
Осы жерден
(1.15)
(1.16)
Δi, Δj, Δk шамаларын дөңгелектеу кезінде туындайтын қосымша есептік кемістікті жою мақсатында автор тордың бір нүктесінде емес, (2 суретті қараңыз) нүктесінен жойылуна пропорционал болатын нүктесін қоршайтын тордың дереу төрт нүктесінде концентрацияның мәнін келесі формула бойынша ауыстыруды ұсынады:
(1.17)
(1.18)
(1.19)
(1.20)
мұндағы
(1.21)
Достарыңызбен бөлісу: |