Журнал вычислительной математики и математической физики, 2022, T. 62, № 7, стр. 1224-1232

Разностная схема метода декомпозиции начально-краевой задачи для сингулярно возмущенного уравнения переноса

Г. И. Шишкин 1*, Л. П. Шишкина 1

1 ИММ УрО РАН
620108 Екатеринбург, ул. С. Ковалевской, 16, Россия

* E-mail: shishkin@imm.uran.ru

Поступила в редакцию 17.12.2021
После доработки 17.12.2021
Принята к публикации 11.02.2022

Полный текст (PDF)

Аннотация

Рассматривается начально-краевая задача для сингулярно возмущенного уравнения переноса. Предлагается новый подход к построению разностной схемы, основанный на специальной декомпозиции решения в виде суммы регулярной и сингулярной компонент решения. Строится разностная схема метода декомпозиции решения, в котором регулярная и сингулярная компоненты решения рассматриваются на равномерных сетках, и устанавливается их $\varepsilon $-равномерная сходимость в равномерной норме с первым порядком скорости сходимости. По сеточным решениям компонент решения строится континуальное решение, аппроксимирующее решение начально-краевой задачи для сингулярно возмущенного уравнения переноса, и устанавливается его $\varepsilon $-равномерная сходимость в равномерной норме с первым порядком скорости сходимости. Предлагаемый подход позволяет в дальнейшем применять технику повышения скорости сходимости сеточных решений на вложенных сетках для построения разностных схем, сходящихся $\varepsilon $-равномерно со вторым порядком скорости сходимости и выше, для начально-краевой задачи для сингулярно возмущенного уравнения переноса. Библ. 11.

Ключевые слова: уравнение переноса, сингулярно возмущенная начально-краевая задача, пограничный слой, стандартная разностная схема, декомпозиция решения, равномерная сетка, $\varepsilon $-равномерная сходимость, равномерная норма, континуальная аппроксимация решения.

1. ВВЕДЕНИЕ

Задачи для регулярных уравнений переноса часто возникают в теоретических исследованиях и в приложениях; методы их решения рассматривались в работах [1]–[6] (см. также библиографию там). Начально-краевая задача для сингулярно возмущенного уравнения переноса с возмущающим (малым) параметром $\varepsilon $, $\varepsilon \in (0,1]$, при “конвективном члене” рассматривалась в [7], [8]; для этой задачи построена разностная схема на кусочно-равномерной сетке (такие сетки известны в литературе как “сетки Шишкина”), сгущающейся в окрестности пограничного слоя, которая сходится $\varepsilon $-равномерно в равномерной норме с первым порядком скорости сходимости. Однако использование таких сеток для построения разностных схем для сингулярно возмущенных задач с порядком скорости сходимости выше первого приводит к достаточно “громоздким” схемам (см., например, [9]).

Таким образом, появляется интерес к использованию более простых равномерных сеток для построения разностных схем для сингулярно возмущенного уравнения переноса, хотя известно, что непосредственное использование стандартных сеточных аппроксимаций сингулярно возмущенных задач (на равномерных сетках) не позволяет получить $\varepsilon $-равномерную сходимость решения.

В настоящей работе для начально-краевой задачи для сингулярно возмущенного уравнения переноса предлагается новый перспективный подход к построению разностной схемы, основанный на специальной декомпозиции решения как суммы регулярной и сингулярной компонент решения. При этом решение дифференциальной задачи рассматривается как сумма решений подзадач для регулярной и сингулярной компонент, рассматриваемых на своих равномерных сетках, каждая из которых сходится $\varepsilon $-равномерно в равномерной норме с первым порядком скорости сходимости. По сеточным решениям разностной схемы метода декомпозиции решения строится континуальная аппроксимация компонент решения и само континуальное решение начально-краевой задачи для сингулярно возмущенного уравнения переноса. Устанавливается $\varepsilon $-равномерная сходимость континуального решения в равномерной норме с первым порядком скорости сходимости.

Работа построена следующим образом. Постановка начально-краевой задачи для сингулярно возмущенного уравнения переноса и цель исследования приводятся в разд. 2. Стандартная разностная схема – схема на основе монотонной аппроксимации начально-краевой задачи – рассматривается в разд. 3. Априорные оценки решения и производных, используемые при построении и обосновании сходимости разностных схем, устанавливаются в разд. 4. В разд. 5 строится разностная схема метода декомпозиции решения, регулярная и сингулярная компоненты которого решаются на равномерных сетках, и устанавливается их $\varepsilon $-равномерная сходимость в равномерной норме с первым порядком скорости сходимости. В разд. 6 по сеточным решениям разностной схемы метода декомпозиции решения строится континуальная аппроксимация решения начально-краевой задачи для сингулярно возмущенного уравнения переноса и устанавливается $\varepsilon $-равномерная сходимость континуального решения с первым порядком скорости сходимости. Выводы приводятся в разд. 7.

2. ПОСТАНОВКА ЗАДАЧИ. ЦЕЛЬ ИССЛЕДОВАНИЯ

На множестве $\bar {G}$:

(2.1a)
$\bar {G} = G \cup S,\quad G = D \times (0,T],\quad S = {{S}_{0}} \cup {{S}^{l}},$
где $D = (0 < x \leqslant d]$, ${{S}_{0}}$ и ${{S}^{l}}$ – нижняя и боковая части границы $S$, и
(2.1б)
${{S}_{0}} = \{ (x,t):0 \leqslant x \leqslant d,\;t = 0\} ,\quad {{S}^{l}} = \{ (x,t):x = 0,\;0 < t \leqslant T\} ,$
причем $d \sim 1$, рассмотрим начально-краевую задачу для сингулярно возмущенного уравнения переноса (из [7], [8])
(2.2)
$\begin{gathered} L{\kern 1pt} u(x,t) = f(x,t),\quad (x,t) \in G; \\ u(x,t) = \varphi (x,t),\quad (x,t) \in S. \\ \end{gathered} $
Здесь
(2.3a)
$L = \varepsilon {\kern 1pt} b(x,t)\frac{\partial }{{\partial x}} + c(x,t) + p(x,t)\frac{\partial }{{\partial t}},\quad (x,t) \in G,$
функции $b(x,t)$, $c(x,t)$, $p(x,t)$, $f(x,t)$ предполагаются достаточно гладкими на $\bar {G}$, функция $\varphi (x,t)$, $(x,t) \in S$, – достаточно гладкая на множествах ${{S}_{0}}$ и $\bar {S}{{{\kern 1pt} }^{l}}$ и непрерывна на $S$, причем,
(2.3б)
$\begin{gathered} {{b}_{0}} \leqslant b(x,t) \leqslant {{b}^{0}},\quad 0 \leqslant c(x,t) \leqslant {{c}^{0}},\quad {{p}_{0}} \leqslant p(x,t) \leqslant {{p}^{0}}, \\ \left| {f(x,t)} \right| \leqslant M,\quad (x,t) \in \overline G ;\quad \left| {\varphi (x,t)} \right| \leqslant M,\quad (x,t) \in S;\quad {{b}_{0}},{{p}_{0}} > 0; \\ \end{gathered} $
параметр $\varepsilon $ принимает произвольные значения из полуинтервала $(0,1]$. (Через M (через m) обозначаем достаточно большие (малые) положительные постоянные, не зависящие от величины параметра ε. В случае сеточных задач эти постоянные не зависят и от шаблонов разностных схем.)

Считаем, что данные задачи (2.2), (2.1) на множестве ${{S}^{c}} = {{S}_{0}} \cap {{\bar {S}}^{l}}$ – угловой точке, удовлетворяют условиям согласования, обеспечивающим требуемую по построениям гладкость решения на $\bar {G}$.

При малых значениях параметра $\varepsilon $ в окрестности множества $S{{{\kern 1pt} }^{l}}$ появляется регулярный пограничный слой – узкая подобласть, примыкающая к множеству $S{{{\kern 1pt} }^{l}}$, ширины порядка $\mathcal{O}(\varepsilon )$, в которой решение задачи изменяется на конечную величину.

В [7], [8] для сингулярно возмущенной начально-краевой задачи для уравнения переноса (2.2), (2.1) с использованием кусочно-равномерной сетки (“сетки Шишкина”) построена монотонная разностная схема, сходящаяся $\varepsilon $-равномерно в равномерной норме со скоростью $\mathcal{O}({{N}^{{ - 1}}} + N_{0}^{{ - 1}})$, где $N$ и ${{N}_{0}}$ – число сеточных интервалов сеток по $x$ и $t$ соответственно. Однако использование таких сеток для повышения скорости сходимости численных решений выше первого порядка вызывает затруднения при построении разностных схем, сходящихся $\varepsilon $-равномерно.

В настоящей работе разрабатывается подход к построению разностной схемы, основанный на декомпозиции решения задачи, которая используется обычно при выводе априорных оценок (см., например, [9], [10]). Решение задачи (2.2), (2.1) представляется в виде следующей декомпозиции:

(2.4)
$u(x,t) = U(x,t) + V(x,t),\quad (x,t) \in \bar {G},$
где $U(x,t)$ и $V(x,t)$ – регулярная и сингулярная компоненты решения соответственно. Рассматриваются сеточная аппроксимация компонент на соответствующих равномерных сетках и построение континуальной аппроксимации решения рассматриваемой задачи.

Таким образом, цель настоящего исследования – для сингулярно возмущенной начально-краевой задачи для уравнения переноса (2.2), (2.1) разработать разностную схему метода декомпозиции решения, т.е. разностную схему на основе специальной декомпозиции решения в виде суммы регулярной и сингулярной компонент решения, рассматриваемых на соответствующих равномерных сетках, и построить континуальное решение, сходящееся к решению задачи для уравнения переноса (2.2), (2.1) $\varepsilon $-равномерно в равномерной норме с первым порядком скорости сходимости.

Такой подход позволит в дальнейшем применять технику повышения скорости сходимости сеточных решений на вложенных сетках (см., например, [9]) для построения разностных схем, сходящихся $\varepsilon $-равномерно со вторым порядком скорости сходимости и выше, для задач типа (2.2), (2.1).

3. СТАНДАРТНАЯ РАЗНОСТНАЯ СХЕМА

В этом разделе мы рассмотрим стандартную разностную схему, строящуюся на основе монотонной сеточной аппроксимации начально-краевой задачи (2.2), (2.1).

3.1. Приведем стандартную разностную схему для уравнения переноса.

На множестве $\bar {G}$ введем прямоугольную сетку

(3.1)
${{\bar {G}}_{h}} = \bar {\omega } \times {{\bar {\omega }}_{0}},$
где $\bar {\omega }$ и ${{\bar {\omega }}_{0}}$ – произвольные, вообще говоря, неравномерные сетки на отрезках $[0,d]$ и $[0,T]$ соответственно. Пусть ${{h}^{i}} = {{x}^{{i + 1}}} - {{x}^{i}}$, ${{x}^{i}},{{x}^{{i + 1}}} \in \bar {\omega }$, $h = {{\max }_{i}}{{h}^{i}}$, и ${{\tau }^{k}} = {{t}^{{k + 1}}} - {{t}^{k}}$, ${{t}^{k}},{{t}^{{k + 1}}} \in {{\bar {\omega }}_{0}}$, $\tau = {{\max }_{k}}{{\tau }^{k}}$. Предполагаем выполненным условие $h \leqslant M{\kern 1pt} {{N}^{{ - 1}}}$, $\tau \leqslant M{\kern 1pt} N_{0}^{{ - 1}}$, где $N + 1$ и ${{N}_{0}} + 1$ – число узлов сеток $\bar {\omega }$ и ${{\bar {\omega }}_{0}}$ соответственно.

