# Inviscid and Viscous Interactions in Subsonic Corner Flows

Authors: Kung-Ming Chung; Po-Hsiung Chang; Keh-Chin Chang
Journal: The Scientific World Journal
 (2013)

---

## Abstract

A flap can be used as a high-lift device, in which a downward deflection results in a gain in lift at a given geometric angle of attack. To characterize the aerodynamic performance of a deflected surface in compressible flows, the present study examines a naturally developed turbulent boundary layer past the convex and concave corners. This investigation involves the analysis of mean and fluctuating pressure distributions. The results obtained indicate strong inviscid-viscous interactions. There are upstream expansion and downstream compression for the convex-corner flows, while the opposite trend is observed for the concave-corner flows. A combined flow similarity parameter, based on the small perturbation theory, is proposed to scale the flow characteristics in both subsonic convex- and concave-corner flows.


---

## Body

## 1. Introduction

Corner flows occur in a wide variety of internal and external aerodynamic problems. Previous studies have been mainly on supersonic and hypersonic speeds [1–4]. In subsonic flow regime, aircraft designs have employed flaps for take-off and landing performance and ailerons for routine turning maneuver. A study by Bolonkin and Gilyard [5] demonstrated that active modification of control surfaces (variable camber wings) potentially could play a role in performance optimization for fighter aircraft and transport aircraft. While cruising, there could be more than 10 percent in maximizing the lift-to-drag ratio, especially for nonstandard flight conditions [6]. Further, a simplified model of a deflected surface comprised convex-corner and concave-corner flows. For a compressible convex-corner flow (or upper deflected surface), there are strong upstream expansion and downstream compression, caused by viscous-inviscid interactions, near the corner. The displacement thickness near the corner is affected by the overlapping region that lies between the viscous sublayer and the main part of the boundary layer [7]. On the lower deflected surface (or concave corner), the flow decelerates upstream of the corner followed by the downstream acceleration. Previous studies [8–10] demonstrated that the lift coefficient increases linearly with the deflection angle.Chung [11] demonstrated that M and η are the two major parameters affecting the type of flow field in compressible convex-corner flows. The hypersonic similarity parameter (Mη) [12] and a similar combined supersonic-hypersonic similarity parameter (1-M2η) [13] were examined for scaling the expansion flows. However, the parameter M2η appears to be more suitable to characterize the flow characteristics, including peak Mach number, interaction region, and amplitude of peak pressure fluctuations [11, 14]. Further, according to the hypothesis that small streamline deflections produce proportionally small change in Mach number and pressure, a hodograph solution for compressible flow past a corner was given by Verhoff et al. [15]. Another flow similarity parameter,  M2η/1-M2, was employed by Chung et al. [16] to categorize flow regimes of compressible convex-corner flows. In compressible concave-corner flows, Chung [17] demonstrated that the characteristics of the flow can be scaled with Mη, in which stronger upstream compression and downstream expansion are observed with increasing M and η. Note that there are only slight variations in surface pressure fluctuation coefficients.Small perturbation theories have been applied to a number of aeronautical problems, in which the flow is characterized by a small deviation of the flow from its original uniform flow. Linearized solutions are useful for explicitly identifying trends and governing parameters. Consider a slender body at hypersonic speeds, the surface pressure coefficient can be written in terms of hypersonic similarity parameter(K) and specific heat ratio (γ), in which [13]
(1)Cp=2η2γK2[(1-γ-12K)2γ/(γ-1)-1]
or
(2)Cpη2=f(K,γ).
As mentioned above, the parameters M2η/1-M2 and Mη were employed to characterize compressible convex- and concave-corner flows, respectively [16]. A common similarity parameter for subsonic corner flows is of interest. Therefore, the present study adopts a combined parameter K*(=M2/η1-M2) as a governing parameter to scale the flow characteristics in both subsonic convex- and concave-corner flows for rapid estimation of the interaction region, peak Mach number, and peak pressure fluctuations.


## 2. Experimental Setup

### 2.1. Transonic Wind Tunnel

