Anderson-et-al-1 (1185923), страница 18

Файл №1185923 Anderson-et-al-1 (Андерсон Д., Таннехилл Дж., Плетчер Р. Вычислительная гидромеханика и теплообмен) 18 страницаAnderson-et-al-1 (1185923) страница 182020-08-25СтудИзба
Просмтор этого файла доступен только зарегистрированным пользователям. Но у нас супер быстрая регистрация: достаточно только электронной почты!

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

е. прн 4г„яп' й/2+ 4г„яп' й/2 ~ (2. Последнее условие выполняется лишь при г, + г„(1/2. Возвращаясь к шагам разностной сетки, получим условие устойчивости в виде сй! [1/(Лх)з + 1/(Лу)з[ (~ 1/2. Это условие аналогично полученному в одномерном случае, хотя эффективный шаг по времени в двумерном случае оказывается меньше, чем в одномерном. В рассмотренном примере результат удалось получить довольно легко, но в большинстве случаев провести анализ устойчивости конечно-разностной схемы решения уравнения в частных производных очень трудно для более чем одной пространственной переменной и времени. Условие устойчивости удается часто определить, лишь вычислив коэффициент перехода для различных значений г„ и г„.

3.6.2. Анализ устойчивости систем уравнений в частных производных Описанный в предыдущем разделе метод Неймана может быть' использован и для анализа устойчивости конечно-разностных схем решения систем уравнений в частных производных. Встречающиеся в газовой динамике и теплопередаче системы уравнений можно обычно записать в виде дн дР— + — =-О, дт дх= где Е и à — векторы, причем Г = Г(Е). Обычно система уравнсний являстся нелинейной, а методы анализа устойчивости й 3,6, устойчивость конечно-рааностных схем ных первого порядка — +с — =0 дн до да дх — + с — =О.

до ди др дх В рассматриваемом случае -П Перепишем систему в матричном виде дс +[А] д. — — О, дЕ мдЕ где [А]= '0 Следовательно, максимальное собственное значение матрицы [А] равно с и конечно-разностная схема устойчива при выполнении обычного условия устойчивости КФЛ ~ с — „~(1. Отметим, что описанный выше метод анализа устойчивости конечно-разностных схем не учитывает влияния граничных условий, хотя и использует матричную запись системы уравнений в частных производных. Как учесть влияние граничных условий, будет показано ниже. Соотношение (3616) показывает, что устойчивость конечноразностной схемы определяется видом матрицы перехода. Перепишем (3.116) в виде е"+'(Й) = [6(Ы, л]" [е'(й)].

(3.120) Теперь можно сформулировать условие устойчивости разностной схемы [К!сЫшуег, Мог1оп, 1967]. Оно сводится к требованию существования такого положительного т, чтобы матрица [6(схг,й)]" была равномерно ограничена при 0 < схт < т, 0 < < иИ < Т для любых й, где Т вЂ” максимальное время. Отсюда следует необходимое условие устойчивости Неймана: [о,(Ы, й)[(~-1+0(Ж) при 0<ар< с (3.121) для всех собственных значений и волновых чисел, где гн — собственные значения матрицы [6(И,Й)]. В рассмотренных примерах мы считали конечно-разностную схему устойчивой, если максимальное собственное значение 4 д. Андерсон н др. том ! Гл. 3.

Основы метода нонечных разностей матрицы перехода по модулю непревосходитединицы. Этоусловие более сильное, чем (3.118). Необходимое условие устойчивости Неймана накладывает требование приемлемого роста локальной величины сЫ, что в действительности возможно для решения многих физических задач. Классическим примером, иллюстрирующим это, является разностная схема решения уравнения теплопроводности с источниковым членом.

Пример 3.6. Пусть мы хотим решить уравнение теплопроводности с источниковым членом ди дзи — =а — +си д~ дкз используя простую явную конечно-разностную схему. Применив анализ Фурье, найдем коэффициент перехода б = 1 — 4г з(пз,((1/2) + сдам. Отсюда следует, что решение данного разностного уравнения может расти по времени и удовлетворяет необходимому условию устойчивости Неймана.

Мы видим, что прн анализе устойчивости разностных схем полезно привлекать физические соображения. Следует помнить, что для системы гиперболических уравнений строгое условие заключается в том, что максимальное собственное значение не должно превосходить по модулю единицу. Это связано с тем, что гиперболические уравнения аналогичны волновому и не имеют экспоненциально растущих по времени решений.

Мы исследовали устойчивость различных конечно-разностных схем, используя метод Неймана. Если мы хотим изучить влияние граничных условий на устойчивость конечно-разностной схемы, то мы должны использовать матричный метод. Это проще всего показать на примере применения метода Лакса к решению одномерного линейного волнового уравнения первого порядка ди1 - 'ди — "'"' — +„с — = О. дЗ " дх Пусть сетка по х состоит из яз точек, а граничные условия периодические, т. е. и» и» (3.122) Применяя для решения этой задачи метод Лакса, получим систему алгебраических уравнений вида п"+' = [Х] и", (3.123) э 3.6. Устойчивость конечно-разностных схем где а т и =]иь и„..., и,„], (3.124) 1 та 2 1 — т Π— О 2 1+а 1 — а О 2 2 О 1+ 2 (3225) (х) = Устойчивость конечно-разностной схемы (3.123) зависит от собственных значений матрицы ]Х].

Так как заданы периодические граничные условия, то в матрице 1Х] существенны лишь элементы на трех отмеченных в соотношении (3.123) диагоналях и два угловых элемента. Такая матрица называется апериодической ]1лтах е1 а1., 1970]. Для матрицы вида (34 26) ат аа О аа а, собственные значения имеют внд 2я 2м йч — — а, + (по+ а,) соз — (] — 1) + 1(ав — ' ав) з1п — (1 — 1), 1=1, лт. В рассматриваемом случае 1+в аь аи О ' ' ао аа а1 аг О О ао 1-а 2 1эа — О 2 Гл. 3. Основы метода конечных разностей и собственные значения имеют вид 2н 2н Л! — — соз — (1 — 1) + Рр з(п — (1 — 1). нз нз (3. 127) !+т ! — т — О 2 2 Π— О !+т 2 ! — ч 2 О О 1 О Собственные значения этой матрицы легко вычисляются.

Они равны А! — — 1, Аз — — О, Аз,ч=:Ь 2 Ч(1 — т)(З+ т). Условие )).((1 приводит к тем же ограничениям на т, что и условие КФЛ: ! т ! ( 1. В рассмотренном примере граничные условия не изменили условия КФЛ устойчивости схемы. Однако граничные условия обычно изменяют условия устойчивости и этого следует ожидать. Итак, конечно-разностная схема устойчива при !т(( 1, т. е. при выполнении условия КФЛ. Это показывает, что матричный анализ устойчивости схемы Лакса привел к тем же результатам, что и метод Неймана для простого волнового уравнения. При периодических граничных условиях матричный метод и метод Неймана приводят по существу к тождественным результатам.

Проиллюстрируем влияние граничных условий и разностной сетки на другом примере. Пусть, как и в предыдущем примере, метод Лакса применяется к решению одномерного линейного волнового уравнения первого порядка. Ограничимся случаем, когда для решения уравнения используются лишь четыре точки по пространственной координате. В первой и четвертой точках зададим граничные условия. Предположим для простоты, что в первой точке величина и постоянна для всех времен, тогда граничное условие можно записать в виде и",+' = и",. Так как мы ищем решение волнового уравнения, то граничное условие в четвертой точке на правом конце не может быть произвольным, а должно быть согласовано с направлением распространения волны.

Зададим его в виде и",+!ч и", т. е. будем определять значение функции и на границе по ее значению во внутреннем узле сетки на предыдущем слое. При таких граничных условиях матрица [Х] имеет вид 1 О О О )0( Задачи Ясно, что матричный метод анализа устойчивости разностных схем учитывает заданные на сетке граничные условия, Это означает, что при матричном методе анализа влияние граничных условий автоматически учитывается. К сожалению, определить аналитически собственные значения матриц при произвольно заданных граничных условиях обычно не удается. В этом разделе проведен анализ устойчивости конечно-разностных схем при помощи метода Фурье (метода Неймана) и матричного метода.

Эти два метода широко применяются для анализа устойчивости разностных схем. Однако разработаны и часто используются и другие методы. В этой связи укажем на работы (Н!г1, 1968] и 1%агш)пд, Нуе(1, 19741. Наиболее всесторонний математический анализ устойчивости конечно-разностных схем с доказательством многих теорем изложен в монографии Рихтмайера и Мортона [Рс!сЫгпуег, Мог1оп, 1967). Задачи ЗЛ. Покажите, что дзи ~ гт,иг > з — = — + О (йх). дхз ( (Ьх)з 3.2. Рассмотрите функцию 1(х) = е". Определите на сетке с шагом Ьх = Од производную 1'(х) при х = 2, используя разности вперед (3.26), пентральные разности (3.28) и трехточечную аппроксимацию второго порядка (3.29).

Сравните полученные результаты с точным значением производной. Повторите вычисления при Лх = 0.2. Правильно ли описывает точность аппроксимации производной порядок погрешности аппроксимации? Проанализируйте полученный результат.

3.3. Проверьте, является ли консервативной конечно-разностная схема для уравнения неразрывности двумерного стационарного движения несжимаемой жидкости (и +и — иг — и ) (ог à — о г Г >) 2Ьх ад О, где и, о — составляющие скорости по асям х и у соответственно, 3.4. Задание то же, что и в задаче 3.3, по для разностной схемы (игь> г "; >!) (ог >,г ог! >) 2Ьх — О. 2ау 3.3. Рассмотрите нелинейное уравнение ди дзи и — =и— дх дра с постоянным козффнцпентом р. (а) Записано ли оно в ливергентпой форме? Если нет, то ма>кете лч вы записать его в днвергсптиой форме? 102 Гл. 3.

Основы метода конечных разностей (Ь) Постройте конечно-разностный аналог этого уравнения, используя интегральный метод. 3.6. Проверьте конечна-разностиую аппроксимацию (3.50) для производной дзи/дхду. 3.7. Проверьте конечно-разностную аппроксимацню (3,40) для производной д'и/дхз. 3.8. Проверьте соотношение (3.79) в табл.

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

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

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