Задачу (2.2), (2.1) аппроксимируем стандартной разностной схемой (см. [6])

(3.2a)
$\Lambda z(x,t) = f(x,t),\;\;(x,t) \in {{G}_{h}};\quad z(x,t) = \varphi (x,t),\quad (x,t) \in {{S}_{h}}.$
Здесь ${{G}_{h}} = G \cap {{\bar {G}}_{h}}$, ${{S}_{h}} = S \cap {{\bar {G}}_{h}}$,
(3.2б)
$\Lambda \equiv \varepsilon {\kern 1pt} b(x,t){\kern 1pt} {{\delta }_{{\bar {x}}}} + c(x,t) + p(x,t){\kern 1pt} {{\delta }_{{\bar {t}}}},$
${{\delta }_{{\bar {x}}}}{\kern 1pt} z(x,t)$ и ${{\delta }_{{\bar {t}}}}{\kern 1pt} z(x,t)$ – первые разностные производные (производные назад) по $x$ и $t$ соответственно.

Схема (3.2), (3.1) монотонна $\varepsilon $-равномерно (определение монотонности разностных схем см., например, в [6]). Для схемы (3.2), (3.1) справедлив сеточный принцип максимума.

Теорема 3.1. Пусть для стандартной разностной схемы (3.2), (3.1) выполняется условие $\Lambda z(x,t) \geqslant 0$, $(x,t) \in {{G}_{h}}$; $z(x,t) \geqslant 0$, $\;(x,t) \in {{S}_{h}}$. Тогда для функции $z(x,t)$ справедлива оценка $z(x,t) \geqslant 0$, $(x,t) \in {{\bar {G}}_{h}}$.

3.2. Для решения стандартной разностной схемы (3.2) на сетке ${{\bar {G}}_{{h{\kern 1pt} (3.1)}}}$, с использованием априорных оценок (4.10) и соотношения для сингулярной компоненты (4.11) получаем оценку

(3.3)
$\left| {u(x,t) - z(x,t)} \right| \leqslant M{\kern 1pt} \left[ {{{{\left( {\varepsilon + {{N}^{{ - 1}}}} \right)}}^{{ - 1}}}{\kern 1pt} {{N}^{{ - 1}}} + N_{0}^{{ - 1}}} \right],\quad (x,t) \in {{\bar {G}}_{{h(3.1)}}}.$

(Запись ${{D}_{{(i.j)}}}$ (${{L}_{{(i.j)}}}$, ${{m}_{{(i.j)}}}$, ${{M}_{{(i.j)}}}$, ${{D}_{{h(i.j)}}})$ означает, что эти множества (операторы, постоянные, сетки) введены в формуле (i, j).) Для решения стандартной разностной схемы (3.2) на сетке $\bar {G}_{h}^{u}$, равномерной по $x$ и $t$,

(3.4)
${{\bar {G}}_{h}} = \bar {G}_{h}^{u} \equiv \bar {\omega } \times {{\bar {\omega }}_{0}},$
имеем следующую оценку, подобную (3.3):
(3.5)
$\left| {u(x,t) - z(x,t)} \right| \leqslant M{\kern 1pt} \left[ {{{{\left( {\varepsilon + {{N}^{{ - 1}}}} \right)}}^{{ - 1}}}{\kern 1pt} {{N}^{{ - 1}}} + N_{0}^{{ - 1}}} \right],\quad (x,t) \in \bar {G}_{{h{\kern 1pt} (3.4)}}^{u},$
где ${{M}_{{h{\kern 1pt} (3.3)}}}$, ${{M}_{{h{\kern 1pt} (3.5)}}} = \mathcal{O}(1)$.

Таким образом, стандартные разностные схемы (3.2), (3.1) и (3.2), (3.4) сходятся при $N,{{N}_{0}} \to \infty $ при фиксированных значениях параметра $\varepsilon $ с первым порядком скорости сходимости.

Теорема 3.2. Пусть для данных начально-краевой задачи (2.2), (2.1) выполняются условия $(2.3)$, а для решения задачи – оценки теоремы 4.1 при $K = 2$. Тогда для решения стандартных разностных схем (3.2), (3.1) и (3.2), (3.4) справедливы оценки (3.3) и (3.5) соответственно.

Замечание 1. Схемы (3.2), (3.1) и (3.2), (3.4) сходятся при неулучшаемом условии

(3.6)
${{N}^{{ - 1}}} = o(\varepsilon ),\quad N_{0}^{{ - 1}} = o(1).$
В соответствии с оценкой (3.5) для сходимости стандартных разностных схем (3.2), (3.1) и (3.2), (3.4) требуется использовать сетки по $x$ с числом узлов, удовлетворяющим условию $N \gg {{\varepsilon }^{{ - 1}}}$, т.е. неограниченно растущим при $\varepsilon \to 0$. Таким образом, стандартные схемы (3.2), (3.1) и (3.2), (3.4) , сходящиеся при фиксированных значениях параметра $\varepsilon $, не сходятся $\varepsilon $-равномерно.

4. АПРИОРНЫЕ ОЦЕНКИ РЕШЕНИЯ И ПРОИЗВОДНЫХ

Приведем ряд априорных оценок решения задачи (2.2), (2.1), используемых при построении и обосновании разностных схем. Вывод оценок подобен выводу оценок производных регулярных и сингулярных компонент решения сингулярно возмущенных эллиптических и параболических уравнений из [9]–[11].

