最近, Alias 与 Buenzli^[1-2] 提出了一种基于细胞和组织生长的数学模型, 组织合成细胞表面演化的法向加速度线性地依赖于平均曲率. 假设二维空间中的组织界面$F(t)$ 用显式的前沿跟踪参数$\gamma(s,t)$来描述, 其中$s$是一个任意的一维参数, $t$是时间. 如果参数化$\gamma(s,t)$正交, 即时间线$t\mapsto\gamma(s,t)$始终垂直于界面, 控制组织边界和细胞密度的演化方程为

其中下标表示偏导数, $v$对应于界面的法向速度, $k$是细胞分泌率(每单位时间每个细胞分泌的新组织的体积), $\rho$是组织合成细胞表面密度, $\ell$为弧长, $d\ell=gds$, $g=|\gamma_{s}|$, ${\vec{t}}=\frac{\gamma_{s}}{|\gamma_{s}|}$ 为 $F(t)$的单位切向量, ${\vec{n}}$ 是垂直于组织基质的单位外法向量, 而$\kappa={\vec{t}}\cdot{\vec{n}}_{l}$是有符号曲率, 当$\kappa>0$时, 组织基质为凸面; 当$k<0$时, 组织基质为凹面. 方程(1.2)中$D\rho_{\ell\ell}$ 表示平行于具有扩散率D的界面的单元扩散, 而$-A\rho$表示以速率A耗散减少的活性细胞数量. 如果分泌速率$k$为常数, 方程(1.2)等价于

其中, $v_{t}$对应于界面的法向加速度$\gamma_{tt}\cdot{\vec{n}}$. 因此, 在右侧的第一项使此流成为双曲曲率流. Alias 与 Buenzli^[3]利用水平集的方法研究了曲率控制增长时细胞和组织的演化, 进一步采用三种不同的数值方法进行数值模拟.

双曲曲率流是曲线或曲面运动的热点研究内容之一. 为了描述氦晶体融化时表现出来的固态原子仅在他们融化并进入液态时才运动的现象,Gurtin与 Podio-Guidugli^[4] 研究了一类满足双曲型偏微分方程的平面曲线演化问题.Rostein等^[5] 进一步分析了在双曲晶体曲率作用下的闭多项式曲线运动,并给出了数值分析.丘^[6]提出振动膜(或曲面运动)以平均曲率向量场为加速度场的双曲几何流方程. Lefloch 与Smoczyk^[7] 提出了一类法向双曲平均曲率流, 证明了曲线在有限时间内产生奇点;同时还证明了在一维图情形下演化方程弱解的整体存在性. 李与王^[8] 利用拟线性双曲方程组的局部解理论研究了文献[7]中的法向双曲平均曲率流Cauchy问题经典解的生命跨度. Lefloch与闫^[27] 证明了（奇异）自相似解的存在性及其 基于Nash-Moser迭代的加权的Sobolev空间的稳定性. 何、孔与刘^[10] 研究了丘^[6]提出的双曲平均曲率流. 孔、刘与王^[11]研究了法向加速度为平均曲率的平面曲线演化问题, 首先证明了平面曲线的保凸性, 其次分析了在初始速度非负时, 平面曲线演化的渐近性态. 孔与王^[12] 根据双曲方程的特征值理论分析了平面曲线在双曲平均曲率流作用下奇性的形成问题. 何、黄与邢^[13]研究了该曲率流的自相似解. 在双曲曲线收缩流的基础上, 王^[14⇓-16] 研究了带有外力场的双曲平均曲率流的演化问题及其奇性分析. 王^[17]对Minkowski 空间中的双曲平均曲率流作用下类空曲线Cauchy问题经典解的生命跨度进行了研究. 毛^[18]研究了强迫双曲平均曲率流, 给出了解的局部存在唯一性, 推广了文献[10,12] 中相应的结果. 毛等^[19]研究了双曲逆平均曲率流, 假设初始紧光滑超曲面是平均凸与星形的, 证明了双曲逆平均曲率流的短时间存在性. 毛等^[20] 研究了一类双曲曲线流, 曲线法向加速度为平均曲率的幂函数, 证明了解的局部存在唯一性. 王^[21]利用拟线性双曲方程组的局部解理论, 研究了一类双曲逆平均曲率流初值问题经典解的生命跨度.

