# Implicit LES for Supersonic Microramp Vortex Generator: New Discoveries and New Mechanisms

Authors: Qin Li; Chaoqun Liu
Journal: Modelling and Simulation in Engineering
 (2011)

---

## Abstract

This paper serves as a summary of our recent work on LES for supersonic MVG. An implicitly implemented large eddy simulation (ILES) by using the fifth-orderWENO scheme is applied to study the flow around the microramp vortex generator (MVG) at Mach 2.5 and Reθ=1440.  A number of new discoveries on the flow around supersonic MVG have been made including spiral points, surface separation topology, source of the momentum deficit, inflection surface, Kelvin-Helmholtz instability, vortex ring generation, ring-shock interaction, 3D recompression shock structure, and influence of MVG decline angles. Most of the new discoveries, which were made in 2009, were confirmed by experiment conducted by the UTA experimental team in 2010.  A new 5-pair-vortex-tube model near the MVG is given based on the ILES observation. The vortex ring-shock interaction is found as the new mechanism of the reduction of the separation zone induced by the shock-boundary layer interaction.


---

## Body

## 1. Introduction

It is well known that for the supersonic ramp jet flow, shock boundary layer interaction (SBLI) can significantly degrade the quality of the flow field by triggering large-scale separation, causing total pressure loss, making the flow unsteady and distorted, and can even make an engine unable to start. In order to improve the “health” of the boundary layer, a series of new devices named as micro vortex generator is designed for flow control, which is with a height approximately 20–40% (more or less) of the boundary layer thickness. Because of the robust structure, the specific microramp vortex generator (MVG) becomes more attractive to the inlet designer. Intensive computational and experimental studies have been made on the MVG recently.Lin indicated [1] that the device like MVG could alleviate the flow distortion in compact ducts to some extent and control boundary layer separation due to the adverse pressure gradients. Similar comments were made in the review by Ashill et al. [2]. The formal and systematic studies about the micro VGs including micro ramp VG can be found in the paper by Anderson et al. [3]. Babinsky et al. [4–7] made a series of experiments on different kinds of mirco VGs and investigated their control effects in details. The mechanism of MVG flow control from his work was described as that a pair of counter-rotating primary streamwise vortices is generated by MVG, which is mainly located within the boundary layer and travel downstream for a considerable distance. Secondary vortices are located underneath the primary ones and even more streamwise vortices could be generated under suitable conditions. Streamwise vortices inside the boundary later will bring low-momentum fluids up from the bottom and high-momentum fluids down to the boundary layer. A striking circular momentum deficit region is observed in the wake behind the MVG. The vortices will keep lifting up slowly, which is thought to be the consequence of the upwash effect of the vortices.Numerical simulations have been made on MVG for comparative study and further design purposes. Ghosh et al. [8] made detailed computations under the experimental conditions given by Babinsky by using RANS, hybrid RANS/LES, and immersed boundary (IB) techniques. Lee et al. [9, 10] also made computations on the micro VGs problems by using Monotone Integrated Large Eddy Simulations (MILES). Basic flow structures like momentum deficit and streamwise vortices were reproduced in the computation. Further studies were also conducted on the improvement of the control effect.Finding physics of MVG for design engineers is definitely needed. RANS, DES, RANS/DES, RANS/LES, and so forth are good engineering tools, but may not be able to reveal the mechanism and get deep understanding of MVG. We need high-order DNS/LES. A powerful tool is the integration of high-order LES and experiment. Recently, an implicit LES for MVG to control the shock-boundary layer interaction around a 24 degree ramp at Mach number of 2.5 and Reynolds number of 1440 has been carried out in 2009 [12]. In this work, several new discoveries have been made and the mechanism of MVG is found quite different from those reported by previous experimental and numerical work in the literature. The new findings include (1) spiral points and flow topology around MVG, (2) new theory of five pairs of vortices near MVG, (3) origin of momentum deficit, (4) inflection points (surface for 3D) and Kelvin-Helmholtz- (K-H-) type instability, (5) vortex ring generation, (6) ring-shock interaction and separation reduction, (7) 3D re-compressed shock structure, and (8) effects of trailing-edge decline  angles.Most of the new discoveries were made by October 2009 and were later confirmed by the UTA experimental work in April 2010 [13]. Based on our observation, the mechanism of MVG to reduce the flow separation is really caused by interaction of shock and vortex rings generated by MVG.The paper is arranged as following. The numerical methods are discussed in Section2. In Section 3, all new discoveries are summarized and analyzed and the new mechanism is described. Finally, conclusions are made in Section 4.


## 2. Numerical Methods

### 2.1. Governing Equations

The governing equations are the nondimensional Navier-Stokes equations in conservative form as follows:(1)∂Q∂t+∂E∂x+∂F∂y+∂G∂z=∂Ev∂x+∂Fv∂y+∂Gv∂z,
where(2)Q=[ρρuρvρwe],E=[ρuρu2+pρuvρuw(e+p)u],F=[ρvρvuρv2+pρvw(e+p)v],G=[ρwρwuρwvρw2+p(e+p)w],Ev=1Re[0τxxτxyτxzuτxx+vτxy+wτxz+qx],Fv=1Re[0τyxτyyτyzuτyx+vτyy+wτyz+qy],Gv=1Re[0τzxτzyτzzuτzx+vτzy+wτzz+qz],e=pγ-1+12ρ(u2+v2+w2),qx=μ(γ-1)M∞2Pr∂T∂x,qy=μ(γ-1)M∞2Pr∂T∂y,qz=μ(γ-1)M∞2Pr∂T∂z,p=1γM∞2ρT,Pr=0.72,τ=μ[43∂u∂x-23(∂v∂y+∂w∂z)∂u∂y+∂v∂x∂u∂z+∂w∂x∂u∂y+∂v∂x43∂v∂y-23(∂w∂z+∂u∂x)∂v∂z+∂w∂y∂u∂z+∂w∂x∂v∂z+∂w∂y43∂w∂z-23(∂u∂x+∂v∂y)]The viscous coefficient is given by Sutherland's equation:(3)μ=T3/21+CT+C,C=110.4T∞.
The nondimensional variables are defined as follows:(4)x=x̃L,y=ỹL,z=z̃L,u=ũU∞,v=ṽU∞,w=w̃U∞,T=T̃T∞,μ=μ̃μ∞,k=k̃k∞,ρ=ρ̃ρ∞,p=p̃ρ∞U∞2,e=ẽρ∞U∞2,
where the variables with “~” are the dimensional counterparts.Considering the following grid transformation,(5)ξ=ξ(x,y,z),η=η(x,y,z),ζ=ζ(x,y,z),
the Navier-Stokes equations can be transformed to the system using generalized coordinates:(6)∂Q̂∂τ+∂Ê∂ξ+∂F̂∂η+∂Ĝ∂ζ=∂Êv∂ξ+∂F̂v∂η+∂Ĝv∂ζ,
where Q̂=J-1Q and(7)Ê=J-1(ξxE+ξyF+ξzG),F̂=J-1(ηxE+ηyF+ηzG),Ĝ=J-1(ζxE+ζyF+ζzG),Êv=J-1(ξxEv+ξyFv+ξzGv),F̂v=J-1(ηxEv+ηyFv+ηzGv),Ĝv=J-1(ζxEv+ζyFv+ζzGv),J-1,ξx, and so forth, are grid metrics, and J-1=det(∂(x,y,z)/∂(ξ,η,ζ)), ξx=J(yηzζ-zηyζ), and so forth.


