数学物理学报, 2023, 43(6): 1744-1758

一类带漂移-扩散项的非线性抛物型方程组的能控性

张志朋,, 张亮,*

武汉理工大学理学院数学系 武汉 430070

Controllability of a Nonlinear Parabolic Systems with Drift-Diffusion Term

Zhang Zhipeng,, Zhang Liang,*

Department of Mathematics, Wuhan University of Technology, Wuhan 430070

通讯作者: *张亮，E-mail: zhangl@whut.edu.cn

收稿日期: 2023-03-9   修回日期: 2023-08-16  

基金资助: 国家自然科学基金(61573012)

Received: 2023-03-9   Revised: 2023-08-16  

Fund supported: NSFC(61573012)

作者简介 About authors

张志朋，E-mail:1005829259@qq.com

摘要

该文研究一类非线性抛物型方程组的能控性问题及其时间最优控制的存在性问题. 该方程具有漂移-扩散项和非线性项 $h(u,v)$, $g(u,v)$. 文中利用半群的 $L^{p}$-$L^{q}$ 估计和近年来发展的最大正则化研究方程组的正则性和控制函数的估计, 利用不动点定理证明了非线性系统的局部精确能控性, 并将该结论应用于证明时间最优控制的存在性.

关键词： 抛物型方程组; 能观性估计; 精确能控性; Kakutani不动点

Abstract

This paper studies the controllability and the existence of time optimal control for a class of nonlinear parabolic equations, which has drift-diffusion term with nonlinear terms $h(u,v)$ and $g(u,v)$. It applies the $L^{p}$-$L^{q}$ estimate of semigroups and the recently developed maximal regularity theory to study the regularity and the cost estimate of control functions. Moreover, this paper establishes the local exact controllability of the nonlinear control system by utilizing the Kakutani's fixed point theorem. Therefore, it is applied to the existence of time optimal controls.

Keywords： Parabolic equations; Observability estimate; Exact controllability; Kakutani's fixed point

