Pressure-based Solver Research Articles

Non-Condensation Turbulence Models with Different Near-Wall Treatments and Solvers Comparative Research for Three-Dimensional Steam Ejectors

Steam ejectors are important energy-saving equipment for solar thermal energy storage; however, a numerical simulation research method has not been agreed upon. This study contributes to a comprehensive selection of turbulence models, near-wall treatments, geometrical modeling (2-D and 3-D), solvers, and models (condensation and ideal-gas) in the RANS equations approach for steam ejectors through validation with experiments globally and locally. The turbulence models studied are k-ε Standard, k-ε RNG, k-ε Realizable, k-ω Standard, k-ω SST, Transition SST, and linear Reynolds Stress. The near-wall treatments assessed are Standard Wall Functions, Non-equilibrium Wall Functions, and Enhanced Wall Treatment. The solvers compared are pressure-based and density-based solvers. The root causes of their distinctions in terms of simulation results, applicable conditions, convergence, and computational cost are explained and compared. The complex phenomena involving shock waves, choking, and vapor condensation captured by different models are discussed. The internal connections of their performance and flow phenomena are analyzed from the mechanism perspective. The originality of this study is that both condensation and 3-D asymmetric effects on the simulation results are considered. The results indicate that the k-ω SST non-equilibrium condensation model coupling the low-Re boundary conditions has the most accurate prediction results, best convergence, and fit for the widest range of working conditions. A 3-D asymmetric condensation model with a density-based solver is recommended for simulating steam ejectors accurately.

Charging of an Air–Rock Bed Thermal Energy Storage under Natural and Forced Convection

An air-rock bed thermal storage system was designed for small-scale powered generation and analyzed with computational fluid dynamics (CFD) using ANSYS-Fluent simulation. An experimental system was constructed to compare and validate the simulation model results. The storage unit is a cylindrical steel container with granite rock pebbles as a storage medium. The CFD simulation used a porous flow model. Transient-state simulations were performed on a 2D axisymmetric model using a pressure-based solver. During charging, heat input that keeps the bottom temperature at 550 °C was applied to raise the storage temperature. Performance analysis was conducted under various porosities, considering natural and forced convection. The natural convection analysis showed insignificant convection contribution after 10 h of charging, as observed in both average air velocity and the temperature profile plots. The temperature distribution profiles at various positions for both convection modes showed good agreement between the simulation and experimental results. Additionally, both cases exhibited similar temperature growth trends, further validating the models. Forced convection reduced the charging time from 60 h to 5 h to store 70 MJ of energy at a porosity of 0.4, compared to natural convection, which stored only 50 MJ in the same time. This extended charging period was attributed to poor natural convective heat transfer, indicating that relying solely on natural convection for thermal energy storage under the given conditions is not practical. Using a small fan to enhance heat transfer, forced convection is a more practical method for charging the system, making it suitable for power generation applications.

Design, Fabrication, and Commissioning of Transonic Linear Cascade for Micro-Shock Wave Analysis

Understanding shock wave behavior in supersonic flow environments is critical for optimizing the aerodynamic performance of turbomachinery components. This study introduces a novel transonic linear cascade design, focusing on advanced blade manufacturing and experimental validation. Blades were 3D-printed using Inconel 625, enabling tight control over the geometry and surface quality, which were verified through extensive dimensional accuracy assessments and surface finish quality checks using coordinate measuring machines (CMMs). Numerical simulations were performed using Ansys CFX with an implicit pressure-based solver and high-order numerical schemes to accurately model the shock wave phenomena. To validate the simulations, experimental tests were conducted using Schlieren visualization, ensuring high fidelity in capturing the shock wave dynamics. A custom-designed test rig was commissioned to replicate the specific requirements of the cascade, enabling stable and repeatable testing conditions. Experiments were conducted at three different inlet pressures (0.7-bar, 0.8-bar, and 0.9-bar gauges) at a constant temperature of 21 °C. Results indicated that the shock wave intensity and position are highly sensitive to the inlet pressure, with higher pressures producing more intense and extensive shock waves. While the numerical simulations aligned broadly with the experimental observations, discrepancies at finer flow scales suggest the need for the further refinement of the computational models to capture detailed flow phenomena accurately.

Multiscale CFD analysis of urban air pollution dome and ventilation enhancement via an urban chimney

Air pollution episodes are critical environmental challenges in urban areas, primarily linked to the lack of ventilation during calm and stable atmospheric conditions, rather than a significant increase in emissions. The prevailing circulation pattern under these conditions is the Urban Heat Island Circulation (UHIC), which governs the characteristics of the urban pollution dome and the distribution of pollutants. This study aims to assess the UHIC's influence on spatial and temporal variations in pollutant concentrations and to investigate the efficiency of an Urban Chimney (UC) in improving air quality. Due to the interaction of micro- and meso-scale flows, a multi-scale analysis is required. For this purpose, a pressure-based CFD solver is adapted for analyzing airflow and pollutant dispersion in a multi-scale environment. A coordinate transformation method is used to account for meso-scale effects, and the resulting source terms are applied to the solver. The species transport equation is transformed to the new coordinate system, and the effects of species diffusion on enthalpy transport are included in the governing equations. The results demonstrate that UHIC reduces the average pollutant concentration and causes pollution peaks near ground level and at the city center. The interaction of micro-scale flow within the chimney with meso-scale UHIC impacts the flow field and pollutant concentration across the entire domain, decreasing the urban plume's vertical velocity by 40 percent and the mixing height at the city center by 20 percent. It also increases the vertical velocity within the UC by 23 percent. The UC can serve as a controlling tool for urban ventilation, and its capacity to remove pollutants increases sixfold when interacting with UHIC, reducing pollutant concentration citywide.

Effect of Process Parameters on the Substrate Surface in Cold Spray Coating

In this work, we used the cold spray process to coat a surface and investigated to understand the influence of various process parameters on the substrate surface temperature. The parameters we varied included the pressure, temperature, particle size, and particle speed of the titanium powder that we used in the cold spray coating process. It is a unique finding in the field of cold spray coating. For the substrate material, we chose steel and simulated the spray geometry using a two-dimensional axisymmetric model. This model employed a k-ɛ turbulence model with an implicit pressure-based solver of second-order precision. Our observations revealed that the surface temperature of the substrate reached its maximum value when the length of the injector was 15 mm. We also found that the most compatible length for the nozzle barrel was equal to the length of the particle injector.

Investigation on the Blade Number of Pico-scale Crossflow Turbine for Low Head by Numerical Method

Pico-scale crossflow turbines (CFT) can be an alternative solution to meet electrical energy needs, especially in remote rural areas. CFT is recommended because of its suitability in low head (< 5 m) conditions and fluctuating discharge conditions. One of the parameters that influences the performance of a CFT is the number of blades of the runner. CFT was discovered in 1903 and is still developing; however, the study of the physical phenomena of flow due to the blade number on the energy conversion process has yet to be comprehensively depicted. Therefore, this study aims to analyze the effect of the blade's number of runners on CFT performance using the computational fluid dynamics (CFD) method. The CFD method can visualize the flow field more detail than analytical and experimental. The CFD method is run with a moving mesh feature (transient) and pressure-based solver with a head condition of 3 m. The blades number studied were 16, 18, 22, 24, 26, and 30. Based on the results, the relationship of the CFT efficiency to blade number is described using a second-order multiple regression polynomial, and runner rotation is parabolic. Based on the performance curve, the CFT with 26 blades has the highest performance for low-head conditions.

Optimization of Flanged Diffuser for Small-Scale Wind Power Applications

The development of renewable and clean energy has become more crucial to societies due to the increasing energy demand and fast depletion of fossil fuels. A state-of-the-art design for an augmented wind turbine has been introduced in the past years to increase the efficiency of compact horizontal axis wind turbines, exceeding the ideal Betz’s limit of the maximum energy captured from the wind. The optimization of the flanged diffuser - so-called diffuser augmented wind turbine DAWT - is investigated numerically using the multi-objective genetic algorithm “MOGA”. A 2D computational model is developed using ICEM CFD and solved by ANSYS Fluent. The Turbulence model selected is shear stress transport K-omega, with a pressure-based solver and a coupled algorithm scheme. The optimization objectives are to maximize the velocity ratio at the shroud throat and minimize shroud form dimensions. 517 design points were solved, and the design dimensions were categorized into four types: compact, small, medium, and large design. The results showed that the diffuser dimensions are the main parameters to increase velocity inside the shroud throat, where a long diffuser with a low converging angle drags more air inside the shroud, reaching in some cases more than double the upwind velocity. While the nozzle and flange are also effective in the different design types. It was found that a super long diffuser with a length ratio of 2.9 LD to throat diameter D is optimal with a diverging angle of 7.6˚, accompanied by a nozzle of ratio 1.2 LN/D and 12.6˚ converging angle and a flange length ratio of 0.6 LF/D. This optimal design increased the velocity ratio by almost 2.5 times.

Validation of High Speed Reactive Flow Solver in OpenFOAM with Detailed Chemistry

An OpenFOAM® based hybrid-central solver called reactingPimpleCentralFoam is validated to compute hydrogen-based detonations. This solver is a pressure-based semi-implicit compressible flow solver based on central-upwind schemes of Kurganov and Tadmor. This solver possesses the features of standard OpenFOAM® solvers namely, rhoCentralFoam, reactingFoam and pimpleFoam. The solver utilizes Kurganov & Tadmor schemes for flux splitting to solve the high-speed compressible regimes with/without hydrodynamic discontinuity. In this work, we present the validation results that were obtained from one-dimensional (1D) and two-dimensional (2D) simulations with detailed chemistry. We consider three different mixtures that fall into the categories of weakly unstable mixture (2H2 +O2 +3.76Ar and 2H2 +O2 +10Ar), and moderately unstable mixture (2H2 +O2 +3.76N2), based on their approximate effective activation energy. We performed the numerical simulations in both laboratory frame of reference (LFR) and shock-attached frame of reference (SFR) for the 1D cases. The 1D simulation results obtained using this solver agree well with the steady-state calculations of Zel’dovich von Neumann Döring (ZND) simulations with an average error below 1% in all cases. For the 2D simulations, circular hot-spots were used to perturb the initially-planar detonations to develop into spatio-temporally unstable detonation front. The convergence is declared when the front does not deviate much from the CJ speed (Chapman-Jouguet) and the regularity of cellular pattern on the numerical smoke foils reaches a steady state. We have verified from our preliminary studies that the SFR-based simulations are computationally cheaper in comparison to the LFR simulations and that the required grid resolution is always lesser in the former than the latter to reach the same level of accuracy (in terms of speed of the detonation front and cell sizes from the numerical smoke foil). We have also verified that at least 24 points per induction zone length (for weakly unstable mixture) and 40 points per induction zone length (for moderately unstable mixture) are required to sufficiently resolve the detonation structures that are independent of grids, boundary and initial conditions. Further reduction in computational cost of approximately 50% is achieved by using non-uniform grids, which converge effectively to the same solutions in comparison to the results from twice the number of grids with uniform resolution.

A new semi-implicit pressure-based solver considering real gas effect

A new semi-implicit pressure-based solver considering the real gas effect is developed and tested. The advantage of the semi-implicit solver is that it can take a large time step, which makes the calculation speed faster than the classical density-based explicit solver. The comparison of the calculation speed between the explicit density-based solver and the new semi-implicit pressure-based solver is conducted by performing two-dimensional numerical simulations of a cryogenic nitrogen planar jet. Moreover, the developed pressure-based solver is validated by comparing the numerical results of a three-dimensional large-eddy simulation of a cryogenic nitrogen jet under the supercritical pressure condition with the experimental result. Results show that the developed pressure-based solver can calculate more than five times faster than the classical explicit density-based solver and accurately predict the jet behavior under the supercritical pressure condition.

Modelling of liquid oxygen and nitrogen injection under flashing conditions

The present numerical investigation of two-phase flashing flows examines the injection of liquid oxygen and liquid nitrogen into near-vacuum conditions prevailing in the upper-stage boosters of rocket engines. The predictive capability of a pressure-based solver and a density-based solver, each employing distinct approaches related to the imposed phase-change rate and thermodynamics closure, has been comparatively evaluated. Regarding the pressure-based solver, the departure from thermodynamic equilibrium during phase-change has been taken into account via the implementation of a bubble-dynamics model employing the Hertz-Knudsen equation. In contrast, the density-based solver relies on the adoption of thermodynamic equilibrium while real-fluid thermodynamic properties are assumed by loading tabulated values to the solver. Each thermodynamic property value was calculated in advance by solving the Helmholtz Equation of State (EoS) for a wide range of density and internal energy conditions. Numerical findings have been compared against experimental data available in the literature. The comparison demonstrates the capability of both methodologies in capturing the evolution of cryogenic flashing flow expansion, phase-change, and spray formation. The salient features identified in the numerical results, i.e., the expansion sphere immediately downstream of the injector exit, the bell-shaped topology of the spray, as well as the dependency of the spray cone angle on superheat, are in agreement with experimental measurements. Especially the density-based approach has been proven highly accurate with respect to the steady expanding flow described by a level of superheat in the range of 3 to 245, while also being independent of any parameter tuning.

Performance of OWC device under Stokes second-order waves

The present work simulates the flow behavior in the system having oscillating water column (OWC) device for wave energy conversion using numerical modelling. To achieve this, the OWC device is structured inside the 2-dimensional numerical wave tank (NWT). Using the ANSYS Fluent software and the multiphase-based CFD technique, the issue is resolved. Unsteady, transient, and pressure-based solvers are included in this work. In this numerical model, the governing equations are based on Reynolds Averaged Navier-Stokes (RANS) equations with the Shear Stress Transport (SST) model, to capture the turbulence effect and the interface of 2-phases are implicitly captured by using volume of fluid (VOF) method. Using an open-channel wave generator, this study generates nonlinear “second-order Stokes waves”. In addition, the impact of nonlinearity on the “free surface description” at specific locations, the vertex average of total pressure near the air-duct, and the vertex average of y-velocity through the air-duct per unit width of the OWC device are investigated in detail. The outcomes indicate that the OWC device will consistently extract energy over an extended period of time.

Investigation of the effects of volume change on flow structure and acoustic in a silencer

This study numerically examined the propellant flow from gunfire using the kω-SST turbulence model and their sound attenuation using the Ffowcs-Williams and Hawkins equations (FW-H). For simulation, a pressure-based solver and 3D axisymmetric geometry were used. The second-order implicit time approach and the second-order upwind scheme spatial discretization were used in the simulation. The maximum exit pressure was 3.748 MPa for the suppressor with a length of 70 mm and diameter of 20 mm. However, when the diameter suppressor increased by 1/6, the maximum exit pressure was reduced to 3.4961 MPa. When the length increased by 1/6, the maximum pressure became 3.3636 MPa. Lastly, when the diameter and length were increased by 1/6, the maximum exit pressure became 3.177 MPa. For this suppressor, 20.835 dB (12.29%) sound pressure level attenuation was achieved with 16.823 MPa (84.115%) overpressure reduction and 484.86 K or 32.32% temperature reduction. Generally, the attenuation increased with the increase in the suppressor’s internal volume.

Pressure-based large-eddy simulation of under-expanded hydrogen jets for engine applications

An assessment of Large-Eddy Simulations (LES) of non-reactive under-expanded hydrogen jets by using a pressure-based algorithm is presented. Such jets feature strong compressible discontinuities often considered to be best dealt with by a density-based solver. The crucial contribution of this work is to evaluate the suitability of the pressure-based solver to correctly describe the flow field of gaseous hydrogen jets for engine applications, despite the presence of shock waves in the under-expanded near-orifice region. Inherently, the paper aims at providing guidance on the corresponding numerical aspects to simulate these flows. Hydrogen jets in an argon atmosphere at three different injection pressures are simulated and the results are compared to experiments in literature. Jet tip penetration and cone angle are the main investigated parameters. A good match is found, confirming the solidity of the proposed model. Different LES sub-grid scale models and discretisation schemes are then investigated in order to find the best approach in terms of accuracy and required computational cost. In particular, it is found that the WALE model coupled with a 4th-order cubic scheme for the convective terms yields the most suitable configuration.

A super-grid approach for LES combustion closure using the Linear Eddy Model

LES–LEM is a simulation approach for turbulent combustion in which the stochastic Linear Eddy Model (LEM) is used for sub-grid mixing and combustion closure in Large-Eddy Simulation (LES). LEM resolves, along a one-dimensional line, all spatial and temporal scales, provides on-the-fly local turbulent flame statistics, captures finite rate chemistry effects and directly incorporates turbulence-chemistry interaction. However, the approach is computationally expensive as it requires advancing an LEM-line in each LES cell. This paper introduces a novel turbulent combustion closure model for LES using LEM to address this issue. It involves coarse-graining the LES mesh to generate a coarse- level ‘super-grid’ comprised of cell-clusters. Each cell-cluster, instead of each LES cell, then contains a single LEM domain. This domain advances the combined advection–reaction–diffusion solution and also provides suitably conditioned statistics for thermochemical scalars such as species mass fractions. Local LES-filtered thermochemical states are then obtained by probability-density-function (PDF) weighted integration of binned conditionally averaged scalars, akin to standard presumed PDF approaches for reactive LES but with physics-based determination of the full thermochemical state for particular values of the conditioning variables. The proposed method is termed ‘super-grid LEM’ or ‘SG-LEM’. The paper describes LEM reaction–diffusion advancement, the LEM representation of turbulent advection, a novel splicing algorithm (a key feature of LES–LEM) formulated for the super-grid approach, a wall treatment, and a thermochemical LES closure procedure. To validate the proposed model, a pressure-based solver was developed using the OpenFOAM library and tested on a premixed ethylene flame stabilised over a backward facing step, a setup for which some DNS data is available. SG-LEM provides high resolution flame structures, temperature and mass fractions suitable for LES thermochemical closure. Additionally, it provides reaction-rate data at the coarse level, a unique feature compared to other mapping-type closure methods. Quantitative comparisons are made between the proposed model and time-averaged DNS data, focussing on velocity, temperature and species mass fraction. Results show good agreement downstream of the step. Furthermore, comparison with an equivalent Partially-Stirred Reactor (PaSR) simulation demonstrates the superior predictive capability of SG-LEM. Additionally, the paper briefly examines the sensitivity of the model to coarse-graining parameters and finally, explores computational efficiency highlighting the substantial speedup achieved when compared to the standard LES–LEM approach with potentially significant speedup relative to PaSR closure for the intensely turbulent regimes of principal interest.

Numerical Study of the Torque and Power of a Hydraulic Turbine with Oscillating Blades

This paper presents results from a physical and numerical study of a new type of axial hydraulic turbine with oscillating blades, which is used to utilize wind waves energy. The pilot studies were conducted on a test bench constructed in one of the labs of the “Department of Hydrodynamics and Hydraulic Machines” in the Technical University of Sofia. The numerical computations were performed with the commercial software package Ansys Fluent 2022. The flow has been modeled with the k-ω (SST) turbulence model, whose main advantage is to resolve the viscous sublayer in over refined meshes. A pressure-based solver was used since the fluid is incompressible and the flow velocity is low. The study investigated several different pitch angles of the blades ranging between 0 and 80 deg at prescribed upstream flow velocities from 0.15 m/s to 2.0 m/s. The dependencies of the torque coefficient, power coefficient, and the optimal tip speed ratio on the flow velocity are presented and discussed.

Simulasi CFD Distribusi Temperatur pada Pengering Biji Kopi dengan Sistem Konveksi Paksa

Drying is a mass transfer and heat transfer process that is strongly influenced by temperature and air velocity. The importance of analyzing the temperature distribution and air velocity is expected to be used as a reference to optimize the drying process to be faster and better. This study aims to analyze the temperature distribution and air velocity that occurs in a coffee bean dryer with a forced convection system using CFD simulation. CFD simulations are carried out using ANSYS Fluent software. with steady conditions, and using a pressure-based solver method. The obtained temperature distribution occurs with an estimated absorber temperature difference with the drying chamber is 31% during the day.

3D simulations and laboratory experiments to evaluate a dynamic airbag valve

Airbag pressure determines the restraint effect during a vehicle crash. The pressure required to restrain the occupant depends on pre-crash detection, collision parameters and the occupant’s mass and position. This work modulated airbag pressure for optimum safety using a novel airbag control valve for cold-gas inflators. This paper evaluates the valve’s stationary and dynamic performances for Helium by 3D flow simulations using a pressure-based solver in ANSYS Fluent® and SAE J2238 laboratory tank tests. The predicted and measured tank pressures for the fully open (stationary) valve were agreed by an average 93.73% with an excellent correlation (correlation coefficient, R = 0.9995). For the first dynamic operation with 10 ms switching time, the results agreed by 92.78% with R = 0.9975. In the second test with 30 ms switching, 83.67% agreement was observed with R = 0.9893. The research concluded that the valve modulates the bag pressure and is implementable in vehicles.

Impact of fin and hybrid nanofluid on hydraulic-thermal performance and entropy generation in a solar collector using a two-phase approach

This research investigates the impact of cylindrical fins and a hybrid nanofluid (HNF) composed of water and copper-zinc oxide on the hydraulic-thermal performance and entropy generation of a solar collector (SC). Numerical simulations are conducted for steady flow utilizing the pressure-based solver. The HNF flow is modeled using the two-phase mixture approach, while the turbulent flow is represented using the standard k-ε turbulence model. The study explores varying inlet Reynolds numbers (Re) from 17,000 to 41,000, volume fractions of nanoparticles (φ) ranging between 1 and 3%, and cylindrical fin heights (λ) of 5, 10, and 15 mm. The results indicate that employing twisted tape in the SC leads to improved thermo-hydraulic performance compared to a SC without twisted tape. Also, the addition of cylindrical fins plays an important role in the flow field because the HNF flow separates after colliding with the fins. This causes the formation of swirling vortices and spatial and temporal fluctuations. The amount of the performance evaluation criterion (PEC) is similar in all cases and decreases with Re. The PEC is much higher than the desired one. Although this trend is decreasing, the results show that in this range of Re, the use of cylindrical fins creates a more reasonable pressure drop (Δp) compared to the enhancement in thermal performance. Finally, the thermal and total entropy generation are reduced and frictional entropy is enhanced by increasing Re. In other words, the increment in the height of the cylindrical fins in the solar system leads to a reduction in the entropy generation.

Hydrodynamics of shear thinning fluid in a square microchannel: a numerical approach

Abstract In this present work, a numerical study was conducted for the formation of a slug bubble for shear thinning non-Newtonian fluid in a cross-junction 2-D square horizontal microchannel. Carboxymethyl cellulose (CMC) of concentration 0.2 (w/w%) percent was used as a continuous phase that shows the shear thinning behavior of non-Newtonian fluid and Nitrogen (N2) was used as the discrete phase. The pressure-based double precision solver was used in ANSYS FLUENT 2021 R2 with the volume of fluid (VOF) method. The finite volume method is applied for the discretization of the continuity and momentum equation. This article also focuses on the fluctuation of static pressure, mechanism of slug, annular, and churn annular flow i.e., obtained by the variation in the inlet velocities. On the other hand, a concept that was applied in this work was also validated with the prior literature data.

Improving performance of H-Type NACA 0021 Darrieus rotor using leading-edge stationary/rotating microcylinders: Numerical studies

In the current paper, a computational study has been performed to simulate a novel concept in which a fixed/rotating microcylinder is installed upstream of the leading edge of the NACA 0021 airfoil, aiming to enhance the performance of the H-type Darrieus rotor. The idea behind implementing a fixed/rotating microcylinder is to enhance the lift to drag ratio due to increasing near-wall momentum through generated vortices transferring the momentum from the outer flow to the near-wall flow region of the airfoil. Parametric analyses for the microcylinder including its location, size, shape (circular, square and rhombus), and rotation were performed. The commercial Computational Fluid Dynamics (CFD) package ANSYS Fluent was used for solving the Unsteady Reynolds Averaged Navier-Stokes (URANS), and turbulence equations. The code has been based on Finite Volume Method (FVM) with the pressure-based solver developed for low-speed incompressible flows. The semi-implicit method for pressure-linked equation (SIMPLE) was used to solve the discretized equations. The URANS was used with adopting the SST k-ω turbulence model to find out the optimum rotor by utilizing the microcylinder model. An optimization methodology using Response Surface Optimization (RSM) based on Kriging method has been first performed to find the optimum size and location of a circular microcylinder. The optimization study indicated an optimum microcylinder diameter (d/C = 0.0085313), at an upstream chordwise distance and normal to it 0.070275 and 0.02303 of the chord length (C), respectively. The results showed also that a small static circular microcylinder of d/C = 0.009 and installed at 0.075C from the leading edge (MC5) is an efficient geometry at high regime of tip-speed ratios (TSR≥2.2) rather than lower regime of TSRs. When these superior set of rotating (5 rad/s) microcylinder (MC5) parameters were adopted, significant enhancement of power coefficient (Cp) could be achieved and a considerable improvement of the maximum power coefficient (Cpmax) up to 120 % at TSR = 3 appeared. A physical analysis of the flow fields using the contours of vorticity, pressure and turbulence energy has been performed to illustrate the effect of rotating microcylinder for improving rotor performance. The vorticity contours showed the ability of rotating microcylinder for generating a strong vortex structure in its wake for enhancing turbulence production around the blades, which delayed the known strong stall of the blades by diminishing the separation bubble size. This generally improved the blades aerodynamic performance and enhancing the lift-to-drag ratio.