Hostname: page-component-78c5997874-mlc7c Total loading time: 0 Render date: 2024-11-05T07:53:52.718Z Has data issue: false hasContentIssue false

Laminar post-stall wakes of tapered swept wings

Published online by Cambridge University Press:  24 November 2023

Jean Hélder Marques Ribeiro*
Affiliation:
Department of Mechanical and Aerospace Engineering, University of California, Los Angeles, CA 90095, USA
Jacob Neal
Affiliation:
Department of Mechanical, Aeronautical, and Nuclear Engineering, Rensselaer Polytechnic Institute, Troy, NY 12180, USA
Anton Burtsev
Affiliation:
Department of Mechanical and Aerospace Engineering, University of Liverpool, Brownlow Hill, Liverpool L69 3GH, UK
Michael Amitay
Affiliation:
Department of Mechanical, Aeronautical, and Nuclear Engineering, Rensselaer Polytechnic Institute, Troy, NY 12180, USA
Vassilios Theofilis
Affiliation:
Department of Mechanical and Aerospace Engineering, University of Liverpool, Brownlow Hill, Liverpool L69 3GH, UK
Kunihiko Taira
Affiliation:
Department of Mechanical and Aerospace Engineering, University of California, Los Angeles, CA 90095, USA
*
Email address for correspondence: [email protected]

Abstract

While tapered swept wings are widely used, the influence of taper on their post-stall wake characteristics remains largely unexplored. To address this issue, we conduct an extensive study using direct numerical simulations to characterize the wing taper and sweep effects on laminar separated wakes. We analyse flows behind NACA 0015 cross-sectional profile wings at post-stall angles of attack $\alpha =14^\circ$$22^\circ$ with taper ratios $\lambda =0.27$$1$, leading-edge sweep angles $0^\circ$$50^\circ$ and semi aspect ratios $sAR =1$ and $2$ at a mean-chord-based Reynolds number of $600$. Tapered wings have smaller tip chord length, which generates a weaker tip vortex, and attenuates inboard downwash. This results in the development of unsteadiness over a large portion of the wingspan at high angles of attack. For tapered wings with backward-swept leading edges, unsteadiness emerges near the wing tip. On the other hand, wings with forward-swept trailing edges are shown to concentrate wake-shedding structures near the wing root. For highly swept untapered wings, the wake is steady, while unsteady shedding vortices appear near the tip for tapered wings with high leading-edge sweep angles. For such wings, larger wake oscillations emerge near the root as the taper ratio decreases. While the combination of taper and sweep increases flow unsteadiness, we find that tapered swept wings have more enhanced aerodynamic performance than untapered and unswept wings, exhibiting higher time-averaged lift and lift-to-drag ratio. The current findings shed light on the fundamental aspects of flow separation over tapered wings in the absence of turbulent flow effects.

Type
JFM Papers
Copyright
© The Author(s), 2023. Published by Cambridge University Press

1. Introduction

Flow separation over aerodynamic lifting bodies has been a subject of research interest for decades, especially for small-scale air vehicles (Mueller Reference Mueller2001; Anderson Reference Anderson2010). To further understand the post-stall wake dynamics, it is important to analyse the influence of wing planform geometry. This characterization is challenging for high-Reynolds-number flows, due to the multiscale nature of the wakes. Nevertheless, for massively separated flows, the large vortex structures observed in higher-Reynolds-number flows are topologically analogous to the core structures in low-Reynolds-number flows (Hunt et al. Reference Hunt, Abel, Peterka and Woo1978; Dallmann Reference Dallmann1988; Délery Reference Délery2001). To examine the fundamental aspects of unsteady three-dimensional ($3$-D) flow separation, we study post-stall flows in the absence of turbulence (Taira & Colonius Reference Taira and Colonius2009; Zhang et al. Reference Zhang, Hayostek, Amitay, Burstev, Theofilis and Taira2020a,Reference Zhang, Hayostek, Amitay, He, Theofilis and Tairab). This characterization has been largely unexplored for low-Reynolds-number flows over tapered wings.

In aircraft design, tapered wings are used to approximate the elliptic aerodynamic loading over the wingspan. Tapered wings are more feasible to manufacture due to their less complex geometry compared with elliptic wings (Prandtl Reference Prandtl1920; McCormick Reference McCormick1995). The usage of tapered wings in aeronautics led to initial studies that explored the wing taper effect, especially for high-Reynolds-number flows (Anderson Reference Anderson1936; Millikan Reference Millikan1936; Irving Reference Irving1937; Soule & Anderson Reference Soule and Anderson1940; Falkner Reference Falkner1950). For the laminar flow regime, the effect of wing taper on the wake dynamics is critical as the local Reynolds number is drastically reduced near the tip. For flows over wings at a chord-based Reynolds number $Re_c = O(10^4)$, taper affects the aerodynamic loading with an increase in the pressure drag (Traub Reference Traub2013; Traub et al. Reference Traub, Botero, Waghela, Callahan and Watson2015). For $Re_c = O(10^3)$, the aerodynamic characteristics are affected significantly by the viscous effects and the influence of wing taper on the wakes remains elusive, especially for massively separated flows.

Post-stall wake dynamics has attracted the attention of aeronautical researchers for many decades. The early efforts to understand post-stall flows over wings were performed over two-dimensional ($2$-D) spanwise homogeneous wings (Abbott & Von Doenhoff Reference Abbott and Von Doenhoff1959; Gaster Reference Gaster1967; Tobak & Peake Reference Tobak and Peake1982). Valuable insights were obtained from $2$-D analysis characterizing the behaviour of the separated laminar boundary layer (Horton Reference Horton1968) and describing the relation between vortex-shedding structures, adverse pressure gradient and shear-layer characteristics (Pauley, Moin & Reynolds Reference Pauley, Moin and Reynolds1990). Moreover, the emergence of wake patterns associated with 3-D separation bubbles, as predicted in topological studies (Hornung & Perry Reference Hornung and Perry1984; Perry & Hornung Reference Perry and Hornung1984), was shown by global linear stability analysis to arise from self-excitation of the laminar separation bubble (Theofilis, Hein & Dallmann Reference Theofilis, Hein and Dallmann2000).

The analysis of $2$-D flows around canonical wings continues providing fundamental insights into the effect of angle of attack and Reynolds number on the wake-shedding structures (Lin & Pauley Reference Lin and Pauley1996; Huang et al. Reference Huang, Wu, Jeng and Chen2001; Yarusevych, Sullivan & Kawall Reference Yarusevych, Sullivan and Kawall2009; Rossi et al. Reference Rossi, Colagrossi, Oger and Le Touzé2018; Durante, Rossi & Colagrossi Reference Durante, Rossi and Colagrossi2020). For separated flows, an increase in Reynolds number and the angle of attack yields a $3$-D flow field even around infinite and spanwise homogeneous wings (Bippes & Turk Reference Bippes and Turk1980; Winkelman & Barlow Reference Winkelman and Barlow1980; Braza, Faghani & Persillon Reference Braza, Faghani and Persillon2001; Schewe Reference Schewe2001; Hoarau et al. Reference Hoarau, Braza, Ventikos, Faghani and Tzabiras2003; Pandi & Mittal Reference Pandi and Mittal2019). In such cases, spanwise fluctuations emerge, producing $3$-D vortices in the wake, as a result of the growth of $3$-D structures associated with secondary linear instability (He et al. Reference He, Gioria, Pérez and Theofilis2017a).

For finite wings, $3$-D wakes result from tip effects, as a strong streamwise vortex is formed rolling up around the wing tip (Winkelman & Barlow Reference Winkelman and Barlow1980; Freymuth, Finaish & Bank Reference Freymuth, Finaish and Bank1987; Toppings & Yarusevych Reference Toppings and Yarusevych2022). While turbulence has an important influence on the $3$-D wake (Pandi & Mittal Reference Pandi and Mittal2019), some of the core global flow structures remain coherent over a broad range of Reynolds numbers, including the quasi-spanwise midspan shedding and the tip vortex (Neal & Amitay Reference Neal and Amitay2023; Pandi & Mittal Reference Pandi and Mittal2023). Tip vortices induce downwash inboard over the wing, which reduces the effective angle of attack near the tip, even suppressing stall formation (Dong, Choi & Mao Reference Dong, Choi and Mao2020; Toppings & Yarusevych Reference Toppings and Yarusevych2021) and the wake shedding for low-aspect-ratio wings (Taira & Colonius Reference Taira and Colonius2009; Zhang et al. Reference Zhang, Hayostek, Amitay, He, Theofilis and Taira2020b). The tip vortex has been extensively studied to reveal its influence on the wake dynamics, aerodynamic forces and pitch moments (Francis & Kennedy Reference Francis and Kennedy1979; Green & Acosta Reference Green and Acosta1991; Devenport et al. Reference Devenport, Rife, Liapis and Follin1996; Pelletier & Mueller Reference Pelletier and Mueller2000; Birch et al. Reference Birch, Lee, Mokhtarian and Kafyeke2004; Torres & Mueller Reference Torres and Mueller2004; Buchholz & Smits Reference Buchholz and Smits2006; Yilmaz & Rockwell Reference Yilmaz and Rockwell2012; Ananda, Sukumar & Selig Reference Ananda, Sukumar and Selig2015; He et al. Reference He, Tendero, Paredes and Theofilis2017b). Beyond understanding the tip vortex formation and evolution, a characterization of its instabilities has enabled the development of control techniques that improve the aerodynamic performance around finite wings (Edstrand et al. Reference Edstrand, Sun, Schmid, Taira and Cattafesta2018; Gursul & Wang Reference Gursul and Wang2018; Navrose, Brion & Jacquin Reference Navrose, Brion and Jacquin2019).

Wing sweep also has a strong influence on post-stall wake dynamics. For laminar flow regimes, a number of experimental and numerical efforts were made to examine the effects of backward and forward wing sweep (Yen & Hsu Reference Yen and Hsu2007; Yen & Huang Reference Yen and Huang2009; Zhang et al. Reference Zhang, Hayostek, Amitay, Burstev, Theofilis and Taira2020a) and identify global modes (Burtsev et al. Reference Burtsev, He, Hayostek, Zhang, Theofilis, Taira and Amitay2022; Ribeiro, Yeh & Taira Reference Ribeiro, Yeh and Taira2023b) that give rise to fundamental global coherent structures of flow separation around swept wings. It is noteworthy that some of the wing-sweep effects on laminar post-stall flows are topologically analogous over a wide range of Reynolds numbers. For instance, for low-Reynolds-number flows over swept wings, some of the core coherent structures emerging in the near wake, such as the ‘ram's horn’ vortex and the canard leading-edge vortices, were also identified in experiments performed at higher Reynolds numbers (Black Reference Black1956; Breitsamter & Laschka Reference Breitsamter and Laschka2001; Neal et al. Reference Neal, Burtsev, Ribeiro, Taira, Theofilis and Amitay2023a). Moreover, the stabilizing effect of the sweep-induced spanwise flow on the wake structures, which significantly impacts stall characteristics, as observed at low-Reynolds-number flows (Zhang et al. Reference Zhang, Hayostek, Amitay, Burstev, Theofilis and Taira2020a; Ribeiro et al. Reference Ribeiro, Yeh, Zhang and Taira2022Reference Ribeiro, Yeh and Taira2023b), is further noticed in both experiments and high-fidelity large-eddy simulations performed at a higher-Reynolds-number regime (Harper & Maki Reference Harper and Maki1964; Visbal & Garmann Reference Visbal and Garmann2019).

