文档介绍：万方数据
型旦耍。Euler型有限体积法(UdX)IF景荒ā:雗ds((^)L()钱旭,东晨,宋松争Euler=tT{UdX()q0=1q=2其中,lAAknAk^AAkqGaussA1986Harten(nonoscillatory)EngquistOsherChakravathy11(essentially猳袷降姆桨负头椒āS捎贓格式在重构时,在多个备选模板中仅选择了一个最优模板而摒弃了其它模板,这样就把很多有价值的信息都丢弃ENO本质无振荡旧。ENOWENO)WENOENO格式利用各个备选模板的凸组合的方式重构,每个模板的权值也就是组合系数的选取依赖于该模WENOENOENOENOWENOENOWENO用弹簧方法在监测函数的指导下拉伸网格节点,其原理简单,操作直观,因此该方法在移动网格领域是一种常用的做法。EulerF=(g)X=(y)(1)AGreenShuOsherTVD狵方法作者简介:钱旭,男,重庆人,硕士生,研究方向:偏微分方程数值解及其应用。314年月湘潭师范学院学报匀豢蒲ОEulerWENO要思想是先对每一个三角形单元构造插值多项式,而在计算交界面的流通量采用两点高斯积分公式以保证格式的整体精TVD狵方法。最后用该格式对一些典型算例进行了数值模拟,结果表明该方法计算速度快,对间断有很好的分辨能力。EulerWENO中图分类号:A文章编号:———锓揽萍即笱Ю硌г海铣ど收稿期:——基金项目:国家重点基础研究发展规划钅—;国家自然科学基金项目Vol3lNo4Dec2009a^.
万方数据
觥季纔,‰口吲警,学,.一疩.,冢一÷玩手玩帧鱰·￡駻’,。构造式捎玫脑际跫!篢:瑈以。p2()=2)加权插值多项式的构造sp2())基于弹簧技术的移动网格个二次多项式名,作光滑探测器礁,构造茗,,,保个相关单元组成的连通区域为酝,把这个连通P()=o+a1+a2yk=l23iP()3(4)3II6P(y)k=l26k一蝣一听D=+一西4W4控制单元缤所示,它周围的龅ピF颍琧,,⋯,记,上的网格平均值为致。要在I瞎乖煲桓鱿咝远嘞钍剑蛑恍柙赾:,琧腥稳个即可。(5)3在与合嗔诘个单元琧陀隿。相邻的#琧腥稳个,再加上烧个单元上的以来确定二区域中的所有共边的三角形单元上的网格平均值之差的平方和相加,记作胚。。令作如下加权二次多项式厶哟弹簧振子方法的基本思想是将网格结构中所有的边均看作倔强系数相等但长度不等的可伸展弹簧,监测函数看作是在流场布置安放这些不同长度弹簧的准则。在监测函数值小的地方弹簧长度自然小,反之亦然。在移动网格中,我们最关心的就是监测函数石,,它指导网格点在运算中移动,引导整个计算区域网格的疏密程度。设节点琘M裰幸桓鼋诘悖窠峁怪杏孝鬚条边与其相连,设其中一条边为PQP())=())+k其中,5删笄肯凳尬>槿ㄖ担粃,为两节点间的距离,口为矢量的单位矢,万为网格归一化因叼∞=算汁三角形网格单元2从以出扯褚41●
万方数据
()=()=(c)+Ce-"I则执行此次操作,否则将∞变为等,执行上述操作直到符合上述限制条件。在数值计算中,监测函数石,数值试验及结果子,一般取;∑为网格中所有边长