Optimization loops for shape and error control - Sophia Antipolis

0 downloads 0 Views 1MB Size Report
Dec 15, 2004 - 64 B ` ! 4 %$'&. Г#". 6 w 4 ! 4( ). Г0". 8D. This statement shows that the ...... with anisotropic adaptation, a convergence order of %Y86`c % 88.
INSTITUT NATIONAL DE RECHERCHE EN INFORMATIQUE ET EN AUTOMATIQUE

Optimization loops for shape and error control: extended lecture notes A. Dervieux, F. Courty, T. Roy, M. Vázquez, B. Koobus

N° 5413 December 15, 2004

ISSN 0249-6399

ISRN INRIA/RR--5413--FR+ENG

Thème NUM

apport de recherche

 

  "!#  $%!&' )(*! +,!- $ !&$. !&'/0 !  !& 132546879:;68=5?A@CBD2CEDFG†)wysutv{‹ˆM|§wAicŠ”)…Qt‚|¬hxAh¹©¸‰{vhxA¦§ˆMtv|§)w£&|¬tve,ˆ)Šƒ„)|§w‰tc•)wAiª sveA)t™‰…Qtv|§h,|¬¯cˆtv|§)w3ˆM¦§š)‰{v|¬tveAhs™ˆMwyŠ"hxy¦¨t‚|¨ª«¦¬iI¤)iI¦‰…Qtv|§h,|¬¯cˆt‚|¬‰w¥…A{‚iI†)wyŠQ|¬tv|§)wy|¬wAšX²dfeAi,ˆM…A…y¦¬|~†ˆMtv|§)w®t‚ˆ¶sveyˆM…Gi ŠAiIsv|¬š‰w…A{‚)“A¦§ihº{‚i¦~ˆt‚iIŠt‚Dsv)wy|§†l“X\‰h»hD|§wA|§h,|¬¯cˆtv|§)w |§s {‚i…G){vtvicм²5´µw ˆDsuic† )wXб…yˆ){utc•\ˆM©½tviI{fsv)h,il{viIh,ˆ){v¾Qs¿)w h,icsue¥ˆ‰ŠAˆM…At‚ˆt‚|¬‰wZ•`£¿ii ÀAˆMh,|§wAit‚eAiD…A{v‰“A¦§ihÁM©€©¸‰{vhxA¦§ˆMtv|§wAšh,icsue¥ˆ‰ŠAˆM…At‚ˆt‚|¬‰w®|§w®tviI{vhslM©¿ˆ)w®)…Atv|§hD|§¯IˆMtv|§)w …y{v‰“A¦¬iIh¶•8…y{vicsuiIw‰t‚|¬wyšDsv)h,iˆM…A…A¦§|~†ˆt‚|¬‰wys¿tv,|§w‰t‚i{‚…X‰¦§ˆMtv|§)wˆMwyбtv5ˆM{vtv|~ˆM¦¼’Â|¬Ã`i{‚iw8tv|~ˆM¦ZÄC}8xyˆMtv|§)wysI²€dfeAi|¬h,…A¦§|¬ª †Iˆt‚|¬‰wys-)©t‚eAi…y{v‰…X8suicжhDitveAQŠO“GiI†)h,i†¦¬icˆM{c°¿“\r¶† ‰h“A|§wA|§wAštv‰š)itveAiI{“GMtve¥|§ŠQicˆ)sl† ‰sutliÃGic†Kt‚|¬¤‰ih,i tvey\Šys †IˆMwO“GiŠQiI¤)iI¦¬‰…Xicж£&eAi{‚itveAi)…Qt‚|¬h,|§¯IˆMtv|§)wO¦§8‰…ys-|§w\¤)‰¦¬¤‰i”“X)tve¥sveyˆM…Gi‰…Qtv|§h,|¬¯cˆtv|§)w>•yt‚±|¬h,…A{‚¤)iˆ)w®ˆMiI{vQŠQr8ª wXˆMh,|§†ŠQicsu|§š)w>•>ÅÆGǙh,icsueOˆ)ŠyˆM…Qt‚|¬¤\|¬t„r)•\tv±†)…Gi£&|¨t‚e t‚eAiŠQ|¬È±† xA¦¬tv|§iIs&)©>t‚eAiÉ Ê5’Œ…A{‚)“A¦§ih|¬w\¤)‰¦¬¤‰iIм² ËCÌ Ç)ÍcÎÏQÐ Ì ÆyÑGÒXӂÔKÕ Ì ÆyÑ ÌÂÖ Å× Ì Ó‹Õ‹Ø;ÍMÆ®Ç)ÔÑ«ÅMØ Ö½Ö Ô Ì Ç Ì ÕÂÆGÍMÑ Ì ÕlÇ Ì Î‹ÍÏQӋÕlÅMÏO٧ڀÛÜ Ýßޙà)áÂâ3ãXä\ÍMÓKÑ Ë ÍMÏQÓ‹Õ Ì ÍMÆ ÝOÏ Ö Ñ³Ø;ÇØ~ÕKå ÎØ Ò Ö Ø½ÆGÅÓKæ±Ý®ÍcÇ ÌÖ½Ö Ø½Æ\ç±ãGؽÐÏ Ö ÅѳØ;ÍÆ"ÅMÆGÇèàXÅ Ö Ø;Ç)ÅMÑéØ;ÍMƥؽÆá Ì ÓvÍÆ`ÅÏQѳØ;ÎKÕê~ë/ì ÏQØÕAí Ì Õ‹Ñ¿Ñ Ì Æy϶îïlÅÓvÎ ÌÖ ÍÆ ̱ð;ñò ë Ö¬Ì Õ™ó8ôMåµó8õ ö ÏQؽÆ"ó8÷‰÷ øA² ùOú8û`ü‚ýþ Ÿÿ)¢ 5ˆM{vtv|~ˆM¦>’l|¬ÃGiI{viIw‰t‚|§ˆ)¦ÄC}8xyˆMtv|§)wysI•`Ê|§wA|¬tviD¦¬iIh,iw8tm¶itveAQм•¼m¶icsue &ŠAˆ)…Qtv|§¤\|¨t„r‰•`mOxA¦¨t‚|¬¦§i¤‰i¦ l…Qª t‚|¬h,|§¯IˆMtv|§)wZ•XÉ¿‰hD…yxQt‚ˆMtv|§)wyˆ)¦Zʦ§xA|~ж’Âr\wyˆ)hD|~†sI• l…Qt‚|¬hˆ)¦>p\eyˆM…Gi’liIsv|¬š‰w

B      "!$#%&(')*)"!(,+.-0/12,3154678 ':9')7?*)')!( @,A9= &" C         ""!$!$#%#%&(&(')')*)*)""!(!(,,+.+.-0-0//1122,,33115546467788 ':':9')')7?7?*)*)')')!(!( @,@,AA99== &"&" DFE V E == '' G>G>""AHAH!I!I'' IW99° $ !% & é–  dfeyi¤­ˆ){v|~ˆM“y¦¬ics  ˆMwyŠ  ŠQicsv†{v|§“Gi,tveAi…A{‚)“y¦¬iIh £&|¬tveA|§w † iI{ut‹ˆM|§wèsv…yˆ‰† iIs)© £ ){‚¾`² Zitxyss‚ˆ­rOtveXˆt  ˆ){viDt‚eAi †)w8tv{‚)¦Q¤ˆ){v|~ˆM“A¦§iIs€ˆ)wyŠ  ˆM{‚iftveyiÂs„t‹ˆt‚i&¤­ˆ){v|~ˆM“y¦¬ics²dfeAir,†ˆ)wt‹ˆM¾‰i-ŠQ|¨Ã`i{‚iw8tC{‚)¦§iIsCˆ‰††){‹ŠQ|¬wyšÂt‚”tveAi…A{v‰“A¦§ih £ i ˆ){vil|¬w\¤‰)¦§¤)iIŠ,|§wZ(° y£ ·yiI¦§Š¤ˆM{‚|§ˆ)“A¦¬ics/‰{¿ˆ)…A…A{‚­À\|§hˆt‚|¬‰wi{‚{v‰{C©¸){  ˆMwyбsveyˆ)…Xi™…yˆM{‹ˆMh,itvi{‚|§¯Iˆt‚|¬‰w,){fˆh,i t‚{v|~† ŠAi ·ywA|§wAš,ˆD…yˆM{vtv|~† xA¦~ˆM{ ¦§\†IˆM¦`h,iIsve ŠQicsv†{v|§…Qt‚|¬‰w±©¸){ ²CdfeAi™|§ŠQicˆD|§s¿tv,h,|§wA|¬h,|§¯i™tveyi Í ö Ì ÎÑéØ½× Ì ‹ÏQÆGÎÑéØ;ÍMÆ ¦§\)¾\|§wA𩏁‰{ $O£&eA|~†‹e¶s‚ˆt‚|§su·yics tveAi Ì ìÏyÅ Ö Ø½Ñ³æ΋ÍÆAÕKѳÓvÅ ؽÆXÑ   $ !'" ² Zi t&xysˆ)s‚suxAh,i™tveXˆt&tvey)i ()ˆ‰† )“y|§ˆ)w   * !,+ ;%‘   + |~s&ˆM¦§£fˆ­r\s |§w\¤)i{vtv|§“A¦§i)² p8t‹ˆM{vtv|§wAšD©¸{v‰ht‚eAi $.- !% -  -  •A£ i”ŠQi ·XwAiˆMw|¨t‚i{‹ˆt‚|¬¤‰i™‰…Qtv|§h,|¬¯cˆt‚|¬‰w …A{‚Q† icŠQxA{‚i™“\rt‹ˆM¾\|§wAš

 !

* -

!

 $ -   $ - & + .  +

´µw¶){‹ŠQiI{ftv,i ȱ† |§iw8t‚¦¬r† )h,…AxAtvitveyi”{vic}‰xy|¬{‚iIŠsurQsutviIh)xAtv…AxQtŠAi{‚|¬¤ˆt‚|¬¤‰iIsf£&|¬tve¶{‚iIsv…Xic†Ktftv,t‚eAi”|¬wy…AxQt‚sI•y£ i Xx sui™tveyiˆ)Š­ƒ„‰|¬w8t&h,i t‚eAQм•A£&eA|~†‹e¶…A{v¤\|~ŠQiIs š‰iwAiI{‚ˆ)¦¬|§¯icŠˆMš){‹ˆMwyš)iÂhxA¦¬tv|§…A¦¬|§i{‹s&|§wtveAi©¸‰{vh)©5ˆ,sv…yˆt‚|§ˆ)¦`·yi¦~ЁM© t‚eAis‚ˆMh,iŠA|¬h,iwXsu|§)wO)©t‚eAi‰{v|§š)|§wyˆ)¦Z…A{‚)“A¦§ih xAwA¾‰£&wys²-Êy){-icˆ)†‹eO|¬tviI{‚ˆMtv|§)w0/`•XtveAi¶Î‹Íå«ÕKÑ«ÅMÑ Ì ){,Å)Ç ö ÍMؽÆyÑ&ÕKѵÅÑ Ì 1-|§s&ˆ,sv)¦§xQtv|§)wM©   1- ! ¸  2 +  $.- 4365.7 +  $.- 8&   '3iŠQiwyMtvit‚eAi ‹ÏQÆGÎѳØ;ÍÆGÅ Ö ç)ÓvÅ)ÇMØ Ì ÆXÑZ“\: r 9;- ! + 9  $! +   $.- @? A 1=-  +  $.- CB éŽ    ´µwt‚eAi”†ˆ)svi)wwAMt‚eAi{-†)wysutv{‹ˆM|§w8tI•8tveyis„t‹ˆtv|§+)wXˆM{‚|¨t„r)©t‚eAi™©¸xAwy†Kt‚|¬‰+ wyˆM¦>£&{v|¬tvics° 9  $ -  1 - >!D"E

