期刊
COMPUTERS & FLUIDS
卷 136, 期 -, 页码 35-47出版社
PERGAMON-ELSEVIER SCIENCE LTD
DOI: 10.1016/j.compfluid.2016.05.026
关键词
Shock/turbulence interaction; Shock-unsteadiness; Turbulent kinetic energy; Turbulent dissipation rate; Conservative formulation; RANS
High-speed turbulent flows with shock waves pose significant challenge in terms of physical modeling and numerical accuracy. In this paper, we develop a k-epsilon turbulence model for canonical shock-turbulence interaction, based on the physical ideas of shock-unsteadiness damping. The shock-unsteadiness k-e. model proposed earlier by Sinha et al. (Phys. Fluids, 15(8), 2003) has model parameters that are functions of upstream Mach number. Our main objective is to eliminate the dependence of the model parameters on upstream flow quantities, and thus overcome difficulties in numerical implementation in complex flow applications. For the turbulent kinetic energy, this is achieved by proposing a new model for the shock-unsteadiness source term that is physically and dimensionally consistent with the earlier formulation. An alternate form of the dissipation rate equation is also used to eliminate upstream dependence of the model parameter in the production term. The proposed model is employed to compute the canonical interaction of a normal shock with homogeneous isotropic turbulence. Large numerical errors at the shock wave are eliminated by employing a transformation of the turbulence variables that result in a conservative form of the governing equations. The new model thus gives physically consistent and numerically accurate results at shock waves. The model predictions for turbulence kinetic energy and its dissipation rate are found to match DNS data over a range of Mach numbers. (C) 2016 Elsevier Ltd. All rights reserved.
作者
我是这篇论文的作者
点击您的名字以认领此论文并将其添加到您的个人资料中。
推荐
暂无数据