### 2.2. Finite Difference Schemes and Boundary Conditions

1. The Fifth-Order WENO Scheme for the Convective Terms [14]
In order to decrease the dissipation of the scheme, the less dissipative Steger-Warming flux splitting method is used in the computation, but not the commonly used dissipative Lax-Friedrich splitting method.2. The Difference Scheme for the Viscous Terms
Considering the conservative form of the governing equations, the traditional fourth-order central scheme is used twice to compute the second-order derivatives for viscous terms.3. The Time Schems
The basic methodology for the temporal terms in the Navier-Stokes equations adopts the explicit third-order TVD-type Runge-Kutta scheme:(8)u(1)=un+ΔtL(un),u(2)=34un+14u(1)+14ΔtL(u(1)),un+1=13un+23u(2)+23ΔtL(u(2)).4. Boundary Conditions
The adiabatic, zero-gradient of pressure and nonslipping conditions are used for the wall as(9)∂T∂n=0,∂p∂n=0,U⃗=0.
To enforce the free stream condition, fixed-value boundary condition with the free parameters is used on the upper boundary. The boundary conditions at the front and back boundary surface in the spanwise direction are given as the mirror-symmetry condition. The outflow boundary conditions are specified as a kind of characteristic-based condition, which can handle the outgoing flow without reflection. The inflow conditions are generated using the following steps.(a)
A turbulent mean profile is obtained from DNS by Liu and Chen [15] for the streamwise velocity (w-velocity) and the distribution is scaled using the local displacement thickness and free stream velocity.(b)
Random fluctuations are added on the primitive variables, that is,u, v, w, p, ρ. The disturbance has the form: εdistbe-(y-yw)2/Δydistb×(random-0.5)/2, where the subscript “distb” means the disturbance, “random” is the random function with the value between 0~1, εdistb equals to 0.1, and Δydistb equals to 2/3δ0.Such inflow conditions are, of course, not the fully developed turbulent flow, but we can consider it as a weakly disturbed inflow while propagating downstream.5. Body-Fitted Grid Generation
The geometry of MVG is shown in Figure1. In order to alleviate the difficulty to grid generation caused by original vertical trailing-edge, a modification is made by declining the edge to 70°. The other geometric parameters in the figure are the same as those given by Holden and Babinsky [5], that is, c=7.2h, α = 24°, and s=7.5h, where h is the height of MVG and s is the distance between the center lines of two adjacent MVGs. So the distance from the center line to the spanwise boundary of the computation domain is 3.75h.According to experiments by Holden and Babinsky [5], the ratio h/δ0of the models has the range from 0.3~1. The appropriate distance from the trailing-edge to the control area is around 19~56h or 8~19δ0. So in this study, the height of MVG h is assumed to be δ0/2 and the horizontal distance from the apex of MVG to the ramp corner is set to be 19.5h or 9.75δ0. The distance from the end of the ramp to the apex is 32.2896h. The distance from the starting point of the domain to the apex of MVG is 17.7775h. The height of the domain is from 10h to 15h and the width of the half domain is 3.75h. The geometric relation of the half of the domain can be seen in Figure 2, where the symmetric plane is the centre plane.The grid number for the whole system is nspanwise×nnormal×nstreamwise = 128 × 192 × 1600. We try to make the grids smooth and orthogonal as much as possible (Figure 3).Figure 1
The geometry of MVG.Figure 2
The schematic of the half-grid system of case (1).Body-fitted grid system.
(a)
Computational domain(b)
The grid section intersecting the middle of MVG(c)
The grids at the corner of MVG


### 2.3. Code Validation

1. Supersonic Ramp
Figure4 gives the instantaneous numerical schlieren image of the central plane, which uses the value of |∇ρ|. For qualitative comparison, an experimental picture of ramp flow at M=2.9 with large Reynolds number is given in Figure 5 obtained from the experiment by Loginov et al. [11]. From both pictures, we can observe that the separation shock wave has a declining angle which is a little lager but nearly same to the ramp angle and is almost aligned with the reflection shock (which is more obvious in the picture of the experiment); complex compressive waves are generated underneath the shock by vortices in the boundary layer.Figure 4
The digital schlieren at central plane.Figure 5
The experimental schlieren by Loginov et al. [11].2. Supersonic MVG
Based on the schlieren picture of the experiment, Babinsky et al. [7] presented the structure of the wave system, that is, the first reflection shock, the expansion wave system, and the recompression shock, as shown in Figure 6. The first shock will bend within the incoming boundary layer due to “the change in Mach number”, and “a second shock wave then turns the flow back to horizontal”. A subtle oblique “λ” structure is observed at the foot of the recompression shock from the picture. The numerical schlieren picture of the time-averaged flow field at the central plane is presented in Figure 7. From the figure, the two shock waves are well described by the computation. The bending of the main shock to the wall can be distinguished; the oblique “λ” structure can be found on the foot of the second shock as well. Table 1 gives a comparison between numerical and experimental results. The measured shock angles from Babinsky's experiment are 26.869° for the first shock and 21.93° for the recompression shock. In contrast, the computational values are 26.988° for the first shock, and 24.656° for the recompression shock. Considering the declining angle of the trailing-edge is 70°, not the 90° as in the experiment, the computational results show a reasonable agreement with experiment and of a high degree of resolution.Table 1
Comparison between IlES and experiment in shock angles.
1st shock2nd shockExperiment26.869°21.93°Computation26.988°24.656°Figure 6
The experimental schlieren by Babinsky et al. [7].Figure 7
The numerical schlieren at central plane.


## 2.1. Governing Equations

