Capturing of Shock Wave of Supersonic Flow over the Bump Channel with TVD, ACM and Jameson Methods

Document Type : Research Paper


1 Corresponding Author, PhD, Department of Mechanical Engineering, K.N.Toosi University, Tehran, Iran

2 Assistant Professor, Department of Mechanical Engineering, Azarbaijan Shahid Madani University, Tabriz, Iran


The phenomenon of capturing of sharp gradients, shock waves, and contact discontinuities in compressible aerodynamic flows is the major concern while numerical methods based on the characteristic variables are favorable solutions. The aim of the present study is to introduce a method based on the characteristic variables (Riemann solution) and control the diffusion term in the classic methods in order to capture the shock waves. An efficient blending procedure based on the density-based algorithm is presented to solve the compressible Euler equations, using a non-orthogonal mesh with collocated finite volume formulation. The fluxes of the convected quantities are approximated by employing the characteristic based total variation diminishing (TVD), artificial compressibility method (ACM) and Jameson methods. Results show that the ACM and TVD methods captured the shock waves with higher resolution than the Jameson method. Moreover, not only quality of shock wave capturing was improved for all flows at the discontinuities by employing ACM and TVD methods, but also the computational time and convergence were improved for the supersonic flows.


Main Subjects

[1] Harten, A., "High Resolution Scheme for Hyperbolic Conservation Laws", Journal of
Computational Physics, Vol. 49, No. 3, pp. 357-393, (1983).
[2] Jameson, A., Schmidt, W., and Turkel, E., "Numerical Solutions of the Euler Equations by
Finite-Volume Methods using Runge–kutta Time-stepping Schemes", AIAA, pp. 81-125,
[3] Harten, A., "The Artificial Compression Method for Computation of Shocks and Contact
Discontinuities. I. Single Conservation Lows", Comm. Pure & Appl., pp. 611-637, (1977).
[4] Harten, A., "The Artificial Compression Method for Computation of Shocks and Contact
Discontinuities: III. Self-Adjusting Hybrid Schemes", Math. Comput, Vol. 32, No. 142, pp.
363-389, (1978).
[5] Yee, H., Warming, R.F., and Harten, A., "Implicit Total Variation Diminishing (TVD)
Schemes for Steady State Calculations", J. Comput. Phys, Vol. 57, No. 2, pp. 327-360,
[6] Yee, H.C., Sandham, N.D., and Djomeri, M.J., "Low Dissipative High-order Shock-
Capturing Methods using Characteristic-Based Filters", J. Comput. Phys, Vol. 150, pp.
199-238, (1999).
[7] Lie, K.A., and Noelle, S., "On the Artificial Compression Method for Second-Order Non
Oscillatory Central Difference Schemes for Systems of Conservation Laws", SIAM J. Sci.
Comput, Vol. 24, No. 4, pp. 1157-1174, (2003).
[8] Yadegari, M., and Abdollahi Jahdi, M.H., "Shock Capturing Method by Numerical
Dissipation Control on Symmetric Airfoil", Journal of Solid and Fluid Mechanics, Vol. 6,
pp. 285-304, (2016).
[9] Lien, F.S., and Leschziner, M.A.," A Pressure–velocity Solution Strategy for Compressible
Flow and Its Application to Shock Boundary Layer Interaction using Second-moment
Turbulence Closure", J. Fluids Eng, Vol. 115, No. 4, pp. 717-725, (1993).
[10] Van Leer, B., "Towards the Ultimate Conservation Difference Scheme. II. Monotonicity
and Conservation Combined in Second-order Scheme", J. Comput. Phys, Vol. 14, No. 4,
pp. 361-370, (1974).
[11] Shyy, W., and Thakur, S., "Controlled Variation Scheme in a Sequential Solver for
Recirculating Flows", Heat Transf, Vol. 25, pp. 273-286, (1994).
[12] Thakur, S., Wright, J., Shyy, W., Liu, J., Ouyang, H., and Vu, T., "Development of
Pressure-based Composite Multigrid Methods for Complex Fluid Flows", Prog,
Aerospace Sci, Vol. 32, No. 4, pp. 313-375, (1996).
[13] Mulder, W.A., and Van Leer, B., "Experiments with Implicit Upwind Methods for the
Euler Equations", J. Comput. Phys, Vol. 59, pp. 232-246, (1985).
[14] Lin, H., and Chieng, C.C., "Characteristic-Based Flux Limiters of an Essentially Third-
Order Flux-Splitting Method for Hyperbolic Conservation Laws", Int. J. Numer. Methods
Fluids, Vol. 13, No. 3, pp. 287-307, (1991).
[15] Issa, R.I., and Javareshkian, M.H., "Application of TVD Schemes in Pressure-based Finitevolume
Methods", In: Proceedings of the Fluids Engineering, (1996).
[16] Issa, R.I., and Javareshkian, M.H., "Pressure-based Compressible Calculation Method
Utilizing Total Variation Diminishing Schemes", AIAA J. Vol. 36, No. 9, pp. 1652-1657,
[17] Kobayashi, M.H., and Pereira, J.C.F., "Characteristic-based Pressure Correction at All
Speeds", AIAA J., Vol. 34, No. 2, pp. 272-280, (1996).
[18] Turkel, E., "Preconditioned Methods for Solving the Incompressible and Low Speed
Compressible Equations", J. Comput. Phys, Vol. 72, pp. 277-298, (1987).
[19] Turkel, E., Radespiel, R., and Kroll, N., "Assessment of Preconditioning Methods for
Multidimensional Aerodynamics", Comp. Fluids, Vol. 26, No. 6, pp. 613-634, (1997).
[20] Guillard, H., and Viozat, C., "On the Behavior of Upwind Schemes in the Low Mach
Number Limit", Comp. Fluids, Vol. 28, No. 1, pp. 63-86, (1999).
[21] Zamzamian, K., and Razavi, S.E., "Multidimensional Upwinding for Incompressible
Flows Based on Characteristics", J. Comput. Phys, Vol. 227, No. 19, pp. 8699-8713,
[22] Godunov, S. K.A., "Difference Scheme for Numerical Computation of Discontinuous
Solutions of Hydrodynamic Equations", MathSbornic, English Translation in
U.Sjointpublications, Vol . 47, pp. 271-306, (1959).
[23] Roe, P. L., "Approximate Riemann Solvers, Parameter Vectors and Difference Schemes",
Journal of Computational Physics, Vol. 43, pp. 357-372, (1981).
[24] Colella, P., and Woodward, P. R., "The Piecewise Parabolic Method (PPM) for Gas
Dynamical Simulations", J. of Comput. Phys, Vol. 54, pp. 174-201, (1984).