High-Fidelity Computational Multi-Physics Laboratory: Projects
Hypersonic Flight Analysis
Hypersonic flight is characterized by a harsh environment which results in extreme temperatures and heat transfer rates. During flight, the flow regime changes from laminar to turbulent, with corresponding changes in drag. Ground tests cannot often reproduce what happens in flight, while flight tests are expensive and difficult to instrument. We are advancing simulation techniques to bridge this gap by exploring the main phenomena whose signatures are clearly evident in experiments.
Accordions
Leveraging experimental data, turbulence models are assessed for their capacity to capture pressure and heat-flux loads in complex transitional and Shock/Boundary Layer Interactions (SBLIs) environments.
Jet Noise
Our goals in this research are to answer fundamental questions on jet noise and its control. For example:
- What is the relation between coherent structures in the flow and noise?
- How does the jet filter turbulent distrubances to generated fluctuations in the near acoustic field?
- What is the nature of hydrodynamic and acoustic fluctuations in the near acoustic field?
- How exactly do active control techiniques based on small perturbations influence noise generation?
We employ experimental data extensively to ensure that our simulations represent reality. The validated simulations are then post-processed with sophisticated statistical techniques to generate enabling insight. Some of our recent work may be viewed with the links below.
Accordions
Large Eddy Simulation of supersonic jets impinging on a flat surface replicate the propulsion system of VTOL aircraft in hover. Our research focus is on understanding the turbulent fountain-flow generated underneath an aircraft and the intense sound produced from aeroacoustic resonance. Observations of the hydrodynamic and acoustic coupling of the two jets informs new models that provide insight and prediction of sound intensity over a range of impinging heights and propulsion configurations.
Simulation-Based Analysis of the Near Field in a Supersonic Jet Controlled by Plasma Actuators
The near field of a Mach 1.3 jet is analyzed under the influence of m = 0 and m = 1 excitation at St = 0.3, 1.0, and 2.0. The analysis, which includes azimuthal mode decomposition of the pressure field, shows that at higher frequencies, the large structure formation is inhibited and peak SPL amplitudes diminish, in a manner that is generally consistent with experimental results.
Scramjet Flowpaths
The hypersonic flight environment provides numerous challenges in achieving sustained combustion with air-breathing propulsion systems. Two particular transient phenomena of interest are the processes of mode-transition and unstart. The objectives of this study are:
- Understand the role of shock/boundary-layer interactions during these transient processes
- Explore the large-scale structures present in these flow fields and their connection to SBLI
- Evaluate the evolution of separated regions and their contributions to such events
Carefully-tailored, time-accurate RANS simulations are computed with validation against the available experimental data.
Accordions
The transition between ramjet and scramjet operation modes in a dual-mode scramjet is not fully understood. The recent flight test of the HIFiRE 2 provided the opportunity to experimentally study this event. In order to supplement this study, computations were performed to understand the influence of inlet distortion, SBLI, and corner separation during the mode-transition event.
The recent flight tests of the X-51 have highlighted the need to better understand the sensitivity of scramjets to unstart events. These typically violent events increase heat and pressure loads on scramjet vehicles and limit control authority of the vehicle. This study is intended to identify precursors of unstart and to understand the sensitivity of SBLI to modulation of heat release as a result of transient fuel-staging.
Shock Interactions
The interaction of shock waves with turbulent boundary layers occurs in all supersonic and hypersonic aircraft. The consequences are usually severe and include loss of control authority in external flows and distortion losses in propulsion systems. Our simulations explore the many individual processes, including three-dimensional separation, unsteadiness and vortical structure generation.
Wing Stall
The phenomena that occur at high angles of attack are extremely complex. In particular, separation and laminar-turbulent transition, accompanied by vortical structure formation are major features. We seek to answer the following questions:
- What are the mechanisms by which three-dimensional phenomena set in for a stalled wing section?
- What is the effect of Reynolds number on the stall phenomenon?
- How do active actuators, based on nano-second pulsed discharges couple to the flow?
- What is the optimal frequency of excitation that can produce the most response from the flow?
- What happens in reverse flight encountered in forward flight of helicopters?
We employ experimental data extensively to ensure that our simulations represent reality. The validated simulations are then post-processed with sophisticated statistical techniques to generate enabling insight. Some of our recent work can be viewed through use of the links below.
Accordions
A Semi-Empirical Model of a Nanosecond Pulsed Plasma Actuator For Flow Control Simulations with LES
A phenomenological model, suitable for coupling to Large Eddy Simlations (LES) of flow control applications, is developed to simulate the effects of a nanosecond dielectric barrier discharge (NS-DBD) actuator. The model considers various surface and volume heating profiles to reproduce the key qualitative and quantitative features of NS-DBD actuators. The model is tuned to match the unique wave shape observed in Schlieren imagery and quantitative wave speed and displacement data. The induced wave structure, comprised of a cylindrical component with a tail, can be reproduced with a spatially varying heating distribution - a Gaussian variation in the direction of the discharge provides accurate results. The instantaneous location of the wave, which propagates at approximately acoustic speed after an initial transient may be matched by considering a volumetric deposition model. Power input estimates in the current model are consistent with those postulated in experimental investigations.
Characterization of Sharp-Edged Airfoils Using LES
Large-Eddy Simulations are performed to characterize the flow over a sharp-edged (retreating) NACA0015 wing section. The goals were to 1) demonstrate the need for three-dimensional simulations, as oppposed to less-expensive two-dimensional simulations; 2) assess the necessary spanwise width to capture the three-dimensional structures and compare to the necessary width for blunt (advancing) airfoils and 3) compare differing angles of attack for three-dimensional simulations to explore flow separation phenomena. Comparisons were made between the simulations as well as experimental results. Two-dimensional simulations are inaccurate because the absence of a third dimension does not allow for spanwise instabilities to develop as are necessary for true breakdown to turbulence. In addition, it was discovered that a spanwise width of 40% of the airfoil chord length is necessary to allow for proper development of spanwise structures, which is twice what is needed for an advancing airfoil. Finally, asymmetric vortex shedding at α=15° was observed and the details of this structure are investigated.
Fluid-Structure Interactions
Fluid-structure interaction (FSI) plays a key role in the design of many engineering systems, e.g. aircrafts, turbines, heat exchangers, bridges etc. Escalating FSI can lead to catastrophic failures such as the infamous first Tacoma Narrows Bridge collapse. FSI are often too complex to solve analytically, and are mostly studies by performing experimental and numerical simulations. The advancement in computational fluid dynamics and computational structural dynamics is used together to simulate FSI problems numerically.
The focus of current FSI research at HFCMPL is high-speed laminar/turbulent boundary layers with/without impinging shock and their interactions with compliant surfaces.
Accordions
Control of transitional shock wave boundary layer interaction using surface morphing
The potential of surface morphing techniques, including passive shock control bumps (SCB) and active surface morphing is explored to control transitional shock wave boundary layer interactions (SWBLI). In addition to reduction of the size of the separation bubble, a key objective is to mitigate the low-frequency unsteadiness that can cause detrimental structural response. To this end, three-dimensional flow simulations are performed using direct numerical simulations (DNS) at Mach 2 and Reynolds number based on inflow boundary layer thickness Re_δin = 996. An incident oblique shock at the shock angle of σ = 35 deg and shock strength of p2/p1= 1.4 impinges on a laminar boundary layer that evolves from a Blasius profile. The boundary layer separation due to the shock impingement leads to G ̈ortler-like instability, where the nominally two-dimensional and steady inflow undergoes flow transition, giving rise to a three-dimensional and unsteady interaction. An aero-structural solver framework is developed and employed to examine surface morphing control. To avoid unrealistic structural deformation in the transient or final states, the structural integrity is concurrently monitored so that the intermediate morphing solutions are restricted to achievable elastic deformation. The results indicate that the transitional SWBLI can be controlled in this manner, to essentially eliminate the three-dimensionality and unsteadiness associated with the G ̈ortler-like vortices. Both passive SCB and active surface morphing reduce separation by modulating the sharp increases in surface pressure at separation and shock-impingement points encountered in uncontrolled SWBLI, without incurring additional loss of the stagnation pressure.
The video shows the effect of active surface morphing on the flow transition exhibited in a shock wave boundary layer interaction, in terms of Q criterion isosurface, completely eliminating the flow separation and transition. The vertical axis is scaled for clarity.
Transitional shock wave boundary layer interaction over a flexible panel
We investigate a full 3-D transitional shock wave boundary layer interaction over a flexible panel by performing high-fidelity direct numerical simulations. A nominally 2-D laminar boundary layer at Mach number 2 interacts with an oblique shock wave with the shock angle of 35 deg and shock strength of 1.8 (p3/p1) in the presence of flexible panel. Simulations are performed for a range of Reynolds numbers (based on panel length), and the coupling between fluid and structure is explored. The flow transition to turbulence occurs at a lower Reynolds number for the flexible panel (Re≈40000) compared to the rigid panel (Re≈70000) simulation. This is manifested through the appearance of unsteadiness and three-dimensionality. The transitional SWBLI exhibits the Görtler instability and low/high frequency unsteadiness, which are characterized interms of Görtler number and wall pressure power spectral density respectively for both the flexible and rigid panels. The principal Reynolds stresses are significantly modified due to the flexible panel, particularly in the near wall region, resulting in overall increased level of turbulence and skin friction coefficient.
The video shows the density gradient magnitude in two planes, displaying unsteady transitional SWBLI over a flexible panel. The vertical axis is scaled by factor 4 for clarity.