Исмаилов Б. Р. Техногенді сипаттағы апаттарды модельдеу және бағалау пәнінен дәрістер жинағы
Атмосферадағы қозғалыс теңдеуін сандық тәсілмен шешу
жүктеу/скачать 1,37 Mb.
|
лекция
| 1. Атмосферадағы қозғалыс теңдеуін сандық тәсілмен шешу
Көп жағдайда қоспалардың атмосферада таралу процесіне модел жасау және шешу кезінде жорамалданғаннан елеулі дәрежеде көп болатын жағдайлар әсер етеді. Мұндай жағдайларда атмосферадағы белсенді қоспаның концентрациясының орташа мәндерін анықтау үшін бастапқы және шекті шарттары берілген (1.1)-(1.2) турбулентті диффузияның жартылай эмпирикалық теңдеуін қолданған мақсатқа сай келеді. Бұл параграфта автор ұсынып отырған қоспаның атмосферадағы қозғалыс теңдеуін және атмосферадағы қоспа диффузиясы теңдеулерін сандық жолмен шешуге негізделген турбуленттік диффузияның жартылай эмпирикалық теңдеуін сандық жолмен шешу тәсілі баяндалады [2]. Бұл (1.1теңдеулерін сандық шешу кезінде қоспа төселетін бет қабат арқылы толық жұтылады деп жорамалдаймыз.Сәйкесінше шекті жағдай шарты (1.1) мына төмендегі шекті жағдай шартына алмастырылады g(t,x,y,z0) = 0 (1.1) Сондай- ақ қоспаның ауада отырып қалу жылдамдығын сипаттайтын Vx коэффициенті турбулентті диффузияның (1.1) теңдеңіне «минус» емес «плюс» таңбасымен кіреді деп санаймыз. Мұндай алмастыру ары қарайғы есептеулерді жеңілдету мақсатында жасалады және координаталар жүйесін «аударып бұрудың» қажеті болмайды. Енді (1.1) теңдеуін [шекті жағдай шарты өзгертілген] сандық шешудің тиімді тәсілін көрсетейік.Бұл жерде Е+3қадамдарыы OX, OY, OZ остері бойында сәйкесінше ∆x, ∆y, ∆z болатын біркелкі (xi, yj, zk), (i, j, k € Z) торымен жабамыз және (1.1) теңдеуінің [шекті жағдай шарты (1.1) өзгертілген] сандық шешуін тек осы тор түйіндеріндегі жеткілікті аз болатын уақыт аралында [t,t + ∆t] іздейміз. [4] сәйкес (1.1) есебін [шекті жағдай шарты (1.1) өзгертілген] жеткілікті аз болатын уақыт аралында [t,t + ∆t] үш есепке бөліктеуге болады: ауыстыру қозғалысы есебі [4]: (.1.2) (1.3) кезінде (1.4) диффузия есебі [10,46]: (1.5) (1.6) кезінде (1.7) таралу және экологиялық мәні бар аймқта қоспаларды әкеліп кіргізу есебі: (1.8) (1.9) (1.2)-(1.4) есебін шешуге атмосферадағы қоспалардың қозғалысы есебін сандық жолмен шешуге ұқсас тәсілді пайдаланамыз. (1.2)-(1.4) үлгісінде Vx = const, Vy = const болғандықтан, t уақыты мезетіндегі қоспа бөлшектерінің (xi, yj, zk) нүктесіндегі жылдамдығы t +∆t уақыты мезетіндегі (xsx, ysy, zsz) нүктесіндегі қоспа бөлшектерінің жылдамдығына тең болады (1.10) Сонда (1.11) Мейлі, v0 нүктесі (xi, yj, zk) нүктесімен, ал v∆v нүктесі (xsx, ysy, zsz) нүктесімен тура келетін болсын. Бұл жерде Vx = |V| cos α, Vy = |V| cos π/2,Vz = |V| cos (π/2-α) болғандықтан, мұнда cos α, cos π/2, cos (π/2-α) – қоспа массасы ауысуының v бағытын анықтап беретін V векторының бағыттаушы косинустары, онда (1.2) өрнегінен ∆t аз уақыт аралығы үшін төмендегілерді аламыз (1.12) (1.13) (1.14) Енді (1.14) өрнегі үшін нақты айырмалану үлгісін жазамыз және онша күрделі емес түрлендірулерден соң мыналарды аламыз , (1.15) (1.16) (1.17) (1.18) (1.19) (1.20) Бұдан келіп шығатыны (1.21) (1.18) теңдеуінен ∆t · |V| / ∆v = 1 қатынасы орындалатыны келіп шығады, сондықтан (1.20) үлгісі ∆v және ∆t бойынша аппроксимациялардың шексіз қатарының үлгісі болып табылады, демек (1.21) үлгісі де Sx, Sy, Sz және ∆t бойынша аппроксимациялардың шексіз қатарының үлгісі болып саналады [6]. Бұл жерде біздер (1.2)-(1.4) есебінің шешімін тек тордың тораптық нүктелеріне қатысты іздеуді шарт еткендіктен (xsx, ysy, zsz) нүктелерінің координаталары мына қатынастар талаптарын қанағаттандырулары тиіс (1.22) мұнда ∆i, ∆j, ∆k – бүтін сандар. Бұдан (1.23) (1.23)-тен келіп шығатыны, ∆i және ∆k есептеу кезінде дөңгелектеу әсерінен қосымша есептеулік дәлсіздіктер туындайды, олардың ықпал ету дәрежесінің жоғарылығы соншалықты ∆x пен ∆z үлкен шамаларында және Vx пен Vz кіші мәндерінде (1.21) формуласы есептеулер жүргізу үшін мүлдем жарамсыз болып қалады. жүктеу/скачать 1,37 Mb. Достарыңызбен бөлісу:
|