PDF (585KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

张志朋, 张亮. 一类带漂移-扩散项的非线性抛物型方程组的能控性[J]. 数学物理学报, 2023, 43(6): 1744-1758

Zhang Zhipeng, Zhang Liang. Controllability of a Nonlinear Parabolic Systems with Drift-Diffusion Term[J]. Acta Mathematica Scientia, 2023, 43(6): 1744-1758

1 引言

设 $\Omega\subset\mathbb{R}^{3}$ 为具有光滑边界 $\partial\Omega$ 的有界区域, $\omega$ 为 $\Omega$ 的非空开子集. 令 $T>0$, 记 $Q_{T}=\Omega\times(0,T)$, $\Sigma_{T}=\partial\Omega\times(0,T)$. 函数空间 $L^{p}(\Omega)$ 和 $L^{p}(Q_{T})\ (1\leq p\leq\infty)$ 的范数分别为 $|\cdot|_{p}$ 考虑非线性抛物型方程组

$\begin{cases}u_{t}=\nabla\cdot(\nabla u-k_{0}u\nabla v)+g(u,v)+1_{\omega}f & \mbox{在 $Q_{T}$中},\\v_{t}=\triangle v-h(u,v) & \mbox{在 $Q_{T}$ 中},\\\partial_{\nu}u=\partial_{\nu}v=0 & \mbox{在 $\Sigma_{T}$ 上},\\u(x,0)=u_{0}(x),\ v(x,0)=v_{0}(x) & x\in\Omega,\end{cases}$

其中 $u_{t}=\partial u/ \partial t$, $v_{t}=\partial v/\partial t$. $\partial_{\nu}=\partial/\partial\nu$ 表示边界 $\partial \Omega$ 的外法向量导数. $1_{\omega}$ 指 $\omega$ 上的特征函数, $f=f(x,t)$ 为控制函数, $u=u(x,t)$ 为细菌密度, $v=v(x,t)$ 为催化剂浓度, $u_{0}(x)$ 和 $v_{0} (x)$ 表示初始值, $k_{0}>0$ 为趋化作用系数. 这里 $h(u,v)$ 和 $g(u,v)$ 都是 $\mathbb{R} \times\mathbb{R}\rightarrow\mathbb{R}$ 的连续可微有界函数, $g(0,0)=h(0,0)=0.$ 为简化符号, 本文省略函数表达式中的 $x$ 和 $t$.

函数对 $(u,v)$ 称为方程组 (1.1) 的弱解, 如果 $u,v\in L^{2}(0,T;H^{1}(\Omega)),$ $u_{t},v_{t}\in L^{2}(0,T;$ $H^{1} (\Omega)^{\ast}),$ $h(u,v),g(u,v)\in L^{2}(Q_{T}),$ $u\nabla v\in L^{2}(Q_{T})^{3},$ 在 $\Sigma_{T}$ 上 $\partial_{\nu}u=0,\partial_{\nu}v=0,$ 且

$\int_{0}^{T}\langle u_{t},\psi\rangle {\rm d}t=-\iint_{Q_{T}}(\nabla u-k_{0}u\nabla v)\nabla\psi {\rm d}x{\rm d}t+\iint_{Q_{T}}[g(u,v)+1_{\omega}f]\psi {\rm d}x{\rm d}t,$
$\int_{0}^{T}\langle v_{t},\phi\rangle {\rm d}t=-\iint_{Q_{T}}\nabla v\nabla\phi {\rm d}x{\rm d}t-\iint_{Q_{T}}h(u,v)\phi {\rm d}x{\rm d}t,$

对所有的 $\phi,\psi\in L^{2}(0,T;H^{1}(\Omega))$ 成立, 这里 $\left\langle \cdot,\cdot\right\rangle$ 表示 $H^{1}(\Omega)^{\ast}$ 与 $H^{1}(\Omega)$ 之间的对偶.

当方程组 (1.1)中控制项 $f=0$ 时可得

$\begin{cases} u_{t}=\nabla\cdot(\nabla v-k_{0}u\nabla v)+g(u,v) \mbox{在 $Q_{T}$ 中},\\ v_{t}=\triangle v-h(u,v) & \mbox{在 $Q_{T}$ 中},\\ \partial_{\nu}u=\partial_{\nu}v=0 & \mbox{在 $\Sigma_{T}$ 上},\\ u(0,x)=u_{0}(x),\ \ v(0,x)=v_{0}(x) & x\in\Omega.\end{cases}$

非线性抛物型偏微分方程组 (1.2) 由文献 [1] 提出, 是较一般形式的 Keller-Segel 方程组, 旨在描述在催化剂的趋化作用下细菌向液体与空气接触面的游动和在细菌重力的作用下液体产生流动这一生物过程. 以 Keller-Segel 抛物方程为代表的具有漂移-扩散项的非线性偏微分方程组是近些年研究的热点问题之一. 该类方程在生物医学上有着重要的应用, 而其特殊的数学结构使得其在偏微分方程理论上的研究具有一定的难度. 其解的全局存在性因结构或初始值的变化出现差异[2]; 解的性态会出现一些比较有趣的现象, 例如会出现爆破或聚集等现象等[3]. 关于方程组 (1.2) 在偏微分方程理论领域的研究比较丰富. 如在有界区域 $\Omega \subset\mathbb{R}^{N}$ $(N=2,3)$ 上, 当 $h(u,v)=-up(v)$时, ${\normalsize p(v)}$ 为连续可微的单调递增函数且满足 $p(0)=0$, 文献 [4] 得到方程组 (1.2) 局部解的存在性. 若 $N=3$, $p\in C^{1}([T]),\ p(0)=0;$ 在 $[T]$ 内 $p>0;$ 在 $(0,T)$ 上 $p^{\prime}>0,\ p^{\prime\prime}\leq0$, 文献 [5] 给出方程组 (1.2) 在带光滑边界的有界凸区域上弱解的存在性. 更多的关于方程组 (1.2) 在偏微分方程和应用领域中的结论可参见文献 [6,7] 等.

本文将讨论非线性控制系统 (1.1) 的能控性及其在时间最优控制的存在性问题中的应用. 设 $(\overline{u},\overline{v})$ 是方程组 (1.2) 关于初值 $(\overline{u}_{0},\overline{v}_{0})$ 的轨迹, 称非线性控制系统 (1.1) 关于轨迹 $(\overline{u},\overline{v})$ 在时刻 $T$ 处局部精确能控, 如果存在 $(\overline{u}_{0},\overline{v}_{0})$ 的邻域 $\vartheta$, 对任一初值 $(u_{0},v_{0})\in\vartheta$, 存在函数 $f$ 使控制系统 (1.1) 解 $(u,v)$ 满足 $(u,v)(x,T)=(\overline{u},\overline{v})(x,T),$ a.e. $x\in\Omega.$ 记集合 $\Pi=\{T>0;$ 存在控制 $f\in L^{\infty}(Q_{T}),$ 使方程组 (1.1) 的解 $(u,v)$ 在时刻 $T$ 处满足 $(u,v)(x,T)=(0,0),$ a.e. $x\in\Omega\}$. 考虑时间最优控制问题

$\mbox{inf}\ \Pi. $

关于非线性抛物型方程组 (1.1) 在控制理论领域的相关成果较少. 据我们所知, 文献 [8] 为第一篇关于此类方程组在控制理论领域的论文, 该文首次研究了经典的 Keller-Segel 方程组 (即 $g(u,v)=0,h(u,v)=-\gamma v+\delta u)$ 的最优控制问题, 并给出最优控制的存在性和必要条件. 在文献 [9] 中, 作者给出抛物-椭圆型 Keller-Segel 方程组 (即 $v_{t} =0,g(u,v)=0,h(u,v)=-\gamma v+\delta u$) 的局部精确能控性; 其后在文献 [10] 中, 作者研究了如下抛物型 Keller-Segel 方程组的局部精确能控性问题

$\begin{cases} u_{t}=\nabla\cdot(\nabla u-\chi u\nabla v)+1_{w}f & \mbox{在 $Q_{T}$中},\,\\ v_{t}=\Delta v-\gamma v+\delta u & \mbox{在 $Q_{T}$ 中},\\ \partial_{\nu}u=\partial_{\nu}v=0 & \mbox{在 $\Sigma_{T}$ 上},\\ u(0,x)=u_{0}(x),\ v(0,x)=v_{0}(x) & x\in\Omega. \end{cases}$

应该指出, 文献 [9,10] 的研究方法不同. 抛物-椭圆型 Keller-Segel 方程组可以视为单个非局部抛物型方程, 因此抛物-椭圆型 Keller-Segel 方程组可转化为单个非线性抛物型方程, 进而其能控性可采用单个非线性抛物型方程的能控性的处理方法; 但是抛物型 Keller-Segel 方程组 (1.4) 与抛物-椭圆型 Keller-Segel 方程组相比有本质的区别, 不能将之转化为单个非线性抛物型方程. 此外, 控制系统 (1.4) 为单个控制问题, 即控制仅施加于其中一个状态方程, 此时控制系统的能控性在建立能观性时存在一定的技术困难. 事实上, 具有单个控制的抛物型方程组的能控性问题目前仍属公开问题, 相关的综述性工作可参见文献 [11]; 关于非线性抛物型偏微分方程的能控性结论可参见文献 [12,13] 等.

本文为文献 [10] 的后续工作, 这里控制系统 (1.1) 具有非线性项 $h(u,v)$ 和 $g(u,v).$ 首先, 从偏微分方程理论角度来看, 方程组 (1.1) 与抛物-椭圆型偏微分方程组有本质的不同, 不属于非局部问题, 因此不能采用文献 [9] 的处理方法. 其次, 方程组 (1.1) 相比于 (1.4) 而言更具一般性, 其中的非线性项 $h(u,v)$ 和 $h(u,v)$ 的存在导致方程组的解的探讨具有一定的难度. 本文没有直接探讨解的全局或局部存在性问题, 方程组 (1.1) 的解可能会在有限时间爆破; 这里从能控性的角度探讨解的性态, 能控性结论不仅说明了方程组 (1.1) 解的存在性, 而且证明在控制的作用下可以保证方程组 (1.1) 的解不会在有限时间发生爆破现象. 再次, 控制系统 (1.1) 为具有单个控制的抛物型方程组的能控性问题, 该类问题目前尚未完全解决, 为研究的热点问题之一. 最后, 关于控制系统 (1.1) 的时间最优控制问题, 目前还没有相关的结果, 本文仅证明了时间最优控制的存在性, 时间最优控制的必要条件是值得继续探讨的问题. 本文拟研究线性偏微分方程组的正则性和零能控性, 运用不动点定理推导非线性抛物型偏微分方程组的局部能控性, 并给出该结论在时间最优控制的存在性问题中的应用. 本文所用的方法可用于研究其他的具有漂移-扩散项的非线性抛物型方程组的能控性.

2 线性抛物方程的正则性

考虑线性抛物型方程组

$\begin{cases}U_{t}=\triangle U-\nabla\cdot(B_{1}\nabla V)-\nabla\cdot(B_{2}U)+a_{1}U+a_{2}V+F & \mbox{在 $Q_{T}$ 中},\\V_{t}=\triangle V-c_{1}U-c_{2}V & \mbox{在 $Q_{T}$ 中},\\\partial_{\nu}U=\partial_{\nu}V=0 & \mbox{在 $\Sigma_{T}$ 上},\\U(x,0)=U_{0}(x),\ V(x,0)=V_{0}(x) & x\in\Omega,\end{cases}$

其中 $a_{i}=a_{i}(x,t),B_{i}=B_{i}(x,t),c_{i}=c_{i} (x,t),i=1,2,F=F(x,t).$

命题 2.1 设 $a_{1}$, $a_{2}$, $B_{1}$, $c_{1},c_{2}\in L^{\infty}(Q_{T})$, $B_{2}\in L^{\infty}(Q_{T})^{3}$, 在 $\Sigma_{T}$ 上 $B_{2}\cdot\nu=0$.

$(1)$ 若 $U_{0},V_{0}\in L^{2}(\Omega)$, $F\in L^{2}(Q_{T})$, 则方程组 $(2.1)$ 有唯一解 $(U,V)\in(L^{2}(0,T;H^{1}(\Omega)) \cap W^{1,2}([T];H^{1}(\Omega)^{\ast}))^{2}$, 且

$\Vert(U,V)\Vert_{(L^{2}(0,T;H^{1}(\Omega))\cap W^{1,2}([T];H^{1} (\Omega)^{\ast}))^{2}}^{2}\leq {\rm e}^{C\varpi T}(|U_{0}|_{2}^{2}+|V_{0}|_{2} ^{2}+\Vert F\Vert_{2}^{2}).$

$(2)$ 当 $2\leq p<\infty$ 时, 若 $F\in L^{p}(Q_{T})$, $U_{0}\in L^{p}(\Omega),\ V_{0}\in W^{2(1-\frac{1}{p}),p}(\Omega), $ 在 $\partial\Omega$ 上 $\partial_{\nu}V_{0}=0$, 则方程组 $(2.1)$ 有唯一解 $(U,V)\in L^{p}(Q_{T})\times W_{p}^{2,1} (Q_{T})$, 且

$\Vert U\Vert_{p}^{p}+\Vert V\Vert_{W_{p}^{2,1}(Q_{T})}^{p}\leq {\rm e}^{C\varpi T}(|U_{0}|_{p}^{p}+\Vert V_{0}\Vert_{W^{2(1-\frac{1}{p}),p}(\Omega)}^{p}+\Vert F\Vert_{p}^{p})$.

$(3)$ 当 $p>5$ 时, 若 $F\in L^{\infty}(Q_{T})$, $U_{0}\in L^{\infty}(\Omega),\ V_{0}\in W^{1,p}(\Omega), $ 在 $\partial\Omega$ 上 $\partial_{\nu}V_{0}=0$, 则方程组 (2.1) 有唯一解 $(U,V)\in L^{\infty}(Q_{T})^{2}$, 且

$\Vert U\Vert_{\infty}+\Vert V\Vert_{\infty}\leq {\rm e}^{C\varpi T}(|U_{0}|_{\infty}+\Vert V_{0}\Vert_{W^{1,p}(\Omega)}+\Vert F\Vert_{\infty}),$

其中 $C=C(\Omega),$

$\varpi =1+\Vert a_{1}\Vert_{\infty}+\Vert a_{2}\Vert_{\infty}+\Vert a_{1}\Vert_{\infty}^{2}+\Vert a_{2}\Vert_{\infty}^{2}+\left\Vert B_{1}\right\Vert _{\infty}+\left\Vert B_{2}\right\Vert _{\infty}\nonumber\\ +\Vert B_{1}\Vert_{\infty}^{2}+\Vert B_{2}\Vert_{\infty}^{2}+\Vert c_{1}|\Vert_{\infty}+\Vert c_{2}\Vert_{\infty}+\Vert c_{1}\Vert_{\infty} ^{2}+\Vert c_{2}\Vert_{\infty}^{2}.$

证 线性方程组 (2.1) 弱解的存在性和不等式 (2.2) 的建立是经典的, 可参见文献 [14]. 故只证不等式 (2.3) 和 (2.4). 方程组 (2.1) 的前两式两端分别乘以 $|U|^{p-2}U$ 和 $|V|^{p-2}V$, 并在 $\Omega$ 上积分可得

$\begin{matrix} \frac{\rm d}{{\rm d}t}|U|_{p}^{p}+\int_{\Omega}|\nabla U|^{2}|U|^{p-2}{\rm d}x\leq\ ; C(1+\Vert B_{1}\Vert_{\infty}^{2}+\Vert B_{2}\Vert_{\infty}^{2}+\Vert a_{1}\Vert _{\infty}+\Vert a_{2}\Vert_{\infty})\nonumber\\ \cdot(|\nabla V|_{p}^{p}+|U|_{p} ^{p}+|V|_{p}^{p})+C|F|_{p}^{p}, \end{matrix}$
$\begin{matrix} \frac{\rm d}{{\rm d}t}|V|_{p}^{p}+\int_{\Omega}|\nabla V|^{2}|V|^{p-2}{\rm d}x\leq\ C(\Vert c_{1}\Vert_{\infty}+\left\Vert c_{2}\right\Vert _{\infty})(|U|_{p}^{p}+|V|_{p}^{p}).\end{matrix}$

对 $|\nabla V|_{p}^{p}$ 关于 $t$ 微分, 再利用 (2.1) 式的第二个方程有

$\frac{\rm d}{{\rm d}t}|\nabla V|_{p}^{p}+\int_{\Omega}|\nabla V|^{p-2}|\triangle V|^{2}{\rm d}x\leq C(\Vert c_{1}\Vert_{\infty}^{2}+\left\Vert c_{2}\right\Vert _{\infty}^{2})(|\nabla V|_{p}^{p}+|U|_{p}^{p}+|V|_{p}^{p}),$

于是联立不等式 (2.6)-(2.8) 并利用 Gronwall 不等式得到

$|U(\cdot,t)|_{p}^{p}+|V(\cdot,t)|_{p}^{p}+|\nabla V(\cdot,t)|_{p}^{p}\leq {\rm e}^{C\varpi T}(|U_{0}|_{p}^{p}+|\nabla V_{0}|_{p}^{p}+\left\Vert F\right\Vert_{p}^{p})$

此外, 对 (2.1) 式的第二式应用发展方程的最大正则性[15]得

$\Vert V_{t}\Vert_{p}^{p}+\Vert\triangle V\Vert_{p}^{p}+\Vert V\Vert_{p}^{p}\leq C(\left\Vert c_{1}\right\Vert _{\infty}\Vert U\Vert_{p}^{p}+\left\Vert c_{2}\right\Vert _{\infty}\Vert V\Vert_{p}^{p}+\Vert V_{0}\Vert_{W^{2(1-\frac{1}{p}),p}(\Omega)}^{p}),$

其中常数 $C$ 与 $T$ 无关. 于是, 不等式 (2.3) 由不等式 (2.9) 和 (2.10) 可得. 下证不等式 (2.4). 令参数 $0<\delta<1$, 设 $A$ 为扇形算子

$Au:=-\triangle u,\ \forall u\in D(A):=\left\{ u\in W^{2,p}(\Omega);\ \left. \partial_{\nu}u\right\vert _{\partial\Omega}=0\right\},$

则 $-(A+\delta)$ 和 $-A$ 在 $L^{p}(\Omega)(1<p<\infty)$ 上分别生成解析半群 $\{{\rm e}^{-t(A+\delta )}\}_{t\geq0}$ 和 $\{{\rm e}^{-tA}\}_{t\geq0}.$ 令 $\alpha>0,$ 设 $D(A^{\alpha})$ 为 $L^{p}(\Omega)$ 中赋以图范数的 Banach 空间, 则有如下嵌入定理[16]

$D((A+\delta)^{\alpha})\hookrightarrow W^{1,p}(\Omega),\ \ \alpha>\frac{1}{2},$
$D((A+\delta)^{\alpha})\hookrightarrow C^{\gamma}(\overline{\Omega}),\ \ 0\leq\gamma<2\alpha-\frac{3}{p}. $

对任意的 $\varepsilon>0,1<p\leq q<\infty,$ 半群 $\{{\rm e}^{-tA}\}_{t\geq0}$ 满足 $L^{p}$-$L^{q}$ 估计

$|(A+\delta)^{\alpha}{\rm e}^{-tA}\nabla\cdot u|_{p}\leq C_{\varepsilon}t^{-\alpha-\frac{1}{2}-\varepsilon}|u|_{p},$
$|(A+\delta)^{\alpha}{\rm e}^{-t(A+\delta)}u|_{q}\leq Ct^{-\frac{3}{2}(\frac{1}{p}-\frac{1}{q})-\alpha}|u|_{p}, \forall u\in L^{p}(\Omega),t>0.$

设 $U_{0}\in C(\overline{\Omega}),V_{0}\in C^{1}(\overline{\Omega})$, $F\in C(\overline{Q}_{T})$. 用半群 $\{{\rm e}^{-tA}\}_{t\geq0}$ 将 (2.1) 式中第一个方程的解写为

$U(\cdot,t)={\rm e}^{-tA}U_{0}+\int_{0}^{t}{\rm e}^{-(t-s)A}[-\nabla\cdot(B_{1}\nabla V)-\nabla\cdot(B_{2}U)+a_{1}U+a_{2}V+F](\cdot,s){\rm d}s.$

上式两端取 $C(\overline{\Omega})$ 范数可得: $\forall t\in\lbrack0,T],$

$\Vert U(\cdot,t)\Vert_{C(\overline{\Omega})} \leq\ \int_{0}^{t}\Vert {\rm e}^{-(t-s)A}[-\nabla\cdot(B_{1}\nabla V)-\nabla\cdot(B_{2}U)+a_{1}U+a_{2}V](\cdot,s)\Vert_{C(\overline{\Omega})}{\rm d}s\nonumber\\ +\Vert {\rm e}^{-tA}U_{0}\Vert_{C(\overline{\Omega})}+\int_{0}^{t}\Vert {\rm e}^{-(t-s)A}F(\cdot,s)\Vert_{C(\overline{\Omega})}{\rm d}s.$

