Numerical simulations of wind farms consisting of innovative wind energy harvesting systems are conducted. The novel wind harvesting system is designed to generate strong lift (vertical force) with lifting-devices. It is demonstrated that the tip-vortices generated by these lifting-devices can substantially enhance wake recovery rates by altering the vertical entrainment process. Specifically, the wake recovery of the novel systems is based on vertical advection processes instead of turbulent mixing. Additionally, the novel wind energy harvesting systems are hypothesized to be feasible without requiring significant technological advancements, as they could be implemented as Multi-Rotor Systems with Lifting-devices (MRSLs), where the lifting-devices consist of large airfoil structures. Wind farms with these novel wind harvesting systems, namely MRSLs, are termed regenerative wind farm, inspired by the concept that the upstream MRSLs actively entrain energy for the downstream ones. With the concept of regenerative wind farming, much higher wind farm capacity factors are anticipated. Specifically, the results indicate that the wind farm efficiencies can be nearly doubled by replacing traditional wind turbines with MRSLs under the tested conditions, and this disruptive advancement can potentially lead to a profound reduction in the cost of future renewable energy.

The airfoil DU91-W2-150 was investigated in the Low Speed Low Turbulence Tunnel at the Delft University of Technology to study unsteady aerodynamics. This experimental study tested the airfoil under a wide range of angles of attack (AoA) from 0° to 310° at three Reynolds numbers ((Formula presented.)) from (Formula presented.) to (Formula presented.). Pressure on the airfoil surface was measured and particle image velocimetry (PIV) measurements were conducted to capture the flow field in the wake. By examining the force coefficient and comparing the wake contours, it shows that an upwind concave surface provides a higher load compared to a convex surface upwind case, highlighting the critical role of surface shape in aerodynamics. When comparing separation at specific locations along the chord for all three (Formula presented.) values, it is observed that as (Formula presented.) increases, separation tends to occur at lower AoA, both for positive stall and negative stall. The examination of the aerodynamic force variation indicates that, during reverse flow, fluctuations are more pronounced compared to forward flow. This is owing to separation occurring at the aerodynamic leading edge (geometric trailing edge) in reverse flow. In terms of vortex shedding frequency, the study found a nearly constant normalized Strouhal number ((Formula presented.)) of 0.16 across various (Formula presented.) and AoA values in fully separated regions, indicating a consistent pattern under these conditions. However, a slight increase in (Formula presented.), between 0.16 and 0.20, was observed for AoA values exceeding 180°, possibly due to the convex curvature of the airfoil in the upwind direction. In conclusion, this research not only corroborates previous findings for small AoA values but also adds new data on the aerodynamic behavior of the DU91-W2-150 airfoil under large AoA values, offering various perspectives on the effects of surface curvature, (Formula presented.), and flow conditions on key aerodynamic parameters.

Accurately determining experimental blade loading distributions is crucial for analyzing rotor performance but challenging due to the limitations of conventional measurement techniques. This paper presents a so-called wake-informed lifting line model that estimates blade loading distributions from phase-locked velocity measurements in the slipstream, eliminating the need for blade instrumentation. The model is evaluated against computational fluid dynamics (CFD) simulations under both attached and separated flow conditions. For the attached flow condition, the model achieves excellent agreement with CFD, with errors in the peak value of thrust distribution below 1%. In the separated flow condition, the model captures radial gradients and the shape of the thrust distribution but exhibits discrepancies in absolute values, with a 10% error in the peak value. These differences arise from the inherent limitations of the potential flow model, the increased significance of drag, and the heightened influence of the spinner’s presence in separated flows. Incorporating profile drag through external polar data improves the model prediction, reducing the error to 4%. The model cannot reliably predict power distributions without external polar data for both attached and separated flows due to the crucial role of drag in the torque direction. The application of the model to experimental flowfield data shows a performance similar to that of the validation case. Therefore, the wake-informed lifting line model offers a promising approach for obtaining experimental blade loading distributions, overcoming the limitations of traditional methods.

This study presents the experimental validation of regenerative wind farms (RGWFs), a novel wind farm concept designed to enhance overall wind farm performance. RGWFs employ Multi-Rotor Systems with Lifting devices (MRSLs), an innovative wind energy harvesters engineered to promote strong vertical energy entrainment, thereby accelerating wake recovery. In the experiments, MRSLs are scaled for wind tunnel testing, with their rotors modeled using porous disks and their lifting devices represented by airfoils. The tested RGWFs comprise up to three by three MRSLs. Flow quantities within RGWFs and aerodynamic loads on MRSLs are measured using volumetric particle tracking velocimetry and strain gauges, respectively. Compared to conventional wind farms, flow analysis indicates that RGWFs significantly enhance vertical energy entrainment, as evidenced by a more than 200 % increase in thrust on second-row-MRSLs. These experimental results, consistent with previous numerical predictions, highlight the promising potential of RGWFs.

