Updates to Multi-Dimensional Flux Reconstruction for Hypersonic Simulations on Tetrahedral Grids

The quality of simulated hypersonic stagnation region heating with tetrahedral meshes is investigated by using an updated three-dimensional, upwind reconstruction algorithm for the inviscid flux vector. An earlier implementation of this algorithm provided improved symmetry characteristics on tetrahedral grids compared to conventional reconstruction methods. The original formulation however displayed quantitative differences in heating and shear that were as large as 25% compared to a benchmark, structured- grid solution. The primary cause of this discrepancy is found to be an inherent inconsistency in the formulation of the flux limiter. The inconsistency is removed by employing a GreenGauss formulation of primitive gradients at nodes to replace the previous Gram-Schmidt algorithm. Current results are now in good agreement with benchmark solutions for two challenge problems: (1) hypersonic flow over a three-dimensional cylindrical section with special attention to the uniformity of the solution in the spanwise direction and (2) hypersonic flow over a three-dimensional sphere. The tetrahedral cells used in the simulation are derived from a structured grid where cell faces are bisected across the diagonal resulting in a consistent pattern of diagonals running in a biased direction across the otherwise symmetric domain. This grid is known to accentuate problems in both shock capturing and stagnation region heating encountered with conventional, quasi-one-dimensional inviscid flux reconstruction algorithms. Therefore the test problems provide a sensitive indicator for algorithmic effects on heating. Additional simulations on a sharp, double cone and the shuttle orbiter are then presented to demonstrate the capabilities of the new algorithm on more geometrically complex flows with tetrahedral grids. These results provide the first indication that pure tetrahedral elements utilizing the updated, three-dimensional, upwind reconstruction algorithm may be used for the simulation of heating and shear in hypersonic flows in upwind, finite volume formulations.

[1]  Bernard Grossman,et al.  Characteristic-based, rotated upwind scheme for the Euler equations , 1992 .

[2]  Pramod K. Subbareddy,et al.  Unstructured grid approaches for accurate aeroheating simulations , 2007 .

[3]  Peter A. Gnoffo,et al.  CFD Validation Studies for Hypersonic Flow Prediction , 2001 .

[4]  Peter A. Gnoffo,et al.  Conservation equations and physical models for hypersonic air flows in thermal and chemical nonequilibrium , 1989 .

[5]  Peter A. Gnoffo,et al.  Computational Aerothermodynamic Simulation Issues on Unstructured Grids , 2004 .

[6]  Peter A. Gnoffo,et al.  Multi-Dimensional, Inviscid Flux Reconstruction for Simulation of Hypersonic Heating on Tetrahedral Grids , 2009 .

[7]  Garrett. Barter,et al.  Shock Capturing with PDE-Based Artificial Viscosity for an Adaptive, Higher-Order Discontinuous Galerkin Finite Element Method , 2008 .

[8]  Richard A. Thompson,et al.  Computational Aerothermodynamic Analysis for the STS-107 Accident Investigation , 2004 .

[9]  Peter A. Gnoffo.  Conservation equations and physical models for hypersonic air flows over the aeroassist flight experiment vehicle , 1989 .

[10]  David A. Stewart,et al.  Surface Catalysis and Characterization of Proposed Candidate TPS for Access-to-Space Vehicles , 1997 .

[11]  Timothy J. Barth,et al.  The design and application of upwind schemes on unstructured meshes , 1989 .

[12]  Peter A. Gnoffo,et al.  Computational Fluid Dynamics Technology for Hypersonic Applications , 2003 .

[13]  Graham V. Candler,et al.  Development of a hybrid unstructured implicit solver for the simulation of reacting flows over complex geometries , 2004 .

[14]  Peter A. Gnoffo,et al.  Simulation of Stagnation Region Heating in Hypersonic Flow on Tetrahedral Grids , 2007 .

[15]  Graham V. Candler,et al.  CFD Validation for Hypersonic Flight: Hypersonic Double-Cone Flow Simulations , 2002 .

[16]  W. K. Anderson,et al.  An implicit upwind algorithm for computing turbulent flows on unstructured grids , 1994 .

[17]  C. Fletcher,et al.  Galerkin Finite-Element Methods , 1984 .

[18]  Graham V. Candler,et al.  Effect of Vibrational Nonequilibrium on Hypersonic Double-Cone Experiments , 2003 .

[19]  J. A. White,et al.  A Psuedo-Temporal Multi-Grid Relaxation Scheme for Solving the Parabolized Navier-Stokes Equations , 1999 .

[20]  Cheatwood F. McNeil,et al.  User''s Manual for the Langley Aerothermodynamic Upwind Relaxation Algorithm (LAURA) , 1996 .

[21]  Peter A. Gnoffo,et al.  Multiblock analysis for Shuttle Orbiter reentry heating from Mach 24 to Mach 12 , 1994 .

[22]  G. Candler,et al.  Data-Parallel Line Relaxation Method for the Navier -Stokes Equations , 1998 .

[23]  Alireza Mazaheri,et al.  Laura Users Manual: 5.1-41601 , 2009 .

[24]  Bernard Grossman,et al.  A rotated upwind scheme for the Euler equations , 1991 .

[25]  Peter A. Gnoffo,et al.  A multiblock analysis for shuttle orbiter re-entry heating from Mach24 to Mach 12 , 1993 .

[26]  D. McRae,et al.  Rotated upwind strategies for solution of the Euler equations , 1994 .