Rotational Reynolds Number Research Articles

Measurement and Analysis of Buoyancy-Induced Heat Transfer in Aero-Engine Compressor Rotors

Abstract The flow inside cavities between corotating compressor disks of aero-engines is driven by buoyancy, with Grashof numbers exceeding 1013. This phenomenon creates a conjugate problem: the Nusselt numbers depend on the radial temperature distribution of the disks, and the disk temperatures depend on the Nusselt numbers. Furthermore, Coriolis forces in the rotating fluid generate cyclonic and anticyclonic circulations inside the cavity. Such flows are three-dimensional, unsteady, and unstable, and it is a challenge to compute and measure the heat transfer from the disks to the axial throughflow in the compressor. In this paper, Nusselt numbers are experimentally determined from measurements of steady-state temperatures on the surfaces of both disks in a rotating cavity of the Bath compressor-cavity rig. The data are collected over a range of engine-representative parameters and are the first results from a new experimental facility specifically designed to investigate buoyancy-induced flow. The radial distributions of disk temperature were collected under carefully controlled thermal boundary conditions appropriate for analysis using a Bayesian model combined with the equations for a circular fin. The Owen-Tang buoyancy model has been used to compare predicted radial distributions of disk temperatures and Nusselt numbers with some of the experimentally determined values, taking account of radiation between the interior surfaces of the cavity. The experiments show that the average Nusselt numbers on the disk increase as the buoyancy forces increase. At high rotational speeds, the temperature rise in the core, created by compressibility effects in the air, attenuates the heat transfer, and there is a critical rotational Reynolds number for which the Nusselt number is a maximum. In the cavity, there is an inner region dominated by forced convection and an outer region dominated by buoyancy-induced flow. The inner region is a mixing region, in which entrained cold throughflow encounters hot flow from the Ekman layers on the disks. Consequently, the Nusselt numbers on the downstream disk in the inner region tend to be higher than those on the upstream disk.

Thermoelectric generation from vented cavities with a rotating conic object and highly conductive CNT nanofluids for renewable energy systems

In the present work, thermoelectric power generation from cavities with ventilation ports is considered by using a rotating conic object and carbon-nanotube particles in the base fluid. Effects of different pertinent parameters such as Reynolds numbers of hot and cold fluid streams (between 200 and 1000), rotational Reynolds number of the conic object (between −400 and 400), size (between 0.05H and 0.25H) and horizontal location (between 0.2H and 0.6H) of the object and nanoparticle volume fractions of nanoparticles (between 0 and 0.02) on fluid flow, interface temperature and generated thermoelectric output power characteristics were studied. It was observed that exit port location of the hot cavity, fluid stream Reynolds number and rotational Reynolds number of the object have significant impacts on the fluid flow, interface temperatures and output power. When lowest and highest fluid stream Reynolds number is compared, 43.75% variations in the output power is obtained. The clockwise rotation of the conic object results in higher thermoelectric power generations as compared to a stationary cone while up to 49.20% increments of the power is attained at the highest rotational speed. The size of the rotating object is influential on the power generation while its horizontal location has slight effects. When nanofluid at the highest solid volume fraction configuration is compared with pure water case, 20% increment in the thermoelectric power is obtained. A predictive model for power estimations with adaptive network-based fuzzy inference system is proposed for input parameters of hot and cold fluid stream Reynolds number and rotational Reynolds number of the conic object which delivers accurate and fast prediction results.

Numerical study of turbulent mixed convection in a rotating inter-disk cavity with axial throughflow of cooling air

