计算物理 ›› 2018, Vol. 35 ›› Issue (1): 22-28.DOI: 10.19596/j.cnki.1001-246x.7573

• • 上一篇    下一篇

非结构多面体二阶局部保界全局重映算法

徐喜华1, 刘娜1,2,3, 陈艺冰1,2   

  1. 1. 北京应用物理与计算数学研究所, 北京 100088;
    2. 计算物理实验室, 北京应用物理与计算数学研究所, 北京 100088;
    3. 中物院高性能数值模拟软件中心, 北京 100088
  • 收稿日期:2016-11-08 修回日期:2017-02-22 出版日期:2018-01-25 发布日期:2018-01-25
  • 作者简介:徐喜华(1986-),男,江苏省张家港,博士后,从事计算流体力学方法研究,E-mail:xihuaxu@126.com
  • 基金资助:
    国家自然科学基金(11701036,11671050,11501043,U1630247,91430218),科学挑战计划(JCKY2016212A502)和国家高技术研究发展计划(2015AA01A304)资助项目

Second-order Local-Bound-Preserving Conservative Remapping on Unstructured Polyhedral Meshes

XU Xihua1, LIU Na1,2,3, CHEN Yibing1,2   

  1. 1. Institute of Applied Physics and Computational Mathematics, Beijing 100088, China;
    2. Laboratory of Computational Physics, IAPCM, Beijing 100088, China;
    3. CAEP Software Center for High Performance Numerical Simulation, Beijing 100088, China
  • Received:2016-11-08 Revised:2017-02-22 Online:2018-01-25 Published:2018-01-25

摘要: 提出一种三维非结构多面体二阶保界全局重映算法.在旧网格上选取模板利用最小二乘构造插值多项式,采用凸包算法计算多面体相交部分,最后使用局部保界修正技术修补重映后的越界量.多项数值实验表明这种格式同时具有高精度、高分辨率和高效率的特点.

关键词: 全局重映, 局部保界算法, 多面体求交, 贡献网格方法

Abstract: We present conservatively remapping cell-centered variables from one mesh to another with second-order accuracy and boundary-preservation. It is generally applicable to any polyhedral source or target mesh. The algorithm consists of four parts:A least square based polynomial reconstruction of physical gradient; an octree-based fast donor-cell searing algorithm; a convex hull algorithm for intersection of polyhedrons and a modifying procedure for local bound preservation. The remapping scheme is scalable, second-order accurate and enjoys bound preservation property. Various benchmark problems demonstrate these properties. Numerical results show that it takes hundreds seconds to remap physical variables on tessellation with hundreds thousands to millions polyhedrons.

Key words: global remap, local bound preservation method, intersection of polyhedrons, donor cell method

中图分类号: