Главная » Просмотр файлов » Н.Н. Калиткин - Численные методы

Н.Н. Калиткин - Численные методы (1133437), страница 92

Файл №1133437 Н.Н. Калиткин - Численные методы (Н.Н. Калиткин - Численные методы) 92 страницаН.Н. Калиткин - Численные методы (1133437) страница 922019-05-12СтудИзба
Просмтор этого файла доступен только зарегистрированным пользователям. Но у нас супер быстрая регистрация: достаточно только электронной почты!

Текст из файла (страница 92)

Поэтому многие вычислители предпочитают использовать безусловно устойчивые неявные схемы. Построим хорошую неявную схему для задачи (1). Возьмем изображенный на рис. 92 шаблон и составим схему с весами при 428 ГИПЕРБОЛИЧЕСКИЕ УРАВНЕНИЯ !ГЛ. ХН1 пространственных производных на разных слоях: —, (у — 2у+ у) = Л (оу + (1 — 2о) у+ оу ~+ ), 1 (12) Лул= !1» Ьл"1 2ул+Ул-1)! чтобы все веса были неотрицательны, следует брать О - 1г«1/». В граничных узлах решение определяется из краевых условий (26). Исследуем построенную схему.

Значения решения на нулевом и первом слоях вычисляют, как и в п, 1, по формулам (3) и (4б). На остальных слоях схема (12) с краевыми условиями (2б) образует относительно у линейную систему уравнений с трехдиагональной матрицей, в которой диагональные элементы преобладают; решение этой системы существует, единственно и вычисляется методом прогонки. Разложением решения по формуле ТейлоРяс 82 ра нетрудно установить, что на решениях с непрерывными четвертыми производными схема (12) аппроксимирует уравнение (!а) с погрешностью 0(т»+Ь») при любом о. Устойчивость проверяется методом разделения переменных.

Делая подстановку (7), получим для множителя роста квадратное уравнение Рд 2сгдуд+ 1 = О> 1 — 2 (1 — 2д) Рд дт . 41! (13) На основании тех же рассуждений, что и в п. 1, можно сделать следующий вывод: устойчивость будет только при комплексно сопряженных корнях, т. е. при !сд ~ «1. Отсюда вытекает условие устойчивости схемы: ( — ~ (! — 4о) 1. (14) Из неравенства (14) видно, что при о)'/л схема (12) безусловно устойчива.

Если а(1)„то схема условно устойчива при ст=- «й (1 — 4о)- ч . Таким образом, при выборе веса 1(, -.= о«11» неявная схема (12) безусловно сходится с точностью 0(т»+Ь»). Замечание !. Схема (12) при о=-О переходит в схему «крест», а условие, устойчивости (14) — в условие Куранта (9).

429 ВОЛНОВОЕ УРАВНЕНИЕ Замечание 2. Обобщим неявную схему (12) на случай задачи с переменной скоростью звука: —,= —,)-~й(х, 1) — — ")+[(х, (), й(х, 1)=с'(х, 1))0, (15) где коэффициенты й(х, 1), 1 (х, 1) переменны и кусочно-непрерывны вместе со своими вторыми производнымн, причем разрывы неподвижны (т. е. лежат на линиях х=сопз1).

