Numerical investigations on interactions between 2D/3D conical shock wave and axisymmetric boundary layer at Ma=2.2

Numerical simulation and analysis are carried out on interactions between a 2D/3D conical shock wave and an axisymmetric boundary layer with reference to the experiment by Kussoy et al., in which the shock was generated by a 15-deg half-angle cone in a tube at 15-deg angle of attack (AOA). Based on the RANS equations and Menter's SST turbulence model, the present study uses the newly developed WENO3-PRM211 scheme and the PHengLEI CFD platform for the computations. First, computations are performed for the 3D interaction corresponding to the conditions of the experiment by Kussoy et al., and these are then extended to cases with AOA = 10-deg and 5-deg. For comparison, 2D axisymmetric counterparts of the 3D interactions are investigated for cones coaxial with the tube and having half-cone angles of 27.35-deg, 24.81-deg, and 20.96-deg. The shock wave structure, vortex structure, variable distributions, and wall separation topology of the interaction are computed. The results show that in 2D/3D interactions, a new Mach reflection-like event occurs and a Mach stem-like structure is generated above the front of the separation bubble, which differs from the model of Babinsky for 2D planar shock wave/boundary layer interaction. A new interaction model is established to describe this behavior. The relationship between the length of the circumferentially unseparated region in the tube and the AOA of the cone indicates the existence of a critical AOA at which the length is zero, and a prediction of this angle is obtained using an empirical fit, which is verified by computation. The occurrence of side overflow in the windward meridional plane is analyzed, and a quantitative knowledge is obtained. To elucidate the characteristics of the 3D interaction, the scale and structure of the vortex and the pressure and friction force distributions are presented and compared with those of the 2D interaction

Improvements to enhance robustness of third-order scale-independent WENO-Z schemes

Although there are many improvements to WENO3-Z that target the achievement of optimal order in the occurrence of the first-order critical point (CP1), they mainly address resolution performance, while the robustness of schemes is of less concern and lacks understanding accordingly. In light of our analysis considering the occurrence of critical points within grid intervals, we theoretically prove that it is impossible for a scale-independent scheme that has the stencil of WENO3-Z to fulfill the above order achievement, and current scale-dependent improvements barely fulfill the job when CP1 occurs at the middle of the grid cell. In order to achieve scale-independent improvements, we devise new smoothness indicators that increase the error order from 2 to 4 when CP1 occurs and perform more stably. Meanwhile, we construct a new global smoothness indicator that increases the error order from 4 to 5 similarly, through which new nonlinear weights with regard to WENO3-Z are derived and new scale-independents improvements, namely WENO-ZES2 and -ZES3, are acquired. Through 1D scalar and Euler tests, as well as 2D computations, in comparison with typical scale-dependent improvement, the following performances of the proposed schemes are demonstrated: The schemes can achieve third-order accuracy at CP1 no matter its location in the stencil, indicate high resolution in resolving flow subtleties, and manifest strong robustness in hypersonic simulations (e.g., the accomplishment of computations on hypersonic half-cylinder flow with Mach numbers reaching 16 and 19, respectively, as well as essentially non-oscillatory solutions of inviscid sharp double cone flow at M=9.59), which contrasts the comparative WENO3-Z improvement