3.1. Computational domain and boundary conditions

Flow past a rectangular wing of span length $b$ and with a NACA 0012 section of chord length $c$ is considered. The Reynolds number, based on the chord of the wing, free-stream speed of the incoming flow and kinematic viscosity of fluid, is $1000$. A schematic of the problem set-up and computational domain is shown in figure 1. We simulate only one half of the span of a rectangular wing to reduce the requirement of computational resources. Symmetry flow conditions are imposed at the mid-span on the plane $ADHE$. The semi-aspect ratio of the wing is defined as $sAR = b/2{\bar {c}}$. The upstream ($ABCD$) and downstream ($EFGH$) boundaries are located at a distance of $L_{xu}$ and $L_{xd}$ from the leading edge of the airfoil. Here $L_y$ is the separation between the lateral boundary faces $CDHG$ and $ABFE$. These dimensions for the present study are $L_{xu} = 4.5c$, $L_{xd} = 15.5c$ and $L_y = 10c$. The spanwise extent of the computational domain is $L_z+b/2$. For each $sAR$, the value of $L_z$ is chosen so that the boundary face $BCGF$ is sufficiently far from the tip of the wing. Toppings & Yarusevych (Reference Toppings and Yarusevych2021) and Toppings, Kurelek & Yarusevych (Reference Toppings, Kurelek and Yarusevych2021) used $L_z = 0.5c$ in their experimental investigations for $1.25 \le sAR \le 2.75$ carried out at $Re = 1.25\times 10^{5}$. Here $L_z$ varies with $sAR$ in the present study to ensure that the computations are unaffected by the location of the lateral boundary. It is $2.75c$ for $sAR = 0.25$ while it is $14.5c$ for $sAR = 7.5$, which is the largest span considered in the study. A no-slip condition for velocity is imposed on the surface of the wing. An inflow with a uniform free-stream speed $U$ along the $x$ direction is prescribed on the upstream boundary face $ABCD$. At the downstream boundary face $EFGG$, corresponding to the outflow boundary, the stress vector is assigned a zero value. The component of velocity normal to the plane and the component of stress vector on the lateral boundary faces $ABFE$ and $CDHG$ are prescribed to zero value. To test the assumption of symmetry of the flow about the mid-span of the wing, a few simulations are carried out for a wing with full span. The details of the set-up as well as the comparison of the results from full- and half-span computations are presented in Appendix C. It is found that the flow features and time-averaged aerodynamic coefficients for the full span and half-span, with symmetry conditions imposed at the mid-span, are identical. Therefore all simulations are carried out using half-span. A few computations have also been carried out for an EEW where an airfoil spans the entire lateral extent of the computational domain. The span length of the wing for these computations is $5c$ and symmetry conditions are imposed on the lateral walls ($ADHE$ and $BCGF$).

Figure 1.

Flow past a finite wing: schematic of the computational domain along with the boundary conditions. One half of the wing span is considered and symmetry conditions are imposed at the plane $ADHE$. It lies at mid-span of the wing and is shaded in the sketch for ease of identification.

3.2. Finite element mesh

The finite element mesh for the EEW is formed by stacking, along the span of the wing, several copies of the 2-D mesh around an airfoil. The 2-D mesh consists of a structured region around the airfoil and downstream of it to resolve the boundary layer, its separation and ensuing wake. The height of the first element lying on the surface of the airfoil is $0.005c$. The mesh outside the structured zone is obtained via Delaunay triangulation. The combination of the structured and unstructured mesh is similar to that described in our earlier work for the Eppler airfoil (Pandi & Mittal Reference Pandi and Mittal2019). Such a mesh enables adequate resolution of the flow structures while keeping the number of unknowns to a reasonable level. A mesh convergence study is first carried out for the 2-D mesh. The flow at $Re=1000$ and $\alpha = 14^{\circ }$ is used as a test case. The details are presented in Appendix A.1. It is found that a mesh with $43\,474$ nodes and $86\,602$ triangular three-noded elements is adequate to resolve this flow. A view of the mesh close to the airfoil is shown in figure 2. This mesh consists of $250$ nodes on the surface of the airfoil. The lift and drag coefficients obtained from present computations with this mesh and those from the literature (Mittal & Tezduyar Reference Mittal and Tezduyar1994; Liu et al. Reference Liu, Li, Zhang, Wang and Liu2012; Kurtulus Reference Kurtulus2015; Meena, Taira & Asai Reference Meena, Taira and Asai2017; Di Ilio et al. Reference Di Ilio, Chiappini, Ubertini, Bella and Succi2018) for various $\alpha$ are presented in Appendix A.1. They are in very good agreement.

Figure 2.

Flow past a finite wing: close-up view of the 2-D mesh $M_{2D}^2$ (details are listed in table 1 of Appendix A.1) for NACA 0012 airfoil in the $xy$ plane at $\alpha = 14^{\circ }$.

The mesh for the EEW is generated by stacking $N_z=160$ copies of the 2-D mesh for the airfoil. A convergence study, by utilizing an additional mesh with twice the number of elements along the span, is presented in Appendix A.2.1. The mesh for the finite wing is formed by stacking copies of two different 2-D meshes. As for the EEW, the mesh on the finite wing is formed by stacking copies of the 2-D mesh for the airfoil. An additional 2-D mesh is generated from the mesh around the airfoil by introducing nodes inside the region occupied by the airfoil. The number of nodes and elements of this mesh are $45\,760$ and $91\,421$, respectively. Copies of this 2-D mesh are stacked beyond the spanwise extent of the finite wing. Unlike in the mesh for EEW, the copies of the 2D mesh are not uniformly spaced along the span, but rather stacked in a manner such that they form a finer mesh towards the wing tip to adequately resolve the wing-tip vortex and the boundary layer at the face of the wing tip. The size of the element at the wing tip is $0.0035c$. The convergence study for the finite wing with $sAR = 5$ for $Re = 1000$ and $\alpha = 14^{\circ }$ is presented in Appendix A.2.2. Three meshes of varying spanwise resolution are utilized. It is found that the mesh with $186$ elements along the span of the wing provides adequate spatial resolution for the range of parameters in this work. Meshes for wings with other $sAR$ are generated to keep the same resolution along the span. We extend the convergence study for the finite wing with $sAR = 1.25$ for $Re = 1000$ and $\alpha = 5^{\circ }$ by utilizing two meshes (see Appendix A.2.3). Both meshes result in almost the same result and it is in excellent agreement with that reported by Edstrand et al. (Reference Edstrand, Sun, Schmid, Taira and Cattafesta2018). The details are presented in Appendix A.2.3.

3.3. Identification of vortex structures

One of the objectives of the present study is to identify the various flow structures and their interactions in the presence of a wing-tip vortex. Epps (Reference Epps2017) and Zhang et al. (Reference Zhang, Liu, Xian and Du2018) presented a brief review of several vortex identification methods. Some of the widely used methods are based on the $Q$ (Hunt, Wray & Moin Reference Hunt, Wray and Moin1988), $\lambda _2$ (Jeong & Hussain Reference Jeong and Hussain1995) and $\varOmega$ criterion (Liu et al. Reference Liu, Wang, Yang and Duan2016). These are described in Appendix B. Our earlier studies (Pandi & Mittal Reference Pandi and Mittal2019; Mittal et al. Reference Mittal, Pandi and Hore2021) utilized $Q$ criterion to visualize the vortex structures. Mittal et al. (Reference Mittal, Pandi and Hore2021) showed that all three methods ($Q$, $\lambda _2$ and $\varOmega$) give identical vortical structures. A similar analysis is carried out in this study, using $Q$ and $\lambda _2$ criterion, for various flows and is presented in Appendix B. It is found that the $Q$ and $\lambda _2$ methods reveal identical flow structures. The $Q$ criterion is used for identification of vortices in this work.

4.1. Flow structures

