A high-order cut-cell method for numerical simulation of hypersonic boundary-layer instability with surface roughness

Laminar-turbulent transition of hypersonic boundary layers can be affected significantly by the existence of surface roughness. Currently many important mechanisms of roughness-induced transition are not well understood. In recent years, direct numerical simulation (DNS) has been extensively applied for investigating instability and transition mechanisms of hypersonic boundary layers. Most of the past DNS studies, however, have been based on body-fitted grids for smooth surfaces without roughness. Due to complex geometry of embedded roughness, the use of body-fitted grids can be very difficult for flow with arbitrary surface roughness. In this paper, we present a new high-order cut-cell method to overcome the natural complexities in grid generation around arbitrary surface of roughness. The new method combines a non-uniform-grid finite-difference method for discrete grid points near the solid boundary and a shock-fitting method for the treatment of the bow shock. The non-uniform-grid finite-difference formulas are expressed in a general explicit form so that they can be applied to different multi-dimensional problems without any modification. The computational accuracy of new algorithms of up to O(h^4) are tested on several one- and two-dimensional elliptic equations in irregular domains. In addition, the new method is applied to the simulation of the receptivity process of a Mach 5.92 flow over a flat plate under the combined effect of an isolated surface roughness element and surface blow and suction. A good agreement is found between the unsteady flow results and those obtained by a Linear Stability Theory (LST).

[1]  W. Shyy,et al.  Computation of Solid-Liquid Phase Fronts in the Sharp Interface Limit on Fixed Grids , 1999 .

[2]  James P. Keener,et al.  Immersed Interface Methods for Neumann and Related Problems in Two and Three Dimensions , 2000, SIAM J. Sci. Comput..

[3]  Paul Fischer,et al.  Roughness-Induced Transient Growth , 2005 .

[4]  Murshed Hossain,et al.  The application of spectral methods in simulating compressible fluid and magnetofluid turbulence , 1993 .

[5]  Gianluca Iaccarino,et al.  Numerical simulation of the effect of a roughness element on high-speed boundary-layer instability , 2008 .

[6]  M. D. Salas,et al.  Spectral methods for the Euler equations. II - Chebyshev methods and shock fitting , 1985 .

[7]  Ronald Fedkiw,et al.  Regular Article: The Ghost Fluid Method for Deflagration and Detonation Discontinuities , 1999 .

[8]  Xiaolin Zhong,et al.  Direct Numerical Simulation and the Theory of Receptivity in a Hypersonic Boundary Layer , 2006 .

[9]  Nikolaus A. Adams,et al.  A High-Resolution Hybrid Compact-ENO Scheme for Shock-Turbulence Interaction Problems , 1996 .

[10]  R. LeVeque,et al.  A comparison of the extended finite element method with the immersed interface method for elliptic equations with discontinuous coefficients and singular sources , 2006 .

[11]  Derek M. Causon,et al.  A Cartesian cut cell method for shallow water flows with moving boundaries , 2001 .

[12]  M. D. Salas,et al.  Spectral methods for the Euler equations , 1983 .

[13]  L. Mack.  Boundary-Layer Linear Stability Theory , 1984 .

[14]  Zhilin Li,et al.  Convergence analysis of the immersed interface method , 1999 .

[15]  J. Quirk.  An alternative to unstructured grids for computing gas dynamic flows around arbitrarily complex two-dimensional bodies , 1994 .

[16]  Xiaolin Zhong,et al.  Numerical Simulation and Experiment Comparison of Leading-Edge Receptivity of A Mach 5.92 Boundary Layer , 2006 .

[17]  S. Berry,et al.  Discrete Roughness Transition for Hypersonic Flight Vehicles , 2007 .

[18]  Xiaolin Zhong,et al.  High-Order Finite-Difference Schemes for Numerical Simulation of Hypersonic Boundary-Layer Transition , 1998 .

[19]  Xiaolin Zhong,et al.  Receptivity of a supersonic boundary layer over a flat plate. Part 2. Receptivity to free-stream sound , 2003, Journal of Fluid Mechanics.

[20]  Andreas Wiegmann,et al.  The Explicit-Jump Immersed Interface Method: Finite Difference Methods for PDEs with Piecewise Smooth Solutions , 2000, SIAM J. Numer. Anal..

[21]  T. A. Zang,et al.  Spectral methods for simulations of transition and turbulence , 1990 .

[22]  Peter McCorquodale,et al.  A Cartesian grid embedded boundary method for the heat equation on irregular domains , 2001 .

[23]  Miguel R. Visbal,et al.  Further development of a Navier-Stokes solution procedure based on higher-order formulas , 1999 .

[24]  Chi-Wang Shu,et al.  Efficient Implementation of Weighted ENO Schemes , 1995 .

[25]  A. D. Gosman,et al.  A new spectral method for calculation of the time-varying area of a laminar flame in homogeneous turbulence , 1991 .