The governing equations are the nondimensional Navier-Stokes equations in conservative form as follows:(1)∂Q∂t+∂E∂x+∂F∂y+∂G∂z=∂Ev∂x+∂Fv∂y+∂Gv∂z,
where(2)Q=[ρρuρvρwe],E=[ρuρu2+pρuvρuw(e+p)u],F=[ρvρvuρv2+pρvw(e+p)v],G=[ρwρwuρwvρw2+p(e+p)w],Ev=1Re[0τxxτxyτxzuτxx+vτxy+wτxz+qx],Fv=1Re[0τyxτyyτyzuτyx+vτyy+wτyz+qy],Gv=1Re[0τzxτzyτzzuτzx+vτzy+wτzz+qz],e=pγ-1+12ρ(u2+v2+w2),qx=μ(γ-1)M∞2Pr∂T∂x,qy=μ(γ-1)M∞2Pr∂T∂y,qz=μ(γ-1)M∞2Pr∂T∂z,p=1γM∞2ρT,Pr=0.72,τ=μ[43∂u∂x-23(∂v∂y+∂w∂z)∂u∂y+∂v∂x∂u∂z+∂w∂x∂u∂y+∂v∂x43∂v∂y-23(∂w∂z+∂u∂x)∂v∂z+∂w∂y∂u∂z+∂w∂x∂v∂z+∂w∂y43∂w∂z-23(∂u∂x+∂v∂y)]The viscous coefficient is given by Sutherland's equation:(3)μ=T3/21+CT+C,C=110.4T∞.
The nondimensional variables are defined as follows:(4)x=x̃L,y=ỹL,z=z̃L,u=ũU∞,v=ṽU∞,w=w̃U∞,T=T̃T∞,μ=μ̃μ∞,k=k̃k∞,ρ=ρ̃ρ∞,p=p̃ρ∞U∞2,e=ẽρ∞U∞2,
where the variables with “~” are the dimensional counterparts.Considering the following grid transformation,(5)ξ=ξ(x,y,z),η=η(x,y,z),ζ=ζ(x,y,z),
the Navier-Stokes equations can be transformed to the system using generalized coordinates:(6)∂Q̂∂τ+∂Ê∂ξ+∂F̂∂η+∂Ĝ∂ζ=∂Êv∂ξ+∂F̂v∂η+∂Ĝv∂ζ,
where Q̂=J-1Q and(7)Ê=J-1(ξxE+ξyF+ξzG),F̂=J-1(ηxE+ηyF+ηzG),Ĝ=J-1(ζxE+ζyF+ζzG),Êv=J-1(ξxEv+ξyFv+ξzGv),F̂v=J-1(ηxEv+ηyFv+ηzGv),Ĝv=J-1(ζxEv+ζyFv+ζzGv),J-1,ξx, and so forth, are grid metrics, and J-1=det(∂(x,y,z)/∂(ξ,η,ζ)), ξx=J(yηzζ-zηyζ), and so forth.


## 2.2. Finite Difference Schemes and Boundary Conditions

1. The Fifth-Order WENO Scheme for the Convective Terms [14]
In order to decrease the dissipation of the scheme, the less dissipative Steger-Warming flux splitting method is used in the computation, but not the commonly used dissipative Lax-Friedrich splitting method.2. The Difference Scheme for the Viscous Terms
Considering the conservative form of the governing equations, the traditional fourth-order central scheme is used twice to compute the second-order derivatives for viscous terms.3. The Time Schems
The basic methodology for the temporal terms in the Navier-Stokes equations adopts the explicit third-order TVD-type Runge-Kutta scheme:(8)u(1)=un+ΔtL(un),u(2)=34un+14u(1)+14ΔtL(u(1)),un+1=13un+23u(2)+23ΔtL(u(2)).4. Boundary Conditions
The adiabatic, zero-gradient of pressure and nonslipping conditions are used for the wall as(9)∂T∂n=0,∂p∂n=0,U⃗=0.
To enforce the free stream condition, fixed-value boundary condition with the free parameters is used on the upper boundary. The boundary conditions at the front and back boundary surface in the spanwise direction are given as the mirror-symmetry condition. The outflow boundary conditions are specified as a kind of characteristic-based condition, which can handle the outgoing flow without reflection. The inflow conditions are generated using the following steps.(a)
A turbulent mean profile is obtained from DNS by Liu and Chen [15] for the streamwise velocity (w-velocity) and the distribution is scaled using the local displacement thickness and free stream velocity.(b)
Random fluctuations are added on the primitive variables, that is,u, v, w, p, ρ. The disturbance has the form: εdistbe-(y-yw)2/Δydistb×(random-0.5)/2, where the subscript “distb” means the disturbance, “random” is the random function with the value between 0~1, εdistb equals to 0.1, and Δydistb equals to 2/3δ0.Such inflow conditions are, of course, not the fully developed turbulent flow, but we can consider it as a weakly disturbed inflow while propagating downstream.5. Body-Fitted Grid Generation
The geometry of MVG is shown in Figure1. In order to alleviate the difficulty to grid generation caused by original vertical trailing-edge, a modification is made by declining the edge to 70°. The other geometric parameters in the figure are the same as those given by Holden and Babinsky [5], that is, c=7.2h, α = 24°, and s=7.5h, where h is the height of MVG and s is the distance between the center lines of two adjacent MVGs. So the distance from the center line to the spanwise boundary of the computation domain is 3.75h.According to experiments by Holden and Babinsky [5], the ratio h/δ0of the models has the range from 0.3~1. The appropriate distance from the trailing-edge to the control area is around 19~56h or 8~19δ0. So in this study, the height of MVG h is assumed to be δ0/2 and the horizontal distance from the apex of MVG to the ramp corner is set to be 19.5h or 9.75δ0. The distance from the end of the ramp to the apex is 32.2896h. The distance from the starting point of the domain to the apex of MVG is 17.7775h. The height of the domain is from 10h to 15h and the width of the half domain is 3.75h. The geometric relation of the half of the domain can be seen in Figure 2, where the symmetric plane is the centre plane.The grid number for the whole system is nspanwise×nnormal×nstreamwise = 128 × 192 × 1600. We try to make the grids smooth and orthogonal as much as possible (Figure 3).Figure 1
The geometry of MVG.Figure 2
The schematic of the half-grid system of case (1).Body-fitted grid system.
(a)
Computational domain(b)
The grid section intersecting the middle of MVG(c)
The grids at the corner of MVG


## 2.3. Code Validation

1. Supersonic Ramp
Figure4 gives the instantaneous numerical schlieren image of the central plane, which uses the value of |∇ρ|. For qualitative comparison, an experimental picture of ramp flow at M=2.9 with large Reynolds number is given in Figure 5 obtained from the experiment by Loginov et al. [11]. From both pictures, we can observe that the separation shock wave has a declining angle which is a little lager but nearly same to the ramp angle and is almost aligned with the reflection shock (which is more obvious in the picture of the experiment); complex compressive waves are generated underneath the shock by vortices in the boundary layer.Figure 4
The digital schlieren at central plane.Figure 5
The experimental schlieren by Loginov et al. [11].2. Supersonic MVG
Based on the schlieren picture of the experiment, Babinsky et al. [7] presented the structure of the wave system, that is, the first reflection shock, the expansion wave system, and the recompression shock, as shown in Figure 6. The first shock will bend within the incoming boundary layer due to “the change in Mach number”, and “a second shock wave then turns the flow back to horizontal”. A subtle oblique “λ” structure is observed at the foot of the recompression shock from the picture. The numerical schlieren picture of the time-averaged flow field at the central plane is presented in Figure 7. From the figure, the two shock waves are well described by the computation. The bending of the main shock to the wall can be distinguished; the oblique “λ” structure can be found on the foot of the second shock as well. Table 1 gives a comparison between numerical and experimental results. The measured shock angles from Babinsky's experiment are 26.869° for the first shock and 21.93° for the recompression shock. In contrast, the computational values are 26.988° for the first shock, and 24.656° for the recompression shock. Considering the declining angle of the trailing-edge is 70°, not the 90° as in the experiment, the computational results show a reasonable agreement with experiment and of a high degree of resolution.Table 1
Comparison between IlES and experiment in shock angles.
1st shock2nd shockExperiment26.869°21.93°Computation26.988°24.656°Figure 6
The experimental schlieren by Babinsky et al. [7].Figure 7
The numerical schlieren at central plane.


