Hostname: page-component-7b7567c8c-tnzcf Total loading time: 0 Render date: 2026-11-15T05:05:33.760Z Has data issue: false hasContentIssue false

Osher Flux with Entropy Fix for Two-Dimensional Euler Equations

Published online by Cambridge University Press:  27 May 2016

Huajun Zhu*
Affiliation:
State Key Laboratory of Aerodynamics, China Aerodynamics Research and Development Center, Mianyang, Sichuan 621000, China
Xiaogang Deng
Affiliation:
National University of Defense Technology, Changsha, Hunan 410073, China
Meiliang Mao
Affiliation:
State Key Laboratory of Aerodynamics, China Aerodynamics Research and Development Center, Mianyang, Sichuan 621000, ChinaComputational Aerodynamics Institute, China Aerodynamics Research and Development Center, Mianyang, Sichuan 621000, China
Huayong Liu
Affiliation:
State Key Laboratory of Aerodynamics, China Aerodynamics Research and Development Center, Mianyang, Sichuan 621000, China
Guohua Tu
Affiliation:
State Key Laboratory of Aerodynamics, China Aerodynamics Research and Development Center, Mianyang, Sichuan 621000, China
*
*Corresponding author. Email:hjzhu@skla.cardc.cn (H. J. Zhu)

Abstract

Core share and HTML view are not available for this content. However, as you have access to this content, a full PDF is available via the 'Save PDF' action button.

We compare in this paper the properties of Osher flux with O-variant and P-variant (Osher-O flux and Osher-P flux) in finite volume methods for the two-dimensional Euler equations and propose an entropy fix technique to improve their robustness. We consider both first-order and second-order reconstructions. For inviscid hypersonic flow past a circular cylinder, we observe different problems for different schemes: A first-order Osher-O scheme on quadrangular grids yields a carbuncle shock, while a first-order Osher-P scheme results in a dislocation shock for high Mach number cases. In addition, a second-order Osher scheme can also yield a carbuncle shock or be unstable. To improve the robustness of these schemes we propose an entropy fix technique, and then present numerical results to show the effectiveness of the proposed method. In addition, the influence of grid aspects ratio, relative shock position to the grid and Mach number on shock stability are tested. Viscous heating problem and double Mach reflection problem are simulated to test the influence of the entropy fix on contact resolution and boundary layer resolution.

Information

Type
Research Article
Copyright
Copyright © Global-Science Press 2016