Wind turbine blades in standstill or parked conditions often experience large angles of attack (AoA), where vortex-induced vibrations (VIV) may occur that increase the risk of structural damage. To better understand the VIV of airfoils at high AoA from an aerodynamic perspective, we conducted experimental investigations into the vortex dynamics of a surging airfoil at a 90∘ incidence undergoing forced vibrations. Experiments were conducted at two reduced frequencies (k) to demonstrate the lock-in effect, where the vortex shedding frequency aligns with the motion frequency. Results indicate distinct vortex shedding behaviors: at higher k value of 0.38, downstream wake vortices form when the airfoil is moving upwind, while upstream vortices emerge during the downwind motion, interacting with the downstream vortices and leading to an outward flow. At lower k value of 0.19, the wake remains directed to the downwind side, regardless of the airfoil’s motion direction. Lock-in is evident in both cases, with one vortex pair shed per cycle at lower k and two pairs at higher k. Furthermore, the study examines the influence of vortex dynamics on unsteady aerodynamic loads. The results show that drag peaks when the airfoil moves upwind near the center position of its trajectory; at higher k, negative drag occurs as the airfoil moves downwind near the center, driven by the interactions among convection, turbulent momentum, pressure, and viscous forces. A reduced-order load estimation model for a flat plate is applied to the experimental data, showing good agreement during the upwind motion of the airfoil, which is the design condition for the original flat plate model. However, during the downwind motion, as the flow condition does not match the original flat plate design condition, the circulatory part of the model is modified to account for the presence of two pairs of vortices in the flow field, yielding improved agreement with the drag values determined from the measured flow field. The findings highlight distinct flow patterns and vortex interactions for the two motion cases, offering insights into their impact on aerodynamic loads.

This study experimentally investigates the dynamic stall of a pitching NACA 643418 airfoil under reverse flow conditions, whereby the geometric trailing edge is located upstream of the geometric leading edge. The investigation focuses on the correlation between surface pressure, aerodynamic forces, and the evolution of the dynamic stall vortex (DSV) and the aerodynamic trailing edge vortex (TEV). A range of mean angles of attack from 5° to 25°, pitching amplitudes from 5° to 15°, and reduced frequencies between 0.05 and 0.21 were tested on an airfoil using a combination of particle image velocimetry (PIV) and surface pressure measurements. The results reveal a distinct dynamic stall mechanism in reverse flow conditions: multiple flow separations occur during both the upstroke and downstroke periods, yet the DSV maintains partial attachment. While the outer layer of the DSV sheds into the flow, its inner core remains anchored and is subsequently reinforced by shear layer feeding. This sustained vorticity injection leads to periodic DSV re-growth, which manifests in the force hysteresis loop as multiple distinct peak regions. In addition, it is also found that the largest difference between conventional and reverse flow dynamic stall lies in the initial stage of the DSV development. For reverse flow dynamic stall cases, the initial flow separation near the leading edge causes a gradual decrease in the pressure coefficient. However, for conventional dynamic stall cases, the laminar separation bubbles that occur before the DSV cause the maximum suction on the airfoil surface. Furthermore, by comparing the dominant modes from Proper Orthogonal Decomposition (POD) analysis, it is found that the type of flow regime depends mainly on the mean angle of attack within the tested range of pitching amplitudes and frequencies. In particular, the dominant vortex determines the flow dynamics. At low mean angles of 5° and 10°, the most energetic flow features are associated with DSV dynamics. At a mean angle of 15°, both DSV and TEV dynamics contribute. At higher mean angles of 20° and 25°, the flow is dominated more by TEV dynamics. Despite these mean-angle-dependent variations in modal dominance, all dominant POD modes share a consistent physical interpretation where they capture the growth or decay of DSV and TEV structures.