## 3. New Discoveries

### 3.1. Discovery I: Spiral Points and Vortex Structure around MVG

A pair of spiral points on the bottom plate and another pair on the sides of MVG (Figure8) have been found by ILES. These spiral points indicate the rising from the wall surface of the secondary vortex pairs. These new findings have been confirmed by UTA oil flow experiment (Figure 9), which shows accumulation blocks of paint. A further comparison is given in Figure 10.Figure 8
Spiral point pairs by ILES in 2009.Figure 9
Oil accumulation blocks by UTA experiment in 2010.Comparison of ILES with UTA experiment.
(a)(b)


### 3.2. Discovery II: New Model of Five Pairs of Vortices near MVG

The traditional vortex structure around MVG is depicted in Figure11 by Babinsky et al. [7], that is, a pair of counter-rotating primary streamwise vortices is generated by MVG, which is mainly located within the boundary layer and travel downstream for a considerable distance. Secondary vortices are located underneath the primary ones and even more streamwsie vortices could be generated under suitable conditions. However, a new model of 5 pairs of vortex tubes near MVG is given by us in Figure 12. The main difference between the two models is about the secondary vortex structures. There are two pairs of rising secondary vortex tubes corresponding to two pairs of spiral points. After the two pairs of secondary vortex tubes rising, a new pair of the secondary tubes is induced by the primary vortex tubes. Of course, there is additional one pair of horseshoe vortices. Recently, in April 2010, the experimentalists in the University of Texas at Arlington (UT Arlington) took a video from the top view to record the process of the MVG oil flow. A distinct pair of oil accumulating points was found in the video (see the experimental snapshot in Figures 9 and 10).Figure 11
Traditional vortex structure around MVG.Figure 12
New 5-pair-vortex-tube model near MVG.


### 3.3. Discovery III: Surface Separation Topology

For comparison, the pictures of the surface oil flow from experiment and limiting streamlines from the computation are given in Figures13 and 14. The topology of the two figures is essentially the same. The separation lines of the horseshoe vortex, the secondary separation lines beside the MVG, and the ones after the MVG are clearly described by computation (Figures 15 and 16). Such lines are the most obvious traces found in experiment due to the deposit of the oil. Figure 17 provides a comparison of computation and experiment in surface topology. The black area in experiment shows no oil flow from MVG top surface, which is corresponding to flow separation in computation. The agreement between computation and experiment is well obtained.Figure 13
The surface oil flow by Babinsky7.Figure 14
The surface limiting streamlines by computation.Figure 15
The separation pattern near the end of MVG.The various views of surface separation pattern (SSP).
(a)
SSP on the side of MVG(b)
SSP backgrounded using the pressure contourComparison of ILES and experiment in surface flow topology.
(a)
Computation:l1/l = 0.57(b)
Experiment:l1/l = 0.545The original secondary vortex tubes induced by the primary vortex tubes on the bottom plate and MVG sides were lifted up from the spiral points and further move downstream around the primary (Figures18 and 19). After that, a pair of new secondary vortex tubes is induced by the primary vortex tubes (Figure 20).Secondary vortex tubes lift up from the bottom plate.
(a)
Secondary vortex lifts up(b)
Moving downstream around the primarySecondary vortex tubes lift up from MVG sides.
(a)
Secondary vortex lifts up(b)
Moving downstream around the primaryFigure 20
New secondary vortex tubes induced by the primary.


### 3.4. Discovery IV: Origin of Momentum Deficit

The momentum deficit is a unique phenomenon which was first observed in the experiment by Ford and Babinsky [6], and later confirmed by the computation of Ghosh et al. [8] and Lee et al. [10], as shown in Figure 21, using the averaged streamwise velocity at different sections. About the mechanism of the deficit, Babinsky thought it is the wake of MVG, and Lee only mentioned “the two vortical tubes merge together to create a larger tube with two counter-rotating vortices inside”. What is the relation of the deficit to the flow structure? What does the low-speed flow come from? The existing explanations did not provide a clear mechanism about the formation of the deficit.Figure 21
Momentum deficit shown in different streamwise directions.In order to investigate the origins of the deficit, a heuristic analysis is made by investigating specific streamlines that are defined by certain cross-section of the momentum deficit. First, a cross-section after the MVG is selected as a reference plane as shown in Figure22, and the momentum deficit on the plane is shown by using the contour of the instantaneous streamwise velocity, which appears in a green circular region. Next, we put some seeds around the boundary of the deficit and within the deficit area. Then we draw the 3D streamlines starting from the seeds in both backward and forward directions. Because the seeds surround or locate within the deficit, the distribution of the backward streamlines qualitatively reflects the origins of the deficit. It can be seen from Figure 22 that (a) all backward streamlines rotate around the primary streamwise vortex, which indicates that the formation of the deficit is caused by primary vortices; (b) all backwards streamlines come from the upper surface of the MVG, which indicates that the main source of the deficit is the shedding of the upstream low-speed boundary layer, but is not the boundary layer flow around the MVG or after the MVG.Origin of the momentum deficit comes from the upper surface flow of MVG entrained by the primary vortex tubes.
(a)
Stream track for momentum deficit origin(b)
Locally enlarged stream tracks


### 3.5. Discovery V: Inflection Surface, K-H Instability, and Vortex Ring Generation by MVG

#### 3.5.1. 2D Kelvin-Helmoholz Instability

2D instability caused by shear layer with inflow disturbances and formation of pairing vortex rings was obtained by our previous calculation (Figure23). It is usually called Kelvin-Helmholtz instability which is described as inviscid instability.2D Kelvin-Helmoholzs instability and vortex ring formation.
(a)(b)(c)(d)


#### 3.5.2. Inflection Surface in 3D Flow behind MVG

