非线性扩散波关于可压缩微极流体方程解的渐近稳定性
Asymptotic Stability of Nonlinear Diffusion Wave Solutions for Compressible Micropolar Fluid Equations
摘要: 本文研究一维可压缩等熵微极流体方程解的大时间行为,我们证明在初始扰动和波的强度适当小的条件下,当时间时,该方程的解收敛于平面扩散波,并得到了相应的衰减速度,本文主要的研究方法为能量方法,结合了反导数法、Cauchy不等式和Young-不等式证明了非线性扩散波关于可压缩微极流体方程解的渐近稳定性。
Abstract: In this paper, we study the large time behavior of the solution of one dimensional compressible isentropic micropolar fluid equation. Under the condition of the initial disturbance and wave intensity being appropriately small, at time , we prove that the solution of the equation converges to plane diffusion wave , and the corresponding decay rate is obtained. The main research method in this paper is the energy method, which combines the inverse derivative method, Cauchy inequality and Young inequality. By using these methods, we prove the asymptotic stability of the solution of the nonlinear diffusion wave to the compressible micropolar fluid equation.
文章引用:高倩. 非线性扩散波关于可压缩微极流体方程解的渐近稳定性[J]. 应用数学进展, 2022, 11(4): 2316-2332. https://doi.org/10.12677/AAM.2022.114244

1. 引言

在本文中我们研究如下具有粘度相关系数的一维可压缩等熵微极流体方程的大时间行为

{ v t u x = 0 u t + p ( v ) x = α u x , t > 0 ω t = ( A ( v ) ω x v ) x 4 μ r ( v ) v ω , (1.1)

给定初值

( v , u , ω ) ( x , 0 ) = ( v 0 , u 0 , ω 0 ) ( x ) , (1.2)

且满足条件

lim x ± ( v 0 , u 0 , ω 0 ) ( x ) = ( v ± , u ± , ω ± ) , (1.3)

此处t和x分别代表时间和空间变量。这里v表示比容,u表示速度, ω 表示微旋转速度, p ( v ) 表示压力, α > 0 表示阻尼系数。

在这篇文章中

p ( v ) = v γ , γ 1 , α 0 ,

并且 A ( v ) , μ r ( v ) 是关于v的光滑的正函数 ( v > 0 ) ,表示微粘度系数。

如果忽略流体的微观结构,即 ω = 0 ,系统(1.1)可视为拉格朗日坐标下带阻尼的等熵欧拉方程,动量方程中带有阻尼项 α u ,可用于模拟多孔介质中的可压缩流。

此时,系统(1.1)等同系统

{ v t u x = 0 u t + p ( v ) x = α u , (1.4)

初值满足

( u ( x , 0 ) , v ( x , 0 ) ) = ( u 0 ( x ) , v 0 ( x ) ) ( u ± , v ± ) , as x ± (1.5)

v v + ,此时方程有波的现象。刘太平 [1] [2] 证明了在初始扰动和波的强度 | u + u | 较小的情况下,当 t 时,

| v + v | < δ , ( v + v ) ,

方程(1.4)的解趋向于非线性扩散波 ( v ¯ , u ¯ ) 。其中 ( v ¯ , u ¯ ) 满足方程

{ v ¯ t u ¯ x = 0 , p ( v ¯ ) x = α u ¯ , x , t > 0 (1.6)

而对于带阻尼的非线性发展方程的柯西问题,已经有许多关于基本波的渐近稳定性的研究,参考文献 [3] [4] [5] [6]。如在一维情况下,刘太平 [5] [6] 先证明了在初始扰动和波的强度满足小性条件的情况下,非线性扩散波关于带阻尼的双曲守恒律方程的渐近稳定性。对于方程(1.4),Nishihara在 [7] [8] 进一步利用能量方法得到了和刘太平相同的结果,并获得在 L 2 L 范数下的收敛速度。对于小性条件下收敛速度其他相关结论,我们参考 [9] [10],与此同时,也有很多作者针对初始扰动较大的情况进行了研究,如Zhao [11] 证明了对于给定的大初始值情况下,带摩擦阻尼的p系统存在唯一的整体光滑解,且当时间 t ± 时,该方程的解渐近趋向于对应的非线性扩散波 ( v ¯ , u ¯ ) ( x , t ) ,并且他还得到了 L p ( 2 p ) 范数下的衰减速率。对于具有非线性阻尼的p系统的柯西问题的其他结果,我们参考 [12] - [23]。

我们的方程与(1.4)相比,增加一个方程,考虑了微旋转速度 ω ,即思考了微观结构该模型的渐近行为,增加了耦合性和研究该模型的难度,这也是本文的创新之处。

本文研究了一维可压缩微极流体模型(1.1)问题解的大时间行为,该模型具有一般的微粘度相关系数。我们将证明,在粘度系数 A ( v ) , μ r ( v ) , ( v > 0 ) 是关于v的光滑的正函数,并且初始扰动和波的强度足够小的情况下,那么方程(1.1)~(1.3)的解整体存在,并且当 t 时,该解收敛到非线性扩散波 ( v ¯ , u ¯ , 0 ) ( x , t ) 。文章主要结果的精确陈述在下面的定理1.1。

定理1.1设 A ( v ) , μ r ( v ) , ( v > 0 ) 是关于v的光滑正函数,设 ( v 0 , u 0 , ω 0 ) ( x ) H 2 × H 2 × H 2 ,则存在适当小的 δ > 0 ,当 | u + u | + | v + v | + ( v v ¯ ) ( x ) 2 + ( u u ¯ ) ( x ) 2 + ω ( x ) 2 δ 时,方程(1.1)~(1.3)的解 ( v , u , ω ) 全局存在,并且满足

v v ¯ m C ( 0 , , H 2 ) , u u ¯ u ^ C ( 0 , , H 2 ) , ω C ( 0 , , H 2 ) ,

v x v ¯ x m x L 2 ( 0 , , H 1 ) , u x u ¯ x u ^ L 2 ( 0 , , H 1 ) , ω x L 2 ( 0 , , H 1 ) ,

当t趋于无穷时,进一步有以下结论:

1 ) x k ( v ( x , t ) v ¯ ( x + x 0 , t ) m ) L 2 = O ( 1 ) δ ( 1 + t ) k 2 , 2 ) x k ( u u ¯ u ^ ) ( x , t ) L 2 = O ( 1 ) δ ( 1 + t ) k 2 , 3 ) x k ω ( x , t ) L 2 = O ( 1 ) δ ( 1 + t ) k 2 , k = 1 , 2 (1.9)

