TOPREFERAT.COM.KZ - Қазақша рефераттар

МАЗМҰНЫ - www.topreferat.com.kz

1 Дискретті теңдеулер

1.1 Ақырлы – айырымдық әдістің негізі

1.2 Аппроксимацияның негізгі әдістері. Тейлордың кесілген қатарымен берілген аппроксимация әдісі

1.3 Айырымдық схема құрудың вариациялық принципі

1.4. Шектік – айырымдық теңдеулер

2 Дербес туындылы теңдеулер үшін айырымдық схемалар

2.1. Есептің қойылу сипаты

2.2. Коэффициенттері тұрақты жылу өткізгіштік теңдеу

2.3 Сызықтық тасымалдау теңдеуін шешудің сандық әдістері

2.4. Консервативтік қасиет

2.5. Ағынға қарсы өзгертілген схеманың бір түрі

2.6. Бір өлшемді сызықтық емес теңдеулерді шешу

Қорытынды

Пайдаланған әдебиеттер тізімі

Жұмыс түрі: Дипломдық жұмыс
Жұмыс көлемі: 70 бет
Пәні: Соңғы қосылған дипломдық жұмыстар

-----------------------------------------------------------------------------------

ДИПЛОМДЫҚ ЖҰМЫСТЫҢ ҚЫСҚАРТЫЛҒАН МӘТІНІ

Мазмұны

1 Дискретті теңдеулер

1.1 Ақырлы – айырымдық әдістің негізі

1.2 Аппроксимацияның негізгі әдістері. Тейлордың кесілген қатарымен берілген аппроксимация
1.3 Айырымдық схема құрудың вариациялық принципі

1.4. Шектік – айырымдық теңдеулер

2 Дербес туындылы теңдеулер үшін айырымдық схемалар

2.1. Есептің қойылу сипаты

2.2. Коэффициенттері тұрақты жылу өткізгіштік теңдеу

2.3 Сызықтық тасымалдау теңдеуін шешудің сандық әдістері

2.4. Консервативтік қасиет

2.5. Ағынға қарсы өзгертілген схеманың бір түрі

2.6. Бір өлшемді сызықтық емес теңдеулерді шешу

Қорытынды

Пайдаланған әдебиеттер тізімі

Кіріспе

Соңғы кезде есептеу математикасына байланысты көп есептер, сызықты емес
Диссертациялық жұмыстың негізгі мақсаты – Сызықтық емес теңдеулерді шешу,
Біз уақыт бойынша айнымалыға және оған сәйкес дифференциялдық операторға
Үзіліссіз дәл шешімдер үшін дененің кез келген облысына сақталу
Анықтама. Энергияның сақталу заңын қанағаттандыратын (интегралдық инварианттар)сандық шешімдер алынатын
Консервативтік қасиеттің маңыздылығын сығымдалушы ортаның үзіліс еместік теңдеуінің мысалында
1 Дискретті теңдеулер

1.1 Ақырлы – айырымдық әдістің негізі

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

кеңістігінің үздіксіз аргументінің
Көп жағдайларда функциясын Фурье қатарына
(1.2)

ң мәндері коэффициенттерінің
(1.3)

Осылайша таңдап, дискретті нүктелерінің
Көп жағдайда төмендегідеу центрлік айырымдар қолданылады.

Центрлік (орталық) айырымдар, осы айырым есептелетіндей нүктеге қатысты симметриялы
Дифференциалдық теңдеулерді аппроксимациялауларды құру жолдарының бірі –туындылардың ақырлы –
(1.4)

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

Алға (вперед) айырымының көмегімен өрнекті төмендегідей түрде жазамыз.

аппроксимация қателігі.

Мұндағы айырымы
-торлық мәнінің орнына туындылар үшін жуық өрнекке
(1.1) қатынасы үшін мұндай процедура төмендегідей өрнекке әкеледі.

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

Тейлордың кесілген қатарымен берілген аппроксимация әдісі

Айталық, кеңістік бойынша
Тейлор қатарын функциясы үшін (1.5) формуласы
,
Егер оның екінші мүшесінен бөліп алсақ, онда (1.6) формула
.
Демек (1.7) формуласы орынды,

.
Интегралдық әдіс. Сызықты Хопфа теңдеуін дивергентті (ажырату) түрде қарастырайық.

Кеңістікте –ден –ге және
.
Жақша ішіндегі өрнектерді интегралдау арқылы төмендегіні (1.9) аламыз:

.
Қалған интегралдарды – кеңістікті интегралды орта туралы теорема, ал
(1.10)

(1.10) қатынасынан орта туралы теораманы қолданудан қалған туындылар
.
(1.11) – ді (1.10) – ға қою және
.
Көпмүшелілік жуықтау. Айырымдық өрнекті алу әдісі – алдымен тордың
фукциясының мәні және
,
ыңғайлы болу үшін
,
береді.

Онда (1.14) формуласы арқылы келесі (1.15) қатынасты аламыз:

,
- нүктесінде бірінші және екінші ретті туындылар сәйкесінше (1.16)
,
түрінде болады.

(1.16) есеп – қитаппен (1.15) формулалары дәлдіктерімен 2 –
1.3 Айырымдық схема құрудың вариациялық принципі

Айырымдық схема құру әдісінің ең көп тараған әдісі баланстың
Торлық облыстың элементар ұяшығы («бокс») үшін баланс теңдеуінен айырымдық
Жылу өткізгіштік теңдеуінде (1.17) вариациондық әдісті қолдануды
,
,
мұндағы бойынша периодтылық
.
Бұдан бөліктеп интегралдаудан (1.19) интегралдық тепе –теңдікті аламыз:

, (1.19)

мұндағы – кіріс параметрлерінің векторы.Торлық облысты бір
,

,

және де және
(1.20)

Мұндағы . Егер кез келген
,
түріндегі жазуға болатын (1.19) тепе – теңдігі бастапқы және
функциялары салыстырмалы түрде (1.17) және (1.19) да барлық амалдардың
.
Айырымдық аппроксимация құру процедурасы келесіден тұрады: Бастапқыда интегралдар квадратуралық
,
(1.19) аппроксимациясының қателігі әлсіз мағынада екі
(1.24)

Мұндағы

.

1.4. Шектік – айырымдық теңдеулер

Шектік – айырымдық жақындау сәйкестенділетін дифференциалдық теңдеулерде туындыны ауыстыру
,
мұндағы – оң тұрақты. (1.25) теңдеуі
аламыз.Бұл өрнектерді (1.25) – ке апарып қойсақ
.
Бұл түзулер адвекция теңдеуінің сипаттамалары болып табылады.Солардың бірі (сурет-1.1)
.

Сурет 1.1 - Сызықты адвекция теңдеуінің характеристикаларының бірі

Бұл схема «ағынға қарсы» деп аталуы мүмкін.Соңғы сөз
Сурет 1.2 - Адвекция теңдеуі жуық шешімін анықтау үшін
айырымдық сұлба

өсімшелері нольге ұмталғанда осы дифференциалдық теңдеуге жуықталатын көптеген (1.27)
,
(1.27) тасымалдау теңдеуіне (1.28) формуласы орындалады

,
диффузия теңдеуіне гидрогазодинамика теңдеуі үшін модельді болатын айырымдық схема
арқылы торабындағы
,
,
.
болғандықтан (1.29) және (1.31) сұлбалары (1.27) – (1.28)
,
(1.27) теңдеуі үшін және
Анықтама: Келесі қабатта функцияның бір ғана мәні болатын әрбір
Жоғарыда қарастырылған схемалар – айқын схема (1.27) теңдеуіне айқындалмаған
,
алуға болады. Жартылай айқындалмаған екі қабатты схемалар класын (1.34)
.
Мұндағы – параметр.

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

болғанда болсын деп талап қояйық.

Келісімділік. Дискретизация процесі нәтижесінде алынған алгебралық теңдеулер жүйесі тор
Орнықтылық. Айырымдық схема орнықты деп аталады, егер уақыттық координатаның
,
сандық шешімнің қателігі болып табылады.

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

1.Белгіленген қосындылы уақытта
2.Егер белгіленген мәндерінде
Бірінші сұраққа жауап сандық шешімнің жинақтылығына байланысты.Егер ұсақталған тордың
Оның (1.37) аппроксимация қателігі тродың ұсақталу кезінде нольге ұмтылады
Сурет 1.3 - Тәуелділік облыс және характеристика орналасуына қатысты

мүмкіндіктер

Координата басын қиятын характеристика А нүктесінен де өтеді.Сондықтан А
Яғни, болғанда.