In order to explore the mechanism of the vortex ring generation, the distributions of averaged streamwise velocity are given in Figure24 along the normal grid lines at the center plane. The streamwise positions of the lines are Lfromapex/h≈3.3, 6.7, 10, and 11, where Lfromapex is the streamwise distance measured from the apex of MVG. The dip of the lines corresponds to the momentum deficit. From the results, it can be clearly seen that there are at least two high-shear layers in the central plane: one is located at the upper edge of the dip and the other is located at the lower edge. Within the shear layer, there is at least one inflection point. In order to demonstrate the existence of the inflection points, the second-order derivative ∂2w/∂y2 (w is the streamwise velocity and y is the normal direction) is calculated along the lines, and the result of the line at Lfromapex/h≈ 6.7 is plotted in Figure 25 as an example. The existence and correspondence of the inflection points at the upper and lower shear layers is illustrated by two dashed lines intersecting the distribution of the streamwise velocity and its second-order derivative.Figure 24
Averaged streamwise velocity at different sections.Inflection points (surface for 3D).
(a)
Averaged streamwise velocity atLfromapex/h≈6.7(b)
∂2u/∂y2 = 0


#### 3.5.3. Vortex Ring Generation

Based on the above analyses, it can be concluded that the existence of the inflection points (surface in 3D) in the shear layer causes the flow instability and generates vortex rollers by K-H instability in a cylindrical coordinate system. Therefore, the mechanism of the vortex ring generation should be K-H instability. The lost of the stability of the shear layer will result in the roll-up of the vortex, which appears in ringlike structure in a 3D view (Figure26). In the Figure 26(b), λ2 is a certain eigenvalue of the stress tensor, and its iso-surface is usually used to describe the vortex surface. The intensity of the upper shear layer appears to be stronger than that of the lower shear layer (Figure 26).Vortex ring generation by MVG due to K-H instability.
(a)
Iso-surface of pressure(b)
Iso-surface ofλ2In Figure27, another qualitative checking about the shear layer and K-H instability is made by using the instantaneous flow field. In the figure, the background at the central plane and the spanwise plane is colored by the value of the streamwise velocity, so that the green regions in two planes represent the momentum deficit. In the central plane, the pressure contours are superimposed on the background cloud-map. The figure shows that the blue circle structures, which indicate the core of the ringlike vortices cut by the plane, are located on the boundary of the deficit circle. Such positions are exactly the same place where the high-shear layer exists.Figure 27
The instantaneous pressure and streamwise velocity contour on different cross-sections.In Figure28, we give the instantaneous numerical schlieren picture at the central plane. From the figure, we can see many vortex rings appear in the circular shapes; after being told the prediction of the vortex rings, the same experimentalists in UT Arlington tried some techniques to validate the discovery. They used techniques of the particle image velocimetry (PIV) and the acetone vapor screen visualization to track the movement of the flow, and specifically the flash of a laser sheet is used to provide the light exposure at the time level of micro seconds. In Figure 29, a typical image at the center plane is presented taken by using PIV and the acetone vapor [13]. It is clearly demonstrated that a chain of vortex rings exist in the flow field after the MVG! And these structures qualitatively resemble those in Figures 27 and 28.Figure 28
The numerical schlieren at the center plane.The laser-sheet flash image at the center plane.
(a)
Using PIV(b)
Using the acetone vapor


### 3.6. Discovery VI: Reduction of Separation Zone by Ring-Shock Interaction

#### 3.6.1. New Mechanism of Separation Zone Reduction by MVG

Traditional concept for MVG is that the MVG produces streamwise vortices inside the boundary later, which will bring low-momentum fluid up from the bottom and high-momentum fluid down to the boundary layer and then the flow separation induced by shock will be reduced due to the higher shear. However, our LES just gives a quite different mechanism which can be described as MVG generates vortex rings and vortex rings destroy the shocks and reduce the boundary layer separation bubble size, which was induced by the separation bubble and reflection shocks.Figure30 shows a graphic about the shock structure around a supersonic ramp, which includes a separation bubble shock and a reflection shock. The boundary layer separation is caused by the strong adverse pressure gradients induced by the shock. Many vortex rings are generated by MVG and convect downstream toward the ramp. The ring first destroys the separation bubble shock and the separation bubble looks broken (Figure 31). The vortex rings are very stable and continue to move downstream to destroy the reflection shock as shown in Figure 32. After many rings continuously move toward the reflection shock, the shock tail disappeared and the separation area is substantially reduced (Figure 33).Figure 30
Shock structure around ramp.Figure 31
Ring-separation bubble shock interaction.Figure 32
Ring-reflection shock interaction.Figure 33
Separation reduction by vortex rings.


#### 3.6.2. Shock-Vortex Ring Interaction

The mechanism why the vortex ring can destroy the shock is believed as that the vortex ring has a low pressure line inside the ring which could destroy the shock. On the other hand, the ring has a high-pressure area outside the ring and near the center of ring, which will change the shock location. The fast rotation of the ring could make the shock surface deformed like a bump.


#### 3.6.3. Stable Structure of the Vortex Ring

The shock is pretty strong and the pressure gradients are pretty large. We use both digital schlieren (Figure34) and λ2 techniques (Figure 35) for 3D visualization and 2D visualization on a central plane in the spanwise direction. The vortex rings are found very stable and keep their forms unchanged after passing the weak separation bubble shock and strong ramp shock.Vortex ring-shock interaction at the central plane (numerical schlieren).
(a)
Before the separation shock(b)
Past the separation shock(c)
Past the ramp shockVortex ring-shock interaction at the central plane (λ2 techniques).
(a)
Before the shock(b)
Hit the shock(c)
Past the shock


### 3.7. Discovery VII: Three-Dimensional Recompressed Shock Structure

In Figure36, the spatial structures of the wave system are given by instantaneous cross-section pressure contours. From the first section in Figure 36(a), the wave system is found to be the main reflection shock above the MVG and the expansion wave emanating from the edge of the MVG. The remaining sections show an arclike structure of the recompression shock wave, which is located behind the first reflection shock. The size of the arc keeps growing while moving downstream. The curve shape of the recompression shocks implies that the function of the recompression shock wave is to make the expansion flow satisfy the “virtual” boundary condition created by the streamwise vortices, at least at the initial stage of the formation of the shock wave. Detailed investigations show that the head and feet of the arc-like shock are separated from each other in the initial stage (see Figure 36(b)). They begin to connect and form a complete arc curve at a certain downstream location, as shown in Figure 36(c).The 3D structure of the recompression shock.
(a)
Section  1(b)
Section  2(c)
Overview


### 3.8. Discovery VIII: Effects of Trailing-Edge Decline Angles

Comparative studies are made on the MVG flow atM=2.5 and Reθ=1440 while the MVGs are of different declining angle, that is, β = 70° and β = 45°. The numerical study finds (see Figure 37) the following(1)
The basic structures are similar in two cases, like the wave system, separation pattern, momentum deficit, and the vortex rings. The effect of declining angle does not bring about the structural difference, while differences exist in the details of the flow.(2)
The smaller declining angle atβ = 45° makes the vortices closer to the wall, which include the initial streamwise vortices and the layer vortex rings. This should be more favorable to flow control. Nevertheless, the quantitative analyses on the time-averaged data at a downstream do not exhibit a significant improvement about the quality of the boundary layer.Comparison of the vortex structure with different decline angles of MVG trailing edge.
(a)
β = 70°(b)
β = 45°


