1626435584-7c6402f545ecf856225d6cf8d21519c9 (844233), страница 83
Текст из файла (страница 83)
Достаточно соответствующим образом записать разностный аналог пространственной производной: 2 脄— оа и,„— х„ з — х, , ( х„„з — ха нл пл-т 1 х„— ха, /' и (х, 1) = ~~ ссеехр~ — й —,(~з)п —, ласт ~ . лсх е=! *) Это нужно в задачах фильтрации нефти в пласте при многолетней эксплуатациии скважин, задачах прогрева слоя печной мерзлоты и ряде других. В этом случае схема по-прежнему сходится в '1 1, с погрешностью 0(йа); однако доказательство этого утверждения значительно сложнее и проводится методом энергетических неравенств (см. 1301). Подведем итоги. Поскольку погрешность почти для всех значений о есть 0(Ьз), то для получения хорошей точности при расчете по схемам (6) надо брать довольно малый шаг Ь.
В этом случае явная схема устойчива при настолько малом т-=.йа/(3гг), что для доведения расчета до заданного момента Т требуется сделать очень много шагов по времени, т. е. выполнить большой объем вычислений. Поэтому явные схемы для решения параболических уравнений почти никогда не употребляются. Обычно для расчетов берут двуслойные неявные безусловно устойчивые схемы. Чаще всего используют симметричную схему или схему повышенной точности, обеспечивающие хорошую точность расчета при не слишком малых шагах т и й.
Чисто неявная схема в случае а=сонэ( редко употребляется из-за невысокой точности (хотя при Й=Й(а) она часто выгодна благбдаря своей монотонности). 3. Асимптотическая устойчивость неявной схемы. Исследуем, при каких условиях схема (6) позволяет рассчитывать задачи с нулевыми краевыми значениями для очень больших промежутков времени е), т. е. каковы условия асимптотической устойчивости схемы. Выход решения параболического уравнения (1) на асимптотику при 1 — ь.со определяется скоростью затухания начальных данных.
Приведенное в главе 1Х разложение решения и(х, г) в ряд Фурье (9.7): 275 ОДНОМЕРНЫЕ УРАВНЕНИЯ показывает, что медленнее всего затухает первая гармоника. Ей соответствует множитель роста лт ~ йлтт р,=ехр( — й —,т)=1 — — ", +О(т'). Чтобы схема (6) была асимптотически устойчивой, ее множители роста (13) не должны превосходить по модулю величины р„ т. е. должно выполняться условие 5!Пт— лай йлтс 2а йлтт — 1-1- — "~1— ат йт лай а5 1 — — ', 1 д-)() — 1; (18) — + 6 5!П5— 4йт 2а здесь мы преобразовали выражение (13) для ра К более удобному виду.
Разумеется, достаточно выполнения этих неравенств с точностью до членов О (т'), потому что наличие таких членов приведет к умножению амплитуд гармоник на величину 11+ О (тт)1нт = = 1+(О(т)„чем при т-5.0 можно пренебречь, даже если ! велико. Нетрудно проверить, что правое неравенство (18) всегда выполняется. В самом деле, р монотонно убывает при увеличении з)п (лдй!(2а)), т.
е. при увеличении д. Поэтому наибольшим является р„которое с учетом малости Ь равно ( Ь)га)5 йл (йт/4йт)+а (лй)2а)5 ат и совпадает с р, с точностью до членов О(т'). Рассмотрим левое неравенство (18). Величина ра минимальна при д=!У' — 1, когда з!Н(л5)й)2а)-1. Подстановка этого значения в левое неравенство (18) после несложных выкладок приводит к условию асимптотической устойчивости л5й~ йт о) — +— 2 4ат 4йт (19) Оно несколько более жестко, чем условие обычной устойчивости (14). Его можно переписать в следующем виде: т — — (2о — 1) т — ( — 7! «= О. 2ат !ай(5 л5й '1 лй7'— аз Г 5 ! лй !51 тт,= —,— ~ (2о — 1) + у (2о — 1)'+( — ) ~. Стоящий слева квадратный трехчлен отрицателен, если т лежит между его корнями: 376 ПАРАБОЛИЧЕСКИЕ УРАВНЕНИЯ 1гл.
х! Один из корней отрицателен, а другой положителен. Поэтому условие асимптотической устойчивости (!9) принимает вид т ~ —, 12о — 1+ 1/ (2о — 1)'+ (пй/а)'1. (20) В частности, симметричная схема асимптотически устойчива не при любом т, а только при т ~ — (о= !Я. (21) Таким образом, схема (б) при любом о формально является лишь асимптотически условно устойчивой. Однако фактически устойчивость условна только при о="'/,+О (й), когда ограничение на шаги принимает вид т.-=.й сопз1. Если же о)",.„то условие (20) требует, чтобы выполнялось неравенство т == == [2а' (2о — 1)/изй) = сопз(, и по существу схема является асимптотически безусловно устойчивой.
Замечан ие. При больших / схемы с о>!/, дают низкую точность. Поэтому для таких расчетов обычно используют схему со=/м 4. Монотонность. Точное решение и(х, () уравнения и,=йи,„ при определенных условиях сохраняет монотонность. Например, если начальные данные монотонны и температура на концах отрезка постоянна, то профиль температуры будет монотонен в любой момент времени. То же будет при постановке задачи Коши на бесконечной прямой. Выясним, сохраняет ли схема (б) монотонность решения. Ограничимся задачей на бесконечной прямой, хотя при использовании результатов надо помнить, что краевые условия тоже влияют на монотонность (если разностное краевое условие не точное, то его неудачное написание может привести к иемонотонности схемы), Для случая о=0 результат почти очевиден. Получающаяся при этом явная схема имеет форму (10.34)! 1!о =1 Ит Ат А2 ' ' -' 7Я ' Необходимым и достаточным условием монотонности является неотрицательность коэффициентов р!.
Видно, что если выполнено условие устойчивости этой схемы 2йт ( й', то коэффициенты иеотрицательны и схема монотонна. В противном случае явная схема немонотонна, 377 ОДНОМЕРНЫЕ УРАВНЕНИЯ При а~О запишем неявную схему (6), полагая у„=0 н выделяя преобладающий член на новом слое: ! — а У.=~2+е — „) ~У.-е+У.ы+ —,(Уп-~+У.-)+ + ~ета 2 а ~ Уп~ (22) Напишем для уп, и уп„аналогичные выражения и подставим их в правую часть (22). При этом появятся другие значения с нового слоя; будем их исключать тем же способом.
Коэффициенты при значениях у„ а на новом слое в правой части будут при этом убывать в геометрической прогрессии. Поэтому в пределе соотношение (22) перейдет в явную схему вида (10.34) с бесконечной суммой, т. е. с бесконечной зоной влияния. Очевидно, если выполнено условие Уп=()пУп+ .~~ Рс (Уп ~+Упн), (24а) где 4ьт 4аьт (О 1 т(А+ ) ( 1 т(Л+т)пг 4апт ()с=()г-е „+, при т у = )/й'+ 4а/гт.
(246) Очевидно, (),зе0 при 1==-:1 и отрицательным может быть только коэффициент ()и. Он неотрицателен„если (2 — а) йп т =4А() — а)п Это условие необходимо и достаточно для монотонности схемы (6); оно несколько слабее ограничения (23). Таким образом, неявные схемы монотонны только при очень малом шаге по времени т а'.
По абсолютно устойчивым неявным схемам расчеты обычно проводят с шагом т й, не гаран- Ап 2А 0-а) ' (23) то все коэффициенты в (22) неотрицательны. Тогда все коэффициенты в соответствующей явной схеме также будут неотрицательны. Следовательно,. неравенство (23) есть достаточное условие монотонности схемы (6). Можно получить необходимое и достаточное условие монотонности, приведя схему (22) после выполнения громоздких выкладок к явной форме (10.34): 378 (гл. хг ПАРАБОЛИЧЕСКИЕ УРАВНЕНИЯ тируюшим монотонности.
Единственное исключение †чис неявная схема с о =- 1, которая монотонна при любых шагах. Напомним, что достаточно гладкое решение на подробных сетках можно хорошо находить и по немонотонным схемам. На грубых же сетках, особенно при разрывных начальных данных, симметричная схема может привести к «разболтке» счета. Чисто неявная схема даже в этих условиях дает плавно меняющееся разностное решение, хотя точность его невысока. 3 а меч а ни е. Монотонные схемы для параболического уравнения могут иметь второй порядок точности по пространству.
Но, как и для уравнения переноса, для параболического уравнения не известно ни одной монотонной схемы, которая имела бы второй порядок точности по времени (хотя никаких теорем о невозможности построения таких схем не доказано). 5. Явные схемы. Явные схемы имеют важное достоинство: они просто записываются и легко программируются на ЭВМ. Поэтому предпринималось много попыток построить для параболического уравнения и,=йи,,-,+) хорошую явную схему. Однако все эти попытки были неудачными.
Например, Ричардсоном была предложена трехслойная схема, использую. шая шаблон рис. 77 с аппроксимацией производных двусторонними разностями: 1 й 2 .' (уп уп) Аз (уп-т 2уп+ упю) +(и' (26) Из симметрии схемы легко усмотреть, что локальнан погрешность ее аппроксимации есть 0(та+Аз). Однако схема Ричардсона непригодна для расчетов, ибо она безусловно неустойчива. В самом деле, используем метод разделения переменных и сделаем подстановку у„=ехр ((ухп),.уп=рду; поскольку схема трехслойна, надо дополнительно положить уп=(1/рд) у„. Тогда для множителя роста получим квадратное уравйенйе 8йт ., ЧА р»+ — -р яп» вЂ” — 1=0, !Р д 2 (27) Рис.
77; один из корней которого нри любом д ~0 по модулю больше единицы на величину 0 (йт)дз). Дюфорт и Франкел в 1953 г. видоизменили схему Ричардгона, заменив в правой части (26) величину уп на (уп + уп)72: ! - й 2т(Уп Уп)=й» (Уп-т Уп Уп+Упдт)+гп (28) Эта схема также явно разрешается относительно уп. Методом разделения переменных нетрудно показать, что она безусловна устойчива. Однако погрешность аппроксимации схемы (28) равна 0(тз+/гз+т»дгз), т. е.
аппроксимация условная, Поэтому сходимость имеет место, только если (т!А) — «.О прн Ь-»0. Фактически, чтобы в расчетах по схеме (28) получить точность 0(йе), надо положить йт — Ь', как в явной схеме (6). Правда, коэффициент пропорциональности а=(йг(Ь») можно брать любым, ибо его величина влияет только на ОДНОМЕРНЫЕ УРАВНЕНИЯ 379 (29а) точность расчета, а не на устойчивость. Позтоиу схема Дюфорта †Френке удобнее явной схемы (6), но ненамного. Плохие качества явных схем обусловлены одним принципиальным ограничением: явная скелга для параболического уравнени может сходиться, только если (т(п) ь-О при й- О.
В самом деле, пусть решение в точке нового слоя выражается через г точек исходного слоя, т. е. через отрезок длиной гп (рис. 78). Тогда оно выражается через отрезок нулевого слоя длиной тгй=г(й!т; этот отрезок будет зоной влияния. Для точного решения зона влияния бесконечна.
![](https://s.studizba.com/z.php?f=/templates/si/i/noavatar.png&w=100&h=100&t=1"){kind=avatar}