The aforementioned studies highlight the importance of the low-Reynolds-number post-stall wake characterization for revealing fundamental aspects of the unsteady 3-D flow separation physics. In fact, the insights obtained from studies of post-stall laminar flows have been important for expanding our knowledge of the stalled flow physics over a wide range of Reynolds numbers. Thus far, however, most studies have not considered wing-taper effects on low-Reynolds-number flows at high angles of attack. Only recently, a combined experimental, numerical and theoretical effort has been initiated towards the understanding of the laminar flow over tapered wings in post-stall flow conditions (Burtsev et al. Reference Burtsev, Theofilis, Ribeiro, Taira, Neal and Amitay2023; Neal et al. Reference Neal, Gares, Amitay, Burtsev, Theofilis, Ribeiro and Taira2023b; Ribeiro et al. Reference Ribeiro, Taira, Neal, Amitay, Burtsev and Theofilis2023a). Effects of taper have been analysed for planforms with tubercles to analyse swimming of whales (Wei, New & Cui Reference Wei, New and Cui2018), for flows over tapered cylinders (Piccirillo & Van Atta Reference Piccirillo and Van Atta1993; Techet, Hover & Triantafyllou Reference Techet, Hover and Triantafyllou1998; Valles, Andersson & Jenssen Reference Valles, Andersson and Jenssen2002) and for separated wakes over tapered plates (Narasimhamurthy, Andersson & Pettersen Reference Narasimhamurthy, Andersson and Pettersen2008). For wing planforms with continuously variable chord length over the wingspan, delta wings have also received substantial attention (Rockwell Reference Rockwell1993; Gursul, Gordnier & Visbal Reference Gursul, Gordnier and Visbal2005; Taira & Colonius Reference Taira and Colonius2009). For laminar post-stall flows, wing taper was studied using trapezoidal plates (Huang et al. Reference Huang, Venning, Thompson and Sheridan2015). Nonetheless, there still is a lack of fundamental studies for understanding the role of taper ratio, and how it interplays with leading-edge (LE) and trailing-edge (TE) sweep-angle effects for massively separated laminar flows.

For laminar separated flows, the combined effect of wing taper and sweep remains elusive. In the present work, we aim to reveal the effects of taper in the laminar wake dynamics and the influence of LE and TE sweep angles on the vortical interactions through a comprehensive campaign of direct numerical simulations of $3$-D flows over finite NACA 0015 wings. We characterize the stalled wakes of wings with backward-swept LE and forward-swept TE, identifying the combined effects of taper and sweep on the post-stall wake dynamics. Our work is organized as follows. In § 2, we present our wing planform geometry definitions and the set-up for direct numerical simulations. In § 3, we offer a detailed analysis and classification of the wake structures, highlighting the effects of taper and sweep on the wakes and aerodynamic forces. Finally, we conclude our study by summarizing our findings in § 4.

2. Problem set-up

We consider laminar flows over tapered wings with a NACA 0015 cross-sectional profile. The spatial coordinates of streamwise, transverse and spanwise directions are denoted by $(x,y,z)$, respectively. The origin is placed at the LE of the wing root, as shown in figure 1. The NACA 0015 profile is defined on the $(x,y)$ plane, which is extruded from the wing root in the spanwise direction to form the $3$-D wing. Wing taper is defined by the taper ratio $\lambda = c_{tip} / c_{root}$, where $c_{tip}$ and $c_{root}$ are tip and root chord lengths, respectively, as shown in figure 1(a). For all wings considered herein, the chord length decreases linearly from root to tip. The non-dimensional mean chord length $c$ at the spanwise location of $z = b/2$ is taken to be the characteristic length used to non-dimensionalize all spatial variables. The mean chord $c$ is fixed and independent of $\lambda$ for all wing planforms studied herein.

Figure 1. Problem set-up for tapered wings. (a) Geometrical parameters shown in a wing planform with $sAR = b/c = 2$, $\alpha = 18^\circ$, $\lambda = 0.27$ and $\varLambda _{LE} = 18.4^\circ$. (b) For a $(\lambda,\varLambda _{LE}) = (0.5,0^\circ )$ wing, we show the computational domain and (c,d) grids with $2$-D planes at $z/c = 1$ and $y/c = -0.5$, respectively.

The semi aspect ratio of the wings is set as $sAR = b/c = 1$ and $2$, where $b$ is the half-span length, as shown in figure 1(d). We consider half-span wing models with symmetry imposed at the root. The angles of attack, $\alpha = 14^\circ$, $18^\circ$ and $22^\circ$, are defined between the airfoil chord line and the streamwise direction. The present wing geometries have sharp TE and straight-cut wing tip. The mean-chord-based Reynolds number is set to $Re_{c} = 600$ and the free-stream Mach number is set to $M_\infty = 0.1$. Taper changes the local Reynolds number $Re_{L_c}$, defined as a function of the spanwise location (Traub et al. Reference Traub, Botero, Waghela, Callahan and Watson2015). For the present study, the difference between $c_{tip}$ and $c_{root}$ accounts for a maximum variation of $60\,\%$ on $Re_{L_c}$ along the span, from $\min (Re_{L_c}) = 250$ to $\max (Re_{L_c}) = 950$ at the lowest taper ratio.

For tapered swept wings, the $3$-D computational set-up is sheared in the chordwise direction and the LE sweep angle is defined between the $z$ direction and the LE. Tapered wings have different LE and TE sweep angles ($\varLambda _{LE}$ and $\varLambda _{TE}$, respectively), as shown in figure 1(a). Note that the wing planform can be specified with two parameters out of the three parameters of taper ratio ($\lambda$), LE and TE sweep. Given that $(c_{tip}+c_{root})/2 = c$, for a chosen $\lambda$, $\varLambda _{LE}$ and $sAR$ we have

(2.1ac)\begin{align} \lambda = \frac{c_{tip}}{c_{root}}, \quad \frac{c_{root}}{c}= \frac{2}{1 + \lambda} \quad \text{and} \quad \varLambda_{TE} = \arctan\left[ -\frac{2}{sAR}\left(\frac{1 - \lambda}{1 + \lambda}\right) + \tan{(\varLambda_{LE})} \right]. \end{align}

In this work, we explore the combined effects of the LE and TE sweep angles on the wake dynamics for LE sweep angles $0 \le \varLambda _{LE} \le 50^\circ$ and taper ratios $0.27 \le \lambda \le 1$. The corresponding TE sweep angles take $-30^\circ \le \varLambda _{TE} \le 50^\circ$. Herein, negative sweep angles indicate a forward sweep, as shown in figure 1(a), while a positive sweep angle represents a backward sweep.

Traditionally in aeronautics, tapered swept wings have wing-sweep angles observed with respect to the quarter-chord line (Anderson Reference Anderson1936Reference Anderson2010; Falkner Reference Falkner1950) denoted by $\varLambda _{c/4}$, as shown in figure 1(a). Anderson (Reference Anderson1999) considered the half-chord sweep angle $\varLambda _{c/2}$, such that aerodynamic load distribution becomes independent of the taper ratio. Straight tapered wings with $\varLambda _{c/4} = 0^\circ$ were studied by Traub et al. (Reference Traub, Botero, Waghela, Callahan and Watson2015). On the other hand, Irving (Reference Irving1937) considered the effect of the LE and TE sweep angles. For the present laminar post-stall wakes, due to the crucial role played by the LE vortex in defining the wake characteristics (Videler, Stamhuis & Povel Reference Videler, Stamhuis and Povel2004; Eldredge & Jones Reference Eldredge and Jones2019), we focus on the distinct effects of $\varLambda _{LE}$ and $\varLambda _{TE}$ in our analysis and describe their influence on the wake dynamics. We note, however, that it is also possible to translate the findings reported herein with respect to the traditional quarter-chord and half-chord sweep angles, $\varLambda _{c/4}$ and $\varLambda _{c/2}$, respectively.

2.1. Direct numerical simulations

We conduct direct numerical simulations with a compressible flow solver, CharLES (Khalighi et al. Reference Khalighi, Ham, Nichols, Lele and Moin2011; Brès et al. Reference Brès, Ham, Nichols and Lele2017), which uses a second-order-accurate finite-volume method in space with a third-order-accurate total-variation diminishing Runge–Kutta scheme for time integration. The computational domain is discretized with a C-type grid with mesh refinement near the wing and in the wake. With the origin at the airfoil LE on the symmetry plane $(x/c, y/c,z/c) = (0,0,0)$, the computational domain extends over $(x/c, y/c, z/c) \in [-20,25] \times [-20,20] \times [0,20]$, which yields a maximum blockage ratio of $0.8\,\%$ for the wing with $\lambda = 0.27$, $sAR = 2$ and $\alpha = 22^\circ$. The computational set-up is shown in figure 1(bd).

We have prescribed a Dirichlet boundary condition of $(\rho, u_x, u_y, u_z, p) = (\rho _\infty, U_\infty, 0, 0, p_\infty )$ at the inlet and far-field boundaries, where $\rho$ is density, $p$ is pressure and $u_x$, $u_y$ and $u_z$ are velocity components in $x$, $y$ and $z$ directions, respectively. The subscript $\infty$ denotes the free-stream values. A symmetry boundary condition is prescribed along the root plane, $z/c = 0$. We have evaluated the applicability of the root-symmetry boundary condition by conducting direct numerical simulations of flows over full wing configurations, without root symmetry, for wings at $\alpha = 22^\circ$ and $\lambda = 0.27$ and $1$. For both wings, we note that the wake exhibits root-concentrated vortex shedding and remains symmetric with respect to the wing root over large computational times.

A no-slip adiabatic boundary condition is set on the airfoil surface. For vortical structures to convect out of the domain, a sponge layer is applied over $x/L_c \in [15,25]$ with the target state being the running time-averaged state over $5$ convective time units (Freund Reference Freund1997). Simulations start from uniform flow and are performed with a constant acoustic Courant–Friedrichs–Lewy (CFL) number of $1$ until transients are washed out of the computational domain. The time to flush out the transients varies depending on the wing planform and angle of attack, generally ranging from $50$ to $300$ convective time units. After the transients are washed out of the domain, flows are simulated with a constant time step defined such that the CFL number is smaller than one. Flow statistics are collected for $100$ to $300$ convective time units, depending on the flow-field characteristics and spectral content to ensure convergence. A detailed discussion on verification is provided in Appendix A.

3. Results

3.1. Overview of tapered wing wakes

In figure 2, we present instantaneous post-stall flows over tapered wings, which exhibit a rich diversity of wake structures through the combined effects of LE and TE sweep. Taper effects on laminar separated flows are entwined with the effects of LE and TE sweep angles. However, by studying straight tapered wings, that is, wings with $\varLambda _{c/2}$ and $\varLambda _{c/4}$ approximately zero, we can distinguish the effects of taper from other geometrical parameters.

Figure 2. Instantaneous flows around tapered wings with $sAR = 2$, $0.27 \le \lambda \le 1$ and $0^\circ \le \varLambda _{LE} \le 30^\circ$, at $\alpha = 18^\circ$ visualized using grey-coloured isosurfaces of $Q = 1$.

For instance, let us explore the flows over wings with $(\lambda,\varLambda _{LE}) = (1,0^\circ )$ and compare them with the wake structures around $(\lambda,\varLambda _{LE}) =(0.27, 10^\circ )$ wings; these flows have $\varLambda _{c/4} = 0^\circ$ and $1.8^\circ$, respectively. For the tapered wing, we note a significant reduction of the tip vortex length caused by the smaller $c_{tip}$. The downstream root shedding, however, exhibits similar hairpin-like structures for both wings. In the near-wake region, for the tapered wing, the spatial spanwise flow fluctuations emerge near the root, over the spanwise vortex on the suction side. Such wake oscillations are absent in the vortical structure that forms over the untapered wing.

We can further explore the distinct taper effects on the wake dynamics by considering wings with $\varLambda _{c/2} \approx 0^\circ$, as shown for the similar flow patterns that develop at the root region for $(\lambda,\varLambda _{LE}) = (1,0^\circ )$ and $(0.27, 18.4^\circ )$ wings. Here, with a lower taper ratio, tip vortices are considerably weakened when compared with the structures near the free end of the untapered wing. The vortical roll structures emerging over the tapered wing appear slanted and aligned with $\varLambda _{LE}$, suggesting that the LE sweep angle plays an important role in defining the behaviour of the near-wake shedding structures.

