谱元法的指数收敛特性

您所在的位置:网站首页 指数方程计算方法 谱元法的指数收敛特性

谱元法的指数收敛特性

2023-06-23 14:45| 来源: 网络整理| 查看: 265

谱元法不仅能像h型有限元法一样通过增加网格数量来处理任意复杂地形,而且能像p型有限元法和谱方法一样通过提高插值阶数来达到非常高的精度,同时具有谱方法的指数收敛特性。本节自编Matlab程序,用二维亥姆霍兹方程验证谱元法的高精度和指数收敛特性,这是出于两方面考虑:1)时间离散采用高阶分裂法,解不可压缩N-S方程时存在速度场的亥姆霍兹方程,如式;2)亥姆霍兹方程通常存在解析解,便于误差对比分析。

考虑二维区域内亥姆霍兹方程:                    

其解析解为:。

为研究谱元法的计算精度及收敛特性,对比h型和p型展开下误差随自由度变化的收敛情况。对于h型展开,取插值阶数P=2,增加单元数改变自由度;对于p型展开,将区域划分为2×2共4个单元,通过提高插值阶数P改变自由度。对比和两种误差形式,其中表示误差绝对值的最大值,表示误差绝对值的和。

h型展开和p型展开下,误差随自由度的变化曲线如图2.3。可以发现,对于h型展开,误差随自由度增加衰减缓慢,最大误差始终在0.01以上。对于p型展开,误差随自由度增加呈指数减小,最大误差低于10-12。

由此可见,h型展开方法,如有限元法,误差随单元数增加的收敛速度十分缓慢,且误差始终较大;p型展开方法,如谱方法和谱元法,误差随着单元数的增加呈指数衰减,误差远小于h型展开方法。因此,当模拟的问题对精度要求很高时,谱元法具有极大的优势。



【本文地址】


今日新闻


推荐新闻


CopyRight 2018-2019 办公设备维修网 版权所有 豫ICP备15022753号-3