尽管目前 HIV/AIDS预防和治疗取得了巨大进展, 但仍有一些低收入和中等收入国家受到 HIV 流行的影响. 据估计, 截止到2020 年底全球有 3760 万 HIV 感染者, 大约 680,000 人死于与AIDS相关的疾病, 其中大多数在东部和南部非洲. 在影响 HIV/AIDS 传播和分布的所有因素中, 人群的扩散被认为是一个重要因素. 这促使关于HIV/AIDS 空间传播动力学的研究工作广泛展开. 例如, Cuadros 等^[2] 利用一种基于空间变量的模型在几个非洲国家创建高分辨率的艾滋病毒感染率估计图. Ayalew 等^[3] 比较了三种空间 HIV 模型, 以估计南非不同地区的艾滋病毒感染率. Wu 和 Zhao^[4] 研究了一类带有抗逆转录病毒治疗的年龄结构空间艾滋病毒流行病模型的动力学行为. Zhao 等^[5] 建立了一类具有三年龄段的 HIV/AIDS 流行病模型, 并进行了动力学和最优控制分析. 其他一些时空结构 HIV 感染模型也已被用来研究 HIV 在宿主体内的传播动力学(参见文献[6⇓⇓⇓-10]).

事实上, 随着时间的推移, HIV 感染者可能会发展成艾滋病或保持在无症状阶段, 这在很大程度上取决于感染者的感染年龄, 即个体感染 HIV 后的时间. 因此有必要将感染年龄纳入空间 HIV/AIDS 模型, 以揭示该疾病传播的动力学机制, 这就是所谓的“年龄-空间结构”. 除了文献[4-5]所做的研究工作之外, 年龄-空间结构也被一些学者在建立其他传染病模型时所考虑. 最近, Chekroun和Kuniya^[11]研究了具有 Neumann 边界条件的年龄空间结构SIR 流行病模型的全局动力学. Yang等^[12]研究了年龄结构空间布鲁氏菌病模型的阈值动力学行为. Liu 等^[13]研究了具有年龄结构和空间扩散的多群体SEIR流行病模型的全局演化行为. Wang等^[14]建立了一类具有年龄-空间结构的 HIV 感染模型并研究了模型解的全局动力学行为. 值得提出的是, 上述这些模型都是考虑了齐次 Neumann 边界条件. 然而, 当空间区域的边界具有对个体生存不利的环境 (例如, 沼泽、沙漠或高山) 时, 模型更适合考虑齐次狄利克雷边界条件 (参见文献[15-16]). 有鉴于此, 具有狄利克雷边界条件的空间模型已被用于研究相关传染病的动力学传播. 例如, Chekroun 和 Kuniya^[17]分析了 Dirichlet 边界条件下具有年龄-空间结构 SIR 反应扩散传染病模型的阈值动力学. Wang等^[18]在齐次 Dirichlet 边界条件下分析了年龄-空间结构口蹄疫模型的时空动力学.

受到上述研究工作的启发, 为了研究人群扩散、感染年龄和 Dirichlet 边界条件环境对人群中 HIV/AIDS 传播的综合影响, 本文建立如下模型

模型的初值条件和齐次 Dirichlet 边界条件如下所示

其中, $\overline{\Gamma}$ 和 $\partial \Gamma$ 分别代表有界连通开集合 $\Gamma\subset {\Bbb R} ^n$ 的闭包和光滑边界, 其中 $n\in{\mathbb N}_+$. 变量 $U_1(x,t),u_2(x,b,t),u_3(x,b,t),U_4(x,t)$ 分别表示在时间 $t$ 和位置$x$ 处感染年龄为 $b$ 的易感个体密度, 初期 HIV 感染个体密度, 慢性无症状感染个体密度和患有 AIDS 病症且性活动减少的个体密度(参见文献[19]). 参数 $d_j(j=1,2,3,4)$ 分别代表这四类人群的个体扩散率. 参数 $\Pi$ 表示易感个体的输入率, 易感个体的自然死亡率用 $\sigma$ 表示. 变量 $u_2$ 和 $u_3$ 因疾病导致的依赖年龄的死亡率分别为 $\delta_2(b)$ 和 $\delta_3(b)$. 变量 $U_4$ 的疾病导致的死亡率为 $\delta_4$. 参数 $\eta_i(b)$ 和 $\gamma_i(b)$ 分别代表变量 $u_i$, $i=2,3$ 中感染个体的转移率和感染率.

