描写气体运动的基本方程是欧拉方程,它由质量、动量和能量三个守恒律组成,它的最大特点和困难在于解中会出现间断现象,冲击波就是一种压缩性的间断。1858年,黎曼紧紧抓住了间断现象这一特点,提出并解决了欧拉方程一种最简单的间断初值问题即初值为含有一个任意间断的阶梯函数,被后人称为黎曼问题。黎曼构造出了它的四类解,它们分别由前、后向疏散波和前、后向冲击波组装而成。并利用相平面分析方法给出了此四类解的判别条件。黎曼的这一工作开创了“微分方程广义解”概念及“相平面分析”方法之先河,具有极大的超前性。黎曼以其敏锐的洞察力和巨大的创造力为非线性双曲型守恒律的数学理论奠定了第一块基石。
对于守恒双曲型方程(1)的分片光滑弱解
∫0∞∫∞-∞[wtu+wxf(u)]dxdt+∫∞-∞w(x,0)φ(x)dx=0 (21)
其中ω(x,t)是具有致密支集的试函数.设在(x,t)平面上的不连续解,沿不连续曲线满足RH关系:
f(uR)-f(uL)=S(uR-uL) (22)
其中dx/dt=S.uR,uL分别为在不连续右方及左方函数值.由Olenik证明过解存在的熵条件
(f(u)-f(uL))/(u-uL)≥S≥(f(u)-f(uR))/(u-uR).(u∈(uL,uR)) (23)
定义1 若式(23)中严格不等式成立时,则称此不连续为激波.
定义2 若式(23)中等号成立时,则称此不连续为接触不连续.
方程式(1)的黎曼问题及其求解,设方程式(1)在标量情形下对于下列初值问题:
u(x,0)=uL,(x<0),u(x,0)=u
R,(0<x) (24)
的解称为方程(1)的
黎曼问题的解.这时在一维标量情形下我们只讨论具有下列形式的解:(见图1)在(x,t)平面上存在OA(f′(
L))及OB(f′(
R))方向.沿OA、OB是激波或接触不连续.而在扇形AOB中为连续解(稀疏波).在X
*OA中u=u
L.在BOX中u=u
R.但在AOB中每一个以O为原点的半放射线上,即在x/t=常数方向上u(x,t)=u(x/t)=常数.这种类型的解称为稀疏波.因此这里有两个问题要解决:1)在AOB中的稀疏波如何求得?2)OA、OB方向如何决定?现分别讨论如下.
图1 黎曼问题的解
Fig.1 The solution of the Riemannian problem
1)在AOB中满足隐式方程f′(u)=x/t的u=u(x/t),当f″≠0时为(1)的解.(即u=(f′)-1(x/t)为(1)的解,其中(f′)-1表示f′的逆函数).
证 由f′(u)=x/t微分之得f″(u)ux=1/t,f″(u)ut=-x/t2.相加得f″(ut+f′ux)=0,故当f″≠0时,u为解.
2)应用几何方法求OA、OB方向.(见图2)
|
图2 不连续方向的确定
Fig.2 The determination of the direction of the
discontinuous
先讨论OA方向如何决定(OB也同理)在(u,f)平面上作图形f=f(u).(uL,f(uL))为图形上一点.我们的目的是求 L使 L( L,f( L)),P(u,f(u)),pL(uL,f(uL))满足熵条件(23).而且PL, L方向为f′( L)=OA方向.(其中u在 L,uL中的任一值)在这样要求下,作法是在pL点作曲线的切线,使在 L点同曲线相切.故得
f( L)-f(uL)=f′( L)( L-uL) (25)
这时 L点切线的斜率等于直线pL L的斜率.这样所得的f′( L)值及 L,即为OA的方向及OA右边的不连续函数值.这里若曲线在切线的下方,可证条件(23)满足.这时为激波或接触不连续,对于式(1)为方程组情形,由于情况复杂,不在此讨论.
|