Three-dimensional supersonic flow over a spike-nosed body of revolution

The unsteady, viscous, supersonic flow over a spike-nosed body of revolution is numerically investigated by solving the Navier-Stokes equations. The time-accurate computations are performed employing an implicit algorithm based on the second-order time-accurate LU-SGS scheme with the incorporation of a subiteration procedure to maintain time accuracy. The characteristics of the flow field for a Mach number of 3.0, Reynolds number of 7.87 x 10(6)/m, and angles of attack of 5 and 10 degrees are described. Self-sustained asymmetric shock wave oscillations were observed in the numerical computations for these angles of attack. The main characteristic of the flow field, as well as its influence on drag coefficient is discussed.

On the development of an unstructured grid solver for inert and reactive high speed flow simulations

An unstructured grid Euler solver for reactive compressible flow applications is presented. The method is implemented in a cell centered, finite volume context for unstructured triangular grids. Three different schemes for spatial discretization are implemented and analyzed. Time march is implemented in a time-split fashion with independent integrators for the flow and chemistry equations. The capability implemented is tested for inert flows in a hypersonic inlet and for inert and reactive supersonic flows over a 2-D wedge. The results of the different schemes are compared with each other and with independent calculations using a structured grid code. The strengths and the possible weaknesses of the proposed methods are discussed.

Solution of aerospace problems using structured and unstructured strategies

Products developed at industries, institutes and research centers are expected to have high level of quality and performance, having a minimum waste, which require efficient and robust tools to numerically simulate stringent project conditions with great reliability. In this context, Computational Fluid Dynamics (CFD) plays an important role and the present work shows two numerical algorithms that are used in the CFD community to solve the Euler and Navier-Stokes equations applied to typical aerospace and aeronautical problems. Particularly, unstructured discretization of the spatial domain has gained special attention by the international community due to its ease in discretizing complex spatial domains. This work has the main objective of illustrating some advantages and disadvantages of numerical algorithms using structured and unstructured spatial discretization of the flow governing equations. Numerical methods include a finite volume formulation and the Euler and Navier-Stokes equations are applied to solve a transonic nozzle problem, a low supersonic airfoil problem and a hypersonic inlet problem. In a structured context, these problems are solved using MacCormack’s implicit algorithm with Steger and Warming’s flux vector splitting technique, while, in an unstructured context, Jameson and Mavriplis’ explicit algorithm is used. Convergence acceleration is obtained using a spatially variable time stepping procedure.

Panel Method Formulation for Oscillating Airfoils in Sonic Flow

The mathematical model for two-dimensional unsteady sonic flow, based on the classical diffusion equation with imaginary coefficient, is presented and discussed. The main purpose is to develop a rigorous formulation in order to bring into light the correspondence between the sonic, supersonic and subsonic panel method theory. Source and doublet integrals are obtained and Laplace transformation demonstrates that, in fact, the source integral is the solution of the doublet integral equation. It is shown that the doublet-only formulation reduces to a Volterra integral equation of the first kind and a numerical method is proposed in order to solve it. To the authors' knowledge this is the first reported solution to the unsteady sonic thin airfoil problem through the use of doublet singularities. Comparisons with the source-only formulation are shown for the problem of a flat plate in combined harmonic heaving and pitching motion.

Direct Simulation Monte Carlo Method Applied to Aerothermodynamics

An axisymmetric supersonic flow of rarefied gas past a finite cylinder was calculated applying the direct simulation Monte Carlo method. The drag force, the coefficients of pressure, of skin friction, and of heat transfer, the fields of density, of temperature, and of velocity were calculated as function of the Reynolds number for a fixed Mach number. The variation of the Reynolds number is related to the variation of the Knudsen number, which characterizes the gas rarefaction. The present results show that all quantities in the transition regime (Knudsen number is about the unity) are significantly different from those in the hydrodynamic regime, when the Knudsen number is small.

Numerical study of wedge supported oblique shock wave-oblique detonation wave transitions

The results of a numerical study of premixed Hydrogen-air flows ignition by an oblique shock wave (OSW) stabilized by a wedge are presented, in situations when initial and boundary conditions are such that transition between the initial OSW and an oblique detonation wave (ODW) is observed. More precisely, the objectives of the paper are: (i) to identify the different possible structures of the transition region that exist between the initial OSW and the resulting ODW and (ii) to evidence the effect on the ODW of an abrupt decrease of the wedge angle in such a way that the final part of the wedge surface becomes parallel to the initial flow. For such a geometrical configuration and for the initial and boundary conditions considered, the overdriven detonation supported by the initial wedge angle is found to relax towards a Chapman-Jouguet detonation in the region where the wedge surface is parallel to the initial flow. Computations are performed using an adaptive, unstructured grid, finite volume computer code previously developed for the sake of the computations of high speed, compressible flows of reactive gas mixtures. Physico-chemical properties are functions of the local mixture composition, temperature and pressure, and they are computed using the CHEMKIN-II subroutines.