Ток функциясы, құйын


Пән: Физика
Жұмыс түрі:  Курстық жұмыс
Тегін:  Антиплагиат
Көлемі: 31 бет
Таңдаулыға:   

РЕФЕРАТ

Бітіру жұмысында «Ток функциясы, құйын» айнымалыларындағы тұтқыр сығылмайтын сұйықтың қозғалысын сипаттайтын теңдеулер қарастырылған. Есептің қойылымында сыртқы күш жоқ болған кездегі қасиеттері тұрақты болатын сығылмайтын ньютондық тұтқыр сұйықтың жазық ағынын сипаттайтын негізгі теңдеулер: қозғалыс мөлшерінің екі теңдеуі (Навье - Стокс теңдеулері) және үзіксіздік теңдеуінен құйын тасымалдау теңдеуі мен ток функциясын анықтаған. Құйын тасымалдау теңдеу тәртібінің көптеген аспектісін зерттеу үшін, модельдік теңдеу ретінде конвективті және диффузиялық мүшелері бар сызықтандырылған бірөлшемді теңдеу мен Бюргерс теңдеуін қарастырылған.

Аталған теңдеулер бойынша әр түрлі сандық сұлбаларды келтіріп, зерттелген. «Ток функциясы, құйын» айнымалылары бойынша сығылмайтын сұйықтың торлық теңдеулеріне әр түрлі сұлбалар бойынша сандық есептеулер жүргізіп, зерттелген.

Көлемі 40 беттен тұратын бітіру жұмысы кіріспе бөлімінен, 3 бөлімнен, қорытынды бөлімнен, пайдалынылған әдебиеттер тізімінен тұрады. Жалпы бітіру жұмыстың құрамына 7 сурет кіреді. Бітіру жұмыстың соңында қосымшалар бар.

Кілт сөздер: сығылмайтын сұйық, конвективті және диффузиялық мүшелері бар сызықтандырылған бірөлшемді теңдеу, құйын тасымал теңдеуі, сұлбалар, Том формуласы, Вудс формуласы, орнықтылық шарт.

АНЫҚТАМАЛАР

Гидродинамика-
гидроаэромеханиканың сығылмайтын сұйықтықтың қозғалысын және оның өзімен шекаралас орналасқан қатты денемен әсерлесуін зерттейтін бөлімі.
Гидродинамика-: Рейнольдс саны-
гидроаэромеханиканың сығылмайтын сұйықтықтың қозғалысын және оның өзімен шекаралас орналасқан қатты денемен әсерлесуін зерттейтін бөлімі.: (ағылшын ғалымы О. Рейнольдстың атымен) - инерциялық күш пен тұтқырлық күш арасындағы қатынасты анықтайтын, тұтқыр сұйықтық пен газ ағысының ұқсастық критерилерінің бірі.
Гидродинамика-: Тұтқырлық -
гидроаэромеханиканың сығылмайтын сұйықтықтың қозғалысын және оның өзімен шекаралас орналасқан қатты денемен әсерлесуін зерттейтін бөлімі.: сұйықтар мен газдардың негізгі қасиеттерінің бірі, сұйыктың қозғалысына кедергі жасайтын ішкі үйкеліс күштері.
Гидродинамика-: Сызықтық теңдеу -
гидроаэромеханиканың сығылмайтын сұйықтықтың қозғалысын және оның өзімен шекаралас орналасқан қатты денемен әсерлесуін зерттейтін бөлімі.: белгісіздері (айнымалы шамалары) 1-дәрежелі болып келетін және белгісіздердің көбейтінділері қатыспайтын теңдеу
Гидродинамика-: Диффузия-
гидроаэромеханиканың сығылмайтын сұйықтықтың қозғалысын және оның өзімен шекаралас орналасқан қатты денемен әсерлесуін зерттейтін бөлімі.: (лат. dіffusіo - таралу, жайылу) - нақтылы дене бөлшектерінің жылулық қозгалыстарга ұшырай отырып, сол дене конңентрациясының селдір аудандарына қарай жылжуы; молекулалардың жылулық қозғалысы салдарынан шеқаралас орналасқан әр түрлі заттардың бір-біріне өту құбылысы. .
Гидродинамика-: Проекция (лат. projectio - алға лақтыру) -
гидроаэромеханиканың сығылмайтын сұйықтықтың қозғалысын және оның өзімен шекаралас орналасқан қатты денемен әсерлесуін зерттейтін бөлімі.: бұйымның бетке, көбінесе жазықтыққа белгілі бір әдіспен тұрғызылған кескіні.
Гидродинамика-: Айырымдық сұлба (the difference circuit) -
гидроаэромеханиканың сығылмайтын сұйықтықтың қозғалысын және оның өзімен шекаралас орналасқан қатты денемен әсерлесуін зерттейтін бөлімі.: дифференциалдық тендеулер мен белгілі бір нүктелердегі туынды функциялар мәндерінің ақырғы саны арқылы жуықталып ұсынылған айырымдық теңдеулер жүйесінің қосымша шарттарының аппроксимациясы.
Гидродинамика-: Итерация (лат. іteratіo - қайталау) -
гидроаэромеханиканың сығылмайтын сұйықтықтың қозғалысын және оның өзімен шекаралас орналасқан қатты денемен әсерлесуін зерттейтін бөлімі.: қандай да бір математикалық амалды қайталап қолдану.
Гидродинамика-: Түйін -
гидроаэромеханиканың сығылмайтын сұйықтықтың қозғалысын және оның өзімен шекаралас орналасқан қатты денемен әсерлесуін зерттейтін бөлімі.: тордағы нүктелер.
Гидродинамика-: Конвекция (лат. convectio - әкелу, жеткізу) -
гидроаэромеханиканың сығылмайтын сұйықтықтың қозғалысын және оның өзімен шекаралас орналасқан қатты денемен әсерлесуін зерттейтін бөлімі.:

1. атмосферада - жер бетіндегі неғұрлым жылыған (тығыздығы кем) ауа массасының немесе ағынының жекелеген бөліктерінің жоғары көтеріліп, онымен бір мезгілде неғұрлым салқын (тығыздау) ауа массасының төмен түсуі.

2. мұхиттағы конвекция- температура немесе тұздылықтың өзгеруі нәтижесінде судың тығыздығы өзгеруінен туындайтын вертикаль қозғалысы.

Белгілеулер

ρ \rho -
Тығыздығық
ρ−\rho -: u , υ u, \ \upsilon -
Тығыздығық: жылдамдық.
ρ−\rho -: p p -
Тығыздығық: қысым.
ρ−\rho -: ω \omega -
Тығыздығық: Құйын
ρ−\rho -: ψ \psi -
Тығыздығық: ток функциясы
ρ−\rho -: L L -
Тығыздығық: өзіндік ұзындық.
ρ−\rho -: C r Cr -
Тығыздығық: Курант саны.
ρ−\rho -: U 0 U_{0} -
Тығыздығық: есептің өзіндік жылдамдығы
ρ−\rho -: R e Re -
Тығыздығық: өлшемсіз параметр, Рейнольдс саны
ρ−\rho -: α \alpha
Тығыздығық: құйын тасымалдау теңдеуіндегі 1 / R e 1/Re өлшемге сәйкес келетін диффузияның жалпыланған коэффициенті
ρ−\rho -: I I -
Тығыздығық: жорамал бірлік
ρ−\rho -: G G -
Тығыздығық: өту көбейткіші

МАЗМҰНЫ

Кіріспе . . . 6
: 1
Кіріспе. . . 6: «ТОК ФУНКЦИЯСЫ, ҚҰЙЫН» АЙНЫМАЛЫЛАРЫНДАҒЫ ТҰТҚЫР СЫҒЫЛМАЙТЫН СҰЙЫҚТЫҢ ҚОЗҒАЛЫСЫН СИПАТТАЙТЫН ТЕҢДЕУЛЕР . . . 8
: 1. 1
Кіріспе. . . 6: Есептің қойылымы . . . 9
: 1. 2
Кіріспе. . . 6: Модельдік есептер . . . 13
: 1. 2. 1
Кіріспе. . . 6:

Тасымалдаудың бірінші модельді теңдеуі:

конвективті және диффузиялық мүшелері бар сызықтандырылған бірөлшемді теңдеу . . . 13

