A discontinuous Galerkin method/HLLC solver for the Euler equations

This paper proposes a fully three‐dimensional non‐linear Euler methodology for solving aerodynamic and acoustic problems in the presence of strong shocks and rarefactions. It uses a discontinuous Galerkin method (DGM) within the element, and a Riemann solver (HLLC) at the boundaries to propagate rarefactions while preserving the entropy condition and capturing shocks with no spurious oscillations. This approach is thought to marry the best aspects of finite element and finite volume methods, achieving conservation while not requiring the solution of a large matrix. Examples in which shock and rarefaction waves are well captured are presented and the propagation of acoustic pulses is well demonstrated. Copyright © 2003 John Wiley & Sons, Ltd.     

一种适用于磁流体力学切向间断的HLLC黎曼算子

磁场的存在使得磁流体力学特征波模不同于流体力学, 因此直接由流体力学HLLC黎曼算子导出的HLLC双中间态在交界的间断面会出现不守恒的问题。通常降级采用HLL单磁场中间态代替HLLC双磁场中间态以实现守恒和计算稳定, 代价是切向间断的模拟精度不足。本文对此进行改进, 在模拟切向间断时仍然保留原有的HLLC双磁场中间态, 同时各守恒量仍然能够满足Toro相容条件; 改进型HLLC算子在间断两侧的磁场分量存在差异, 因此能够更精确还原切向间断面。基于数值测试, 包括一维激波管和切向间断的时变模拟, 以及地球磁层三维数值模拟, 将模拟结果进行对比, 结果表明: 相比于已发展的HLLC算子, 改进型HLLC算子对切向间断具有更好的捕捉精度, 能够达到或接近耗时更多的HLLD算子的模拟精度。     

Energy relaxation method for chemical non‐equilibrium flow computations

In this paper, an algorithm for chemical non‐equilibrium hypersonic flow is developed based on the concept of energy relaxation method (ERM). The new system of equations obtained are studied using finite volume method with Harten–Lax–van Leer scheme for contact (HLLC). The original HLLC method is modified here to account for additional species and split energy equations. Higher order spatial accuracy is achieved using MUSCL reconstruction of the flow variables with van Albada limiter. The thermal equilibrium is considered for the analysis and the species data are generated using polynomial correlations. The single temperature model of Dunn and Kang is used for chemical relaxation. The computed results for a flow field over a hemispherical cylinder at Mach number of 16.34 obtained using the present solver are found to be promising and computationally (25%) more efficient. The present solver captures physically correct solution as the entropy conditions are satisfied automatically during the computations. Copyright © 2007 John Wiley & Sons, Ltd.     

用HLLC方法处理运动边界

针对含有运动边界非结构网格点的显式有限体积方法，研究了如何使用HLLC方法处理运动边界条件，计算了活塞拉伸、压缩问题，并和解析解进行了比较．作为应用示例，比较了静止圆柱绕流和静止流场中运动圆柱的绕流情形。结果表明：应用HLLC方法处理运动边界是可行的。     

管系内一维可压缩流动数值建模与特性分析

针对复杂管系内可压缩流体,基于有限体积法,采用HLLC(Harten-Lax-van Leer Contact)格式和黎曼求解器构建了有限控制体数值离散方法,引入虚拟节点用于连接有限控制体,借助虚拟节点给出控制体之间数值通量的计算格式,发展了一种管道内一维流动数值建模方法.针对含有分支管路的管系,在管道连接部位构建了分...     

一种简单的精确捕捉接触间断的黎曼求解器

基于Godunov型数值格式的有限体积法是求解双曲型守恒律系统的主流方法,其中用来计算界面数值通量的黎曼求解器在很大程度上决定了数值格式在计算中的表现。单波的Rusanov求解器和双波的HLL求解器具有简单、高效和鲁棒性好等优点,但是在捕捉接触间断时耗散太大。全波的HLLC格式能够精确捕捉接触间断,但是在计算中出现的激波不稳定现象限制了其在高马赫数流动问题中的应用。本文利用双曲正切函数和五阶WENO格式来重构界面两侧的密度值,并且结合边界变差下降算法来减小Rusanov格式耗散项中的密度差,从而提高格式对于接触间断的分辨率。研究表明,相比于全波的HLLC求解器,本文构造的黎曼求解器不仅具有更高的接触分辨率,而且还具有更好的激波稳定性。     

一维含化学反应流体力学方程组的ALE间断有限元方法

研究一维含化学反应流体力学方程组的数值模拟方法.结合理想气体状态方程并利用HLLC解法器在各个单元边界处的数值通量,给出ALE间断有限元方法.高阶计算时,使用TVD斜率限制器对数值解可能产生的非物理振荡进行抑制.结果表明：该算法能够保持物理量的守恒性和高精度,并能够清晰地捕捉爆轰波的结构特征.     

HLLC scheme with novel wave‐speed estimators appropriate for two‐dimensional shallow‐water flow on erodible bed

This paper presents a first‐order HLLC (Harten‐Lax‐Van Leer with contact discontinuities) scheme to solve the Saint‐Venant shallow‐water equations, including morphological evolution of the bed by erosion and deposition of sediments. The Exner equation is used to model the morphological evolution of the bed, while a closure equation is needed to evaluate the rate of sediment transport. The system of Saint‐Venant–Exner equations is solved in a fully coupled way using a finite‐volume technique and a HLLC solver for the fluxes, with a novel wave‐speed estimator adapted to the Exner equation. Wave speeds are usually derived by computing the eigenvalues of the full system, which is highly time‐consuming when no analytical expression is available. In this paper, an eigenvalue analysis of the full system is conducted, leading to simple but still accurate wave‐speed estimators. The new numerical scheme is then tested in three different situations: (1) a circular dam‐break flow over movable bed, (2) an one‐dimensional bed aggradation problem simulated on a 2D unstructured mesh and (3) the case of a dam‐break flow in an erodible channel with a sudden enlargement, for which experimental measurements are available. Copyright © 2010 John Wiley & Sons, Ltd.     

基于HLL-HLLC的高阶WENO格式及其应用研究

HLL-HLLC格式能够克服HLLC在强激波附近的激波不稳定现象,并且保持了HLLC的低耗散特性,是一种适合更大马赫数范围的近似黎曼求解器。本文从RANS方程出发,将HLL-HLLC近似黎曼求解器结合五阶WENO重构,实现了对无粘通量的高阶离散;同时,采用完全守恒形式的四阶中心差分格式处理粘性项,建立了RANS 方程的高阶数值求解格式。通过对四个经典算例,钝头体、 ONERA M6机翼、DLR F6-WB翼身组合体和DLR F6-WBNP复杂外形的数值模拟,考察了两种WENO改进格式在复杂流场中的表现,研究了高阶格式的收敛特性;给出了在复杂流动中 WENO自由参数的推荐值,以增强求解的收敛性。算例结果表明,本文构造的高阶格式鲁棒性好,能够显著改善激波位置和激波强度,捕获更丰富的流场细节,满足复杂工程应用需求。