Numerical Simulation of Shock Wave Reflecting Patterns for Different Flow Conditions

  • Choi, Sung-Yoon (Department of Aerospace Engineering Pusan National University) ;
  • Oh, Se-Jong (Department of Aerospace Engineering Pusan National University)
  • Published : 2002.05.30


The numerical experiment has been conducted to investigate the unsteady shock wave reflecting phenomena. The cell-vertex finite-volume, Roe's upwind flux difference splitting method with unstructured grid is implemented to solve unsteady Euler equations. The $4^{th}$-order Runge-Kutta method is applied for time integration. A linear reconstruction of the flux vector using the least-square method is applied to obtain the $2^{nd}$-order accuracy for the spatial derivatives. For a better resolution of the shock wave and slipline, the dynamic grid adaptation technique is adopted. The new concept of grid adaptation technique, which is much simpler than that of conventional techniques, is introduced for the current study. Three error indicators (divergence and curl of velocity, and gradient of density) are used for the grid adaptation procedure. Considering the quality of the solution and the numerical efficiency, the grid adaptation procedure was updated up to $2^{nd}$ level at every 20 time steps. For the convenience of comparison with other experimental and analytical results, the case of interaction between the straight incoming shock wave and a sharp wedge is simulated for various flow conditions. The numerical results show good agreement with other experimental and analytical results, in the shock wave reflecting structure, slipline, and the trajectory of the triple points. Some critical cases show disagreement with the analytical results, but these cases also have been proven to show hysteresis phenomena.


  1. Ben-Dar, Gabi, Shock Wave Reflection Phenomena, Springer-Verlag Inc., New York, 1992.
  2. Baum, J. D. and Lohner, R., "Numerical Simulation of Shock-Box Interaction Using an Adaptive Finite Element Scheme," AIAA Journal, Vol. 32, No.4, 1994, pp. 682-692.
  3. Choi, W. S., "Numerical analysis of steady, unsteady compressible flow fields using finite volume TVD schemes," Ph. D. Thesis, Postech, 1994.
  4. Yang, H., "Impingement of a Moving Shock Wave into a Cavity," Master Thesis, KAIST, 1998.
  5. Lohner, R.and Baum, J. D., "Numerical Simulation of shock Interaction with Complex Geometry Three-Dimensional Structures Using A New Adaptive H-Refinement Scheme on Unstructured Grids," AlAA paper 90-0700, 1990.
  6. Ivanov, M. S., Gimelshein, S. F. and Beylich, A. E., "Hysteresis effect in stationary reflection of shock waves," Phys. Fluids, Vol. 7, No.4, 1995, pp. 685-687.
  7. Anderson, J. D., Jr ., Modem Compressible Flow, McGraw-Hill Publishing Company, 1990.
  8. Roe, P. L., "Characteristic Based Schemes for the Euler Equations," Annual Review of Fluid Mechanics, Vol. 18, 1986, pp. 337-365.
  9. Fink, N. T., "Upwind Scheme for Solving the Euler Equations on Unstructured Tetrahedral Meshes," AIAA Journal, Vol. 30, No. 1, 1992, pp. 70-77.
  10. Anderson, W. K. and Lhaus, D., "An Implicit Upwind Algorithm Computing Turbulent Flows In Unstructured Grid," Computers and Fluids, Vol. 23, No. 1, 1994, pp. 1-21.
  11. Barth. T. J. and Jespersen, D. C, "The Design and Application of Upwind Scheme on Unstructured Meshes," AIAA paper 89-0366, 1989.
  12. Kozic, M. S. and Sredojevic, D., "Development of unstructured dynamic grids for solving unsteady two-dimensional Euler equations," The Aeronautical Journal, Vol. 102, No. 1014, 1998, pp. 195-200.
  13. Lohner R., "AN ADAPTIVE FINITE ELEMENT SCHEME FOR TRANSIENT PROBLEMS IN CFD," Computer Methods in Applied Mechanics and Engineering, Vol. 61, 1987, pp. 323-338.
  14. Pizadeh, S. Z.,"An Adaptive Unstructured Grid Method by Grid Subdivision, Local Remeshing, and Grid Movement," AIAA paper 99-3255, 1999.
  15. Ivanov, M. S., et. aI, "Numerical Analysis of Shock Wave Reflection Transition in Steady Flows," AIAA Journal, Vol. 36, No. 11, 1998, pp. 2079-2086.
  16. Cho, K. H., "The Optimal shape design using adjoint method on unstructured meshes," Master Thesis, Pusan National University, 2000.
  17. Kim, S. H., "Numerical Simulation of Shock_Induced Combustion on Adaptive Unstructured Grid," Master Thesis, Pusan National University, 2001.