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

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

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



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






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


Математические заметки, 2024, том 116, выпуск 4, страницы 636–640
DOI: https://doi.org/10.4213/mzm14368
(Mi mzm14368)
 

Краткие сообщения

О существовании волнового фронта в задаче Коши для уравнения Гуртина–Пипкина

И. В. Романовa, А. С. Шамаевb

a Национальный исследовательский университет "Высшая школа экономики", г. Москва
b Институт проблем механики им. А. Ю. Ишлинского Российской академии наук, г. Москва
Список литературы:
Ключевые слова: уравнение Гуртина–Пипкина, волновой фронт.
Финансовая поддержка Номер гранта
Министерство науки и высшего образования Российской Федерации 075-15-2022-284
Работа выполнена при поддержке Министерства науки и высшего образования РФ (грант № 075-15-2022-284).
Поступило: 29.04.2024
Дата публикации: 11.10.2024
Англоязычная версия:
Mathematical Notes, 2024, Volume 116, Issue 4, Pages 862–866
DOI: https://doi.org/10.1134/S0001434624090414
Реферативные базы данных:
Тип публикации: Статья

1. Введение

Рассмотрим задачу Коши для одномерного уравнения Гуртина–Пипкина:

$$ \begin{equation} \theta_t(t,x)-\int^t_0K(t-s)\theta_{xx}(s,x)\,ds=0,\qquad x\in\mathbb{R},\quad t>0, \end{equation} \tag{1} $$
$$ \begin{equation} \theta|_{t=0}=\theta_0(x), \end{equation} \tag{2} $$
$$ \begin{equation} K(t)=\sum_{j=1}^{m}c_je^{-\gamma_jt},\qquad m\geqslant 2, \end{equation} \tag{3} $$
где $c_j$, $\gamma_j$ – заданные положительные постоянные, причем $\gamma_1<\gamma_2<\cdots<\gamma_m$. Ставится вопрос: имеется ли в данной задаче фронт распространения колебаний? Другими словами, является ли решение задачи (1)(3) финитной по $x$ функцией для каждого фиксированного конечного момента времени, если начальная функция $\theta_0$ финитна? В данной работе будет доказано, что это верно.

Если $m=1$, то уравнение (1) сводится к телеграфному уравнению, для которого известно, что фронт существует (см. [1]).

Постоянно в центре внимания был известный парадокс – скорость распространения тепла в классической модели теплопроводности бесконечна, что противоречит основным физическим принципам. Использование моделей теплопроводности с интегральным последствием в определенном смысле снимает это противоречие – скорость распространения тепла в этой модели становится конечной. Уравнение вида (1) (см. [2]) описывает не только тепловые процессы, но также и распространение звуковых волн в суспензиях, эмульсиях, взвесях и пр. (см. [3]). Поэтому вопрос о существовании фронта и скорости распространения фронта в уравнении (1) важен для целого ряда прикладных задач.

2. Существование фронта

Теорема 2.1. Пусть $\theta_0$ – достаточно гладкая и финитная функция. Тогда соответствующее ей решение задачи (1)(3) является финитной по $x$ функцией для каждого фиксированного конечного момента времени $t$.

Доказательство. Для доказательства воспользуемся хорошо известной теоремой Пэли–Винера (см. [4; с. 26, теорема X]): следующие два класса целых функций тождественны: Пусть $\widehat{\theta}(t,z)$ – преобразование Фурье от решения уравнения (1). Покажем сначала, что при каждом фиксированном $t>0$ это целая функция переменного $z$, затем, что она принадлежит классу (а).

Сделаем преобразование Фурье по переменной $x$ от обеих частей уравнения (1), затем проинтегрируем это уравнение по переменной $t$:

$$ \begin{equation} \widehat{\theta}(t,z)+z^2\int^t_0Q(t-s)\widehat{\theta}(s,z)\,ds= \widehat{\theta}_0(z),\qquad z\in\mathbb{C},\quad t>0, \end{equation} \tag{5} $$
$$ \begin{equation} Q(s)=\int^s_0K(\tau)\,d\tau. \end{equation} \tag{6} $$
Так как функция $Q$ непрерывна на $[0,+\infty)$, то для любого отрезка $[0,T]$ найдется $M>0$ такое, что
$$ \begin{equation*} |Q(s)|\leqslant M,\qquad s\in [0,T]. \end{equation*} \notag $$
Уравнение (5) представляет собой уравнение Вольтерры 2-го рода с параметром $z^2$. Свойства решения такого уравнения хорошо известны. Само решение можно представить в виде
$$ \begin{equation} \widehat{\theta}(t,z)=\widehat{\theta}_0(z) \sum_{k=0}^{+\infty}(-z)^{2k}(\mathbb{Q}^k 1)(t),\qquad t\in [0,T], \end{equation} \tag{7} $$
где $\mathbb{Q}^k$ – степень оператора $\mathbb{Q}$, который, в свою очередь, определяется следующим образом:
$$ \begin{equation*} (\mathbb{Q}g)(t)=\int^t_0Q(t-s)g(s)\,ds. \end{equation*} \notag $$
Кроме того,
$$ \begin{equation*} |(\mathbb{Q}^kg)(t)|\leqslant \|g\|_{C[0,T]}\frac{(Mt)^k}{k!}\,,\qquad t\in[0,T],\quad k=0,1,\dots\,. \end{equation*} \notag $$
Так как $\theta_0$ – финитная функция, то $\widehat{\theta}_0$ является целой. Следовательно, учитывая (7), получаем, что $\widehat{\theta}(t,\,\cdot\,)$ – целая для любого $t>0$.

Покажем теперь, что при каждом фиксированном $t>0$ имеет место представление $\widehat{\theta}(t,z)=O(e^{A|z|})$. Для этого сделаем преобразование Фурье по переменной $x$ от обеих частей уравнения (1), а по переменной $t$ – преобразование Лапласа. Тогда получим

$$ \begin{equation} \widetilde{\widehat{\theta}}(\lambda,z)= \frac{\widehat{\theta}_0(z)}{\lambda+z^2\widetilde{K}(\lambda)}\,, \end{equation} \tag{8} $$
где
$$ \begin{equation*} \widetilde{K}(\lambda)=\sum_{j=1}^{m}\frac{c_j}{\lambda+\gamma_j}\,. \end{equation*} \notag $$
Рассмотрим уравнение относительно неизвестного $\lambda\in\mathbb{C}$:
$$ \begin{equation} \frac{\lambda}{z^2}+\sum_{j=1}^{m}\frac{c_j}{\lambda+\gamma_j}=0. \end{equation} \tag{9} $$
В уравнении (9) комплексное число $z$ выступает в роли параметра. Так как (9) сводится к алгебраическому уравнению порядка $m+1$, то оно имеет $m+1$ корень. Покажем, что для достаточно больших по модулю значений $z$ эти корни попарно различны. Для этого будем использовать методы, близкие к тем, что применялись в [5] для похожей задачи (там число $z^2$ было вещественным). Введем обозначения
$$ \begin{equation*} \varphi(\lambda)=\frac{\lambda}{z^2}\,,\qquad \psi(\lambda)=\sum_{j=1}^{m}\frac{c_j}{\lambda+\gamma_j}\,. \end{equation*} \notag $$
Если $\lambda$ сделать вещественным переменным, то функция $\psi$ будет иметь $m-1$ попарно различных вещественных нулей $\mu_1,\mu_2,\dots,\mu_{m-1}$ таких, что
$$ \begin{equation} -\gamma_{j+1}<\mu_j<-\gamma_j,\qquad j=1,2,\dots,m-1. \end{equation} \tag{10} $$
На комплексной плоскости в окрестности каждой точки $\mu_j$ возьмем круг $D_j$ достаточно малого радиуса таким образом, чтобы в него не попали точки $-\gamma_j$, $j=1,2,\dots,m$. Тогда на каждом таком круге функция комплексного переменного $\psi$ будет голоморфной. Выбирая параметр $z_0\in\mathbb{C}$ достаточно большим по модулю, получим, что для любого $z$, для которого верно $|z|>|z_0|$, на границе каждого круга $D_j$ выполнено неравенство
$$ \begin{equation} |\psi(\lambda)|>|\varphi(\lambda)|. \end{equation} \tag{11} $$
Следовательно, по теореме Руше уравнение (9) имеет $m-1$ изолированных корней $\lambda_j(z)$, $j=1,\dots,m-1$.

Рассмотрим теперь круг $D$ на комплексной плоскости сколь угодно большого радиуса и с центром в нуле. Пусть все полюсы функции $\psi$ попадают внутрь $D$. Снова существует такое число $z_0$, что для любого $z$ такого, что $|z|>|z_0|$, на границе круга $D$ выполняется неравенство (11). По известной теореме из комплексного анализа имеем