## 3.1. Discovery I: Spiral Points and Vortex Structure around MVG

A pair of spiral points on the bottom plate and another pair on the sides of MVG (Figure8) have been found by ILES. These spiral points indicate the rising from the wall surface of the secondary vortex pairs. These new findings have been confirmed by UTA oil flow experiment (Figure 9), which shows accumulation blocks of paint. A further comparison is given in Figure 10.Figure 8
Spiral point pairs by ILES in 2009.Figure 9
Oil accumulation blocks by UTA experiment in 2010.Comparison of ILES with UTA experiment.
(a)(b)


## 3.2. Discovery II: New Model of Five Pairs of Vortices near MVG

The traditional vortex structure around MVG is depicted in Figure11 by Babinsky et al. [7], that is, a pair of counter-rotating primary streamwise vortices is generated by MVG, which is mainly located within the boundary layer and travel downstream for a considerable distance. Secondary vortices are located underneath the primary ones and even more streamwsie vortices could be generated under suitable conditions. However, a new model of 5 pairs of vortex tubes near MVG is given by us in Figure 12. The main difference between the two models is about the secondary vortex structures. There are two pairs of rising secondary vortex tubes corresponding to two pairs of spiral points. After the two pairs of secondary vortex tubes rising, a new pair of the secondary tubes is induced by the primary vortex tubes. Of course, there is additional one pair of horseshoe vortices. Recently, in April 2010, the experimentalists in the University of Texas at Arlington (UT Arlington) took a video from the top view to record the process of the MVG oil flow. A distinct pair of oil accumulating points was found in the video (see the experimental snapshot in Figures 9 and 10).Figure 11
Traditional vortex structure around MVG.Figure 12
New 5-pair-vortex-tube model near MVG.


## 3.3. Discovery III: Surface Separation Topology

For comparison, the pictures of the surface oil flow from experiment and limiting streamlines from the computation are given in Figures13 and 14. The topology of the two figures is essentially the same. The separation lines of the horseshoe vortex, the secondary separation lines beside the MVG, and the ones after the MVG are clearly described by computation (Figures 15 and 16). Such lines are the most obvious traces found in experiment due to the deposit of the oil. Figure 17 provides a comparison of computation and experiment in surface topology. The black area in experiment shows no oil flow from MVG top surface, which is corresponding to flow separation in computation. The agreement between computation and experiment is well obtained.Figure 13
The surface oil flow by Babinsky7.Figure 14
The surface limiting streamlines by computation.Figure 15
The separation pattern near the end of MVG.The various views of surface separation pattern (SSP).
(a)
SSP on the side of MVG(b)
SSP backgrounded using the pressure contourComparison of ILES and experiment in surface flow topology.
(a)
Computation:l1/l = 0.57(b)
Experiment:l1/l = 0.545The original secondary vortex tubes induced by the primary vortex tubes on the bottom plate and MVG sides were lifted up from the spiral points and further move downstream around the primary (Figures18 and 19). After that, a pair of new secondary vortex tubes is induced by the primary vortex tubes (Figure 20).Secondary vortex tubes lift up from the bottom plate.
(a)
Secondary vortex lifts up(b)
Moving downstream around the primarySecondary vortex tubes lift up from MVG sides.
(a)
Secondary vortex lifts up(b)
Moving downstream around the primaryFigure 20
New secondary vortex tubes induced by the primary.


## 3.4. Discovery IV: Origin of Momentum Deficit

The momentum deficit is a unique phenomenon which was first observed in the experiment by Ford and Babinsky [6], and later confirmed by the computation of Ghosh et al. [8] and Lee et al. [10], as shown in Figure 21, using the averaged streamwise velocity at different sections. About the mechanism of the deficit, Babinsky thought it is the wake of MVG, and Lee only mentioned “the two vortical tubes merge together to create a larger tube with two counter-rotating vortices inside”. What is the relation of the deficit to the flow structure? What does the low-speed flow come from? The existing explanations did not provide a clear mechanism about the formation of the deficit.Figure 21
Momentum deficit shown in different streamwise directions.In order to investigate the origins of the deficit, a heuristic analysis is made by investigating specific streamlines that are defined by certain cross-section of the momentum deficit. First, a cross-section after the MVG is selected as a reference plane as shown in Figure22, and the momentum deficit on the plane is shown by using the contour of the instantaneous streamwise velocity, which appears in a green circular region. Next, we put some seeds around the boundary of the deficit and within the deficit area. Then we draw the 3D streamlines starting from the seeds in both backward and forward directions. Because the seeds surround or locate within the deficit, the distribution of the backward streamlines qualitatively reflects the origins of the deficit. It can be seen from Figure 22 that (a) all backward streamlines rotate around the primary streamwise vortex, which indicates that the formation of the deficit is caused by primary vortices; (b) all backwards streamlines come from the upper surface of the MVG, which indicates that the main source of the deficit is the shedding of the upstream low-speed boundary layer, but is not the boundary layer flow around the MVG or after the MVG.Origin of the momentum deficit comes from the upper surface flow of MVG entrained by the primary vortex tubes.
(a)
Stream track for momentum deficit origin(b)
Locally enlarged stream tracks


## 3.5. Discovery V: Inflection Surface, K-H Instability, and Vortex Ring Generation by MVG

### 3.5.1. 2D Kelvin-Helmoholz Instability

2D instability caused by shear layer with inflow disturbances and formation of pairing vortex rings was obtained by our previous calculation (Figure23). It is usually called Kelvin-Helmholtz instability which is described as inviscid instability.2D Kelvin-Helmoholzs instability and vortex ring formation.
(a)(b)(c)(d)


### 3.5.2. Inflection Surface in 3D Flow behind MVG

In order to explore the mechanism of the vortex ring generation, the distributions of averaged streamwise velocity are given in Figure24 along the normal grid lines at the center plane. The streamwise positions of the lines are Lfromapex/h≈3.3, 6.7, 10, and 11, where Lfromapex is the streamwise distance measured from the apex of MVG. The dip of the lines corresponds to the momentum deficit. From the results, it can be clearly seen that there are at least two high-shear layers in the central plane: one is located at the upper edge of the dip and the other is located at the lower edge. Within the shear layer, there is at least one inflection point. In order to demonstrate the existence of the inflection points, the second-order derivative ∂2w/∂y2 (w is the streamwise velocity and y is the normal direction) is calculated along the lines, and the result of the line at Lfromapex/h≈ 6.7 is plotted in Figure 25 as an example. The existence and correspondence of the inflection points at the upper and lower shear layers is illustrated by two dashed lines intersecting the distribution of the streamwise velocity and its second-order derivative.Figure 24
Averaged streamwise velocity at different sections.Inflection points (surface for 3D).
(a)
Averaged streamwise velocity atLfromapex/h≈6.7(b)
∂2u/∂y2 = 0