For tapered wings, the backward-swept LE effect can be observed by fixing $\varLambda _{TE} = 0^\circ$ while the LE is swept backwards with $\varLambda _{LE} = 18.4^\circ$ and $30^\circ$ for $\lambda = 0.5$ and $0.27$, respectively. For such wings, taper shifts the wake-shedding structures closer to the wing tip. An opposite effect is shown in the top row of figure 2, for flows over forward-swept TE wings. These planforms have fixed $\varLambda _{LE} = 0^\circ$, while $\varLambda _{TE} = -18.4^\circ$ and $-30^\circ$ for $\lambda = 0.5$ and $0.27$, respectively. For these cases, we observe that taper reduces the tip vortex length and changes the topology of the root shedding structures. Let us further study the taper effect for highly swept wings, shown in the bottom row of figure 2, with fixed $\varLambda _{LE} = 30^\circ$, while $\varLambda _{TE} = 13.7^\circ$ and $0^\circ$ for $\lambda = 0.5$ and $0.27$, respectively. Here, taper increases the amplitude of wake oscillations. We further detail the discussions of the effects of taper, LE and TE sweep in § 3.3.

The variety of wake structures that appear around tapered wings, as seen in figure 2, calls for a proper characterization of the wake dynamics that associates its behaviour with the wing planform geometry. The above discussions suggest that taper affects the location where unsteadiness emerges and the characteristics of the vortical structures. In the following section, we provide a map that characterizes the wakes of tapered wings.

3.2. Wake classification and aerodynamic forces

We now classify the flow patterns with respect to the wing geometry. Our criterion is based on the examination of the flow characteristics downstream of the airfoil on a $2$-D plane at $x/c = 4$, where we identify the spatial location of maximum time-averaged $\bar {Q}$ and the maximum fluctuating component of $Q^\prime = Q - \bar {Q}$, where $Q$ is the second invariant of the velocity gradient tensor used to identify the vortical structures (Hunt, Wray & Moin Reference Hunt, Wray and Moin1988; Jeong & Hussain Reference Jeong and Hussain1995). Maximum $\bar {Q}$ and $Q^\prime$ located between $0 \le z/(c\ sAR) < 0.5$ are labelled root-dominant, while points with maximum $\bar {Q}$ or $Q^\prime$ between $0.5 \le z/(c\ sAR) \le 1$ are named tip-dominant. We consider the flow as steady when the maximum fluctuating value of $Q^\prime$ is smaller than $0.1$ at $x/c=4$. Using the root and tip locations of $\bar {Q}$ and $Q^\prime$, we classify their wakes into three unsteady and two steady regimes, as shown in figure 3, where the steady–unsteady threshold (black dotted line) is computed via biharmonic spline interpolation. We further verify our classification criterion by carefully inspecting the flow fields. Instantaneous flow fields for all tapered wings shown in figure 3 are provided in Appendix B using isosurfaces of $Q = 1$ coloured by streamwise velocity $u_x$.

Figure 3. Classification of laminar flows over tapered wings into five distinct wake patterns (ae) shown for $sAR = 2$ wings visualized with time-averaged $\bar {Q} = 1$ in grey and instantaneous $Q^\prime = 0.2$ coloured by $u_x^\prime$. Classification map for (fh) $sAR = 1$ and (ik) $sAR = 2$ coloured by $\overline {C_L/C_D}$. Black dashed lines mark transition from steady to unsteady flows. (a) $(\lambda, \varLambda_{LE}, \alpha) = (0.27, 20^{\circ}, 22^{\circ})$, (b) $(\lambda, \varLambda_{LE}, \alpha) = (0.27, 50^{\circ}, 22^{\circ})$, (c) $(\lambda, \varLambda_{LE}, \alpha) = (1, 10^{\circ}, 22^{\circ})$, (d) $(\lambda, \varLambda_{LE}, \alpha) = (0.5, 50^{\circ}, 22^{\circ})$, (e) $(\lambda, \varLambda_{LE}, \alpha) = (1, 30^{\circ}, 14^{\circ})$, (f) $sAR = 1, \alpha = 14^{\circ}$, (g) $sAR = 1, \alpha = 18^{\circ}$, (h) $sAR = 1, \alpha = 22^{\circ}$, (i) $sAR = 2, \alpha = 14^{\circ}$, (j) $sAR = 2, \alpha = 18^{\circ}$ and (k) $sAR = 2, \alpha = 22^{\circ}$.

The first flow regime (${\boldsymbol {\triangle }}$) is composed of tapered wing wakes that have both maximum $\bar {Q}$ and $Q^\prime$ found near the root region. Such wakes appear for tapered wings with low LE sweep angles. For such wings, the tip vortex tends to be short in length and the taper and forward-swept TE effects concentrate shedding at the wing root, as shown in figure 3(a) for $(\lambda,\varLambda _{LE}) = (0.27, 20^\circ )$ at $\alpha = 22^\circ$. The second flow regime of unsteady wakes (${\boldsymbol {\diamond }}$) occurs when both maximum $\bar {Q}$ and $Q^\prime$ are found over the tip region. Such wakes are observed around tapered wings over a broad range of $\lambda$ values, being present for wings with high LE sweep angles. The flow over such wings often exhibits hairpin-like vortices downstream in the wake aligned with the wing tip, as shown in figure 3(b) for $(\lambda,\varLambda _{LE}) = (0.27, 50^\circ )$ at $\alpha = 22^\circ$.

The third flow regime of unsteady wakes (${\boldsymbol {\triangleleft }}$) around tapered wings presents maximum $\bar {Q}$ at the wing tip with maximum $Q^\prime$ at the root. This wake characteristic is often present for slightly tapered and swept wings, that is, wings with high $\lambda$ and low LE sweep angles. Such wings exhibit a distinct tip vortex formation, at the location of the maximum $\bar {Q}$, and wake shedding near the root. On some occasions, the tip vortex exhibits weak unsteady flow oscillations, as shown in figure 3(c) for $(\lambda,\varLambda _{LE}) = (1, 10^\circ )$ at $\alpha = 22^\circ$, while the most energetic vortices are generally observed over the root region.

There are two distinct flow regimes of steady wakes shown herein, as seen in figure 3(d,e). The first one (${\boldsymbol {\triangledown }}$) is comprised of wakes with a steady streamwise vortex that develops into the wake. Such flows are mainly exhibited around highly swept $sAR = 2$ wings with high and moderate taper ratios, $\lambda \ge 0.5$, as shown in figure 3(d) for $(\lambda,\varLambda _{LE}) = (0.5, 50^\circ )$ at $\alpha = 22^\circ$. The second steady wakes regime (${\boldsymbol {\square }}$) is comprised of flows with no significant wake structures, with maximum $\bar {Q} \le 0.1$ in the wake and are commonly observed for $sAR = 1$ wings, as shown in figure 3(e) for a $(\lambda,\varLambda _{LE}) = (1, 30^\circ )$ wing at $\alpha = 14^\circ$, and for $sAR = 2$ wings at lower angles of attack, low taper ratios and high LE sweep angles.

In figure 3(fk), we present the classification for all wings studied herein. For $sAR = 1$ wings, whose classification is shown in figure 3(fh), there are fewer changes in wake class, when compared with $sAR = 2$ wings. For the higher-aspect-ratio wings, with a fixed $\varLambda _{LE}$, we often notice two or three distinct classes of wake behaviour as $\lambda$ changes. On the other hand, for $sAR = 1$ wings, the same class for all $\lambda$ is observed regularly. Taper effects become increasingly important for $sAR = 2$ wings to alter their wake characteristics, as shown in figure 3(ik), not only affecting the steady–unsteady wake behaviour, but also producing distinct wakes as a function of $\lambda$. For $sAR = 2$ wings with high LE sweep angles, the transition from steady to unsteady wakes is dependent on $\lambda$. Generally, untapered wings with high LE sweep angle wakes remain steady, while unsteady flow structures emerge in the wakes of tapered swept wings.

The combination LE sweep angle increase and taper ratio decrease is shown to promote wake unsteadiness. In addition, tapered wings with high LE sweep angles exhibit enhanced aerodynamic performance when compared with untapered and unswept wings, as shown in figure 3(fk). To visualize this trend, each symbol associated with a wake class is coloured by the time-averaged lift-to-drag ratio, $\overline {C_L/C_D}$, showing that the higher lift-to-drag ratio coefficients, for all $sAR$ and $\alpha$ combination, appear for the wings with lower $\lambda$ and higher backward-swept LE. Here, the aerodynamic forces are reported with lift and drag coefficients defined as

(3.1a,b) \begin{equation} C_L = \frac{F_y}{\dfrac{1}{2}\rho U_\infty^2 b c} \quad \text{and} \quad C_D = \frac{F_x}{\dfrac{1}{2}\rho U_\infty^2 b c} ,\end{equation}

where $F_x$ and $F_y$ are the $x$ and $y$ components of the force on the wing, respectively. Furthermore, we study the aerodynamic loads through the time-averaged $\overline {C_L}$, $\overline {C_D}$ and $\overline {C_L/C_D}$ for selected $sAR = 2$ wings, as shown in figure 4, to reveal the influence of taper and sweep on the aerodynamic forces. The blue symbols present the aerodynamic loads for tapered wings with unswept LE and forward-swept TE. The red symbols show the results for tapered wings with backward-swept LE and unswept TE, while the yellow symbols represent tapered wings with $\varLambda _{LE} = 40^\circ$.

Figure 4. Time-averaged lift, drag and lift-to-drag coefficients, $\overline {C_L}$, $\overline {C_D}$ and $\overline {C_D}$, respectively, for $sAR = 2$ tapered wings with $0.27 \le \lambda \le 1$ at (ac) $14^\circ \le \alpha \le 22^\circ$. Blue downward-pointing triangles, wings with unswept LE and forward-swept TE; red upward-pointing triangles, wings with backward-swept LE and unswept TE; yellow diamonds, wings with $\varLambda _{LE} = 40^\circ$. (a) $\alpha = 14^\circ$, (b) $\alpha = 18^\circ$ and (c) $\alpha = 22^\circ$.

The effects of wing taper on the wakes and aerodynamic forces are strongly dependent on the combination of taper and sweep angle. Let us start from the untapered and unswept wings, marked by blue downward-pointing triangles at $\lambda = 1$. The flow fields around these wings are characterized by root shedding and a strong tip vortex, as seen in figure 2. While keeping the LE unswept, the TE is forward-swept for low $\lambda$. In the wake, such taper produces a concentration of both steady and fluctuating wake structures near the root, as shown in figure 3(ik). A root-concentrated wake with a small tip vortex substantially decreases $\overline {C_L}$, $\overline {C_D}$ and $\overline {C_L/C_D}$ with $\lambda$. It is noteworthy that at the same $\lambda$, tapered wings having backward-swept LE and unswept TE exhibit a higher $\overline {C_L}$ and $\overline {C_L/C_D}$, as shown in figure 4. For such wings, we recall that taper shifts wake structures towards the tip region, as shown in figure 2.

Backward-swept LE enhances the aerodynamic efficiency of tapered wings in post-stall laminar flow conditions. This LE-sweep-induced improvement in aerodynamic loads also occurs for tapered wings with high LE sweep angles. For instance, untapered swept wings present lower values of $\overline {C_L}$ for all angles of attack, as shown by the diamond-shaped yellow symbols for $\lambda = 1$ in figure 4, while significantly reducing wake oscillations, as shown in figures 2 and 3(df).

For $sAR = 2$ wings with high LE sweep at high incidence $\alpha \ge 18^\circ$, taper causes a change in the wake regime, as shown in figure 3(ik). For such wings, wake shedding emerges near the wing tip. The change in wake flow regime for tapered wings at high incidence causes an increase in $\overline {C_L/C_D}$. At lower incidence, $\alpha = 14^\circ$, the wake regime remains steady without noticeable vortices for tapered wings and the aerodynamic forces remain fairly constant for all $\lambda$, as shown in figure 4(a). Let us now describe in detail the taper effects on the wake characteristics in the following section and relate the wakes to the lift distribution over the wing.