The experiments were conducted at a transonic wind tunnel of blowdown type, located at Aerospace Science and Technology Research Center, National Cheng Kung University, Taiwan. Major components of the facility include compressors, air dryers, cooling water system, storage tanks, and the tunnel. The dew point of high-pressure air through the air dryers is maintained at −40°C under normal operation conditions. Operating Mach number ranges from 0.2 to 1.4, and simulated Reynolds number is up to 20 million per meter. The test section is 600 × 600 mm and 1500 mm long. In the present study, the test section was assembled with solid sidewalls and perforated top/bottom walls. The freestream Mach numbers were 0.33 and 0.64 ± 0.01, and the stagnation pressure and temperature were 172 ± 0.5 kPa and room temperature, respectively. For the data acquisition system, the LeCroy waveform recorders were used. A host computer with CATALYST software controlled the setup of LeCroy waveform recorders through a LeCroy 8901A interface. All input channels were triggered simultaneously.


### 2.2. Test Model

The test model consisted of a flat plate and an interchangeable instrumentation plate. The test model was 150 mm wide and 600 mm long, which was supported by a single sting mounted on the bottom wall of the test section, as shown in Figure1. The concave corner with 5-, 7-, 10-, and 15-deg angles or the convex corner with 5-, 10-, 13-, and 15-deg angles was located at 500 mm from the leading edge of the flat plate. Along the centerline of each instrumentation plate, 19 pressure taps (6 mm apart and 2.5 mm in diameter) were drilled perpendicular to the test surface. All the pressure transducers were flush-mounted to the test surface. The side fences of the instrumentation plate were also installed to prevent cross-flow. A study by Miau et al. [18] indicated that the transition of the boundary layer under the present test condition is close to the leading edge of the flat plate, indicating a turbulent boundary layer at the measurement locations. The boundary layer thickness, δ, at 25 mm upstream of the corner was approximately 7.0 mm.Figure 1
Experimental setup.


### 2.3. Experimental Techniques

The Kulite (Model XCS-093-25A, B screen) pressure transducers, which were powered by a TES Model 6102 power supply at 15.0 volts, were employed for pressure measurements. Their outer diameter and sensing element are 2.36 mm and 0.97 mm, respectively. Note that the perforated screen of the pressure transducers might limit the frequency response to only 50 kHz [19]. To improve the signal-to-noise ratio, external amplifiers (Ecreon Model E713) were also employed. With a gain of 20, the roll-off frequency is about 140 kHz. The typical sampling period is 5 μs (200 kHz). Each data record possesses 131,072 data points for statistical analysis. The data were divided into 32 blocks. The mean values of each block (4,096 data points) were calculated. Variations of the blocks are estimated to be 0.43 and 0.13 percent for the mean surface pressure coefficient, Cp, and the fluctuating pressure coefficient, Cσp, respectively, which were taken as uncertainty of the experimental data.


## 2.1. Transonic Wind Tunnel

The experiments were conducted at a transonic wind tunnel of blowdown type, located at Aerospace Science and Technology Research Center, National Cheng Kung University, Taiwan. Major components of the facility include compressors, air dryers, cooling water system, storage tanks, and the tunnel. The dew point of high-pressure air through the air dryers is maintained at −40°C under normal operation conditions. Operating Mach number ranges from 0.2 to 1.4, and simulated Reynolds number is up to 20 million per meter. The test section is 600 × 600 mm and 1500 mm long. In the present study, the test section was assembled with solid sidewalls and perforated top/bottom walls. The freestream Mach numbers were 0.33 and 0.64 ± 0.01, and the stagnation pressure and temperature were 172 ± 0.5 kPa and room temperature, respectively. For the data acquisition system, the LeCroy waveform recorders were used. A host computer with CATALYST software controlled the setup of LeCroy waveform recorders through a LeCroy 8901A interface. All input channels were triggered simultaneously.


## 2.2. Test Model

The test model consisted of a flat plate and an interchangeable instrumentation plate. The test model was 150 mm wide and 600 mm long, which was supported by a single sting mounted on the bottom wall of the test section, as shown in Figure1. The concave corner with 5-, 7-, 10-, and 15-deg angles or the convex corner with 5-, 10-, 13-, and 15-deg angles was located at 500 mm from the leading edge of the flat plate. Along the centerline of each instrumentation plate, 19 pressure taps (6 mm apart and 2.5 mm in diameter) were drilled perpendicular to the test surface. All the pressure transducers were flush-mounted to the test surface. The side fences of the instrumentation plate were also installed to prevent cross-flow. A study by Miau et al. [18] indicated that the transition of the boundary layer under the present test condition is close to the leading edge of the flat plate, indicating a turbulent boundary layer at the measurement locations. The boundary layer thickness, δ, at 25 mm upstream of the corner was approximately 7.0 mm.Figure 1
Experimental setup.