4.1. Для сингулярно возмущенной начально-краевой задачи (2.2), (2.1) справедлив принцип максимума, подобный принципу максимума для начально-краевых задач для сингулярно возмущенных параболических уравнений [9].

Теорема 4.1. Пусть для данных начально-краевой задачи для уравнения переноса (2.2), (2.1) выполняется условие

$Lu(x,t) \geqslant 0,\quad (x,t) \in G;\quad u(x,t) \geqslant 0,\quad (x,t) \in S.$
Тогда для функции $u(x,t)$ справедлива оценка $u(x,t) \geqslant 0$, $(x,t) \in \bar {G}$.

4.2. Приведем “стандартные” априорные оценки для решения сингулярно возмущенного уравнения переноса (2.2), (2.1). При не слишком малых значениях параметра $\varepsilon $ исследуемая задача подобна задаче для регулярного уравнения переноса, рассматриваемой в [1], [2].

1. Применяя технику мажорантных функций (подобную приведенной в [9], [10] для параболического уравнения конвекции-диффузии), находим оценку решения начально-краевой задачи (2.2), (2.1):

(4.1)
$\left| {u(x,t)} \right| \leqslant M,\quad (x,t) \in \bar {G}.$

2. При исследовании производных решения начально-краевой задачи считаем, что коэффициенты и правая часть уравнения являются достаточно гладкими на $\bar {G}$, начальная и граничная функции – достаточно гладкие на боковой и нижней частях границы $S$, производные $\frac{{{{\partial }^{k}}}}{{\partial {{x}^{k}}}}\varphi (x,0)$, $\frac{{{{\partial }^{{{{k}_{0}}}}}}}{{\partial {{t}^{{{{k}_{0}}}}}}}\varphi (0,t)$ ограничены при $k,{{k}_{0}} \leqslant 2$, причем в угловой точке $(x,t) = (0,0)$ выполняются условия согласования, обеспечивающие требуемую непрерывность решения начально-краевой задачи на $\bar {G}$ и его производных по $x$ и $t$ до второго порядка. Пусть выполняется условие

(4.2)
$u \in {{C}^{{k,{{k}_{0}}}}}(\bar {G}),\quad k + {{k}_{0}} \leqslant 2.$

3. При оценке производных решения задачи (2.2), (2.1) перейдем к переменным $\xi ,\;t$, где $\xi = {{\varepsilon }^{{ - 1}}}x$. В новых переменных имеем следующую начально-краевую задачу:

(4.3)
$\begin{gathered} \tilde {L}\tilde {u}(\xi ,t) \equiv \left\{ {\tilde {b}(\xi ,t)\frac{\partial }{{\partial \xi }} + \tilde {c}(\xi ,t) + \tilde {p}(\xi ,t)\frac{\partial }{{\partial t}}} \right\}\tilde {u}(\xi ,t) = \tilde {f}(\xi ,t),\quad (\xi ,t) \in \tilde {G}; \\ \tilde {u}(\xi ,t) = \tilde {\varphi }(\xi ,t),\quad (\xi ,t) \in \tilde {S}. \\ \end{gathered} $

Теперь уже для регулярной задачи (4.3) в переменных $\xi ,{\kern 1pt} t$ в случае условий согласования второго порядка в переменных $x,\;t$ на множестве ${{\tilde {S}}^{c}}$, где

(4.4)
$\begin{gathered} \frac{\partial }{{\partial t}}u(0,0) = \frac{d}{{dt}}{{\varphi }^{l}}(0),\quad \frac{{{{\partial }^{2}}}}{{\partial {{t}^{2}}}}u(0,0) = \frac{{{{d}^{2}}}}{{d{{t}^{2}}}}{{\varphi }^{l}}(0); \\ \frac{{{{\partial }^{k}}}}{{\partial {{t}^{k}}}}{{\varphi }^{l}}(0) = \mathop {\lim }\limits_{t \to 0} \frac{{{{\partial }^{k}}}}{{\partial {{t}^{k}}}}\varphi (0,t),\quad k = 1,2, \\ \end{gathered} $
обеспечивающих включение (4.2), находим оценку производных
(4.5)
$\left| {\frac{{{{\partial }^{{k + {{k}_{0}}}}}}}{{\partial {{\xi }^{k}}\partial {{t}^{{{{k}_{0}}}}}}}\tilde {u}(\xi ,t)} \right| \leqslant M,\quad (\xi ,t) \in \bar {\tilde {G}},\quad k + {{k}_{0}} \leqslant 2.$
В переменных $x,{\kern 1pt} t$ получаем оценку

(4.6)
$\left| {\frac{{{{\partial }^{{k + {{k}_{0}}}}}}}{{\partial {{x}^{k}}\partial {{t}^{{{{k}_{0}}}}}}}u(x,t)} \right| \leqslant M{{\varepsilon }^{{ - k}}},\quad (x,t) \in \bar {G},\quad k + {{k}_{0}} \leqslant 2.$

Справедлива следующая

Теорема 4.2. Пусть для данных начально-краевой задачи для уравнения переноса (2.2), (2.1) выполняется условие $b,c,p,f \in {{C}^{{k,{{k}_{0}}}}}(\bar {G})$, $k + {{k}_{0}} \leqslant 2$, и пусть на множестве ${{\tilde {S}}^{c}}$ выполняется условие согласования второго порядка (4.4). Тогда для решения начально-краевой задачи $u(x,{\kern 1pt} t)$ и его производных справедливы оценки (4.1), (4.6).

4.3. Приведем априорные оценки решений и производных для начально-краевой задачи для сингулярно возмущенного уравнения переноса (2.2), (2.1), используемые в разд. 5 при построении разностных схем метода декомпозиции решения и обосновании их сходимости.

1. Рассмотрим оценки регулярной и сингулярной компонент решения задачи (2.2), (2.1) из декомпозиции (2.4).