Осылайша, бұл (1.30) сұлбасының жинақталуы үшін қажетті шарттары болып
Кейбір айырымдық теңдеулер үшін шешім анықтаудың қандайда бір кезеңінде
Эквиваленттік туралы Лакс теоремасы. Егер шекаралық шарттармен берілген нақты
Максимум принципі. Дәл шешімнің шектелген екендігі белгілі, онда сандық
(1.38)

Жинақтылықтың қажеттілік шартымен сәйкес келетін шартында
,
уақыттық деңгейде максимумын қабылдайтын нүктеде бұл
Бұл теңсіздік сандық шешімнің шектелгендігін дәлелдейді. Яғни,
Энергетикалық әдіс. Бұл әдіс өте көп қолданылады.Оны сызықты емес
.
Циклдік шекаралық шарттар орындалатын деп есептейміз,мысалға
Сонда (1.41) өрнегі келіп шығады

.
Енді Шварц теңсіздігін қолданайық және (1.41)-
.
(1.41) және (1.42) қолдана отырып және егер
Сонымен, (1.38) схемасының орнықты болуы үшін
Бастапқы берілгендері бойынша орнықтылық. Тасымалдау теңдеуі үшін (1.43) Коши
,
есебі үшін (1.44) айырымдық теңдеуін қарастырамыз.

, (1.44)

(1.44) айырымдық теңдеуін бастапқы берілгендері бойынша орнықтылығын (1.45) өрнегі
.
Мұндағы -тен тәуелсіз және
2 Дербес туындылы теңдеулер үшін айырымдық схемалар

2.1 Есептің қойылу сипаты

Кез келген есепті шешкенде оған кіру мәліметтерін алғашқы, ізделінуші
Әрбір есеп үшін белгілі,бірдей сұрақтар қойылады: есептің шешуі бола
Адамар шарты. Математикалық есептің жалпы жазылуын операторлық теңдеу түрінде
(2.1)

мұндағы u және f ізделінуші және әлдеқандай U
(2.1) теңдеуінің есебін шешу Адамар шарты бойынша дұрыс қойылған
1) Кез келген үшін
2) U - да бір ғана шешімі болады .

3) шешімі үзіліссіз түрде - қа тәуелді
Егер осы айтылған талаптардың біреуі болмаса (2.1) есебі
Дұрыс қойылған есептің мысалы, интегралдау есебі, ал сыпайы емес
Мысал. Интегралдау есебі f(x) функциясы берілген; (2.2) интегралын тап:

(2.2)

Функция – ты
, , .

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

Сонымен интегралды квадратуралық формулалар арқылы есептеу дұрыс қойылған есеп.

Мысал. Дифференциалдау есебі. Жуық шамамен берілген u(x) функциясын дифференциалдау
болсын, мұндағы барынша үлкен. Онда
Сонымен, С кеңістігінде u(x) функциясының аз өзгерісіне
2.2 Коэффициенттері тұрақты жылу өткізгіштік теңдеу

1. Бастапқы есеп. Бір өлшемді жағдайда диффузияның
(2.3.)

Мұндағы х нүктесіндегі t уақыт
Егер тұрақты болса, онда (2.3.) теңдеуін
(2.4)

