基于混合重构高阶间断伽辽金方法的二维层流和湍流数值模拟

doi:10.12066/j.issn.1007-2861.2242

上海大学学报(自然科学版) ›› 2022, Vol. 28 ›› Issue (4): 619-631.doi: 10.12066/j.issn.1007-2861.2242

基于混合重构高阶间断伽辽金方法的二维层流和湍流数值模拟

熊为¹, 张卫国², 丁珏¹, 杨小权¹(), 翁培奋¹

1.上海大学力学与工程科学学院, 上海 200444
2.中国空气动力研究与发展中心旋翼空气动力学重点实验室, 四川绵阳 621000

收稿日期:2020-02-13 出版日期:2022-08-30 发布日期:2022-08-29
通讯作者: 杨小权 E-mail:quanshui@shu.edu.cn
作者简介:杨小权(1983—), 男, 副教授, 博士生导师, 博士, 研究方向为计算流体力学、空气动力学等. E-mail: quanshui@shu.edu.cn
基金资助:
国家自然科学基金资助项目(11702329);空气动力学国家重点实验室开放课题资助项目(SKLA201801);旋翼空气动力学重点实验室开放课题资助项目(RAL20180403)

Methods for two-dimensional laminar and turbulent numerical simulations

XIONG Wei¹, ZHANG Weiguo², DING Jue¹, YANG Xiaoquan¹(), WENG Peifen¹

1. School of Mechanics and Engineering Science, Shanghai University, Shanghai 200444, China
2. State Key Laboratory of Aerodynamics, China Aerodynamics Research and Development Center, Mianyang 621000, Sichuan, China

Received:2020-02-13 Online:2022-08-30 Published:2022-08-29
Contact: YANG Xiaoquan E-mail:quanshui@shu.edu.cn

摘要/Abstract

摘要：

为了提高间断伽辽金方法的计算效率, 解决 least-squares 重构方法无法满足 2-exact 的缺陷, 发展了基于 recovery 重构和 least-squares 重构相结合的三阶混合重构方法, 用于求解可压缩层流和湍流流动. 将 Navier-Stokes 方程和修正的一方程 Negative Spalart-Allmaras 模型方程耦合成为系统方程, 采用三阶重构间断伽辽金方法进行求解. 时间推进采用基于半解析精确 Jacobian 矩阵的上-下对称高斯赛德尔格式(lower-upper symmetric Gauss-Seidel scheme, LU-SGS)预处理广义极小剩余(generalized minimal residual, GMRES)方法和四阶隐式 Runge-Kutta 方法; 空间对流项离散采用 Haten-Lax-van Leer 接触(Haten-Lax-van Leer contact, HLLC) 格式; 黏性项离散采用第二Bassi-Rebay (second Bassi-Rebay, BR2)格式, 并对 BR2 局部和全局提升算子开展三阶重构, 达到提高计算精度的目的. 通过典型算例验证了发展 rDGP₁P₂方法的准确性和计算效率. 研究结果表明: 重构的 rDGP₁P₂方法不仅具有较高的计算精度, 而且还具有较高的计算效率.

关键词: 高阶精度, 间断伽辽金方法, 混合重构, 第二Bassi-Rebay (second Bassi-Rebay, BR2)格式

Abstract:

Third-order hybrid reconstructed methods based on the combination of least-squares recovery and reconstruction were developed for solving compressible laminar and turbulent flows to improve the calculation efficiency of the discontinuous Galerkin methods and overcome the deficiency of the least-squares reconstruction in satisfying the two-exact property. Navier-Stokes equations and the modified equation of the negative Spalart-Allmaras model were coupled as the equation system and solved using the developed third-order reconstructed discontinuous Galerkin methods. The lower-upper symmetric Gauss-Seidel preconditioning generalized minimal residual method based on the semianalytical exact Jacobian matrix and the fourth-order implicit Runge-Kutta method were implemented in temporal advance. Harten-Lax-van Leer contact and second Bassi-Rebay (BR2) schemes were adopted to calculate the inviscid and viscous terms, respectively. Third-order reconstruction was applied to solve the local and global lifting operators of the BR2 scheme, improving the calculation accuracy. The benchmarks were selected to verify the accuracy and calculation efficiency of the developed rDGP₁P₂ method. The research results indicate that the developed reconstructed rDGP₁P₂ methods have high computational accuracy and efficiency.

Key words: high-order accuracy, discontinuous Galerkin method, hybrid reconstruction, second Bassi-Rebay (BR2) scheme

中图分类号:

V221.3

