石国权, 燕振国, 朱华君, 等. 非结构高阶CPR方法的子单元限制中的问题单元侦测[J]. 空气动力学学报, 2023, 41(2): 38−52. doi: 10.7638/kqdlxxb-2021.0406
引用本文: 石国权, 燕振国, 朱华君, 等. 非结构高阶CPR方法的子单元限制中的问题单元侦测[J]. 空气动力学学报, 2023, 41(2): 38−52. doi: 10.7638/kqdlxxb-2021.0406
SHI G Q, YAN Z G, ZHU H J, et al. Troubled cell detection in subcell limiting for high-order CPR method on unstructured meshes[J]. Acta Aerodynamica Sinica, 2023, 41(2): 38−52. doi: 10.7638/kqdlxxb-2021.0406
Citation: SHI G Q, YAN Z G, ZHU H J, et al. Troubled cell detection in subcell limiting for high-order CPR method on unstructured meshes[J]. Acta Aerodynamica Sinica, 2023, 41(2): 38−52. doi: 10.7638/kqdlxxb-2021.0406

非结构高阶CPR方法的子单元限制中的问题单元侦测

Troubled cell detection in subcell limiting for high-order CPR method on unstructured meshes

  • 摘要: 重构修正(correction procedure via reconstruction, CPR)是一种计算高效、模板紧致、适用于非结构网格的高阶方法,但是无法处理流场含有较强间断的问题。针对这个问题,前期基于二阶非等距非线性加权插值的激波捕捉格式以及非结构四边形网格,发展了高阶CPR方法的一种子单元限制技术。该技术首先通过侦测因子确定出问题单元,然后将问题单元基于CPR求解点划分成非等距子单元,并在子单元中采用二阶激波捕捉格式进行计算。这种子单元限制技术保证了主单元和子单元的求解点重合,可以十分方便地处理主单元和子单元之间的数据交换。问题单元的侦测是该子单元限制技术的关键步骤,侦测出的问题单元的数量以及分布情况都将会影响数值结果。基于此,本文同时考虑了单元侦测和分维侦测两种问题单元侦测方式。通过激波管问题、Shu-Osher问题、双马赫反射问题、激波-旋涡干扰等典型算例对比分析了TVB、MDHE和JST侦测因子对数值格式分辨率的影响。研究结果表明:分维侦测能够减少问题单元的区域,进而提高分辨率。

     

    Abstract: The correction procedure via reconstruction (CPR) method is an efficient and compact method, suitable for unstructured meshes, but it cannot solve the problem with strong discontinuities in flow fields. Based on the shock-capturing scheme with second-order non-uniform nonlinear weighted interpolation, we previously developed a subcell limiting technique for the high-order CPR method on two-dimensional unstructured quadrilateral meshes. Within this technique, the troubled cells are first divided into non-uniform subcells, then a second-order shock-capturing scheme is adopted for the calculation. As the solution points are overlaid between the main cells and the subcells, it is very convenient to handle the data exchange in the two layers. The detection of troubled cells is a crucial step in the subcell limiting technique, and the number and distribution of troubled cells can both affect the numerical results. Based on this, the present study considers two methods of troubled cell detecting, i.e. detecting by cell and detecting by dimension. This paper compares and analyzes the influence of different indicators, i.e. TVB, MDH, and JST, on the resolution of numerical schemes through typical benchmark cases such as the shock tube problem, the Shu-Osher problem, the double Mach reflection problem, and the shock-vortex interaction. The results show that detecting by dimension can reduce the number of troubled cells, thereby improving the resolution, which has a certain scientific value.

     

/

返回文章
返回