Wake Region Research Articles

A new 3D asymmetric double-Gaussian wake analytical model for horizontal-axis wind turbines

A new three-dimensional analytical model named 3DADG that is able to describe asymmetric evolution of wind turbine wake is proposed in this paper. Due to the combined effects of wind wheel rotation and wind shear effects wake attenuation shows asymmetric properties. However, existing wake analytical models fail to describe such unique properties. To solve this problem, spatial dependencies of the wake expansion rate and deficit coefficient are revealed to successfully present horizontal attenuations of the wake varying from asymmetric-double-Gaussian in the near and transition wake regions, and then to single-Gaussian in the far wake region. The wind shear effect is also considered in the model to include asymmetric vertical wake speed attenuation. This model is derived from the laws of mass and momentum conservation. It is verified through comparisons with the results of wind tunnel “blind test” and LES simulation. The model proposed is capable to show spatial evolution of wind turbine wake from a three-dimensional perspective. Sensitivities of thrust coefficient and turbulence intensity to the model are also investigated. It not only promotes development of wake analytical model but also provides a useful tool for power prediction and wake control of wind farms under wake effects.

Retardation mechanisms and modeling of fatigue crack growth of a high-strength steel after single overload

This study investigated the retardation mechanisms and modeling methodology for the fatigue crack growth (FCG) of a high-strength Q500qE steel after a single overload. Different stress ratios (R = 0.1, 0.3, and 0.5) and overload ratios (OLR = 1.5, 2.0, and 2.5) were considered as key parameters. The digital image correlation (DIC) testing technique was used to capture the deformation at the crack tip during the tests. The crack closure response in the crack wake region was investigated using a virtual extensometer. The results show that, under an identical stress ratio, the fatigue life was increased as the OLR increased. The FCG was completely arrested when OLR reached 2.5 when R ≥0.3. Under an identical stress ratio, the residual strain before the crack tip and the residual deformation in the crack wake region increased with the rise of the overload ratio after a single overload. It’s found that the FCG was retarded or completely arrested due to the constraint of the crack tip plastic zone caused by single overload. Finally, a modified Wheeler model was proposed to predict the FCG rate based on the measured plastic zone size at the crack tip, and the predicted results agree well with the experimental findings.

Enhanced flow variability and morphological changes through individual wood placements on a gravel bed

To design and evaluate river restoration projects with wood, it is crucial to understand how the placement of wood will alter flow and bed morphology. This information is still missing. To address this, flume experiments were conducted with a single log placed on a mobile bed with a uniform grain size. The effect of log diameter, log orientation to the flow, and log length on flow variability and bed morphology was studied for various initial flow conditions. Even a single log increased flow variability, herein defined as spatial variability in both velocity and turbulent kinetic energy. For each log placed at the channel center, the logs generated a symmetric scour pattern alongside and downstream of the log due to increased bed shear and turbulence generated by the log. A steady, adjusted cross-section was reached at the cessation of erosion, which occurred when the bed stress decreased below the critical shear velocity. Logs with larger diameter tended to erode more deeply, so that the final, adjusted cross-section was similar for different log dimensions. Changes in bed topography affected wake characteristics. Logs with small changes in bed topography generated longer zones of reduced velocity, which can increase the residence time of organic matter or nutrients in the wake region. Logs with big changes in bed topography created the deepest and longest scour holes, which may locally reduce stream temperature and improve fish habitat during summer.

Influence of magnetohydrodynamics configuration on aerothermodynamics during Martian reentry