下面估计 (2.12) 式中各项. 注意 $-A$ 在 $C(\overline{\Omega})$ 上生成压缩的解析半群, 则对 $0\leq s<t\leq T$ 有

$\Vert {\rm e}^{-tA}U_{0}\Vert_{C(\overline{\Omega})}\leq\Vert U_{0}\Vert _{C(\overline{\Omega})},\ \ \Vert {\rm e}^{-(t-s)A}F(\cdot,s)\Vert_{C(\overline {\Omega})}\leq\Vert F(\cdot,s)\Vert_{C(\overline{\Omega})}.$

取 $\varepsilon$ 和 $\alpha$ 满足 $ 0<\varepsilon<\frac{p-5}{2p},\frac{3}{2p}<\alpha<\frac{1}{2}-\frac{1}{p}-\varepsilon,$ 由此根据嵌入 (2.15) 式和不等式 (2.15a) 并应用 Hölder 不等式有

$ \int_{0}^{t}\left\Vert {\rm e}^{-(t-s)A}\nabla\cdot(B_{1}\nabla V)(\cdot,s)\right\Vert_{C(\overline{\Omega})}{\rm d}s \leq C\int_{0}^{t}\left\vert (A+\delta)^{\alpha}{\rm e}^{-(t-s)A}\nabla \cdot(B_{1}\nabla V)(\cdot,s)\right\vert _{p}{\rm d}s\\ \leq C\int_{0}^{t}(t-s)^{-\alpha-\frac{1}{2}-\varepsilon}\left\vert (B_{1}\nabla V)(\cdot,s)\right\vert _{p}{\rm d}s\\ \leq C\left\Vert B_{1}\right\Vert _{\infty}\left\Vert \nabla V\right\Vert _{p}T^{\frac{1}{2}-\alpha-\varepsilon-\frac{1}{p}},\forall t\in\lbrack0,T],$

和

$ \int_{0}^{t}\left\Vert {\rm e}^{-A(t-s)}(a_{1}U)(\cdot,s)\right\Vert _{C(\overline{\Omega})}{\rm d}s \leq C\int_{0}^{t}|(A+\delta)^{\alpha}{\rm e}^{-(t-s)(A+\delta)}{\rm e}^{\delta (t-s)}(a_{1}U)(\cdot,s)|_{p}{\rm d}s\\ \leq C\int_{0}^{T}(t-s)^{-\alpha}{\rm e}^{\delta(t-s)}|(a_{1}U)(\cdot,s)|_{p}{\rm d}s\\ \leq C{\rm e}^{\delta T}\left\Vert a_{1}\right\Vert _{\infty}\left\Vert U\right\Vert _{p}T^{-\alpha+1}.$

类似估计 (2.16) 式中其余各项并联立 (2.17) 式可得

$\left\Vert U(\cdot,t)\right\Vert _{C(\overline{\Omega})}\leq C\varpi {\rm e}^{T}(\left\Vert U\right\Vert _{p}+\left\Vert V\right\Vert _{p}+\left\Vert \nabla V\right\Vert _{p}),\forall t\in\lbrack0,T],$

由此式以及 (2.9) 式得

$\left\Vert U(\cdot,t)\right\Vert _{C(\overline{\Omega})}\leq {\rm e}^{C\varpi T}(|U_{0}|_{p}+|\nabla V_{0}|_{p}+\left\Vert F\right\Vert _{p}),\forall t\in\lbrack0,T].$

应用半群 $\{{\rm e}^{-t(A+\delta)}\}_{\geq0}$ 于 (2.1) 式中第二个方程有表达式

$V(\cdot,t)={\rm e}^{-t(A+\delta)}V_{0}+\int_{0}^{t}{\rm e}^{-(A+\delta)(t-s)}[-c_{1}U-(c_{2}-\delta)V)(\cdot,s){\rm d}s.$

上式两端取 $W^{1,p}(\Omega)$ 范数, 则对任意的 $0\leq t\leq T,$

$\Vert V(\cdot,t)\Vert_{W^{1,p}(\Omega)} \leq\Vert {\rm e}^{-t(A+\delta)}V_{0}\Vert_{W^{1,p}(\Omega)} \! +\!\int_{0}^{t}\Vert {\rm e}^{-(t-s)(A+\delta)}[-c_{1}U-(c_{2}-\delta )V](\cdot,s)\Vert_{W^{1,p}(\Omega)}{\rm d}s$

下面估计 (2.19) 式中各项. 首先, 类似于 (2.3) 式进行能量估计可得

$\Vert {\rm e}^{-tA}V_{0}\Vert_{W^{1,p}(\Omega)}\leq {\rm e}^{CT}\left\vert \nabla V_{0}\right\vert _{p}.$

其次, 根据嵌入 (2.12) 和不等式 (2.15) 知, 对 $t\in \lbrack0,T],$

$ \int_{0}^{t}\Vert {\rm e}^{-(t-s)(A+\delta)}(c_{1}U)(\cdot,s)\Vert_{W^{1,p}(\Omega)}{\rm d}s \leq C\int_{0}^{T}|(A+\delta)^{\alpha}{\rm e}^{-(t-s)(A+\delta)}(c_{1}U)(\cdot,s)|_{p}{\rm d}s\\ \leq C\int_{0}^{T}(t-s)^{-\alpha}|(c_{1}U)(\cdot,s)|_{p}{\rm d}s\\ \leq C\left\Vert c_{1}\right\Vert _{\infty}|U|_{p}T^{-\alpha+1}.$

类似的方式估计 (2.19) 式中其余各项, 可以得到

$\Vert V(\cdot,t)\Vert_{W^{1,p}(\Omega)}\leq {\rm e}^{C\varpi T}(|U_{0}|_{p}+\left\vert \nabla V_{0}\right\vert _{p}+\left\Vert F\right\Vert_{p}),\forall t\in\lbrack0,T].$

联立不等式 (2.18), (2.21), 以及 Sobolev 嵌入 $W^{1,p}(\Omega)\hookrightarrow C(\overline{\Omega})\ (p>3)$, 有

$\left\Vert U\right\Vert _{\infty}+\Vert V\Vert_{\infty}\leq {\rm e}^{C\varpi T}(|U_{0}|_{\infty}+\Vert V_{0}\Vert_{W^{1,p}(\Omega)}+\Vert F\Vert_{\infty}).$

即证当 $U_{0}\in C(\overline{\Omega}),V_{0}\in C^{1}(\overline{\Omega})$ 和 $F\in C(\overline{Q}_{T})$ 时, 不等式 (2.4) 成立. 对 $U_{0}\in L^{\infty}(\Omega),V_{0}\in W^{1,p}(\Omega),$ $F\in L^{\infty}(Q_{T})$ 的情形, 可取光滑函数列 $\{(U_{0n},V_{0n})\}_{n=1}^{\infty}$ 和 $\{F_{n}\}_{n=1}^{\infty}$分别在 $L^{2}(\Omega)^{2}$ 和 $\ L^{2}(Q_{T})$ 中收敛于 $(U_{0},V_{0})$ 和 $F$. 由于方程组 (2.1) 相应于 $(U_{0n},V_{0n},F_{n})$ 的解$(U_{n},V_{n})$ 满足不等式 (2.4), 故可取 $\{(U_{n},V_{n})\}_{n=1}^{\infty}$ 的子列使其在 $L^{2}(Q_{T})$ 中收敛于 $(U,V),$ 从而 $(U,V)$ 为 (2.1) 式相应于$(U_{0},V_{0},F)$ 的弱解, 且满足不等式 (2.4).

3 线性抛物方程组的能控性

考虑线性方程组刻画的线性控制系统

$\begin{cases} U_{t}=\triangle U-\nabla\cdot(B_{1}\nabla V)-\nabla\cdot(B_{2}U)+a_{1}U+a_{2}V+1_{\omega}f & \mbox{在 $Q_{T} $ 中},\\ V_{t}=\triangle V-c_{1}U-c_{2}V & \mbox{在 $Q_{T} $ 中},\\ \partial_{\nu}U=\partial_{\nu}V=0 & \mbox{在 $\sum_{T}$ 上},\\ U(x,0)=U_{0}(x),\ \ V(x,0)=V_{0}(x) & x\in\Omega, \end{cases}$

其中系数 $a_{1},a_{2},c_{1},c_{2},B_{1}\in L^{\infty}(Q_{T}),B_{2}\in L^{\infty}(Q_{T})^{3};$ 在 $\Sigma_{T}$ 上 $B_{2}\cdot\nu=0;$ $c_{1}$ 满足

$c_{1}(x,t)\geq\sigma_{0}>0,\forall(x,t)\in Q_{T}.$

注 3.1 这里关于系数 $c_{1}$ 的条件 (3.2) 是合理的. 事实上, 如果 $c_{1}=0,$ 那么由于关于 $V$ 的状态方程无控制作用, $V$ 的轨迹无法在给定的时刻 $T$ 到达 $0,$ 亦即零能控不存在. 条件 3.2 可弱化为: 存在非空开子集 $\omega^{\prime}\subset\omega,$ 使得 $ c_{1}(x,t)\geq\sigma_{0}>0\text{ }(\text{或}\ c_{1}(x,t)\leq-\sigma _{0})\text{a.e. }(x,t)\in\omega^{\prime}\times(0,T).$

控制系统 3.1 的共轭线性方程组为

$\begin{cases} -\eta_{t}=\triangle\eta+B_{2}\cdot\nabla\eta+a_{1}\eta-c_{1}\zeta \mbox{在 $Q_{T}$ 中},\\ -\zeta_{t}=\triangle\zeta-\nabla\cdot(B_{1}\nabla\eta)+a_{2}\eta-c_{2}\zeta \mbox{在 $ Q_{T}$ 中},\\ \partial_{\nu}\eta=\partial_{\nu}\zeta=0 & \mbox{在 $\Sigma_{T}$ 上},\\ \eta(x,T)=\eta^{T}(x),\ \zeta(x,T)=\zeta^{T}(x) & x\in\Omega, \end{cases}$

其中 $\eta^{T}$, $\zeta^{T}\in L^{2}(\Omega)$. 本节将建立线性控制系统 3.1在单个控制 $f$下的能控性. 为此, 需建立合适的Carleman估计并由此推导所需的能观性估计. 设 $\omega_{0}$ 是 $\omega$ 的非空开子集, $\overline{\omega_{0}}\subset\omega,$ 则存在函数[12] $\beta\in C^{2}(\overline{\Omega})$, 满足 $\beta(x)>0,\ \forall x\in\Omega$; $\beta|_{\partial\Omega}=0$; $|\nabla\beta(x)|>0$, $\forall x\in\overline{\Omega\setminus\omega_{0}}$. 有如下 Carleman 不等式[10]

引理 3.1 设 $g_{i}=g_{i}(x,t)\in L^{2}(Q_{T}),i=0,1,\cdots,3,$ $z=z(x,t)$ 为方程

$\begin{cases} -z_{t}-\triangle z=g_{0}+\sum\limits_{i=1}^{N}\dfrac{\partial g_{i}}{\partial x_{i}} & \mbox{在 $Q_{T}$ 中},\\ \partial_{\nu}z=0 & \mbox{在 $\Sigma_{T}$ 上},\\ z(\cdot,T)=z_{T}\in L^{2}(\Omega) & \quad \end{cases}$

的弱解, 则对任意的实数 $d,$ 存在常数 $C=C(\Omega,\omega_{0},d)$ 和常数 $\lambda_{0}=\lambda_{0}(\Omega,\omega_{0},d)>1$, 当 $\lambda\geq\lambda_{0}$ 和 $s\geq\gamma(\lambda)(T+T^{2})$ 时有

