Формула Алексея Юрьевича Виноградова для начала вычислений по методу прогонки Годунова для краевых условий любой сложности

Формулы были опубликованы в «Журнале вычислительной математики и математической физики» кажется в 1993 году (точно только то, что опубликовано до 1996 года). Формула заново восстановлена в ходе рассуждений в феврале 2006 года. 1. Краткое изложение основных матрично-векторных понятий в их классическом виде.

Составлено так, чтобы было понятно выпускникам вузов. В матричном виде система линейных дифференциальных уравнений (в том числе «жёстких») записывается так: Y(x) ’=A(x)· Y(x) + F(x), где Y ( x ) - вектор-столбец искомых функций, Y ( x ) ’ - вектор-столбец производных искомых функций, A ( x ) - квадратная матрица коэффициентов, F ( x ) – вектор внешних воздействий на систему.

Условия на левом крае записываются в виде: L ·Y(0) = L, где Y (0) - вектор-столбец значений функций Y ( x ) на левом крае x =0 , L - вектор-столбец «правой части» краевых условий левого края, L - прямоугольная горизонтальная матрица коэффициентов краевых условий левого края.

Аналогично записываются условия на правом крае (не будем их здесь записывать). 2. Про «жесткие» краевые задачи. При моделировании пространственных систем при помощи дифференциальных уравнений они иногда оказываются «жёсткими». Это, например, задачи типа расчёта на прочность тонкостенных оболочек в ракето и самолёто-строении, в кораблестроении, в трубопроводах, баках, прочие задачи для тонких и изогнутых конструкций из металла, пластика или композиционного материала. Для решения таких краевых задач с «жёсткими» дифференциальными уравнениями обычно применяют специальные приёмы-методы. «Жёсткие» краевые задачи можно решать методом Алексея Юрьевича Виноградова. Этому методу не свойственны никакие проблемы, какие есть у метода Годунова.

Познакомиться с «методом переноса краевых условий» Алексея Юрьевича Виноградова можно на страничке www . AlexeiVinogradov . narod . ru . 3. Про проблему метода прогонки С.К.Годунова.

Краевые задачи для «жёстких» дифференциальных уравнений решаются различными методами: Годунова, Абрамова, Гельфанда-Локуциевского, Алексея Юрьевича Виноградова.

Сейчас наиболее распространён метод прогонки С.К.Годунова (появившийся в 1961 году в журнале «Успехи математических наук»). На применении метода прогонки С.К.Годунова к разным задачам в разных областях науки и техники было сделано множество кандидатских и докторских диссертаций.

Рассмотрим проблему метода прогонки С.К.Годунова.

Предположим, что решается задача об оболочке ракеты. Это тонкостенная труба. Тогда система дифференциальных уравнений будет состоять из 8-ми уравнений и матрица А( x ) коэффициентов системы дифференциальных уравнений будет иметь размерность 8х8, а векторы Y ( x ) , Y ( x ) ’, F ( x ) будут иметь размерность 8х1. Соответственно матрицы краевых условий будут прямоугольными горизонтальными с размерностью 4х8. Тогда в методе прогонки С.К.Годунова для такой задачи решение ищется в следующем виде: Y ( x ) = Y 1( x )с1 + Y 2( x )с2 + Y 3( x )с3 + Y 4( x )с4 + Y *( x ) или Y ( x ) = Y матрица( x ) · с + Y *( x ), где векторы Y 1( x ), Y 2( x ), Y 3( x ), Y 4( x ) – это линейно-независимые вектора-решения однородной системы дифференциальных уравнений (системы, где F ( x )=0 ), а вектор Y *( x ) – это вектор частного решения неоднородной системы дифференциальных уравнений, а с1, с2, с3, с4 - это константы, которые надо вычислить. Здесь Y матрица( x )= | Y 1( x ), Y 2( x ), Y 3( x ), Y 4( x ) |, а с это вектор | с1,с2,с3,с4| . Но вообще то решение для такой краевой задачи с размерностью 8 (вне рамок метода прогонки С.К.Годунова) может состоять не из 4 линейно-независимых векторов Yi ( x ) , а полностью из всех 8 линейно-независимых векторов-решений однородной системы дифференциальных уравнений: Y ( x ) = Y 1( x )с1 + Y 2( x )с2 + Y 3( x )с3 + Y 4( x )с4 + Y 5( x ) c 5 + Y 6( x ) c 6 + Y 7( x ) c 7 + y 8( x ) + Y *( x ). И как раз трудность и проблема метода прогонки С.К.Годунова и состоит в том, что решение ищется только с половиной возможных констант (с1, с2, с3, с4) и проблема в том, что такое решение с половиной констант (с1, с2, с3, с4) должно удовлетворять условиям на левом крае (стартовом для прогонки) при всех возможных вариантах констант (с1, с2, с3, с4), чтобы потом найти эти константы (с1, с2, с3, с4) из условий на правом крае. То есть в методе прогонки С.К.Годунова есть проблема нахождения таких начальных значений Y 1(0), Y 2(0), Y 3(0), Y 4(0), Y *(0) векторов Y 1( x ), Y 2( x ), Y 3( x ), Y 4( x ), Y *( x ), чтобы можно было начать прогонку с левого края х=0, то есть чтобы удовлетворялись условия L · Y (0) = L на левом крае х=0 при любых значениях констант (с1, с2, с3, с4) . Обычно эта трудность «преодолевается» тем, что краевые условия рассматриваются самые простейшие и-или специально подгоняется форма записи физических параметров, через которые записывается система дифференциальных уравнений, чтобы начальные значения Y 1(0), Y 2(0), Y 3(0), Y 4(0), Y *(0) можно было угадать. То есть задачи со сложными краевыми условиями так решать нельзя: например задачи с упругими условиями на краях, не говоря уже о более сложных возможных механических формах опор. 4. Формула Алексея Юрьевича Виноградова для метода прогонки С.К.Годунова.

Выполним построчное ортонормирование матричного уравнения краевых условий на левом крае: L ·Y(0) = L, где матрица L прямоугольная и горизонтальная размерности 4х8. В результате получим эквивалентное уравнение краевых условий на левом крае, но уже с прямоугольной горизонтальной матрицей L орт (размерности 4х8), у которой будут 4 ортонормированные строки: L орт · Y (0) = L орт, где в результате ортонормирования вектор L преобразован в вектор L орт . Как выполнять построчное ортонормирование систем линейных алгебраических уравнений можно посмотреть в книгах по численным методам.

Дополним прямоугольную горизонтальную матрицу L орт до квадратной матрицы U : | L орт | U = |--------| | R |, где матрица R размерности тоже 4х8 должна достраивать матрицу L орт до невырожденной квадратной матрицы U размерности 8х8. В качестве строк матрицы R можно взять те краевые условия, то есть выражения тех физических параметров, которые не входят в параметры левого края или линейно независимы с ними. Это вполне возможно, так как у краевых задач столько независимых физических параметров какова размерность задачи, то есть в данном случае их 8 штук и если 4 заданы на левом крае, то ещё 4 есть где взять.

Завершим ортонормирование построенной матрицы U , то есть выполним построчное ортонормирование и получим матрицу U орт размерности 8х8 с ортонормированными строками: | L орт | U орт = |--------| | R орт |. Можем записать, что Y матрица(0) = ( R орт)транспонированная = ( R орт)тр. (1) Тогда подставив в формулу метода прогонки С.К.Годунова получим: Y (0) = Y матрица(0) · с + Y *(0), или Y (0) = ( R орт)тр · с + Y *(0). Подставим эту последнюю формулу в краевые условия L орт · Y (0) = L орт и получим: L орт · ( ( R орт)тр · с + Y *(0) ) = L орт.