This paper investigates the role of magnetohydrodynamics (MHD) on the aerothermodynamics (ATD) of a representative entry vehicle while flying into the Martian atmosphere. By strategically placing a flight-ready superconducting magnet at varied positions in the Schiaparelli reentry capsule of the ExoMars mission, we discern its impact on essential flow properties. The primary consequence of MHD during atmospheric entry is the generation of the Lorentz force, which increases the shock standoff distance resulting in a reduction of the heat flux on the spacecraft by pushing high-energy plasma particles away. Through different magnet configurations, three distinct cases are formed to comprehensively understand the effects and implications of each setup. The study is performed using the COOLFluiD MHD for EnTries, an in-house ATD solver. For case 1, the magnet's placement behind the ExoMars forebody at the stagnation point reduces the heat flux. In case 2, the magnet's relocation to the shoulder region explores its potential to mitigate communication blackouts by influencing the wake region's flow. However, this positioning also induces shock bending, leading to variations in post-shock species mass fractions and heat flux spikes in the post-shock region. Case 3, involving an additional magnet where the shock bends in case 1, showcases a consistent increase in shock standoff distance across the forebody, providing a longer relaxation zone for species equilibration. Our findings highlight that while the strength of the applied magnetic field is crucial, the magnet's size is equally pivotal in determining ATD behavior. Case 3 emerges as the most promising configuration, consistently reducing heat flux across the forebody and maintaining it in the afterbody. This study underscores the potential of multi-magnet configurations as next-generation MHD heat shields for Martian atmospheric entry, emphasizing the criticality of magnet placement and configuration in enabling future MHD-enhanced deep space exploration missions.

Plasma Actuator Separation Control Investigated with Spectral Proper Orthogonal Decomposition

A single dielectric-barrier discharge plasma actuator is an active flow control device that imparts momentum to the fluid through ion acceleration using electromagnetic forces and has been used to suppress flow separation. This paper studies flow over an airfoil and how adding a single dielectric-barrier discharge actuator influences flow characteristics through numerical modeling. Using the spectral proper orthogonal decomposition and large-eddy simulation, flow instabilities are analyzed at their different temporal and spatial scales in the wake region. This study aims to explore the viability of spectral proper orthogonal decomposition for separation control and correlate the decomposed flow modes to the different actuation modes.

Experimental Study of Wake Evolution under Vertical Staggered Arrangement of Wind Turbines of Different Sizes

During the expansion of a wind farm, the strategic placement of wind turbines can significantly improve wind energy utilization. This study investigates the evolution of wake turbulence in a wind farm after introducing smaller wind turbines within the gaps between larger ones, focusing on aspects such as wind speed, turbulence intensity, and turbulence integral length scale. The flow field conditions are described using parameters like turbulence critical length and power spectral density, as determined through wind tunnel experiments. In these experiments, a single large wind turbine model and nine smaller wind turbine models were used to create a small wind farm unit, and pressure distribution behind the wind turbines was measured under various operating conditions. The results indicate that downstream wind speed deficits intensify as the number of small wind turbines in operation increases. The impact of these smaller turbines varies with height, with a relatively minor effect on the upper blade tip and increasingly adverse effects as you move from the upper blade tip to the lower blade tip. Through an analysis of power spectral density, the contribution of vortex motion to wake turbulence kinetic energy is further quantified. In the far wake region, the number of small wind turbines has a relatively small impact on wind speed fluctuations.

Parallel Electron Beams at Io: Numerical Simulations of the Dense Plasma Wake

AbstractIn 1995, the Galileo spacecraft traversed the wake of Io at ∼900 km altitude. The instruments onboard detected intense bi‐directional field‐aligned electron beams (∼140 eV–150 keV), embedded in a dense, cold and slow plasma wake (Nel ∼ 35,000 cm−3, Ti < 10 eV, V < 3 km/s). Similar electron beams were also detected along subsequent Galileo flybys. Using numerical simulations, we show that these electron beams are responsible for the formation of Io's dense plasma wake. We prescribe the composition of Io's atmosphere in S, O, SO and SO2, compute the atmospheric ionization by the beams with a parameterization adapted from study of auroral electrons at Earth, the plasma flow into Io's atmosphere with a Magneto‐Hydro‐Dynamic code, and the ion composition and temperature with a multi‐species physical chemistry code. Results reveal contrasting chemistries between the upstream and wake regions, leading to different ion compositions. The upstream chemistry is driven by the torus thermal electrons at 5 eV with SO2+ becoming the dominant ion because of electron‐impact ionization of the SO2 atmosphere. The wake chemistry is driven by the high‐energy electrons in the beams with S+ and SO+ becoming the dominant ions produced by dissociative‐ionization of SO2. We show that the wake ion composition is highly sensitive to the atmospheric composition. Juno, in its extended mission, will traverse Io's wake and determine its ion composition, which, compared with our numerical simulations will enable us to infer the detailed composition of the atmosphere.

Quantification of three-dimensional added turbulence intensity for the horizontal-axis wind turbine considering the wake anisotropy

The accurate evaluation of wake turbulence intensity (TI) is of great significance for the layout optimization of wind farms and research on control strategies. However, the existing engineering wake models still have many shortcomings in evaluating the TI in the wake of horizontal-axis wind turbines (HAWTs), especially in the near wake region. In response to this, this article proposed an analytical model based on the double-Gaussian ellipse function that can describe the three-dimensional added TI distribution in the entire wake region. Firstly, the anisotropy of wake expansion is taken into account in both the horizontal and vertical directions, assuming that the added TI distribution downstream of the HAWT is a double-Gaussian elliptical shape. Secondly, taking into account the self-similarity characteristics of flow TI, the maximum added TI and its position in the entire wake region are evaluated. In addition, considering the influence of the ground effect, the vertical turbulence distribution was corrected. Finally, the proposed added turbulence model was validated and relative error analysis was conducted using large eddy simulation (LES) data, wind field measurement data, and wind tunnel experimental data. The results show that the model had good prediction accuracy in the entire wake region, and its prediction performance was significantly improved compared to traditional models, especially in the near wake region, with the lowest even relative error of 3.66%. This model has significant advantages such as low computational cost, wide applicability, and high accuracy. It can reduce energy losses in wind farms, improve wind energy utilization efficiency, and has great potential for widespread application in large-scale wind farms.

Numerical Simulation of Cavitation Control around a Circular Cylinder Using Porous Surface by Volume Penalized Method

In this work, we conducted a numerical study on the cavitation flow around a circular cylinder with Re=200 and σ=1, through the implementation of a porous coating. The primary objective addressed the effectiveness of utilizing a porous surface to control cavitation. We analyzed the cavitation dynamics around the cylinder and the hydrodynamic performance at different permeability levels of the porous surfaces (K=10−12−10−10). The flow was governed by the density-based homogeneous mixture model, and the volume penalization method was used to deal with the porous layer. A high-order compact numerical method was adopted for the simulation of the cavitating flow through solving the preconditioned multiphase equations. The hydrodynamic findings demonstrated that the fluctuations in the lift coefficient decreased when the porous layer was applied. However, it is not possible to precisely express an opinion about drag because the drag coefficient may vary, either increasing or decreasing, depending on the permeability within a constant thickness of the porous layer. The results revealed that the application of a porous layer led to the effective suppression of cavitation vortex shedding. In addition, a reduction of the shedding frequency was obtained, which was accompanied by thinner and elongated vortices in the wake region of the cylinder. With the proper porous layer, the inception of the cavitation on the cylinder was suppressed, and the amplitude of pressure pulsations due to the cavitation shedding mechanism was mitigated.

A new exploration on passive control of transonic flow over a backward-facing step

PurposeThis paper aims to study passive control techniques for transonic flow over a backward-facing step (BFS) using square-lobed trailing edges. The study investigates the efficacy of upward and downward lobe patterns, different lobe widths and deflection angles on flow separation, aiming for a deeper understanding of the flow physics behind the passive flow control system.Design/methodology/approachLarge Eddy Simulation and Reynolds-averaged Navier–Stokes were used to evaluate the results of the study. The research explores the impact of upward and downward patterns of lobes on flow separation through the effects of different lobe widths and deflection angles. Numerical methods are used to analyse the behaviour of transonic flow over BFS and compared it to existing experimental results.FindingsThe square-lobed trailing edges significantly enhance the reduction of mean reattachment length by up to 80%. At Ma = 0.8, the up-downward configuration demonstrates increased effectiveness in reducing the root mean square of pressure fluctuations at a proximity of 5-step height in the wake region, with a reduction of 50%, while the flat-downward configuration proves to be more efficient in reducing the root mean square of pressure fluctuations at a proximity of 1-step height in the near wake region, achieving a reduction of 71%. Furthermore, the study shows that the up-downward configuration triggers early spanwise velocity fluctuations, whereas the standalone flat-downward configuration displays less intense crosswise velocity fluctuations within the wake region.Practical implicationsThe findings demonstrate the effectiveness of square-lobed trailing edges as passive control techniques, showing significant implications for improving efficiency, performance and safety of the design in aerospace and industrial systems.Originality/valueThis paper demonstrates that the square-lobed trailing edges are effective in reducing the mean reattachment length and pressure fluctuations in transonic conditions. The study evaluates the efficacy of different configurations, deflection angles and lobe widths on flow and provides insights into the flow physics of passive flow control systems.

Global dynamics and topology of two-phase mixing layer flow through simultaneous gas and liquid velocity measurements

This study reports the first time-resolved particle image velocimetry characterization of a planar two-phase mixing layer flow, whose velocity field is measured simultaneously in gas and liquid streams. Two parallel air and water flows meet downstream of a splitter plate, giving rise to an initially spanwise invariant configuration. The aim is to elucidate further the mechanisms leading to the flow breakup in gas-assisted atomization. The complete experimental characterization of the velocity field represents a database that could be used in data-driven reduced-order models to investigate the global behaviour of the flow system. After the analysis of a selected reference case, a parametric study of the flow behaviour is performed by varying the liquid ( $Re_l$ ) and gas ( $Re_g$ ) Reynolds numbers, and as a consequence also the gas-to-liquid dynamic pressure ratio ( $M$ ), shedding light on both time-averaged (mean) and unsteady velocity fields. In the reference case, it is shown that the mean flow exhibits a wake region just downstream of the splitter plate, followed by the development of a mixing layer. By increasing both $Re_l$ and $Re_g$ , the streamwise extent of the wake decreases and eventually vanishes, the flow resulting in a pure mixing layer regime. The spectral analysis of the normal-to-flow velocity fluctuations outlines different flow regimes by variation of the governing parameters, giving more insights into the global characteristics of the flow field. As a major result, it is found that at high $Re_g$ and $M$ values, the velocity fluctuations are characterized by low-frequency temporal oscillations synchronized in several locations within the flow field, which suggest the presence of a global mode of instability. The proper orthogonal decomposition of velocity fluctuations, performed in both gas and liquid phases, reveals finally that the synchronized oscillations are associated with a low-frequency dominant flapping mode of the gas–liquid interface. Higher-order modes correspond to interfacial wave structures travelling with the so-called Dimotakis velocity. For lower gas Reynolds numbers, the leading modes describe higher frequency fingers shedding at the interface.

Control of separated flow at low Reynolds number around NACA0012 airfoil by boundary layer suction

The separated flow at low Reynolds number around the NACA0012 airfoil is numerically studied by large-eddy simulation. Strategies of boundary layer suction to control flow separation are investigated. A method of using two-zone suctions, near the leading edge and near the trailing edge, are calculated. Based on verification with direct numerical simulation (DNS) and experimental data, the results of the lift and the drag, the vortices, and the strength of near-field pressure fluctuations, are checked. The results show that the two-zone suctions can supress flow transition and separation, thereby increase the lift and reduce the drag. The shedding of vortices is weakened, and the near-field pressure fluctuations are attenuated. For comparison with the two-zone suctions, the strategies of suction near the leading edge only and suction near the trailing edge only are also studied. It is found that suction near the leading edge only may suppress transition and delay separation when the suction zone is large enough, but the flow property deteriorates due to shedding vortices in the wake. The suction near the trailing edge only may improve the flow performance by reducing the size of the vortices in the rear section of the airfoil and in the wake region, but it has little effect on the separation bubble and transition.

Aerodynamic analysis of building with zigzag-patterned façade: Insights from large eddy simulation

High-rise buildings have diverse shapes and configurations to harmonize with the architectural, engineering, environmental, economic, and functional considerations, and aspirations. However, the architectural shape and configuration plays a vital role in determining the aerodynamic characteristics of the surrounding flow field and wind forces acting on the building. This paper numerically investigates the aerodynamic characteristics of the zigzag-patterned building by adopting the three-dimensional large eddy simulation (LES). The flow around a square cylinder in high Reynolds number is used to valid the numerical method. The results show that the zigzag pattern has a little effect on the mean streamwise velocity but significantly influences the fluctuation of streamwise and vertical velocities in the wake region. It is observed that the zigzag pattern on the leeward wall results in a narrower wake region, while on the windward wall, leads to a thinner and longer shear layers and shorter recirculation region. As for the pedestrian level wind, the flow velocity within the separation area is noticeably reduced when incorporating a zigzag pattern on the windward wall. Moreover, the fluctuating pressure coefficient values at the sidewalls of a building with a zigzag windward wall is significantly smaller than that of a smooth configuration building due to the façade pattern's impact in diminishing the curvature of the nearby separation at the leading edge.

