A computational study of vorticity reconnection, associated with the breaking and reconnection of vortex lines, during vortex cutting by a blade is reported. A series of Navier–Stokes simulations of vortex cutting with different values of the vortex strength are described, and the different phases in the vortex cutting process are compared to those of the more traditional vortex tube reconnection process. Each of the three phases of vortex tube reconnection described by Melander & Hussain ( Phys. Fluids A, vol. 1(4), 1989, pp. 633–635) are found to have counterparts in the vortex cutting problem, although we also point out numerous differences in the detailed mechanics by which these phases are achieved. Of particular importance in the vortex cutting process is the presence of vorticity generation from the blade surface within the reconnection region and the presence of strong vortex stretching due to the ambient flow about the blade leading edge. A simple exact Navier–Stokes solution is presented that describes the process by which incident vorticity is stretched and carried towards the surface by the ambient flow, and then interacts with and is eventually annihilated by diffusive interaction with vorticity generated at the surface. The model combines a Hiemenz straining flow, a Burgers vortex sheet and a Stokes first problem boundary layer, resulting in a nonlinear ordinary differential equation and a partial differential equation in two scaled time and distance variables that must be solved numerically. The simple model predictions exhibit qualitative agreement with the full numerical simulation results for vorticity annihilation near the leading-edge stagnation point during vortex cutting.
more »
« less
Transient lift force on a blade during cutting of a vortex with non-zero axial flow
The problem of orthogonal penetration of a blade into the core of a vortex with non-zero axial flow was studied using a combination of scaling theory, a heuristic plug-flow model and full Navier–Stokes simulations. The particular focus of this paper was to understand the mechanics of the transient lift force that occurs during the initial penetration of the blade leading edge into the vortex core, and the relationship of this transient force to the steady-state lift force that develops due to the difference in vortex core radius over the blade surface. The three modelling approaches all lead to the conclusion that the maximum value of the lift coefficient for the transient blade penetration force is proportional to the impact parameter and inversely proportional to the axial flow parameter. This observation is used to develop a simple expression that collapses the predictions of the full Navier–Stokes simulations for lift coefficient over a wide range of parameter values.
more »
« less
- Award ID(s):
- 1144388
- NSF-PAR ID:
- 10343147
- Date Published:
- Journal Name:
- Journal of Fluid Mechanics
- Volume:
- 819
- ISSN:
- 0022-1120
- Page Range / eLocation ID:
- 258 to 284
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
more » « less
In this paper, we investigate the three-dimensional nature of dynamic stall. Conducting the investigation, the flow around a harmonically pitching National Advisory Committee for Aeronautics (NACA) 0012 airfoil is numerically simulated using Unsteady-Reynolds-Averaged Navier–Stokes (URANS) and multiple detached eddy simulation (DES) solvers: the Delayed-DES (DDES) and the Improved-DDES (IDDES). Two- and three-dimensional simulations are performed for each solver, and the results are compared against experimental measurements in the literature. The results showed that three-dimensional simulations surpass two-dimensional ones in capturing the stages of dynamic stall and predicting the lift coefficient values, with a distinguished performance of the DES solvers over the URANS ones. For instance, the IDDES simulations, as an inherently three-dimensional solver, predicted the necessary cascaded amalgamation process of vortices to form the adequate strength of the dynamic stall vortex. This vortex size and timing provided accurate and sufficient suction that resulted in identical matching of the numerical and experimental lift coefficients at the peak value. Hence, the hypothesis that dynamic stall has a three-dimensional nature is supported by the superiority of the three-dimensional simulation in all aspects. In conclusion, it is found that dynamic stall is intrinsically a three-dimensional phenomenon.
-
null (Ed.)Hill's vortex is a classical solution of the incompressible Euler equations which consists of an axisymmetric spherical region of constant vorticity matched to an irrotational external flow. This solution has been shown to be a member of a one-parameter family of steady vortex rings and as such is commonly used as a simple analytic model for a vortex ring. Here, we model the decay of a Hill's vortex in a weakly rotating flow due to the radiation of inertial waves. We derive analytic results for the modification of the vortex structure by rotational effects and the generated wave field using an asymptotic approach where the rotation rate, or inverse Rossby number, is taken to be small. Using this model, we predict the decay of the vortex speed and radius by combining the flux of vortex energy to the wave field with the conservation of peak vorticity. We test our results against numerical simulations of the full axisymmetric Navier–Stokes equations.more » « less
-
We evaluate two leading-edge-based dynamic stall-onset criteria (namely, the maximum magnitudes of the leading-edge suction parameter and the boundary enstrophy flux) for mixed and trailing-edge stall. These criteria have been shown to successfully predict the onset of leading-edge stall at Reynolds numbers of O(10^5), where the leading-edge suction drops abruptly. However, for mixed/trailing-edge stall, leading-edge suction tends to persist even when there is significant trailing-edge reversed flow and stall is underway, necessitating further investigation into the effectiveness of these criteria. Using wall-resolved large-eddy simulations and the unsteady Reynolds-averaged Navier–Stokes method, we simulate one leading-edge stall and three mixed/trailing-edge stall cases at Reynolds numbers of 200,000 and 300,000. We contrast the progression of flow features such as trailing-edge separation and vortex formation across different stall types and evaluate the stall-onset criteria relative to critical points in the flow. We find that the criteria nearly coincide with the instance of leading-edge suction collapse and are reached in advance of dynamic stall vortex formation and lift stall for all four cases. We conclude that the two criteria effectively signal dynamic stall onset in cases where the dynamic stall vortex plays a prominent role.more » « less
-
Lifting line theory describes the cumulative effect of shed vorticity from finite span lifting surfaces. In this work, the theory is reformulated to improve the accuracy of the actuator line model (ALM). This model is a computational tool used to represent lifting surfaces, such as wind-turbine blades in computational fluid dynamics. In ALM, blade segments are represented by means of a Gaussian body force distribution with a prescribed kernel size. Prior analysis has shown that a representation of the blade using an optimal kernel width $\unicode[STIX]{x1D716}^{opt}$ of approximately one quarter of the chord size results in accurate predictions of the velocity field and loads along the blades. Also, simulations have shown that use of the optimal kernel size yields accurate representation of the tip-vortex size and the associated downwash resulting in accurate predictions of the tip losses. In this work, we address the issue of how to represent the effects of finite span wings and tip vortices when using Gaussian body forces with a kernel size larger than the optimal value. This question is relevant in the context of coarse-scale large-eddy simulations that cannot afford the fine resolutions required to resolve the optimal kernel size. For this purpose, we present a filtered lifting line theory for a Gaussian force distribution. Based on the streamwise component of the vorticity transport equation, we develop an analytical model for the induced velocity resulting from the spanwise changes in lift force for an arbitrary kernel scale. The results are used to derive a subfilter-scale velocity model that is used to correct the velocity along the blade when using kernel sizes larger than $\unicode[STIX]{x1D716}^{opt}$ . Tests are performed in large-eddy simulation of flow over fixed wings with constant and elliptic chord distributions using various kernel sizes. Results show that by using the proposed subfilter velocity model, kernel-size independent predictions of lift coefficient and total lift forces agree with those obtained with the optimal kernel size.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1017/jfm.2017.188