3.3. Wake characteristics

3.3.1. Tapered wings with unswept LE and forward-swept TE

Let us take a closer look at the effect of wing taper for unswept LE wings with forward-swept TE, as it allows us to isolate the $\varLambda _{TE}$ effect on the wake dynamics. For tapered wings with $\lambda = 0.27, 0.5, 0.7$ and $1$, the planforms we study in this section have $\varLambda _{TE} = -30^\circ, -18.4^\circ, -10^\circ$ and $0^\circ$, respectively. The negative $\varLambda _{TE}$ indicates forward sweep. The LE is fixed with $\varLambda _{LE} = 0^\circ$. For such wings, taper has a negative impact on the aerodynamic performance, while concentrating the unsteady shedding to a narrow region near the root, and significantly reducing the tip vortex strength, as shown in figure 5(a).

Figure 5. (a) Isosurfaces of flow fields around tapered wings with $sAR = 2$, $\varLambda _{LE} = 0^\circ$, $\lambda = 0.27$ and $1$, $\alpha = 14^\circ$ and $22^\circ$. Time-averaged $\bar {Q} = 1$ isosurface is shown in grey. Instantaneous $Q^\prime = 0.2$ isosurface is shown coloured by $u_x^\prime$. (b) Spanwise distribution of $\| \boldsymbol {\textit {u}}^\prime \|_2$ for different $\lambda$ for $\varLambda _{LE} = 0^\circ$ wings. (c) Spatial–temporal (top) and PSD (bottom) of ${u}_y$ distribution over the spanwise direction from probes located at $(x,y)/c = (3,-0.5)$ for the $\lambda = 0.27$ and $1$ tapered wings at $\alpha = 22^\circ$ shown above. (d) Sectional lift distribution over wingspan for tapered wings at $\alpha = 22^\circ$.

Tapered wings have a smaller $c_{tip}$, which weakens the tip vortices and decreases their length, alleviating the inboard downwash over the wing. Such tip vortex attenuation and the aforementioned concentration of shedding over the root region occur for all angles of attack shown herein. The influence of the incidence angle appears on the formation of secondary vortices near the wing tip. For wings at high incidence angle, a secondary tip vortex is known to emerge from the LE, as shown in figure 6 (DeVoria & Mohseni Reference DeVoria and Mohseni2017; Zhang et al. Reference Zhang, Hayostek, Amitay, He, Theofilis and Taira2020b). For the tapered wings with forward-swept TE at $\alpha = 22^\circ$, there is also another core vortex that emerges near the wing tip from the TE. This structure is seen over the vortex sheet rolling up the TE as a slanted vortex pointing towards the root for the lower taper ratio, as visualized in figure 6.

Figure 6. Vortical structures emerging near the wing tip around tapered wings with unswept LE and forward-swept TE. Top view of time-averaged flow fields around $sAR = 2$ tapered wings at $\alpha = 14^\circ$ and $22^\circ$, with $\varLambda _{LE} = 0^\circ$, visualized using grey-coloured time-averaged isosurfaces of $\bar {Q} = 2$.

To gain further insights into the characteristics of wake vortices, we study the unsteady flow behaviour over the wingspan using probe measurements of velocity fluctuations over $x/c \in [3,4]$, $y/c \in [-1.5,0.5]$. The $x/c$ location is arbitrary and does not affect significantly the results. The $y/c$ range encompasses the region where vortical structures appear. Over this region, we probe the norm of the root mean square of the velocity, $\|\boldsymbol {\textit {u}}^\prime \|_2$. This measurement represents the spanwise distribution of flow unsteadiness, as shown in figure 5(b).

By examining at the spanwise $\|\boldsymbol {\textit {u}}^\prime \|_2$ distribution in figure 5(b) for untapered wings (blue), we notice that the flow unsteadiness peaks at $z/c \approx 0.5$ and decays towards the wing tip for both angles of attack. For tapered wings, the spanwise $\|\boldsymbol {\textit {u}}^\prime \|_2$ curves are independent of the taper ratio for $\lambda \le 0.7$. For such wings, taper yields an attenuation of the $\|\boldsymbol {\textit {u}}^\prime \|_2$ peak. The peak of $\|\boldsymbol {\textit {u}}^\prime \|_2$ also moves towards $z/c \approx 0$, showing a concentration of unsteadiness towards the wing root for tapered wings.

Next, we analyse the spatial–temporal distribution of $u_y$ from probes located at $(x,y)/c = (3,-0.5)$ over the spanwise direction, to investigate how wing taper affects the shedding behaviour. Herein, temporal frequency is characterized through the Strouhal number defined as $St = f (c \sin \alpha / U_\infty )$, where $f$ is the frequency. For comparison, the wake spectrum for the flow over an untapered wing is shown on the left of figure 5(c). For this wing, there is a narrow peak of oscillations at $St \approx 0.14$. The wake spectrum is clean with a vortex-shedding pattern comprised of spanwise-dominated vorticity near the root, forming hairpin vortices and a steady streamwise vortex at the wing tip. For the tapered wing, the spectrum is broadband as a result of the mixing of streamwise and spanwise vortices near the wing root. Even though the wake exhibits more mixing, the spanwise structures remain dominant, being related to the power spectral density (PSD) peak at $St \approx 0.13$. We note that the PSD peak occurs at a lower $St$ than that observed for the untapered wing, as the core unsteady structures that populate the downstream wake arise from the root region of the wing, where the chord length is large.

As the post-stall wakes around tapered wings with forward-swept TE concentrate shedding near the root, they also alter the sectional load distribution as the near-wake vortices play an important role in generating lift and drag over the wing. While untapered wings exhibit a peak in sectional lift near the tip region, we note that the root contribution to lift is higher for tapered wings, as shown in figure 5(d). Such load distribution is generally positive for flight stability (Anderson Reference Anderson2010). For laminar post-stall flows over wings, the emergence of near-wake vortices closer to the wing surface can provide added lift (Lee et al. Reference Lee, Hsieh, Chang and Chu2012; Zhang & Taira Reference Zhang and Taira2022). As shown in figure 5(a), there are fewer large near-wake structures over the wing for the lower $\lambda$. This is a possible reason of the decrease in $C_l$ over the entire wingspan experienced by the $\lambda = 0.27$ wing.

3.3.2. Tapered wings with backward-swept LE and unswept TE

Next, let us analyse the taper effects of wings with backward-swept LE and fixed unswept TE, to understand and separate the effects of $\varLambda _{LE}$ on the global wake. For such wings with $\lambda = 0.27, 0.5, 0.7$ and $1$, the planforms have $\varLambda _{LE} = 30^\circ, 18.4^\circ, 10^\circ$ and $0^\circ$, respectively. The positive $\varLambda _{LE}$ indicates backward sweep. The TE is fixed with $\varLambda _{TE} = 0^\circ$. For such wings, taper yields an opposite effect on the wake characteristics, when compared with those discussed in § 3.3.1. As shown in figure 4, such wake pattern results in a better aerodynamic performance for wings with the same $\lambda$ but distinct LE and TE sweep angles. Herein, taper shifts the unsteadiness region towards the wing tip, as shown in figure 7(a).

Figure 7. (a) Isosurfaces of flow fields around tapered wings with $sAR = 2$, $\varLambda _{TE} = 0^\circ$, $\lambda = 0.27$ and $0.7$, $\alpha = 14^\circ$ and $22^\circ$. Time-averaged $\bar {Q} = 1$ isosurface is shown in grey. Instantaneous $Q^\prime = 0.2$ isosurface is shown coloured by $u_x^\prime$. (b) Spanwise distribution of $\| \boldsymbol {\textit {u}}^\prime \|_2$ for different $\lambda$ for $\varLambda _{LE} = 0^\circ$ wings. (c) Spatial–temporal (top) and PSD (bottom) of ${u}_y$ distribution over the spanwise direction from probes located at $(x,y)/c = (3,-0.5)$ for the $\lambda = 0.27$ and $0.7$ tapered wings at $\alpha = 22^\circ$ shown above. (d) Sectional lift distribution over wingspan for tapered wings at $\alpha = 22^\circ$.

Concurrently, the tip vortex weakens for tapered wings with the shortened $c_{tip}$, which alleviates the inboard downwash near the tip, similar to what was observed for the wings in § 3.3.1. This increases the effective angle of attack near the tip and allows for the flow to detach from the wing surface and form wake-shedding structures near $z/c \approx 1$, as shown in figure 7(a). We quantify the effect of wing taper on flow unsteadiness through the winsgpan distribution of $\|\boldsymbol {\textit {u}}^\prime \|_2$, as shown in figure 7(b). For both angles of attack, taper affects the wake-shedding distribution over the wingspan. For $\lambda = 0.27$ (purple), at $\alpha = 22^\circ$, the peak of $\|\boldsymbol {\textit {u}}^\prime \|_2$ appears near the quarter-span at $z/c \approx 1.25$, with a gradual transition towards $z/c \approx 0.5$ from $\lambda =0.27$ to $1$.

As seen in figure 7(b), tapered wings with backward-swept LE and unswept TE exhibit unsteadiness over a larger spanwise length than untapered wings. For instance, let us observe the spanwise $\|\boldsymbol {\textit {u}}^\prime \|_2$ distribution for wings at $\alpha = 22^\circ$. For the untapered wing (blue), $\boldsymbol {\textit {u}}^\prime \ge 0.02$ over $0 \le z/c \le 1$, which is the region where significant unsteady wake structures appear. Now, for the tapered wing with $\lambda = 0.27$, $\boldsymbol {\textit {u}}^\prime \ge 0.02$ over $0 \le z/c \le 1.6$, hence large unsteady structures can be observed over a larger spanwise portion of the wake.

The spatial–temporal distribution of the transverse velocity $u_y$ over the spanwise direction also shows that the wake of backward-swept LE and unswept TE tapered wings exhibits $3$-D vortical structures that result in a broadband wake spectrum, as shown in figure 7(c). The wake, however, is mainly dominated by large quasi-$2$-D spanwise aligned vortex rolls observed for all taper ratios. For $\lambda = 0.27$, as unsteadiness appears over a larger portion of the wingspan, the stronger shedding structures are hairpin-like vortices that appear between $0.5 \le z/c \le 1.5$, as shown on the right of figure 7(c).

The sectional lift distribution, shown in figure 7(d), reveal that tapered wings with backward-swept LE and unswept TE significantly increase the root contribution and reduce the influence of the near-tip region on the overall lift. The increase in root contribution to lift results from the shifting of the separation bubble towards the tip. This shifting causes a pair of near-wake vortices to emerge over the wing surface at the root region, as shown in figure 7(a). Using force element analysis (Chang Reference Chang1992), we reveal the near-wake structures that contribute to lift. In this approach, the volume force elements are identified by the dot product of the Lamb vector $\boldsymbol {u} \times \boldsymbol {\omega }$ and an auxiliary potential $\boldsymbol {\nabla } \phi _i$ (details in Appendix C). Force element analysis shows a vortex pair emerging near the root, which increases the local contribution to the total lift over the wing.

In figure 8, force elements further show that vortical structures with major contribution to lift appear over the separation bubble. Here, this region is illustrated by a black solid line contour at $\overline {u_x} = 0$ on the $2$-D slices at $z/c = 0.1$ and $0.3$. In particular, the emergence of the near-root vortex pair is persistent for wings with backward-swept LE as similar structures have been identified for backward-swept untapered wings by Zhang & Taira (Reference Zhang and Taira2022). These structures are absent for unswept LE wings, both tapered and untapered, as shown in figure 8(b). In fact, force elements over tapered unswept LE wings show that the lift elements emerging over the wing are much smaller than those over tapered wings with backward-swept LE. As the separation bubble near the root becomes larger over tapered wings with unswept LE, the wake structures are shifted far from the wing, reducing their contribution to the total lift.