This study investigates the impact of zigzag tape on the aerodynamic performance and wake characteristics of the Delft Vertical Axis Wind Turbine (VAWT). The primary aim is to understand how the zigzag tape affects blade loads and the resulting aerodynamic wake. A comprehensive analysis was conducted using the Actuator Line Model (ALM) with airfoil characteristics measured in the wind tunnel at the Technical University of Denmark (DTU). Additionally, a 2-D CFD analysis with k-ω SST and γ-Re_θ turbulence models were employed to evaluate the influence of laminar transition phenomena on rotor characteristics. Results indicate that while the zigzag tape linearizes the lift coefficient characteristic, it leads to a notable reduction in aerodynamic efficiency due to increased drag and decreased lift below the critical angle of attack. The simulations were performed at a tip-speed ratio (TSR) of 4.5 to avoid a dynamic stall, as this operating condition ensures that the rotor blades remain below the static stall threshold and large offshore VAWTs are designed to operate near their maximum aerodynamic efficiency (C_P) for the majority of their operational time. The aerodynamic wake behind the rotor also shows significant changes, with the zigzag tape promoting asymmetry and affecting the wake recovery distance. The study's findings highlight the importance of considering surface contamination effects, represented by zigzag tape, in evaluating VAWT performance and wake behavior, offering valuable insights for wind turbine design and optimization.

The aerodynamics of the multi-rotor system with lifting-devices (MRSL), an innovative concept of wind energy harvesting machine, is preliminary investigated using Large Eddy Simulation (LES) with actuator techniques. In the current setup, turbulent inflow conditions are considered, but inflow wind shear is excluded. Consistent with previous studies, the results demonstrate faster wake recovery of the MRSL compared to its conventional counterpart, namely the wind turbine system without the lifting-devices. Additionally, a set of high-fidelity simulations further reveals that the enhanced wake recovery is robust under both laminar and turbulent inflow conditions, remaining largely unaffected by variations in the ambient turbulence level. The present work provides proof-of-concept evidence that the effectiveness of MRSLs is not significantly hindered by ambient turbulence, motivating future research to evaluate their performance within a realistic atmospheric boundary layer.

Wake losses are a significant source of inefficiency in wind farm arrays, hindering the development of high-energy-density wind farms offshore. Studies have demonstrated the potential of vertical-axis wind turbines (VAWTs) to achieve high-energy-density configurations, due to their increased rate of wake recovery compared with their horizontal-axis counterparts. Recent works have demonstrated a wake control technique for VAWTs that utilizes blade pitch to accelerate the wake recovery, hereinafter referred to as the “vortex-generator” method. The present work is an experimental investigation of the wake topology using this control technique for the novel X-Rotor VAWT. The time-averaged wake topology of the X-Rotor has been measured by stereoscopic particle-image velocimetry at three fixed-pitch conditions of the top blades, namely pitch-in, pitch-out, and a baseline case with no pitch applied. The results demonstrate the wake recovery mechanism linked to the streamwise vorticity system of the rotor and the mechanisms that lead to a streamwise momentum recovery, where the pitched-in case injects high-momentum flow from above the rotor while ejecting the wake from the sides. In contrast, the pitched-out case operates in a mirrored fashion, with high-momentum flow injected into the wake from the sides while low-momentum flow is ejected out axially above the rotor. These modes of operation demonstrate a significant increase in the available power for hypothetical downstream turbines, reaching as high as a factor of 2.2 two rotor diameters downstream compared with the baseline case. The pitched-in case exhibits a higher rate of momentum recovery in the wake, compared with the pitch-out configuration.

This study investigates the potential of regenerative wind farming using multirotor systems equipped with paired multirotor-sized wings, termed atmospheric boundary layer control (ABL-control) devices, positioned in the near-wake region of the multirotor. These ABL-control devices generate vortical flow structures that enhance vertical momentum flux from the flow above the wind farm into the wind farm flow, thereby accelerating the wake recovery process. This work presents numerical assessments of a single multirotor system equipped with various ABL-control configurations. The wind flow is modeled using steady-state Reynolds-averaged Navier-Stokes (RANS) computations, with the multirotor and ABL-control devices represented by three-dimensional actuator surface models based on momentum theory. Force coefficient data for the actuator surface models, as well as validation data for the numerical computations, were obtained from a scaled model at TU Delft's Open Jet Facility. The performance of the ABL-control devices was evaluated by analyzing the net momentum entrained from the flow above the wind farm and the total pressure and power available in the wake. The results indicate that, when the ABL-control strategy is employed, vertical momentum flux may become the dominant mechanism for wake recovery. In configurations with two or four ABL-control wings, the total wind power in the wake recovers to 95 % of the free-stream value at positions as early as x/D≈6 downstream of the multirotor system, representing a recovery rate that is approximately an order of magnitude faster than that observed in the baseline wake without ABL-control capabilities. It should be noted, however, that this study employs a simplified numerical setup to provide a proof of concept, and the current findings are not yet directly applicable to real-world scenarios.