: 1. 2. 2
Кіріспе. . . 6: Тасымалдаудың екінші модельді теңдеуі: Бюргерс теңдеу . . . 14
: 1. 3
Кіріспе. . . 6: Гидродинамиканың толық есебінің барлық шешу процедурасының жалпы түрі . . . 16
: 1. 4
Кіріспе. . . 6: Ток функциясы мен құйын үшін шекаралық мәндер . . . 17
: 1. 5
Кіріспе. . . 6:

Ток функциясы мен құйын үшін шекаралық шарттардың

қойылымы . . . 20

: 1. 6
Кіріспе. . . 6: Құйынның шекаралық мәні үшін тестілік есеп . . . 22
: 2
Кіріспе. . . 6: САНДЫҚ СҰЛБАЛАРҒА ШОЛУ . . . 25
: 2. 1
Кіріспе. . . 6:

Конвекция-диффузия сызықтық теңдеуі үшін

айырымдық сұлбалар . . . 25

: 2. 2
Кіріспе. . . 6: Қуалау әдісі . . . 30
: 3
Кіріспе. . . 6: САНДЫҚ ЕСЕПТЕУЛЕРДІҢ НӘТИЖЕЛЕРІ . . . 34
:
Кіріспе. . . 6: Қорытынды . . . 38
:
Кіріспе. . . 6: Пайдаланылған әдебиеттер . . . 39
:
Кіріспе. . . 6:

Кіріспе

Жұмыстың жалпы сипаттамасы . Дипломдық жұмысым «ток функциясы, құйын» айнымалыларындағы тұтқыр сығылмайтын сұйықтың қозғалысын сипаттайтын теңдеулер үшін айырымдық сұлбаларды құру мен зерттеуге арналған.

Зерттеу нысанасы - «ток функциясы, құйын» айнымалыларындағы бірөлшемді тұтқыр сығылмайтын сұйықтың қозғалысын сипаттайтын теңдеу мен модельдік теңдеулер .

Мәселенің (проблема) өзектілігі. Техника және жаратылыстанудың дамуында сұйық қозғалысысының заңын зерттеу әрқашан маңызды роль ойнаған. Бұл саладағы зерттелу авиация, кеме жасау, жылу энергетикасы, атом энергетикасы, геофизика және т. б. қажеттіліктеріне жағдай жасап отырады. Осы соңғы он жылда сұйық қозғалысына байланысты зерттелу сферасы мен осы сферадағы құбылыс қолданысы айтарлықтай кеңейді. Бұл сфераға техниканың алдыңғы қатарлы бағыттарымен ( химиялық технология, металлургия, мұнай өндірісі және т. б. ) қоса негізгі жаратылыстану ғылымдары

( биология, атмосфера және мұхит физикасы және т. б. ) кіреді.

Сұйық динамикасын зерттеу барысында шыққан әр түрлі болып келетін есептер теориялық жолмен немесе мұқият та ұқыпты қойылған физикалық эксперимент көмегімен зерттелуі мүмкін. Көптеген жағдайларда сұйық ағысы кезінде өз орны болатын құбылысты модельдеу лабораториялық және табиғи жағдайында аса қиындатылған. Осы сұйық қозғалысынын зерттеу сияқты бағыттардағы физикалық эксперименттер әдетте техникалық жағынан ауыр, қиын және қымбатқа шығады. Онымен қоса, тәжірибелі өлшеулердің берілімі жалпы жағдайда шектеулі сипатқа ие болып келеді. Сол себептен гидродинамикалық зерттеулерде математикалық модельдеу маңызды роль атқарады.

Сығылмайтын сұйықтың дифференциалдық теңдеулер жүйесінің сандық шешілуіне көптеген монографиялар мен ғылыми мақалалар арналған, оның ішінде қазақ зерттеушілерінің құнды жұмыстары баршылық. Есептеу гидродинамика саласында жұмыс атқаратын мамандар үшін тұтқыр сығылмайтын сұйықтың Навье-Стокс теңдеуінің шешіміне тиімді сандық алгоритм құру аса қызығушылық тудырады. Бұл жүйені «ток функциясы, құйын» айнымалылары бойынша қарастыратын болсақ, есептеулік және теориялық қиындықтар туады.

