Perfectly matched layer as an absorbing boundary condition for the linearized Euler equations in
上传者:秦龙|上传时间:2015-05-07|密次下载
Perfectly matched layer as an absorbing boundary condition for the linearized Euler equations in
JOURNALOFCOMPUTATIONALPHYSICS
ARTICLENO.144,213–234(1998)CP985997
PerfectlyMatchedLayerasanAbsorbing
BoundaryConditionfortheLinearizedEuler
EquationsinOpenandDuctedDomains
ChristopherK.W.Tam,?LaurentAuriault,?andFrancescoCambuli?
?DepartmentofMathematics,FloridaStateUniversity,Tallahassee,Florida32306-4510;and?Dipartimento
diIngegneriaMeccanica,UniversitadegliStudidiCagliari,Piazzad’Armi,09123,Cagliari,Italy
E-mail:tam@math.fsu.edu
ReceivedOctober30,1997;revisedApril6,1998
Recently,perfectlymatchedlayer(PML)asanabsorbingboundarycondition
hasfoundwidespreadapplications.Theideawas?rstintroducedbyBerengerfor
electromagneticwavescomputations.Inthispaper,itisshownthatthePMLequations
forthelinearizedEulerequationssupportunstablesolutionswhenthemean?owhasa
componentnormaltothelayer.Tosuppresssuchunstablesolutionssoastorenderthe
PMLconceptusefulforthisclassofproblems,itisproposedthatarti?cialselective
dampingtermsbeaddedtothediscretizedPMLequations.Itisdemonstratedthat
withaproperchoiceofarti?cialmeshReynoldsnumber,thePMLequationscan
bemadestable.Numericalexamplesareprovidedtoillustratethatthestabilized
PMLperformswellasanabsorbingboundarycondition.Inaductedenvironment,
thewavemodesaredispersive.Itwillbeshownthatinthepresenceofamean
?owthegroupvelocityandphasevelocityofthesemodescanhaveoppositesigns.
ThisresultsinabandoftransmittedwavesinthePMLtobespatiallyamplifying
insteadofevanescent.Thusinacon?nedenvironment,PMLmaynotbesuitableas
c1998AcademicPressanabsorbingboundaryconditionunlessthereisnomean?ow.??
1.INTRODUCTION
Recently,Berenger[1,2]succeededinformulatinganabsorbingboundaryconditionforcomputationalelectromagneticsthathastheunusualcharacteristicthatwhenanoutgoingdisturbanceimpingesontheinterfacebetweenthecomputationdomainandtheabsorbinglayersurroundingit,nowaveisre?ectedbackintothecomputationdomain.Inotherwords,alltheoutgoingdisturbancesaretransmittedintotheabsorbinglayerwheretheyaredampedout.Suchalayerhascometobeknownasaperfectlymatchedlayer(PML).
Sinceitsinitialdevelopment,PMLhasfoundwidespreadapplicationsinelasticwavepropagation[3],computationalaeroacoustics,andmanyotherareas.Hu[4]wasthe?rstto
213
0021-9991/98$25.00
c1998byAcademicPressCopyright??
Allrightsofreproductioninanyformreserved.
214TAM,AURIAULT,ANDCAMBULI
applyPMLtoaeroacousticsproblemsgovernedbythelinearizedEulerequations,linearizedoverauniformmean?ow.Hehassinceextendedhisworktononuniformmean?owandforthefullynonlinearEulerequations[5].FurtherapplicationsofPMLtoacousticsproblemsincludingwavemodesinductscanbefoundinthemostrecentworksofHuandco-workers[6,7].Inthesereferences,examplesareprovidedthatindicatethathighqualitynumericalsolutionscouldbefoundwithPMLusedasradiationorout?owboundaryconditions.
Inopenunboundeddomains,acousticwavesarenondispersiveandpropagatewiththespeedofsoundrelativetothelocalmean?ow.Insideaduct,thesituationiscompletelydifferent.Acousticwavesarerepeatedlyre?ectedbackbythecon?ningwalls.Forductswithparallelwalls,thecontinuousre?ectionoftheacousticwavesbythewallleadstotheformationofcoherentwavepatternscalledductmodes[8,9].Unliketheopendomain,ductmodesaredispersivewithphaseandgroupvelocitiesvarywithaxialwavenumber.Becauseofthedispersivenatureoftheductmodesmanyradiationboundaryconditionsthatworkwellinopendomainsareknowntobeinappropriateforductedenvironments.Forthisreason,Tam[10]inarecentreviewonnumericalboundaryconditionsforcomputationalaeroacousticssuggestedthattheboundaryconditionforaductedenvironmentberegardedasacategoryofitsown.
Therearethreeprimaryobjectivesinthiswork.First,weintendtoshowthatinthepresenceofamean?ownormaltoaPML,thestandardPMLequationsofthelinearizedEulerequationssupportunstablesolutions.EarlierTam[10]pointedoutthatthePMLequationswithmean?owhaveunstablesolutions.However,hedidnotshowthattheexistenceofinstabilitiesisduetothemean?owcomponentnormaltothelayer.Theoriginandcharacteristicsoftheseinstabilitiesareinvestigatedandanalyzed.Itisinterestingtomentionthatinhisearliestwork,Hu[4]reportedthathiscomputationencounterednumericalinstability.Butbyapplyingnumerical?ltering,hewasabletoobtainstablesolutions.Inlightofour?nding,webelievethatwhatHuencounteredwasnotinstabilityofhisnumericalschemebutthathisnumericalsolutioninadvertentlyexcitedtheintrinsicunstablesolutionofthePMLequations.NotdirectlyrelatedtotheinstabilityofthePMLequations,AbarbanelandGottlieb[11]recentlyanalyzedtheelectromagneticPMLequations.Theyconcludedthattheequationsareonlyweaklywell-posed.
Second,wewillshowthattheinstabilityisnotverystrong,namely,thegrowthratesaresmall.Alsotheinstabilitiesarecon?nedprimarilytoshortwaves.Itis,therefore,possibletosuppresstheinstabilitiesbytheadditionofarti?cialselectivedampingterms[12]tothediscretizedPMLequations.Itisimportanttopointoutthatarti?cialselectivedampingeliminatesmainlytheshortwavesandhasnegligibleeffectonthelongorthephysicalwaves.ThustheadditionofthesedampingtermsdoesnoteffecttheperfectlymatchedconditionsofthePML.
Third,wewillshowthataperfectlymatchedlayermaynotbesuitableasanabsorbingboundaryconditionforwavesinaducted?owenvironment.Themajordifferencebetweenacousticwavesinanopendomainandacousticwavesinsideaductisthatinanunboundedregionacousticwavesarenondispersivewhereasductmodesaredispersive.Itwillbeshownthatinthepresenceofamean?owthegroupandphasevelocityoftheductmodescanhaveoppositesigns.Becauseofthis,abandoftransmittedwaveswillactuallygrowspatiallyinsteadofbeingdampedinthePML.Inotherwords,thePMLequationsdonotdampthesewavemodesasabsorbingboundaryconditionoughttodo.Theexceptioniswhenthereisnomean?owintheduct.Inthisspecialcase,allthetransmittedwavesarespatiallydamped.
PMLFORLINEARIZEDEULEREQUATIONS215
InSection2,theuseofPMLforopendomainproblemsisdiscussed.ThestabilityofthePMLgoverningequationsisinvestigated.Itwillbeshownthattheadditionofdamp-ingtermstoformthePMLequationscanactuallycausethevorticityandacousticwavemodestobecomeunstable.ThesplittingofthevariablesinformulatingthePMLequationsleadstoahigherordersystemofequations.Thishighersystemsupportsextrasolutions.Theseextraorspurioussolutionsarefoundtobecomeunstablewhenthedampingcoef?cientislarge.NumericalexamplesareprovidedtoillustratethespreadoftheunstablesolutionfromthePMLbackintotheinteriorofthecomputationdomain.
InSection3,theeffectoftheadditionofarti?cialselectivedampingtermstothedis-cretizedPMLequationsisinvestigated.ItisshownthatwithanappropriatechoiceofmeshReynoldsnumber,theunstablesolutionsofthePMLequationscanbesuppressed.Numericalexamplesaregiventodemonstratetheeffectivenessofthemodi?edPMLasaradiation/out?owboundarycondition.
Section4dealswiththetheoryandapplicationofPMLtoductedinternal?owproblems.AneigenvalueanalysisiscarriedouttoshowtheexistenceofabandoffrequencyforwhichthePMLexertsnodampingontheacousticductmodes.ThesewavemodesactuallywouldgrowinamplitudeastheypropagatethroughthePML.Numericalresultsareprovidedtoillustratetheexistenceofthiskindofamplifyingductedacousticmodes.
2.OPENDOMAINPROBLEMS
LetusconsidertheuseofPMLasabsorbingboundaryconditionforthesolutionofthelinearizedEulerequations(linearizedoverauniformmean?ow)inatwo-dimensionalopendomainasshowninFig.1.Wewilluse??x=??y(themeshsize)asthelengthscale,a0x2(thesoundspeed)asthevelocityscale,
内容需要下载文档才能查看??asthetimescale,andρ0a0(whereρ0isthea0
FIG.1.TwodimensionalcomputationdomainwithPerfectlyMatchedLayersasboundaries.
216TAM,AURIAULT,ANDCAMBULI
meandensity)asthepressurescale.ThedimensionlessgoverningequationsinthePMLareformedbysplittingthelinearizedEulerequationsaccordingtothespatialderivatives.Anabsorptiontermisaddedtoeachoftheequationswithspatialderivativeinthedirectionnormaltothelayer.Forexample,forthePMLontherightboundaryofFig.1(notatthecorners)thegoverningequationsare[4]
???u1+σu1+Mx(u1+u2)+(p1+p2)=0?u2?+My(u1+u2)=0?v1?+σv1+Mx(v1+v2)=0?v2??+My(v1+v2)+(p1+p2)=0???p1+σp1+Mx(p1+p2)+(u1+u2)=0?p2??+My(p1+p2)+(v1+v2)=0,whereMxandMyarethemean?owMachnumbersinthexandydirections.σistheabsorptioncoef?cient.
Supposewelookforsolutionswith(x,y,t)dependenceintheformexp[i(αx+βy?ωt)].Itiseasyto?ndfrom(1)thatthedispersionrelationsofthePMLequationsare
??βMyαMx?1???2α2β2??=0βMyαMx?=0.(2)(3)(1)1?
Inthelimitσ→0,(2)and(3)becomethewell-knowndispersionrelationsoftheacousticandthevorticitywavesofthelinearizedEulerequations.
2.1.MeanFlowParalleltoPML
Dispersionrelations(2)and(3)behaveverydifferentlydependingonwhetherthereisanymean?ownormaltothePML.Whenthemean?owisparalleltothelayer,i.e.,Mx=0,thesolutionsarestable.Thisiseasytoseefrom(3)forthevorticitywave.Physically,ifthemean?owisparalleltothePML,thevorticitywavesinthecomputationdomain,beingconvectedbythemean?ow,cannotenterthelayerandhencewouldnotleadtounstablesolution.
ToshowthatforMx=0allthesolutionsof(2)arestable,asimplemappingwillsuf?ce.Rewrite(2)intheform
α2ω2=β2.F≡(ω?βMy)?22(4)
Figure2showstheimageoftheupper-halfω-planeintheFplane.Theupper-halfω-planeismappedintotheentireFplaneexceptfortheslitADC.Butsinceβ2isrealandpositive,
PMLFORLINEARIZEDEULEREQUATIONS
内容需要下载文档才能查看217
FIG.2.Theimageoftheupperhalfω-planeintheF-plane.
forsubsonicmean?owthepointβ2liesoutsidetheimage.Thusnovalueofωintheupper-halfω-planewouldsatisfyEq.(2)indicatingthatthereisnounstablesolution.
2.2.UnstableSolutionsofthePMLEquations
ForMx=0,thePMLequationssupportunstablesolutions.Itistobenotedthat,unliketheoriginaldispersionrelationoftheacousticwaves,Eq.(2)isaquadricequationinω.Ithastwoextrarootsinadditiontothetwomodi?edacousticmodes.Forsmallσ,thetwospuriousrootsaredampedbutoneofthemodi?edacousticrootsisunstable.Forlargerσ,numericalsolutionsindicatethatoneofthespuriousrootsbecomesunstable.Inanycase,theequationsplittingprocedureandtheadditionofanabsorptionterm,botharevitaltothesuppressionofre?ectionsattheinterfacebetweenthecomputationdomainandthePML,inadvertently,leadtoinstabilities.
Forsmallσ,therootsof(2)and(3)canbefoundbyperturbation.Let
(a)(a)(a)+σω1+σ2ω2+···ω(a)=ω0
(v)(v)(v)+σω1+σ2ω2+···,ω(v)=ω0(5)(6)
wheretherootsof(2)and(3)aredesignatedbyasuperscripta(foracousticwaves)andv
下载文档
热门试卷
- 2016年四川省内江市中考化学试卷
- 广西钦州市高新区2017届高三11月月考政治试卷
- 浙江省湖州市2016-2017学年高一上学期期中考试政治试卷
- 浙江省湖州市2016-2017学年高二上学期期中考试政治试卷
- 辽宁省铁岭市协作体2017届高三上学期第三次联考政治试卷
- 广西钦州市钦州港区2016-2017学年高二11月月考政治试卷
- 广西钦州市钦州港区2017届高三11月月考政治试卷
- 广西钦州市钦州港区2016-2017学年高一11月月考政治试卷
- 广西钦州市高新区2016-2017学年高二11月月考政治试卷
- 广西钦州市高新区2016-2017学年高一11月月考政治试卷
- 山东省滨州市三校2017届第一学期阶段测试初三英语试题
- 四川省成都七中2017届高三一诊模拟考试文科综合试卷
- 2017届普通高等学校招生全国统一考试模拟试题(附答案)
- 重庆市永川中学高2017级上期12月月考语文试题
- 江西宜春三中2017届高三第一学期第二次月考文科综合试题
- 内蒙古赤峰二中2017届高三上学期第三次月考英语试题
- 2017年六年级(上)数学期末考试卷
- 2017人教版小学英语三年级上期末笔试题
- 江苏省常州西藏民族中学2016-2017学年九年级思想品德第一学期第二次阶段测试试卷
- 重庆市九龙坡区七校2016-2017学年上期八年级素质测查(二)语文学科试题卷
- 江苏省无锡市钱桥中学2016年12月八年级语文阶段性测试卷
- 江苏省无锡市钱桥中学2016-2017学年七年级英语12月阶段检测试卷
- 山东省邹城市第八中学2016-2017学年八年级12月物理第4章试题(无答案)
- 【人教版】河北省2015-2016学年度九年级上期末语文试题卷(附答案)
- 四川省简阳市阳安中学2016年12月高二月考英语试卷
- 四川省成都龙泉中学高三上学期2016年12月月考试题文科综合能力测试
- 安徽省滁州中学2016—2017学年度第一学期12月月考高三英语试卷
- 山东省武城县第二中学2016.12高一年级上学期第二次月考历史试题(必修一第四、五单元)
- 福建省四地六校联考2016-2017学年上学期第三次月考高三化学试卷
- 甘肃省武威第二十三中学2016—2017学年度八年级第一学期12月月考生物试卷
网友关注
- 苏州园林
- 如何对待逆反孩子
- 驴小弟变石头
- 陕西艺术幼儿园音乐活动教研心得(3.14)
- 2016年四川省内江市中考化学试卷
- 端午节放假安全教育材料
- 学生服使用单位履行质量义务情况专项检查记录表(幼儿园)
- 端午节安全教育
- 洋县理光复印土管局大门北:幼儿园校安全工作日志
- 广东省深圳市西丽幼儿园分园装修工程可行性研究报告-广州中撰咨询
- 留守儿童之家读书活动安排意见
- 登鲁的故事
- 六年级下册6《半截蜡烛》优质教案
- 2017人教版七下语文一单元检测题
- 小学古诗词归类整理
- 2010-2011学年北京市海淀区高三第二学期期末练习(语文)答案
- 广东省佛山市均安镇星槎幼儿园工程可行性研究报告-广州中撰咨询
- 2011-2012学年吉林省长春外国语学校初二上学期第二次月考语文试卷答案
- 给父母的一份廉洁家书
- 兰德公司“T项目”研究及其对我国学前教育的启示
- 2016----2017年度小班名画欣赏
- 王春亮民间推拿传艺文化
- 人教版初中语文七年级(下)第7课:《土地的誓言》课件
- 黄道婆
- 关于幼儿园转让协议
- 语段练习卷
- 父母永不放弃自我成长,是孩子最好的教育!
- 幼儿教育心得
- 小古文诵读100篇[1]
- 母亲节所思
网友关注视频
- 【部编】人教版语文七年级下册《泊秦淮》优质课教学视频+PPT课件+教案,广东省
- 冀教版小学数学二年级下册第二单元《有余数除法的整理与复习》
- 《空中课堂》二年级下册 数学第一单元第1课时
- 沪教版牛津小学英语(深圳用) 四年级下册 Unit 8
- 外研版英语三起6年级下册(14版)Module3 Unit1
- 第4章 幂函数、指数函数和对数函数(下)_六 指数方程和对数方程_4.7 简单的指数方程_第一课时(沪教版高一下册)_T1566237
- 二年级下册数学第二课
- 沪教版牛津小学英语(深圳用) 五年级下册 Unit 12
- 沪教版八年级下册数学练习册21.3(2)分式方程P15
- 苏科版数学八年级下册9.2《中心对称和中心对称图形》
- 精品·同步课程 历史 八年级 上册 第15集 近代科学技术与思想文化
- 北师大版八年级物理下册 第六章 常见的光学仪器(二)探究凸透镜成像的规律
- 冀教版英语五年级下册第二课课程解读
- 冀教版小学数学二年级下册第二单元《租船问题》
- 冀教版小学数学二年级下册第二单元《有余数除法的竖式计算》
- 冀教版英语三年级下册第二课
- 30.3 由不共线三点的坐标确定二次函数_第一课时(市一等奖)(冀教版九年级下册)_T144342
- 沪教版牛津小学英语(深圳用) 四年级下册 Unit 3
- 沪教版八年级下册数学练习册21.3(3)分式方程P17
- 第五单元 民族艺术的瑰宝_15. 多姿多彩的民族服饰_第二课时(市一等奖)(岭南版六年级上册)_T129830
- 沪教版牛津小学英语(深圳用)五年级下册 Unit 1
- 河南省名校课堂七年级下册英语第一课(2020年2月10日)
- 【部编】人教版语文七年级下册《泊秦淮》优质课教学视频+PPT课件+教案,湖北省
- 苏科版数学七年级下册7.2《探索平行线的性质》
- 沪教版牛津小学英语(深圳用) 六年级下册 Unit 7
- 【部编】人教版语文七年级下册《逢入京使》优质课教学视频+PPT课件+教案,安徽省
- 第五单元 民族艺术的瑰宝_16. 形形色色的民族乐器_第一课时(岭南版六年级上册)_T1406126
- 化学九年级下册全册同步 人教版 第25集 生活中常见的盐(二)
- 冀教版小学英语四年级下册Lesson2授课视频
- 【部编】人教版语文七年级下册《泊秦淮》优质课教学视频+PPT课件+教案,天津市
精品推荐
- 2016-2017学年高一语文人教版必修一+模块学业水平检测试题(含答案)
- 广西钦州市高新区2017届高三11月月考政治试卷
- 浙江省湖州市2016-2017学年高一上学期期中考试政治试卷
- 浙江省湖州市2016-2017学年高二上学期期中考试政治试卷
- 辽宁省铁岭市协作体2017届高三上学期第三次联考政治试卷
- 广西钦州市钦州港区2016-2017学年高二11月月考政治试卷
- 广西钦州市钦州港区2017届高三11月月考政治试卷
- 广西钦州市钦州港区2016-2017学年高一11月月考政治试卷
- 广西钦州市高新区2016-2017学年高二11月月考政治试卷
- 广西钦州市高新区2016-2017学年高一11月月考政治试卷
分类导航
- 互联网
- 电脑基础知识
- 计算机软件及应用
- 计算机硬件及网络
- 计算机应用/办公自动化
- .NET
- 数据结构与算法
- Java
- SEO
- C/C++资料
- linux/Unix相关
- 手机开发
- UML理论/建模
- 并行计算/云计算
- 嵌入式开发
- windows相关
- 软件工程
- 管理信息系统
- 开发文档
- 图形图像
- 网络与通信
- 网络信息安全
- 电子支付
- Labview
- matlab
- 网络资源
- Python
- Delphi/Perl
- 评测
- Flash/Flex
- CSS/Script
- 计算机原理
- PHP资料
- 数据挖掘与模式识别
- Web服务
- 数据库
- Visual Basic
- 电子商务
- 服务器
- 搜索引擎优化
- 存储
- 架构
- 行业软件
- 人工智能
- 计算机辅助设计
- 多媒体
- 软件测试
- 计算机硬件与维护
- 网站策划/UE
- 网页设计/UI
- 网吧管理