ТООЦОН БОДОХ МАТЕМАТИКИЙН ХӨГЖЛИЙН ШИНЭ ЧИГ ХАНДЛАГА
Академич Т. Жанлав ШУА, Математик Тоон Технологийн Хүрээлэн
Оршил
Байгалийн шинжлэл (физик, хими, биологи) болон инженерийн шинжлэх ухаанд тавигдах бодлогууд, үзэгдэл процессын төлөв байдал, хөдлөл зүйг гүн гүнзгий судалъя гэвэл эхлээд түүнийг математикийн хэлэн дээр томьёолох шаардлагатай бөгөөд энэ нь ихэнхдээ ердийн болон тухайн уламжлалт дифференциал тэгшитгэлийн хэлбэрээр бичигдэх нь олонтой. Уг үзэгдэл процесс дахь компонентуудын харилцан үйлчлэлийг тооцвол шугаман бус тухайн уламжлалт тэгшитгэл (PDE)-д хүргэнэ.
Ийм шугаман бус тэгшитгэл, түүний системийн шийд нь аналитик байдлаар олдох нь маш ховор учраас шийдийг ямарваа аргаар ойролцоогоор олох асуудал өргөн тавигддаг. Ийм ойролцоо бодох арга, алгоритмыг байгуулах, түүний алдаа, нийлэлтийг тогтоох асуудлыг судалдаг математикийн салбар нь тооцон бодох математик юм.
Орчин үеийн компьютерын хөгжил нь тооцон бодох математикийн хөгжилд онцгой чухал ач холбогдолтой. Одон орон судлаачдад өндөр хүчин чадалтай дуран авиа (телескоп) ямар чухал байдаг, биологичдод өндөр хүчин чадалтай микроскоп ямар чухал байдгийн адилаар тооцон бодох математикт орчин үеийн хүчтэй компьютер тийм чухал үүрэг гүйцэтгэдэг юм. Тухайн судлагдаж байгаа объект, үзэгдэл процессын математик загварыг гарган түүнийг бодох арга, алгоритм боловсрогдсоны дараах шат нь түүний дагуу код программ бичиж компьютерын тоон туршилт, тооцоо судалгаа хийдэг бөгөөд туршилтын үр дүн хир зэрэг нарийвчлалтай, хир үнэмшилтэй гарч байгаа эсэхийг өөр аргаар гаргасныг, өөр багаж техник хэрэглэн гаргаж авсан туршилтын үр дүнтэй жиших замаар шинжилж зайлшгүй тохиолдолд математик загвараа эргэн харах, зохих засвар, коррекц оруулах шаардлагатай болдог цикл процесс юм.
Тухайлбал шугаман бус тухайн уламжлалт тэгшитгэлүүдийн системийг тоон аргаар бодоход асар их хурд, санах ой бүхий компьютер, түүний кластер шаардлагатай. Ялангуяа системийн тэгшитгэлийн тоо ихсэх дутам эдгээр шаардлагууд улам хурцаар тавигддаг болно. Тооцон бодох математикийн чиглэлд дэлхийн олон орны эрдэмтэд идэвхтэй, эрчимтэй ажиллаж өөрсдийн хувь нэмрээ оруулж байгааг тэмдэглэх хэрэгтэй. Манай улсад ч энэ чиглэлийн шавь сургууль хөгжиж олон залуу эрдэмтэд ажиллаж байгаа билээ. Тэдгээрийн бүтээлийн тоймыг Монгол улсын шинжлэх ухаан 76-р ботиос үзэж болно. Энэ чиглэлээр хэвлэгдэн гарч байгаа томоохон сэтгүүлүүдээс дурдвал:
Нэрс | Impact Factor |
Acta Numerica (Cambridge University Press) | 7.417 |
Journal of Approximation Theory (Elsevier) | 1.022 |
Numerical Methods for Partial Differential Equations (Wiley) | 1.633 |
SIAM, Journal on Numerical Analysis | |
Numerical Linear Algebra with Applications (Wiley) | 1.298 |
Journal of Computational and Applied Mathematics (Elsevier) | 1.883 |
Applied Mathematics and Computation (Elsevier) | 3.092 |
Numerical Algorithms (Springer) | 2.417 |
Computational and Applied Mathematics (Springer) | 1.26 |
Calcolo (Springer) | 1.981 |
BIT, Numerical mathematics | 1.670 |
Computational Mathematics and Mathematical Physics (Springer) | 0.774 |
International Journal of Computer Mathematics (Taylor France) | 1.196 |
Applied Numerical Mathematics (Elsevier) | 1.678 |
Advances in Computational Mathematics (Springer) | 1.638 |
(Энэ жагсаалт гүйцэд биш болохыг анхаарна биз ээ!)
Шугаман бус системийг бодох итерацийн аргууд, түүнд тавигдах шаардлагууд
Шугаман бус дифференциал тэгшитгэлийг ойролцоо бодоход дискретчлэл, шугамчлалын хоёр шатыг дамжих шаардлагатай. Төгсгөлөг ялгаврын болон элементийн аргууд, сплайн, вейвлет схемүүдийн тусламжтайгаар дискретчлэх процесс хийгдсэний дүнд тухайн тэгшитгэл нь төгсгөлөг хэмжээст огторгуйд тодорхойлогдсон ялгаварт тэгшитгэлүүдийн системд шилжинэ. Ийм шугаман бус системийг итерацийн аргуудаар алхам тутамд шугаман тэгшитгэлүүдийн систем бодох руу шилжүүлнэ. Бодолтын дүнд олсон ойролцоо шийд нь хэр зэрэг нарийвчлалтайгаар олдох нь хэрэглэсэн дискретчилэлийн аргын алдаа болон тогтворжилт, шугамчилсан аргын нийлэлтийн чанараас шууд хамаардаг. Иймээс шийдийг өндөр нарийвчлалтайгаар олох шаардлагатай тохиолдолд (шийдийн салаалалт, бифуркацид хүргэх параметрийн критик утгыг маш нарийн тооцоолох, цөмийн реакторын тооцоо, сансрын хөлгийн траекторын судалгаа тооцоо зэрэг) найдвартай, өндөр нарийвчлал, түргэн нийлэлтэй аргуудыг сонгох нь зүйтэй. Бодолтыг компьютерын тусламжтай гүйцэтгэх учраас аль болох бага санах ой, багтаамж, хугацаа шаардах арга алгоритм сонгох нь чухал. Энэ нь дискретчлэлийн алхам (h)-г хэтэрхий бага авахаас зайлсхийх доод босгыг
тавихад хүргэдэг. Тухайн тэгшитгэлд тодорхой физик агуулга бүхий параметрүүдийн тооцох нь элбэг бөгөөд ийм үед уг системийг параметерүүдын хувирах муж дахь төлөөллийн утга бүр дээр бодох учраас дээрх шаардлага бүр ч хурцаар тавигддаг. Бид энд жишээ болгож турблент урсгалын үрэлтийн коэффициент олох бодлого авч үзье. Турболент (хуйлрал) урсгалын үрэлтийн коэффициентын тооцоолох асуудал иргэний болон үйлдвэрлэл, гидравлик, химийн инженерчлэлийн практикт тавигддаг. Хий, шингэн, нефть тараах сүлжээний гол хэсэг нь хоолой бөгөөд түүгээр дамжих хуйларсан урсгалын шинж чанар нь уг хоолойн доторх гадаргуугийн гөлгөр байх болон арзгар байдлаас ихээхэн хамаардаг. Ийм урсгалын үрэлтийн хүчин зүйл ( -коэффициент)-ийг аль болох өндөр нарийвчлалтай тодорхойлох шаардлагатай байдаг. Энэ коэффициентыг тооцоолох тэгшитгэл нь![](blob:https://niitlel.mn/ebb9e646-bb07-473d-8362-7594f42633b4)
трансцендент тэгшитгэл юм. Үүний шийдийг элементар функцүүдийн тусламжтайгаар илэрхийлж чадахгүй учраас өндөр эрэмбийн нийлэлттэй итерацийн аргуудаар олж болно.
(1) тэгшитгэлд Рейнольдсын тоо болон харьцангуй арзгар байдлыг тодорхойлох параметр гэсэн хоёр параметр орсон байна. Эдгээрийн практикт хэрэглэгдэх утгын муж нь
![](https://niitlel.mn/wp-content/uploads/2020/11/image.png)
гэсэн тэгш өнцөгт муж юм. Энэ мужийн цэг бүр дээр (1) тэгшитгэлийг бодох шаардлагатай учраас өндөр эрэмбийн нийлэлттэй итерацийн арга хэрэглэх нь үр дүнтэй. Ердийн (1) тэгшитгэлийг ганц, хоёрхон удаа бодох байсан бол ямар ч итерацийн арга хэрэглэж шийдийг олох нь чанарын ялгаагүй юм. Гэтэл (1) бодлогын тавилын муж маш том учраас тухайн тэгш өнцөгтийг жишээлбэл хувааж хуваалтын жижиг тэгш өнцөгтүүдээс нэг нэг цэг аваад (ийм цэгийн тоо 65536) түүний хувьд (1) тэгшитгэлийг бодно гэвэл (1) -ийг 65536 удаа бодох шаардлагатай болж байна. Тэгэхээр энд өндөр эрэмбийн нийлэлтэй итерацийн арга хэрэглэх нь машины ресурс, бодолтын цаг хугацааг хэмнэх чухал ач холбогдолтой.
Символ тооцоолол, хэрэглээ
Цэвэр математикт, тухайлбал алгебрт компьютерын тусламжтайгаар теоремын баталгаа хийдэг арга өргөн хэрэглэгдэх болсон. Бодолтыг алдаагүйгээр гүйцэтгэх символ-тооцоолол, компьютер-алгебр, символ дифференциалчлалын программууд (Matlab, Mathematica, Maple) шинэ чиглэл, дэвшлийг авчирч байна. Тооцон бодох математикт ойролцоо аргуудын нийлэлтийн теоремыг символ-тооцоолол ашиглан баталдаг шинэ чиглэл сүүлийн арваад жил эрчимтэй хөгжиж байна. Энэ чиглэлийн давуу тал нь нүсэр гар тооцооллыг халж завсрын бодолт, хувиргалтууд, хялбарчлалуудыг ямар ч алдаагүйгээр гүйцэтгэх замаар тухайн аргын алдааны тэгшитгэлийг аналитик байдлаар гаргадаг. Улмаар нийлэлтийн эрэмбийг хамгийн их утгад хүргэдэг параметрүүдийн сонголтыг олох явдал юм. Жишээлбэл J. R. Sharma, H. Arora нарын [4] ажилд
хоёр алхамт итерацийн арга 5-р эрэмбийн нийлэлтэй байхаар параметрүүдийн сонголтыг олсон теоремын баталгааг символ-тооцооллын тусламжтай яаж гүйцэтгэхийг авч үзье. -р дөхөлтийн алдааг гэж тэмдэглэе. Үүний хувьд
биелж байвал тухайн аргыг эрэмбийн хурдтай нийлж байна гэдэг. -уудын Тейлорын задаргааг эхний тэгшитгэлд ашиглах замаар
үнэлгээг гаргаж авна. Түүнчлэн -уудын Тейлорын задаргааг символ-тооцооллын бодолтууд хийснээр дараах илэрхийллийн
задаргааг олж түүнийгээ 2-р шатанд ашиглах замаар дөхөлтийн алдааны тэгшитгэлийг бичвэл
Нийлэлтийн эрэмбэ хамгийн их (манайд 5-р эрэмбийн) байх параметрүүдийн сонголт нь
буюу
болно. Өөрөөр хэлбэл энэхүү сонголтод (2) итерац 5-р эрэмбийн нийлэлтэй байна гэсэн теорем батлагдлаа.
Тухайн PDE-ийн хувьд түүний шийдийн онцлог, гол гол чанаруудыг хадгалсан тийм ойролцоо ялгаварт схемийг ашиглах нь чухал юм. Энгийн заншсан ялгаварт схемүүдийн ихэнх нь ийм “хадгалах” чанаргүй учраас сүүлийн жилүүдэд стандарт бус байгууламжтай схемүүд ашиглах нь онол, практикийн хувьд чухал ач холбогдолтой, ирээдүйтэй шинэ чиглэл болон хөгжиж байна. Энэ стандарт бус аргын хэрэглээний талбар өргөн бөгөөд түүний зарим төлөөллүүдийг энд тоймлоё. 1. Бюргерсийн тэгшитгэл Байгалийн шинжлэх ухаан, инженерийн ухааны олон үзэгдэл процессууд шугаман бус дифференциал тэгшитгэл, түүний системээр загварчлагддаг. Хий, шингэний динамик, дулаан тархах процесс, шугаман бус диссипатив (dissipative) систем дэх долгион тархах процессыг судлахад Бюргерсийн тэгшитгэл (Burgers’ Equation) өргөн ашиглагддаг.
(3) нь шугаман бус гишүүн, зуурамтгайн коэффициент бүхий диффузийн (diffusion) гишүүнийг агуулдгаараа Навье-Стоксын тэгшитгэлтэй төстэй, түүний хялбаршуулсан хэлбэр гэж үздэг. (3) нь аналитик шийд нь мэдэгддэг цөөхөн шугаман бус PDE-үүдийн нэг учраас PDE тэгшитгэлүүдийг бодох тоон аргуудыг үнэлэх тест тэгшитгэл болдгоороо чухал ач холбогдолтой. (3) нь
Коүл-Хопфийн орлуулгаар (Cole-Hopf transform)![](blob:https://niitlel.mn/59628aa4-0891-482c-8512-7d8cce2a7c3d)
нэг хэмжээст дулааны тэгшитгэлд шилждэг. Диффузийн тэгшитгэлийн
анхны нөхцөл хангасан шийд нь
юм.(4), (6) -ийг ашиглаад (3) тэгшитгэлийн шийдийг олох боломжтой. (3)-ийн
анхны болон захын нөхцөлийг хангасан шийд нь![](blob:https://niitlel.mn/15c8d9ce-5e63-459b-8381-eb62cf443c44)
![](blob:https://niitlel.mn/e0b6ee6e-423f-4c5d-91d8-ac42eb07db67)
захын нөхцөл хангах гүйгч долгион шийдийг
хэлбэртэй хайя. Энд нь долгионы хурд юм. (8)-ийг (3)-д орлуулахад
гэсэн ердийн дифференциал тэгшитгэл гарах бөгөөд үүнийг мужаар интегралчилахад
болж эндээс
Өөрөөр хэлбэл цохилтын долгионы хурд нь зуурамтгай орчны параметр - ээс биш харин зөвхөн захын утгаас хамаарч байна. (9)-аас
мөрдөх бөгөөд үүнийг интегралчилахад
Зураг 1![](blob:https://niitlel.mn/034a0233-60cc-4dad-9e54-7ba8095012e5)
багасах дутам долгионы фронт нь улам огцом уналттай болж байна. Нэг загвар өөр өөр үзэгдэл процессыг тайлбарлахад хэрэглэгддэгийн жишээ нь Бюргерсийн тэгшитгэл юм. (3) тэгшитгэл нь популяцийн динамик, зүйлс (species)-ийн тархалтыг судлахад мөн ашиглагдана. Үүнд (3)-ийн шугаман бус гишүүн нь нягтралаас хамааралтай шилжилтийг тооцох, харин баруун тал нь зүйлсийн санамсаргүй хөдөлгөөнийг заах жишээтэй.
2. Реакци-диффузийн тэгшитгэлүүдийн систем
Байгалийн шинжлэх ухааны олон үзэгдлийн математик загвар нь
![](https://niitlel.mn/wp-content/uploads/2020/11/image-1.png)
реакци-диффузийн тэгшитгэлүүдийн (reaction-diffusion equation) системээр илэрхийлэгддэг. Энд нь хамаарах хувьсагч, өөрөөр хэлбэл эерэг утгатай байх физик хувьсагчууд юм.
нь харгалзан урсгалын нэвчилтийн болон хурдны коэффициент, нь харилцан үйлчлэл урвалдуулах гишүүд бөгөөд сөрөг бус байна. Ө.х
![](https://niitlel.mn/wp-content/uploads/2020/11/image-3.png)
Биологи, химийн олон үзэгдлүүдийн загварууд (11)-д багтана. Бид энд жишээ болгон гоц халдварт өвчний (эпидеми) математик (SIR) загварыг авч үзье.
![](https://niitlel.mn/wp-content/uploads/2020/11/image-2.png)
N(t)=S(t)+I(t)+R(t) нь t хугацаан дахь хүний тоо
S нь өвчинд мэдрэмтгий хүний тоо,
I нь халдвар авсан хүний тоо,
R нь халдвар аваад эдгэрсэн хүний тоо.
N мэдэгдэх учраас (12) -н эхний хоёр тэгшитгэлийг авч үзэх нь хангалттай. Уг систем нь
өвчлөөгүй үеийн тэнцвэрийн цэг,
өвчин тархсан үеийн тэнцвэрийн цэг.
Ийм хоёр тэнцвэрийн цэгтэй
![](https://niitlel.mn/wp-content/uploads/2020/11/image-4.png)
Үүнд нь үржилийн хурд. Халдварт өвчний тархалтын динамик нь дээрх тэнцвэрийн цэгүүдийн тогтворжилтыг судлах уруу шилждэг. (12)-ийг
![](https://niitlel.mn/wp-content/uploads/2020/11/image-5.png)
анхны болон захын нөхцөлтэйгөөр авч үздэг. (12) дахь параметрүүдийн физик утга нь:
нь харгалзан үхлийн хурд (1/μ нь амьдрах магадлал), өвчлөөд эдгэрэх хурд (1/ν нь халдвартай байх үе), шилжилтийн коэффициент (мэдрэмтгийгээс өвчлөлд шилжих) юм. Хүмүүс холилдсон хөдөлгөөнтэй байх учраас S, I нь хугацаа болон огторгуйн хувьсагч z-ээс хамаарсан функцүүд юм.
![](https://niitlel.mn/wp-content/uploads/2020/11/image-6.png)
системээр загварчлагдана. Энд нь урвалын хурд (давтамж) нь эерэг тогтмолууд, 2-р тэгшитгэл нь диффузийн процессыг тайлбарлана. Реактор гадна орчноос тусгаарлагдсан гэвэл
![](https://niitlel.mn/wp-content/uploads/2020/11/image-8.png)
захын нөхцөл тавигдана. Энд 2 янзын процесс явагдах боломжтой. Хий дэх автомат шаталт болж болно. Энэ нь реакторыг тэсрэхэд хүргэх аюултай. Хоёрдох нь процесс хэвийн явагдах болно. Тэгвэл процесс хэвийн явагдах нөхцөлийг ямар нөхцөлөөр хангах вэ? гэдэг асуудал чухал юм. Үүнийг хариулахын тулд гэвэл стационар тохиолдол болж дээрх систем нь нэг тэгшитгэл
![](https://niitlel.mn/wp-content/uploads/2020/11/image-7.png)
захын нөхцөл тавигдана. Энд 2 янзын процесс явагдах боломжтой. Хий дэх автомат шаталт болж болно. Энэ нь реакторыг тэсрэхэд хүргэх аюултай. Хоёрдох нь процесс хэвийн явагдах болно. Тэгвэл процесс хэвийн явагдах нөхцөлийг ямар нөхцөлөөр хангах вэ? гэдэг асуудал чухал юм. Үүнийг хариулахын тулд гэвэл стационар тохиолдол болж дээрх систем нь нэг тэгшитгэл
![](https://niitlel.mn/wp-content/uploads/2020/11/image-7.png)
болно. M_μ=√μ M, 0∈M (μ нь реакторын хэмжээсийг заах параметр) гэсэн масштаб оруулах замаарбодлого тавигдана. Химийн инженерчлэлийн практикт f(u)=A"exp"(Bu), A,B эерэг тогтмолууд байдлаар ихэвчлэн авч үздэг байна. Боломжит температурын дээд заагийг
гэвэл
байх шийдийг сонирхъё. Цилиндр хэлбэрийн реактор авч үзэхэд
(
-реакторын хэмжээсний критик утга ) үед шийдгүй,
үед 2 шийдтэй, үүний их температуртай шийд нь тогтворгүй (зураг 2-ийг харна уу).
Зураг 2
![](https://niitlel.mn/wp-content/uploads/2020/11/image-9.png)
Эсрэгээр сфер хэлбэрийн реактор авахад n шийдтэй байх
утгууд олдоно. Энэ нь реакторын гүдгэр шугаман бус чанарыг ихэсгэхэд шийдийн төвөгтэй төлвийг үүсгэж болохыг харуулж байна. (13) хэлбэрийн бодлогыг Gelfand-Bratu problem гэж нэрлэдэг.
Сүүлийн жилүүдэд амьдралын эрэлт хэрэгцээ, практик шаардлагаас улбаалан бутархай уламжлал бүхий PDE тэгшитгэл, түүний шийд ганц байх чанар, ийм тэгшитгэлийг ойролцоо бодох талаар судалгаа маш эрчимтэй явагдаж байгааг тэмдэглэх нь зүйтэй.
Сэтгэгдэл үлдээх