F%

)ˆ wyŠtveAisvi t-)©t‚eA{viIiiI}8xyˆt‚|¬‰wysfy•GFˆ)wyŠtveAi†)wysutv{‹ˆM|§w8t&ŠAi ·ywA|¬tv|§)w¶|§w萙|~s tveyi)H”ˆM{‚xysveQªIHxyeAwQªµd>xy†‹¾‰i{ HJH”d  svi tM©&)…Qt‚|¬hˆ)¦¬|¬t„r"†)wyŠQ|¬tv|§)wXs² Âss„t‹ˆt‚iIŠè|§wètveAi |¬w8tv{‚QŠQxy† tv|§)wZ•£¿i£&|§¦§¦ sviieA£ºhxy†‹eè“X)tve …Gi{‚¤ˆ)sv|¬¤‰i±ˆMwyŠ xXsui©¸xA¦Z|~s tveAi†)wy†i…Qt&M©tveyiˆ)Š­ƒ„‰|¬w8t&…A{‚)“A¦§ih¶² ,L  .MONQP »   K PRP Á,LTS ´µwtvey|§s5sviI†Kt‚|¬‰wZ•M£ i ŠAiIs‚† {‚|¬“Gi/t‚eA|~s©¸){‚h,ˆ)¦¬|~svhŒˆM…A…A¦§|§iIŠtvb l…Qtv|§hˆM¦yp\eyˆ)…Xif’liIsv|¬š‰w©¸‰{5ˆl…yˆ){ut‚|§†xA¦§ˆ){ˆ)i{‚QŠQr\wyˆMh,|~†s …y{v‰“A¦¬iIh¶² UWVYX

