Научный журнал
Современные наукоемкие технологии
ISSN 1812-7320
"Перечень" ВАК
ИФ РИНЦ = 0,899

ВАРЬИРУЕМОЕ КУСОЧНО-ИНТЕРПОЛЯЦИОННОЕ РЕШЕНИЕ ЗАДАЧИ КОШИ ДЛЯ УРАВНЕНИЯ ПЕРЕНОСА С ИТЕРАЦИОННЫМ УТОЧНЕНИЕМ

Ромм Я.Е. 1 Джанунц Г.А. 1
1 Таганрогский институт имени А.П. Чехова (филиал) ФГБОУ ВО «РГЭУ (РИНХ)»
Кусочно-интерполяционное приближение решения задачи Коши для уравнения переноса строится в каждом элементе разбиения прямоугольной области на подобласти. Интерполяционным полиномом Ньютона от двух переменных, представленным в форме алгебраического полинома с числовыми коэффициентами, приближается выражение частной производной по времени. Интеграл от полинома подставляется на место зависимой переменной. Процесс циклически повторяется при фиксированной степени полинома по аналогии с последовательными приближениями Пикара. Представлены доказательства сходимости кусочно-интерполяционного приближения и итерационного уточнения. Даны оценки скорости сходимости. Показано, что сконструированное приближение равномерно сходится к решению с ростом числа подобластей, при фиксированном количестве подобластей – с ростом числа пикаровских итераций. Приближение является равномерно-непрерывным в каждой подобласти. Программно реализован выбор степени полинома и числа итераций для наиболее точного приближения с наименьшей временной сложностью. Представлены результаты численных экспериментов. Согласно эксперименту в прямоугольнике единичной высоты абсолютная погрешность приближения в условиях гладкости решения составляет 10–19–10-18. Формальные аналоги метода строятся для некоторых разновидностей уравнений в частных производных, интегро-дифференциальных и интегральных уравнений.
задача Коши для уравнения переноса
интерполяционный полином Ньютона от двух переменных
кусочно-интерполяционная аппроксимация функций
кусочно-интерполяционное решение задачи Коши
двухмерный аналог последовательных приближений Пикара
1. Жамбалова Д.Б., Черный С.Г. Метод интерполяционного профиля решения уравнений переноса // Вестник НГУ. Серия: Информ. техн. 2012. Т. 10. № 1. С. 33–54.
2. Попов И.В. Построение разностной схемы повышенного порядка аппроксимации для нелинейного уравнения переноса с использованием адаптивной искусственной вязкости // Препринты ИПМ им. М.В. Келдыша. 2017. № 68. 21 с.
3. Рогов Б.В., Михайловская М.Н. Монотонная высокоточная компактная схема бегущего счета для квазилинейных уравнений гиперболического типа // Матем. моделирование. 2011. Т. 23. № 12. С. 65–78.
4. Самарский А.А., Вабищевич П.Н. Разностные схемы для уравнения переноса // Дифференц. ур-ния. 1998. Т. 34. № 12. С. 1675–1685.
5. Джанунц Г.А., Ромм Я.Е. Варьируемое кусочно-интерполяционное решение задачи Коши для обыкновенных дифференциальных уравнений с итерационным уточнением // Ж. вычисл. матем. и матем. физ. 2017. Т. 57. № 10. С. 1641–1660.
6. Ромм Я.Е., Джанунц Г.А. Кусочно-интерполяционное решение задачи Коши для уравнения переноса. Таганрог, 2019. Деп. в ВИНИТИ 20.08.19, № 68. 40 с.
7. Демидович Б.П., Марон И.А. Основы вычислительной математики. СПб.: Лань, 2016. 672 с.
8. Ромм Я.Е. Локализация и устойчивое вычисление нулей многочлена на основе сортировки. II // Кибернетика и системный анализ. 2007. № 2. С. 161–174.
9. Gasca M., Sauer T. On the history of multivariate polynomial interpolation. J. Comput. and Appl. Math. 2000. V. 122. Р. 23–35.
10. Kurganov A., Tadmor E. New high-resolution central schemes for nonlinear conservation laws and convection-diffusion equations. J. Comput. Phys. 2000. V. 160. Р. 241–282.
11. Якимов А.С. Аналитический метод решения краевых задач. 2-е изд., доп. Томск: Изд-во Том. ун-та, 2011. 199 с.
12. B?lb?l B., Sezer M. Taylor polynomial solution of hyperbolic type partial differential equations with constant coefficients. Int. J. Computer Math. 2011, V. 88. № 3. Р. 533–544.
13. Nazir T, Abbas M., Yaseen M., Lai Sh. (Reviewing Editor) Numerical solution of second-order hyperbolic telegraph equation via new cubic trigonometric B-splines approach. Cogent Mathematics. 2017. V. 4. № 1.
14. Kumar R., Choudhary A., Baskar S. Modified cubic B-spline quasi-interpolation numerical scheme for hyperbolic conservation laws. Applicable Analysis. 2018.
15. Ромм Я.Е., Джанунц Г.А. Компьютерный метод варьируемой кусочно-полиномиальной аппроксимации функций и решений обыкновенных дифференциальных уравнений // Кибернетика и системный анализ. 2013. № 3. С. 169–189.

Для уравнений в частных производных снижение погрешности приближенного решения является актуальной задачей. В частности, для уравнения переноса эта проблема обсуждается в [1–3]. Глубокому исследованию проблемы посвящены работы [4], где отмечается принципиальное значение уравнения переноса для численного решения нестационарных задач механики сплошной среды и задач газовой динамики. Границы погрешности существующих методов в прямоугольной области небольшого размера, в условиях гладкости решения, как правило, находятся в диапазоне 10–11–10-8. Для снижения погрешности в статье предлагается кусочно-интерполяционное решение задачи Коши для уравнения переноса с итерационным уточнением. Прототипом является представленный в [5] метод для случая обыкновенных дифференциальных уравнений. В случае уравнений в частных производных метод строится на основе интерполяционного полинома Ньютона от двух переменных, преобразуемого к виду алгебраического полинома с числовыми коэффициентами. Такая форма полинома позволяет естественным образом строить последовательные приближения для уточнения решения. В случае модельной одномерной задачи Коши для линейного уравнения переноса доказывается равномерная сходимость метода, оценивается скорость сходимости, для квазилинейного уравнения показана возможность аналогичных оценок. В [6] приведена программа и описан численный эксперимент, где рассматриваемая задача в прямоугольнике единичной высоты решается с погрешностью 10–19–10-18, ниже в статье приводятся результаты расширенного эксперимента. Схема построения метода формально допускает аналоги для разновидностей уравнений в частных производных, интегро-дифференциальных и интегральных уравнений.

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

1. От разностных методов разрабатываемый метод отличается построением на основе интерполяционного полинома Ньютона для двух переменных с равноотстоящими по направлению осей декартовых координат узлами. Полином варьируемой степени строится в каждой прямоугольной подобласти, на которые делится исходная прямоугольная область, метод является кусочно-интерполяционным.

2. От интерполяционных аналогов метод отличается тем, что полином имеет форму алгебраического полинома с числовыми коэффициентами. Коэффициенты конструктивно вычисляются на основе восстановления по корням полинома с помощью формул, отличных от формул Виета.

3. На основе алгебраической формы выполняется табличное восстановление первообразной от полинома, интерполирующего правую часть уравнения. Первообразная подставляется в правую часть на место искомой переменной. Такие подстановки итерируются, в результате достигается итерационное уточнение приближенного решения, выполняемое по аналогии с последовательными приближениями Пикара.

