Simulation of a Compression Cone with a Freestream Impulsive Entropy Spot in Mach 6 Freestream


Entropy spot diagram 
Compression cone 
Objectives
 Development of computer code for simulating the interaction of a freestream hotspot with the bow shock and the boundary layer in high order of accuracy
 Simulation of the hotspot perturbed shock layer on compression cone under the effect of a freestream hotspot
 Analysis and investigation on the receptivity and instability growth of hotspot perturbed boundary layer based on the simulation results
 Comparison of numerical results with those obtained by Professor
Schneider’s group in the Purdue Mach6 quiet tunnel
Mean Flow for Compression Cone Cases
Freestream Flow Conditions
 M = 6.0
 T = 52.8 K
 P = 610.8 Pa
 Twall = 300 K
 Re_u = 1.026e7 1/m


Geometrical parameters
 R_flared = 3.0 m
 cone halfangle = 2.0 deg
 cone length = 0.45 m
 grid resolution = 3720 x 240
 R=0.001 m


Objectives
 Our goal is to utilize numerical simulation to study the role of roughness on the stability of a hypersonic boundary layer, mainly focus on roughness effect on transition
 In the current project, we consider 2D roughness on a hypersonic flat plate at Mach 5.92
 The roughness is modeled as a hump and is treated by a high order cut cell method
 Pure mode S (slow acoustic wave), mode F (fast acoustic wave) perturbation at 100KHz, and wall normal velocity pulse which has a frequency range 1MHz are imposed into the mean flow separately
 The evolution of disturbance on the wall due to different roughness height and locations are studied using FFT


Discrete roughness element 
Four separate computational cases with one roughness element 
FFT results for 100 khz slow acoustic disturbance (Mode S) with different roughness heights


Roughness located upstream of synchronization point 
Roughness located downstream of synchronization point 
 Upstream of synchronization point, roughness amplifies disturbance
 The overall amplification of disturbance depends roughness height
 However, when roughness is placed downstream of synchronization
point, it damps disturbance with strength related to roughness height
2D roughness movies
Mode S disturbance
Free stream waves receptivity to nonlinear breakdown over Mach 5.5 circular cone
Objectives
 develop a direct numerical simulation code to simulate hypersonic boundary layer flow from laminar to nonlinear breakdown in transition over blunt circular cone
 build a linear receptivity database for different types of freestream waves with a selected frequency specturm
 construct the inflow condition for subsequent nonlinear region using receptivity database and conduct the breakdown simulation
Linear receptivity simulations were conducted to obtain freestream waves database. The database were used as entrance conditions for subsequent nonlinear breakdown simulations.
3D nonlinear breakdown simulation: pressure disturbances contour at wall surface
Numerical Study of WaveEnergy Extraction Device with Electroactive Polymer
Objectives
 Collaborate with the Soft Material Research Laboratory of Prof. Pei (UCLA) to develop technology using electroactive polymer to extract and convert wave energy to electricity
 Implement a finite volume method with freesurface tracking done by Volume of Fluid method to simulate the dynamic interaction between energycontaining incident wave and the energy extracting device
 In preliminary numerical study, the targeted strain of 20 to 30 percent was achieved, suggesting dielectric elastomer may be practical for largescale power generation
Change from high to low capacitance state raises the potential of the charge of the film
Schematic of numerical wave tank
Effects of Graphite Ablation Induced Outgassing on Hypersonic Boundary Layer Stability
Objectives
 Study effects of graphite ablation on the receptivity process and linear growth of Mack's first and second modes in a high speed boundary layer
 Surface chemistry model to approximate ablation includes oxidation, recombination of atomic oxygen and sublimation of C_3, C_2 and C
 Develop nonequilibrium LST code for ablation gas model and compare stability results to DNS
Direct Numerical Simulation Code Validation
11 species ablation code validation to Keenan's computations. Left: surface mass fractions. Right: surface mass flux.
Surface chemistry model comparison
Direct Numerical Simulation Results


Instantaneous snapshot of oxygen (O2) and translationrotation temperature perturbation with 550 khz freestream fast acoustic forcing waves 
Wall pressure perturbation amplitude for real gas simulation with ablation (left) and ideal gas simulation (right). Real gas flow is unstable while ideal gas is stable. Significant difference is shown between the two gas models.
Motivations
 Porous coating significantly stabilizes Mack’s second mode (acoustics) whereas it moderately destabilizes the first mode (TS wave)
 In previous studies, porous coating covers either the entire flat plate or the surface around half the cone circumference
 No work is reported on the effect of porous coating location on the firstmode destabilization, and how to increase stabilization efficiency of porous coating


Felt metal coating  Regular porous coating 
 The two types of porous coatings both destabilize mode S in Mack's first mode region and stabilize it in Mack's second mode region.
 Pressure perturbation amplitude decreases as the phase angle of admittance decreases.
 The destabilization/stabilization behavior of porous coating is affected by the phase angle of admittance.
 For pure mode S propagating downstream, the overall effect of porous coating is destabilizing.
 At approximately the same porosity, regular porous coating is weaker in firstmode destabilization and secondmode stabilization than feltmetal porous coating.
 For regular porous coating, porosity decreasing leads to even weaker firstmode destabilization and secondmode stabilization.
Conclusions  There are three approaches to increase the stabilization efficiency of porous coating
 Put porous coating downstream of the synchronization point
 Use regular porous coating if one must put porous coating along the whole surface
 Design new porous coating of smaller admittance phase angle to attenuate first mode destabilization
Objectives
 To conduct extensive DNS studies on strong shock and turbulence interactions for perfect gas flow with mean Mach numbers ranging from 2 to 30
 To validate our new 3D highorder shockfitting code for nonequilibrium flow
 To conduct DNS studies on strong shock and turbulence interactions for nonequilibrium flow
Numerical Analysis


Inflow: Temporal DNS  DNS with shock fitting 
Variation of various turbulence statistics in simulation of decaying isotropic turbulence.
Conclusions
 Increasing shockstrength reduces the shock deformation
 The trend of maximum values of streamwise vorticity fluctuations versus shock Mach number reverses at M1 = 2.8 from increase to decrease
 The trend of maximum values of Reynolds stress R11 versus shock Mach number reverses at M1 = 8.8 from decrease to increase
Objectives
 Develop and validate new highorder shockfitting numerical method for hypersonic nonequilibrium flow simulations
 Conduct numerical simulation studies on transient growth of reentry vehicles, including the effects of small/finite surface roughness, the interaction of surface roughness and freestram disturbances
 Update and validation of the code with more advanced 5species models and the more realistic 11species air models are ongoing
 The simulation of the steady flow over 9degree halfangle cone shows that real gas effects have significant impacts on flow fields
Two Temperature Model and Governing Equations
 Translation and rotation energy modes are in equilibrium at translation temperature (T)
 Vibration and electronic energy modes are in equilibrium at vibration temperature (Tv)


Two temperature comparison along the stagnation line  Species densities along the stagnation line 
Motivations
 Receptivity studies have focused more on freestream disturbances than on wall disturbances.
 Wall disturbances, together with freestream disturbances, are the main disturbances that hypersonic vehicles experience under the real flight conditions.
 Blowing–suction is not only the most sensitive wall disturbance for hypersonic boundary layers, but is also widely used to control the boundarylayer transition.
Mode S is strongly excited only when the blowingsuction actuator is upstream of the synchronization point. When the blowingsuction actuator is downstream of the synchronization point (Case 7), there is very little excitation of mode S, despite the fact the actuator is still located within the unstable region of mode S.
Conclusions
 The synchronization point of mode F and mode S played an important role in the excitation of mode S.
 To excite strong mode S at a specific frequency, it is necessary to place the blowing–suction actuator upstream of the corresponding synchronization point.
 New control strategy of hypersoinc boundarylayer transition can be developed based on such results.