## 2.3. Experimental Techniques

The Kulite (Model XCS-093-25A, B screen) pressure transducers, which were powered by a TES Model 6102 power supply at 15.0 volts, were employed for pressure measurements. Their outer diameter and sensing element are 2.36 mm and 0.97 mm, respectively. Note that the perforated screen of the pressure transducers might limit the frequency response to only 50 kHz [19]. To improve the signal-to-noise ratio, external amplifiers (Ecreon Model E713) were also employed. With a gain of 20, the roll-off frequency is about 140 kHz. The typical sampling period is 5 μs (200 kHz). Each data record possesses 131,072 data points for statistical analysis. The data were divided into 32 blocks. The mean values of each block (4,096 data points) were calculated. Variations of the blocks are estimated to be 0.43 and 0.13 percent for the mean surface pressure coefficient, Cp, and the fluctuating pressure coefficient, Cσp, respectively, which were taken as uncertainty of the experimental data.


## 3. Results and Discussions

### 3.1. Mean and Fluctuating Surface Pressure Distributions

Distributions of the mean surface pressure coefficientCp at M=0.64 are shown in Figure 2, where x*(=x/δ) is the normalized streamwise distance. The origin of the x coordinate is set at the corner. The solid symbol corresponds to convex-corner flows, while the hollow symbol denotes concave-corner flows. It is known that viscous-inviscid interactions in subsonic corner flows affect displacement thickness (or effective local wall surface) near the corner apex [7]. Thus, as can be seen, convex-corner flows accelerate gradually upstream of the corner followed by stronger expansion and then downstream compression. The minimum pressure coefficient is observed near the corner apex. There are stronger upstream expansion and downstream recompression at η = 15°. It is also noted that the level of Cp tends to an equilibrium value at further downstream locations and decreases when η increases. For concave-corner flows at M=0.64, the pressure distributions show a similar shape. The flow decelerates upstream of the corner followed by expansion. The amplitude of Cp appears to increase and decrease linearly within the upstream and downstream influence regions. At further downstream locations, there is a more positive Cp at η = 15°. Moreover, the interaction region tends to expand in both upstream and downstream directions for both test cases, indicating stronger viscous-viscid interactions when η increases. It is also noted that both upstream and downstream pressure gradients are less significant for convex-corner flows. Further, Mach number distributions at M=0.64 are shown in Figure 3. For concave-corner flows, mild variations in streamwise Mach number are observed. The flow expands suddenly and reaches a peak Mach number, Mp, near the convex corner, approaching sonic condition at η = 15°.Figure 2
Distributions of pressure coefficient,M=0.64. Hallow symbol: concave corner. Solid symbol: convex corner.Figure 3
Mach number distributions,M=0.64. Hallow symbol: concave corner. Solid symbol: convex corner.The distributions of normalized surface pressure fluctuations for convex-corner flows atM=0.64 are shown in Figure 4. σp/pw corresponds to the local variation of surface pressure fluctuations. It can be seen that σp/pw increases upstream of the corner and reaches the maximum at immediately downstream location. The rise in σp/pw corresponds to the initial pressure rise of mean surface pressure or downstream adverse pressure gradient, as shown in Figure 2. The value of peak pressure fluctuations is higher when η increases. At further downstream locations, the amplitude of σp/pw approaches an equilibrium level. Note that there are minor variations in σp/pw for subsonic flow over a concave corner.Figure 4
Distributions of surface pressure fluctuations,M=0.64.


### 3.2. Flow Similarity

