本文研究如下 Cauchy 问题解的大时间行为
方程 (1.1) 是以下单个守恒律方程的一个特殊情形
我们知道,即使初值非常光滑,方程 (1.3) 的特征线也可能会在有限时间内相交而形成激波,从而使得解产生间断, 因而,在一般情况下问题 (1.3)(1.2) 整体经典解是不存在的. 因此,我们自然考虑问题 (1.3)(1.2)的弱解即在分布(或广义函数)意义下的解. 一个函数 u(x,t) 称之为问题 (1.3)(1.2)的熵解, 如果它是弱解,并且 对几乎所有的时间t>0,u(x,t) 满足下面的熵条件
对齐次问题 (1.4)(1.2),即 g(u)=0的情形,其问题解的大时间性态已经被很多人研究,参见文献[4, 10, 13]等. 对初值是周期的且 g(u)=0 的情形,Lax 在文献[12]中证明了问题 (1.4) (1.2) 的解以O(t−1) 的速率收敛到初始值 u0 在一个周期内的平均值. 类似的 结果也可以在 文献[1, 9]中找到.
类似于 齐次守恒律方程的情形,非齐次问题 (1.3)(1.2) 解的大时间行为也有一定的物理意义. 例如, 当源项g(u) 满足 ug(u)<0,其中u>0 或 u<0,则随着时间的增长整个系统中流体的质量随着时间的无限增大而趋向于零. 此时g(u) 表示耗散项,即对流体的质量消耗. 在1990 年,Lyberopoulos 在文献[14]中 研究了特殊的源项,即 g(u)=u的情形,并在初值 u0 是周期的以及
情形1 解收敛于行波并且波速为 f′(a2m+1), 其中a2m+1 是 g(u) 的一个零点;
情形2 解收敛于 g(u)的一个零点.
在本文中,我们将继续研究以下的问题: 对给定的一个初值u0,如何去判断解是收敛于一个行波或者收敛于常数. 本文中,我们假设
(H1)~ 初值 u0(x) 在区间[0,L]上全变差有界并且是周期的, 即对某些常数 L>0,有
(H2)~ g(u)∈C(R) 并且具有有限个零点: a1<a2<⋯<a2k+1,并且
g(u)的图像如图 1 所示.
在这里不等式(1.11)可以确保问题 (1.1)(1.2) 弱解的存在性,而且 在本文中 我们假设 f 是关于u 的偶函数 使得 u 保持与 初值具有同样的性质. 因而,在这里我们只考虑 f(u)=u22 的特殊情形.然而,条件(1.10) 和 (1.11) 都是很强的,我们希望在今后的工作中减弱.
本文的主要结果如下:
定理1.1 假设条件(H1) 和 (H2) 成立,并且 u(x,t) 问题 (1.1)(1.2)的熵解.则 当 t→+∞ 时,u(x,t) 几乎处处收敛于行波解 φ(x−a2i+1t) 并且函数 φ 可以被显式的表达出来.
本文的主要内容安排如下. 在第二节,我们简单回顾一下广义特征线的定 义及其相关性质. 第三节 主要分析解渐近性态并在最后我们给出 定理 1.1 的证明. 最后,在第四节里 我们举出了两个关于定理 1.1 具体例子.
首先,让我们简要的回顾一下广义特征线理论的一些基本事实, 这些结果可在文献 []找到. 这种方法可以直接应用到非齐次方程 (1.3),其证明类似于齐次情形. 因此,在这里我们省略具体的证明.
在区间[a,b]上的一条Lipschitz曲线称为方程 (1.3) 的特征线, 如果对几乎所有的 t∈[a,b],在分布意义下满足方程
根据具有非齐次间断项的常微分方程解的理论[8],对任意的 (ˉx,ˉt)∈R×(0,∞),至少存在一条定义在最大存在区间(s,ˉt],s≥0上的后向广义特征线 ξ(t;ˉx,ˉt),使得ξ(ˉt;ˉx,ˉt)=ˉx在(x,t). 平面上所有经过点(ˉx,ˉt) 的后向广义特征线都限制在最小和最大后向广义特征线所张成的 一个锥形区域之中. 我们用ξ−(t;ˉx,ˉt), ξ+(t;ˉx,ˉt) 分别表示最小和最大后向广义特征线, 如图 2所示.
在 (H1)和(H2)的假设下,问题 (1.3)(1.2)的特征线具有如下性质.
引理2.1 设 ξ(t):[a,b]⟶R 和ζ(t):[a,b]→R,0≤a<b<∞, 是 Lipschitz 连续曲线. 则,对几乎所有的 σ,τ,a≤σ<τ≤b,
由这个引理,我们可以推出
引理2.2 设 ξ(t):[a,b]⟶R 是方程 (1.3)的特征线. 则对几乎所有的 t∈[a,b],
注2.1 引理 2.2 表明,广义特征线要么沿着真正特征线速度传播,要么沿着激波速度传播.
引理2.3 设ξ(t):[a,b]⟶R 是真正的特征线,则 存在一个函数 v:[a,b]⟶R, 使得对几乎所有的 t∈(a,b),(ξ,v) 满足
引理2.4 最小和最大的特征线 ξ−(t;ˉx,ˉt),ξ+(t;ˉx,ˉt) 都是真正的特征线, 并且它们是分别以在 t=ˉt 时值
由(H2) 可以推出 g(a2i+1)=0,i=0,1⋯,k, 并且对某个固定的 i,(1.10)式成立. 不失一般性,我们只需考虑其中一个零点 a2i+1 即可, 因为其它的零点可以作类似的处理.
引理3.1 在定理1.1假设条件下,问题(1.1) (1.2)的解 u(x,t) 也是周期的,即
证 利用本文的假设和弱解的定义,我们可以直接得出结论.
设 [a,a+L] 为在 t=0 时的一个固定周期区间,且 ξ0(t)=a2i+1t+a. 记
引理3.2 Xa(t) 是非空的并且是紧的.
证 首先,我们用反证法证明 Xa(t)≠∅. 若不然,则存在 ˉt 使得 Xa(ˉt)=∅. 由 (3.3) 式,存在 ˉx∈(ξ0(ˉ(t)),ξ0(ˉ(t))+L),使得 u(ˉx+,ˉt)>a2i+1, u(ˉx−,ˉt)<a2i+1 或 u(ˉx+,ˉt)<a2i+1<u(ˉx−,ˉt). 由熵条件 (1.5) 第一种情况不成立. 若第二种情况成立, 我们推出在 (ˉx,ˉx+L) 上一定存在 g(u) 的零点. 若不然,因为 u(x,t) 是周期的,我们有 u(ˉx+,ˉt)<a2i+1,u((ˉx+L)−,ˉt)>a2i+1. 因此,一定存在一个 ˆx∈(ˉx,ˉx+L) 使得 u(ˆx+,ˉt)>a2i+1, u(ˆx−,ˉt)<a2i+1,而这与 (1.5) 式相矛盾. 紧性的证明直接通过 Xa(t) 的定义即可. 引理证毕.
引理3.3 对任意的 0<s<t,Xa(t)⊆Xa(s).
证 对任意的 x∈Xa(t),存在经过点 (x,t) 最小和最大特征线 ξ(s) 满足 ξ±(t)=x,u(x+,t)=u(x−,t)=a2i+1. 我们断言,对所有的 s∈(0,t),
令 E(s)=(v(s)−a2i+1)2,那么
记 Xa(∞)=⋂t>0Xa(t),则我们有如下引理.
引理3.4 Xa(∞) 非空并且是紧的.
证 选取正的单调递增序列 {tk}∞k=1, 则由引理 3.3 可知 Xa(tk)⊆⋯⊆Xa(t1). 因此
定义集合 Ht 如下
对任意 t>0,记
引理3.5 ℜ(u(⋅,t)) 是闭区间并且存在 T0>0,使得当 t>T0 时,ℜ(u(⋅,t))⊂(a2i,a2i+2).
引理3.6 若对任意 x∈Sk(t),ξ(s;x,t) 是经过 (x,t) 的最小或最大特征线. 那么对任意 ε>0, 存在 T(ε)>0 使得当 t>T(ε) 时, ξ(0;x,t)∈[ak,ak+ε)∪(bk−ε,bk].
证 若引理不正确,则对任意 T0>0,存在 t1>T0 和 x1∈k(t),以及定义在 (0,t1) 上的最小或最大特征线 ξ(s;x1,t1),使得 ξ(0;x1,t1)∈(ak+ε, bk−ε). 对此 t1,存在 t2>t1+1 和 x2⊂Sk(t) 使得ξ(0;x2,t2)∈(ak+ε, bk−ε). 由此继续下去,我们可以得到一个序列 {ξ(s;xn,tn)}n∈N.
记 ξn(s)=ξ(s;xn,tn),由引理 2.3 和 2.4,存在函数 vn(s) 使得
另一方面,由于 ˙u(ξ(ˉt),ˉt)=g(u(ξ(ˉt),ˉt))>0,我们得到
对任意固定的 t>0 以及 ε>0,令
引理3.7 对由引理 3.6 得出的 T(ε)>0,若 t>T(ε), 则 Sk(t)=S−k,ε(t)∪S+k,ε(t) 并且存在一条曲线 χk(t) 使得 S−k,ε(t)=[ξk(t),χk(t)], S+k,ε(t)=[χk(t),ηk(t)].
证 由引理 3.6,对 t>T(ε),显然有Sk(t)=S−k(t)⋃S+k(t).
对固定的 ˉt,我们考虑区间 [ξk(ˉt),ηk(ˉt)],其中ξk(ˉt)≡a2i+1ˉt+ak, ηk(ˉt)≡a2i+1ˉt+bk.
假设 a0k(ˉt)=ξk(ˉt),b0k(ˉt)=ηk(ˉt),把 区间 [a0k(ˉt),b0k(ˉt)] 分成两等份, 分别记为 [a0k(ˉt),c0k(ˉt)] 和 [c0k(ˉt),b0k(ˉt)].
若 ξ−(0;c0k(ˉt),ˉt)∈[ak,ak+ε] 且 ξ+(0;c0k(ˉt),ˉt)∈[bk−ε,bk],则令 χk(ˉt)=c0k(ˉt).
若 ξ±(0;c0k(ˉt),ˉt)∈[ak,ak+ε],记 a1k(ˉt)=c0k(ˉt),b1k(ˉt)=b0k(ˉt). 若 ξ±(0;c0k(ˉt),ˉt)∈[bk−ε,bk],记 a1k(ˉt)=a0k(ˉt),b1k(ˉt)=c0k(ˉt).
由归纳法. 如果我们已经得到了闭区间 [ank(ˉt),bnk(ˉt)],令 cnk(ˉt)=(bnk(ˉt)−ank(ˉt))/2. 若 ξ−(0;cnk(ˉt),ˉt)∈[ak, ak+ε] 且 ξ+(0;cnk(ˉt),ˉt)∈[bk−ε,bk],则记χk(ˉt)=cnk(ˉt). 若 ξ±(0;cnk(ˉt),ˉt)∈[ak, ak+ε],记 an+1k(ˉt)=cnk(ˉt),bn+1k(ˉt)=bnk(ˉt). 若 ξ±(0;cnk(ˉt),ˉt)∈[bk−ε,bk],记 an+1k(ˉt)=ank(ˉt),bn+1k(ˉt)=cnk(ˉt).
由此,我们得出一闭区间序列
引理3.8 由引理 3.7 定义的曲线 χk(t) 在 [T(ε),∞) 上是 Lipschitz 连续的广义特征线,它实际上是一个激波,并且满足
证 取 Δt>0 且 t,t+Δt∈[T(ε),∞),则由 引理 3.7 可推出 ξ−(0;χk(t),t), ξ−(0;χk(t+Δt),t+Δt)∈[ak,ak+ε],且 ξ+(0;χk(t),t), ξ+(0;χk(t+Δt),t+Δt)∈[bk,bk−ε]. 因此,
由于
记
引理3.9 对 引理 3.6 中 给定的 T(ε) ,如果 t>T(ε),则当 t→∞ 时,
证 任取 x∈[ξk(t),χk(t)],由 引理 3.6, 对任意固定的 t0>0,当 t→∞ 时
定理 1.1 的 证明 首先,定义函数 φ 如下
最后,我们将要证明 φ(x) 是问题 (1.1) (1.2) 的行波解. 事实上,作变换 (x,t)=(y+a2i+1t,t) 和 U(y,t)≐u(y+a2i+1t,t). 则方程 (1.1) 化为
在本节中,给出关于定理 1.1两个具体的例子.
例1 考虑以下Cauchy 问题的解 ut+(u22)x=u13, u0(x)={1,x∈[−32,0),−1, x∈[0,32], u0(x)=u0(x+3). 直接计算可知,方程的解表示如下 u(x,t)={(53(x+32)+1)35, −32≤x<25(23t+1)52−1910,(23t+1)32, 25(23t+1)52−1910≤x<0,−(23t+1)32, 0<x<1910−25(23t+1)52,(53(x−23)+1)35, 1910−25(23t+1)52≤x≤32, 且 u(x,t)=u(x+3,t). 由直接验证知,当 t→∞ 时,u(x,t) 趋向于 φ(x)={(53(x+32)+1)35, x∈[−32,0),(53(x−32)+1)35,x∈(0,32].相应的解曲线 u(x,t) 在 (x,t) 平面上如图 3所示.
例2 考虑以下Cauchy 问题的解