The flow past an EEW is investigated at various angles of attack. The $Q (=0.1)$ isosurface for the fully developed unsteady flow at a time instant corresponding to the peak value of lift coefficient is shown in figure 3 for various angles of attack. The vortices are coloured with the spanwise component of vorticity ($\omega _z = \pm 2$). The flow stays steady for $\alpha \le 7^{\circ }$ and becomes unsteady thereafter via primary instability of the wake leading to vortex shedding. A similar observation was made by Kurtulus (Reference Kurtulus2015) and Di Ilio et al. (Reference Di Ilio, Chiappini, Ubertini, Bella and Succi2018) via a 2-D calculation. The shed vortices are parallel to the axis of the wing. They develop spanwise undulations causing the wake to transition from a two- to three-dimensional state via a mode C instability (Zhang et al. Reference Zhang, Fey, Noack, Konig and Eckelmann1995; Yildirim, Rindt & van Steenhoven Reference Yildirim, Rindt and van Steenhoven2013) for $\alpha > 12^{\circ }$. Pandi & Mittal (Reference Pandi and Mittal2019); Deng et al. (Reference Deng, Sun and Shao2017) reported a similar transition, but at different $Re$ and $\alpha$, for the Eppler 61 and NACA 0015 sections. As observed by Pandi & Mittal (Reference Pandi and Mittal2019) for an Eppler 61 airfoil, each wave of the mode C instability, along the span, evolves to a hairpin vortex (see the close-up view in figure 3c). The two limbs of each hairpin vortex are associated with a streamwise vorticity of opposite sign. Pandi & Mittal (Reference Pandi and Mittal2019) presented detailed features of the mode C instability in terms of the $RT$ symmetry ($R$ denotes reflection about the wake axis and $T$ denotes translation in time) and time periodicity of the three-dimensional flow. The flow in the present study has the same characteristics. It exhibits neither odd-$RT$ nor even-$RT$ symmetry and has a time period of $2T$, where $T$ corresponds to the time period of primary vortex shedding. The flow at $\alpha = 14^{\circ }$ is regular and periodic along the span. The spanwise wavelength ($\lambda _z$) estimated from the streamwise flow structure is $0.33c$. Hoarau et al. (Reference Hoarau, Braza, Ventikos, Faghani and Tzabiras2003) proposed the length scale for an airfoil, placed at an angle of attack $\alpha$, to be $d=c \sin \alpha$. The spanwise wavelength of the flow structures, based on this scaling, is $\lambda _z \sim 1.36d$ and is close to the value reported by Pandi & Mittal (Reference Pandi and Mittal2019) ($1.12d \le \lambda _z \le 1.45d$), Deng et al. (Reference Deng, Sun and Shao2017) ($1.35d \le \lambda _z \le 1.62d$) and He et al. (Reference He, Gioria, Perez and Theofilis2017) ($1.63d \le \lambda _z \le 1.95d$) for other airfoils. Pandi & Mittal (Reference Pandi and Mittal2019) presented characteristics of various modes of instabilities that cause transition to three dimensionality in flow past airfoils and a circular cylinder, from past studies. It was shown that the onset Reynolds number and wavelength of the instability for the circular cylinder and airfoils are comparable if one uses a diameter as the length scale for the cylinder and $d = c\sin \alpha$ for the airfoil.

Figure 3.

Flow past an EEW at $Re = 1000$: $Q (= 0.1)$ isosurface for an instantaneous flow coloured with the spanwise component of vorticity ($\omega _z = \pm 2$) for $\alpha =$ (a) $8^{\circ }$, (b) $12^{\circ }$ and (c) $14^{\circ }$. Also shown in (c) is a close-up view of the hairpin vortices. The spanwise extent of the domain is $5c$.

The difference in pressure between the upper and lower surfaces of the wing results in modification of the flow near the wing tip of a finite wing, causing the formation of a wing-tip vortex. This streamwise vortex interacts with the vortex shedding to further modify the flow. We compare the flow past a finite wing of $sAR = 5$ with that past an EEW to understand the differences and their effect on aerodynamic coefficients at various $\alpha$. The $Q (=0.1)$ isosurface of instantaneous flow coloured with the spanwise component of vorticity ($\omega _z = \pm 2$) for an $sAR = 5$ wing at various angles of attack is shown in figure 4. The wing-tip vortex suppresses the vortex shedding near the tip of the wing. In addition, it weakens the vortex shedding over the bulk of the span compared with that for the EEW. As a result, at $\alpha = 8^{\circ }$, the vortex shedding for the finite wing is quite weak compared with that of the EEW (see figures 3a and 4a).

Figure 4.

Flow past the $sAR = 5$ wing at $Re = 1000$: $Q (=0.1)$ isosurface for an instantaneous flow coloured with the spanwise component of vorticity ($\omega _z = \pm 2$) for $\alpha =$ (a) $8^{\circ }$, (b) $10^{\circ }$, (c) $12^{\circ }$ and (d) $14^{\circ }$. The wing-tip vortex, central and end cells, dislocations, vortex splitting, vortex linkages and arch vortices are identified in the images. The inset in (d) shows the flow mirrored about the mid-span of the wing to visualize the arch vortices.

According to Helmholtz's theorem (Batchelor Reference Batchelor1967), a vortex must extend to the boundaries of the fluid or form a closed path. The spanwise vortices of an EEW extend to the lateral boundaries of the computational domain (see figure 3). However, owing to the wing-tip vortex, linkages form between vortices of opposite polarity near the tip of a finite wing. They are clearly seen in figure 4(b–d) for $\alpha \ge 10^{\circ }$. Figure 5 shows the linkages via the $Q$ isosurface as well as a few vortex lines that pass through the core of an adjacent pair of counter-rotating vortices and connect to form a closed loop. We note that the linkages are not clear for $\alpha =8^{\circ }$ from the $Q$ isosurface (figure 4a). However, vortex lines in figure 5(a) clearly show the linkages even at this $\alpha$. Similar linkages have been reported for a finite wing (Zhang et al. Reference Zhang, Hayostek, Amitay, He, Theofilis and Taira2020) and flow past a circular cylinder near the side wall (Mittal et al. Reference Mittal, Pandi and Hore2021). Another flow feature observed at $\alpha = 14^{\circ }$ is the arch vortices in the mid-span region and at streamwise locations far downstream of the wing. These appear more clearly in images that show the flow for the full span of the wing. The inset in figure 4(d) shows a close-up view of the arch vortices for the full span. This image is constructed by mirroring the flow about the mid-span of the wing. A study to assess the effect of the symmetry boundary condition at the mid-span to constrain the flow is presented in Appendix C. The flow for the full span of the wing is computed for certain cases and compared with those for the half-span. It reveals that computations with half-span and full span lead to the same results. Arch vortices for $sAR = 1$ and $4$ are shown in the appendix. They are very similar to those reported in earlier studies for flow past a non-stationary wing (Visbal, Yilmaz & Rockwell Reference Visbal, Yilmaz and Rockwell2013; Rockwood et al. Reference Rockwood, Medina, Garmann and Visbal2019; Visbal & Garmann Reference Visbal and Garmann2019).

Figure 5.

Flow past the $sAR = 5$ wing at $Re = 1000$: $Q (=0.1)$ isosurface for an instantaneous flow coloured with the spanwise component of vorticity for $\alpha =$ (a) $8^{\circ }$ and (b) $10^{\circ }$. A few vortex lines that pass through the core of the adjacent pair of counter-rotating vortices connect to form a closed loop.