Сондықтан сығылмайтын сұйық теңдеулерін шешудегі айырымдылық әдістер теориясының ары қарай дамуы есептеу математикасында өзекті мәселелердің бірі болып табылады.

Жұмысымның мақсаты «ток функциясы, құйын» айнымалыларындағы сығылмайтын сұйықтың торлық теңдеулері үшін ең тиімді итерациялық сұлбаларды зерттеу болып табылады.

Жұмыстағы көрсетілген мақсатқа жету үшін келесі міндеттер қойылды:

  1. сығылмайтын сұйықтың торлық теңдеулері үшін итерациялық

сұлбалар құру.

  1. айқын емес итерациялық сұлбалардың орнықтылығы мен дәлдікке зерттеу.
  2. сандық есептеулер жүргізу және алынған нәтижелерге талдау

жасау.

1 «ТОК ФУНКЦИЯСЫ, ҚҰЙЫН» АЙНЫМАЛЫЛАРЫНДАҒЫ ТҰТҚЫР СЫҒЫЛМАЙТЫН СҰЙЫҚТЫҢ ҚОЗҒАЛЫСЫН СИПАТТАЙТЫН ТЕҢДЕУЛЕР

Тұтқыр сығылмайтын сұйықтың ағынын сипаттайтын Навье - Стокстың теңдеуі көптеген жылдар ағынында жеке туындылы теңдеулерді шешу мәселесімен айналысатын зерттеушілер мен сандық талдау аймағында жұмыс атқаратын мамандардың назарын аударып келеді. Бұндай қызығушылыққа қарамастан, үш кеңістік айнымалылар жағдайындағы стационарлы емес Навье - Стокс теңдеуінің шешімінің болуы мен жалғыздығы туралы сұрақ әлі де ашық қалуда.

Бұл теңдеулердің сандық шешу жағдайы өте ауыр сипатқа ие. Мұндағы мәселе, бір түрдегі есепті жақсы шығаруда көрінген сандық әдіс басқа түрдегі есепті шығаруында тиімсіз болып келеді. Өлшемнің аз болып келуі тордағы қадамның ұсақталуына байланысты проблемаға тірейді.

Тұтқыр сұйық ағыны туралы есептің шешімі сызықты емес теңдеулердің шешілуін қажет етеді, және де оның сызықты еместігі теңдеудің сол жақ бөлігіндегі тұрған, жылдамдықтың конвективті бөлігін сипаттайтын инерциондық мүшесінде тұр. Бұл мүшені оған жақындатылған сызықты өрнекке ауыстыратын болсақ теңдеудің сызықтандырылуына алып келеді. Бұндай сызықтандырудың қарапайым мысалы ретінде Стокстың классикалық есебі бола алады.

Навье - Стокс теңдеулер жүйесін интегралдау кезінде нақты шешімді алу үшін шекаралық шарт қолданылуы тиіс, егер қозғалыс стационарлы емес болған кезде бастапқы шарттар қолданылуы керек. Тұтқыр сұйықтың қозғалысы туралы жеке түрдегі есептер үшін шешімнің бар болуы мен жалғыздығы туралы теоремалардың қатал дәлелдеулері бар [1] . Бұл теоремалар өзінің жалпы математикалық мазмұнын қоса дифференциалдық теңдеулерге қойылған шекаралық және бастапқы шарттардың қандай болу керек екендігін көрсететіні үшін маңызды.

Қазіргі уақытта гидродинамика есептеріне негізделген сандық әдістерде бағыттардың қатары анықталған. Оның інінде ерекше орын алатын әдістер: ақырлы айырымдық әдістер, «ірі бөлшектер», ақырлы элементтердің, интегралдық қатынастардың, торлық-варияциялықтың және т. б. әдістер болыа табылады. Тұтқыр сұйық динамикасына байланысты айтатын болсақ, онда бұл жерде ақырлы айырымдық әдісін қолдану аса сәттілік алып келуде. Бұл әдіс өзінің қарапайымдылығы мен негізінің жан жақты болып келуімен ерекшеленеді және де нәтижесінің жоғарғы деңгейде дәл шығуын қамтамасыз етеді. Бұл әдіс әр түрлі шекаралық және бастапқы шарттары бар сызықты емес теңдеулері мен сызықты да теңдеулерін сандық шешу үшін қолданылады.