Регулярная компонента решения – функция $U(x,t)$ – определяется соотношением

(4.7)
$U(x,t) = {{U}^{e}}(x,t),\quad (x,t) \in \bar {G},$
где $U(x,t)$ есть сужение на $\bar {G}$ функции ${{U}^{e}}(x,t)$, $(x,t) \in {{\bar {G}}^{e}}$, являющейся решением следующей “расширенной” задачи:
(4.8a)
$\begin{gathered} {{L}^{e}}{{U}^{e}}(x,t) = {{f}^{e}}(x,t),\quad (x,t) \in {{G}^{e}}; \\ {{U}^{e}}(x,t) = {{\varphi }^{e}}(x,t),\quad (x,t) \in {{S}^{e}}. \\ \end{gathered} $
Здесь множество
(4.8б)
$\bar {G}{{{\kern 1pt} }^{e}} = [ - d{\kern 1pt} *,d] \times [0,T],\quad d{\kern 1pt} * > 0,$
есть продолжение $\bar {G}$ за боковую границу ${{S}^{l}}$, оператор ${{L}^{e}}$ и функции ${{f}^{e}}(x,t)$ и ${{\varphi }^{e}}(x,t)$ – расширения оператора $L$ на множество $\bar {G}{{{\kern 1pt} }^{e}}$ и функций $f(x,t)$ и $\varphi (x,t)$ на множества $\bar {G}{{{\kern 1pt} }^{e}}$ и ${{S}^{e}}$ соответственно с сохранением свойств данных задачи (2.2), (2.1). Предполагаем выполненным условие
$d{\kern 1pt} * \geqslant T{\kern 1pt} \mathop {\max }\limits_{{{{\bar {G}}}^{e}}} \{ {{b}^{e}}(x,t)[{{p}^{e}}{{]}^{{ - 1}}}(x,t)\} .$
Сингулярная компонента – функция $V(x,t)$ – решение следующей задачи:

(4.9)
$\begin{gathered} L{\kern 1pt} V(x,t) = 0,\quad (x,t) \in G; \\ V(x,t) = \varphi (x,t) - {{U}^{e}}(x,t),\quad (x,t) \in {{S}^{l}}; \\ V(x,t) = 0,\quad (x,t) \in {{S}_{0}}. \\ \end{gathered} $

2. Для компонент ${{U}^{e}}(x,t)$ и $V(x,t)$, $(x,t) \in \bar {G}$, справедливы оценки

(4.10a)
$\left| {\frac{{{{\partial }^{{k + {{k}_{0}}}}}}}{{\partial {{x}^{k}}\partial {{t}^{{{{k}_{0}}}}}}}{{U}^{e}}(x,t)} \right| \leqslant M,\quad (x,t) \in {{\bar {G}}^{e}};$
(4.10б)
$\left| {\frac{{{{\partial }^{{k + {{k}_{0}}}}}}}{{\partial {{x}^{k}}\partial {{t}^{{{{k}_{0}}}}}}}V(x,t)} \right| \leqslant M{{\varepsilon }^{{ - k}}},\quad (x,t) \in \bar {G},\quad k + {{k}_{0}} \leqslant 2,$
причем,
(4.11a)
$V(x,t) = 0,\quad (x,t) \in \bar {G}\quad {\text{при}}\quad x \geqslant {{M}_{1}}\varepsilon t,$
где ${{M}_{1}}$ – произвольная постоянная, удовлетворяющая условию

(4.11б)
${{M}_{1}} \geqslant {\text{ma}}{{{\text{x}}}_{{\bar {G}}}}{\kern 1pt} [{\kern 1pt} b(x,t){\kern 1pt} {{p}^{{ - 1}}}(x,t){\kern 1pt} ].$

Справедлива следующая

Теорема 4.3. Пусть для данных начально-краевой задачи для сингулярно возмущенного уравнения переноса (2.2), (2.1) выполняется условие $b,c,p,f \in {{C}^{{k,{{k}_{0}}}}}(\bar {G})$, $k + {{k}_{0}} \leqslant 2$, и пусть на множестве ${{\tilde {S}}^{c}}$ выполняется условие согласования второго порядка (4.4). Тогда для решения начально-краевой задачи и его компонент из представления (2.4) справедливы оценки (4.10), а также соотношение (4.11) для сингулярной компоненты $V(x,t)$.

5. РАЗНОСТНАЯ СХЕМА МЕТОДА ДЕКОМПОЗИЦИИ РЕШЕНИЯ

Построим разностную схему метода декомпозиции решения, рассматривая сеточные аппроксимации регулярной и сингулярной компонент решения из декомпозиции (2.4) на своих равномерных сетках.

Заметим, что в силу оценки (4.10) регулярная компонента $U(x,t)$ и ее производные по $x$ и $t$ ограничены $\varepsilon $-равномерно; производные по $t$ сингулярной компоненты $V(x,t)$ ограничены $\varepsilon $-равномерно, производные $k$-го порядка по $x$ растут как $\mathcal{O}\left( {{{\varepsilon }^{{ - k}}}} \right)$.

5.1. Построим сеточную аппроксимацию регулярной компоненты ${{U}^{e}}(x,t)$, $(x,t) \in {{\bar {G}}^{e}}$. Заметим, что функция $U(x,t)$, $(x,t) \in \bar {G}$, связана с функцией ${{U}^{e}}(x,t)$, $(x,t) \in {{\bar {G}}^{e}}$, соотношением (4.7).

На множестве $\bar {G}_{{(4.8)}}^{e}$ строим равномерную сетку

(5.1a)
$\bar {G}_{h}^{e} = {{\bar {\omega }}^{e}} \times {{\bar {\omega }}_{0}},$
где ${{\bar {\omega }}^{e}}$ и ${{\bar {\omega }}_{0}}$ – равномерные сетки по $x$ и $t$. Шаги этих сеток определяются соотношениями
(5.1б)
${{h}^{e}} = (d + d{\kern 1pt} *){\kern 1pt} {{N}^{{ - 1}}},\quad \tau = T{\kern 1pt} N_{0}^{{ - 1}},$
где $N$ и ${{N}_{0}}$ – число интервалов разбиения множеств $[ - d{\kern 1pt} *,d]$ и $[0,T]$ соответственно; $d{\kern 1pt} * = d_{{(4.8)}}^{*}$.

Задачу (4.8) аппроксимируем разностной схемой

(5.2)
${{\Lambda }^{e}}z_{U}^{e}(x,t) = {{f}^{e}}(x,t),\quad (x,t) \in G_{h}^{e};\quad z_{U}^{e}(x,t) = {{\varphi }^{e}}(x,t),\quad (x,t) \in S_{h}^{e},$
где оператор ${{\Lambda }^{e}}$ есть оператор ${{\Lambda }_{{(3.2)}}}$, но на множестве ${{\bar {G}}^{e}}$.

Разностная схема (5.2), (5.1) является монотонной. С использованием априорных оценок регулярной компоненты решения “расширенной задачи” (4.8) устанавливается оценка

(5.3)
$\left| {{{U}^{e}}(x,t) - z_{U}^{e}(x,t)} \right| \leqslant M{\kern 1pt} [{\kern 1pt} {{N}^{{ - 1}}} + N_{0}^{{ - 1}}{\kern 1pt} ],\quad (x,t) \in \bar {G}_{h}^{e}.$
Сеточное решение $z_{U}^{e}(x,t)$ сходится $\varepsilon $-равномерно с первым порядком скорости сходимости по $x$ и $t$.

Справедлива следующая

Теорема 5.1. Пусть для регулярной компоненты ${{U}^{e}}(x,t)$, $(x,t) \in {{\bar {G}}^{e}}$, из декомпозиции (2.4) решения начально-краевой задачи для сингулярно возмущенного уравнения переноса (2.2), (2.1) выполняется априорная оценка (4.10a). Тогда решение разностной схемы (5.2), (5.1) сходится $\varepsilon $-равномерно в равномерной норме с первым порядком скорости сходимости по $x$ и $t$ с оценкой (5.3).

5.2. Построим сеточную аппроксимацию сингулярной компоненты $V(x,t)$, $(x,t) \in \bar {G}$.

На множестве

(5.4)
${{\bar {G}}_{V}} = [0,{\kern 1pt} \min ({{M}_{1}}\varepsilon T,d)] \times [0,T]$
построим равномерную сетку
(5.5)
${{\bar {G}}_{{Vh}}} = \{ {{\bar {\omega }}_{V}} \times {{\bar {\omega }}_{0}}\} \cap {{\bar {G}}_{V}},$
где ${{\bar {\omega }}_{V}}$ и ${{\bar {\omega }}_{0}}$ – равномерные сетки на множествах $[0,\min ({{M}_{1}}\varepsilon T,d)]$ и $[0,T]$ соответственно с шагами ${{h}_{V}} = \min ({{M}_{1}}{\kern 1pt} \varepsilon {\kern 1pt} T,{\kern 1pt} d){\kern 1pt} {{N}^{{ - 1}}}$ и $\tau = T{\kern 1pt} N_{0}^{{ - 1}}$ по $x$ и $t$ соответственно; $\tau = {{\tau }_{{(5.1)}}}$, ${{\bar {\omega }}_{{0(5.5)}}} = {{\bar {\omega }}_{{0(5.1)}}}$.

Задачу (4.9) аппроксимируем разностной схемой

(5.6)
$\begin{gathered} \Lambda {\kern 1pt} {{z}_{V}}(x,t) = 0,\quad (x,t) \in {{G}_{{Vh}}}; \\ {{z}_{V}}(x,t) = {{\varphi }^{e}}(x,t) - z_{U}^{e}(x,t),\quad (x,t) \in {{S}_{{Vh}}}, \\ \end{gathered} $
где $z_{U}^{e}(x,t) = z_{{U(5.2)}}^{e}(x,t)$, ${{\bar {G}}_{{Vh}}} = {{G}_{{Vh}}} \cup {{S}_{{Vh}}}$.

Разностная схема (5.6), (5.5) является монотонной.

С использованием априорных оценок сингулярной компоненты решения задачи (4.9) устанавливается оценка

(5.7)
$\left| {V(x,t) - {{z}_{V}}(x,t)} \right| \leqslant M[{{N}^{{ - 1}}} + N_{0}^{{ - 1}}],\quad (x,t) \in {{\bar {G}}_{{Vh}}}.$
Сеточное решение ${{z}_{V}}(x,t)$ сходится $\varepsilon $-равномерно с первым порядком скорости сходимости по $x$ и $t$.

Справедлива следующая

Теорема 5.2. Пусть для сингулярной компоненты $V(x,t)$, $(x,t) \in {{\bar {G}}_{V}}$, из декомпозиции (2.4) решения начально-краевой задачи для сингулярно возмущенного уравнения переноса (2.2), (2.1) выполняется априорная оценка (4.10б). Тогда решение разностной схемы (5.6), (5.5) сходится $\varepsilon $-равномерно; для решения разностной схемы справедлива оценка (5.7).

6. КОНТИНУАЛЬНАЯ АППРОКСИМАЦИЯ РЕШЕНИЯ СИНГУЛЯРНО ВОЗМУЩЕННОГО УРАВНЕНИЯ ПЕРЕНОСА