The present study extends the idea of the vertical-axis wind turbine (VAWT) “vortex generator mode” for wake recovery on a wind farm scale, working towards the concept of “regenerative wind farming”, where upstream turbines entrain vertical momentum for those downstream. An experimental wind tunnel demonstration of the regenerative wind farming concept for an array of nine H-type VAWTs arranged in a 3×3 grid layout is performed. Volumetric particle tracking velocimetry measures the wake within the simulated wind farm while using two vortex generator modes achieved through a fixed blade pitch. The results demonstrate the strong dependence of the wake topology of a VAWT on the streamwise vorticity system, which can be effectively modified by pitching the blades and subsequently changing the load distribution of the different quadrants of a VAWT. An increase in momentum entrainment in the wake is observed for both vortex generator modes of operation, highlighting the potential of achieving regenerative wind farming. The derived available power within the farm increases by factors of 6.4 and 2.1 for the pitch-in and pitch-out cases compared to the baseline case, respectively, considering potential rotors directly in line with those upwind.

This paper studies the dynamic stall characteristics of thick flatback and nonflatback wind turbine airfoils. Two airfoils with a maximum thickness of (Formula presented.) were studied, with trailing edge thicknesses of (Formula presented.) and (Formula presented.), respectively. The static and dynamic experimental measurements were performed in the wind tunnel using surface pressure measurements for clean and tripped airfoils at the Reynolds number of (Formula presented.) and dynamic reduced frequency ranging from (Formula presented.) to (Formula presented.). The effects of the trailing edge gap, roughness, mean angle of attack, and reduced frequency on the dynamic stall characteristics of the airfoils were investigated. The results show that increasing the trailing edge gap delays the onset of dynamic stall. However, the lift loss after the onset of dynamic stall is for the flatback airfoil higher than the sharp trailing edge airfoil. Moreover, the flatback airfoil show higher lift overshoot compared to the sharp trailing edge airfoil in the dynamic stall condition. Increasing the reduced frequency affects the dynamic behavior both airfoils differently.

Conventional propellers operating at negative thrust conditions, even at 0 deg angle of attack, are characterized by flow separation and significantly different noise emissions than at positive thrust conditions. Operating the propeller at nonzero angles of attack at negative thrust conditions can further impact aerodynamic performance and far-field noise emission. This paper studies these effects using lattice-Boltzmann very large eddy simulations coupled with the Ffowcs Williams and Hawkings analogy. At positive thrust, operation at 10 deg angle of attack increases thrust along the freestream direction by approximately 3% compared to operation at 0 deg angle of attack, while efficiency remains constant. Conversely, the negative thrust condition shows approximately a 7% decrease in thrust magnitude and a 10% reduction in regenerated power. In this condition, the positively cambered blade sections exhibit dynamic stall, resulting in broadband fluctuations of up to 10% of the mean loading near the blade tip. The nonzero angle of attack induces opposite variations in absolute blade loading between positive and negative thrust conditions, resulting in opposite changes in the noise directivity. At positive thrust, noise increases in the region from which the propeller is tilted away (i.e., below the propeller at a positive angle of attack), while the opposite occurs at negative thrust. The varying blade loading over the azimuth results in destructive interference between loading and thickness noise for the negative thrust case at the 10 deg angle of attack. These findings highlight the crucial role of considering nonzero angles of attack in propeller design and optimization analyses.

This study examined the effect of vortex generators on the dynamic stall characteristics of thick wind turbine airfoils with a relative thickness of 35% and trailing edge thickness of 10% and 2%. The experiments were conducted in the TU Delft LTT wind tunnel at a Reynolds number of Re=1×10⁶ and dynamic reduced frequency ranging from 0.032 to 0.096. The study investigated the impact of various factors on the dynamic stall characteristics of the airfoils, including the vortex generator's chord position, trailing edge gap, roughness, mean angle of attack, and reduced frequency. The study found that vortex generators delay dynamic stall for thick airfoils by stabilizing the flow during the upstroke phase. However, this can increase the maximum lift overshoot, particularly with flatback airfoils, resulting in a higher drop in lift during dynamic stall. This can potentially increase the dynamic loads on a wind turbine blade due to stall-induced vibrations. The study noted a significant difference in dynamic stall behavior between flatback and non-flatback airfoils. Overall, this research provides valuable insights into the dynamic stall and flow physics characteristics of thick wind turbine airfoils using vortex generators, aiding in more accurate rotor blade design.

