基于重构的修正过程(correction procedure via reconstruction,CPR)方法是一种适用于非结构网格、紧致高效的高阶方法,但在离散非线性对流项时容易因为混淆误差的累积而出现数值不稳定。本文研究了基于LG(Legendre-Gauss)点分裂形式的...基于重构的修正过程(correction procedure via reconstruction,CPR)方法是一种适用于非结构网格、紧致高效的高阶方法,但在离散非线性对流项时容易因为混淆误差的累积而出现数值不稳定。本文研究了基于LG(Legendre-Gauss)点分裂形式的CPR方法在欠解析流动中的稳定性,并将其与子单元限制技术结合,求解含激波的欠解析流动问题。结果表明:经过边界通量修正的LG点分裂形式CPR格式能够满足离散守恒律,且与子单元限制相结合后守恒律依然成立;在无激波的欠解析流动问题中,相比于散度形式CPR格式,分裂形式明显提高了计算的稳定性,而且比使用LGL(Legendre-Gauss-Lobatto)点的分裂形式CPR格式具有更小的数值误差;在含激波的欠解析流动求解中,相比基于LGL点的间断伽辽金谱元法子单元限制策略,基于LG点的分裂形式CPR格式子单元限制策略具有更高的分辨率和更小的振荡。展开更多
文摘基于重构的修正过程(correction procedure via reconstruction,CPR)方法是一种适用于非结构网格、紧致高效的高阶方法,但在离散非线性对流项时容易因为混淆误差的累积而出现数值不稳定。本文研究了基于LG(Legendre-Gauss)点分裂形式的CPR方法在欠解析流动中的稳定性,并将其与子单元限制技术结合,求解含激波的欠解析流动问题。结果表明:经过边界通量修正的LG点分裂形式CPR格式能够满足离散守恒律,且与子单元限制相结合后守恒律依然成立;在无激波的欠解析流动问题中,相比于散度形式CPR格式,分裂形式明显提高了计算的稳定性,而且比使用LGL(Legendre-Gauss-Lobatto)点的分裂形式CPR格式具有更小的数值误差;在含激波的欠解析流动求解中,相比基于LGL点的间断伽辽金谱元法子单元限制策略,基于LG点的分裂形式CPR格式子单元限制策略具有更高的分辨率和更小的振荡。
