A General Object Oriented Framework for DiscretiCing Nonlinear

 !" #$
% &'()*,+-./ 012
..35476
.89* 4;:<8=3*&
>@?BA-CEDGF5HJIKAA-C
LMONQPQRTSVUXW-YZ\[B]_^a`bU\c'd]fegPQRb`b]hSj i SfkBlmRnPQ]oNTY(Rp^qk9r<PsRbt j dvu
wsxzy|{a}E~€‚ƒ…„(„‡†ˆŠ‰(‹ŒmˆŽ*qŒmˆŽ†q‘’Xƒ‚“m†”\•(„vv†–(ƒv„—’,‰…
>˜F?BA™D›šœž™Ÿ?BF-™A
LMONQPQRTSVUXW-YZ\[B]_^a`bU\c'd]fegPQRb`b]hSj i SfkBlmRnPQ]oNTY(Rp^qk9r<PsRbt j dvu
wsxzy|{a}E~€‰‰Š“*„‹ŒmˆŽ(‡qŒmˆŽm†\‘’Xƒ(“m†”q•„v…†q(ƒ„‚—’z‰v
œ=™mFKFKCfšœ Cf™m¡K¢
LMONQPQRTSVUXW-YZ\[B]_^a`bU\c'd]fegPQRb`b]hSj i SfkBlmRnPQ]oNTY(Rp^qk9r<PsRbt j dvu
wsxzy|{a}E~€‰“(“m†Š£‹ŒmˆŽ(‡qŒmˆŽm†\‘’Xƒ(“m†”q•„v…†q(ƒ„‚—’z‰v
¤<¥G¦ ™AЧ©¨ª¢ ¥¨ «*¬¥ AF
LMONQPQRTSVUXW-YZ\[B]_^a`bU\c'd]fegPQRb`b]hSj i SfkBlmRnPQ]oNTY(Rp^qk9r<PsRbt j dvu
wsxzy|{a}E~€„­*„ŠŽ‘…‹ŒmˆŽ*qŒmˆŽ†q‘’Xƒ‚“m†”\•(„vv†–(ƒv„—’,‰…
® DGA-C ¥¯ ¡K¢ ¦ ™mA)°›™A
LMONQPQRTSVUXW-YZ\[B]_^a`bU\c'd]fegPQRb`b]hSj i SfkBlmRnPQ]oNTY(Rp^qk9r<PsRbt j dvu
wsxzy|{a}E~€Œmˆ„‡†­Š‹ŒmˆŽ(‡qŒmˆŽm†\‘’Xƒ(“Ÿ†(”a•„…v†q‚ƒ…„‚—±’,‰v
1
2
²
$
.
³
´
·
µ
¸
¶
¹
µ
f
´
¼
º
›
»
p
½
¸
¶
.
µ
V
¾
¾
¹
³
¿
$
À
.
³
À
Â
¶
Â
Á
$
À
z
»
¹
µ
´
p
»
¼
Ã
.
³
Â
¶
T
Ä
f
Å
Æ
Á
.
³
È
À
T
Ç
V
´
.
³
T
É
Æ
¶
p
»
B
Ê
¾
…
Ë
›
»
T
Ì
p
»
f
´
Â
Á
¼
Ã
»
¹
µ
È
À
¹
µ
$
É
f
¾
f
Å
h
´
.
µ
p
½
Í
Å
Î
¿
.
³
f
´
Ê
T
Ä
¸
¶
¹
µ
f
Å
Æ
Á
.
³
À
Åf»XÏÐ×nµ¹ÄÐÅ'µ¹ÅfÁÆÁƾ³.À$Éоzµ.؟¾V»XÙbÌ9Åhµ¹³.´fÀÅfÁÂÀTË*º2´fÁ¿ÎÅf´VÁÂÀT³.ÊÚº±ÅfÅfÏ$ÏT»BÁƾ³.´f´VÁ»oºÑ¿ÎÁ³.Àдfµ¹Ê¶(ÄTÊB¶¸µ¹³QÅfÌTÁƳ.»pÀa¶…Ûmµ.¾Íµµ·¾V»p¾VÅ·ÒQ¾f³¹ÅV¿»pÊÓÁÂÀnÅV³¹»p¿v´Ü¿OÔ$µ.´V½X¾f»JÅ'½p³.¶¸´Vµ.ÌT¾V»p¾V´»X¾·Çе¹µ¹´f´VÅf»<Á¸µ¹ÌT¶»pÌs´fÁÖÁÂÕqü»X»p̘´V»pÀbÅfÅfÏÐÁ¸µ¹µ¹Å ¶
¾V»Xµ¹½V×n¶Â¶ÆÏ$Äгz»p˞µ¹ÊBÅfÁƵ»X³.¾zÀ$ÇTØn¾´Vß_³.ÀKË*ÉT¶ÆÇÐÁ»pÅfµ¹ÏÊÝ´fÅfÅfÁÂÁÆÏ$À$½pÄT»GÌT¶¸»pµ¹ãaÇq´z³Q»pÛ¼½zÀ$Åfµ¹ÌTÏ$¶v»p»‡Àn䵹ŊÉ$ÁƾVÁ¾VÊ3½XÅf³.´hÇTÀbµ.¶ÆÅf½p»pÁÂŊÊBÀn¿àÄ$»p´h³.Àbµ¹Ä$ÊBÅhµ¹¾›»pÅfËvåÁƳ.³.µ¹ÀK´f¶Æ»pá-³¹´f¿$áQÁÆÃѾ(Áµ¹ÀÈü´f»pÁƵ¹´f³.À$ÒÄ$¾h¾aËvµ¹ÅVÅf»p»pÞç¶ÂÊ3¶$æÂ¾fÇqèpÄTé_³.ÁÂØÅV´h»Xêµ¹ÌK¶TÏ$µ¹¿â»›À$³.´ŠÌÍÁÂÊ3Ìs¾fÇTÁÆÇоh¶Æµ.½p»p½p´VÊB»pÁ¸µ¹Åf»p¶QÁÂÀnÞXÌsÅhÁÂÀTÁƵ¹¾VºÍÅf½pÁÆ´V³.»p»pÀ=üÅf³.ÁÂÞ,³¹¶ÂÄTµ¹¿*ÅfÅfÅfÁÆÁƳ.³.Ï$ÀÀ »
ÇTî·´Vï±³.ǖð·³Ñ¾Vñò»XÌæôós¿ÎÛ$´hµ¹õ.ÊBé_Ø »p˅³.´fáÁƾ-ÌT³.À$»GË*ÁÂÅfÏTÁÂÀÈÅfÏ$»äÁƾfÅf´fÁÂÉTÄTÅV»X̞µ¹À$ÌëmÀTÁÖÔлXÌÈìŸÄTÊB»p´fÁƽX¾Kí‚ÀbÃbÁ´V³.ÀTÊB»pÀbÅ
ÿ
ö÷±ø‚ù…úû9üýŸø(þúG÷
"!#$!&%(')*+-,.***/%0,1233%'('(43%135.67'(98
*:;<35235>=?'"3)@12ACBDFE.CGD10)'7HI4JK%20> MLONOPRQTSVUW%!X5>0>ZY)%
%*'[%!V'1\5R4)!#W@.'] 10'X4JK%20> <'^352W',.363%%!V.\T'"3)11'D
*%)3520>'S4_.%+'"\4352?,)%@*0`0,*04\%$a@bD2%c%>5?@.''%!d35D43)!#
e)fhgZiZijlknmnopcqhrcmnsZoItugZv6pwox#yXzvZz{x&m|o"k|z}gZy~|Z
€g kVtuzpwg6vZ‚XgZv6pWƒZjk|x|„…shgZvZ‚^gZv6pwo)k„†jlvm|k|‡„mˆl‰ŠI‹ŒŽuƒ‹I
f gZiZijlknmnopcqhrcmnsZoI‘‡lv6pwoxnx#mnz}~’mngZv6‚<tu‡pwo"vZ“”CZ
€g knm#m|o"y[qo)kn‚^gZv6pwo)kV„†jvmnk…‡l„m•w–"“˜——™Z •l‰š\›6
 h
œ
!"$%'# &)(% +*,.-/!1032
,b24<%!352W.3T')3)11W5.
76
'^1')!&%
8
54
7',bFE.
,,2*3% S R%%wJ.')1\5]T%8
9
,2W**%Z;<'<1'7%! 2')%'D@d %2R352M4)!#S
%![a 10T@b?%!,%6Y)'W;<35
4
6350
'?
<8
:8
1'"3)˜*/
+;
_ )@12S D*352%15
·ÿ
52'cR "!#
,5 'h* ,,*h63%'?6?352
;<Ta%!,,*h63%'D5'?3%
9
)%1>!&%I35D1* M0>3520>3h*0%*'[
@bD!&%12 S
'c
43*V,."
;?,,*
')31K3M!&%
@dT%wJ) >%00>%
5'[%@'))J%>'1'"3'^%M98
E.2 *')%(?%0>0%743)!#^!&%/5X 10)h*K')\50'+2M3%(,)%wJKXW[0,*0>4\%
!&%+352D43*d,)'[%!u35c'\50>'S2GD%;<F35'"\2
(')'363%12'2
ÿ
,.,d<;[
'))@b
>=
5
5
6
@?
A
%K*.=?'% 1%12'
8
*
'363:
>?
A
CB
lJ%*123%
ÿ FE
1.3%'
Kl3%
1')5
6
F3' 3'%
')),2%
LON]P]Q
ÿ
%
!&%0
35
IA
5
";?,2,*G4
D?
')),2%O%!V35
%*12
5$0
A
*
'1*F3'
,2%0')S '\520S
CB
3520>3'
6350
D9
76
'0>,2*
)8
20>,2*S
_.@1)-;<5\5
*0>%*'(
HE
)
5?!&%0>*
?'ŽlJ%
6
':3%3527LONOPRQ
JE%LKS
35'
1.3%'S V%C')0,*F!
5h')R%!W
T'3:^' ,,*h@2*S K%
0>,*04\63%.* 3*'X%!V35
;<F35O'%0>
')31K3]%!
)$,'4 !&%1')%
63*0%K2*'
a'3:D!&%X'):7E')<%
'$20%')3
)!&%)0T1*˜63%7')%W35.Ž35X*%K*=?')%4 1%1'
%
8
%A0%]%0>,*3lJ%*123%
352]2'l3:3%
5],%6Y)'
*
nÿ
%! 35c0>350
œ .'0'J
')3 !&%<$* 0,*04\%>%!V35
35? <'%a;[(%12*2D352(=?')%4 1%1'
5]@')3l
18
35DJ)'
%43 S
MNPORQTSVU þ W(ýúG÷·ø(þ¹÷üJúüXWZY\[] S _ù ^þ÷ U þ W(ýmù S ø(þ`a[Bø(þ,úG÷
bdc3bDcfehgjiklnmfiopqisrTtduDvjwxlgjudyzrX{/w|}lgjudyk
IAD@b%1+')1')'%;<F35> '"31~>
8 %!uE.)')[%+'8')0'^%!d35<!&%0L ∂ U(t, ·) = L[U(t, ·)](·)
t
;<5 35 '),.3* %,b3%X'<lE.]@:8
€ œ$
L[V] = S(V) − ∇ · F(V) − A(V)∇V
;<5
'+ '%0><'12\@*X!&123% ',.
S‚A
I <!&%12'52c% 35c',c'))3:h68
V : Rd → R p
V
% |S SF;^D%'"3)1^
ƒ8
Vh
B
1.%
€ ' ;<5
œ
,.)3F3%$43%
WE.F3
%R35T**'
@12?0
:8
Lh
%43%]%5DE.l8…0')%* !&1%
,* @
2'l3c%,d3%
,.)3'(5
ϕA(V)∇V
∆j
„8
3'"^!&123%'
10T@b+%!d**'
'c$0>'12
Wh
∆j
%wJ^35c%0>
6
0M1*F8
;<5\5 '
Ω ⊂ Rd
)c')0>%%5T!&1%2'
j
%'M5(Y)120>,A%!
0h')1
5,.5
}9
S
…K 9
F(V) · ∇ϕdx −
∆j
<8
%%2'"J3FJ,%1l
5 1'<;[7FJ
XZ
j
ϕA(V)∇Vdx +
∆j
ϕ^
F(V) · ndσ +
∂∆j
XZ
j
-l3%A%!X35>%)0
V
(ϕA(V)∇V)Φ
Φ
;[>43),l(35
∆j
j
j
,d'$%
:8
52'(%,b3%D'W%'"3123A@
lE.R- |S
(ϕA(V)∇V)Φ
Z
XZ
L[V]ϕdx =
S(V)ϕdx +
'W35TJ*1$%!+35
q9
S
ϕ
S‚IA<'')10<5.
(∆j )j=1,...,N
V, ϕ
'
XZ
^
ϕA(V)
:8
E.l8…0')%*K!&123%',.
2a@ >433
%ZJ(352**'
XZ
[V]n
'),.
,2'%
]@bM'% 1%12'W%wJD357**/43"!#'SuGD d;[M!&%0>** a433
J0
Ω
;<5
;<5\5>0
∂∆j
*+6M**+ "!#'S
18
,d%4(%![')%4 1
5a3*%3'
hJ
€ 
^
ϕA(V)[V]ndσ
@dl;^
352
^
ϕA(V)[V]n
wJJ*12
35**^4)!#'ŽnS}S
T
% ( & % x
C
2 % % % ~6
10)h* .1 K'
5.
9
R
!&%W5
CB
W h ϕh
5Ω
%'"3123%
~6
%!
,d S _.%(352
= 6
‡ÿ
!&123%
5'))3$!&123%
Lh
,d'a%35
~6
.1 AhJ3
6
Wh := Lh [Vh ]
€ wJ !&%<**b3J')*<1!&1'1l')3% ' ϕh ∈ IW;<h5S\5
€
35(J*1'
ϕh A(V)
!&%0
,2,% K0>3%
43)!#(%!
5'))3:h%
6
35
'
;<535
J%*Ja!&%( 20,*T55(%)
'8
,)%K21V;[a0,*% A
'),.˜*W%,d3%
FJ63J'(%!
=
2%0>,b%'F3%A%!
L
A
~6
.1 C* %'')3'>%!?5
5
.1
hJ3
S
5.'O
%(%K8…*2h)3'c
V
,2%@*0
<0>'%*J[@'
∂∆j
bdc"!ac$#Fp ypqi!|dv&%('/kl~p*) kfu+ rTtdudv3wxl~g3udyZrX{/w|}lgjudyk
_.%
1')aO
%!W5R0h'1)
@d%5R')'<%!/5
56
%hA5
^
F(V)
„6
'0>,2*'"h')R352
%!
'(E.2C')%
!&1l3%>%+!&%+ 20,*c')%0c0%)c')%,5')h63
6
nS 352
@b;[
,.)S
54<5.]'(%!
%')J3FJ
8†_.)2\52' 1
q9
%')"J63J
1.*'X5
76
0%
%0,*
!&%)0R
352$%K8…%')"J63J
S lT1'T'')10>35.6M;[a5.hJa'),.˜*%,b363%)'
Ls
;<5
'1\5 35
s = 1, . . . , S
Vs−1
V0 = V
WS = V
12'2]5>,2%6Y)3%C%,d3%
S …!Xh\5A%!^5
L = Π S ◦ LS ◦ · · · ◦ L 1
ΠS (VS , . . . , V1 , V ) := VS
%,b3%'
'+?E.)')%2I%,b3%Ž%! 35 !&%)0
'1'')a@d%wJ 5$;[c %2')3)1'
Ls
%,b3%'
'<@b!&%) ]%0T@2W5')7%,b3%'X3% lE.
S
Lh,i
Lh = ΠS ◦ Lh,S ◦ · · · ◦ Lh,1
!&%
0>,,2
=
[1'0%'"3363
L
Vs := Ws × Vs−1
%
35',,)%4\5
q
5ch')D%! ‚ÿ
')0,*
-,
')h*+2Jl3%
,
CB
u1')%
œ 13%
€ 
L[u] = −∂x f (u) − K(u)∂x (k(u)) − A(u)∂x a(u) .
. F %43' O'%A%2 )J3FJT0 € œ'A
352 JJ>30
8D35'M*' 3% ,2%@*0'M
lE.35 0h')1$% 5 @d%12.38 …K
|SqIAT352l!&%$!&%)0 **<8O);<)3T5'(CB 1.3% 5
!&%**%w;<2M!&%0
1'1~6K*˜38>!&12%'

w = (w , w )
9
5'
B
1.%- d)'
;^%
',bl3'
!&%0
5a%8%')J3FJ>,%1lT'7%K8…*h
1
1,1
1,2
w1,1 = ∂x k(u) ,
w1,2 = a(u) ,
w2 = −∂x (f (u) − K(u)w1,1 ) − A(u)∂x w1,2 .
=?lE.M352?%,b3%)'
%!
ÿ
L
%2''"3
L1 [u] = (w1 , u)
%!Ž%,d3%'<%!V35
/
6
a')%0>Dh')'X*%')1D*˜%'+FJ
_6
%12 S H 0>,*'<)W;^%68,5.'
5.hJ
.%w;
2%.*/30'(%!^5$!&%0
9
%,b3%S
L2 [w1 , u] = (w2 , w1 , u)
'O!&%)0RS
.
(U))
0%K*'
ƒ8
Q(A−1
52'(0'XO% @bWl3*
2,b4M'
3*
)$!&30;[%
qÿ
A
**,23X%,b3%)'^5.hJ?3%7@dc\ 4%
:8
8
O33%a5 % .0'S O352'
5
5.*
%')1l3%A%!^,)%@*0
ƒ8
:8
!&% 20,*c@
;^W)JDX7%0,d%'%
Q
'
Š‚ÿ ÿ
38
5W!&%)0
%,d3%'
18
O%8*h R2
'@bO')%
A
''c;[
'%0$*
B 1?5
8 h-@b
!#^@12<)
−1 S …!^5')>hJ*˜@*
35
A ')@bAA@b%ZhJS 35'W;X„8 5 ,,)%4\5-h-@b
,2,*O3% %'"312W2'l3M%,d3%'<!&%xB 17%0>,2*3R' ')30>'D%!IlJ%*123%RCB 1.3% '
;[ ;<**d0%'"3363 LEK% KS
*1
35
bdc10>c32 g4) p65ng3k87+ipl~g49+|}lgjudy
EK%!#<;^75.hJ$
;<5\5
9V%
*h2'a3%
38
CB
'0> 80,*F^2'l3:3%]%! 3527lJ%*123%R
8
' ')0
'%*J>35'7' K'"30
:
∂t U(t, ·) = L[U(t, ·)](·)
!&% 35%K > 'OE.2
S
d
Uh (t) = Lh [Uh (t)]
Uh (t)
dt
')13@*>')%*JV|S SF/ c12l8 (123R0352%KAh @ba,2,* S
%!TU(=WHI'
%!DU(=WHI'$
1.3%
;
=<
K
!"$%'# &)(% +*,.-/!1032
18
=?,d%
5T')3@* O)'"3%'c0,b%'
2,2*
0,*F(0352%K S V%
76
%M
352
'30>
30>
76
5'
3%
:8
@ R5T',3*/%,d3%D%Th
%wJ%0>30l8…'),C)'"3%'
C8
l\2O5>0
6
13 A%! K,*
0352%K'
!&%
l;<(35
B
J%*1K3%
1'$5?
6
!&%
K,*F7')\50'T;<52*
%K8…')3)l3J)0>'7
9
%!? 2,2*0,*F('%*FJ''T')%0>l30>'T35@b'"
')13@*%0T@.%
5' %2
.9
,,%\5 S V%C\52J
1.%'<1')$;[%>%,d3%'
€K 
d
U(t, ·) = Lexpl [U(t, ·)](·) + Limpl [U(t, ·)](·)
dt
;<5
Limpl [U(t, ·)](·)
%0T@'W**V35
18
q9
')3@* ])')3a30'S
<
52%)'),d%22'l3
%,b3%'[c%2')3)15?'0><0>^'[%1K3*'%7!#>!&%5c0l8…')l3:h%
6
T')0F80>,2* <198 (12)\M05%K'X1')
!&%
S
K,*F
,2,%4\5O!&%
OO0,*F
Lexpl,h
;Limpl,h
?
];^$%*17%1?2'1'')%R;<35 '0>,* 620,*S^%2'2D5TlJ%*123% CB13%
;<5(352(',.˜*d%,b363%X' JO@:8
 S!9 5(E.)')<%)X)0>' *ha%>M0l8…')3,')3)l3%
€
5%(%!^5 ':
;<5)h' 352>-,u12'% 3)0
43)%K21'7'"\@*F18 )')3%%!
h
5 %
S
9
5
)
&
!
%
)
M
$
'
2
,
*
)
3
T
%
V
!
3
5
(
&
!
%
)
0
h2
Lexpl [u] = −∂x f (u) − A(u)∂x a(u) ,
€ 

€ Limpl [u] = ∂x (K(u)∂x k(u)) .
' ,,)%,)˜SCA
I I1')[X')0*d2%0>,b%'F3%(' @b!&%E.2 L [u] = (w , u), L [w , u] =
expl,1
1
expl,2
1
2
@:8
(w , w , u)
L
[u] = (ŵ , u), L
[ŵ , u] = (ŵ , ŵ , u)
2
1
impl,1
1
impl,2
9
w1
w2
ŵ1
ŵ2
=
=
=
=
1
2
1
a(u) ,
−∂x f (u) − A(u)∂x w1 ,
∂x k(u) ,
∂x (K(u)ŵ1 ) .
5'<*h'X3%
9
d
U(t, ·) = Π2 [Lexpl,2 [Lexpl,1 [U(t, ·)]]](·) + Π2 [Limpl,2 [Limpl,1 [U(t, ·)]]](·) .
dt
:8
A
‚I
5c'0>,2*'"I30>98'):h63%T'IFJ>@ $352c!&%";X@ 4;X)aHI1*+0>l35% S
'"3,
IA
∆t
352'
*'X%
F35
0l8
U(tn+1 , ·) − U(tn , ·)
= Π2 [Lexpl,2 [Lexpl,1 [U(tn , ·)]]](·) + Π2 [Limpl,2 [Limpl,1 [U(tn+1 , ·)]]](·) .
∆t
%*12
5.
5
'')43*'),'<!&%X352
œ Sc);<)3(5(',.63˜* %,b363%<'c$'120
Sc2%0>,b%'
KSc12'
35
@b%35O'),.3* %,b3%'
=?'%43 1%1'
K2Sc12'7O6K,*F0>,*
Gc%3W5.<'),'c%
57]!&%1
?
A
*
%!/;[%%,d3%'
4%
<
<198 (123U
=?H
5
Lexpl , Limpl
E')c%)<%,b3%'
O,,)%4\5O3%
;^%>(,2%@*0
h]@b
2'l3:3%'12''
)
S
%!V35(!&%)0
%2')3)1c'
'%*FJ<%2J
5(!&%**%w;<
œ S
%,d3%'S
5
~6
'%*123%0S
,b24D;<52'<35(!&%**%w;<7'))3:h%
<8
0>,2*04)h** >!V')13@*( 10)* .1 K'?)
hJ*˜@*S
')3,'
T
% ( & % x
C
2 % % % !/c3bDc32u
9 iip y l~gj|Dv e vjuauaok
N [3ù…ø(þýþL[Bø(þ÷
5
.%w;
ùvú S ýmø W
%![;^"8…'04
6
0> 31)'
3E ‚A
%0>,2* KF %!Ž5
B
0>%* ')5.**%w;
:
0>,2** $lE.
<8
,5'(OG?hJ†8 % 'X
6 18
!&)
38
@d%1. '
9
5A0
18
1*
8
)$%! 3A0%K2*
1.63%'^!&%[M%OGc;
.%w;
%˜
,2,%4\5
:8
@ O'%*FJ
')3)a!&)%0
.1
;<5
5
'<1')];<5\5R*h'
CB
O 13%]!&%X35
4J%*FJ'O352-2FJ3%
%!M')13@*
%0>,2')'@*M;[%
38 59
!&W@d%1. S /%$)12?5
3%a>]0>%* ;<5)
,235]%!V35
.%w;
h
0>%*'a!&%35 ).*
38
>9/%0%'"3335-*
!&)%uS HIJJ >'0>,*cJ)'%'X2*12cJ a%0>,*30'^4J%*FJKM@b%5
FJ63J'R
%K8…%')"J63JA,2%1'S
<8
,2,%4\5R'@b],)J%12'*
9
B
18
57@.*
<8
%*
ƒ8
%* a,)')4W')%0>(,)<%!V352
1.43F3'D
Sτ
*12W;^%
'
357,K35
]')10
%!Ž35
.%w;
h
0%* 2FJ
,
R5M0%04310
%'')3O%!3)0>'
0%K2*
'%2]%

€

S
hu Sg
'c5M)J2
5>2 d47)*J4
30'ca35(!&%**%w;<2
@b
']!&%]12'25
∂t U + ∇ · F(U) − Sg (U) − Sτ (U, ∇U, ∇2 U) = 0 .
!&%)$3%]wJF IA
;^
5
S
')')''S
U = (h, hu)T ,
1
F = (hu, huuT + gz h2 I)T ,
2
Sg = (0, gz ∇bh) ,
Sτ = (0, −µh∂x2 u1 + sin(φint )gz sign(∂y u1 )h∂y h + . . . , . . . ) .
c
)
!&%
gz
'5W0>F31c%! 3hJ\63%.* *3% %'"\4'D0%K2*75
I
43.* !&l3%a%!V35
B
5Mc'0>,*Ž3'"!&%0>3%(;^hMl;<F3+
U
+B
]')c%!ŽF3%.* 143'
V1

b
1*˜
1.%
8
2'@d'^35?@da3%,b%*% 2
€ 
0
)˜*nS
8
φint , µ
B
'Vc' ')0
%!2E')/%)V
1.%2'
€ 
V1 − ∇ux = 0 ,
∂t U + ∇ · F(U) − Sg (U) − S2 (U, V1 ) = 0
;<5
35
l;
')%1
)0
'
FJ
€ 
'
S2 (U, V1 ) = (0, −µh∂x V1,1 + gz sin(φint )sign(V1,2 )h∂y h + . . . , . . . ) .
Gc%3
5a%K8%')JJa)0>'T
S2
;<5\5)a*h
B
A35a2FJ3FJ'M
':3%O,,)%4\5]'))@bR@d%wJS
e vju gjyXu i!udwk p ogj|
352'('12@'3%
%'W;[%8…,5.')0M1*%0>,b% 0>
ÿ
!/c"!ac32 u |dk5p
5
;^
(
)75T'0T1*3%R%!;^(2
'0M1*˜%
%!/;^c]4' .%w;
DŽ;^a!&%12'T%
8
%* .%w;
!&1* **'
'0>,2*E-' ')0
%Z;
]
,d%)%1'<010
;<5
.%w;`
5T'%*
€ ')(!&% 7 60>,*
+;
,,*h63%R
35274'c,25.'
!&%75
,d%%1'?0>˜2S ?,,2*3%'
€
S Sd@%3.63%
& 
!&%7%0,''@*;^%68,5.'
%0>,b%4W3',b%"?;<5\5-'D5T@.')'?!&%D352
2];^3
35)%15
1
9
0> S
S
.%w;

*1
%D5
0M1*F8
5$0%K2*V%')'(4'
8
„658u
*'%$3527'),d'<5 )% d%
!"$%'# &)(% +*,.-/!1032
;^3 '%0T')
CB
€ 0>*ƒ8O%2'0 '"'3' 2%!0>%0%410 
4'
9
5.hJ$3%@b$0>%** S
13%' )WFJO'^!&%**%w;<'SK_ )%0
2]5
4'
€ i = g
,5.'
;[
5
%wJ
@.*˜c!&%^352?;X63^,5.')
€ i = w
€ œ:
∂t (φsi ) + ∇ · (ui ) = qi (si ),
ui = −λi (si )K∇pi ,
;<F3535(!&%**%w;<2T*%'1)
I
'
,5''
5>,**˜
12
%w;<]J)˜@*'
V2C35,5.'
8
pi
76
!&%.* .%w;`3$%![5
F3%
')
qi
S
K
2,b24
352(,''12
sg
„8
*F;X 'W5.hJ
5
9
35
sg > 0
B
1.%
!&%
ui φ
8
18
,b)0>@*
18
pg
')335]E%'
%!
35
?
;
O2
H2
(c , c , c
9
B
5?
H2 o
35
!&%
qi
),)'43'('%1'W
;X63?R4'<,5.')M'c0%**5%125
'313%
%![5
4'
:8
@
]5
*%@.*u')%1)(30
,5.')
q(s) = qg (s) + qw (1 − s)
€œ 
∇ · u = q(sg ).
1')
B
43%-
*3%
€œ
1.63%
€ S œ 9 5  13'% ;[ ,2*˜ p;<w35@:8 352pgR−**pc,2(13a−%s,bg)3%
'%1',5.'
;[
5.hJ
3%
CB
%'OC32',b%)
R T ;[75wJ
€œK 
1.3% S^_ %
∂t (φsg c) + ∇ · (ug c) − ∇ · (φsg Dg ∇c) = rg (sg , sg c).
a@b?)%,,b
8
@ 1')
ÿ
35D*˜%
c :=
€œ 
,c )
1.%!&%5?'),d'
pc (sw )
,.30:h%
A(pg ) = q(sg ) − ∇ · (λw (1 − sg )K∇pc (1 − sg )).
',b'
~9 5)@:8;^ ''120>235.6X35(,258')'c)W'12\5]5.X;[
pg
5TB 1.63%'W
€ œ‚ 1' 352T*˜% € œœ' ;^MD!&%
u, ui

A(pg ) := −∇ · ((λw (1 − sg ) + λg (sg ))K∇pg )
_.%a\5
235(,2')'1)'c%!
 S
€%'!>
λi(si) 2%3'c5
35 ,5.'''),dJ*<8S
12\5 '36313%
%!Ž570210
35
si
S
u = u g + uw
'Dl30
,5''
18
,b%%' %!X35>0>120
%0>,212\%2'?;^\5%%'$5
%!35?4' ,5.'
S K120>01,
€œ 
;X7'313%uS
0>%* %!
32'%O@b;[
J˜@*'?!&%(%1
qE
'313%2'X%!35
%3'75
35-,5.' 2,b]%
,5.')
'%h***%@.*dJ*%F
5M
S
JR@ 5
i
?5
'MO!&1%-%!X5
%3'D5M@'%*123
'Sd_.%<'"\2
;?'
J*%F3'
,2')'1)/J
5' !&1l3%R'c!&%c 20>,2*
ÿ !A
€ œœ$
sw + sg = 1,
pw = pg − pc (sw ) = pg − pc (1 − sg ).
!A
F35a35'X0>%*d352
5
*3%'
H2 O
R
O2
H2
9 5l!&% 5 J3%%! 35 %43363%'1'3% c := (cH2 , ccO2 , c+H2cO )T+S c 5+'ac 0=%1*
2%3'T35',b)'%C3'%$2'@ 5a0>)%'%, u1')%C%!<35a'),d'S
'M
Dg
rg
!A
')%1 ') a)0
76
<8
=="?
%!/1!A %w;<'D'
R
76
;<5\5O!&% 0>,*(,'4'?%2'
8
%!/5
GD h;[(;<)W35(@d%wJ(2'@dR0%K2*d'<T' ')0
,2,* ]35
0>l35% Su=?%2'%u;[$2-12,;<F35
U := (φsg , φsg c)T
V1 , V2
-,
',b'S
%!VE')<%
5T!&%**%w;<>'
CB
%![
8
B
1.%'^%)^%
9
1.3%'S
5$J%
WJ%'c%!V30>,b%3 a22RJ*1'S
V1,1 + K∇pc (1 − U1 /φ) = 0,
V1,2 − ∇(U2 /U1 ) = 0,
V2 + K∇{A−1 B(V1,1 , U1 /φ) } = 0,
∂t U + ∇ · (F(U, V1 , V2 )) − S(U) = 0.
€œ 
€œ 
€œ 
S
T
9
5
% ( & % x
C
2 % % % E.'"?,''D'
8
lE.
@ 5
9
5.R')S
5
8 B 1.%' € œ . 5 p '<0,*F3*<8
€ œ ;< 5.352 7')% @:]
g
'(352 54
pg = A−1 (B(V1,1 , V1,1 ))
B(x, y) = q(y) − ∇ · (λw (1 − y)x)
*˜'"<,.')'c%')')'<%! CB13%
€ œ  ;<F35352/1~6!&1l3% F
:8 B
!&%0R@ O
1.%'
F(U, V1 , V2 ) :=
!A
]5
'%12 ' a30'
λg (U1 /φ)V2
λg (U1 /φ)V2 U2 /U1 − U1 Dg V1,2
_.% 5D')0T12*˜%T%! J
%'/'/E.2$
B
GDwJ"8 3% 'c
qg (U1 /φ)
rg (U/φ)
| t/g3pqi %@l~udp k rL{@w|}l~g3udyk
!/c10>csp | 7l~gjtdp
œ
.
jE A
B
%Z;<'4'I1@<%0M@1')%'435<J?GDwJ"8 3% '
9
21' S
1.63%'
S
S(U) :=
jE A
h]@b
52X0
''/T0%04310
;<))3]'
€œ 
€ :
∂ρ
+ ∇ · (ρu) = 0,
∂t
ÿ
C5
∂ρu
+ ∇ · (ρuuT + pI − τ ) = 0.
∂t
B
0%043120
5(!&%0
€ ‚ 1.%
τ
%'M5>'"3)''732'%
2
τ = µ ∇u + (∇u) − (∇ · u)I
3
ÿ
;<5
µ
%3'<352
]% 3%
j
'W
J'%' a%!/5
%$35
\50>h*u',b'D
;<5
18
J
T
.1R
3E ‚A
<8
43*
1
)')1*
5
ÿ
Pnspecies
j=1
!&% h\5R',b
'
=?1
3%5
!&)%0
CB
T5' ÿ
λ
CB
,
8 B
9
',b'(2 d1'% S
9
;<5
Dj ∇wj
S
i
-,
Di
1.%O0T12')<@b
-,
5',b'
'(
%')3%%![ u1')%
')\.)R) O
;<5
35
G?hJ†8 % 'D

%4J)'i%'M21$3%]\50*I%'S
−ρDi ∇wi − ρ
I
18
12
;<5\5C!&%MGc;
3%˜
1.3%
€ 'œ 
-,
')%S
CB
1.3%',."3˜*u u4˜* 1.3%!&%^5
€ 
,
5T'%1T3)0
d1')%
u1')%
'
ṁi
%.63'<!&%0
A '
0%*1'R_Ž
%'"\4(;<5\5
-,
h
j
:
-,
5>5)0>% .0$J@*'
CB
%![
'
~6
:8
:8
T
6
∆Hi
jT = −λ∇T
1.3%'W2'1'')-5) F3%.* K,)'')%-*˜
S cV5
h*'7*h;
p = ρRT
'
€ 
.1 T$'/J$@ $_.%12 'Ž*w;'
352W30,b3631)S
j =
@bT\5%',bi8
')1,,*043R;<F35]]F3%.*b3)0R2*hT%
'/352X30,b3631)^ u1')%
8
'),d'
*h;
%!^'),d'(;<35 u435*, %!!&%0>%
5hT6D%123%]% 4D
%2?%*%'M5$' ')0
8
.1'75'
nspecies
X
∂ρe
+ ∇ · (u(ρe + p)) = −∇ · (jT + τ u +
∆Hi ji ).
∂t
i=1
'X5
∂ρwi
+ ∇ · (ρwi u) + ∇ · ji = ṁi , i ∈ {1, nspecies}
∂t
~6
1.Z8
%'))J%7,."Ž%!235[%0,'')@*
1' /;<535
!"$%'# &)(% +*,.-/!1032
h*/4'(%')34
ÿ
@b$35
%!+35
ƒ6
6
0 312M'( K,)'')C'
1FJ'3* 4'<%2')\ S
18
…!?5O',bE.]5h6>,.
'),dE
R = R∗
8
%\* ρe
'
J
cv
Pnspecies
j=1
wj ;<5
Wj
''')10>%A@bO2,b4
R∗ := 8.136J/mol
%!c352O0,d31
'
35
lE.F3%
γ := cp /cv
18
2-5
CB
%!^')33SŽU?J3** ;[a)J7]')7%!X 13%'
€ HK 
}B
!&%(35,."3˜* d)43*
*%') 1')
% lE.a352a,'')1
(ρ, ρu, ρe, ρw)T
F35]$'0,*?3')!&%0>3%a;^(E.** a)FJ c$' ')0
I
,
'(1') V'(;[**I'(35
R = c p − cv
ƒ8
8
p = p(U)
1.43'
352a*˜%
B
%!VE.'"c% 1.%2'
CB
13%
U=
p S
T = Rρ
2
V1 − µ ∇u + (∇u) − (∇ · u)I = 0,
3
ρR
V2 + λ∇
= 0,
Xnspec p(U)
V3 + ρD∇w + ρ
Di ∇wi = 0,
T
€ 
€ 
i=1
;<5
35
∂t U + ∇ · F2 (U, V1 , V2 , V3 ) = S2 (U)
56
l;
.1 OJl3%
O'%12
F2
F2 (U, V1 , V2 , V3 ) =
30
S2
FJ]'
ρu, ρuuT + p(U)I + V1 ,
nspecies
X
(ρe + p(U))u + V2 + V1 u +
∆Hi V3i , ρwuT + V3
i=1
S2 (U) =
Ž5
€ HK 
nspecies
X
1 2
1
p + ρu +
∆Hi ρwi
ρe =
γ−1
2
i=1
;<5
T
0, 0, 0, ṁ(U)
T
T
,
.
gj{/wgjo J| udwiGe vju k g l |Dk5p |Dy dp
IA %' 35>8.0'M%!<]%0,)'')@*>* B 1K8…J,b%1 .%w; 12)%2 ,5.')>3')%'
;<5\5$'/%wJ)@:7
8 5X%0>,2')'@*^G?hJ†83E %‚A'"8=<D%"3;[ '38K'"30 œœ |SH9 5'/'IT63'%
%!u35?G?hJ†83
E %‚A 'X'8 '"30 ;<5aaF3%.*25) %)3)0 35.\HA'[35D ,d^%! ')1)!#
!/c c
')%R<,5.'
!A
;<5(352
E%‚A'
12
@d%1.)'c43%
%w;<'?W352
35
,.30X
:8
38
18
2')
ρ(x, t) > 0
I
CB
JK8…†8 C*'[ 13%
*b0,d31W%
**%w;
18
2O5WJ*%F
18
λ>0
A
' u(x, t) τ
%w;<
'X35
1')1.* GDwJ"8
%')34+T5[!&123%
6
%! ')\63D;<52D35?0,d31c'+E Ka%$7%K8
!&%
qÿ
,5.')(33'F3%'SKGD%W35.6
35
'X2%%
8
35'X0%K*b5. h'[5(,5.''S 35'X0%K*b5(,5.''<Wl30@
5(J*1(%!V352
5<' K'"30
€ 
∂t ρ +
∇ · (ρu)
= 0,
T
∂t (ρu) + ∇ · (ρuu ) + ∇p(ρ) = ∇ · τ + λρ∇∆ρ
'Ja@
'"\4c@d*%w;
%14S
3'% 'ŽM35^'l3%T@b!&%)5X,**˜)
p = p(ρ)
9
18
)!&)'c35
'))JJW!&%0
ƒ8 € *%w;
2') O%*
*')%

@dX;<)
3%
18
18 <B
2') ?J,d%12.552' W*
%')JJ
'):(5(35]%
!&%)0
3)0
1
@1K[F 312>%1KI3%
3%5
;<535
B
a%)X%$c123(')1*'X;<5O% FJK*u*%K*b
S
'M œ Z
b
€
30
O%%K8
@d<@b)3
∇p(ρ)
1*@)10
%2E.13%2'
T
)
% ( & % x
C
2 % % % ,'4S
;<5
κ
:8
' 2E.
€ 
∂t ρ +
∇ · (ρu)
= 0,
∂t (ρu) + ∇ · (ρuuT ) + ρ∇κ = ∇ · τ
@ a5
*3%
κ = κ(ρ, ∆ρ) := −λ∆ρ + µ(ρ),
µ(ρ) :=
Zρ
€ 
p0 (s)
ds.
s
0
9 5 2J4\ %!W5 %%2'"J3FJ !&%)0 '35.>35]J@* κ ,,bh' 76K,**<8 5
CB13% S
' CB1*%R')%0>%')34T6M]')33CB 1*@)10
')33
nS}SR')33 ;<5
€
κ
∂t ρ = 0
2
S K
E
%]5')l3:h%%![35%%'))JJ>)!&%)0T1*˜63% %!^5>CB 1.3%-*h'

u=0
.631**ƒ>
8 %>$;[** 8@.*˜')\50S
€ :
€ œ$
Vρ − ∇ρ = 0,
Vu − ∇u = 0,
κ − µ(ρ) + λ∇ · Vρ = 0,
€ 4 
∂t ρ +
∇ · (ρu)
= 0,
∂t (ρu) + ∇ · (ρuuT ) + ρ∇κ − ∇ · τ (Vu ) = 0.
ö X] S S ±÷ ø [3ø(þúG÷ úø QTSVU þW(ýmù S ø S ú]¹ü<ø(þ,ú÷
hN
9
5'('3%-h*'?;<5
EK12@'3%aKS œ
')1@'%'
35
76
ƒ8
22Sa K,*˜
04\63% S
0xc3bDc 5ngjk 7i!pl~p
9
*:h%R%!35$0>3520>3h*V0>%*/!&%0
'))@b'+!&1.043*.%2,K3'^%! LONOPRQ
2S
352T;<*
e>wy 7 l~g3udyk |dyo m
50,*04\%R%!^O2'l3
%1
'%
=="?
*J4+!&%I352
'):h63%
pqi!|}l~u i!k
S _ù [3ø(ú›ù
,
B
0>%* %!^],.)3*Ž2 d43˜* %,23'
œ
-LONOPRQTS
2'3:<;<52'
5?0>,*l8
1.%'(@.'A%-5$!&%*F8
*%Z;<$* %2,K?%!V!&1l3%O',.' !&12%'<R%,d3%'235.6?<%O!&1l3%'S
!#"%$&(')*',+&-+/.0.21#&#*'+/&34(5"67&#8.91(&#*',+&#
;
!&1%>',.
R := KnR
2S S
V
%1[%,2
'^7'^%! 0>,,'+!&%0`5?2%0
D := KdD
3%735D3
V := {u : KdD → KnR }
c
KD
%'?5M%0>E.* <8
2R',bJ* S
)
;
9V%!&1"35
)8
35737E.*
d, n
352M0>'%'c%! 352M%0>
',b! a5(!&1%'),.bF3%.*b,%,d)3' R@b
%<% @d*%2 %$35
Cm
'))3X!&12% ',
18
KR
EK%@b%*JXE,.
Vh
H
;<F35$E.F3X0>'%
ƒ8
2 .S}S5W!&1l3%
mS
m
'+('1@2'[%!dW!&1l3%'),.<;<F35T5
,)%,b) ?35. 5 !&1l3%')2E.7*%K** c%(35I*04'
e
<8
%!35I12)* 2c%0,12\63%.*
!"$%'# &)(% +*,.-/!1032
œ
)
T
gÿ
S …!
%3'W5$l!&)$*0>4(%!
ê
*%K*u@.'
!&1%O')
Vê
%]5
l!&
e
*04
35T0>,,2
Fe
35)%15
ê
Fe : ê → e
d;[$E.2T35
Vê := {ϕ1 , ·, ϕdim(Vê ) }.
9
5
')
!&12%O',
'X5]FJ]'
Vh
n
X
g(ue,ϕ) ϕ ◦ Fe−1 , for all e ∈ T .
Vh := uh ∈ V : uh |e := ue :=
ϕ∈Vê
IA]**
ÿ
*%K*[!&123%
',.
Ve := span{ϕ ◦ Fe−1 : ϕ ∈ Vê }
35M'W%!+*%h*'(%!I!&)%0]S
DOFe := {ue,ϕ , ϕ ∈ Vê }
2
,)%,b)'Ž%!.5X')^!&1l3%M',35!&1l3%T',.
5
*%K*
IA
0, · · · , m}
)
'
€ DOFe 
%!?!&%0
S
'120>0>:5.MR'>!&123%C'),.
T
.357@.'
!&1l3%
'l3'
*%@.*d'D%!/!&)2%0RS
;
35R*%@*
'
!&1l3%
Vê
uh ∈ V h
'
ê
ƒ8
'M%* E.2'$
D0
2
:8
@ A]!&123%C'),.
3570
,,2
g
!&)%0
')$%!<*%h*I!&1l3%'
€ DOFe  S
*%h*d)'D%!V!&%0
@b;[
,2,
g
DOF := {ui : i =
%!?!&%0
02
!&%D** )!&)T*043'
*%K*d!&123%O,%wJ2'D''<3%$5
%D3%a%),d%37*%@.*
5.'V3%W,2%wJ
'
Vh
*%K*[!&12% u e ∈ Ve
ue
V
V
*%h* 3%
;<5)O
!#"%$&(')*',+&-+/.0+ 4 " * + ! '& +/&#' 5"% " .21#&( '+/&
;
')
9
%,d3%
L h : Vh → W h .
8
5@ ;^
Lpre
'<0>,,T35c'<%O')l3
Lh
„8
$8
')1,,b%'[35.I?'))3^%,d3%Ž0> T*;^ ' @bX%0,b%'$43%(D*%@*KU(,b3%
.')c%!/*%h*d%,b363%)'
Le
.R$*%@.*u%,d3%
Lpre
Le
Lpost
!&1l3%'S}2S
:
:
:
Lpost
.nS}S
Vh → {Ve , e ∈ T },
Ve → We , for all e ∈ T ,
{We , e ∈ T } → Wh ,
Lh = Lpost ◦ diag{Le , e ∈ T } ◦ Lpre .
c
diag{Le , e ∈ T }
GD%3$35.6
2
;<F35
ƒ6
'<%*u0>3)
„8
2 ]*%K*d;^ ,)%ZJ35
Lh
:8
%0,b%']@ 5
352'(lE.%%!^O')l3
43)'
Le
S
%,b3%W;[>@*T3%O%0T@2$%,d3%'
L2pre ◦ L1post = Id
nS}S
L1h
L2h ◦ L1h = L2post ◦ diag{L2e ◦ L1e , e ∈ T } ◦ L1pre .
;
&"6*.95" , "6 .9+ #'5*"% " .21#&(*',+/&( &#+ 4 "6* +/
%)3%
357@'"3D')),2%]%! 2'l37!&12%'c2