4. Отличия дают положительное преимущество относительно известных методов приближенного решения уравнения переноса, позволяющее достигать сравнительного повышения точности приближения на несколько десятичных порядков.

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

Кусочно-интерполяционное вычисление функций двух переменных

В декартовой системе координат UXT кусочно-интерполяционное приближение действительной функции u = u(x, t) двух действительных переменных в прямоугольной области

romm01.wmf (1)

строится следующим образом. Область (1) разбивается на прямоугольные подобласти Gij с пересекающимися границами:

romm02.wmf, (2)

romm03.wmf,

romm04.wmf, romm05.wmf, romm06.wmf. (3)

Интерполяционный полином Ньютона от двух переменных в подобласти (3) использует romm07.wmf узлов с треугольным расположением:

romm08.wmf (4)

При расположении узлов (4) функция интерполируется в нижней треугольной части Gij, в верхней части фактически нужна экстраполяция, что учитывается в дальнейшем. Пусть произвольно задана граница ε абсолютной погрешности приближения функции u(x, t). В Gij интерполяционный полином Ньютона romm09.wmf строится с равноотстоящими на шаги hx, ht по направлениям координат узлами romm10.wmf, где

romm11.wmf, romm12.wmf, romm13.wmf, romm14.wmf, romm15.wmf, (5)

romm16.wmf, romm17.wmf. (6)

Искомый полином примет вид [7]

romm18.wmf, (7)

где romm19.wmf – конечные разности k-го порядка. Степень полинома n выбирается одинаковой для всех Gij и минимальной при условии

romm20.wmf romm21.wmf, romm22.wmf. (8)

Проверка точности приближения (8) в каждой подобласти выполняется на множестве точек romm23.wmf, romm24.wmf, где γ ≥ 3 – параметр. Минимальность n обеспечивается следующим образом. Построение и проверка точности начинаются с n = 1 и kx = 0, kt = 0 во всех Gij из (3). При нарушении неравенства (8) хоть в одной проверочной точке какой-либо подобласти значения kx и kt увеличиваются на единицу, и проверка начинается сначала. Процесс продолжается до нарушения априори заданных границ romm25.wmf, romm26.wmf. Если в результате заданная точность не достигнута, то полагается kx = 0, kt = 0, степень n увеличивается на единицу и проверка возобновляется сначала для этой степени в тех же границах kx и kt. Описанный процесс циклически воспроизводится до нарушения априори заданной границы romm27.wmf. В качестве искомого фиксируется наименьшее n, при котором неравенство (8) выполняется одновременно во всех проверочных точках всех romm28.wmf подобластей, в соответствии с этим фиксируются текущие значения kx и kt. Предполагается, что romm29.wmf, но значения kx и kt абстрактно не ограничиваются. Полином (7) с видоизменением формул Виета [8] эквивалентно преобразуется к виду алгебраического полинома с числовыми коэффициентами

romm30.wmf. (9)

Значение (9) вычисляется по аналогии со схемой Горнера:

romm31.wmf

Для вычисления u(x, t), romm32.wmf, дешифрируются индексы Gij: romm33.wmf, romm34.wmf, int – целая часть числа, romm35.wmf, romm36.wmf, romm37.wmf, romm38.wmf. Индексы определяют адрес массива коэффициентов полинома (9) в памяти компьютера, программа вычисления приводится в [6]. С целью оценки погрешности полином (7) рассматривается в эквивалентной форме

romm39.wmf, (10)

узлы интерполяции определяются из (5), (6). Сначала (10) рассматривается в области G, при этом не используются индексы подобласти. Остаточный член интерполяции в этом случае можно представить в виде [9]

romm40.wmf, (11)

где romm41.wmf – узлы интерполяции, romm42.wmf, romm43.wmf – некоторая точка внутри G. Если не оговорено иное, внутри области (подобласти) производные понимаются в обычном смысле, на границе – как односторонние производные по направлению изнутри к границе. Для дальнейшего предполагается существование, непрерывность, и, следовательно, ограниченность в замкнутой области G всех частных производных до порядка 2n + 1 включительно. Согласно (10) достаточно было бы ввести такое предположение относительно порядка n + 1. Смысл формального завышения порядка гладкости заключается в следующем. Вследствие треугольности расположения узлов (4) остаточный член (11) фактически относится только к нижней треугольной части, а не ко всей прямоугольной области G. Для оценки погрешности именно во всей области G преобразуется расположение узлов. Интерполяционный полином строится в описанном прямоугольном треугольнике, катеты которого продолжают нижнюю и левую стороны прямоугольника G, длина каждого катета вдвое больше продолжаемой стороны. Гипотенуза треугольника, охватывающего все узлы, пройдет через правую вершину G. В описанном треугольнике узлы (4) полностью охватят область G, остаточный член (11) будет распространяться на всю эту область. Аналогичное преобразование выполняется для каждой подобласти Gij. Вследствие преобразования расстояния hx, ht между узлами в (5), при неизменности степени полинома n, окажутся вдвое больше первоначально предполагавшихся в (5) значений. Чтобы не возросла погрешность, расстояния будут сохранены такими, как они были даны в (5), тогда степень интерполирующего полинома станет равной 2n. Это соответствует количеству romm44.wmf узлов в их преобразованном расположении. В дальнейших оценках и в соотношениях вида (7)–(11) n всегда будет заменяться на 2n. С данными изменениями имеет место соотношение

romm45.wmf. (12)

Неравенство (12) получается следующим образом. Если для точки (x, t) из (11) в слагаемом под знаком Σ (c заменой n на 2n) выполняется romm46.wmf, romm47.wmf, то

romm48.wmf.

Поэтому

romm49.wmf.

Отсюда

romm50.wmf,

или,

romm51.wmf

С учетом romm52.wmf получится

romm53.wmf.

С другой стороны, если для (x, t) в рассматриваемом слагаемом с индексом i0, romm54.wmf, под знаком Σ выполняется romm55.wmf, romm56.wmf, но romm57.wmf, то romm58.wmf. Поэтому romm59.wmf, следовательно, romm60.wmf. Если под знаком Σ в (11) с рассматриваемым видоизменением romm61.wmf, но для i0, romm62.wmf, выполняется romm63.wmf, то romm64.wmf. Поэтому romm65.wmf, следовательно, romm66.wmf. В результате

romm67.wmf.

Тем более,

romm68.wmf. (13)

Неравенство (13) сохраняется в случае i = 0 и выполняется romm69.wmf. Из (11) (при замене n на 2n)

romm70.wmf. (14)

В рассматриваемых предположениях

romm71.wmf. (15)

Подстановка в (14) правых частей из (13) и C0 из (15) влечет (12). Всюду ниже предполагается, что размеры G позволяют считать расстояния между узлами меньшими единицы при всех рассматриваемых n. Узлы являются равноотстоящими вдоль направлений осей, поэтому найдутся h, q, p, такие, что

romm72.wmf. (16)

Из (12) и (16) romm73.wmf. Очевидно, romm74.wmf. В обозначении romm75.wmf,

romm76.wmf. (17)

Пусть теперь по этой же схеме в каждой подобласти Gij построен интерполяционный полином (10) с заменой показателя степени n на 2n. Тогда расстояния между проекциями узлов на оси координат уменьшатся соответственно в обратной пропорции romm77.wmf и romm78.wmf, где kx, kt из (3). В (12) тех же пропорциях соответственно уменьшатся hx, ht, поэтому romm79.wmf. Если выбрать romm80.wmf, то h в (17) обратно пропорционально 2k, в результате

romm82.wmf. (18)

В (18) C, h из (17), h < 1, i, j из (2). Таким образом, имеет место

Лемма 1. Пусть функция u(x, t) определена в G из (1), где у нее существуют и непрерывны все частные производные до порядка 2n + 1 включительно. Тогда при условии разбиения G на 22k подобластей (2), (3) в случае romm83.wmf, кусочно-интерполяционное приближение данной функции в G с помощью полиномов вида (7), взятых в степени 2n, может быть выполнено с абсолютной погрешностью (18), где шаги интерполяции hx, ht из (5) связаны с h < 1 из (17) соотношениями (16).

Рассматриваемое приближение инвариантно относительно i, j из (2), (3), поэтому в левой части (18) можно взять максимум по всем подобластям:

romm84.wmf. (19)

Отсюда вытекает

Теорема 1. В условиях леммы 1 кусочно-интерполяционное приближение равномерно сходится к функции u(x, t) в области G при k → ∞ со скоростью сходимости (19).

В каждой подобласти полином (7) может быть преобразован к виду (9) со степенью 2n без изменения оценок (18), (19). В дальнейшем используются соотношения

romm85.wmf, romm86.wmf, romm87.wmf, (20)

romm88.wmf, romm89.wmf, (21)

romm90.wmf. (22)

Значения констант и параметров определяются по ходу изложения, hx, ht в (20)–(22) пропорциональны h согласно (16).

Исходные предположения

Вначале рассматривается задача Коши для линейного уравнения

romm91.wmf (23)

где a(x, t), f(x, t) – заданные функции, рассматриваемые в полуплоскости romm92.wmf, φ(x) – заданная функция romm93.wmf [4]. Для построения метода приближенного решения выбраны прямоугольная область задания этих функций и одна из ее границ, определяемые непосредственно ниже. Ввиду применения кусочной интерполяции с оценками (14), (19) для анализа сходимости используются следующие ограничения.

I. Приближенное решение задачи (23) строится в области G из (1), объединяющей подобласти Gij из (2), (3) при romm94.wmf, областью определения φ(x) служит основание G на оси OX. Значения a и b не конкретизируются, вместе с тем всюду ниже c = 0, d = T, romm95.wmf. Предполагается, что u(x, t) принадлежит области romm96.wmf, где romm97.wmf – постоянная, значение которой может быть задано произвольно, при необходимости оно конкретно оговаривается.

II. Предполагается, что в области romm98.wmf, при romm99.wmf, существуют и непрерывны все частные производные u(x, t) до порядка 2n + 1 включительно romm100.wmf, с таким же порядком непрерывно дифференцируемы a(x, t) и f(x, t), функция φ(x) определена и непрерывно дифференцируема 2n + 1 раз на отрезке romm101.wmf.

III. Предполагается, что в подобласти romm102.wmf значение Δ может быть произвольным в границах romm103.wmf.

В предположениях I, II обусловлено существование и единственность, а также устойчивость решения задачи (23) относительно возмущения начальных данных [4].

В Gij приближенное решение строится с помощью интерполяционного полинома с итерационным уточнением, последовательно по j выполняется переход от Gij к Gi(j+1). За начальные условия в Gi(j+1) принимается приближение из Gij на смежной с Gi(j+1) границе. В G00 узловые значения интерполяции задает функция φ(x). Конкретно интерполируется romm104.wmf, что дает приближение romm105.wmf. Интеграл по времени от интерполяционного полинома принимается за приближение решения, которое подставляется в выражение romm106.wmf из (23). Полученное приближение romm107.wmf снова интерполируется, и описанный процесс повторяется с использованием (20)–(22). На практике итерации выполняются до искомой точности приближения, абстрактно их количество предполагается неограниченным. Процесс воспроизводится в каждой подобласти до полного прохода области G.

Итерационное уточнение с учетом остаточных членов интерполяции. В Gij

romm108.wmf. (24)

В дальнейшем обозначения корректируются с целью отличия от использованных при описании интерполяции функции. К нижним индексам полинома romm109.wmf из (20) добавлен индекс 2k в соответствии с числом подобластей 22k, аналогичная индексация применяется к другим выражениям. Полином romm110.wmf интерполирует не u(x, t), а подынтегральную функцию правой части (24), что отмечается слитным индексом ∂u, остаточный член интерполяции romm111.wmf оценивается из (19), –

romm112.wmf, (25)

где z, w из (20), romm113.wmf – приближение решения на предыдущей итерации. Аналогично (20), romm114.wmf преобразуется к виду (9). Решение на текущей итерации приближается интегралом от romm115.wmf, определяемым аналогично (22). Полученный полином обозначается romm116.wmf – по степени 2n, числу подобластей 22k и приближению u(x, t) путем интегрирования по t полинома romm117.wmf. Итерационный процесс примет вид

romm118.wmf,

romm119.wmf,

где romm120.wmf, romm121.wmf – остаточный член от приближения u(x, t) полиномом romm122.wmf на отдельно взятой итерации. С учетом (25) romm123.wmf, согласно предположению II romm124.wmf. Без явного выражения romm125.wmf этот же процесс запишется в виде

romm126.wmf (26)

где romm127.wmf, коэффициенты romm128.wmf связаны с коэффициентами romm129.wmf согласно (21), остаточный член оценивается из соотношения

romm130.wmf. (27)

Пусть прямоугольник Gij произвольно фиксирован и уравнение (24) рассматривается в нем с теми же начальными условиями на границе, с которыми выполняются итерации (26). В этом случае решение не является точным, что отмечается чертой сверху, аналогично отмечается приближающий его полином, из (26) следует

romm131.wmf, (28)

где использовано romm132.wmf, romm133.wmf, romm134.wmf, romm135.wmf, romm136.wmf

В предположениях I, II u(x, t) удовлетворяет соотношению

romm137.wmf

romm138.wmf, (29)

в частности (29) выполняется в случае Δx = Δ. В самом деле, по теореме о среднем, romm139.wmf, romm140.wmf. Повторное применение теоремы влечет romm141.wmf, romm142.wmf, romm143.wmf. Отсюда romm144.wmf. Область G замкнута, поэтому romm145.wmf. Очевидно, romm146.wmf достаточно взять любое Δ = const, romm147.wmf, чтобы выполнялось (29). Значение Δ в (29) можно произвольно уменьшить.

Для исследования сходимости (28) рассматривается вспомогательная задача

romm148.wmf, (30)

где romm149.wmf, f(x, t) из (23), Δ произвольно выбрано в соответствии с (29) и зафиксировано. Задача (30) является эквивалентным преобразованием (23), в области G решения задач совпадают и одновременно устойчивы к возмущению начальных данных. В предположениях I, II romm150.wmf удовлетворяет условию Липшица относительно uΔ. В условиях (30) выполнено romm151.wmf, romm152.wmf, romm153.wmf. В Gij решение задачи (30) представимо в виде

romm154.wmf, (31)

где romm155.wmf. Приближающий romm156.wmf полином будет обозначаться romm157.wmf. Полином, приближающий romm158.wmf, обозначается romm159.wmf, где z и w из (20). Значения остаточных членов интерполяции изменятся, но их обозначения сохраняются, что не должно приводить к недоразумениям. В этих обозначениях решение (31) приближает последовательность

romm160.wmf, (32)

где romm161.wmf, ввиду неточности начальных данных решение отмечается чертой. Соотношение (31) перейдет в соотношение

romm162.wmf, (33)

где romm163.wmf, romm164.wmf romm165.wmf.