По сеточным решениям компонент разностной схемы метода декомпозиции решения строим континуальную аппроксимацию компонент решения и самого решения начально-краевой задачи для сингулярно возмущенного уравнения переноса (2.2), (2.1) и устанавливаем $\varepsilon $-равномерную сходимость в равномерной норме с первым порядком скорости сходимости как компонент, так и самого решения.

6.1. Рассмотрим аппроксимацию регулярной компоненты $U(x,t)$, $(x,t) \in \bar {G}$, из декомпозиции.

Пусть $z_{U}^{e}(x,t)$, $(x,t) \in \bar {G}_{h}^{e}$, – решение разностной схемы метода декомпозиции решения (5.2), (5.1). По значениям функции $z_{U}^{e}(x,t)$ в узлах сетки

(6.1)
$\bar {G}_{h}^{U} = \bar {G}_{h}^{e} \cap \bar {G}$
строим интерполянт
(6.2)
$\tilde {U}(x,t),\quad (x,t) \in \bar {G}.$
В узлах элементарной ячейки, определяемой узлами
(6.3)
$({{x}_{i}},{\kern 1pt} {{t}_{j}}),\;({{x}_{{i + 1}}},{\kern 1pt} {{t}_{j}}),\;({{x}_{i}},{\kern 1pt} {{t}_{{j + 1}}}),\;({{x}_{{i + 1}}},{\kern 1pt} {{t}_{{j + 1}}}),$
полагаем
(6.4)
$\tilde {U}(x,t) = z_{U}^{e}(x,t).$
На элементарных прямоугольниках, определяемых узлами (6.3), интерполянт (6.2) строится линейным по $x$ и $t$.

Приведенную технику построения интерполянта $\tilde {U}(x,t)$, $(x,t) \in \bar {G}$, назовем континуальной аппроксимацией регулярной компоненты $U(x,t)$ решения начально-краевой задачи для сингулярно возмущенного уравнения переноса.

С использованием априорных оценок регулярной компоненты решения устанавливается следующая оценка для интерполянта $\tilde {U}(x,t)$, подобная оценке (5.3):

(6.5)
${\text{|}}U(x,t) - \tilde {U}(x,t){\text{|}} \leqslant M[{{N}^{{ - 1}}} + N_{0}^{{ - 1}}],\quad (x,t) \in \bar {G}.$

Таким образом, интерполянт $\tilde {U}(x,t)$, $(x,t) \in \bar {G}$, регулярной компоненты решения сходится $\varepsilon $-равномерно с первым порядком скорости сходимости по $x$ и $t$.

Справедлива следующая

Теорема 6.1. Пусть для регулярной компоненты ${{U}^{e}}(x,t)$, $(x,t) \in {{\bar {G}}^{e}}$, решения начально-краевой задачи для сингулярно возмущенного уравнения переноса (2.2), (2.1) выполняется априорная оценка (4.10a). Тогда континуальная аппроксимация регулярной компоненты решения – интерполянт $\tilde {U}(x,t)$, $(x,t) \in \bar {G}$, – сходится $\varepsilon $-равномерно в равномерной норме; для функции $\tilde {U}(x,t)$ справедлива оценка (6.5).

6.2. Рассмотрим аппроксимацию сингулярной компоненты $V(x,t)$, $(x,t) \in \bar {G}$.

Пусть ${{z}_{V}}(x,t)$, $(x,t) \in {{\bar {G}}_{{Vh}}}$, есть решение разностной схемы метода декомпозиции решения (5.6), (5.5).

По значениям функции ${{z}_{V}}(x,t)$ в узлах сетки ${{\bar {G}}_{{Vh}}}$ строим интерполянт

(6.6)
$\tilde {V}(x,t),\quad (x,t) \in {{\bar {G}}_{V}}.$
В узлах элементарной ячейки, определяемой узлами
(6.7)
$({{x}_{i}},{\kern 1pt} {{t}_{j}}),\;({{x}_{{i + 1}}},{\kern 1pt} {{t}_{j}}),\;({{x}_{i}},{\kern 1pt} {{t}_{{j + 1}}}),\;({{x}_{{i + 1}}},{\kern 1pt} {{t}_{{j + 1}}}),$
полагаем
(6.8)
$\tilde {V}(x,t) = z_{V}^{e}(x,t).$
На элементарных прямоугольниках, определяемых узлами (6.7), интерполянт (6.6) строится линейным по $x$ и $t$.

Приведенную технику построения интерполянта $\tilde {V}(x,t)$, $(x,t) \in {{\bar {G}}_{V}}$, назовем континуальной аппроксимацией сингулярной компоненты $V(x,t)$ решения начально-краевой задачи для сингулярно возмущенного уравнения переноса.

С использованием априорных оценок сингулярной компоненты решения устанавливается следующая оценка для интерполянта $\tilde {V}(x,t)$, подобная оценке (5.7):

(6.9)
${\text{|}}V(x,t) - \tilde {V}(x,t){\text{|}} \leqslant M[{{N}^{{ - 1}}} + N_{0}^{{ - 1}}],\quad (x,t) \in {{\bar {G}}_{V}}.$

Таким образом, интерполянт $\tilde {V}(x,t)$, $(x,t) \in {{\bar {G}}_{V}}$, сингулярной компоненты решения сходится $\varepsilon $-равномерно с первым порядком скорости сходимости по $x$ и $t$.

Справедлива следующая

Теорема 6.2. Пусть для сингулярной компоненты $V(x,t)$, $(x,t) \in {{\bar {G}}_{V}}$, решения начально-краевой задачи для сингулярно возмущенного уравнения переноса (2.2), (2.1) выполняется априорная оценка (4.10б). Тогда континуальная аппроксимация сингулярной компоненты решения – интерполянт $\tilde {V}(x,t)$, $(x,t) \in {{\bar {G}}_{V}}$, – сходится $\varepsilon $-равномерно в равномерной норме; для функции $\tilde {V}(x,t)$ справедлива оценка (6.9).

Далее будет удобно определить интерполянт $\tilde {V}(x,t)$ на всем множестве $\bar {G}$, полагая $\tilde {V}(x,t) = 0$ при $(x,t) \in \bar {G}{{\backslash }}{{\bar {G}}_{V}}$.

6.3. Построим континуальное решение начально-краевой задачи для сингулярно возмущенного уравнения переноса (2.2), (2.1) – функцию $\tilde {u}(x,t)$, $(x,t) \in \bar {G}$, – как сумму интерполянтов:

(6.10)
$\tilde {u}(x,t) = \tilde {U}(x,t) + \tilde {V}(x,t),\quad (x,t) \in \bar {G}.$

С учетом оценок (6.5) и (6.9) имеем

(6.11)
${\text{|}}u(x,t) - \tilde {u}(x,t){\text{|}} \leqslant M[{{N}^{{ - 1}}} + N_{0}^{{ - 1}}{\kern 1pt} ],\quad (x,t) \in \bar {G}.$
Континуальное решение $\tilde {u}(x,t)$, $(x,t) \in \bar {G}$, сходится к решению $u(x,t)$, $(x,t) \in \bar {G}$, начально-краевой задачи для сингулярно возмущенного уравнения переноса (2.2), (2.1) $\varepsilon $-равномерно с первым порядком скорости сходимости.

Теорема 6.3. Пусть для регулярной и сингулярной компонент из декомпозиции (2.4) решения начально-краевой задачи для сингулярно возмущенного уравнения переноса (2.2), (2.1) выполняются априорные оценки (4.10). Тогда континуальное решение $\tilde {u}(x,t)$, $(x,t) \in \bar {G}$, сходится к решению $u(x,t)$, $(x,t) \in \bar {G}$, начально-краевой задачи для сингулярно возмущенного уравнения переноса (2.2), (2.1) $\varepsilon $-равномерно в равномерной норме с оценкой (6.11).

7. ВЫВОДЫ

1. Рассмотрена начально-краевая задача для сингулярно возмущенного уравнения переноса с возмущающим параметром $\varepsilon $ при пространственной производной, где параметр $\varepsilon $ принимает произвольные значения из полуинтервала $(0,1]$. В этой задаче при малых значениях параметра $\varepsilon $ в окрестности боковой границы $S{{{\kern 1pt} }^{l}}$ (через которую характеристики входят в область) появляется пограничный слой.

2. Показано, что разностная схема на основе монотонных сеточных аппроксимаций задачи на сетке с произвольным распределеним узлов, а также на равномерной сетке, в случае ее сходимости в равномерной норме, не сходится $\varepsilon $-равномерно.

3. Разработана монотонная разностная схема метода декомпозиции решения, в которой регулярная и сингулярная компоненты решения вычисляются на своих равномерных сетках (см. построения в разд. 5 и утверждения теорем 5.1 и 5.2 об $\varepsilon $-равномерной сходимости сеточных решений регулярной и сингулярной компонент решения).

4. По сеточным аппроксимациям регулярной и сингулярной компонент решения в разд. 6 построена континуальная аппроксимация компонент решения и само континуальное решение начально-краевой задачи для сингулярно возмущенного уравнения переноса, и установлена $\varepsilon $-равномерная сходимость в равномерной норме с первым порядком скорости сходимости как компонент, так и самого континуального решения.

5. Разработанный в настоящей работе подход к построению разностной схемы метода декомпозиции решения позволит применить технику повышения скорости сходимости сеточных решений на вложенных сетках (см., например, [9]) для построения робастных разностных схем, сходящихся со вторым порядком скорости сходимости и выше, для начально-краевой задачи для сингулярно возмущенного уравнения переноса.

Список литературы

  1. Калиткин Н.Н. Численные методы. М.: Наука, 1978. 512 с.

  2. Калиткин Н.Н., Корякин П.В. Численные методы. Методы математической физики. М.: Изд-кий центр “Академия”, 2013. 304 с.

  3. Марчук Г.И. Методы вычислительной математики. М.: Наука, 1989. 608 с.

  4. Марчук Г.И., Шайдуров В.В. Повышение точности решений разностных схем. М.: Наука, 1979. 320 с.

  5. Самарский А.А. Введение в численные методы. М.: Наука, 1982. 272 с.

  6. Самарский А.А. Теория разностных схем. М.: Наука, 1989. 616 с.

  7. Шишкин Г.И. Разностная схема для начально-краевой задачи для сингулярно возмущенного уравнения переноса // Ж. вычисл. матем. и матем. физ. 2017. Т. 57. № 11. С. 1824–1830.

  8. Shishkina L.P., Shishkin G.I. Development and Numerical Study of Robust Difference Schemes for a Singularly Perturbed Transport Equation // Finite Difference Methods: Theory and Applications, FDM 2018, I. Dimov, I. Farago, and L. Vulkov (eds.), Lecture Notes in Computer Science (Springer, Cham, 2019). P. 476–483.

  9. Shishkin G.I., Shishkina L.P. Difference Methods for Singular Perturbation Problems. V. 140 of Chapman & Hall/CRC Monographs and Surveys in Pure and Applied Mathematics. Boca Raton: CRC Press, 2009. 408 p.

  10. Шишкин Г.И. Сеточные аппроксимации сингулярно возмущенных эллиптических и параболических уравнений. Екатеринбург: УрО РАН, 1992. 233 с.

  11. Miller J.J.H., O’Riordan E., Shishkin G.I. Fitted numerical methods for singular perturbation problems. Error estimates in the maximum norm for linear problems in one and two dimensions. Revised Ed. Singapore: World Sci., 2012. 176 p.

Дополнительные материалы отсутствуют.

Инструменты

Журнал вычислительной математики и математической физики