1626435584-7c6402f545ecf856225d6cf8d21519c9 (844233), страница 97
Текст из файла (страница 97)
В первом способе линеаризация выполняется так, как описано выше, а к давлению добавляется вязкий член, взятый с предыдущей итерации: Ы Р(Р )+))ор (о +) о» ) и) ш и) и) *и) (77) Это означает, что вязкость включена в итерационный процесс методом последовательных приближений.
Такой способ прост„но ухудшает сходимость итераций: уменьшает скорость сходимости и усиливает ограничение на шаг т, хотя не слишком сильно. Второй способ — полная линеаризация — сложнее, но надежнее. Линеаризируя уравнение (71а), учтем зависимость а не только от р, но и непосредственно от о через вязкость (?1г): При этом вместо (72) и (74) получаются более громоздкие выражения, которые мы не приводим. Однако такой процесс является чисто ньютоновским и хорошо сходится. Не и вот ер ми чески й ел уча й требует включения в итерационный процесс уравнения энергии (71д), что часто делают способом двухкруговых шпераций (последовательных прогонок).
Сначала считаем энергию е„(или температуру) известной во всех точках нового слоя. Тогда в каждой точке р (е„, р„) =р„(р„), т. е. применимы формулы изотермического случая (74), (76); по ним проводят первый малый круг итераций. Когда эти итерации сойдутся, полученные значения г, о, р подставляют в уравнение энергии (71д). Неизвестными в нем остаются значения е; их можно определить, линеаризируя уравнение (71д) с учетом зависимости р(е): о»»19»ы бе»»1+ (1 + о»»тз» Р»9») бе» о»9»-т бе»-) = » е» + 4 (о»+) о») 4 (о»+1 о»)+ в" гн + — [О»(у»-) +у» ) — о»мЯ» +у»+))1~ 2т 2)» ) дв )» (79) Итерации (79) образуют второй малый круг.
На каждой итерации трехточечное уравнение (79) решается прогонкой. Найденные значения е„передают в уравнения (74), (76) н снова проводят первый малый круг итераций и т. д, Это взаимное согласование уравнений импульса и энергии составляет большой круг итераций. 4ЗО (гл. Хп! ГИПЕРБОЛИЧЕСКИЕ УРАВИЕЬ!ИЯ Обычно считают нормалы!ым, если малые круги сходятся за 3 — 5 итераций, а большой круг — за 2 — 3. Большее число итераций указывает на целесообразность уменьшения шага т. Замечание. Можно провести итерации в один круг, если полностью линеаризировать систему (71), считая )т = р (е, р), Однако при этом получаются существенно более громоздкие уравнения в вариациях, для решения которых надо применять матричную прогонку (см, дополнение к 130)).
Устойчивость. Методом разделения переменных в линейном приближении можно показать, что схема (71) безусловно устойчива. Таким образом, шаг т ограничивает только условие сходимости итераций при решении нелинейной системы (71). Аппроксимация и сходимость. Схема (71) не симметрична по ( и поэтому даже на гладких течениях имеет аппроксимацию О (т+/тз). Тем самым, на гладких течениях схема <крест» может оказаться более точной.
Однако при расчете течений с ударными волнами и другими особенностями неявная схема дает сушественио лучшие результаты, чем схема «крест». Поэтому она широко применяется в практике вычислений„ особенно в кбольших задачах». Сходимость схемы (71) строго не доказана, но многократно проверена на сложных задачах-тестах с известными точными решениями. б. О других схемах. Схемы (бб) и (7!) являются одпородиыми. Имеется много близких к ним алгоритмов, отличающихся деталями написания отдельных членов разностных схем или другой организацией итерационных процессов решения нелинейных разностных уравнений.
Из них следует отметить лолногглью коньервотивныг схемы, в которых автоматически выполняются разностные законы сохранения не только массы, импульса и полной энергии, но также законы сохранения энтропии и внутренней энергии. В настоящее время построены полностью консервативныс схемы лля задач одномерной газодинамики в лагранжевых и зйлеровых переменных, задач магнитной газодинамики и двумерных газодинамических течений (подробнее см. в 134)).
Есть иначе построенные однородные схемы. Из них отметиы схему распада разрыла. Она составлена так, что в акустичсскол! приближении *) переходит в явную схему бегущего счета для инвариантов (ЗЗ), обладающую хорошей аппроксимационной вязкостью. Благодаря этому схема позволяет рассчитывать любые разрывы без введения псевдовязкости. В акустическом приближении схема распада разрыва монотонна; в газо- динамике на сильных ударных волнах возможна немонотонность, хотя фактически она невелика. Схема имеет аппроксимацию О (т+(!), поэтому для расчета гладких течений она невыгодна.
Но фронгы ударных волн она воспроизводит хорошо, с малым сглагкиванием. Схема распада разрыва — явная и имеет ограничение на шаг типа От ~ Лг, где 0 — скорость ударной волны. Это ограничение, а также громоздкость схемы препятствуют широкому ее применению. *) Если р, р, е лишь слабо колеблются около равновесных значений, то уравнения газодинамики переходит в уравнения акустики (см., наприь!ер, 1401), 451 ЗАДАЧИ ЗАДАЧИ 1. Составить схему «крест» для задачи (1) при неравномерных сетках по х и 1 и исследовать аппроксимацию схемы. 2. Найти невязку схемы (!2). 3.
Для волнового уравнения (1) составить схему с весом и на шаблоне рис. 101 и провести исследование этой схемы; показать, ! ио пРи а -. «7» схема бе»Условно неУстойчива, пРи о -= :=" ~,'« — безусловно устойчива и при о ) т(а обладает аппроксимацнонной вязкостью. 4. Установить аппроксимацию схемы (19). 5. Проверить исследование устойчивости схемы (!9), данное в 9 1, п.
3. 1-б 6. Доказать, что схема (26) имеет аппроксимацию Рис. 101. О (тз+ А«). 7. Составить схему типа «крест» для задачи (!8), приписывая значения рз Узлам сетки, а з — центРам Ячеек к„ьм, 1 +у. написать ДлЯ нее начаЛьэ ные данные точности 0 (т«+Аз) 8. Провести полное исследование схемы (ЗЗ). 9. Рассмотреть, как в схеме (33) вычисляется разностное решение в граничных узлах. 10. Вывести формулы циклической прогонки для случая матрицы, изобра»кенной на рис.
95. 11. Исследовать устойчивость многомерной схемы с весами (40). 12. Провести полную лннеаризацию сисгемы (7(а) — (71г) для случая изотермической газодинамики с учетом вязкости и свести задачу к решению трехточечного уравнения относительно бо. Помимо однородных схем сушествуют схемы с явным выделением особенностей, в которых точно прослеживается движение всех сильных и слабых разрывов. Одна такая схема предложена и подробно описана в (87). Но такие схемы очень сложны, и их применяюг только в тех случаях, когда требуется особенно высокая точность расчета.
ГЛАВ А Х1Ч ИНТЕГРАЛЬНЫЕ УРАВНЕНИЯ В главе Х1Ч рассмотрены простейшие методы решения интегральных уравнений. Корректно поставленным задачам посвшпен й !. В нем изложены некоторые типичные постановки задач и даны методы нх решения: разностный метод и некоторые приближенные методы. В 1 2 рассмотрены некорректно поставленные задачи для линейных интегральных уравнений первого рода. Изложена теория построения регуляризирушщих алгоритмов по А. Н. Тихонову. Для некоторых некорректных задач, возникших в предыдущих главах, даны алгоритмы решения, доведенные до схем численного расчета. $1.
Корректно поставленные задачи 1. Постановки задач. Интегральным называют уравнение, в котором неизвестная функция и (х) стоит под знаком интеграла. Одномерное нелинейное интегральное уравнение имеет вид ь ')К(х, $, и(1)) д$=г" (х, и(х)), а~х~й, (1) а где ядро К(х, $, и) и правая часть г (х, и) — заданные функции. К интегральным уравнениям приводят многие физические задачи. Так, задача восстановления переданного радиосигнала и (1) по принятому сигналу ) (1) сводится к решению интегрального уравнения типа свертки: ~ К (1 — т) и (т) ат =- ) (1), (2) о где ядро К(й) зависит от свойств приемной аппаратуры и среды, через которую проходит сигнал.
Заметим, что даже для задач, записанных в терминах уравнений в частных производных, первичной обычно является формулировка в виде интегральных законов сохранения, т. е. интегральных уравнений. В предыдущих главах такие формулировки использовались, например, для построения консервативных разностных схем. 5 н КОРРЕКТНО ПОСТАВЛЕННЫЕ ЗАДАЧИ Интегральные уравнения в некоторых отношениях удобнее дифференциальных. Во-первых, интегральное уравнение содержит в себе полную постановку задачи. Например, интегральное уравнение х (х)=и,+~(Я, (В))а$ Р) х~ эквивалентно задаче Коши для дифференциального уравнения аи (х) =Г" (х, и), и(х) =и,.
(4) Тем самым, для уравнения (3) не требуется задавать никаких дополнительных условий, начальных или граничных (см. также задачу 1). Во-вторых, в интегральных уравнениях переход от одной переменной ко многим является естественным. Так, многомерным аналогом (1) является уравнение ~ К(х, $, и(й)) г(4=г (х, и(х)), о (5) х=(х„х„..., х,) я б(х), отличающееся от (1) только тем, что интегрирование проводится по многомерной области 6. Поскольку оба уравнения не требуют дополнительных условий и полностью определяют задачу, аналогия является полной.
![](https://s.studizba.com/z.php?f=/templates/si/i/noavatar.png&w=100&h=100&t=1"){kind=avatar}
