耿艳辉;秦国良;王阳;贺唯;

• 1:西安交通大学流体机械研究所

摘要：

提出了时空耦合谱元方法,并将其用于带第一类边界条件的非齐次一维、二维、三维波动方程的求解。分别采用四边形、六面体和超六面体作为计算单元,在每个单元内采用Chebyshev多项式的极值点作为Lagrange插值节点,并且探讨了区域剖分方式对计算精度的影响。时空耦合谱元法能够得到精度很高的数值结果,并且其色散随时间推移是稳定的;当总网格节点数相同时,不同的网格剖分方式所得数值误差不同,当空间方向Chebyshev多项式的阶数较高和时间方向Chebyshev多项式的阶数较低时,得到的数值精度较高;在总节点数相同的情况下,与时间全域方式相比,逐时间子区域方式计算所需要的时间更经济,两种方式可以得到相同的精度。结果表明:时空耦合谱元方法使时空方向精度相匹配,可以提高整体精度;空间方向的Chebyshev多项式对数值精度起主要影响作用;时间子区域方式的采用可以扩大问题的计算区域。

关键词：

Abstract: A space-time coupled spectral element method based on Chebyshev polynomials is presented for solving timedependent wave equations.Acoustic propagation problems in 1 + 1,2 + 1,3 + 1 dimensions with the dirichlet boundary conditions are simulated via space-time spectral element method using quadrilateral,hexahedral and tesseractic elements respectively.Space-time coupled spectral element method can obtain high-order precision and the dispersion is stable over time.With the same total number of grid nodes,higher numerical precision is obtained if the higher-order Chebyshev polynomials in space directions and lower-order Chebyshev polynomials in time direction are adopted.When space-time coupled spectral element method is used,time subdomain-by-subdomain approach is more economical than time domain approach.

Keywords:

基金项目: 国家自然科学基金资助项目(51076123)

通讯作者:

Email:

参考文献：

• 1 Patera A T.A spectral element method for fluid dynamics: laminar flow in a channel expansion.Journal of Computational Physics,1984:54(3):468-488
• 2 Dauksher W.Emery A F.Accuracy in modeling the acoustic wave equation with Chebyshev spectral finite elements. Finite elements in analysis and design,1997:26(2):115- 128
• 3 Rong Z J,Xu C J.Numerical approximation of acoustic waves by spectral element methods.Applied numerical mathematics,2008:58(7):999-1016
• 4 Zhu C Y,Qin G L,Zhang J Z.Implicit Chebyshev spectral element method for acoustics wave equations.Finite elements in analysis and design,2011:47(2):184-194
• 5 Zampieri E,Pavarino L F.Approximation of acoustic waves by explicit Newmark's schemes and spectral element methods. Journal of computational and applied mathematics, 2006:185(2):308-325
• 6 Gopalakrishnan S,Doyle J F.Spectral super-elements for wave propagation in structures with local non-uniformities. Computer methods in applied mechanics and engineering, 1995:121(1-4):77-90
• 7张荣欣,秦国良.用切比雪夫谱元方法求解均匀流场中的声传播问题.西安交通大学学报,2009;43(7):120-124
• 8耿艳辉,秦国良.Chebyshev谱元方法求解含吸收边界的二维有均匀来流的声传播问题.西安交通大学学报,2012;46(3): 100-106
• 9 Seriani G.Double-grid Chebyshev spectral elements for acoustic wave modeling.Wave Motion,2004:39(4):351- 360
• 10 French D A,Peterson T E.A continuous space-time finite element method for the wave equation.Mathematics of computation,1996:65(214):491-506
• 11 Anderson M,Kimn J H.A numerical approach to spacetime finite element for the wave equation.Journal of Computational Physics,2007:226(1):466-476
• 12 Thompson L L,He D T.Adaptive space-time finite element methods for the wave equation on unbounded domains. Computer methods in applied mechanics and engineering, 2005:194(18-20):1947-2000
• 13 Kurdi M H,Beran P S.Spectral element method in time for rapidly actuated systems.Journal of Computational Physics,2008:227(3):1809-1835
• 14 Yoseph P B,Moses E,Zrahia U,Yarin A.Space-time spectral element methods for one-dimensional nonlinear advection-diffusion problems.Journal of Computational Physics,1995:119(1):32-74
• 15 Tang J G,Ma H P.A legendre spectral method in time for first-order hyperbolic equations.Applied numerical mathematics, 2007:57(1):1-11
• 16 Maerschaick B D,Gerrisma M I.The use of Chebyshev polynomials in the space-time least-squares spectral element method.Numerical algorithms,2005:38(1-3):173- 196