Optimization loops for shape and error control: extended lecture notes

0 downloads 0 Views 1MB Size Report
Dec 15, 2004 - I INTR O DU С TI O N. Multidisciplinary modelling is an important ingredient for the solution of an engineering problems of utmost importance in ...
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  ° "



?   $