The presence of a convex corner in a subsonic uniform flow results in expansion near the corner, corresponding to the increment in displacement thickness by viscous-inviscid interaction. The minimum pressure coefficientsCp,min are associated with freestream Mach number and deflection angle. In Figure 5, Cp,min/η2 is plotted against K* for all test cases. It is noted that η is in radian in the following analyses. Although the data are slightly scattered, it can be seen that Cp,min/η2 decreases linearly with K* in convex-corner flows. For concave-corner flows at M=0.64, there is mild adverse pressure gradient upstream of the corner. Cp,max/η2 can also be scaled with K*, as shown in Figure 6, in which Cp,max/η2 increases linearly with K*, but not for the test case of K* = 6.1. Cp,max/η2 appears to be a quadratic function of K* for concave-corner flows. To elaborate further the correlation of flow expansion/compression near the corner with K*, Mp2/η2 is plotted against K* for all test cases. As shown in Figure 7, Mp2/η2 is a quadratic function of K* for both test cases. Mp2/η2 increases significantly when K* increases. However, there are only mild variations in Mp2/η2 with K* for concave-corner flows, demonstrating more significant viscous-inviscid interactions in subsonic convex-corner flow than in subsonic concave-corner flows.Figure 5
Minimum pressure coefficient, convex-corner flows.Figure 6
Maximum pressure coefficient, concave-corner flows.Figure 7
Peak Mach number.Interaction region, including upstream influenceXu and downstream influence Xd, can be employed to highlight viscous-viscid interactions in subsonic convex- and concave-corner flows. Upstream influence, Xu*(=Xu/δ), can be determined as the intercept of the tangent to the maximum pressure gradient with the undisturbed surface pressure (or Cp = 0). Downstream influence, Xd*(=Xd/δ), represents the distance for a disturbed boundary layer back to an equilibrium state and can be estimated from the peak pressure near the corner to the intersection of the tangent through the downstream pressure data with the approximately equilibrium downstream pressure [20]. At M=0.64, variations of Xu* and Xd* with K* are shown in Figure 8. As can be seen, Xu* and Xd* decrease linearly with K*for concave-corner flows. For convex-corner flows, there is a shorter interaction region, corresponding to higher upstream and downstream pressure gradients near the corner. Xu* and Xd* appear to be a quadratic function of K*.Figure 8
Upstream and downstream influences.Pressure fluctuations are coupled with global flow unsteadiness. In general, the shear layer structures in the buffer region are responsible for the generation of high-amplitude wall pressure peaks [21]. Laganelli et al. [22] examined wall pressure fluctuations in the attached boundary layer flow. They noted that σp/pw is proportional to M2. In Figure 9, peak pressure fluctuations, (σp/pw)max/η2, are plotted against K* for convex-corner flows. Note that there is only slight variation in surface pressure fluctuations for concave-corner flows. It can be seen that (σp/pw)max/η2 increases linearly with K*. More intense pressure fluctuations are associated with higher peak Mach number or flow expansion near the convex-corner apex.Figure 9
Peak pressure fluctuations, convex-corner flows.


## 3.1. Mean and Fluctuating Surface Pressure Distributions

Distributions of the mean surface pressure coefficientCp at M=0.64 are shown in Figure 2, where x*(=x/δ) is the normalized streamwise distance. The origin of the x coordinate is set at the corner. The solid symbol corresponds to convex-corner flows, while the hollow symbol denotes concave-corner flows. It is known that viscous-inviscid interactions in subsonic corner flows affect displacement thickness (or effective local wall surface) near the corner apex [7]. Thus, as can be seen, convex-corner flows accelerate gradually upstream of the corner followed by stronger expansion and then downstream compression. The minimum pressure coefficient is observed near the corner apex. There are stronger upstream expansion and downstream recompression at η = 15°. It is also noted that the level of Cp tends to an equilibrium value at further downstream locations and decreases when η increases. For concave-corner flows at M=0.64, the pressure distributions show a similar shape. The flow decelerates upstream of the corner followed by expansion. The amplitude of Cp appears to increase and decrease linearly within the upstream and downstream influence regions. At further downstream locations, there is a more positive Cp at η = 15°. Moreover, the interaction region tends to expand in both upstream and downstream directions for both test cases, indicating stronger viscous-viscid interactions when η increases. It is also noted that both upstream and downstream pressure gradients are less significant for convex-corner flows. Further, Mach number distributions at M=0.64 are shown in Figure 3. For concave-corner flows, mild variations in streamwise Mach number are observed. The flow expands suddenly and reaches a peak Mach number, Mp, near the convex corner, approaching sonic condition at η = 15°.Figure 2
Distributions of pressure coefficient,M=0.64. Hallow symbol: concave corner. Solid symbol: convex corner.Figure 3
Mach number distributions,M=0.64. Hallow symbol: concave corner. Solid symbol: convex corner.The distributions of normalized surface pressure fluctuations for convex-corner flows atM=0.64 are shown in Figure 4. σp/pw corresponds to the local variation of surface pressure fluctuations. It can be seen that σp/pw increases upstream of the corner and reaches the maximum at immediately downstream location. The rise in σp/pw corresponds to the initial pressure rise of mean surface pressure or downstream adverse pressure gradient, as shown in Figure 2. The value of peak pressure fluctuations is higher when η increases. At further downstream locations, the amplitude of σp/pw approaches an equilibrium level. Note that there are minor variations in σp/pw for subsonic flow over a concave corner.Figure 4
Distributions of surface pressure fluctuations,M=0.64.