Validation and Assessment of a Hybrid VOF-Lagrangian Numerical Methodology for Turbulent Liquid Fuel Jets in High-Speed Crossflow

Abstract For gas turbine combustors, injecting liquid fuel in crossflow can achieve superior mixing characteristics to improve performance and reduce emissions. Robust numerical design tools are needed to accelerate the development of low-emissions technologies. Atomization modeling is often facilitated by using the Lagrangian approach rather than the more complex and computationally expensive Volume-of-Fluid (VOF) approach. However, most Lagrangian breakup models rely on empirical constants that must be fully calibrated for jets in crossflow and representative conditions. A hybrid approach could deliver a better compromise between accuracy and computational cost. This study proposes and validates a novel numerical methodology for coupling the VOF and Lagrangian approaches using a stochastic breakup model with Adaptive Mesh Refinements for turbulent liquid fuel jets in crossflow under more representative gas turbine conditions. The predictions were validated in the near and far-field regions using datasets at high pressures (1–8 bar), Weber numbers (720–1172), and momentum flux ratios (6–33). The predictive capabilities and computational cost were also compared to the Lagrangian approach coupled with large eddy simulations (LES) and with the unsteady Reynolds-Averaged Navier–Stokes (RANS) methodology previously developed and validated by the authors. The effect of different VOF-Lagrangian transition criteria on the computational cost was also assessed and recommendations were provided for further improvements. The overall LES predictions were significantly improved by the proposed hybrid methodology. Although it tends to underpredict the spray trajectory and Sauter Mean Diameter compared to the URANS methodology, it better captures the diameter in the wake region and the droplet velocities.

Predicting the upper bound of two-dimensional flow regimes of symmetric objects through two-dimensional computations

The onset of secondary wake instability is generally predicted via experiments, linear stability analysis, and three-dimensional direct numerical simulations. The current work stems from an open question that is very intriguing and fundamental: Can the upper bound of a two-dimensional flow be predicted purely on the basis of two-dimensional computational results? It is found that spatial distribution of a field variable, i.e., streamwise velocity in the vortex formation region, aids in determining the upper limit of a two-dimensional flow regime of a symmetric object. The vortex formation length attains its least value at the second critical Reynolds number. In addition, streamwise extents of mean wake and vortex formation region along wake axis become the same. Under this circumstance, the streamwise velocity at the terminal point of vortex formation region is such that its mean value vanishes while intensity of fluctuations or corresponding Reynolds stress becomes the maximum. The predicted values of critical Reynolds numbers for circular, square, and diamond cross sections exhibit excellent agreement with the results available in the literature.

Influence of the suspension gap on the wake characteristics of a 600 km/h superconducting maglev train

Although the suspension gap between a maglev train and its track is narrow, it has a significant impact on the aerodynamic performance of the train. This study investigated the effect of the suspension gap on the wake flow of a superconducting maglev train traveling at a speed of 600 km/h by comparing the aerodynamics and flow field structure of a tail car (TC) with suspension gaps of 8, 10, 12, and 14 cm. The results showed that with increasing suspension gap, the vortex structures in the wake region move closer to the TC, the airflow from the suspension gap inhibits the high-pressure region of the wake, and the pressure and velocity fluctuations significantly increase near the nose tip of the TC. Moreover, the suspension gap was found to have a significant impact on lift and almost no effect on drag. Increasing the suspension gap leads to a balance in the lift between the head car (HC) and the TC; as the gap continues to increase, the lift ratio of the HC and TC rapidly decreases.

Turbulent Characteristics of a Submerged Reef under Various Current and Submergence Conditions

Submerged Reefs (SRs) are a kind of artificial fish habitat that can protect coasts and maintain ecological biodiversity. In this study, the flow field of the SR is simulated by solving a Reynolds-averaged Navier–Stokes equation closed with the Realizable k-ε model based on the finite volume method. The turbulent characteristics of SRs under different inflow velocities and submergences in the vicinity of the SR are analyzed. The wake vorticities are the primary turbulent pattern within and around the SR. The back wake and vorticity are chosen as critical indicators to quantitatively assess the hydrodynamic characteristics induced by the SR. The results show: (1) as the main flow passes through the SR, the upwelling is produced in front of the SR and a large-scale wake region is formed behind the SR which contains a clockwise vortex; (2) the length of the wake region formed behind the SR is positively and linearly correlated with both the inflow velocity and submergence; (3) the dipole-type vorticity patterns are induced within the compartment of the SR, where the area and average value of high vorticity have a positive correlation with the flow velocity and a negative correlation with the submergence, respectively.