Figure 8. Lift elements around $\lambda = 0.27$ tapered wings at $\alpha = 22^\circ$ at the time of maximum lift. (a) Backward-swept LE with unswept TE and (b) unswept LE with forward-swept TE . Perspective view shown with isosurfaces of $Q = 1$ coloured by $(\boldsymbol {u} \times \boldsymbol {\omega }) \boldsymbol{\cdot } \boldsymbol {\nabla } \phi _y$. Slices at selected $z/c$ locations with isocontours of lift elements and a black solid line over $\overline {u_x} = 0$.

The tip vortex is also substantially affected by wing taper, becoming smaller than the tip vortex around untapered wings, as shown in figure 9(a). Tip vortices have high importance in terms of the aerodynamic characteristics of the wing (Francis & Kennedy Reference Francis and Kennedy1979; Green & Acosta Reference Green and Acosta1991; Devenport et al. Reference Devenport, Rife, Liapis and Follin1996; Birch et al. Reference Birch, Lee, Mokhtarian and Kafyeke2004; Taira & Colonius Reference Taira and Colonius2009; Dong et al. Reference Dong, Choi and Mao2020; Zhang et al. Reference Zhang, Hayostek, Amitay, He, Theofilis and Taira2020b; Toppings & Yarusevych Reference Toppings and Yarusevych2021Reference Toppings and Yarusevych2022) and, in the case of tapered wings, due to the small $c_{tip}$, tip vortices are attenuated, as a result of the reduced pressure differences between upper and lower sides of the wing near the tip. Beyond that, even for wings with the same $\lambda$, the tip vortex behaviour can be shifted in the $x$ direction depending on how the wing is tapered, whether it has a backward-swept LE or a forward-swept TE, as shown in figure 9(b,c). For this reason, we must also analyse the tip vortex at a distance $x/c$ from the LE at the wing tip, which is identified herein as $x_{tip,LE}$.

Figure 9. Streamwise circulation $|\varGamma _x|$ of tip vortex around (a) an untapered wing and tapered wings with $\lambda = 0.5$ with (b) unswept LE and forward-swept TE and (c) backward-swept LE and unswept LE at $\alpha = 22^\circ$. Flow field visualized with grey-coloured isosurfaces of $\bar {\omega }_x = -2$ and $2$-D slices with isolines of $\bar {\omega }_x$ at specific $(x - x_{tip,LE})/c$ locations. The magnitude of $|\varGamma _x|$ computed for the isocontour of (d) $\overline {\omega _x} = -2$ for tapered swept wings with different planform configurations.

We can observe how wing taper affects the strength of the tip vortex by analysing $\overline {\omega _x}$ near the tip, as shown in figure 9(ac). The isosurfaces of $\overline {\omega _x}$ and the contour lines at representative $(x - x_{tip,LE})/c$ locations show the decay of vorticity magnitude for tapered wings with $\lambda = 0.5$. However, the effect of taper is not the same for both wings, even though they share the same taper ratio. This difference can be quantified as we compute the streamwise circulation $\varGamma = \int _{C} = \boldsymbol {\textit {u}} \boldsymbol{\cdot } \text {d}\boldsymbol {\textit {l}}$. Here, $C$ is the isocontour of $\overline {\omega _x} = -2$, as shown in figure 9(d). The choice of $\overline {\omega _x}$ level is carefully chosen to isolate the tip vortex.

The tip vortex diffuses downstream of the wing, which makes the $|\varGamma _x|$ profiles decay slowly (Edstrand et al. Reference Edstrand, Sun, Schmid, Taira and Cattafesta2018; Zhang et al. Reference Zhang, Hayostek, Amitay, He, Theofilis and Taira2020b). In general, for tapered wings the reduction in $c_{tip}$ is the main cause of the tip vortex weakening; thus the $|\varGamma _x|$ circulation decays with $\lambda$ at any distance from the wing tip. The circulation $|\varGamma _x |$ further reveals how different types of wing taper can affect the strength of the tip vortex, as shown in figure 9(d). At any given streamwise distance from the wing tip at the LE, $(x - x_{tip,LE})/c$, for $\lambda = 0.5$, the streamwise circulation decay is similar for both tapered wings. For $\lambda = 0.27$, the tip vortex strength decays considerably depending on the LE and TE sweep angles. For tapered wings with backward-swept LE and unswept TE, the root shedding is shifted towards the tip region and reduces the tip vortex strength at any distance from the TE, when compared with the wing with same $\lambda$ and forward-swept TE. For forward-swept TE wings, as vortex shedding is concentrated near the wing root, it has a minor influence on the tip vortex.

3.3.3. Tapered wings with high LE sweep angles

The final class of tapered wing shapes we consider are those with high LE sweep. For the swept wings discussed herein, with $\varLambda _{LE} > 30^\circ$, wake oscillations are strongly attenuated. For laminar flows over untapered wings with high sweep angles at moderate angles of attack, the wake becomes steady, while at high angles of attack, unsteadiness may develop in the wing tip region (Zhang et al. Reference Zhang, Hayostek, Amitay, Burstev, Theofilis and Taira2020a; Ribeiro et al. Reference Ribeiro, Yeh and Taira2023b). For highly swept and tapered wings, the flow exhibits wake shedding for small $\lambda$, as shown in figure 10(a).

Figure 10. (a) Isosurfaces of flow fields around tapered swept wings with $sAR = 2$, $\varLambda _{LE} = 40^\circ$, $\lambda = 0.27$ and $0.7$, $\alpha = 18^\circ$ and $22^\circ$. Time-averaged $\bar {Q} = 1$ isosurface is shown in grey. Instantaneous $Q^\prime = 0.2$ isosurface is shown coloured by $u_x^\prime$. (b) Spanwise distribution of $\| \boldsymbol {\textit {u}}^\prime \|_2$ for different $\lambda$ for $\varLambda _{LE} = 0^\circ$ wings. (c) Spatial–temporal (top) and PSD (bottom) of ${u}_y$ distribution over the spanwise direction from probes located at $(x,y)/c = (3,-0.5)$ for the $\lambda = 0.27$ and $0.7$ tapered wings at $\alpha = 22^\circ$ shown above. (d) Sectional lift distribution over wingspan for tapered wings at $\alpha = 22^\circ$.

Here, we analyse wings with a fixed $\varLambda _{LE} = 40^\circ$, while the TE is swept with $\varLambda _{TE} = 10^\circ, 21.6^\circ, 30^\circ$ and $40^\circ$ for $\lambda = 0.27, 0.5, 0.7$ and $1$, respectively. The onset of shedding for highly swept and tapered wings results from the distinct effects of $\varLambda _{LE}$ and $\varLambda _{TE}$. For the present tapered swept wings, the vortical structures emerging from the TE promote unsteadiness in the wake near the wing tip. For lower taper ratios, wings have a low $\varLambda _{TE}$, which causes wake oscillation to appear and become large towards the root. Such effects show that while high $\varLambda _{LE}$ has the effect of stabilizing wake oscillations for untapered wings, the combination of wing taper and sweep can promote wake unsteadiness.

For instance, at $\alpha = 18^\circ$ the wake is steady for $\lambda = 0.7$ with long steady streamwise vortices developing from both LE and TE. At $\lambda = 0.5$, unsteadiness appears with vortex rolls at the wing tip, with wake shedding appearing for $\lambda = 0.27$. We quantify the wing taper effect in figure 10(b). For instance, for the wings with $\lambda \ge 0.7$ at $\alpha = 18^\circ$, the flow is steady and $\|\boldsymbol {\textit {u}}^\prime \|_2$ is negligible in the wake. At $\alpha = 22^\circ$, $\|\boldsymbol {\textit {u}}^\prime \|_2$ is small for untapered wings, increasing considerably in magnitude and spanwise length as the taper ratio decreases. For highly swept tapered wings, the flow fluctuations are exhibited at the tip, further appearing over the midspan for the lower taper ratios.

The unsteady vortices exhibited in the wakes of tapered wings with high LE sweep angles behave as vortex-shedding structures, as shown by the probed $u_y$ in the wake in figure 10(c). For $\lambda = 0.7$, the vortices appear as a consistent flow oscillation departing from the wing tip. For $\lambda = 0.27$, the wake is dominated by spanwise-aligned roll structures that occupy a large portion of the wingspan. As these structures develop from the wingspan region near the wing tip, which has a reduced chord length, their frequency $St \approx 0.15$ is slightly higher than the shedding frequency of untapered wings.

The sectional lift distributions around tapered wings with high LE sweep, seen in figure 7(d), show how wing taper substantially increases the root contribution to the overall lift. Backward-swept wings have a higher contribution of lift from the root region due to the pair of vortical structures near the midspan that attach closer to the wing surface (Zhang & Taira Reference Zhang and Taira2022) similar to the structures shown in figure 8(a). For tapered wings, the combined effect of the midspan vortex pair and the development of shedding structures near the wing results in a considerable increase of the overall lift over the entire wingspan for the $\lambda = 0.27$ wing.

4. Conclusions

We have examined the influence of taper and sweep on the dynamics of wake structures for finite NACA 0015 wings with straight-cut tip at a Reynolds number of $600$ and a Mach number $0.1$. For this study, we performed an extensive campaign of direct numerical simulations of flows over half-span wings with symmetry boundary condition imposed at the wing root. The present numerical study spans over a wide parameter space with angles of attack between $14^\circ \le \alpha \le 22^\circ$, aspect ratios $sAR = 1$ and $2$, LE sweep angles $0^\circ \le \varLambda _{LE} \le 50^\circ$ and taper ratios between $0.27 \le \lambda \le 1$. This parameter space was chosen to characterize the effects of wing taper as well as the LE and TE sweep angles on the wake dynamics.

Through direct numerical simulations, we observe that the flow over unswept and untapered wings forms a strong tip vortex, which interacts with the spanwise vortex detaching from the wing surface at the root region. This flow yields a 3-D and unsteady wake for all angles of attack considered herein. Untapered and swept wings are observed to advect the shedding region towards the wing tip for lower angles of sweep. At higher sweep angles, the wake oscillations are attenuated, yielding a steady wake around wings at lower angles of attack.

Wing taper has a strong influence on the wake dynamics. For tapered wings, the LE and TE are not parallel and have a distinct influence on the flow structures within the stalled flow region. For tapered wings with unswept LE and forward-swept TE, taper concentrates shedding structures towards the wing root and yields a broadband spectral content downstream in the wake as a result of increased mixing in that region. Beyond the unsteady wake shedding, the tip vortex is heavily affected by wing taper, reducing its length considerably for tapered wings, as the chord length decreases towards the tip.

For tapered wings with backward-swept LE and unswept TE, the spanwise length where wake unsteadiness is observed increases as shedding is promoted over a larger portion of the wingspan. For this type of tapered wing planform, in contrast to the forward-swept wing effect, the peak of wake unsteadiness moves towards the wing tip region for lower taper ratios. Moreover, for wings with high LE sweep, although the flow is steady for $\lambda = 1$, taper causes wake unsteadiness to appear. The wake oscillations develop near the wing tip for moderate taper ratios. For low $\lambda$, wings with high LE sweep angles exhibit strong wake-shedding structures occupying a large portion of the wingspan.

Through the detailed analysis of the wake structures, we also provide a map that classifies the wake behaviour of tapered wings associating the behaviour with the wing planform geometry and the angle of attack. The map provides a unique description of the overall flow physics of the wakes around tapered wings and reveals, for each semi aspect ratio and angle of attack, how the steady–unsteady flow behaviour is related to the wing taper and LE sweep angle. The present study shows the effect of taper, as well as the effects of LE and TE sweep and evaluates its impact on the formation of the wake structures.

