Fletcher-1-rus (1185917), страница 56
Текст из файла (страница 56)
7. Одномерное уравнение диффузии довольно велик, особенно на мелкой сетке. Это наводит на мысль о том, что вклады в решение различных собственных векторов, отраженные в формуле (7.60), будут затухать с различными скоростями. Если для решения системы (7.42) используется явная схема, то ограничение на величину шага по времени будет определяться значением (йеЛа~ ., Однако время, в течение которого нужно проводить интегрирование, будет определяться значением !йеЛа(;а. Отношение (йеЛа1 „/!йеЛа~ м называется коэффициентом жесткости. Если этот коэффициент очень велик, порядка 10" — !О', то система обыкновенных дифференциальных уравнений называется жесткой.
Ясно, что для решения жестких систем желательно иметь алгоритмы, несвязанные ограничениями на шаг по времени. Безусловная устойчивость подобного рода гарантируется, если численный метод обладает А-устойчивостью (РаЫ- йц!з(, 1963), т. е. если область его абсолютной устойчивости Рис. 7АО. Уетоинивоеть ти"а А(н)' включает в себя всю левую полуплоскость йеЛ б((0. Но понятие А-устойчивости содержит весьма высокое требование. Ни одна из явных схем, будь то линейная многошаговая или схема Рунге — Кутты, не является А-устойчивой.
Наивысший возможный порядок из числа А-устойчивых схем имеет неявная линейная многошаговая схема второго порядка. Родственное, хотя и менее суровое понятие, связанное с устойчивостью, — это устойчивость типа А(я) (%16!вид, 1967], представление о которой иллюстрируется на рис. 7.10. Метод обладает устойчивостью типа А(сс), если можно построить бесконечный клин, симметричный относительно отрицательной вещественной полуоси и имеющий угол полураствора а, внутри которого метод является абсолютно устойчивым. Если все собственные значения матрицы А вещественны, что имеет место для уравнения (7.41), то численный метод, обладающий устойчивостью типа А(0), будет эффективным, даже если система жесткая. При построении решения задач об установившемся течении численный алгоритм нередко основан иа эквивалентной псевдонестационарной процедуре (см.
$ 6.4). Чтобы достичь ста- $5.7. Заключение 321 ционарного решения наиболее эффективным путем, строятся неявные схемы, основанные на известных свойствах, касаюшихся устойчивости эквивалентных алгоритмов для обыкновенных дифференциальных уравнений 1Веат, %агш1пд, 19791. Неявные алгоритмы такого рода (гл. 8 — 10, 14, 17 и 18) оказываются обычно эквивалентными линейным многошаговым методам низкого порядка. Процессы приспособления вычислительных алгоритмов к специфическим свойствам векторных и параллельных процессоров [Ог1еда, Уо(81, 1985) возродили интерес к явным методам, даже с целью построения стационарных решений на основе псевдонестационарных формулировок. Такие формулировки часто основаны на использовании схем Рунге — Кутты с маршевым продвижением по времени в сочетании с построениями метода прямых (гл.
14 и 18). $7.5. Заключение Для одномерного уравнения диффузии были построены явные схемы, имеющие порядок вплоть до четвертого. Повышенные порядки точности достигались за счет взаимной компенсации ошибок усечения, связанных с дискретизацией пространственного члена дтг/дхт на различных временнйх слоях а и члена с производной по времени ду/д1 в различных пространственных точках 1.
На практике при этом требуется сделать такой специальный подбор свободных параметров, как в соотношении (7.15), чтобы члены низкого порядка в выражении для ошибки усечения стали тождественно равными нулю. При использовании явных схем для получения устойчивых решений обычно возникает необходимость выбора И, не превышающего некоторого максимального значения (посредством ограничения на параметр з).
Как правило, применение явных схем повышенного порядка связано с более суровыми ограничениями на М, обусловленными устойчивостью. В случае схемы Дюфорта — Франкела необходимо ограничивать величину М, чтобы получить приемлемо точные решения, несмотря на то что устойчивые решения могут быть получены при очень больших значениях М. Неявные схемы часто позволяют добиться безусловной устойчивости ценой добавочных вычислительных усилий. До тех пор пока неявную схему можно ограничить трехдиагоиальной формой, суммарный вычислительный труд, вкладываемый в иее, будет примерно вдвое превышать труд по решению на основе эквивалентной явной схемы.
Включение дополнительных членов на неявном (и + 1)-м временнбм слое позволяет 2! к. Флетчер, т. ! Гл.?. Одномерное уравнение диффузии 322 проявить больше гибкости в построении схем повышенного порядка. Однако, как свидетельствуют результаты (табл.
7.6), для получения значительного улучшения точности может оказаться необходимым вводить мелкую сетку. Если задача, описываемая уравнением диффузии, требует задания граничного условия для производной, то возникают дополнительные ошибки, связанные с дискретизацией граничного условия. При решении задач, определяемых параболическим дифференциальным уравнением, весьма важно, чтобы дискретизация граничных условий с производными осуществлялась с помощью формул такого же порядка точности, как и внутри области. При помощи дискретизации одного лишь пространственного члена исходное дифференциальное уравнение в частных производных (ДУЧП) преобразуется к полудискретной форме, сводящейся к 'системе обыкновенных дифференциальных уравнений (ОДУ) по времени.
В результате различные алгоритмы, разработанные специально для ОДУ, становятся пригодными для решения ДУЧП. Подчеркивается, однако, что эквивалентная система ОДУ содержит ошибки, обусловленные пространственной дискретизацией. Это обстоятельство несколько уменьшает стремление к использованию алгоритмов решения ОДУ, имеющих очень высокий порядок. Обращение к полудискретной форме оказывается удобным, когда для дискретизации пространственных членов применяется спектральный метод 5 5.6, п.
1 5.3.3 и 1 7.!.6). $7.6. Задачи Явные методы (% 7.1) 7.1. Выведите (а) формулу (7.!4), (Ь) уравнение (7,17). 7.2. Проведите приближенный подсчет числа операций для схем ВВЦП, Дюфорта — Франкела и трехслойной схемы четвертого порядка (т = О) при з = 0.30, затем сравните их вычислительную эффективность. 7.3. Постройте двухслойную пятиточечную схему, имеющую четвертый порядок точности. В окрестности границ понадобятся несимметричные пити- точечные формулы. С помощью анализа по Нейману проверьте устойчивость данной схемы. Получите численные результаты (если возможно), подтверждающие необходимость ограничений, связанных с устойчивостью и со скоростью сходиности.
7.4. Рассмотрите уравнение (7.13) с отрицательными значениями у, поставив целью убедиться в возможности построении устойчивых и точных рещений. В случае у = — 0.3 может ли схема (7.13) быть видоизменена так, чтобы она походила на схему Дюфорта — Франкела, т. е. за счет изменения р-комплекса в дискретном представлении пространственной производнойу Если да, то установите, может ли быть построена устойчивая н имеющая высокую точность схема Дюфорта — Франкела? 9 7.6. Задачи Неявные методы (9 7.2) 7.5.
Выведите формулу (7.30). 7.6. Примените анализ устойчивости по Нейману к системе (7.31) с параметром (), выраженным по формуле (7.30), чтобы в плоскости у, з определить границу области устойчивого поведения для случаев: (а) В = О, (Ь) б = 1/12, (с) б = 1/6. 7.7. При б = 1/6, у = О, 6 = 1 и з = 1(6 схема (7.31) становится явной. Сравните полученные для этого случая решения с теми, которые получаются с применением схемы ВВЦП. Исследуйте другие комбинации параметров, превращающие схему (7.3!) в явную. 7.8.
С помощью программы 01Р1М получите решения для случаев МЕ = = 4, 5 при у = 0.5, Сравните точность и скорость сходимости с данными, приводимыми в табл. 7.6. 7.9. Проведите приближенный подсчет числа операций для случаев МЕ= = 1, 4 и у = 0 в табл. 7.6 и сравните вычислительную эффективность при з =!.0 с эффективностью схемы ВВЦП при з = 0.5 н схемы Дюфорта— Франкела при з = 0.3. Граничные н начальные условия (9 7.3) 7.!О.
Видоизмените программу П1РЕХ, чтобы авеста и учесть в ней: (а) граничное условие Неймана (7.32) с первым порядком точности, (Ь) граничное условие Неймана (7.34) со вторым порядком точности, и воспроизведите результаты, приводимые в табл. 7.8. 7.11. Видоизмените программу П!Р)М, чтобы ввести и учесть в ней! (а) граничное условие Неймана (7.32) с первым порядком точности, (Ь) граничное условие Неймана (7.34) со вторым порядком точности, и получите результаты, приведенные в табл.
7.10. При ЯМАХ = 6 сравните распределение ошибки с тем, которое приводится в табл. 7.9. 7.12. Для схемы Дюфорта — Фраикела, начиная с точного решения (подпрограмма ЕХТРА) при ! = 2.00, введите схему ВВЦП на более мелкой сетке по л и 1 для получения второго слоя начальных данных. При з = 0.3 и з = 0.5 сравните точность решения по схеме Дюфорта — Франкела с ис,пользованием схемы ВВЦП на начальном участке в момент времени = 9.00 с решением, в котором два слоя начальных данных задаются точным решением при ! = 2.00 и ! = 2.00 — йй Глава 8 Многомерное уравнение диффузии 5 8.1.
Двумерное уравнение диффузии Уравнение диффузии в двух измерениях имеет вид дТ д1Т дтТ вЂ” — а — — а — г — — О. д1 "дх~ "ду (8.1) Для области, показанной на рис. 8.1, граничные условия Ди- рихле записываются в виде Т (О, у, 1) =а(у, г), Т(1, у, г) = д(у, г), Т (х, О, 1) = с (х, г), Т(х, 1, 1) =Н(х, ~), (8.2) а начальные условия — в виде Т (х, у, 0) = То (х, у). (8 8) Материал, изложенный в гл.
![](https://s.studizba.com/z.php?f=/templates/si/i/noavatar.png&w=100&h=100&t=1"){kind=avatar}
