Search NASA⌕ Search

SEARCH · Search NASA

Results for “Wake Vortex”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

A vortex sheet based analytical model of the curled wake behind yawed wind turbines

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.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Numerical Modeling of a Two-Stage Ocean Current Turbine

The Equinox Ocean Turbines (EQOT) current energy converter has a unique design with power generation in two small-diameter turbines attached to the tips of a large-diameter passive rotor. This configuration offers some key advantages for capturing ocean currents. With no centrally placed generator, almost no reaction torque is required at the nacelle of the main large-diameter rotor, and the small-diameter tip turbine generators operate at a higher speed and lower torque. The physics that determine the performance and loads on the turbine are also unique. The interactions of the flow field between the two stages and the general architecture of the system cannot be captured with traditional mid-fidelity modeling tools. For design iterations and large sets of load cases, it is important to have mid-fidelity models that can capture the important phenomenon with enough accuracy to identify global trends. This work uses a limited set of high-fidelity computational fluid dynamics (CFD) simulations to help inform the selection of and construction of a custom mid-fidelity model. Mid-fidelity modeling approaches were verified by comparing key turbine performance quantities to those found with the CFD model. Hydrodynamic interactions of the two-stage rotor were identified through high-fidelity CFD modeling. This highlighted the impact of the main rotor tip vortex and wake on the secondary rotor apparent inflow. This results in a relative flow rotation and sharp deficit, that change the optimal secondary rotor rotation speed and adds unsteadiness to the blade loading respectively. Multiple mid-fidelity approaches were evaluated for their ability to capture these effects. A simple approximation of the combined-stage performance based on single-stage BEM provides a reasonable rough prediction, especially near the peak TSR values, with some larger discrepancy at higher TSRs. Predicting the combined-stage performance based on single-stage CFD data improves this prediction across the TSR range. Although the combined-stage modeling in OLAF was not successful in this stage of the project, it showed promise as a mid-fidelity method, assuming the parameters can be tuned to account for the significant differences in time and length scales between the main and secondary rotors. This may be addressed through code changes in future work. A significant finding from the OLAF work was the agreement between the vortex core radius values found independently via a parameter space search and via CFD. The technique of using single-stage secondary rotor BEM, with a custom inflow taken from single-stage main rotor CFD or OLAF, provides an efficient method to capture one-way coupled flow interactions. This method provided generally good predictions of the impact of the flow rotation on the secondary rotor but struggled to accurately predict the peaks of the unsteady load progression. Future work could include some superposition of a tuned main rotor trailing edge viscous wake into the custom inflow to better predict this interaction.

16 TIDAL AND WAVE POWER↗
DOE OSTI2023

Investigation of thermal hydraulic behavior of the High Temperature Test Facility's lower plenum via large eddy simulation

A high-fidelity computational fluid dynamics (CFD) analysis was performed using the Large Eddy Simulation (LES) model for the lower plenum of the High–Temperature Test Facility (HTTF), a ¼ scale test facility of the modular high temperature gas-cooled reactor (MHTGR) managed by Oregon State University. In most next–generation nuclear reactors, thermal stress due to thermal striping is one of the risks to be curiously considered. This is also true for HTGRs, especially since the exhaust helium gas temperature is high. In order to evaluate these risks and performance, organizations in the United States led by the OECD NEA are conducting a thermal hydraulic code benchmark for HTGR, and the test facility used for this benchmark is HTTF. HTTF can perform experiments in both normal and accident situations and provide high-quality experimental data. However, it is difficult to provide sufficient data for benchmarking through experiments, and there is a problem with the reliability of CFD analysis results based on Reynolds–averaged Navier–Stokes to analyze thermal hydraulic behavior without verification. To solve this problem, high-fidelity 3-D CFD analysis was performed using the LES model for HTTF. It was also verified that the LES model can properly simulate this jet mixing phenomenon via a unit cell test that provides experimental information. As a result of CFD analysis, the lower the dependency of the sub-grid scale model, the closer to the actual analysis result. In the case of unit cell test CFD analysis and HTTF CFD analysis, the volume-averaged sub-grid scale model dependency was calculated to be 13.0% and 9.16%, respectively. As a result of HTTF analysis, quantitative data of the fluid inside the HTTF lower plenum was provided in this paper. As a result of qualitative analysis, the temperature was highest at the center of the lower plenum, while the temperature fluctuation was highest near the edge of the lower plenum wall. The power spectral density of temperature was analyzed via fast Fourier transform (FFT) for specific points on the center and side of the lower plenum. FFT results did not reveal specific frequency-dominant temperature fluctuations in the center part. It was confirmed that the temperature power spectral density (PSD) at the top increased from the center to the wake. The vortex was visualized using the well-known scalar Q-criterion, and as a result, the closer to the outlet duct, the greater the influence of the mainstream, so that the inflow jet vortex was dissipated and mixed at the top of the lower plenum. Additionally, FFT analysis was performed on the support structure near the corner of the lower plenum with large temperature fluctuations, and as a result, it was confirmed that the temperature fluctuation of the flow did not have a significant effect near the corner wall. In addition, the vortices generated from the lower plenum to the outlet duct were identified in this paper. It is considered that the quantitative and qualitative results presented in this paper will serve as reference data for the benchmark.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Investigation of grid-based vorticity-velocity large eddy simulation off-body solvers for application to overset CFD