An advanced three-dimensional analytical model for wind turbine near and far wake predictions

Most existing analytical wake models have been designed to represent the velocity distribution of the far wake region behind a wind turbine, but fail to reproduce the trend of the near wake. In order to meet the demand for both near and far wake velocity predictions, a concise three-dimensional dual-cosine shape (3D-COU) model is proposed. Moreover, by taking into account the wind shear and ground effects in the vertical direction, this model is capable of describing the anisotropic property of the 3-D wake field and is much closer to the reality. The proposed model is then calibrated and validated through a series of test cases, including different scales of turbines that are operating under a wide range of inflow conditions, and meanwhile the reference data resource is diverse. Overall, the results demonstrate that the 3D-COU model can accurately predict the variations of the wake velocity in both the vertical and horizontal directions throughout the whole wake region. Given its good accuracy, simplicity and applicability, this model provides the potential to effectively address practical problems in wind energy projects.

Particle-resolved study of the onset of turbulence

The transition from laminar to turbulent flow is an immensely important topic that is still being studied. Here we show that complex plasmas, i.e., microparticles immersed in a low temperature plasma, make it possible to study the particle-resolved onset of turbulence under the influence of damping, a feat not possible with conventional systems. We performed three-dimensional (3D) molecular dynamics (MD) simulations of complex plasmas flowing past an obstacle and observed 3D turbulence in the wake and forewake region of this obstacle. We found that we could reliably trigger the onset of turbulence by changing key parameters such as the flow speed and particle charge, which can be controlled in experiments, and show that the transition to turbulence follows the conventional pathway involving the intermittent emergence of turbulent puffs. The power spectra for fully developed turbulence in our simulations followed the −5/3 power law of Kolmogorovian turbulence in both time and space. We demonstrate that turbulence in simulations with damping occurs after the formation of shock fronts, such as bow shocks and Mach cones. By reducing the strength of damping in the simulations, we could trigger a transition to turbulence in an undamped system. This work opens the pathway to detailed experimental and simulation studies of the onset of turbulence on the level of the carriers of the turbulent interactions, i.e., the microparticles. Published by the American Physical Society 2024

Effect of typical arch structure on slipstream and wake flow of 600 km/h maglev train

PurposeThe aerodynamic differences between the head car (HC) and tail car (TC) of a high-speed maglev train are significant, resulting in control difficulties and safety challenges in operation. The arch structure has a significant effect on the improvement of the aerodynamic lift of the HC and TC of the maglev train. Therefore, this study aims to investigate the effect of a streamlined arch structure on the aerodynamic performance of a 600 km/h maglev train.Design/methodology/approachThree typical streamlined arch structures for maglev trains are selected, i.e. single-arch, double-arch and triple-arch maglev trains. The vortex structure, pressure of train surface, boundary layer, slipstream and aerodynamic forces of the maglev trains with different arch structures are compared by adopting improved delayed detached eddy simulation numerical calculation method. The effects of the arch structures on the aerodynamic performance of the maglev train are analyzed.FindingsThe dynamic topological structure of the wake flow shows that a change in arch structure can reduce the vortex size in the wake region; the vortex size with double-arch and triple-arch maglev trains is reduced by 15.9% and 23%, respectively, compared with a single-arch maglev train. The peak slipstream decreases with an increase in arch structures; double-arch and triple-arch maglev trains reduce it by 8.89% and 16.67%, respectively, compared with a single-arch maglev train. The aerodynamic force indicates that arch structures improve the lift imbalance between the HC and TC of a maglev train; double-arch and triple-arch maglev trains improve it by 22.4% and 36.8%, respectively, compared to a single-arch maglev train.Originality/valueThis study compares the effects of a streamlined arch structure on a maglev train and its surrounding flow field. The results of the study provide data support for the design and safe operation of high-speed maglev trains.