• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    Numerical Study on Flow and Heat Transfer Characteristics Inside Rotating Rib-Roughened Pipe with Axial Throughflow

    2023-11-22 09:11:20,,,

    ,,,

    College of Energy and Power, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, P.R.China

    Abstract: A numerical study was performed on the flow and heat transfer characteristics inside a rotating ribroughened pipe with an axial throughflow, under a fixed axial throughflow Reynolds number Rex 6 400 and a series of rotational Reynolds number Reω ranging from 0 to 3.77×104.Three rib configurations were taken into consideration,including ring rib, discrete-ring rib and longitudinal rib, wherein the ring-type ribs were distributed along the axial direction and the longitudinal-type ribs were distributed along the circumferential direction.The results show that the flow field inside a rotating pipe is affected by both axial through flow and rotation-induced azimuthal flow.Under small rotational Reynolds numbers, the axial through flow is dominant such that the ring-type rib plays significant heat transfer enhancement role.Whereas under high rotational Reynolds numbers, the rotation-induced azimuthal flow is dominant such that the longitudinal-type rib is superior to the ring-type rib on heat transfer enhancement.In the view of a comprehensive performance factor that considers both heat transfer enhancement and flow loss together, the discrete-ring rib is identified to be the best among the current three rib configurations under small rotational Reynolds numbers.However, under high rotational Reynolds numbers, the ring-type ribs show no advantage over smooth pipe, and even lead to a reduction in comprehensive performance when compared to the smooth pipe.For the highspeed rotating pipe, the longitudinal-type rib is a promising rib configuration as applied to the heat transfer enhancement.

    Key words:rotating pipe; rib configuration; axial throughflow; flow and heat transfer; comprehensive performance factor

    0 Introduction

    The hydrodynamic gas bearing utilizes the gas dynamic pressure compression effect to provide rigid rotor bearing capacity.As it owns some inherent advantages, such as self-acting, no requirements of external pressurization and lubricant supply, and low-cost maintenance, the hydrodynamic gas bearing has been widely applied in the airborne electricmachine system and turbomachineary[1].To satisfy the technological advancement of airborne electricmachine system and turbomachineary towards more compact size and greater efficiency, the rotational speed needs further prompted such that the thermal effects on the performance of foil journal bearings will significantly behave in the practical applications, due to the strong viscous-shearing inside the thin hydrodynamic gas film[2-4].Apparently the thermal management becomes one of the quite challenging subjects in the ultra-high speed journal bearings.

    Considering that the natural cooling scheme to remove the heat generation inside the hydrodynamic gas film is crucially limited, the use of additional forced cooling is a significant necessity for avoiding impervious degradation and even thermal damage of ultra-high speed journal bearings[5].During past decades, a lot of efforts were devoted to the thermal management of hydrodynamic foil journal bearings.Refs.[6-7] proposed two cooling schemes for the radial gas foil bearings, such as the axial throughflow cooling scheme by forcing the cooling air through the film layer axially and air injection cooling scheme by supplying the cooling air into the film layer radially from feeding holes.Ref.[8] proposed a hybrid cooling scheme, wherein two cooling air flow paths were designed, referring as the outer cooling flow that flows through the gas film and foil gap, and the inner cooling flow that flows through the hollow shaft.Ref.[9] performed laboratory tests to identify the aerodynamics of foil bearing under hot surroundings with the use of axial cooling flow.Ref.[10] performed an investigation on the design of the cooling flow passage for the gas foil bearings used in the micro gas turbines.Ref.[11] performed an experimental and theoretical investigation on the thermal behavior of a hybrid bump-metal foil bearing.Different thermal managements with cooling air flows in hollow rotor and bearing substructure were considered and the heat-carrying ratios between two cooling air flows.From aforementioned works, it was exactly sure that the outer cooling flow through the gas film and foil gap played a dominant role on the heat removal of the heat generation inside hydrodynamic gas film.However, the pressure drop inside the outer cooling flow channel was dramatically increased with the increase of operational speed,such that the throughflow capacity was seriously affected in the high-speed operational situations.Therefore, the inner cooling inside a hollow shaft was essential for controlling the temperature rise of rotor.

    In the viewing of fundamental flow dynamics,the flow inside a hollow shaft is reasonably modeled as a classical internal flow inside a rotating pipe with an axial throughflow.With respect to the stationary situation, the flow structure inside a rotating pipe is significantly altered by the centrifugal force caused by the rotation, generally leading to a convective heat transfer enhancement but also a stronger pressure drop[12-15].For the rotating pipes, the installedrib turbulator is still an effective means for enhancing the convective heat transfer coefficient, following the same principle as that in the stationary pipes[16-18].To the best of our knowledge, the transverse rib or circular ring-type rib is the main configuration as applied to the stationary pipes, wherein multiple ring-type ribs are distributed along the axial direction to introduce the secondary flow and cause the flow separation and reattachment[19-21].Discrete ribs have the ability to interrupt the continuous development of the buoyancy layer near the wall,which can reduce the negative effect of buoyancy and enhance heat transfer to a certain extent[22-24].However, for the circular pipe that rotates about its axis, as the internal flow field is featured by a helicoidal flow structure caused by the superposing of rotation-induced zimuthal movement on the axial through flow, whether the ring-type rib plays the same role on heat transfer enhancement as that in the stationary situation is doubted.In particular,when the rotational speed is high enough the rotation-induced Zimuthal movement is dominate in the internal flow, the flow disturbance role of transverse ribs to the internal helical flow in the rotating situation is conjectured to be changed with respect to the stationary situation.Therefore the rib configuration has significant effects on the internal flow and heat transfer of a rotating pipe with an axial throughflow, which is still lack of deeper identification.

    The objective of the present investigation is to assess the effects of rib configuration on the flow and heat transfer inside a rotating rib-roughened pipe with an axial throughflow, by using a numerical methodology.Three rib configurations, including ring-type rib, discrete ring-type rib and longitudinal-type rib, are taken into consideration under a fixed axial-throughflow Reynolds number and a series of rotational Reynolds numbers.From current study, the influences of rotational speed or rotational Reynolds number on the pressure drop and convective heat transfer coefficient are illustrated.In particular, reasonable rib configurations are proposed in accordance with the rotational conditions,in the viewing of the comprehensive performance factor.

    1 Physical Model and Simplifications

    A typical bump-type hydrodynamic foil journal bearing consists of shaft, foil assembly and bearing sleeve, as schematically shown in Fig.1.

    Fig.1 Schematic diagram of hydrodynamic gas bearing structure

    During the operation, a self-acting offset of the rotor from the global bearing center is formed, with an eccentricity between the shaft and the bearing sleeve.The hydrodynamic pressure effect of the lubricating gas itself is built up in an eccentric circular gap around the circumference of the bearing, to provide the bearing capacity.Meanwhile, the heat generation inside the hydrodynamic gas film occurs due to the strong viscous heating effect.Considering the current study concerns the internal flow and heat transfer inside a rotating hollow shaft, a researched physical model is taken, as schematically displayed in Fig.2, wherein the supporting section in hydrodynamic foil journal bearing (L=25 mm) is simplified as the heating section by imposing a fixed heat load.In this physical model, the main geometric parameters are denoted as the follows: the radius of shaftR=10.95 mm, the inner diameter of rotating tubed=15 mm, the thickness of tubeb=3.45 mm.The front section relative to the heating section has a length ofLe=20 mm and rear section has a length ofLc=45 mm.In order to enhance the heat transfer capacity inside the rotating tube, the rib-roughened turbulators are installed.Three rib configurations are considered, referring as ring-type rib(RT)(Fig.2(a)), discrete ring-type rib(DRT) (Fig.2(b)) and longitudinal-type rib(LT) (Fig.2(c)).Seven rows of ring-type ribs or discrete ring-type ribs are distributed along the axial direction, and twelve longitudinal-type ribs are distributed in the circumference direction.It has been verified that the rib size is the best comprehensive effect in each rib type.All of the ribs have the same rib height ofe=1 mm and rib width ofδ=1 m.For the ring-type rib layout, the interval between adjacent ribs is designed asΔ=2.5 mm.For the discrete ring-type rib layout, the intersection angle corresponding to the circumference interval is designed asθ1=30°.For the longitudinal-type rib layout, the intersection angle between adjacent ribs is designed asθ2=24°.

    Fig.2 Schematic diagram of rib configurations inside rotating hollow rotor

    2 Description of Computational Methodology

    2.1 Numerical modeling

    Considering the internal flow in a pipe that rotates about its axis with a rotating angular velocityω, the governing equations for momentum and energy transports are modified in the relative velocity or rotating reference frame, wherein the relative velocity is used[25]

    wherevris the relative velocity andτthe viscous stress.The momentum equation contains two additional terms, the Coriolis acceleration term (2ω×vr) and the centripetal acceleration term (ω×ω×r).Since local temperature changes and local pressure changes will cause density changes, resulting in buoyancy, the centrifugal buoyancy termF=ρ0βΔT(ω×ω×r) is introduced, andρ0is the reference density.The energy equation uses the relative total enthalpyHrand (τ·vr) is the energy transport due to viscous dissipation.ρandkare the density and thermal conductivity of air, respectively.

    The corresponding boundary conditions are summarized as the follows:

    (1) For the rotating pipe, a rotating angular velocityωis given, and the rotational Reynolds numberReωis defined as

    whereμisthe dynamic viscosity of air.

    (2) The inlet of axial throughflow is set as velocity inlet by giving the mass flowrateQm.The inlet temperature is set as 300 K.The axial throughflow inlet Reynolds numberRexis defined as

    whereUinis the bulk-averaged velocity at tube inlet.

    (3) The outlet is set as a pressure out, by settingpout=101 325 Pa.

    (4) On the heating section, a constant heat fluxqis imposed.Whereas for the other solid surface, the adiabatic thermal condition is adopted.On the solid surface, no slip velocity is adopted.

    In the current study, the mass flowrate of axial throughflow is set asQm=5 kg/h and the corresponding axial-throughflow inlet Reynolds numberRexis about 6 400.The rotating angular velocityωis varied between 0 r/min and 1×105r/min, and the corresponding rotational Reynolds numberReωis varied from 0 to 3.77×104.The heat flow of the rotor surface under different conditions is obtained by numerical calculation, and the average value is selected as the thermal boundary condition of this study.The heat flux imposing on the heating section is set asq=6 950 W/m2.During the fluid-solid coupling, the thermal conductivity of solid wall is set as

    16.3 W/(m·K).

    2.2 Computational methodology

    The current numerical simulations are performed by using the commercial CFD solver ANSYS Fluent, wherein 3-D steady-state flow and temperature fields are solved.The working fluid is treated as an ideal gas on account of the compressible effect.Its thermal physical parameters are set as the ideal-gas based density.Kinetic-theory based specific heat and thermal conductivity.Sutherland law based viscosity.The second order upwind scheme is used for the spatial discretization of convective terms in the conservation equations.In the momentum equation, the PRESTO scheme is used for the discretization of pressure gradient term.With regard to the pressure-velocity coupling, the SIMPLE algorithm is adopted.The convergence is regarded to be achieved when all normalized residuals of concerned variables fell below 1.0×10-5.

    Computational grids are generated by using ICEM software.The structured meshes are adopted in the whole computational domain.The first layer of computational meshes is located at near-wall zone with a space of 0.02 mm normally to the solid surface, andy+is 1, which meets the requirement of the enhanced wall function.Then the computational meshes are stretched away from the viscous wall using a stretching ratio of 1.1.To make sure the computational result is independence of the computational mesh, a grid independency test is carried out in advance.As demonstrated in Fig.3 (l/Dis the ratio of the length between the position on the heating section and the inlet to the diameter of the channel), when the grid number is beyond 2.9 million the computed circumference-averaged temperature distribution at the heating section is not sensitive to the grid number.Thus a final grid number is selected around 3 million.

    Fig.3 Grid independence test

    To justify the current computational methodology on modelling the flow field inside a rotating pipe, an example validation is performed in advance against the experimental test by Imao et al[12].Beside, the selection of reasonable two-equation turbulence model for the the turbulence closure is also made from the comparison.As demonstrated in Fig.4 (r/Ris the ratio of the radial position to the radius of the pipe), it is confirmed that all the selected two-equation turbulence models provide nearly the same prediction on the azimuthal velocity profile in radial direction, andW/ωris the ratio of the tangential velocity to the tangential velocity of the rotating pipe.Regarding the axial velocity profile in radial direction (U/Uinis the ratio of the axial velocity to the inlet velocity), it is found both RNGk-εturbulence model and Realizablek-εturbulence model could provide more correct prediction with respect to the experimental data.Referring previous researches relating to the numerical simulations of rotating flow[26-28], RNGk-εturbulence model is selected in current study.

    Fig.4 Validation of turbulence models

    3 Results and Discussion

    3.1 Flow fields

    In a rotating tube, the Coriolis force and centrifugal force caused by the rotation make the internal flow complicated when compared to the stationary situation.Fig.5 shows the local three-dimensional flow trajectories inside a smooth pipe(SP), colored by the temperature value.Generally, the flow field inside a rotating pipe suffers to dual roles of forced axial through flow and rotation-induced azimuthal flow.When the pipe rotates at a small rotating speed or rotational Reynolds number (such asReω=807), the azimuthal flow induced by pipe rotation makes the flow trajectory deflect slightly in the circumferential direction.At a big radial location where the local flow zone is close to the rotating wall, the flow trajectory deflection behaves more strongly because of radially growing centrifugal forces.In a total, the forced axial through flow is dominant in the situations of small rotating speeds.However, when the pipe rotates at a high rotating speed or rotational Reynolds number (such asReω=3.77×104), the rotation-induced azimuthal flow is dominate on the internal flow field, making the flow trajectories take on a significant spiral flow feature,because of the superposing of axial through flow.When compared to the small rotational Reynolds number situation, the axial flow in the entrance section is concentrated to the central zone of the rotating pipe, indicating that the inlet flow of forced axial flow is locally blocked at the big radial location because of the rotation effect.Along the axial direction, the spiral flow trajectory extends outward in the radial direction gradually.Meanwhile, the fluid temperature in the near-wall zone of heating section is obviously higher than that in the small rotational Reynolds number situation, indicating that the convective heat transfer capacity is raised in the high rotational Reynolds number situation.

    Fig.5 Flow trajectories inside smooth pipe under stationary and rotating states

    The rotation-induced azimuthal flow changes the static pressure distribution inside the pipe, as demonstrated in Fig.6.In the stationary situation,the internal flow along the flow direction is classically characterized by a relatively regular favourablepressure flow.With the increase of rotational Reynolds number, the static pressure in the near-wall zone is increased because the near-wall fluid suffers the stronger Coriolis force and centrifugal force.Meanwhile, the static pressure in the central zone of a rotating pipe is decreased on the otherwise.Especially, when the pipe rotates at a high rotational Reynolds number, a low-pressure zone occurs inside the rotating pipe, occupying nearly entire axial length.This formation of low-pressure zone not only affects the internal flow inside the pipe, but also affects the outflow at the outlet.Due to this lowpressure zone, the external ambient air would even be suctioned into the central zone of a rotating pipe from the outlet plane.

    Fig.6 Pressure contours in side smooth pipe

    As the internal flow field is dramatically changed by the rotational effect, the roles of ribroughened turbulators on the flow disturbance are conjectured to be changed in accordance with the rotational Reynolds number.The local streamlines in the vicinity of rib-roughened zone are displayed in Figs.7—9, corresponding to different rib configurations.In these figures, the streamlines are colored by the dimensionless velocity, or the velocity ratio of locally absolute velocity to the bulk-averaged velocity at tube inlet (Uin).

    Fig.7 Local streamlines inside RT

    Fig.8 Local streamlines inside DRT

    Fig.9 Local streamlines inside LT

    For the RT, a large-scale vortex structure is formed in the interval region after the axial through flow passes through the ribs in the stationary situation, as displayed in Fig.7(a).The transverse rib plays a promising flow disturbance role on the internal flow.As the rotational Reynolds number increases, the flow velocity in the high radial location is increased.The rotation-induced azimuthal flow becomes dominate gradually, making the size of vortex structure between ribs reduction radially, as displayed in Fig.7(b).In particular, under a high rotational Reynolds number, the rotation-induced azimuthal flow is dominant significantly, the absolute velocity near the wall surface is obviously greater than the bulk-averaged velocity of the axial throughflow, the flow disturbance role of ring-type rib on the helicoidal flow is apparently weakened with respective to the stationary situation, as demonstrated in Fig.7(c).Also, it is conjectured that the heat transfer enhancement role of ring-type rib layout would be reduced in the high rotational Reynolds number situation.

    For the DRT, the flow disturbance role on the internal flow is nearly the same as the ring-type rib layout.Because there is a gap between discrete ring ribs in the circumferential direction, the local flow disturbance also appears in the vicinity of this gap in the stationary situation, as displayed in Fig.8(a).In the rotating situations, the discrete ring ribs could introduce a flow disturbance role to the rotation-induced azimuthal flow, as displayed in Fig.8(b) and Fig.8(c).

    For the LT, the flow disturbance role of longitudinal rib on the axial flow is very weak in the stationary situation, as displayed in Fig.9(a).However, in the rotating situations, the longitudinal rib layout shows its superiority on disturbing the local flow in the rib-roughened region, to the rotation-induced azimuthal flow.Under a small rotational Reynolds number, the rotation-induced azimuthal velocity is small, such that the vortex structure is composed by the low-velocity fluid, as displayed in Fig.9(b).Under a high rotational Reynolds number, the rotationinduced azimuthal velocity is strong, such that the vortex structure is composed by the high-velocity fluid, as displayed in Fig.9(c).

    From the above-mentioned analysis, it is confirmed that the rotation has a significant effect on the internal flow field, altering the flow disturbance roles of rib layout on the internal flow.As a consequence, the heat transfer enhancement role of rib layout would be changed in accordance with the rotational Reynolds number.

    3.2 Flow and heat transfer performances

    Fig.10 displays the temperature contours in the solid region at the middle plane of heating section,in the stationary situation.It is seen that temperature of heating section could be effectively reduced by utilizing the inner rib-roughened configurations,when compared to the SP at the same axial throughflow mass flow rate.By comparing three different rib configurations, the ring-type ribs, either in continuous (RT) or discrete (DRT), provide better cooling roles than the LT.With respect to the SP,the RT could reduce the peak temperature of heating section about 64K.Only about 34 K reduction is achieved by the LT.

    Fig.10 Temperature contours on middle-axial plane of heating section under Reω=0

    Fig.11 displays the temperature contours in the solid region at the middle plane of heating section,in the rotating situation atReω=3.77×104.When compared to the stationary situation, the peak temperature is reduced obviously because the strong helicoidal flow enhances the convective heat transfer capacity.Taking the SP as an example, the peak temperature is reduced about 70 K in the high rotational Reynolds number situation relative to the stationary situation.With respect to the SP, the rotation effect on internal convective heat transfer of ribroughened tube is weakened.In particular, for the RT, the peak temperature underReω=3.77×104has only about 16 K reduction in relative to the stationary situation.When compared to the SP at the same rotational Reynolds number, only about 9 K reduction in the peak temperature is achieved by the RT.Apparently, in the high rotational Reynolds number situation, the LT presents the best convective heat transfer enhancement owing to its strongest flow disturbance role to the rotation-induced azimuthal flow.

    Fig.11 Temperature contours on middle-axial plane of heating section under Reω=3.77×104

    Fig.12 displays the peak temperature variations of heating section with the rotational Reynolds number.It is seen that the peak temperature is decreased with the increase of rotational Reynolds number generally.For the SP and LT, the peak temperature is reduced rapidly with the rotational Reynolds number whenReωis increased from 0 to 1×104.While for the RT and DRT, the peak temperature is not sensitive to the rotational Reynolds number in a wide range.It is noteworthy that when the rotational Reynolds number is less than 2 400, the RT and DRT provide a better heat transfer enhancement than the LT.Whereas under high rotational Reynolds numbers, the LT is more pronounced for reducing the peak temperature of heating section.

    Fig.12 Effects of rotational Reynolds number and rib type on peak temperature of rotor

    The rib-roughened configuration inside a rotating tube enhances the convective heat transfer but also leads to a serious pressure loss.In order to assess the comprehensive effects of ribbed configuration, a comprehensive performance factor is introduced,which is defined as

    whereNuis the area-averaged Nusselt number andfthe friction factor.The subscript “0” represents the baseline case for the comparison.Here the corresponding smooth tube under the same rotational Reynolds number is selected as baseline case.

    The area-averaged Nusselt number over the heating section and the friction factor across the entire pipe are defined as

    wherehis the area-averaged convective heat transfer coefficient over the inner surface of heating section;λandρare the thermal conductivity and the density of the air, respectively; (pin-pout) is the pressure drop of the internal flow across the lengthLof the pipe.

    For the SP, the effects of rotation on the areaaveraged Nusselt numberNu0and friction factorf0are displayed in Fig.13.In general, both the convective heat transfer and the flow loss are increased monotonically with the rotational Reynolds number.However, the influencing tendency of rotational Reynolds number onNu0is somewhat different from that onf0.Regarding the convective heat transfer capacity, as seen in Fig.13(a),Nu0is increased rapidly along the rotational Reynolds number in a nearly linear mode whenReωis increased from 0 to 1×104.WhenReωis beyond 1×104, the increase rate ofNu0alongReωbecomes degraded.Whereas regarding the friction factor, as seen in Fig.13(b),f0is weakly increased along the rotational Reynolds number whenReωis less than 0 to 1×104.WhenReωis beyond 2×104, the friction factor is increased sharply with the rotational Reynolds number.Identified from Fig.13, under a high rotational Reynolds number ofReω=3.77×104, the Nusselt number and friction factor are about two times and five times of the corresponding values in the stationary situation, respectively.

    Fig.13 Effects of rotational Reynolds number on Nu0 and f0 for smooth pipe

    Fig.14 shows the variations of average Nusselt numberNuand friction factorfwith the rotational Reynolds number, for the rib-roughened pipes.It is confirmed that the effects of rotation on the Nusselt number and friction factor of rib-roughened pipe are similar to the smooth pipe.However, the variation tendencies of convective heat transfer and friction with rotational Reynolds number are somewhat distinct for different rib-roughened configurations, because of different flow disturbance roles of rib layout on the internal flow.Taking the SP as a baseline case for comparison, Fig.15 presents directly the variation patterns of convective heat transfer enhancement ratio (——-Nu/Nu0), friction factor ratio(f/f0) and comprehensive performance factor(PEC) along with the the rotational Reynolds number, for three different rib-roughened configurations.With respect to the baseline SP, all the ribroughened pipes show their heat transfer enhancement role but also increase the flow loss.Their influencing roles are tightly dependent on the rotational Reynolds number, showing a non-monotonic variation.

    Fig.14 Effects of rotational Reynolds number on Nu and f for rib-roughened pipes

    Fig.15 Comparison between different rib configurations in view of comprehensive factor

    Under small rotational Reynolds numbers(such asReω<2 400), the RT and DRT provide bigger convective heat transfer enhancement ratio but also bigger friction factor ratio than the LT.In the viewing of the comprehensive performance factor that considers both heat transfer enhancement and flow loss together, the DRT is identified to be the best among the current three rib configurations.

    Under high rotational Reynolds numbers(such asReω>2×104), the LT provide bigger convective heat transfer enhancement ratio but also bigger friction factor ratio than the RT.It is noteworthy the RT show no advantage over SP, and even lead to a reduction in the comprehensive performance factor when compared to the smooth pipe.Therefore, for the high-speed rotating tube, the LT is suggested to be a more promising rib configuration as applied to the heat transfer enhancement.

    4 Conclusions

    A numerical study is conducted to illustrate the flow and heat transfer characteristics inside a rotating rib-roughened pipe with an axial throughflow,under a fixed axial-throughflow Reynolds numberRexof 6 400 and a series of rotational Reynolds numberReωranging from 0 to 3.77×104.Three rib configurations, including RT, DRT and LT, are taken into consideration.From current study, the main conclusions are summarized as the follows:

    (1) The flow field inside a rotating pipe is of high complexity, due to the dual influencing aspects from the axial through flow and the rotation-induced azimuthal flow.The dominant flow feature and subsequently the flow of rib-roughened flow disturbance role is tightly dependent on the competition between axial through flow and rotation-induced azimuthal flow.

    (2) Under small rotational Reynolds numbers,the axial through flow is dominant such that the RT plays a significant flow disturbance role and subsequently the strong heat transfer enhancement.Whereas under high rotational Reynolds numbers,the rotation-induced azimuthal flow is dominant such that the LT is superior to the RT on the flow disturbance and subsequently heat transfer enhancement.

    (3) Evaluated on the comprehensive performance factor that considers both heat transfer enhancement and flow loss together, the DRT is identified to be the best among the current three rib configurations under small rotational Reynolds numbers.However, under high rotational Reynolds numbers, the LT layout is a promising rib configuration as applied to the heat transfer enhancement.

    国产精品亚洲av一区麻豆| 母亲3免费完整高清在线观看| 亚洲色图 男人天堂 中文字幕| 欧美另类亚洲清纯唯美| 亚洲成av片中文字幕在线观看| 国产成人欧美在线观看 | 亚洲成国产人片在线观看| 高清视频免费观看一区二区| videos熟女内射| 999久久久国产精品视频| 777久久人妻少妇嫩草av网站| 变态另类成人亚洲欧美熟女 | 久久久国产成人免费| 国产精品自产拍在线观看55亚洲 | 狂野欧美激情性xxxx| 精品一品国产午夜福利视频| 国产精品美女特级片免费视频播放器 | 久久国产精品影院| 校园春色视频在线观看| 亚洲五月婷婷丁香| 老鸭窝网址在线观看| 中亚洲国语对白在线视频| 18禁黄网站禁片午夜丰满| 久久久久久亚洲精品国产蜜桃av| 欧美乱妇无乱码| 久久人妻福利社区极品人妻图片| 丁香欧美五月| 亚洲午夜精品一区,二区,三区| 亚洲片人在线观看| 一本大道久久a久久精品| 热99国产精品久久久久久7| 亚洲专区字幕在线| 成年人午夜在线观看视频| 91精品国产国语对白视频| 亚洲精品美女久久久久99蜜臀| 精品视频人人做人人爽| 黄色 视频免费看| 欧美激情极品国产一区二区三区| 一级,二级,三级黄色视频| 日韩制服丝袜自拍偷拍| 熟女少妇亚洲综合色aaa.| 亚洲人成电影观看| 国产高清激情床上av| 久久天堂一区二区三区四区| 两个人免费观看高清视频| 亚洲专区中文字幕在线| 女人高潮潮喷娇喘18禁视频| 国产精品.久久久| 香蕉国产在线看| 人妻久久中文字幕网| 美女福利国产在线| 成人亚洲精品一区在线观看| 宅男免费午夜| 欧美激情极品国产一区二区三区| 成在线人永久免费视频| 国产一区二区激情短视频| 亚洲精品自拍成人| 国产精品偷伦视频观看了| 女人精品久久久久毛片| 久久国产精品人妻蜜桃| 男人操女人黄网站| 欧美日韩成人在线一区二区| 丝袜人妻中文字幕| 女人精品久久久久毛片| 亚洲黑人精品在线| 久久人妻福利社区极品人妻图片| 十八禁人妻一区二区| 日韩有码中文字幕| 欧美日韩福利视频一区二区| 免费不卡黄色视频| 校园春色视频在线观看| 一进一出抽搐动态| tocl精华| av福利片在线| 久久久精品国产亚洲av高清涩受| 99国产精品一区二区蜜桃av | 女性生殖器流出的白浆| 国产高清激情床上av| 99re在线观看精品视频| 黄色女人牲交| 一级a爱片免费观看的视频| 中文亚洲av片在线观看爽 | 身体一侧抽搐| 999久久久国产精品视频| 国产在线精品亚洲第一网站| 久久午夜综合久久蜜桃| 久久热在线av| 精品国内亚洲2022精品成人 | 亚洲成a人片在线一区二区| 在线观看午夜福利视频| 国产精品免费一区二区三区在线 | 一区福利在线观看| 亚洲片人在线观看| a级片在线免费高清观看视频| 午夜福利乱码中文字幕| 两个人免费观看高清视频| 免费在线观看完整版高清| 国产蜜桃级精品一区二区三区 | 国产免费av片在线观看野外av| 国产不卡av网站在线观看| 1024香蕉在线观看| 精品国产美女av久久久久小说| 成年动漫av网址| 婷婷成人精品国产| 日韩欧美国产一区二区入口| av线在线观看网站| 丝瓜视频免费看黄片| 精品少妇一区二区三区视频日本电影| 色老头精品视频在线观看| 欧美不卡视频在线免费观看 | 免费黄频网站在线观看国产| 国产亚洲av高清不卡| 欧美乱码精品一区二区三区| 一本综合久久免费| 大片电影免费在线观看免费| svipshipincom国产片| 搡老乐熟女国产| 欧美黄色淫秽网站| 亚洲黑人精品在线| а√天堂www在线а√下载 | 欧美日韩亚洲高清精品| 国产在线一区二区三区精| 成人国产一区最新在线观看| 国产亚洲欧美精品永久| 大型av网站在线播放| 叶爱在线成人免费视频播放| ponron亚洲| av网站免费在线观看视频| 亚洲九九香蕉| 国产精品成人在线| 成人18禁高潮啪啪吃奶动态图| 丰满饥渴人妻一区二区三| www.999成人在线观看| 亚洲精品自拍成人| 日本wwww免费看| bbb黄色大片| 一边摸一边抽搐一进一小说 | 色在线成人网| 午夜福利乱码中文字幕| 久久中文字幕人妻熟女| 99国产精品免费福利视频| 国产亚洲精品久久久久5区| 香蕉国产在线看| 丝袜美足系列| 午夜福利乱码中文字幕| 美女 人体艺术 gogo| 丁香六月欧美| 又黄又粗又硬又大视频| 男男h啪啪无遮挡| 免费在线观看黄色视频的| 麻豆av在线久日| 国产麻豆69| 国产亚洲精品久久久久5区| 亚洲五月婷婷丁香| 久久国产精品人妻蜜桃| 精品久久久久久,| 母亲3免费完整高清在线观看| 免费久久久久久久精品成人欧美视频| 成人国语在线视频| 国产精品九九99| 纯流量卡能插随身wifi吗| 黄网站色视频无遮挡免费观看| 国产日韩欧美亚洲二区| 亚洲熟妇中文字幕五十中出 | 纯流量卡能插随身wifi吗| 国产精品.久久久| 欧美一级毛片孕妇| 丰满饥渴人妻一区二区三| 91麻豆精品激情在线观看国产 | 欧美激情极品国产一区二区三区| 乱人伦中国视频| 欧美+亚洲+日韩+国产| 伊人久久大香线蕉亚洲五| 男女床上黄色一级片免费看| 黄频高清免费视频| 一级,二级,三级黄色视频| 操美女的视频在线观看| 黑人巨大精品欧美一区二区mp4| 岛国毛片在线播放| 成熟少妇高潮喷水视频| 午夜精品国产一区二区电影| 国产一区二区三区综合在线观看| 男人操女人黄网站| 人人澡人人妻人| 亚洲精品在线美女| 婷婷成人精品国产| 欧美成人免费av一区二区三区 | 1024视频免费在线观看| 成人手机av| 两人在一起打扑克的视频| 国产成人免费观看mmmm| 精品国产亚洲在线| 天天操日日干夜夜撸| 久久ye,这里只有精品| 老司机在亚洲福利影院| 久久精品aⅴ一区二区三区四区| 亚洲av电影在线进入| 91av网站免费观看| 最近最新中文字幕大全免费视频| 91国产中文字幕| 日韩一卡2卡3卡4卡2021年| 日本一区二区免费在线视频| 99热国产这里只有精品6| 黑人猛操日本美女一级片| 咕卡用的链子| 三上悠亚av全集在线观看| 日本黄色日本黄色录像| 久久香蕉国产精品| 国产黄色免费在线视频| 国产激情久久老熟女| √禁漫天堂资源中文www| 国产无遮挡羞羞视频在线观看| 成人亚洲精品一区在线观看| 久久久久久久久免费视频了| 一区二区日韩欧美中文字幕| av片东京热男人的天堂| 下体分泌物呈黄色| 少妇裸体淫交视频免费看高清 | 色综合欧美亚洲国产小说| 老司机靠b影院| 亚洲av第一区精品v没综合| 高潮久久久久久久久久久不卡| 欧美精品高潮呻吟av久久| 制服诱惑二区| 国产亚洲精品一区二区www | 免费少妇av软件| 波多野结衣一区麻豆| 成人精品一区二区免费| 男女床上黄色一级片免费看| 国产精品 欧美亚洲| 日韩免费高清中文字幕av| 一本一本久久a久久精品综合妖精| 精品一区二区三区四区五区乱码| 国产成+人综合+亚洲专区| 国产一区在线观看成人免费| 视频区欧美日本亚洲| 人成视频在线观看免费观看| 久久久国产成人免费| 亚洲免费av在线视频| 色尼玛亚洲综合影院| 人妻久久中文字幕网| 亚洲av成人一区二区三| 岛国在线观看网站| 精品国内亚洲2022精品成人 | 天堂√8在线中文| 中文字幕色久视频| av免费在线观看网站| www.999成人在线观看| 亚洲国产欧美一区二区综合| 久久久久国内视频| 亚洲欧美激情综合另类| 久久人人爽av亚洲精品天堂| 欧美人与性动交α欧美软件| 亚洲成人免费电影在线观看| 久久狼人影院| a级毛片在线看网站| 国产精品一区二区免费欧美| 波多野结衣av一区二区av| 欧美色视频一区免费| 欧美精品啪啪一区二区三区| 香蕉丝袜av| 窝窝影院91人妻| 极品少妇高潮喷水抽搐| 午夜福利在线免费观看网站| 99久久国产精品久久久| 亚洲成a人片在线一区二区| 在线观看午夜福利视频| 99精国产麻豆久久婷婷| 精品人妻熟女毛片av久久网站| 日韩人妻精品一区2区三区| 久久久久精品国产欧美久久久| 国产野战对白在线观看| 极品人妻少妇av视频| 国产99白浆流出| 人成视频在线观看免费观看| 国产aⅴ精品一区二区三区波| 日本黄色视频三级网站网址 | 深夜精品福利| 91av网站免费观看| 免费看十八禁软件| 久久草成人影院| 不卡av一区二区三区| 18禁观看日本| 天天添夜夜摸| 狂野欧美激情性xxxx| 日韩中文字幕欧美一区二区| 久久精品成人免费网站| 亚洲成人手机| 啦啦啦免费观看视频1| a在线观看视频网站| 王馨瑶露胸无遮挡在线观看| 91字幕亚洲| 999精品在线视频| 午夜免费成人在线视频| 精品人妻1区二区| 欧美老熟妇乱子伦牲交| av超薄肉色丝袜交足视频| 久久影院123| 黄色成人免费大全| 国产午夜精品久久久久久| 人成视频在线观看免费观看| 亚洲va日本ⅴa欧美va伊人久久| 捣出白浆h1v1| 人人妻人人澡人人看| 搡老熟女国产l中国老女人| 男人舔女人的私密视频| 成年动漫av网址| 18禁国产床啪视频网站| 亚洲精品美女久久久久99蜜臀| 一级片免费观看大全| 亚洲第一欧美日韩一区二区三区| 国产免费男女视频| 国产精品影院久久| 中文欧美无线码| 成年版毛片免费区| 最近最新中文字幕大全免费视频| 国产精品乱码一区二三区的特点 | 日本一区二区免费在线视频| www.熟女人妻精品国产| av不卡在线播放| 欧美日韩黄片免| 天堂√8在线中文| 国产成人免费观看mmmm| 王馨瑶露胸无遮挡在线观看| 日韩成人在线观看一区二区三区| 欧美精品人与动牲交sv欧美| 丰满迷人的少妇在线观看| 美女 人体艺术 gogo| 久久久国产成人精品二区 | 法律面前人人平等表现在哪些方面| 又黄又爽又免费观看的视频| 欧美 日韩 精品 国产| 精品卡一卡二卡四卡免费| 嫩草影视91久久| 亚洲少妇的诱惑av| 少妇裸体淫交视频免费看高清 | 18禁黄网站禁片午夜丰满| 国产人伦9x9x在线观看| 欧美日韩精品网址| 亚洲熟妇中文字幕五十中出 | av超薄肉色丝袜交足视频| 正在播放国产对白刺激| 亚洲国产精品sss在线观看 | 日韩大码丰满熟妇| 看免费av毛片| 成年人免费黄色播放视频| 免费av中文字幕在线| 一a级毛片在线观看| tocl精华| 欧美黑人欧美精品刺激| 免费看十八禁软件| 国产精品 欧美亚洲| 一二三四在线观看免费中文在| 国产精品国产av在线观看| 色精品久久人妻99蜜桃| 制服人妻中文乱码| 妹子高潮喷水视频| 亚洲av日韩精品久久久久久密| 国产亚洲欧美精品永久| 国产精品香港三级国产av潘金莲| 欧洲精品卡2卡3卡4卡5卡区| 在线永久观看黄色视频| 大码成人一级视频| 亚洲专区中文字幕在线| 中亚洲国语对白在线视频| 精品少妇久久久久久888优播| 美女视频免费永久观看网站| 精品一品国产午夜福利视频| 亚洲精品美女久久av网站| 1024视频免费在线观看| 亚洲 欧美一区二区三区| 久久久国产欧美日韩av| 中文字幕另类日韩欧美亚洲嫩草| 欧美久久黑人一区二区| 每晚都被弄得嗷嗷叫到高潮| 国产精品成人在线| 欧美日韩亚洲高清精品| 精品国产一区二区久久| 成人av一区二区三区在线看| 狠狠狠狠99中文字幕| 久久这里只有精品19| 黑人操中国人逼视频| 老司机在亚洲福利影院| 午夜免费成人在线视频| 国产精品自产拍在线观看55亚洲 | 99国产精品99久久久久| 亚洲成人免费电影在线观看| 亚洲欧美一区二区三区黑人| 丁香六月欧美| 大陆偷拍与自拍| 免费一级毛片在线播放高清视频 | 日韩精品免费视频一区二区三区| 视频区图区小说| 国产精品乱码一区二三区的特点 | 欧美日韩亚洲综合一区二区三区_| 麻豆av在线久日| 热re99久久国产66热| 国产男靠女视频免费网站| 狠狠狠狠99中文字幕| 久久ye,这里只有精品| 国产在线观看jvid| 久久精品国产99精品国产亚洲性色 | 国产乱人伦免费视频| 一区福利在线观看| 日韩人妻精品一区2区三区| 中文字幕人妻熟女乱码| 夫妻午夜视频| 亚洲欧美激情在线| 精品人妻熟女毛片av久久网站| av片东京热男人的天堂| 国产亚洲精品久久久久久毛片 | 妹子高潮喷水视频| 国产欧美日韩综合在线一区二区| 精品亚洲成a人片在线观看| 99精品久久久久人妻精品| 天天躁夜夜躁狠狠躁躁| 后天国语完整版免费观看| av免费在线观看网站| 老汉色av国产亚洲站长工具| 久久香蕉激情| 热re99久久国产66热| 久久久久久免费高清国产稀缺| 老汉色av国产亚洲站长工具| x7x7x7水蜜桃| 黄色成人免费大全| 国产欧美日韩精品亚洲av| www.精华液| 久久人人97超碰香蕉20202| 国产精华一区二区三区| 亚洲国产欧美日韩在线播放| 性少妇av在线| 精品一品国产午夜福利视频| www.精华液| 亚洲国产看品久久| 久久久国产成人精品二区 | 两个人看的免费小视频| 丝袜在线中文字幕| 亚洲av电影在线进入| 很黄的视频免费| 亚洲五月色婷婷综合| 视频在线观看一区二区三区| 成人国语在线视频| 国产不卡一卡二| 亚洲欧洲精品一区二区精品久久久| 国产一卡二卡三卡精品| 高潮久久久久久久久久久不卡| 国产亚洲欧美在线一区二区| 超碰成人久久| 这个男人来自地球电影免费观看| 不卡一级毛片| 色婷婷av一区二区三区视频| avwww免费| av一本久久久久| 国产不卡av网站在线观看| 怎么达到女性高潮| 欧美日韩亚洲综合一区二区三区_| 亚洲国产精品合色在线| 99riav亚洲国产免费| 国产亚洲av高清不卡| 亚洲av日韩在线播放| 亚洲欧美一区二区三区久久| 国产精品国产av在线观看| 90打野战视频偷拍视频| 狂野欧美激情性xxxx| 欧美日韩黄片免| 亚洲七黄色美女视频| 真人做人爱边吃奶动态| 欧美黑人欧美精品刺激| av天堂久久9| 国产av一区二区精品久久| 国产男女内射视频| 1024香蕉在线观看| 99国产精品99久久久久| 久久99一区二区三区| 水蜜桃什么品种好| 男人的好看免费观看在线视频 | 黑人欧美特级aaaaaa片| 三上悠亚av全集在线观看| 亚洲精品粉嫩美女一区| 大片电影免费在线观看免费| 身体一侧抽搐| 天天操日日干夜夜撸| 一本综合久久免费| 久久国产精品大桥未久av| 久久精品国产99精品国产亚洲性色 | 成年人免费黄色播放视频| 韩国av一区二区三区四区| av网站免费在线观看视频| 国产亚洲精品久久久久久毛片 | 国产精品香港三级国产av潘金莲| 我的亚洲天堂| 捣出白浆h1v1| 亚洲 欧美一区二区三区| 色在线成人网| 久久香蕉国产精品| 国产精品永久免费网站| tube8黄色片| 一区二区三区精品91| 国产精品乱码一区二三区的特点 | 女警被强在线播放| 脱女人内裤的视频| 中文字幕精品免费在线观看视频| 高清视频免费观看一区二区| 成人国语在线视频| 久久久精品区二区三区| 男女之事视频高清在线观看| 一级毛片女人18水好多| 国产高清激情床上av| 国产成人欧美在线观看 | 中文字幕制服av| 欧美乱妇无乱码| 午夜91福利影院| 国产精品二区激情视频| 国产区一区二久久| 欧美成人午夜精品| 这个男人来自地球电影免费观看| 不卡一级毛片| 狂野欧美激情性xxxx| 免费在线观看完整版高清| 91九色精品人成在线观看| 性色av乱码一区二区三区2| 国产激情欧美一区二区| 国产免费现黄频在线看| 老司机亚洲免费影院| 老鸭窝网址在线观看| 亚洲五月色婷婷综合| 如日韩欧美国产精品一区二区三区| 免费日韩欧美在线观看| 久久天堂一区二区三区四区| 美女 人体艺术 gogo| 国产精品亚洲一级av第二区| 精品乱码久久久久久99久播| 99精品久久久久人妻精品| 亚洲欧美激情在线| 国内毛片毛片毛片毛片毛片| 午夜免费鲁丝| 国产精品久久久人人做人人爽| 国产国语露脸激情在线看| 欧美 亚洲 国产 日韩一| 国产av精品麻豆| 看免费av毛片| 国产亚洲精品久久久久5区| 黄色成人免费大全| 极品人妻少妇av视频| 一进一出好大好爽视频| 亚洲av成人一区二区三| 国产精品久久视频播放| 老司机在亚洲福利影院| a级毛片在线看网站| 在线十欧美十亚洲十日本专区| 91av网站免费观看| 国产成人免费观看mmmm| 国产无遮挡羞羞视频在线观看| 18禁美女被吸乳视频| 欧美久久黑人一区二区| 午夜福利乱码中文字幕| 日本vs欧美在线观看视频| 欧美亚洲日本最大视频资源| 日韩大码丰满熟妇| 少妇 在线观看| 久久久久久久精品吃奶| 日韩免费av在线播放| 国产精品欧美亚洲77777| 国产欧美亚洲国产| 欧美精品亚洲一区二区| 亚洲欧美激情综合另类| 精品乱码久久久久久99久播| 看免费av毛片| 国产精品.久久久| 欧美乱码精品一区二区三区| 一个人免费在线观看的高清视频| 久久国产亚洲av麻豆专区| 久热这里只有精品99| 热re99久久精品国产66热6| 超色免费av| 身体一侧抽搐| 亚洲一卡2卡3卡4卡5卡精品中文| 国产无遮挡羞羞视频在线观看| 别揉我奶头~嗯~啊~动态视频| 日韩欧美国产一区二区入口| 成人影院久久| 久久人妻熟女aⅴ| 欧美日韩亚洲高清精品| 多毛熟女@视频| 露出奶头的视频| 日韩中文字幕欧美一区二区| 久久久国产成人免费| av天堂在线播放| 国产亚洲欧美精品永久| 一进一出好大好爽视频| 99精品在免费线老司机午夜| 亚洲,欧美精品.| 国产片内射在线| 午夜福利,免费看| 最新美女视频免费是黄的| 久久中文看片网| 精品欧美一区二区三区在线| 精品国产一区二区三区四区第35| 91麻豆av在线| 欧美黑人精品巨大| 人人妻人人添人人爽欧美一区卜| 久久久久精品人妻al黑| 少妇猛男粗大的猛烈进出视频| 精品国产一区二区三区久久久樱花| 欧美精品啪啪一区二区三区| 日韩大码丰满熟妇| 欧美+亚洲+日韩+国产| 可以免费在线观看a视频的电影网站| 亚洲国产看品久久| 伦理电影免费视频| 人妻一区二区av| 黄色毛片三级朝国网站| 99国产综合亚洲精品| 欧美日韩av久久| 久久这里只有精品19|