双曲型方程组

1、(1)其中以及wx,0=w0xWix,tW2x,tWm(x,t)-二：x二f1wf2Wfm(w)flWi，W2，,Wmf2Wi，W2，,WmfmWi，W2，,Wmw0都是m维向量，分别叫做守恒变量、通量和初值。双曲型守恒律方程组1.拟线性双曲型守恒律方程初值问题拟线性双曲型守恒律方程组初值问题的一般形式为；wfw=0:t.x在初值问题(1)中，要求守恒变量具有紧支集，即:集合S=(xw(x,t)#。的闭包为有界集。从而存在常数R>0,使得：当x>R时,w(x,t)=0；要求通量满足归零条件，即：当|w|t0时,(2)ux,0=x例2:一维气体力学方程组初值问题(3)wx,0i=w0

2、x例1:单个守恒律方程初值问题其中（对于完全气体,利用关系式1E二-11-2p-u2,其中的¥是比热比）fi:uPu2pu<u(E+P)J2.守恒性w1www3w2上+U-1)W3w1w2Wi|P：uE2w1w22w13-7w；丁3-1w322w1w?w3Wi上积分，得任取A<-R,B>R,将方程（1）在区间lA,BB；wAftBjfwdxdx=0A.:xdBr-.-/awdx'f|wB,t-f|wA,t”=0但根据A、B的选取，有f-w（A,t）=0,f-w(B,t)=0所以有dbwdx=0dtA令At-co,Bt",就得到dx,tdx=0(4)

3、dwdt-二此式对任意的t>0都成立。这表明,f：w（x,t）dx是一个与时间无关的常数,或者说，这个量是一个守恒量。对任意的t>0,再将（4）时在区间10,t上积分，得.wx,dxd>；=0也就是"bo-bowx,tdx-iwx,0dx=0从而有wx,tdx=wx，0dx二-bew-OO0xdx.wx,tdx=-bew-oO0xdx-t03.非守恒形式定义通量的Jacobi矩阵A；fAw=wvf1外.%:w1.:w2-：wm开2，开2.:w1.:w2.：wma二+.a二九：fm；：fm“1:w2-：wm则方程组（1）可写成.:t：w：x:wftA卫二0二x(5)称

4、为方程组（1）的非守恒形式。例3:对于一维气体力学方程组（3),-13f"2w22w39-3-u23-u-13EE-12-1u3-u-3u2u:；-2引入音速0=收,则p11二一1EP：2E1所以尸三二,不十u2,于是上述Jacobi矩阵最终可写成；-120103-2u23-u二2十4、2V-1123-22一cu-1u4,双曲型方程组对于方程（5）,以及相应的守恒形式（1）,定义：若矩阵A的所有特征值都是实数，并且矩阵A是可对角化的，则称方程（5）是双曲型的。（注）关于“矩阵A是可对角化的“，有以下几种等价的描述：矩阵A是可对角化的：矩阵A可以通过相似变换，变换成一个对角矩阵D;存在

5、可逆矩阵R，使得R，AR=D是对角矩阵；矩阵A存在m个线性无关的特征向量。事实上，以上述m个线性无关的特征向量为列，就得到矩阵R。例4:对于单个守恒律方程（2）,矩阵A成为标量a=a（u）=f'（u）,所以其特征值也是a5.不变性考虑因变量的变换，设Wi(V1,V2,，Vm)；W2(M,V2,，)w=w(v)=,Wm(Vi,V2,，Vm)则方程组(5)可写成w：v:f二w二vc=0v二t：w二v二x定义因变量变换的Jacobi矩阵'电加>>>>>当'cVi5V20Vm二2刖2二2P=6VlCV2CVmcv9+I-孤孤£Wm<

6、CViCV2CVm/则上式就是：vZPAP=0.trx三vft如果矩阵P可逆，就是/vPAP一二0.x记B=P,AP,就是vv八B二0,t二x经过因变量变换，双曲型方程中的矩阵A变成了矩阵B。上述推导表明,B与A相似。如果A能够相似于对角矩阵，则B也相似于对角矩阵。也就是说，双曲型方程组经过因变量变换还是双曲型方程或者说,方程组的双曲性质在因变量变换下具有不变性。例5:对于一维气体力学方程组（3）,如果用原变量表示，可写成cPcPu:t；:xfuftfuu一；:x二0；x1fp-=0p:：x或写成就是.:pFpr2fucu：c2=0.:txxufp)2u+0a山01PfP)色u=0exuP0B