$$ \begin{equation*} N(\psi)-P(\psi)=N(\psi+\varphi)-P(\psi+\varphi), \end{equation*} \notag $$
где $N$ и $P$ – количество нулей и полюсов соответственно внутри круга $D$. Из последнего равенства получаем, что $N(\psi+\varphi)=m-1$. Это означает, что оставшиеся два корня уравнения (9) стремятся к бесконечности по модулю при $|z|\to +\infty$. Обозначим эти два корня $\lambda_{m}(z)$ и $\lambda_{m+1}(z)$ соответственно. Найдем порядок их асимптотического роста. Перепишем (9) в виде
$$ \begin{equation*} \lambda=-\frac{z^2}{\lambda}\sum_{j=1}^{m} \frac{c_j}{1+\gamma_j/\lambda}\,. \end{equation*} \notag $$
Из последнего получим
$$ \begin{equation} \lambda^2=-z^2\sum_{j=1}^{m}c_j\biggl(1-\frac{\gamma_j}{\lambda}+ \frac{\gamma_j^2}{\lambda^2}+ o\biggl(\frac{1}{\lambda^2}\biggr)\biggr). \end{equation} \tag{12} $$
Из равенства (12) следует, что главный член асимптотики для модулей растущих корней имеет порядок $|z|$.

Покажем теперь, что $\lambda_{m}(z)$ и $\lambda_{m+1}(z)$ различны для достаточно больших по модулю значений $z$. Перепишем (9) в виде

$$ \begin{equation*} \frac{\lambda}{z^2}+\frac{S}{\lambda}+ \sum_{j=1}^{m}\frac{c_j}{\lambda+\gamma_j}-\frac{S}{\lambda}=0, \end{equation*} \notag $$
где
$$ \begin{equation*} S=\sum_{j=1}^{m}c_j. \end{equation*} \notag $$
Обозначим
$$ \begin{equation*} \varphi_1(\lambda)=\frac{\lambda}{z^2}+\frac{S}{\lambda}\,,\qquad \psi_1(\lambda)=\sum_{j=1}^{m}\frac{c_j}{\lambda+\gamma_j}- \frac{S}{\lambda}\,. \end{equation*} \notag $$
Пусть $|z|$ достаточно велико. У функции $\varphi_1$ имеется два нуля: $\pm iz\sqrt{S}$ . Рассмотрим сперва нуль $iz\sqrt{S}$ . Возьмем круг $G_{\varepsilon|z|}$ радиуса $\varepsilon |z|$ с центром в точке $iz\sqrt{S}$ , $0<\varepsilon<\sqrt{S}/2$. Предположим, что $|z|$ выбрано так, чтобы полюсы функции $\psi_1$ не попали в замыкание круга $G_{\varepsilon|z|}$. Пусть $\lambda\in\partial G_{\varepsilon|z|}$, т.е. $|\lambda-iz\sqrt{S}\,|=\varepsilon|z|$, тогда
$$ \begin{equation*} \begin{aligned} \, |\varphi_1(\lambda)|&=\frac{|\lambda-iz\sqrt{S}\,|\,|\lambda+ iz\sqrt{S}\,|}{|z^2|\,|\lambda|}= \frac{\varepsilon|z|\,|\lambda+iz\sqrt{S}\,|}{|z^2|\,|\lambda|}= \frac{\varepsilon|(\lambda-iz\sqrt{S}\,)+ 2iz\sqrt{S}\,|}{|z|\,|\lambda|} \\ &\geqslant\frac{\varepsilon(2|z|\sqrt{S}- \varepsilon|z|)}{|z|\,|\lambda|}= \frac{\varepsilon(2\sqrt{S}-\varepsilon)}{|\lambda|}\,. \end{aligned} \end{equation*} \notag $$
С другой стороны,
$$ \begin{equation*} \psi_1(\lambda)=\frac{1}{\lambda}\biggl(\sum_{j=1}^{m}\frac{c_j}{1+\gamma_j/\lambda}-S\biggr)= \frac{1}{\lambda}\biggl(S-\frac{\Gamma}{\lambda}-S+ o\biggl(\frac{1}{\lambda}\biggr)\biggr)= -\frac{\Gamma}{\lambda^2}+o\biggl(\frac{1}{\lambda^2}\biggr), \end{equation*} \notag $$
где
$$ \begin{equation*} \Gamma=\sum_{j=1}^{m}\gamma_j. \end{equation*} \notag $$
Из этого следует, что при растущем $|z|$ за счет соответствующего роста $|\lambda|\in\partial G_{\varepsilon|z|}$ можно добиться выполнения неравенства
$$ \begin{equation*} |\varphi_1(\lambda)|>|\psi_1(\lambda)|. \end{equation*} \notag $$
Следовательно, у функции $\varphi_1+\psi_1$ в круге $G_{\varepsilon|z|}$ также один нуль. Аналогично рассматривается случай для $-iz\sqrt{S}$ . Это и означает, что $\lambda_{m}(z)$ и $\lambda_{m+1}(z)$ различны для достаточно больших по модулю значений $z$ и не сближаются при $|z|\to+\infty$.

Таким образом, для достаточно больших значений $z$ из (8) следует представление

$$ \begin{equation} \begin{aligned} \, \widetilde{\widehat{\theta}}(\lambda,z)&=\widehat{\theta}_0(z) \sum_{j=1}^{m+1}\frac{B_j(z)}{\lambda-\lambda_j(z)}\,, \\ \nonumber B_j(z)&=\frac{\prod_{k=1}^{m}(\lambda_j(z)-\gamma_k)} {\prod_{k=1,k\neq j}^{m+1}(\lambda_j(z)-\lambda_k(z))}\,. \end{aligned} \end{equation} \tag{13} $$
Из (13) получим
$$ \begin{equation} \widehat{\theta}(t,z)=\widehat{\theta}_0(z) \sum_{j=1}^{m+1}B_j(z)e^{\lambda_j(z)t}. \end{equation} \tag{14} $$
Используя (14), а также в силу описанных выше свойств корней $\lambda_j(z)$ и выбора $\widehat{\theta}_0(z)$, получаем представление
$$ \begin{equation} \widehat{\theta}(t,z)=O(e^{A|z|}),\qquad |z|\to+\infty. \end{equation} \tag{15} $$

Докажем теперь, что $\widehat{\theta}(t,\,\cdot\,)\in L_2(\mathbb{R})$. Для этого в уравнении (9) положим $z=x\in\mathbb{R}$. В этом случае для всех достаточно больших по модулю $x$ соответствующие корни $\lambda_j(x)$, $j=1,\dots,m+1$, имеют следующие асимптотики ($|x|\to+\infty$):

$$ \begin{equation*} \lambda_j(x)=\mu_j+o(1),\qquad j=1,\dots,m-1, \end{equation*} \notag $$
$$ \begin{equation*} \lambda_j(x)=\pm i\sqrt{S}\,|x|-\frac{1}{2\sqrt{S}}\sum_{j=1}^{m} c_j\gamma_j+o(1),\qquad j=m,m+1. \end{equation*} \notag $$
Истинность асимптотик для $j=1,\dots,m-1$ была установлена на стр. 3, представления для $j=m,m+1$ следуют из (12). Вопрос об асимптотиках для похожей ситуации (в более общем случае) рассматривается, например, в работе [5]. Из всего этого и из того, что $\widehat{\theta}_0\in L_2(\mathbb{R})$, получим требуемое.

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

1. В. И. Смирнов, Курс высшей математики, 2, БХВ-Петербург, 2008  mathscinet
2. M. E. Gurtin, A. C. Pipkin, Arch. Rational Mech. Anal., 31:2 (1968), 113–126  crossref  mathscinet
3. Э. Санчес-Паленсия, Неоднородные среды и теория колебаний, Мир, М., 1984  mathscinet
4. Н. Винер, Р. Пэли, Преобразование Фурье в комплексной плоскости, Наука, М., 1964  mathscinet
5. В. В. Власов, Н. А. Раутиан, А. С. Шамаев, Уравнения в частных производных, СМФН, 39, РУДН, М., 2011, 36–65  mathnet  mathscinet

Образец цитирования: И. В. Романов, А. С. Шамаев, “О существовании волнового фронта в задаче Коши для уравнения Гуртина–Пипкина”, Матем. заметки, 116:4 (2024), 636–640; Math. Notes, 116:4 (2024), 862–866
Цитирование в формате AMSBIB
\RBibitem{RomSha24}
\by И.~В.~Романов, А.~С.~Шамаев
\paper О~существовании волнового фронта в~задаче~Коши для
уравнения Гуртина--Пипкина
\jour Матем. заметки
\yr 2024
\vol 116
\issue 4
\pages 636--640
\mathnet{http://mi.mathnet.ru/mzm14368}
\crossref{https://doi.org/10.4213/mzm14368}
\transl
\jour Math. Notes
\yr 2024
\vol 116
\issue 4
\pages 862--866
\crossref{https://doi.org/10.1134/S0001434624090414}
\scopus{https://www.scopus.com/record/display.url?origin=inward&eid=2-s2.0-85213313826}
Образцы ссылок на эту страницу:
  • https://www.mathnet.ru/rus/mzm14368
  • https://doi.org/10.4213/mzm14368
  • https://www.mathnet.ru/rus/mzm/v116/i4/p636
  • Citing articles in Google Scholar: Russian citations, English citations
    Related articles in Google Scholar: Russian articles, English articles
    Математические заметки Mathematical Notes
     
      Обратная связь:
     Пользовательское соглашение  Регистрация посетителей портала  Логотипы © Математический институт им. В. А. Стеклова РАН, 2025