Generation and Absorption of Periodic Waves Traveling on a Uniform Current in a Fully Nonlinear BEM-based Numerical Wave Tank

Accurate and efficient numerical wave generation and absorption of two-dimensional nonlinear periodic waves traveling on a steady, uniform current were carried out in a potential, fully nonlinear numerical wave tank. The solver is based on the Βoundary Εlement Μethod (ΒΕΜ) with linear singularity distributions and plane elements and on the mixed Eulerian–Lagrangian formulation of the free surface equations. Wave generation is implemented along the inflow boundary by imposing the stream function wave solution, while wave absorption at both end-boundaries is effectively treated by introducing absorbing layers. On the absorbing beach side, the outflow boundary condition is modified to ensure that the solution accurately satisfies the dispersion relation of the generated waves. The modification involves a free-parameter that depends on the mass flux through the domain and is determined through a feedback error-correction loop. The developed method provides accurate time domain wave solutions for shallow, intermediate, and deep water depths of high wave steepness (wave heights up to 80% of the maximum value) that remain stable for 150 wave periods. This also holds in case a coplanar or opposing uniform current of velocity up to 20% of the wave celerity interacts with the wave

Investigating the Level of Fidelity of an Actuator Line Model in Predicting Loads and Deflections of Rotating Blades under Uniform Free-Stream Flow

In this paper, the accuracy of an in-house Actuator Line (AL) model is tested on aeroelastic simulations of a Wind Turbine (WT) rotor and a helicopter Main Rotor (MR) under uniform free-stream flow. For the scope of aeroelastic analyses, the AL model is coupled with an in-house multibody dynamics code in which the blades are modeled as beams. The advantage from the introduction of CFD analysis in rotorcraft aeroelasticity is related to its capability to account in detail for the interaction of the rotor wake with the boundary layer developed on the surrounding bodies. This has proven to be of great importance in order to accurately estimate the aerodynamic forces and thus the corresponding structural loads and deflections of the blades. In wind turbine applications, a good example of the above is the rotor/ground interaction. In helicopter configurations, the interaction of MR with the ground or the fuselage and the interaction of tail rotor with the duct in fenestron configurations are typical examples. Furthermore, CFD aerodynamic analysis is an obvious modeling option in which the above mentioned asset can be combined with the consideration of the mutual interaction of the rotor with the ambient turbulence. A WT rotor operating inside the atmospheric boundary layer under turbulent free-stream flow is such a case. In the paper, AL results are compared against Blade Element Momentum (BEM) and Lifting Line (LL) model results in the case of the WT, whereas LL and measured data are considered in the helicopter cases. Blade loads and deflections are mainly compared as azimuthal variations. In the helicopter MR cases, where comparison is made against experimental data, harmonic analysis of structural loads is shown as well. Overall, AL proves to be as reliable as LL in the canonical cases addressed in this paper in terms of loads and deflections predictions. Therefore, it can be trusted in more complex flow conditions where viscous effects are pronounced

Assessment of a Hybrid Eulerian–Lagrangian CFD Solver for Wind Turbine Applications and Comparison with the New MEXICO Experiment

In this paper, the hybrid Lagrangian–Eulerian solver HoPFlow is presented and evaluated against wind tunnel measurements from the New MEXICO experiment. In the paper, the distinct solvers that assemble the HoPFlow solver are presented, alongside with details on their mutual coupling and interaction. The Eulerian solver, MaPFlow, solves the compressible Navier–Stokes equations under a cell-centered finite-volume discretization scheme, while the Lagrangian solver uses numerical particles that carry mass, pressure, dilatation and vorticity as flow markers in order to represent the flow-field by following their trajectories. The velocity field is calculated with the use of the decomposition theorem introduced by Helmholtz. Computational performance is enhanced by utilizing the particle mesh (PM) methodology in order to solve the Poisson equations for the scalar potential ϕ and the stream function ψ→. The hybrid solver is tested in 3-D unsteady simulations concerning the axial flow around the wind turbine (WT) model rotor tested in the New MEXICO experimental campaign. Simulation results are presented as integrated rotor loads, radial distribution of aerodynamic forces and moments and pressure distributions at various span-wise positions along the rotor blades. Comparison is made against experimental data and computational results produced by the pure Eulerian solver. A total of 5 PM nodes per chord length of the blade section at 75% have been found to be sufficient to predict the loading at the tip region of the blade with great accuracy. Discrepancies with respect to measurements, observed at the root and middle sections of the blade, are attributed to the omission of the spinner geometry in the simulations

Simulation of oscillating trailing edge flaps on wind turbine blades using ranging fidelity tools

Abstract Recent research developments have indicated that substantial reduction of both the fatigue and ultimate loads can be achieved by adopting trailing edge (TE) flap control strategies. Their aeroelastic tools employ blade element momentum (BEM) aerodynamic models enhanced with a sectional 2D treatment of the TE flap, neglecting the 3D effect of the trailed vorticity in the vicinity of the moving flap. In the present paper, a cross comparison of the BEM‐based models used in the aeroelastic analysis tools against higher fidelity, free‐wake lifting line, and fully resolved CFD models is performed, with the aim to highlight limitations of the first. A second level of comparison assesses the differences among tools of the same level of fidelity from different research groups. Moreover, a number of engineering‐based correction models that are used in conjunction with BEM and account for the complex 3D trailed vorticity effect are assessed. Simulations of a stiff rotor configuration of the DTU 10 MW Reference Wind Turbine are performed for a prescribed, harmonic TE flap motion, and aerodynamic loads are compared at the sectional and rotor‐integrated level. For the studied stiff rotor with the chosen flaps configuration, the results of the code‐to‐code comparisons indicate that low‐fidelity BEM tools consistently predict 1P variations of the rotor thrust due to the TE flap motion, but fail to reproduce the details of the load distributions especially in the vicinity of the flap section. BEM‐based corrected models, which account for 3D‐induced velocity effects, provide load distribution predictions closer to higher fidelity free‐wake and CFD models