The paper presents the results of eddy-resolving numerical simulation of mixed convection developing in a rapidly rotating inter-disk cavity with axial throughflow of cooling air. Under conditions of the experiments available in the literature for the case of the cavity heating from the side of the disk surfaces, the computations were performed using two methods: Implicit LES and zonal RANS/ILES, with and without taking into account the gravity effects. In case of the zonal RANS/ILES approach, the axial pipe flow and the zone of turbulent mixing of throughflow with particles of heated gas were simulated on the base of the RANS equations added by a turbulence model, whereas the major part of the buoyancy-controlled cavity flow was computed with the ILES method. The ratio of the cavity width to its outer radius was 0.13. The computations were carried out at the rotational Reynolds number equal to 2·105, if evaluated with the outer radius; the axial flow Reynolds number was 2·104. A moderate size grid used in the computations (of about 5·105 cells) provided a good resolution of quasi-laminar Ekman layers. Data on dynamics of unsteady vortex structures forming in the cavity, on characteristics of local heat transfer, as well as on effects of the gravity force action are reported. The results of computations of the time-averaged local Nusselt number on the disks surfaces are compared with measurement data; an acceptable agreement has been obtained in case of the zonal RANS/ILES approach.

Heat Transfer Characteristics outside the Condenser of a Rotating Heat Pipe Grinding Wheel with a Lateral Air Impinging Jet

This study numerically analyzed the heat transfer characteristics outside the condenser of a rotating heat pipe grinding wheel (RHP-GW). The goal of this investigation is to determine the optimal structure and parameters for the condenser section of RHP-GW. Different fin height (f=0–7 mm), rotational Reynolds number (Rer=1602–6408) and jet Reynolds number (Rej=42 379–108 302) were analyzed under input heat flux of 4000 W/m2. A fully developed flow was imposed at the outlet of the nozzles. Results showed that the optimal heat transfer rate was obtained for fin height of 5 mm, which improved the average Nusselt number by 84% compared to the structure without fins. A critical Rej for each Rer that the impinging jet can reach the condenser section was found. The critical Rej value increases with Rer, which is in the range from 42 379 to 61 215 and 61 215 to 80 050 for Rer=6408 and Rer=9610, respectively.

Mixed convection between rotating sphere and concentric cubical enclosure

This paper presents numerical investigations of the mix convection between a rotating inner sphere and a concentric cubical enclosure, using the recently developed immersed boundary-simplified lattice Boltzmann method. The validity of the method has been established through benchmark tests, and a grid independence study is also carried out to ensure the accuracy of the conveyed results. Various factors that may influence the mixed convection system, such as the rotational direction, Rayleigh number, and rotational Reynolds number, are considered in the present study. Three representative rotational axes, namely, the vertical, the horizontal, and the diagonal axes, are selected. The Rayleigh number spans from 104 to 106, which covers the transition range from a conduction-dominated system to a convection-dominated one. Moreover, the rotational Reynolds number varies from 0 to 300. Distinct flow patterns, global heat transfer behavior, and the heat transfer rates on cubic walls are studied to reveal the characteristics of this problem. It is found that the rotationality of the inner sphere would increase the global heat transfer rate of the system and the rotating motion would stimulate heat transfer in the radial direction of the rotational axis while suppressing the thermal expansion in the axial direction.

Mixed convection flow of hybrid nanofluid through a vented enclosure with an inner rotating cylinder

A numerical solution for the influence of an inner adiabatic rotating cylinder inside a vented cavity on mixed convection of hybrid nanofluid is presented in this study. The governing equations of mixed convection flow for an incompressible Newtonian hybrid nanofluid are assumed to be two-dimensional, steady, and laminar. These equations are solved numerically by using the finite volume technique. In this work, (Al2O3 –Cu/ water) is introduced as a hybrid nanofluid. The influences of relevant parameters including nanoparticle concentrations (0 ≤ φ1 ≤ 0.02 and 0 ≤ φ2 ≤ 0.02), cylinder radiuses (0.1 ≤ R ≤ 0.3), cylinder locations (0.25 ≤ δ ≤ 0.75), angular rotational speeds (−5 ≤ Ω ≤ 5), Grashof numbers (103 ≤ Gr ≤ 105) and Reynolds numbers (50 ≤ Re ≤ 500) are examined. The obtained data are exhibited by means of different qualitative and quantitative figures. The results show that the energy transport of hybrid nanofluid enhances with an increase in solid particle concentration, but it is associated with an increased pressure drop. With respect to the stationary cylinder, it is reported that the counter-clockwise cylinder rotation increases convective heat transfer whereas the clockwise direction has a reverse effect. In addition, for the counter-clockwise rotation of the cylinder, the heat transfer enhances when the position of the cylinder approaches the hot wall. The maximum heat transfer enhancement over the stationary cylinder with 21% occurs in case of R = 0.3 and δ = 0.5 under Gr = 104, φ1 = φ2 = 0.02, Re = 100 and Ω =5. The current study can provide useful guidelines to the designers of rotary heat exchangers.

Influence of the hydrodynamics and the reaction-rate formulation in modeling infinitely fast irreversible reactions in a turbulent un-baffled chemical reactor

Numerical results are presented from large eddy simulations (LES) of a turbulent flow in an un-baffled stirred vortex reactor. The open source code TrioCFD is employed with a discontinuous front-tracking algorithm to detect the free surface separating the two fluids in the reactor. The influence of increasing the rotation Reynolds number on the flow hydrodynamics is investigated; a very good agreement is reported between the LES and both theoretical/experimental data. Infinitely fast and irreversible chemical reactions taking place in the liquid phase are modeled and different reaction models are examined. It is figured out that the chemical reactions are highly influenced by both hydrodynamics and the reaction model. Two macro-mixing zones are identified in the reactor where the first is in the center due to the downward-oriented forced vortex, while the second is in the free helicoidal vortex that forms about 96% of the liquid-phase volume. Numerical results show that only reaction models based on an Eddy-Dissipation Concept (EDC) are capable to reproduce the effect of turbulence in enhancing mixing, unlike classical simple-rate law reaction models that hold no information concerning the micro-mixing. Such observations demonstrate the necessity of using the EDC-based models as far as turbulence and mixing enhancement are two important phenomena encountered in the precipitation process of un-baffled mixing reactors. This study is a first mandatory step preceding the description of crystals formation due to precipitation which takes place in a wide range of energy related applications or in waste management.

Tuning the Morphology of Protein-Inorganic Calcium-Phosphate Supraparticles via Directed Assembly.

Understanding the phenomena that govern complex interfacial and directed assemblies is essential for both control and scale-up of particle syntheses. The present work describes an effort to understand, control, and tune the formation of protein-inorganic calcium-phosphate supraparticles that are produced at an oscillating air-water interface created by end-over-end rotation of the synthesis solution. Supraparticles were synthesized under an array of different conditions that varied reagent concentration, the presence of additives, tube size, and rotational speed. Paired with a fluid mechanics model of the end-over-end rotation and dimensional analysis, the sensitivity of the synthesis to physicochemical and mechanical parameters was determined. Surface tension and bubble formation were found to be important criteria for changing the size distribution of supraparticles. Thresholds for the values of the Froude, Iribarren, and rotational Reynolds numbers were identified for narrowing particle size distribution. These results both guide the specific protein-inorganic supraparticle synthesis described here and inform future manipulation and scale-up of other complex interfacial colloidal assemblies.

The cooling performance of three-row compound angle holes on the suction surface of a rotating turbine blade