4 ) ( v v ¯ m ) ( t ) L = O ( 1 ) δ ( 1 + t ) 1 4 , 5 ) ( u u ¯ u ^ ) ( t ) L = O ( 1 ) δ ( 1 + t ) 1 4 , 6 ) ω ( t ) L = O ( 1 ) δ ( 1 + t ) 1 4 . (1.10)

其中 x 0 满足

( v ( x , 0 ) v ¯ ( x + x 0 , t ) ) d x = u + u α , (1.11)

并且

u ^ u e α t + x m t ( η , t ) d η , m 0 ( x ) d x = 1 , m ( x , t ) u + u α m 0 ( x , t ) e α t . (1.12)

本文的安排如下:在第二部分,我们给出一些在证明过程中所需相关引理;在第三部分,我们对问题(1.1)~(1.3)的扰动方程进行基本能量估计和时间加权估计,从而证明定理1.1。

标记在文章中,我们用 C , O ( 1 ) 表示正常数。对于函数空间来说, L P = L P ( n ) ( 1 p ) 表示一般的Lebesgue空间,其范数为

f L p = ( | f ( x ) | P d x ) 1 p , (1.13)

| f | L = sup , < x < | f ( x ) | , 1 p < , (1.14)

并且积分区域R将会被省略,而不产生混淆概念的后果,用 H m ( m 0 ) 用来表示m阶Sobolev空间,其范数为

f = ( k = 0 m x k f 2 ) 1 2 , = 0 = L 2 . (1.15)

2. 准备工作

非线性扩散方程,其中一自相似解 ( v ¯ , u ¯ ) ( x 1 + t ) 满足

{ v ¯ t u ¯ x = 0 , p ( v ¯ ) x = α u ¯ , x , t > 0 (2.1)

命题2.1 (参考 [24] )对于上述该自相似解,当 x ± 时,有 v ( ± , t ) = v ± ,让 ζ = x 1 + t ,则有以下性质:

1 ) | v ¯ ( ζ ) v | ζ > 0 + | v ¯ ( ζ ) v + | ζ < 0 C | v + v | e μ ζ 2 , (2.2)

2 ) | x r t l v ¯ | C | v + v | ( 1 + t ) r + 2 l 2 e μ ζ 2 1 + t , r + l 1 , r , l 0 , (2.3)

3 ) | x r t l u ¯ | C | v + v | ( 1 + t ) r + 2 l + 1 2 e μ ζ 2 1 + t , r , l 0. (2.4)

其中C是大于0的常数, r , l 是正整数,该引理在下面第三部分的证明中起了一定的作用。

3. 定理1.1的证明

为了证明定理1.1,我们首先构造扰动方程

ϕ ( x , t ) = v ( x , t ) v ¯ ( x + x 0 , t ) m ( x , t ) ,

ψ ( x , t ) = u ( x , t ) u ¯ ( x + x 0 , t ) u ^ ( x , t ) ,

ω ( x , t ) = ω ( x , t ) ,

其中

u ^ u e α t + x m t ( η , t ) d η , m 0 ( x ) d x = 1 , m ( x , t ) u + u α m 0 ( x , t ) e α t ,

我们得到以下扰动方程

{ ϕ t ψ x = 0 ψ t + [ p ( ϕ + v ¯ + m ) p ( v ¯ ) ] x + u ¯ t + α ψ = 0 ω t = ( A ( v ) ω x v ) x 4 μ r ( v ) v ω , (3.1)

初值满足

( ϕ ( t , x ) , ψ ( t , x ) , ω ( t , x ) ) | t = 0 = ( ϕ 0 ( x ) , ψ 0 ( x ) , ω 0 ( x ) ) = ( v 0 ( t ) v ¯ 0 ( t ) m 0 ( t ) , u 0 ( t ) u ¯ 0 ( t ) u ^ 0 ( t ) , ω 0 ( x ) ) (3.2)

我们令

y ( x , t ) = x ϕ ( η , t ) d η ,

注意到

y ( , t ) = ϕ ( η , t ) d η = 0 = y ( , t ) ,

那么根据(3.1)1,可以得到

y x = ϕ , (3.3)

y t = x ϕ t ( η , t ) d η = x ψ x ( η , t ) d η = ψ , (3.4)

则扰动方程(3.1)转换为以下形式

{ y t t + [ p ( y x + v ¯ + m ) p ( v ¯ ) ] x + α y t 1 α p ( v ¯ ) x t = 0 , ω t = ( A ( y x + v ¯ + m ) ω x y x + v ¯ + m ) x 4 μ r ( y x + v ¯ + m ) ( y x + v ¯ + m ) ω , (3.5)

初值为

( y , ω ) ( x , 0 ) = ( y 0 ( x ) , ω 0 ( x ) ) , y t ( x , 0 ) = y 1 , (3.6)

由于定理1.1的条件以及(3.2)知

y ( x , 0 ) H 2 ( ) , y t ( x , 0 ) H 1 ( ) , ω ( x , 0 ) H 2 ( ) , (3.7)

通过以上分析,我们知道,只要证明问题(3.5)~(3.6)的解的整体存在性,即可证明问题(1.1)~(1.2)的解的整体存在性,因此问题转化为如何证明(3.5)~(3.6)解的整体存在性并得到相应的衰减估计。

定义初值问题(3.5)~(3.6)的解空间如下:

X ( 0 , T ) = { y L ( 0 , T ; H 3 ) , ω L ( 0 , T ; H 2 ) L 2 ( 0 , T ; H 1 ) , y x L 2 ( 0 , T ; H 2 ) }

定理3.1 (整体存在性)设初值满足(3.7),则存在与 u ± , v ± 无关的 δ ˜ ,使得当

| u + u | + | v + v | + N 0 δ ˜ (3.8)

时,初值问题(3.5)~(3.6)的解 ( y , ω ) X [ 0 , ] 中全局存在且唯一。

根据偏微分方程的经典理论,初值问题(3.5)~(3.6)的解局部存在且唯一。

定理3.2 (局部存在性)假设 ( y 0 , ω 0 ) H 3 × H 2 ,且 N 0 ε 0 2 ,则存在正常数 T 0 ,其中

N 0 = y 0 H 2 + ω 0 H 2 + y 1 H 1 , (3.9)

使得问题(3.5)~(3.6)在 [ 0 , T 0 ] 上有唯一解 ( y , ω ) X ( 0 , T 0 ) ,且满足

( y , ω ) H 2 2 C 1 ( ( y 0 , ω 0 ) H 2 2 + y t ( 0 ) H 1 2 ) , (3.10)

为了研究解的全局存在性和解的大时间渐近行为,我们只需要证明解的一致先验估计。为此,我们做如下先验假设:

先验假设3.1对于时间 T > 0 ,设 ( y , ω ) 是问题(3.5)~(3.6)问题的解,并且满足

N ( T ) = : k = 0 1 ( 1 + k ) k t k y ( , t ) 2 + k = 0 2 ( 1 + k ) k x k y ( , t ) 2 + k = 0 2 ( 1 + k ) k k ω ( , t ) 2 C ε 2 , (3.11)

其中C是大于0的常数, ε 是依赖于初值和波的强度的某个充分小的正常数。

根据先验假设3.1以及Sobolev引理,有

y L C y L 2 1 2 y x L 2 1 2 C δ ( 1 + t ) 1 4 , y x L C y x L 2 1 2 y x x L 2 1 2 C δ ( 1 + t ) 3 4 , (3.12)

同样地,

ω L C δ ( 1 + t ) 1 4 , ω x L C δ ( 1 + t ) 3 4 .

在先验假设3.1下,可证明如下先验估计:

定理3.3 (先验估计)设 ( y , ω ) X ( 0 , T ) 是初值问题(3.5)~(3.6)的解,则存在与 u ± , v ± 无关的 δ 2 ,使得 | u + u | + | v + v | + N 0 δ 2 时,下面估计对 t [ 0 , T ] 成立

N ( t ) 2 + 0 t ( y x H 1 2 + y t H 1 2 + ω x H 1 2 ) d t C 2 ( N 0 2 + δ 2 ) . (3.13)

其中 N ( t ) = ( y , ω ) H 2 + y t H 1

下面我们将通过对 ( y , ω ) 进行低阶、高阶以及时间衰减估计来证明(3.13)式。

证明过程由下面几个引理构成。

引理3.1设 ( y , ω ) X ( 0 , T ) 是初值问题(3.5)~(3.6)的解,那么有以下估计

( y 2 + y x 2 + y t 2 ) ( x , T ) d x + 0 T ( y t 2 + y x 2 ) d x d t C N 0 2 , (3.14)

证明:

在(3.5)1两端同时乘以y,并在 [ 0 , T ] × R 上进行积分,得到

α y 2 2 d x | 0 T + y y t d x | 0 T 0 T y t 2 d x d t 0 T y x [ p ( y x + v ¯ + m ) p ( v ¯ ) ] d x d t + 0 T 1 α y x p ( v ¯ ) t d x d t = 0 (3.15)

由于 p < 0 ,当 δ 足够小时,利用Cauchy-Schwarz不等式,由上式可得

y 2 ( x , T ) d x + 0 T y x 2 d x d t = O ( 1 ) [ δ 2 + y t 2 ( x , T ) d x + 0 T y t 2 d x d t ] (3.16)

在(3.5)1两端同时乘以 y t ,并在 [ 0 , T ] × R 进行积分,得到

1 2 y t 2 d x | 0 T + 0 T α y t 2 d x d t 0 T y x t [ p ( y x + v ¯ + m ) p ( v ¯ ) ] d x d t 0 T 1 α p ( v ¯ ) x t y t d x d t = 0. (3.17)

q ( x , t ) 0 y x [ p ( η + v ¯ ( x , t ) ) p ( v ¯ ( x , t ) ) ] d η ,

q t y x t [ p ( y x + v ¯ ) p ( v ¯ ) ] + 0 y x [ p ( η + v ¯ ) p ( v ¯ ) ] v ¯ t d η ,

由于 p < 0 ,则存在 C ˜ 1 , C ˜ 2 > 0 ,使得

C ˜ 2 ( y x ) 2 q C ˜ 1 ( y x ) 2

从而

0 T y x t [ p ( y x + v ¯ ) p ( v ¯ ) ] d x d t = q d x | 0 T + 0 T 0 y x [ p ( η + v ¯ ) p ( v ¯ ) ] | v ¯ t | d η d x d t = q d x | 0 T + 0 T O ( 1 ) y x 2 | v ¯ t | d x d t C ˜ 1 y x 2 ( x , T ) d x C ˜ 2 y x 2 ( x , 0 ) d x + O ( 1 ) δ 0 T y x 2 d x d t , (3.18)

根据(3.13),(3.17)满足

( y x 2 + y t 2 ) ( x , T ) d x + 0 T y t 2 d x d t = O ( 1 ) δ 2 + O ( 1 ) δ 0 T y x 2 d x d t , (3.19)

则结合等式(3.16)得到基本估计

( y 2 + y x 2 + y t 2 ) ( x , T ) d x + 0 T ( y t 2 + y x 2 ) d x d t C N 0 2 , (3.20)

引理3.2设 ( y , ω ) X ( 0 , T ) 是初值问题(3.5)~(3.6)的解,那么有以下估计

ω 2 ( x , T ) d x + 0 T ω 2 d x d t + 0 T ω x 2 d x d t C ω 0 2 d x , (3.21)

证明:

将方程(3.5)2乘上 ω ( x , t ) ,得到

ω ω t = ω ( A ( y x + v ¯ + m ) ω x y x + v ¯ + m ) x 4 μ r ( y x + v ¯ + m ) ( y x + v ¯ + m ) ω 2 , (3.22)

将(3.22)关于 ( t , x ) [ 0 , T ] × R 上积分,我们得到

1 2 ω 2 ( x , T ) d x + 0 T A ( y x + v ¯ + m ) y x + v ¯ + m ω x 2 d x d t + 4 0 T μ r ( y x + v ¯ + m ) ( y x + v ¯ + m ) ω 2 d x d t = 1 2 ω 0 2 d x , (3.23)

根据先验假设(3.10)以及引理2.1,存在一个大于0的常数 M ,使得

1 M v = y x + v ¯ + m M ,

那么在方程(3.23)中,我们有

0 T A ( y x + v ¯ + m ) y x + v ¯ + m ω x 2 d x d t C 0 T ω x 2 d x d t , (3.24)

4 0 T μ r ( y x + v ¯ + m ) ( y x + v ¯ + m ) ω 2 d x d t C 0 T ω 2 d x d t , (3.25)

根据(3.24)、(3.25),并结合(3.23),整理得到

R ω 2 ( x , T ) d x + 0 T ω 2 d x d t + 0 T ω x 2 d x d t C ω 0 2 d x , (3.26)

则有

ω 2 ( x , T ) d x + 0 T ω 2 d x d t + 0 T ω x 2 d x d t C N 0 2 , (3.27)

现在我们来证明 ω x 的高阶估计。

引理3.3设 ( y , ω ) X ( 0 , T ) 是初值问题(3.5)~(3.6)的解,那么有以下估计

ω x 2 ( x , T ) d x + 0 T ω x 2 d x d t + 0 T ω x x 2 d x d t C ω 0 x 2 d x , (3.28)

证明:

将方程(3.5)2对x求导,再乘上 ω x

ω x ω x t = ω x [ A ( y x + v ¯ + m ) ω x y x + v ¯ + m ] x x 4 ω x [ μ r ( y x + v ¯ + m ) ( y x + v ¯ + m ) ω ] x , (3.29)

将(3.29)关于 ( t , x ) [ 0 , T ] × R 上积分得到,有

1 2 ω x 2 ( x , T ) d x + 0 T ω x x [ A ( y x + v ¯ + m ) y x + v ¯ + m ω x ] x d x d t + 4 0 T ω x [ μ r ( y x + v ¯ + m ) ( y x + v ¯ + m ) ω ] x d x d t = 1 2 ω 0 x 2 d x , (3.30)

其中

0 T ω x x [ A ( y x + v ¯ + m ) y x + v ¯ + m ω x ] x d x d t = 0 T [ A ( y x + v ¯ + m ) y x + v ¯ + m ] x ω x ω x x + A ( y x + v ¯ + m ) y x + v ¯ + m ω x x 2 d x d t = 0 T [ ( y x x + v ¯ x + m x ) A ( y x + v ¯ + m ) y x + v ¯ + m ( y x x + v ¯ x + m x ) A ( y x + v ¯ + m ) ( y x + v ¯ + m ) 2 ] ω x ω x x d x d t + 0 T A ( y x + v ¯ + m ) y x + v ¯ + m ω x x 2 d x d t C 0 T ω x x 2 d x d t C ( ε + δ ) 0 T ω x x 2 d x d t C ( ε + δ ) 0 T ω x 2 d x d t , (3.31)

并且

4 0 T ω x [ μ r ( y x + v ¯ + m ) ( y x + v ¯ + m ) ω ] x d x d t 4 0 T ω x [ μ r ( y x + v ¯ + m ) ω ] x d x d t C 0 T μ r ( y x + v ¯ + m ) ω x 2 + μ r ( y x + v ¯ + m ) ( y x x + v ¯ x + m x ) ω ω x d x d t C ( ε + δ ) 0 T ω x 2 d x d t + C 0 T ω 2 d x d t , (3.32)

把(3.31)、(3.32),代入(3.30),取 ε , δ 适当小,可得

ω x 2 ( x , T ) d x + 0 T ω x 2 d x d t + 0 T ω x x 2 d x d t C ω 0 x 2 d x d t C ( ε + δ ) ω 2 d x d t , (3.33)

应用柯西不等式,并利用(3.21)有

ω x 2 ( x , T ) d x + 0 T ω x 2 d x d t + 0 T ω x x 2 d x d t C ω 0 x 2 d x + ω 0 2 d x , (3.34)

由初值条件以及引理3.2则有

ω x 2 ( x , T ) d x + 0 T ω x 2 d x d t + 0 T ω x x 2 d x d t C N 0 2 , (3.35)

引理3.4设 ( y , ω ) X ( 0 , T ) 是初值问题(3.5)~(3.6)的解,那么有以下估计

y x 2 ( x , T ) d x + y x t 2 ( x , T ) d x + y x x 2 ( x , T ) d x + 0 T y x t 2 d x d t + 0 T y x x 2 d x d t C δ 2 , (3.36)

证明:

对(3.5)1有关x求导,同时乘以 y x 得到

y x t t y x + [ p ( y x + v ¯ + m ) p ( v ¯ ) ] x x y x + α y x y x t 1 α p ( v ¯ ) x x t y x = 0

并在 [ 0 , T ] × R 上进行积分,得到

y x t y x d x | 0 T 0 T y x t 2 d x d t + 0 T [ p ( y x + v ¯ + m ) p ( v ¯ ) ] x x y x d x d t + α 2 y x 2 ( x , T ) d x α 2 y x 2 ( x , 0 ) d x 0 T 1 α p ( v ¯ ) x x t y x d x d t = 0 , (3.37)

其中

0 T [ p ( y x + v ¯ + m ) p ( v ¯ ) ] x x y x d x d t = 0 T p ( y x + v ¯ + m ) x y x x d x d t + 0 T p ( v ¯ ) x y x x d x d t = 0 T p ( y x + v ¯ + m ) ( y x x + v ¯ x + m x ) y x x d x d t + 0 T p ( v ¯ ) x y x x d x d t C 0 T ( y x x + v ¯ x + m x ) y x x d x d t + 0 T p ( v ¯ ) x y x x d x d t C 0 T y x x 2 d x d t , (3.38)

0 T 1 α p ( v ¯ ) x x t y x d x d t = 1 α p ( v ¯ ) x x y x t d x | 0 T + 1 α 0 T p ( v ¯ ) x x y x t d x d t C p ( v ¯ ) x x y x t d x | 0 T + C 0 T [ p ( v ¯ ) v ¯ x 2 + p ( v ¯ ) v ¯ x x ] y x t d x d t C p ( v ¯ ) x x y x t d x | 0 T + C 0 T y x t 2 d x d t , (3.39)

同样地对(3.5)1有关x求导并乘以 y x t ,在 [ 0 , T ] × R 上进行积分,得到

1 2 y x t 2 ( x , t ) d x | 0 T + 0 T [ p ( y x + v ¯ + m ) p ( v ¯ ) ] x x y x t d x d t + 0 T α y x t 2 d x d t 0 T 1 α p ( v ¯ ) x x t y x t d x d t = 0 , (3.40)

其中

0 T [ p ( y x + v ¯ + m ) p ( v ¯ ) ] x x y x t d x d t = 0 T p ( y x + v ¯ + m ) x x y x t d x d t 0 T p ( v ¯ ) x x y x t d x d t C 0 T ( y x x + v ¯ x + m x ) y x x t d x d t 0 T p ( v ¯ ) x x y x t d x d t C y x x 2 d x | 0 T 0 T ( v ¯ x x + m x x ) y x t d x d t 0 T p ( v ¯ ) x x y x t d x d t C y x x 2 d x | 0 T + C 0 T y x t 2 d x d t , (3.41)

0 T 1 α p ( v ¯ ) x x t y x t d x d t C 0 T y x t 2 d x d t , (3.42)

将(3.38)~(3.39)代入(3.37),同样地将(3.41)~(3.42)代入(3.40),结合得到

y x t 2 ( x , T ) d x + y x x 2 ( x , T ) d x + 0 T y x t 2 d x d t y x t 2 ( x , 0 ) d x + y x x 2 ( x , 0 ) d x + y x 2 ( x , 0 ) d x , (3.43)

则得到以下了基本衰减估计:

y x 2 ( x , T ) d x + y x t 2 ( x , T ) d x + y x x 2 ( x , T ) d x + 0 T y x t 2 d x d t + 0 T y x x 2 d x d t C N 0 2 , (3.44)

接下来,我们进行有关时间衰减估计的证明。

引理3.5设 ( y , ω ) X ( 0 , T ) 是初值问题(3.5)~(3.6)的解,那么有以下估计

( 1 + T ) ( y x 2 + y t 2 ) ( x , T ) d x C δ 2 , (3.45)

证明:

在(3.5)1两端同时乘以 ( 1 + t ) y t ,并在 [ 0 , T ] × R 上进行积分,得到

0 T ( 1 + t ) y t y t t d x d t + 0 T ( 1 + t ) y t [ p ( v ¯ ) y x ] x d x d t + 0 T α ( 1 + t ) y t 2 d x d t = 0 T F ( 1 + t ) y t d x d t (3.46)

这里

F = 1 α p ( v ¯ ) x t + [ p ( v ¯ ) y x ] x [ p ( y x + v ¯ + m ) p ( v ¯ ) ] x

将(3.46)整合得到

1 2 ( 1 + T ) y t 2 ( x , T ) d x + 1 2 ( 1 + T ) p ( v ¯ ) y x 2 ( x , T ) d x + α 0 T ( 1 + t ) y t 2 d x d t = α 0 T p ( v ¯ ) y x 2 d x d t + 1 + t 2 p ( v ¯ ) y x 2 v ¯ t d x + α 0 T ( 1 + t ) y t F d x d t C δ 2 + α 0 T ( 1 + t ) y t F d x d t , (3.47)

其中

α 0 T ( 1 + t ) y t F d x d t = 0 T ( 1 + t ) y t p ( v ¯ ) x t d x d t α 0 T ( 1 + t ) y t [ p ( y x + v ¯ + m ) p ( v ¯ ) p ( v ¯ ) y x ] x d x d t , (3.48)

在上式中

0 T ( 1 + t ) y t p ( v ¯ ) x t d x d t = 0 T ( 1 + t ) y t [ p ( v ¯ ) v ¯ x t + p ( v ¯ ) v ¯ t v ¯ x ] d x d t = 0 T ( 1 + t ) y t p ( v ¯ ) v ¯ x t d x d t + 0 T ( 1 + t ) y t p ( v ¯ ) v ¯ t v ¯ x d x d t C δ 0 T ( 1 + t ) 1 2 y t e μ ζ 2 1 + t d x d t + C δ 2 0 T ( 1 + t ) 1 2 y t e 2 μ ζ 2 1 + t d x d t C δ 0 T ( 1 + t ) y t 2 d t + C δ 2 , (3.49)

α 0 T ( 1 + t ) y t [ p ( y x + v ¯ + m ) p ( v ¯ ) p ( v ¯ ) y x ] x d x d t = α 0 T ( 1 + t ) y t [ p ( v ¯ ) m + o ( 1 ) ( y x + m ) 2 ] x d x d t = α 0 T ( 1 + t ) y t p ( v ¯ ) v ¯ x m d x d t + α 0 T ( 1 + t ) y t p ( v ¯ ) m x d x d t + C 0 T ( 1 + t ) y t [ ( y x + m ) 2 ] x d x d t C 0 T ( 1 + t ) y t 2 m d x d t + C 0 T ( 1 + t ) y x t ( y x + m ) 2 d x d t C 0 T ( 1 + t ) y t 2 d t + C 0 T y x t 2 d x d t , (3.50)

同样地将(3.48)~(3.50)代入(3.47),则得到以下了基本衰减估计:

( 1 + T ) ( y x 2 + y t 2 ) ( x , T ) d x C N 0 2 , (3.51)

引理3.6设 ( y , ω ) X ( 0 , T ) 是初值问题(3.5)~(3.6)的解,那么有以下估计

( 1 + T ) ω 2 ( x , t ) d x + 0 T ( 1 + t ) ω x 2 d x d t C ω 0 x 2 d x , (3.52)

证明:

对方程(3.5)2两边同时乘以 ( 1 + t ) ω ,并在 [ 0 , T ] × R 上积分,得到

( 1 + T ) ω 2 ( x , T ) d x + 0 T ( 1 + t ) A ( y x + v ¯ + m ) ω x 2 y x + v ¯ + m d x d t + 4 0 T ( 1 + t ) μ r ( y x + v ¯ + m ) ( y x + v ¯ + m ) ω 2 d x d t = 1 2 ω 0 2 d x + 1 2 0 T ω 2 d x d t , (3.53)

由于 y x + v ¯ + m 有正下界,且 A ( v ) , μ r ( v ) 为大于0的光滑函数,所以有

0 T ( 1 + t ) A ( y x + v ¯ + m ) ω x 2 y x + v ¯ + m d x d t C 0 T ( 1 + t ) ω x 2 d x d t , (3.54)

4 0 T ( 1 + t ) μ r ( y x + v ¯ + m ) ( y x + v ¯ + m ) ω 2 d x d t C 0 T ( 1 + t ) ω 2 d x d t , (3.55)

把(3.54)、(3.55),代入至(3.53),整理得到

( 1 + T ) ω 2 ( x , T ) d x + 0 T ( 1 + t ) ω x 2 d x d t + 0 T ( 1 + t ) ω 2 d x d t C ω 0 2 d x + C 0 T ω 2 d x d t , (3.56)

由初值条件及引理3.2得到

( 1 + T ) ω 2 ( x , T ) d x + 0 T ( 1 + t ) ω 2 d x d t + 0 T ( 1 + t ) ω x 2 d x d t C ω 0 2 d x , (3.57)

引理3.7设 ( y , ω ) X ( 0 , T ) 是初值问题(3.5)~(3.6)的解,那么有以下估计

( 1 + T ) ω x 2 ( x , T ) d x + 0 T ( 1 + t ) ω x x 2 d x d t C ω 0 x 2 d x . (3.58)

证明:

将方程(3.5)2对x求导,两边同时再乘以 ( 1 + t ) ω x ,得到

( 1 + t ) ω x ω x t = ( 1 + t ) ω x [ A ( y x + v ¯ + m ) ω x y x + v ¯ + m ] x x 4 ( 1 + t ) ω x [ μ r ( y x + v ¯ + m ) ( y x + v ¯ + m ) ω ] x , (3.59)

将(3.59)关于 ( t , x ) [ 0 , T ] × R 上分部积分,得到

1 2 ( 1 + T ) ω x 2 ( x , T ) d x + 0 T ( 1 + t ) ω x x [ A ( y x + v ¯ + m ) ω x y x + v ¯ + m ] x d x d t + 4 0 T ( 1 + t ) ω x [ μ r ( y x + v ¯ + m ) ( y x + v ¯ + m ) ω ] x d x d t = 1 2 ω 0 x 2 d x + 1 2 0 T ω x 2 d x d t , (3.60)

其中

0 T ( 1 + t ) ω x x [ A ( y x + v ¯ + m ) y x + v ¯ + m ω x ] x d x d t = 0 T ( 1 + t ) [ A ( y x + v ¯ + m ) y x + v ¯ + m ] ω x x 2 + ( 1 + t ) [ A ( y x + v ¯ + m ) y x + v ¯ + m ] x ω x ω x x d x d t = 0 T ( 1 + t ) [ A ( y x + v ¯ + m ) y x + v ¯ + m ] ω x x 2 d x d t + 0 T R ( 1 + t ) ( y x x + v ¯ x + m x ) A ( y x + v ¯ + m ) y x + v ¯ + m ω x ω x x d x d t 0 T ( 1 + t ) ( y x x + v ¯ x + m x ) A ( y x + v ¯ + m ) ( y x + v ¯ + m ) 2 ω x ω x x d x d t C 0 T ( 1 + t ) ω x x 2 d x d t C ( ε + δ ) 0 T ( 1 + t ) ω x 2 d x d t C ( ε + δ ) 0 T ( 1 + t ) ω x x 2 d x d t , C δ 0 T ( 1 + t ) ω x x 2 d x d t C ( ε + δ ) 0 T ( 1 + t ) ω x 2 d x d t . (3.61)

并且

4 0 T ( 1 + t ) ω x [ μ r ( y x + v ¯ + m ) ( y x + v ¯ + m ) ω ] x d x d t 4 0 T ( 1 + t ) ω x [ μ r ( y x + v ¯ + m ) ω ] x d x d t C 0 T ( 1 + t ) μ r ( y x + v ¯ + m ) ω x 2 d x d t + 0 T ( 1 + t ) ( y x x + v ¯ x + m x ) μ r ( y x + v ¯ + m ) ω ω x d x d t C 0 T ( 1 + t ) ω x 2 d x d t C ( ε + δ ) 0 T ( 1 + t ) ω x 2 d x d t C 0 T ω 2 d x d t C δ 0 T ( 1 + t ) ω x 2 d x d t C 0 T ω 2 d x d t , (3.62)

把(3.61)、(3.62),代入至(3.60),整理得到

( 1 + T ) ω x 2 ( x , T ) d x + 0 T ( 1 + t ) ω x x 2 d x d t C ω 0 x 2 d x + C 0 T ( 1 + t ) ω x 2 d x d t , (3.63)

结合公式(3.28)和(3.41),我们有

( 1 + T ) ω x 2 ( x , T ) d x + 0 T ( 1 + t ) ω x x 2 d x d t C N 0 2 , (3.64)

引理3.8 设 ( y , ω ) X ( 0 , T ) 是初值问题(3.5)~(3.6)的解,那么有以下估计

( 1 + T ) 2 ω x x 2 ( x , T ) d x + 0 T ( 1 + t ) 2 ω x x x 2 d x d t C ω 0 x x 2 d x . (3.65)

证明:

将方程(3.5)2对x求二次导,再乘上 ( 1 + t ) 2 ω x x ,则有

( 1 + t ) 2 ω x x ω x x t = ( 1 + t ) 2 ω x x ( A ( y x + v ¯ + m ) ω x y x + v ¯ + m ) x x x 4 ( 1 + t ) 2 ω x x [ μ r ( y x + v ¯ + m ) ( y x + v ¯ + m ) ω ] x x , (3.66)

将(3.66)对 ( t , x ) [ 0 , T ] × R 上积分得到,有

1 2 ( 1 + T ) 2 ω x x 2 ( x , T ) d x + 0 T ( 1 + t ) 2 ω x x x [ A ( y x + v ¯ + m ) ω x y x + v ¯ + m ] x x d x d t + 4 0 T ( 1 + t ) 2 ω x x [ μ r ( y x + v ¯ + m ) ( y x + v ¯ + m ) ω ] x x d x d t = 1 2 ω 0 x x 2 d x + 0 T ( 1 + t ) ω x x 2 d x d t , (3.67)

其中

0 T ( 1 + t ) 2 ω x x x [ A ( y x + v ¯ + m ) ω x y x + v ¯ + m ] x x d x d t = 0 T ( 1 + t ) 2 ω x x x [ A ( y x + v ¯ + m ) y x + v ¯ + m ω x x x + 2 ( A ( y x + v ¯ + m ) y x + v ¯ + m ) x ω x x + ( A ( y x + v ¯ + m ) y x + v ¯ + m ) x x ω x ] d x d t = 0 T A ( y x + v ¯ + m ) y x + v ¯ + m ( 1 + t ) 2 ω x x x 2 d x d t + 2 0 T ( A ( y x + v ¯ + m ) y x + v ¯ + m ) x ( 1 + t ) 2 ω x x ω x x x d x d t + 0 T ( A ( y x + v ¯ + m ) y x + v ¯ + m ) x x ( 1 + t ) 2 ω x ω x x x d x d t C 0 T ( 1 + t ) 2 ω x x x 2 d x d t + 0 T ( 1 + t ) 2 ω x x ω x x x d x d t + 0 T ( 1 + t ) 2 ω x ω x x x d x d t , (3.68)

并且

4 0 T ( 1 + t ) 2 ω x x [ μ r ( y x + v ¯ + m ) ( y x + v ¯ + m ) ω ] x x d x d t C 0 T ( 1 + t ) 2 ω x x [ μ r ( v ¯ ) ( y x + v ¯ + m ) ω ] x x d x d t C 0 T ( 1 + t ) 2 ω x x [ ( y x + v ¯ + m ) x x ω + 2 ( y x + v ¯ + m ) x ω x + ( y x + v ¯ + m ) ω x x ] d x d t C 0 T ( 1 + t ) 2 ω x x 2 d x d t + C 0 T ( 1 + t ) 2 ω x ω x x d x d t + C 0 T ( 1 + t ) 2 ω ω x x d x d t , (3.69)

则将上式整理,得到

1 2 ( 1 + T ) 2 ω x x 2 ( x , T ) d x + C 0 T ( 1 + t ) 2 ω x x x 2 d x d t C ω 0 x x 2 d x + H , (3.70)

其中

H = C 0 T ( 1 + t ) 2 ω x 2 d x d t C 0 T ( 1 + t ) 2 ω ω x x d x d t 0 T ( 1 + t ) 2 ω x x ω x x x d x d t 0 T ( 1 + t ) 2 ω x ω x x x d x d t C 0 T ( 1 + t ) 2 ω x x x 2 d x d t + C 0 T ( 1 + t ) ω x x 2 d x d t + C 0 T ω x 2 d x d t + C 0 T ω 2 d x d t C 0 T ( 1 + t ) 2 ω x x x 2 d x d t + C δ 2 , (3.71)

根据先验估计以及上述引理,有

( 1 + T ) 2 ω x x 2 ( x , T ) d x + 0 T ( 1 + t ) 2 ω x x x 2 d x d t C ε 2 . (3.72)

下面我们开始定理3.1的证明。

δ = min { ε 2 C 1 + 1 , ε ( 2 C 1 + 1 ) ( 2 C 2 + 1 ) }

则在定理1.1条件下, N ( t ) δ ε 2 ε 0 2

由局部存在性定理3.2,在 [ 0 , T 0 ] 上,令 N ( t ) = ( y , ω ) H 2 + y t H 1 ,问题(3.5)~(3.6)的解存在,且满足

N ( t ) C 1 , N 0 C 1 δ ε 2 C 1 + 1 < ε , (3.73)

从而满足先验假设3.1,由先验估计定理2.2知,在 [ 0 , T 0 ] 上,

N ( t ) C 2 N 0 C 2 δ ε 2 C 1 + 1 < ε < ε 0 , (3.74)

特别地,上式对 T = T 0 也成立。以 T = T 0 为初始,局部存在性定理,可延拓到 [ T 0 , 2 T 0 ] ,依次进行下去,可以将解延拓到 T = +

根据引理3.1至引理3.8我们已经证明先验估计(3.13)成立,由解的局部存在性定理和先验估计(3.13),我们通过以上可以证明(3.5)~(3.6)的解整体存在,进一步证明了(3.1)~(3.2)的解整体存在,并满足如下估计:

( u u ¯ u ^ ) ( , t ) L 2 + ( v v ¯ m ) ( , t ) L 2 + ω ( , t ) L 2 = O ( 1 ) δ ( t + 1 ) 1 2 , (3.75)

对于衰减估计(3.49),我们也可以推广到更高阶导数,如此我们有

x k ( u u ¯ u ^ ) ( , t ) L 2 + x k ( v v ¯ m ) ( , t ) L 2 + x k ω ( , t ) L 2 = O ( 1 ) δ ( t + 1 ) k 2 , k = 1 , 2 (3.76)

基于上式,从Sobolev引理得到逐点估计

sup t [ 0 , T ] ( u u ¯ u ^ L + v v ¯ m L + ω L ) = O ( 1 ) δ ( t + 1 ) 1 4 . (3.77)

从而定理1.1得到证明。

参考文献

[1] Hsiao, L. and Liu, T.-P. (1992) Convergence to Nonlinear Diffusion Waves for Solutions of a System of Hyperbolic Conservation Laws with Damping. Communications in Mathematical Physics, 143, 599-605.
https://doi.org/10.1007/BF02099268
[2] Hsiao, L. and Liu, T.-P. (1993) Nonlinear Diffusive Phenomena of Nonlinear Hyperbolic Systems. Chinese Annals of Mathematics, 14, 465-480.
[3] Hsiao, L and Li, T.T. (1983) Global Smooth Solution of Cauchy Problems for a Class of Quasilinear Hyperbolic Systems. Chinese Annals of Mathematics, Series B, 4, 109-115.
[4] Li, T.T., Qin, T.H. (1985) Global Smooth Solutions for a Class of Quasilinear Hyperbolic Systems with Dissipative Terms. Chinese Annals of Mathematics, Series B, 6, 199-210.
[5] Nishida, T. (1978) Nonlinear Hyperbolic Equations and Related Topics in Fluid Dynamics. Publications Mathematiques D’Orsay 78.02, Department de mathématique, Paris-Sud.
[6] Zheng, Y.S. (1996) Global Smooth Solution to the Adiabatic Gas Dynamics System with Dissipation Terms. Chinese Annals of Mathematics, Series A, 17, 155-162.
[7] Nishihara, K. (1996) Convergence Rates to Nonlinear Diffusion Waves for Solutions of System of Hyperbolic Conservation Laws with Damping. Journal of Differential Equations, 131, 171-188.
https://doi.org/10.1006/jdeq.1996.0159
[8] Nishihara, K. (1997) Asymptotic Behavior of Solutions of Quasilinear Hyperbolic Equations with Linear Damping. Journal of Differential Equations, 137, 384-395.
https://doi.org/10.1006/jdeq.1997.3268
[9] Nishihara, K., Wang, W.K. and Yang, T. (2000) Lp-Convergence Rate to Nonlinear Diffusion Waves for p-System with Damping. Journal of Differential Equations, 161, 191-218.
https://doi.org/10.1006/jdeq.1999.3703
[10] Zhu, C.J. and Jiang, M.N. (2006) Lp-Decay Rates to Nonlinear Diffusion Waves for p-System with Nonlinear Damping. Science in China Series A, 49, 721-739.
https://doi.org/10.1007/s11425-006-0721-5
[11] Zhao, H.J. (2000) Asymptotic Behaviors of Solutions of Quasilinear Hyperbolic Equations with Linear Damping, II. Journal of Differential Equations, 167, 467-494.
https://doi.org/10.1006/jdeq.2000.3793
[12] Jiang, M. and Zhu, C.J. (2009) Convergence to Strong Nonlinear Diffusion Waves for Solutions to p-System with Damping On quadrant. Journal of Differential Equations, 246, 50-77.
https://doi.org/10.1016/j.jde.2008.03.033
[13] Dafermos, C. (1995) A System of Hyperbolic Conservation Laws with Frictional Damping. In: Casey, J. and Crochet, M.J., Eds., Theoretical, Experimental, and Numerical Contributions to the Mechanics of Fluids and Solids, Birkhäuser, Basel, 294-307.
https://doi.org/10.1007/978-3-0348-9229-2_16
[14] Hsiao, L. and Tang, S.Q. (1995) Construction and Qualitative Behavior of the Solution of the Perturbed Riemann Problem for the System of One-Dimensional Isentropic Flow with Damping. Journal of Differential Equations, 123, 480-503.
https://doi.org/10.1006/jdeq.1995.1170
[15] Huang, F.-M. and Pan, R.-H. (2003) Convergence Rate for Compressible Euler Equations with Damping and Vacuum. Archive for Rational Mechanics and Analysis, 166, 359-376.
https://doi.org/10.1007/s00205-002-0234-5
[16] Li, H.-L. and Saxton, K. (2003) Asymptotic Behavior of Solutions to Quasilinear Hyperbolic Equations with Nonlinear Damping. Quarterly of Applied Mathematics, 61, 295-313.
https://doi.org/10.1090/qam/1976371
[17] Pan, R.-H. (2006) Darcy’s Law as Long-Time Limit of Adiabatic Porous Media Flow. Journal of Differential Equations, 220, 121-146.
https://doi.org/10.1016/j.jde.2004.10.013
[18] Yu, S.H. (1999) Zero-Dissipation Limit of Solutions with Shocks for Systems of Hyperbolic Conservation Laws. Archive for Rational Mechanics and Analysis, 146, 275-370.
https://doi.org/10.1007/s002050050143
[19] Zhu, C.J. (2003) Convergence Rates to Nonlinear Diffusion Waves for Weak Entropy Solutions to P-System with Damping. Science in China Series A: Mathematics, 46, 563-575.
https://doi.org/10.1007/BF02884028
[20] Nishida, T. (1978) Nonlinear Hyperbolic Equations and Related Topics in Fluid Dynamics. Université de Paris-Sud, Département de Mathématique, Orsay.
[21] Duyn, C.T. and Van Peletier, L.A. (1977) A Class of Similarity Solutions of the Nonlinear Diffusion Equation. Nonlinear Analysis: Theory, Methods & Applications, 1, 223-233.
https://doi.org/10.1016/0362-546X(77)90032-3
[22] Matsumura, A. (1977) Global Existence and Asymptotics of the Solutions of the Second-Order Quasilinear Hyperbolic Equation with the First-Order Dissipation. Publications of the Research Institute for Mathema, 13, 349-379.
https://doi.org/10.2977/prims/1195189813
[23] Matsumura, A. and Mei, M. (1999) Convergence to Travelling Fronts of Solutions of the P-System with Viscosity in the Presence of a Boundary. Arch. Ration. Archive for Rational Mechanics and Analysis, 146, 1-22.
https://doi.org/10.1007/s002050050134
[24] Huang, F.-M., Mei, M. and Wang, Y. (2011) Large Time Behavior of Solutions to N-Dimensional Bipolar Hydrodynamic Models for Semiconductors. SIAM Journal on Mathematical Analysis, 43, 1595-1630.
https://doi.org/10.1137/100810228