$\iint_{Q_{T}}[\lambda^{2}(s\varphi)^{1+d}|\nabla z|^{2}+\lambda ^{4}(s\varphi)^{3+d}|z|^{2}]{\rm e}^{2s\alpha}{\rm d}x{\rm d}t\nonumber\\ \leq \ C\iint_{Q_{T}}(s\varphi)^{d}{\rm e}^{2s\alpha}|g_{0}|^{2}{\rm d}x{\rm d}t+C\sum_{i=1}^{3}\iint_{Q_{T}}\lambda^{2}(s\varphi)^{2+d}{\rm e}^{2s\alpha}|g_{i} |^{2}{\rm d}x{\rm d}t\nonumber\\ +C\iint_{\omega_{0}\times(0,T)}\lambda^{4}(s\varphi)^{3+d}{\rm e}^{2s\alpha }|z|^{2}{\rm d}x{\rm d}t, $

其中

$\varphi=\frac{{\rm e}^{\lambda\beta}}{t(T-t)},\ \alpha=\frac {{\rm e}^{\lambda\beta}-{\rm e}^{2\lambda\Vert\beta\Vert_{C(\overline{\Omega})}}} {t(T-t)},\ \gamma(\lambda)={\rm e}^{2\lambda\Vert\beta\Vert_{C(\overline{\Omega})}}. $

下面建立适合于方程组 (3.3) 的Carleman估计.

引理 3.2 设共轭方程组 (3.3) 中系数 $c_{1}$ 满足条件 (3.2), 则存在正常数 $\lambda_{1}=C_{0}(\Omega,\omega,\omega_{0})\varpi$, 满足 $\gamma(\lambda_{1})\geq\lambda_{1}>1$, 对于任意的 $\lambda \geq\lambda_{1}$, $s\geq\gamma(\lambda)(T+T^{2})$ 和 $\eta^{T}$, $\zeta^{T}\in L^{2}(\Omega)$, 方程组 (3.3) 的解 $(\eta,\zeta)$ 满足

$\Phi_{1}(\eta)+\Phi_{2}(\zeta)\leq C\iint_{\omega\times(0.T)}\lambda ^{17}(s\varphi)^{9}|\eta|^{2}{\rm e}^{2s\alpha}{\rm d}x{\rm d}t, $

其中常数 $C=C(\Omega,\omega,\omega_{0})$, 常数 $\varpi$ 由 (2.5) 给出,

$\Phi_{1}(\eta)=\iint_{Q_{T}}[\lambda^{2}(s\varphi)^{3}|\nabla\eta|^{2}+\lambda^{4}(s\varphi)^{5}|\eta|^{2}]{\rm e}^{2s\alpha}{\rm d}x{\rm d}t,$
$\Phi_{2}(\zeta)=\iint_{Q_{T}}[\lambda^{2}(s\varphi)|\nabla\zeta|^{2}+\lambda^{4}(s\varphi)^{3}|\zeta|^{2}]{\rm e}^{2s\alpha}{\rm d}x{\rm d}t.$

证 以下常数 $C_{i}$ 表示不同的依赖于 $\Omega,\omega$ 和 $\omega_{0}$ 的常数. 对方程组 (3.3) 的前两式分别应用引理 3.1, 其中参数 $d$ 分别取 $2$ 和 $0,$ 由此存在常数 $\widetilde{\lambda}_{1}$ 满足 $\gamma(\widetilde{\lambda}_{1})\geq\widetilde{\lambda}_{1}=C(\Omega,\omega_{0})\varpi>1,$ 使对所有的 $\lambda>\widetilde{\lambda}_{1}$ 和 $s\geq \gamma(\lambda)(T+T^{2})$, 方程组 (3.3)的解 $(\eta,\zeta)$ 满足不等式

$\Phi_{1}(\eta)+\Phi_{2}(\zeta)\leq C_{1}\iint_{\omega_{0}\times(0,T)}\lambda^{4}[(s\varphi)^{5}|\eta|^{2}+(s\varphi)^{3}|\zeta|^{2}]{\rm e}^{2s\alpha }{\rm d}x{\rm d}t.$

下面希望消除 (3.7) 式右端中关于 $\zeta$ 的项. 为此, 采用乘子方法. 设截断函数$\xi\in C_{0}^{\infty}(\Omega)$, 在 $\omega_{0}$ 上 $\xi=1$; 在 $\Omega\backslash\overline {\omega}$ 中 $\xi=0$; 在 $\Omega$ 中 $0\leq\xi\leq1$, 并且满足 $\xi^{-1/2}\triangle\xi\in L^{\infty}(\Omega)$ 和 $\xi^{-1/2}\nabla\xi\in L^{\infty}(\Omega)^{3}.$ 函数 $\xi$ 的存在性易证. 事实上, 可取截断函数$\psi\in C^\infty_0({\Omega})$, 在 $\omega_0$ 上 $\psi=1$; 在 $\Omega\setminus\overline{\omega}$ 上, $\psi=0$; 在 $\Omega$ 上 $0\leq\psi\leq 1$. 若取 $\xi=\psi^4$, 则 $\xi$ 即满足所需条件.

令 $\Theta=\lambda^{4}(s\varphi)^{3}{\rm e}^{2s\alpha}.$ 根据方程组 (3.3) 中第一式可知 $c_{1}\zeta=\eta_{t}+\triangle\eta+B_{2}\cdot\nabla\eta+a_{1}\eta$, 于是两边同乘 $\Theta\xi\zeta$, 并再 $Q_{T}$ 上积分可得

$\iint_{Q_{T}}c_{1}\Theta\xi|\zeta|^{2}{\rm d}x{\rm d}t=\iint_{Q_{T}}\Theta\xi\zeta \cdot\lbrack\eta_{t}+\triangle\eta+B_{2}\cdot\nabla\eta+a_{1}\eta]{\rm d}x{\rm d}t,$

上式应用分部积分并根据 (3.1) 中第二式有

$\iint_{Q_{T}}c_{1}\lambda^{4}(s\varphi)^{3}|\zeta|^{2}{\rm e}^{2s\alpha}\xi {\rm d}x{\rm d}t=J_{1}+J_{2}+J_{3},$

其中

$J_{1} =\iint_{Q_{T}}\zeta\lbrack-\Theta_{t}\xi\eta-\nabla(\Theta\xi)\cdot\nabla\eta+\Theta\xi B_{2}\cdot\nabla\eta+a_{1}\Theta\xi\eta]{\rm d}x{\rm d}t,\\ J_{2} =-\iint_{Q_{T}}\nabla\zeta\cdot\lbrack\nabla(\Theta\xi)\eta+2\Theta\xi\nabla\eta]{\rm d}x{\rm d}t,\\ J_{3} =\iint_{Q_{T}}[B_{1}\eta\nabla(\Theta\xi)\cdot\nabla\eta+\Theta\xi B_{1}|\nabla\eta|^{2}+a_{2}\Theta\xi|\eta|^{2}]{\rm d}x{\rm d}t.$

下面估计 $J_{1}$, $J_{2}$ 和 $J_{3}.$ 已知 $s$ 足够大时 $s\varphi\geq1,$ 且有估计

$|\Theta_{t}|\leq\lambda^{4}(s\varphi)^{5}{\rm e}^{2s\alpha},\ |\nabla(\Theta \xi)|\leq C\lambda^{4}(\xi^{\frac{1}{2}}s^{3}\varphi^{3}+\xi\lambda s^{4}\varphi^{4}){\rm e}^{2s\alpha},$

则利用 Cauchy 不等式可知对任意的正实数 $\varepsilon_{1},$

$J_{1} \leq\frac{\varepsilon_{1}}{2}\Phi_{2}(\zeta)+\frac{C_{2}\varpi }{2\varepsilon_{1}}\iint_{Q_{T}}[\lambda^{6}(s\varphi)^{5}|\nabla\eta |^{2}+\lambda^{4}(s\varphi)^{7}|\eta|^{2}]{\rm e}^{2s\alpha}\xi {\rm d}x{\rm d}t,$
$J_{2} \leq\frac{\varepsilon_{1}}{2}\Phi_{2}(\zeta)+\frac{C_{2}\varpi }{2\varepsilon_{1}}\iint_{Q_{T}}[\lambda^{7}(s\varphi)^{5}|\nabla\eta |^{2}+\lambda^{8}(s\varphi)^{7}|\eta|^{2}]{\rm e}^{2s\alpha}\xi {\rm d}x{\rm d}t,$
$J_{3} \leq\varepsilon_{1}\Phi_{1}(\eta)+\frac{C_{2}\varpi}{\varepsilon_{1} }\iint_{Q_{T}}[\lambda^{4}(s\varphi)^{3}|\nabla\eta|^{2}+\lambda^{8}(s\varphi)^{5}|\eta|^{2}]{\rm e}^{2s\alpha}\xi {\rm d}x{\rm d}t.$

下面将近一步消除(3.9)-(3.11) 式中关于 $|\nabla\eta|$ 的项. 为此, 令 $\Psi=\lambda^{7}(s\varphi)^{5}{\rm e}^{2s\alpha},$ 再次根据方程组 (3.3) 中第一式可知 $-\triangle\eta=\eta_{t}+B_{2}\cdot\nabla\eta+a_{1}\eta-c_{1}\zeta$, 此式乘以 $\Psi\xi\eta$ 并在 $Q_{T}$ 上积分可得

$\iint_{Q_{T}}\lambda^{7}(s\varphi)^{5}{\rm e}^{2s\alpha}\xi|\nabla\eta|^{2}{\rm d}x{\rm d}t\\ =\ \iint_{Q_{T}}[-\frac{1}{2}\Psi_{t}\xi|\eta|^{2}+\Psi\xi\eta(B_{2}\cdot\nabla\eta+a_{1}\eta-c_{1}\zeta)-\nabla\eta\cdot\nabla(\Psi\xi)\eta]{\rm d}x{\rm d}t.$

注意

$|\Psi_{t}|\leq\lambda^{7}(s\varphi)^{7}{\rm e}^{2s\alpha},\ |\nabla(\Psi\xi)|\leq C\lambda^{7}[\xi^{\frac{1}{2}}(s\varphi)^{5}+\xi\lambda(s\varphi )^{6}]{\rm e}^{2s\alpha},$

故利用 Cauchy 不等式对上面的等式做估计得到

$\iint_{Q_{T}}\lambda^{7}(s\varphi)^{5}{\rm e}^{2s\alpha}\xi|\nabla\eta|^{2}{\rm d}x{\rm d}t \leq\varepsilon_{2}[\Phi_{1}(\eta)+\Phi_{2}(\zeta)] +\frac{C_{3}\varpi}{\varepsilon_{2}}\iint_{Q_{T}}\lambda^{14}(s\varphi )^{9}{\rm e}^{2s\alpha}|\eta|^{2}\xi {\rm d}x{\rm d}t.$

因为系数 $c_{1}$ 满足条件 (3.2), 所以

$\iint_{\omega_{0}\times(0,T)}\lambda^{4}(s\varphi)^{3}|\zeta|^{2}{\rm e}^{2s\alpha }{\rm d}x{\rm d}t\leq\frac{1}{\sigma_{0}}\iint_{Q_{T}}c_{1}\Theta\xi|\zeta|^{2}{\rm d}x{\rm d}t,$

根据此 (3.8)-(3.12) 式, 并取

$\varepsilon_{1}=\frac{\sigma_{0}}{4C_{1}},\varepsilon_{2}=\frac{\sigma_{0}^{2}}{32C_{1}^{2}C_{2}\varpi},$

于是有不等式

$\Phi_{1}(\eta)+\Phi_{2}(\zeta)\leq C_{4}\varpi^{3}\iint_{\omega\times (0,T)}\lambda^{14}(s\varphi)^{9}{\rm e}^{2s\alpha}|\eta|^{2}{\rm d}x{\rm d}t,$

最后取 $\lambda_{1}=C(\Omega,\omega,\omega_{0})\varpi \geq\widetilde{\lambda}_{1}$ 足够大, 故当 $\lambda \geq\lambda_{1}$, $s\geq\gamma(\lambda)(T+T^{2})$ 时, 由 (3.13) 式即可得所需不等式 (3.6).

引理 3.3 设共轭方程组 (3.3) 中系数 $c_{1}$ 满足条件(3.2), 则存在正常数 $\lambda$ 和 $s$, 使对任意的 $T>0,\ \eta^{T}$, $\zeta^{T}\in L^{2}(\Omega)$, 方程组 (3.1) 的解 $(\eta,\zeta)$ 满足

$\int_{\Omega}[|\eta(\cdot,0)|^{2}+|\zeta(\cdot,0)|^{2}]{\rm d}x\leq {\rm e}^{C(\varpi T+\frac{1}{T}+1)}\iint_{\omega\times(0,T)}\lambda^{17}(s\varphi)^{9}{\rm e}^{2s\alpha}|\eta|^{2}{\rm d}x{\rm d}t. $

特别地, 当 $\lambda$ 和 $s$ 足够大时,

$\int_{\Omega}[|\eta(\cdot,0)|^{2}+|\zeta(\cdot,0)|^{2}]{\rm d}x\leq {\rm e}^{C(\varpi T+\frac{1}{T}+1)}\iint_{\omega\times(0,T)}{\rm e}^{\kappa\alpha}|\eta|^{2}{\rm d}x{\rm d}t,$

其中参数 $\kappa<2s,$ $\varpi$由 (2.5) 式给出.

证 首先根据能量估计可得对任意的 $t\in \lbrack0,T],$

$\int_{\Omega}|\eta(\cdot,0)|^{2}+|\zeta(\cdot,0)|^{2}{\rm d}x\leq {\rm e}^{C\varpi T} \int_{\Omega}|\eta(\cdot,t)|^{2}+|\zeta(\cdot,t)|^{2}{\rm d}x.$

上式两端取 $[T/4,3T/4]$ 上的积分, 有

$\int_{\Omega}|\eta(\cdot,0)|^{2}+|\zeta(\cdot,0)|^{2}{\rm d}x\leq\frac{2} {T}{\rm e}^{C\varpi T}\iint_{\Omega\times\lbrack\frac{T}{4},\frac{3T}{4}]} (|\eta|^{2}+|\zeta|^{2}){\rm d}x{\rm d}t.$

因为在 $\Omega\times\left[ T/4,3T/4\right] $ 上 $(s\varphi)^{-5}{\rm e}^{-2s\alpha}\leq {\rm e}^{Cs/T^{2}},\ (s\varphi)^{-3}{\rm e}^{-2s\alpha}\leq {\rm e}^{Cs/T^{2}},$ 所以上式表明

$\int_{\Omega}|\eta(\cdot,0)|^{2}+|\zeta(\cdot,0)|^{2}{\rm d}x\leq\frac{2} {T}{\rm e}^{C\varpi T+\frac{Cs}{T^{2}}}[\Phi_{1}(\eta)+\Phi_{2}(\zeta)].$

最后取 $\lambda=\lambda(\Omega,\omega,\varpi)>1$ 和 $s=\gamma(\lambda)(T+T^{2})$ 足够大使得

$\frac{2}{T}{\rm e}^{C\varpi T+\frac{Cs}{T^{2}}}\leq {\rm e}^{C(\varpi T+\frac{1}{T}+1)},$

故由引理 3.2 中 (3.6) 式以及 (3.16) 式可得 (3.14) 式; 进一步, 若取 $\lambda$ 和 $s$ 足够大时也有 (3.15) 式.

下面建立线性控制系统 (3.1) 的零能控性.

定理 3.1 设方程组 (3.1) 的系数 $c_{1}$ 满足条件 (3.2), 则对任意的初值 $(U_{0},V_{0})\in L^{2}(\Omega)^{2}$, 存在控制函数 $f\in L^{\infty}(Q_{T})$, 使得方程组 (3.1) 的解 $(U,V)$ 满足 $U,V\in L^{2}([0,T);H^{1} (\Omega))\cap C([0,T);L^{2}(\Omega),$ $U(x,T)=0, V(x,T)=0,\text{ a.e. }(x,t)\in Q_{T},$ 且控制函数 $f$ 满足

$\Vert f\Vert_{\infty}\leq {\rm e}^{C(\varpi T+\frac{1}{T}+1)}(|U_{0}|_{2}+|V_{0}|_{2}), $

这里的 $C$ 为仅依赖于 $\Omega$ 和 $\omega$ 的正常数, $\varpi$ 由 (2.5) 式给出.

证 设参数 $s$ 和 $\kappa<2s$ 由引理 3.3 给出. 对于任意给定的实数 $\varepsilon>0$, 考虑最优控制问题

$\text{min}\{J_{\varepsilon}(f),\ f\in L^{2}(Q_{T})\},$

其中 $J_{\varepsilon}(f)=\frac{1}{2}\Vert1_{\omega}f{\rm e}^{-\frac{1}{2}\kappa\alpha}\Vert_{2}^{2}+\frac{1}{2\varepsilon}\left\Vert (U,V)(\cdot,T)\right\Vert_{L^{2}(\Omega)^{2}}^{2},$ 函数对 $(U,V)$ 为线性方程组 (3.1) 相应于 $f$ 的解. 易知该最优控制问题存在最优解 $(f_{\varepsilon},(U_{\varepsilon},V_{\varepsilon}))$, 其中控制函数 $f_{\varepsilon}$ 有表达式

$f_{\varepsilon}=1_{\omega}\eta_{\varepsilon}{\rm e}^{\kappa\alpha},$

函数对 $(U_{\varepsilon},V_{\varepsilon})$ 为方程组 (3.1) 相应于 $f=f_{\varepsilon}$ 的解, $(\eta_{\varepsilon},\zeta_{\varepsilon})$ 为共轭方程组 (3.3)相应于

$(\eta_{\varepsilon},\zeta_{\varepsilon})(\cdot,T)=-\frac{1}{\varepsilon}(U_{\varepsilon},V_{\varepsilon})(\cdot,T)$

的解. 根据线性方程组 (3.1) 和共轭方程组 (3.3) 之间的对偶关系以及引理 3.3 可知

$\Vert1_{\omega}f{\rm e}^{-\frac{1}{2}\kappa\alpha}\Vert_{2}^{2}+\frac{1}{\varepsilon}\left\Vert (U_{\varepsilon},V_{\varepsilon})(\cdot,T)\right\Vert_{L^{2}(\Omega)^{2}}^{2}\leq {\rm e}^{C(\varpi T+\frac{1}{T}+1)}\left\Vert(U_{0},V_{0})\right\Vert _{L^{2}(\Omega)^{2}}^{2}. $

注意控制函数 $f_{\varepsilon}$ 的表达式 (3.18), 则可知 $f_{\varepsilon}$ 在 $L^{2}(Q_{T})$ 中有估计

$\Vert f_{\varepsilon}\Vert_{2}\leq {\rm e}^{C(\varpi T+\frac{1}{T}+1)}\left\Vert (U_{0},V_{0})\right\Vert _{L^{2}(\Omega)^{2}}.$

下证控制函数 $f_{\varepsilon}$ 在 $L^{\infty}(Q_{T})$ 中的估计. 令 $1<\mu<\kappa$, 并取单调递增的有限序列 $\{\mu_{j}\}_{0}^{\ell+1}$ 使得 $0<\mu_{j}<\mu_{\ell+1}=\mu.$ 令 $\alpha_{0}:=\alpha_{0}(t)=\min\{\alpha(x,t);x\in\overline{\Omega}\}$, 则

$\alpha_{0}\leq\alpha\leq\alpha_{0}[1+\gamma(\lambda)^{-\frac{1}{2}}]^{-1}<0,$

其中 $\gamma(\lambda)$ 和 $\alpha=\alpha(x,t)$ 由 (3.5) 式给出. 对每个 $i$ $(i=0,1,\cdots,\ell,\ell+1)$, 定义

$\widehat{\eta}_{i}(x,t)={\rm e}^{\mu_{i}\alpha_{0}}\eta_{\varepsilon}(x,T-t),\ \widehat{\zeta}_{i}(x,t)={\rm e}^{\mu_{i}\alpha_{0}}\zeta_{\varepsilon}(x,T-t),$

则函数 $ (\widehat{\eta}_{i},\ \widehat{\zeta}_{i})$ 为方程组

$\begin{matrix}\partial_{t}\widehat{\eta}_{i}-\triangle\widehat{\eta}_{i}=\widehat{B}_{2}\cdot\nabla\widehat{\eta}_{i}+\widehat{a}_{1}\widehat{\eta}_{i}-\widehat{c}_{1}\widehat{\zeta}_{i}+\mu_{i}\alpha_{0}^{\prime}\widehat{\eta}_{i} & \mbox{在 $Q_{T}$ 中},\\ \partial_{t}\widehat{\zeta}_{i}-\triangle\widehat{\zeta}_{i}=-\nabla \cdot(\widehat{B}_{1}\nabla\widehat{\eta}_{i})+\widehat{a}_{2}\widehat{\eta}_{i}-\widehat{c}_{2}\widehat{\zeta}_{i}+\mu_{i}\alpha_{0}^{\prime}\widehat{\zeta}_{i} \mbox{在 $Q_{T}$ 中},\\ \partial_{\nu}\widehat{\eta}_{i}=\partial_{\nu}\widehat{\zeta}_{i}=0 \mbox{在 $\Sigma_{T}$ 上},\\ \widehat{\eta}_{i}(x,0)=0,\ \ \widehat{\zeta}_{i}(x,0)=0 & x\in\Omega. \end{matrix}$

的解, 其中系数 $\widehat{a}_{1}$ 定义为 $\widehat{a}_{1}(x,t)=a_{1}(x,T-t),$ 其余系数 $\widehat{a}_{2},\widehat{c}_{1}, \widehat{c}_{2},\widehat{B}_{1}$ 和 $\widehat{B}_{2}$ 用类似方式给出. 利用半群 $\{{\rm e}^{-tA}\}_{t\geq0}$ 将方程组 (3.22) 的解表示为

$\widehat{\eta}_{i}(\cdot,t)=\int_{0}^{t}{\rm e}^{-(t-s)A}[\widehat{B}_{2}\cdot \nabla\widehat{\eta}_{i}+\widehat{a}_{1}\widehat{\eta}_{i}-\widehat{c}_{1}\widehat{\zeta}_{i}+\mu_{i}\alpha_{0}^{\prime}\widehat{\eta}_{i}](\cdot,s){\rm d}s, $
$\widehat{\zeta}_{i}(\cdot,t)=\int_{0}^{t}{\rm e}^{-(t-s)A}[-\nabla\cdot (\widehat{B}_{1}\nabla\widehat{\eta}_{i})+\widehat{a}_{2}\widehat{\eta}_{i}-\widehat{c}_{2}\widehat{\zeta}_{i}+\mu_{i}\alpha_{0}^{\prime}\widehat{\zeta}_{i}](\cdot,s){\rm d}s.$

取单调递增的有限序列 $\{p_{i}\}_{0}^{\ell}$ 使得 $p_{0}=2,$ $p_{\ell}>5/2,$ 且 $0<1/p_{i-1}-1/p_{i}<1/10.$ 应用半群的 $L^{p_{i-1}}$-$L^{p_{i}}$ 估计[15]有

$\int_{0}^{t}\left\vert {\rm e}^{-(t-s)A}(\widehat{a}_{1}\widehat{\eta}_{i})(\cdot,s)\right\vert _{p_{i}}{\rm d}s\leq\left\Vert a_{1}\right\Vert _{\infty}\int_{0}^{t}m(t-s)^{-\frac{3}{2}(\frac{1}{p_{i-1}}-\frac{1}{p_{i}})}|\widehat{\eta}_{i}(\cdot,s)|_{p_{i-1}}{\rm d}s,$

其中 $m(t)=\min\{1,t\}.$ 对上式应用 Young 卷积不等式可得

$\left\Vert \int_{0}^{t}{\rm e}^{-(t-s)A}\widehat{a}_{1}\widehat{\eta}_{i}(\cdot,s){\rm d}s\right\Vert _{p_{i}}\leq\left\Vert a_{1}\right\Vert _{\infty}{\rm e}^{CT}\left\Vert \widehat{\eta}_{i}\right\Vert _{p_{i-1}}.$

用类似的方式估计 (3.23) 和 (3.24) 式中其余各项得到

$\Vert(\widehat{\eta}_{i},\widehat{\zeta}_{i})\Vert_{p_{i}}\leq {\rm e}^{C\varpi T}[\Vert(\widehat{\eta}_{i},\widehat{\zeta}_{i})\Vert_{p_{i-1}}+\Vert \nabla\widehat{\eta}_{i}\Vert_{p_{i-1}}+\Vert\mu_{i}\alpha_{0}^{\prime}(\widehat{\eta}_{i},\widehat{\zeta}_{i})\Vert_{p_{i-1}}].$

一方面, 根据命题 2.1

$\Vert(\widehat{\eta}_{i},\widehat{\zeta}_{i})\Vert_{p_{i-1}}+\Vert \nabla\widehat{\eta}_{i}\Vert_{p_{i-1}}\leq {\rm e}^{C\varpi T}\Vert\mu_{i} \alpha_{0}^{\prime}(\widehat{\eta}_{i},\widehat{\zeta}_{i})\Vert_{p_{i-1}},$

另一方面, 根据 $\widehat{\eta}_{i}$ 和 $\widehat{\zeta}_{i}$ 的表达式有估计

$\Vert\mu_{i}\alpha_{0}^{\prime}(\widehat{\eta}_{i},\widehat{\zeta}_{i})\Vert_{p_{i-1}}\leq CT\Vert(\widehat{\eta}_{i-1},\widehat{\zeta}_{i-1})\Vert_{p_{i-1}},$

由此根据 (3.25)-(3.27) 式得到迭代不等式

$\Vert(\widehat{\eta}_{i},\widehat{\zeta}_{i})\Vert_{p_{i}}\leq {\rm e}^{C\varpi T}\Vert(\widehat{\eta}_{i-1},\widehat{\zeta}_{i-1})\Vert_{p_{i-1}},$
$\Vert(\widehat{\eta}_{\ell},\widehat{\zeta}_{\ell})\Vert_{p_{\ell}}\leq {\rm e}^{C\varpi T}\Vert(\widehat{\eta}_{0},\widehat{\zeta}_{0})\Vert_{2}\leq {\rm e}^{C(\varpi T+\frac{1}{T}+1)}\left\Vert (U_{0},V_{0})\right\Vert _{L^{2}(\Omega)^{2}}$

此外, 对 (3.22) 式的第一个方程应用命题 2.1, 可得

$\Vert\widehat{\eta}_{\ell+1}\Vert_{W_{p_{\ell}}^{2,1}(Q_{T})}\leq {\rm e}^{C\varpi T}\Vert(\widehat{\eta}_{\ell},\widehat{\zeta}_{\ell})\Vert_{p_{\ell}}.$

因为 $p_{\ell}>5/2$ 有嵌入 $W_{p_{\ell}}^{2,1}(Q_{T})\subset C(\overline{Q}_{T}),$ 所以

$\Vert\widehat{\eta}_{\ell+1}\Vert_{\infty}\leq {\rm e}^{C(1+T+\frac{1}{T})}\Vert\widehat{\eta}_{\ell+1}\Vert_{W_{p_{\ell}}^{2,1}(Q_{T})}\leq {\rm e}^{C(\varpi T+\frac{1}{T}+1)}\left\Vert (U_{0},V_{0})\right\Vert _{L^{2}(\Omega)^{2}}.$

注意 $\widehat{\eta}_{\ell+1}$ 和 $f_{\varepsilon}$ 的定义以及不等式 (3.21)

$\left\Vert f_{\varepsilon}\right\Vert _{\infty}=\left\Vert {\rm e}^{\kappa\alpha -\mu\alpha_{0}}\widehat{\eta}_{\ell+1}\right\Vert _{\infty}\leq\Vert {\rm e}^{[\kappa-\mu-\mu\gamma(\lambda)^{-1/2}]\alpha}\widehat{\eta}_{\ell+1}\Vert_{\infty}.$

故可在引理 3.3 中取参数 $\lambda$ 足够大使得 $\kappa-\mu-\mu\gamma(\lambda)^{-1/2}>0$, 从而根据 (3.29) 和 (3.30) 式可得 $f_{\varepsilon}$的$L^{\infty}$ 估计式

$\left\Vert f_{\varepsilon}\right\Vert _{\infty}\leq {\rm e}^{C(\varpi T+\frac{1} {T}+1)}\left\Vert (U_{0},V_{0})\right\Vert _{L^{2}(\Omega)^{2}}.$

最后, 由于序列 $\{f_{\varepsilon}\}$ 在 $L^{\infty}(Q_{T})$ 中有界, 故存在子列 $\{f_{\varepsilon_{k}}\}_{k=1}^{\infty}$, 当 $\varepsilon_{k}\rightarrow0$ 时, $f_{\varepsilon_{k}}$ 在 $L^{2}(Q_{T})$ 中弱收敛于 $f\in L^{\infty}(Q_{T}).$ 设 $(U_{\varepsilon},V_{\varepsilon})$ 是线性方程组 (3.1) 相对于 $f_{\varepsilon}$ 的解. 根据命题 2.1, 序列 $\{U_{\varepsilon}\}$ 和 $\{V_{\varepsilon}\}$ 在 $L^{2}(0,T;H^{1}(\Omega))\cap W^{1,2}([T];H^{1}(\Omega)^{\ast})$ 中有界, 从而也存在子列 $\{U_{\varepsilon_{k}}\}_{k=1}^{\infty}$ 和 $\{V_{\varepsilon_{k}}\}_{k=1}^{\infty}$ 使得 $U_{\varepsilon_{k}}$, $V_{\varepsilon_{k}}$ 分别在 $L^{2}(0,T;H^{1}(\Omega))$ 和 $W^{1,2}([T];H^{1}(\Omega)^{\ast})$ 中弱收敛于 $U$, $V$, 且 $U$ 和 $V$ 属于 $L^{2}([0,T);H^{1}(\Omega))\cap C([T],L^{2}(\Omega))$. 取 $\varepsilon_{k}\rightarrow0$ 时的极限可知 $(U,V)$ 是方程组 (3.1) 相对于 $f\in L^{\infty}(Q_{T})$ 的解. 在式 (3.19) 中取 $\varepsilon_{k}\rightarrow0$ 的极限可知 $(U,V)(x,T)=(0,0)$, a.e. $x\in\Omega;$ 即证方程组 (3.1) 的零能控性. 定理 3.1 得证.

4 非线性控制系统的能控性及其应用

本节将证明非线性控制系统(1.1)的局部精确能控性, 并探讨其在时间最优控制问题中的应用.

定理 4.1 设函数 $h(u,v)$ 满足

$h_{u}^{\prime}\geq\sigma_{0}>0. $

令 $p>5$, $(\overline{u},\overline{v})$ 为方程组 (1.2) 相对于 $(\overline{u}_{0},\overline{v}_{0})$ 的轨迹, 则存在不依赖于 $T$ 的正常数 $\tau_{0}$, 使对任意的 $(\overline{u}_{0},\overline{v}_{0})$, 当

$|u_{0}-\overline{u}_{0}|_{\infty}+\Vert v_{0}-\overline{v}_{0}\Vert _{W^{2(1-\frac{1}{p}),p}(\Omega)}\leq {\rm e}^{-\tau_{0}(1+T+\frac{1}{T})}$

时, 存在 $f\in L^{\infty}(Q_{T})$, 使得 (1.1) 的解 $(u,v)$ 满足 $(u,v)(x,T)=(\overline{u},\overline{v})(x,T)$, a.e. $x\in\Omega$.

证 令 $u=\overline{u}+U,$ $v=\overline{v}+V,\ u_{0}=\overline{u}_{0}+U_{0},\ v_{0}=\overline{v}_{0}+V_{0}$, 则方程组 (1.1) 可写为

$\begin{cases} U_{t}=\triangle U-k_{0}\nabla\cdot(\overline{u}\nabla V)-k_{0}\nabla \cdot(U\nabla(\overline{v}+V) & \quad\\ +g(\overline{u}+U,\overline{v}+V)-g(\overline{u},\overline{v})+1_{w}f & \mbox{在 $Q_{T}$ 中},\\ V_{t}=\triangle V-h(\overline{u}+U,\overline{v}+V)+h(\overline{u},\overline {v}) \mbox{在 $Q_{T}$ 中},\\ \partial_{\nu}U=\partial_{\nu}V=0 & \mbox{在 $\Sigma_{T}$ 上},\\ U(x,0)=U_{0}(x)\ \ V(x,0)=V_{0}(x) & x\in\Omega, \end{cases}$

因此方程组 (1.1) 关于轨迹 $(\overline{u},\overline{v})$ 的精确能控性等价于方程组 (4.3) 的零能控性. 令 $X=L^{2}(Q_{T})\times L^{2}(0,T;H^{1}(\Omega)),$ 取其子集

$K=\{(\alpha,\beta)\in X\ ;\Vert\alpha\Vert_{\infty}+\Vert\beta\Vert _{L^{\infty}(0,T;W^{1,\infty}(\Omega))}\leq1\}.$

对任意的 $(\alpha,\beta)\in K$, 有如下的线性化方程组

$\begin{cases} U_{t}=\triangle U-\nabla(B_{1}\nabla V)-\nabla\cdot(B_{2}^{\beta} U)+a_{1}^{\alpha\beta}U+a_{2}^{\alpha\beta}V+1_{w}f & \mbox{在 $Q_{T}$ 中},\\ V_{t}=\triangle V-c_{1}^{\alpha\beta}U-c_{2}^{\alpha\beta}V \mbox{在 $Q_{T}$ 中},\\ \partial_{\nu}U=\partial_{\nu}V=0 & \mbox{在 $\Sigma_{T} $ 上},\\ U(x,0)=U_{0}(x),\ \ V(x,0)=V_{0}(x) & x\in\Omega,\end{cases}$

其中 $B_{1}=k_{0}\overline{u}\in L^{\infty}(Q_{T}),\ B_{2}:=B_{2}^{\beta}=k_{0}\nabla(\overline{v}+\beta)\in L^{\infty}(Q_{T})^{3}, $ 在 $\Sigma_{T}$ 上 $B_{2}\cdot\nu=0,$

$a_{1} :=a_{1}^{\alpha\beta}={\normalsize \int_{0}^{1}g_{u}^{\prime }(\overline{u}+\lambda\alpha,\overline{v}+\lambda\beta){\rm d}\lambda\in L^{\infty}(Q_{T}),}\\ a_{2} :=a_{2}^{\alpha\beta}=\int_{0}^{1}g_{v}^{\prime} ({\normalsize \overline{u}+\lambda\alpha,\overline{v}+\lambda\beta){\rm d}\lambda\in L^{\infty}(Q_{T}),} $

系数 $c_{1}:=c_{1}^{\alpha\beta},c_{2}:=c_{2}^{\alpha\beta}$ 类似 $a_{1}^{\alpha\beta},a_{2}^{\alpha\beta}$ 予以定义. 根据假设条件 (4.1) 可知 $c_{1}\geq\sigma_{0},$ 从而满足定理 3.1 的条件, 因此线性方程组 (4.4) 具有零能控性, 亦即对任意的 $(U_{0},V_{0})\in L^{2}(\Omega)^{2}$ 有 $f\in L^{\infty}(Q_{T})$ 使方程组 (4.4) 的解 $(U,V)\in L^{2}(0,T;H^{1}(\Omega))\cap C([T];L^{2}(\Omega))$ 满足 $(U,V)(x,T)=(0,0)$, a.e. $x\in\Omega,$ 且 $f$ 满足

$\Vert f\Vert_{\infty}\leq {\rm e}^{c_{0}(1+T+\frac{1}{T})}\left\Vert (U_{0},V_{0})\right\Vert _{L^{2}(\Omega)^{2}},$

这里以及下述中 $C$ 表示不依赖于 $T$ 的正常数. 定义多值映射 $\Lambda:K\rightarrow2^{L^{2}(Q_{T})},$

$\Lambda(\alpha,\beta)=\left\{\begin{array} [c]{c} (U,V)\in(L^{2}(0,T;H^{1}(\Omega))\cap C([T];L^{2}(\Omega))^{2};\ \exists\ f\ \text{满足 (4.5) 式,}\\ \text{使}\ (U,V)\ \text{为方程组}\ (4.4)\ \text{的解, 且}\ (U,V)(x,T)=(0,0),\ \text{a.e}\ x\in\Omega.\end{array} \right\}$

下证 $\Lambda$ 满足 Kakutani 不动点定理[17]. 首先, 根据定理 3.1 可知对每个 $(\alpha,\beta)\in K,$ $\Lambda(\alpha,\beta)$ 非空; 此外, 线性方程组的线性结构保证 $\Lambda(\alpha,\beta)$ 为凸集. 应用命题 2.1 可知 $(U,V)$ 满足

$ \Vert(U,V)\Vert_{(L^{2}([0,T);H^{1}(\Omega))\cap W^{1,2}([T];H^{1} (\Omega)^{\ast}))^{2}}+\left\Vert V\right\Vert _{L^{2}(0,T;H^{2}(\Omega))} \leq {\rm e}^{C(1+T)}[|U_{0}|_{2}+\left\Vert V_{0}\right\Vert _{H^{1}(\Omega)}].$

由此可知 $\Lambda(\alpha,\beta)$ 为 $[L^{2}(0,T;H^{1}(\Omega))\times L^{2}(0,T;H^{2}(\Omega))]\cap W^{1,2}([T];H^{1}(\Omega)^{\ast})^{2}$ 中的有界集, 于是 Aubin-Lions 引理[17] 保证 $\Lambda(\alpha,\beta)$ 是 $X$ 中的紧集.

其次, 证明 $\Lambda$ 上半连续. 为此, 取 $K$ 的序列 $\{(\alpha_{n},\beta_{n})\}_{n=1}^{\infty}$ 使其当 $n\rightarrow\infty$ 时在 $X$ 中强收敛于 $(\alpha,\beta)$. 对每个 $n,$ 取 $(U_{n},V_{n})\in\Lambda(\alpha_{n},\beta_{n})$, 则由 $\Lambda(\alpha_{n},\beta_{n})$ 的定义, 存在 $f_{n}$ 使 $(U_{n},V_{n})$ 为以下方程组的解

$\begin{cases} \partial_{t}U_{n}=\triangle U_{n}-\nabla\cdot(B_{1}\nabla V_{n})-\nabla \cdot(B_{2n}U_{n}) +a_{1n}U_{n}+a_{2n}V_{n}+1_{w}f_{n} & \mbox{在 $Q_{T}$ 中},\\ \partial_{t}V_{n}=\triangle V_{n}-c_{1n}U_{n}-c_{2n}V_{n} \mbox{在 $Q_{T}$ 中},\\ \partial_{\nu}U_{n}=\partial_{\nu}V_{n}=0 & \mbox{在 $\Sigma_{T}$ 上},\\ U_{n}(x,0)=U_{0}(x),\ \ V_{n}(x,0)=V_{0}(x) x\in\Omega,\end{cases}$

且 $(U_{n},V_{n})(x,T)=(0,0)$, a.e. $x\in\Omega$. 这里 $f_{n}$ 和 $(U_{n},V_{n})$ 分别满足不等式 (4.5) 和 (4.6), 此即表明存在子序列 $\{f_{k}\}_{k=1}^{\infty}$ 和 $\{(U_{k},V_{k})\}_{k=1}^{\infty}$ 使当 $k\rightarrow\infty$ 时, 函数 $f_{k}$ 在 $L^{\infty}(Q_{T})$ 中弱*收敛而在 $L^{2}(Q_{T})$ 中弱收敛于 $f;$ $(U_{k},V_{k})$ 在 $L^{2}(0,T;H^{1}(\Omega))^2$ 和 $W^{1,2}([T];H^{1}(\Omega)^{\ast})^2$ 弱收敛于 $(U,V).$ 于是根据 Aubin-Lions 引理, $(U_{k},V_{k})$ 在 $L^{2}(Q_{T})^{2}$ 中强收敛于 $(U,V).$ 令 $k\rightarrow\infty,$ 则 $(U,V)$ 为方程组 (4.4) 相应于 $f$ 的弱解. 下证 $(U,V)\in\Lambda(\alpha,\beta)$. 令 $Y_{k}=U_{k}-U$, $Z_{k}=V_{k}-V$. 设 $(Y_{k},Z_{k})$ 为以下方程的解

$\begin{cases} \partial_{t}Y_{k}=\triangle Y_{k}-\nabla\cdot(B_{1}\nabla Z_{k})-\nabla \cdot(B_{2k}Y_{k})+a_{1k}Y_{k}+a_{2k}Z_{k}+G_{k} \mbox{在 $Q_{T}$ 中},\\ \partial_{t}Z_{k}=\triangle Z_{k}-c_{1n}Y_{k}-c_{2k}Z_{k}+H_{k} \mbox{在 $Q_{T}$ 中},\\ \partial_{t}Y_{k}=\partial_{\nu}Z_{k}=0 & \mbox{在 $\Sigma_{T}$ 上},\\ Y_{k}(x,0)=0\ \ Z_{k}(x,0)=0 & x\in\Omega, \end{cases}$

其中 $G_{k} =-\nabla\cdot\lbrack(B_{2k}-B_{2})U]+(a_{2k}-a_{2})V+(a_{1k}-a_{1})U+1_{\omega}(f_{k}-f),$ $H_{k} =-(c_{1k}-c_{1})U-(c_{2k}-c_{2})V.$ 根据能量估计可得

$|Y_{k}(\cdot,t)|_{2}^{2}+|Z_{k}(\cdot,t)|_{2}^{2}\\ \leq\ {\rm e}^{C\varpi T}[\left\Vert U\right\Vert _{\infty}^{2}+\left\Vert V\right\Vert _{\infty}^{2}]\iint_{Q_{T}}\bigg[ |B_{2k}-B_{2}|^{2}+\sum_{i=1}^{2}\left( |a_{ik}-a_{i}|^{2}+|c_{ik}-c_{i}|^{2}\right) \bigg] {\rm d}x{\rm d}t\\ +\iint_{Q_{T}}1_{\omega}(f_{k}-f)Y_{k}{\rm d}x{\rm d}t.$

这里用到 $(U,V)\in L^{\infty}(Q_{T})$. 由于 $\alpha_{k}$ 和 $\beta_{k}$ 分别在 $L^{2}(Q_{T})$ 和 $L^{2}(0,T;H^{1}(\Omega))$ 中强收敛于 $\alpha$ 和 $\beta,$ $f_{k}$ 在 $L^{2}(Q_{T})$ 中弱收敛于$f,$ 从而 $B_{2k},a_{ik}$ 和 $c_{ik}$ 均在 $L^{2}(Q_{T})$ 中收敛于 $B_{2},a_{i}$ 和 $c_{i}(i=1,2),$因此上式右端当 $k\rightarrow\infty$ 时趋于$0,$ 亦即在 $L^{2}(\Omega)^{2}$ 中 $(U_{k},V_{k})(\cdot,t)\rightarrow(U,V)(\cdot,t).$ 又 $(U_{k},V_{k})(\cdot,T)=(0,0),$ 故可知 $(U,V)(x,T)=0$ a.e. $x\in\Omega,$ 即证 $(U,V)\in\Lambda(\alpha,\beta),$ 即证 $\Lambda$ 上半连续.

最后, 证明 $\Lambda(K)\subset K.$ 对任意的 $(U,V)\in\Lambda(K)$, 因为当 $p>5$ 时, 有嵌入 $W_{p}^{2,1}(Q_{T})\subset L^{\infty}(0,T;W^{2(1-\frac{1}{p}),p}(\Omega))\subset L^{\infty}(0,T;W^{1.\infty}(\Omega)),$ 所以根据命题 2.1 和控制函数 $f$ 的估计式 (4.5)

$\left\Vert U\right\Vert _{\infty}+\left\Vert V\right\Vert _{L^{\infty }(0,T;W^{1,\infty}(\Omega))} \leq\Vert U\Vert_{\infty}+\left\Vert V\right\Vert _{\infty}+C\left\Vert V\right\Vert _{W_{p}^{2,1}(Q_{T})}\\ \leq {\rm e}^{\tau_{0}(1+T+\frac{1}{T})}(|U_{0}|_{\infty}+\Vert V_{0}\Vert_{W^{2(1-\frac{1}{p}),p}(\Omega)}),$

显然, 当 $(U_{0},V_{0})=(u_{0}-\overline{u}_{0}, v_{0}-\overline{v}_{0})$ 满足不等式 (4.2) 时, $(U,V)\in K,$ 亦即 $\Lambda(K)\subset K.$ 综上所述, 多值映射 $\Lambda$ 在 $K\subset X$ 中满足 Kakutani 不动点定理, 即当 $(u_{0},v_{0})$ 满足 (4.2) 式时, 有不动点 $(U,V)\in\Lambda(U,V),$ 从而 $(U,V)$ 为非线性方程组 (4.3) 或 $(u,v)$ 为非线性方程组 (1.1) 的解, 且 $u(x,T)=\overline{u}(x,T)$, $v(x,T)=\overline {v}(x,T)$, a.e. $x\in\Omega$.

定理 4.2 设 $p,\tau_{0}$ 如定理 4.1 所述, 如果非零初值 $({u}_{0},{v}_{0})\in L^{\infty}(\Omega)\times W^{2(1-\frac{1}{p}),p}(\Omega)$ 满足

$|u_{0}|_{\infty}+\Vert v_{0}\Vert_{W^{2(1-\frac{1}{p}),p}(\Omega)}\leq {\rm e}^{-\tau_{0}(1+T+\frac{1}{T})},$

则问题 (1.3) 至少存在一个时间最优控制.

证 首先, 证明集合 $\Pi$ 非空. 当 $T=1$ 时, ${\rm e}^{-\tau_{0}(1+T+\frac{1}{T})}$ 可取得最大值 ${\rm e}^{-3\tau_{0}}$. 于是当 $(u_{0},v_{0})$ 满足不等式$|u_{0}|_{\infty}+\Vert v_{0}\Vert_{W^{2(1-\frac{1}{p}),p}(\Omega)}\leq {\rm e}^{-3\tau_{0}}$时, 根据定理 4.1, 存在 $f\in L^{\infty}(Q_{T})$, 使得方程组 1.1 的解 $(u,v)$ 在 $T=1$ 时满足 $(u,v)(x,1)=(0,0)$, a.e. $x\in\Omega,$ 从而表明 $T=1\in$ $\Pi,$ 即证集合 $\Pi$ 非空. 其次, 证明时间最优控制的存在性. 令 $T^{\ast}=\inf\Pi$, 并在 $\Pi$ 中取一列不增的数列 $\{T_{n}\}_{n=1}^{\infty}$ 使 $n\rightarrow\infty$ 时, $T_{n}\rightarrow T^{\ast}$. 对 $\forall n$, 由定理 4.1 可知存在控制函数 $f_{n}\in L^{\infty}(Q_{T_{n}})$, 使得方程 (1.1) 的解 $(u_{n},v_{n})\in(L^{2}(0,T_{n};H^{1}(\Omega))\cap C([T_{n}];L^{2}(\Omega)))^{2}$ 满足 $(u_{n},v_{n})(x,T_{n})=(0,0),$ a.e. $x\in\Omega.$ 设 $T>T_{1}$, 将 $u_{n}$, $v_{n}$ 和 $f_{n}$ 在 $[T_{n}]$ 外零延拓

$(\widetilde{u}_{n},\widetilde{v}_{n}) = \begin{cases} (u_{n},v_{n}), & \mbox{在 $\Omega\times(0,T_{n})$ 中},\\ (0,0), & \mbox{在 $\Omega\times\lbrack T_{n},T)$ 中}; \end{cases} \widetilde{f}_{n} = \begin{cases} f_{n}, &\mbox{在 $\Omega\times(0,T_{n})$ 中},\\ 0, &\mbox{在 $\Omega\times\lbrack T_{n},T)$ 中}, \end{cases} $

则 $(\widetilde{u}_{n},\widetilde{v}_{n},\widetilde{f}_{n})$ 满足

$\begin{cases} \partial_{t}\widetilde{u}_{n}=\nabla\cdot(\nabla\widetilde{u}_{n}-k_{0}\widetilde{u}_{n}\nabla\widetilde{v}_{n})+g(\widetilde{u}_{n},\widetilde{v}_{n})+1_{w}\widetilde{f}_{n} & \quad\mbox{在 $Q_{T}$ 中},\\ \partial_{t}\widetilde{v}_{n}=\triangle\widetilde{v}_{n}-h(\widetilde{u}_{n},\widetilde{v}_{n}) & \quad\mbox{在 $Q_{T}$ 中},\\ \partial_{\nu}\widetilde{u}_{n}=\partial_{\nu}\widetilde{v}_{n}=0 \quad\mbox{在 $\Sigma_{T}$ 上},\\ \widetilde{u}_{n}(x,0)=u_{0}(x)\ \ \widetilde{v}_{n}(x,0)=v_{0}(x) \quad x\in\Omega. \end{cases}$

注意定理 4.1 保证了控制函数 $f_{n}$ 满足

$\left\Vert f_{n}\right\Vert _{L^{\infty}(0,T_{n};L^{\infty}(\Omega))}\leq {\rm e}^{-\tau_{0}(1+T_{n}+\frac{1}{T_{n}})}\left\Vert (U_{0},V_{0})\right\Vert _{L^{2}(\Omega)^{2}}\leq {\rm e}^{-3\tau_{0}}\left\Vert (U_{0},V_{0})\right\Vert _{L^{2}(\Omega)^{2}}.$

由此 $\{\widetilde{f}_{n}\}_{n=1}^{\infty}$ 在 $L^{\infty}(Q_{T})$ 中有界; 且根据定理 4.1 的证明可知序列

$\left\Vert (\widetilde{u}_{n},\widetilde{v}_{n})\right\Vert _{L^{\infty}(Q_{T})}+\left\Vert \widetilde{v}_{n}\right\Vert _{L^{\infty}(0,T;W^{1,\infty }(\Omega))}\leq\rho\text{, }$

其中 $\rho$ 为有限的正实数, 于是$\left\Vert g(\widetilde{u}_{n},\widetilde{v}_{n})\right\Vert _{L^{\infty}(Q_{T})}+\left\Vert h(\widetilde{u}_{n},\widetilde{v}_{n})\right\Vert _{L^{\infty}(Q_{T})}\leq C_{\rho},$ 这里 $C_{\rho}$ 为依赖于 $\rho$ 的正实数. 根据能量估计可知

$\left\Vert (\widetilde{u}_{n},\widetilde{v}_{n})\right\Vert _{(L^{2}(0,T;H^{1}(\Omega))\cap W^{1,2}([T];H^{1}(\Omega)^{\ast}))^{2}}+\left\Vert \widetilde{v}_{n}\right\Vert _{L^{2}(0,T;H^{2}(\Omega))}\leq {\rm e}^{C_{\rho}T},$

由此根据 Aubin-Lions 引理可知 $(\widetilde{u}_{n},\widetilde{v}_{n})$ 在 $L^{2}(Q_{T})^{2}$ 中强收敛于 $(u^{\ast},v^{\ast}),$ 于是存在子列 (仍用相同记号) $\{(\widetilde{u}_{n},\widetilde{v}_{n})\}_{n=1}^{\infty}$ 使当 $n\rightarrow\infty$ 时, $(\widetilde{u}_{n},\widetilde{v}_{n})(x,t)\rightarrow(u^{\ast},v^{\ast})(x,t),$ a.e. $(x,t)\in Q_{T}.$ 因此当 $n\rightarrow\infty$ 时下列子序列均在 $L^{2}(Q_{T})$ 中弱收敛

$\widetilde{u}_{n}\nabla\widetilde{v}_{n}\rightarrow u^{\ast}\nabla v^{\ast}\text{;}g(\widetilde{u}_{n},\widetilde{v}_{n})\rightarrow g(u^{\ast},v^{\ast});h(\widetilde{u}_{n},\widetilde{v}_{n})\rightarrow h(u^{\ast},v^{\ast});f_{n}\rightarrow f^{\ast}.$

于是, 在 (4.10) 式中令 $n\rightarrow\infty$ 可知 $(u^{\ast},v^{\ast})\in(L^{2}(0,T;H^{1}(\Omega))\cap W^{1,2}([T];H^{1}(\Omega)^{\ast}))^{2}$ 是方程组 (1.1) 相应于 $f^{\ast}\in L^{\infty}(Q_{T})$ 的弱解$.$ 最后, 注意到 $\widetilde{u}_{n}(\cdot,T_{n})=0,$

$|u^{\ast}(\cdot,T^{\ast})|_{2}=|u^{\ast}(\cdot,T^{\ast})-u^{\ast}(\cdot,T_{n})|_{2}+|u^{\ast}(\cdot,T_{n})-\widetilde{u}_{n}(\cdot,T_{n})|_{2}\rightarrow 0.$

可得 $u^{\ast}(T^{\ast},\cdot)=0$; 同理, $v^{\ast}(\cdot,T^{\ast})=0,$ 由此 $T^{\ast}\in\Pi,$ 即证时间最优控制的存在性.

参考文献

Schaaf R.

Stationary solutions of chemotaxis systems

Trans Amer Math Soc, 1985, 292(2): 531-556

DOI:10.1090/tran/1985-292-02      URL     [本文引用: 1]

Biler P.

Local and global solvability of some parabolic systems modelling chemotaxis

Adv Math Sci Appl, 1998, 8: 715-743

[本文引用: 1]

Gajewski G, Zacharias K.

Global behavier of a reaction-diffusion system modelling chemotaxis

Math Nathr, 1998, 195(1): 77-114

[本文引用: 1]

Lorz A.

Coupled chemotaxis fluid model

Math Models Methods Appl Sci, 2010, 20(6): 987-1004

DOI:10.1142/S0218202510004507      URL     [本文引用: 1]

We consider a model system for the collective behavior of oxygen-driven swimming bacteria in an aquatic fluid. In certain parameter regimes, such suspensions of bacteria feature large-scale convection patterns as a result of the hydrodynamic interaction between bacteria. The presented model consist of a parabolic–parabolic chemotaxis system for the oxygen concentration and the bacteria density coupled to an incompressible Stokes equation for the fluid driven by a gravitational force of the heavier bacteria.

Winkler M.

Global weak solutions in a three-dimensional chemotaxis-Navier-Stokes system

Ann Inst H Poincaré Anal Non Linéaire, 2016, 33(5): 1329-1352

DOI:10.4171/aihpc      URL     [本文引用: 1]

Arumugam G, Tyagi J.

Keller-Segel chemotaxis models: A review

Acta Appl Math, 2021, 171: Article number 6

[本文引用: 1]

Zhong X, Jiang S.

Globally bounded in-time solutions to a parabolic-elliptic system modeling chemotaxis

Acta Mathematica Scientia, 2007, 27B(2): 421-429

[本文引用: 1]

Ryu S U, Yagi A.

Optimal control of Keller-Segel equations

J Math Anal Appl, 2001, 256: 45-66

DOI:10.1006/jmaa.2000.7254      URL     [本文引用: 1]

Guo B Z, Zhang L.

Local null controllability for a chemotaxis system of parabolic-elliptic type

Systems Control Lett, 2014, 65: 106-111

DOI:10.1016/j.sysconle.2013.10.010      URL     [本文引用: 3]

Guo B Z, Zhang L.

Local exact controllability to positive trajectory for parabolic system of chemotaxis

Math Control Relat Fields, 2016, 6(1): 143-165

DOI:10.3934/mcrf      URL     [本文引用: 4]

Ammar-Khodja F, Benabdallah A, Gonzalez-Burgos M, et al.

Recent results on the controllability of linear coupled parabolic problems: A survey

Math Control Relat Fields, 2011, 1: 267-306

DOI:10.3934/mcrf.2011.1.267      URL     [本文引用: 1]

Fursikov A, Imanuvilov O Yu. Controllability of Evolution Equations. Seoul: Seoul National University, 1996

[本文引用: 2]

Barbu V. Controllability and Stabilization of Parabolic Equations. Berlin: Springer, 2018

[本文引用: 1]

Evans L C. Partial Differential Equations. Providence RI: American Mathematical Society, 2010

[本文引用: 1]

Arendt W. Semigroups and evolution equations:Functional calculus, regularity and kernel estimates. Handbook of Differential Equations: Evolutionary Equations, 2002, 1: 1-85

[本文引用: 2]

Henry D. Geometric Theory of Semilinear Parabolic Equations. Berlin: Springer, 1981

[本文引用: 1]

Barbu V. Analysis and Control of Nonlinear Infinite Dimensional Systems. Bosten: Academic Press, 1993

[本文引用: 2]

/

〈 〉

版权所有 © 《数学物理学报》杂志社
地址： 武汉市71010号信箱(430071)　电话: 027-87199206(中文版)　027-87199087(英文版)
E-mail: actams@apm.ac.cn
本系统由北京玛格泰克科技发展有限公司设计开发