Experiments were carried out to study the effects of compound angle, hole arrangement, and blowing ratio on the film cooling performance of multiple rows of holes on the suction surface of a turbine blade. The turbine worked at rotational speed of 600 rpm corresponding to the rotational Reynolds number of 5.36 × 105. Three rows of cylindrical holes arranged in line or in stagger were drilled on the rotor blade suction surface at the streamwise location of 12.4%, 17.8%, and 23.2%, respectively. Three compound angles, with the same streamwise angle of 45° but different lateral deflection angles of 45°, 0°, and −45°, were studied. The film cooling effectiveness was obtained using pressure sensitive paint (PSP) technique with average blowing ratios varied from 0.5 to 2.0. The results showed that the application of compound angle changes the jet direction in the near-hole region and makes the film spread laterally. Compared with the film cooling without compound angle, using positive and negative compound angle can improve overall average film cooling effectiveness by about 20% and 25%, respectively. The effects of the secondary flow also can be weakened. A stagger film trajectory arrangement can achieve more uniform film coverage with higher overall film cooling effectiveness. The film trajectory arrangement of a positive compound angle injection is determined by the combined effect of hole arrangement and blowing ratio. While, the film trajectory arrangement of a negative compound angle injection is almost the same as the hole arrangement and nearly does not change with the blowing ratio.

An Investigation on Forced Convection in Clearance of Corotating Disks With Diverging Laminar Flow

Abstract In this work, forced convection in the laminar outflow between two disks is predicted. The disks can be corotating or contrarotating with equal or with different angular velocities. A published closed form solution is used to model the forced convection. The simplified energy equation is discretized using a finite difference method. The discretized equation is solved using the tridiagnal matrix algorithm (TDMA). Results are obtained for a heat flux as well as a specified temperature boundary condition and are presented for various flow parameters such as throughflow Reynolds number, rotational Reynolds number, gap ratio and speed ratio. The parameters except the speed ratio affect the temperature field strongly. The computed Nusselt number is in good agreement with the published experimental data. The coupling of closed form solution and energy equation predicts temperature field accurately. Additionally, the solution procedure is simplified considerably and the solution is obtained with lesser computational resources.

Numerical investigation of flow and heat transfer between concentric cylinders with slit wall

This paper investigates the flow field distribution and heat transfer characteristics of turbulent flow between concentric cylinders with slit wall in the outer cylinder. The slit wall cross sections are triangle, rectangle and trapezoid respectively. The governing equations are established by the PANS model which is modified from the standard k - ε model. The rotating Reynolds number ranges from 2336 to 4682 and axial Reynolds number ranges from 161 to 645, the influence of slit wall cross section on the formation and intensity of small vortex in slit region is analyzed. The results show that the slit wall geometry has a significant influence on the flow field distribution and its heat transfer process, the main vortex in annular region and small scale vortex in slit region has a tendency to merge together in the model with trapezoid cross section, resulting in more intense mixing of fluid with different temperature, therefore, the slit model with trapezoid cross section shows the best heat transfer performance among three models. The present analysis results may be applied for the structural design and optimization of many related rotating machinery.

Empirical Modelling of Hydrodynamic Effects on Starch Nanoparticles Precipitation in a Spinning Disc Reactor.

Empirical correlations have been developed to relate experimentally determined starch nanoparticle size obtained in a solvent–antisolvent precipitation process with key hydrodynamic parameters of a spinning disc reactor (SDR). Three different combinations of dimensionless groups including a conventional Reynolds number (Re), rotational Reynolds number (Reω) and Rossby number (Ro) have been applied in individual models for two disc surfaces (smooth and grooved) to represent operating variables affecting film flow such as liquid flowrate and disc rotational speed, whilst initial supersaturation (S) has been included to represent varying antisolvent concentrations. Model 1 featuring a combination of Re, Reω and S shows good agreement with the experimental data for both the grooved and smooth discs. For the grooved disc, Re has a greater impact on particle size, whereas Reω is more influential on the smooth disc surface, the difference likely being due to the passive mixing induced by the grooves irrespective of the magnitude of the disc speed. Supersaturation has little impact on particle size within the limited initial supersaturation range studied. Model 2 which characterises both flow rate and disc rotational speed through Ro alone and combined with Re was less accurate in predicting particle size due to several inherent limitations.