曹与沃^[22]研究了一类凸超曲面的保法向双曲曲率流, 进一步在文献[23]中, 沃等人提出并研究了双曲仿射曲线流, 在此基础上, 王^[24] 对带有耗散项的双曲仿射曲线流进行了分析和研究. Notz^[25]研究了Riemannian流形中闭超曲面的运动几何方程, 在球形曲面情况下该方程可以描述一个移动肥皂泡的非理想化模型. 闫^[26] 引进了一个拟线性退化的双曲方程描述在球对称位势中闭曲面的运动, 基于一个新的Nash-Morse迭代方法证明了这种Cauchy问题大时间的存在性和具有小初值的稳定性.

本文主要对$\gamma(s,t)$为严格凸闭曲线的情况进行研究, 特别地, 当$D=A=0$时, 演化方程方程(1.1)、方程(1.3)的Cauchy问题, 即

主要结果如下.

${\bf定理1.1}$ 假设$\gamma_{0}$为光滑闭曲线, 曲率$k_{0}>0$, $v_{0}(s)\in C^{1}(S^{1})$, 则初值问题(1.4)局部解的生命跨度的下界$\underline{\delta}$为

其中, ${{\varphi}}=({\varphi}_{1},\cdots,{\varphi}_{5})^{T}=\big(\frac{\tilde{v_{0}}-\tilde{v_{0}}'}{h''+h}, \frac{-\tilde{v_{0}}-\tilde{v_{0}}'}{h''+h},h, \tilde{v_{0}},h'\big)^{T}.$

论文结构如下. 在第2节, 根据严格凸曲线的支撑函数, 曲线演化方程(1.4)被约化成曲线支撑函数满足的一个双曲型Monge-Ampére方程, 进一步给出其Riemann不变量满足的拟线性双曲方程组, 得到解的局部存在唯一性定理.在第3节, 利用双曲拟线性双曲方程组的局部解理论, 给出定理1.1的证明.

假设严格凸闭曲线$\gamma(s,t)$满足(1.4), $\theta$是曲线$\gamma$的单位外法向角, 则单位外法向量${\vec{n}}$和单位切向量${\vec{t}}$分别为

由Frenet公式可得$\frac{\partial\theta}{\partial \ell}=\kappa$. 令

其中$t(\theta,\tau)=\tau$.则曲线$\gamma$ 的支撑函数为

求$S$关于$\theta$的一阶偏导数

则曲线$\tilde{\gamma}$可表示为

易得曲线$\tilde{\gamma}$的平均曲率为

由

得支撑函数$S$关于$\tau$的一阶偏导数满足

及

由$S_{\theta}=(\widetilde{\gamma},{\vec{t}})$ 知

又由

可得

从而曲线的支撑函数$S$满足方程

即

因此方程(1.4)可简化为

其中$h$是$\gamma_{0}$的支撑函数, 且初始速度$\widetilde{v_{0}}>0$. 易知$S\equiv h(\theta)$是(2.3)的静态解, 此时, $\widetilde{v_{0}}(\theta)=0$. 本文我们只考虑$S_{\tau}>0$的情况, 该方程为双曲型Monge-Ampére 方程. 事实上, 由

可知

则

下文中将$\tau$记为$t$, 设 $S(\theta,t)\in C^{3}\tilde{D}$是方程(2.3)的解. 则易知函数

是方程(2.3)的Rimeann不变量. 设 $p=S_{t}$, $q=S_{\theta}$, 则方程(2.3) 转化成Riemann不变量形式下的拟线性双曲方程组的Cauchy问题

对应初值为

若令${u}=(r,s,S,p,q)^{T}$, 其中$u_{1}=r,u_{2}=s,u_{3}=S,u_{4}=p,u_{5}=q$.

则方程(2.5)及初值(2.6)可改写为

首先假设对任意向量${u}=(u_{1},\cdots,u_{5})^{T}\in {\Bbb R} ^{5}$, 定义

对于向量函数${u}={u}(t,\theta)=(u_{1}(t,\theta),\cdots,u_{5}(t,\theta))^{T}$ 的模$\|{u}\|$定义为

而如果向量函数${u}(t,\theta)$是连续可微的, 还可定义$\|{u}\|_{1}$为

类似地, 可定义任一有限的函数集$G$的模 $\|G\|, \|G\|_{1}$.

由(2.6)式可知, ${\varphi}=(\varphi_{1},\cdots,\varphi_{5})^{T}$ 是$D$ 上有界的,且${\varphi}\in C^{1}(D)$, 则存在正常数$M$, 使得$\|{\varphi}\|<M$. 记$R(\delta_{0})=\{(t,\theta)|(t,\theta)\in [\delta_{0}]\times [2\pi]\}$, 并将其扩展区域定义为$E(\delta_{0})=\{(t,\theta,{u})\}|(t,\theta)\in R(\delta_{0}),|{u}|\leq M\}$. 进一步假设$\lambda_{i},\frac{\partial\lambda_{i}}{\partial \theta},\frac{\partial\lambda_{i}}{\partial u_{j}},\mu_{i}, \frac{\partial\mu_{i}}{\partial \theta},\frac{\partial\mu_{i}}{\partial u_{j}}$是$E(\delta_{0})$上的连续函数, 那么由标准的双曲型偏微分方程理论, 存在常数$0<\delta_{*}\leq\delta_{0}$, 使得方程组(2.7)的Cauchy问题在$R(\delta_{*})$上存在唯一的连续可微解$u(t,\theta)$, 即定理2.1.

