Shock Wave Boundary Layer Interaction and Passive Flow Control Devices
Shock Wave Boundary Layer Interaction and Passive Flow Control Devices
Submitted By
iii
DAYANANDA SAGAR COLLEGE OF ENGINEERING
DEPARTMENT OF AERONAUTICAL ENGINEERING
Kumaraswamy Layout, Bengaluru-560078
2)
iv
DECLARATION
CHETHAN S
(1DS14AE017)
D SHOBITHA
(1DS14AE063)
v
ACKNOWLEDGEMENT
Before introducing our thesis work, we would like to thank the people without whom
the success of this thesis would have been only a dream.
It is with great pleasure, we extend our gratitude and thanks to Dr. C P S Prakash,
Principal, Dayananda Sagar College of Engineering, for his encouragement throughout the
project.
We express our deep sense of gratitude and indebtedness to Dr. A Arokkiaswamy,
Professor and HOD, Department of Aeronautical Engineering, for his valuable guidance,
continuous assistance and in the critical appraisal of the thesis and also for providing the
facilities required for the completion of this project work.
We are very thankful to Dr. S B Verma, Deputy Head, Experimental Aerodynamics
Division, CSIR-NAL for his technical support for the present project.
We are very grateful to Mr. C Manisankar, Scientist, Experimental Aerodynamic
Division, CSIR-NAL for being a great source of knowledge and help during design and
analysis stage and guidance in spite of his busy schedule.
Our sincere thanks to Dr. T K Pradeepa, project coordinator for his encouragement
and guidance.
We express our heartfelt thanks to all our family members and friends and all those
who have directly or indirectly helped us during our course.
vi
ABSTRACT
vii
TABLE OF CONTENTS
INNER TITLE PAGE………………………………………………….………………...i
CERTIFICATE………………………………..…………………………………………ii
DECLARATION................................................................................................................ v
ACKNOWLEDGEMENT ................................................................................................ vi
ABSTRACT ......................................................................................................................vii
CHAPTER-1 ....................................................................................................................... 1
INTRODUCTION ............................................................................................................... 1
CHAPTER-2 ....................................................................................................................... 3
viii
CHAPTER-3 ....................................................................................................................... 9
CHAPTER-4 ..................................................................................................................... 11
METHODOLOGY ............................................................................................................ 11
CHAPTER-5 ..................................................................................................................... 18
5.2.3. POST-PROCESSOR.................................................................................... 20
ix
5.7.1. MESHING ................................................................................................... 26
CHAPTER-6 ..................................................................................................................... 31
CHAPTER-7 ..................................................................................................................... 46
REFERENCES .................................................................................................................. 47
x
LIST OF TABLES
RVGs
xi
LIST OF FIGURES
xii
TABLE OF ABBREVIATIONS
xiii
TABLE OF NOMENCLATURE
xiv
ANALYSIS OF SWBLI AND ASSESSMENT OF PASSIVE FLOW CONTROL DEVICE
CHAPTER-1
INTRODUCTION
1.1. OBJECTIVES
The main objectives of the present project work are following:
• To carry out a validation case of shock wave boundary layer interaction over a flat
plate.
• To conduct viscous 2-D and 3-D numerical flow simulation (CFD analysis) of a
SWBLI over a flat plate.
• To perform grid independent study for a typical case of free stream Mach number
of 2.05.
• To investigate the flow characteristics of a rectangular vane type micro vortex
generators.
• To study the capabilities of the rectangular vane type micro vortex generators for
controlling the flow separation induced by SWBLI.
• To study the effect of variation of freestream Mach number on the recirculation
zone thickness.
• The 2D and 3D geometries were generated using ANSYS 15.0. A structured mesh
has been generated using Hexa-Dominant method to capture the flow field.
• A numerical analysis has been carried out in CFD – Fluent. At inlet of the fluid
domain, initial static pressure, static temperature and velocity are given as
boundary conditions.
• CFD is used to solve the defined problem, the K-ω SST model is opted for
turbulence model.
CHAPTER-2
2.1. INTRODUCTION
The prominent problem faced by high speed flow intakes or on lifting surfaces is the
complex phenomenon of shock wave boundary layer interactions (SWBLI) which leads to
boundary layer separation and adverse pressure gradients. The goal is to minimize the total
pressure losses. Therefore, a fundamental understanding of this physical phenomenon is
important for the development of improved techniques to control their counterproductive
effects.
Fig.2.1 shows the schematic of the incident shock induced separation when
undergoing regular reflection. The impingement of the incident shock wave(C1) causes a
large pressure difference across it hence induces the boundary layer to separate due to the
consequent effect of severe adverse pressure gradient. The separation occurs well ahead the
impingement location, Ximp. The upstream propagation of the pressure gradient within the
subsonic part of the turbulent boundary-layer (TBL) induces compression waves in the
supersonic part of the TBL, which coalesce to the reflected shock (C2). This reflected shock
(C2) intersects the incident shock (C1) at ‘I’ and it forms transmitted shocks (C3 and C4).
The shock (C4) penetrates into the separated shear layer and is reflected as an
expansion fan which causes the flow to reattach. This flow turning – reattachment process
results in the formation of compression waves that merge together to form a reattachment
shock (C5). Finally, the boundary layer reaches equilibrium downstream of the
SWBLI[1][2].
If the strength of the induced shock (C1) reaches the detachment limit for a given
freestream Mach number, the intersection point ‘I’ gets replaced with a normal shock which
is known as a Mach Stem. In both the causes the separation formed acts as an obstacle to
the incoming flow.
Fig 2. 2 Static Pressure along Plate Length showing First and Second Pressure
rise[3].
Fig.2.2. shows the separation point is denoted by S and further downstream exists a
recirculating ‘bubble’. This area is bounded by a layer called dividing streamline (S) that
separates the recirculating region from the incoming flow. The dividing streamline
originates from point S and ends at the reattachment point R. Strong mixing takes place in
the detached shear layer emanating from S and causes mechanical energy transfer to occur
from outer high-speed flows into the separated region. As a result, the velocity on the
dividing streamline (S) increases gradually until it approaches the reattachment point (R)
where deceleration starts[3].
Department of Aeronautical Engineering, DSCE, Bengaluru-78 Page 4
ANALYSIS OF SWBLI AND ASSESSMENT OF PASSIVE FLOW CONTROL DEVICE
affect not only the propulsion system but also endanger the structural strength of any aero-
vehicle[3].
2.3.2. SUCTION
Suction is applied locally inside the interaction region or in the immediate vicinity. It
is done through drilling holes perpendicularly to the surface. The other way of implementing
this method is at a certain distance upstream (usually specified in scale of the boundary layer
thickness, δ) of the shock impingement. This will produce a fuller velocity profile and more
robust boundary layer towards separation.
body is in motion relative to the air, the VG creates a vortex, which, by removing some part
of the slow-moving boundary layer in contact with the airfoil surface, delays local flow
separation and aerodynamic stalling, thereby improving the effectiveness of wings and
control surfaces, such as flaps, elevators, ailerons, and rudders.
Vortex generators are most often used to delay flow separation. Conventional VGs
which have the height of the boundary layer thickness, δ have been used for decades to
control the separation by increasing the near-wall momentum by transferring the high
momentum flow from the outer region into the lower momentum region, at the wall. This
type of flow control, which is usually characterized as passive flow control, so doesn’t
require the addition of extra energy to the system nor it needs any actuators. These are
becoming more favorable despite active devices have their own attractions. More researches
proved that the micro VGs, termed after having the height less than the boundary layer
thickness, are able to outperform the performance of the conventional.
CHAPTER-3
LITERATURE REVIEW
During the literature review on design of various Micro-Vortex Generators and its
effects on Shockwave Boundary Layer Interactions in various high-speed flows, it is
observed that various studies and analysis
S. B. Verma and C. Manisankar[6] investigated the performance of several Micro-
Vortex Generators of height 30% and a ramp-vane of height 50% of the local boundary layer
thickness(δ) of no-control case which was implemented at 10δ upstream of separation
location.
Shashi.B. Verma and Manisankar[2] conducted experiments by splitting the
trapezoidal control configuration into rectangular vane type control devices and found that
a local pressure relieving effect is created and also found out that the height of the Vortex
Generators is significant to control the extent of separation.
S. B. Verma and A. Hadjadj[7] focused on the flow control mechanisms by
implementing various VG configurations which effect the flow unsteadiness, bow-shock
strengths and reduce the overall sound pressure levels by 4.5dB.
Amjad Ali Pasha[8] studied the effect of deflection angle, the Reynolds Number, the
wall temperature and the Mach Number on the recirculation zone in hypersonic SWTBLI
flows using CFD as a tool for compression ramp in which it was observed that by keeping
the above-mentioned parameters constant, a decrease in freestream Mach number results in
larger length of the recirculation zone.
Mohd R. Saad, Hossein Zare-Behtash, Azam Che-Idris and Konstantinos
Kontis[9] stated that the height of the Micro-Vortex Generators should be in the range 30%
- 90% of the local boundary layer thickness, which are embedded inside the boundary layer,
which is intended to reduce parasitic drag relative to the conventional full-size VG.
Francis K. Acquaye[10] evaluated various turbulence models for three different
supersonic flow conditions and compared WA, SA, and SST k-ω models and found out that
SST k-ꞷ model made accurate predictions of the flow separation in SWBLI.
Mohd R. Saad[3] investigated the micro ramps in detail using various flow diagnostic
methods. The oil-flow and oil-dot visualization enabled technique the visualization of the
detail structures. A significant alteration in the separation region and minimization in the
reversed flow region was also observed.
Frank K. Lu, Qin Li, Yusi Shih, Adam J. Pierce and Chaoqun Liu[11] performed
an assessment of various micro vortex generator configurations and studied that these MVG
adds momentum to the boundary layer which acts as a beneficiary factor in reducing the
shock induced separation.
Thomas Herges, Erik Kroeker, Greg Elliott and Craig Dutton[12] investigated
micro ramp array at Mach 1.4 and its effects of on normal SWBLI. The height of the vortex
generators was scaled to 40% of the incoming boundary layer thickness. The incompressible
displacement thickness was decreased for the SBLI with the micro-ramp array relative to
the no-array case.
The survey of literature indicates that specific and focused research work has been
carried out all over the world on the flow physics in the complex phenomena SWBLI. It
includes study of various parameters affecting the boundary layer thickness, separation
location and recirculation zone height. Most efforts to determine the optimum height of the
micro vortex generators.
CHAPTER-4
METHODOLOGY
𝛾
𝑃𝜊
= [1 + (𝛾 − 1)/2 ∗ 𝑀2 ](𝛾−1) ()
𝑃1
𝑇𝜊
= [1 + (𝛾 − 1)/2 ∗ 𝑀2 ] ()
𝑇1
M = Mach Number
𝑃𝜊 = Total Pressure
𝑇𝜊 = Total Temperature
𝑇1 = Static Temperature
𝛾 = Specific Heat Ratio
𝜌2 (𝛾+1)∗𝑀1 2 ∗(𝑠𝑖𝑛2 𝛽)
= 2+(𝛾−1)∗𝑀 2 ∗(𝑠𝑖𝑛2 𝛽 ) ()
𝜌1 1
4.5. 3D MODELING
3-D modeling is done using ANSYS 15.0. A flat plate of length 340mm is made and
extruded 110 mm. A fluid domain of volume 840x360x500 mm is created. Wedge of
154x160x12.5 with a turning angle of 14˚ is made. Fig. 4.2. shows the 3-D geometry divided
into several blocks.
Fig.4.4. shows the arrangement of RVG on the surface of the Flat Plate at a distance
of 5.5δ upstream of the Separation line. The RVGs produce counterrotating vortices that add
momentum to TBL.
CHAPTER-5
Computational fluid dynamics (CFD) is the science of predicting fluid flow, heat
transfer, mass transfer, chemical reactions, and related phenomena by solving the
mathematical equations which govern these processes using a numerical process. Computers
are used to perform the calculations required to simulate the interaction of liquids and gases
with surfaces defined by boundary conditions. With high-speed super computers, better
solutions can be achieved. Ongoing research yields software that improves the accuracy and
speed of complex simulation scenarios such as transonic or turbulent flows. Initial
experimental validation of such software is performed using a wind tunnel with the final
validation coming in full-scale testing, e.g. flight tests. The result of CFD analyses is relevant
engineering data used in, Conceptual studies of new designs, Detailed product development,
Troubleshooting, Redesign.
Fluid dynamics is the science of fluid motion. The study of the fluid flow can be
possible in three various ways as 1) Experimental 2) Theoretical and 3) Numerically. The
numerical approach is called Computational fluid dynamics. CFD uses numerical methods
and algorithms to solve and algorithms to solve and analyze problems that involve fluid
flows by using computers.
5.2.1. PRE-PROCESSOR
Pre-processor consists of the input of the flow problem to a CFD program by means
of an operator friendly interface and the subsequent transformation of this input into a form
suitable for use by the solver. The region of fluid to be analyzed is called the computational
domain and it is made up of a number of discrete elements called the mesh (or grid). After
the mesh generation, to define the properties of fluid and to specify appropriate boundary
condition.
5.2.2. SOLVER
Solver calculates the solution of the CFD problem by solving the governing equations.
The equations governing the fluid motion are Partial Differential Equations (PDE), made up
of combinations of the flow variables (e.g. velocity and pressure) and the derivatives of these
variables. Computers cannot directly produce a solution of it. Hence the PDEs must be
transformed into algebraic equations. This process is known as numerical discretization.
There are four methods for it as 1) Finite difference method 2) Finite element method and
3) Finite volume method and 4) Spectral method. The finite difference method and the finite
volume method both produce solutions to the numerical equations at a given point based on
the values of neighboring points, whereas the finite element method produces equations for
each element independently of all other elements.
5.2.3. POST-PROCESSOR
It is used to visualize and quantitatively process the results from the solver part. In a
CFD package, the analyzed flow phenomena can be presented in vector plots or contour
plots to display the trends of velocity, pressure, kinetic energy and other properties of the
flow.
about the solution than for either FDM or FVM. In FEM, residuals are used to measure the
errors.
CONTINUITY EQUATION
Continuity equation it is based on the principle of conservation of mass. Net mass flow
out of control volume = time rate of decrease of mass inside control volume mass
conservation equation.
𝐷𝜌 𝜕𝑈𝑖
+𝜌 =0 ()
𝐷𝑡 𝜕𝑥𝑖
MOMENTUM EQUATION
ENERGY EQUATION
Energy equation it is based on the principle that total energy is conserved. Total energy
entering control volume = total energy leaving control volume energy equation.
𝜕𝑇 𝜕𝑇 𝜕𝑈𝑖 𝜕2 𝑇 𝜕𝑈𝑗
𝜌𝑐
⏟ 𝜇 𝜕𝑡 + 𝜌𝑐 𝑈 = −𝑃 + 𝜆 2 − 𝜏𝑖𝑗 (5.4)
⏟ 𝜇 𝑖 𝜕𝑥𝑖 ⏟ 𝜕𝑥𝑖 ⏟ 𝜕𝑥 ⏟ 𝜕𝑥𝑖
i ii iii iv v
where 𝑢̅𝑖 and 𝑢𝑖′ are the mean and fluctuating velocity components (i = 1,2,3).
Substituting expressions of this form for the flow variables into the instantaneous
continuity and momentum equations and taking a time (or ensemble) average (and dropping
the overbar on the mean velocity, 𝑢̅ ) yields the ensemble-averaged momentum equations.
They can be written in Cartesian tensor form as[13]:
(5.7)
𝜕𝑘 𝜕𝑘 ̅𝑖
𝜕𝑢 𝜕 𝜈𝜏 𝜕𝑘
+ 𝑢̅ = 𝜏𝑖𝑗 −𝜖+ [(𝜈 + ) ] (5.8)
𝜕𝑡 𝜕𝑥𝑗 𝜕𝑥𝑗 𝜕𝑥𝑗 𝜎𝜅 𝜕𝑥𝑗
𝜕𝜖 𝜕𝜖 𝜖 ̅𝑖
𝜕𝑢 𝜖2 𝜕 𝜈𝜏 𝜕𝜖
+ 𝑢̅𝑗 = 𝐶𝜖1 𝜏𝑖𝑗 − 𝐶𝜖2 ∗ + [(𝜈 + ) ] (5.9)
𝜕𝑡 𝜕𝑥𝑗 𝜅 𝜕𝑥𝑗 𝜅 𝜕𝑥𝑗 𝜎𝜅 𝜕𝑥𝑗
𝜕𝑘 𝜕𝑘 ̅′ 𝑖
𝜕𝑢 𝜕 𝜕𝑘
+ 𝑢̅𝑗 = 𝜏𝑖𝑗 − 𝛽 ∗ 𝑘𝜔 + [(𝜈 + 𝜎 ∗ 𝑘/ 𝜔) ] (5.10)
𝜕𝑡 𝜕𝑥𝑗 𝜕𝑥𝑗 𝜕𝑥𝑗 𝜕𝑥𝑗
𝜕𝜔 𝜕𝜔 𝜔 ̅′ 𝑖
𝜕𝑢 𝜎𝑑 𝜕𝑘 𝜕𝜔 𝜕 𝑘 𝜕𝜔
+ 𝑢̅𝑗 = 𝛼 𝜏𝑖𝑗 − 𝛽 ∗ 𝜔2 + + [(𝜈 + 𝜎 ∗ ) ]
𝜕𝑡 𝜕𝑥𝑗 𝜅 𝜕𝑥𝑗 𝜔 𝜕𝑥𝑗 𝜕𝑥𝑗 𝜕𝑥𝑗 𝜔 𝜕𝑥𝑗
(5.11)
5.7.1. MESHING
The finite volume method (FVM) is a common approach used in CFD codes, as it has
an advantage in memory usage and solution speed, especially for large problems, high
Reynolds number turbulent flows. In the finite volume method, the governing partial
differential equations are recast in a conservative form, and then solved over discrete control
volumes. This discretization guarantees the conservation of fluxes through a particular
control volume.
The required mesh size is defined and various steps are performed to produce a mesh
of high quality. Grid angles are exceptionally good, mesh sizes transition smoothly, and high
aspect-ratio elements are generated in the near-wall regions to resolve these regions
efficiently and capture boundary layer and recirculation zone accurately.
The total geometry shown in the fig. 5.2 is divided into several blocks of different
sizes. The block size is smaller between the wedge and the flat plate where the whole
shockwave boundary layer interaction phenomena occurs and larger blocks are made in the
rest of the domain. The whole domain is divided into a total of 49 blocks out of which 47
are cuboids and 2 are with triangular face.
All the cuboids are given edge sizing with respect to block dimensions i.e., 2/3th of
the block edge length. The remaining 2 triangular faced blocks are given body sizing of
1mm and 3mm which are near to the wedge and the plate respectively.
Fig. 5.2. Meshing pattern is seen with 25,58,000 grid elements and 28,73,000 nodes.
5.7.2. PREPROCESSING
After finishing the meshing process, ANSYS-Fluent boundary conditions are defined.
5.7.3. ANALYSIS
The file is solved in CFD – Fluent for the input boundary conditions as given. This
solver uses RANS equations to solve for required parameters. Residual criteria are given as
10−8. Typical convergence residual graph is shown.
involves performing the simulation on two or more successively finer grids. The term grid
convergence study is equivalent to the commonly used term grid refinement study.
The analysis is carried out for grid elements of 1058000, 2558000 and 5201000 of
coarse, medium and fine mesh.
CHAPTER-6
CFD analysis has been carried out for 2-D and 3-D – with and without flow control
devices. Results of each value has been studied and plotted. It will further be discussed in
detail.
Fig. 6.3. shows the zoomed in velocity vector plot to view the recirculation zone where
it is observed that velocity reduces to subsonic range at the near wall flow regime due to the
low momentum boundary layer.
Fig. 6.4. Normalized Pressure Ratio plot for CFD Results and Experimental Data
Fig.6.6. shows that the Maximum Pressure Ratio increases with the Mach number and
the Separation point shifts downstream.
Fig. 6.6. Pressure Ratio vs Plate Length for Various Freestream Mach Numbers
Fig. 6.7. Variation of Recirculation zone height with respect to Mach Number
Fig.6.7. shows that the recirculation zone height increases with the increase in Mach
Number.
Fig. 6.8. Variation in Mach Stem Height with respect to Mach Number
Fig.6.11. shows the Velocity Contour along with the small recirculation zone at the
bottom in blue colour. The Vector Contour of zoomed recirculation zone is shown in
Fig.6.12. and vector direction shows direction of flow.
Fig. 6.13. shows the pressure ratio plot for coarse, medium and fine mesh. And it is
observed that the locations of pressure rise are appropriate and constant in the case of
medium and fine mesh, so the medium mesh case is considered for the analysis to cut down
the computation time.
𝑷𝒔
Fig. 6.13. Variation of Vs Plate Length for different mesh qualities
𝑷∞
Table 6.3 Percentage Error for 3-D Results for various Mesh Qualities without
RVG
ACTUAL VALUES COARSE MEDIUM FINE
𝑀2 = 1.53 1.778 % 1.699 % 1.653 %
Fig. 6.14. Pressure Contour with RVG for freestream Mach Number 2.05
Fig 6.15 shows the Pressure Contour on plate surface and corrugation patterns in
the separation line after employing RVGs.
(a) Height=0.01mm
Fig. 6.16. Velocity Contours for Various Offset Planes from Plate surface
Fig. 6.17. Velocity Contour with RVG for freestream Mach Number 2.05
Fig. 6.18. Recirculation Zone with RVG for Freestream Mach Number 2.05
After implementing the RVGs, it has been observed that the recirculation zone has
been pushed downstream due to the interaction of vortices with the separation line.
Fig. 6.20. Pressure Ratio plots for RVG and No-RVG case
Fig 6.20 shows comparison between Pressure Ratio plots for RVG and No- RVG
case. Reduction in first and second pressure ratio can be seen which a direct measure of
reduction in Shock- Foot unsteadiness.
δ 4.5 3.6 20 %
CHAPTER-7
7.1. CONCLUSION
Numerical investigations were carried out to study the shockwave boundary layer
interactions on a flat plate of 2-D and 3-D cases for a freestream Mach number of 2.05.
Wedge of a flow turning angle 14˚ was used to generate the incident shock wave over the
flat plate. The effect of variation in freestream Mach number on the recirculation zone
thickness is studied. The analysis was carried out for with and without RVG cases. Also,
assessment of the effectiveness of the Micro Rectangular Vane Type Vortex Generators in
controlling shock induced separation in flow of freestream Mach number 2.05 is studied.
The RVG height spanned 37.78% of the local BLT of no-RVG case. Array of RVGs is
placed at 5.5δ upstream of the separation location of no-RVG case. Following are the
conclusions drawn from the present study:
• The parametric study revealed that the Maximum Pressure Ratio increases with
the Mach number and the Separation point shifts downstream.
• The Recirculation Zone height increases with Mach Number because the shock
deflection angle reduces leading to reduction in shock strength.
• The Mach steam length decreases with increasing Mach Number and eventually
disappears.
• The separation point is shifted 2.796% (1.34δ) downstream with respect to the
separation point for no control case.
REFERENCES