Предполагается, что на этих разрывах выполняются условия сопряжения [и) = О, [яи,) = О. Выберем по 1 равномерную сетку, а по х — специальную неравномерную сетку (у которой все точки разрыва коэффициентов являются узлами). Построим аналог наилучшей консервативной схемы (11.34), используя во всех пространственных операторах значения я(х, 1) со среднего слоя: —, (у — 2у+у) =Л [ау+ (1 — 2о) у+ну)+<р, (1) 1 [ (1 ) Улп! (О Ул (б и (1 )Уп (О Уп-! (!)~ л л 1 1= (т-» )т (тм Лп = хлк! хл~ йп = (йп+ йл — !)~ 2 и.!.! ! 1 1' !1к !т ! кл+Ч, 1 <рл= — !(1 ~ ) (~, 1) !(~. т л ! !."! 'л и (16) к к о(х, 1)=)и!($, 1)сЦ, г(х, 1)=)1($, 1)Щ, (17) о а Известно, (см. [30)), что при сделанных предположениях (и достаточно гладких начальных и граничных данных) эта схема равномерно сходится со скоростью 0(т'+шахйк), если выполясно условие устойчивости (14).

Из схемы (16) нетрудно получить схемы для гладких и для постоянных коэффициентов на произвольных неравномерных по х сетках. В случае )1 = сопз1 и равномерной сетки схема (16) совпадает со схемой (12). 3. Двуслойная акустическая схема. Уравнение второго порядка (1а) можно заменить эквивалентной ему парой уравнений первого порядка. Для этого введем потенциалы скоростей и правой части: 430 ГИПЕРБОЛИЧЕСКИЕ УРАВНЕНИЯ 1ГЛ Х111 Функции и(х, 1), о(х, 1) удовлетворяют системе уравнений акустики и1 = Ол, О1 = Сги, + Р (Х, 1).

(18а) Начальные условия (!б) с учетом (17) примут вид к и (х, 0) = р,д (х), о (х, 0) = ) р,(с) с(а, 0 (18б) а граничные условия (1в) останутся без изменения: и (О, Г) = р, (1), и (а, 1) = р, (1). (18в) нередко оказывается более удобной для чем волновое уравнение (1); в частности, она позволяет построить двуслойные разностные схемы, допускающие неравномерную сетку по Г. Неявная схема. Будем рассматривать в узлах неравномерной пространственной сетки величины у, = —.и(хл, 1), О=п=аХ, а в серединах интервалов — ве- 1), О~п~Л1 — 1.

Возьмем шаблон, изоби составим на нем схему с весами Задача акустики (18) численного решения, л л УП ЕП апик хп гп л 1/К гп Рис. 93 личины гл о(хл+ и раженный иа рис. 93, г.)= — „-[о,(у,.„— у,)+(1 — о,) (у.,— у.)1+1р, л Ул) [ол(г л г л-1) + (1 02) (гл гл .1Н~ л 1 Ь Х +1 Х й л (6 +й -1) (19а) (19б) у =яе'л л 1Лк гл=рЕ, ул=рул, г,=ргл, (20) подразумевается, что 0-=. о; == 1, 0-.=.

д, ==. 1. Исследуем эту схему, подробно останавливаясь только на наиболее важных деталях и для простоты ограничиваясь равномерной по х сеткой. Схема (19) составлена симметрично по переменной х, если положить 1рл =- г (х, ч л, 1); поэтому нетрудно сообразить, что она имеет аппроксимацию 0(т+Ь'). Если взять о,=о,=", и 1р„=г (х„+ и, 1„+т/2), то схема становится вдобавок симметрйчной по времени и приобретает аппроксимацию 0(т'+/1').

У с т о й ч и в о с т ь исследуем методом разделения переменных, рассматривая возмущения функций у и г в виде гармоник ВОЛНОВОЕ УРАВНЕНИЕ 43! с одной и той же частотой и множителем роста, но с разными амйлитудами я и р. Подставляя (20) в (19) и полагая»р,=О, получим для амплитуд систему линейных однородных уравнений -- (р — 1) а — —. (а,р + 1 — а,) (1 — е-'» ) р = О, ! ! и (21) — „(а,р + 1 — о;) (е!»А — 1) а — — (р — 1) !) = О.

(25) и становится явной. В самом деле, величины 2, явно выражаются из уравнения (2ба) через значения величин на исходном слое. Чтобы она имела нетривиальное решение, ее определитель должен обращаться в нуль. Это дает квадратное уравнение для нахождения множителей роста р: ер' — 2!Ар+ у = 0; (22а) е = 1 + 2у»а,а, ) 1, р = 1 — у» (а, + ໠— 2а1а»), т = 1+ 2у» (1 — а,) (1 — ал) ~ 1, у» — — 2 ~ — з1п -2 — ) ) О. ! лт .

»А !» (22б) Оба корня уравнения (22а) меньше единицы по модулю тогда и только тогда, если т~е, 2(р)(е+т. (23) Первое из этих неравенств очевидно, поскольку по теореме Виета р„'р» — — у/е; второе доказывается несложными, но громоздкими выкладками. Неравенства (23) выполняются для всех гармоник только в том случае, если а,+а,~1, ( — „— ) (2ат — 1)(2а,— 1)~ — 1, (24) что является условием устойчивости схемы (19).

При выполнении этого условия схема сходится со скоростью, соответствующей порядку аппроксимации. Из неравенств (24) вытекает, что если а, ~ '/, и а, — '/, то схема (19) безусловно устойчива. Если а, + а,'=- 1, но один из весов меньше '/,, то схема условно устойчива при ст Ь/)' (2а, — 1) (1 — 2а,). Если а,+ае«„-1, то схема безусловно неустойчива.

Рассмотрим два частных случая схемы (19). Явная схема. Положим а,=О, а,=1; тогда схема (19) принимает вид т (Ел Ел) = А (!/л+1 ул)+Гл+1у,', 0~~!! -/!/ — 1, (2ба) — (р,— р)= —,(г,— Ел т), 1(п~/!/ — 1, (266) 1гл, х»и ГИПЕРБОЛИЧЕСКИЕ УРАВНЕНИЯ После того, как вычислены все значения гл, можно найти ул также по явным формулам (266). Из (26) следует, что схема (26) устойчива при выполнении условия Куранта ст ( й. Заметим, что схема (26) является схемой типа «крест».

В самом деле, будем считать, что величина г„=с(хл+ и, 1 ) сдвинута на полшага по х, а величина у„=и(х„1„4 и) — на полшага по 1 относительно узлов сетки (рис'. 94). Тогда этой схеме соответствует шаблон из двух крестов, показанных на рисунке жирными линиями. Зададим согласованные с этим шаблоном граничные данные: р„=)»з((~еч,) р~~=)»«((меч,) (26в) и начальные данные, уточненные аналогично (4б), где надо вместо т взять т(2: Уй»=Р»(Х»)+ й )Х»(Хл)+ В ~С,Ц~» +1(Х» )1 кл-~-Ч~ (26г) е."= ~ р (~) $.

к» Тогда схема (26) при выполнении условия Куранта ст(й сходится со скоростью О (т»+Ь»). Симметрич н а я схема. Положим о,=а»='7»; тогда схема (19) является безусловно устойчивой и сходится со скоростью О (т»+ Ь»). Эта схема двуслойна, поэтому она позволяет произвольно менять шаг т в ходе расчета, обеспечивая при этом точность 0(тахт,"„). Кроме того, поскольку значения у«и г' соответствуют моменту 1= О, начальные данные для расчета берут непосредственно из постановки задачи (18): Й» "л»- Ч, дл = )х» (хл), ел = ~ р» (5) с(«, (27) Рис.

94. без каких-либо сдвигов по времени; такая аппроксимация начальных условий является точной. Однако при любых значениях весов п„о» (если только один из них не равен нулю) схема (19) неявна. Рассмотрим, как целесообразно вычислять разностное решение в этом случае. Определим г„из уравнения (19а) и напишем аналогичное выражение для г„». Подставляя эти вь|ражения в (19б) и полагая ВОЛНОВОЕ УРАВНЕНИЕ для простоты Ь = сопз(, получим —, (у — ул)= „ь ол(о~(ул+1 — 2у +у -~)+ 1 + (1 — О,) (Ул„— 2Ул + Ул,)1+ — (гл — гл,) + т + и Оа (Ч'л Ч~л-~) Это линейная система относительно неизвестных у, имеющая трехдиагональную матрицу с преобладанием диагональных элементов; ее решение легко вычисляется прогонкой.

Найдя у, нетрудно определить е по явным формулам (19а). Таким образом, симметричная схема (19) приводит к несложному' вычислительному алгоритму, безусловно устойчива и имеет хорошую точность. Она является одной из лучших схем для расчета задач акустики. По аналогии с ней строятся надежные однородные схемы расчета газодинамических и других сложных задач. 3 а м е ч а н и е 1. Разиостное решение схемы (19) можно, вообще говоря, вычислять методом последовательных приближенйй: ев! = .+ — "1о (уп"+ — уо')+(! — ) (у.

— у.)1+ р., (29) ул =уп+ ~ !Ол(ел ел — 1)+(1 ол) (ел ел-1)] Однако это эквивалентно применению последовательных прибли- жений к решению системы (28), когда в левой ее части берется у"'", а в правой — у"'. Этот метод, записанный в форме ум'» = = Су~ю+ф сходится при !~С ~~ ~1. Выбирая одну из норм: ~~с~~= .*(~~С„,) =-пф,п„ л получим для сходимости итераций условие типа Куранта: ст~ Ь 2 !' В1лл (30) Поэтому метод последовательных приближений невыгодно применять к вычислению разностного решения безусловно устойчивых схем. Замечание 2. Для задач с разрывными или недостаточно гладкими решениями нередко используют чисто неявную схему (19) при а = а, = 1, поскольку она подавляет лразболтку» счета. Однако на достаточно гладких решениях эта схема существенно уступает по точности симметричной схеме.

(гл. хги ГИПЕРБОЛИЧЕСКИЕ УРАВНЕНИЯ 4. Инварианты. Рассмотрим запись системы уравнений акустики через инвариантеп (31) г (х, () = о — си, з (х, () = о+ си. Умножая первое из уравнений (18а) на с, прибавим его ко второму уравнению (18а) и вычтем. Получим систему уравнений, которым удовлетворяют инварианты: г,+сг„=Р(х, У), зт — сз,=Р(х, )). (32а) Из соотношений (18б) нетрудно получать для инвариантов началь- ные условия: к э г(х, О) =~рз($) с$ — ср (х), з(х, О) =~р Я) ГК+ср,(х), (32б) о о а из соотношений (18в) — краевые условия: (з г)э-е= 2 рэ (~) ( )э-а= 2 ра (). с С (32в) — (г„— г„) + л (г.а — г„,) = Є— (з„— з„) — -,— (з„,, — з„) = Р„. (33) 1 с 1 с Она действительно является схемой бегущего счета, и организация вычислений здесь почти такая же, как для одномерного уравнения переноса, Нетрудно показать, что при выполнении условия ст ==6 эта схема устойчива, монотонна и равномерно сходится с порядком тОЧНОСтИ О (т+ гт) На ДВаждЫ НЕПРЕРЫВНО ДнффЕРЕНЦИРУЕМЫХ решениях.

Характеристики

Тип файла
PDF-файл
Размер
21,94 Mb
Тип материала
Высшее учебное заведение

Список файлов книги

Свежие статьи
Популярно сейчас
Почему делать на заказ в разы дороже, чем купить готовую учебную работу на СтудИзбе? Наши учебные работы продаются каждый год, тогда как большинство заказов выполняются с нуля. Найдите подходящий учебный материал на СтудИзбе!
Ответы на популярные вопросы
Да! Наши авторы собирают и выкладывают те работы, которые сдаются в Вашем учебном заведении ежегодно и уже проверены преподавателями.
Да! У нас любой человек может выложить любую учебную работу и зарабатывать на её продажах! Но каждый учебный материал публикуется только после тщательной проверки администрацией.
Вернём деньги! А если быть более точными, то автору даётся немного времени на исправление, а если не исправит или выйдет время, то вернём деньги в полном объёме!
Да! На равне с готовыми студенческими работами у нас продаются услуги. Цены на услуги видны сразу, то есть Вам нужно только указать параметры и сразу можно оплачивать.
Отзывы студентов
Ставлю 10/10
Все нравится, очень удобный сайт, помогает в учебе. Кроме этого, можно заработать самому, выставляя готовые учебные материалы на продажу здесь. Рейтинги и отзывы на преподавателей очень помогают сориентироваться в начале нового семестра. Спасибо за такую функцию. Ставлю максимальную оценку.
Лучшая платформа для успешной сдачи сессии
Познакомился со СтудИзбой благодаря своему другу, очень нравится интерфейс, количество доступных файлов, цена, в общем, все прекрасно. Даже сам продаю какие-то свои работы.
Студизба ван лав ❤
Очень офигенный сайт для студентов. Много полезных учебных материалов. Пользуюсь студизбой с октября 2021 года. Серьёзных нареканий нет. Хотелось бы, что бы ввели подписочную модель и сделали материалы дешевле 300 рублей в рамках подписки бесплатными.
Отличный сайт
Лично меня всё устраивает - и покупка, и продажа; и цены, и возможность предпросмотра куска файла, и обилие бесплатных файлов (в подборках по авторам, читай, ВУЗам и факультетам). Есть определённые баги, но всё решаемо, да и администраторы реагируют в течение суток.
Маленький отзыв о большом помощнике!
Студизба спасает в те моменты, когда сроки горят, а работ накопилось достаточно. Довольно удобный сайт с простой навигацией и огромным количеством материалов.
Студ. Изба как крупнейший сборник работ для студентов
Тут дофига бывает всего полезного. Печально, что бывают предметы по которым даже одного бесплатного решения нет, но это скорее вопрос к студентам. В остальном всё здорово.
Спасательный островок
Если уже не успеваешь разобраться или застрял на каком-то задание поможет тебе быстро и недорого решить твою проблему.
Всё и так отлично
Всё очень удобно. Особенно круто, что есть система бонусов и можно выводить остатки денег. Очень много качественных бесплатных файлов.
Отзыв о системе "Студизба"
Отличная платформа для распространения работ, востребованных студентами. Хорошо налаженная и качественная работа сайта, огромная база заданий и аудитория.
Отличный помощник
Отличный сайт с кучей полезных файлов, позволяющий найти много методичек / учебников / отзывов о вузах и преподователях.
Отлично помогает студентам в любой момент для решения трудных и незамедлительных задач
Хотелось бы больше конкретной информации о преподавателях. А так в принципе хороший сайт, всегда им пользуюсь и ни разу не было желания прекратить. Хороший сайт для помощи студентам, удобный и приятный интерфейс. Из недостатков можно выделить только отсутствия небольшого количества файлов.
Спасибо за шикарный сайт
Великолепный сайт на котором студент за не большие деньги может найти помощь с дз, проектами курсовыми, лабораторными, а также узнать отзывы на преподавателей и бесплатно скачать пособия.
Популярные преподаватели
Добавляйте материалы
и зарабатывайте!
Продажи идут автоматически
6417
Авторов
на СтудИзбе
307
Средний доход
с одного платного файла
Обучение Подробнее