${\bf定理2.1}$ (局部存在唯一性) 假设$\gamma_{0}$是严格凸闭曲线, $v_{0}(s)\in C^{1}(S^{1})$, 则存在$T>0$ 与一簇严格闭凸曲线$\gamma(\cdot,t)$, 使得当$t\in [0,T)$时, $\gamma(\cdot,t)$满足(1.4)式.

在本节中, 根据拟线性双曲方程组的Cauchy问题的局部解理论研究局部解的生命跨度.

第一步, 线性化拟线性双曲方程组(2.7). 为了进行解的第二先验估计, 选取适当的正常数$M_{1}$, 使得

对于任意${\nu}\in\sum(\delta|M,M_{1})$, 引入下列函数集合

根据一阶线性偏微分方程组

定义函数集$G^{*}=\left\{\lambda_{l},\mu_{l}, \frac{\partial\lambda_{l}}{\partial \theta},\frac{\partial\lambda_{l}} {\partial u_{j}},\frac{\partial\mu_{l}}{\partial \theta},\frac{\partial\mu_{l}}{\partial u_{j}},l,j=1,\cdots,5\right\}$, 其中$\frac{\partial\lambda_{l}}{\partial \theta}=\frac{\partial\mu_{l}}{\partial \theta}=0$$(l=1,\cdots,5,)$,有

则由

可知, $\|G^*\|$为

令函数集合$G_{2}[{\nu}]=\{\tilde{\lambda}_{l}(t,\theta),\frac{\partial \tilde{\lambda}_{l}}{\partial \theta}(t,\theta,),\tilde{\mu}_{l}(t,\theta),\frac{\partial \tilde{\mu}_{l}}{\partial \theta}(t,\theta)\}$, 则由

可知

从而, $\|G_{2}[{\nu}]\|$为

首先, 证明Cauchy问题(2.7)解的第一限先验估计式. 定义

为沿第$l$族特征线的微商$\frac{d}{d_{l}\varsigma}$. 则由第一相关积分,

其中$u_{k}^{0}(t,\theta)=\varphi_{k}(x_{k}(0;t,\theta)),k=1,\cdots,5.$ 当$0\leq t\leq\delta$时

下证Cauchy问题(2.7)解的第二先验估计式, 即解的一阶导数的上界. 记

第二相关积分为

其中

则在$R(\delta)$上, 可得

令

则成立

由(3.5)式, 当$0\leq t\leq\delta$时

由Gronwall不等式

从而有

由

将$w(t)\leq\left[\|\dot{{\varphi}}\|+\left(4M+1\right)M_{1}\delta\right]e^{M_{1}t}$ 代入上式,得

从而有

于是由$\|{u}\|_{1}$定义可知

假设$\delta\leq\frac{1}{M_{1}}$, 由Taylor展开式

代入(3.8)式, 并简化$\delta^{2}$项,

下面将给出Cauchy问题(2.7)的第三估计式, 即估计解的一阶偏导数的连续性模. 设$\psi(t,\theta)$是定义在区域$R$上的函数, $\psi(t,\theta)$在$R$上的连续性模定义为非负值函数

类似地, 可分别定义向量函数${\psi}=(\psi_{1}(t,\theta),\cdots,\psi_{n}(t,\theta))^{T}$ 或一般的有限个函数的集合$\Psi=\{\psi_{i}\}$ 的连续性模为

第$i$条特征曲线$x_i=x_{i}(\varsigma;t,\theta)$满足

有