### 3.5.3. Vortex Ring Generation

Based on the above analyses, it can be concluded that the existence of the inflection points (surface in 3D) in the shear layer causes the flow instability and generates vortex rollers by K-H instability in a cylindrical coordinate system. Therefore, the mechanism of the vortex ring generation should be K-H instability. The lost of the stability of the shear layer will result in the roll-up of the vortex, which appears in ringlike structure in a 3D view (Figure26). In the Figure 26(b), λ2 is a certain eigenvalue of the stress tensor, and its iso-surface is usually used to describe the vortex surface. The intensity of the upper shear layer appears to be stronger than that of the lower shear layer (Figure 26).Vortex ring generation by MVG due to K-H instability.
(a)
Iso-surface of pressure(b)
Iso-surface ofλ2In Figure27, another qualitative checking about the shear layer and K-H instability is made by using the instantaneous flow field. In the figure, the background at the central plane and the spanwise plane is colored by the value of the streamwise velocity, so that the green regions in two planes represent the momentum deficit. In the central plane, the pressure contours are superimposed on the background cloud-map. The figure shows that the blue circle structures, which indicate the core of the ringlike vortices cut by the plane, are located on the boundary of the deficit circle. Such positions are exactly the same place where the high-shear layer exists.Figure 27
The instantaneous pressure and streamwise velocity contour on different cross-sections.In Figure28, we give the instantaneous numerical schlieren picture at the central plane. From the figure, we can see many vortex rings appear in the circular shapes; after being told the prediction of the vortex rings, the same experimentalists in UT Arlington tried some techniques to validate the discovery. They used techniques of the particle image velocimetry (PIV) and the acetone vapor screen visualization to track the movement of the flow, and specifically the flash of a laser sheet is used to provide the light exposure at the time level of micro seconds. In Figure 29, a typical image at the center plane is presented taken by using PIV and the acetone vapor [13]. It is clearly demonstrated that a chain of vortex rings exist in the flow field after the MVG! And these structures qualitatively resemble those in Figures 27 and 28.Figure 28
The numerical schlieren at the center plane.The laser-sheet flash image at the center plane.
(a)
Using PIV(b)
Using the acetone vapor


## 3.5.1. 2D Kelvin-Helmoholz Instability

2D instability caused by shear layer with inflow disturbances and formation of pairing vortex rings was obtained by our previous calculation (Figure23). It is usually called Kelvin-Helmholtz instability which is described as inviscid instability.2D Kelvin-Helmoholzs instability and vortex ring formation.
(a)(b)(c)(d)


## 3.5.2. Inflection Surface in 3D Flow behind MVG

In order to explore the mechanism of the vortex ring generation, the distributions of averaged streamwise velocity are given in Figure24 along the normal grid lines at the center plane. The streamwise positions of the lines are Lfromapex/h≈3.3, 6.7, 10, and 11, where Lfromapex is the streamwise distance measured from the apex of MVG. The dip of the lines corresponds to the momentum deficit. From the results, it can be clearly seen that there are at least two high-shear layers in the central plane: one is located at the upper edge of the dip and the other is located at the lower edge. Within the shear layer, there is at least one inflection point. In order to demonstrate the existence of the inflection points, the second-order derivative ∂2w/∂y2 (w is the streamwise velocity and y is the normal direction) is calculated along the lines, and the result of the line at Lfromapex/h≈ 6.7 is plotted in Figure 25 as an example. The existence and correspondence of the inflection points at the upper and lower shear layers is illustrated by two dashed lines intersecting the distribution of the streamwise velocity and its second-order derivative.Figure 24
Averaged streamwise velocity at different sections.Inflection points (surface for 3D).
(a)
Averaged streamwise velocity atLfromapex/h≈6.7(b)
∂2u/∂y2 = 0


## 3.5.3. Vortex Ring Generation

Based on the above analyses, it can be concluded that the existence of the inflection points (surface in 3D) in the shear layer causes the flow instability and generates vortex rollers by K-H instability in a cylindrical coordinate system. Therefore, the mechanism of the vortex ring generation should be K-H instability. The lost of the stability of the shear layer will result in the roll-up of the vortex, which appears in ringlike structure in a 3D view (Figure26). In the Figure 26(b), λ2 is a certain eigenvalue of the stress tensor, and its iso-surface is usually used to describe the vortex surface. The intensity of the upper shear layer appears to be stronger than that of the lower shear layer (Figure 26).Vortex ring generation by MVG due to K-H instability.
(a)
Iso-surface of pressure(b)
Iso-surface ofλ2In Figure27, another qualitative checking about the shear layer and K-H instability is made by using the instantaneous flow field. In the figure, the background at the central plane and the spanwise plane is colored by the value of the streamwise velocity, so that the green regions in two planes represent the momentum deficit. In the central plane, the pressure contours are superimposed on the background cloud-map. The figure shows that the blue circle structures, which indicate the core of the ringlike vortices cut by the plane, are located on the boundary of the deficit circle. Such positions are exactly the same place where the high-shear layer exists.Figure 27
The instantaneous pressure and streamwise velocity contour on different cross-sections.In Figure28, we give the instantaneous numerical schlieren picture at the central plane. From the figure, we can see many vortex rings appear in the circular shapes; after being told the prediction of the vortex rings, the same experimentalists in UT Arlington tried some techniques to validate the discovery. They used techniques of the particle image velocimetry (PIV) and the acetone vapor screen visualization to track the movement of the flow, and specifically the flash of a laser sheet is used to provide the light exposure at the time level of micro seconds. In Figure 29, a typical image at the center plane is presented taken by using PIV and the acetone vapor [13]. It is clearly demonstrated that a chain of vortex rings exist in the flow field after the MVG! And these structures qualitatively resemble those in Figures 27 and 28.Figure 28
The numerical schlieren at the center plane.The laser-sheet flash image at the center plane.
(a)
Using PIV(b)
Using the acetone vapor


## 3.6. Discovery VI: Reduction of Separation Zone by Ring-Shock Interaction

### 3.6.1. New Mechanism of Separation Zone Reduction by MVG

