一类二维单个守恒律方程的Riemann问题
2021-08-30
来源:年旅网
2010年12月 Dec..2010 应用数学与计算数学学报 C0MM.ON APPL.MATH.AND C0MPUT 第24卷第2期 Vo1.24 NO.2 一类二维单个守恒律方程的Riemann问题 张泓知 ,。 盛万成。 摘要本文研究一类二维单个守恒律方程的Riemann问题.用广义特征分析方法研究 了这类方程,给出了基本波的分类,解决了初值为两片常数的二维Riemann问题,给出 了Riemann解. 关键词二维Riemann问题,单个守恒律方程,广义特征分析方法,基本波,激波, 疏散波 The TWO Dimensional Riemann Problems for a Class of Scalar Conservation Laws Zhang Hongzhi , Shang Wancheng。十 Abstract This paper is conserned with the two dimensional Riemann problems of a class of SCalar conservation laws.By use of the generalied characteristic analysis method, the classification of elementary waves is given,and the Riemann solutions for the two pieces constant initial data are constructed. Keywords Two.dimension Riemann problem,SCalar conservation laws,generalized characteristic analy’sis method,elementary wave,shock wave,rarefaction wave. 1 引 ..口 -l 一 考虑二维单个守恒律方程 ut+_厂( ) +9( “) =0 及其Riemann问题 “(0, ,Y)=uo(O) (o≤0≤2丌) 其中0=arctan 是x, )平面中的极角,uo(O)是分段常数函数,l厂( ),g(札)满足下 列条件: (1),(u),9( )∈C。( ); (2)f ( )≠0,夕 ( )≠0; 收稿日期:2007年12月3日. 1上海应用数学与力学研究所,上海,200072;Shanghai Institute of Applied Mathematics and Mechanics Shanghai 200072,China 2.海南大学信息学院应用数学系,海南,570228;Department of Applied Mathematics,College of Infor mation Science&Technology Hainan University,Hainan 570228,China 3.上海大学数学系,上海200444;Department of Mathematics,Shanghai University,Shanghai 200444 China 十通讯作者Corresponding author 2 应用数学与计算数学学报 24卷 (3)(H~ 。)( ) >0. 条件(1)一(3)表明单个守恒律方程是真iE_x维的,且其基本曲线非凸,tt=u0是它 的一个拐点.为简单起见,又不失一般性,假设: m ㈤ ( 咖,,(u)>。,乱( ) >。. (H) 文[3—4]中给出了方程(1.1)的Cauchy问题的存在性和唯一性定理,这表明了 Riemann问题(1.1)和(1.2)仅依赖于( ,叩)=(x/t,v/t)的解唯一.因此,我们可以考虑 (∈, r/)平面中的解.此时,(1.1)变成 ∈u∈一qu +_厂(n)∈+夕(Ⅱ)叩=0, 而(1.2)得到在无穷远点的边界值 l imtg/∈ f2+rl2一十 (∈,叩)∈ 2, (1.3) (∈,叩)=u0( ), ∈[0,2丌】 .(1.4) uo(O)的间断在(∈,叩)平面的无穷远点的邻域将产生含有基本波:稀疏波和激波(包含 半接触间断)的波.因此,构造(1.3)和(1.4)的Riemann解,实质上是研究基本波的相 互作用并最终在一起匹配. [2,8—12]等在不同的条件下,得到T.维单个守恒律方程的Riemann问题的解. 我们将以f(uJ=U2/2,g(u)=U4/4为例,在假设(H)下讨论基本波的分类及Riemann 问题. 2广义特征分析 2.1特征线及其奇性曲线 对于光滑解,方程(1.3)可以写为 (, (u)一f)u∈+(gl(札)~叩)u =0 因此,(1.3)的特征方程是 tf 一叩一 ) 对该系统积分,我们得到特征线 } 其中 和C是任意常数.以上特征线在(∈,叼)一平面上的投影也同样称为特征线这些 特征线显然都是半射线,其终点是(2.2)的积分曲线的奇点(∈叩):(,,( ),9,(u))上. .,2期 张泓知,等:一类二维单个守恒律方程的Riemann问题 3 这些奇点对应到下列形式的特征线 1f =, :-9 ,((Ⅱ)t )£ (2.4) 由于Riemann数据(1.2)在原点处的奇性,每一条这些线都是一个特征面的边界. 当u取Ul时的每一条特征线都是指向(∈,叩)一平面中的奇性点(∈,卵)=(.厂 (乱1),gt( )) 1(见图1).(2.3)式所有的奇性点的集合是一条参数曲线 : c……。。 5 同文[9,11],称其为基本曲线或特征线的奇性曲线,记为I1(“):叩= (∈).它具有以下 性质. 引理1在条件(H)下,特征线的奇性曲线r(u):叩=叩r( ),是单调递增的,并 且(∈,卵)=(, (0),gl(0))全A0是它唯一的拐点. 证明根据条件(H)和r(u)的定义,我们有塞=笋 >0,即r(札)单调, 髻:( ) du f-G ; ;南 'tIt>0=0., (2.6) <0. 由此引理,在r(u)和tt轴之间存在一个1—1映射,而一条特征线上的“可以定 义为这些特征线与奇性曲线r(u)的交点.因此,在不引起异义的情况下,可以把I、(札) 看成为“轴.如果某区域内所有的特征线直到奇性曲线c(u)才相交,那么我们就称其 为疏散波,记为R. 2.2间断及其奇性曲线 考虑(2.1)与(2.2)的间断解.令s:r'/=叩(∈)为解的一个间断,其上点P(∈o,叩o)两 侧u的极限值分别为U一和u+,那么按照Rankine—Hugoniot(R-H)关系,有: dT]: ∈(2.一7) 其中 。 J +'+, 一t上一 = i 二 上一 f上一 ,’ 夕 +,u+, 一= 一 , (2.8) J_~钆一 以下,点(∈,叩)=( , )记为A+,一等等.根据R-H关系(2.7),在点P(∈o,卵0) 的切线,必然过点 +,一.因此,我们规定s的方向为其切线向量指向 +,一,或者简单 的说,s指向点 +.一. 4 应用数学与计算数学学报 24卷 g- . 与特征线的奇性曲线类似,(2.7)的积分曲线同样有它自己的奇性点(∈,叩)=( +, , ).如果固定S一侧的状态u一,而使另一侧的状态札变化,则所有的奇性点构成 一了一个以u为参数的曲线: ,u一 U- (一。C<“<+。。) (2.9) .这条曲线称为间断线的奇性曲线,记为Fs(u, 一)或者卵=叩r。(f).它同样有与特征线 的奇性曲线相类似的性质. 引理2对于任意固定的“ ∈ ,奇性曲线rs(u, ):叩=叩r。(∈)有如下性质: (I)Fs(u,“1)单调递增; (2)rs(“,n )与r( )在r(“ )=(.厂 (J“ ),I9 (“1))处相切; (3)除了在 1=, ( “1)点外, . l>77 (∈), s ∈>, (0) 1<叩r( ), <,,(0) (4)V“l< 2,除了在0点相交外,Fs(u,721)与Fs(u,%t2)在Fs(u2,U1)=rs(ul, 2)= ( 胁, ) Am处有唯一的交点.且 , ( )< 。<, (u2) (5)V 1<“2,r( 1)位于Ps(u,U2)在点PS(Ul,U2)处的切线上;F(u2)位于Fs(u,U1)在 点rs(u2,“1)处的切线上. 证明 drh,s一: d∈ = d“/ d“ 2/ : .厂 (“)一 . 二 ! . f2.10) d∈。 (, )一 ( 一_, (n) . , )一咒 /)『 . (2.il) (1)根据条件(H),有 ( (u)一, , )(u一札 )>0,(9- (¨)一_q 、 )(u—u )>0, 所以, >0・ 一 ,(2)由(2.8),有lim "11— U1 , ( 1), linl U—}“1 一 :9 (u1).于是有 Fs(u,11,1)lu= ̄l:r(札)ku= ̄l, 且 lira …=lim …一 = … 2期 张泓知,等:一类二维单个守恒律方程的Ri。ma 问题 5 (3)由(2.11), -厂,,(u) /,9,,(“,) 夕 (,f上) d∈0 l=“一“ …lJ.i1m111 , U=Ul .7 一 _厂 ( ) =兰 (丽 而/ )J (2.13) 从而 (警一 ) _一一 1而1( ) < ul>0。 ̄. 再结合(2)的结果, 可知当,“ >0时,在∈:, (,“ )点的某个邻域(∈ :‘厂,(札 )除 外)内,『/r (∈)> (∈).从而有叩r ( )> (∈), >,,(0).事实上,若不成立,则必然 存在 =, ( 2)(不失一般性,假设 > ),使得叩r )=叩r( ).从而有 , 2): 。 , 9I 2)=9 , , (2.14) 其中,屯2是Fs(u,u1)对应于交点( 2,叩r( ))的参数.从而 l< < 2.另一方面, 我们有, d叩 l ,d叩 l 。 即 l厂 。)一 。 ≤ . (2.15)5J 由(2.14)和(2.15)式,我们可得 9 ( 2)一gl(,“2) 9 ( 2) .厂 ( 2)一.厂 ( 2) , (札2)‘ 根据Cauchy中值定理, ∈( “2, 2), s.t. 夕 (面)/,g (“2) 而这与(笋 ) >0矛盾; 同理可以证明卯 ( )<叩r(∈),V∈<, (0). (4)由Fs(u, )的性质(1),可知Fs(u," )与 一轴之间存在一个1—1映射, 而 值是由S的一侧上的一个特征线给出的.这就意味着本结论成立. (5)由下列式子,本结论显然成立: = 1( “2. ) , ( )~ 。 ’,i:1,2Fs 6 应用数学与计算数学学报 2.3熵条件 不失一般性,以下我们均假设it一<it+.为了获得唯一的分片光滑解,我们需 要Rankine—Hugoniot条件(2.7)和一个熵条件.根据Kruzkov的存在性和唯一性定 理,文[1,7]中的熵条件应为: H(k) k Vk∈( 一,u+). (2.16) 其中,日(“)= ,(札)+ 9(札),且n(e):( , )( 。+//'。一1)足S在P点的法向,其方 向由u一指向 z+.我们同样可以用几何形式来重写Rankine—Hugoniot条件(2.7)和熵 条件(2.16): if(P)・Prs(u+, 一)=0, if(P)・Fs(“,乱一)rs(札,乱+)≤0, (2.17) (2.18) (2.19) (P)・Ps(u,it+)rs( ,札一)≥0, 关系(2.17)一(2.19)意味着线段PFs( ,u一)是间断线s在尸点处的切线,Vu∈ ( ,it+),向量 (P)与向量Fs(u,— )rs( ,“一)之间的夹角是钝角,且Vu∈(it一,“+), 向量 (尸)与向量 ‘n 之间的夹角是锐角.特别地,当(2.18)或者(2.19) 式中的等式成立,则必为一个半接触间断. 如果间断线S:叩=叩(∈)满足Rankine—Hugoniot条件(2.7)和熵条件(2.16),则我 们称其为激波(包括接触间断,记为J).所以我们有以下的几何准则. 引理3在条件(H)下,一个连接u 和it2的间断线是激波当且仅当 (1)间断线的方向指向点Fs(u1,it2); (2)间断线两侧的特征线不指向外. 最后,我们把一个疏散波和一个激波分别划分为两类: R一兄土若 u与特征线的方向组成一个左手系(有手系); S=S士若法向向量与S的切向向量组成一个右手系(左手系).当它包含接触间 ,,断时,激波记为JS±或者 士. R=R士与S—s土统称为基本波. 3基本波的分类 当我们求解系统(1.3)的边值问题(1.4)时,假设 是初始间断线,在(∈,叩)一平 面无穷远点的某个临域内从无穷远点沿着 方向投射了一个波( 或S,记为 ), 1(it2)分别是“在 左(右)侧的极限.根据R士,s士, 土的定义及他们在引理1和 引理2中的性质,针对l厂(u)=?t2/2,_9(u)=214/4,我们可以把初始间断分为以下几类。 (∈,叼)一平面被分为4个角域Qi(札 ,仳2)(也简记为Qi),角度分别为Q (i=l,…,4),满 足E =27r,且Q 以射线 和 +1为边界(g5= 1).而下列每个图中 与 都互 z=1 相平行. 2期 张泓知,等:一类二维单个守恒律方程的Riemann问题 7 1.乱1<0< 2:根据点Fs(u1, 2)与直线F(U1)F(u2)的位置关系分为两个小类: 1)点Fs(“1,札2)在直线F( )r(札2)的下方,此时 的方向为FS( ̄tl,“2)F(札2), 的方 向为 , 在点F( 。)处与特征线的奇性曲线r(u)相切, 为r( )向间断 的奇性曲线Fs(-“,“ )引的一条切线(切点其实是两条间断的奇性曲线不同于Fs(it1,r“z)的 另一个交点,并且也不是r(“,)与Fs( ,U1)的公共切点F(u1),在这里f(u)=u /2,g(饥)= 乱 /4时,切点其实就是原点): ● L , 1 = ,l●●●●●-,、●●●●【/ Ql / / % // |Q‘ } ,一 / + % 一 2 / Q3 rfu 。 图1 ∈【21, ∈ 2, ∈Q3, (3.1) ∈Q4, 其中, 代表“ ∈( , z)在Q2或者Q4中依次满足: l: ‘厂 (u )一 , ’, i=1,2 2)点rs(u ,u。)在直线 的上方,此时 的方向为 , 的 方向为 ,%在点r( )处与特征线的奇性曲线r(u)相切, 为F(钆。)向 间断的奇性曲线rs( ,r“z)引的一条切线(在这里切点其实就是原点): JL吖 1 Q1 / / 3 Q3 图2 8 应用数学与计算数学学报 :JI + R , f , ∈【21 ∈Q。 ∈ 3 (3.2) 【R + s 其中“ 与情形(a)中的含义相同. 2. 。<0< :根据点I"s( ,n2)与直线 ∈Q 的位置关系分为两个小类: ,-I●I●●●J‘l●●●●1)点Fs(“ .u )在直线 曲线r(札)相切, 丽的下方,此时 在点r( )处与特征线的奇性 【 为r(-“z)向间断的奇性曲线Fs(u, z)引的一条切线(在这里切点 其实就是原点).%的方向为r( 1)rs(Ⅱ1, “z), +他+砘硫 的方向为I1( )rs( 1, ): 一, 一 一+ R Ir/ , Q1 / l 24 i/ 【22/ r 3 3 Q3 ∈Ql ∈ 2 ∈Q3 ∈ 4 (3.3) 2)点I1s(札 , z)在直线 的上方,此时 在点r( 。)处与特征线的奇性 曲线r(u)相切, 为r(“ )向间断的奇性曲线Fs(u, )引的一条切线(在这里切点 其实就是原点). 的方向为 , 的方向为 : . , 1 【2 / Q4 |k / Q3 3 4 2期 张泓知,等:一类二维单个守恒律方程的Riemann问题 9 ∈Q1, ∈Q2, ∈【23, (3.4) ∈ft4 3.it1,it2≥0或者 “1,札2≤0.这两种情况事实上就是无拐点的奇异曲线中的情形, 它们在[2,7,9,11—12]中已经详细讨论过了,在这里我们略过. 4 具有两片常状态的Riemann解的结构 在这一节里,我们将根据引理3指出的几何标准以及关于初始间断(3.1)~(3.4)的 分类,解出具有两片常状态Ri黔 + o= 一emann数据的Riemann问题(1.1),这其中包含了基本 波的相互作用. 假设初始平面被两条从原点出发的射线,划分成了两块,在每一块上分别是Rie— mann数据的一个常状态,表达式如下: 札。c臼 ={ :;茎 , 。 , 。≤ < z≤2丌 (4.1) 其中it1 it2是常数。我们只需要考虑 2<0<it1这种情形,而u1<0<扎2这种情形 是类似的。其它的情形,如:“ ,it。≥0以及u ,U2≤0就是特征线的奇性曲线无拐点 的情形,它们已经在[2,7,9,11—12]中讨论过了,本文略过. 根据第三节的关于基本波的分类讨论,我们可以知道,针对仳z<0<u 这种情 形,由Fs(ul,71'2)在直线F(u1)F(u2)上方和下方以及01,02所代表的两条初始间断线 的方向,可以具体划分为32种小的情形.由于本文长度所限,现从其中挑出几种典型 的情况予以讨论,其它的详细过程略过. 1.R十2+R ,及 +s .这是最简单的两种情况.其中R十2与R 在特征线的 奇性曲线r(r“)上相交,而s击与s 都指向同一个终点rs(ul,71'2)=( l 9i2). 2.基本情形 : 十 +2.我们假设初始间断线如图5所示,因此一个激波和一 个疏散波: =kx \u 2 \ U "it1 ‘\ \~ / 乱= f上1 乱 1 图5 10 应用数学与计算数学学报 24卷 {慧 , (“+) ( (∈一U )2≤“, ≤U1) 从无穷远点把两片状态Ul和U2分割开来,其中k是常数. 和R 相交于点 P=(‘厂 (u ), +k(f (“z)一 )),然后激波将会刺穿疏散波,从而构成一条新的间断 线叩=r『s(∈),从点P起满足如下Rankine—Hugoniot条件: = , )‘(uz 1)1 (4.2) 【 ∈ ,,(1 。)=夕 1+ (, 2)一 ). 因此,(4.2)的积分曲线取如下参数形式(71,为参数) f (“), U2≤u≤Ul, 叼s(“) u exp+(9 + (, (扎z)一.,i2))exp(一 。(s, )ds), 其中, ( ̄s。c ,u d )9,s, 。cs,乱1)ds、 )= , (4・3) (4・4) 与f12]中的证明类似,可以很容易地得到积分曲线有一个最小值的点札="/-tO(由r/s(uo)= 9 所确定的).当 ≤ o时,它是单调递减的,而当U≥UO时单调递增,并且总 是在特征线的奇性曲线r(u)的上方直至在点F(U1)=(I厂 ( ),9 (札1))相切结束.即: Vs(u)>9 (札)(U2<“<l/1)且T/S(U1)=夕 ( 1). 方程组(4.2)表明,间断线叩=叩s(∈)的切线总是指向间断的奇性曲线rs(u,U1) (U2≤“≤n1).且有: 孥=d∈0 “ 南一乱1 f ( )f一 >。一. 、。5 因此,(4.2)的积分曲线是凸的. 进一步,我们很容易就可以知道,除了F(u1)点是特征线从状态 =u 相切到间 断线之外,曲线两边的所有的特征线都是进入的,这表明间断线在点F(U1)处消失了. 因此,间断线叩=W(u)是容许的 3.基本情形B:R 十S .为了简单起见,我们假设初始间断曲线如图6所示, 则一个疏散波和一个激波: 2期 张泓知,等:一类二维单个守恒律方程的Riemann问题 11 _“2 u2 厂 — I J 钆2 \Y=CX 图6 f R :叩 = gl(u)+c(∈一.厂 (u)) ( 2≤u≤ 1), 【sf2:叼=g l+ (∈一 1), 从无穷远点把两片状态 “ 和11,z分割开来,其中c,k是常数.这两个波在点 尸1= 91(1“1)一g2/1+kf;1一cf (U1)k gt(u1)一c911 + C 一C (处开始相交,并构成一个新的间断曲线叩一叩s( ).与[1,7,12]类似,可以知道,激波将 首先刺穿疏散波R (U2< ≤“1),然后从点P2(∈一, ( ))处发出一个包络疏散波. 这个包络疏散波以一个新的半接触间断与疏散波R (722< ≤ ).这个新的半接触 间断是单调递减的、凹的、与上方接触,并且在P3( :f (五))点处与r( )相切结束, 与疏散波R 在F(u)上相连.其中, = (u1,U2)且 = (u1, “2). 4.其它的情形.通过把之前的结果进行组合,容易得到剩下的28种情况下的结 果: (1)点rS(ul,U2)在直线r( ̄1)r(u2)之上时 f + +JS22, f R +R + _2, 【s +JS ̄-,+R , 【R + s五+R+2, r R + s +R十2, r S珐+R +R十2, J R + s +R + s J s 砬+R21+R +JS22 l R云+t厂 +s , l sJ;+R + , L R云+r,s +JS1.+R十2 L sJ;+t?2-1+JSf,+R l2 应用数学与计算数学学报 24卷 (2)点rs('ul,7.12)在直线r(“1)r( 之下时: Jf +J s +R_S2, 1.+R 1【 +sf2, 5茜+R +JS22 f Js++R +R+2, J JSL+R +JS ̄-.+R I f, +R + , L!,s去+ + + f 五+Js5+R十2, j R + 5 +Js ̄+R . Jf R +R十R 2, 1【 R +s . R + +JsV.+ +Js22 I R +JS+x+s12, 【兄云+ s +月 + s . 由于本文的长度所限,具体的讨论在此略过 参考文献 [1]Chang T.&Hsiao L.Riemann Problem and Interaction of Waves in Gas Dynamics[J].Pitman Itlonographs,Surveys Pure App1.Math.,1989,41. f 2 IChen G.Q.,Li D.&Tan D.Structure of Rienlann solutions for 2-Dimensional 8Calar conser— ration laws[a、 Di乳【"antial Equations 1996 127:124—147 [3]Conway E&Smoler J Global solution of the Cauchy problem for quasilinear first—order equations in several variables[J].Co:ulm.Pure Apple.Math.1966,XlX:95—105. [4】Guckenheimer J.Shocks and rarefactions in two space dimensions[J1.Arch.Rational Mech. Aha1.,1975 59:281—291. [5]Hiao L.&Kilingenberg C.The structure of the solution for the two—dimensional Riemann problem.Preprint,Heidelberg University,1986. [6】Kruzkov s.N.First order quasilinear equations with several independent variables[J].Math. USSR Sb..1970.10:271—243. [7]Li J.,Zhang T.&Yang S.Two—dimensional Riemann Problem in Gas Dynamics/J].Pitman Monographs.Surveys Pure AppI.Math.,1999,98. I8l Lindquist W.B.The scalar Riemann problem in two spatial dimensions:Piecewise smoothness of solutions[J].sjAM J.Math.Ana1..1986.17:l178—1197. 9] Sheng Wancheng.The Two—dimensional Riemann problem for scalar conservation laws[J] D Equs..2002.38:534—559. 10】Wagner D.The Riemann problem in two space dimentions for a single conservation law[J SIAM  ̄lath.Anal,,1983,38:534—559. 11 Zhang P.&Zhang T.Genralized characteristic analysis and Guckenheimer structure[J].Dif- ferential Equations,1999,152:409—430. [12]Zhang T.&Zheng Y.Two—dimensional Riemann problem for a single conservation law[J] Trans.Amer Math.Soc.,1989,312:589—619.