мұндағы - жылуөткізгіштік коэффицинті. Жалпы шектемей-ақ
Біз бірінші шеттік есепті мына облыста қарастырайық (кейде бастапқы
(2.5)

шартымен және шекаралық шарттың диффузия теңдеуінің
,
.

Максимум принципінің күшімен (2.5) есебінің шешімі үшін төменлегі (2.7)
(2.7)

Біртекті шектік шартымен (2.8) түріндегі біртекті жылу өткізгіштік теңдеуді
(2.8)

Бұл есептің шешуін айнымалыларды айыру әдісімен (2.9) түрде табамыз:

(2.9)

мұндағы және -
,

Сонымен қатар скалярлық көбейтіндісі (2.10) теңдеулер жүйесі түрінде беріледі

,
Тегінде барлық дербес шешімдер (гармоникалар)
(2.11)

коэффициенттерін табады. (2.8) және (2.10) өрнектері (2.12) формуласын береді

(2.12)

өйткені

.

Сонымен (2.8) есебінің шешуіне (2.13) орынды:

,
бұл (2.7) есебінің алғашқы мәліметтері бойынша (
(2.14)

Процестің бұл стадиясы жүйелі режим деп аталады.

Енді жылу өткізгіштік теңдеуін шешетін әлдебір өте маңызды айырымдық
(2.15)

Уақыт бойынша бірінші ретті және х бойынша екінші
мұндағы .
Қарастырып отырған жағдайда модификацияланған теңдеудің түрі (2.17) түрінде болады:

(2.17)

болғанда аппроксимацияның қатесі тең болатындығын атап
2.1 сурет - Ауысу коэффициентіндегі айқын схемаға айқын әдіс
Аппроксимация қатесінің өрнегінде тақ ретті туындылардың кірмейтіндігін ескеру де
.
Оның үстіне ауысу коэффициентінің жорамал бөлігі нольге тең және
, немесе

(2.19)

мұндағы .

Демек, жылу өткізгіштік теңдеуінің дәл шешімінің амплитудасы уақыт бойынша
2.2 сурет - Айқын схемадағы тәуелсіз болмау нүктелерінің зонасы

2.2-суреттен, жылу өткізгіштік теңдеуін шешудің қарапайым айқын әдісі r=
Қарапайым айқын әдісті пайдаланғанда жылуөткізгіштік теңдеуі алғашқы мәліметтер
2.3 Сызықтық тасымалдау теңдеуін шешудің сандық әдістері

(2.20) сызықтық тасымалдау теңдеуін қарастырайық

,
с>0, с=const болғанда сипаттамалық теңдеудің түрі мынадай болады:
(2.20) теңдеуінің шекті – айырымдық шешімдерінің негізгі қателерді мыналар:

- дисперсиялық толқындар тудыратын фазалық: жекелеген гармоникалар физикалық жылдамдықтардан
Бұл тегістеуші әсер немесе схемалық есептеу тұтқырлығын тудырады.
- осцилляциялық немесе «есептеу шуы», кеңістіктік туындыны аппроксимаялуадың сандық
- амплитудалық, гармониканың шамасының өсуіне қарай орнықсыздығына немесе азаюына
Уақыт бойынша айқын емес аппроксимациялауда (субстанцияның тасымалдау жылдамдығы с-ны
(2.21)

мұндағы , егер с≥0 және кері жағдайда
Нейман орнықтылығының шартынан схема мына жағдайда орнықты болады:

, мұндағы – Курант саны

(2.21) - ге болғанда және
Жеңіл-желпі түрлендіруден кейін мынау шығады:

Соңғы теңдікті t және х бойынша дифференциялдап және –
Бұларды қоссақ

Осы жолмен табамыз:

Жоғарыдағы келтірілген барлық теңдіктерден (2.22) теңдеуі шығады:

(2.22)

Кейбір әдебиеттерде (2.22) теңдеуін айырымдық схеманың дифференциалдық жуықтауы деп
(2.22) диференциалдық жуықтауының оң жағы берілген дербес туындылы теңдеудің
.

Берілген дербес туындылы теңдеудің дәл шешімін алуға ықпал ететін
Біз қарастырып отырған жағдайдағы аппроксимация қатесінің өрнегіндегі бас мүше
Айырымдық схеманың басқа физикалық қасиетке жақынын дисперсия деп атайды.
Үлкен градиентті облыстарда қоспаларды тасымалдау теңдеулерін шешкенде «ара»
Сонымен тасымалдау теңдеуін шешкенде, алдыңғы параграфтардағы қарастырылған (аппроксимация,орнықтылық, консервативтік)
Анықтама. Біз схеманы монотонды деп атаймыз,егер ол барлық монотонды
Бірінші ретті схемалар үшін диссипативтік қатені азайтудың басқа бір
осы приципке негізделген.

(2.20) теңдеуінің шешімінің барлық схемаларын шешілген түрде (2.23) теңдігі
.
мұндағы j бойынша қосынды егер схема айқын емес болса,
Анықтама. (2.23) схемасы (2.20) теңдеуін шешу үшін Р индексіне
Монотондық схемалар жалған максимумдар мен минимумдар (жалған осциляциялар) бере
Теорема. (2.20) - ны аппроксимациялайтын және индексі 2 ≥р
Дәлелдеу. Алғы мәндер ретінде екінші дәрежелі көпмүшені алайық:

,

мұндағы һ – тордың адымы.

Тордың нүктелерінде

.

Қарапайымдылық үшін былайша ұйғарайық:

.

(2.20) - ның осының алғы мәндер бойынша дәл шешімі
.

Дегенмен, біз зерттелуші схеманың индексі р ≥2 деп ұйғардық,
.

Егер схема монотонды болса, онда барлық Sj-i
.

Бірақта бұл бұлай емес, мысалға i=1 үшін біз
,

бұл rБұл вариантта «ағынға қарсы» схемасы формалдық түрде екінші ретті
2.3 сурет - Адвенция теңдеуінің сандық шешімі

1 – дәл шешімі; 2 – ағынға қарсы айырымдардың
теңдеуін аппроксимациялдаушы әртүрлі схемалар үшін схемалық жасанды тұтқырлықтың коэффициенттері

Кесте 1 – Схемалық тұтқырлық формуласы коэффициенттері

Схема Формуланың жазылуы Схемалық тұтқырлық формуласы Аппроксимация Монотондық
1 2 3 4 5

Ағынға қарсы болғанда

Айқын
Айқын емес
Кранка-Никольсондыкі
Орталық айырымдардың айқын схемасы
Кесте 1 – дің жалғасы

Орталық айырымдардың айқын емес схемасы
Кранк-Никольсондікі
Лакстыкі
Такачтыкі

әлсізмонотонды емес

Лейтаныкі, Лакс-Вендрофатыкі, Мореттидыкі, Маккормаканыкі
Дәлдігі бірінші реттен жоғары монотонды емес схема үшін саны
Шекті-айырымдық схеманы жақсартатын бағыттардың бірі фильтр ретінде арнайы диффузиялық
Егер болса, (2.24) теңдіктері орындалады

(2.24)

Егер болса, (2.25) теңдіктері орындалады

(2.25)

Мұндағы - коэффициенттер,
Дәлдігі екінші ретті монотонды схема алу үшін
Осы кемшіліктерді есептеу функциясының градиенттері үлкен болатын және осциляциондық
Мұндағы ;

коэффицинтін бұл түрде таңдау, осцилляция болатын облыстарда ғана локалдық
Қолданбалы есептерде, Борис және Бук ұсынған FCT ағындардың клоррекция
1 Конвективтік және диффузиондық мәнін есептеу
(2.26)

мұндағы - оптимизацияға қажетті тасымалдаудың және
2 (2.27) теңдігі дөрекі антидиффузиондық ағындарын
(2.27)

мысалға сызықтық түрде анықталатын антидиффузиялық өлшемсіз коэффициент:

Шарт бойынша Ф түзетілген ағымдарын есептеу (2.28) теңдігінде көрсетілген

,
мұндағы , .

3 Жаңа мәнін есептеу бірақ
(2.29)

Такач схемасын қарастырайық:

(2.30)

α=0 болған жағдада кәдімгі екінші ретті Лакс-Вендроф схемасын
2.4 сурет - Қоспаны тасымалдау теңдеуінің сандық шешімі (1
2.5 сурет - Лакс-Вендрофа схемасын құру үшін қолданылатын кеңістікті-уақыттық
Қорытындысында Лакс және Вендроф ұынған схеманы талдайық, дәлірек айтқанда
(2.31)

Осы аралық мәндердің көмегімен және кеңістік бойынша орталық –
(2.32)

(2.31) - дағы аралық мәндері ( 2.32) - ге
(2.33)

( 2.33) – де соңғы
Оның үстіне бұл мүше енді 2
Бұл проблемаларды диссипативтік схемаларды пайдалану арқылы жұмсартуға болатындығын, әсіресе
((2.32) схемасының аппроксимация қатесі O[(∆x)2]+O[(∆t)2] болады. Демек, бұл схема
,

қоямыз. Сонда

(2.34)

Демек,

(2.35)

Дегенмен

, .

Түпкілікті алатынмыз:

.
Квадрат жақшаның ішіндегі өрнек екі квадраттардың қосындысы және ол
Бұл теңсіздік Курант орнықтылығының критеримен тағы да дәл келіп
(2.37)

ұзындығы екі есе үлкен толқынға (4∆x) R∆x= /2
.
Жалпы жағдайда мына нәтижеге байланысты

барлық қисықтары
Демек, толқынның ұзындығы 2 -
Ұзындығы 2 және 4
бұдан берілген периодта (белгіленген n∆t шамасында) жалпы демпфировандығы жуық
Бірақ, тегінде - ны, жақсы дәлдікті
2.6 сурет - Ланс-Вендроф схемасы үшін үлкейту көбейткіштерінің тәуелділігі.
Ланс-Вендроф схемасы гидродинамикалық үлгілерде кеңінен қолданылады. Бұл туралы Рихтмайердің
Лакс-Вендроф схемасының диссипативтігі қолайсыздықтар тудырмайды, егер бұл
Өзімізге қажетті жағдайда әсер
2.4 Консервативтік қасиет

Біз уақыт бойынша айнымалыға және оған сәйкес дифференциялдық операторға
Үзіліссіз дәл шешімдер үшін дененің кез келген облысына сақталу
Анықтама. Энергияның сақталу заңын қанағаттандыратын (интегралдық инварианттар)сандық шешімдер алынатын
Тасымалдау теңдеуін қарастырайық:

мұндағы - набла операторы.
Осы теңдеуді әлдеқандай кеңістік R облысы бойынша интегралдайық:
(2.40)

Остроградский – Гаус формуласын қолданып мынаны аламыз

(2.41)

мұндағы ∂R – R облысының шекарасы, u –
.
Сонда (2.39) теңдеуді мынадай түрде болады

(2.43)

(2.43) теңдеуді R облысындағы u шамасының жылдамдықтар жиынтығының
Консервативтікті талап ету шекті – айырымдық схемада интегралдық қатынастарды
(2.39) теңдеуінен алынатын тұтқыр емес сұйықтың шекті жағдайы үшін
.
Егер басқа жағынан u шамасын массаның тығыздығы ретінде түсіндіретін
,
(мұнда қарапайымдылық үшін жоғары n индексі жазылмалды). Енді бірөлшемді
,

бұл (2.43) теңдуіндегі интегралына сәйкес келеді:

(2.46)

Бірінші бөлімді қосу былайша орындалады:

+( )
+( )
+( )
+ ( )
+ . . . . . . .. .
+( )
+( )
+( )
+( )
= ( ) +( )
Сонда ( 5.45) теңдеуі мына түрде болады:

(2.47)

Берілген теңдеу R облысындағы
.

Алдыңғы мысалдағы схеманы пайдаланып, яғни уақыт бойынша алға айырымы
.
Сонда (2.48) - ге сәйкес келетін қосындының түрі
.

Қайтадан қосатын болсақ

=

=
+
+
+
+ . . . . . . . .
Бұдан, дифференциялдық теңдеудің мұндай формасындағы көршілес ячейкалардың қыры арқылы
(2.49)

болатын дербес жағдайдан басқасы. Яғни бұл жағдайда шекті –
Бірінші жағдайда консервативтік, шекті –айырымдық өрнекті қорыту кезінде
Консервативтік форманы пайдаланғанда шамасының центрі i
Консервативтік қасиеттің маңыздылығын сығымдалушы ортаның үзіліс еместік теңдеуінің мысалында
Консервативтік еместік әсердің мысалы.Айырымдық схема үзіліссіз ортаның
Осы текті мысалдардың біреуін баян етейік. Қарапайым үлгі –
, 0 бұл мына шартты қанағаттандырсын

(2.51)

Жылуөткізгіштік коэффициентінің k(x) нүктесі үзіліс нүктесі.
(2.52)

Бұл шарттар шеттік шарттармен бірге жалғыз шешімді бөліп шығарады
(2.53)

Қолданыстағы дифференциалдық теңдеуді дивергенттік емес түрде жазамыз:

0 ≤ x ≤ 1 облысында бір қалыпты тор
= xn + θh,
Дифференциалдық есепті айырымдық схемамен айырбастаймыз және дивергенттік емес
,

Айырымдық теңдеуді мына түрге түрлендіреміз:

,

мұндағы ai =ki - (ki+1 –
Құрылған айырымдық схема консервативті болмайтындығыны белгілі. Шын мәнінде
Айырымдық теңдеулер жүйесінің шешуі айқын түрде табылады. Торды
, 0 .
және x = нүктесіндегі қосымша
q шамасы - - тен +
2.7 – суретте k1 = 10, k2 = 1,
2.7 сурет- Консервативтік схема

2.5 Ағынға қарсы өзгертілген схеманың бір түрі

Ағынға қарсы айырымдық яғни донорлық ячейкалы (қуысты) схемада кеңістіктік
, мұндағы

(2.55)

(орталандырудың басқа да жолдары болуы мүмкін).

- дың мәндері былайша алынады:

Бұл схеманың консервативті және транспортивті болатындығын тексеру оңай. Бұл
Ескерту. Егер шамасын ячейканың қабырғаларында орта
Ағынға қарсы айырымдық бірінші схемамен салыстырғанда қарастырып отырған схемаға
тұрақты болатын жағдайды қарастырамыз, яғни , ал
немесе

(2.56)

бұл конвекцияны сипаттайтын мүше үшін екінші ретті дәлдікті береді.
2.6 Бір өлшемді сызықтық емес теңдеулерді шешу

Гидродинамиканың теңдеуі сияқты қарапайым сызықтық емес теңдеуді қарастырайық. Бұл
.
Бұл теңдеудің сол жағындағы бірінші және екінші қосылғыш оларға
.
Мұны Эйлер теңдеуі үшін модельдік ретінде қарастыруға болады, ол
(2.58) теңдеуін сызықтық емес толқындық теңдеу ретінде де түсіндіруге
(2.59)

Жалпы жағдайда u белгісізі де және F(u) функциясы да
(2.60)

мұндағы, жалпы жағдайда А=А (и) – Якоби матрицасы
Анықтама. (2.59) теңдеуінің тегіс шешімі деп, и функциясы
(2.59) теңдеуінің шешімдері шекті уақыт аралығында егер алғы шарттар
,
Бұл сияқты шектік шешімдердің үзілістегі шарттары мына Энтропия теңсіздігі
,

барлық ( және
Анықтама. (2.59) теңдеуінің, u функциясы үзілісті болуы мүмкін
(2.61)

және . Тегіс шешім әр уақытта онымен
Егер үзілістер энтропия шарты деп аталатын теңсіздікті қанағаттандырса, онда
Коши есебінің жалпылама мағынасындағы (2.61) жалғыз шешімін анықтау үшін
Интегралдық мағынада (2.61) теңдігінің (2.59) теңдеуі үшін Коши есебінің
энтропия теңсіздігін қанағаттандырушы функция мынау болады:

егер және егер
түзуі бұл теңдеу үшін үзіліс сызығы болады. Газогидродинамиканың осы
Бұл шешім өзінің құрылымы бойынша газогидродинамикадағы жеңілдету толқындары деп
Үлгілік (2.59) теңдеуінің жалпыланған шешімін құру, біріңғай орта теңдеуіндегі
Сұйықтың қозғалысын зерттегенде тұтқырлықтың және жылуөткізгіштіктің есептелуі, бірінші ретті
Үлгілік (2.59) теңдеуі үшін тұтқырлыққа сәйкестендіретін мүше кіргізу (2.59)
,
мұндағы – кіші параметр.

Анықтама.(2.59) теңдеуінің жалпыланған шешуін ''тұтқырлықты'' енгізудің көмегімен анықтау арнайы
Оның үстіне егер , онда (2.57) теңдеуінің
2.8 сурет - Бюргерс теңдеуі үшін үзілістің таралуы

Бюргерстің тұтқыр емес теңдеуіне қайта оралайын және осы теңдеудің
– кез келген үзіліссіз және үзіліссіз туындысы болатын функция
D – кез келеген (x,t) жазықтығындағы тікбұрышты облыс. Бұл
, .
Егер u және F функциялары үзіліссіз және үзіліссіз бірінші
Анықтама. Кез келген функциясы үшін (2.62)
(х, t) жазықтығындағы тікбұрышты D облысы, u
(2.63)

2.9 сурет - Үзіліс сызығындағы облыстың орналасуы

Соңғы интеграл D1 және D2 облыс
(2.62) - не сәйкес (2.63) – ға кіретін
.
Соңғы қатынас Бюргерс теңдеуі және
,

(штрихпен x бойынша диференциалдау белгіленген),
немесе .

Түпкілікті

(2.65)

яғни, үзілістің жылдамдығы оның сол жақтағы және оң жақтағы
Толқынды сейілту дыбыстан жоғары ағындарда екпінді толқындардан
,
(x, t) жазықтығындағы сипаттың түрі
Толқынның сейілу үлестірімі туралы
u=0, x≤0;
2.10 сурет - Толқын сейілуінің алғы шарты

2.11 сурет - Орталық толқын сейілуі жағдайындағы сипаттар

Сонымен, берілген бастапқы үлестерім ені уақытқа байланысты сызықты түрде
Біз дыбыстан жоғары газодинамикалық ағындарда жиі кездесетін екі
Алғы шарты болатын (2.60) теңдеуін қарастырайық.
(2.67)

коэффициенті нүктесінде тұрақтандырсақ, пайда болатын коэффициенттері
Сондықтан (2.67) схемасының орнықтылығынын күтуге болады. Егер диффузиялық есеп
Бірақта егер дифференциялдық есептің үзілісті шешімі болса, онда интегралдық
Сондықтан айырымдық схеманы құрғанда интегралдық сақталу заңын, ең болмағанда
(2.68)

- ге ұмтылғанда бізді қызықтыратын үзіліссіз жуық шешімді аламыз.

Тұтқырлық еленбейтін жағдайдағы (2.60) теңдеуін шешу үшін әлдеқандай схеманы
Лакспен және Вендрофпен үзілісті шешімдер үзілісті арнайы қарастырмай ақ
Екінші ретті дәлдікті әдіс, тұтқыр емес сұйықтық теңдеуін шешу
(2.69)

және туындылары қолданыстағы (2.69) теңдеуді
(2.70)

Егер осы өрнектерді (2.70) қойсақ және туындыларды орталықтық айырымдармен
(2.71)

мұндағы төмендегідей қатынаспен анықталады:

Мұндай схеманың орнықтылығы деп ұйғару арқылы
(2.72)

Лакс – Вендрофтың схемасының ені адымдық варианты алғашқы схемадан
,
(2.74)

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

(2.75)

(2.76)

,
Мұндағы, - еркін Т-ға байланысты функциялар.

(2.75) - (2.76) теңдеуді және шеттік шарты айқын емес
Бұл жағдайда температураға байланысты,
Айырымдық теңдеулер ішкі нүктелер үшін мына түрде болады:

,
ал шекаралық нүктелер үшін

,
,
мұндағы , сонымен бұл торлық функция үздіксіз
Жылу өткізгіштік коэффициенті ағымдардың торлық функцияларын
,

.

Екі айырымдық схеманы қарастырайық, -ге байланысты, бірінші
Квазисызықты теңдеуде коэффициенттері
Қарастырылған сызықты емес есепті (2.78) - (2.80) шешуге екі
. (2.81)

Мұндағы итерацияның нөмірі (2.81) теңдеуі уақыт
Нөлдік жуықтауды температураның мәндерін бұрынғы адымнан
.
Қорытынды

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

Андерсон Д., Таннехилл Дж., Плетчер Р. Вычислительная гидродинамика и
Бабушка Н., Витасен Э., Прагер М. Численные процессы решения
Бахвалов Н.С. Численные методы. – М.: Наука, 1973. –
Волков Е.А. Численные методы. – М.: Наука, 1982. –
Дьяченко В.Ф. Основные понятия вычислительной математики. – М.: Наука,
Годунов С.К., Рябенький В.С. Разностные схемы. – М.: Наука,
Калиткин Н.Н. Численные методы. – М.: Наука, 1978. –
Курант Р. Уравнения с частными производными. – М.: Мир,
Марчук Г.И. Методы вычислительной математики. - .: Наука. –
10. Марчук Г.И. Численные методы. – М.: Наука, 1989.
Пасканов В.М., Полежаев В.И., Чудов Л.А. Численное
Патанкар С.Численные методы решения задач теплообмена и динамики жидкости.
Пейре Р., Тейлор Т.д. Вычислительные методы в задач механики
Роуч П. Вычислительная гидродинамика. – М.: Мир, 1980. –
Самарский А.А. Теория разностных схем. – М.: Наука, 1983.
Самарский А.А., Николаев Е.С. Методы решения сеточных уравнений. –
Самарский А.А., Михаилов А.П. Математическое моделирование. Идеи. Методы. Примеры.
Шокин Ю.И. Метод дифференциального приближения. Новосибирск: Наука, 1979. –
Флетчер К. Вычеслительные методы в динамике жидкостей. – М.:
Хемминг Р.В. Численные методы. – М.: Наука, 1972. –
Яненко Н.Н. Метод дробных шагов решения многомерных задач математической
Годунов С.К. Уравнения математической физики. – Наука, 1979. –
Ковеня В.М., Яненко Н.Н. Метод расщепления в задачах газовой
Марчук Г.И. Методы расщпления. – М.: Наука, 1998. –
Тарасевич Ю.Ю. Математическое и компьютерное моделирование. Вводный курс.М.: Физматгиз.
Рихтмайер Р., Мортон К. Разностные методы решения краевых задач.-М.: