Расщепления метод

- сеточный метод решения нестационарных задач со многими пространственными переменными, в к-ром переход от заданного временного слоя tn к новому слою [img: http://localhost:8080/file/041871-30.jpg] осуществляется за счет последовательного решения сеточных аналогов родственных нестационарных задач е меньшим числом пространственных переменных (см. [1] - [4]). Часто в этом классе методов могут быть найдены такие, что 1) весь переход от сеточного слоя в момент времени tn к новому сеточному слою [img: http://localhost:8080/file/041871-31.jpg] является достаточно простым и может быть осуществлен при затрате О(N)арифметич. действий, где N - число узлов пространственной сетки; 2) гарантируется абсолютная устойчивость метода; 3) гарантируется наличие приемлемой точности метода (наличие аппроксимации в том или ином смысле). Р. м. довольно широко применяются при практич. решении многомерных задач математич. физики, связанных, напр., с линейными и нелинейными системами параболического, гиперболического или смешанного типа (см. [1] - [8]). Обычно для задачи с рпространственными переменными переход от tn к [img: http://localhost:8080/file/041871-32.jpg] в Р. м. производится с использованием рвспомогательных (дробных) шагов: [img: http://localhost:8080/file/041871-33.jpg] (1) где [img: http://localhost:8080/file/041871-34.jpg] - матрица, соответствующая разностной аппроксимации нек-рого дифференциального оператора, содержащего производные только по xs (одномерного дифференциального оператора), а правые части (1) легко вычислимы. При соответствующей нумерации неизвестных, связанной с выбором направления xs, матрицы [img: http://localhost:8080/file/041871-35.jpg] становятся обычно диагональными и решение систем (1) при каждом s сводится к многократному решению одномерных разностных систем но направлению xs. Поэтому часто Р. м. наз. также или переменных направлений методом или дробных шагов методом. Одним из типичных примеров в случае уравнения [img: http://localhost:8080/file/041871-36.jpg] с начальным условием [img: http://localhost:8080/file/041871-37.jpg] и краевым условием [img: http://localhost:8080/file/041871-38.jpg], где [img: http://localhost:8080/file/041871-39.jpg], Г- граница W, может служить следующий метод, построенный на квадратной сетке с шагом h: [img: http://localhost:8080/file/041871-40.jpg] (2) где [img: http://localhost:8080/file/041871-41.jpg] [img: http://localhost:8080/file/041871-42.jpg] - простейшие разностные аппроксимации для Wn- совокуп- [img: http://localhost:8080/file/041871-43.jpg] ность внутренних узлов [img: http://localhost:8080/file/041871-44.jpg], [img: http://localhost:8080/file/041871-45.jpg] Имеются два альтернативных подхода к теории Р. м. В одном из них промежуточные шаги ни в чем существенном не отличаются от целых шагов, и сами разностные уравнения на дробных шагах и граничные условия для них, подобно методу (2), устроены одинаково, и можно ожидать, что [img: http://localhost:8080/file/041871-46.jpg] и [img: http://localhost:8080/file/041871-47.jpg] будут служить аппроксимациями для решения исходной задачи в моменты времени [img: http://localhost:8080/file/041871-48.jpg] и [img: http://localhost:8080/file/041871-49.jpg]. Этот подход основан на использовании понятия составной схемы и суммарной аппроксимации (см. [2]). Схемы такого типа часто наз. л о к а л ь н о о д н о м е р н ы м и с х е м а м и или а д д и т и в н ы м и с х е м а м и; их можно также трактовать как обычные разностные схемы для нек-рого уравнения с сильно осциллирующими по времени коэффициентами, решение к-рого должно быть близко к решению исходной задачи (см. [1]-[4]). Достоинства этого подхода в его простоте и общности, напр, обобщения метода (2) возможны и для случая криволинейных областей W и более общих задач. Точность же получаемых на этом пути методов обычно не очень высока. Известны и иногда успешно применяются варианты Р. м., в к-рых расщепление производится не по пространственным переменным, а по физич. процессам (см. [5]). Второй подход в плане анализа устойчивости и сходимости исключает какие-либо дробные шаги из рассмотрения. Сама разностная схема и аппроксимация трактуются традиционным образом. Необычность разностной схемы проявляется лишь в том, что на верхнем слое схемы появляется необычный разностный оператор. Напр., вместо метода (2) рассматривается метод [img: http://localhost:8080/file/041871-50.jpg] (3) где [img: http://localhost:8080/file/041871-51.jpg], Е - тождественный оператор. Такие операторы Аобычно наз. р а с щ е п л яю щ и м и с я или ф а к т о р и з о в а н н ы м и о п е р а т о р а м и. Дробные шаги связываются лишь с методом решения возникающих систем и для одной и той же схемы (3) могут быть введены различными способами, граничные условия для них должны выбираться в зависимости от этого. Сами схемы типа (3) можно трактовать как обычные схемы с весом для e-уравнения, напр., вида [img: http://localhost:8080/file/041871-52.jpg] решения к-рого отличаются от решения исходной задачи на О(e) (см. [4]). В случае области W, составленной из прямоугольников, матрицы возникающих систем в методах типа (3) уже не нредставимы в виде произведения "одномерных" матриц. Все же решения подобных систем могут быть найдены при затрате O(N). арифметич. действий (см. [4]), операторы подобного типа наз. р а с ш и р е н н о р а с щ е п л я ю щ и м и с я о п ер а т о р а м и. При исследовании устойчивости и сходимости схем с расщепляющимися и расширенно расщепляющимися операторами большую роль играет метод энергетич. неравенств (см. [2], [4], [6] - (8]).