An interesting aspect of vortex shedding on a wing of finite span is the formation of cells. This has been studied in detail for bluff body flows (Gerrard Reference Gerrard1978; Eisenlohr & Eckelmann Reference Eisenlohr and Eckelmann1989; Williamson Reference Williamson1989; Behara & Mittal Reference Behara and Mittal2010; Mittal et al. Reference Mittal, Pandi and Hore2021). The variation of frequency of vortex shedding along the span is used to identify the number of cellular structures in the wake (Williamson Reference Williamson1989; Behara & Mittal Reference Behara and Mittal2010; Mittal et al. Reference Mittal, Pandi and Hore2021). The frequency of vortex shedding is constant within a cell. Similar to the flow of an EEW, the axis of the primary vortices shed from a wing of $sAR = 5$ and $\alpha = 8^{\circ }$ is parallel to the span of the wing and the frequency of vortex shedding is the same at all spanwise locations, resulting in a single cell across the span. With an increase in angle of attack, the primary vortices are no longer parallel to the span of the wing. Furthermore, the vortex shedding occurs in two cells for $\alpha \ge 10^{\circ }$. We refer to the cell near the mid-span as the central cell and the one towards the wing tip as the end cell. The frequency of vortex shedding changes across the boundary of adjacent cells. A larger number of vortices, per unit time, are shed along one part of the span compared with other. Periodically in time, the extra vortex in the cell with a larger number of vortices connects with the neighbouring vortices in the adjacent cell leading to the formation of ‘knots’ (Eisenlohr & Eckelmann Reference Eisenlohr and Eckelmann1989) or ‘dislocations’ (Williamson Reference Williamson1989; Zhang et al. Reference Zhang, Hayostek, Amitay, He, Theofilis and Taira2020; Mittal et al. Reference Mittal, Pandi and Hore2021).

Dislocations can be of several kinds (Mittal et al. Reference Mittal, Pandi and Hore2021). For example, a pair of vortices of the same polarity in one cell may connect with a single vortex of the same polarity in an adjoining cell to form a fork-type ($D_f$) dislocation. A modification to the fork-type dislocation was reported by Mittal et al. (Reference Mittal, Pandi and Hore2021) wherein vortices of opposite polarity form an additional linkage along with those in a fork-type dislocation. This is referred to as the connected fork-type dislocation. Another variant is the mixed-type dislocation that is of a fork type at one instant and a connected fork type at another instant (Mittal et al. Reference Mittal, Pandi and Hore2021). Vortex dislocations are marked in figure 4(b–d). Mixed-type dislocations are observed at $\alpha = 10^{\circ }$ while they are of fork type at $\alpha = 14^{\circ }$. The mixed-type dislocations have been observed earlier in bluff body flows. However, they are being reported for a finite wing for the first time. An interesting feature of the dislocations at $\alpha = 12^{\circ }$ is that the vortices split and reconnect in the far wake forming fork-type and reverse fork-type dislocations. As a result, a single cell in the near wake degenerates to two cells in the far wake. We note that this phenomenon has not been reported earlier in any flow.

4.2. Review of lifting line theory

We briefly review the lifting line theory for a finite wing and propose a small modification to apply it to relatively low $Re$ flows. Consider a rectangular wing of span $b$ with no twist and the same airfoil section of chord length $c$ all along the span (see figure 1), placed in a uniform flow at an angle of attack, $\alpha$. We refer to $\alpha$ as the geometric angle of attack. The airfoil is assumed to be symmetric. The lifting characteristic of the wing is modelled by an infinite number of horseshoe vortices. The bound vortex of each horseshoe vortex lies along the quarter-chord line of the wing forming a ‘lifting line.’ The downwash induced by the vortices, estimated using the Biot–Savart law, results in tilting of the local lift vector at each spanwise station by $\alpha _i(z)$, referred to as the induced angle of attack. The effective angle of attack experienced by the local airfoil at each spanwise section is therefore $\alpha _{e}(z)=\alpha -\alpha _i(z)$. Assuming that the flow along the span is not significant, the local lift coefficient at each span location is estimated from the Kutta–Joukowski theorem as $C_l(z) = 2 \varGamma (z) / U c$, where $\varGamma (z)$ is the circulation at each spanwise section (Bertin Reference Bertin2002; Anderson Reference Anderson2017). Here $C_l(z)$ is also related to the effective angle of attack as $C_l(z) = \alpha _{e}(z)a_o$, where $a_o$ is the lift curve slope of the airfoil section. Equating the local lift coefficient obtained from these two approaches results in the following equation for circulation:

(4.1)\begin{equation} \alpha = \frac{{2 \varGamma(z)} }{{a_o U c}} + \frac{1 }{{4 {\rm \pi}U}} \int_{{-}b/2}^{b/2} \frac{{{\rm d}\varGamma}/{{\rm d}z} }{{z_o-z}} \,{\rm d}z. \end{equation}

Using a transformation $z = -b/2\cos \theta$, $0 \le \theta \le {\rm \pi}$, the spanwise variation of circulation is approximated using a Fourier sine series with $N$ coefficients: $\varGamma (\theta ) = 2 b U \sum _{1}^{N} {A_n~{\rm sin}~n\theta }$. Substitution of this approximation of $\varGamma$ in (4.1) and rearranging the terms results in the following equation that can be utilized to estimate the Fourier coefficients for $\varGamma (\theta )$:

(4.2)\begin{equation} \alpha = \frac{{2 b} }{{a_oc}} \sum _{1}^{N} {A_n\sin n\theta} + \sum_{1}^{N} nA_n\frac{\sin n\theta}{\sin\theta}. \end{equation}

A collocation method is used to generate $N$ equations at the same number of discrete spanwise locations. The collocation points are distributed all along the span. However, their density increases towards the wing tip. The resulting Fourier coefficients are utilized to estimate the spanwise distribution of the coefficient of lift ($C_l(\theta ) = 4 b U c \sum _{1}^{N} {A_n\sin n\theta }$), the total lift coefficient for the finite wing ($C_L = A_1 {\rm \pi}b c$) and the total induced drag coefficient (Bertin Reference Bertin2002; Anderson Reference Anderson2017):

(4.3)\begin{equation} C_{Di} = \frac{{C_L}^2 }{{{\rm \pi} AR}} \left(1 + \sum_{2}^{N} n \left(\frac{A_n}{A_1}\right)^2 \right). \end{equation}

In the classical lifting line theory, the slope of the variation of the local lift coefficient with angle of attack is assumed to be $2{\rm \pi} \,{\rm rad}^{-1}$ for a thin airfoil (Bertin Reference Bertin2002; Anderson Reference Anderson2017). Unlike at large $Re$, the lift curve slope ($a_o$) for the airfoil section at $Re = 1000$ is much lower. We estimate $a_o$ from direct numerical simulations at various $\alpha$ for an EEW, and utilize this estimated value in (4.1) and (4.2) to obtain the coefficients of a finite wing.

4.3. Force coefficients

Figure 6 shows the variation of time-averaged aerodynamic force coefficients and aerodynamic efficiency, with $\alpha$, for an EEW and $sAR = 5$ wing at $Re = 1000$. As expected, the lift coefficient of a finite wing is lower than that for the EEW at each $\alpha$. It is inline with lifting line theory (Prandtl Reference Prandtl1921) wherein the finite wing, unlike an EEW, experiences a reduction in the effective angle of attack at each spanwise section. The prediction from the lifting line theory, with a small modification to account for the low $Re$ flow, is also shown in figure 6(a). Compared with flows at high $Re$ (Winkelmann et al. Reference Winkelmann, Barlow, Saini, Anderson and Jones1980; Sheldahl & Klimas Reference Sheldahl and Klimas1981), the lift curve slope at $Re = 1000$ is much lower. Consistent with the observations of Taira & Colonius (Reference Taira and Colonius2009) and Zhang et al. (Reference Zhang, Hayostek, Amitay, He, Theofilis and Taira2020), unlike at high $Re$, as indicated by the $C_L-\alpha$ curve, the flow does not exhibit stall at low $Re$. The drag coefficient shown in figure 6(b) from the lifting line theory is the sum of $\overline {C_D}$ for the EEW and the estimate of the induced drag coefficient for the finite wing. The induced drag is due to the tilting of the local lift vector. Also shown in the figure is the time-averaged drag coefficient from the computations for the viscous flow for the EEW as well as the $sAR=5$ finite wing. Contrary to the prediction from the lifting line theory, the present results show that the drag coefficient for the finite wing is lower than that for the EEW. Computations have also been carried out to obtain the steady flow by leaving out the unsteady terms in the governing equations. It is found that the drag coefficient for the EEW and $sAR=5$ wing are quite comparable; the values for the EEW are marginally higher. Furthermore, they are both lower than those for the time-averaged flow. The higher difference in the drag coefficient for the unsteady flow, between the finite wing and EEW, is primarily because of reduced vortex shedding activity in the finite wing owing to the wing-tip vortex. We investigate this further in a later section in the paper. The loss in lift for the finite wing outweighs the decrease in drag leading to its overall decrease in aerodynamic efficiency compared with the EEW, as shown in figure 6(c).