Traditional concept for MVG is that the MVG produces streamwise vortices inside the boundary later, which will bring low-momentum fluid up from the bottom and high-momentum fluid down to the boundary layer and then the flow separation induced by shock will be reduced due to the higher shear. However, our LES just gives a quite different mechanism which can be described as MVG generates vortex rings and vortex rings destroy the shocks and reduce the boundary layer separation bubble size, which was induced by the separation bubble and reflection shocks.Figure30 shows a graphic about the shock structure around a supersonic ramp, which includes a separation bubble shock and a reflection shock. The boundary layer separation is caused by the strong adverse pressure gradients induced by the shock. Many vortex rings are generated by MVG and convect downstream toward the ramp. The ring first destroys the separation bubble shock and the separation bubble looks broken (Figure 31). The vortex rings are very stable and continue to move downstream to destroy the reflection shock as shown in Figure 32. After many rings continuously move toward the reflection shock, the shock tail disappeared and the separation area is substantially reduced (Figure 33).Figure 30
Shock structure around ramp.Figure 31
Ring-separation bubble shock interaction.Figure 32
Ring-reflection shock interaction.Figure 33
Separation reduction by vortex rings.


### 3.6.2. Shock-Vortex Ring Interaction

The mechanism why the vortex ring can destroy the shock is believed as that the vortex ring has a low pressure line inside the ring which could destroy the shock. On the other hand, the ring has a high-pressure area outside the ring and near the center of ring, which will change the shock location. The fast rotation of the ring could make the shock surface deformed like a bump.


### 3.6.3. Stable Structure of the Vortex Ring

The shock is pretty strong and the pressure gradients are pretty large. We use both digital schlieren (Figure34) and λ2 techniques (Figure 35) for 3D visualization and 2D visualization on a central plane in the spanwise direction. The vortex rings are found very stable and keep their forms unchanged after passing the weak separation bubble shock and strong ramp shock.Vortex ring-shock interaction at the central plane (numerical schlieren).
(a)
Before the separation shock(b)
Past the separation shock(c)
Past the ramp shockVortex ring-shock interaction at the central plane (λ2 techniques).
(a)
Before the shock(b)
Hit the shock(c)
Past the shock


## 3.6.1. New Mechanism of Separation Zone Reduction by MVG

Traditional concept for MVG is that the MVG produces streamwise vortices inside the boundary later, which will bring low-momentum fluid up from the bottom and high-momentum fluid down to the boundary layer and then the flow separation induced by shock will be reduced due to the higher shear. However, our LES just gives a quite different mechanism which can be described as MVG generates vortex rings and vortex rings destroy the shocks and reduce the boundary layer separation bubble size, which was induced by the separation bubble and reflection shocks.Figure30 shows a graphic about the shock structure around a supersonic ramp, which includes a separation bubble shock and a reflection shock. The boundary layer separation is caused by the strong adverse pressure gradients induced by the shock. Many vortex rings are generated by MVG and convect downstream toward the ramp. The ring first destroys the separation bubble shock and the separation bubble looks broken (Figure 31). The vortex rings are very stable and continue to move downstream to destroy the reflection shock as shown in Figure 32. After many rings continuously move toward the reflection shock, the shock tail disappeared and the separation area is substantially reduced (Figure 33).Figure 30
Shock structure around ramp.Figure 31
Ring-separation bubble shock interaction.Figure 32
Ring-reflection shock interaction.Figure 33
Separation reduction by vortex rings.


## 3.6.2. Shock-Vortex Ring Interaction

The mechanism why the vortex ring can destroy the shock is believed as that the vortex ring has a low pressure line inside the ring which could destroy the shock. On the other hand, the ring has a high-pressure area outside the ring and near the center of ring, which will change the shock location. The fast rotation of the ring could make the shock surface deformed like a bump.


## 3.6.3. Stable Structure of the Vortex Ring

The shock is pretty strong and the pressure gradients are pretty large. We use both digital schlieren (Figure34) and λ2 techniques (Figure 35) for 3D visualization and 2D visualization on a central plane in the spanwise direction. The vortex rings are found very stable and keep their forms unchanged after passing the weak separation bubble shock and strong ramp shock.Vortex ring-shock interaction at the central plane (numerical schlieren).
(a)
Before the separation shock(b)
Past the separation shock(c)
Past the ramp shockVortex ring-shock interaction at the central plane (λ2 techniques).
(a)
Before the shock(b)
Hit the shock(c)
Past the shock


## 3.7. Discovery VII: Three-Dimensional Recompressed Shock Structure

In Figure36, the spatial structures of the wave system are given by instantaneous cross-section pressure contours. From the first section in Figure 36(a), the wave system is found to be the main reflection shock above the MVG and the expansion wave emanating from the edge of the MVG. The remaining sections show an arclike structure of the recompression shock wave, which is located behind the first reflection shock. The size of the arc keeps growing while moving downstream. The curve shape of the recompression shocks implies that the function of the recompression shock wave is to make the expansion flow satisfy the “virtual” boundary condition created by the streamwise vortices, at least at the initial stage of the formation of the shock wave. Detailed investigations show that the head and feet of the arc-like shock are separated from each other in the initial stage (see Figure 36(b)). They begin to connect and form a complete arc curve at a certain downstream location, as shown in Figure 36(c).The 3D structure of the recompression shock.
(a)
Section  1(b)
Section  2(c)
Overview


## 3.8. Discovery VIII: Effects of Trailing-Edge Decline Angles

Comparative studies are made on the MVG flow atM=2.5 and Reθ=1440 while the MVGs are of different declining angle, that is, β = 70° and β = 45°. The numerical study finds (see Figure 37) the following(1)
The basic structures are similar in two cases, like the wave system, separation pattern, momentum deficit, and the vortex rings. The effect of declining angle does not bring about the structural difference, while differences exist in the details of the flow.(2)
The smaller declining angle atβ = 45° makes the vortices closer to the wall, which include the initial streamwise vortices and the layer vortex rings. This should be more favorable to flow control. Nevertheless, the quantitative analyses on the time-averaged data at a downstream do not exhibit a significant improvement about the quality of the boundary layer.Comparison of the vortex structure with different decline angles of MVG trailing edge.
(a)
β = 70°(b)
β = 45°


## 4. Conclusions

The following new discoveries and new mechanisms have been made by implicit LES since 2009. Most of these new findings were confirmed later by UTA experiment in 2010.(1)
Detailed flow structures (topology) around MVG and spiral points are found. Then the new five pairs of vortex tubes model around MVG are given.(2)
Origin of the momentum deficit is found, not only caused by stream wise vortices through mushroom, but also, mainly, viscosity of MVG and wall surface.(3)
Inflection points (surface in 3D) inside the deficit area are found.(4)
Kelvin-Helmholtz-type instability is caused by the momentum deficit.(5)
Vortex rings are generated by K-H-type instability after MVG.(6)
Shock-ring interaction is found.(7)
Mechanism of MVG for separation reduction is due to the ring shock interaction, not only mixing as traditional explanation on VG or MVG.(8)
Ring is the only form of vortex inside the fluid field and never breaks down.(9)
Structure of recompressed shock is found.(10)
Influence of the MVG trailing edge decline angles is studied.All these new findings may be very important to regular VG, subsonic flow, flow control, flow transition, and turbulence studies.



---

*Source: 934982-2011-04-03.xml*