Anderson-et-al-1 (1185923), страница 16
Текст из файла (страница 16)
3. Основы метода конечных разностей либо по формуле ад ив+ ан „ран ("и "в). Можно использовать и среднее арифметическое этих двух зна- чений: ир =(ир + ир )/2. 3. Построим конечио-разностную схему, пригодную для решения уравнения в частных производных на нерегулярной сетке. Последнее проще всего сделать, используя интегральный метод на неравномерной сетке. Проинтегрируем дифференциальное уравнение в окрестности точки (хо,уо), выбрав верхние и нижние ! )+1 пределы интегрирования отстоящими на полшага сетки от этой точки. Тогда в соответствии с приведенными на рис. 3.9 обо!-!,/ ;+! значениями получим Ух+да+/2 «о+ах+/2 уо Ду-/2 «о Дх /2 хо + ах+/2 Уо+ да+/2 «о ох-/2 Уо-да /2 Рис. 3«к Обозначения, используемые ири решении уравнений на нерегулярной сетке. (3.97) Внутренние интегралы могут быть вычислены аналитически: но+ ЛУ+/2 ~ — (хо + й о У) (хо й У)( с(У + Уо-ДУ /2 хо+ах+/2 + ~ ~ — (х, уо+ — ) — — (х, уо — — )~с/х=О.
хо-Ьх /2 Остальные интегралы можно вычислить, используя теорему о среднем значении: (3.98) з 3.5. Применение нерегулярных сеток 85 Если для аппроксимации производных в последнем соотношении воспользоваться центральными разностями, как это было сделано в п. 3.4.3, то получим следуюгций конечно-разностный аналог уравнения Лапласа: (и.. — и.
и. — и. ( ее!! и! ,1 е-!,1) йх +ах х Ьх Ьх Используя это выражение для точек, прилежащих к нерегулярной границе (рис. 3.8), получим следующие конечно-разностные аналоги вторых п)роизводных: дти ( 2 (ис ир ир ил '! дхе !н Ьх !! + и) т, мох Ьх (3. 100) д'и / 2 (ив ив ин ив1 дуе !в Ьу (! + 5) ~ аду Ьу /' Разностную схему (3.99) можно построить, применив также разложение функций в ряд Тейлора.
При использовании нерегулярных сеток построение разностных схем методом разложения функций в ряд Тейлора требует куда больше усилий, чем в случае равномерных сеток, тогда как применение интегрального метода на равномерной и неравномерной сетках требует примерно одного и того же объема аналитических вычислений. Для определения погрешности аппроксимации и проверки условия согласованности разностной схемы (3.99) необходимо все величины разложить в ряд Тейлора в окрестности узла (!,/).
Мы предлагаем читателю проделать это в качестве упражнения. В качестве предупреждения напомним, что при аппроксимации вторых производных на равномерной сетке второй порядок погрешности аппроксимации получался вследствие сокращения членов рядов Тейлора для узлов с большим и меньшим номерами. На неравномерной сетке такое сокращение нс происходит. Если расчет проводится при примерно одинаковом числе узлов разностной сетки, следует ожидать, что метод 3 учета неравномерности сетки, вызванной формой границы, должен быть наиболее точным, так как уравнение в частных производных аппроксимируется в каждой внутренней точке (в отличие от метода 2) и положение границы не меняется (как в методе 1).
Если на нерегулярной границе задана производная искомой функции (задача Неймана), то и в этом случае можно построить необходимую разностную схему, хотя алгебраические соотношения оказываются более сложными, чем в рассмотренном случае. бб Гл.
3. Основы метода конечных разностей Ряд простых примеров, иллюстрирующих методы построения разностных схем при задании на нерегулярной границе производной искомой функции, можно найти в работах 1гогзу!(зе, 'ттазотч, 1960;,)ашез е! а!., 19б7; А1!еп, 1954]. 3.5.2. нерегулярные сетки, не сзазанные с формой границы При расчете некоторых газодинамических течений для обеспечения требуемой точности часто приходится использовать сетки с мелким шагом вблизи твердой границы (стенки) или вблизи ударных волн, т.
е. в тех областях течения, где градиенты параметров велики. Для экономии времени и памяти ЭВМ целесообразно использовать более грубую сетку вдали от этих критических областей. Следовательно, для решения задач следует использовать нерегулярные или переменные сетки. Их можно строить различными методами, укажем по крайней мере два из них. 1. Можно применить преобразование координат, переводящее неравномерную сетку в физических координатах в равномерную, при этом изменяется вид уравнения в частных производных. Этот метод мы подробно опишем несколько ниже. 2. Так же как в методе 3, рассмотренном выше, можно построить конечно-разностную схему, пригодную для решения уравнения в частных производных на неравномерной сетке.
Этот подход ничем не отличается от ранее описанного, просто неравномерность сетки определяется не формой границы, а особенностями решения уравнения. В качестве примера укажем на разностную схему решения уравнения Лапласа (3.99). 3.3Л. Заклатчительные замечания Цель данного раздела состояла лишь в том, чтобы в общих чертах показать некоторые проблемы, возникающие в задачах с нерегулярными границами при использовании неравномерных сеток, а также описать методы их решения. Наше описание этого вопроса ни в коем случае не претендует на полноту. Дальнейшее его изучение быстро приводит к специальным задачам. Хорошая педагогика предполагает, что мы сначала продвинемся вперед и рассмотрим весь «лес», а уж потом вернемся к изучению отдельного дерева.
В дальнейшем мы рассмотрим некоторые вопросы использования неравномерных сеток в связи с решением конкретных задач гидродинамики н теплообмена. 87 5 3.6. Устойчивость конечно-разностных схем й 3.6. Устойчивость конечно-разностных схем Не любая согласованная конечно-разностная схема обеспечивает сходимость решения конечно-разностных уравнений к решению уравнения в частных производных. В соответствии с теоремой Лакса об эквивалентности (см. п.
3.3.4) такую сходимость решения обеспечивает лишь устойчивая разностная схема. Ватой главе мы подробно рассмотрим теорию устойчивости конечноразностных схем. Понятие устойчивости конечно-разностных схем аналогично понятию устойчивости в теории управляемых систем. Передаточная функция в теории управляемых систем играет ту же роль, что и разностный оператор в вычислительной математике.
Рис. 330. Схематическое изображение коисчво-разностной схемы как управляемой системы. Рассмотрим маршевую задачу. Пусть на п-м шаге по времени начальные значения известны, а на (и+ 1)-м шаге по времени значения этих величин надо определить. Газностный оператор можно тогда интерпретировать как «черный ящик» с некоторой передаточной функцией. Схематически такая интерпретация разностного оператора показана на рис. 3.10. Устойчивость такой управляемой системы определяется преобразованием исходных данных, проводимым «черным ящиком». Из теории управления известно, что управляемая система работает устойчиво тогда и только тогда, когда все нули характеристического полинома этой системы расположены в левой полуплоскости.
Если это не так, то входной сигнал усиливастся неправильно, а выходной сигнал бесполсзсн, так как он неограниченно растет. Теория устойчивости разиостных схем изучает способ, которым разностный оператор преобразует начальные значения величин в их значения на следующем шаге по времени, и является центральной проблемой анализа устойчивости. Начнем изучение устойчивости конечно-разностных схем с анализа простой явной схемы для уравнения теплопроводности о+! и ит — ит о Решив это уравнение относительно ич+', получим ! (3.101) аа Гл.
3. Основы метода конечных разностей Обозначим буквой Р точное решение этого разностного уравнения, т. е, решение, которое мы получили бы на ЭВМ при отсутствии ошибок округления, а буквой Л/ — решение, полученное на реальной ЭВМ. Если А — аналитическое решение исходного уравнения в частных производных, то можно записать Погрешность аппроксимации = А — Р, Погрешность округления = /1/ — Р. Устойчивость конечно-разностной схемы определяется изменением погрешности в процессе вычисления.
О'Брайен и др. [О'Вг!еп е! а!., 1950) предложили следующую классификацию устойчивости разностных схем: 1. Если полная погрешность округления растет (не растет), то разностная схема называется сильно неустойчивой (устойчивой). 2. Если отдельная погрешность округления растет (не растет), то разностная схема называется слабо неустойчивой (усгойчивой). Обычно изучают лишь слабую устойчивость, так как для ее анализа можно использовать метод разложения решения в ряд Фурье, называемый в вычислительной математике методом Неймана.
При этом предполагают, что если выполнено условие слабой устойчивости, то выполнено и условие сильной устойчивости. З.вп. Метод Фурье или метод Неймана Гассмотрим разностное уравнение (3.101). Если е — погрешность округления, то численное решение разностного уравнения можно представить в виде (3. 102) Л/= Р+ е. Так как численное решение должно удовлетворять разностному уравнению, то, подставляя (3.102) в (3.101), получим //»т/ + е»+1 /Э» е» Л //» + е» З/Э» Зе~ + /1~ + 1 / / / /е! /е/ / / 1 — 1 1 — / а/ Ьлт Точное решение Р удовлетворяет разностному уравнению (3.101), поэтому и погрешность е удовлетворяет тому же уравнению; (3.103) й 3.6. Устойчивость конечно-разностных схем 89 Так как точное решение разностного уравнения 1) и погрешность округления е удовлетворяют одному и тому же уравнению, то и растут по времени они одинаково.
![](https://s.studizba.com/z.php?f=/templates/si/i/noavatar.png&w=100&h=100&t=1"){kind=avatar}
