Computation of shock/boundary-layer interactions in bump channels with transport-type turbulence models

In this paper, an explicit time-marching finite-volume scheme has been used together with a number of convergence acceleration techniques such as the multigrid strategy. Two types of turbulence models, a Johnson-King (J-K) mo del and a two-layer k-epsilon/k-l model, have been incorporated and modifie d to model internal compressible flows with multiple walls. Some modificati ons have been made of the inner layer viscosity formulations of the J-K mod el in order to improve its predictive capability for flow separation. Parti ally implicit treatment of the transport-type equations of turbulence in th e models is adopted, because the source terms in these equations can cause numerical stiffness when there are flow separation, sharp gradients and hig h cell-aspect ratio near the solid wall. A two-dimensional are-bump flow in vestigated experimentally by Liu and Squire [X. Liu and L.C. Squire, Intera ction on curved surface at transonic speed, in: Turbulent Shear/Shock Wave Interactions, IUTAM Symposium Palaiseau 1985 (Springer, Berlin Heidelbeg, 1 985) 93-104.] was calculated using the J-K model with satisfactory agreemen t with the corresponding measurement. Although efficient and accurate, it i s found that the J-K model lacks the theoretical generality to be extended to model three-dimensional (3D) complex internal flows with multiple walls. Therefore, a two-layer k-epsilon model is employed for 3D how computation. Various measures are adopted to ensure stable and convergent numerical sol ution. A three-dimensional transonic channel flow with multiple shock/bound ary layer interactions was studied with the aforementioned two-layer model and numerical methods. The results are compared with experimental measureme nts [J. Cahen, V. Couaillier, J. Delery and T. Pot, Validation of Navier-St okes code using a k-epsilon turbulence model applied to a three-dimensional transonic tunnel, AIAA paper AIAA-93-0293, AIAA, 1993] and numerical resul ts obtained by using a Low-Reynolds-Number (LRN) k-epsilon model [VUB/FFA, Turbulence Models in EURANUS and the 3D Delery bump, Technical Report SNWP3 .3/01,VUB, Pleinlaan 2, 1050 Brussels, Belgium and FFA, P.O. Box 11021, S-1 61 11 Bromma, Sweden, 1993]. Compared with other (LRN) two-equation models, the two-layer model implemented is promising in modeling very complex 3D i nternal flows in terms of efficiency, robustness and accuracy. The two-laye r model permits uniform distribution of flow properties to be specified as initial condition which makes the simulation easier to be carried out. (C) 1999 Elsevier Science S.A. All rights reserved.