Convergence of an Upstream Finite Volume Scheme for a Nonlinear ...

3 downloads 0 Views 222KB Size Report
convergence of the solution given by an upstream finite volume scheme towards ... scheme for the time discretisation and a finite volume scheme with upstream.
       !"#$ &%'(")#*,+-/.01* 2 &435(""0  7689 :9(";=?)0 @; ACBD @;E0 )(; &F  G. HJILK@MLNPO-QSRUT$VXW$Y[ZIL\NP]^]U_a`cT$b dfegN"hLij1ISklT$VnmSRUhLe ow-pJx1qsy{rsz{t$|2u2}2v"~r axf€{x1z{2x}P‚y{ƒf€{xsz„‚y†…Xz„‚‡;ˆ9‡X‰ z{2xˆa‡;ˆPŠ‹‚ˆ2x$…X€"~"Œ[xs€{[‡$Š‹‚ƒxsŽ"|…Xz„‚‡ˆ’‘[“J”9•"–˜—a™˜šJ›œ™‘[ž{žlŸ¡ 1‚ˆ¢ £/¤¥¢ £¦¨§J©ª‚z{ «›1‚¬;‚yµxsˆ … ‚« ˆ"‚‚¬;z„xs^…­ˆŠ®ˆ2ƒs‡‡ˆ2ˆ2}a}Pxs‚z„ƒf‚€{‡;x$ˆ#…Xy„‘¯‚ˆ ™{«°‚±† ‰˜žJ|2ˆaŸ²ƒsz„‘a‚‡;³ˆ ™{° ž¨‰˜€{‚‡ˆ·¢ £@¢ £D¤§z{©G‡ 2¢ xs£ €{xl¸J¹@š´ˆ‚yµxs…/º"Œ"‰¶Š‹|a‚ƒnˆ2‚ƒszz„‚»¼‡ˆ|"Š‚xs‰˜€{€µ‡;y{·¡ƒn2¢xs£@· ¤gx@¥‚y ¢ £|2¦y{xs} z{‡ ‰˜‡¢ £l€®¤lz{ay{xl|aƒ†z„‚· z{xl…X}"z ‚y{}Pƒs€{‚¬xsz„š‚y†Ÿ²…$z„‚  ‡ˆ…Xˆ2‡$} ‰ z{ƒs2‡;x@ˆ"¬xsxsŽ"€ |…Xxsz„ˆa‚ƒf‡;xlˆ¯§L‡X‰J…$ˆaz{2}x@…ly{‡$z{€„Š‚|a^…Xz„ˆ ‚‡;« ˆ |"Š^…X‚€@¬;xs·½ˆ xs"y{ ~‰˜…$‡ˆ €ª|2z{2Œax@y{z{y{€{Œx­…X…Xz„·!^…sŠJ¿ˆ"}"‚‚z{y{xlƒs€{¬xs‡Xz„‚Š‚y†|2…$· z„‚x‡ˆSy{ƒn¸l2¾ xs· ˆa}2xµxsz{€ ‡$…l© |a…X€{y{}2|yª…­ŠLz{2y{z†x@…$|2PˆP‹Š‹‚‚Ž"z~ |2xgƒs‡;xfˆaˆ"}"z{€{‚z„‡‚Œ"‡~ ˆS§œ©ª©®x­xG…XÀ Œ2€{y{‡$‡$¬Š‚|axµz„z{‚‡;aˆ x « z{‡ z{2xl« xsŽ"|…Xz„‚‡ˆS¸ Á   -Ò)0l@; w-xlƒs‡ˆ2y„‚}axf€laxs€{xlz{2x@‰˜‡$ЋЂ‡X©ª‚ˆ « ˆ2‡ˆ"Š‹‚ˆax­…X€"~"Œ[xs€{[‡$Š‹‚ƒ1xfŽ"|[…$z„‚‡ˆ‚ˆ z„©®‡ y{Œ[…$ƒsx}P‚· xfˆay„‚‡;ˆays§L©ª‚z{1‚ˆ"‚z„^…­Š®ƒs‡ˆ2}P‚z„‚‡;ˆ¯Ä Å ‘[“{™Æ¯±{Ç{žS”#}P‚¬œ™˜š™Æ¯±Ç{ž{›œ™‘¯™Æ¯±{ÇȞ{ž{ž,Ÿ9 "± ÆÊÉË¢ £l¤±²Ç‚ˆ¢ £¦¨± ™{Í$ž ‘¯™ÆS±† ž4Ÿ9‘2³P™Æžn±ÌÆÊÉË¢ £l¤± Ï Ï © axf€{xl‘ “ }axfˆa‡;z{xsy@z{2x@}2Ïxs€„‚¬;…Xz„‚¬;x ‡X‰¼‘1©ª‚z{€{xsy{Œ[xsƒsz z{‡ dž§"}"‚¬ Ÿ ΠςÐS¤ Ñ Ò §2© 2xs€{x Ò }axfˆa‡;z{xsy z{axlŒÑ …X€{z„^…sŠJ}2xs€„‚¬;…Xz„‚¬x©l¸ €s¸ zs¸ z{…X2ˆ2x1}‹y{ӝ|2z{#Œ[ÛUƒsܐ‡Ý · Œ‡ˆ2xsˆ"zsÑ §µÆš¨™§ÆS‡X±„ÇȉԞ ÆÕŸ¡É9çè¢ £/Éé¤-Öª¢ £gš×§J© ‚y2… xs€{‰˜x |2ˆa° ƒsz„}a‚‡;xfˆˆa‡;‰˜z{€{xs‡;yg·Øz{2¢ £@x¤xs|2¥-ƒnŠ‹‚¢ £}P^¦Ù…$ˆ9z{‡ˆa‡;¢ £l€{·ê¤§‡‡Xˆ0‰l¢ ƒ†£ÔŠ^…Xë-y{yÖڛ#… §y{|a‚¬ƒ†xs#ˆ#z{ˆ2[‡…$ˆ2z }a}"xs‚ƒf¬[€{šDx$…Xy„Ÿ'‚ˆ   ‰˜· |2xsˆay{ƒfz„‚‡…$ˆˆa}E‡X‰íî“UÞàƒ†ß;z{Š^á2…Xâcx@y{ãPylƒfä;ڇá âLˆ2å y{z†‰¶æ €{…X‡ˆ"·¡z z„¢‚£D· æ x@z{‡ y{z{¢ xf£ Œ§™˜…Xz{ˆ22}-x ‘ xs³ ˆ2… xs€†« …­‚Š‹¬‚y†xsæ …$ˆæ z„‚[‡"‡;ˆ´|az{ˆ2‡}2xs¬;}…X€„‰¶^…X|a"ˆ2Š‚ƒsxz„‚z„‡‚ˆ· x@©ª‚y{z{z{xsŒ2ƒsy¨‡;‚· y Œy{z{…X€†ƒs…­z  y{"|2z„Œa‰˜Œ[‡€{‡;©ï€{zs…X¸ €{« }ìSžnxs¸ zÙ쯍xfz¨x|2…gy z{}P€„^‚…Xy{ƒsˆ €{« xf|"z„Š^‚…Xy{« x€ xs©ªŽ"xn| …X"z„‚z„‡;‚ˆ ˆ´™{Í$‰˜‡;ž€l|az{y„a‚ˆ x « y{z{Œ2…XxGz„^xf…sº"Š ŒP}"Š‹‚‚ƒny{‚ƒszԀ{xs»¼z„‚|"y†Š‚…Xxsz„€ ‚‡;y{ˆœƒ†a¸xs« ·½ð 2xGx ‰˜‡;}P€¨‚y{z{ƒf2€{x@xsz{z„x ‚· |2xµˆa}PÀ"‚ˆ2y{‡$ƒf© €{xsˆaz„‚y½y†…X…Xz„€{‚x ‡;ˆz{2…$xˆa}¬;…­…GŠ‚|2¿xsˆPy ‚z{‘xlògñ ¬;« §J‡Xó׊‚|2· Éxl¢ ôîy{ƒn§c2õöxs· É9xµ©ª§ ‚z{‚1¬;xs|2ˆ#Œ2y{"z{~€{x$z{…Xa· x ˆ"‡X‰ï|2· |´« xs‚€„ˆ‚ƒ$z{…s« 2Š@x y{ƒ†z{2€„^xs…X· ˆ xŠ‚x-™÷õ7ž ™˜}aÉøxfy{ƒs…$€„zg‚[z„xs‚· }éx½[ÇxnŠ‚‡$©µŸ¡žîó¯Öí[z{¸2x1¢„ˆ-¬;…s‡Š‚€{|a}2xExs€‘ z{ñò!‡‚}ay xfxsy{ƒsº"€„Œ[‚[xfƒsx z{xsz{}2x z{‡y{ƒ†[axxf· …$xˆ×§J…$y{‡;Œa· Œ2€{x‡$ºˆa‚‡;·z†…X…Xz„z„‚‚‡;‡ˆ9ˆ2y‡X…$‰ €{z{x12€{x«xfŽ"· |Px$‚€{…$xsˆ}S¸¬;…­w´Š‚|2xx }2z{€„xs^ˆa…X‡;ˆ z{xgŠ‚xî"õ?~EƯ§"° ùg~ ŸDü Ï « ™™Æõ-“ ž{žnù §–Jz{2ŸýxgÍ;y{±fƒ­÷P…­±†Š^…Xþg€lz{Œ22€{x‡"ñ y„}2‚|a}aƒsxszµyl‡X‡$‰®‰®Æ ÿg§SÉ…$¢ ˆa£ } ¤ ©ïò ‚z{Ý Ï-z{ù axlɈ2¢ £‡€{¤ ·î¸ï…­ì¯Šxsz{zµ‡õúy„‚}aÉx½§ü Ï©ª™xlõ-}2žnxs§aˆa‡;‡;|az{zxg© "…X€{~E}û¨z{™‡îõ´õ?žª¸ïz{2ì¯x xsz …$€{¡x$…g}2‡$xs‰ˆ2z{‡az{xx z{š 2xg° y{«ñ xsz  ‡$‰ "xs§S} ©G« xs2y@xs€{‡$x ‰š z{ax° · ñ xs}2y{xsˆ2Ö2‡‰˜z{‡xs€yl1z{2É x  · x­§ …Xà ™ˆ-"žª¬…s‚y@Š‚|az{x½2x@‡X‰®Š‚xsš¨ˆ ™{« °‚z{±{´Ç žn‡$°‰z{ñ ax‡ˆ-xs} z{« a x xaxs §} x ñ a…î§L|2õ ˆ" ‚ñ z@Ý ¦ ˆ2‡™˜€{€{·xsy{…sŒSЏ¬õ xsƒs ñ z{݇; € ž z{Ï z{ ‡ a x y{|2|2Œ2ƒny{z{€{z{x$…X…X· z ™ ˜€{ñ xsŸ²y{ŒSü ¸Ï ™}aõ ‡$© ñ Ý ˆa y{ž{z{žn€{§ x­…X…X·ˆ2} žªñz{€„^…Xˆ å« Š‚Ïx½Ÿ z{‡  ñ …X™˜z €{xsz„y{‚· ŒS¸ xlÇ ñ §"˜¸ x¸Ï y{Ÿ |aƒ†z{ ñ žn…X¸¨zªñ ð z{22xsxg€{xlˆ"|2xsº· ‚y{xsz{€„yԁ‚ƒ­–¼…­ŠÉ « y{ƒ†a͐±fxs÷"·½±†xgþ l€{xsy{Ž"|2|"ƒ†‚ €{xsz{y…X…$zˆg…XŒ2Ÿ9Œ2€{ü ‡Xº™õ-‚·î C ñ …XÝ ¦z{xž@¬;™˜€{…­xsŠ‚|2y{Œ¯x§¸ ò  Ý ò  Ý $‘!ñ  §"‡X‰¼‘1‡ˆE…Xˆ"~ xf} « "x   …Xzªz„‚·½x@Ç ñ Ö2y„‚ˆ2ƒsx$› #   P$§ $‘!ñ  ‚y z†…XÀxfˆ1z{‡[xlz{2x@|2Œay{z{€{x­…X·¡¬;…­Š‚|2x§a˜¸ x¸ $‘%ñ  Ÿ9‘[&òñ ' å ¸  w-xl·…s~1z{2xsˆ}2xn¿ˆaÑ xz{axl|2Œ2y{z{€{x$…X·¡¿[ˆ"‚z{x¬‡XŠ‚|2· xy{ƒn2xs· x"~Ä Ï Ï () ‘ òñ ¦ 0‘ ñò û¨™õ-žœÎ ‚ Ï ¯ Ð Ñ Ý ò ” , à   ™ ü  ™ ´ õ { ž ž š ° ± 1 õʁ‚ˆÌ3± 1 ó ‚ˆ¢ ôg± ñ ›œ.™ $‘ ñ /.0 Û ò Þ žJŸÙ "2 í * ™÷ž )+ ³ Í 7Ñ 698;:=@?BAC:=DF>9E G;‘ >ò9HIŸ DCJLKMû¨DC:=™>.õ-8;8N5ž >'4 OQPRò >.E ‘ SN³ T7™?BÆ[DU>.ž{V"•"Ư>@8W± D XRTDCYZ>.E V"TDC:=[\N>]A^O]_F`baF`bc&HIDLJLKdDL:=>@8;8N>fed>^G;>^1 ghOFiIõÊ ?.Tb8N‚ˆÌ j]> ± 698;:=@?BAC:=DF>9E G;>HITb'OP >.E SNTb?.D >@V">.8;D X T'DLYf>@E V"TDL:=['\N>^A]O^l9monnba_;OFp7caFnne!YNT7Vqf>@E ?Bred>]G;>^ghOFiI?.Tb8Nj^> ¤ s