Large wind turbines face more intricate atmospheric conditions with turbulent coherent structures sized similarly to the rotor diameter, posing loading challenges. The present study assesses twelve distinct wind fields using the Large Eddy Simulations (LES) and International Electrotechnical Commission (IEC) Kaimal model scaled to their LES counterpart. The hub height wind speed in the different cases was set to 8.5 m/s (below-rated), 11.5 m/s (at-rated), and 14.5 m/s (above-rated). In a previous study, it was found that the unscaled IEC model-based wind field is conservative and scaled IEC model-based wind fields were found to yield different loads than upon use of LES-based wind fields in different atmospheric stability conditions. The present study aims to understand these differences. Utilizing Spectral Proper Orthogonal Decomposition (SPOD), the original wind fields were decomposed and reconstructed to study the influence of large and small coherent structures represented by their distinct frequencies. SPOD analysis was complemented by wind field spectral analysis considering atmospheric surface layer height, integral length scales, and co-coherence estimates. Integral length scales in the scaled IEC Kaimal model were found to be half of those in unstable atmosphere LES wind fields. The aero-elastic impact on the IEA 22 MW reference wind turbine with a 280 m rotor diameter was evaluated. The analysis reveals that large coherent structures, particularly low-frequency (≤0.06 Hz) ones, significantly impact wind turbine loads, contingent upon atmospheric stratification. Compared to the scaled IEC Kaimal model wind field, the maximum tower fore–aft bending moment and the maximum blade root flap-wise bending moment were found to be higher, for example, by 10% and 5% respectively in an unstable atmosphere during below-rated wind turbine operation. In the same scenario, standard deviation of the tower fore–aft bending moment was found to be higher by up to 50% while standard deviation of the blade root flap-wise bending moment was found to be lower by up to 25%. These findings underscore the critical importance of accurately modeling atmospheric turbulence and its coherent structures for more reliable design and operation of large wind turbines.

With the growing trend towards larger wind turbine rotor diameters, the impact of wind shear on rotor performance and loads becomes increasingly significant. Atmospheric stability strongly influences wind shear, leading to higher wind shear under stable atmospheric conditions. In this study, the aeroelastic performance of the IEA 22 MW rotor is assessed under inflow conditions generated by different methods. Inflow conditions were generated using turbulence models specified in the IEC Standards and also by Large Eddy simulations. Standalone OpenFAST simulations were conducted with the respective inflow conditions. It was found that at rated and above-rated wind speeds, the time-averaged wind turbine design loads were higher in stable atmospheric conditions, in comparison to the IEC NTM inflow conditions, while the opposite held for below-rated wind speeds. Specifically, the time-averaged root flapwise bending moment and rotor thrust were found to be higher by up to 7% in stable atmospheres. However, maximum design and fatigue loads were considerably higher in the IEC NTM case due to elevated turbulence levels. Compared to the IEC NTM case, the damage equivalent root flapwise bending moment was found to be 30% to 70% lower in the different scenarios.

The aim of this study was to assess the accuracy of predicting the aerodynamic loads and investigate the aerodynamic wake characteristics of a vertical axis wind turbine (VAWT) rotor using a simplified two-dimensional numerical rotor model and an advanced numerical approach – the Scale Adaptive Simulation (SAS) coupled with the four-equation γ – Re_θ turbulence model. The challenge for this approach lies in the operating conditions of the rotor, the blade pitch angles, and the very small geometric dimensions of the rotor. The rotor, with a diameter of 0.3 m, operates at a low tip speed ratio of 2.5 and an extremely low blade Reynolds number of approximately 22.000, whereas the pitch angles, β, are:-10, 0, and 10 degrees. Validation was conducted based on high-fidelity measurements obtained using the PIV technique at TU Delft. The obtained results of rotor loads and velocity profiles are surprisingly reliable for cases of β = 0° and β = –10°. However, the 2-D model is too imprecise to estimate both aerodynamic loads and velocity fields accurately.