%,b3%)'c@d%wJ ;^7E.2
5
!&%**%w;<2T43)!#(*'')'
œ S
9
h
!"#!#$%'&)(#*!#$%'&
52'(*''?%)'),d%'(3%352T!&1%
KD = 2%0
Ž5>
"#!$%
d =
+
E.*
O32
Šÿ
TS
'),.
"
…
(#*!#$%
KR =
(#*#
n=
+uS
+,&-(#*#
'W,.0>l3:
+
:i8
@
5$%0>
Ž'7;[**['7352>0>'%'
E.*
%!^5
T
% ( & % x
C
2 % % % S
;
h
18
_.1%>'[,.30>l3):a@ 35c K,bW%! 5D!&12%a'),.
!9
€  %S V%>lJ*1.63
KS
:8
i
œœ
T!&1l3% 2352
& #5uH
,2%wJ
a]12'<5(J*1#bS
9
q
!&%**%w;<T0352%KO'
J*1.'X5(!&1%]c,b%4
#$
h
!%'&
bÿ
V
&
52' *˜')' %)',b%' %?35^')l3^!&12%$'),.
5(!&1l3%]'),.
T =
!
,2%wJ'(
18

,)%wJK
$ $
*%K*u 10M@d
<12' 35
#$
$>%O352
@')3l<E.23% S
€@ 
5
&
aÿ
10T@bI%!.35
*# <1'<@.'(!&1l3%]')
K2S h
; ')>!&1l3%C'7,.30:A;<5A518,b
Ðÿ
aF[@b*%'+3%2S
%![F3'W*%K*/!&12%'
ue
##$
&
S
$ ')l3
€

% #
;
$ #$
h
%,d3%
C
'<!&%**%w;<'
18
$
cl31)'c5
$ ,b
Le
"
#$ #! '&
Lh
352
5c*%h* UW,b3%)'
"!#
&
e
%#
u )
) ( gjygjy
18
18
%!Ž43
5$*˜')'
%! !&%0
0>,,
g
;<5
]%1
S
i
€
3
q9
S V%]')'
18
Vh =
5$*%h*/!&1l3%'
18
$
i
35C ,b
ƒ8
%# 0xc"!ac
(#*#
9
S
 Saÿ]23%a352W18K,b