利用特征线法, 将系统 (1.1)-(1.2) 改写为积分反应扩散方程. 为此, 令 $b-t=a$, $u_i(x,b,t)=u_i(x,t+a,t)$, $\eta_i(t)=\eta_i(t+a),\delta_i(t)=\delta_i(t+a)$, $i=2,3$. 类似文献[11-12,20-21] 中的推导过程, 我们可以得到 $u_i(x,b,t)$ 表达式如下

其中 $\Omega_i(b)=\exp\{-\int_0^b(\sigma+\eta_i(s)+\delta_i(s))ds\}$, $G_i(x,y,t)$ 是关于拉普拉斯算子 $d_i\Delta u_i$ 的格林函数, $i=2,3$, $G_i(x,y,t)$ 函数的性质参加文献[22-23].

令 $U_i(x,t)=u_i(x,0,t)$ 并将 (1.3) 带入到系统 (1.1) 第四个方程中可得

由于 $U_3(x,t)$ 和 $U_4(x,t)$ 被 $U_2(x,t)$ 所确定, 因此我们接下来研究如下积分反应扩散方程模型, 其中 $x\in\Gamma, t>0$,有

其中

接下来, 我们致力于研究系统 (1.5) 的动力学行为.

为了研究系统 (1.5) 的适定性,我们首先做出如下假设.

${\bf定义 2.1}$ 我们假设

(1) 模型参数 $\Pi,\sigma,\delta_4$, $d_j\ (j=1,2,3,4)$ 都是正的;

(2) $\gamma_i(\cdot)\in L^1_+({\Bbb R} _+)\cap L^{\infty}_+({\Bbb R} _+)$. 此外, 存在 $0<b_1<b_2$ 使得 $\gamma_i(b)>0, b\in (b_1,b_2)$;

(3) $\eta_i(\cdot),\delta_i(\cdot)\in L^{\infty}_+({\Bbb R} _+)$, $m^{\infty}=ess\cdot\sup\limits_{b\in{\Bbb R} _+}m(b)<\infty$, 其中 $m(\cdot)=\gamma_i(\cdot),\eta_i(\cdot),\delta_i(\cdot),$$i=2,3$.

我们令空间 ${\mathbb M}\triangleq C({\Bbb R},\overline{\Gamma})$ 具有范数 $\| \psi\| _{{\mathbb M}}\triangleq\sup\limits_{x\in\Gamma}|\psi(x)|,\psi\in{\mathbb M}$. 定义空间 ${\mathbb Q}\triangleq BC({\Bbb R} _+,{\mathbb M})$, 其范数为 $\| \phi\| _{{\mathbb Q}}\triangleq\sup\limits_{b\ge 0}\| \phi(b)\| _{{\mathbb M}}=\sup\limits_{(x,b)\in\Gamma\times{\Bbb R} _+}|\phi(x,b)|$, $\phi\in{\mathbb Q}$. 定义 ${\mathbb M}^0\triangleq\{\psi\in{\mathbb M}|\psi(x)=0,$$x\in\partial\Gamma\}$, ${\mathbb Q}^0\triangleq\{\phi\in{\mathbb Q}|\phi(x,\cdot)=0,$$x\in\partial \Gamma\}$ 分别是空间 ${\mathbb M}$ 和 ${\mathbb Q}$ 的正锥. 在空间 ${\mathbb Q}^0$ 中定义线性算子 $\Phi$

根据格林函数 $G_i(x,y,t)$ 性质以及假设2.1, 易知 $\Phi({\mathbb Q}^0)\subset {\mathbb Q}^0_+$, 其中 ${\mathbb Q}^0_+$ 是空间 ${\mathbb Q}^0$ 的正锥. 此外, 由 $\sigma+(\Phi U_2)(x,t)+{\cal F}_1+{\cal F}_2$ 的有解性和连续性, 再根据强极大值原理^[24] 和假设 2.1, 我们可得如下结论.

${\bf引理 2.1}$ 如果 $(U_1,U_2)$ 是系统 (1.5) 具有 $(\psi_1,\psi_2)\in{\mathbb M}^0\times{\mathbb Q}^0$ 初值条件的解, 那么存在 $t_0>0$ 使得 $U_1(x,t)>0,U_2(x,t)\ge 0$, $(x,t)\in\Gamma\times(0,t_0]$.

进一步地, 我们可得系统(1.5)解的全局存在性和一致有界性的结论.

${\bf引理 2.2}$ 如果假设 2.1 成立且 $(\psi_1,\psi_2,\psi_3)\in{\mathbb M}^0\times{\mathbb Q}^0\times{\mathbb Q}^0$, 那么系统 (1.5)有正解 $(U_1,U_2)$, $(x,t)\in\Gamma\times[0,\infty)$. 此外, $U_1$ 和 $U_2$ 都是一致有界的, 即存在两个正数 ${\cal B}_1={\cal B}_1(\psi_1,\psi_2,\psi_3),$${\cal B}_2={\cal B}_2(\psi_1,\psi_2,\psi_3)$, $(\psi_1,\psi_2,\psi_3)\in{\mathbb M}^0_+\times{\mathbb Q}^0_+\times{\mathbb Q}^0_+$, $(x,t)\in\Gamma\times(0,+\infty)$, 有 $0<U_1(x,t)\le{\cal B}_1,0<U_2(x,t)\le{\cal B}_2$ 成立.

${\bf证}$ 从系统 (1.5) 第一个方程中可得

其中 ${\cal F}_3(x,t)=e^{-\sigma t}\int_{\Gamma}G_1(x,y,t)\psi_1(y){\rm d}y$. 于是, $U_2(x,t)$ 可以被表示为

由此, 我们可知当(2.1)式存在解 $U_2(x,t)$, 那么 $U_1(x,t)$ 就可以被确定.

令 ${\cal F}^{\infty}_k:=\sup\limits_{(x,t)\in\Omega\times{\Bbb R} _+}{\cal F}_k(x,t),k=1,2,3$. 由于 $(\psi_1,\psi_2,\psi_3)\in{\mathbb M}^0_+\times{\mathbb Q}^0_+\times{\mathbb Q}^0_+$ 以及格林函数 $G_i(x,y,t)\ (i=2,3)$ 的性质可知 ${\cal F}^{\infty}_k$ 是有界的. 设 ${\cal B}>{\cal F}^{\infty}_3({\cal F}^{\infty}_1+{\cal F}^{\infty}_2)$ 是足够大的数以及

注意到 $f_1(x,0)={\cal F}_3(x,0),f_2(x,0)={\cal F}_1(x,0)+{\cal F}_2(x,0)$, 我们可选取足够小的 $T>0$ 使得

对任意固定的 $T$, 我们令 ${\mathbb M}_T:=C({\mathbb M},[T])$ 具有范数

定义

以及 ${\mathbb M}^0_{T,{\cal B}}$ 中的线性算子 ${\cal H}$, $(x,t)\in\bar{\Gamma}\times(0,T]$,有

令 $\psi\in{\mathbb M}^0_{T,{\cal B}}$, 我们可得 $({\cal H}\psi)(x,\cdot)=0,\ x\in\partial \Gamma$,有

这意味着 ${\cal H}\psi\in{\mathbb M}^0_{T,{\cal B}}$. 由于 ${\cal H}({\mathbb M}^0_{T,{\cal B}})\subset{\mathbb M}^0_{T,{\cal B}}$ 我们可知 ${\cal H}$ 在 ${\mathbb M}^0_{T,{\cal B}}$ 中是严格收缩的. 为了延拓解的存在区域, 我们还需证明解 $(U_1,U_2)$ 永远不会爆破. 从系统 (1.5) 第一个方程中, 可知 $U_1(x,t)\le \Pi/\sigma(1-e^{-\sigma t})+e^{-\sigma t}\| \psi_1\| _{{\mathbb M}}<+\infty$. 于是, 存在 ${\cal B}_1:=\Pi/\sigma+\| \psi_1\| _{{\mathbb M}}$ 使得 $0<U_1(x,t)\le {\cal B}_1$.

假设 $U_2(x,t)$ 会爆破, 则存在 $\tilde{t}>0, \tilde{x}\in\Gamma$ 使得 $\lim\limits_{t\to\tilde{t}^{-}}\partial_tU_1(\tilde{x},t)=-\infty$, 这意味着 $U_1(\tilde{x},t)$ 在 $\tilde{t}$ 邻域中会变为负的. 这显然与 $U_1(x,t)$ 的正性相矛盾. 因此, $U_2(x,t)$ 不会爆破. 于是, 存在 ${\cal B}_2$ 使得 $0<U_2(x,t)\le{\cal B}_2$. 引理证毕.

${\bf定理 2.1}$ 如果假设条件 2.1 成立, 那么系统 (1.5) 存在唯一的具有初始条件 $(\psi_1,\psi_2,\psi_3)\in{\mathbb M}^0_+\times{\mathbb Q}^0_+\times{\mathbb Q}^0_+$ 的有界解 $(U_1,U_2)$, $(x,t)\in\Gamma\times[0,+\infty)$.

${\bf证}$ 我们已经在引理 2.2 中证明了系统 (1.5) 经典解的全局存在性. 现在我们利用反证法证明解的唯一性. 假设系统 (1.5) 有两个解 $(\widetilde{U}_1,\widetilde{U}_2)$ 和 $(\widehat{U}_1,\widehat{U}_2)$. 从引理 2.2 中我们可知存在 ${\cal B}^+>0$ 使得 $0\le\widetilde{U}_2,\widehat{U}_2\le {\cal B}_+$, $(x,t)\in\Gamma\times(0,+\infty)$. 注意到 $\widetilde{U}_2={\cal H}\widetilde{U}_2, \widehat{U}_2={\cal H}\widehat{U}_2$. 令 $U^*_2=\widetilde{U}_2-\widehat{U}_2$, 可得

于是利用 Gronwall 不等式, 可得 $\| U^*_2(b)\| _{{\mathbb M}}=0,t>0$. 这意味着 $\widetilde{U}_2=\widehat{U}_2$. 类似地, 可得 $U_1(x,t)$ 的唯一性. 定理证毕.

此节中, 我们致力于推导出系统 (1.5) 的基本再生数泛函表达式. 根据命题 4.1 和引理 4.1^[17], 可知方程

有唯一的非负解 $U_{10}\in{\mathbb M}^0_+\setminus\{0\}$, 其在 $\Gamma$ 中是严格正的且满足 $\limsup\limits_{t\to+\infty}U_1(x,t)\le U_{10}(x)$, $x\in\overline{\Gamma}$. 系统 (1.5) 中 $U_2(x,t)$ 在无病平衡态 $E_0(U_{10}(x),0)$ 处的线性化方程

根据文献[25], 我们定义再生算子 ${\cal R}:{\mathbb M}^0\to{\mathbb M}^0$

其中 $\psi\in{\mathbb M}^0, x\in\overline{\Gamma}$. 由文献[11-12]中的结论可知系统 (1.5) 的基本再生数 $R_0$ 定义为 $R_0=r({\cal R})$, 其中 $r({\cal R})$ 是算子 ${\cal R}$ 的谱半径. 此外, 根据文献[17,引理 4.2,命题 4.2], 我们有如下结论.

${\bf定理 3.1}$ 若算子 ${\cal R}$ 如(3.1)式所定义, 则 ${\cal R}$ 是强正的, 有界的和紧的. 此外, $R_0$ 是 ${\cal R}$ 唯一简单的特征值, 其具有正特征函数 $\zeta_0\in{\mathbb M}^0_+\setminus\{0\}$.

此节中, 我们讨论系统 (1.5) 无病平衡态 $E_0(U_{10}(x),0)$ 的全局吸引性. 为此, 我们首先给出如下引理.

${\bf引理 4.1}$ 如果 $(\psi_1,\psi_2,\psi_3)\in{\cal D}$, $R_0<1$, 那么有

其中 $\epsilon>0$ 是足够大的常数, ${\cal D}$如下所示

${\bf证}$ 根据引理 2.2 易知 $0\le U_1(x,t)\le U_{10}(x)$. 从 $U_2(x,0)$ 的方程中我们可得

假设存在 $t_2>0$ 和 $x_2\in\Gamma$ 使得 $U_2(x,t)<\epsilon\zeta_0(x)$, $t\in(0,t_2)$ 以及 $U_2(x,t_2+\rho)>\epsilon\zeta_0(x)$ 当 $\rho>0$ 足够小. 从而可得

引理证毕.

${\bf定理 4.1}$ 假设条件 2.1 成立, $\epsilon$ 如引理 4.1 中所示. 如果 $R_0<1$ 和 $(\psi_1,\psi_2,\psi_3)\in{\cal D}$, 那么 $E_0$ 是全局吸引的, 即 $\lim\limits_{t\to\infty}\| U_1(\cdot,t)-U_{10}(x)\| _{{\mathbb M}}=0, \lim\limits_{t\to\infty}\| U_2(\cdot,t)\| _{{\mathbb M}}=0$.

${\bf证}$ 设 $U^{\infty}_2(x)=\limsup\limits_{t\to\infty}U_2(x,t),x\in\overline{\Gamma}$. 显然地, $U^{\infty}_2(x)\le\epsilon\zeta_0(x)$. 于是对任意的 $x\in\overline{\Gamma}$, 我们有

这意味着 $\lim\limits_{t\to\infty}\| U_2(\cdot,t)\| _{{\mathbb M}}=0$. 于是, 对任意的 $\varrho>0$, 存在 $T>0$, $\| U_2(\cdot,t)\| _{{\mathbb M}}\le \varrho$, $t\ge T$, 使得

继而可得

其中 $U^{\varrho}_{10}(x)=(\Pi-\varrho)\int_0^{\infty}e^{-\sigma b} \int_{\Omega}G_1(x,y,b){\rm d}y{\rm d}b$. 由于 $\varrho>0$, 可知 $\lim\limits_{t\to\infty}\| U_1(\cdot,t)-U_{10}(\cdot)\| _{{\mathbb M}}=0$. 定理证毕.

此节中, 我们证明当 $R_0>1$ 时, HIV 感染在人群中会持久存在的. 我们首先给出如下结论.

${\bf引理 5.1}$ 如果假设 2.1 成立, $(\psi_1,\psi_2,\psi_3)\in{\mathbb M}^0_+\times{\mathbb Q}^0_+\times{\mathbb Q}^0_+$, 那么存在一个紧的连续解半流 $\Sigma(t,\psi_1,\psi_2,\psi_3)=(U_1(\cdot,t),U_2(\cdot,t)) \in{\mathbb M}^0_+\times{\mathbb Q}^0_+,t>0$.

${\bf证}$ 根据文献[11,引理 5.1], 可得 $\Sigma(t,\psi_1,\psi_2,\psi_3)$的存在性. 接下来, 我们证明半流的紧性. 设 ${\cal B}\subset{\mathbb M}^0_+\times{\mathbb Q}^0_+\times{\mathbb Q}^0_+$ 是一个有界集, 定义

选取 ${\psi^n=(\psi^n_1,\psi^n_2,\psi^n_3)}^{\infty}_{n=1} \in{\cal B}$, 则存在 $C>0$ 使得

我们用 $(U^n_1,U^n_2)$ 表示具有初始条件 $(\psi^n_1,\psi^n_2,\psi^n_3)$ 的解. 接下来, 我们证明对固定的 $t>0,$${(U^n_1(\cdot,t)U^n_2 (\cdot,t)}^{\infty}_{n=1}$ 存在一个收敛子序列. 根据文献[11,引理 5.1], 对于 $(x,t)\in\overline{\Gamma}\times(0,+\infty),$$k=1,2$, 我们有 $0\le U^n_k(x,t)\le \widetilde{B}k$. 因此, ${U^n_k(\cdot,t)}^{\infty}_{n=1}$是一致有界的. 基于 $G_1(x,y,t)$ 的性质, 可知存在 $\varsigma>0$, 使得对于任意 $\xi>0$, 当 $|x_1-x_2|<\varsigma$ 且 $x_1,x_2\in\overline{\Gamma}$ 时, 有

因此, 对于 $|x_1-x_2|<\varsigma$, 有

继而可知 ${U^n_1(\cdot,t)}^{\infty}_{n=1}$ 具有等度连续性. 同样地, ${U^n_2(\cdot,t)}^{\infty}_{n=1}$ 也具有等度连续性. 综上, ${U^n_k(\cdot,t)}^{\infty}_{n=1}, k=1,2$ 有收敛子序列, 从而半流的紧性得证. 引理证毕.

${\bf引理 5.2}$ 如果假设 2.1 成立以及 $R_0>1$, 那么存在 $\xi_1>0$ 使得

其中

${\bf证}$ 对足够小的 $\xi_1>0$ 和足够大的 $g>0$, 根据 Kerin-Rutman 定理可知 $r(\widetilde{{\cal R}})>1$, 其中

这里 $\widetilde{U}_{10}(x)=(\Pi-\xi_1)\int_0^ge^{-\sigma b}\int_{\Gamma}G_1(x,y,b){\rm d}y{\rm d}b$, $\Omega_i(\lambda)=\exp\{-\int_0^b(\sigma+\eta_i(s)+\delta_i(s)+\lambda){\rm d}s\},$$i=2,3$. Set $\widetilde{R}^+_0=r(\widetilde{{\cal R}}^+)$, 其中

从而可得 $R^+_0=r(\widetilde{{\cal R}}^+)=r(\widetilde{{\cal R}})>1$. 进而可知 $R^+_0\zeta^+_0(\cdot)=\widetilde{{\cal R}}^+\zeta^+_0$, 其中 $\zeta^+_0\in{\mathbb M}^0_+\setminus\{0\}$. 假设引理 5.2 结论不成立, 那么存在足够大的 $T^*>0$ 使得 $U_2(x,t)\le\xi_1$, $t\ge T^*$. 于是, 可得 $U_1(x,t)=(\Pi-\xi_1)\int_0^{t-T^*}e^{-\sigma b}\int_{\Gamma}G_1(x,y,b){\rm d}y{\rm d}b$. 对 $(x,t)\in\overline{\Gamma}\times(0,+\infty)$, 我们有

设 ${\cal L}(U_2)(\lambda)=\int_0^{\infty}e^{-\lambda t}U_2(x,t){\rm d}t$ 并对式 (5.1) 两边同乘以 $\zeta^+_0(x)$ 可得

这与 $R^+_0>1$ 矛盾. 引理证毕.

${\bf定理 5.1}$ 如果假设 2.1 成立以及 $R_0>1$, 那么系统 (1.5) 是一致持久的, 即存在 $\xi=\xi(x,b)>0$ 使得 $\liminf\limits_{t\to\infty}u_2(x,b,t)>\xi$, $(\psi_1,\psi_2,\psi_3)\in {\cal D}_0$. 此外, 系统 (1.5) 至少存在一个空间依赖的正平衡态 $E^+=(U^+_1(x),U^+_2(x))$.

${\bf证}$ 类似文献[11,命题5.3] 中的证明过程, 我们可以断言存在 $\xi_0>0$ 使得

设 $\kappa(t)=(U_1(\cdot,t),u_2(\cdot,\cdot,t))$ 是算子 $\Sigma$ 使得 $\Sigma(\kappa(s),t)=\kappa(t_s)$, $(s,t)\in{\Bbb R} \times(0,+\infty)$ 的总轨迹. 那么有

对 $b>0,x\in\Gamma$, 定义$\rho_k:{\mathbb M}^0_+\times{\mathbb Q}^0_+\times{\mathbb Q}^0_+\to{\Bbb R} ^+, k=1,2$的算子如下

由于 $\lim\limits_{t\to\infty}\| U_2(\cdot,t)\| {{\mathbb M}}>\xi_0$, 因此$(\psi_1,\psi_2,\psi_3)\in{\cal D}0$的一致 $\rho_1-$持久性(见文献[26,定义3.1]) 成立. 根据文献[26,推论4.22], 可假设$\kappa(t)$ 是具有预紧缩值域的总轨迹, 且$\inf{t\in{\Bbb R} }\rho_1(\kappa(t))=\inf{t\to{\Bbb R} }\| U_2(\cdot,t)\| {{\mathbb M}}>0$, 因而可得 $\rho_2(\kappa(0))>0$. 根据文献[26,推论4.22],可知存在$\xi=\xi(x,b)>0$ 使得 $\inf{t\to+\infty}\rho_2(\kappa(t))>\xi$. 再根据 (5.2)式, 可知, 对于 $b>0$, 有 $\liminf\limits_{t\to+\infty}u_2(x,b,t)>\xi$.

接下来, 我们证明系统 (1.5) 存在空间依赖的正稳态 $E^+$. 为此, 我们首先定义 ${\mathbb M}^0_+$ 上的线性算子

然后,正稳态满足方程

从而可得

其中 Itô是关于 $d_1\Delta$ is $\{\chi_t\}$ 的扩散, $P^x$ 表示 $\chi_t$ 的概率分布的期望, $s\Gamma:=\inf\{t>0:\chi_t\notin\Gamma\}$. 定义非线性算子 $L^1_+(\Gamma)$

以及从 $L^1_+(\Gamma_{\varpi})$ 到 $L^1_+(\Gamma)$ 的算子

其中 $\phi\in L^1_+(\Gamma_{\varpi}),\ \bar\omega>0$. 定义如下从 $L^1_+(\Gamma_{\varpi})$ 到 $L^1_+(\Gamma_{\varpi})$ 的算子

其中 ${\Gamma_{\varpi}}$ 是 $\Gamma$ 的一组适当的子域, 满足 $d(\Gamma_{\varpi},\partial \Gamma)>0$. 对于 $\varpi_1<\varpi_2$, 有 $\Gamma_{\varpi_1}\subset\Gamma_{\varpi_2}$ 且 $\lim\limits_{\varpi\to 0}\Big|\Gamma\setminus\Gamma_{\varpi}\Big|{\Gamma{\varpi}}=0$, 其中 $\Big|{\Gamma{\varpi}}$ 表示域在 $\Gamma_{\varpi}$ 上的约束. 注意到, 当 ${\cal A}{\varpi}$ 在 $L^1+(\Gamma_{\varpi})\setminus{0}$ 中具有正定点 $\phi^+{\varpi}$ 时, 则 ${\cal A}$ 有一个正定点.

一方面, 由于 $R_0=r({\cal R})>1$, 可知显然存在一个足够小的常数 $\varpi>0$ 使得 $r({\cal R}{\varpi})>1$. 对于固定的 $\varpi$, 可知 ${\cal A}{\varpi}(L^1+(\gamma_{\varpi})\subset L^1_+(\Gamma_{\varpi})$ 和 ${\cal A}{\varpi}(0)=0$. 另一方面, 根据文献[17,引理 7.1], 可知 ${\cal A}'{\varpi}(0)={\cal R}{\varpi}$, ${\cal A}'{\varpi}(\infty)=0$ 相对于 $L^1_+(\Gamma_{\varpi})$. 此外, ${\cal A}'{\varpi}(\infty)$ 的谱是 0, 并位于以 $0$ 为中心的半径小于 $1$ 的圆内. 事实上, ${\cal A}'{\varpi}(0)={\cal R}{\varpi}$ 有一个正特征函数 $\vartheta_{\varpi}\in L^1_+(\Gamma_{\varpi})\setminus{0}$ 使得 ${\cal A}'(0)\vartheta_{\varpi}={\cal R}{\varpi}\vartheta{\varpi}=r({\cal R}{\varpi})\vartheta{\varpi}$. 根据Krein-Rutman定理, 由于 $r({\cal R}{\varpi})>1$, ${\cal A}'(0)$ 没有任何特征函数对应特征值1. 此外, ${\cal A}{\varpi}$ 的紧性得证 (详见文献[17,引理7.1]). 应用 Krasnoselskii 不动点定理, 可得 ${\cal A}{\varpi}$ 至少有一个非平凡不动点 $\phi^+{\varpi}\in L^1_+(\Gamma_{\varpi})\setminus{0}$, 这意味着系统 (1.5) 至少有一个空间依赖的正稳定态解 $E^+=(U^+_1,U^+_2)$. 证毕.

在本节中, 我们根据中国50岁及以上感染艾滋病人群的报告病例, 对系统 (1.5) 进行数值模拟. 与文献[17]类似, 我们考虑矩形域 $\Gamma=(0,w_1)\times(0,w_2)\subset {\Bbb R} ^2, w_1=m,w_2=1/m, m>0$. 因此, 格林函数 $G_i(x,y,t)$ 形式为

其中 $x_1,y_1\in[w_1]$, $x_2,y_2\in[w_2], i=2,3$. 根据文献[5-6], 我们假设模型 (1.5) 参数值与初值如下

其中 $i=2,3$, $U_2(0)=2318,U_3(0)=1128$, $\Gamma_0=(0.2w_1,0.8w_1)\times(0.2w_2,0.8w_2)\subset \Gamma$.

利用文献[27,3.1.2 节]中的 Fredholm 离散化方法, 我们可以得到 $R_0$ 随迭代次数 $m\in(0,40]$ 的变化图.当$\gamma_1=5.7, \gamma_2=6.8$ 和 $\gamma_1=7.6, \gamma_2=9$ 时 $R_0$ 的值变化分别在图1(a)和(b)中给出. 这些图表明: 基本再生数 $R_0$ 随着 $m$ 的增加先增加后减少, 并且 $R_0$的值在图 1(a)和(b)中分别在 $m=10$ 时达到最大值 $0.927$ 和 $1.492$. 在图 2 中, 我们可以看到 当 $R_0<1$ 时$U_2(x,t)$ 密度分布曲面随着时间的推移逐渐收敛于零平面, 这与定理 4.2 的结果一致. 图 3 表明在 $R_0>1$ 时 $U_2(x,t)$ 密度分布曲面随着时间的推移远离零平面, 验证了定理 5.1 的结果. 图 4 比较了扩散率 $d_1$, $d_2$ 和 $d_3$ 对 $R_0$ 的影响. 从中我们可以看到, 当这些扩散率变小时 $R_0$ 变大. 此外,$R_0$ 对$d_1$ 最敏感, 对 $d_3$ 最不敏感.

在本文中, 我们建立了一类具有齐次 Dirichlet 边界条件的年龄空间结构 HIV/AIDS 传染病模型. 这种边界条件给动力学分析带来了挑战. 首先, 由于系统的无病稳态是非常数的, 即使模型参数是空间无关的, 这与齐次Neumann 边界条件下的空间同质模型的结论不同. 其次, 由于反应扩散方程和积分方程的耦合系统的耗散性不能直接得到, 我们证明了系统的一致持久性之后, 不能直接获得依赖于空间的正稳态解的存在性. 为了克服这些困难, 我们首先研究了系统的适定性, 并导出了基本再生数 $R_0$ 的泛函表达式. 然后我们得到了无病稳态的全局吸引性, 并证明了系统的一致持久性. 最后, 通过使用强 Fréchet 导数、紧致性理论和 Krasnoselskii 不动点定理证明了空间依赖的正稳态的存在性. 数值模拟部分, 我们在二维空间区域中通过数值模拟验证了文中理论结果. 值得指出的是, 本文所用的方法也可推广到研究其他具有年龄空间结构和齐次 Dirichlet 边界条件的传染病模型中.