Lastly, we show how the wing taper affects the aerodynamic forces over the wing. We show that wings with the same taper ratio may present distinct overall lift and aerodynamic performance, as these characteristics are also influenced by the LE sweep of the wing. Our findings show that the combination of wing taper and high LE sweep can considerably improve lift and the aerodynamic performance of the wing in laminar post-stall flows conditions. The present insights gained into the effect of wing taper in the absence of turbulence serve as a stepping stone for future efforts that aim to study, interpret and control higher-Reynolds-number post-stall flows over tapered wings.

Acknowledgements

The first author thanks T.R. Ricciardi and K. Zhang for enlightening discussions. Computational resources were provided by the High Performance Computing Modernization Program at the US Department of Defense and the Texas Advanced Computing Center.

Funding

We acknowledge the support from the US Air Force Office of Scientific Research (programme manager: Dr G. Abate; grant: FA9550-21-1-0174).

Declaration of interest

The authors report no conflict of interest.

Appendix A. Grid verification

We verify the convergence of grid resolution for the numerical results using a wing with $(sAR,\alpha,\varLambda _{LE},\lambda ) = (2,22^\circ,40^\circ,0.27)$. This planform combines a high LE sweep angle and the lowest taper ratio considered in the present study. Herein, we report the aerodynamic forces through their lift coefficients $C_L$. Two meshes are used for verification: a medium and a refined mesh. The medium mesh refinement is the one used throughout the present work. This mesh has $80$ grid points on both pressure and suction sides of the wing and $48$ grid points along the wingspan, with a total of approximately $3.1 \times 10^6$ control volumes. The refined mesh has $120$ grid points on pressure and suction sides, with $64$ grid points along the wingspan, resulting in approximately $4.3 \times 10^6$ control volumes in total. For the refined mesh we have increased the temporal resolution by setting the CFL number to $0.5$. The quality of our medium mesh is assessed through the forces exerted over the wing and the instantaneous vortical elements as shown in figure 11.

Figure 11. Lift coefficient over time and instantaneous flow-field structures visualized using $Q = 1$ coloured by instantaneous velocity component $u_x$ at the lift peak for the two sets of meshes used for grid verification for the wing with $(sAR,\alpha,\varLambda _{LE},\lambda ) = (2,22^\circ,40^\circ,0.27)$.

Appendix B. A portfolio of flow fields around tapered wings

In this appendix, we provide flow-field visualizations of the wake structures around all tapered wings considered in the present study. Flows around $sAR = 1$ wings at $\alpha = 14^\circ$, $18^\circ$ and $22^\circ$ are shown in figures 12, 13 and 14, respectively. Similarly, flows around $sAR = 2$ wings at $\alpha = 14^\circ$, $18^\circ$ and $22^\circ$ are shown in figures 15, 16 and 17, respectively. All flows are visualized using isosurfaces of $Q = 1$, coloured by the streamwise velocity $u_x$.

Figure 12. Instantaneous flow fields around wings of $sAR = 1$ at $\alpha = 14^\circ$ visualized using isosurfaces of $Q = 1$ coloured by streamwise velocity $u_x$.

Figure 13. Instantaneous flow fields around wings of $sAR = 1$ at $\alpha = 18^\circ$ visualized using isosurfaces of $Q = 1$ coloured by streamwise velocity $u_x$.

Figure 14. Instantaneous flow fields around wings of $sAR = 1$ at $\alpha = 22^\circ$ visualized using isosurfaces of $Q = 1$ coloured by streamwise velocity $u_x$.

Figure 15. Instantaneous flow fields around wings of $sAR = 2$ at $\alpha = 14^\circ$ visualized using isosurfaces of $Q = 1$ coloured by streamwise velocity $u_x$.

Figure 16. Instantaneous flow fields around wings of $sAR = 2$ at $\alpha = 18^\circ$ visualized using isosurfaces of $Q = 1$ coloured by streamwise velocity $u_x$.

Figure 17. Instantaneous flow fields around wings of $sAR = 2$ at $\alpha = 22^\circ$ visualized using isosurfaces of $Q = 1$ coloured by streamwise velocity $u_x$.

Appendix C. Force element analysis

Force element analysis (Chang Reference Chang1992) is a method used to identify vortical structures emerging near the wake that play a significant role in exerting aerodynamic loads over the wing. We note that force element theory is similar in spirit to other force element methods derived through a variational approach (Quartapelle & Napolitano Reference Quartapelle and Napolitano1983), as well as to vortex force maps (Li, Zhao & Graham Reference Li, Zhao and Graham2020) and force partition approaches (Menon & Mittal Reference Menon and Mittal2021). The present force element analysis was previously employed to analyse incompressible flows over wings in multiple configurations (Lee et al. Reference Lee, Hsieh, Chang and Chu2012; Zhang et al. Reference Zhang, Hayostek, Amitay, He, Theofilis and Taira2020b; Zhang, Shah & Bilgen Reference Zhang, Shah and Bilgen2022; Zhang & Taira Reference Zhang and Taira2022).

Initially, we define an auxiliary potential with a specific boundary condition of $-\boldsymbol {\textit {n}} \boldsymbol{\cdot } \boldsymbol {\nabla } \phi _i = \boldsymbol {\textit {n}} \boldsymbol{\cdot } \boldsymbol {\textit {e}}_i$ set on the wing surface, where $\phi$ is the auxiliary potential, $\boldsymbol {\textit {n}}$ is the unit wall normal vector and $\boldsymbol {\textit {e}}_i$ is the unit vector in the $i$th direction. The inner product of the Navier–Stokes equations with $\boldsymbol {\nabla } \phi$ and performing an integral over the fluid domain, the forces exerted in the $i$th direction may be expressed as

(C1)\begin{equation} F_i = \int_V \boldsymbol{\omega} \times \boldsymbol{\textit{u}} \boldsymbol{\cdot } \boldsymbol{\nabla} \phi_i \,\text{d}V + \frac{1}{Re} \int_S \boldsymbol{\omega} \times \boldsymbol{\textit{n}} \boldsymbol{\cdot } (\boldsymbol{\nabla} \phi_i + \boldsymbol{\textit{e}}_\boldsymbol{\textit{i}}) \,\text{d}S,\end{equation}

where the first integral represents the volume force elements and the second integral term is comprised of the surface force elements. For $Re_c = 600$ flows, the volume elements are responsible for the major contribution to the total force over the wing. We note that the auxiliary potential velocity field $\phi _i$ decays rapidly far from the surface. For this reason, the structures that have a higher contribution to lift are the ones that emerge near the surface.

To visualize the vortical structures associated with lift generation, one can take the Hadamard product of $\boldsymbol {\nabla } \phi _i$ and the Lamb vector ($\boldsymbol {\omega } \times \boldsymbol {\textit {u}}$), as shown in figure 8. The resulting $(\boldsymbol {u} \times \boldsymbol {\omega }) \boldsymbol{\cdot } \boldsymbol {\nabla } \phi _y$ variable is often called the lift element. Lastly, we recall that the force element theory used herein considers incompressible Navier–Stokes equations. Nevertheless, important insights into the flow structures may still be obtained for weakly compressible flows, such as the ones considered herein and in previous studies as well (Ribeiro et al. Reference Ribeiro, Yeh, Zhang and Taira2022).

Footnotes

Present address: Faculdade de Engenharia Mecânica, Universidade Estadual de Campinas, Campinas, SP 13083-860, Brazil.

Present address: Department of Aerospace Engineering and Engineering Mechanics, The University of Texas at Austin, Austin, TX 78712, USA.

§

Present address: Faculty of Aerospace Engineering, Technion - Israel Institute of Technology, Haifa 32000, Israel.

References