Accurately predicting unsteady wakes and vortex-dominated flows is essential to a wide range of engineering applications, including aircraft, rotorcraft, shipboard operations, bio-inspired unsteady flight and propulsion, wind turbines, and urban flows. While current CFD software can model the complete flow field and wake system, the computational costs incurred in high Reynolds number unsteady turbulent flow simulations often remain prohibitive for routine engineering use, particularly for applications involving moving components. Prior work has demonstrated that by adopting a vorticity-velocity formulation in a grid-based off-body flow solver (VorTran-M and VorTran-M2) one can lower these costs by several orders of magnitude when compared to conventional approaches. This paper describes the extensions made to VorTran-M2 to support turbulent flows, and associated benchmarking activity to assess its performance for problems involving strong stretching and diffusion processes, whose competing contributions to the vorticity field are core drivers of turbulent flow evolution. Predictions are presented for: (i) the Kida-Pelz problem whose inviscid form is of mathematical interest due to its apparent formation of singular flow in finite time; and (ii) the Taylor Green vortex arrangement, which has been extensively studied as a fundamental simulation challenge in the turbulent modeling community. Here, the results are used to evaluate the overall predictive ability and performance of two sub-grid scale models incorporated into VorTran-M2. Results indicate that the computational cost savings seen previously for inviscid and convection dominated problems extend to turbulent flow simulations supporting the viability of VorTran-M2 as a low cost means for accurately modeling the far-field and background flow, particularly when long duration vorticity evolution is of interest.

42 ENGINEERING↗
DOE OSTI2022

A time‐varying formulation of the curled wake model within the FAST.Farm framework

Abstract In this article, we present a time‐varying formulation of the curled wake model that we implemented in FAST.Farm. The curled wake model, originally developed for steady‐state conditions, is used to produce realistic wake profiles behind a wind turbine in yawed (or skewed) conditions. We begin by introducing the key elements of the FAST.Farm framework. Then, after briefly summarizing the original wake dynamics formulation of FAST.Farm based on a polar wake profile, we present the new time‐varying formulation of the curled wake model, compare the two, and highlight the differences with the original curled wake model. After discussing some implementation details, we present different applications with increasing levels of complexity: single turbine with uniform and turbulent inflow, fixed and transient yaw, and multiple turbines. We verify our results using the original FAST.Farm implementation and large‐eddy simulations. The results with the new curled wake model are improved compared to the original implementation, as they include cross‐flow velocities and wake asymmetry. Yet, large‐eddy simulation results show a more pronounced lateral convection of the wake and a stronger concentration of vorticity at the top vortex. The new curled wake implementation in FAST.Farm should enable the calculation of not only generator power but also wind turbine structural loads for applications involving intentional or unintentional skewed flow and wind‐farm control involving wake steering.

17 WIND ENERGY↗
DOE OSTI2021

Experimental Measurements of the Wake of a Sphere at Subcritical Reynolds Numbers

