1. 中国空气动力研究与发展中心 空气动力学国家重点实验室,绵阳,621000
2. 中国空气动力研究与发展中心,绵阳,621000
3. 北京航空航天大学 能源与动力工程学院,北京,100191
刘旭亮*(1986-),男,助理研究员,研究方向:计算流体力学. E-mail:liuxuliang@cardc.cn
纸质出版:2023
移动端阅览
刘旭亮, 范召林, 张树海, 等. 一类新型自适应反扩散近似Riemann求解器及其应用[J]. 空气动力学学报, 2023, 41(4): 52−63.
LIU Xuliang, FAN Zhaolin, ZHANG Shuhai, et al. A new type of approximate Riemann solver with adaptive anti-diffusion and its application[J]. 2023, 41(4): 52-63.
刘旭亮, 范召林, 张树海, 等. 一类新型自适应反扩散近似Riemann求解器及其应用[J]. 空气动力学学报, 2023, 41(4): 52−63. DOI: 10.7638/kqdlxxb-2022.0009
LIU Xuliang, FAN Zhaolin, ZHANG Shuhai, et al. A new type of approximate Riemann solver with adaptive anti-diffusion and its application[J]. 2023, 41(4): 52-63. DOI: 10.7638/kqdlxxb-2022.0009.
对于包含激波、剪切层等复杂结构的流动问题,为了精确模拟剪切层等精细结构,且保证激波计算的稳定性,必须采用低耗散且强鲁棒的数值通量方法。传统的HLL近似Riemann求解器的耗散性较大,Roe、HLLEM和HLLC等近似Riemann求解器在计算某些含有强激波的物理问题时会出现非物理解,容易导致不稳定。针对这一问题,本文在Riemann求解器中通过合理设计反扩散矩阵,发展了一类具有自适应反扩散的新型Riemann求解器,并将其应用到高阶加权紧致格式,实现了高阶精度求解。通过典型数值算例验证了新型方法的计算精度和稳定性,结果表明本文提出的新型自适应反扩散Riemann求解器克服了传统Riemann求解器的缺陷,既能准确识别剪切层等精细结构,又能保证激波解的稳定性。
Constructing numerical fluxes using Riemann solvers is one of the most important steps for solving hyperbolic systems of conservation laws. However
it is still of great challenge for existing Riemann solvers to simulate high-speed flows with shock waves. The Riemann solver based on the HLL approximation has a high dissipation
while those based on the Roe
HLLEM
and HLLC approximation may result in non-physical solutions and numerical instability when simulating complex flows with strong shock waves. Because of this
a new type of Riemann solver with adaptive anti-diffusion is developed by a dedicated design of the anti-diffusion matrix. Due to its low dissipation
the proposed Riemann solver yields high-accuracy solutions when applied to high-order weighted compact finite difference schemes. The standard numerical experiments verify computational accuracy and stability for the new method
which show that the new Riemann solver overcomes the shortcomings of the traditional approximate Riemann solver. The new Riemann solver can accurately identify the fine structures such as shear layers
and maintain the numerical stability for the computation of shock wave.
LAX P D. Weak solutions of nonlinear hyperbolic equations and their numerical computation[J]. Communications on Pure and Applied Mathematics, 1954, 7(1): 159-193. DOI:10.1002/cpa.3160070112
STEGER J L, WARMING R F. Flux vector splitting of the inviscid gasdynamic equations with application to finite-difference methods[J]. Journal of Computational Physics, 1981, 40(2): 263-293. DOI:10.1016/0021-9991(81)90210-2
BRAM VAN LEER. Flux-vector splitting for the Euler equations[C]//8th International Conference on Numerical Methods in Fluid Dynamics, 1982. Lecture Notes in Physics, 170: 507-512. doi: 10.1007/3-540-11948-5_66
GODUNOV S K. A finite difference method for the computation of discontinuous solutions of the equations of fluid dynamics[R]. Mathematicheskii Sbornik 1959: 271–306.
ROE P L. Approximate Riemann solvers, parameter vectors, and difference schemes[J]. Journal of Computational Physics, 1981, 43(2): 357-372. DOI:10.1016/0021-9991(81)90128-5
HARTEN A, LAX P D, VAN LEER B. On upstream differencing and Godunov-type schemes for hyperbolic conservation laws[M]// SIAM Review, 1983, 25(1): 35–61. doi: 10.1137/1025002
EINFELDT B. On Godunov-type methods for gas dynamics[J]. SIAM Journal on Numerical Analysis, 1988, 25(2): 294-318. DOI:10.1137/0725021
EINFELDT B, MUNZ C D, ROE P L, et al. On Godunov-type methods near low densities[J]. Journal of Computational Physics, 1991, 92(2): 273-295. DOI:10.1016/0021-9991(91)90211-3
TORO E F, SPRUCE M, SPEARES W. Restoration of the contact surface in the HLL-Riemann solver[J]. Shock Waves, 1994, 4(1): 25-34. DOI:10.1007/BF01414629
OSHER S, CHAKRAVARTHY S. Upwind schemes and boundary conditions with applications to Euler equations in general geometries[J]. Journal of Computational Physics, 1983, 50(3): 447-481. DOI:10.1016/0021-9991(83)90106-7
LIOU M S, STEFFEN C J Jr. A new flux splitting scheme[J]. Journal of Computational Physics, 1993, 107(1): 23-39. DOI:10.1006/jcph.1993.1122
PEERY K, IMLAY S. Blunt-body flow simulations[C]//24th Joint Propulsion Conference, Boston, MA, USA. Reston, Virigina: AIAA, 1988. AIAA 88-2924. doi: 10.2514/6.1988-2904
HARTEN A. High resolution schemes for hyperbolic conservation laws[J]. Journal of Computational Physics, 1983, 49(3): 357-393. DOI:10.1016/0021-9991(83)90136-5
QUIRK J J. A contribution to the great Riemann solver debate[J]. International Journal for Numerical Methods in Fluids, 1994, 18(6): 555-574. DOI:10.1002/fld.1650180603
KIM S D, LEE B J, LEE H J, et al. Realization of contact resolving approximate Riemann solvers for strong shock and expansion flows[J]. International Journal for Numerical Methods in Fluids, 2009, 62(4): 1107–1133. DOI:10.1002/fld.2057
HUANG K B, WU H, YU H, et al. Cures for numerical shock instability in HLLC solver[J]. International Journal for Numerical Methods in Fluids, 2011, 65(9): 1026-1038. DOI:10.1002/fld.2217
FU L, GAO Z H, ZHANG X D. Construction and application research of HLL-HLLC scheme[J]. Acta Aerodynamica Sinica, 2014, 32(1): 116-122, 135. (in Chinese)
BALSARA D S. Multidimensional HLLE Riemann solver: application to Euler and magnetohydrodynamic flows[J]. Journal of Computational Physics, 2010, 229(6): 1970-1993. DOI:10.1016/j.jcp.2009.11.018
BALSARA D S, DUMBSER M, ABGRALL R. Multidimensional HLLC Riemann solver for unstructured meshes - with application to Euler and MHD flows[J]. Journal of Computational Physics, 2014, 261: 172-208. DOI:10.1016/j.jcp.2013.12.029
REN Y X. A robust shock-capturing scheme based on rotated Riemann solvers[J]. Computers Fluids, 2003, 32(10): 1379-1403. DOI:10.1016/S0045-7930(02)00114-7
NISHIKAWA H, KITAMURA K. Very simple, carbuncle-free, boundary-layer-resolving, rotated-hybrid Riemann solvers[J]. Journal of Computational Physics, 2008, 227(4): 2560-2581. DOI:10.1016/j.jcp.2007.11.003
DUMBSER M, BALSARA D S. A new efficient formulation of the HLLEM Riemann solver for general conservative and non-conservative hyperbolic systems[J]. Journal of Computational Physics, 2016, 304: 275-319. DOI:10.1016/j.jcp.2015.10.014
XIE W J, LI W, LI H, et al. On numerical instabilities of Godunov-type schemes for strong shocks[J]. Journal of Computational Physics, 2017, 350: 607-637. DOI:10.1016/j.jcp.2017.08.063
LIU X L, ZHANG S H, ZHANG H X, et al. A new class of central compact schemes with spectral-like resolution I: Linear schemes[J]. Journal of Computational Physics, 2013, 248: 235-256. . Doi: 10.1016/j.jcp.2013.04.014
LIU X L, ZHANG S H, ZHANG H X, et al. A new class of central compact schemes with spectral-like resolution II: hybrid weighted nonlinear schemes[J]. Journal of Computational Physics, 2015, 284: 133-154. DOI:10.1016/j.jcp.2014.12.027
BATTEN P, CLARKE N, LAMBERT C, et al. On the choice of wave speeds for the HLLC Riemann solver[J]. SIAM Journal on Scientific Computing, 1997, 18(6): 1553-1570. DOI:10.1137/s1064827593260140
DAVIS S F. Simplified second-order Godunov-type methods[J]. SIAM Journal on Scientific and Statistical Computing, 1988, 9(3): 445-473. DOI:10.1137/0909030
TORO E F. Riemann solvers and numerical methods for fluid dynamics[M]. Berlin, Heidelberg: Springer Berlin Heidelberg, 2009. doi: 10.1007/b79761
SANDERS R, MORANO E, DRUGUET M C. Multidimensional dissipation for upwind schemes: stability and applications to gas dynamics[J]. Journal of Computational Physics, 1998, 145(2): 511-537. DOI:10.1006/jcph.1998.6047
DENG X G, ZHANG H X. Developing high-order weighted compact nonlinear schemes[J]. Journal of Computational Physics, 2000, 165(1): 22-44. DOI:10.1006/jcph.2000.6594
TU G H, ZHAO X H, MAO M L, et al. Evaluation of Euler fluxes by a high-order CFD scheme: shock instability[J]. International Journal of Computational Fluid Dynamics, 2014, 28(5): 171-186. DOI:10.1080/10618562.2014.911847
SHU C W, OSHER S. Efficient implementation of essentially non-oscillatory shock-capturing schemes[J]. Journal of Computational Physics, 1988, 77(2): 439-471. DOI:10.1016/0021-9991(88)90177-5
LIU X L, ZHANG S H, ZHANG H X. Development of high-order weighted compact schemes with various difference methods[J]. Computers Fluids, 2016, 136: 114-131. DOI:10.1016/j.compfluid.2016.06.003
PIROZZOLI S. On the spectral properties of shock-capturing schemes[J]. Journal of Computational Physics, 2006, 219(2): 489-497. DOI:10.1016/j.jcp.2006.07.009
TU G H, DENG X G, MAO M L. Spectral property comparison of fifth-order nonlinear WCNS and WENO difference schemes[J]. Acta Aerodynamica Sinica, 2012, 30(6): 709-712. (in Chinese) doi: 10.7638/kqdlxxb-2014.0115
MAO M L, YAN Z G, LIU H Y, et al. Study of quasi-linear spectral analysis method of high-order weighted nonlinear schemes[J]. Acta Aerodynamica Sinica, 2015, 33(1): 1-9. (in Chinese)
JIANG G S, SHU C W. Efficient implementation of weighted ENO schemes[J]. Journal of Computational Physics, 1996, 126(1): 202-228. DOI:10.1006/jcph.1996.0130
TITAREV V A, TORO E F. Finite-volume WENO schemes for three-dimensional conservation laws[J]. Journal of Computational Physics, 2004, 201(1): 238-260. DOI:10.1016/j.jcp.2004.05.015
WOODWARD P, COLELLA P. The numerical simulation of two-dimensional fluid flow with strong shocks[J]. Journal of Computational Physics, 1984, 54(1): 115-173. DOI:10.1016/0021-9991(84)90142-6
0
浏览量
0
下载量
0
CSCD
关联资源
相关文章
相关作者
相关机构
京公网安备11010802024621