Ï Ï Ï w-z„‚· x x€{xsÇ ƒ$…sŠ‹[Š xz{[}2…$xnz ¿[ˆ2š xsñ }#° "ò ~#Ý Ä }a‘xfˆañ ‡;™Æ[z{xsžlygŸz{2‘[x ò!ñ · ‹x­‰Ô…Xƈ¬;É酭Š‚õ|2x§œy{‡X‡‰ šz{™{…X° ±ÈzsÇ §œñ ‰¶žn‡° €ò … Ý ‚‡¬;ˆ-xsˆy„‚}2· xxfy{ü  ™õ-…$žnˆa¸ }?ìS…1xszlz„‚z{· 2xx y{…$z{ŒaxsŒ2Œ#€{‡$í[º§J‚·z{2…Xxz{x…XŒ2y{Œa‡XŠ‚€{|2‡$z„º‚‚‡·ˆ’…$z{…Xx z « y{‡XŠ‚|2z„‚‡ñ ˆ‚y@}2xn¿ˆaxf}xs¬xs€{~"© 2xs€{x"~Ä ‘ Ý ™Æ¯±{ÇȞ¨Ÿ9‘ ñ ™Æ[žn±9‰˜‡;€9ǪÉ  óSí ±s™ó½”’Í­ž{í ‚± Æ1É¢ £ ¤ ° ™þž ð xsˆ"az{x½€{‡;…­Œ"‚~´· ©ª‡Xx­‰…Xz{À 2xy{‡$ŒaŠ‚|a€{xsz„‚y{‡;xfˆéˆ"zl‘´Œ[z{…$‡´Œ[xs™{€GÍ­žn‚yµ§œ© z{‡2xsy{aˆ´‡$z{©D2x z{2z„xg‚· ƒfx½‡ˆ"y{¬z{xfxs€Œ« xs…$ˆ2ˆaƒs}-x z{™U2‚ˆ-x … y{ƒsŒ‡;…Xˆ"ƒs¬xxsy{ˆ"z{‚xsxsŒˆ"zlz{‡‡ Œ‡Xz{‡Š‚‡ «s~2xs€{ž/‡a‡$¸‰®ðJ‘ ‡Ý z{Pz{‡$‚yl© Œ2…$€{|a}a€{yµŒ‡z{y{axx§J™˜©®|ax ˆ"‚y{Ž"[|2…sx$Š‹ž Š ˆ2xnŠ‚xsxsxs·½}xsy{ˆ"‡z · · xµxsˆaz{‡;2ˆa‡"}2}2x y†« ž1xsˆa…$ˆaxf€{}éxsy{… ƒsxsƒnˆ2Š^ƒs…$xy{y„‚…$ƒ$y{…­y{ŠÔ|a· y{z†Œ2…Xz„"‚‹‡Š‹‚ˆ2z„y@~’‡;…Xˆy{y{z{|2a· x Œa· z„xs‚‡;y{ˆ9™˜‡¬;ˆxs€{z{~ 2·x |az„‚ƒ†·½1x Š‹‚À;y{z{xlxsz{ŒSa ¸ ‡;«y{ x@‡|2€{y{x xs}Œ2€{©ïxs‚ƒnz{‚"y{‚x†ˆ Š‚~z{§ª2Š‚xGxs z ‰¶€† ¨…X· ± 9xs©®[‡x €{À z©ª‡$‡ ‰¯¿¿ˆ"º"‚xsz{} x Œ[‡;y„‚z„‚¬x €{x$…sŠJˆ"|2·gxs€{ys§[©ªxlƒf‡ˆ2y„‚}axs€l·½xsy{axfyÙy{|2ƒnz{…Xz z{2xs€{xlxsº‚y{z{y EÉ¢ £ ¦ §œ™˜z{2x y{Œ…Xƒsxy{z{xsŒ "žªy†…Xz„‚y„‰¶~‚ˆ « Ä    ¤  û¨™õ´ž   ¤ ± 1 õ É0± ™ "ž   à™ "ž  ¯± 1 îÉ ° w-x…Xy{y{|2· xgz{2x@‰˜‡$ЋЂ‡X©ª‚ˆ « y{z†…$P‹Š‹‚z~´ƒs‡;ˆa}"‚z„‚‡ˆ-‡ˆz{axz„‚· xly{z{xsŒ-Ä Ý í  F÷ à ™ "ž û¨™æ õ´š ñ ñ ° ñ ž æ  ™{Í  "žn3± 1 

É 3± 1 óÉ ¢ ôg± ™ ž © , ë a™˜xf¢ €{£ x ¤ žª‡;é-€  Ô, ‚ë yµ™˜¢ … £ « ¤ ‚¬¥xsˆ¢ £lˆ"¦¨|ažn·§"‰˜[‡;xs€l€s§ Í  !./Ÿ 1y{0 |aŒ¸®ôG™› ‡;z{# ™xg$žnz{±/…Xzs!É §©  a"xsˆ#± "%z{$¶ažnx§2…Xˆ2ˆ2‡} ˆ2}a"’x 'Ÿ xsˆ2&sxs‘a€{xs(³ y{&*ƒs)xsˆ2¸¨ƒsw-xx@…$y{}2y{|axsˆ2·½‡Œaz{xlz„‚‡;"ˆa!~ y&X™ +° "& ž ë …Xz{2€{xlxly†ˆ2…$‡z„€{‚y„· ¿[xf‡;}¯ˆ § « z{2x « Š‚‡…­ŠJŠ‹‚ˆ2x$…$€lƒs‡ˆ2}P‚z„‚‡;ˆ´xsz„©ªxfxs!ˆ …Xˆ2}íĪí  2÷ ç   ™{Í  "ž¨‚· ŒPŠ‹‚xfylƒs‡ˆ2}P‚z„‚‡;ˆ™ ;žn¸ w-354xl68z{a7:xs9ˆ6€{‡$= ¬xlz{2x@‰˜‡XЋЂ‡$©ª‚ˆ « z{2xs‡€{xs·Ä Ñ ?

Suggest Documents