skip to main content


Title: Filtered lifting line theory and application to the actuator line model
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
Award ID(s):
1635430
NSF-PAR ID:
10106934
Author(s) / Creator(s):
;
Date Published:
Journal Name:
Journal of Fluid Mechanics
Volume:
863
ISSN:
0022-1120
Page Range / eLocation ID:
269 to 292
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    The filtered lifting line theory is an analytical approach used to solve the equations of flow subjected to body forces with a Gaussian distribution, such as used in the actuator line model. In the original formulation, the changes in chord length along the blade were assumed to be small. This assumption can lead to errors in the induced velocities predicted by the theory compared to full solutions of the equations. In this work, we revisit the original derivation and provide a more general formulation that can account for significant changes in chord along the blade. The revised formulation can be applied to wings with significant changes in chord along the span, such as wind turbine blades.

     
    more » « less
  2. Laboratory experiments were performed on a geometrically scaled vertical-axis wind turbine model over an unprecedented range of Reynolds numbers, including and exceeding those of the full-scale turbine. The study was performed in the high-pressure environment of the Princeton High Reynolds number Test Facility (HRTF). Utilizing highly compressed air as the working fluid enabled extremely high Reynolds numbers while still maintaining dynamic similarity by matching the tip speed ratio (defined as the ratio of tip velocity to free stream, $\unicode[STIX]{x1D706}=\unicode[STIX]{x1D714}R/U$ ) and Mach number (defined at the turbine tip, $Ma=\unicode[STIX]{x1D714}R/a$ ). Preliminary comparisons are made with measurements from the full-scale field turbine. Peak power for both the field data and experiments resides around $\unicode[STIX]{x1D706}=1$ . In addition, a systematic investigation of trends with Reynolds number was performed in the laboratory, which revealed details about the asymptotic behaviour. It was shown that the parameter that characterizes invariance in the power coefficient was the Reynolds number based on blade chord conditions ( $Re_{c}$ ). The power coefficient reaches its asymptotic value when $Re_{c}>1.5\times 10^{6}$ , which is higher than what the field turbine experiences. The asymptotic power curve is found, which is invariant to further increases in Reynolds number. 
    more » « less
  3. We study the emergence of precessing vortex core (PVC) oscillations in a swirling jet experiment. We vary the swirl intensity while keeping the net mass flow rate fixed using a radial-entry swirler with movable blades upstream of the jet exit. The swirl intensity is quantified in terms of a swirl number $S$ . Time-resolved velocity measurements in a radial–axial plane anchored at the jet exit for various $S$ values are obtained using stereoscopic particle image velocimetry. Spectral proper orthogonal decomposition and spatial cross-spectral analysis reveal the simultaneous emergence of a bubble-type vortex breakdown and a strong helical limit-cycle oscillation in the flow for $S>S_{c}$ where $S_{c}=0.61$ . The oscillation frequency, $f_{PVC}$ , and the square of the flow oscillation amplitudes vary linearly with $S-S_{c}$ . A solution for the coherent unsteady field accurate up to $O(\unicode[STIX]{x1D716}^{3})$ ( $\unicode[STIX]{x1D716}\sim O((S-S_{c})^{1/2})$ ) is determined from the nonlinear Navier–Stokes equations, using the method of multiple scales. We show that onset of bubble type vortex breakdown at $S_{c}$ , results in a marginally stable, helical linear global hydrodynamic mode. This results in the stable limit-cycle precession of the breakdown bubble. The variation of $f_{LC}$ with $S-S_{c}$ is determined from the Stuart–Landau equation associated with the PVC. Reasonable agreement with the corresponding experimental result is observed, despite the highly turbulent nature of the flow in the present experiment. Further, amplitude saturation results from the time-averaged distortion imposed on the flow by the PVC, suggesting that linear stability analysis may predict PVC characteristics for $S>S_{c}$ . 
    more » « less
  4. We experimentally investigate the rise velocity of finite-sized bubbles in turbulence with a high energy dissipation rate of $\unicode[STIX]{x1D716}\gtrsim 0.5~\text{m}^{2}~\text{s}^{-3}$ . In contrast to a 30–40 % reduction in rise velocity previously reported in weak turbulence (the Weber number ( $We$ ) is much smaller than the Eötvös number ( $Eo$ ); $We\ll 1 more » « less
  5. Geometric, kinematic and dynamic properties of focusing deep-water surface gravity wave packets are examined in a simplified model with the intent of deriving a wave breaking threshold parameter. The model is based on the spatial modified nonlinear Schrödinger equation of Dysthe ( Proc. R. Soc. Lond.  A, vol. 369 (1736), 1979, pp. 105–114). The evolution of initially narrow-banded and weakly nonlinear chirped Gaussian wave packets are examined, by means of a trial function and a variational procedure, yielding analytic solutions describing the approximate evolution of the packet width, amplitude, asymmetry and phase during focusing. A model for the maximum free surface gradient, as a function of $\unicode[STIX]{x1D716}$ and $\unicode[STIX]{x1D6E5}$ , for $\unicode[STIX]{x1D716}$ the linear prediction of the maximum slope at focusing and $\unicode[STIX]{x1D6E5}$ the non-dimensional packet bandwidth, is proposed and numerically examined, indicating a quasi-self-similarity of these focusing events. The equations of motion for the fully nonlinear potential flow equations are then integrated to further investigate these predictions. It is found that a model of this form can characterize the bulk partitioning of $\unicode[STIX]{x1D716}-\unicode[STIX]{x1D6E5}$ phase space, between non-breaking and breaking waves, serving as a breaking criterion. Application of this result to better understanding air–sea interaction processes is discussed. 
    more » « less