Сызықты айырымдық сұлбалардың теориялары, яғни айырымдық сұлбалардың аппроксимация, рет, орнықтылық, жинақтылық сияқты қастиеттерін зерттеу А. А. Самарский, А. В Гулин [2], А. А. Самарский, Е. С. Николаев [3], А. А. Самарский [4], А. А Самарский., Р. Д. Лазаров, В. Л. Макаров [5], В. М Ковеня [6] және де тағы басқа өз еліміздегі және шетел елдеріндегі есептеуіш математиктердің монографиярында кеңінен, тереңінен жазылған.

Навье-Стокс теңдеуін шешуде Н. Т. Данаев, Ш. С Смагулов [15] үлкен еңбегін сіңдірген. Бұл еңбекте қатты өзгермелі коэффициенттері бар Навье-Стокс теңдеуі үшін тиімді сандық әдістер қарастырылған.

Екі айнымалы болған кездегі Навье-Стокс теңдеуін шешу тәсілінің бірі ток функциясы мен құйын айнымалыларын қолдануында бекітіледі және де сығылмайтын сұйық үшін екіөлшемді теңдеуді шешудің таралған әдістерінің бірі болып саналады.

Сандық әдістердің көптеген бөлігі ток функциясы мен құйын айнымалыларын қолданғандағы теңдеулер жүйесі үшін құракстырылған. Бұл әдістердің жалпы кемшілігі қатты беттегі құйын үшін шекаралық шарттарды қолдану болып табылады, ол есептің физикалық қойылымда жоқ.

Стационарлы емес күйде жағдай нашарлай түседі, ψ \psi -дан уақыт бойынша алынған туынды қосылады. Онымен қоса, ток функциясы үшін екі шекаралық мән беріледі, ал құйын шекаралық мән жоқ болып келеді. Бұның бәрі арнайы сандық әдістерді құруды талап етеді. Ереже бойынша бұл әдістер құйын үшін жасанды шекаралық шарт енгізуге негізделген. (Том, Вудс формуласы және т. б) [22], [23] .

Шаршы кавернадағы тұтқыр сығылмайтын сұйық ағынына сандық зерттеу [21] мақалада қарастырылған. Және де есептің қойылымы бір белгісіз функция-ток функциясы алынған, ол үшін төртінші ретті теңдеу жазылған.

Бұл мақаланың аса бір ерекшелігі кавернаның жоғарғы беті ақырындап нөлден кей шекті мәнге дейін өсетін айнымалы жыламдықтан тұратыны ([21] - те жоғарғы бетінің шекті мәні 1 тең) .

1. 1 ЕСЕПТІҢ ҚОЙЫЛЫМЫ

Сыртқы күш жоқ болған кездегі қасиеттері тұрақты болатын сығылмайтын ньютондық тұтқыр сұйықтың жазық ағынын сипаттайтын негізгі теңдеулер қозғалыс мөлшерінің екі теңдеуі (Навье - Стокс теңдеулері) және үзіксіздік теңдеуі болып табылады және келесі түрде жазылады:

u t + u u x + υ u y = 1 ρ P x + ν ( 2 u x 2 + 2 u y 2 ) , ( 1. 1 . 1 ) \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \frac{\partial u}{\partial t} + u\frac{\partial u}{\partial x} + \upsilon\frac{\partial u}{\partial y} = - \frac{1}{\rho}\frac{\partial P}{\partial x} + \nu\left( \frac{\partial^{2}u}{\partial x^{2}} + \frac{\partial^{2}u}{\partial y^{2}} \right), \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ (1. 1. 1)

υ t + u υ x + υ υ y = 1 ρ P y + ν ( 2 υ x 2 + 2 υ y 2 ) , ( 1. 1 . 2 ) \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \frac{\partial\upsilon}{\partial t} + u\frac{\partial\upsilon}{\partial x} + \upsilon\frac{\partial\upsilon}{\partial y} = - \frac{1}{\rho}\frac{\partial P}{\partial y} + \nu\left( \frac{\partial^{2}\upsilon}{\partial x^{2}} + \frac{\partial^{2}\upsilon}{\partial y^{2}} \right), \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ (1. 1. 2)

u x + υ y = 0 ( 1. 1 . 3 ) \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \frac{\partial u}{\partial x} + \frac{\partial\upsilon}{\partial y} = 0\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ (1. 1. 3)

Теңдеулер 𝐕 \mathbf{V} жылдамдығының u u , υ \ \upsilon векторлары мен P қысымнан тұратын физикалық айнымалылар үшін жазылған. Сұйықтың қасиеттері тығыздық

ρ = c o n s t \rho = const пен тұтқырлық кинематикалық коэффициентімен ν = c o n s t \nu = const сипатталады. Бұл теңдеулер келесі физикалық заңдарға негізделген: (1. 1. 1), (1. 1. 2) теңдеулер қозғалыс мөлшерінің векторлық теңдеуінің F = m a F = ma (Ньютонның екінші заңы) проекциясы болып табылады, онымен қоса тұтқыр күштер жанама кернеу үшін сызықты ньютондық заң бойынша деформация жылдамдығымен байланысқан, ал (1. 1. 3) теңдеуі массаның сақталу заңын білдіреді. Келтірілген теңдеулер координатаның эйлерлік жүйесінде жазылған, дәлірек айтқанда қозғалмайтын жүйеде. Сұйық сол қозғамайтын жүйеге салыстырмалы қозғалады. Тікелей осы теңдеулерді сандық шешуге болатынына қарамастан, өте жақсы нәтижелерді құйын мен ток функция үшін сандық шешуде алуға болады.

(1. 1. 1), (1. 1. 2) теңдеулерінің біріншісін у бойынша, ал екінші теңдеуді х бойынша дифференциалдап, шыққан нәтижені бір бірінен алып, қысымды шығырып тастауға болады. Құйынды төмендегідей анықтап

ω = u y υ x , ( 1. 1 . 4 ) \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \omega = \frac{\partial u}{\partial y} - \frac{\partial\upsilon}{\partial x}\, \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ (1. 1. 4)

параболалық типті құйынды тасымалдау теңдеуін аламыз:

ω t = u ω x υ ω y + ν ( 2 ω x 2 + 2 ω y 2 ) = 𝐕 ( ω ) + ν 2 ω , ( 1. 1 . 5 ) \ \ \ \ \ \ \ \ \ \ \ \ \frac{\partial\omega}{\partial t} = - u\frac{\partial\omega}{\partial x} - \upsilon\frac{\partial\omega}{\partial y} + \nu\left( \frac{\partial^{2}\omega}{\partial x^{2}} + \frac{\partial^{2}\omega}{\partial y^{2}} \right) = - \mathbf{V} \bullet (\nabla\omega) + \nu\nabla^{2}\omega, \ \ \ \ \ \ \ \ \ \ \ \ \ (1. 1. 5)

Субстанциональды туынды қолданып, бұл теңдеуді келесі түрде көрсетуге болады:

D ω d t = ν 2 ω ( 1. 1 . 6 ) \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \frac{D\omega}{dt} = \nu\nabla^{2}\omega\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ (1. 1. 6)

Ток функциясын ψ \psi келесі арақатынастармен анықтап

ψ y = u , ψ x = υ , ( 1. 1 . 7 ) \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \frac{\partial\psi}{\partial y} = u, \ \ \ \ \ \ \ \frac{\partial\psi}{\partial x} = - \upsilon, \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ (1. 1. 7)

(1. 1. 4) теңдеуді эллиптикалық типті Пуассон теңдеуі ретінде жазуға болады:

2 ψ = ω . ( 1. 1 . 8 ) \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \nabla^{2}\psi = \omega. \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ (1. 1. 8)

Құйынды тасымалдау теңдеу (1. 1. 5) құрамына тұрақты емес мүше ω \partial\omega / t \partial t , конвективті мүшелер u ω u\partial\omega / x \partial x және υ ω \upsilon\partial\omega / y \partial y және тұтқыр диффузиямен байланысты мүше ν 2 ω \nu\nabla^{2}\omega кіреді. Бұл теңдеу конвективті мүщелері үшін сызықты емес, өйткені (1. 1. 7) және (1. 1. 8) күшінен u u мен υ \upsilon өздігімен ω \omega тәуелді айнымалының функциялары түрінде көрсетіледі. Ол уақыт бойынша параболалық болып табылады, сондықтан оған бастапқы шарттары болатын есеп қойылады. Ол есепте шешім кейбір бастапқы берілімнен әр қадам бойынша «жылжиды».

Ток функциясының (1. 1. 8) теңдеуі эллиптикалық болып табылады, сондықтан оған итерациялық әдіспен шығарылатын, шекаралық шарттары берілген есеп қойылады. Көптеген практикалық есептерде шешімнің уақыт бойынша шығу нәтижесі емес, тек тұрақты шешімі ғана қызықтырады; бұл жағдайда (1. 1. 5) теңдеуінің сол жақ бөлігіне ω \partial\omega / t \partial t =0 қойып, бір тәуелсіз айнымалы - уақытты алып тастаймыз. Негізі ереже бойынша аналитикалық зерттеу жүргізген уақытта дәл солай істейміз; сондықтан осы есептеу гидродинамикасымен жұмыс істеп көрмеген адамдар гидродинамиканың тіпті стационарлы есептерінің көптеген(тек барлығы емес) тиімді сандық әдістерімен шешілуі стационарлы емес теңдеулерді интегралдауға негізделетініне, ал стационарлы шешім (егер ол болса) стационарлы емес теңдеулерінің шешілу шегі уақыт бойынша асимптоталы болып келетініне таңқалады.

Тағы бір есекеретін жағдай, құйынды тасымалдау теңдеуі (1. 1. 5) басқа да көптеген үрдістерге модельдік сипаттама жасау қызметін атқарады.

Әдетте математиктер жеке туындылы (сызықты) дифференциалдық теңдеулердің келесі типтегі классификациясына қанағаттанады: параболалық, эллиптикалық және гиперболалық. Бұндай классификацияда құйын тасымадау теңдеуі мен диффузия теңдеуінің ω \partial\omega / t \partial t = α 2 ω \alpha\partial^{2}\omega / x 2 \ \partial x^{2}\ арасында айырмашылық жасалмайды, бірақ, төменде көрсетіп отырғанымдай, (1. 1. 5) теңдеуде бірінші ретті туындының болуы, оның диффузия теңдеуіне қарағанда сапалы өте жақсы екендігін көрсетеді, және де конвективті мүше сандық шешу кезінде аса маңызды рөл ойнайды. Өкінішке орай, көрсетілген екі мүше үшін әртүрлі сандық сұлбалар ең тиімдірек болып шығуы мүмкін.

Теңдеудің консервативті формасы

Үзіксіздік теңдеуін (1. 1. 3)

u x + υ y = 0 \frac{\partial u}{\partial x} + \frac{\partial\upsilon}{\partial y} = 0

толық жылдамдық векторы арқылы келесі түрде жазуымызға болады:

𝐕 = 0 \nabla\mathbf{\bullet V} = 0 (1. 1. 9)

( 𝐕 ω ) \nabla \bullet \left( \mathbf{V}\omega \right) қарастырайық. Векторлық алгебрада келесі тепе - теңдік белгілі

( 𝐕 ω ) = 𝐕 ( ω ) + ω ( 𝐕 ) = 𝐕 ( ω ) . \nabla \bullet \left( \mathbf{V}\omega \right) = \mathbf{V} \bullet (\nabla\omega) + \omega\left( \nabla \bullet \mathbf{V} \right) = \mathbf{V} \bullet (\nabla\omega) .

Сонымен құйын тасымалдау теңдеуінің консервативті формасын алу үшін (1. 1. 5) теңдеуде 𝐕 ( ω ) \mathbf{V} \bullet (\nabla\omega) мүшесін ( 𝐕 ω ) \nabla \bullet \left( \mathbf{V}\omega \right) осы мүшеге ауыстыру керек, нәтижесінде

ω t = 𝐕 ( ω ) + ν 2 ω = ( u ω ) x ( υ ω ) y + ν ( 2 ω x 2 + 2 ω y 2 ) . ( 1. 1 . 10 ) \ \ \ \ \ \frac{\partial\omega}{\partial t} = - \mathbf{V} \bullet (\nabla\omega) + \nu\nabla^{2}\omega = - \frac{\partial(u\omega) }{\partial x} - \frac{\partial(\upsilon\omega) }{\partial y} + \nu\left( \frac{\partial^{2}\omega}{\partial x^{2}} + \frac{\partial^{2}\omega}{\partial y^{2}} \right) . \ \ \ \ \ \ \ \ \ \ (1. 1. 10)

осындай түрге келеді.

Өлшемсіз айнымалылар теңдеулері

Менің дипломдық жұмысымда қолданылған өлшемсіз айнымалылармен келетін теңдеулер жүйесі барлық жерінде конвективті масштабқа L / U 0 L/U_{0} негізделеді, мұндағы L \ L -өзіндік ұзындық, ал U 0 U_{0} -есептің өзіндік жылдамдығы; мысалы, егер L L -қанатты пішіндегі хорда ұзындығы және U 0 U_{0} -жүгірмелі ағынның жылдамдығы, онда L / U 0 L/U_{0} -уақыт, осы уақыт аралығында жүгірмелі ағынның бөлшегі бүкіл пішіннен(профиль) өтеді. Келесі өлшемсіз шамалар енгіземіз:

u ¯ = u U 0 , υ ¯ = υ U 0 , x ¯ = x L , y ¯ = y L , ω ¯ = ω U 0 / L , t ¯ = t L / U 0 ( 1. 1 . 11 ) \overline{u} = \frac{u}{U_{0}}, \ \ \overline{\upsilon} = \frac{\upsilon}{U_{0}}, \ \ \overline{x} = \frac{x}{L}, \ \ \overline{y} = \ \frac{y}{L}\, \ \ \ \ \ \overline{\omega} = \frac{\omega}{U_{0}/L}, \ \ \ \ \overline{\ \ t} = \frac{t}{L/U_{0}}\ \ \ \ (1. 1. 11) \

Осыдан (1. 1. 10) және (1. 1. 8) теңдеулер келесі түрге келеді

ω t = ( 𝐕 ω ) + 1 R e 2 ω , ( 1. 1 . 12 ) \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \frac{\ \partial\omega}{\partial t} = - \nabla \bullet \left( \mathbf{V}\omega \right) + \ \frac{1}{Re}\nabla^{2}\omega, \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ (1. 1. 12)

2 ψ = ω ( 1. 1 . 13 ) {\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \nabla}^{2}\psi = \omega\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ (1. 1. 13)

мұндағы R e \ Re - өлшемсіз параметр, Рейнольдс саны,

R e = U 0 L ν . ( 1. 1 . 14 ) \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ Re = \frac{U_{0}L}{\nu}. \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ (1. 1. 14)

... жалғасы

Сіз бұл жұмысты біздің қосымшамыз арқылы толығымен тегін көре аласыз.
Ұқсас жұмыстар
Еркін ағыстың сипатына сыртқы акустикалық әсердің ықпалы
МЕХАНИЗМДЕРДІҢ ЕСЕПТЕУ КӨРСЕТКІШТЕРІ
Жазық турбуленттік ағыншаның бастапқы бөлігін есептеу
Ірі құйындар әдісімен пішіндеу
«Медеу бөгетінің суағытқыштары» ТУ абж үшін ОРС-серверді Masterscada құралдарымен жобалау және баптау
Абсолюттік нөлге жақын жағдайдағы материалдардағы ғажайыптар
Лабораториялық кешен сипаттамасы базалық сынақтарды жүргізу басқармасы
Магнетиктердің түрлері және қасиеттері
АТМОСФЕРАНЫҢ ШЕКАРАЛЫҚ ҚАБАТЫНЫҢ САНДЫҚ ҮЛГІЛЕРІ
Метеорологилық элементтерге анықтама, метеоөлшеуіш құралдар
Пәндер



Реферат Курстық жұмыс Диплом Материал Диссертация Практика Презентация Сабақ жоспары Мақал-мәтелдер 1‑10 бет 11‑20 бет 21‑30 бет 31‑60 бет 61+ бет Негізгі Бет саны Қосымша Іздеу Ештеңе табылмады :( Соңғы қаралған жұмыстар Қаралған жұмыстар табылмады Тапсырыс Антиплагиат Қаралған жұмыстар kz