Всюду ниже значение Δ в (30)–(33) рассматривается как параметр, выбор которого в границах условия (29) позволит оценить сходимость последовательности (28) на основе оценки сходимости последовательности (32).

Из (32), (33)

romm166.wmf

Замечание 1. С учетом постоянного числа подобластей 22k, k = const, не умаляя общности, можно считать, что к romm167.wmf для всех romm168.wmf из области romm169.wmf применимо условие Липшица при некотором достаточно большом значении romm170.wmf. Последовательность romm171.wmf не выводит из romm172.wmf при таком определении romm173.wmf по следующим причинам. При r = 1 полином romm174.wmf отличается от romm175.wmf не более, чем на остаточный член однократного интерполирования romm176.wmf, под интегралом можно применить условие Липшица при этом r: romm177.wmf, romm178.wmf. Тогда romm179.wmf, r = 1, при некотором romm181.wmf, таком, что romm182.wmf. Дальнейшие рассуждения строятся по индукции. Именно, с повторением данного приема для r ≥ 1 ниже доказывается сходимость romm183.wmf к romm184.wmf при r → ∞. В частности, полином данной последовательности не будет выводить из romm185.wmf romm186.wmf, поэтому к romm187.wmf применимо условие Липшица. В качестве константы Липшица romm188.wmf в дальнейшем принят ее максимум по всем подобластям из (1)–(3).

Таким образом, можно предположить, что для некоторого r ≥ 1 выполнено –

romm189.wmf (34)

В (34) и ниже romm190.wmf, M = const, romm191.wmf – отмеченное значение константы. Отсюда

romm192.wmf

и

romm193.wmf. (35)

Сначала рассматривается случай, когда в (35) погрешность однократного интерполирования не превосходит погрешности r – 1 итераций с некоторым постоянным коэффициентом, точнее,

romm194.wmf (36)

В (36) значение Q можно произвольно увеличить, и для romm195.wmf из (29) оно выбирается так, чтобы

romm196.wmf. (37)

В дальнейшем исследуется случай, когда (36) нарушается: для некоторого r0 > r окажется выполненным соотношение

romm197.wmf, (38)

при этом в силу (38) и (37) нарушение (36) влечет romm198.wmf.

Пусть сначала (36) не нарушено, в этом случае согласно (35)

romm199.wmf (39)

Не умаляя общности, можно считать Δ ≤ 1, тогда

romm200.wmf.

Если обозначить romm201.wmf, то неравенство примет вид

romm202.wmf (40)

где romm203.wmf, romm204.wmf – функция на смежной с Gij границе Gi(j–1), romm205.wmf, romm206.wmf – номер заключительной итерации, выполнявшейся в Gi(j–1). Из (40) и (27) romm207.wmf, или, romm208.wmf, здесь и ниже romm209.wmf. Тогда romm210.wmf. По индукции romm211.wmf Согласно предположению I, с учетом размеров Gij, верно неравенство romm212.wmf, в результате

romm213.wmf, (41)

и

romm214.wmf (42)

В (42) romm215.wmf, поэтому romm216.wmf. Отсюда следует, что (36) необходимо окажется нарушенным, –

romm217.wmf: romm218.wmf, romm219.wmf, romm220.wmf. (43)

Пусть (43) и соответственно (38) впервые выполнились при romm221.wmf. Для индекса r – 1 еще сохранялось (36), не выполнялось (38) и оставалось верным (39), поэтому при ? = 0

romm222.wmf (44)

и

romm223.wmf. (45)

Если предположить, что (44), (45) выполняются для некоторого ? ≥ 0, то для оценки погрешности (r + ? + 1)-й итерации можно воспользоваться соотношением (35) при соответственных индексах:

romm224.wmf.

С учетом (45) и (44)

romm225.wmf,

или, в прежнем обозначении,

romm226.wmf. (46)

Согласно (41) правая часть (46) не превосходит romm227.wmf при L = r. С учетом (43)

romm228.wmf.

Отсюда romm229.wmf, и соотношения (44), (45) сохраняются при замене ? на ? + 1. В силу индукции и согласно (37)

romm230.wmf (47)

Из изложенного вытекает

Лемма 2. Пусть в (1)–(3) зафиксировано 22k подобластей и произвольно выбрана подобласть Gij. При условии, что romm231.wmf из (33) рассматривается в Gij, romm232.wmf, с теми же начальными условиями на границе с Gi(j–1), с которыми выполняются итерации (32), romm233.wmf, romm234.wmf, последовательность (32) равномерно romm235.wmf сходится к romm236.wmf. До тех пор, пока не нарушается соотношение (36), скорость сходимости оценивается из (42). В любом случае имеет место (47).

Попутно показано, что последовательность romm237.wmf не выводит из области romm238.wmf, как и предполагалось в замечании 1.

Рассматриваемое приближение непрерывно, следовательно, равномерно непрерывно в замкнутой подобласти Gij. Поскольку для (33) и для (32) в Gij выполнено romm239.wmf, где romm240.wmf – номер заключительной итерации, выполнявшейся в Gi(j–1), то непрерывность и равномерная непрерывность приближения (32) сохраняется при переходе от Gi(j–1) к Gij и, таким образом, во всем временном слое romm241.wmf для каждого i = const.

Отсюда вытекает

Следствие 1. Кусочно-интерполяционное приближение решения задачи (30) с итерационным уточнением в каждой подобласти Gij при любом количестве итераций является равномерно-непрерывной функцией. В частности, при k = 0 приближение равномерно-непрерывно в G, при k > 0 оно кусочно-непрерывно в G и сохраняет равномерную непрерывность в слое подобластей romm242.wmf для каждого i = const.

Для romm243.wmf (47) сохранится, если взять romm244.wmf, априори в (29) взять соответственное значение Δ и в (47) указать зависящее от них r:

romm245.wmf (48)

Из того, что на выходе из Gij полином romm246.wmf задает начальные условия romm247.wmf, romm248.wmf, в Gi(j+1) , и при этом выполняется (48), вытекает оценка изменения начальных условий при переходе от одной подобласти к другой.

Следствие 2. В условиях леммы 2 romm250.wmf найдется romm251.wmf, такое что romm252.wmf из (29), romm253.wmf, максимальное отклонение последовательности (32) в Gij от romm254.wmf оценивается из (48). Этого же значения не превысит максимальное отклонение начальных условий в Gi(j+1) от тех начальных условий romm255.wmf, которые задавались бы при переходе из Gij в Gi(j+1) точным решением, взятым с начальными условиями из Gij в виде romm256.wmf.

Устойчивость решения задачи (30) сохраняется в каждой подобласти и означает, что romm257.wmf romm258.wmf, такое, что если romm259.wmf, romm260.wmf, то romm261.wmf romm262.wmf. В качестве невозмущенного решения в Gij рассматривается romm263.wmf. Его возмущение – точное (не получающееся вследствие приближения) решение в этой подобласти romm264.wmf, соответственное возмущенным начальным условиям romm265.wmf. По следствию 2 максимальное отклонение начальных данных при переходе от Gij к Gi(j+1) оценивается из (48). Аналогично, при переходе от Gi(j–1) к Gij, где это достигается за счет выбора Δ и r в Gi(j–1). Отсюда для romm266.wmf найдется δ > 0, и такие Δ, r для последовательности (32) в Gi(j–1), что выполняются соотношения:

romm267.wmf (49)

и

romm268.wmf, romm269.wmf,

поэтому, в силу устойчивости,

romm270.wmf. (50)