This work experimentally investigated the flow phenomena and vortex structures in the wake of a sphere located in a water loop at Reynolds numbers of Re = 850, 1,250, and 1,700. Velocity fields in the wake region were obtained by applying the time-resolved stereoscopic particle image velocimetry (TR-SPIV) technique. From the acquired TR-SPIV velocity vector fields, the statistical values of mean and fluctuating velocities were computed. Spectral analysis, two-point velocity–velocity cross-correlation, proper orthogonal decomposition (POD) and vortex identification analyses were also performed. The velocity fields show a recirculation region that decreases in length with an increase of Reynolds numbers. The power spectra from the spectral analysis had peaks corresponding to a Strouhal number of St = 0.2, which is a value commonly found in the literature studies of flow over a sphere. The two-point cross-correlation analysis revealed elliptical structures in the wake, with estimated integral length scales ranging between 12% and 63% of the sphere diameter. The POD analysis revealed the statistically dominant flow structures that captured the most flow kinetic energy. It is seen that the flow kinetic energy captured in the smaller scale flow structures increased as Reynolds number increased. The POD modes contained smaller structure as the Reynolds number increased and as mode order increased. In addition, spectral analysis performed on the POD temporal coefficients revealed peaks corresponding to St = 0.2, similar to the spectral analysis on the fluctuating velocity. The ability of POD to produce low-order reconstructions of the flow was also utilized to facilitate vortex identification analysis, which identified average vortex sizes of 0.41D for Re1, 0.33D for Re2, and 0.32D for Re3.

Engineering↗
DOE OSTI2024

Simulation of the Multi-Wake Evolution of Two Sandia National Labs/National Rotor Testbed Turbines Operating in a Tandem Layout

The future of wind power systems deployment is in the form of wind farms comprised of scores of such large turbines, most likely at offshore locations. Individual turbines have grown in span from a few tens of meters to today’s large turbines with rotor diameters that dwarf even the largest commercial aircraft. These massive dynamical systems present unique challenges at scales unparalleled in prior applications of wind science research. Fundamental to this effort is the understanding of the wind turbine wake and its evolution. Furthermore, the optimization of the entire wind farm depends on the evolution of the wakes of different turbines and their interactions within the wind farm. In this article, we use the capabilities of the Common ODE Framework (CODEF) model for the analysis of the effects of wake–rotor and wake-to-wake interactions between two turbines situated in a tandem layout fully and partially aligned with the incoming wind. These experiments were conducted in the context of a research project supported by the National Rotor Testbed (NRT) program of Sandia National Labs (SNL). Results are presented for a layout which emulates the turbine interspace and relative turbine emplacement found at SNL’s Scaled Wind Technologies Facility (SWiFT), located in Lubbock, Texas. The evolution of the twin-wake interaction generates a very rich series of secondary transitions in the vortex structure of the combined wake. These ultimately affect the wake’s axial velocity patterns, altering the position, number, intensity, and shape of localized velocity-deficit zones in the wake’s cross-section. This complex distribution of axial velocity patterns has the capacity to substantially affect the power output, peak loads, fatigue damage, and aeroelastic stability of turbines located in subsequent rows downstream on the farm.

Baruah, Apurva (ORCID:0000000252354068)↗
DOE OSTI2021

Mechanisms of Dynamic Near-Wake Modulation of a Utility-Scale Wind Turbine

The current study uses large eddy simulations to investigate the transient response of a utility-scale wind turbine wake to dynamic changes in atmospheric and operational conditions, as observed in previous field-scale measurements. Most wind turbine wake investigations assume quasi-steady conditions, but real wind turbines operate in a highly stochastic atmosphere, and their operation (e.g. blade pitch, yaw angle) changes constantly in response. Furthermore, dynamic control strategies have been recently proposed to optimize wind farm power generation and longevity. Therefore, improved understanding of dynamic wake behaviours is essential. First, changes in blade pitch are investigated and the wake expansion response is found to display hysteresis as a result of flow inertia. The time scales of the wake response to different pitch rates are quantified. Next, changes in wind direction with different time scales are explored. Under short time scales, the wake deflection is in the opposite direction of that observed under quasi-steady conditions. Finally, yaw changes are implemented at different rates, and the maximum inverse wake deflection and time scale are quantified, showing a clear dependence on yaw rate. To gain further physical understanding of the mechanism behind the inverse wake deflection, the streamwise vorticity in different parts of the wake is quantified. The results of this study provide guidance for the design of advanced wake flow control algorithms. The lag in wake response observed for both blade pitch and yaw changes shows that proposed dynamic control strategies must implement turbine operational changes with a time scale of the order of the rotor time scale or slower.

flow-structure interactions↗
DOE OSTI2022

Numerical investigation of the influence of shear and thermal stratification on the wind turbine tip‐vortex stability

Summary The interaction between wind turbine wakes and atmospheric turbulence is characterised by complex dynamics. In this study, two major components of the atmospheric boundary layer dynamics have been isolated, namely, the mean velocity profile shear and the thermal stratification, to examine their impact on the near‐wake development by undertaking a series of highly resolved large‐eddy simulations. Subsequently, instantaneous flow fields are extracted from the simulations and used to conduct Fourier analysis and proper orthogonal decomposition (POD) and compute the mean kinetic energy fluxes by different POD modes to better understand the tip‐vortex instability mechanisms. Our findings indicate that shear can significantly affect the breakup of the wind turbine tip‐vortices and the shape and stable length of the wake, whereas thermal stratification seems to only have limited contribution to the spatial development of the near‐wake field. Finally, our analysis shows that the applied perturbation frequency determines the tip‐vortex breakup location as it controls the onset of the mutual inductance instability.

17 WIND ENERGY↗
DOE OSTI2020

The curled wake model: equivalence of shed vorticity models

In this work, we explore the spanwise velocities in the wake of a yawed wind turbine. In the curled wake model, these motions are induced by a collection of vortices shed from the rotor plane. The direction of the vorticity generated by yaw is aligned with the main flow (streamwise) direction. The streamwise vorticity induces velocities in the spanwise directions. These are the motions responsible for creating the curled wake mechanism. In this work, we explore a more accurate formulation for this mechanism, using a vortex cylinder. Under certain assumptions, the new and original curled wake models yield the same mathematical formulation. Also, both models predict an elliptic distribution of vortex strength, where the main difference is the location of the vortices.

17 WIND ENERGY↗
DOE OSTI2022

Implications of Shear and Thermal Stratification on Wind Turbine Tip-Vortex Stability

The interaction between wind turbines in a wind farm through their wakes is a phenomenon that has been studied for decades and is still relevant today. Turbines clustered together in arrays will often operate in the wake of other upstream turbines which may lead to significant power losses and fatigue loads. For modern large-scale wind turbines, the mean shear velocity profile and thermal stratification are major components of the atmospheric boundary layer so it is important to understand their impact on near-wake development. Additionally, veer is present due to the rotation of the Earth. The impact of shear, thermal stratification and veer on the stable wake length of turbines with a dynamic control strategy is studied numerically in this work using a suite of highly resolved large-eddy simulations. Instantaneous flow fields are extracted from the simulations and used to conduct proper orthogonal decomposition (POD) and compute the mean kinetic energy fluxes by different POD modes to better understand the tip-vortex instability mechanisms. Our findings show that the dynamic pitch control scheme is able to shorten the stable wake length to about 1.5R in uniform flow. Shear can significantly affect the break up of wind turbine tip-vortices as well as the shape and stable length of the wake, whereas thermal stratification seems to only have limited contribution to the spatial development of the near-wake field. Veer causes the wake boundary to skew but has a limited impact on the wake length.

length thermal stratification↗
DOE OSTI2020

Assessing the blockage effect of wind turbines and wind farms using an analytical vortex model

Wind farm blockage effects are currently neglected in the prediction of wind farm energy yield, typically leading to an overestimation of the production. This work presents a novel method to assess wind farm production, while accounting for blockage effects. We apply a vortex model, based on a cylindrical wake, to assess induction effects. We present variations of the model to account for finite wake length, finite tip-speed ratios, and the proximity to the ground. The results are applied to single rotors in aligned and yawed conditions and to different wind farm layouts. We provide far-field approximations for faster estimates of the velocity field. Further, this article includes a new methodology to couple the induction model to engineering wake models, such as the ones present in the FLOw Redirection and Induction in Steady State (FLORIS). We compare the results to actuator disk simulations for various operating conditions of a single turbine and different wind farm layouts. We found that the mean relative error of the model in the induction zone is typically around 0.2% compared with actuator disk simulations. The computational time of the velocity field using the analytical vortex model is three orders of magnitude less than the one obtained with the actuator disk simulation.

17 WIND ENERGY↗
DOE OSTI2022

Do ambient shear and thermal stratification impact wind turbine tip-vortex breakdown?

Modern wind turbines experience uneven inflow conditions across the rotor, due to the ambient flow’s shear and thermal stratification. Such conditions alter the shape and length of turbine wakes and thus impact the loads and power generation of downstream turbines. To this end, understanding the spatial evolution of the individual wakes under different atmospheric conditions is key to controlling and optimising turbine arrays. With this numerical study we aim to obtain a better understanding of the fundamental physics governing the near-wake dynamics of wind turbines under shear and thermal stability, by examining their tip-vortex breakup mechanisms. Our approach considers scale-resolving simulations of a single turbine wake under a linear shear profile as well as the application of harmonic tip perturbations to trigger flow instabilities. For the subsequent analysis we use the proper orthogonal decomposition (POD) method to extract coherent structures from the flow, and we also calculate mean kinetic energy fluxes to quantify each coherent structure’s contribution to wake recovery. The wake’s helical spiral is found to hinder wake recovery for all studied ambient flow conditions, whereas the mutual inductance instability has positive MKE flux leading to an enhanced wake recovery. Finally, the ambient shear has the largest impact on the local MKE flux with respect to downstream location by changing the shape of the curve and location of extrema, whereas thermal stratification has only a minimal impact on the magnitude of the near-wake local MKE flux distribution.

17 WIND ENERGY↗
DOE OSTI2025

Freestream turbulence effects on unsteady wind turbine loads and wakes: An IDDES study

We investigate numerically the effects of freestream turbulence on the unsteady aerodynamics and wakes of the National Renewable Energy Laboratory Phase VI wind turbine rotor for increasing wind speed. Turbulence is modeled using the Improved Delayed Detached-Eddy Simulation (IDDES) method. As a first step, a detailed mesh resolution study is conducted with the decaying freestream turbulence model at turbulence intensity of 0.5%. Our blade-resolved IDDES simulations show that grid-independent average torque and thrust results can be achieved with relatively coarse meshes, whereas dramatically higher mesh resolution is required for grid-independent results for power spectral densities of thrust force, especially in the deep-stall regime. Comparing the loads with the Shear-Stress Transport model demonstrates the superiority of IDDES in predicting massively separated flows. The aerodynamic performance and wake predictions with the decaying freestream turbulence model are compared with the synthetic freestream turbulence model. Both models predict nearly the same loads, spectral energy content, and wake characteristics. The properties of both the near- and far-wake regions are then examined. Furthermore, we show that separated boundary layers accelerate turbulent mixing and entrainment of the external flow, which results in faster wake recovery. The effect of increasing turbulence intensity to 6% is investigated using the synthetic freestream turbulence model. In contrast with the fully attached boundary layer, higher freestream turbulence in deep stall does not significantly affect the loads and vortex-shedding characteristics. However, the turbulent mixing in the wake is enhanced, which further hastens the recovery of the self-similar velocity profile. In general, increasing the wind speed at high turbulence intensity shifts the recovery farther upstream and increases the wake width.

17 WIND ENERGY↗
DOE OSTI2021

Experimental and computational investigation of interactive flow induced oscillations of two tandem rough cylinders at 3×10 4 ≤Re≤1.2×10 5

