Байгалийн шинжлэл (физик, хими, биологи) болон
инженерийн шинжлэх ухаанд тавигдах бодлогууд, үзэгдэл процессын төлөв байдал,
хөдлөл зүйг гүн гүнзгий судалъя гэвэл эхлээд түүнийг математикийн хэлэн дээр
томьёолох шаардлагатай бөгөөд энэ нь ихэнхдээ ердийн болон тухайн уламжлалт
дифференциал тэгшитгэлийн хэлбэрээр бичигдэх нь олонтой. Уг үзэгдэл процесс
дахь компонентуудын харилцан үйлчлэлийг тооцвол шугаман бус тухайн уламжлалт
тэгшитгэл (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)-г хэтэрхий бага авахаас зайлсхийх доод босгыг
тавихад хүргэдэг. Тухайн тэгшитгэлд тодорхой физик агуулга бүхий параметрүүдийн тооцох нь элбэг бөгөөд ийм үед уг системийг параметерүүдын хувирах муж дахь төлөөллийн утга бүр дээр бодох учраас дээрх шаардлага бүр ч хурцаар тавигддаг. Бид энд жишээ болгож турблент урсгалын үрэлтийн коэффициент олох бодлого авч үзье. Турболент (хуйлрал) урсгалын үрэлтийн коэффициентын тооцоолох асуудал иргэний болон үйлдвэрлэл, гидравлик, химийн инженерчлэлийн практикт тавигддаг. Хий, шингэн, нефть тараах сүлжээний гол хэсэг нь хоолой бөгөөд түүгээр дамжих хуйларсан урсгалын шинж чанар нь уг хоолойн доторх гадаргуугийн гөлгөр байх болон арзгар байдлаас ихээхэн хамаардаг. Ийм урсгалын үрэлтийн хүчин зүйл (
трансцендент тэгшитгэл юм. Үүний шийдийг элементар функцүүдийн тусламжтайгаар илэрхийлж чадахгүй учраас өндөр эрэмбийн нийлэлттэй итерацийн аргуудаар олж болно.
(1) тэгшитгэлд Рейнольдсын тоо
гэсэн тэгш өнцөгт муж юм. Энэ мужийн цэг бүр дээр (1) тэгшитгэлийг бодох шаардлагатай учраас
өндөр эрэмбийн нийлэлттэй итерацийн арга хэрэглэх нь үр дүнтэй. Ердийн (1) тэгшитгэлийг ганц, хоёрхон удаа бодох
байсан бол ямар ч итерацийн арга хэрэглэж шийдийг олох нь чанарын ялгаагүй юм.
Гэтэл (1) бодлогын тавилын муж маш том учраас
тухайн тэгш өнцөгтийг жишээлбэл
Цэвэр математикт, тухайлбал алгебрт компьютерын тусламжтайгаар теоремын баталгаа хийдэг арга өргөн хэрэглэгдэх болсон. Бодолтыг алдаагүйгээр гүйцэтгэх символ-тооцоолол, компьютер-алгебр, символ дифференциалчлалын программууд (Matlab, Mathematica, Maple) шинэ чиглэл, дэвшлийг авчирч байна. Тооцон бодох математикт ойролцоо аргуудын нийлэлтийн теоремыг символ-тооцоолол ашиглан баталдаг шинэ чиглэл сүүлийн арваад жил эрчимтэй хөгжиж байна. Энэ чиглэлийн давуу тал нь нүсэр гар тооцооллыг халж завсрын бодолт, хувиргалтууд, хялбарчлалуудыг ямар ч алдаагүйгээр гүйцэтгэх замаар тухайн аргын алдааны тэгшитгэлийг аналитик байдлаар гаргадаг. Улмаар нийлэлтийн эрэмбийг хамгийн их утгад хүргэдэг параметрүүдийн сонголтыг олох явдал юм. Жишээлбэл J. R. Sharma, H. Arora нарын [4] ажилд
хоёр алхамт итерацийн арга
5-р эрэмбийн нийлэлтэй байхаар
биелж байвал
тухайн аргыг
үнэлгээг
гаргаж авна. Түүнчлэн
задаргааг
олж түүнийгээ 2-р шатанд ашиглах замаар
Нийлэлтийн эрэмбэ хамгийн их (манайд 5-р эрэмбийн) байх параметрүүдийн сонголт нь
буюу
болно.
Өөрөөр хэлбэл энэхүү сонголтод (2) итерац 5-р эрэмбийн нийлэлтэй байна гэсэн теорем батлагдлаа.
Тухайн
PDE-ийн хувьд түүний шийдийн онцлог, гол гол чанаруудыг хадгалсан тийм ойролцоо
ялгаварт схемийг ашиглах нь чухал юм. Энгийн заншсан ялгаварт схемүүдийн ихэнх
нь ийм “хадгалах” чанаргүй учраас сүүлийн жилүүдэд стандарт бус байгууламжтай
схемүүд ашиглах нь онол, практикийн хувьд чухал ач холбогдолтой, ирээдүйтэй
шинэ чиглэл болон хөгжиж байна. Энэ стандарт бус аргын хэрэглээний талбар өргөн
бөгөөд түүний зарим төлөөллүүдийг энд тоймлоё.
Байгалийн шинжлэх ухаан, инженерийн ухааны олон
үзэгдэл процессууд шугаман бус дифференциал тэгшитгэл, түүний системээр
загварчлагддаг. Хий, шингэний динамик, дулаан тархах процесс, шугаман бус
диссипатив (dissipative) систем дэх долгион тархах процессыг судлахад
Бюргерсийн тэгшитгэл (Burgers’ Equation) өргөн ашиглагддаг.
(3) нь шугаман бус
Коүл-Хопфийн
орлуулгаар (Cole-Hopf transform)
нэг хэмжээст
дулааны тэгшитгэлд шилждэг. Диффузийн тэгшитгэлийн
анхны нөхцөл хангасан шийд нь
юм.
(4), (6) -ийг ашиглаад (3) тэгшитгэлийн шийдийг олох боломжтой. (3)-ийн
анхны болон захын нөхцөлийг хангасан шийд нь
захын нөхцөл хангах гүйгч долгион шийдийг
хэлбэртэй
хайя. Энд
гэсэн ердийн
дифференциал тэгшитгэл гарах бөгөөд үүнийг
болж эндээс
Өөрөөр
хэлбэл цохилтын долгионы хурд нь зуурамтгай орчны параметр
мөрдөх бөгөөд үүнийг интегралчилахад
Зураг 1
багасах дутам долгионы фронт нь улам огцом
уналттай болж байна. Нэг загвар өөр өөр үзэгдэл процессыг тайлбарлахад
хэрэглэгддэгийн жишээ нь Бюргерсийн тэгшитгэл юм. (3) тэгшитгэл нь популяцийн динамик, зүйлс (species)-ийн
тархалтыг судлахад мөн ашиглагдана. Үүнд (3)-ийн шугаман бус гишүүн нь нягтралаас хамааралтай
шилжилтийг тооцох, харин баруун тал нь зүйлсийн санамсаргүй хөдөлгөөнийг заах
жишээтэй.Байгалийн шинжлэх ухааны олон үзэгдлийн математик загвар нь
реакци-диффузийн
тэгшитгэлүүдийн (reaction-diffusion equation) системээр илэрхийлэгддэг. Энд 
нь харгалзан урсгалын нэвчилтийн болон хурдны
коэффициент,
Биологи,
химийн олон үзэгдлүүдийн загварууд (11)-д багтана. Бид энд жишээ болгон гоц халдварт өвчний (эпидеми) математик
(SIR) загварыг авч үзье.
N(t)=S(t)+I(t)+R(t) нь t хугацаан дахь хүний тоо
S нь өвчинд мэдрэмтгий хүний тоо,
I нь халдвар авсан хүний тоо,
R нь халдвар аваад эдгэрсэн хүний тоо.
N мэдэгдэх учраас (12) -н эхний хоёр тэгшитгэлийг авч үзэх нь хангалттай. Уг систем нь
өвчлөөгүй үеийн тэнцвэрийн
цэг,
өвчин тархсан үеийн тэнцвэрийн
цэг.
Ийм хоёр тэнцвэрийн цэгтэй
Үүнд 
нь үржилийн
хурд. Халдварт өвчний тархалтын динамик нь дээрх тэнцвэрийн цэгүүдийн
тогтворжилтыг судлах уруу шилждэг. (12)-ийг
анхны болон захын нөхцөлтэйгөөр авч үздэг. (12) дахь параметрүүдийн физик утга нь:
нь харгалзан үхлийн хурд (1/μ нь амьдрах магадлал), өвчлөөд эдгэрэх хурд (1/ν нь халдвартай байх үе), шилжилтийн коэффициент (мэдрэмтгийгээс өвчлөлд шилжих) юм. Хүмүүс холилдсон хөдөлгөөнтэй байх учраас S, I нь хугацаа болон огторгуйн хувьсагч z-ээс хамаарсан функцүүд юм.
системээр
загварчлагдана. Энд
захын нөхцөл
тавигдана. Энд 2 янзын процесс явагдах боломжтой. Хий дэх автомат шаталт болж
болно. Энэ нь реакторыг тэсрэхэд хүргэх аюултай. Хоёрдох нь процесс хэвийн
явагдах болно. Тэгвэл процесс хэвийн явагдах нөхцөлийг ямар нөхцөлөөр хангах
вэ? гэдэг асуудал чухал юм. Үүнийг хариулахын тулд
захын нөхцөл
тавигдана. Энд 2 янзын процесс явагдах боломжтой. Хий дэх автомат шаталт болж
болно. Энэ нь реакторыг тэсрэхэд хүргэх аюултай. Хоёрдох нь процесс хэвийн
явагдах болно. Тэгвэл процесс хэвийн явагдах нөхцөлийг ямар нөхцөлөөр хангах
вэ? гэдэг асуудал чухал юм. Үүнийг хариулахын тулд
гэвэл 
байх шийдийг сонирхъё. Цилиндр хэлбэрийн реактор авч үзэхэд 
(
-реакторын хэмжээсний критик утга ) үед шийдгүй, 
үед 2 шийдтэй, үүний их температуртай шийд нь тогтворгүй (зураг 2-ийг харна уу).Зураг 2
Эсрэгээр сфер хэлбэрийн
реактор авахад n шийдтэй 
байх 
утгууд
олдоно. Энэ нь реакторын гүдгэр шугаман бус чанарыг ихэсгэхэд шийдийн төвөгтэй
төлвийг үүсгэж болохыг харуулж байна. (13) хэлбэрийн бодлогыг Gelfand-Bratu problem гэж нэрлэдэг.
Сүүлийн жилүүдэд амьдралын эрэлт хэрэгцээ, практик шаардлагаас улбаалан бутархай уламжлал бүхий PDE тэгшитгэл, түүний шийд ганц байх чанар, ийм тэгшитгэлийг ойролцоо бодох талаар судалгаа маш эрчимтэй явагдаж байгааг тэмдэглэх нь зүйтэй.
Тооцон бодох
математик, тоон аргууд байгалийн болон инженерийн олон асуудлуудыг шийдвэрлэх,
судлахад өргөн хэрэглэгдэж байгаагийн нэг үзүүлэлт нь олон улсын хэмжээнд
тогтмол хэвлэгдэн олны хүртээл болж байгаа сэтгүүлүүдийн жагсаалтаас бэлхнээ
харж болно.
|
Нэрс |
Impact Factor |
|
Journal of
Computational Physics |
2.845 |
|
Journal of
Computational Biology |
1.191 |
|
Computational
Geometry |
0.343 |
|
Journal of
Computational Finance |
0.758 |
|
Journal of
Computational Chemistry |
3.224 |
|
International
Journal for Numerical Methods in Biomedical Engineering |
2.082 |
|
International
Journal for Numerical Methods in Engineering |
2.746 |
|
International
Journal of Computational Methods in Engineering Science and Mechanics |
0.76 |
|
International
Journal for Numerical Methods in Fluids |
1.631 |
•
Тухайн
юм уу эсвэл түүнтэй ойролцоо шинжлэх ухаануудын болон математикийн гүнзгий
мэдлэг нь объект, үзэгдлийн математик загварыг босгоход илүү чухлаар шаардагдах
учраас тал талын мэргэжилтнүүдийн нягт хамтын ажиллагаа, интеграци, идэвх
санаачилга чухал ач холбогдолтой болж байна.
•
Академийн
хүрээлэнгүүд, их дээд сургуулиудын хамтын ажиллагааг өргөжүүлэх, хамтарсан
семинар уулзалт хийх, үүнд магистр, доктор оюутнуудыг идэвхтэй оролцуулах, бие
биесээ сонсож асуудлаа дэвшүүлж, энэ чиглэлд хүний нөөцөө дайчлах шаардлага
хурцаар тавигдаж байна.
•
Судалгааны
орчин үеийн математик аппарат, онолуудыг хэрэглэх, практикт нэвтрүүлэх асуудал
бас чухал байна. Жишээлбэл вейвлет анализын хэрэглээ маш өргөн юм. Дүрс, дуу,
дохиог боловсруулах (шахах, сэргээх, хадгалах, алсад дамжуулах), эмнэлгийн
практик нэвтрүүлэх, элдэв аппарат, техник хэрэгслийн оношилгоо тохиргоо хийх
зэрэгт вейвлет анализ өргөн ашиглагддаг. Энэ чиглэлээр олон улсад ашиглагдаж
байгаа бэлэн пакет-программыг интернэтээс татан өөр өөрийн судалгаанд ашиглаж
сурах нь шинэ дэвшлийг авчирна гэдэг нь эргэлзээгүй юм.
Ашигласан ном
[1] Mickens
Ronald E, Nonstandard Finite Difference Models of Differential Equations, World
Scientific, 1993.
[2] P. Praks,
D. Brki´c, Advanced Iterative Procedures for Solving the Implicit Colebrook Equation for Fluid
Flow Friction, Advances in Civil Engineering Vol.18, Article ID5452034, 18 pages.
[3] E.
Zeidler, Nonlinear functional analysis and its applications, I: Fixed-Point
Theorems, Springer-Verlag New York, 1986.
[4] J.R.
Sharma, H. Arora, Improved Newton-like methods for solving systems of nonlinear equations, SeMA
Journal 74 (2), 147-163.
[5] S.V.
Petrovskii, B.L. Li, Exactly Solvable Models of Biological Invasion, 2006.
[6] Тооцон бодох математик, Монгол улсын шинжлэх ухааны цуврал, 76-р боть, Улаанбаатар 2009.