Communications on Applied Mathematics and Computation ›› 2024, Vol. 6 ›› Issue (1): 533-549.doi: 10.1007/s42967-023-00271-z

• ORIGINAL PAPERS • Previous Articles     Next Articles

Verification and Validation of High-Resolution Inviscid and Viscous Conical Nozzle Flows

Luciano K. Araki1, Rafael B. de R. Borges2, Nicholas Dicati P. da Silva3, Chi-Wang Shu4   

  1. 1. Department of Mechanical Engineering, Federal University of Paraná, Curitiba, 81531-980, Paraná, Brazil;
    2. Mathematics and Statistics Institute, Rio de Janeiro State University, Rio de Janeiro, 20550-900, Rio de Janeiro, Brazil;
    3. Department of Mechanical Engineering, State University of Maringá, Maringá, 87020-900, Paraná, Brazil;
    4. Division of Applied Mathematics, Brown University, Providence, RI, 02912, USA
  • Received:2023-01-30 Revised:2023-03-04 Published:2024-04-16
  • Contact: Nicholas Dicati P. da Silva,E-mail:ndicati@gmail.com;Luciano K. Araki,E-mail:lucaraki@ufpr.br;Rafael B. de R. Borges,E-mail:rafael.borges@ime.uerj.br;Chi-Wang Shu,E-mail:chi-wang_shu@brown.edu E-mail:ndicati@gmail.com;lucaraki@ufpr.br;rafael.borges@ime.uerj.br;chi-wang_shu@brown.edu
  • Supported by:
    Research of Chi-Wang Shu is partially supported by the AFOSR grant FA9550-20-1-0055 and the NSF grant DMS-2010107.

Abstract: Capturing elaborated flow structures and phenomena is required for well-solved numerical flows. The finite difference methods allow simple discretization of mesh and model equations. However, they need simpler meshes, e.g., rectangular. The inverse Lax-Wendroff (ILW) procedure can handle complex geometries for rectangular meshes. High-resolution and high-order methods can capture elaborated flow structures and phenomena. They also have strong mathematical and physical backgrounds, such as positivity-preserving, jump conditions, and wave propagation concepts. We perceive an effort toward direct numerical simulation, for instance, regarding weighted essentially non-oscillatory (WENO) schemes. Thus, we propose to solve a challenging engineering application without turbulence models. We aim to verify and validate recent high-resolution and high-order methods. To check the solver accuracy, we solved vortex and Couette flows. Then, we solved inviscid and viscous nozzle flows for a conical profile. We employed the finite difference method, positivity-preserving Lax-Friedrichs splitting, high-resolution viscous terms discretization, fifth-order multi-resolution WENO, ILW, and third-order strong stability preserving Runge-Kutta. We showed the solver is high-order and captured elaborated flow structures and phenomena. One can see oblique shocks in both nozzle flows. In the viscous flow, we also captured a free-shock separation, recirculation, entrainment region, Mach disk, and the diamond-shaped pattern of nozzle flows.

Key words: High-resolution, Compressible, Navier-Stokes, Free-shock separation, Nozzle flow