The interactive flow induced oscillations (FIO) of two adjacent, elastically mounted, rigid, tandem, locally-rough cylinders in transverse flow direction are analyzed utilizing two-dimensional Unsteady Reynolds-Averaged Navier-Stokes (2-D URANS) equations and verified experimentally in the proximity-wake region. Three sets of experiments and simulations (K = 600N/m, ζ = 0.14; K = 600N/m, ζ = 0.26; K = 1,200N/m, ζ = 0.26; K is spring stiffness, ζ is damping ratio) of two tandem cylinders with turbulence stimulation are tested and simulated for Reynolds number in the range of 30,000 = Re ≤ 120,000. The reduced velocity range is U* = 2.48–14.22, the mass ratio is m* = 1.343, and the center-to-center in-flow spacing to diameter ratio is d/D = 2.57. The characteristics of amplitude response, frequency response, lift force, and interactive wake patterns are presented and discussed. The trends of the amplitude and frequency responses from numerical simulations are in good agreement with experimental results. The main conclusions of Reynolds number effect on interactive flow induced oscillations are: (1) Five significant flow patterns between two tandem locally-rough cylinders for different Reynolds numbers are observed through analyzing the complex but stable interactions between vortices and cylinders. (2) In the initial and upper VIV branches, the downstream cylinder's FIO is seriously interfered by the wake of the upstream cylinder. (3) The downstream cylinder is strongly impinged by the vortices shed from the upstream cylinder resulting in nearly 180O out-of-phase oscillations in transition from VIV to galloping and in-phase oscillations in galloping.

42 ENGINEERING↗
DOE OSTI2020

Wind farm blockage effects: comparison of different engineering models

The work presents four engineering methods to estimate the induction zone in front of a wind turbine and account for the wind farm blockage effect. The methods comprise the vortex cylinder model, vortex dipole model, self-similar model, and wake projection model. The majority of the models presented account for yaw misalignments and ground effect. Actuator disk simulations are used to verify the individual models. The performance of each model is evaluated both in terms of precision and computational time. The induction models are coupled to wake models within the FLOw Redirection and Induction in Steady State framework to provide the full velocity field within a wind farm. Sample wind farm computations are presented, and the impact of including induction effects into wind farm performance predictions is reported. The different codes are publicly available online.

17 WIND ENERGY↗
DOE OSTI2024

Three-Dimensional Aerodynamics and Vortex-Shedding Characteristics of Wind Turbine Airfoils over 360-Degree Angles of Attack

In this work, we present the first three-dimensional (3D) computational investigation of wind turbine airfoils over 360° angles of attack to predict unsteady aerodynamic loads and vortex-shedding characteristics. To this end, static–airfoil simulations are performed for the FFA-W3 airfoil family at a Reynolds number of 107 with the Improved Delayed Detached Eddy Simulation turbulence model. Aerodynamic forces reveal that the onset of boundary-layer instabilities and flow separation does not necessarily coincide with the onset of stall. In addition, a comparison with two-dimensional simulation data and flat plate theory extension of airfoil polars, suggest that, in the deep stall regime, 3D effects remain critical for predicting both the unsteady loads and the vortex-shedding dynamics. For all airfoils, the vortex-shedding frequencies are found to be inversely proportional to the wake width. In the case of slender airfoils, the frequencies are nearly independent of the airfoil thickness, and their corresponding Strouhal number St is approximately 0.15. Based on the calculated St, the potential for shedding frequencies to coincide with the natural frequencies of the International Energy Agency 15 MW reference wind turbine blades is investigated. The analysis shows that vortex-induced vibrations occur primarily at angles of attack of around ±90° for all airfoils.

17 WIND ENERGY↗
DOE OSTI2022

Fluid structure interaction of a subaqueous pendulum: Analyzing the effect of wake correction via large eddy simulations

The dynamic behavior of a subaqueous cylindrical pendulum and corresponding flow dynamics are investigated. The objectives were twofold: (i) to examine whether the two-dimensional model equations sufficiently capture the three-dimensional dynamics and (ii) to investigate the emerging three-dimensional vortical flow structures. Large eddy simulations with two-way coupling fluid structure interaction were carried out using the immersed boundary method to simulate the motion of the pendulum and its interactions with the initially stagnant water. The resulting pendulum motion is compared against measured data obtained in a series of experimental tests to validate the simulation results and the model equations with and without wake corrections. An analysis of the flow vorticity revealed the development of a vortex ring during the first swing and the formation of tip vortices. The evolution of the vortex rings emerging from the motion of the subaqueous cylindrical pendulum was visualized using Q-criteria showing a reasonable agreement with vortical structures observed in the experiment using particle imaging velocimetry. The hydrodynamic moments acting on the simulated pendulum and the moments calculated from the model equations are analyzed. Using the insights from these numerical simulations, a modification of the wake correction is proposed to enhance the accuracy of the rate of decay and period. The transient effect of coherent flow on pendulum dynamics, especially the added mass effect, is discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