%!$352C!&1l3%'
(*#
3'R%0>
5ngjk87+i!pl~p
3%
18,b
*W%!K
(#*#.S
"7
l31)'V5+'1*F32
m pqi| l~u i!k+*-,.( k:l~i|7l0/‚y l~pqi+| 7+p vj|Dk5k5p k
ÿ=?'5%D43!&% **1%w%;<1'? 3*52?A*'')='^="!&? %[%'0T36@3: '))'@b)]?%3,b5 ,3)%)J'^%>1' 352'D!&3%%)0 'D)B (1,2>')!&4%3 5S =

9 5 ; *'')'c]LO€ N]P]Q H
&
1# 21 1!$ 3 4
h
œ S
;D'
'))@b
LK [LK−1
%!
')l3%
2
1
Li
gÿ
1
œ c5
D9
[· · · L [L [U ]] · · · ]]
'1@%,b3%
S
i
%,b3%
L
'R%0>,b%')
5'*''')J'O'
76
S …(%43)%*'c35T K1K3%
\2|SS(Ž%0M@' 35^'%*123%T%! 352
9
35R@.')
0>l35%
>3%$35'<
€  1 3
!$ Bc'†8lE.
9
S
52
*''
*˜')'
A
43%
*''
,.'')'
\ '<;^%
Li
'
(Vi−1 , . . . , V1 )
>')1@d*˜')' ;<5\5
L[U] =
''0T@*'?5
;<F35T5
.3
,%@*0
U
'?3%>%wJ)