Пусть romm271.wmf, где r и Δ1 из (48), romm272.wmf и Δ0 из (49) . Из (48) и (50)

romm273.wmf. (51)

Поскольку число подобластей 22k зафиксировано, можно выбрать наименьшее по всем i, j значение romm274.wmf, и, при этом значении romm275.wmf, наибольший по всем i, j номер rmax, которые обеспечат выполнение (48), следствия 2, соотношений (49), а также (50), (51) одновременно во всех подобластях Gij:

romm276.wmf, (52)

где romm277.wmf, romm278.wmf. Суммой отклонений во всех подобластях слоя romm279.wmf, i = const, можно оценить максимальное отклонение кусочно-интерполяционного приближения с итерационным уточнением от точного невозмущенного решения задачи (30) во всем этом слое:

romm280.wmf,

или

romm281.wmf.

Согласно выбору максимума в (52), использованному в этих оценках, они сохраняются для любого слоя romm282.wmf с постоянным romm283.wmf из области G. Отсюда, а также из того, что по построению вспомогательной задачи (30) при любом Δ из (29) romm284.wmf, где u(x, t) – решение уравнения (23), вытекает

Лемма 3. В условиях леммы 2 кусочно-интерполяционное приближение с итерационным уточнением (32) решения задачи (30) равномерно сходится в области G к решению u(x, t) задачи (23). При этом

romm286.wmf (53)

Относительно равномерной и кусочной непрерывности последовательных приближений без изменений повторяется утверждение следствия 1, с той оговоркой, что оно относится не только к решению задачи (30), но также и к решению задачи (23).

Замечание 2. Условия и утверждение леммы 3 сохранятся без существенных изменений, если кусочно-интерполяционное приближение решения задачи (30) с итерационным уточнением применить к решению задачи

romm287.wmf, (54)

в виде

romm288.wmf. (55)

Аналогично предыдущему, доказывается сходимость (55) к решению uΔ(x, t) задачи (54), а также приближение этого решения с точностью до romm289.wmf. Из того, что romm290.wmf, где u(x, t) – решение задачи (23), с учетом (29), следует romm291.wmf. При выборе romm292.wmf получится romm293.wmf. В результате (55) будет приближать u(x, t) с точностью до ε.

Тем не менее сходимость (32) к решению задачи (30) означает сходимость непосредственно к решению задачи (23), что используется ниже.

Сходимость основной последовательности к решению линейного уравнения переноса

Для перехода от (32) к (28), в условиях леммы 3 и с учетом (53), предварительно оценивается разность

romm294.wmf. (56)

С применением теоремы о среднем,

romm295.wmf (57)

Повторное применение теоремы влечет

romm296.wmf (58)

По построению, romm297.wmf, где z, w из (20), отсюда

romm300.wmf. (59)

Правая часть (59) ограничена в силу следующих причин. Последовательность (32) сходится к решению u(x, t) задачи (23), и выполнено (53). Функция u(x, t) ограничена в G (в GΔ), в частности в Gij, поэтому при romm301.wmf значения romm302.wmf также ограничены в Gij. Поскольку romm303.wmf, то и вся последовательность romm304.wmf ограничена в Gij romm305.wmf В G, GΔ и Gij ограничена функция romm306.wmf. Как следствие, при замене romm307.wmf, u(x, t) в выражении romm309.wmf из (30) на полиномы romm310.wmf, romm311.wmf последовательность данных выражений останется ограниченной в Gij. Элементы именно этой последовательности интерполируются полиномами romm312.wmf:

romm313.wmf.

В процессе итерационного уточнения полином romm314.wmf не меняет степень и расстояние между узлами, при этом множество всех узловых значений оказывается ограниченным в продолжение всего итерационного процесса. Структура полинома аналогична (7), таким образом, ограничены его коэффициенты. В то же время коэффициенты стоят перед не меняющимися частями полинома, которые непрерывны в Gij. Отсюда полином romm315.wmf ограничен в Gij. Взятие производной от полинома с рассматриваемыми свойствами сохраняет ограниченность коэффициентов и непрерывность выражения полинома, представляющего производную. В результате romm316.wmf. С учетом romm317.wmf, и k = const, romm318.wmf. Подстановка C0 в (59), затем в (58), с учетом (57) влечет romm319.wmf. Отсюда romm320.wmf при априорном выборе romm321.wmf будет выполнено

romm322.wmf (60)

Пусть наряду с (32) рассматривается последовательность

romm323.wmf (61)

Из (33) и (61), при равенстве начальных условий в Gij,

romm324.wmf (62)

В обозначении romm325.wmf (62) эквивалентно

romm326.wmf

romm327.wmf из (56). Применение условия Липшица к romm328.wmf влечет

romm329.wmf (63)

Пусть romm330.wmf. Согласно (60) romm331.wmf. Если априори выбрать romm332.wmf так, чтобы romm333.wmf, то при romm334.wmf из (60), взятых для данного ε,

romm335.wmf. (64)

В предположении romm336.wmf из (63) следует romm337.wmf

romm338.wmf

и

romm339.wmf. (65)

Для (65) с точностью до обозначений и постоянных множителей воспроизводятся все рассуждения, преобразования и соотношения от (35) до (53) включительно. С такой оговоркой рассматриваемые преобразования инвариантны относительно Δ, выбранного в (60) и, соответственно, в (64). В (64) предполагалось romm340.wmf, romm341.wmf, или, romm342.wmf. С этим ограничением получится, –

romm343.wmf (66)

где romm344.wmf из (61), romm345.wmf – решение задачи (30). Поскольку romm346.wmf, последовательность (61) с оценкой (66) приближает решение задачи (23). По сравнению с (40)–(53) изменится коэффициент в аналоге (40). Именно, при выполнении оценки εk0 с учетом (64) получится romm347.wmf. Соответственно, romm348.wmf потребуется заменить на romm349.wmf. Если обозначить romm350.wmf, то вид дальнейших оценок сохранится. Таким образом, имеет место

Лемма 4. В условиях леммы 2 кусочно-интерполяционное приближение с итерационным уточнением вида (61) решения задачи (30) равномерно сходится в области G к решению u(x, t) задачи (23). При этом romm351.wmf выполняется соотношение (66), где romm352.wmf. Относительно равномерной и кусочной непрерывности последовательных приближений (61) без изменений повторяется утверждение леммы 3.

В (61), в силу леммы 4, для Δ из (66) romm353.wmf, r → ∞, и

romm354.wmf romm355.wmf: romm356.wmf romm357.wmf. (67)

С другой стороны, из (23), (30) и (29) следует

romm358.wmf:romm359.wmf. (68)

Пусть рассматривается аналог последовательности (61), получаемый заменой в (61) romm360.wmf на romm361.wmf, при этом разность romm362.wmf присоединяется к остаточному члену romm363.wmf и в сумме образует новый остаточный член romm364.wmf. Из (67), (68) romm365.wmf при условии romm366.wmf и romm367.wmf. В предположении romm368.wmf, получится romm369.wmf. В результате,

romm370.wmf, (69)

где

romm371.wmf, romm372.wmf, romm373.wmf. (70)

Из (33) и (69), при условии равенства начальных условий в Gij,

romm374.wmf

или

romm375.wmf, (71)

где romm376.wmf – новый остаточный член. При этом согласно (70) и (68), с учетом предположения romm377.wmf,

romm378.wmf, romm379.wmf, romm380.wmf, romm381.wmf:

romm382.wmf romm383.wmf. (72)

Тождественное преобразование (71) с учетом (56) влечет

romm384.wmf

или

romm385.wmf