熊为, 张卫国, 丁珏, 杨小权, 翁培奋. 基于混合重构高阶间断伽辽金方法的二维层流和湍流数值模拟[J]. 上海大学学报(自然科学版), 2022, 28(4): 619-631.

XIONG Wei, ZHANG Weiguo, DING Jue, YANG Xiaoquan, WENG Peifen. Methods for two-dimensional laminar and turbulent numerical simulations[J]. Journal of Shanghai University（Natural Science Edition）, 2022, 28(4): 619-631.

图/表 11

图1

图2

图3

表1

图4

图5

图6

图7

图8

图9

图10

参考文献 14

[1]	Dumbser M, Balsara D S, Toro E F, et al. A unified framework for the construction of one-step finite volume and discontinuous Galerkin schemes on unstructured meshes[J]. Journal of Computational Physics, 2008, 227(18): 8209-8253. doi: 10.1016/j.jcp.2008.05.025
[2]	Dumbser M. Arbitrary high order PNPM schemes on unstructured meshes for the compressible Navier-Stokes equations[J]. Computers & Fluids, 2010, 39(1): 60-76. doi: 10.1016/j.compfluid.2009.07.003
[3]	Van Leer B, Lo M, Van Raalte M. A discontinuous Galerkin method for diffusion based on recovery[C]// AIAA Computational Fluid Dynamics Conference. 2007: 4083.
[4]	Van Raalte M, Van Leer B. Bilinear forms for the recovery-based discontinuous Galerkin method for diffusion[J]. Communications in Computational Physics, 2009, 5(2): 683-693.
[5]	Luo H, Luo L, Nourgaliev R, et al. A reconstructed discontinuous Galerkin method for the compressible Navier-Stokes equations on arbitrary grids[J]. Journal of Computational Physics, 2010, 229(19): 6961-6978. doi: 10.1016/j.jcp.2010.05.033
[6]	Xia Y, Luo H, Frisbey M, et al. A set of parallel, implicit methods for a reconstructed discontinuous Galerkin method for compressible flows on 3D hybrid grids[J]. Computers & Fluids, 2014, 98: 134-151. doi: 10.1016/j.compfluid.2014.01.023
[7]	Cheng J, Liu T G, LUO H. A hybrid reconstructed discontinuous Galerkin method for compressible flows on arbitrary grids[J]. Computers & Fluids, 2016, 139: 68-79. doi: 10.1016/j.compfluid.2016.04.001
[8]	Cheng J, Yang X Q, Liu T G, et al. A direct discontinuous Galerkin method for the compressible Navier-Stokes equations on arbitrary grids[J]// Journal of Computational Physis, 2016, 327: 484-502.
[9]	杨小权, 杨爱明, 孙刚. 一种强耦合Spalart-Allmaras湍流模型的RANS方程的高效数值计算方法[J]. 航空学报, 2013, 34(9): 2007-2018. doi: 10.7527/S1000-6893.2013.0056
[10]	Yang X Q, Cheng J, Liu X, et al. A reconstructed direct discontinuous Galerkin method for compressible turbulent flows on hybrid grids[J]. Computers & Fluids, 2018, 168: 216-231. doi: 10.1016/j.compfluid.2018.04.011
[11]	Bassi F, Crivellini A, Rebay S, et al. Discontinuous Galerkin solution of the Reynolds-averaged Navier-Stokes and $k$-$\Omega $ turbulence model equations[J]. Computers & Fluids, 2005, 34(4/5): 507-540. doi: 10.1016/j.compfluid.2003.08.004
[12]	Yang X Q, Cheng J, Luo H. A robust implicit direct discontinuous Galerkin method for simulating the compressible laminar and turbulent flows[J]. AIAA Journal, 2019, 57(3): 1113-1132. doi: 10.2514/1.J057172
[13]	Batten P, Leschziner M A, Goldberg U C. Average-state Jacobians and implicit methods for compressible viscous and turbulent flows[J]. Journal of Computational Physics, 1997, 137(1): 38-78. doi: 10.1006/jcph.1997.5793
[14]	Kristensen M R, Jrgensen J B, Thomsen P G, et al. An ESDIRK method with sensitivity analysis capabilities[J]. Computers & Chemical Engineering, 2004, 28(12): 2695-2707. doi: 10.1016/j.compchemeng.2004.08.004

基于混合重构高阶间断伽辽金方法的二维层流和湍流数值模拟

Methods for two-dimensional laminar and turbulent numerical simulations

RichHTML

PDF

可视化

摘要/Abstract

引用本文

使用本文

图/表 11

参考文献 14

相关文章 0

编辑推荐

Metrics

本文评价