[26]  Nikolaus A. Adams,et al.  Subharmonic transition to turbulence in a flat-plate boundary layer at Mach number 4.5 , 1996, Journal of Fluid Mechanics.

[27]  Ponnampalam Balakumar,et al.  Transition in a Supersonic Boundary-Layer Due to Roughness and Acoustic Disturbances , 2003 .

[28]  Kazufumi Ito,et al.  Three-dimensional elliptic solvers for interface problems and applications , 2003 .

[29]  S. Schneider.  Summary of Hypersonic Boundary-Layer Transition Experiments on Blunt Bodies with Roughness , 2008 .

[30]  Hermann F. Fasel,et al.  Direct numerical simulation of transition in supersonic boundary layers: oblique breakdown , 1993 .

[31]  William S. Saric,et al.  BOUNDARY-LAYER RECEPTIVITY TO FREESTREAM DISTURBANCES , 2002 .

[32]  Alexander N. Shiplyuk,et al.  Hypersonic flow stability experiments , 2002 .

[33]  Randall J. LeVeque,et al.  An Immersed Interface Method for Incompressible Navier-Stokes Equations , 2003, SIAM J. Sci. Comput..

[34]  A. Demetriades,et al.  Hypersonic viscous flow over a slender cone. III - Laminar instability and transition , 1974 .

[35]  Xiaolin Zhong,et al.  Receptivity of a supersonic boundary layer over a flat plate. Part 1. Wave structures and interactions , 2003, Journal of Fluid Mechanics.

[36]  R. J. LeVeque,et al.  ICASE I STABLE BOUNDARY CONDITIONS FOR CARTESIAN GRID CALCULATIONS , 1990 .

[37]  R. Fedkiw,et al.  A fourth order accurate discretization for the Laplace and heat equations on arbitrary domains, with applications to the Stefan problem , 2005 .

[38]  Krishnan Mahesh,et al.  High order finite difference schemes with good spectral resolution , 1997 .

[39]  Yao H. Chu,et al.  A new type of TVD schemes for computations of high speed flows , 1993 .

[40]  Anthony Demetriades.  Laminar Boundary Layer Stability Measurements at Mach 7 Including Wall Temperature Effects , 1977 .

[41]  C. Peskin.  The immersed boundary method , 2002, Acta Numerica.

[42]  S. Lele.  Compact finite difference schemes with spectral-like resolution , 1992 .

[43]  Timothy P. Chartier,et al.  New Geometric Immersed Interface Multigrid Solvers , 2004, SIAM J. Sci. Comput..

[44]  A. Shiplyuk,et al.  Leading-edge receptivity of a hypersonic boundary layer on a flat plate , 2001, Journal of Fluid Mechanics.

[45]  M. Malik.  Numerical methods for hypersonic boundary layer stability , 1990 .

[46]  E. Koschmieder.  Taylor vortices between eccentric cylinders , 1976 .

[47]  Zhilin Li,et al.  Immersed interface methods for moving interface problems , 1997, Numerical Algorithms.

[48]  Xiaolin Zhong,et al.  Leading-edge receptivity to free-stream disturbance waves for hypersonic flow over a parabola , 2001, Journal of Fluid Mechanics.

[49]  Randall J. LeVeque,et al.  Immersed Interface Methods for Stokes Flow with Elastic Boundaries or Surface Tension , 1997, SIAM J. Sci. Comput..

[50]  A. H. Boudreau.  A Correlation of Artificially Induced Boundary-Layer Transition Data from Blunt Slender Cones at Hypersonic Speeds , 1979 .

[51]  Kazufumi Ito,et al.  Solving a Nonlinear Problem in Magneto-Rheological Fluids Using the Immersed Interface Method , 2003, J. Sci. Comput..

[52]  C. David Pruett,et al.  Spatial direct numerical simulation of high-speed boundary-layer flows Part II: Transition on a cone in Mach 8 flow , 1995 .

[53]  James P. Keener,et al.  An Immersed Interface Method for Solving Anisotropic Elliptic Boundary Value Problems in Three Dimensions , 2003, SIAM J. Sci. Comput..

[54]  H. S. Udaykumar,et al.  A Sharp Interface Cartesian Grid Methodfor Simulating Flows with ComplexMoving Boundaries , 2001 .

[55]  Zhilin Li.  A Fast Iterative Algorithm for Elliptic Interface Problems , 1998 .

[56]  P. Colella,et al.  A Cartesian Grid Embedded Boundary Method for Poisson's Equation on Irregular Domains , 1998 .

[57]  P. Moin,et al.  Direct numerical simulation of transition and turbulence in a spatially evolving boundary layer , 1991 .

[58]  T. Poinsot.  Boundary conditions for direct simulations of compressible viscous flows , 1992 .