## 3.2. Flow Similarity

The presence of a convex corner in a subsonic uniform flow results in expansion near the corner, corresponding to the increment in displacement thickness by viscous-inviscid interaction. The minimum pressure coefficientsCp,min are associated with freestream Mach number and deflection angle. In Figure 5, Cp,min/η2 is plotted against K* for all test cases. It is noted that η is in radian in the following analyses. Although the data are slightly scattered, it can be seen that Cp,min/η2 decreases linearly with K* in convex-corner flows. For concave-corner flows at M=0.64, there is mild adverse pressure gradient upstream of the corner. Cp,max/η2 can also be scaled with K*, as shown in Figure 6, in which Cp,max/η2 increases linearly with K*, but not for the test case of K* = 6.1. Cp,max/η2 appears to be a quadratic function of K* for concave-corner flows. To elaborate further the correlation of flow expansion/compression near the corner with K*, Mp2/η2 is plotted against K* for all test cases. As shown in Figure 7, Mp2/η2 is a quadratic function of K* for both test cases. Mp2/η2 increases significantly when K* increases. However, there are only mild variations in Mp2/η2 with K* for concave-corner flows, demonstrating more significant viscous-inviscid interactions in subsonic convex-corner flow than in subsonic concave-corner flows.Figure 5
Minimum pressure coefficient, convex-corner flows.Figure 6
Maximum pressure coefficient, concave-corner flows.Figure 7
Peak Mach number.Interaction region, including upstream influenceXu and downstream influence Xd, can be employed to highlight viscous-viscid interactions in subsonic convex- and concave-corner flows. Upstream influence, Xu*(=Xu/δ), can be determined as the intercept of the tangent to the maximum pressure gradient with the undisturbed surface pressure (or Cp = 0). Downstream influence, Xd*(=Xd/δ), represents the distance for a disturbed boundary layer back to an equilibrium state and can be estimated from the peak pressure near the corner to the intersection of the tangent through the downstream pressure data with the approximately equilibrium downstream pressure [20]. At M=0.64, variations of Xu* and Xd* with K* are shown in Figure 8. As can be seen, Xu* and Xd* decrease linearly with K*for concave-corner flows. For convex-corner flows, there is a shorter interaction region, corresponding to higher upstream and downstream pressure gradients near the corner. Xu* and Xd* appear to be a quadratic function of K*.Figure 8
Upstream and downstream influences.Pressure fluctuations are coupled with global flow unsteadiness. In general, the shear layer structures in the buffer region are responsible for the generation of high-amplitude wall pressure peaks [21]. Laganelli et al. [22] examined wall pressure fluctuations in the attached boundary layer flow. They noted that σp/pw is proportional to M2. In Figure 9, peak pressure fluctuations, (σp/pw)max/η2, are plotted against K* for convex-corner flows. Note that there is only slight variation in surface pressure fluctuations for concave-corner flows. It can be seen that (σp/pw)max/η2 increases linearly with K*. More intense pressure fluctuations are associated with higher peak Mach number or flow expansion near the convex-corner apex.Figure 9
Peak pressure fluctuations, convex-corner flows.


## 4. Conclusions

This paper investigates the flow characteristics of subsonic convex- and concave-corner flows. Mean and fluctuating pressures are presented. A combined flow similarity parameter,K*, based on the small perturbation theory, is employed as a governing parameter to identify the trends in flow expansion/compression, interaction region, and peak pressure fluctuations. For convex-corner flows, variations in minimum pressure coefficient and peak pressure fluctuations can be scaled linearly with K*, but not for upstream and downstream influences. Upstream compression in concave-corner flows appears to be a quadratic function of K*. The results imply that there are small streamline deflections near corner apex, producing proportionally small change in Mach number and pressure. The present results can also be used for a quick estimation of aerodynamic characteristics of a deflected surface in subsonic flows.



---

*Source: 940862-2013-06-27.xml*