1. 引言
对于偏微分方程组在航空航天、燃烧与爆破理论、大气物理等诸多方面的研究一直都是人们关注的重点,有着十分广泛的应用。而其自身在发展中也极大的丰富了数学研究的理论,充分推动了数学学科理论的发展。理想流体,即忽略流体中的粘性因素以及热传导因素,常在气体动力学的研究中出现。而描述理想流体进行运动的具体方程为包括质量守恒、能量守恒、动量守恒方程在内的欧拉方程。因此本文中所研究的退化双曲欧拉方程具有十分重要的物理意义,是研究流体力学的重要内容。黎曼问题作为双曲守恒律方程组的具有标度不变量的初值的柯西初值问题,是分片的具有常状态的初值问题,其初始状态的间断线汇聚在原点处成为几何奇异点。一般黎曼问题所取形式为有限个常状态,取该形式的意义在于可将该方程组进行自相似变换。针对欧拉方程二维黎曼问题的研究由 [1] 开创,利用广义特征分析法和数值实验提出了关于问题解的多个猜想。非线性退化双曲守恒律方程组在初值充分光滑或充分小的情况下,该方程组问题的解随时间的发展可能会出现间断点,即该问题的局部解易证明。而该问题解的间断现象存在于各种物理现象中,比如在流体力学中与激波的产生与传播相对应。对于带有初值的二维黎曼问题,研究其整体光滑解时,可利用特征分解的方法将局部解进行延拓从而得到问题的全局解。对含跨音速激波的可压缩流问题的研究已具有较为成熟的理论方法,例如Chen和Zhang在 [2] 中定义了二维守恒律方程组中的双曲性、凸性、中心简单波等,对二维欧拉方程黎曼问题的研究发展提供了重要准备。Zhang和Zheng在 [3] 中研究了某个单个二维守恒律方程黎曼问题,从中得到五类解。针对多方气体二维绝热欧拉方程黎曼问题,其中该问题的初值为四片常数,分布在四个象限中,文章中将问题归于初始间断线上发出的四道平面基本波之间的相互作用。利用广义特征分析法得到解的性质。具体方法参考 [4] [5] [6] [7] [8] 。在文章 [9] 中针对一类二维黎曼问题利用特征分解法及数值模拟法给出了十二类分析与猜想。在超音速解的范围中,Dai和Zhang [10] 利用特征分解法得到静止的气体在向真空扩散的过程中存在超音速解。Yang和Zhang在 [11] 中证明了压差方程在两个平面稀疏波的相互作用下所存在的整体经典解,将 [10] 中Dai和Zhang的结果进行了拓展延伸。
Song和Zheng [12] 为了得到在某一固定条件下退化到声速线上的解的存在性,在二维压力梯度系统中构造了解的半双曲片这样的新型波。由二维超音速流在绕开拐角向真空区域扩散所引起的多方气体向真空扩散问题,该问题与一中心稀疏波和一平面稀疏波相互作用有关。在数学研究中,该问题即为某一二维可压缩欧拉方程的一个Goursat问题,且该方程组具有自相似性。Feng和Xiao [13] 证明了一类压强定律的二维可压缩欧拉方程半双曲片解的存在性。本文与 [14] 中考虑欧拉方程Goursat型边值型问题利用特征分解法及自举法的思想二维黎曼问题在真空边界处的整体解的方法不同,主要根据特征角建立估计从而证明了存在解且解具有整体光滑性。
本文主要研究一个二维欧拉方程组的整体光滑解的存在性。具体分为四个部分:第二节建立压力P的特征分解;第三节利用解的边界值估计得到局部解的存在性;第四节主要进行解的梯度估计;第五节利用先验估计结合三四节得到退化双曲方程组Goursat问题整体解的存在性。
2. 提出问题
我们首先引入以下二维可压缩欧拉系统
(1)
其中P为压力,
为速度。将系统(1)进行自相似变换,令
,有
(2)
且(2)可简化为
(3)
(3)的正负特征线定义为如下方程的特征曲线
(4)
进一步地,类似于 [14] [15] 中的方法,我们定义特征角的概念:
(5)
且引入
,
。则有
(6)
证:根据特征角概念(4),得到
(7)
则有
(8)
即证明结论。
为了证明发生在两个连接不连续点以及两个稀疏波的相互作用的半双曲片,我们令
为三个实数且有
。令
,定义
,思考如下二维稀疏波
:
(9)
建立一个半双曲片区域:
在二维稀疏波
中,给出一个正特征线
,作切线的延长线BD变为恒定状态。其中
两点为声速点。同时给出严格凸的负特征线
,其中点C为声速点。该曲线三角形ABC为退化双曲区域,在该区域内求解,
为声速线。
3. 特征分解
下面给出关于P的特征分解。根据 [16] 中的方法,首先对特征角
进行限制,令
(10)
且给出方向导数
(11)
将压力P代入以上方向导数公式,利用P的方向导数
表示出特征角
的一阶方向导数。
3.1. 性质1
给出
以及
的一阶方向导数:
(12)
证:对(8)求导结合方向导数定义(11)可得
(13)
同理,根据以上方法可以得到
(14)
对(8)求导后结合(13)有
(15)
同理
(16)
得到(15)后,根据(13)中
,
之间的关系可得
,
的关系式:
(17)
同理可得
(18)
根据流角,马赫角与特征角之间的关系,得到
(19)
(20)
3.2. 性质2
计算压力P的二阶导数,得到P的特征分解如下:
(21)
证:首先根据方向导数定义可得
(22)
利用(11),结合(16)有
(23)
且有
(24)
则
(25)
同理
(26)
由(6)可以得到
与
之间的关系,根据
有解可知
有解。令
,有
(27)
且
。因此,
(28)
因此对方程组(2)及(10)的解可转化为
上的解。
4. 边界值估计及局部存在性
以下做P特征分解的先验估计:首先给出边界
,
上的条件。在B点的位置靠近A点的范围内,圆弧
弧长不超过1/4圆。假设对于凸负特征线在点C处的倾斜角
有
。
且有如下边界值条件:
4.1. 推论1
(P的二阶方程的齐次形式)
(29)
根据 [17] 中的证明可以得到针对退化双曲方程组(1)及边界值条件,有区域
上存在方程组的唯一
-解,其中
唯一依赖于边界
,
处
的
-范数。
4.2. 引理1
(不变三角形)对于任意
-解,有
(30)
证:由方程组(3) (10)的局部解可1得存在
使得在区域
中存在解。对于
,区域
中有
,
。
由于
(31)
由于
,结合(32)可得
。
根据(17)得
。由于
,结合(19)有
。
则结合(18) (19) (20),在区域
中有
(32)
5. 梯度估计
在第四节当中,首先利用连续性方法得到关于压力P的
的一致边界,其次利用压力P的特征分解,给出其常微分方程形式,从而建立梯度估计。
5.1. 引理2
(
的一致边界)在半双曲片区域ABC中,
的最大值是有界的,且
(33)
证:根据引理1结论
,
,则可知证明(34)即证明
(34)
只需证明对任意
,有
。当在B点附近区域
时,结论显然成立。否则一定有
上的一点
,在
内
,
。在点
处有
或
。不妨设在
处,满足
,
。由点
处开始作负特征线,结合(33)可知
。根据连续性可得在
的某邻域
内,有
成立。因此在该邻域内必存在点
,使得
成立,与假设矛盾。根据
的任意性可知引理2成立。
5.2. 引理3
(一阶梯度估计)若在区域
(
)内存在
-连续解,存在常数C且与
无关,有
(35)
证:根据引理2可得
(36)
结合(10)有
(37)
同理,有
。
6. 全局解
6.1. 引理4
若在区域
(
)内存在
-解,存在常数C且与
无关,有
(38)
证:结合(18)及(19)有
低阶项,
低阶项。 (39)
则可得
,
有界。对于
,利用(19)可得
低阶项。 (40)
因此
有界,同理有
有界可证。
利用(22)及(18) (19)有
(41)
同理
(42)
通过计算有
(43)
同理
。以上证明过程
同理成立。
根据所得方程组的局部存在性定理及先验估计,利用有限覆盖定理将局部解进行延拓从而得到如下定理。
6.2. 定理2
退化双曲方程组Goursat问题(1) (10) (11)在半双曲片区域ABC内存在整体光滑解。