Numerical Analysis for the Assessment of Factors Influencing the Breakdown of Swirl Flow in a Cylinder Driven by a Rotating End Wall

Abstract Finite element solution for the classical problem of swirl flow in a cylinder with a rotating lid has been used to study the characteristic features of the stream-tube and identify the factors contributing to axial vortex breakdown. An increase of rotational Reynolds number has been found to result in (i) a decrease of total flow rate; (ii) an increase of flow rate through the boundary layer over the stationary walls; (iii) an increase of the throat area of the stream-tube, with the upstream axial vortex flow in some cases having a deficit in momentum flux needed to overcome the pressure and viscous forces; and (iv) an increase of distance for the axial flow to sustain deceleration in the diverging passage. Based on the analysis, it is hypothesized that “flow with particles in axial vortex motion having a deficit of momentum flux for axial flow when subjecting to a fluctuating radial force undergoes axial vortex breakdown.” This explanation has been able to justify the disappearance of vortex breakdown at larger Re of laminar regime and the absence of vortex breakdown in small aspect ratio cylinders. We report novel results pertaining to total flow rate and its distribution within the vessel. The momentum flux of axial vortex, a main determinant of bubble breakdown, is found to be governed by the total flow rate, distribution of flow through the boundary layers, and the Reynolds number. The proposed hypothesis has been verified by analyzing two cases, one involving a passive and the other involving an active mechanism for regulating the axial momentum.

Experimental Study of Heat Transfer Through a Rotating Heat Exchanger.(Dept.M)

This paper presents an experimental study of laminar forced convection heat transfer for hot water flows through annular tube. The test section consisted of two annular spaces the inner of which for the tested hot water while the outer for the cooling water flow. The inner tube of the annulus is rotated with rotational speeds which varied from 350 to 1580 r.p.m. to give rotational Reynolds number in the range (200 ≤ Rew ≤ 1750). The outer tube of the annulus have a constant diameter of 54.5mm but the diameter of inner tube of the annulus is equal to 21, 26.5 and 33.5 mm so as to give diameter ratios of 0.62, 0.487 and 0.386. Hot water flowed axially through the annular space with velocities that ranged from (0.0042 to 0.52m/s), to give axial Reynolds numbers in the range (130 ≤ Re ≤ 2300), to cover the laminar flow regime. The cold water that flowed through the outer annulus was cooled by a refrigerating circuit. In the present work the effects of radius ratio, axial Reynolds number and rotational Reynolds number on the heat transfer were investigated. The results show that. The rate of heat transfer of the rotational annular tubes is higher than that of the stationary one of the same radius ratio and axial Reynolds numbers. An increase of as much as 200% in the heat transfer is observed for annular tube of radius ratio equal to 0.62, Reynolds number of 1951 rotational Reynolds number of 1744 and rotational speed of 1580 r.p.m respectively, a decreases in heat transfer as much as 36% in the heat transfer is reported for annular tube of radius ratio (B) of 0.386, Reynolds number of 1100, and rotational Reynolds number of 205 and rotational speed of 350 r.p.m respectively. Two correlations were established between Nusselt number and axial Reynolds number. Rotational Reynolds number and radius ratio of the annulus for both stationary and rotational annular tubes.

Numerical investigation on the flow and heat transfer for a millimeter-scale rotor-stator system with ultrahigh rotational speed

Micro turbine engine (MTE) is a promising device for the future energy application. However, the literature to date is limited in the research on the millimeter-scale rotor-stator system with the ultrahigh rotational speed in MTE. This work investigates the flow and heat transfer for a rotor-stator system with the radius of 4 mm and the rotational speed from 0.2 × 106 to 1.0 × 106 rpm. The mechanisms of centrifugal force and centrifugal buoyancy force are analyzed at the rotational Reynolds number within 0.2 × 105 and 1.0 × 105 and the rotor excess temperature of 25 and 275 K. The effect of geometric configuration including the spacing and opening is also discussed. The results demonstrate that the difference between the millimeter-scale and conventional rotor-stator system lies in not the flow structure, but the heat transfer, which is due to the viscous dissipation induced by the ultrahigh rotational speed. With the increasing rotational Reynolds number, the larger centrifugal force enhances the radial flow and raises the heat transfer. Meanwhile, the increasing viscous dissipation results in the heat transfer deterioration at high radius, especially at low rotor excess temperature. The large centrifugal buoyancy force induced by the high rotor excess temperature reduces the heat transfer by suppressing the radial flow, though it weakens the viscous dissipation. The spacing and opening determine the flow structure, the increase of which enhance the average heat transfer especially at low gap ratios. These findings can provide the reference for the design of millimeter-scale rotor-stator system in MTE.

Effects of a rotating tube bundle on the hydrothermal performance for forced convection in a vented cavity with Ag–MgO/water hybrid and CNT–water nanofluids

In this study, effects of rotating tube bundle on the hydrothermal performance for the forced convection in a vented cavity are numerically studied with a mixture of Ag and MgO nanoparticles suspended in water, forming a hybrid nanofluid and CNT–water nanofluid. The finite volume method is used for the numerical simulations. The simulation study is performed for various values of Reynolds number (between 100 and 1000), rotational Reynolds number (between − 100 and 100), size (between 0.01H and 0.12H), vertical (between 0.3H and 0.65H) and horizontal locations (between 0.3H and 0.65H) of the cylinder in the tube bundle and nanoparticle volume fraction (between 0 and 0.02). An experimental correlation was used for the description of the Ag–MgO/water hybrid nanofluid. Results indicated that the rotational effects of the tube bundle contribute to the hydrothermal performance. The average Nusselt number enhancements are 8.95% and 40% when CNT–water is used for solid volume fraction of 0.02 at the highest speed, while these values are 31.5% and 37.7% for hybrid nanofluid. The pressure coefficients are higher for negative value of $${\text {Re}}_{\upomega }$$ and hybrid nanofluid. The size of the tube bundle is effective at the highest size, while the increment amount is 16.8% both for CNT–water and for hybrid nanofluids. The center location of the tube bundle resulted in heat transfer enhancement which is in the range of 11% and 18%. As the solid volume fraction of nanoparticles becomes higher, discrepancies between the average Nu values of CNT–water nanofluid and Ag–MgO/water hybrid nanofluid become higher which is also observed for the pressure coefficient. It is advantageous to use CNT particles since higher heat transfer coefficients with lower values of pressure coefficient are obtained.

The cooling performance of multiple rows of film holes on the suction surface of a turbine blade under rotating conditions

In this paper, the experimental results on the cooling performance of multiple rows of holes on the suction surface of a turbine blade were indicated. The effects of hole arrangement, blowing ratio, and rotation on the film cooling performance were analyzed. Four different film hole arrangements were studied. Each blade was modeled with three rows of cylindrical holes with streamwise angle of 45°. The first and the second injection rows were fixed at the position of 12.4% and 17.8% of the streamwise length of the suction surface, respectively. The third injection row was located at streamwise location of either 23.2% or 34%. To investigate the rotation effect, experiments were performed at three rotational speeds of 400 rpm, 600 rpm, and 800 rpm corresponding to the rotational Reynolds number of 3.57 × 105, 5.36 × 105, and 7.15 × 105, respectively. Average blowing ratios of 0.5, 0.75, 1.0, 1.25, 1.5, and 2.0 were tested. Adiabatic film cooling effectiveness was obtained using pressure sensitive paint (PSP) technique. The results showed that inward deflection for the trajectories complicates the superposition principle. Staggered rows of holes can achieve more uniform film coverage and higher film cooling effectiveness. For the case with small array spacing between each adjacent rows, the film cooling area is too concentrated and can achieve better film cooling performance only at a small blowing ratio. The increasing rotational speed slightly weakens the tendency of inward film deflection and increases the film cooling effectiveness at a low blowing ratio.

