Метод Алексея Юрьевича Виноградова для решения краевых задач

Пожалуйста, пишите на адрес AlexeiVinogradov @ yandex . ru о том, где может найти применение метод А.Ю.Виноградова: например, в каких областях он может быть сопоставлен с методом С.К.Годунова помимо тонкостенных оболочек J . Хорошо было бы составить на этой страничке перечень задач разных областей знаний, где нужно преодолевать трудности неустойчивого компьютерного численного решения краевых задач: то есть хорошо было бы собрать координаты специалистов (и описания их задач), которые сейчас пользуются методом С.К.Годунова потому, что им можно предложить попробовать простенький и элегантно-эффективный метод А.Ю.Виноградова J . Пишите: если Вы захотите, то Ваши сообщения будут выложены на прилинкованных страничках J или будут даны ссылки на Ваши странички или специализированные сайты J . Дополнительно хотелось бы узнать в каких вузах-университетах читаются соответствующие численные методы, чтобы предложить студентам и аспирантам пробовать проводить расчёты методом А.Ю.Виноградова (в сравнении, например, с методами Годунова или Абрамова или Гельфанда-Локуциевского и т.п.) J . 1. Далее идёт краткое изложение основных матрично-векторных понятий в их классическом виде вне рамок метода Алексея Юрьевича Виноградова, а как это принято для любых краевых задач.

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

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

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

Многообразие всех возможных линейно независимых решений системы дифференциальных уравнений в матричном виде выглядит как квадратная матрица, которая называется матрицей Коши или интегралом Коши и для обозначения можно использовать просто букву К. Суть (объяснение) матрицы К: она состоит из векторов типа Y , но эти векторы типа Y имеют принципиально различные варианты поведения, то есть они являются «линейно независимыми»: K = [ Y1 Y2 Y3 Y4 Y5 Y6 Y7 Y8 ] При помощи матрицы Коши (то есть матрицы К) можно (на основе дифференциальных уравнений) установить связи между левым краем задачи и правым краем задачи: K(x_left x_ right). То есть вычислив некоторым из известных способов матрицу К можно записать, что: Yright = K(x_right x_left) · Yleft Далее для краевых задач решение находится привнесением краевых условий в совместное рассмотрение с полученной матрицей К: L ·Yleft = L R ·Yright = R Yright = K(right left) · Yleft Обычно в краевых задачах традиционно требуется найти либо искомый вектор левого края Yleft либо искомый вектор правого края Yright , чтобы привести краевую задачу как бы к задаче не краевой, а к задаче Коши (не к матрице Коши, а задаче Коши – к задаче с начальными, а не краевыми условиями). То есть чтобы можно было искать решение задачи в какой-либо точке по формуле решения задачи Коши: Y ( x ) = K ( x 0) · Y (0) то есть Y ( x ) = K ( x left ) · Yleft . Это можно делать следующим образом. Можно подставить Yright = K ( right left ) · Yleft в R · Yright = R и тогда можем записать: L · Yleft = L R · K(right left) · Yleft = R, откуда можно получить искомый вектор левого края Yleft и далее можно вычислять решения в произвольной точке по формуле: Y ( x ) = K ( x left ) · Yleft . Или иначе можем записать - не Yright = K ( right left ) · Yleft , а Yleft = K ( left right ) · Yright и можем тогда записать: L ·Yleft = L R ·Yright = R Yleft = K(left right) · Yright и подставив одно в другое можем записать: L · K(left right) · Yright = L R · Yright = R, откуда можно получить искомый вектор правого края Yright и далее можно вычислять решения в произвольной точке по формуле: Y ( x ) = K ( x right ) · Yright . Приблизительно 10 лет назад Алексеем Юрьевичем Виноградовым было сформулировано, что можно делать иначе, был предложен новый простенький в своей элегантности ( J ) метод, который излагается далее. Было показано, что можно не приводить краевую задачу к задаче Коши, а можно вычислять решение краевой задачи в произвольной рассматриваемой точке при помощи переноса в эту точку краевых условий . 2. Метод «переноса краевых условий в произвольную точку» А.Ю.Виноградова можно сформулировать следующим образом.

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

Поэтому можно не только традиционно интегрировать слева на право и обратно от края и до края, а можно как это предложено А.Ю.Виноградовым интегрировать от некоторой внутренней точки x интервала краевой задачи в разные стороны к краям ( lef и right ) интервала краевой задачи: Yleft = K(left x) · Y(x) Yright = K(right x) · Y(x) Это означает, что можно вычислить не одну матрицу Коши на всю длину участка рассматриваемой задачи K ( right left ), а можно вычислить независимо отдельно две матрицы Коши.

Вычисляться будут две матрицы Коши, но это не дольше, так как каждая из двух матриц вычисляется не на всём участке задачи, а на своём отдельном участке общей задачи: K ( left x ) и K ( right x ). Тогда при совместном рассмотрении краевых условий: L ·Yleft = L R ·Yright = R и матриц Коши K ( left x ) и K ( right x ) можно записать: L · K ( left x ) · Y ( x ) = L R · K(right x) · Y(x) = R Отсюда получаем систему обыкновенных линейных алгебраических уравнений с квадратной матрицей коэффициентов для вычисления вектора искомого решения Y ( x ) краевой задачи в любой точке x : | L · K ( left x ) | | L | |--------------------------| · Y ( x ) = |----| | R · K ( right x ) | | R | Для нахождения решения задачи в окрестности точки x (то есть при x + или x - ) то есть для нахождения решений задачи в окрестности найденного вектора Y ( x ) можно использовать формулы начальной задачи, то есть формулы задачи Коши: Y(x+ ) = K( (x+ ) x ) · Y(x) Y(x- ) = K( (x- ) x ) · Y(x) хотя это не столько формулы задачи Коши сколько вычислительная суть матрицы Коши. Было предложено для вычисления матрицы Коши использовать не методы численного интегрирования типа Рунге-Кутта для вычисления векторов составляющих матрицу Коши, а было предложено использовать теорию матриц и вычислять матрицу Коши как матричную экспоненту, которая вычисляется при помощи матричного ряда в виде K(x)=exp(Ax). Можно посмотреть формулу матричного ряда exp ( Ax ) в «Теории матриц» Гантмахера. Как потом было показано - применение матричной экспоненты exp ( Ax ) значительно ускоряет вычисления по сравнению с методами типа Рунге-Кутта в рамках метода. 3. Этот же предложенный метод как детально показано Алексеем Юрьевичем Виноградовым подходит и для успешного решения краевых задач с «жёсткими» дифференциальными уравнениями. При моделировании систем при помощи дифференциальных уравнений они иногда оказываются «жёсткими». Это задачи типа расчёта на прочность тонкостенных оболочек в ракето и самолёто-строении, в кораблестроении, в трубопроводах, баках и прочих конструкциях из тонкого металла, пластика или композиционного материала. Для решения таких краевых задач с «жёсткими» дифференциальными уравнениями обычно применяют специальные приёмы-методы.

Краевые задачи для «жёстких» дифференциальных уравнений решаются различными методами.

Сейчас наиболее распространён метод С.К.Годунова (появившийся в 1961 году в журнале «Успехи математических наук»). На применении метода С.К.Годунова к разным задачам в разных областях науки и техники было сделано множество кандидатских и докторских диссертаций. В методе С.К.Годунова (с 1961 года) применяется ортонормирование. В методе С.К.Годунова ортонормированию подвергается только половина Y 1( x ) Y 2( x ) Y 3( x ) Y 4( x ) из возможных для вычисления независимых векторов-столбцов Y 1( x ) Y 2( x ) Y 3( x ) Y 4( x ) Y 5( x ) Y 6( x ) Y 7( x ) Y 8( x ) , то есть векторов-функций-решений Y ( x ) , которые могли бы входить в матрицу Коши. Для того чтобы можно было вычислять только половину, то есть Y 1( x ) Y 2( x ) Y 3( x ) Y 4( x ) из возможных векторов-решений приходиться их подгонять под то, чтобы они при любом своем виде не противоречили краевым условиям на левом крае . Подогнать под левый край можно изменением формы записи дифференциальных уравнений так, чтобы условия на краях записывались очень несложно. Но даже при такой подгонке нельзя решать задачи с более-менее не примитивными условиями на краях потому, что не получается «угадать» так, чтобы искать только половину из возможных векторов-решений. Хотя можно преодолеть «угадывание» свойственное методу С.К.Годунова специальными математическими приёмами (которые были предложены тоже Алексеем Виноградовым в «Журнале вычислительной математики и математической физики» Академии наук), которые были сформулированы тоже приблизительно 10 лет назад чуть ранее появления метода «переноса краевых условий в рассматриваемую точку». В методе С.К.Годунова от левого края векторы-решения Y 1( x ) Y 2( x ) Y 3( x ) Y 4( x ) с константами-коэффициэнтами влияния Y 1( x )с1 Y 2( x )с2 Y 3( x )с3 Y 4( x )с4 переносятся с ортонормированием на правый край. Далее на правом крае вычисляются неизвестные коэффициенты (коэффициенты влияния) с1 с2 с3 с4 этих векторов в составе решения задачи. Далее в методе С.К.Годунова выполняется перенос векторов Y 1( x ) Y 2( x ) Y 3( x ) Y 4( x ) совместно с вычисленными коэффициентами их влияния с1 с2 с3 с4 с правого края на левый для последовательного вычисления решений задачи Y ( x ) во всех точках рассматриваемого участка с правого края по левый. В методе «переноса краевых условий в произвольную точку» А.Ю.Виноградова тоже может применяться ортонормирование. Но применяться может иначе - ортонормированию подвергаются переносимые краевые условия. Далее показано как применять построчное ортонормирование прямоугольных горизонтальных матриц краевых условий (комбинированных с матрицами Коши) метода А.Ю.Виноградова.

Прямоугольные горизонтальные матрицы L и R краевых условий переносятся матрицами Коши в рассматриваемую точку интервала интегрирования и построчному ортонормированию подвергаются матричные линейные алгебраические уравнения переносимых краевых условий. У ортонормируемых матричных уравнений – всегда прямоугольные горизонтальные матрицы коэффициентов.