7、=0u0Pc2uvvcB=0-t：x对于一维气体力学方程组3）,如果将守恒变量w换成原变量v,例2中的Jacobi矩阵A变成了这里的矩阵Bo根据双曲性方程组的不变性，这两个矩阵是相似的，从而有相同的特征值。显然，计算矩阵B的特征值更容易些。事实上，由KI-B=,：-u000PU'-u-c1-U'-u“c'-uc0这样，我们很容易就求出了矩阵B,同时也是矩阵A的三个特征值%=uc,九2=u,%=u+c6 .特征线特征线是双曲型方程的重要性质，我们先通过一个线化的模型方程对此展开讨论。考虑对流方程初值问题O=UXa+=O,audx考虑x-t平面上由万程一dta确定的一族平行

8、直线(因为这里的a为常数)。沿着这族直线中的任意一条,x=x(t)成为t的函数，所以方程的解u=u(x(t),t)也是t的函数。所以有du桃udx札u=tt+=+adt抖xdt抖xu=常数也就是说，沿着这族曲线中的每一条,程的特征线。从而有以下结论：方程的解都是常数。这样的曲线称为该方对于双曲型方程，沿着方程的特征线，方程的解为常数如图，对x-t平面(上半平面,t>0)上的任何一个点p(x,t),设过P点的特征线与x轴相交于Q(xo,0)点，根据上面的性质，就有u(x,t)=%=%=中)。x=at+bTx0=a?0b=b=x-at所以u(x,t)=(x-at)这样，利用特征线的性质，我们

9、实际上已经得到了初值问题的解。对于标量守恒律方程工?f(u)抖记a(u)=f),可有非守恒形式桃u八+a(u)=0抖()x此时，上面的所有推导依然成立。沿一条特征线，由于u为常数，所以a(u)也是常数，因此特征线仍是直线。但对于不同的特征线，u的变化导致特征线斜率a(u)的变化，所以此时的特征线不再是一族平行直线，这是非线性方程与线性方程之间的重要区止匕外，虽然初值问题的解仍可写成u(x,t)=(x-a(u)t)但由于右端隐含了u,上式并没有给出显式的解。7 .Riemann不变重再考虑拟线性双曲型守恒律方程组。直接考虑非守恒形式：wftA-X根据方程组的双曲性，有R"AR=D。其中

10、，矩阵D是对角矩阵，其对角线元素就是矩阵A的特征值九j(j=1,2,L,m);矩阵R是可逆矩阵，其列向量行向量rj就是对应的右特征向量，即Ar）=%口；还可以证明，逆矩阵R的LT就是对应的左特征向量，即有iTa=%iTo于是，用iT左乘方程可将方程组解耦，成为m个独立的方程式=0(j=1,2,L,m)此时，引入变量玛，满足dRjfTdw。则上式成为.:Rj；Rj十%=0(j=1,2,L,m).t.x于是，根据上一节的特征线理论，可定义（方程组的）第j族特征线Cj：-jjftj而沿着其中的每一条特征线，都有Rj=常数o因此将Rj称为（方程组的）第j个Riemann不变量。（注）由lTA=九jlT

11、,取转置，得ATlj=4lj,所以矩阵A的左特征向量实际上就是矩阵AT的右特征向量。在因变量变换w=w（v）下，考虑行向量qT=lT=lTP，则而根据双二v曲型方程组的不变性，B=P，AP，有qTB=l:PB=lTPPAP=l:AP=jl：P=jlTP=jq所以qT恰好就是矩阵B的左特征向量。止匕时,TT-WTTdRj=ljdw=lj-dv=ljPdv=qdv这表明,Riemann不变量是双曲型方程组在因变量变换下的不变量。10例6:对于一维气体力学方程组,设qT=«0小）是矩阵对特征值%2%2：c2的左特征向量，则)q=0o(BT-uI,2=0,有方程组q1q22q=_cq3qi=

12、一cv不171cv二，0，一从而dR2=q；dvcvVil.一dP+dp=Cvd(lnp)-'d(lnP)=Cvdpp)所以R2=cvIn-丫=S就是烯Jo对特征值K13=u不c,由（B,3T-(u+c)I)为=0,有方程组-c-Cq3二一In与Plq1q2g3J0、04，解得q=0,q3=+4T3=0,1,£11T,1,dR13=q3dv=du一dp，：cdp=p2dc2对绝热等嫡流动，有羽=二上，所以p一1c2dc=dc-1从而dR,3,f丁2dc=du+cJ最终得到Fc8 .Riemann问题考虑一个截面积不变的无限长直管，在直管内的x=0处有一金属膜片。在膜片左侧（x

13、<0）的管道内充满了一种气体，其密度为P1,压力为p1;而在膜片右侧（x>0）的管道内则充满了另一种气体，其密度为P2,压力为p2。设在t=0时刻，金属膜片被突然打破，考虑直管内两种气体在破膜后的流动。这一问题称为激波管问题。可以证明，激波管问题存在自相似解，即：所有的流动参数都是组合自变量u=uG),p=p。，p=Pf)对于激波管问题的解，有如下的描述（解的波系结构）12在t>0时刻，将产生两个波，分别在两种气体中传播。第一种气体中的波记作L,向左传播；第二种气体中的波记作R,向右传播。根据两种气体的状态不同，波L和R,既可以是激波（强间断），也可以是中心稀疏波（弱间断）。

14、第一种气体经过波L,达到新的常数状态p=p3,p=P3;而第二种气体经过波R,也达到新的常数状态p=p4,p=p4。由于是无粘流，两种气体不会掺混，两者之间始终保持着一个接触面，称为接触间断。但是这个接触间断不是静止不动的，而是移动的（也就是会传播的），这就形成了第三个波，记作M。纵上所述，激波管问题的解在x-t平面内被上述三个波分成了四个区域，每个区域内的解都是常数。这些常数解的区域由各种波分隔开，如下图所示。Ox激波管问题可以用一维气体力学方程组间断初值问题:w:fw0二二x:00：x二(6)二txwLwx,0=wR来描述，称为Riemann问题。139 .弱解作为一个典型的例子，Riem

15、ann问题表明,双曲型方程组的解并不总是连续的、光滑的，可以产生间断。这样的间断解显然不能用传统的方式来描述，需要引进弱解的概念。记x-t平面的上半平面（t30）为区域W,考虑函数集合v蛹C1(),suppv有界任取v?V,则由suppv的有界性，一定存在T>0和有限区间b，使得v在矩形区域Sb备T的三条边界x=a、x=b和t=T以及边界以外的区域上等于零。用v乘以方程组（1）的两边，并在矩形区域薪b'箱上积分，得f,Cvdxdt=0x+分部积分14-v-+fzdxdt+?-fvdt+wvdx=0这里，G表示矩形区域T的边界。在组成G的四条边上，v只在t=0这一条边上有可能不等于

16、零，而沿着这条边，dt=0,所以上式中的边界积分变成蝌fvdt+wvdx=bbw(x,0)v(x,0)dx=?w0(x)v(x,0)dx从而有执w+融抖v-fzdxdt+x+b?w0(x)v(x,0)dx=0?w+触抖v-fzdxdt=b?w0(x)v(x,0)dx现在令b?vq?fzdxdt=?w0(x)v(x,0)dxx.'''''根据以上推导,可以给出双曲型守恒律方程组初值问题弱解的定义。定义（弱解）:若存在分片连续的函数w（x,t）及相应的通量函数f(u),使得vq?-?w杆+f-ddxdt=?w0(x)v(x,0)dx,x-?则称w为初值问题（

17、1）的弱解15激波关系式不失一般性，假设初值函数w0(x)在x=x0处有一个间断。这使得初值问题(1)的解w在区域W上形成了一条间断线C:x=x(t)(t30)这条间断线将矩形区域薪b'露T分成了W.和WR两个部分。abx在这种情形下，我们重复上面的推导。&f+dvdxdt=0Wa啊x+分部积分+fV-zdxdtx+?-槽Lfvdt+wvdxv+-zdxdt+x+?-禧Rfvdt+wvdx=016v-4-±xdt+x+-fvdt+wvdx+Clx0aw(x,0)v(x,0)dxH涉#V丘vf,质+fzdxdt+-fvdt+wvdx+W?>x+CV吸crb?w(x

18、,0)v(x,0)dxv.jdxdt+x+-fvdt+CLwvdx+x0aw0(x)v(x,0)dxv=dxdt+x+-fvdt+wvdx+CR'xw0x0(x)v(x,0)dx=0式中Cl和Cr分别表示沿间断线c两侧的积分。显然，按照逆时针方向积分时，cl和cr方向相反。令cl?c,cr?c,上式可合并写成v.ddxdt+x+(x)v(x,0)dx；f；vdt+:w；vdx=0其中(w)=wR-wL，f=fR-fL表示物理量在跨过间断线c时的跳跃量。再令a?,bT?，就有J?射持"丫、?/、/c、-蝴-?wIT+f岁出+?w0(x)v(x，0)dxvdt+iw)vdx=017如果w是初值问题（1）的弱