Abbott, I.H. & Von Doenhoff, A.E. 1959 Theory of Wing Sections: Including a Summary of Airfoil Data. Dover.Google Scholar
Ananda, G.K., Sukumar, P.P. & Selig, M.S. 2015 Measured aerodynamic characteristics of wings at low Reynolds numbers. Aerosp. Sci. Technol. 42, 392406.CrossRefGoogle Scholar
Anderson, R.F. 1936 Determination of the characteristics of tapered wings. Tech. Rep. NACA-TR-572. NASA Langley.Google Scholar
Anderson, J.D. 1999 Aircraft Performance and Design. WCB/McGraw-Hill.Google Scholar
Anderson, J.D. 2010 Fundamentals of Aerodynamics. McGraw-Hill.Google Scholar
Bippes, H. & Turk, M. 1980 Windkanalmessungen in einem Rechteckflügel bei anliegender und abgelöster Strömung. Tech. Rep. 251-80 A 18. DFVLR Forschungsbericht IB.Google Scholar
Birch, D., Lee, T., Mokhtarian, F. & Kafyeke, F. 2004 Structure and induced drag of a tip vortex. J. Aircraft 41 (5), 11381145.CrossRefGoogle Scholar
Black, J. 1956 Flow studies of the leading edge stall on a swept-back wing at high incidence. Aeronaut. J. 60 (541), 5160.CrossRefGoogle Scholar
Braza, M., Faghani, D. & Persillon, H. 2001 Successive stages and the role of natural vortex dislocations in three-dimensional wake transition. J. Fluid Mech. 439, 141.CrossRefGoogle Scholar
Breitsamter, C. & Laschka, B. 2001 Vortical flowfield structure at forward swept-wing configurations. J. Aircraft 38 (2), 193207.CrossRefGoogle Scholar
Brès, G.A., Ham, F.E., Nichols, J.W. & Lele, S.K. 2017 Unstructured large-eddy simulations of supersonic jets. AIAA J. 55 (4), 11641184.CrossRefGoogle Scholar
Buchholz, J.H.J. & Smits, A.J. 2006 On the evolution of the wake structure produced by a low-aspect-ratio pitching panel. J. Fluid Mech. 546, 433443.CrossRefGoogle Scholar
Burtsev, A., He, W., Hayostek, S., Zhang, K., Theofilis, V., Taira, K. & Amitay, M. 2022 Linear modal instabilities around post-stall swept finite wings at low Reynolds numbers. J. Fluid Mech. 944, A6.CrossRefGoogle Scholar
Burtsev, A., Theofilis, V., Ribeiro, J.H.M., Taira, K., Neal, J.M. & Amitay, M. 2023 Wake dynamics of tapered wings. Part III: triglobal linear stability analysis. AIAA Paper 2023-2299.CrossRefGoogle Scholar
Chang, C.-C. 1992 Potential flow and forces for incompressible viscous flow. Proc. R. Soc. Lond. A 437 (1901), 517525.Google Scholar
Dallmann, U.C. 1988 Three-dimensional vortex structures and vorticity topology. Fluid Dyn. Res. 3 (1–4), 183189.CrossRefGoogle Scholar
Délery, J.M. 2001 Robert Legendre and Henri Werlé: toward the elucidation of three-dimensional separation. Annu. Rev. Fluid Mech. 33, 128.CrossRefGoogle Scholar
Devenport, W.J., Rife, M.C., Liapis, S.I. & Follin, G.J. 1996 The structure and development of a wing-tip vortex. J. Fluid Mech. 312, 67106.CrossRefGoogle Scholar
DeVoria, A.C. & Mohseni, K. 2017 On the mechanism of high-incidence lift generation for steadily translating low-aspect-ratio wings. J. Fluid Mech. 813, 110126.CrossRefGoogle Scholar
Dong, L., Choi, K.-S. & Mao, X. 2020 Interplay of the leading-edge vortex and the tip vortex of a low-aspect-ratio thin wing. Exp. Fluids 61 (9), 115.CrossRefGoogle Scholar
Durante, D., Rossi, E. & Colagrossi, A. 2020 Bifurcations and chaos transition of the flow over an airfoil at low Reynolds number varying the angle of attack. Commun. Nonlinear Sci. Numer. Simul. 89, 105285.CrossRefGoogle Scholar
Edstrand, A.M., Sun, Y., Schmid, P.J., Taira, K. & Cattafesta, L.N. 2018 Active attenuation of a trailing vortex inspired by a parabolized stability analysis. J. Fluid Mech. 855, R2.CrossRefGoogle Scholar
Eldredge, J.D. & Jones, A.R. 2019 Leading-edge vortices: mechanics and modeling. Annu. Rev. Fluid Mech. 51, 75104.CrossRefGoogle Scholar
Falkner, V.M. 1950 Sweepback and wing taper: the effect of sweepback on the lift, drag, and aerodynamic centre of a tapered wing. Aircraft Engng Aerosp. Technol. 22 (10), 296300.CrossRefGoogle Scholar
Francis, M.S. & Kennedy, D.A. 1979 Formation of a trailing vortex. J. Aircraft 16 (3), 148154.CrossRefGoogle Scholar
Freund, J.B. 1997 Proposed inflow/outflow boundary condition for direct computation of aerodynamic sound. AIAA J. 35 (4), 740742.CrossRefGoogle Scholar
Freymuth, P., Finaish, F. & Bank, W. 1987 Further visualization of combined wing tip and starting vortex systems. AIAA J. 25 (9), 11531159.CrossRefGoogle Scholar
Gaster, M. 1967 The structure and behaviour of laminar separation bubbles. AGARD CP-4.Google Scholar
Green, S.I. & Acosta, A.J. 1991 Unsteady flow in trailing vortices. J. Fluid Mech. 227, 107134.CrossRefGoogle Scholar
Gursul, I., Gordnier, R. & Visbal, M. 2005 Unsteady aerodynamics of nonslender delta wings. Prog. Aerosp. Sci. 41 (7), 515557.CrossRefGoogle Scholar
Gursul, I. & Wang, Z. 2018 Flow control of tip/edge vortices. AIAA J. 56 (5), 17311749.CrossRefGoogle Scholar
Harper, C.W. & Maki, R.L. 1964 A review of the stall characteristics of swept wings. Tech. Rep. NASA/TN D-2373. NASA, Washington, DC.Google Scholar
He, W., Gioria, R.S., Pérez, J.M. & Theofilis, V. 2017 a Linear instability of low Reynolds number massively separated flow around three NACA airfoils. J. Fluid Mech. 811, 701741.CrossRefGoogle Scholar
He, W., Tendero, J.A., Paredes, P. & Theofilis, V. 2017 b Linear instability in the wake of an elliptic wing. Theor. Comput. Fluid Dyn. 31, 483504.CrossRefGoogle Scholar
Hoarau, Y., Braza, M., Ventikos, Y., Faghani, D. & Tzabiras, G. 2003 Organized modes and the three-dimensional transition to turbulence in the incompressible flow around a NACA 0012 wing. J. Fluid Mech. 496, 6372.CrossRefGoogle Scholar
Hornung, H.G. & Perry, A.E. 1984 Some aspects of three-dimensional separation. Part I. Streamsurface bifurcations. Z. Flugwiss. Weltraumforsch. 8, 7787.Google Scholar
Horton, H.P. 1968 Laminar separation bubbles in two and three dimensional incompressible flow. PhD thesis, Queen Mary University of London.Google Scholar
Huang, Y., Venning, J., Thompson, M.C. & Sheridan, J. 2015 Vortex separation and interaction in the wake of inclined trapezoidal plates. J. Fluid Mech. 771, 341369.CrossRefGoogle Scholar
Huang, R.F., Wu, J.Y., Jeng, J.H. & Chen, R.C. 2001 Surface flow and vortex shedding of an impulsively started wing. J. Fluid Mech. 441, 265292.CrossRefGoogle Scholar
Hunt, J.C.R., Abel, C.J., Peterka, J.A. & Woo, H. 1978 Kinematic studies of the flows around free or surface-mounted obstacles; applying topology to flow visualization. J. Fluid Mech. 86 (1), 179200.CrossRefGoogle Scholar
Hunt, J.C.R., Wray, A.A. & Moin, P. 1988 Eddies, streams, and convergence zones in turbulent flows. In Proceedings of the Summer Program, Center for Turbulence Research, Stanford, CA, pp. 193–208.Google Scholar
Irving, H.B. 1937 Some notes on tapered wings: the effects of brake flaps, taper and ‘sweep’ on stalling. Aircraft Engng Aerosp. Technol. 9 (2), 3136.CrossRefGoogle Scholar
Jeong, J. & Hussain, F. 1995 On the identification of a vortex. J. Fluid Mech. 285, 6994.CrossRefGoogle Scholar
Khalighi, Y., Ham, F., Nichols, J., Lele, S.K. & Moin, P. 2011 Unstructured large eddy simulation for prediction of noise issued from turbulent jets in various configurations. AIAA Paper 2011–2886.CrossRefGoogle Scholar
Lee, J.-J., Hsieh, C.-T., Chang, C.-C. & Chu, C.-C. 2012 Vorticity forces on an impulsively started finite plate. J. Fluid Mech. 694, 464492.CrossRefGoogle Scholar
Li, J., Zhao, X. & Graham, M. 2020 Vortex force maps for three-dimensional unsteady flows with application to a delta wing. J. Fluid Mech. 900, A36.CrossRefGoogle Scholar
Lin, J.C.M. & Pauley, L.L. 1996 Low-Reynolds-number separation on an airfoil. AIAA J. 34 (8), 15701577.CrossRefGoogle Scholar
McCormick, B.W. 1995 Aerodynamics, Aeronautics, and Flight Mechanics. John Wiley & Sons.Google Scholar
Menon, K. & Mittal, R. 2021 Quantitative analysis of the kinematics and induced aerodynamic loading of individual vortices in vortex-dominated flows: a computation and data-driven approach. J. Comput. Phys. 443, 110515.CrossRefGoogle Scholar
Millikan, C.B. 1936 On the stalling of highly tapered wings. J. Aerosp. Sci. 3 (5), 145150.Google Scholar
Mueller, T.J. 2001 Fixed and Flapping Wing Aerodynamics for Micro Air Vehicle Applications. American Institute of Aeronautics and Astronautics, Inc.CrossRefGoogle Scholar
Narasimhamurthy, V.D., Andersson, H.I. & Pettersen, B. 2008 Cellular vortex shedding in the wake of a tapered plate. J. Fluid Mech. 617, 355379.CrossRefGoogle Scholar
Navrose, , , Brion, V. & Jacquin, L. 2019 Transient growth in the near wake region of the flow past a finite span wing. J. Fluid Mech. 866, 399430.CrossRefGoogle Scholar
Neal, J.M. & Amitay, M. 2023 Three-dimensional separation over unswept cantilevered wings at a moderate Reynolds number. Phys. Rev. Fluids 8, 014703.CrossRefGoogle Scholar
Neal, J., Burtsev, A., Ribeiro, J.H.M., Taira, K., Theofilis, V. & Amitay, M. 2023 a Similarities in massive separation across Reynolds numbers for swept and tapered finite span wings. Preprint, arXiv:2308.12442.Google Scholar
Neal, J.M., Gares, B., Amitay, M., Burtsev, A., Theofilis, V., Ribeiro, J.H.M. & Taira, K. 2023 b Wake dynamics of tapered wings. Part II: an experimental study. AIAA Paper 2023-2298.CrossRefGoogle Scholar
Pandi, J.S.S. & Mittal, S. 2019 Wake transitions and laminar separation bubble in the flow past an Eppler 61 airfoil. Phys. Fluids 31 (11), 114102.CrossRefGoogle Scholar
Pandi, J.S.S. & Mittal, S. 2023 Streamwise vortices, cellular shedding and force coefficients on finite wing at low Reynolds number. J. Fluid Mech. 958, A10.CrossRefGoogle Scholar
Pauley, L.L., Moin, P. & Reynolds, W.C. 1990 The structure of two-dimensional separation. J. Fluid Mech. 220, 397411.CrossRefGoogle Scholar
Pelletier, A. & Mueller, T.J. 2000 Low Reynolds number aerodynamics of low-aspect-ratio, thin/flat/cambered-plate wings. J. Aircraft 37 (5), 825832.CrossRefGoogle Scholar
Perry, A.E. & Hornung, H.G. 1984 Some aspects of three-dimensional separation. Part II. Vortex skeletons. Z. Flugwiss. Weltraumforsch. 8, 155160.Google Scholar
Piccirillo, P.S. & Van Atta, C.W. 1993 An experimental study of vortex shedding behind linearly tapered cylinders at low Reynolds number. J. Fluid Mech. 246, 163195.CrossRefGoogle Scholar
Prandtl, L. 1920 Theory of lifting surfaces. Tech. Rep. NACA/TN-9. NASA, Washington, DC.Google Scholar
Quartapelle, L. & Napolitano, M. 1983 Force and moment in incompressible flows. AIAA J. 21 (6), 911913.CrossRefGoogle Scholar
Ribeiro, J.H.M., Taira, K., Neal, J.M., Amitay, M., Burtsev, A. & Theofilis, V. 2023 a Wake dynamics of tapered wings. Part I: a computational study. AIAA Paper 2023-2297.CrossRefGoogle Scholar
Ribeiro, J.H.M., Yeh, C.-A. & Taira, K. 2023 b Triglobal resolvent analysis of swept-wing wakes. J. Fluid Mech. 954, A42.CrossRefGoogle Scholar
Ribeiro, J.H.M., Yeh, C.-A., Zhang, K. & Taira, K. 2022 Wing sweep effects on laminar separated flows. J. Fluid Mech. 950, A23.CrossRefGoogle Scholar
Rockwell, D. 1993 Three-dimensional flow structure on delta wings at high angle-of-attack-experimental concepts and issues. AIAA Paper 1993-0550.CrossRefGoogle Scholar
Rossi, E., Colagrossi, A., Oger, G. & Le Touzé, D. 2018 Multiple bifurcations of the flow over stalled airfoils when changing the Reynolds number. J. Fluid Mech. 846, 356391.CrossRefGoogle Scholar
Schewe, G. 2001 Reynolds-number effects in flow around more-or-less bluff bodies. J. Wind Engng Ind. Aerodyn. 8, 12671289.CrossRefGoogle Scholar
Soule, H.A. & Anderson, R.F. 1940 Design charts relating to the stalling of tapered wings. Tech. Rep. NACA-TR-703. NASA.Google Scholar
Taira, K. & Colonius, T. 2009 Three-dimensional flows around low-aspect-ratio flat-plate wings at low Reynolds numbers. J. Fluid Mech. 623, 187207.CrossRefGoogle Scholar
Techet, A.H., Hover, F.S. & Triantafyllou, M.S. 1998 Vortical patterns behind a tapered cylinder oscillating transversely to a uniform flow. J. Fluid Mech. 363, 7996.CrossRefGoogle Scholar
Theofilis, V., Hein, S. & Dallmann, U.C. 2000 On the origins of unsteadiness and three-dimensionality in a laminar separation bubble. Phil. Trans. R. Soc. Lond. A 358, 3229–324.CrossRefGoogle Scholar
Tobak, M. & Peake, D.J. 1982 Topology of three-dimensional separated flows. Annu. Rev. Fluid Mech. 14 (1), 6185.CrossRefGoogle Scholar
Toppings, C.E. & Yarusevych, S. 2021 Structure and dynamics of a laminar separation bubble near a wingtip. J. Fluid Mech. 929, A39.CrossRefGoogle Scholar
Toppings, C.E. & Yarusevych, S. 2022 Structure and dynamics of a laminar separation bubble near a wing root: towards reconstructing the complete lsb topology on a finite wing. J. Fluid Mech. 944, A14.CrossRefGoogle Scholar
Torres, G.E. & Mueller, T.J. 2004 Low-aspect-ratio aerodynamics at low Reynolds numbers. AIAA J. 42 (5), 865873.CrossRefGoogle Scholar
Traub, L.W. 2013 Aerodynamic impact of aspect ratio at low Reynolds number. J. Aircraft 50 (2), 626634.CrossRefGoogle Scholar
Traub, L.W., Botero, E., Waghela, R., Callahan, R. & Watson, A. 2015 Effect of taper ratio at low Reynolds number. J. Aircraft 52 (3), 734747.CrossRefGoogle Scholar
Valles, B., Andersson, H.I. & Jenssen, C.B. 2002 Oblique vortex shedding behind tapered cylinders. J. Fluids Struct. 16 (4), 453463.CrossRefGoogle Scholar
Videler, J.J., Stamhuis, E.J. & Povel, G.D.E. 2004 Leading-edge vortex lifts swifts. Science 306 (5703), 19601962.CrossRefGoogle ScholarPubMed
Visbal, M.R. & Garmann, D.J. 2019 Effect of sweep on dynamic stall of a pitching finite-aspect-ratio wing. AIAA J. 57 (8), 32743289.CrossRefGoogle Scholar
Wei, Z., New, T.H. & Cui, Y.D. 2018 Aerodynamic performance and surface flow structures of leading-edge tubercled tapered swept-back wings. AIAA J. 56 (1), 423431.CrossRefGoogle Scholar
Winkelman, A.E. & Barlow, J.B. 1980 Flowfield model for a rectangular planform wing beyond stall. AIAA J. 18 (8), 10061008.CrossRefGoogle Scholar
Yarusevych, S., Sullivan, P.E. & Kawall, J.G. 2009 On vortex shedding from an airfoil in low-Reynolds-number flows. J. Fluid Mech. 632, 245271.CrossRefGoogle Scholar
Yen, S.-C. & Hsu, C.M. 2007 Flow patterns and wake structure of a swept-back wing. AIAA J. 45 (1), 228236.CrossRefGoogle Scholar
Yen, S.-C. & Huang, L.-C. 2009 Flow patterns and aerodynamic performance of unswept and swept-back wings. J. Fluids Engng 131 (11), 111101.CrossRefGoogle Scholar
Yilmaz, T.O. & Rockwell, D. 2012 Flow structure on finite-span wings due to pitch-up motion. J. Fluid Mech. 691, 518545.CrossRefGoogle Scholar
Zhang, K., Hayostek, S., Amitay, M., Burstev, A., Theofilis, V. & Taira, K. 2020 a Laminar separated flows over finite-aspect-ratio swept wings. J. Fluid Mech. 905, R1.CrossRefGoogle Scholar
Zhang, K., Hayostek, S., Amitay, M., He, W., Theofilis, V. & Taira, K. 2020 b On the formation of three-dimensional separated flows over wings under tip effects. J. Fluid Mech. 895, A9.CrossRefGoogle Scholar
Zhang, K., Shah, B. & Bilgen, O. 2022 Low-Reynolds-number aerodynamic characteristics of airfoils with piezocomposite trailing surfaces. AIAA J. 60 (4), 27012706.CrossRefGoogle Scholar
Zhang, K. & Taira, K. 2022 Laminar vortex dynamics around forward-swept wings. Phys. Rev. Fluids 7 (2), 024704.CrossRefGoogle Scholar
Figure 0

