Головна Природознавство
СЕТОЧНО-ХАРАКТЕРИСТИЧНІ ЧИСЕЛЬНІ МЕТОДИ
|
|
|||||
ПЕРЕДМОВА ДО ДРУГОГО ВИДАННЯВидана в 1988 р наукова монографія «сіткових-характеристичні чисельні методи» давно стала бібліографічною рідкістю, а в деяких пропозиціях про продаж оголошується навіть «антикварної». На її основі в Московському фізико-технічному інституті одним з авторів всі ці роки читається семестровий курс лекцій з обчислювальним практикумом «Нелінійні обчислювальні процеси» для ряду спеціалізацій факультетів управління і прикладної математики і аерофізікі і космічних досліджень. Відеолекції з даного курсу доступні також в інтернеті 1і мають багатотисячну аудиторію. У МФТІ ці лекції доповнюють початковий (річний) курс обчислювальної математики і фундаментальний (також річний) курс теорії різницевих схем. У зв'язку з цим представляється корисним перевидання зазначеної книги, доповнене коротким оглядом вельми істотних досягнень за минулі без малого 30 років в розвитку обчислювальної математики, в тому числі «сеточно-характеристичного» напряму, пов'язаного з чисельним рішенням нелінійних систем рівнянь гіперболічного типу і ряду інших підходів . Рівняння і системи гіперболічного типу складають значну частину математичних моделей, що використовуються для вирішення найрізноманітніших прикладних задач. З безлічі можливих прикладів відзначимо рівняння Ейлера для нестаціонарних і надзвукових стаціонарних течій стиснення газу, нестаціонарні рівняння Максвелла і різні їх комбінації з рівняннями Ейлера, зокрема рівняння магнітогазо- динаміки і ряду інших математичних моделей фізики плазми, нестаціонарні рівняння механіки деформівних твердих тіл для ряду реологічних моделей і геометрій (як в багатовимірних моделях суцільного середовища, так і в ряді теорій стрижневих і оболонкових конструкцій), нестаціонарні одно- і д вумерние рівняння течії нестисливої рідини в полі сили тяжіння в наближенні дрібної води, нестаціонарні рівняння течії нестисливої рідини в пружно деформуються трубах (зокрема, описують пульсові течії крові в великих кровоносних судинах і повітря в трахейно-бронхіальному тракті), динамічні рівняння переносу енергії в електроенергетичних і радіотехнічних системах, інтенсивного автодорожнього руху в мегапо- 1 URL: http://fvl.fizteh.ru/courses/kholodov/; http://univertv.ru/video/ matematika / matematicheskie_metody_i_modeli /? mark = sciencel; http://www.intuit. ru / goods_store / video / 5663. лисиць, інтенсивних інформаційних потоків в телекомунікаційних і комп'ютерних мережах і багато інших. Ще більш широкою є область математичних моделей, в яких гіперболічна частина оператора є суттєвою складовою, багато в чому визначальною властивості шуканих рішень і використовуваних чисельних методів (наприклад, рівняння Нав'є-Стокса та ін.). З моменту появи перших чисельних методів для рівнянь гіперболічного типу їх число досягло важко піддається оцінці величини і безперервно поповнюється новими розробками. З одного боку, це є свідченням важливості даного класу задач в різних додатках, однак така велика кількість свідчить також про те, що універсального, що задовольняє всім вимогам користувачами вимогам, методу, мабуть, не існує (принаймні, він до цих пір не знайдене) . Огляди але даної проблеми (існують навіть огляди оглядових робіт) і монографії включають посилання на багато сотень оригінальних робіт, що не дозволяє провести тут навіть побіжний їх аналіз. Тому обмежимося аналізом лише найбільш значущих ідей, запропонованих при розробках різних методів і підходів, відсилаючи читача за більш детальним описом, наприклад, до монографій [92, 66, 57] і численним оригінальним роботам, посилання на які в них наведені. Відмінною рисою рівнянь гіперболічного типу є кінцева швидкість розповсюдження збурень (хвиль) в області інтегрування і існування характеристичних різноманіть - характеристичних ліній і поверхонь (обмежують області залежності і впливу рішень), на яких розмірність вихідної системи рівнянь в приватних похідних (число незалежних змінних) можна зменшити на одиницю. Основною тут є робота [188], в якій, хоча і в окремому випадку рівнянь баротропного газу, були детально розглянуті характеристичні властивості цих рівнянь, в тому числі вперше введені в розгляд їх перші інтеграли, названі згодом інваріантами Рімана. З інших робіт за характеристичними властивостями рівнянь газової динаміки і більш загальних систем рівнянь гіперболічного типу в якості прикладу відзначимо: [94, 1611, а також монографії та підручники [83, 59,112, 92, 20, 23, 171, 58] та ін. Перші чисельні методи для рівнянь і систем гіперболічного типу були засновані на використанні характеристичних властивостей цих рівнянь, як для апроксимації самих рівнянь, так і для побудови різницевої сітки ( прямі методи характеристик , перший з них був запропонований ще в кінці XIX століття Массо [176]) . По суті це були перші різницеві схеми на неструктурованих, адаптуються до вирішення сітках, настільки популярні і інтенсивно розробляються в даний час. Характеристичні властивості гіперболічних рівнянь активно використовувалися також для доказу збіжності чисельного рішення до точному - [156] - відоме умова стійкості Куранта-Фрідріхса-Леві і число Куранта), в тому числі в нелінійному випадку [157]. Подібні чисельні методи (для двовимірних задач газової динаміки - [37, 145 [и др., Для тривимірних - [154, 153, 180, 46, 39, 64, 11] та ін.) Були придатні для побудови тільки безперервних, а часто ще й досить гладких рішень і вимагали явного виділення особливостей (точок, ліній або поверхонь розривів в рішеннях) з постановкою на розривах відповідних граничних умов, що сильно ускладнювало розрахунковий алгоритм в цілому, особливо в багатовимірному випадку і при виникненні нових розривів, коли потрібно їх ідентифікація і введення в розрахункову схему. Детальний огляд цього напрямку можна знайти, наприклад, в оглядовій роботі [+1461. Подальшим розвитком характеристичного напрямки були зворотні , або сеточно-характеристичні методи , коли стали використовувати, як в звичайних сіткових методах, регулярну разностную сітку, але апроксимувались НЕ вихідні рівняння, а спеціальні їх лінійні комбінації - умови спільності уздовж характеристичних напрямків з тієї чи іншої інтерполяцією параметрів в точках перетину характеристик (в багатовимірному випадку ліній перетину характеристичних і координатних поверхонь, біхарактерістік і т.п.) з площиною , На якій дані вже відомі (початкової або раніше розрахованої): [157, 167] - схема Карлсона [61, 3, 65, 54, 56] і ряд інших відомих чисельних методів. Як і прямі методи характеристик, ці методи не були консервативними, і при їх використанні потрібно явне виділення сильних розривів. Такі підходи, як і раніше необхідні для розрахунку граничних точок в задачах з більш складними, ніж найпростіші (періодичні і т.п.), граничними умовами. Іншою особливістю рівнянь і систем гіперболічного типу є можливість появи розривних рішень, в нелінійному випадку навіть при гладких (в тому числі аналітичних) крайових умовах - відома градиентная катастрофа. Для лінійних рівнянь і систем джерелом розривів можуть бути розривні функції в початкових або граничних умовах, а також неузгодженість початкових і граничних умов. В альтернативних методів характеристик «звичайних» різницевих схемах (схемах на регулярних, що не характеристичних сітках і без приведення вихідних рівнянь до характеристическому увазі), в тому числі здатних відтворювати розривні чисельні рішення без явного виділення розривів (наприклад, методи [172, 19] та ін .), там, де вдасться виділити особливості, також краще використовувати явне виділення розривів, що забезпечує можливість використання невеликого числа сіткових вузлів для отримання прийнятної точності чисельного вирішене я. Коли продуктивність ЕОМ була порівняно невеликою, такі методи з явним виділенням розривів були цілком природними, а часто і єдино можливими, так як це дозволяло навіть на слабкій обчислювальної базі вирішувати досить складні класи задач. Як приклад можна вказати відомі в обчислювальної газодинаміки і механіки деформівних тіл сіткові методи [5, 21, 32, 52, 63, 142, 147, 66] та ін. З ростом продуктивності ЕОМ все більшого поширення набули так звані методи наскрізного рахунку (однорідні різницеві схеми), в яких розривні чисельні рішення розраховувалися без явного виділення розривів. В таких методах для рівнянь і систем гіперболічного типу, крім загальних для всіх рівнянь з приватними похідними питань апроксимації і стійкості різницевих схем (див., Наприклад, [96, 24, 97, 187] та ін.), Не менш важливу роль відіграють такі властивості чисельних методів, як консервативність і монотонність різницевих схем. Консервативні схеми (першими такими методами, мабуть, були [183, 172, 31. 33, 1111 та ін.) Засновані на апроксимації дівергентних рівнянь, записаних в інтегральної формі, і забезпечують виконання дискретних аналогів відповідних законів збереження, що і дозволяє чисельно знаходити розривні рішення без явного виділення розривів. У разі дуже сильних розривів як додатковий елемент в консервативних різницевих схемах використовують точне ( автомодельного) рішення задачі про розпад довільного розриву, якщо його вдається знайти (метод [19] і його численні наступні зміни - [50,134] і ін.). У деяких випадках виникає необхідність наділення різницевих схем аналогом деяких додаткових властивостей вихідних рівнянь. Наприклад, в рівняннях Ейлера поряд зі збереженням повної енергії важливо забезпечити правильний баланс окремих її видів при істотному - на порядки - їх відміну. Приклади таких повністю консервативних різницевих схем запропоновані в роботі [84, 9], а подальші узагальнення ( « вариационно-різницеві » схеми, метод опорних операторів) - в [26, 76 | і ін. Певну проблему для обчислювачів представляють ситуації, коли в вигляді законів збереження записується лише частина рівнянь, як, наприклад, в одіоскоростпих мпоготемпера- туріих моделях плазми ( частково консервативні системи - див. [57]), нелінійні гіперболічні рівняння і системи з розривними матрицями Якобі - см. [110,40] та ін., виконання в чисельному методі аналогів групових властивостей диференціальних рівнянні [150] та ін. Монотонні схеми (за іншою термінологією - схеми з позитивною аппроксимацией, мажорантності схеми) в нестрогому визначенні - це схеми, які переводять монотонний профіль рішення на якомусь часовому шарі в також монотонний профіль рішення на наступному (прилеглому) часовому шарі, як це має місце в точній вирішенні, по крайней мере, для лінійного рівняння переносу або для інваріантів Рімана в разі систем рівнянь. Найбільш сильно відмінності в численних рішеннях при розрахунках по монотонним і немонотонним схемами проявляються на грубих просторових сітках і при великих часах розрахунку. Існують різні визначення і критерії монотонності різницевих схем, а також способи регуляризації (монотонізаціі) розривних чисельних рішень, початок яким покладено введенням в немонотонні схеми з порядком апроксимації вище першого штучної в'язкості [183, 99. 60] та ін., « Згладжуванням » ( фільтрацією ) високочастотних осциляцій в численних рішеннях [163, 144] та ін. Для рівнянь і систем гіперболічного типу поняття про монотонних схемах вперше введено К. Фрідріха [162]. З інших критеріїв монотонності слід зазначити критерії [19, 174,164]. Критерій Фрідріхса є найбільш універсальним і сформульований як неотрицательность скалярних або неотрицательная визначеність матричних коефіцієнтів різницевих схем. Цей же критерій застосуємо для параболічних і еліптичних рівнянь і систем - для еліптичних систем у разі попарно комутуючих матричних коефіцієнтів [123, 118, 120]. У разі лінійного рівняння переносу (найпростішого з гіперболічних рівнянь), в роботі [191 для явних двошарових сіткових шаблонів доведена відома теорема Годунова про неможливість побудови лінійних монотонних схем з порядком апроксимації вище першого. В роботі [122] теорема С.К. Годунова про монотонних але Фрідріхса різницевих схемах була узагальнена на випадок довільних (в тому числі багатошарових і неявних) сіткових шаблонів. До таких монотонним по Фрідріхса схемами належать вже обговорювалися вище методи першого порядку апроксимації: [157, 167, 61, 3, 65, 541 та ін. Ряд монотонних по Фрідріхса різницевих схем першого порядку апроксимації одночасно має властивість консервативності. Зокрема, такими властивостями володіють схема [172], а також схема [19] і її різні лінеаризації, першими з яких були явні і неявні консервативні схеми з роботи [122], придатні для загальних систем рівнянь гіперболічного типу. Однак в разі розривів дуже великої інтенсивності лінеаризовані схеми вимагають сильного подрібнення кроку інтегрування по гіперболічної змінної. Так як в цьому випадку тимчасової крок інтегрування обмежується тепер швидкістю руху поверхні розриву (ударної або магнітозвукових хвилі і т.п.), яка істотно більше швидкості поширення малих лінійних збурень. Пізніше були побудовані багато інших лінеаризовані схеми Годунова, з яких слід зазначити схему [189] метод, заснований на наближеному вирішенні задачі про розпад довільного розриву (линеаризованной задачі Рімана), і використовує при побудові рішень цього завдання тільки розриви. Цей метод точно описує співвідношення на одиночних розривах і швидкості їх переміщення, проте хвилі розрідження через це може перетворювати в нефі- зічяие ударні хвилі розрідження. Цей дефект був усунений в роботі [101]. Спочатку цей метод був побудований для рівнянь газової динаміки і згодом поширений на МГД-рівняння, рівняння мілкої води та ін. Метод [185] також заснований на наближеному вирішенні задачі про розпад довільного розриву, але використовує при побудові рішень цього завдання тільки хвилі стиснення і розрідження , тому в розрахунках ударні хвилі замінює ріманово хвилями стиснення. Він також має реалізації для рівнянь газової динаміки, МГД-рівнянь, рівнянь мілкої води. Обмеження, які випливають з теореми Годунова, вдалося подолати переходом до нелінійним різницевим схемам. До таких нелінійним монотонним схемш відносяться гібридні схеми (схеми зі змінним порядком апроксимації - першим поблизу розривів і більш високим в області гладкого рішення - [115, 27, 81, 103] та ін.), Різні варіанти методу корекції потоків (по суті також гібридні схеми) - [152, 174, 70] та ін. і подальший їх розвиток - TVD схеми [164], різні «ква- зімонотонние » схеми [17, 25] та ін. На цьому шляху вдалося істотно підвищити реальну точність чисельних методів і дозволу розривів в численних рішеннях, спираючись на помічений ще Р.Н. Федоренко «експериментальний» факт, що полягає в тому, що в таких схемах ширина області розмазування розривів така ж, як в базових схемах високого порядку апроксимації (істотно менше, ніж у базових монотонних схемах першого порядку апроксимації), але при цьому відсутні нефізичні осциляції в рішеннях . Величезне число існуючих такого роду різницевих схем відрізняються вибором різних сіткових шаблонів, базових різницевих схем (монотонних першого порядку апроксимації поблизу розривів і немонотонність з високим порядком апроксимації в області гладкого рішення) і перемикачів між цими базовими схемами. Оригінальний підхід був запропонований в роботі [8], в ньому гарне дозвіл в розривних численних рішеннях досягається використанням двох (хоча можна використовувати і більше їх число) монотонних схем першого порядку апроксимації. Найбільш загальний підхід до побудови подібних схем з використанням лінійних просторів невизначених коефіцієнтів запропонований в [122, 121] і широко використовується для класифікації і аналізу монотонних по Фрідріхса і близьких до них найменш коливальних різницевих схем. Хоча в разі розривних рішень говорити про схемах наскрізного рахунку з порядком апроксимації вище першого не цілком коректно (під цим розуміють, природно, порядок апроксимації схеми па рішеннях відповідної гладкості, а потім використовують їх для наскрізного рахунку негладких, в тому числі розривних рішень), спроби побудови таких схем шляхом комбінації (в залежності від поведінки шуканого чисельного рішення) двох або більше схем з порядком апроксимації у кожної з них вище першого виявилися цілком успішними [53, 6], і узагальнення цієї схеми а загальні системи гіперболічного типу 116] та ін., привели до створення різницевих схем високого дозволу (UNO-, ENO-, WENO- схем [165] та ін.). Опис численних прикладів побудови цих схем читач може знайти, наприклад, в оглядовій монографії [57], а найбільш загальний підхід до побудови цього класу схем з використанням аналізу різницевих схем в просторах сіткових функцій в [170, 132, 128, 133] і в ряді інших робіт. Характеристичні властивості гіперболічних рівнянь зберігаються і в наслідках вихідних рівнянь, одержуваних дифференцированием вихідних рівнянь ( продовжених або розширених системах рівнянь): [158, 921 і ін. Такі продовжені системи використовуються в явній ([29] та ін.) Або неявній формі для побудови ефективних чисельних методів схем високого порядку апроксимації на нерас- Ширяєв сіткових шаблонах, в тому числі монотонних різницевих схем подібного типу ([132] та ін.). Іншим, більш поширеним прийомом для побудови схем високого порядку апроксимації без явного використання багатоточкових сіткових шаблонів (по аналогії з двоточковими методами Рунге-Кутта для систем звичайних диференціальних рівнянь), є методи типу предиктор-коректор (схеми другого порядку апроксимації [173, 175, 193 ] і ін., схеми третього порядку [191, 95, 194] та ін.), компактні різницеві схеми ([ІЗ] і ін.) і ряд інших підходів. Використання продовжених систем дозволяє також позбутися від виникає в ряді випадків немонотонності « чисельних похідних» в монотонних за рішенням різницевих схемах. Ще один підхід боротьби з зазначеним явищем немонотонності «чисельних похідних» в розривних рішеннях (а іноді і в гладких, але з різкою зміною похідних, рішеннях) запропонований, наприклад, в [77] ( «схеми з сильною монотонністю»). Для багатовимірних рівнянь і систем гіперболічного типу найбільш активно використовуються підходи, пов'язані з розщепленням вихідного багатовимірного оператора на одномірні [98, 52, 28, 49, 691 та ін., Спочатку запропоновані для параболічних рівнянь [186, 160, 149]. Деякі підходи до побудови різницевих схем на неструктурованих сітках , які використовуються для вирішення завдань в складних, в тому числі багатозв'язних, областях інтегрування, на адаптуються до шуканого чисельному рішенню різницевих сітках , а також для побудови самих генераторів сіток читач може знайти в роботах [22, 142 , 127, 190, 184, 87, 10, 166, 178, 105, 18, 38] та ін. Слід зазначити, що в процесі еволюції уявлень про те, якими властивостями повинен володіти ефективний чисельний метод для рівнянь гіперболічного типу вже вдалося подолати оч ень багато труднощів і поєднати раніше здавалися погано сумісні одна з одною властивості (наприклад, явность різницевих схем і її абсолютна стійкість, високий порядок апроксимації і монотонність, консервативність різницевих схем, побудованих з використанням характеристичних властивостей рівнянь і ін.). Це дозволяє сподіватися, що процес вдосконалення чисельних методів розв'язання рівнянь і систем гіперболічного типу і надалі буде продовжений. Поряд із зазначеними вище підходами є багато інших розробок, в яких гіперболічність систем рівнянь в явній формі не використовується, а присутні інші оригінальні підходи до чисельного вирішення таких рівнянь. Не маючи можливості зупинитися на них в цьому короткому огляді, наведемо в якості прикладів один з перших « бессеточная» методів - метод вільних точок [34], метод частинок в осередках [117], метод великих часток [141, 7], метод гладких частинок [179, 901, кінетичні узгоджені різницеві схеми 11351, бікомпактних схеми 1911 і ін. У першому виданні монографії також були представлені результати досліджень з використанням розроблених авторами сеточіо-харак- тсрістіческіх методів в області надзвукового обтікання літальних апаратів, механіки деформованого твердого тіла і фізики плазми (глави 5-7 відповідно). За минулі роки ці наукові напрямки також отримали розвиток в роботах авторів монографії, їх колег та учнів. Можливість приведення вихідних гіперболічних систем рівнянь до еквівалентним системам для умов спільності уздовж характеристик (в тому числі в багатовимірному випадку уздовж біхарактерістік або ліній перетину характеристичних і координатних поверхонь) дозволяє істотно спростити перехід від побудови різницевих схем для скалярних рівнянь до гіперболічних систем зі збереженням відповідних властивостей в скалярних аналогах багатовимірних різницевих схем. В області аерогазодинаміки літальних апаратів після першого видання досліджувалося надзвукове обтікання тіл, що деформуються | 75 |, надзвукове обтікання системи близько розташованих тел з урахуванням їх інтерференції [2], активно розвиваються зараз в роботах Ф.А. Максимова [67] та ін., Багатошарових течій стратифікованою по щільності нестисливої рідини в наближенні дрібної води 114,151 і ін. Деякі з результатів досліджень гінерзвукового обтікання затуплених тел сеточно-характеристичними методами увійшли в фундаментальну монографію (по суті енциклопедію гіперзвукових течій) [62]. У напрямку механіки деформівних твердих тіл слід зазначити роботи групи І.Б. Петрова з моделювання різноманітних динамічних задач з використанням сеточно-характеристичних методів, в тому числі високоточних версій па неструктурованих сітках, розроблених в [106, 104, 143, 114] та ін. Розглядалися 2D- і ЗБ-зіткнення тіл з багатошаровими деформуються перешкодами [78 , 80, 79] та ін., а також високошвидкісних зіткнень з урахуванням пружно-пластичних деформацій і різних моделей руйнування [40, 90] та ін .; впливу потужних потоків енергії на деформуються тіла 175] та ін .; проходження сейсмічних хвиль в різних, в тому числі тріщинуватих геологічних середовищах, а також їх і інших силових впливів на споруди і елементи складних конструкцій [47, 48, 139, 71, 140, 82, тисячі триста вісімдесят один і ін. Поряд з динамічними розглядалися стаціонарні навантаження різних конструкцій, в тому числі термопружних деформацій графена і аналогічних двовимірних матеріалів зі складною многосвязной областю інтегрування рівнянь Дюгамеля-Неймана [124]. У напрямку фізики плазми з використанням рівнянь магіі- тогазодінамікі розраховувалися: дво- і тривимірні задачі обтікання сонячним вітром магнітосфери Землі [130]; акустико-гравітаційні хвилі, що виникають в іоносфері від природних і антропогенних джерел збурень в приземному шарі атмосфери [55]. У розвиток роботи [30] по моделюванню експериментів типу «Морська зірка» досліджувалася еволюція в іоносфері на довгі часи високо енергетичних плазмових утворень і взаємодія декількох плазмових згустків [136, 116, 182, 181] та ін .; вхід і рух в атмосфері великих космічних тіл, включаючи взаємодію з поверхнею Землі і подальшу еволюцію генеруються плазмових утворень [125]; а також плазмових течій в дивертора токамаков [168]. Були отримані істотні результати в ряді інших наукових напрямків. Перш за все, це обчислювальні моделі деяких фізіологічних процесів в людському організмі. У офтальмохирургии проводилося моделювання іригаційне-аспіраційної техніки факоемульсифікації при екстракапсулярної екстракції катаракти [109, 129] і моделювання етапу руйнування кришталика в процесі ультразвукової та лазерної екстракції катаракти [12, 89]. Був вперше розроблений комплекс глобальних динамічних і квазістаціонарних моделей зовнішнього дихання і кровообігу з урахуванням їх зв'язності і перенесення речовин (в тому числі включає чотирикамерну модель серця) [119, 35] і проведено дослідження багатьох нормальних і патологічних процесів в дихально-кровоносній системі людини [151 , 177, 131, 107, 108, 155, 192, 159] та ін. Оскільки в медичних дослідженнях часто доводиться мати справу з великими обсягами сильно коливальних даних, перспективним напрямком є створення більш досконалої их, ніж традиційне середньоквадратичне наближення, методів апроксимації таблично заданих функцій. Як приклад можна привести підхід, заснований на доповнення умов середньоквадратичного наближення умовами безперервності на кордонах ділянок розбиття «довгих» таблиць і умовою мінімальності різниці найбільшого і найменшого значення функції на кожній ділянці розбиття [169]. Розроблено ряд мережевих обчислювальних моделей на основі нелінійних гіперболічних систем рівнянь на графах: інтенсивне вуличний рух автотранспорту в мегаполісах [73, 137, 85, 88] та ін .; руслових течії в великих річкових басейнах з використанням рівнянь мілкої води 174]; великі електроенергетичні мережі [72], інформаційні потоки в комп'ютерних і телекомунікаційних мережах [102] та ін., поширення домішок в великих вентиляційних мережах [86]. В рамках сеточно-характеристичного підходу також розроблений метод розв'язання обернених задач для нелінійних систем гіперболічного твань [1 |, що є обчислювальної версією управління рішеннями за допомогою граничних умов, активно досліджувалися раніше в рамках різних лінійних моделей [13, 93, 41, 43, 44, 42, 45] та ін. Нижче наведено кілька прикладів розрахунку з використанням сеточно-характеристичних методів ряду завдань, що моделюються нелінійними системами гіперболічного типу (в тому числі па графах), а також 20-параболічними і еліптичними рівняннями зі складними (в тому числі багатозв'язними) областями інтегрування на неструктурованих хаотичних сітках ( рис. Н.1-п.З). На рис. II. 1 представлені деякі результати розрахунку з використанням високоточної консервативної версії сеточно-характеристичного методу в рамках МГД рівнянь одним із завдань еволюції збурень в атмосфері і в іоносферній плазмі, викликаних рухом Челябінського метеороида по пологій траєкторії з імітацією його миттєвого руйнування [125]. Зліва показана спостерігалася в реальності світиться частина обуреної Челябинским болідом атмосфери з [4] (фото А. Ахметвалеева). На правих панелях показані розраховані розподілу швидкості (верхня панель), питомої внутрішньої енергії (середня панель) і магнітного поля (нижня панель) в момент взаємодії з поверхнею Землі ударної хвилі, що відокремилася від швидко гальмуються залишків метеороида. Форма високотемпературної області (на середньої панелі) цілком порівнянна зі світиться областю на лівій частині малюнка. ![]() Мал. 11.1 На рис. П.2 наведені результати чисельного моделювання переносу кисню в артеріальній системі людини з використанням 20-розподілу моделі кровообігу і конвективно-дифузійного переносу речовин кров'ю. Параметри несучої фази (нуля тисків і нуля швидкостей) визначалися на основі рівнянь потенційного течії нестисливої рідини в великих судинах і фільтрації в «пористої середовищі» (органи і тканини). Незважаючи на досить грубу двовимірну модель, характер перебігу відтворюється цілком адекватно реальності. Зокрема, час перенесення кисню від лівої легені (праве належало віддаленим) до мочки лівого вуха в розрахунках склало 4-4.5 з (при експериментальних його значеннях - 3.5 с) [35]. ![]() Мал. 11.2. Приклад розрахунку динаміки концентрації кисню в організмі людини (3 подиху) На рис. П.З показані результати розрахунків сеточно-характеристичним методом з використанням розробленої в [73, 137] мережевий обчислювальної моделі добового автомобільного руху на одній з ділянок автотраси 1-80 в районі затоки Сан-Франциско (рис. П.З а). На рис. П.З 6 для однієї з ділянок автостради наведені регулярно оновлювані з датчиків дані потік-щільність в різні моменти часу доби, що передують моделювати діб (хрестики) і закладає в якості останнього розрахункову гідравлічну модель рівняння стану їх апроксимація (регулярно оновлювана фундаментальна діаграма крива). На рис. П.З в вгорі для цієї ж ділянки приведені в залежності від часу виміряні значення автомобільного потоку в моделюються добу (хрестики) і дані розрахунку (крива). ![]() Мал. з За минулі роки пішли з життя мої вчителі - видатні російські вчені: багаторічний ректор МФТІ і директор Інституту автоматизації проектування РАН академік О.М. Білоцерківський (керівник обох моїх дисертацій, відповідальний редактор першого видання монографії і співавтор багатьох спільних робіт), співавтор багатьох робіт але сеточно-характеристичним методам і першого (наукового) видання монографії, творець першого в Дагестані політехнічного інституту, керівник Дагестанського наукового центру РАН, директор Інституту геотермії РАН К.М. Магомедов (керівник моєї дипломної роботи в лабораторії ЦНДІ машинобудування, на щастя, досі живий професор В.В. Луньов - мій науковий «дідусь»). Всім, з ким пощастило працювати довгі роки, низький уклін і подяка. Автор вдячний також редакції видавництва і Поліні Сеі- чило за допомогу в оформленні книги. А.С. Холодов, м Волоколамськ, серпень 2016 р ЛІТЕРАТУРА
В.А. Харламов, І.А. Трубецкая // Астрономічний вісник. 2013. Т. 47. №4. С. 116.
С. 144-155.
вирахував. математики і мат. фізики. 2011. Т. 51. №2. С. 282-302.
A. В. Гасніков, І.І. Морозов, B. Н. Тарасов // Праці МФТІ. 2010. Т. 2. №4. С. 152-162. 74.0 чисельному моделюванні деяких завдань взаємодії літосфери, гідросфери та атмосфери Землі / О.М. Білоцерківський, М.О. Васильєв, А.Б. Ведерников, В.П. Димні- ков, Б.В. Замишляєв, Б.Ю. Кри- санів, II.В. Ковшов, В.Є. Куніцин, Е.А. Молоков, АЛО. Ренін, Н.А. Сидоренкова, Е.Л. Стуніц- кий, Я.А. Холодов, А.С. Холодов // Фрагменти історії та досягнення ІАП РАН. 1986-2011. М .: Політ Джонатана, 2011. С. 14-71. 75. Про чисельному рішенні пов'язаних завдань надзвукового обтікання деформованих оболонок кінцевої товщини / П.М. Коро- тин, І.Б. Петров, В.Б. Пирогов, О.А. Пиркова, А.С. Холодов // Журі, вирахував. математики і мат. фізики. 1988. Т. 27. № 8. С. 1233- +1243. 76. Операційні різницеві схеми / А.А. Самарський, В.Ф. Тішкін, А.П. Фаворський, М.Ю. Шашков // Д іфферен ци ал ь н и е у равн е і і я. 1981. Т. 17. №7. С. 1317-1321.
С. 58-72.
A. С. Холодов, І.В. Цибулін // Мат. моделювання. 2016. Т. 28. № 8. С. 65-81. 87. Різницеві схеми на нерегулярних сітках / А.А. Самарський, Ю.А. Колдоба, В.Ф. Повещенко, B. Ф. Тішкін, А.П. Фаворський. Мінськ: Критерій, 1996.. 88. Розробка, калібрування і верифікація моделі руху трафіку в міських умовах. Частина I. Комп'ютерні дослідження і моделювання / А.Є. Алексєєнко, Я.А. Холодов, А.С. Холодов, А.І. Горєва, М.О. Васильєв, Ю.В. Чехович, В.Д. Мішин, В.М. Старожілец // Комп'ютерні дослідження і моделювання. 2015. Т. 7. № 6. С. 1185-1203.
В.І. Кондауров, І.Б. Петров, О.С. Холодов // Мат. моделювання. 1990. Т. 2. № І. С. 10-29.
С. 42-52. 107. Симаков С.С., Євдокимов А.В., Холодов А.С. Методи розрахунку глобального кровотоку в організмі людини з використання гетеро ге інших ви ч і з л і тел ьн и х моделей // Медицина в дзеркалі інформатики / За ред. О.М. Білоцерківського, А.С. Холодова. М .: Наука. 2008. С. 76-110.
ПО. Тихонов А.Н., Самарський АЛ. Про розривних рішеннях квазілінійних рівнянь першого порядку // Докл. АН СРСР. 1959. Т. 124. № 3. С. 529-532.
С. 104-113. 119. Холодов А.С. Деякі динамічні моделі зовнішнього дихання і кровообігу з урахуванням їх зв'язності і перенесення речовин // Комп'ютерні моделі і прогрес медицини / За ред. О.М. Білоцерківського, А.С. Холодова. М .: Наука, 2001.. С. 127-163.
С. 1601-1620.
С. 3-23.
A. В. Фаворська, І.Б. Петров, B. І. Голубєв, II.І. Хохлов // Мат. моделювання. 2015. Т. 27. № 12. С. 109-120. 140. Чисельне моделювання сеточно-характеристичним методом впливу землетрусів на спорудження / A. В. Фаворська, І.Б. Петров, B. І. Голубєв, II.І. Хохлов // Мат. моделювання. 2015. Т. 27. № 12. С. 109-120.
C. К. Годунов, А.В. Забродін, М.Я. Іванов, О.М. Край до і ін. М .: Наука, 1976.
Adaptive Grid and Grid-free Methods for the Edge Plasma Fluid Simulations // Journal Plasma Physics. 1999. Vol. 61. №5. P.701-722.
Cratering // 4 th Aerodynamic Testing Conference. 1969. P. 69-354.
SS Siniakov // Russian Journal of Numerical Analysis and Mathematical Modelling. 2015. Vol. 30. № 5. P. 269-276.
|
<< | ЗМІСТ | >> |
---|