30>,*3?,.30>l3)'
;<5\5R,2%wJ'<5
%h*
!&%',b˜*:J)'%'
%!I35M,lJK%1'W,'')'
,lJK%1'I,.'''
52Tl31.*V%0>,1K\%]%!+,.''D'c0,*0>43]
5
%!
&'h8%,d3%+'+E
'+FJ S /%M,,* T352D'?%,b3%5
' 35+%,b3%
€  2 ' l3(!&1 l3'%;?,,2*
%!V18,d
*'&
&
:8
@
%!/5W!&1l3%'<F3'
:9
"!a%!Ž
*%h*b!&1l3%
!"
'R,.30>l3):
bS
%O352M!&%**%w;<
['[*')%7,.30:;<F35>5D ,b
5$43F3'
q
18
%0,1233%.*
%!^35>')l3!&1%A',.
F3%
A%
Vê
##$
5(!&%**%w;<T05%]' ,%wJ2
€ 
%!/5
5 1'bX%)'),d%'D%$35
*

*%@.*
18 S9
43
18
235( ,b
Vh ⊂ V
')
Vê =
%!I5M) S e
5M ,dM%!+35T@'M!&12%
:8
i
S …I'/,.30>l3):T@ M5^ ,bX%!
h
O')1\5O35.
F33%c!&%c5$')'W%!I5M '
05%'<
€ 
V =
%
[@d*%'
18
,d4D K,b'S
S
!"$%'# &)(% +*,.-/!1032
œ 5 ,lJ%1'<,.''<0,*0>432
S
€ @  11 11 d9 ' 8K,bF 35%!/M]
Li−1
'1@218,ba%!1 R% -h')>352 l31.*,''7' 35 E.'"
,'' %!V18K,b
1 i S
h
1 '"J'W'<]2R0>
A X%5 *'" %!/,.')''S
h
i
;X+)1 30Z352c'"3)11<%! 352c,'')'+'0M@*'[(*!A$*')S=?1 % 5X30>,2*˜X@.')
0>,2*0433% Z35<')1l31)D%!d5'+*!A *'"+h@bclJ*1.63%0>,2*98…0Z'% 35.6
lJT
8 ,''}A %w;<' 5X26 l 18 ,d^%! 'Ž,)2W,.')''+$5.'Ž'' %?352^18 ,d' E
@:a
8 50]S
#
#
*
#
21"
1 "
S
!$
h
1!$ 3 +
&
1# !$c')3'
%!^,''
Li
;<55
i
+ 1
i
!&% ')V%!KI*˜''' ;<5\5
.3O%@2\
%00>%Ah')'T;<5\5
8
*
+ ')1*7!&%0
56
!&)%0
352
5
0,*0>4u352lJ*1.63%W%!35+%,b3%
,.')'S IEA*˜'''7 K')(!&%
3')!&%0>3%A%!c5a%wJ3**[,2%@*0
B
%T' ')0]S+_.%M')3 ',bE.0>,2*0433%'(3%
€ 
;<35
.1 R!&12%
'%1273)0
hJ*@*DOLONOPRQTS …!/5(,%@*0
5O1'
1
hEO5
h...i 21!$ +
1 ! $ + KS
A
\5
43%
RE.'"
1.63%'$%! 352a!&%0
R!&%D')%*J2>*V**,2(,)%@*0'
FJK
E< %
1-
!&%0
%(%!V35'
%ZJ)2
h63%'
352]05%
nÿ
$
'"\2'Ž!&%ID'+%! 1'"82E.
$
!&%Ž35X!&12%'
*
:8
S
lJ*1.a
<5.O%K'<%
%w;<*˜')''@
S
0%'"<h')'5
?
–Fÿ
76
,lJK%1'
*''
18
%'')'c%!/,2%@*0$8'),dE? ,d
S F
UW,d3%
*˜''' ;<5\5$2'@dX5
A
B
M
LONOPRQ
1.%
'M)12
,.')'M;<5*>5a0%>43)h63
€ œ$ % S
,
l31.*K,%@*0RS lE.%2'<]O0,*0>4\%
c5
d%)Ž3%
.,K 35<=D'%4312%1'
;<F3aR'0,*>,)%@*0
2\*'7)a'5*
!&)%0
E%A*˜''7!&%
35a12'$C352a,E
*˜')''S
0xc10>c
9V%
u )
) ( gjygjy
**1'"335%w;
'
')5%w;<
'3%
CB
.9
5)S
m pqi| l~u i!k+* rx|
5R*˜')''a'))@b
76
œ Ž12'2
')%
)
v3p /)
vjp*) p y l~|}l~gjuDy
2S{-]1') ^A'0>,2*E
0,*0>4\%
€ œ'
E F3a-,u''3%A!&%M'0,*F18S 9 5')0,*FE.-!&%0
.
5>%> 20>,2*
13%Oh'
0,*0>43'
]')0,*FE.
J')%A%!^5
gÿ
,2%@*0
E
%!
57%)M%!/35
]35
'0>,*'"c!&%0] %,b3%<;[%1*
1 ,d%43)'b'CB1243**ƒ8OJ*13(3520 ;<5R352M.3 2 35 '1*F3'c!&%0
16'WO)312O35 '1*F<!&%0 35 *˜')<%S9 5 '%*123%0,*0>4352 '
%4\
>*'"D%!
,
,lJK%1'
*%')
1S
0,*0>4\%'<5M*˜')'
5
€ 
∂t U (x, t) = −∂x · (aU (x, t) + ∂x U (x, t)).
9
5
5ngjk87+i!pl~p
C
%$35'X!&%0]K;<F35O;^% d)'
>9
œ S
527'<%*%')%,b3%
B
')
<8
)54<**T% '
>9
S
1 6'?'
527*')X%! 2!&%)0
5
ƒ8
1 * b5M''0T@* %!Ž5
*˜')' @12c3527*')<,.''?')"J'?'<35
%,b363%<'*!†S^%8
&1>O352M25)73%$5
%)l?,12<\$3%h\5
8
0>.O@
1
35
!A
*˜')'<3')*F!†S
18
0>,2*04O'c$* R*')XJ*1@*Wc%0,*l8|30c'%$35.6c ,b
%h@dW)a%wJX!&)%0
'))3(!&123%'<%!V35
.6\
1 6'S9
,lJK%1' O%!V35
)'12*3'
38
18
52'X'^')'3 'D35W ,bW%!
-1T3% 1S
J)'X!&)%0
T
% ( & % x
C
2 % % % œ ;D'X%*',b%')@*1835 1 7*'''X3%7,%wJ2?')%3c!&%[ )0>˜<')1*'
<8
9 5?!&%**%w;< *') '5%w;<'<T')0,*FE.0,*0>4\%
;<5\5O?%* 1');<535?,.'')'S
9
%!V352' %
*''S
!"#%$$&('*)+,.-/102"#%$$4365879;:=<>)67?@A9-02"#%$$B&C'D-E"<F)G5H&IAHH8J
"#%$$4&AHKHML
IN%OI*P*RQ
KSNTLU@8AHVH7G'XWY&Z'[-E"<\)]5H&^AH_HZQ=Q@"AH_H76';`Ya;bc
Vd]efg3]5"79:V<F)G7?%@A9-4hU-H=:V<X7A:<)678ikj6@-!c
lmonpqVrsmpt"uq;v8wm
Vd]efxyd_Sz{&Z'|-E"<\)]5H&^AHVH}Q~Q€,)"+A,‚k'„ƒ"5/!-7:Ki€j@-s},)8+A,‚k'„ƒ85/!-7:yi€j@-Rc
Vd]efxyd_Sz{&Z'|-E"<\)]5H&^AHVH}Q~Q…^-=†:[‚U'*ƒ5/Y-7:[iZj@-ˆ‡)9A,‚U'#ƒ85/Y-7:yi€j@-Cc
lmonq;pBu‰Šp‰
{v8
‹wq;po‰ŒŠm
Vd]efg&(A<\'V%ŽGS]$=},)8+A,‚U'#ƒ85/Y-7:yi€j@- 0x‡)9_A,.‚}'*ƒ85/Y-7:yiZj@-J
∗
i}):~A,‚U'#ƒ85/Y-7:yi€j@-%c
Vd]efg&(A<\'VRhU-H=:V<X7A:<)67ij"@- 0‡)9_A,.‚}'*ƒ5/-7:KiZj@-Jˆ…I-†:[‚}'*ƒ85/Y-7;:yiCj@-(c
∗
IN%OI*P*RQ
ˆm
uu;u‰6u_Š"u_
u_‘’”“GuVŠ4uŠmpqp4m
u"•–]—M˜tplm
oŠp6uV‹‹p8
Š
q
Š"u;q;‰™q;nŠVrm8qp
pl{r6Vr6qpl{q_‰4uŠompqp!—›š=‰ˆ’qu‰‹wŠGuuŠ{u;
‹v
p
‰Š’œŠ
ŠˆŠ
;Šˆm
o
q
ЉqŠVrmqp
pl{r"VrGqplnq;pŠp;vŠ
w_“6u;4u;m
’wŠquŠolŠ{u_‰ŠIž
&%AH~H"Ÿ6&Z'[-E6<)]5H&^AH_H_ s@AHHC¡1Q
@"'|-E8<)65H&(AHH¢£Ÿ8@AHH%¡C0£¤-VHy:¢¥Ÿ=¦6¡¥L
@AHHZ§8<~7C<D:"<[AC,y<[¨]-Ÿ¡
b
© P#ªNI¬«U&%AH~HŸ;¡L
e("*{¤-VHy:¢­c
¤-_H[:V¢®W™¦]c
b
mmu
Š"u;q‰¯q;mp
Šqp°u;
’s’"uVp;Šl{lŠ{8u_‰ŠR—±6ug“
_l³²
Š
;Šom
ptˆ
Šqmp
Šqp™u=‰ˆŠˆ´v‰‰8
;‰’wŠr"Vr6qpl
nq;pŠp;vŠœ_Š
p‹vr"‰Š4’;Š%žu;m
’wŠquŠ’"u_p;Šl³—
© ŽGP\e4Ž6ª__Ž6ª(Ÿ¡MŸŽGS]$=},)8+A,‚U'#ƒ85/Y-7:yi€j@-8 A'#ƒz0µhU-H=:V<X7A:V<F)G7iUj6@-( ¤6-HV:"¡
L
o˜t
;v
Š{Šˆmpt"u;qv8m
ˆn"uVpŠ
@"'|-E8<)65H&(AHH¢­§_@AHHGŸA8'#ƒR¡c
¶vGu;’ˆvmˆŠw
p‹v;r"‰ŠB8u;Š4q;‰8uŠ"u‰‹™q;nˆŠˆ’
Š
q;Š
8u‰’
npqrŠw´·¸;mp
Šqpœ
;‰’wŠp_vŠ8wq;nŠˆmpt"u;qv8m
]ž
i}):~A,‚k'„ƒ85/!-7:yi€j@-:V)6:_AG,‚Y'yƒ]ŸX !A'„ƒ£0M@"'D-E"<F)G5H&IAHVH¢¹§,)G9AG,\‚'yƒ!Ÿ¡¡8c
g´qwq_r8mvŠ
Š6uq‰"p
Š’ŠqoŠ6um
w;t"uVpŠv
{nv‰8Š"u;q;‰€ž
9)/!@5_:|-(Ÿ:)8:VA6,‚³'[ƒx02¤6-HV:(¡c
b
p"uV‹‹p{Šw
q
Š"uq;‰°q;nŠVrmqp
pl{r"_r6qpl
© P#ªNI © ŽGP\eº<~7^<|:8<yAC,y<|¨-YŸ¡¥L"b
(ªP © 8_Q
nv‰8;Š"u;q‰8{pˆ
pˆ‰q;Šˆm
pŠwq;nˆŠˆmvusu‰Špn
4vŠrv";Šo
o
u4l
wqŠpm
Y"#%$$&k02"#%$$&8&U02"#%$$g&&&J
f[ªP*"S(e»8D$$4&%AH~HIc
u;ˆ
’npqrsmp’"u‰‹4m
€—­q;‰lˆ’6uVnnp‰8BŠqoŠ
o
mmu
Š"uq‰¯q;mp
Šqp™u;Š
ŠŠv;q;nˆŠ6unv‰8;Š"uq;‰°
v8Šˆm
gŠqv;{uŠq_“‰ŠVrmqp
pl”’"u;pŠ4nv‰8;Š"u;q‰¯uV‰8;Š
’sq;no
‹6uVt‰sq_‰
© ŽGP\e¼@AHVH6Ÿ_ŽGS$~w€,)"+A,‚k'„ƒ"5/!-7:Ki€j@-6 wA8'#ƒk¡xL
Ž6]ªVVŽ6ª%Ÿ_¡Ÿ8A'#ƒ½0 ¤-_H[:¢¥¡c
∗
b
¾Švp‰
o8u;Šwq;nŠp_vŠq;nˆ
mpt"u;q;v"wm
4u‰8;v’"u‰‹”ŠGuoq_‰
…I-†:y‚}'*ƒ5/-7;:[iZj@-w,)G9AG,\‚'yƒŸ_¡¥L
ª_;NªS¼…I-†:[‚}'*ƒ85/Y-7;:yiCj@-8Ÿ_¤-H[:=¢¿0M@"'D-E"<F)G5H&IAHVH¢¹§,)G9AG,\‚'yƒ!Ÿ¡¡8c
b
(ªP © 8_Q
w
Šv
4qVr"mvŠ
ŠGuq‰"
p’‹
Š’BŠq
ˆ;v8
© P#ªNI © ŽGP\eB9)/À@5_:|-IŸVŽGS$~i}):KA,.‚k'„ƒ"5/!-7:Ki€j@-; A'„ƒÁ0
hU-H:<*7A:V<F)G7ikj"@-( ¤G-H:¡W™¦c
&Z'[-E"<\)]5H&^AH_HV s@"'|-E8<F)G5H&IAHH¢¹c
ª_Ž8V;_ezQ
Vrmqp
pl°;Šqp
‹³—2Âv8
r8v8;Š4u‰6uVŠ6u_
urTuŠo“‰‰’’—
xN~OGÃhU-H:<*7A:V<F)G7ij"@¤-VHy:¢­c
∗
bc
5
0,*0> 33%
%!(5
05%K
'l!&)
%
A')1@d*˜')'a%2>%
**%w;
!&%
K
œ
!"$%'# &)(% +*,.-/!1032
-,
:4
1 ?5'Ž352
u 0,*0> 33%'Z3*%7%W35X12' }'+ S R%%wJ 2C%0>,*98…30>7%2')3)1'7!&)%0
-,
351')T-!&)%0
76
!&%**%w;<2$*'"3$'5%w;<'?R 20>,2*
5XE.F3X u
9
5
0>l35%A%')')'
#
%! ')1\5
35
' $(0,*0> 33%
35.^;^%
:8
R K,''?@ a;<)
1
)J
%! -3%
\5C)%l31)'
76
~6
lJ*1.%M%!dW120>* 1
*
%wJ
35
%]5
h**S Gc%3$5.
1O'*F!
*''S
E.23%'2;<5\55(1'
·‰p"uBu=r8m8_rG‰Š
ŠGuq‰wq;no
n6uVp;Š{qp’p ´_Vr"
pq;_r{rv8Šmpqt"u_’w‰vr"p"uwn;vŒwnv‰8;Š"uq;‰{‹
!"#%$$&('*)+,.-/102"#%$$4365879;:=<>)67?@A9-02"#%$$B&C'D-E"<F)G5H&IAHH8J
"#%$$w38h}&6A_H|H{Q­IN%OI*P*B&]AHKH_}&%'*)+,.-/ 0 3G5879;:<>)"7?%@A9-£0µ&C'|-E8<F)G5H&IAHHZJ L
IN%OI*P*RQ
Vd]ef &]AHyH;€&%'*)+, -/ 0 365879;:=<>)67?@A9-¬0µ&C'D-E"<F)G5H&IAHHCJ %AHK-ViZj@-^c
Vd]efxyd_Sz %AHK-ViZj@-}Q=QVi}):KA,.‚k'„ƒ"5/!-7:Ki€j@-1‚€'Xƒ5/-;7_:|iCj@;-Uc
Vd]efg3]5"79:V<F)G7?%@A9-4hU-H=:V<X7A:<)678ikj6@-!c
KSNTLU@8AHVH7G'XWk&]AHKHV€&%'*)8+,.-/10 3G579;:<>)"7?@A9;-¬0µ&C'|-E8<F)G5H&IAHHRJÁQ„Q=@"AH_H76'bc
IN%OI*P*RQ
q;‰8;Špv8;Šqp¼u‰6u_Š"u_
u_‘{
ˆ
;‰’s;p
ŠsŠqp
‹
3h}&6A_H|HŸG&Z'[-E6<)]5H&^AH_H_ {@8'#-;E6&(AH=HY0V&%'*)+, -/£ @'*)+À0PDS% @_…Y¡ Q
@'„)8+,-/±Ÿ@'„)8+k¡C0 %AHK-ViZj@-6Ÿ@"'D-E]&(AHH(¡C0¥ƒ'~<X¤GH<X¨-6¢xŸ=@_…k¡L"b
¶v6u;’ovmŠVr8mqp
pl¯;Šqp
‹wnq;pm
w;Š
u"ž
© P#ªNI © ŽGP\eº<~7^<|:8<yAC,y<|¨-!ŸV¡¥L
IP„$ JI¤-_H[:¢€W"S
hU-H=:V<X7A:V<F)G7iUj6@-kŸƒ'~<X¤GH<\¨-6¢ &%'*)+,.-/gQ=Q'|A7ƒ-k¡c
−
∗
b
(ªP © 8_Q
[d_Szˆ&%'*)8+,.-/ Q=Q )7;Hyi^j@;-”' ,„58†¥0 Yc
qrmvŠmpqt"uV’wŠw
Šv
wt
;v
Š"u;q‰™q;nˆŠq;mp
Šqp4q;‰‹"u‰‹4Šq
ŠGum
Z—2
p‹su;g
ˆu;Š{q;‰Š
8u‰6uV‰‹”Šw
p‹v;r"‰ŠB’
Š
“
Šp_vŠ8
© P#ªNI © ŽGP\eB9)/À@5_:|-IŸVŽGS$~‚}'Xƒ5/³-;7;:|i^j@;- A'#ƒz0µhU-H=:V<X7A:V<F)G7iUj6@-I w¤6-HV:¡ L
¤6-H:Z§_9G,X-"A]'(Ÿ¡c
e]Ž6NGO%F WYa§ eGŽ6N]O]Ÿ¤6-H:C§;H<*¨8-³Ÿ¡¡c
?G:A:=- k-9:=)G'}&%'*)+,.-/gQ=QK¤)/A<Ã7%J†6,G0~†G'Gc
†%, „¦ ~W ¦^§Dc
†6' „¦ ~W a(§Dc
!§;9G,X-8A'(Ÿ¡c
f[Ž%ªŸ"PDS(”<~WR¦cD<|€¤6-H:Z§H<X¨-^Ÿ|¡=cÃ`` <¡xL
PKfŸ6<|Z¤G-H:Z§_H<„¨-CŸ[¡ a¡
@'„)8+,-/ §DƒŸ<]0<y`³a−0F†8'I0~†G,G0~A'#ƒÀ0' ;,D56†À¡c
GD$
' ;,D56†Á§;9G,\-"A]'(Ÿ¡8c
W' ;,D56†Ác
−
~W}Ÿ I¡c
¤6-H:Z−§VH;-:GŸ<]0 Y¡c
WÀ' ;,D56†Àc
b
b
&%'*)+, -/4 ½@'„)"+,-/µc
PDS({ƒ'~<X¤GH<\¨-6¢±c
bc
5
$ 'c!&%XM;<((%!V,)%@*0
1!$
$',bl3J
)
9
9 52 1 *''<0>,*043'
:8 1 0,*˜X)104hS
;<5\5Ch**'W5
10>h* .156 %
*''S
lE.>@ $F3'
'1*F3O1,d3$J%
A
43)h63^2\*'
352>0>,*043%(%!X'),dE>,.'''S
9
6
+9
21!$
;<F35OF3%.* ,d1')% )0
!&%0
C
B 13%  S2GD%
ƒ8 ,.)(35212' 0T12') ;<F3S +8 *)3*72!&%0>% @d%12(35 %5€ (,'')'
'^')'338..0>*<$
8 ;<5\5%0,d%43'^%! 35W120> [18,dW%4\ 352W*J4^2!&%)0 63% S
9V%'0>,*!)8a5T1'd352$,%@2*0 *˜')'('('( %,d3%?F3')*! @8 Y)1')?!&%);^235$**
57 (*'"32$'5%w;<'W5%w;
5M,)%@*0
h
@bM,12?43%%S
5('*˜[*h+3'),d%" ,)%@*0
5.(352'('W5
%
5
57*''
%*
*˜'"<,.')'S
mpq;_r°
qp‹
‰6uV‘sŠˆm
ˆ
;‰’mpqt"u_’Šnv‰6ŠGuq‰8
"u_Šlq;n
;‰
qmp
Šqp4“6u_s;
;‰ov8’su=‰
gŠ8uyr8C•"u‰Š‹p
Šqp!—
"#%$$4&%'*)+, -/»L
(ªP © 8;Q
pq;_r°
nq;pw’"uVnnv8u;q;‰}•m
g’Œs;;m t6ž
"#%$$™&Z'|)]+G,F-/hR< =5H<)]7gL
!
0>,*043'
T
% ( & % x
C
2 % % % œ
IN%OI*P*RQ
KSN L¤)/!A<>7WYa;bc
KSN LZ'|A7ƒ-%WYa;bc
_d]ef ?G:A:=- k-9:=)G';€'DA7ƒ-€J )7H[iCj@;-Zc
_d]efh³<HV98'|-:-"35679:_<\)]7U%¤)/A<Ã7Y0 '#A7ƒ-ZJ{hU-H=:V<X7A:V<F)G7iUj6@-Àc
&Z'|)]+G,F-/hR< =5H<)]7ZŸe]Ž6NGO%Fo-@HC¡1Q
)G- Ÿ-@H]¡L"b
Gu r"Šq’œuKrm8Vr"‰ŠsЉv;r6p"u
{n;vŒwq;nŠ’"u_nnv8uq;‰œŠp_r{’ ŒŒv
F!"#%$$o‚}'\ƒ5/³-;7_:|iCj@;-VJ
© ŽGP\esƒŸPDS(s-G,X- ~:z0PDS%s-8'~<Xƒ :À0
?6:A:=- k-9:=)6'(¤)/A<Ã7%J†6,G0V?G:A:=- k-9:=)G'¤)/YA<Ã7Jg†8'I0
ŽGS$~‚€'Xƒ5/-;7_:|iCj@;- A'„ƒz0 )7;H[iCj@- ¤6-H:¡xL
)7;H[iCj@-B56,c
,.-/Y-7:%}¦GJQFQXƒ8-:"ŸA'Dƒ8¡ J³-E6AG,„5"A:V-^Ÿ"-",X- y: 0~†6,G05G,¡c
−
¤6-H:[W )G56,c
∗
b
(ªP © 8_Q
e]Ž6NGO%F
)G- %c
bc
pq;_r°
nq;pwp
‰"_mqpŠ€•m
g’Œs_
;v%•_t6ž
"#%$$B&^'#)6+8,Ã-/Ài³'DA7H_@)6'„:L ’Œ{V
;v%•_t]ž
IN%OI*P*RQ
KSN L¤)/!A<>7WYa;bc
KSN LZ'|A7ƒ-%WYa;bc
_d]ef ?G:A:=- k-9:=)G';€'DA7ƒ-€J )7H[iCj@;-Zc
_d]efh³<HV98'|-:-"35679:_<\)]7U%¤)/A<Ã7Y0 '#A7ƒ-ZJ{hU-H=:V<X7A:V<F)G7iUj6@-Àc
&^'#)6+8,Ã-/ÀiR'DA87H_@)6'„:IŸe]Ž6NGO%FAZ¡Q
)G- Ÿ_A¡L"b
š=r8m"_rG‰Š
ŠGuq‰q;nˆŠˆ‰v;r"p6u
™;vm“6u‰’^ž n;vŒ
F!"#%$$o‚}'\ƒ5/³-;7_:|iCj@;-VJ
© ŽGP\esƒŸPDS(s-G,X- ~:z0PDS%s-8'~<Xƒ :À0
?6:A:=- k-9:=)6'(¤)/A<Ã7%J†6,G0V?G:A:=- k-9:=)G'¤)/YA<Ã7Jg†8'I0
ŽGS$~‚€'Xƒ5/-;7_:|iCj@;- A'„ƒz0 )7;H[iCj@- ¤6-H:¡xL
)7;H[iCj@-58,60_58'"c
,.-/Y-7:%}¦GJQFQXƒ8-:"ŸA'Dƒ8¡ J³-E6AG,„5"A:V-^Ÿ"-",X- y: 0~†6,G05G,¡c
,.-/Y-7:%}¦GJQFQXƒ8-:"ŸA'Dƒ8¡ − J³-E6AG,„5"A:V-^Ÿ-"'~<\ƒ 8:!0y†8'I0V58'8¡c
−
PKfŸ )"- GJ€¦_¡
¤6-H:[W )]56,c
∗
GD$
¤6-H:[W )]58'6c
!
)7;H[iCj@-E"'Gc
‚%:|U¤(< V @]A8HHÀJ½Q„Q|ƒ8-:"ŸA'Dƒ¡
¤8-H~: WCE6'Gc
∗
−
J-E6AG,„5"A:V-IŸ-"'~<\ƒ 8:Y0y†8'(0VE"'"¡c
−
b
(ªP © 8_Q
e]Ž6NGO%F
)G- %c
bc
IN%OI„P*Q
p‹v;r"‰ŠI´;Š"uV‰
Š"uq‰œŠlmnq;p{qmp
Šqp
Vd]efh³<HV98'|-:-"35679:_<\)]7£RaG0^aGJs€,)"+A,‚k'„ƒ"5/!-7:Ki€j@-c
Vd]ef¿},)8+A,‚U'#ƒ85/Y-7:yi€j@-”R,F)]+8A6,#hk-H:_<„78A:_<)]7ijG@-kc
q;‰8;Špv8;Šqp!—
sq;‰8;Špv8;Šqp™q;noŠmpq;Vr®u‰6u_Š"u_
u_‘4
ˆmpq;_r”
;‰’m
ˆq; ;Š
&%'*)+, -/­Ÿ=e]Ž"N]O%F AZ0|e]Ž6NGO%F-@H0=PDS(1…k¡1Q
…¢€ŸK…^¡%0
@'„)8+,-/½a€Ÿ-@H(¡C0 @'„)8+,-/ UŸ_A]¡0
@AHH;a^Ÿ@AH~H~¦ 0V@'„)"+,-/½a³0#…^¡0
@AHH ZŸ@AH~Ha 0V@'„)"+,-/ Y0#…C¡xL
b
ˆmpq;Vrœmpqt"uV’”uŠˆq_“‰s
mm8u
Š6uq;‰¼q;mp
ŠqpB
;‰’
;Šo
g
´v‰
qmp
Šqpœu=‰4uŠˆq_“‰p"uV‹Š—
© ŽGP\e4Ž6ª__Ž6ª(Ÿ¡MŸŽGS]$=},)8+A,‚U'#ƒ85/Y-7:yi€j@-8 A'#ƒz0›R,F)]+8A6,„h-H:_<„78A:_<)]7ij6@8-8 {¤6-HV:¡xL
@AHH ZŸA'„ƒ½0V¤G-H:I¡c
b
Â8u_‘ˆq;nˆŠˆ‹p"uV’
PDS(4H<*¨-RŸV¡¥L
ª_;NªS…¢c
b
(ªP © 8;Q
Vd]ef&IAHHV?G:A8'D:€€,)"+A,‚k'„ƒ"5/!-7:Ki€j@-J&IAHHV?G:=A"'D:~ikj"@-Cc
œ
!"$%'# &)(% +*,.-/!1032
Vd]ef 38h}&6AVH[H|R&^'|)G+",>-/£h< ~5H_<F)G7 0¿hR<\HV9'[-:-635679:_<\)]7£³a0Za(J 0Á&(AH=H?":=A'#:KiUj8@-GJ
&]AHKH=aiCj@-kc
Vd]ef 38h}&6AVH[H|k&^'#)6+8,Ã-/ÀiR'DA87H_@)6'„:³0¿hR<\HV9'[-:-635679:_<\)]7£³a0Za(J 0!&]AHyHa;iZj@-J
&]AHKH iCj@-kc
PDS(›…¢c
&Z'[)]+6,F-/hR< =5H<)]7”@'„)8+,-/½a³c
&^'#)6+8,Ã-/ÀiR'DA87H_@)6'„:4@'„)8+,-/ Yc
&IAHHV?G:A8'D:=iUj"@-@AHHV¦c
&]AHKH=aiCj@-”@AHH;aUc
&]AHKH iCj@-”@AHH c
KSNTLU¤I< ;@]A6HHW&]AHKH=aiCj@-!Q=Q@AHH7"'Z0V:'KAG78H@"A8H;HKW&]AHKH iCj@-!Q=Q@8AHVH76'bc
!
bc
I12)32
+9
.6\2S
7** 3%352M%')')'M%!X'"\43335>,%@*0
5
J*13%
6
ƒ8
*˜''7
%!+352$%,d3%?E.** R'('(')0,*M'W;<F32
˜*:>5%)l
% &
$ @b?')5?2 *')S D) 35D2!&%0>3%@d%12[5D0l8…')3,,
R23%.*b*''
€
% !$

'\50>W'[5*a
%c,' R35' ,.,bS
PDS(/A<Ã7ŸPDS(A'„ƒ9 0 Gª
A"'DƒEZ¡xL
∗ ∗ vrp{q;nˆ‹p"uV’ˆ‰q’
PDS( …VW³A8:_)%<VŸA8'|ƒE ~a ;¡6c o‰
e]Ž6NGO%F iWRA:V) ŸA"'DƒE
¡"c o;‰’wŠ6uKr"
e]Ž6NGO%Fˆ-@H[WRA:V) ŸA"'DƒE ¡"c t"uquVŠl
e]Ž6NGO%FAyWRA:V) ŸA"'DƒE ¡"c ˆ
’tŠ"u;q;‰œtquVŠl
p
Š"uq;‰Bq;nˆŠˆmpq;8Vrœ
;‰’Šsqppm8q‰’6u_‰‹wq;mp
Šqp
&%'*)+, -/ @'„)8+ZŸA}0y-@Hk0„…U¡c
p
Š"uq;‰Bq;nˆŠ’
Š
ˆt;Šqp
&%'*)+, -/ Q=Q;},)8+A,‚k'„ƒ85/!-7:yi€j@-5;¦%ŸK…}¡C0 5GŸK…U¡c
p
Š"uq;‰Bq;nˆŠwŠ6u=r68Šmm]uV‰‹wq;mp
Šqp°V‰q;Šw’n"u‰’pž
38)6' A8'#¤ 58,Ã-'€&%'*)+, -/ J' Ÿ@'*)+ 0[¦€§[¡"c
wš_‰Gu;Š6u;
8u;‘
ŠGuq‰q;nˆŠ’
Š
<#76<\:hRA:_A}ŸV5U0¹5;¦€¡c
"uKr"wqq;m
%PX\Ÿ' §K:<\/À-8Ÿ#¡V^i¡xL
' Ÿ5€0y5I¡c
H:¤zQ~Q9-6''zRˆ' §K:=<\/À-6Ÿ|¡£½BH:¤ Q=Q-7¤%,c
b
¸vŠmvŠ{p_vŠ
)]5:;@65: C-H56,F:Ÿ_5^¡c
b
" '(
"!#$&%'(
)+*(,.-/-/0214335/*657*$33598;:"<=8?>@3359*(AB,DCFE2*G>"8IHF3J8LKL*GM†8INO8?>P3Q*$RS*GTUMV1;MV3J*$5W:9*GMQX$>QE2YZ*$:[E]\^H;_9Y/MQ*$X$3JE28?\a`Fbdc/bfeg8?>
:/*'hg\/Ei3JE]8L\j8.<#33 59*k,lCL*G>J,.T?*G:jCL,.0]_/*$M ,.\/:+\m_/57*G>PE]XG,L0ong_;p;*GMR#*7>Q*'<=*G>q3J8+rtsu`uvbxwyE]57*z:/E]MPXG>Q*'3JE :l,.3QE]8?\aE2M
:/8L\/*6Ym1,.-/-/0]E2Xl,.3QE]8?\k8L<oE25k-90]E]X$E23S{_/\/T?9* 8INq_93Q3J,z5k*'33 5/8F:/MGb
| MO3J*$MP3}XG,LMQ*6<=8?>33 5/*57*$33 598;:R#* \X 5/8;8LMQ* ,L\[E2\/E23QE~,L0ZXG8?\9h/T?_/>J,3JE]8L\^E2\^3€RS84MQ-/,LXG* :/E]57*G\/MPE]8?\9M<=8L>
35g,.33 5/*#MP8?0]_;3JE]8L\ 8L<g3 5/*@AB,lC;E2*G†> 8VHF3Q8?KL*$"M 8INO8?>P3Q*$RS*GTqMV1;MV3J*$5‚E]M#ƒG„9…†'‡ 8€Km\/8DR\ˆbw 5/E]M‰5k*G,L\/M‰3 5/*#MP8?0]_;3JE]8L\
E2MBKm\/8DR\[3Q84*'p;E2MP3q,.\/:[XG,L\[YŠ*(X$8?5k-9_93J*$:^C?*$>P1[,LXGX$_/>J,3J*G0i1"Y1[,7:/E2‹Z*G>P*G\3}MP\X 5/*G57*LbZeE2T?_/>Q*Œs M 5/8DRM
35/*q>Q*$MQ_/0i3y8L<x33 5/E2M@XG8L5k-/_;3,.3QE]8?\ˆbFw 5/*q02*$<3y-g,.>P3y8L<x3 5/*ŽhgT?_9>Q*qM 598uRMy3 5g,.3@3 5/*q*'p;-Š*GX$3Q*G:8?>Q:9*G>y8L<&3 5/*
57*$3 5/8;:^RE]0]0YŠ*>Q*G,L\X 5/*G:‘‘33 5/*(>QE]T 5m3BhgT?_/>P*:/*$5k8?\9MP3J>Q,.3J*$M}33 5/*(*$’7XGE2*G\/X'1“8L<”33 59*5k*'33 5/8F:“Ebd*.d
b 5/E2T 5/*$>
8?>P:/*G>57*$3 5/8;:9My0]*l,.:"3Q8Œ5k8?>P*q*$’7XGE2*G\3O5k*'33 5/8F:/M •-9>Q8DCFE]:/*$:3 5g,.3y33 59*MP8?0]_;3JE]8L\E]MyMQ578F8L33 5f–b
y—@˜"!#
rtsGvz™ šŽ›Žœo&žoŸZ ‰¡x¢¤£}›(š¦¥V§šB¨f©‰ ª•IsD«?«?¬?–w 5/*"0]8FXl,.0#:/E2MQXG8L\m3QE]\m_/8?_9M7­,L02*G>PKFE2\®57*$33598;:®<=8?>3QE]57*l8
:/*$-Z*$\/:/*G\3BXG8?\CL*GX$3QE]8?\8€:/Ei‹‘_/MPE]8?\"MV1;MV3J*$5kM$baH;¯ |}°²± bgAO_/5k*$>Gb | \g,L0³bi´µ;¶f`c?c·u¸9`c¹?´9b
!"!#%$'&()*+,-#.&0/%*1%,-#+
2
s
3547698;:8%<=>@?
>@47AB= h C =DE>@F9> L G =D@DEHIDKJ!L@?9=NM1OQPIRESO747T9=>U47T9VW4XRBPYLE=ZL@?9=Z=\[W]%=RLE=V^HIDEV1=D_`PYT9V L G =D@DEHID C =DE>@F9>
acbed G L@47fg=8
rd`vi}
h šx™kjmlneopj ¢¤£;q šsr(¡Zœtlužtux£›(š;v@¢mwcx
uŠ¡”£eyzšez|{~œt} /žxœo¡x¢S£5€7š‚Kuf Ÿƒj& suf¡¤£sq še„©s{FŸ5uf¡5w‰¥
u ¡”£q šOy m…Z†‡ •`L·?·Lc– w‰8DR#,L>Q:/MŽ,ˆO\/E2hg*$:“eg>J,.5k*'R#8?>PK4<=8?>ŽH;XGE2*G\3JEihgXЉ#8?57-/_93JE2\/Tgbg¯€\Œ‹ >Q8FXLb
f
8L<3 /*Œsuµu3 "¯€\3J*G>P\g,.3QE]8?\g,.0e#
‰ 8?\9<=*$>Q*G\9XG*8?\ŽB 8?54,.E]\.Ž *$XG8?57-Z8LMQE23QE]8?\ ° *$3 98;:ˆb
r ´.vŠU
 š™B ¡x¡s€o £eU štru‡‘ ¢sŠu ¡”£cŒy šcz|{~œc} / žxœo¡x¢S£cq š’(„ ©s;{ ŸƒŠu ¡ƒwfu ¡ •³`.·?·µ?– | \[*$’7X$E]*G\3BE]57-/0]*$5k*$\
3,3JE]8L\k8L<,.\,L:g,L-;3JE2CL*},L\9:U-g,L>Q,L0]02*G0gT?>PE]:4E2\^NB
ˆ Ay“ b9w*GX ‘b9>Q*$-ˆbigH;_/Y95kEi3Q3J*$:k3Q8/¶%‹ >P8;X$*G*G:9E]\/T?MS8L<
w /*U`.\/: {_/MQMPE~,L\ V­q*G>Q5k,L\ | :9C.,L\/X$*G: {*GMQ*G,L>QX ®)+8?>PKFM /8?-®8L\”#
‰ 8?57-/_93,3JE]8L\g,L0”H;X$E]*$\/XG*,.\/:
• E2T .‰‹ *$>P<=8?>P54,.\/XG*Š#‰ 8?57-/_93QE]\/T//HF3J_;3Q3JT?,L>P3G ° ,L>QX jslc OsD¹9b
B
r cLv—U
– š;rNjc{.q˜jmul œ@£e}h še™tuƒš›F{ œo&©S£5Bš š‚q ‰¡;j+n •VsD«?«?µ?–Ž *$hg\9E23JE28?\k,L\/:4RS*l,LK7MP3J,LY/E20]E23€1(8L<&\/8L\/XG8?\
MQ*$>PC.,.3QE2C?*6-/>P8;:9_/X$3QMGb ± b ° ,.3 ‘b‡”
‹ _9>Q*GM | -9-/0³biœ2 c/¶/c¬?´D¸9µ.c¬9b
rdµviz
y š+™tu+iu‡ož su •VsD«L«?·–kAO_/57*G>QE2Xl,L0 ° *$39Ÿ 8;:/Mx<=8?>m#‰ 8?\/MP*G>VCL,3JE]8L\x ,lRMGb¡& *GX'3J_/>P*GME2\ ° ,.39Ÿ *G5k,.3JE2XGMG
¢ E2>QK9Ÿe,L£ _/MP*G>$/hg>QMV3B*$:/\ˆb
#
r ¹.vi}
h št œt{¤{¥…œœt} {¥{¤fu ¡ •³`L·?·.cm– | M\9Ÿ 8;XK“XG,L-93Q_/>QE2\/TUR@,lCL*k-/>P8?-g,LT?,.3JE28?\“57*$39Ÿ 8;:a<=8?>}:9>P1 ,L\9: n/_/E]:¤¦
MJ,3J_/>Q,.3J*$:^T?>Q,L\m_/0~,.>#ng8DRMGb ± b§#
‰ 8?5k-9_93lb‹Ÿ 1FMGbi®sD«?«•Isu–'¶ˆsDµL·u¸Šs2 c/b
r¨2 vi}
h šO™kjmlneop‰j ¢¤£™ škyZo!©tofu ª ¡eu •³`L·?·.cm–KB E]MQX$8?\3JE2\F_98?_/MT,.0]*G>PKmE]\j5k*'3/Ÿ 8F:/M <=8L> 3€RS8€¦ -+/Ÿ ,LMQ*4n/8uR E2\
-Š8?>Q8?_9My5k*$:/E~,9b;w‰*$Xˆ
Ÿ bŠ{O*$-ˆbŠ`L·?·.c¸9`L¬9;¯ ) { •³H9ec¢ ´?µL«–'+Bˆ \/EiC?*G>PMQEi3,.£ 3kO• *GE]:9*G0]YŠ*G>PT/b
r ¬.viz
y š’«Ku‡{¥…Zo¬w •IsD«?«1L2 – ° _/0i3JE2-+gŸ ,LMP*e028uR ,L\/:aw>J,L\/MP-Z8L>P3Z‹ >P8;X$*GMQMP*GM6E]\3/Ÿ *zH;_/Y/MP_/>P<,LX$*L¶ | XG8?\¥¦
3J>PE]Y/_;3JE]8L\43J8Œ3/
Ÿ *6578;:9*G0]E2\/T(8.<tŸ 1F:/>Q8?MV1FMP3J*$5kM$bZH;-/>QE2\/T?*$>Gb
r «.viU
z šƒz­ } ©¢S£ƒq š‚„ ©sF{ Ÿ5fu ¡5wfu ¡”£ƒZ® š&¨f&©‰ s…N&j &©sŠu ¡”£§zy š;¯eo¬u§w{¥fu ¡”£5zy š5zi{’œt} / žxœo¡x¢®•`L·?·?´?– | :911¦
\g,L57E]Xy3€RS8€¦ -%g
Ÿ ,LMQ*Bng8DR 578F:/*G0f8L<x-/>P8L3J8L\4*'p;XgŸ ,.\/T?*q57*G5Y/>Q,L\/*<=_/*G0ZXG*$0]0]M$b°”‹ >P*G-/>PE]\3y‰ H%‰ |q°“°
{*G-Š8?>P3O·?´I¦ ·1;
2 /H;_/Y/57E23P3J*$:"3Q8z3/Ÿ * `.\/:±“’Žˆ {`²Z‹’“ | A³‹’”“ em‰ emq² {Uˆ ° ‡& _ G*$>Q\ˆb
rtsD·.vŒ
€ š`@h œto ¢sul œc”n £`4r š’iu%´ ¢e‰j ¢sneu •`L·?·9sD–w`/Ÿ *G8?>P*$3QE]Xl,.0S,.\/: AO_/57*G>QE2Xl,L0~#‰ 8L5Y/_/MV3JE]8L\ˆb{b wqbs”“ :¥¦
R#,L>Q:/M$g¯€\/XLb
rts?sGvi4
r šc( ¢m‘mŠu ¡eul œo¢¤£~U– š’q sš‡jm‘c‘mŠu ¢¤£m4 š§ª©su‡u‡{¤fu ¡ •Vsl«?«?¬–Ž E2‹Z_/MP*(E]\3Q*G>P<,.XG*(57*$3\/Ÿ 8F:/MŽE2\^ng_/E2:
57*GXg
Ÿ ,L\9E]XGM$b | \/\ˆbf{*$CŠbfe02_/E]: ° *$XˆŸ bia´?·9¶ˆsl´?«u¸ŠsD¹?µ;b
rtsu`vi4
r šZro¬fu ©s{ •³`L·?·?µ?– ) *$0]0Yg,L0],L\/X$*G: :/E]MPXG8?\3JE2\m_/8?_/M4­,.0]*G>PKmE]\ MQX/Ÿ *$5k*$M4<=8?>z39Ÿ *[AB,DCFE2*G>µI¦ HF3J8LKL*GME¦
Nq8?>P3Q*$RS*GT7*Ym
¶ _g,.3QE]8?\9MGb±”‹ >P8;X$*G*$:/E]\/TLMO8.<t•`·g‹ `L·L·Lc/¤¤· 8?K.8WgŸ ,L5k,—‹ _/Y/0]E2M\/Ÿ *$>QM$g¯€\/XLb
1
35
1
35
35
K8
\5
5
8
\5
5
)5
\5
T8
8
35
: