Discontinuous Galerkin and the Crouzeix-Raviart ... - FEniCS Project

¢¡¤£ ¡¦¥ § § ¨©§ § £ ¥ § £ ¥ §
PREPRINT 2000-09
Discontinuous Galerkin and the Crouzeix-Raviart
Element: Application to Elasticity
Peter Hansbo and Mats G. Larson
Chalmers Finite Element Center
CHALMERS UNIVERSITY OF TECHNOLOGY
Göteborg Sweden 2001

¡ ¢ £¥¤ ¦¨§© ¦ ¦¨£¥¦¨¤ ¦¨ ¦¨ ¦¨§
Preprint 2000-09
Discontinuous Galerkin and the
Crouzeix-Raviart Element: Application to
Elasticity
Peter Hansbo and Mats G. Larson
Chalmers Finite Element Center
Chalmers University of Technology
SE–412 96 Göteborg Sweden
Göteborg, January 2001

Discontinuous Galerkin and the Crouzeix-Raviart Element: Application to Elasticity
Peter Hansbo and Mats G. Larson
NO 2000-09
ISSN 1404–4382
Chalmers Finite Element Center
Chalmers University of Technology
SE–412 96 Göteborg
Sweden
Telephone: +46 (0)31 772 1000
Fax: +46 (0)31 772 3595
¡ ¡ £¢¥¤¡¦¨§©¢¦¢
Printed in Sweden
Chalmers University of Technology
Göteborg, Sweden 2001

© ¥§¨ ¡!#"$&%¢¡'(
¢¡¤£¦¥¨§©¡©§¨£¡©
§ 3£4¡¤¥¡'(5
)*)+,)©.-)/0/*¡¤¥1(0¡2§¨©
ª:«:¬

¿¾¡¿Â.¿³\º$·)É&»º‹Ã˳i¹Pµi¿·)·$»º$¼J´]º$¿'µÎ´]ºº‹Ã ¿¨Â¨¡µ Ä“mš›
ê ¿À$¿9‹0ìœne" p
ž8sr}]ì`nŸ"¦× t su}€vx ¡·Œº‹Ã ¿0·Xº‹À$ŕ¡³º$¿'³i·{¹VÀ‡É&–º‹Ã¢¸'¹rÂ.Ä8¹V³i¿³\º$·
Ô
ž^s4}Qìœn" p
}€vx “¦{|s4}@©“|¢¥}kª t
s«vx
ƒ¤„Ry¥p
¦C§¨t
Ø•†¬p
×~t su}®vx É&»º$ísr}¯p
Ô4­s4}
»Ç°NpŽ±àŕ³iµ0­ks4}Gp³²–Ç°1´ pŽ±JØ
§
ŕ³iµ˜ ŕÀ$¿¨ÏV–æV¿'³Î¾¡¹Ž´]µi·'ئŕ³iµ – ¡·)º‹Ã ¿.¹V¼ º¦É}ŕÀ{µä¼i³i–º0³i¹VÀ$Â.ŕ¾ º$¹>“m ¨¢Ð ¼iÀ{º$Ãi¿'À{Â.¹rÀ$¿rصäŕ³iµ‰
Ú ŕÂ0· ¿¸¹V³i·{º‹ŕ³Žº‹·Ø3·$´]º‹»·{ÇaÆP¡³iÏ$²¹¸¤‰x¯¸³‰
ŕÀ{¿#º‹Ãi¿
¸¤‰
ŕ³iµ¹²¹¸º»¸½¼¨›R²¦³äº$¿'À$¨· ¹rÇ&º‹Ã ¿
¨¹ µi¼ ¾¡¼i·&¹rÇ—¿¾žŕ·Xº‹»¸'»º¦Ærؾ.Ø

¾?À|©iì{ì § ‡—À|"'ì § "Œ ‹kìb9"€"|+
‡¬

¨ ÕR" ì%<
É&»º$Ã
?r¾P?Žº$Ãi¿(¾¡¿³iÏrº‹Ã¢¹rǵ¾¯=E)Û¡·Œµi¿'J³i¿'µÊŽÆ
¾P?r" Ç‚¹VÀo¾ ?!‡I?% Z ?¤£D©iìi

¡ £¢¦ ¤¤ ~¦¥ ¦¥N ¢J§¨§© ~ ¤|¨¢¢\§¦¥©^¢ ¨ §¦ §¥¤(¢£Š¢\¢J ¤
Œ…‡l@Œ…‡¦xkŒ ‰ R^,-#3? A
]^0=+T._A áο'JÀ$·Xº‡³i¹rº‹¿)º$ÃJ´]ºŒº‹Ã ¿(Ç‚¹r¾¡¾¡¹¤É&»³iÏ1¡³Žær¿'À$·{¿¿·{º$¡Â.´]º‹¿)Ãi¹r¾¡µi·
)
Ó ¿‰ÅPºxØO¼i·$»³iÏ º‹ÃJ´Qº _ –
§

é à ¿.ÄiÀ{¹ ¹]ÇM¹rÇì%< ¨»§ Õ‘")¸'ŕ³ÍÊ

£¢¦ ¦¥N ¢J§¨§© ~ ¤|¨¢¢\§¦¥©^¢ §¥¤(¢£Š¢\¢J ¤ ¨ §¦
ŕ³iµ¢º‹Ã ¿'³Û¼i·$¿)º‹Ã ¿&#}´]¼i¸ ÃŽÆ\È-+ ¸ ÃŽÉ}´]À$ß ¡³ ¿'Ü\¼Jŕ¾»»º¦Æ¨º‹¹ÏV¿º
~¦¥ ¤¤ ¡
nŸ" [ ?DC&s £
¡³i»³iÏÛì%< ¨ < è "}´]³iµÎì%< ¨

"$#9#5©$5 ,-20£$#n¦)ì%< ¨ 5‘" T.!

£¢¦ ¦¥N ¢J§¨§© ~ ¤|¨¢¢\§¦¥©^¢ §¥¤(¢£Š¢\¢J ¤ ¨ §¦
IP£IO¬I8ªE° IB
~¦¥ ¤¤ ¡
ö£¢ G3L @Lm? lo`2jygtŽ… ? `}ez`'bdp{ldegjmliqrpX`Qsmyƒb‚„*¢Q`2ezbdp4pXygpX|}pX`'bJ|}p~boc2j¤t3h:ezbdcŒt¤eufov~jm`'boez`xw2jwQfPp{yzp{|}pX`'bof{… B _ ?AH ò]L‰@ wQ|Mp{l L
¡
`Qs‰y L ö{ô•‚ö$ô¥¤'óm˜§¦ ¨ ó ¦ ¦mõ L
?
ó¨¢ KML ?*L‹C s‰nmp{l{… ez`2ezbdp

¢ ± ¥ Ē¬I¤£¦¥¨§ Å ŕ¸mº ·{¹V¾¡¼ º‹¡¹r³ ¨
± ¥ Ē¬I©¥ §‡³i·{º‹ŕÊi¾¡¿·{¹V¾¡¼ º$¡¹V³ ìˆŔ"‘¸'¹VÀ{À$¿·$Ä8¹V³iµi»³iÏÛº$¹¸ à ¹ ¹V·{¡³iχ¦xU¸9¸ ‰ Ø¦ŕ³ µ
¢
º$¡¹V³Ç‚¹rÀ˜‡x/p»‰ ¨
ì‚Ê"M·{º‹ŕÊi¾¡¿(·{¹V¾¡¼
¡ ¡£¢¥¤¦¢¥§©¨¦¥¤¥§ !

¡ £¢¦ ¤¤ ~¦¥ ¦¥N ¢J§¨§© ~ ¤|¨¢¢\§¦¥©^¢ ¨ §¦ §¥¤(¢£Š¢\¢J ¤
¢ ± ¥ Ē¬I¡ ¥ + ¿¸'¹V³iµ¹VÀ$µi¿ÀR¸¹V³ŽæV¿'À{ÏV¿'³ ¸'¿Œ¡³(º‹Ã ¿˜<
t
¹]Ç µ ¿'ÏVÀ{¿'¿'·¹rÇ&Ç‚À{¿'¿'µ ¹V »·1¸'¹V¨ÄJŕÀ$ŕÊi¾»¿.º‹¹Ë´à¸¹V³ Ç‚¹VÀ{Â.»³iÏ,'J³ »º‹¿
º$¹Ìº‹Ãi¿³\¼iÂ1Ê8¿'À
µ ¨
¿¾¡¿'¨¿'³Žº&¨¿º$Ãi¹
¢ ± ¥ Ē¬I£¢ ¥ Ó‡¹V³i¸'¹r³ Ç‚¹VÀ$¨¡³ ÏÛìdŔ"}ŕ³ µ¸¹V³ Ç‚¹VÀ{Â.»³iÏ¢ìˆÊ¦"Œ·{¹V¾¡¼ º‹¡¹r³i·}Ç‚¹VÀzÀGp‚²\› è Õ ¨
-7
-7.5
Present method
Conforming linear FEM
-8
Log(L 2
-error)
-8.5
-9
-9.5
-10
-10.5
-5.5 -5 -4.5 -4 -3.5 -3 -2.5 -2 -1.5
Log(1/NDOF 1/2 )
ìimo"
³i¹rÀ$ Äé Ãi¿M¿'À$À{¹VÀ¦À$¿'¾¡´]º‹–æV¿
•»s˜ @ j`Qv~j`¤€jmlo|}ez`Q”
• † ß

¢ ± ¥ Ē¬I¡ ¥ Ó‡¹V³i¸¹V³ Ç‚¹VÀ{Â.»³iÏÛìdŔ"}´]³iµä삾¡¹P¸ Ár¿µ"}¸'¹V³ Ç‚¹rÀ$¨¡³iÏÛì‚Ê"M·$¹V¾»¼ º‹»¹V³i·}Ç‚¹VÀzÀGp¨²\› è ÕVÕrÕ ¨
- ¡£¢¥¤¦¢¥§©¨¦¥¤¥§ !
•»s˜ @ j`Qv~j`¤€jmlo|}ez`Q”
• † ß

ümî ÿrï ‹ð þ¡ ïaðzð£}ðXþQï0ð~ï 2ú$ü‹ý~ú‰û—ï Qð þ2ý¦ïžð îü ø0îù8ð£ýQî ï –ú þQý
©
jmc]s‰` = jV|Œs‰` ò
© ƒï © ümî ÿrï ‹ð"!Oí#$ &%Vú '‹ðXû
¢¡£¡¢¡¥¤¦¡
`Qt¤p{lof N jm”m” ?
Öï © ümîoÿVï ðkû¦ûdú‰û*)Jú‰þQï¡ûdú +~ú û,!OíEî zðXû@? þ Œð~ï 2ú$ü‹ý$~ú û þ£. ú ¦ÿrûoð¦ý ý /2 gð î‰þrü‘þrð î‰û Öø þ£. ú ¦ÿrûdð~ý ý /2 zð
¢¡£¡¢¡¥¤¦¡£<
ý ï /.£ ï»ø"2XøBA ï»ý. Qð3ý ŒðXï Qú{ü
ð£zî
pXbdp{l = s‰`Qf † j&s‰`Qt H s b f KMLQN smlofdjm`
6
© í uý. úD¤þ]ï þ¢QúDxý>:îzðXû@? þE0ð~ï 2ú$üF~ú‰û*ï 2ðGzî ïaðHQî‰ï –ú‰þ
¢¡£¡¢¡¥¤¦¡C
pXbdp{l = s‰`Qf † j&s‰`Qt H s b f KMLQN smlofdjm`
6
)Oú‰þ2ýXð~ûJ î ï –ú þ,Gkûdú ÿ2ðXû¦ï –ð¦ý¨~ú‰ûEï Qð:)Oú‰þQï þ3QúDxý¦î‰þrüRí uý. ú þ]ï þ¢QúDxý>:îzðXû@? þ"0ð~ï 5
¢¡£¡¢¡¥¤¦¡£I
í uý. ú þQï þ3Qú=xýL>Eî zðXû@? þ î þVüï QðM)iûdúD+N$ð£'O /PAîD –î û¦ïQ zð ŒðXþQïSR © ÿmÿ /. î ï –ú þ ïaú
¢¡£¡¢¡¥¤¦¡£K
zî‹ý¦ï /. ï¡ø
¥&TZY¦fŒV8b ]œa]C[gVÛfœV V¥a[.¥&V¥a[@V./ ^VäŠ]œa[^b
¢¡£¡¢¡¥¤¦¡¨§
; egv n s‰l t C p{ld”'f‚bol
j|
ú$ü‹ý
s bof KML]N s‰l fdjm` sm`Qt ?*LQò jm`]smf @ egnxygsfdfdjm`
H
6 pXbdp{l = s‰`Qf † j&s‰`Qt H s b f KMLQN smlofdjm`
é Ãi¿·$¿0ÄiÀ$¿ÄiÀ$»³Žº‹·&¸xŕ³ÛÊ8¿0¹rÊ º ŕ»³i¿'µÇ‚À{¹VÂ
¡ ¡ ¢ ¤¡¦¨§©¢ ¦ ¢ UT ¤ ¤§5V¨Ẅ