谱元方法求解波动方程及影响其数值精度的相关因素.pdfVIP

第42卷 第1期 西 安 交 通 大 学 学 报 Vol.42 Xo1 2008年1月 JOURNALOFVANJIAOTONGUNIVERSITY Jan.2008 谱元方法求解波动方程及影响其数值精度的相关因素 朱昌允，秦国良，徐忠 (西安交通大学能源与动力工程学院，710049，西安) 摘要:为探讨波动方程的高精度数值模拟，采用Chebyshev谱元方法结合隐式Newmark时间积 分方法求解波动方程.求解一个具体算例验证了数值方法的可行性，讨论了时间步长、Newmark 因子以及计算区域的网格剖分方式对数值精度的影响.结果表明:和差分法相比，谱元方法求解波 动方程具有所用网格节点少，数值精度高的特点;数值误差随时间步长减小而减小;在满足稳定性 要求的前提下，数值误差随着Newmark因子的减小而减小;当总网格节点数相同时，不同的网格 剖分方式所得数值误差不同.所述方法和结论可用于模拟声波在空气中的传播. 关键词:波动方程;谱元方法;时间积分方法;气动声学 中图分类号:0353 文献标志码:A 文章编号:0253-987X(2008)01-0056-04 ImplicitSpectralElementMethodforWaveEquationandSome FactorsInfluencingNumericalAccuracy ZHUChangyun,QINGuoliang,XUZhong (SchoolofEnergyandPowerEngineering,XianJiaotongUniversity,Xian710049,China) Abstract;Toinvestigatethenumericalschemewithhighorderofaccuracyforthesimulationof waveequations,ChebyshevspectralelementmethodcombinedwithimplicitNewmarktimeinte- gralmethodisadoptedforsimulatingwaveequations.Thensomefactorsaffectingthenumerical accuracyarediscussedindetail，suchastimesteph,Newmarkfactorandthesubdivisionstyle forcomputationaldomain.Theconclusionsindicatethatthespectralelementmethodhashigher ordernumericalaccuracythandifferencemethodforsimulationofwaveequations，smalltimestep induceslessnumericalerror,smallerNewmarkfactorinducessmallernumericalerror,anddiffer- entmeshstyleinducesdifferentnumericalerrorwhilethegeneralgridmembergetsequaltoeach other.Theproposedmethodenablestosimulateacousticwavepropagationinair. Keywords;waveequation;spectralelementmethod;timeintegralmethod;aeroacoustics 谱元方法是结合谱方法和有限元法各自特点的 流体和声的相互作用等.目前面临的主要问题[3]之 一种数值方法，它既具有谱方法的高精度和收敛特 一是传统数值方法的精度不能完全满足计算气动声 性，又具有较好的几何区域适应性，20