где romm386.wmf. Ввиду romm387.wmf, в предположении romm388.wmf, будет выполняться romm389.wmf. Аналогично предыдущему, предположение относительно ε повлечет изменение ограничений Δ и r. Именно, с учетом (72),

romm390.wmf, romm391.wmf,

romm392.wmf,

romm393.wmf, такие, что romm394.wmf выполняется

romm395.wmf,

и, в тех же условиях,

romm396.wmf.

Остается повторить проделанные ранее рассуждения, чтобы с точностью до обозначений, значений констант и постоянных множителей вывести аналог соотношения (66), в котором romm397.wmf из (23), romm398.wmf из (69) с точностью до обозначения совпадает с romm399.wmf из (28). Отсюда рассуждения, проделанные с данными видоизменениями для romm400.wmf из (69) и romm401.wmf из (30), сохраняются для romm402.wmf из (28) и u(x, t) из (23). Сохраняются также описанные выше ограничения значений констант. Однако они дополнительно скорректируются при введении аналога соотношений (36), (37) путем выбора Q для нового остаточного члена romm403.wmf. Чтобы не усложнять обозначений, ниже данные ограничения неявно подразумеваются, но не детализируются. С этой оговоркой, в условиях леммы 2, выполняется соотношение

romm404.wmf (73)

где u(x, t) – решение уравнения (23), romm405.wmf из (28); по построению romm406.wmf зависит от ε, Δ и romm407.wmf Δ из (29); romm409.wmf существует в качестве параметра условий сходимости (69), на основе которых доказывается сходимость (28).

Имеет место

Теорема 2. Пусть в области G из (1)–(3) зафиксировано 22k подобластей. Пусть в каждой подобласти Gij уравнение (24) рассматривается с теми же начальными условиями на границе с Gi(j–1), с которыми выполняются итерации (28), при этом они сохраняются с изменением номера итерации: romm411.wmf, romm412.wmf, romm413.wmf, romm414.wmf. Тогда кусочно-интерполяционное приближение с итерационным уточнением (28) решения задачи (23) равномерно сходится в G к решению u(x, t) с оценкой (73). Относительно равномерной и кусочной непрерывности последовательных приближений (28) без принципиальных изменений формулируется аналог утверждения леммы 4, данного относительно приближений (61).

Приближение частных производных

Пусть выполнены условия леммы 2 и теоремы 2. Тогда

romm415.wmf,

где

romm416.wmf,

romm417.wmf,

romm418.wmf.

Слагаемое a1 оценивается из (29), для оценки a3 можно повторить рассуждения, проделанные для (56)–(60). При этом Δ = const можно считать столь малым, что romm419.wmf и romm420.wmf. Кроме того, romm421.wmf. В (73) для такого Δ ничто не исключает romm422.wmf, при необходимости можно указать соответственные Δε и romm423.wmf. Тогда romm424.wmf. В неравенствах можно перейти к максимуму по romm425.wmf. В результате

romm426.wmf. (74)

Производная romm427.wmf приближается с оценкой, зависящей от размера подобласти. Из (23) и (69) romm428.wmf, где учитывается (70) и c = 2C. С подстановкой (74), romm429.wmf. Здесь и в (74) вместо ε0 можно взять romm430.wmf, соответственно скорректировав romm431.wmf и romm432.wmf. Тогда

romm433.wmf

Случай квазилинейного уравнения переноса

Пусть рассматривается задача

romm434.wmf, (75)

где принимается, что f(x, t) – заданная функция в области romm435.wmf, φ(x) – заданная функция romm436.wmf, a(u, x, t) – функция, заданная в области romm437.wmf, romm438.wmf – некоторая постоянная, выбор которой оговаривается аналогично условию I. Изложенный выше метод рассматривается в предположениях I–III, с тем изменением, что a(u, x, t) задана в области romm439.wmf, и относительно (75) приняты все предположения, сделанные относительно (23). Из (75)

romm440.wmf.

Аналогично предыдущему, определяются последовательности

romm441.wmf,

romm442.wmf,

где z, w из (20). Отсюда,

romm443.wmf.

Рассматривается вспомогательная задача

romm444.wmf (76)

где

romm445.wmf, romm446.wmf, romm447.wmf из (75) romm448.wmf, romm449.wmf произвольно выбрано в соответствии с (29). В Gij при romm450.wmf (75) преобразуется к виду

romm451.wmf, (77)

где romm452.wmf. Строится последовательность

romm453.wmf (78)

Пусть в (77) и (78) начальные условия одинаковы: romm454.wmf, romm455.wmf, romm456.wmf, решение romm457.wmf отмечается чертой. Тогда

romm458.wmf (79)

Выражение под знаком интеграла преобразуется в сумму romm459.wmf, где

romm460.wmf,

romm461.wmf.

В A1 можно добавить и вычесть romm462.wmf. Отсюда

romm463.wmf

Учитывая аналог замечания 1, можно применить условие Липшица к функции romm464.wmf:

romm465.wmf

romm466.wmf.

В результате romm467.wmf, где romm468.wmf, romm469.wmf, romm470.wmf. Отсюда

romm471.wmf.

В A2 можно добавить и вычесть romm472.wmf. Тогда

romm473.wmf

С применением условия Липшица,

romm474.wmf

Как и выше, romm475.wmf. Отсюда и из (79) romm476.wmf. Ввиду произвольности romm477.wmf,

romm478.wmf. (80)

На основе (80) можно рассматривать преобразования, аналогичные выполненным для линейного уравнения, включая (36), (40), а также переход от (30) к (23), в данном случае – от (76) к (75).

Если в (75) вместо f(x, t) рассматривать f(u, x, t), то в (79) в качестве слагаемого под знаком интеграла появится разность romm479.wmf, модуль которой оценивается по условию Липшица. Это изменит значение константы в (80), но не исключает оценку сходимости метода.

Численный эксперимент

Изложенный метод реализован программно, код программы дан в [6]. Реализация содержит отклонения от исходного описания. Именно, не использован переход к описанному треугольнику области G. Интерполирование выполняется на основе (7) с переводом в форму (9) без удвоения степени полинома. Аналог (28) romm480.wmf, romm481.wmf, служит для итерационного уточнения, коэффициенты преобразуются с учетом (20)–(22). Интерполяция с итерационным уточнением выполняется в G, кроме того, в области, получающейся сдвигом G вдоль OX влево (Gleft) на половину длины катета треугольника, образованного узлами интерполяции. Этот же процесс выполняется в области, полученной аналогичным сдвигом G вправо (Gright). Для приближения решения в G выбираются полиномы из третьей слева направо четверти области Gleft (ниже они обозначаются romm482.wmf) и первой слева направо четверти области Gright (ниже romm483.wmf), а также из промежуточной между ними половины области G (ниже romm484.wmf). Этим исключаются наименее точные приближения вдоль сторон треугольника, образованного узлами интерполяции. Полный процесс приближения со сдвигами вдоль оси OX циклически повторяется со сдвигом по вертикали вдоль оси OT. Окончательные полиномиальные приближения выбираются из десятой части, отсчитанной от нижнего основания вертикально сдвигаемой прямоугольной области («нижний прямоугольный слой»), что устойчиво повышает точность приближения. В процессе смещения вдоль OT полученные на верхней границе текущего нижнего прямоугольного слоя приближения решения принимаются за начальные условия для следующего за ним временного слоя. Контроль точности привязан к текущему нижнему прямоугольному слою. Степень полинома выбирается на основе минимизации аналога невязки (ниже просто невязки) специального вида. Невязка

romm485.wmf

определяется из соотношений

romm486.wmf,

romm487.wmf,

и вычисляется в текущем (k-м) нижнем прямоугольном слое. Здесь xi – соответственные границы сдвигов вдоль OX, tk, romm488.wmf – временные границы k-го нижнего прямоугольного слоя. Геометрически невязка представляет разность объемов параллелепипедов, ограниченных криволинейными поверхностями, которые задаются кусочно-интерполяционными приближениями решения u(x, t) посредством полиномов степени n и n – 1 соответственно. Сходная форма невязки используется при выборе числа итераций r, для каждой разновидности полиномов romm489.wmf, romm490.wmf и romm491.wmf это число выбирается отдельно. Для romm492.wmf вычисляется romm493.wmf, romm494.wmf, где romm495.wmf, romm496.wmf. Аналогичные вычисления выполняются для romm497.wmf и romm498.wmf. Фиксируются значения n и r, при которых Sn и Ir принимают наименьшие значения.

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

Пример 1. Задача Коши

romm499.wmf, (81)

точное решение romm500.wmf, рассматривается как модель для сравнения численных методов [10]. Ее приближенное решение ниже дано в romm501.wmf. В табл. 1 приводятся значения абсолютной погрешности кусочно-интерполяционного решения с итерационным уточнением (здесь и ниже – реализация в Delphi, тип данных extended) в 10 равномерно распределенных точках на отрезке romm502.wmf при значениях t = 0,1 и t = 1,0 (степень интерполяционного полинома n = 13, число итераций – r = 50, область G не делится на подобласти). В первом столбце таблицы – время решения на персональном компьютере (для сравниваемых методов время на том же компьютере дано в табл. 3, 4).

Таблица 1

Абсолютная погрешность приближенного решения задачи (81)

11 s

hx ≈ 0.08

ht ≈ 0.04

t = 0.1

romm503.wmf

romm504.wmf

romm505.wmf

romm506.wmf

romm507.wmf

t = 1.0

romm508.wmf

romm509.wmf

romm510.wmf

romm511.wmf

romm512.wmf

Пример 2. Задача Коши

romm513.wmf (82)

имеет решение romm514.wmf. Приближенное решение построено в квадрате G из примера 1. Табл. 2 содержит абсолютную погрешность предложенного метода в точках, распределенных как в примере 1 (n = 13, r = 25,число подобластей Gij равно 4).

Таблица 2

Абсолютная погрешность приближенного решения задачи (82)

23 s

hx ≈ 0.04

ht ≈ 0.02

t = 0.1

5.4×10-19

0.0

2.2×10-19

4.3×10-19

0.0

t = 1.0

1.8×10-18

3.3×10-19

9.8×10-19

6.5×10-19

5.6×10-18

В системах компьютерной математики решение рассматриваемых уравнений опирается на граничные условия. Для единства условий сравнения ниже даны решения начально-краевых задач в системах MathCAD, Maрle и предложенным методом без его принципиальных изменений. В табл. 3, 4 строка (1) соответствует предложенному методу, (2) – функции pdesolve программы MathCAD (выбор разностного метода отвечает наилучшему приближению), (3) – pdsolve программы Maрle. Соотношение шагов удовлетворяет условию сходимости метода сеток для уравнений гиперболического типа.

Пример 3. Рассматривается начально-краевая задача

romm515.wmf, romm516.wmf, romm517.wmf. (83)

Погрешность ее решения в области из примера 1 сравнивается в 10 равномерно распределенных точках на отрезке romm518.wmf при значении t = 1. В предложенном методе n = 13, r = 20, область G не делится на подобласти.

Таблица 3

Абсолютная погрешность приближенного решения задачи (83)

(1)

4 s

romm519.wmf

romm520.wmf

romm521.wmf

romm522.wmf

romm523.wmf

romm524.wmf

romm525.wmf

(2)

14 s

romm526.wmf

romm527.wmf

romm528.wmf

romm529.wmf

romm530.wmf

romm531.wmf

romm532.wmf

(3)

4 min 27 s

romm534.wmf

romm535.wmf

romm536.wmf

romm537.wmf

romm538.wmf

romm539.wmf

romm540.wmf

Пример 4. В [11] рассматривается нелинейная краевая задача

romm541.wmf, romm542.wmf, romm543.wmf. (84)

Если решение взято в виде romm544.wmf, romm545.wmf, то источник E принимает вид

romm546.wmf.

В качестве опорных согласно [11] взяты значения: m = 0, romm547.wmf, romm548.wmf и b = 1. В табл. 4 приводится погрешность приближенного решения при различных значениях параметров k и A3 (первая строка таблицы) в 10 равномерно распределенных точках из примера 3. Остальные обозначения соответствуют табл. 3. В предложенном методе n = 13, r = 20, r = 35, r = 40 соответственно параметрам данных задач, деление G на подобласти не выполнялось.

Таблица 4

Абсолютная погрешность приближенного решения задачи (84)

k = 1, A3 = 0

(1)

4 s

romm549.wmf

romm550.wmf

romm551.wmf

romm552.wmf

romm553.wmf

romm554.wmf

romm555.wmf

(2)

24 s

romm556.wmf

romm557.wmf

romm558.wmf

romm559.wmf

romm560.wmf

romm561.wmf

romm562.wmf

(3)

4 min 37 s

romm564.wmf

romm565.wmf

romm566.wmf

romm567.wmf

romm568.wmf

romm569.wmf

romm570.wmf

k = 1, A3 = –1

(1)

10 s

romm571.wmf

romm572.wmf

romm573.wmf

romm574.wmf

romm575.wmf

romm576.wmf

romm577.wmf

(2)

25 s

romm578.wmf

romm579.wmf

romm580.wmf

romm581.wmf

romm582.wmf

romm583.wmf

romm584.wmf

(3)

4 min 44 s

romm586.wmf

romm587.wmf

romm588.wmf

romm589.wmf

romm590.wmf

romm591.wmf

romm592.wmf

k = 2, A3 = –1

(1)

16 s

romm593.wmf

romm594.wmf

romm595.wmf

romm596.wmf

romm597.wmf

romm598.wmf

romm599.wmf

(2)

1 min 47 s

romm601.wmf

romm602.wmf

romm603.wmf

romm604.wmf

romm605.wmf

romm606.wmf

romm607.wmf

(3)

51 min 14 s

romm609.wmf

romm610.wmf

romm611.wmf

romm612.wmf

romm613.wmf

romm614.wmf

romm615.wmf

В [11] приводится относительная погрешность 2,21 % при t = 1. Предложенный метод в этом случае дает romm616.wmf.

Пример 5. Рассматривается задача Коши

romm617.wmf,

romm618.wmf

где romm619.wmf. Точное решение имеет вид [2]:

romm620.wmf

Область приближения

romm621.wmf.

В точке с абсциссой x = 21, соответствующей пику волны при t = 1, абсолютная погрешность рекурсивной 5-точечной разностной схемы MathCAD составляет romm622.wmf, предложенный метод в этой точке даст romm623.wmf. Кусочная непрерывность приближения позволяет сравнительно адекватно отразить характер изменения решения в окрестности амплитуды. В начале и в конце волны погрешность можно снизить, окружив эти точки узкой областью с большим количеством подобластей. Так, в точке x = 11,2, соответствующей окрестности начала волны при t = 1, абсолютная погрешность составит romm624.wmf.