Figure 6.

Flow past an EEW and $sAR = 5$ wing at $Re = 1000$: variation of the time averaged (a) $\overline {C_L}$, (b) $\overline {C_D}$ and (c) $\overline {C_L}/\overline {C_D}$ with angle of attack. The data obtained using lifting line theory (LLT) (Prandtl Reference Prandtl1921) are plotted in each image. The best linear fit for the lift coefficient of an EEW is shown in (a) via the dashed green line.

Figure 7(a) shows the spanwise variation of the local lift coefficient for the time-averaged flow for different $\alpha$. Also shown in the figure is the distribution from the lifting line theory (Prandtl Reference Prandtl1921). As per the lifting line theory, the local lift coefficient on a wing with a rectangular planform is maximum at the wing root and decreases monotonically towards the wing tip. The lift distribution from the computations at $Re=1000$ is in good agreement with that from the lifting line theory for $\alpha \le 8^{\circ }$. However, a non-monotonic variation is revealed for $\alpha \ge 10^{\circ }$. A local maximum is seen near the wing-tip region for $\alpha \ge 12^{\circ }$. A similar observation was made by Zhang et al. (Reference Zhang, Hayostek, Amitay, He, Theofilis and Taira2020) for flow at large $\alpha$ at $Re = 400$. We attribute the non-monotonic distribution of $\overline {C_l}$ to streamwise vortices that alter the distribution of circulation along the span. This is explained later in § 7 via a discrete vortex model. The spanwise variation of $\overline {C_d}$ is shown in figure 7(b). Also plotted is the prediction from the lifting line theory. It is obtained by adding the induced drag coefficient at each section from the tilting of the local lift vector to the drag coefficient for the EEW. The figure brings out the limitation of the lifting line theory in prediction of drag at low $Re$.

Figure 7.

Flow past the $sAR = 5$ wing at $Re = 1000$: spanwise variation of the time-averaged sectional force coefficient (a) $C_l$ and (b) $C_d$ at various $\alpha$. The sectional force coefficient distribution from the simulation is compared with theoretical data obtained using lifting line theory (LLT).

5.1. Flow structures

We investigate the effect of $sAR$ at two representative angles of attack: $\alpha = 8^{\circ }$ and $14^{\circ }$. Vortex shedding is quite weak at $\alpha =8^{\circ }$ for the EEW. The wing-tip vortex on a finite wing weakens it further by reducing the suction on the upper surface of the wing, thereby reducing the adverse pressure gradient. As a result, the flow for $sAR \le 4$ at $\alpha = 8^{\circ }$ is devoid of vortex shedding. Vortex shedding is quite strong at $\alpha =14^{\circ }$ and its interaction with the wing-tip vortex results in interesting flow features. To further understand the interplay between the vortex shedding and wing-tip vortex, we explore the strength and structure of the wing-tip vortex later in the paper. The $Q$ isosurface of instantaneous flow for various $sAR$ is presented in figure 8. Vortex shedding is suppressed in the region close to the wing tip; the flow achieves a steady state for $sAR \le 0.25$ (figure 8a) while shedding is observed away from the wing tip for larger $sAR$. Vortex shedding for $sAR=0.5$ occurs via hairpin vortex structures (figure 8b). These hairpin vortices are different than those that appear in the wake of the EEW as a result of the mode C instability (figure 3c). The hairpin vortices resulting from the mode C instability, also reported by Pandi & Mittal (Reference Pandi and Mittal2019), have a significantly smaller spanwise wavelength compared with those for an $sAR=0.5$ wing. Similar hairpins are seen for a finite wing, away from the wing tip, albeit for high $sAR$ (see inset of figure 8g). The interaction of the hairpin vortices with the wing-tip vortex for $sAR=0.5$ as it convects downstream results in undulations and the formation of linkages. The structure of the linkage is highlighted in the inset in figure 8(b) that shows a close-up view of the flow mirrored about the mid-span of the wing. The flow for $sAR = 0.5$ was also computed for the full span of the wing. It is found that the flow retains symmetry about the mid-span and is identical to the flow for the half-span. Vortex linkages for $sAR \ge 2$ have a very different structure compared with that for $sAR=0.5$. As shown in figures 5 and 8(d–g), they form between spanwise vortices of opposite polarity and do not link them with the wing-tip vortices. The linkages for $sAR=1$ have attributes that are common to a low and high $sAR$ wing. They are shown in the close-up view in figure 8(c) and marked via broken lines and form in the near wake to link spanwise vortices of opposite polarity. As in the case of $sAR = 0.5$, they form an additional linkage with the wing-tip vortex as they convect downstream. Arch vortices form in the far wake near the plane of symmetry for $sAR \ge 1$. Unlike a hairpin vortex whose limbs are aligned streamwise (Zhang et al. Reference Zhang, Hayostek, Amitay, He, Theofilis and Taira2020), the limbs of an arch vortex are in the cross-flow direction (Visbal et al. Reference Visbal, Yilmaz and Rockwell2013; Rockwood et al. Reference Rockwood, Medina, Garmann and Visbal2019; Visbal & Garmann Reference Visbal and Garmann2019). The wake consists of one cell for $sAR \le 2$. Two cells are observed for $sAR \ge 3$ with a fork-type dislocation (described later in the paper) at their boundary. The vortex dislocation forms at the same spanwise location, at different time instants, for $sAR = 3$. However, it is generated at different spanwise locations for $sAR \ge 4$.

Figure 8.

Flow past a finite wing at $\alpha =14^{\circ }$ and $Re = 1000$: $Q (=0.1)$ isosurface for an instantaneous flow coloured with the spanwise component of vorticity ($\omega _z = \pm 2$) for $sAR =$ (a) $0.25$, (b) $0.5$, (c) $1$, (d) $2$, (e) $3$, ( f) $4$ and (g) $7.5$. Various flow features such as central and end cells, dislocations, hairpin vortices, vortex linkages and arch vortices are identified in the images. The insets in (b) that show hairpin vortex and linkage for $sAR=0.5$ are with the flow mirrored about the mid-span of the wing.

5.2. Force coefficients and comparison with lifting line theory

The variation of the time-averaged lift and drag coefficient with $sAR$ is plotted in figures 9(a) and 9(b), respectively. For reference, the coefficients for the EEW along with the estimates from lifting line theory for the finite wing are also shown. We note that the lifting line theory cannot be expected to yield good predictions for low sAR wings owing to the intrinsic three dimensionality of the flow. The lift coefficient increases with an increase in $sAR$ and approaches the value for an EEW. The drag coefficient, predicted by the viscous computations, shows a very interesting trend with an increase in $sAR$. As per the lifting line theory, it decreases with an increase in $sAR$ due to a decrease in induced drag. However, the viscous computations for the finite wing reveal a non-monotonic variation of the time-averaged drag coefficient with an increase in $sAR$ (see figure 9b). It decreases with an increase in $sAR$ for $sAR \le 0.5$ and increases thereafter. Furthermore, for all $sAR$, it is smaller than that for the EEW. To explore this curious behaviour, we study the pressure and viscous contributions. It is found that the lift is primarily due to the pressure distribution and has very little contribution from viscous force (not shown here). The contributions from both the components is significant for drag and are shown in figure 9(c). Viscous force is the dominant contributor to drag for low $sAR (=0.25)$. Its contribution decreases while the pressure contribution increases with an increase in $sAR$. This explains the non-monotonic variation of $\overline {C_D}$ with $sAR$. The steady increase in the pressure contribution with $sAR$ causes the net drag coefficient to increase for $sAR > 0.5$, unlike the prediction from lifting line theory. We note that the wing-tip vortex reduces the suction on the upper surface near the wing tip. This leads to a reduced adverse pressure gradient and, thereby, reduced flow separation compared with an EEW. The fraction of wing span affected by the wing-tip vortex reduces with an increase in $sAR$, thereby causing a reduction in its effect. We explore the effect of a wing-tip vortex on vortex shedding. Figure 10(a) shows the spanwise variation of $\overline {v'v'}$ at ($x/c=2.74$, $y/c=0.051$) for various $sAR$ along with that for an EEW. The Reynolds stress is significantly smaller near the wing tip for all $sAR$, indicating suppression of vortex shedding at those locations due to the wing-tip vortex. The net unsteadiness across the span, for each $sAR$, is estimated via spanwise integration of the Reynolds stress as ${\widetilde {v'v'}}=({\int _{0}^{b/2} \overline {v'v'}(z) \,{\rm d}z})/({\int _{0}^{b/2} \,{\rm d}z})$. Its variation with $sAR$ is shown in figure 10(b). We note that $\widetilde {v'v'}$ increases with an increase in $sAR$ and approaches the value for an EEW, signifying the increasing strength of vortex shedding with an increase in $sAR$. The trend is similar to the variation of drag coefficient with $sAR$ (figure 9b).

Figure 9.

Flow past a finite wing at $\alpha =14^{\circ }$ and $Re = 1000$: variation of the time averaged (a) $C_L$ and (b) $C_D$ with $sAR$ along with the prediction from lifting line theory for $C_L$ and $C_D$. Also shown in (c) is the variation of time averaged pressure ($C_{Dp}$) and the viscous component ($C_{Dv}$) of drag coefficient with $sAR$. The sectional $\overline{C_l}$ for steady and time-averaged flows is compared in (d).

Figure 10.

Flow past a finite wing at $\alpha = 14^{\circ }$ and $Re = 1000$: spanwise variation of (a) $\overline {v'v'}$ sampled at ($x/c=2.74$, $y/c=0.051$) for various $sAR$. The net unsteadiness in the wake ($\widetilde {v'v'}$), estimated via spanwise integration of $\overline {v'v'}$, is plotted in (b). Also shown in (a,b) is the data for the EEW.

The spanwise variation of sectional lift coefficient for the time-averaged flow for $\alpha =14^{\circ }$ is shown in figure 9(d) for several $sAR$. As is observed for $sAR=5$ (see figure 7a), peak $\overline {C_l}$ is not at the root but close to wing tip for all $sAR$. The non-monotonic spanwise variation in $\overline {C_l}$ is due to the streamwise vortices that are formed near the wing-tip region. Images of the flow demonstrating the same are shown later in the paper. The peak $\overline {C_l}$ increases with an increase in $sAR$ for $sAR > 0.5$; $\overline {C_l}$ at the wing root also increases with an increase in $sAR$. Figure 9(d) shows the $\overline {C_l}$ distribution for the steady flow, computed by dropping the unsteady terms in the governing equations. We note that, compared with the steady flow, $\overline {C_l}$ is higher for the unsteady flow indicating that vortex shedding enhances $\overline {C_l}$. The higher unsteadiness for the EEW compared with the finite wing is also the reason for larger $\overline {C_D}$ for the former, as shown in figure 6(b). A further observation is that $\overline {C_l}$ for the steady flow at the wing root is quite similar for all the $sAR$ shown in the figure. Therefore, a larger increase of $\overline {C_l}$ for higher $sAR$ for the unsteady flow is indicative of a stronger vortex shedding. This is consistent with the observation from figure 8. Contrary to the prediction from lifting line theory, the peak $\overline {C_l}$ for the steady flow is also not at the wing root but at a location close to that for the unsteady flow. We will address this in a later section.

6.1. Frequency of vortex shedding

The number of cells in the wake can be identified from the frequency of vortex shedding measured at probes located at various locations along the span (Williamson Reference Williamson1989; Behara & Mittal Reference Behara and Mittal2010; Mittal et al. Reference Mittal, Pandi and Hore2021). The frequency is constant within a cell. The varying frequencies from multiple cells leave an imprint in the time histories of aerodynamic forces integrated along the span. The time histories of lift coefficient for various $sAR$ along with their frequency spectra are presented in figure 11. For reference, the plots for the EEW are shown as well. The frequency is non-dimensionalized with the free-stream speed and a length scale corresponding to $d =c\sin \alpha$ (Hoarau et al. Reference Hoarau, Braza, Ventikos, Faghani and Tzabiras2003). The flow is steady for $sAR = 0.25$ and the associated lift coefficient is invariant with time (not shown in the figure). Single cell shedding for $sAR \le 2$ and two cell shedding is observed for $sAR \ge 3$ (see figure 8). The difference in the frequencies of central and end cells lead to low frequency modulation in the time histories (figure 11a). The Strouhal number for the two frequencies are referred to as $St_C$ and $St_E$, respectively, and are marked in figure 11(b). The frequency of vortex shedding has also been estimated from the time variation of the cross-flow component of velocity at various span locations in the near wake at ($x/c=1.25, y/c=-0.068$). The frequencies so measured are identical to those marked in figure 11(b). The flow is explored in more detail later in the paper. The fluctuation in lift coefficient is largest for the EEW and decreases with a decrease in $sAR$.

Figure 11.

Flow past a finite wing at $\alpha =14^{\circ }$ and $Re = 1000$: (a) time histories of $C_L$ for the EEW and finite wing of various $sAR$, and (b) their frequency spectra. The Strouhal number corresponding to the dominant frequency for each spectrum is marked in the figure.

The variation of $St_C$ and $St_E$ with $sAR$ is shown in figure 12(a). Here $St_C$ increases with an increase in $sAR$ and approaches the $St$ for the EEW. The vortex shedding frequency in the end cell is smaller than in the central cell. The difference in $St_C$ and $St_E$ is related to the beat frequency in the time signals of force coefficients and the probes lying at the boundary of the cells (Williamson Reference Williamson1989; Mittal et al. Reference Mittal, Pandi and Hore2021). It increases with an increase in $sAR$ up to $sAR=5$ beyond which it saturates to a constant value. Both $St_C$ and $St_E$ reach close to their respective saturated values at $sAR=5$. The variation of $St_C$ and $St_E$ with $\alpha$ is explored for the $sAR = 5$ wing; see figure 12(b). Also plotted is the $St$ for the EEW. Vortex shedding is very weak at $\alpha =8^{\circ }$ for the finite wing; the flow is close to steady near the wing tip (see figure 4a). Nevertheless, the Strouhal number is close to that for the EEW at this $\alpha$. The vortex shedding at $\alpha =10^{\circ }$ is relatively stronger and occurs in two cells (see figure 4b). The frequency of vortex shedding in the central cell, near the root, is close to that for the EEW; it is smaller in the end cell. A very interesting phenomenon is observed at $\alpha =12^{\circ }$. The near wake exhibits shedding in a single cell (figure 4c). However, vortex splitting downstream in the wake leads to two cells of vortex shedding in the far wake. Therefore, probes located in the near wake ($x/c = 1.25, y/c = -0.068$) show a single cell with the same frequency of shedding along the span ($St_C=0.159$), but indicate two cell shedding in the far wake ($x/c = 11.62$, $y/c = -0.026$) ($St_E=0.159, St_C=0.145$). We note that the $St$ shown in figure 12(b) for $\alpha = 12^{\circ }$ is in the near wake where the shedding is in a single cell. The central cell experiences a decrease in shedding frequency in the downstream wake. Unlike at other $\alpha$, $St_C$ is smaller than $St_E$. The frequency of occurrence of dislocation, the difference in $St_E$ and $St_C$, increases with an increase in $\alpha$.