Z@[8\%]^_^=`Yabt }

Iwy€p?{w=€ | I}~>t‚nƒ.„:… †=‡… }ˆ ‰

4Š 

ãXä\ÅKÒ Ì ÅÆ`Ç Ì ÓKÓvÍÓ΋ÍMÆyѳÓvÍ Ö Ž

£&eyi{‚i tu • t } ˆ)wyŠ t ˆM{‚i† ‰wys„t‹ˆMw8t‚s&tveyˆMtl…A{‚iIs‚† {‚|¬“Gitveyi£ i|§š)e8t‚s&)©5tveA{‚iisvxA“Qªµ† {‚|¬tvi{‚|~ˆA•A{‚i¦~ˆt‚iIŠtv¦¬|¬©½tI•¼ŠQ{‹ˆMš ˆ)wyŠ¥sv)wy|§†“X\‰h ih,|§s‚sv|¬‰wZ²dfeAiIr¥ˆM{‚iiI¤­ˆ)¦¬xXˆtvicŠ®“\r¥sv)¦§¤\|¬wAštveAi‘‰’ÄCxA¦¬iI{”iI}8xyˆMtv|§)wys™ˆM{‚)xAwXŠOtveyiDš‰i‰hDitv{‚r sv…GiI† |¬·yicŠè“8r"tveyisueXˆM…Gi±…yˆM{‹ˆMh,itvi{c² C{‚ˆ‰†Ktv|~†ˆ)¦¬¦§r)•tveyi±)“ysvi{‚¤ˆt‚|¬‰w3“G­À q r |§sˆO…yˆM{vtM©ft‚eAi† ‰h,…AxQt‚ˆMtv|§)wyˆ)¦ ŠA)hˆM|§w¶…A¦§ˆ‰† icŠ “Gi¦§£ tveAiˆM|§{v…A¦~ˆMwyi)²C´«t‹s-xA…A…Gi{l“X‰xAwyŠAˆ){vr |§sˆ…A¦~ˆMwAi†¦¬8sui”“Gi¦§£ t‚eAiˆ)|¬{‹† {‹ˆ©½tc² ©€† )xy{‚svi)•A|§w …y{vicsuiIwy† iM©¿sveAQ†‹¾Qs•XtveAiD¤ˆM¦§xAiDM©Ct‚eAi,pGo ’™Ä tvi{‚h |~sÂ|§wQ·ywA|¬tvi,|§w®t‚eAiD†)w8tv|§w\xA)xXsl†Iˆ)svi)•`“AxQt™·ywA|¬tvi ;ˆ)wyŠ®h,iIsve ŠAi…GiwyŠQiIw8t  |§w t‚eAiŠQ|~sv†{vitvi”†ˆ‰sui‰²€dfeAi† ‰w‰t‚{v‰¦`…yˆM{‹ˆMh,i t‚i{  £&|§¦¬¦Z“Gi”{‚i¦~ˆtvicбt‚DtveAisveyˆM…Gi  M©tveAiˆ)|¬{‹† {‹ˆ©½tc²                                                                          

Sonic Boom Downwards Emission Control Box (near field)

Tail Wave

Bow Wave

Sonic Boom N−Wave (far field) p

Ê|§š)xA{‚i®‰°,dfeyi sv)wA|~†“G8‰h¶²3p\¾‰i t‚†‹e )©-wAiIˆ){ˆ)wyŠè©;ˆM{·yi¦~Š sveAQ†‹¾"£fˆ­¤)i…Xˆtut‚i{‚wys)©Âˆ¥svxA…Gi{‹su‰wA|§†ˆM|§{‚†{‚ˆM©½tI² p o ’™Ä s-† ‰w8tv{‚)¦`“X­À q8rf² G ¸c¡MŸ ú ž  ú ÿ >Ÿ þ œ  ú

dfeAi¿ŠQ|§s‚† {‚i t‚i¿É Ê5’ h,QŠQi¦‰xysuicsZˆMw”xy…8£&|§wyŠĀxA¦§i{su‰¦¬¤‰i{Zˆ)…A…A¦§r8|§wAšft‚ÂˆMwxAwysutv{‚xy† tvxA{‚iIŠltvicsvsvi¦~ˆt‚|¬‰w”† )h,…G‰sviIŠ )©¼tvitv{‹ˆMeAicŠQ{‚ˆy²dfeyi™sveyˆM…GiÂ|~sf†‹eyˆMwyš)iIŠ,“\r,h,¤\|¬wAštveAi™wAQŠQics/‰wtveAi™“G)xAwyŠyˆM{‚rD)©¼tveAi™h,iIsve ˆ)¦¬‰wAšwA‰{vhˆM¦~s/tv t‚eyˆt™“X‰xAwyŠAˆ){vr‰²ÂdfeAiI|¬{™ŠQ|§sv…A¦~ˆ)†ih,iw8tl|~sÂt‚ˆM¾‰iwO|§w8tv ˆ‰††)xAw8t“\r¶ˆtv{‹ˆMwXsu…A|§{‹ˆtv|§)wO†)wyŠQ|¬tv|§)w¥|¬w®‰{‚ŠAi{&t‚ ˆ­¤)‰|§Š †‰sutv¦§r{viIh,iIsveA|¬wyš‰sIF² *š‰{‚ˆ‰ŠQ|§iw8t )©>t‚eAiŠQ|~sv†{vitvi™©¸xAwy† tv|§)wyˆ)¦¼|§s&†)h,…AxQt‚iIŠ£&|¨t‚etveAi”eAiI¦¬…M©5ˆMwOˆ)Šƒ„)|§w‰t-svrQs„t‚ih¶• “yxA|¬¦¬t-xysv|¬wyšt‚eAbi -xAtv)hˆMtv|~†”’Â|¬ÃGiI{viIw8tv|~ˆtv|§)w tv\‰ ¦ 

   `.•  –‰*Ž ³² 



UWVU

>e_\n[8dCczdC^_cai a h]cGa dfe

Om \ŠAi{‚w ·XwA|¨t‚i ªµŠQ|§hDiIwysv|¬‰wyˆM¦)…Qt‚|¬h,|§¯Iˆt‚|¬‰w¹h,i tvey\Šys®{‚i¦§r8|§wA𠁉wŒˆ)Š­ƒ„‰|¬w8t‚s"ˆM{‚i |§s‚svxAiIŠ ©¸{‚)h tveAi p\iI}8xAiIw‰t‚|§ˆ)¦  xXˆ)ŠQ{‹ˆt‚|§†C{v‰š){‹ˆMh,h,|¬wAšh,itveAQŠQ‰¦¬‰š)r)² Œ…G)…Axy¦§ˆ){l…A{‚Mt‚Mt„r\…Gi|§sÂtveyio r\{‚ŠQª"lh,Mƒ„‰¾\xAwOˆM¦§š)){‚|¬tveAh¶•Xsvii ‘8Ž*«² dfey|§sˆM¦§š)){‚|¬tveAh#|¬w |¬t‚s“yˆ)sv|§†D©¸‰{vh#©¸){ˆ)s‚suxAh,icstveyˆMtt‚eAi±{‚iIsv)¦§xQtv|§)wèM©ft‚eAiŠQ|¬ÃGiI{viIw8t¦¬|§wAicˆM{‚|§s‚ˆt‚|¬‰wys”M©-sut‚ˆMtvi svrQs„t‚ihs ;jÂi£ftv‰w |¬tviI{‚ˆMtv|§)w M©sut‚ˆMtviˆMwyŠsv)¦§xQtv|§)w M©ˆ‰Š­ƒ„)|§w8t ¿ˆM{‚ilwyMt&i ÀQ…yˆ)wysu|§¤)i‰² £¿iI¤)iI{I•‰tveA|~s&ˆ)s‚svxAh,…Qtv|§)w |~s,wA)t,¤ˆ)¦¬|~Š |§w l…Qtv|§hˆM¦-p\eyˆ)…Xi®’liIsv|¬š‰wZ² dfey|§sD©;ˆ‰†Kteyˆ)s,¦§iIˆ‰Š sv)h,iOˆMxQt‚eA){‹stv ˆtvt‚ˆ)†‹¾èt‚eAi¶…A{‚)“A¦§ih xXsu|§wAš þ  ú)ü   þ ž ª³‰{C…A{‚)š){‚iIs‚sv|¬¤‰i)•M‰{/su|§hxy¦¨t‹ˆMwAiI)xysuªˆM¦§š)‰{v|¬tveyh,s ;svii-©¸‰{/iÀQˆ)h,…A¦¬i  8— ³• §  K•\£&eA|~†‹e±ˆM{‚i&“yˆ)sviIŠ,|§wt‚eAi ©¸‰¦¬¦§£&|§wAšt„£ ,…X‰|¬w8t‹s° i HJH”dÁsvrQs„t‚ih ˆt tveAièsvˆ)h,iOtv|§h,i"|§w iIˆ‰†‹e )…Atv|§hD|§¯IˆMtv|§)w " Ê|¬{‹s„tc•lˆ‰ŠQ¤ˆMwy†iOtveAi"tvey{viIi¥iI}8xyˆMtv|§)wysM©”t‚eA@ |¨t‚i{‹ˆt‚|¬‰wZ²

Ì ÓK×cØ Ì Ï Ë ÍÏQÓKѳæ  ÛlÍæ  àIìÏ Ì   ÍIÍ KÏ\Õ F

j-iÀ8tc•M|§wys„t‚iIˆ‰ŠM©Xxysu|§wA𙦬|§wAicˆM{š)iwyi{‚|§†¿su‰¦¬¤‰i{‹s•cxysvi ŠQ|~s‚† |§…A¦¬|§wAiª«sv…Xic† |¬·X†¿|¨t‚i{‹ˆtv|§¤)i‰•hˆ­r8“Gi/wy)wA¦§|¬wyiIˆM{5sv)¦§¤)i{‹s ½©¸‰{-i ÀAˆMh,…A¦§i)•Q…XsuiIxyŠQ,xAwysutvicˆ)ŠQrsv)¦§¤)iI{‚s¿©¸‰{-ʦ§xA|~ŠOm¶iI†‹eyˆ)wA|~†s K²    ú)ü   þ ž ü Ÿ >Iž ü ú  þ ß   þ Ÿ ½ž 

dfeyis„t‹ˆtvi”ic}8xyˆt‚|¬‰w¶sv)¦§¤)iI{ftv“Gi”|¬w8t‚iš){‹ˆt‚iIŠ|§w t‚eAi…A{‚)…G‰sviIŠ )…Atv|§hD|§¯iI{-|~s-ˆ)s‚suxyhDicŠ tv“GiˆMw¶|¬tvi{‹ˆt‚|¬¤‰i™·AÀQiIŠ …G)|§w8t&)wAi‰²€´«t-{viI¦¬|§iIsf‰wˆÒyÕ Ì ÏyÇ)ÍMå Ì Ñ«ÍMƱs„t‚i…¶)w¤ˆM{‚|§ˆ)“A¦¬i  ° "



?   $

Suggest Documents