Постоянно в центре внимания был известный парадокс – скорость распространения тепла в классической модели теплопроводности бесконечна, что противоречит основным физическим принципам. Использование моделей теплопроводности с интегральным последствием в определенном смысле снимает это противоречие – скорость распространения тепла в этой модели становится конечной. Уравнение вида (1) (см. [2]) описывает не только тепловые процессы, но также и распространение звуковых волн в суспензиях, эмульсиях, взвесях и пр. (см. [3]). Поэтому вопрос о существовании фронта и скорости распространения фронта в уравнении (1) важен для целого ряда прикладных задач.
Доказательство. Для доказательства воспользуемся хорошо известной теоремой Пэли–Винера (см. [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})$, получим требуемое.