Figure 12.

Flow past a finite wing at $Re = 1000$: variation of non-dimensional vortex shedding frequency in the central cell ($St_C$) and end cell ($St_E$) with (a) $sAR$ for $\alpha = 14^{\circ }$ and (b) $\alpha$ for $sAR = 5$. The $St$ for EEW is also marked. The frequency is estimated from probes placed along the span in the near wake at ($x/c = 1.25$, $y/c = -0.068$).

The dislocations in two cell vortex shedding are of a varied kind depending on the combination of $sAR$ and $\alpha$. They are of mixed type ($D_{f-fc}$) (Mittal et al. Reference Mittal, Pandi and Hore2021), fork type ($D_f$) (Williamson Reference Williamson1989; Mittal et al. Reference Mittal, Pandi and Hore2021) and a mix of fork ($D_f$) and reverse fork type ($D_{rf}$) reported for the first time in the present study. Figure 13 shows a summary of the various states of flow in the $sAR\unicode{x2013}\alpha$ plane. The classification is based on the type of dislocations and number of cells in the flow. The parameters for which the flow achieves a steady state is also shown. The boundaries of the various regimes are indicative extents where a similar flow structure is expected. The hatched region in the figure shows the regime for which the dislocation is generated at various spanwise locations. These features are described in detail in the following sections.

Figure 13.

Flow past a finite wing: classification of the flow in the $sAR\unicode{x2013}\alpha$ plane on the basis of the number of cells, type of dislocation and its formation along the span. The boundaries of the various regions are not exact, but indicative in nature and are based on the cases for which computations have been carried out and marked by hollow symbols. The cases for which a full-span simulation is initiated from the solution from half-span, by reflecting it about the mid-span, are highlighted using solid symbols. The dislocations, for two cell shedding, are generated at the same spanwise location except in the hatched zone where they are formed at different spanwise locations at different time instants. Here $D_f$, $D_{fc}$ and $D_{rf}$ indicate fork-type, connected fork-type and reverse fork-type dislocations.

6.2. Mixed-type dislocation ($D_{f-fc}$) generated at the same spanwise location: $sAR = 5$, $\alpha = 10^{\circ }$

Dislocations occur at the boundary of cells with a different vortex shedding frequency. We present an example of flow wherein the dislocation is always generated at the same spanwise location. We estimate $\widetilde {St}_v$ from each time period of the time history of the cross-flow component of velocity at various spanwise locations in the near wake ($x/c = 1.25$, $y/c = -0.0068$). We note that $\widetilde {St}_v$ is the unsteady vortex shedding frequency estimated from each time period while that shown in figure 12 is from the fast Fourier transform of the entire time history of the unsteady signal. The spatio-temporal variation of $\widetilde {St}_v$ is presented in figure 14(a). Vortex dislocations can be identified from a spatial change in the local vortex shedding frequency, $\widetilde {St}_v$ (Mittal et al. Reference Mittal, Pandi and Hore2021). Behara & Mittal (Reference Behara and Mittal2010) further observed that the appearance of vortex dislocations coincides with local minima in $\widetilde {St}_v$ followed by its sharp rise. Three dislocations, $D_1$, $D_2$ and $D_3$, are marked in figure 14(a). They appear at $t = 175.23$, $196.23$ and $217.33$, respectively, at the same spanwise location. The time period of formation of the dislocations is ${\sim }21$ time units and corresponds to a frequency of $0.0476$. The frequency of appearance of dislocations can also be estimated from the beat frequency in the time signal near the junction of the cells (Williamson Reference Williamson1989; Mittal et al. Reference Mittal, Pandi and Hore2021). The beat frequency estimated from the time history of the cross-flow component of velocity is $0.0458$, and is in good agreement with the estimate from the variation of $\widetilde {St}_v$.

Figure 14.

Flow past the $sAR = 5$ wing at $\alpha =10^{\circ }$ and $Re = 1000$: (a) spatio-temporal variation of $\widetilde {St}_v$ obtained using a probe placed at $x/c = 1.25$ and $y/c = -0.068$. The $Q (=0.1)$ isosurface for an instantaneous flow coloured with the spanwise component of vorticity ($\omega _z = \pm 2$) at two time instants is presented in (b,c). Fork-type ($D_f$) and connected fork-type ($D_{fc}$) dislocation is marked in (b,c). Results are shown for (b) $t = 178.386$, (c) $t = 180.886$.

Images of the isosurface of $Q$ coloured with the spanwise component of vorticity are presented in figure 14(b,c) at two time instants. The dislocation is of a fork type at $t = 178.386$ (figure 14b). At $t = 180.886$ (figure 14c), a linkage between vortices of opposite polarity results in the formation of a ring-type vortex structure in addition to a fork-type dislocation. Mittal et al. (Reference Mittal, Pandi and Hore2021) referred to this as the mixed-type dislocation in their investigation of flow past a cylinder in the presence of a side wall. To the best of our knowledge, this type of dislocation is being reported for the first time in the wake of wings.

6.3. Fork-type ($D_f$) dislocation generated at varying spanwise locations: $sAR =5$, $\alpha = 14^{\circ }$

Dislocations are formed at varying spanwise locations for $sAR \ge 4$ at $\alpha = 14^{\circ }$ (see figure 13). The spatio-temporal variation of $\widetilde {St}_v$ for $sAR = 5$ is shown in figure 15(a). Dislocations formed at various time instants are marked in the figure as $D_1-D_6$. Each occurrence is at a different spanwise location between $z/c=0.89$ to $1.9$. The time period between successive occurrence also varies. The average dislocation frequency is estimated to be $0.091$. It is in good agreement with the beat frequency ($=0.092$) estimated from the time signal of the cross-flow component of velocity in the near wake. The instantaneous flow is shown in figure 15(b–e). The first two images show dislocation $D_1$ as it convects downstream, while the other two images show $D_3$. Both $D_1$ and $D_3$ are of a fork type ($D_f$) and appear at different spanwise locations resulting in a time-varying span length of the end and central cells. However, once generated, each dislocation convects along the same spanwise location. Unlike in the present study, the dislocation at lower $Re (= 400)$ undergoes a spanwise movement as it convects (Zhang et al. Reference Zhang, Hayostek, Amitay, He, Theofilis and Taira2020).

Figure 15.

Flow past the $sAR = 5$ wing at $\alpha = 14^{\circ }$ and $Re = 1000$: (a) spatio-temporal variation of $\widetilde {St}_v$ obtained using a probe placed at $x/c = 1.25$ and $y/c = -0.068$. The $Q (=0.1)$ isosurface for an instantaneous flow coloured with the spanwise component of vorticity ($\omega _z = \pm 2$) at various time instants is presented in (b–e). Also marked are the dislocation $D_1$ in (b–c) and $D_3$ in (d–e). Results are shown for (b) $t = 60.176$, (c) $t = 61.842$, (d) $t = 77.994$, (e) $t = 79.298$.

6.4. Fork ($D_f$) and reverse fork-type ($D_{rf}$) dislocations due to vortex splitting and reconnection: $sAR = 5$, $\alpha = 12^{\circ }$

