Current Research Topics

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

Yuet Huang

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 Mach-6 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 half-angle = 2.0 deg
  • cone length = 0.45 m
  • grid resolution = 3720 x 240
  • R=0.001 m

Compression cone entropy spot movies

Zone 12 pressure disturbance
Zone 12 entropy disturbance

Numerical Study of Hypersonic Boundary Layer Instability Affected by 2D Surface Roughness

Danny Fong, Le Duan, and Xiaowen Wang

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

Jia Lei

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 Wave-Energy Extraction Device with Electroactive Polymer

Kin Wai Leung

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 free-surface tracking done by Volume of Fluid method to simulate the dynamic interaction between energy-containing 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 large-scale 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

Clifton Mortensen

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 translation-rotation 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.

Efficient Control of Hypersonic Boundary Layer Transition Using Surface Porous Coating

Xiaowen Wang

Motivations

  • Porous coating significantly stabilizes Mack’s second mode (acoustics) whereas it moderately destabilizes the first mode (T-S 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 first-mode destabilization, and how to increase stabilization efficiency of porous coating
Felt metal coatingRegular porous coating

Numerical Analysis on A Mach 5.92 Flat-Plate Boundary Layer

The hypersonic boundary layer is destabilized or stabilized by feltmetal coating when it is upstream or downstream of the synchronization point.

  • 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 first-mode destabilization and second-mode stabilization than felt-metal porous coating.
  • For regular porous coating, porosity decreasing leads to even weaker first-mode destabilization and second-mode 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

Direct Numerical Simulation of Strong Shock and Turbulence Interaction Using Shock-Fitting Methods

Xiaowen Wang and Pradeep Rawat

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 3-D high-order shock-fitting code for non-equilibrium flow
  • To conduct DNS studies on strong shock and turbulence interactions for non-equilibrium flow

Numerical Analysis

Inflow: Temporal DNSDNS with shock fitting


Variation of various turbulence statistics in simulation of decaying isotropic turbulence.

Conclusions

  • Increasing shock-strength 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

Development and validation of new high order numerical method for hypersonic nonequilibrium flow simulations

Xiaowen Wang and Akshay Prakash

Objectives

  • Develop and validate new high-order shock-fitting numerical method for hypersonic nonequilibrium flow simulations
  • Conduct numerical simulation studies on transient growth of re-entry 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 5-species models and the more realistic 11-species air models are ongoing
  • The simulation of the steady flow over 9-degree half-angle 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 lineSpecies densities along the stagnation line

Receptivity and Stability of a Hypersonic Boundary Layer to Wall Blowing–Suction

Xiaowen Wang and Yanbao Ma

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 boundary-layer transition.

Mode S is strongly excited only when the blowing-suction actuator is upstream of the synchronization point. When the blowing-suction 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 boundary-layer transition can be developed based on such results.