从而, 对于$i=1,\cdots,5$, 我们有

因此

由

记$\chi(\eta)=\max\{\rho(\eta|\mu_{i}),\rho(\eta|\lambda_{i})\}$, 则

综上可得

上两式相加

令

由第二相关积分可以估计$\rho(\tau,\eta|{w}), \rho(\tau,\eta|{v})$. 令

当$\bar{t}\leq\delta$时

由$\bar{\mu}_{i}=\frac{\partial\tilde{\mu}_{i}}{\partial \theta}-\frac{\partial\tilde{\lambda}_{i}}{\partial \theta}w_{i}, i=1,\cdots,5.$ 得

即

由连续模的性质

可得

因此

由Gronwall不等式知

下一步构造非线性算子的一致有界性. 由拟线性双曲方程组的局部解理论, 在区域$R(\delta)$上存在唯一的连续可微解${u}=(u_{1}(t,\theta),\cdots, u_{5}(t,\theta))^{T}$. 记

则算子$T$的不动点就是Cauchy问题(2.7)的解, 只须证明$\delta$足够小时,算子$T$有唯一的不动点即可.

首先证明存在$\delta_{1}=\delta_{1}(\|{\varphi}\|,\|{\varphi}\|_{1})>0$ 足够小, 算子$T:\sum(\delta|M,M_{1})\to \sum(\delta|M,M_{1})$, 其中$\delta\leq\delta_{1}$. 在$R(\delta)$上, 由${u}=T{v}$可知

取$M=2\|{\varphi}\|$,当

成立

由(3.12)式得

取

记

则可取

于是, 取

则算子$T$是一个从函数集合$\sum(\delta|M,M_{1})$ 到自身的变换, 其中$\delta\leq\delta_{1}$.

为了算子$T$的等度连续性, 我们考虑$\delta_{2}$的选取. 令

由连续模的定义与$M_{1}$的构造可知

由(3.16)式及

可得

令

则可取

令 $C_{0}(\|{\varphi}\|)=12(1+2\|{\varphi}\|)^{2}(8\|{\varphi}\|+1)^2$, $C_{1}(\|{\varphi}\|)=12(1+2\|{\varphi}\|)^{2}$, 取

那么当$0\leq\delta\leq\delta_{2}$ 时, 可得

与

从而取 $\delta_{3}=\min\{\delta_{1},\delta_{2}\}$, 定义

则算子$T$是一个将函数集合$\sum(\delta_{3})$ 映到自身的算子. 由于在$C^{1}$ 模下, $\sum(\delta_{3})$集合一致有界且等度连续, 所以$\sum(\delta_{3})$是准紧的. 而

则

从而$\sum(\delta_{3})$在$C^{1}[R(\delta)]$中闭, 即为紧的, 因此$\sum(\delta_{3})$是完备的. 由$\|{\nu}\|\leq\|{\nu}\|_{1}$可知$C^{1}$模下的收敛性可以推出$C^{0}$ 模下的收敛, 从而$\sum(\delta_{3})$在$C^{0}$ 模下也是完备的. 因此, 只需证明存在$\delta_{*}\leq\delta_{3}$, 使得$T$为$R(\delta_{*})$上的压缩映射, 则不动点问题有解.

令 ${\nu}^{1},{\nu}^{2}\in\sum(\delta_{3})$, 则${u}^{1}=T{\nu}^{1},{u}^{2}=T{\nu}^{2}\in\sum(\delta_{3})$. 设${\nu}^{*}={\nu}^{1}-{\nu}^{2}$, ${u}^{*}={u}^{1}-{u}^{2}$, 且

则有

其中$\tilde{h}_{i}=-(\tilde{\lambda}_{i}^{1}-\tilde{\lambda}_{i}^{2})\frac{\partial u^{2}_{i}}{\partial \theta} +(\tilde{\mu}^{1}_{i}-\tilde{\mu}^{2}_{i})$, $i=1,\cdots,5$. 在$R(\delta)$上, 计算可得

由第一先验估计, 可得

于是取

则有

那么当选取$\delta_{*}=\min\{\delta_{3},\delta_{4}\}$ 时, 算子$T$ 是$R(\delta_{*})$上的压缩映射.

综上所述,可知 $\delta_{*}=\min\{\delta_{3},\delta_{4}\}$, 又因为$\delta_{3}<\delta_{4}$. 因此,Cauchy问题的局部解的存在区间的下界$\delta_{*}$为

定理1.1得证.