Unlike for other flows, $\alpha =12^{\circ }$ exhibits formation of dislocations in the far wake ($x/c \sim 8$). The spatio-temporal variation of $\widetilde {St}_v$ for a probe in the near wake shows a single cell shedding while it shows the appearance of dislocations rather close to the wing root in the far wake. Figure 16(a) shows that the vortex shedding frequency at the wing root shows periodic variations with time while it is approximately constant with time, over the outer span. The instantaneous flow presented in figure 16(b) shows that the vortices are released from the wing at an oblique angle in the wake. They tend to align themselves parallel to the wing, especially near the wing root, as they convect downstream. In the process they undergo splitting and reconnection leading to a dislocation in the far wake. A reverse fork-type dislocation ($D_{rf}$), on the far right, is marked in the image. The dislocation forms when two vortices of the same polarity in the end cell connect to a vortex of the same polarity in the adjoining cell. Some of the vortices marked in the near wake in this panel are identified in figure 16(c) at a later time as they convect downstream. While the dislocation seen in figure 16(b) is washed out of the frame at a later time, a new dislocation appears in figure 16(c) as a result of splitting and reconnection of vortices as shown. Also, a fork-type dislocation ($D_{f}$) forms in the later part of the wake. The convection of the vortices and the dislocations lead to the state shown in figure 16(d). Dislocation $D_{f}$ is washed out, while $D_{rf}$ moves in the far wake.

Figure 16.

Flow past the $sAR = 5$ wing at $\alpha = 12^{\circ }$ and $Re = 1000$: (a) spatio-temporal variation of $\widetilde {St}_v$ obtained using a probe placed at $x/c = 1.25$ and $y/c = -0.068$. The $Q (=0.05)$ isosurface for an instantaneous flow coloured with the spanwise component of vorticity ($\omega _z = \pm 1$) at various time instants is presented in (b–d). The vortex splitting and reconnections are illustrated using coloured solid lines. The filled black circles show the inter-connection between the vortices. Here $D_f$ and $D_{rf}$ refer to fork-type and reverse fork-type dislocations, respectively. Results are shown for (b) $t = 189.61$, (c) $t = 192.61$,(d) $t = 200.61$.

6.5. Role of cellular shedding in spanwise variation of r.m.s. of aerodynamic coefficients

We explore the effect of cellular shedding observed for $sAR \ge 3$ (see figure 8) on the unsteadiness of flow via spanwise variation of the r.m.s. of the lift and drag coefficient and $\overline {u'v'}$ component of Reynolds stress in the near wake. The same is shown in figure 17 along with the flow for $sAR=2$ and $4$ associated with a single cell and two cells, respectively. The r.m.s. coefficients and $\overline {u'v'}$ are maximum close to the wing root for $sAR=2$. However, the unsteadiness in the flow is closer to the wing tip than the wing root for $sAR=4$. For example, $C_{l_{rms}}$ achieves a maxima at $z/c\sim 2.6$. We note the nearly parallel shedding for $sAR=2$ and varying oblique angle of the vortices in the near wake for $sAR=4$. The spanwise minima in the unsteadiness for $sAR=4$ occurs at the location where the oblique angle of the shed vortices is maximum. The unsteadiness near the wing tip is very weak for both $sAR$.

Figure 17.

Flow past a $\alpha = 14^{\circ }$ wing at $Re = 1000$: (a,b) spanwise variation of sectional r.m.s. of lift and drag coefficient, (c,d) $Q (=0.1)$ isosurface for an instantaneous flow coloured with the spanwise component of vorticity ($\omega _z = \pm 2$) and (e, f) $\overline {u'v'}$ component of Reynolds stress on the $y\unicode{x2013}z$ plane at $x/c = 3.5$. The quantities are for $sAR = 2$ (left) and $sAR = 4$ (right).

7.1. Do streamwise vortices exist in a steady flow?

Figure 9(d) shows that, unlike the prediction from lifting line theory presented in figure 7(a), the sectional lift coefficient for the steady flow also shows a peak closer to the wing tip and not at the wing root. A natural question that arises is: ‘Is the steady flow also associated with streamwise vortices other than wing-tip vortices?’ The $Q$ isosurface, coloured with the streamwise component of vorticity, is shown in figure 20 for the steady flow for the $sAR=5$ wing at various $\alpha$. The flow is obtained by dropping the unsteady terms from the flow equations. The time-averaged and steady flows are very similar for $\alpha \le 8^{\circ }$ (not shown here) as the vortex shedding is relatively weak; the wing-tip vortices are the only prominent streamwise vortices. Additional streamwise vortices appear near the wing tip at larger $\alpha$. As seen from figure 20, they increase in number, originate from a larger length of the span and become stronger with an increase in $\alpha$. At $\alpha =14^{\circ }$ they virtually occupy the entire span (figure 20c). Interestingly, for a wing with larger $sAR$ (for example, $sAR=7.5$, shown in figure 20d), the spanwise extent of the generation of streamwise vortices is similar to that for $sAR=5$ (figure 20c).

Figure 20.

Flow past the $sAR = 5$ wing at $Re = 1000$: $Q (=0.01)$ isosurface for a steady flow coloured with the streamwise component of vorticity ($\omega _x = \pm 2$) for $\alpha =$ (a) $10^{\circ }$, (b) $12^{\circ }$, (c) $14^{\circ }$. Also presented in (d) is the flow of the $sAR = 7.5$ wing at $\alpha = 14^{\circ }$.

7.2. Effect of streamwise vortices on spanwise lift distribution

We attempt to explain the local maximum observed in the spanwise distribution of the sectional lift coefficient by using a discrete model of the streamwise vortices. Figure 21(a,b) shows the cores of streamwise vortices corresponding to the steady and time-averaged flows for $\alpha = 14^{\circ }$ and the $sAR = 5$ wing. Also shown in this schematic is the bound vortex on the wing. The circulation of each streamwise vortex is estimated by carrying out line integration of the velocity field along a closed curve, enclosing the vortex, in the $y\unicode{x2013}z$ plane at $x/c=3.5$. To enable a good estimate, the closed curve is chosen to be as large as possible while ensuring that it does not include adjacent vortices. The circulation is non-dimensionalized using the free-stream speed and chord of the airfoil. The streamwise component of vorticity in the $y\unicode{x2013}z$ plane at $x/c = 3.5$ is presented in figure 21(d,e) for the steady and time-averaged flows, respectively. A streamwise vortex with a component of its vorticity pointing in the positive $x$ direction is considered positive and is marked in red in the schematics shown in figure 21(a,b). The circulation of each vortex ‘$i$’ is denoted by $\varGamma _i$ with the index increasing from the wing tip towards the mid-span. The indices for all the major streamwise vortices are marked in figure 21(a–e) and their respective $\varGamma _i$ plotted in figure 21(f, g); $\varGamma _1$, for example, denotes the strength of the wing-tip vortex. Also shown in figure 21(f, g) is the circulation, $\varGamma _z$, along the bound vortex. It is estimated at various spanwise locations via line integration of the velocity field on a closed curve in the $x\unicode{x2013}y$ plane that encloses the airfoil section. We utilize a circle whose centre is located at the mid-chord of the airfoil. The value of $\varGamma _z$ estimated on various circular curves with radius $0.51c-1c$ is very close. The sign convention is indicated in figure 21(c). We recall that $\varGamma _z$ is related to the sectional lift coefficient via the Kutta–Joukowski theorem for 2-D flows (Bertin Reference Bertin2002; Anderson Reference Anderson2017). A reasonable match is seen between $\varGamma _z$ and $C_l/2$ for both the time-averaged and steady flows except in the region of the wing tip, where, expectedly, the flow is highly three dimensional.

Figure 21.

Flow past the $sAR = 5$ wing at $\alpha = 14^{\circ }$ and $Re = 1000$: schematic showing the bound vortex and core of streamwise vortices including the wing-tip vortex for a (a) steady and (b) time-averaged flow. The wing surface is shown in a grey colour and the bound vortex in a solid black line. The sign convention for the circulation around the bound vortex is also shown in (c). Streamwise component of vorticity on the $y\unicode{x2013}z$ plane at $x/c = 3.5$ is shown for asteady flow in (d) and a time-averaged flow in (e). Shown in ( f,g) is the spanwise variation of circulation ($\varGamma _z$) in the $x\unicode{x2013}y$ plane enclosing the airfoil section along with the sectional lift coefficient ($C_l$) for the steady and time-averaged flows, respectively. Also shown in the same plot is the circulation of streamwise vortices $\varGamma _i$ and the strength of the bound vortex ($\varGamma _b$).

Using the discrete vortex model, as shown in figure 21(a,b), we estimate the strength of the bound vortex at a spanwise location $z_i$ by the expression $\varGamma _b(z_i) = \sum _{j=1}^{i}(-\varGamma _j)$. The same is shown in figure 21( f,g). For both, steady and time-averaged flows, $\varGamma _b(z_i)$ shows qualitatively the same behaviour as $\varGamma _z$ and $C_l/2$ confirming that the local peak in the spanwise variation of the sectional lift coefficient is due to streamwise vortices. Figure 21(a,b) shows that the streamwise vortices that are fairly evenly spread over the span of the wing in a steady flow relocate away from the root of the wing for the time-averaged flow, leading to an increased interaction and vortex strength. In general, the strength of vortices in the time-averaged flow is larger than in the steady flow. For example, $\varGamma _1$ and $\varGamma _2$ for the time-averaged flow are larger compared with the steady flow, even though their locations are quite comparable.

8.1. Strength of wing-tip vortex

We explore the dependence of the strength of the wing-tip vortex with angle of attack and aspect ratio of the wing. The circulation is estimated by carrying out integration along a closed curve lying in the $y\unicode{x2013}z$ plane at $x/c=3.5$. The effect of the spatial extent of the curve on $\varGamma$ is assessed for the time-averaged flow at $\alpha =14^{\circ }$ and presented in Appendix D. Isovorticity curves ($\overline {\omega _x}=C$) are utilized to carry out the integration. It is found that a curve corresponding to a value of $|C| \le 0.1$ provides a good estimate of $\varGamma$. The same has been utilized in this study. Figure 22(a) shows the variation of circulation, corresponding to the steady and time-averaged flows for the $sAR=5$ wing, with angle of attack. Here $\varGamma$ increases with an increase in $\alpha$. The variation is linear for the steady flow. However, the slope of the curve changes at $\alpha =8^{\circ }$; it is larger for the variations beyond $\alpha =8^{\circ }$. We recall that the streamwise vortices, other than the wing-tip vortices, become prominent beyond $\alpha =8^{\circ }$ (see figure 20). We attribute this change in slope of $\varGamma -\alpha$ variation to the appearance of additional streamwise vortices across the wing span. For the time-averaged flow, the linear variation up to $\alpha =8^{\circ }$ is followed by a nonlinear increase of $\varGamma$ with $\alpha$. In this regime the strength of the wing-tip vortex for the time-averaged flow is higher than that of the steady flow. This is also seen from figure 21( f,g) for $\alpha =14^{\circ }$. The effect of the aspect ratio of the wing on the strength of the wing-tip vortex is presented in figure 22(b) for $\alpha = 8^{\circ }$ and $14^{\circ }$. For $\alpha =8^{\circ }$, the vortex strength increases with an increase in $sAR$ up to $sAR=2$ and does not show any significant change beyond it. The variation in $\varGamma$ is more significant at $\alpha =14^{\circ }$. It increases with $sAR$ up to $sAR=5$ and does not show any appreciable increase beyond it. For each $sAR$, the strength of the wing-tip vortex for the time-averaged flow is larger than that for the steady flow.

Figure 22.

Flow past a wing at $Re = 1000$: variation of strength of the wing-tip vortex, $\varGamma$, with (a) $\alpha$ for $sAR = 5$ and (b) $sAR$ for $\alpha = 8^{\circ }$ and $14^{\circ }$. Here $\varGamma$ is non-dimensionalized using a free-stream speed and chord of the airfoil.

8.2. Model of wing-tip vortex

We utilize two models to estimate the radius of the wing-tip vortex: Rankine vortex (Rankine Reference Rankine1877) and Lamb–Oseen vortex (Saffman Reference Saffman1995; Jacquin et al. Reference Jacquin, Fabre, Sipp and Coustols2005). The core of a Rankine vortex is modelled by a forced vortex of strength $\varGamma$ and radius $R$ while the outer zone is approximated by a free vortex. The Lamb–Oseen vortex is an exact solution of the Navier–Stokes equation. The velocity across the vortex decreases exponentially beyond its core of radius $R$. The tangential velocity distributions for the Rankine and Lamb–Oseen vortex are, respectively,

(8.1a,b)\begin{gather} V_{\theta_{fo}} = \frac{\varGamma r}{2{\rm \pi} R^2}; V_{\theta_{fr}} = \frac{\varGamma}{2 {\rm \pi}r}; \end{gather}

(8.2)\begin{gather} V_{\theta_{LO}} = \frac{\varGamma}{2{\rm \pi} r}(1-e^{{-}r^2/R^2}). \end{gather}

A schematic of the two models is shown in figure 23(a). The models are applied to the time-averaged flow past the $sAR = 5$ wing at $\alpha = 14^{\circ }$. Figure 23(b,c) shows the velocity distribution in the $y\unicode{x2013}z$ plane along lines parallel to the $y$ and $z$ axes and passing through the core of the wing-tip vortex at two streamwise locations $x/c = 3.5$ and $10$, respectively. Also shown in the figures, with broken lines, is the velocity distribution from the models given by (8.1) and (8.2). The value of $\varGamma$ is estimated from the computations (see figure 22). The value of $R$ is estimated from each of the two models as shown in figure 23(b,c) for the best approximation to the velocity profiles from direct time integration. Here $R$ is estimated using this procedure for various streamwise locations. The streamwise variation of $R$ alongside the best fit is shown in figure 23(d). We note that both the models yield virtually the same value of $R$. It increases linearly with an increase in $x/c$. A similar variation is also seen for the steady flow (not shown here). It is found that $R$ is larger for the time-averaged flow compared with the steady flow. Figure 24(a,b) shows the variation of $R$, at $x/c = 3.5$, with $\alpha$ and $sAR$, respectively. Here $R$ increases with an increase in the angle of attack. For a fixed $\alpha$, $R$ increases with an increase in $sAR$ up to $sAR=2$, beyond which it does not change with any further increase in $sAR$.

Figure 23.

Flow past a finite wing for $Re = 1000$: schematic of the Rankine (forced-free) vortex and Lamb–Oseen vortex model where $R$, $\varGamma$, $o'$, $V_{\theta _{fo}}$, $V_{\theta _{fr}}$ and $V_{\theta _{LO}}$ represents the radius of the vortex, strength of the wing-tip vortex, core of the vortex, velocity for the Rankine (forced-free) vortex and Lamb–Oseen vortex, respectively. The variation of velocity distribution across the wing-tip vortex for a time-averaged flow of an $sAR = 5$ wing at $\alpha = 14^{\circ }$ is shown in (b) $x/c = 3.5$ and (c) $x/c = 10$. Also shown in broken lines are the curve obtained using the theoretical expressions. The variation of $R/c$, obtained using the Rankine and Lamb–Oseen vortex, with $x/c$ along with best linear fit is presented in (d).

Figure 24.

Flow past a wing at $Re = 1000$: variation of $R/c$ estimated at $x/c = 3.5$ with (a) $\alpha$ for an $sAR = 5$ wing and (b) $sAR$ at $\alpha = 8^{\circ }$ and $14^{\circ }$.