Figure 1. Problem set-up for tapered wings. (a) Geometrical parameters shown in a wing planform with $sAR = b/c = 2$, $\alpha = 18^\circ$, $\lambda = 0.27$ and $\varLambda _{LE} = 18.4^\circ$. (b) For a $(\lambda,\varLambda _{LE}) = (0.5,0^\circ )$ wing, we show the computational domain and (c,d) grids with $2$-D planes at $z/c = 1$ and $y/c = -0.5$, respectively.

Figure 1

Figure 2. Instantaneous flows around tapered wings with $sAR = 2$, $0.27 \le \lambda \le 1$ and $0^\circ \le \varLambda _{LE} \le 30^\circ$, at $\alpha = 18^\circ$ visualized using grey-coloured isosurfaces of $Q = 1$.

Figure 2

Figure 3. Classification of laminar flows over tapered wings into five distinct wake patterns (ae) shown for $sAR = 2$ wings visualized with time-averaged $\bar {Q} = 1$ in grey and instantaneous $Q^\prime = 0.2$ coloured by $u_x^\prime$. Classification map for (fh) $sAR = 1$ and (ik) $sAR = 2$ coloured by $\overline {C_L/C_D}$. Black dashed lines mark transition from steady to unsteady flows. (a) $(\lambda, \varLambda_{LE}, \alpha) = (0.27, 20^{\circ}, 22^{\circ})$, (b) $(\lambda, \varLambda_{LE}, \alpha) = (0.27, 50^{\circ}, 22^{\circ})$, (c) $(\lambda, \varLambda_{LE}, \alpha) = (1, 10^{\circ}, 22^{\circ})$, (d) $(\lambda, \varLambda_{LE}, \alpha) = (0.5, 50^{\circ}, 22^{\circ})$, (e) $(\lambda, \varLambda_{LE}, \alpha) = (1, 30^{\circ}, 14^{\circ})$, (f) $sAR = 1, \alpha = 14^{\circ}$, (g) $sAR = 1, \alpha = 18^{\circ}$, (h) $sAR = 1, \alpha = 22^{\circ}$, (i) $sAR = 2, \alpha = 14^{\circ}$, (j) $sAR = 2, \alpha = 18^{\circ}$ and (k) $sAR = 2, \alpha = 22^{\circ}$.

Figure 3

Figure 4. Time-averaged lift, drag and lift-to-drag coefficients, $\overline {C_L}$, $\overline {C_D}$ and $\overline {C_D}$, respectively, for $sAR = 2$ tapered wings with $0.27 \le \lambda \le 1$ at (ac) $14^\circ \le \alpha \le 22^\circ$. Blue downward-pointing triangles, wings with unswept LE and forward-swept TE; red upward-pointing triangles, wings with backward-swept LE and unswept TE; yellow diamonds, wings with $\varLambda _{LE} = 40^\circ$. (a) $\alpha = 14^\circ$, (b) $\alpha = 18^\circ$ and (c) $\alpha = 22^\circ$.

Figure 4

Figure 5. (a) Isosurfaces of flow fields around tapered wings with $sAR = 2$, $\varLambda _{LE} = 0^\circ$, $\lambda = 0.27$ and $1$, $\alpha = 14^\circ$ and $22^\circ$. Time-averaged $\bar {Q} = 1$ isosurface is shown in grey. Instantaneous $Q^\prime = 0.2$ isosurface is shown coloured by $u_x^\prime$. (b) Spanwise distribution of $\| \boldsymbol {\textit {u}}^\prime \|_2$ for different $\lambda$ for $\varLambda _{LE} = 0^\circ$ wings. (c) Spatial–temporal (top) and PSD (bottom) of ${u}_y$ distribution over the spanwise direction from probes located at $(x,y)/c = (3,-0.5)$ for the $\lambda = 0.27$ and $1$ tapered wings at $\alpha = 22^\circ$ shown above. (d) Sectional lift distribution over wingspan for tapered wings at $\alpha = 22^\circ$.

Figure 5

Figure 6. Vortical structures emerging near the wing tip around tapered wings with unswept LE and forward-swept TE. Top view of time-averaged flow fields around $sAR = 2$ tapered wings at $\alpha = 14^\circ$ and $22^\circ$, with $\varLambda _{LE} = 0^\circ$, visualized using grey-coloured time-averaged isosurfaces of $\bar {Q} = 2$.

Figure 6

Figure 7. (a) Isosurfaces of flow fields around tapered wings with $sAR = 2$, $\varLambda _{TE} = 0^\circ$, $\lambda = 0.27$ and $0.7$, $\alpha = 14^\circ$ and $22^\circ$. Time-averaged $\bar {Q} = 1$ isosurface is shown in grey. Instantaneous $Q^\prime = 0.2$ isosurface is shown coloured by $u_x^\prime$. (b) Spanwise distribution of $\| \boldsymbol {\textit {u}}^\prime \|_2$ for different $\lambda$ for $\varLambda _{LE} = 0^\circ$ wings. (c) Spatial–temporal (top) and PSD (bottom) of ${u}_y$ distribution over the spanwise direction from probes located at $(x,y)/c = (3,-0.5)$ for the $\lambda = 0.27$ and $0.7$ tapered wings at $\alpha = 22^\circ$ shown above. (d) Sectional lift distribution over wingspan for tapered wings at $\alpha = 22^\circ$.

Figure 7

Figure 8. Lift elements around $\lambda = 0.27$ tapered wings at $\alpha = 22^\circ$ at the time of maximum lift. (a) Backward-swept LE with unswept TE and (b) unswept LE with forward-swept TE . Perspective view shown with isosurfaces of $Q = 1$ coloured by $(\boldsymbol {u} \times \boldsymbol {\omega }) \boldsymbol{\cdot } \boldsymbol {\nabla } \phi _y$. Slices at selected $z/c$ locations with isocontours of lift elements and a black solid line over $\overline {u_x} = 0$.

Figure 8

Figure 9. Streamwise circulation $|\varGamma _x|$ of tip vortex around (a) an untapered wing and tapered wings with $\lambda = 0.5$ with (b) unswept LE and forward-swept TE and (c) backward-swept LE and unswept LE at $\alpha = 22^\circ$. Flow field visualized with grey-coloured isosurfaces of $\bar {\omega }_x = -2$ and $2$-D slices with isolines of $\bar {\omega }_x$ at specific $(x - x_{tip,LE})/c$ locations. The magnitude of $|\varGamma _x|$ computed for the isocontour of (d) $\overline {\omega _x} = -2$ for tapered swept wings with different planform configurations.

Figure 9

Figure 10. (a) Isosurfaces of flow fields around tapered swept wings with $sAR = 2$, $\varLambda _{LE} = 40^\circ$, $\lambda = 0.27$ and $0.7$, $\alpha = 18^\circ$ and $22^\circ$. Time-averaged $\bar {Q} = 1$ isosurface is shown in grey. Instantaneous $Q^\prime = 0.2$ isosurface is shown coloured by $u_x^\prime$. (b) Spanwise distribution of $\| \boldsymbol {\textit {u}}^\prime \|_2$ for different $\lambda$ for $\varLambda _{LE} = 0^\circ$ wings. (c) Spatial–temporal (top) and PSD (bottom) of ${u}_y$ distribution over the spanwise direction from probes located at $(x,y)/c = (3,-0.5)$ for the $\lambda = 0.27$ and $0.7$ tapered wings at $\alpha = 22^\circ$ shown above. (d) Sectional lift distribution over wingspan for tapered wings at $\alpha = 22^\circ$.

Figure 10

Figure 11. Lift coefficient over time and instantaneous flow-field structures visualized using $Q = 1$ coloured by instantaneous velocity component $u_x$ at the lift peak for the two sets of meshes used for grid verification for the wing with $(sAR,\alpha,\varLambda _{LE},\lambda ) = (2,22^\circ,40^\circ,0.27)$.

Figure 11

Figure 12. Instantaneous flow fields around wings of $sAR = 1$ at $\alpha = 14^\circ$ visualized using isosurfaces of $Q = 1$ coloured by streamwise velocity $u_x$.

Figure 12

Figure 13. Instantaneous flow fields around wings of $sAR = 1$ at $\alpha = 18^\circ$ visualized using isosurfaces of $Q = 1$ coloured by streamwise velocity $u_x$.

Figure 13

Figure 14. Instantaneous flow fields around wings of $sAR = 1$ at $\alpha = 22^\circ$ visualized using isosurfaces of $Q = 1$ coloured by streamwise velocity $u_x$.

Figure 14

Figure 15. Instantaneous flow fields around wings of $sAR = 2$ at $\alpha = 14^\circ$ visualized using isosurfaces of $Q = 1$ coloured by streamwise velocity $u_x$.

Figure 15

Figure 16. Instantaneous flow fields around wings of $sAR = 2$ at $\alpha = 18^\circ$ visualized using isosurfaces of $Q = 1$ coloured by streamwise velocity $u_x$.

Figure 16

Figure 17. Instantaneous flow fields around wings of $sAR = 2$ at $\alpha = 22^\circ$ visualized using isosurfaces of $Q = 1$ coloured by streamwise velocity $u_x$.