В аналогичных исследованиях точность данных численных экспериментов, как правило, не превосходят точность в строках 2, 3 табл. 3, 4. Так, в [12] строится решение гиперболического уравнения второго порядка с помощью полиномов Тейлора от двух переменных. Абсолютная погрешность этого метода в квадрате из примера 1 при t = 0 составляет romm625.wmf, но с ростом t растет до romm626.wmf. В [13] эта же задача решается на основе кубической тригонометрической сплайн-интерполяции. В romm627.wmf абсолютная погрешность приближения не ниже romm628.wmf. В [14] строится кусочно-интерполяционное решение уравнения переноса с помощью кубической B-сплайн квазиинтерполяции, ее абсолютная погрешность в квадрате из примера 1 не ниже romm629.wmf при t = 1.

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

Замечание 3. В табл. 1–4 время решения задачи предложенным методом дано без использования программного выбора степени полинома и числа итераций. С программным выбором n и r время существенно возрастает, наиболее долго решается задача (82) – romm630.wmf. В целом с таким выбором время решения занимает промежуточное значение между MathCAD и Maрle в границах погрешности, приведенной в таблицах.

Алгоритм выбора параметров можно выполнить параллельно по всем сравниваемым значениям невязки Sn одновременно для всех степеней полиномов n (в силу взаимной независимости сравниваемых значений), аналогично, по всем Ir для всех значений r. По параллельно найденным значениям без труда определяются параметры минимизации погрешности при наименьшей временной сложности. Вместе с тем вычислительный алгоритм решения рассматриваемой задачи Коши является взаимно независимым по временным слоям подобластей romm631.wmf, i = const, поэтому решение распараллеливается по всем номерам слоев romm632.wmf, ускоряя последовательный алгоритм пропорционально 2k. Распараллеливание допускает также алгоритм вычисления коэффициентов и значения интерполяционного полинома Ньютона. Для полинома от одной переменной временная сложность оценивается в [15], для полинома от двух переменных – в [6].

Об аналогах метода для уравнений других видов

Интерполяционным полиномом (7), преобразованным к виду (9), приближаются частные производные решения уравнения гиперболического типа второго порядка и выше, итерационное уточнение в этом случае выполнимо с помощью первообразной соответственной кратности. Для линейного гиперболического уравнения второго порядка такой подход реализован в [6]. Формально рассмотренный метод переносится на уравнения порядка выше первого на основе их преобразования в систему уравнений первого порядка. Пусть, например, рассматривается задача Коши для уравнения гиперболического типа

romm633.wmf,

где a, b, c, d, g, f – заданные функции независимых переменных x и t, ab > 0. Пусть приближенное решение ищется в прямоугольнике G из (1)–(3) при начальных условиях

romm634.wmf,

где φ и ψ – заданные функции romm635.wmf. Замена переменных влечет

romm636.wmf,

romm637.wmf.

В каждой подобласти значения y и z из первых двух уравнений можно интерполировать полиномами (7) с преобразованием к виду (9). Для решения romm638.wmf в качестве приближений частных производных подставляются соответственные частные производные преобразованных полиномов. Кроме того, используются полиномиальные приближения y, z. Процесс можно циклически повторять до достижения искомой точности. Вопрос о сходимости при данном подходе требует отдельного исследования.

Аналогично итерационному уточнению можно строить последовательность приближений решения интегро-дифференциальных уравнений. Пусть, например, рассматривается задача Коши для линейного уравнения с интегральным оператором типа Вольтерры

romm639.wmf. (85)

Предполагается, что A(x) и f(x) непрерывны при romm640.wmf, romm641.wmf непрерывно при romm642.wmf. В качестве нулевого приближения принимается romm643.wmf. С подстановкой romm644.wmf подынтегральную функцию можно интерполировать полиномом (7) с преобразованием к (9): romm645.wmf, где romm646.wmf, romm647.wmf (здесь и ниже не индексируется подобласть и соответственные ей выражения). Тогда romm648.wmf. Если romm649.wmf, то romm650.wmf. Правую часть приближения romm651.wmf можно интерполировать полиномом Ньютона от одной переменной с равноотстоящими узлами, преобразовав его согласно [5] к виду алгебраического полинома с числовыми коэффициентами, romm652.wmf, romm653.wmf, h – шаг интерполяции. В результате romm654.wmf. Первообразная от romm655.wmf принимается за первое приближение: romm656.wmf, или, romm657.wmf. Циклическое повторение описанного процесса с рекуррентными подстановками влечет romm658.wmf (знак romm659.wmf ставится, поскольку не учитываются остаточные члены интерполяции), отсюда romm660.wmf. Сходимость метода требует отдельного исследования.

Задача (85) сводится к интегральному уравнению Вольтерры второго рода

romm661.wmf,

где K(x, s) и f(x) – непрерывные функции при romm662.wmf. Метод последовательных приближений для этого уравнения сходится к единственному непрерывному решению. Если в качестве нулевого приближения взять romm663.wmf, то romm664.wmf. Подынтегральную функцию с подстановкой y1 можно интерполировать полиномом (7) и преобразовать его к виду (9): romm665.wmf, romm666.wmf, romm667.wmf. Отсюда romm668.wmf. Если romm669.wmf, то romm670.wmf. Правую часть приближения y(x) можно интерполировать полиномом Ньютона от одной переменной, преобразовав его к виду romm671.wmf. Тогда romm672.wmf, romm673.wmf, и в качестве второго приближения принимается romm674.wmf. В продолжение процесса, с рекуррентными подстановками, получится romm675.wmf, откуда romm676.wmf. Сходимость в данном случае можно анализировать, сопоставляя yr(x) со сходящейся последовательностью приближений, каждое из которых имеет точное аналитическое выражение.

Для неоднородного уравнения Фредгольма 2-го рода romm677.wmf метод последовательных приближений сходится лишь при условии romm678.wmf, где romm679.wmf. Поэтому аналог рассмотренной схемы допустимо строить только при выполнении этого условия.

Заключение

Изложен метод кусочно-интерполяционного решения задачи Коши для уравнения переноса с итерационным уточнением. Решение и его частные производные интерполируются полиномами Ньютона от двух переменных, преобразуемыми к виду алгебраических полиномов с числовыми коэффициентами. На этой основе строится последовательное уточнение решения, которое сходно с двумерным аналогом последовательных приближений Пикара. В случае линейного уравнения переноса метод равномерно сходится к решению в прямоугольной области, строятся приближения частных производных. Приближение решения равномерно-непрерывно в данной области, кусочно-непрерывно – при ее делении на подобласти. Согласно эксперименту решение линейного и квазилинейного уравнений переноса в условиях гладкости приближается с абсолютной погрешностью 10-19–10-18. Метод отличается структурой и сравнительно малой погрешностью, что позволяет детализировать численное моделирование процесса переноса волны.


Библиографическая ссылка

Ромм Я.Е., Джанунц Г.А. ВАРЬИРУЕМОЕ КУСОЧНО-ИНТЕРПОЛЯЦИОННОЕ РЕШЕНИЕ ЗАДАЧИ КОШИ ДЛЯ УРАВНЕНИЯ ПЕРЕНОСА С ИТЕРАЦИОННЫМ УТОЧНЕНИЕМ // Современные наукоемкие технологии. – 2020. – № 1. – С. 21-46;
URL: https://top-technologies.ru/ru/article/view?id=37897 (дата обращения: 17.08.2022).

Предлагаем вашему вниманию журналы, издающиеся в издательстве «Академия Естествознания»
(Высокий импакт-фактор РИНЦ, тематика журналов охватывает все научные направления)

«Фундаментальные исследования» список ВАК ИФ РИНЦ = 1.074