Validation of Under-Resolved Numerical Simulations of the PDC Exhaust Flow Based on High Speed Schlieren

Owing to their high thermodynamic efficiency, pulsating combustion cycles have become an attractive option for future gas turbine designs. Yet, their potential gains should not be outweighed by losses due to unsteady pressure wave interactions between engine components. Consequently, the geometric engine design moves into focus. Ideally, one would quickly test several different principal layouts with respect to their qualitative behavior, select the most promising variants and then move on to detailed optimization. Computational fluid dynamics (CFD) appears as the methodology of choice for such preparatory testing. Yet, the inevitable geometric complexity of such engines makes fully resolved CFD an arduous and expensive task necessitating computations on top high-performance hardware, even with modern adaptive mesh refinement in place. In the present work we look at the initial flow field of a shock generated by a pulse detonation combustor (PDC) which leaves the combustion chamber and enters the plenum. We provide first indicators, however, that overall mechanical loads, represented by large-scale means of, e.g., mass, energy, and momentum fluxes can be well estimated on the basis of rather coarsely resolved CFD calculations. Comparing high-resolution simulations of the exit of a strong shock from a combustion tube with experimental schlieren photographs, we first establish validity of fully resolved CFD. Next we compare several integral quantities representative of overall mechanical loads with a sequence of successively coarser grid simulations, thereby corroborating that the “quick and dirty” coarse-grained simulations indeed allow for good order of magnitude estimates.

[1]  E. Toro Riemann Solvers and Numerical Methods for Fluid Dynamics , 1997 .

[2]  R. E. Carlson,et al.  Monotone Piecewise Cubic Interpolation , 1980 .

[3]  Bram van Leer,et al.  On the Relation Between the Upwind-Differencing Schemes of Godunov, Engquist–Osher and Roe , 1984 .

[4]  Wolfgang Schröder,et al.  A cut-cell method for sharp moving boundaries in Cartesian grids , 2013 .

[5]  P. Colella,et al.  Local adaptive mesh refinement for shock hydrodynamics , 1989 .

[6]  G. Strang On the Construction and Comparison of Difference Schemes , 1968 .

[7]  N. Nikiforakis,et al.  Well-balanced compressible cut-cell simulation of atmospheric flow , 2009, Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[8]  M. Meinke,et al.  A flexible level-set approach for tracking multiple interacting interfaces in embedded boundary methods , 2014 .

[9]  Julio Soria,et al.  An assessment of high-power light-emitting diodes for high frame rate schlieren imaging , 2012 .

[10]  Marsha J. Berger,et al.  A Simplified h-box Method for Embedded Boundary Grids , 2012, SIAM J. Sci. Comput..

[11]  C. O. Paschereit,et al.  A compact shock-focusing geometry for detonation initiation: Experiments and adjoint-based variational data assimilation , 2017 .

[12]  B. W. Skews,et al.  The shape of a diffracting shock wave , 1967, Journal of Fluid Mechanics.

[13]  Phillip Colella,et al.  A Cartesian grid embedded boundary method for hyperbolic conservation laws , 2006 .

[14]  Nikolaus A. Adams,et al.  A conservative interface method for compressible flows , 2006, J. Comput. Phys..