|
Эта публикация цитируется в 7 научных статьях (всего в 7 статьях)
Определение коэффициента теплопередачи в математических моделях тепломассопереноса
С. Г. Пятков, В. А. Баранчук Югорский государственный университет, г. Ханты-Мансийск
Аннотация:
В работе рассматривается вопрос о корректности в пространствах Соболева
обратных задач об определении коэффициента теплопередачи, входящего
в граничное условие типа Робина для уравнений конвекции–диффузии.
Доказана теорема существования и единственности решений.
Библиография: 37 названий.
Ключевые слова:
обратная задача, коэффициент теплопередачи, параболическое уравнение,
тепломассоперенос, диффузия.
Поступило: 03.05.2022
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 |
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 |
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 |
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 |
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 |
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 |
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 |
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 |
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 |
16. |
T. M. Onyango, D. B. Ingham, D. Lesnic, “Restoring boundary conditions in heat conduction”, J. Eng. Math., 62:1 (2008), 85–101 |
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 |
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 |
19. |
B. Jin, X. Lu, “Numerical identification of a Robin coefficient in parabolic problems”, Math. Comp., 81:279 (2012), 1369–1398 |
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 |
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 |
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 |
23. |
A. Rösch, “Stability estimates for the identification of nonlinear heat transfer laws”, Inverse Problems, 12:5 (1996), 743–756 |
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 |
25. |
С. А. Колесник, В. Ф. Формалев, Е. L. Кузнецова, “О граничной обратной задаче теплопроводности по восстановлению тепловых потоков к границам анизотропных тел”, ТВТ, 53:1 (2015), 72–77 |
26. |
A. S. A. Alghamdi, “Inverse Estimation of Boundary Heat Flux for Heat Conduction Model”, JKAU. Eng. Sci., 21:1 (2010), 73–95 |
27. |
А. Б. Костин, А. И. Прилепко, “О некоторых задачах восстановления граничного условия для параболического уравнения. II”, Дифференц. уравнения, 32:11 (1996), 1519–1528 |
28. |
А. Б. Костин, А. И. Прилепко, “О некоторых задачах восстановления граничного условия для параболического уравнения. I”, Дифференц. уравнения, 32:1 (1996), 107–116 |
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 |
30. |
Х. Трибель, Теория интерполяции, функциональные пространства, дифференциальные операторы, Мир, М., 1980 |
31. |
H. Amann, “Compact embeddings of vector-valued Sobolev and Besov spaces”, Glas. Mat. Ser. III, 35:1 (2000), 161–177 |
32. |
О. А. Ладыженская, Н. Н. Уральцева, В. А. Солонников, Линейные и квазилинейные уравнения параболического типа, Наука, М., 1967 |
33. |
В. А. Белоногов, С. Г. Пятков, “О разрешимости задач сопряжения с условиями типа неидеального контакта”, Изв. вузов. Матем., 2020, № 7, 18–32 |
34. |
В. П. Михайлов, Дифференциальные уравнения в частных производных, Наука, М., 1976 |
35. |
М. А. Вержбицкий, С. Г. Пятков, “O некоторых обратных задачах об определении граничных режимов”, Матем. заметки СВФУ, 23:2 (2016), 3–18 |
36. |
С. М. Никольский, Приближение функций многих переменных и теоремы вложения, Наука, М., 1977 |
37. |
P. Grisvard, “Equations differentielles abstraites”, Ann. Scuola Norm. Sup. Pisa Cl. Sci. (4), 2:3 (1969), 311–395 |
Образец цитирования:
С. Г. Пятков, В. А. Баранчук, “Определение коэффициента теплопередачи в математических моделях тепломассопереноса”, Матем. заметки, 113:1 (2023), 90–108; Math. Notes, 113:1 (2023), 93–108
Образцы ссылок на эту страницу:
https://www.mathnet.ru/rus/mzm13573https://doi.org/10.4213/mzm13573 https://www.mathnet.ru/rus/mzm/v113/i1/p90
|
Статистика просмотров: |
Страница аннотации: | 201 | PDF полного текста: | 22 | HTML русской версии: | 144 | Список литературы: | 43 | Первая страница: | 18 |
|