Investigation on the Radial Heat-Transfer of Pumping Secondary Flow in Semiclosed Rotating Disk Cavity

Abstract This study investigates the mean flow and radial heat-transfer behaviors in semiclosed rotating disk cavity within the canned reactor coolant pump. The flow in the semiclosed cavity contains the Stewartson type flow at inner region and the Batchelor type flow at outer region. The heat is radially transported from the outer rim of the semiclosed disk cavity to discharge-hole through the nondirect discharge (ND) portion of the superimposed flow from inlet. The effects of rotating Reynolds numbers, cavity aspect ratio and radial location of discharge-hole on the discharge ratio, pumping mass flow rate, local wall shear stress and radial heat-transfer coefficient are examined in the semiclosed rotating cavity flow, respectively. Based on the radial heat transfer behaviors of pumping secondary flow, an equivalent thermal network is proposed and validated by experiments, which can effectively predict the radial temperature distribution from the discharge hole to periphery with the viscous-heating and nonisothermal effects.

Flow and heat transfer in a rotating cavity with de-swirl nozzles: An LES study

Large eddy simulation was applied for studying the flow and heat transfer characteristics in the rotating cavity with de-swirl nozzles. Simulations were performed at the rotational Reynolds number of 1.4 × 105 and nondimensional flow rate of 4040. Compared with vertical nozzle, de-swirl nozzle generates lower swirl ratio, pressure loss and Nusselt number. Distributions of radial velocity show that de-swirl nozzle generates lower velocity in the Ekman boundary layer, but higher velocity in the core region. In the corner region of shroud and jet, the ‘dead zone’ of coherent vortices exists for vertical nozzle. However, for de-swirl nozzle, the ‘dead zone’ disappears due to the impact of inclined injection, and accordingly, the scale of time-averaged counter-rotating vortex pair is reduced effectively. In the middle and lower sections of cavity, plenty of vortex stripes are formed both for vertical and de-swirl nozzles; however, in the Ekman layer, turbulent kinetic energy and Reynolds normal stress for vertical nozzle is much higher than de-swirl nozzle. The distributions of turbulent energy spectrum and energy integral length scale show that, for de-swirl nozzle, large-scale eddies take the dominant role in middle section, while small-scale eddies dominate the flow fields in upper and lower sections of cavity.

Performance assessment of a thermoelectric module by using rotating circular cylinders and nanofluids in the channel flow for renewable energy applications

In the present study, a novel thermal flow system configuration is offered with rotating cylinders and nanofluids in channel to increase the powers generated in the TEG module mounted between the channels. The 3D coupled heat transfer, fluid flow and electric field equations are solved with finite element method. The numerical study is performed for different Reynolds number (between 350 and 1000), rotational Reynolds number (between −4000 and 4000), number (between 1 and 5), size (between 0 and 0.2H), location of the circular cylinders and solid CNT nanoparticle volume fraction (between 0 and 0.02). It is observed that as the cylinders are used in the channel with rotation, power enhancements of TEG module are achieved. When reference case of no-cylinder is compared, power increments of 56.5% and 27.5% are achieved for rotational Reynolds number of −3000 at Reynolds number of 350 and 1000. Higher number of cylinders provide higher power enhancements in the module and with five cylinders rotating at the highest speed, increments up to 47% are obtained at Reynolds number of 1000. Size and location of the cylinders are also found o be effective in altering the power features of the TEG module. Inclusion of nanoparticle is found to be efficient in the TEG power enhancement when performance of the system is low. Thus, slight improvement is achieved with nanoparticles when system with five cylinder at the highest speed and highest Reynolds number is considered. A reduced order model which is based on proper orthogonal decomposition methods is used to predict the temperature variations in the coupled domains and to effectively calculate the power generated in the TEG module.