Математические заметки
RUS  ENG    ЖУРНАЛЫ   ПЕРСОНАЛИИ   ОРГАНИЗАЦИИ   КОНФЕРЕНЦИИ   СЕМИНАРЫ   ВИДЕОТЕКА   ПАКЕТ AMSBIB  
Общая информация
Последний выпуск
Скоро в журнале
Архив
Импакт-фактор
Правила для авторов
Лицензионный договор
Загрузить рукопись

Поиск публикаций
Поиск ссылок

RSS
Последний выпуск
Текущие выпуски
Архивные выпуски
Что такое RSS



Матем. заметки:
Год:
Том:
Выпуск:
Страница:
Найти






Персональный вход:
Логин:
Пароль:
Запомнить пароль
Войти
Забыли пароль?
Регистрация


Математические заметки, 2023, том 113, выпуск 1, страницы 90–108
DOI: https://doi.org/10.4213/mzm13573
(Mi mzm13573)
 

Эта публикация цитируется в 9 научных статьях (всего в 9 статьях)

Определение коэффициента теплопередачи в математических моделях тепломассопереноса

С. Г. Пятков, В. А. Баранчук

Югорский государственный университет, г. Ханты-Мансийск
Список литературы:
Аннотация: В работе рассматривается вопрос о корректности в пространствах Соболева обратных задач об определении коэффициента теплопередачи, входящего в граничное условие типа Робина для уравнений конвекции–диффузии. Доказана теорема существования и единственности решений.
Библиография: 37 названий.
Ключевые слова: обратная задача, коэффициент теплопередачи, параболическое уравнение, тепломассоперенос, диффузия.
Финансовая поддержка Номер гранта
Российский научный фонд 22-11-20031
Работа выполнена при финансовой поддержке Российского научного фонда и правительства Ханты-Мансийского автономного округа-ЮГРЫ (грант № 22-11-20031).
Поступило: 03.05.2022
Англоязычная версия:
Mathematical Notes, 2023, Volume 113, Issue 1, Pages 93–108
DOI: https://doi.org/10.1134/S0001434623010108
Реферативные базы данных:
Тип публикации: Статья
УДК: 517.95

1. Введение

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

$$ \begin{equation} Mu=u_t-Lu=f(t,x),\qquad (t,x)\in Q=(0,T)\times G, \end{equation} \tag{1.1} $$
где $Lu=\sum_{i,j=1}^na_{ij}(t,x)u_{x_i x_j} +\sum_{i=1}^n a_i(t,x)u_{x_i}+a_0(t,x)u$, $G\in\mathbb R^n$ – ограниченная область с границей $\Gamma$. Уравнение (1.1) дополняется начально-краевыми условиями:
$$ \begin{equation} Bu|_S=g(t,x)\quad (S=(0,T)\times\Gamma),\qquad u|_{t=0}=u_0(x), \end{equation} \tag{1.2} $$
где $Bu=\sum_{i,j=1}^n a_{ij}(t,x)u_{x_i}\nu_j +(\beta(t,x)+\beta_0(t,x))u$, ($\nu$ – внешняя единичная нормаль к $\Gamma$) и условиями переопределения
$$ \begin{equation} u(t,b_i)=\psi_i(t),\qquad i=1,2,\dots,r, \end{equation} \tag{1.3} $$
где $b_i\in\Gamma$, $\{b_i\}_{i=1}^r$ – некоторый набор точек. Задача состоит в нахождении решения уравнения (1.1), удовлетворяющего условиям (1.2), (1.3), и неизвестной функции $\beta(t,x)=\sum_{j=1}^r\beta_i(t)\Phi_i(t,x)$. где функции $\Phi_i$ заданы, а функции $\beta_i(t)$ считаются неизвестными.

Обратные задачи о нахождения неизвестных граничных режимов, в частности, задачи конвективного теплообмена, являются классическими. Они возникают в самых различных задачах математической физики: описание различных процессов тепломассопереноса, диффузии, фильтрации, экологии и т.п. (см. [1]–[4]).

В настоящее время имеется большое количество работ, посвященных численному решению задач (1.1)(1.3) в различных постановках, возникающих в приложениях точки $\{b_i\}$ в (1.3) могут быть как внутренними [5]–[13], так и граничными точками [14]–[16] области $G$ в последних работах постановки совпадают с вышеприведенной. В стационарном случае такая постановка рассматривалась в [17]. В работе [16] рассматривается параболическая система и ищутся постоянные коэффициенты теплопередачи (получена теорема единственности решений и описан численный метод). Коэффициент теплопередачи, зависящий от времени, по значениям решения в некотором наборе внутренних точек численно определяется в [7]–[9], [11]–[13]. В работах [5], [6] рассматривается система $2$-х одномерных параболических уравнений и задача численного определения постоянных коэффициентов теплопередачи (граничные условия матричные) и некоторых других параметров по данным замеров температуры внутри области. В работах [7], [18] коэффициент теплопередачи ищется зависящим от $x$, а данные – данные Дирихле на части боковой поверхности цилиндра. По данным такого же типа коэффициент теплопроводности, зависящий от всех переменных, численно определяется в [19]. Имеется ряд работ, посвященных определению коэффициента теплопередачи в нелинейном граничном условии вида

$$ \begin{equation*} \frac{\partial u}{\partial N}+\rho(t)\varphi(u)=g, \end{equation*} \notag $$
где функция $\rho$ считается неизвестной (см. [20]). Отметим работы (см. библиографию в [21]–[23]), где восстанавливается функция вида $\varphi(t,x,u)$ (иногда не зависящая от независимых переменных) в граничном условии Робина вида
$$ \begin{equation*} \frac{\partial u}{\partial N}+\varphi(t,x,u)=g \end{equation*} \notag $$
или близком условии. В этих работах, используются интегральные условия переопределения различного вида и в некоторых случаях получены теоремы существования и единственности решений таких задач локально по времени.

Большое количество работ посвящено также близкой задаче об определении потоков через границу или ее часть, с использованием тех же условий переопределения (1.3), которая возникает при линеаризации задачи (см. [4], [24]–[26] и библиографию в этих работах). Основной метод построения приближенного решения – сведение задачи к задаче управления и минимизация соответствующего квадратичного функционала.

Теоретических результатов, посвященных задаче (1.1)(1.3), крайне мало. Единственная работа, посвященная задаче (1.1)(1.3) в многомерном случае, есть работа [27] (см. также [28]), где в случае $Mu=u_t-\Delta u$ и $r=1$ была показана теорема существования и единственности классических решений задачи об определении потока и теорема единственности в задаче об определении коэффициента теплопередачи. Теорем существования решений задачи (1.1)(1.3) по всей видимости нет. Доказательство единственности в [27] основано на сведении задачи к интегральному уравнению Абеля, за счет выделения главной части функции Грина. По-видимому, перенести результаты с помощью этого метода на случай общих параболических уравнений и систем вряд ли возможно. В данной работе мы опишем другой подход к исследованию задачи (1.1)(1.3) об определении коэффициента теплопередачи $g(t,x)$, основанный на сведении задачи к интегральной уравнению Вольтерры второго рода; в этом случае решение находится методом последовательных приближений. На этой основе можно построить и новый метод численного решения задачи. В пространствах Соболева мы получаем теорему существования и единственности решений.

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

Пусть $E$ – банахово пространство. Обозначения для пространств Лебега $L_p(G;E)$, Соболева $W_p^s(G;E)$, $W_p^s(Q;E)$ и Гёльдера $C^\alpha(\overline G;E)$ $(\alpha\geqslant 0)$ стандартные (см. [29]–[31]). Если $E=\mathbb R$ или $E=\mathbb R^n$, то символ $E$ не используем и обозначения упрощаются, пишем просто $W_p^s(G)$ и т.д. Все рассматриваемые пространства и коэффициенты уравнения (1.1) мы считаем вещественными. Под нормой вектора понимаем сумму норм координат. Для данного интервала $J=(0,T)$ положим

$$ \begin{equation*} W_p^{s,r}(Q)=W_p^s(J;L_p(G))\cap L_p(J;W_p^r(G)). \end{equation*} \notag $$
Соответственно,
$$ \begin{equation*} W_p^{s,r}(S)=W_p^s(J;L_p(\Gamma))\cap L_p(J;W_p^r(\Gamma)). \end{equation*} \notag $$
Пусть
$$ \begin{equation*} (u,v)=\int_Gu(x)v(x)\,dx. \end{equation*} \notag $$

Говорим, что граница $\Gamma$ данной области $G$ принадлежит классу $C^s$, $s\geqslant 1$ (см. определение в [32; гл. 1]), если для любой точки $x_0\in\Gamma$ найдутся окрестность $Y_{x_0}$ (координатная окрестность) этой точки, и система координат $y$ (локальная система координат), полученная с помощью поворота и переноса начала координат из исходной, такая, что ось $y_n$ направлена по внутренней нормали в $\Gamma$ в точке $x_0$ и уравнение части границы $Y_{x_0}\cap\Gamma$ имеет вид $y_n=\gamma(y')$, $\gamma(0)=0$, $|y'|<\delta$, $y'=(y_1,\dots,y_{n-1})$, причем $\gamma\in C^s(\overline{B'_\delta})$ ($B'_\delta=\{y':|y'|<\delta\}$) и

$$ \begin{equation*} \begin{aligned} \, G\cap Y_{x_0}&=\{y:|y'|<\delta,\,0<y_n-\gamma(y')<\delta_1\}, \\ (\mathbb R^n\setminus G)\cap Y_{x_0}&=\{y:|y'|<\delta,-\delta_1<y_n-\gamma(y')<0\}. \end{aligned} \end{equation*} \notag $$
Числа $\delta$, $\delta_1$ для области $G$ фиксированы, причем без ограничения общности считаем, что $\delta_1>(M+1)\delta$, где $M$ постоянная Липшица функции $\gamma$. Обозначим, такой параметр $\delta$ через $\delta_\Gamma$ (он определен неоднозначно). Пусть $\Gamma_{x_0}=Y_{x_0}\cap\Gamma$.

Фиксируем параметр $\delta_\Gamma$. Пусть $B_\delta(b_i)$ – шар радиуса $\delta$ с центром в точке $b_i$. Параметр $\delta\in(0,\delta_\Gamma)$ назовем допустимым, если $\overline{B_\delta(b_i)}\cap\overline{B_\delta(b_j)}=\varnothing$ для $i\ne j$, $i,j=1,2,\dots,r$. Введем обозначения:

$$ \begin{equation*} \begin{gathered} \, Q^\tau=(0,\tau)\times G,\qquad G_\delta=\cup_i B_\delta(b_i)\cap G, \\ \Gamma_\delta=\Gamma\cap G_\delta,\qquad S_\delta=(0,T)\times\Gamma_\delta,\qquad S^\tau=(0,\tau)\times\Gamma,\quad i\leqslant r. \end{gathered} \end{equation*} \notag $$

Далее, мы считаем, что допустимый параметр $\delta>0$ зафиксирован. Рассматривая задачу (1)–(3), предполагаем, что

$$ \begin{equation} \Gamma\in C^2,\qquad \Gamma_\delta\in C^3. \end{equation} \tag{2.1} $$

Мы будем использовать в пространстве $W_p^s(0,\beta; E)$ ($s\in (0,1)$, $\beta>0$, $E$ – банахово пространство) норму

$$ \begin{equation*} \|q(t)\|_{W_p^s(0,\beta;E)} =(\|q\|_{L_p(0,\beta;E)}^p+\langle q\rangle_{s,\beta}^p)^{1/p},\qquad \langle q\rangle_{s,\beta}^p =\int_0^\beta\int_0^\beta\frac{\|q(t_1)-q(t_2)\|_E^p}{|t_1-t_2|^{1+sp}}\,dt_1\,dt_2. \end{equation*} \notag $$
Если $E=\mathbb R$, то мы получим обычное пространство $W_p^s(0,\beta)$. При $s\in(0,1)$ положим
$$ \begin{equation*} \widetilde W_p^s(0,\beta;E) =\{q\in W_p^s(0,\beta;E):t^{-s}q(t)\in L_p(0,\beta;E)\}. \end{equation*} \notag $$
Наделим это пространство нормой
$$ \begin{equation*} \|q(t)\|_{\widetilde W_p^s(0,\beta;E)}^p =\biggl\|\frac{q}{t^s}\biggr\|_{L_p(0,\beta;E)}^p+\langle q\rangle_{s,\beta}^p. \end{equation*} \notag $$
Если $s>1/p$ и $q\in\widetilde W_p^s(0,\beta;E)$, то $q(0)=0$ и эта норма и обычная норма $\|\cdot\|_{W_p^s(\alpha,\beta;E)}$ для функций $q(t)$ таких, что $q(0)=0$ эквивалентны (см. [30; п. 3.2.6, лемма 1]). Пространства $\widetilde W_p^s(0,\beta;L_p(G))$ и $\widetilde W_p^{s,2s}(Q^\beta)=\widetilde W_p^s(0,\beta;L_p(G)) \cap L_p(0,\beta;W_p^{2s}(G))$ при $s\ne 1/p$ состоят из функций $v(t,x)$ из $W_p^s(0,\beta;L_p(G))$ и $W_p^{s,2s}(Q^\beta)$, соответственно, таких, что $v(0,x)=0$ при $s>1/p$. Нормы ${\|\cdot\|}_{\widetilde W_p^{s,2s}(Q^\beta)}$, ${\|\cdot\|}_{\widetilde W_p^s(0,\beta;L_p(G))}$ определяются естественным образом, например,
$$ \begin{equation*} \|u\|_{\widetilde W_p^{s,2s}(Q^{\beta})} =(\|u\|_{\widetilde W_p^s(0,\beta;L_p(G))}^p +\|u\|_{L_p(0,{\beta};W_p^{2s}(G))}^p)^{1/p}. \end{equation*} \notag $$
Аналогично определяем пространства $\widetilde W_p^s(0,\beta;L_p(\Gamma))$, $\widetilde W_p^{s,2s}(S^\beta)$. Далее, мы считаем, что допустимый параметр $p>n+2$ зафиксирован. Следующая лемма известна (см. [33; леммы 1–4]).

Лемма 1. Существует постоянная $C$, не зависящая от $\tau\in(0,T]$, такая, что

$$ \begin{equation*} \|v\|_{\widetilde W_p^{s_1,2s_1}(S^\tau)} \leqslant C\|v\|_{W_p^{1,2}(Q^\tau)},\qquad \biggl\|\frac{\partial v}{\partial\nu}\biggr\|_{\widetilde W_p^{s_0,2s_0}(S^\tau)} \leqslant C\|v\|_{W_p^{1,2}(Q^\tau)}, \end{equation*} \notag $$
для всех $v\in W_p^{1,2}(Q^\tau)$ таких, что $v(x,0)=0$. Здесь $s_1=1-1/(2p)$ и $s_0=1/2-1/(2p)$.

Лемма 2. Пусть $s\in((n+2)/2p,1)$. Тогда, если $q\in\widetilde W_p^{s,2s}(Q^\tau)$ и $v\in W_p^{s,2s}(Q^\tau)$, то $qv\in\widetilde W_p^{s,2s}(Q^\tau)$ и справедлива оценка

$$ \begin{equation*} \|qv\|_{\widetilde W_p^{s,2s}(Q^\tau)} \leqslant c_0\|q\|_{\widetilde W_p^{s,2s}(Q^\tau)} (\|v\|_{W_p^{s,2s}(Q^\tau)}+\|v\|_{L_\infty(Q^\tau)}). \end{equation*} \notag $$
Если $v\in W_p^{s,2s}(Q)$, то последнее неравенство можно переписать в виде
$$ \begin{equation*} \|qv\|_{\widetilde W_p^{s,2s}(Q^\tau)} \leqslant c_1\|q\|_{\widetilde W_p^{s,2s}(Q^\tau)}\|v\|_{W_p^{s,2s}(Q)}, \end{equation*} \notag $$
а если $v\in\widetilde W_p^{s,2s}(Q^\tau)$, то в виде
$$ \begin{equation*} \|qv\|_{\widetilde W_p^{s,2s}(Q^\tau)} \leqslant c_2\|q\|_{\widetilde W_p^{s,2s}(Q^\tau)} \|v\|_{\widetilde W_p^{s,2s}(Q^\tau)}. \end{equation*} \notag $$
где постоянные $c_i$, $i=0,1,2$, не зависят от $q$, $v$ и $\tau\in(0,T]$. Если функция $v(t,x)$ строго отделена от нуля в $Q^\tau$, т.е. $\delta_0=\inf_{(t,x)\in Q^\tau}|v(t,x)|>0$, $q\in\widetilde W_p^{s,2s}(Q^\tau)$ и $v\in W_p^{s,2s}(Q^\tau)$, то $q/v\in\widetilde W_p^{s,2s}(Q^\tau)$ и
$$ \begin{equation*} \begin{aligned} \, \biggl\|\frac{q}{v}\biggr\|_{\widetilde W_p^s(0,\tau)} &\leqslant c_0\|q\|_{\widetilde W_p^{s,2s}(Q^\tau)} (\|v\|_{W_p^{s,2s}(Q^\tau)}+\|v\|_{L_\infty(Q^\tau)}), \\ \biggl\|\frac{q}{v}\biggr\|_{\widetilde W_p^{s,2s}(Q^\tau)} &\leqslant c_0\|q\|_{\widetilde W_p^{s,2s}(Q^\tau)}\|v\|_{W_p^{s,2s}(Q)}, \end{aligned} \end{equation*} \notag $$
где постоянная $c_0$ не зависит от функции $q$ и $\tau$. Множество $Q^\tau$ в этих утверждениях может быть заменено на $S^\tau$ (при этом считаем, что $s\in((n+1)/2p,1))$. В случае если одна из функций зависит только от одной переменной $t$, норма в пространствах $\widetilde W_p^{s,2s}(Q^\tau)$ или $W_p^{s,2s}(Q^\tau)$ в этих неравенствах заменяется на норму в $\widetilde W_p^s(0,\tau)$.

Доказательство. Доказательство основано на определении нормы и фактически оно доказано в доказательстве леммы 1 в [35]. Поэтому мы его опустим.

Замечание 1. Условие $s\in((n+2)/2p,1)$ гарантирует включение $W_p^{s.2s}(Q)\subset C(\overline Q)$ (см. § 6.3 и теорему 1 в разделе замечания (с. 424) в [36]).

Лемма 3. Пусть $G\subset\mathbb R^n$ – ограниченная область с границей класса $C^m$ или $G=\mathbb R^n$, $E$ – банахово пространство. Тогда

$$ \begin{equation*} \begin{gathered} \, (W_p^{s_1},W_p^{s_2})_{\theta,q}=B_{p,q}^s,\qquad (B_{p,q_0}^{s_1},B_{p,q_1}^{s_2})_{\theta,q}=B_{p,q}^s, \\ \frac{1-\theta}{q_0}+\frac{\theta}{q_1}=\frac{1}{q}\,,\qquad (1-\theta)s_1+\theta s_2=s, \\ B_{p,q}^s\subset B_{p,q_1}^s,\qquad q_1>q,\qquad B_{p_1,q}^{s_1}\subset B_{p,q}^s, \\ s_1-\frac{n}{p_1}=s-\frac{n}{p}\,,\qquad B_{p,p}^s=W_p^s,\qquad s\ne[s], \end{gathered} \end{equation*} \notag $$
где $q_i,p,q\in[1,\infty]$, $s_i,s\in[-m,m]$, $i=1,2$, $s_1\ne s_2$,
$$ \begin{equation*} W_p^s\subset B_{\infty,\infty}^{s-n/p} \subset C^{s-n/p}(\overline G;E),\qquad 1\leqslant p<\infty,\qquad s>\frac{n}{p}\,. \end{equation*} \notag $$

Здесь для сокращения обозначений под символами $W_p^s$, $B_{q,p}^s$ и т.д. понимаем пространства Соболева и Бесова функций, определенных на области $G$ со значениями в банаховом пространстве $E$, т.е. пространства $W_p^s(G;E)$, $B_{q,p}^s(G;E)$ и т.д.

Для областей с бесконечно дифференцируемой границей утверждение леммы содержится в [31] (см. § 2, 3, 4, соотношения (3.1)(3.9), следствие 4.3). Однако, как известно (см. замечение в конце работы [31]) утверждения могут быть доказаны и для областей с конечной гладкостью границы, для тех областей, заданные функции в которых допускают продолжение на все $\mathbb R^n$ с сохранением классов. В данном случае можно использовать стандартный метод продолжения Хестенса (см. утверждение и доказательство теоремы 4.2.2 в [30]).

Оператор $L$ считается эллиптическим, т.е. для некоторой постоянной $\delta_0>0$ выполнено неравенство

$$ \begin{equation*} \sum_{i,j=1}^na_{ij}\xi_i\xi_j\geqslant\delta_0|\xi|^2\qquad \forall\,\xi\in\mathbb R^n\quad \forall\,(t,x)\in Q. \end{equation*} \notag $$

Приведем условия на исходные данные. Считаем, что выполнены условия

$$ \begin{equation} a_i\in L_p(Q),\qquad a_{kl}\in C(\overline Q),\qquad \beta_0,a_{kl}|_\Gamma\in W^{s_0,2s_0}_p(S), \end{equation} \tag{2.2} $$
$$ \begin{equation} a_i\in L_\infty(0,T;W_p^1(G_\delta)),\qquad a_{kl}\in L_\infty(0,T;W_\infty^1(G_\delta)), \end{equation} \tag{2.3} $$
$$ \begin{equation} \beta\in W^{s_0,2s_0}_p(S),\qquad g(0,x)=B(x,0,\partial_x) u_0|_\Gamma, \end{equation} \tag{2.4} $$
где $i=0,1\dots,n$, $k,l=1,\dots,n$. Построим функции $\varphi_i(x)\in C_0^\infty(\mathbb R^n)$ такие, что $\varphi_i(x)=1$ в $B_{\delta/2}(b_i)$ и $\varphi_i(x)=0$ в $\mathbb R^n\setminus B_{3\delta/4}(b_i)$; положим
$$ \begin{equation*} \varphi(x)=\sum_{i=1}^r\varphi_i(x). \end{equation*} \notag $$

Мы используем выпрямление границы. Это преобразование

$$ \begin{equation*} z_n=y_n-\gamma(y'),\qquad z'=y', \end{equation*} \notag $$
где $y$ – локальная система координат в точке $b_i$. При выполнении условия (2.1) оно и обратное к нему $y_n=z_n+\gamma(z')$, $y'=z'$ принадлежат классу $C^3$ (т.е. $y=y(z)\in C^3(\overline{B'_\delta})$). Тоже самое утверждение имеет место и для преобразований $x=x(y(z))=x^i(z)$.

Пусть $U=\{z:|z'|<\delta,\,0<z_n<\delta_1\}$. Положим $Q_0^\tau=(0,\tau)\times U$, $Q_0=(0,T)\times U$ и $S_0^\tau=(0,\tau)\times B'_\delta$, $S_0=(0,T)\times B'_\delta$.

Считаем, что

$$ \begin{equation} u_0(x)\in W_p^{2-2/p}(G),\qquad f\in L_p(Q),\qquad g\in W_p^{s_0,2s_0}(S), \end{equation} \tag{2.5} $$
Пусть $Y_{b_i}$ – координатная окрестность точки $b_i\in\Gamma$. Выпрямим границу и перейдем к системе координат $z=(z',z_n)$. Мы также предполагаем, что
$$ \begin{equation} \nabla_{z'}\beta(x^i(z',0) \in W_p^{s_0,2s_0}(S_0), \end{equation} \tag{2.6} $$
$$ \begin{equation} \nabla_{z'}\varphi_ig(t,x^i(z',0)) \in W_p^{s_0,2s_0}(S_0),\qquad \nabla_{z'}\varphi_if(t,x^i(z))\in L_p(Q_0),\qquad i=1,2,\dots,r, \end{equation} \tag{2.7} $$
и что
$$ \begin{equation} \begin{aligned} \, \nabla_{z'}\varphi_iu_0(x^i(z)) &\in W_p^{2-2/p}(U), \nonumber \\ \nabla_{z'}a_{kl}(t,x^i(z',0)), \nabla_{z'}\beta_0(t,x^i(z',0)) &\in W_p^{s_0,2s_0}(S_0),\qquad k,l=1,2\dots,n,\quad i\leqslant r. \end{aligned} \end{equation} \tag{2.8} $$
Отметим, что условия (2.6)(2.8) не зависят от введенной локальной системы координат $y$ и системы координат $z$. Например, последнее включение в условии (2.8) эквивалентно тому, что для любого гладкого (класса $C^2$) касательного к $\Gamma$ векторного поля $\tau=(\tau_1(x),\dots,\tau_n(x))$
$$ \begin{equation*} \frac{\partial a_{ij}}{\partial\tau} \in W_p^{s_0,2s_0}((0,T)\times\Gamma_\delta),\qquad i,j=1,2\dots,n. \end{equation*} \notag $$
Последнее эквивалентно тому, что $\partial a_{ij}/\partial\tau\in W_p^{s_0,2s_0}((0,T)\times(B_\delta(b_i)\cap\Gamma))$ для всех $i$. В свою очередь $\partial/\partial z_i=\partial/\partial\tau$, $i\leqslant n-1$, для некоторого такого векторного поля $\tau$. Приведем вспомогательный результат. Условия (2.6)(2.8) могут быть переписаны и в терминах переменной $x$, но удобнее оставить их в вышеприведенной форме.

Теорема 1. Пусть выполнены условия (2.2), (2.4), (2.5) и $\Gamma\in C^2$. Тогда существует единственное решение задачи (1.1), (1.2) такое, что $u\in W_p^{1,2}(Q)$, причем справедлива оценка

$$ \begin{equation} \|u\|_{W_p^{1,2}(Q)}\leqslant C_0(\|u_0\|_{W_p^{2-2/p}(G)} +\|f\|_{L_p(Q)}+\|g\|_{W_p^{s_0,2s_0}(S)}). \end{equation} \tag{2.9} $$
Если $u_0\equiv 0$, то для любого $\tau\in(0,T]$ существует единственное решение $u\in W_p^{1,2}(Q^\tau)$ задачи (1.1), (1.2), удовлетворяющее оценке
$$ \begin{equation} \|u\|_{W_p^{1,2}(Q^\tau)} \leqslant C_1(\|f\|_{L_p(Q^\tau)}+\|g\|_{\widetilde W_p^{s_0,2s_0}(S^\tau)}), \end{equation} \tag{2.10} $$
где постоянная $C_1$ не зависит от $\tau$.

Доказательство. Вначале считаем, что $\tau=T$. Утверждение о разрешимости задачи (1.1), (1.2) из класса $u\in W_p^{1,2}(Q)$ вытекает из теоремы 2.1 [29]. В лучае произвольного $\tau\in(0,T]$ доказательство осуществляется по той же схеме, что и доказательство теоремы 2 в [33]. Поэтому мы его не приводим.

В формулировке следующей теоремы в каждой точке $b_i$ и соответствующей окрестности $Y_{b_i}$ мы будем использовать преобразование выпрямления границы и локальную систему координат $z$.

Теорема 2. Пусть выполнены условия (2.1)(2.8). Тогда на промежутке $(0,\tau)$ существует единственное решение $u$ задачи (1.1), (1.2) такое, что $u\in W_p^{1,2}(Q^\tau)$, причем $\nabla_{z'}\varphi_i u(x^i(z))\in W_p^{1,2}(Q_0^\tau)$, $i=1,\dots,r$. Если $u_0\equiv 0$, то решение удовлетворяет оценке

$$ \begin{equation} \begin{aligned} \, \sum_{i=1}^r\|\nabla_{z'}\varphi_iu(t,x^i(z))\|_{W_p^{1,2}(Q_0^\tau)} \leqslant C_1\biggl(&\|g\|_{\widetilde W_p^{s_0,2s_0}(S^\tau)}+\|f\|_{L_p(Q^\tau)} +\sum_{i=1}^r(\|\nabla_{z'}\varphi_if\|_{L_p(Q_0^\tau)} \nonumber \\ &\qquad{} +\|\nabla_{z'}\varphi_ig(x^i(z',0))\|_{\widetilde W_p^{s_0,2s_0}(S_0^\tau)})\biggr), \end{aligned} \end{equation} \tag{2.11} $$
где постоянная $C_1$ не зависит от $\tau\in(0,T]$ и $f$, $g$, $u_0$.

Доказательство. Доказательство использует метод конечных разностей (см. лемму 4.6 гл. 2 в [3] и теорему 4, п. 4, § 3, гл. 3 в [34]). Рассуждения совпадают с теми которые приведены в доказательстве теоремы 4, п. 3, § 2, гл. 4 в [34]. Поэтому мы остановимся только на основных моментах.

Возьмем точку $b_i\in\Gamma$. Умножая уравнение на $\varphi_i$, для $v=\varphi_iu$, имеем

$$ \begin{equation*} Mv=v_t-Lv=\varphi_if+[\varphi_i,L]u=\widetilde f,\qquad v|_{t=0}=\varphi_iu_0(x),\qquad Bv|_S=\varphi_ig-[\varphi_i,B]u, \end{equation*} \notag $$
где
$$ \begin{equation*} \begin{aligned} \, [\varphi_i,L]u&=\varphi_iLu-L(\varphi_iu) =-2\sum_{l,k=1}^na_{lk}u_{x_k}\varphi_{ix_l} -\sum_{l,k=1}^na_{lk}u\varphi_{ix_{l}x_k} -\sum_{k=1}^na_k\varphi_{ix_k}u, \\ [\varphi_i,B]u&=-\sum_{k,l=1}^na_{kl}\nu_k\varphi_{ix_l}u. \end{aligned} \end{equation*} \notag $$
Перейдем к локальной системе координат $y$ и обозначим полученные коэффициенты через $\widetilde a_{ij}$, $\widetilde a_i$. Оператор $B$ будет иметь тот же вид (изменятся фактические координаты вектора $\nu$). Выпрямим границу преобразованием $z_n=y_n-\gamma(y')$, $z'=y'$ и перейдем к новым координатам $z$ в уравнении. Получим задачу
$$ \begin{equation} \begin{aligned} \, Mv &=v_t-\widetilde Lv=\widetilde f,\qquad v|_{t=0}=v_0(z)=\varphi_i(x^i(z))u_0(x^i(z)), \\ \widetilde Bv|_{z_n=0} &=\varphi_ig(t,x^i(z',0))-[\varphi_i,B]u(t,x^i(z',0)) =\widetilde g(t,z'), \end{aligned} \end{equation} \tag{2.12} $$
где $\widetilde L$, $\widetilde B$ – операторы $L$, $B$, записанные в системе координат $z$. Обозначим коэффициенты $\widetilde L$ через $c_{kl}$, $c_k$, оператор $\widetilde B$ запишется в виде
$$ \begin{equation*} \widetilde Bu=\sum_{k=1}^n\gamma_k(t,z')u_{z_k}+(\beta+\beta_0)u,\qquad \text{где}\quad \gamma_n=\sqrt{1+|\nabla\gamma|^2} \sum_{k,l=1}^n\widetilde a_{k,l}\nu_k\nu_l|_{z_n=0}. \end{equation*} \notag $$
Уравнение рассматривается для $z\in U$. Пусть $\Delta_jv(z)=(v(z+e_j\eta)-v(z))/\eta$ ($e_j$ – $j$-й координатный вектор), где $|\eta|<\delta/8$ и $j\leqslant n-1$. Тогда функция $w=\Delta_j v$ есть решение задачи
$$ \begin{equation} \begin{gathered} \, w_t-\widetilde L(t,z,D)w =-[\widetilde L,\Delta_j]v+\Delta_j\widetilde f=\widetilde f_0, \\ \widetilde Bw|_{z_n=0} =[\widetilde B,\Delta_j]v+\Delta_j\widetilde g=\widetilde g_0,\qquad w|_{t=0}=\Delta_jv_0=\widetilde u_{01}, \end{gathered} \end{equation} \tag{2.13} $$
где
$$ \begin{equation*} \begin{aligned} \, [\widetilde L,\Delta_j]v &=-\sum_{k,l=1}^n\Delta_jc_{kl}(t,z)v_{z_{k}z_{l}}(t,z+e_j \eta) \\ &\qquad -\sum_{k=1}^n\Delta_jc_{k}(t,z)v_{z_{k}}(t,z+e_j \eta)-\Delta_jc_0(t,z)v(t,z+e_j \eta), \\ [\widetilde{B},\Delta_j]v &=-\sum_{k=1}^n\Delta_j\gamma_{k}v_{z_{k}}(t,z'+e_j\eta,0)+ \Delta_j(\beta+\beta_0)(t,z'+e_j\eta,0). \end{aligned} \end{equation*} \notag $$
Эти равенства получаются, если мы применим оператор $\Delta_j$ в (2.12). Вернемся к переменным $x$ и продолжим все функции в (2.13) нулем вне $Y_{b_i}\cap G$. Тогда функция $w\in W_p^{1,2}(Q)$ есть решение задачи (1.1), (1.2) с некоторыми новыми правыми частями в граничном условии и уравнении, т.е.
$$ \begin{equation} Mw=w_t-Lw=\widetilde f_0,\quad (t,x)\in Q,\qquad w|_{t=0}=\widetilde u_{01},\quad Bw|_S=\widetilde g_0. \end{equation} \tag{2.14} $$
Ссылаясь на теорему 1, заключаем, что справедлива оценка
$$ \begin{equation} \|w\|_{W_p^{1,2}(Q^\tau)} \leqslant C_1(\|\widetilde f_0\|_{L_p(Q^\tau)} +\|\widetilde g_0\|_{\widetilde W_p^{s_0,2s_0}(S^\tau)} +\|\widetilde u_{01}\|_{W_p^{2-2/p}(G)}), \end{equation} \tag{2.15} $$
В случае $u_0=0$ дополнительно имеем, что
$$ \begin{equation} \|w\|_{W_p^{1,2}(Q^\tau)} \leqslant C_2(\|\widetilde f_0\|_{L_p(Q^\tau)} +\|\widetilde g_0\|_{\widetilde W_p^{s_0,2s_0}(S^\tau)}), \end{equation} \tag{2.16} $$
где постоянная $C_2$ уже не зависит от параметра $\tau$. Далее необходимо показать, что нормы правых частей оцениваются постоянными, не зависящими от параметра $\eta\in(0,\delta/8)$. Оценки более или менее стандартные. Мы используем условия на данные, теоремы вложения и леммы 1, 2, а также представления вида
$$ \begin{equation*} \Delta_jv_0=\frac{1}{\eta} \int_{z_j}^{z_j+\eta}v_{0z_j}(z_1,\dots,z_{j-1},\xi,z_{j+1},\dots,z_n)\,d\xi =\int_0^{1}v_{0z_j}(z+e_j\tau\eta)\,d\tau. \end{equation*} \notag $$

Полученные оценки и лемма 4.11 в [32] или теорема 4, п. 4, § 3, гл. 3 в [34]), гарантируют утверждение теоремы.

3. Основные результаты

Пусть $p\in(n+2,\infty)$. Приведем условия на данные. Считаем, что функции $\Phi_i(t,x)$ при всех $i$ обладают свойствами

$$ \begin{equation} \Phi_i\in W_p^{s_0,2s_0}(S),\quad \nabla_{z'}\Phi_i(t,x^j(z',0))\in W_p^{s_0,2s_0}(S_0),\qquad i,j=1,2,\dots,r. \end{equation} \tag{3.1} $$
Пусть $\Phi(t)$ – матрица с элементами $\phi_{ij}=\Phi_j(t,b_i)$, $i,j=1,2,\dots,r$. В силу теорем вложения $\Phi_i(t,b_j)\in C^{1/2-(n+2)/2p}([0,T])$ (см. § 6.3 и теорему 1 (раздел замечания, с. 424) в [36]). Дополнительные условия на данные имеют вид
$$ \begin{equation} |\psi_i(t)|\geqslant\delta_1,\qquad |{\det\Phi}|\geqslant\delta_1>0\quad \forall\,t\in[0,T],\qquad \psi_i\in W_p^{s_1}(0,T),\qquad u_0(b_i)=\psi_i(0), \end{equation} \tag{3.2} $$
где $\delta_1$ – некоторая положительная постоянная и $i=1,2,\dots,r$.

Однако это не все условия, гарантирующие разрешимость задачи. Возьмем первое из равенств (1.2) в точке $(0,b_j)$. Имеем

$$ \begin{equation} Bu_0(b_j)=\frac{\partial u_0(b_j)}{\partial N}+\beta(0,b_j)u_0(b_j) =g(0,b_j)-\beta_0(0,b_j)u_0(b_i),\qquad j=1,\dots,r. \end{equation} \tag{3.3} $$
Из этой системы в силу (3.2) однозначно определяются величины $\beta_i(0)$ как решение системы
$$ \begin{equation} \sum_{i=1}^r\beta_i(0)\Phi_i(0,b_j) =\frac{1}{u_0(b_j)}\biggl(g(0,b_j)-\frac{\partial u_0(b_j)}{\partial N} -\beta_0(0,b_j)u_0(b_i)\biggr),\qquad j=1,\dots,r. \end{equation} \tag{3.4} $$
Тогда, если решение задачи (1.1)(1.3) существует, то выполнено равенство
$$ \begin{equation} \frac{\partial u_0(x)}{\partial N}+(\beta(0,x)+\beta_0(0,x))u_0(x)=g(0,x)\qquad \forall\,x\in\Gamma, \end{equation} \tag{3.5} $$
где постоянные $\beta_j(0)$ – решение системы (3.4). Это равенство – дополнительное условие согласования данных. Положим $\beta^0=\sum_{i=1}^r\beta_i(0)\Phi_i(t,x)$. Считая, что условия теорем 1, 2 выполнены, построим решение $w_0$ уравнения (1.1), удовлетворяющее условиям
$$ \begin{equation} w_0(0,x)=u_0(x), \end{equation} \tag{3.6} $$
$$ \begin{equation} \frac{\partial w_0(t,x)}{\partial N}+(\beta^0+\beta_0)w_0(t,x)=g,\qquad (t,x)\in S. \end{equation} \tag{3.7} $$
Отметим, что в силу условия (3.1) и леммы 2, $\beta^0\in W_p^{s_0,2s_0}(S)$, $\nabla_{z'}\beta^0(x^{j}(z',0))\in W_p^{s_0,2s_0}(S_0)$ для всех $j$ и таким образом, условия теорем 1, 2 будут выполнены. Сделаем замену $u=v+w_0$ в (1.1)(1.3). Обозначим $\alpha=\beta-\beta^0$, $\alpha_0=\beta^0+\beta_0$. Функция $v$ есть решение обратной задачи
$$ \begin{equation} Mv=v_t-Lv=0,\qquad (x,t)\in Q=G\times(0,T), \end{equation} \tag{3.8} $$
$$ \begin{equation} v|_{t=0}=0,\qquad \frac{\partial v}{\partial N}+\alpha_0v=-\alpha v-\alpha w_0, \end{equation} \tag{3.9} $$
$$ \begin{equation} v(t,b_i)=\widetilde\psi_i(t)=\psi_i(t)-w_0(t,b_i). \end{equation} \tag{3.10} $$
Имеем, что $\alpha=\sum_{i=1}^r\alpha_i\Phi_i(t,x)$, $\alpha_i=\beta_i(t)-\beta_i(0)$. Положим $\vec\alpha=(\alpha_1,\dots,\alpha_r)$. Сформулируем основной результат.

Теорема 3. Пусть выполнены условия (2.1)(2.3), (2.5), (2.7), (2.8), (3.1), (3.2), (3.5). Тогда на некотором промежутке $[0,\tau_0]$ существует единственное решение задачи (1.1)(1.3) такое, что $u\in W_p^{1,2}(Q^\tau)$, $\beta_i(t)\in W_p^{s_0}(0,\tau_0)$, $i=1,2,\dots,r$, причем $\nabla_{z'}\varphi_i u(x^i(z))\in W_p^{1,2}(Q_0^\tau)$, $i=1,\dots,r$.

Доказательство. Достаточно доказать утверждение для вспомогательной задачи (3.8)(3.10). Фиксируем $R_0>0$ (эту величину мы определим позже) и предположим, что

$$ \begin{equation*} \vec\alpha\in B_{R_0}=\{\vec\alpha\in\widetilde W_p^{s_0}(0,\tau): \|\vec\alpha\|_{\widetilde W_p^{s_0}(0,\tau)}\leqslant R_0\}. \end{equation*} \notag $$
Фиксируя $\vec\alpha\in B_{R_0}$ и решая задачу (3.8), (3.9), мы тем самым построим отображение $\vec\alpha\to v(\vec\alpha)$. Кроме этого отображения, нам понадобится еще одно отображение. Пусть $v$ – решение задачи (3.8), (3.9) из класса, описанного в формулировке теоремы. Фиксируя $i$ и умножая уравнение (3.8) на $\varphi_i$, имеем
$$ \begin{equation} Mw=w_t-Lw=[\varphi_i,L]v=f_0,\qquad w=\varphi_iv,\quad w\bigr\vert_{t=0}=0, \end{equation} \tag{3.11} $$
где
$$ \begin{equation*} [\varphi_i,L]v=\varphi_iLv-L(\varphi_iv) =-2\sum_{l,k=1}^na_{lk}v_{x_k}\varphi_{ix_l} -\sum_{l,k=1}^na_{lk}\varphi_{ix_lx_k}v -\sum_{k=1}^na_k\varphi_{ix_k}v. \end{equation*} \notag $$
Сделав замену переменных $x=x^i(z)$, перепишем (3.11) в виде
$$ \begin{equation} w_t-c_{nn}(t,z)w_{z_nz_n}=\sum_{l+k<2n}c_{kl}w_{z_kz_l} +\sum_{k=1}^nc_kw_{z_k}+c_0w+f_0=f_1,\qquad z\in U. \end{equation} \tag{3.12} $$
Отметим, что коэффициент $c_{nn}>0$ для всех $t$, $z$. В силу свойств решения $v$ и условий на коэффициенты, имеем что $f_1\in L_p(Q_0^\tau)$, $\nabla_{z'}f_1\in L_p(Q_0^\tau)$ и, более того, $f_1(t,z',z_n)\in C^\alpha(\overline{B'_\delta};L_p((0,\tau)\times(0,\delta_1)))$ с $\alpha\leqslant 1-(n-1)/p$ (см. лемму 3), после может быть изменения на множестве меры ноль. Рассмотрим уравнения
$$ \begin{equation} w_{it}(t,z_n)-c_{nn}(t,0,z_n)w_{iz_nz_n}=f_1(t,0,z_n)\qquad i=1,2,\dots,r,\qquad z_n\in (0,\delta_1), \end{equation} \tag{3.13} $$
Дополним уравнение (3.13) начальными и краевыми условиями
$$ \begin{equation} w_i(0,z_n)=0,\quad w_i|_{z_n=0}=\widetilde\psi_i(t),\quad w_i|_{z_n=\delta_1}=0,\qquad i\leqslant r. \end{equation} \tag{3.14} $$
Пусть $v(\vec\alpha)$ – решение задачи (3.8), (3.9); построим функции $w_i$ как решение задачи (3.13), (3.14). Очевидно, что если $v(\vec\alpha)$ есть решение задачи (3.8)(3.10), то имеем $w_i(t,z_n)=\varphi_iv(t,x^i(0,z_n))$. Таким образом, каждому $\vec\alpha$ отвечает функция $v$ и набор функций $w_i$, $i=1,2,\dots,r$. Перепишем граничное условие (3.9) в окрестности $Y_{b_i}$ в виде
$$ \begin{equation*} \sum_{j=1}^n\gamma_j(t,z')v_{z_j}(t,x^i(z',0)) +(\alpha_0v+\alpha v+\alpha w_0)(t,x^i(z',0))=0. \end{equation*} \notag $$
Отметим, что если функции $v$, $\vec\alpha$ есть решение обратной задачи, то $\omega_{iz_n}(t,0)=v_{z_n}(t,x^i(t,0,0))$ при $i\leqslant r$. Полагая $z'=0$ и используя (3.10), запишем равенства
$$ \begin{equation} \gamma_n(t,0)\omega_{iz_n}(t,0) +\sum_{j=1}^{n-1}\gamma_j(t,0)v_{z_j}(t,x^i(z',0)) +\alpha_0(t,b_i)\widetilde{\psi}_i(t)+\alpha(t,b_i)\psi_i(t)=0, \end{equation} \tag{3.15} $$
которые также можно переписать в виде
$$ \begin{equation} \alpha(t,b_i)=\frac{1}{\psi_i}\biggl(-\gamma_n(t,0)\omega_{iz_n}(t,0) -\sum_{j=1}^{n-1}\gamma_j(t,0)v_{z_j}(t,x^i(z',0)) -\alpha_0\widetilde\psi_i(t)\biggr)=F_i, \end{equation} \tag{3.16} $$
где $i=1,2,\dots,r$. Отметим, что в силу условия (3.2) знаменатели в этих равенствах строго отделены от нуля. Это и есть искомая система уравнений для нахождения координат вектора $\vec\alpha$. Она также может быть переписана в виде
$$ \begin{equation} \vec\alpha=\Phi^{-1}\vec F(\vec\alpha)=R(\vec\alpha), \end{equation} \tag{3.17} $$
где координаты вектора $F$ определены равенствами (3.16). Отметим, что лемма 2 гарантирует оценку
$$ \begin{equation} \|\Phi^{-1}\vec F\|_{\widetilde W_p^{s_0}(0,\tau)} \leqslant c\|\vec F\|_{\widetilde W_p^{s_0}(0,\tau)}. \end{equation} \tag{3.18} $$
Покажем, что оператор $R(\vec\alpha)$ является сжимающим в некотором шаре $B_{R_0}$ и переводит его в себя. Пусть $\vec\alpha\in B_{R_0}$ и $v$ – соответствующие решения задачи (3.8), (3.9). Из теоремы 1 и леммы 2 имеем оценку
$$ \begin{equation} \begin{aligned} \, \|v\|_{W_p^{1,2}(Q^\tau)} &\leqslant c_0\|\alpha (v+w_0)\|_{\widetilde W_p^{s_0,2s_0}(S^\tau)} \nonumber \\ &\leqslant c_1\|\alpha\|_{\widetilde W_p^{s_0,2s_0}(S^\tau)} (\|v\|_{\widetilde W_p^{s_0,2s_0}(S^\tau)} +\|w_0\|_{W_p^{s_0,2s_0}(S)}), \end{aligned} \end{equation} \tag{3.19} $$
где постоянные $c_0$, $c_1$ не зависит от $\tau$. В силу леммы 2 и условий (3.1) выполнено $\alpha\in\widetilde W_p^{s_0,2s_0}(S^\tau)$, $\nabla_{z'}\alpha(t,x^i(z',0))\in\widetilde W_p^{s_0,2s_0}(S_0^\tau)$, $i=1,\dots,r$, и имеем оценки
$$ \begin{equation} \|\alpha\|_{\widetilde W_p^{s_0,2s_0}(S^\tau)} \leqslant c_2\sum_{i=1}^r\|\alpha_i\|_{\widetilde W_p^{s_0}(0,\tau)} \|\Phi_i\|_{W_p^{s_0,2s_0}(S)} \leqslant c_3\|\vec\alpha\|_{\widetilde W_p^{s_0}(0,\tau)}\leqslant c_3R_0, \end{equation} \tag{3.20} $$
$$ \begin{equation} \sum_{i=1}^r\|\nabla_{z'}\alpha(t,x^i(z',0))\|_{\widetilde W_p^{s_0,2s_0}(S_0^\tau)} \leqslant c_4\|\vec\alpha\|_{\widetilde W_p^{s_0}(0,\tau)}\leqslant c_4R_0, \end{equation} \tag{3.21} $$
где постоянные $c_3$, $c_4$ не зависит от $\tau$. В норму $\|v\|_{\widetilde W_p^{s_0,2s_0}(S^\tau)}$ входят два слагаемых $\|v\|_{L_p(\Gamma;\widetilde W_p^{s_0}(0,\tau))}$ и $\|v\|_{L_p(0,\tau;\widetilde W_p^{2s_0}(\Gamma))}$ Оценим первое слагаемое. Имеем
$$ \begin{equation} \|v\|_{L_p(\Gamma;\widetilde W_p^{s_0}(0,\tau))} \leqslant c\tau^{1/2}\|v\|_{W_p^{1,2}(Q^\tau)}. \end{equation} \tag{3.22} $$
Действительно, по определения нормы, очевидно, что $\|v\|_{\widetilde W_p^{s_0}(0,\tau)} \leqslant\tau^{1/2}\|v\|_{\widetilde W_p^{s_1}(0,\tau)}$. Отсюда вытекает, что (лемма 1) $\|v\|_{L_p(\Gamma;\widetilde W_p^{s_0}(0,\tau))} \leqslant C\tau^{1/2}\|v\|_{W_p^{1,2}(Q^\tau)}$ и таким образом неравенство (3.22) имеет место. Далее имеем
$$ \begin{equation} \begin{aligned} \, \|v\|_{L_p(0,\tau; W_p^{1-1/p}(\Gamma))} &\leqslant C_1\|v\|_{L_p(0,\tau;W_p^1(G))} \nonumber \\ &\leqslant c_5\|v\|_{L_p(Q^\tau)}^{1/2}\|v\|_{W_p^{1,2}(Q^\tau)}^{1/2} \leqslant c_6\tau^{1/2}\|v\|_{W_p^{1,2}(Q^\tau)}. \end{aligned} \end{equation} \tag{3.23} $$
Здесь мы воспользовались теоремой о следах (см. [30; теорема 4.7.1]) и соответствующей оценкой $\|u\|_{W_p^{1-1/p}(\Gamma)}\leqslant c_1\|u\|_{W_p^1(G)}$, интерполяционным неравенством
$$ \begin{equation*} \|u\|_{W_p^1(G)}\leqslant C\|u\|_{L_p(G)}^{1/2}\|u\|_{W_p^2(G)}^{1/2} \end{equation*} \notag $$
(следствие теоремы 2 пункта 4.3.1 в [30]) и неравенством
$$ \begin{equation*} \|u\|_{L_p(0,\phi)}\leqslant\phi\|u_t\|_{L_p(0,\phi)}\qquad (u(0)=0), \end{equation*} \notag $$
вытекающим из формулы Ньютона–Лейбница. Здесь все постоянные не зависят от $\tau$. Из неравенств (3.19)(3.23) вытекает неравенство
$$ \begin{equation} \|v\|_{W_p^{1,2}(Q^\tau)}\leqslant c_7R_0 \tau^{1/2}\|v\|_{W_p^{1,2}(Q^\tau)} +c_8R_0\|w_0\|_{W_p^{s_0,2s_0}(S)}. \end{equation} \tag{3.24} $$
Выбрав $\tau_1$ такое, что $c_7R_0\tau_1^{1/2}\leqslant 1/2$, из (3.24) получим при $\tau\leqslant\tau_1$, что
$$ \begin{equation} \|v\|_{W_p^{1,2}(Q^\tau)} \leqslant 2c_8R_0\|w_0\|_{W_p^{s_0,2s_0}(S)}. \end{equation} \tag{3.25} $$
Здесь все постоянные не зависят от $\tau$. Сейчас мы запишем оценки из теоремы 2. Как вытекает из этой теоремы, имеет место оценка
$$ \begin{equation*} \begin{aligned} \, &\sum_{i=1}^r\|\nabla_{z'}\varphi_iv(t,x^i(z))\|_{W_p^{1,2}(Q_0^\tau)} \\ &\qquad \leqslant c_9\biggl(\|\alpha(v+w_0)\|_{\widetilde W_p^{s_0,2s_0}(S^\tau)} +\sum_{i=1}^r\|\nabla_{z'} \varphi_i\alpha(v+w_0)\|_{\widetilde W_p^{s_0,2s_0}(S_0^\tau)}\biggr), \end{aligned} \end{equation*} \notag $$
где постоянная $C_1$ не зависит от $\tau\in(0,T]$. Повторяя предыдущие рассуждения с использованием оценки (3.25), придем к неравенству вида
$$ \begin{equation*} \begin{aligned} \, \sum_{i=1}^r\|\nabla_{z'}\varphi_i v(t,x^i(z))\|_{W_p^{1,2}(Q_0^\tau)} &\leqslant c_{10}\tau^{1/2}\sum_{i=1}^r\|\nabla_{z'}\varphi_i v(t,x^i(z))\|_{W_p^{1,2}(Q_0^\tau)} \\ &\qquad{}+c_{11}\biggl(\sum_{i=1}^r \|\nabla_{z'}\varphi_i w_0\|_{W_p^{s_0,2s_0}(S_0^\tau)} +\|w_0\|_{W_p^{s_0,2s_0}(S)}\biggr), \end{aligned} \end{equation*} \notag $$
где постоянные $c_i$ не зависят от $\tau$ и зависят от величины $R_0$. Выберем $\tau_2\leqslant\tau_1$ такое, что $c_{10}\tau_2^{1/2}\leqslant 1/2$. Тогда при $\tau\leqslant\tau_2$ будет выполнено неравенство
$$ \begin{equation} \begin{aligned} \, &\sum_{i=1}^r\|\nabla_{z'}\varphi_iv(t,x^i(z))\|_{W_p^{1,2}(Q_0^\tau)} \nonumber \\ &\qquad{} \leqslant c_{12}\biggl(\|w_0\|_{W_p^{s_0,2s_0}(S)} +\sum_{i=1}^r\|\nabla_{z'}\varphi_iw_0(t,x^i(z',0))\|_{W_p^{s_0,2s_0}(S_0)}\biggr). \end{aligned} \end{equation} \tag{3.26} $$
Пусть $\vec\alpha_i=(\alpha_{i1},\alpha_{i2},\dots,\alpha_{ir})\in B_{R_0}$, $i=1,2$, и $v_i$ – соответствующие решения задачи (3.8), (3.9), где функция $\alpha$ заменяется на соответствующие функции
$$ \begin{equation*} \alpha^j=\sum_{i=1}^r\alpha_{ji}\Phi_i,\qquad j=1,2. \end{equation*} \notag $$
Тогда разности $v_1-v_2=\widetilde\omega$, $\widetilde\alpha=\alpha^1-\alpha^2$ есть решение задачи
$$ \begin{equation*} \begin{gathered} \, Mv=\widetilde\omega_t-L\widetilde\omega=0,\qquad (x,t)\in Q=G\times(0,T), \\ \frac{\partial\widetilde\omega}{\partial N}+\alpha_0\widetilde\omega =-\frac{(\alpha^1+\alpha^2)}{2}\,\widetilde\omega -\frac{\widetilde\alpha}{2}(v_1+v_2)-\widetilde\alpha w_0. \end{gathered} \end{equation*} \notag $$
В силу леммы 2 и (3.1) $\widetilde\alpha,\alpha^j\in\widetilde W_p^{s_0,2s_0}(S^\tau)$, $\nabla_{z'}\widetilde\alpha(t,x^i(z',0)), \nabla_{z'}\alpha^j(t,x^i(z',0))\in\widetilde W_p^{s_0,2s_0}(S_0^\tau)$ ($i=1,\dots,r$) и имеем оценки
$$ \begin{equation} \|\widetilde\alpha\|_{\widetilde W_p^{s_0,2s_0}(S^\tau)} \leqslant c_0\sum_{i=1}^r\|\alpha_{1i}-\alpha_{2i}\|_{\widetilde W_p^{s_0}(0,\tau)} \|\Phi_i\|_{W_p^{s_0,2s_0}(S)} \leqslant c_1\|\vec\alpha_1-\vec\alpha_2\|_{\widetilde W_p^{s_0}(0,\tau)}, \end{equation} \tag{3.27} $$
$$ \begin{equation} \|\alpha^1+\alpha^2\|_{\widetilde W_p^{s_0,2s_0}(S^\tau)} \leqslant c_1(\|\vec\alpha_1\|_{\widetilde W_p^{s_0}(0,\tau)} +\|\vec\alpha_2\|_{\widetilde W_p^{s_0}(0,\tau)})\leqslant 2c_1 R_0. \end{equation} \tag{3.28} $$
Аналогично имеем
$$ \begin{equation} \sum_{i=1}^r \|\nabla_{z'}\widetilde\alpha(t,x^i(z',0))\|_{\widetilde W_p^{s_0,2s_0}(S_0^\tau)} \leqslant c_2\|\vec\alpha_1-\vec\alpha_2\|_{\widetilde W_p^{s_0}(0,\tau)}, \end{equation} \tag{3.29} $$
$$ \begin{equation} \sum_{i=1}^r \|\nabla_{z'}(\alpha^1+\alpha^2)(t,x^i(z',0)\|_{\widetilde W_p^{s_0,2s_0}(S_0^\tau)} \leqslant 2c_2R_0, \end{equation} \tag{3.30} $$
где постоянная $c_2$ не зависит от $\tau$.

Повторяя рассуждения использованные при получении оценки (3.25), получим при $\tau\leqslant\tau_1$, что справедливо неравенство

$$ \begin{equation} \|\widetilde\omega\|_{W_p^{1,2}(Q^\tau)} \leqslant c_3\|\widetilde\alpha\|_{\widetilde W_p^{s_0}(0,\tau)}, \end{equation} \tag{3.31} $$
где $c_3$ не зависит от $\tau$. Аналогично при $\tau\leqslant\tau_2$ имеем оценку
$$ \begin{equation} \sum_{i=1}^r \|\nabla_{z'}\varphi_i\widetilde\omega(t,x^i(z))\|_{W_p^{1,2}(Q_0^\tau)} \leqslant c_4\biggl(\|\widetilde\alpha\|_{\widetilde W_p^{s_0,2s_0}(S^\tau)} +\sum_{i=1}^r \|\nabla_{z'}\widetilde{\alpha}(x^i(z',0))\|_{\widetilde W_p^{s_0,2s_0}(S_0^\tau)}\biggr). \end{equation} \tag{3.32} $$
Используя (3.27)(3.30), неравенства (3.31), (3.32) перепишем в виде
$$ \begin{equation} \|\widetilde\omega\|_{W_p^{1,2}(Q^\tau)}+\sum_{i=1}^r\|\nabla_{z'}\varphi_i \widetilde\omega(t,x^i(z))\|_{W_p^{1,2}(Q_0^\tau)} \leqslant c_5\|\vec\alpha_1-\vec\alpha_2\|_{\widetilde W_p^{s_0}(0,\tau)}. \end{equation} \tag{3.33} $$
Далее, изучим свойства отображения $\vec\alpha\to(w_1,w_2,\dots,w_r)$. Приведем оценку для правой части в (3.13) в $L_p((0,\tau)\times(0,\delta_1))$. Имеем
$$ \begin{equation*} \|f_1(t,0,z_n)\|_{L_p((0,\tau)\times(0,\delta_1))} \leqslant c_{6}\|f_1(t,z',z_n)\|_{W_p^s(B_{\delta}';L_p((0,\tau)\times(0,\delta_1)))}=J \end{equation*} \notag $$
при $s>(n-1)/p$ в силу леммы 3. Далее снова используем лемму 3, а точнее неравенства, вытекающие из соответствующих интерполяционных равенств. Имеем
$$ \begin{equation} J\leqslant c_7\|f_1(t,z)\|_{W_p^1(B'_\delta;L_p((0,\tau)\times(0,\delta_1)))}^\theta \|f_1(t,z)\|_{W_p^{-1}(B'_\delta;L_p((0,\tau)\times(0,\delta_1)))}^{1-\theta},\qquad 2\theta-1=s. \end{equation} \tag{3.34} $$
Исходя из определения $f_1$ и условий на коэффициенты, имеем
$$ \begin{equation} \|f_1\|_{W_p^{-1}(B_{\delta}';L_p((0,\tau)\times(0,\delta_1)))} \leqslant c\|v\|_{L_p(0,\tau;W^{1}_p(U))} \leqslant c_8\tau^{1/2}\|v\|_{W_p^{1,2}(Q^\tau)}, \end{equation} \tag{3.35} $$
где постоянная $c_1$ не зависит от $\tau$. Последняя оценка получается, если мы применим интерполяционное неравенство
$$ \begin{equation*} \|v\|_{L_p(0,\tau;W^{1}_p(U))} \leqslant c_9\|v\|_{L_p(0,\tau;W^2_p(U))}^{1/2}\|v\|_{L_p(0,\tau;L_p(U))}^{1/2} \end{equation*} \notag $$
и оценку
$$ \begin{equation*} \|v\|_{L_p(0,\tau;L_p(U))}\leqslant\tau\|v_{t}\|_{L_p(0,\tau;L_p(U)}, \end{equation*} \notag $$
вытекающую из формулы Ньютона–Лейбница, а затем оценим полученные нормы через норму в $W_p^{1,2}(Q^\tau)$. Оценки (3.34), (3.35) влекут, что
$$ \begin{equation} \|f_1(t,0,z_n)\|_{L_p((0,\tau)\times(0,\delta_1)))} \leqslant c_{10}\tau^{(1-\theta)/2}(\|\nabla_{z'}\varphi_iv(t,z)\|_{W_p^{1,2}(Q_0^\tau)} +\|v\|_{W_p^{1,2}(Q^\tau)}), \end{equation} \tag{3.36} $$
где $c_{10}$ – постоянная, не зависящая от $\tau$. Как вытекает из неравенств (3.25), (3.26), неравенство (3.36) можно переписать в виде
$$ \begin{equation*} \begin{aligned} \, &\|f_1(t,0,z_n)\|_{L_p((0,\tau)\times(0,\delta_1)))} \\ &\qquad\leqslant c_3\tau^{(1-\theta)/2}\biggl(\|w_0\|_{W_p^{s_0,2s_0}(S_0)} +\sum_{i=1}^r\|\nabla_{z'} \varphi_iw_0(x^i(z))\|_{W_p^{s_0,2s_0}((0,T)\times B_{\delta}')}\biggr). \end{aligned} \end{equation*} \notag $$
Пусть, как и ранее, $\vec\alpha_i=(\alpha_{i1},\alpha_{i2},\dots,\alpha_{ir})\in B_{R_0}$, $i=1,2$, и $v_i$ – соответствующие решения задачи (3.8), (3.9), где функция $\alpha$ заменяется на соответствующие функции $\alpha^j=\sum_{i=1}^r\alpha_{ji}\Phi_i$, $j=1,2$. Пусть $w_i^j$, $j=1,2$, решения задач (3.13), (3.14) c новыми правыми частями, где вместо $v$ стоят функции $v_i$ и $w^0=\varphi_i\widetilde\omega$. Тогда разности $k_i=w_i^1-w_i^2$ есть решения задач
$$ \begin{equation} k_{it}-c_{nn}(t,0,z_n)k_{i z_nz_n} =\sum_{i+j<2n}c_{ij}\omega^0_{z_iz_j} +\sum_{i=1}^nc_i\omega_{z_i}^0+c_0\omega^0 +[\varphi_i,L]\widetilde\omega|_{z'=0}=\widetilde f_i, \end{equation} \tag{3.37} $$
$$ \begin{equation} k_i|_{t=0} =0,\qquad k_i|_{z_n=0}=0,\qquad k_i|_{z_n=\delta_1}=0,\qquad i\leqslant r. \end{equation} \tag{3.38} $$
Из известных свойств параболических задач (см., например, [2]) имеем оценку
$$ \begin{equation} \sum_{i=1}^r\|k_i\|_{W_p^{1,2}((0,\tau)\times(0,\delta_1))} \leqslant\sum_{i=1}^r\|\widetilde f_i\|_{L_p((0,\tau)\times(0,\delta_1))}. \end{equation} \tag{3.39} $$
Используя аналог оценки (3.36) для оценки правой части, получим
$$ \begin{equation*} \sum_{i=1}^r\|k_i\|_{W_p^{1,2}((0,\tau)\times(0,\delta_1))} \leqslant c_2\tau^{(1-\theta)/2}\biggl(\|\widetilde\omega\|_{W_p^{1,2}(Q^\tau)} +\sum_{i=1}^r\|\nabla_{z'}\varphi_i\widetilde\omega(t,x^i(z))\|_{W_p^{1,2}(Q_0^\tau)}\biggr). \end{equation*} \notag $$
В частности, отсюда и из (3.33) вытекает неравенство
$$ \begin{equation} \sum_{i=1}^r\|k_{iz_n}(t,0)\|_{\widetilde W_p^{s_0}(0,\tau)} \leqslant c_4\tau^{(1-\theta)/2}\|\vec\alpha_1-\vec\alpha_2\|_{\widetilde W_p^{s_0}(0,\tau)}. \end{equation} \tag{3.40} $$
Возьмем $\vec\alpha=0$. Тогда в силе единственности решений задачи (3.8), (3.9) $v=v(\vec\alpha)=0$, в этом случае координаты вектора $\vec F(0)$ запишутся в виде
$$ \begin{equation*} F_i(0)=\frac{1}{\psi_i}(-\gamma_n(t,0)\omega_{iz_n}(t,0) -\alpha_0\widetilde\psi_i(t)),\qquad i=1,2,\dots,r, \end{equation*} \notag $$
где функция $\omega_i$ есть решение задачи
$$ \begin{equation*} w_{it}(t,z_n)-c_{nn}(t,0,z_n)w_{iz_nz_n}=0,\qquad w_i(0,z_n)=0,\quad w_i|_{z_n=0}=\widetilde\psi_i(t),\quad w_i|_{z_n=\delta_1}=0. \end{equation*} \notag $$

Положим $R_0=2\|\Phi^{-1}\vec F(0)\|_{\widetilde W_p^{s_0}(0,T)}$. В этом случае при $\tau\leqslant\tau_2$ будут справедливы все вышеприведенные оценки. Постоянная $\tau_2$ зависит от $R_0$ и норм данных задачи. Получим оценки, считая, что $\vec\alpha_i\in B_{R_0}$, $i=1,2$. Из (3.18) имеем, что

$$ \begin{equation*} \|R(\vec\alpha_1)-R(\vec\alpha_2)\|_{\widetilde W_p^{s_0}(0,\tau)} \leqslant c\sum_{i=1}^r\|F_i(\vec\alpha_1)-F_i(\vec\alpha_2)\|_{\widetilde W_p^{s_0}(0,\tau)}. \end{equation*} \notag $$
Как и ранее, обозначим соответствующие решения задачи (3.8), (3.9) через $v_i$, $i=1,2$, положим $\widetilde w=v_1-v_2$ обозначим через $\omega_j^i$, $i=1,2$, решения задач (3.13), (3.14), положим $k_i=\omega_i^1-\omega_i^2$. Рассмотрим первое слагаемое в координате $F_i(\vec\alpha_1)-F_i(\vec\alpha_2)$. Оно записывается в виде
$$ \begin{equation} J_1=-\gamma_n(\omega_{iz_n}^1-\frac{\omega_{iz_n}^2}{\psi_i(t)}\,, \end{equation} \tag{3.41} $$
где $\gamma_n=\sqrt{1+|\nabla\gamma|^2} \sum_{k,l=1}^n\widetilde a_{k,l}(y^i(z))\nu_k\nu_l|_{z_n=0}$. Здесь
$$ \begin{equation*} \nu_k=-\frac{\gamma_{z_k}(z')}{\sqrt{1+|\nabla\gamma|^2}}\,,\quad k<n,\qquad \nu_n=\frac1{\sqrt{1+|\nabla\gamma|^2}}\,, \end{equation*} \notag $$
и $\widetilde a_{k,l}$ – старшие коэффициенты оператора $L$ записанного в локальной системе координат $y$. В силу леммы 2 и неравенства (3.40) имеем
$$ \begin{equation*} \|J_1\|_{\widetilde W_p^{s_0}(0,\tau)} \leqslant c_1\|\omega_{iz_n}^1-\omega_{iz_n}^2\|_{\widetilde W_p^{s_0}(0,\tau)} \leqslant c_2\tau^{\beta_1}\|\vec\alpha_1-\vec\alpha_2\|_{\widetilde W_p^s(0,\tau)}, \end{equation*} \notag $$
где все постоянные не зависят от $\tau$ и $\beta_1$ – положительная постоянная. Рассмотрим второе слагаемое
$$ \begin{equation*} J_2=-\sum_{j=1}^{n-1}\gamma_j(t,0)\frac{\widetilde w_{z_j}(t,x^i(0))}{\psi_i(t)}\,. \end{equation*} \notag $$
В силу леммы 2 имеем
$$ \begin{equation*} \|J_2\|_{\widetilde W_p^{s_0}(0,\tau)} \leqslant c_2\sum_{i=1}^r(\|\widetilde w(t,x^i(0))\|_{\widetilde W_p^{s_0}(0,\tau)} +\|\nabla_{z'}\widetilde w(t,x^i(0))\|_{\widetilde W_p^{s_0}(0,\tau)}). \end{equation*} \notag $$
Каждое из слагаемых, входящих в эту сумму, оценивается одинаково. Отметим, что условие $v(t,x^i(z))\in\widetilde W_p^{s_1,2s_1}(S_0^\tau)$ влечет, что $\nabla_{z'}v\in\widetilde W_p^{s_0,2s_0}(S_0^\tau)$ и справедлива оценка (см. [37; лемма 7.2])
$$ \begin{equation*} \|\nabla_{z'}v\|_{\widetilde W_p^{s_0,2s_0}(S_0^\tau)} +\|v\|_{\widetilde W_p^{s_0,2s_0}(S_0^\tau)} \leqslant c\|v\|_{\widetilde W_p^{s_1,2s_1}(S_0^\tau)}, \end{equation*} \notag $$
где с помощью замены переменных легко убедиться, что постоянная $c$ не зависит от $\tau$. В частности, отсюда вытекает оценка (см. лемму 1)
$$ \begin{equation} \|v\|_{W_p^{1}(B'_{\delta};\widetilde W_p^{s_0}(0,\tau))} \leqslant c\|v\|_{\widetilde W_p^{s_1,2s_1}(S_0^\tau)} \leqslant c_1\|v\|_{W_p^{1,2}(Q^\tau)}. \end{equation} \tag{3.42} $$
Рассмотрим, например, последнее слагаемое в оценке для $\|J_2\|$. Имеем (см. лемму 3)
$$ \begin{equation*} \|\nabla_{z'}(\widetilde w(t,x^i(z)))|_{z=0}\|_{\widetilde W_p^{s_0}(0,\tau)} \leqslant c_2\|\nabla_{z'} \varphi_i(\widetilde w(t,x^i(z',0)))\|_{W_p^s(B_{\delta}';\widetilde W_p^{s_0}(0,\tau))}, \end{equation*} \notag $$
при $s>(n-1)/p$ (возьмем $s\in(n-1)/p,1)$). Далее, используя интерполяционное неравенство (следствие леммы 3) и (3.42) получим
$$ \begin{equation*} \begin{aligned} \, &\|\nabla_{z'}(\widetilde w(t,x^i(z)))|_{z=0}\|_{\widetilde W_p^{s_0}(0,\tau)} \leqslant c_2\|\nabla_{z'}(\widetilde w(t,x^i(z',0)))\|_{W_p^s(B'_\delta;\widetilde W_p^{s_0}(0,\tau))} \\ &\qquad\leqslant c_3\|\nabla_{z'}\varphi_i(\widetilde w(t,x^i(z',0)))\|_{W_p^1(B'_\delta;\widetilde W_p^{s_0}(0,\tau))}^\theta \|\nabla_{z'}\varphi_i(\widetilde w(t,x^i(z',0)))\|_{W_p^{-1}(B'_\delta;\widetilde W_p^{s_0}(0,\tau))}^{1-\theta} \\ &\qquad\leqslant c_4\|\nabla_{z'}\varphi_i(\widetilde w(t,x^i(z)))\|_{W_p^{1,2}(Q_0^\tau)}^\theta \|\varphi_i(\widetilde w(t,x^i(z',0)))\|_{L_p(B'_\delta;\widetilde W_p^{s_0}(0,\tau))}^{1-\theta} \\ &\qquad\leqslant c_5\|\nabla_{z'}\varphi_i(\widetilde w(t,x^i(z',0)))\|_{W_p^{1,2}(Q_0^\tau)}^\theta \tau^{(1-\theta)/2}\|\varphi_i(\widetilde w(t,x^i(z',0)))\|_{L_p(B'_\delta;\widetilde W_p^{s_1}(0,\tau))}^{1-\theta}. \end{aligned} \end{equation*} \notag $$
Ссылаясь на лемму 1 и используя (3.33), получим оценку
$$ \begin{equation*} \|\nabla_{z'}(\widetilde w(t,x^i(z)))|_{z=0}\|_{\widetilde W_p^{s_0}(0,\tau)} \leqslant c_6\tau^{\beta_2}\|\vec\alpha_1-\vec\alpha_2\|_{\widetilde W_p^{s_0}(0,\tau)}. \end{equation*} \notag $$
Аналогично оцениваются оставшиеся слагаемые в $\|J_2\|$, и можно сказать, что
$$ \begin{equation} \|J_2\|_{\widetilde W_p^{s_0}(0,\tau)} \leqslant c_7\tau^{\beta_2}\|\vec\alpha_1-\vec\alpha_2\|_{\widetilde W_p^{s_0}(0,\tau)}. \end{equation} \tag{3.43} $$
для некоторого $\beta_2>0$ и не зависящей от $\tau$ постоянной $c_7$. Окончательная оценка, как вытекает из (3.41), (3.43), имеет вид
$$ \begin{equation*} \|R(\vec\alpha_1)-R(\vec\alpha_2)\|_{\widetilde W_p^{s_0}(0,\tau)} \leqslant c_8\tau^{\beta_0}\|\vec\alpha_1-\vec\alpha_2\|_{\widetilde W_p^{s_0}(0,\tau)}, \end{equation*} \notag $$
где показатель $\beta_0$ минимальный из полученных в доказательстве, и постоянная $c_{10}$ не зависит от $\tau$. Возьмем в качестве $\tau_0$ число со свойством $\tau_0\leqslant\tau_2$, $c_{8}\tau_0^{\beta_0}\leqslant 1/2$. В этом случае оператор $R$ переводит шар $B_{R_0}$ в себя и является в нем сжимающим. Следовательно уравнение (3.17) имеет решение $\vec\alpha\in \widetilde W_p^{s_0}(0,\tau_0)$. Используя найденное решение $\vec\alpha$, найдем решение $v$ задачи (3.8), (3.9).

Покажем, что у нас выполнены условия переопределения (3.10). Возьмем равенства (3.9), записанные в системе координат $z$ и взятые в точке $t$, $x^i(0)$ ($x^i(0)=b_i$), и вычтем их из соответствующих равенств (3.15), получим

$$ \begin{equation} \gamma_n(t,0)(w_{iz_n}(t,0)-v_{z_n}(t,x^i(0))) +\alpha_0(w_i(t,0)-v(t,x^i(0)))+\alpha(w_i(t,0)-v(t,x^i(0)))=0, \end{equation} \tag{3.44} $$
где $i=1,2,\dots,r$. Функция $w_{0i}=\varphi_iv$ удовлетворяет уравнению (3.12). Возьмем в этом уравнении $z'=0$ и вычтем его из равенства (3.13). Получим равенства
$$ \begin{equation} (w_{it}(t,z_n)-w_{0it}(t,x^i(0,z_n))) -c_{nn}(t,0,z_n)(w_{iz_nz_n}-w_{0iz_nz_n}(t,x^i(0,z_n)))=0,\qquad i\leqslant r. \end{equation} \tag{3.45} $$
Функции $w_i(t,z_n)-w_{0i}(t,x^i(0,z_n))$ удовлетворяют уравнениям (3.45) начальному условию $(w_i(t,z_n)-w_{0i}(t,x^i(0,z_n)))|_{t=0}=0$ и в силу (3.44) граничным условиям
$$ \begin{equation*} \begin{gathered} \, \begin{aligned} \, &\gamma_n(t,0)(w_i(t,z_n)-w_{0i}(t,x^i(0))) +\alpha_0(w_i(t,z_n)-w_{0i}(t,x^i(0))) \\ &\qquad{} +\alpha(w_i(t,0)-w_{0i}(t,x^i(0)))=0,\qquad \end{aligned} \\ w_i(t,z_n)-w_{0i}(t,x^i(0,z_n))|_{z_n=\delta_1}=0,\quad i=1,2,\dots,r, \end{gathered} \end{equation*} \notag $$
в силу единственности решений смешанной начально-краевой задачи $w_i(t,z_n)=w_{0i}(t,x^i(0,z_n))$. Следовательно, выполнены равенства
$$ \begin{equation*} w_{0i}(t,x^i(0))=v(t,x^i(0))=\widetilde\psi_i\qquad \forall\,i. \end{equation*} \notag $$

Поскольку локально по времени задача сводится к уравнению со сжимающим оператором, то утверждение о единственности решений здесь очевидно.

Замечание 2. Утверждение теоремы 3 легко обобщается на случай систем параболических уравнений. В этом случае роль коэффициента $\beta$ играет соответствующая матрица, куда входят неизвестные функции, зависящие от времени. В самом общем случае элементы матрицы имеют вид $\beta_{ij}=\sum_{k=1}^{r_{ij}}\alpha_{ij}^k(t)\Phi_{ij}^k(x)$, где неизвестными являются функции $\alpha_{ij}^k(t)$. Схема рассуждений та же самая, однако, формулировка условий и результатов становится существенно более громоздкой.

СПИСОК ЦИТИРОВАННОЙ ЛИТЕРАТУРЫ

1. О. М. Алифанов, Обратные задачи в исследовании сложного теплообмена, Янус-К, М., 2009
2. В. Н. Ткаченко, Математическое моделирование, идентификация и управление технологическими процессами тепловой обработки материалов, Наукова думка, Киев, 2008  mathscinet
3. M. V. Glagolev, A. F. Sabrekov, “Determination of gas exchange on the border between ecosystem and atmosphere: inverse modeling”, Mat. Biolog. Bioinform., 7:11 (2012), 81–101  crossref
4. А. И. Бородулин, Б Д. Десятков, Г. А. Махов, С. Р. Сарманаев, “Определение эмиссии болотного метана по измеренным значениям его концентрации в приземном слое атмосферы”, Метеорология и гидрология, 1997, № 1, 66–74
5. L. V. Dantas, H. R. B. Orlande, R. M. Cotta, “An inverse problem of parameter estimation for heat and mass transfer in capillary porous media”, Int. J. Heat and Mass Transfer, 46:9 (2003), 1587–1599  crossref
6. J. Jr. Lugon, A. J. S. Neto, “An inverse problem of parameter estimation in simultaneous heat and mass transfer in a onedimensional porous medium”, Proceedings of COBEM 2003. 17-th International Congress on Mechanical Engineering, ABCM, San-Paolo, 2003; https://abcm.org.br/anais/cobem/2003/html/pdf/COB03-1158.pdf
7. K. Cao, D. Lesnic, M. J. Colaco, “Determination of thermal conductivity of inhomogeneous orthotropic materials from temperature measurements”, Inverse Probl. Sci. Eng., 27:10 (2018), 1372–1398  crossref  mathscinet
8. L. A. B. Varan, H. R. B. Orlande, F. L. V. Vianna, “Estimation of the convective heat transfer coefficient in pipelines with the Markov chain Monte-Carlo method”, Blucher Mech. Eng. Proc., 1:1 (2014), 1214–1225
9. A. M. Osman, J. V. Beck, “Nonlinear Inverse Problem for the Estimation of Time-and-Space-Dependent Heat-Transfer Coefficients”, J. Thermophysics, 3:2 (2003), 146–152
10. M. J. Colac, H. R. B. Orlande, “Inverse natural convection problem of simultaneous estimation of two boundary heat fluxes in irregular cavities”, Int. J. Heat and Mass Transfer, 47:6 (2004), 1201–1215  crossref
11. F. Avallone, C. S. Greco, D. Ekelschot, “2D Inverse Heat Transfer Measurements by IR Thermography in Hypersonic Flows”, Proceedings of the 11-th International Conference on Quantitative InfraRed Thermography, Naples, 2012, 1–13
12. S. D. Farahani, F. Kowsary, M. Ashjaee, “Experimental estimation heat flux and heat transfer coefficient by using inverse methods”, Sci. Iranica B, 3:4 (2016), 1777–1786  crossref
13. J. Su, G. F. Hewitt, “Inverse heat conduction problem of estimating time-varying heat transfer coefficient”, Numer. Heat Transfer A, 45 (2004), 777–789  crossref
14. D. N. Hao, R. X. Thanh, D. Lesnic, “Determination of the heat transfer coefficients in transient heat conduction”, Inverse Problems, 29:9 (2013), 095020  crossref  mathscinet
15. J. D. Lee, I. Tanabe, K. Takada, “Identification of the heat transfer coefficient on machine tool surface by inverse analysis”, JSME Int. J. Ser. C, 42:4 (1999), 1056–1060  crossref
16. T. M. Onyango, D. B. Ingham, D. Lesnic, “Restoring boundary conditions in heat conduction”, J. Eng. Math., 62:1 (2008), 85–101  crossref  mathscinet
17. S. Wang, L. Zhang, X. Sun, H. Jial, “Solution to two-dimensional steady inverse heat transfer problems with interior heat source based on the conjugate gradient method”, Math. Probl. Eng., 2017 (2017), 2861342  mathscinet
18. J. Sladek, V. Sladek, P. H. Wen, Y. C. Hon, “The Inverse problem of determining heat transfer coefficients by the meshless local Petrov–Galerkin method”, CMES Comput. Model. Eng. Sci., 48:2 (2009), 191–218  mathscinet
19. B. Jin, X. Lu, “Numerical identification of a Robin coefficient in parabolic problems”, Math. Comp., 81:279 (2012), 1369–1398  crossref  mathscinet
20. W. B. Da Silva, J. C. S. Dutra, C. E. P. Kopperschimidt, D. Lesnic, R. G. Aykroyd, “Sequential particle filter estimation of a time-dependent heat transfer coefficient in a multidimensional nonlinear inverse heat conduction problem”, Appl. Math. Model., 89:1 (2012), 654–668  mathscinet
21. D. N. Hao, B. V. Huong, P. X. Thanh, D. Lesnic, “Identification of nonlinear heat transfer laws from boundary observations”, Appl. Anal., 94:9 (2015), 1784–1799  crossref  mathscinet
22. M. Slodicka, R. Van Keer, “Determination of a Robin coefficient in semilinear parabolic problems by means of boundary measurements”, Inverse Problems, 18:1 (2002), 139–152  crossref  mathscinet
23. A. Rösch, “Stability estimates for the identification of nonlinear heat transfer laws”, Inverse Problems, 12:5 (1996), 743–756  crossref  mathscinet
24. D. Knupp, L. A S. Abreu, “Explicit boundary heat flux reconstruction employing temperature measurements regularized via truncated eigenfunction expansions”, Int. Commun. in Heat and Mass Transfer, 78 (2016), 241–252  crossref
25. С. А. Колесник, В. Ф. Формалев, Е. L. Кузнецова, “О граничной обратной задаче теплопроводности по восстановлению тепловых потоков к границам анизотропных тел”, ТВТ, 53:1 (2015), 72–77  mathnet  crossref
26. A. S. A. Alghamdi, “Inverse Estimation of Boundary Heat Flux for Heat Conduction Model”, JKAU. Eng. Sci., 21:1 (2010), 73–95  crossref
27. А. Б. Костин, А. И. Прилепко, “О некоторых задачах восстановления граничного условия для параболического уравнения. II”, Дифференц. уравнения, 32:11 (1996), 1519–1528  mathnet  mathscinet
28. А. Б. Костин, А. И. Прилепко, “О некоторых задачах восстановления граничного условия для параболического уравнения. I”, Дифференц. уравнения, 32:1 (1996), 107–116  mathnet  mathscinet
29. R. Denk, M. Huber, J. Prüss, “Optimal $L_p-L_q$-estimates for parabolic boundary value problems with inhomogeneous data”, Math. Z., 257 (2007), 193–224  crossref  mathscinet
30. Х. Трибель, Теория интерполяции, функциональные пространства, дифференциальные операторы, Мир, М., 1980  mathscinet  zmath
31. H. Amann, “Compact embeddings of vector-valued Sobolev and Besov spaces”, Glas. Mat. Ser. III, 35:1 (2000), 161–177  mathscinet
32. О. А. Ладыженская, Н. Н. Уральцева, В. А. Солонников, Линейные и квазилинейные уравнения параболического типа, Наука, М., 1967  mathscinet
33. В. А. Белоногов, С. Г. Пятков, “О разрешимости задач сопряжения с условиями типа неидеального контакта”, Изв. вузов. Матем., 2020, № 7, 18–32  mathnet  crossref
34. В. П. Михайлов, Дифференциальные уравнения в частных производных, Наука, М., 1976  mathscinet  zmath
35. М. А. Вержбицкий, С. Г. Пятков, “O некоторых обратных задачах об определении граничных режимов”, Матем. заметки СВФУ, 23:2 (2016), 3–18
36. С. М. Никольский, Приближение функций многих переменных и теоремы вложения, Наука, М., 1977  mathscinet
37. P. Grisvard, “Equations differentielles abstraites”, Ann. Scuola Norm. Sup. Pisa Cl. Sci. (4), 2:3 (1969), 311–395  mathscinet

Образец цитирования: С. Г. Пятков, В. А. Баранчук, “Определение коэффициента теплопередачи в математических моделях тепломассопереноса”, Матем. заметки, 113:1 (2023), 90–108; Math. Notes, 113:1 (2023), 93–108
Цитирование в формате AMSBIB
\RBibitem{PyaBar23}
\by С.~Г.~Пятков, В.~А.~Баранчук
\paper Определение коэффициента теплопередачи в~математических моделях тепломассопереноса
\jour Матем. заметки
\yr 2023
\vol 113
\issue 1
\pages 90--108
\mathnet{http://mi.mathnet.ru/mzm13573}
\crossref{https://doi.org/10.4213/mzm13573}
\mathscinet{http://mathscinet.ams.org/mathscinet-getitem?mr=4563351}
\transl
\jour Math. Notes
\yr 2023
\vol 113
\issue 1
\pages 93--108
\crossref{https://doi.org/10.1134/S0001434623010108}
\scopus{https://www.scopus.com/record/display.url?origin=inward&eid=2-s2.0-85149995814}
Образцы ссылок на эту страницу:
  • https://www.mathnet.ru/rus/mzm13573
  • https://doi.org/10.4213/mzm13573
  • https://www.mathnet.ru/rus/mzm/v113/i1/p90
  • Эта публикация цитируется в следующих 9 статьяx:
    Citing articles in Google Scholar: Russian citations, English citations
    Related articles in Google Scholar: Russian articles, English articles
    Математические заметки Mathematical Notes
    Статистика просмотров:
    Страница аннотации:208
    PDF полного текста:24
    HTML русской версии:150
    Список литературы:44
    Первая страница:18
     
      Обратная связь:
     Пользовательское соглашение  Регистрация посетителей портала  Логотипы © Математический институт им. В. А. Стеклова РАН, 2024