五阶WENO有限差分法在线性双曲守恒律方程中的应用.pdfVIP

第26卷第1期 广西科技大学学报 V01．26No．1 2015年3月 JOURNAL0FGUANGXIUNIVERSl7IY0FSCIENCEAND TECHNOLOGY Mar．2015 文章编号2095．7335(2015)01-0090．06 五阶WENO有限差分法在 线性双曲守恒律方程中的应用 汤淑芳，林贤坤+，覃柏英，秦文东 (广西科技大学 汽车与交通学院，广西柳州 545006) 摘 性双曲守恒律方程中的应用．经过经典数值算例的验证。结果表明五阶WENO有限差分法可实现线性双曲守恒律方 程高精度、高分辨率和本质无振荡的求解，也可实现流体力学中运动界面高精度、高分辨率的追踪． 关键词：WEN0；数值计算；双曲守恒律方程；界面追踪 中图分类号：0241；0354 文献标志码：A 0引言 在科学与工程问题中，普遍存在一类非常重要的偏微分方程，即双曲守恒律方程．空气动力学、流体动 力学、爆炸力学等许多力学问题的求解都与之相关．即使初始条件充分光滑，双曲型守恒律方程的解仍然 可能出现间断．因此，为了能高效地捕捉到间断，并避免间断附近的数值振荡。高阶精度、高分辨率的差分 方法，一直以来都是计算流体力学领域致力研究的热点问题． Galerkin 自1973年Reed和HiU提出间断Galerkin方法(Discontinuous V耐ation 方程求解的理论分析和应用得到迅速发展．而1983年H献en提出高分辨率、无振荡的TVD(Total andminimum Bounded)‘21，NOC (Non—OscillatoryCentral)㈨NND Finite Dissipation)[4]ENO(EssentialN0n—Oscillatory)[5-71，以及SV(Spectral Essential Non—Oscillatory)‘93等差分格式已相继提出和发展，并在科学与工程领域得到应用． 个差值模板的非线性凸组合重构，在相同模板上，WENO格式在光滑区域达到高阶精度，在间断区域相当 于ENO格式，从而使其在格式上更加有效、通量上更加光滑、并且有更加稳定的收敛解．呈现出本质无振荡 的特性．由于WENO格式具有高精度，高分辨率的良好性质，因此特别适合激波湍流问题，航空声学问题， 界面追踪问题等间断问题的处理阻也]，双曲守恒律方程的求解． 本文首先探讨五阶WENO有限差分法应用于求解双曲型偏微分方程的思想与方法，然后通过一维、 二维和三维的相关数值算例验证该方法可实现线性双曲守恒律方程高精度、高分辨率和本质无振荡的求 解，并说明该方法可实现流体力学中运动界面高精度、高分辨率的追踪． 1五阶WENO有限差分法 收稿日期：2014—09—22 基金项目：国家自然科学基金项目11272057)资助． 4通信作者：林贤坤，副教授，博士，研究方向：计算流体力学，E—mail：1in)【l【0209@163．com． 万方数据 第l期 汤淑芳等：五阶WEN0有限差分法在线性双曲守恒律方程中的应用 91 戈“。尼的左右状态值t彪Z尼，可利用五阶wENo有限差分法计算+，其具体思想与方法如下： 函数／在子模板上关于边界菇“。尼的左状态值为： 簏’=争知-+告，蔗’一知t+》扛-，虎’=扣缸t一每 (1) 由此可分别构成一个凸组合计算边界戈许抛五阶精度的左状态值只‰即 庀陀：∑∞崩：∞碟’Ⅷ。患’Ⅷ藏’ (2) 其中，姊是