As the demand for renewable energy increases, wind turbine rotors will become larger with slender blades. Vortex Generators (VGs) are used for passive flow control to avoid flow separation and reduce unsteady loading on the thick root section of slender blades due to their simplicity, inexpensiveness, and the ability to retrofit them to blades. Aerodynamic load calculations for VGs involve long experimental campaigns or resource intensive CFD calculations. Due to their inherent time-intensive nature aeroelastic optimisation and design tools prefer to use simplified but accurate analysis tools for aerodynamic load calculations of clean airfoils. One such class of tools are viscous-inviscid interaction solvers that use Integral Boundary Layer (IBL) methods for viscous calculations in the boundary layer coupled with an inviscid solver for the rest of the domain. The most popular example of this tool is XFOIL. An engineering model for VGs in IBL methods has previously been developed and implemented in XFOIL. In this research, the model parameters are tuned for implementation in another tool RFOIL based on XFOIL. RFOIL has been developed for accurate and robust analysis of wind turbine airfoils with improvements for thick airfoils and rotational corrections. The aerodynamic performance predicted the VG model in both XFOIL and RFOIL is then validated with an extensive database of airfoil data consisting of airfoils between 21% to 60% thickness, as well as Reynolds numbers between 1 million to 14 million, equipped with and without VGs. Finally, the computation time for the VG model is compared with that for a clean airfoil analysis in the same viscous-inviscid interaction solvers RFOIL and XFOIL. The investigation provides an overview of the usability of the engineering model in airfoil design methodologies in the wind turbine industry.

The topic of vortex-induced vibrations on a wind turbine blade has recently gained much attention due to its growing size and flexibility. To address this concern, a wind tunnel test was conducted to study the forced plunging and surging motion of a NACA0021 airfoil at 90° angle of attack. Results indicate that vortex lock-in occurred for a motion amplitude of one chord length even for a small frequency ratio (between motion frequency and static Strouhal frequency) of 0.39. Analysis of the drag coefficient, derived from the phase-averaged Particle Image Velocimetry (PIV) data, shows that a plunging airfoil experiences higher average loading than a surging airfoil, which is deemed to be more harmful considering the higher loading in the crossflow-direction due to the variation of effective angle of attack.

Hybrid Eulerian–Lagrangian solvers have gained increasing attention in the field of external aerodynamics, particularly when dealing with strong body–vortex interactions. This approach effectively combines the strengths of the Eulerian component, which accurately resolves boundary layer phenomena, and the Lagrangian component, which efficiently evolves the wake downstream. This study builds on our team's previous work by enhancing the capabilities of a two-dimensional hybrid Eulerian–Lagrangian solver. We aim to upgrade our solver which was initially designed for static cases, to now also simulate cases involving moving objects. To ensure the reliability and applicability of a new solver, it is essential to validate its performance in complex cases. Here, the solver is validated across the case of a traveling cylinder and the case of a rotating cylinder in two different rotational speeds at low Reynolds numbers. In the realm of Eulerian solvers, such as OpenFOAM (utilized for the Eulerian component of this hybrid approach), traditional techniques include the use of morphing meshes, overset meshes, and Arbitrary Mesh Interfaces (AMI) to model body motion. The proposed methodology involves extending the Eulerian mesh up to a short distance from the solid boundary and moving it entirely as a solid entity. Then the Lagrangian solver is responsible for calculating the updated boundary conditions, thereby completing the hybrid solver's functionality. This approach is very similar to the overset mesh technique. However, unlike the traditional method where an Eulerian mesh moves on top of a static one, our method involves the motion of an Eulerian mesh over a Lagrangian grid. We compared the results from our hybrid solver with those from a purely Eulerian solver, specifically OpenFOAM. The comparison demonstrates that our solver can replicate OpenFOAM's results with high accuracy. Another interesting point highlighted in this study is the presence of high-frequency oscillations in the body forces in hybrid solvers that incorporate the redistribution of Lagrangian particles and do not utilize surface elements such as vortex panels, specifically when dealing with dynamic mesh simulations. When the Eulerian mesh travels on top of the Lagrangian grid of particles, the positions of the particles with respect to the Eulerian mesh continuously change. This results in a constant shift of particles near the solid body, where the highest vorticity is observed. Particles that are close to the solid boundary at one time step may find themselves inside the boundary at the next time step, leading to their removal. This pattern continuously changes during the simulation, causing fluctuations in the boundary conditions of the Eulerian solver and manifesting as oscillations in the forces acting on the body. It is shown that this issue can be alleviated either by increasing the spatial resolution of the Lagrangian solver or by synchronizing the movement of the Lagrangian grid with the motion of the Eulerian mesh. The results of the study make the solver trustworthy and pave the way for more demanding external aerodynamic simulations.