袁礼
- 作品数：27 被引量：48H指数：5
- 供职机构：中国科学院数学与系统科学研究院计算数学与科学工程计算研究所科学与工程计算国家重点实验室更多>>
- 发文基金：国家重点基础研究发展计划国家自然科学基金博士科研启动基金更多>>
- 相关领域：理学航空宇航科学技术天文地球建筑科学更多>>
- 一种健壮的TV通量分裂格式
- 2019年
- 随着计算流体力学的快速发展,设计精确、高效并且健壮的数值格式变得尤为重要。Toro等 [8] 提出的TV通量分裂格式表现出简单、高效和精确分辨接触间断等优点,但是在计算一些多维算例时会出现数值激波不稳定现象。两波近似的HLL格式在计算中非常高效和健壮,但是不能分辨接触间断大大地限制了其应用。本文对TV通量分裂格式进行稳定性分析,据此提出一种混合格式来消除TV格式的数值激波不稳定性。数值试验表明,本文构造的混合格式不仅保留了原始TV格式的优点,而且具有更好的健壮性,在计算二维问题时不会出现数值激波不稳定现象。
- 胡立军袁礼
- 关键词：欧拉方程
- 一种治愈强激波数值不稳定性的混合方法被引量：4
- 2015年
- HLLC(Harten-Lax-Leer-contact)格式是一种高分辨率格式,能够准确捕捉激波、接触间断和稀疏波.但是使用HLLC格式计算多维问题时,在强激波附近会出现激波不稳定现象.FORCE(first-order centred)格式在强激波附近表现出很好的稳定性,并且其数值耗散比HLL(HartenLax-Leer)格式小.分析了HLLC格式和FORCE格式在特定流动条件下的稳定性,构造了HLLCFORCE混合格式并且进一步结合开关函数来消除HLLC格式的激波不稳定现象.数值试验表明新构造的混合格式不仅能够消除HLLC格式的激波不稳定现象,还最大程度地保留HLLC格式高分辨率的优点.
- 胡立军袁礼
- 关键词：稳定性分析
- 颗粒流Savage-Hutter方程的有限体积法及其GPU计算
- Savage-Hutter( SH)方程是描述Coulomb颗粒材料在地表面上流动的薄层近似方程，是滑坡、泥石流等地质灾害研究中的一个重要的力学模型。本文对SH方程的一维和二维形式，利用Godunov有限体积法进行积分离...
- 翟健袁礼
- 关键词：力学模型有限体积法
- 基于PHG平台的非结构四面体网格欧拉方程间断有限元并行求解器
- 2021年
- 针对计算流体力学对高性能计算的需求,基于三维并行自适应有限元程序开发平台PHG(Parallel Hierarchical Grid)开发了在非结构四面体网格上求解可压缩流欧拉方程的间断有限元法并行求解器(Libdgphg库).该求解器以C++函数库的形式实现数值方法中各项功能.实施了模态基一次间断有限元,采用低耗散的MLP(Multi-dimensional Limiting Process)限制器来抑制间断附近的数值振荡.由于MLP限制器需要所有与当前单元共享顶点的邻近单元的信息,模板较宽,这给程序设计带来一定的困难.我们通过引入辅助向量收集共享顶点的所有单元中的最大、最小单元积分平均值,并归属到单元数据结构上,从而利用PHG内在的通信机制实现MPI分区间的信息交换.通过几个数值算例测试了Libdgphg库的数值结果以及并行性能.算例表明:该求解器能得到理论精度阶和较高分辨率,同时有良好的并行性能,在千核测试中可达到60%以上的并行效率,可用于流体问题的大规模计算.
- 杜玉龙徐凯文赵昆磊袁礼
- 关键词：双曲守恒律间断有限元法
- 求解含源项守恒律组的间断有限元方法及在反应流计算中的应用
- <正>由Cockburn和Shu等人构建和发展的龙格库塔间断Galerkin有限元方法(Runge-Kutta DiscontinuousGalerkin Method,RKDG)是一种求解含间断现象问题的方法。RKDG...
- 张磊袁礼
- 文献传递
- 求解可压缩流的二维通量分裂格式被引量：2
- 2020年
- 传统的一维通量分裂格式在计算界面数值通量时,只考虑网格界面法向的波系。采用传统的TV格式分别求解对流通量和压力通量。通过求解考虑了横向波系影响的角点数值通量来构造一种真正二维的TV通量分裂格式。在计算一维数值算例时,该格式与传统的TV格式具有相同的数值通量计算公式,因此其保留了传统的TV格式精确捕捉接触间断和膨胀激波的优点。在计算二维算例时,该格式比传统的TV格式具有更高的分辨率;在计算二维强激波问题时,消除了传统TV格式的非物理现象,表现出更好的鲁棒性;此外,该格式大大提高了稳定性CFL数,从而具有更高的计算效率。因此,本文方法是一种精确、高效并且具有强鲁棒性的数值方法,在可压缩流的数值模拟中具有广阔的应用前景。
- 胡立军翟健袁礼
- 关键词：欧拉方程可压缩流鲁棒性
- 移动网格法在模拟预混合氢氧气体燃烧中的应用
- 将移动网格法应用于化学反应流的数值模拟。整个计算方法由流场计算和网格移动两个独立的部分组成。在流场计算中应用二阶迎风TVD格式和LU-SGS格式求解轴对称多组分反应气体的
- 袁礼
- 关键词：移动网格化学反应流
- 文献传递
- 应用非交错NOC格式求解两相浅水流方程的多GPU计算
- 给出了求解两相浅水流动方程的非时空交错无震荡差分(NOC)格式的应用,并利用多GPU对计算进行加速.NOC格式无需求解近似Riemann问题.基于这一优势,应用一个NOC的简单变化形式,来求解流动模型的数值解.计算中发现...
- 翟健袁礼刘伟
- 用改进的LevelSet方法计算可压缩双介质流动
- 用带有虚拟流体（GhostFluid）修正的LevelSet方法计算了一维可压缩双介质流动。把描述流动的Euler方程和描述流体界面运动的LevelSet方程耦合起来,得到一个整体的双曲型守恒律系统,从而可应用高分辨率差...
- 张镭袁礼
- 文献传递
- 求解二维Euler方程有限单元边插值的降维重构算法被引量：1
- 2019年
- 数值求解二维Euler方程的有限体积法(如k-exact,WENO重构、紧致重构等),无一例外地要进行耗时的网格单元上的二维重构.然而这些二维重构最后仅用于确定网格单元边界上高斯积分点处的解值,单元上二维重构似乎并非必需的.因此,文章提出用网格边上的一维重构来取代有限体积法中网格单元上的二维重构,分别在一致矩形网格和非结构三角形网格上发展了基于网格边重构的求解二维Euler方程的新方法,称为降维重构算法.数值算例表明该算法可以计算有强激波的无黏流动问题,且有较高的计算效率.
- 王亚辉刘伟袁礼袁礼
- 关键词：矩形网格三角形网格EULER方程守恒律
