Exploiting active subspaces to quantify uncertainty in the numerical simulation of the HyShot II scramjet

Abstract We present a computational analysis of the reactive flow in a hypersonic scramjet engine with focus on effects of uncertainties in the operating conditions. We employ a novel methodology based on active subspaces to characterize the effects of the input uncertainty on the scramjet performance. The active subspace identifies one-dimensional structure in the map from simulation inputs to quantity of interest that allows us to reparameterize the operating conditions; instead of seven physical parameters, we can use a single derived active variable. This dimension reduction enables otherwise infeasible uncertainty quantification, considering the simulation cost of roughly 9500 CPU-hours per run. For two values of the fuel injection rate, we use a total of 68 simulations to (i) identify the parameters that contribute the most to the variation in the output quantity of interest, (ii) estimate upper and lower bounds on the quantity of interest, (iii) classify sets of operating conditions as safe or unsafe corresponding to a threshold on the output quantity of interest, and (iv) estimate a cumulative distribution function for the quantity of interest.

[1]  M. Urner.  Scattered Data Approximation , 2016 .

[2]  Datta V. Gaitonde,et al.  Exploratory Simulations of the HIFiRE 2 Scramjet Flowpath , 2012 .

[3]  P. C. Gehlen,et al.  Computer Experiments , 1996 .

[4]  F. Billig,et al.  Prediction of Precombustion Wall Pressure Distributions in Scramjet Engines , 1973 .

[5]  G. Iaccarino,et al.  Risk assessment of scramjet unstart using adjoint-based sampling methods , 2012 .

[6]  Johan Steelant,et al.  Ground Testing of the HyShot Supersonic Combustion Flight Experiment in HEG and Comparison with Flight Data (AIAA 2004–3345) , 2004 .

[7]  Gianluca Iaccarino,et al.  Characterizing the operability limits of the HyShot II scramjet through RANS simulations , 2011 .

[8]  Klaus Hannemann,et al.  Investigation of Unsteady/Quasi-Steady Scramjet Behavior using High-Speed Visualization Techniques , 2012 .

[9]  Klaus Hannemann,et al.  Experimental and Numerical Investigation of the HyShot II Flight Experiment , 2013 .

[10]  Vincent Terrapon.  Flamelet-based model for supersonic combustion , 2009 .

[11]  Anastasios S. Lyrintzis,et al.  Assessment of Turbulent Shock-Boundary Layer Interaction Computations Using the OVERFLOW Code , 2007 .

[12]  Karen Willcox,et al.  A Survey of Projection-Based Model Reduction Methods for Parametric Dynamical Systems , 2015, SIAM Rev..

[13]  Russell R. Boyce,et al.  The HyShot Scramjet Flight Experiment - Flight Data and CFD Calculations Compared , 2003 .

[14]  G. Iaccarino,et al.  Non-intrusive low-rank separated approximation of high-dimensional stochastic models , 2012, 1210.1532.

[15]  R. H. Moore,et al.  Regression Graphics: Ideas for Studying Regressions Through Graphics , 1998, Technometrics.

[16]  Jeroen A. S. Witteveen,et al.  Uncertainty quantification for laminar-turbulent transition prediction in RANS turbomachinery applications , 2011 .

[17]  Johan Steelant,et al.  Ground Testing of the HyShot II Scramjet Configuration in HEG , 2008 .

[18]  H. Bungartz,et al.  Sparse grids , 2004, Acta Numerica.

[19]  Chih-Peng Chen,et al.  Numerical Investigation of Supersonic Combustion of the HyShot II in the Shock Tunnel , 2011 .

[20]  F. Menter.  ZONAL TWO EQUATION k-w TURBULENCE MODELS FOR AERODYNAMIC FLOWS , 1993 .

[21]  Ralph C. Smith,et al.  Uncertainty Quantification: Theory, Implementation, and Applications , 2013 .

[22]  Klaus Hannemann,et al.  Incipient thermal choking and stable shock-train formation in the heat-release region of a scramjet combustor. Part I: Shock-tunnel experiments , 2015 .

[23]  Sebastian Karl,et al.  Incipient thermal choking and stable shock-train formation in the heat-release region of a scramjet combustor. Part II: Large eddy simulations , 2015 .

[24]  Donald R. Jones,et al.  A Taxonomy of Global Optimization Methods Based on Response Surfaces , 2001, J. Glob. Optim..

[25]  Gianluca Iaccarino,et al.  Modeling of structural uncertainties in Reynolds-averaged Navier-Stokes closures , 2013 .

[26]  Jeroen A. S. Witteveen,et al.  Assessment of Uncertainties in Modeling of Laminar to Turbulent Transition for Transonic Flows , 2013 .

[27]  Johan Steelant,et al.  Ground testing of the HyShot supersonic combustion flight experiment in HEG , 2005 .

[28]  Dongbin Xiu,et al.  The Wiener-Askey Polynomial Chaos for Stochastic Differential Equations , 2002, SIAM J. Sci. Comput..

[29]  R. Ghanem,et al.  Stochastic Finite Elements: A Spectral Approach , 1990 .

[30]  Klaus Hannemann,et al.  A Computational Study of the HySHOT II Combustor Performance , 2013 .

[31]  Anthony Donald Gardner.  HyShot scramjet testing in the HEG , 2007 .

[32]  Christine A. Shoemaker,et al.  ORBIT: Optimization by Radial Basis Function Interpolation in Trust-Regions , 2008, SIAM J. Sci. Comput..

[33]  Johan Steelant,et al.  Methodology of a combined ground based testing and numerical modelling analysis of supersonic combustion flow paths , 2010 .

[34]  Shuangzhe Liu,et al.  Global Sensitivity Analysis: The Primer by Andrea Saltelli, Marco Ratto, Terry Andres, Francesca Campolongo, Jessica Cariboni, Debora Gatelli, Michaela Saisana, Stefano Tarantola , 2008 .

[35]  Florian R. Menter,et al.  Transition Modelling for General Purpose CFD Codes , 2006 .

[36]  Gianluca Iaccarino,et al.  Quantification of margins and uncertainties using an active subspace method for approximating bounds , 2014 .

[37]  E. Jaynes.  Information Theory and Statistical Mechanics , 1957 .

[38]  Michael K. Smart,et al.  Flight Data Analysis of the HyShot 2 Scramjet Flight Experiment , 2006 .

[39]  Heinz Pitsch,et al.  Full-system RANS of the HyShot II scramjet Part 2: Reactive cases , 2010 .

[40]  Ronald K. Hanson,et al.  An improved H2/O2 mechanism based on recent shock tube/laser absorption measurements , 2011 .

[41]  Jeffrey Larson,et al.  Derivative-Free Optimization of Expensive Functions with Computational Error Using Weighted Regression , 2013, SIAM J. Optim..

[42]  Dongbin Xiu,et al.  High-Order Collocation Methods for Differential Equations with Random Inputs , 2005, SIAM J. Sci. Comput..

[43]  Andreas Mack,et al.  CFD Analysis of the HyShot II Scramjet Experiments in the HEG Shock Tunnel , 2008 .

[44]  Klaus Hannemann,et al.  CFD Analysis of Unstart Characteristics of the HyShot II Scramjet Configuration in the HEG Shock Tunnel , 2011 .

[45]  Fabio Nobile,et al.  Analysis of Discrete $$L^2$$L2 Projection on Polynomial Spaces with Random Evaluations , 2014, Found. Comput. Math..

[46]  Malik Beshir Malik,et al.  Applied Linear Regression , 2005, Technometrics.

[47]  Klaus Hannemann,et al.  An Experimental Investigation of Steady and Unsteady Combustion Phenomena in the HyShot II Combustor , 2011 .

[48]  Eric R. Ziegel,et al.  The Elements of Statistical Learning , 2003, Technometrics.

[49]  Saltelli Andrea,et al.  Global Sensitivity Analysis: The Primer , 2008 .

[50]  John L. Nazareth,et al.  Introduction to derivative-free optimization , 2010, Math. Comput..

[51]  Heinz Pitsch,et al.  An Ecient Flamelet-based Combustion Model for Supersonic Flows , 2011 .

[52]  D. Gleich,et al.  Computing active subspaces with Monte Carlo , 2014, 1408.0545.

[53]  Qiqi Wang,et al.  Erratum: Active Subspace Methods in Theory and Practice: Applications to Kriging Surfaces , 2013, SIAM J. Sci. Comput..

[54]  Jiyuan Tu,et al.  Computational Fluid Dynamics: A Practical Approach , 2007 .

[55]  Gianluca Iaccarino,et al.  Full-system RANS of the HyShot II scramjet Part 1: Numerics and non-reactive simulations , 2010 .

[56]  M. Kenward,et al.  An Introduction to the Bootstrap , 2007 .

[57]  Matthew Doolan,et al.  Analysis of Transient Thermal Choking Processes in a Model Scramjet Engine , 1998 .

[58]  Michael K. Smart,et al.  Flight Data Analysis of HyShot 2 , 2005 .

[59]  Houman Owhadi,et al.  A non-adapted sparse approximation of PDEs with stochastic inputs , 2010, J. Comput. Phys..

[60]  Christer Fureby,et al.  CFD Analysis of the HyShot II scramjet combustor , 2011 .

[61]  Heinz Pitsch,et al.  Reynolds-Averaged Navier-Stokes Simulations of the HyShot II Scramjet , 2012 .

[62]  Klaus Hannemann,et al.  Experimental investigation of different scramjet hydrogen injection systems , 2008 .

[63]  J. Wagner,et al.  Experimental Investigation of Unstart in an Inlet/Isolator Model in Mach 5 Flow , 2009 .

[64]  Carl E. Rasmussen,et al.  Gaussian processes for machine learning , 2005, Adaptive computation and machine learning.

[65]  Eric R. Ziegel,et al.  Handbook of Statistics 13: Design and Analysis of Experiments , 2000 .

[66]  Ker-Chau Li,et al.  Regression Analysis Under Link Violation , 1989 .

[67]  Peter A. Jacobs,et al.  Fuel Stagnation Temperature Effects on Mixing with Supersonic Combustion Flows , 1997 .

[68]  F. Menter.  Two-equation eddy-viscosity turbulence models for engineering applications , 1994 .

[69]  Heinz Pitsch,et al.  An efficient flamelet-based combustion model for compressible flows , 2015 .

[70]  Krishnan Mahesh,et al.  Modeling shock unsteadiness in shock/turbulence interaction , 2003 .