An official website of the United States government
This paper presents a state-variable formulation to model and simulate the 2D unsteady aerodynamics of an airfoil undergoing arbitrary motion kinematics. The model builds upon a large-angle unsteady aerodynamic formulation in which the airfoil is represented using a lumped vortex element (LVE) model. The airfoil is divided into several panels, with a bound vortex placed on each panel. At any time instant, the bound-vortex strengths are determined by employing zero-normal-flow conditions at the control points located on each panel. The vorticity shed from the trailing edge of the airfoil is modeled using discrete vortices that move freely in the flow field. The required state variables are first identified, and all the time derivative terms of the state variables are then derived to form the final state-variable representation. Trailing-edge vortex shedding is incorporated using the Kelvin condition. The final state variable equation can be solved as an ordinary differential equation using any standard ODE-solving algorithm. Three case studies are presented here to evaluate the predictions of the model. In the cases considered here, the airfoil undergoes various unsteady plunge motions. The aerodynamic load history and the wake patterns are compared against the results from the low-order model developed by Narsipur et al. [1] in previous research. The comparison shows that the current state-variable formulation captures the unsteady flow characteristics and the aerodynamic load in good agreement with the reference results.
more »
« less
The vortex-entrainment sheet in an inviscid fluid: theory and separation at a sharp edge
In this paper a model for viscous boundary and shear layers in three dimensions is introduced and termed a vortex-entrainment sheet. The vorticity in the layer is accounted for by a conventional vortex sheet. The mass and momentum in the layer are represented by a two-dimensional surface having its own internal tangential flow. Namely, the sheet has a mass density per-unit-area making it dynamically distinct from the surrounding outer fluid and allowing the sheet to support a pressure jump. The mechanism of entrainment is represented by a discontinuity in the normal component of the velocity across the sheet. The velocity field induced by the vortex-entrainment sheet is given by a generalized Birkhoff–Rott equation with a complex sheet strength. The model was applied to the case of separation at a sharp edge. No supplementary Kutta condition in the form of a singularity removal is required as the flow remains bounded through an appropriate balance of normal momentum with the pressure jump across the sheet. A pressure jump at the edge results in the generation of new vorticity. The shedding angle is dictated by the normal impulse of the intrinsic flow inside the bound sheets as they merge to form the free sheet. When there is zero entrainment everywhere the model reduces to the conventional vortex sheet with no mass. Consequently, the pressure jump must be zero and the shedding angle must be tangential so that the sheet simply convects off the wedge face. Lastly, the vortex-entrainment sheet model is demonstrated on several example problems.
more »
« less
- Award ID(s):
- 1805776
- PAR ID:
- 10097583
- Date Published:
- Journal Name:
- Journal of Fluid Mechanics
- Volume:
- 866
- ISSN:
- 0022-1120
- Page Range / eLocation ID:
- 660 to 688
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
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
-
Recent studies reveal that at large friction Reynolds number delta^+ the outer, inertially-dominated region of the turbulent boundary layer is composed of large scale zones of uniform momentum segregated by narrow fissures of concentrated vorticity. Experiments show that, when scaled by the boundary layer thickness, the fissure thickness is O(1/sqrt{delta^+}), while the dimensional jump in streamwise velocity across each fissure scales in proportion to the friction velocity u_tau. A simple model that exploits these essential elements of the turbulent boundary layer structure at large delta^+ is developed. First, a master wall-normal profile of streamwise velocity is constructed by placing a discrete number of fissures across the boundary layer. The number of fissures and their wall-normal locations follow scalings informed by analysis of the mean momentum equation. The fissures are then randomly displaced in the wall-normal direction, exchanging momentum as they move, to create an instantaneous velocity profile. This process is repeated to generate ensembles of streamwise velocity profiles from which statistical moments are computed. The modelled statistical moments are shown to agree remarkably well with those acquired from direct numerical simulations of turbulent channel flow at large delta^+. In particular, the model robustly reproduces the empirically observed sub-Gaussian behaviour for the skewness and kurtosis profiles over a large range of input parameters.more » « less
-
Motivated by the need for compact descriptions of the evolution of non-classical wakes behind yawed wind turbines, we develop an analytical model to predict the shape of curled wakes. Interest in such modelling arises due to the potential of wake steering as a strategy for mitigating power reduction and unsteady loading of downstream turbines in wind farms. We first estimate the distribution of the shed vorticity at the wake edge due to both yaw offset and rotating blades. By considering the wake edge as an ideally thin vortex sheet, we describe its evolution in time moving with the flow. Vortex sheet equations are solved using a power series expansion method, and an approximate solution for the wake shape is obtained. The vortex sheet time evolution is then mapped into a spatial evolution by using a convection velocity. Apart from the wake shape, the lateral deflection of the wake including ground effects is modelled. Our results show that there exists a universal solution for the shape of curled wakes if suitable dimensionless variables are employed. For the case of turbulent boundary layer inflow, the decay of vortex sheet circulation due to turbulent diffusion is included. Finally, we modify the Gaussian wake model by incorporating the predicted shape and deflection of the curled wake, so that we can calculate the wake profiles behind yawed turbines. Model predictions are validated against large-eddy simulations and laboratory experiments for turbines with various operating conditions.more » « less
-
This work presents an exact solution of Euler's incompressible equations in the context of a bidirectional vortex evolving inside a conically shaped cyclonic chamber. The corresponding helical flowfield is modeled under inviscid conditions assuming constant angular momentum. By leveraging the axisymmetric nature of the problem, a steady-state solution of the generalized Beltramian type is obtained directly from first principles, namely, from the Bragg–Hawthorne equation in spherical coordinates. The resulting stream function representation enables us to fully describe the ensuing swirl-dominated motion including its fundamental flow characteristics. After identifying an isolated singularity that appears at a cone divergence half-angle of 63.43°, two piecewise formulations are provided that correspond to either fluid injection or extraction at the top section of the conical cyclone. In this process, analytical expressions are readily retrieved for the three velocity components, vorticity, and pressure. Other essential flow indicators, such as the theoretically preferred mantle orientation, the empirically favored locus of zero vertical velocity, the maximum polar and axial velocities, the crossflow velocity, and other such terms, are systematically deduced. Results are validated using limiting process verifications and comparisons to both numerical and experimental measurements. The subtle differences between the present model and a strictly Beltramian flowfield are also highlighted and discussed. The conically cyclonic configuration considered here is relevant to propulsive devices, such as vortex-fired liquid rocket engines with tapered walls; meteorological phenomena, such as tornadoes, dust devils, and fire whirls; and industrial contraptions, such as cyclonic flow separators, collectors, centrifuges, boilers, vacuum cleaners, cement grinders, and so on.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1017/jfm.2019.134
