• <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免费在线看不卡| 欧美xxⅹ黑人| 人人妻人人澡人人爽人人夜夜 | 亚洲图色成人| 久99久视频精品免费| 亚洲精品久久午夜乱码| 亚洲av电影在线观看一区二区三区 | 国产成人a∨麻豆精品| 色综合亚洲欧美另类图片| 国内少妇人妻偷人精品xxx网站| 成人亚洲欧美一区二区av| 不卡视频在线观看欧美| 中国美白少妇内射xxxbb| 日韩一区二区视频免费看| 国产在视频线在精品| 国产亚洲一区二区精品| 日韩强制内射视频| 在线免费观看的www视频| 菩萨蛮人人尽说江南好唐韦庄| 七月丁香在线播放| 人人妻人人澡欧美一区二区| 亚洲激情五月婷婷啪啪| 夜夜看夜夜爽夜夜摸| 一级毛片黄色毛片免费观看视频| 亚洲第一区二区三区不卡| 岛国毛片在线播放| 亚洲,欧美,日韩| 国产成人午夜福利电影在线观看| 最近中文字幕高清免费大全6| 国产精品爽爽va在线观看网站| 久久久久久伊人网av| 亚洲精品,欧美精品| 亚洲精品,欧美精品| 婷婷色综合大香蕉| 边亲边吃奶的免费视频| 天堂中文最新版在线下载 | 国产淫片久久久久久久久| 国产高清国产精品国产三级 | 97精品久久久久久久久久精品| 毛片女人毛片| 好男人在线观看高清免费视频| 偷拍熟女少妇极品色| 色综合站精品国产| 欧美精品国产亚洲| 有码 亚洲区| 欧美精品国产亚洲| 国产成年人精品一区二区| 国产 亚洲一区二区三区 | 免费看a级黄色片| 亚洲欧美精品自产自拍| 国产成人freesex在线| 岛国毛片在线播放| 精品不卡国产一区二区三区| 国产白丝娇喘喷水9色精品| 亚洲成色77777| 国产精品一区二区三区四区免费观看| 一个人看视频在线观看www免费| 中文字幕人妻熟人妻熟丝袜美| 亚洲国产成人一精品久久久| 国产精品一二三区在线看| 国产老妇伦熟女老妇高清| 亚洲av一区综合| 国产精品爽爽va在线观看网站| 天美传媒精品一区二区| 搡老妇女老女人老熟妇| 亚洲国产精品成人综合色| 国产在视频线精品| 能在线免费看毛片的网站| 亚洲av电影在线观看一区二区三区 | 熟女电影av网| 99视频精品全部免费 在线| 老司机影院毛片| 久久久久久久久久人人人人人人| 少妇熟女aⅴ在线视频| 国产毛片a区久久久久| av在线亚洲专区| 人妻制服诱惑在线中文字幕| 你懂的网址亚洲精品在线观看| 我的老师免费观看完整版| 男的添女的下面高潮视频| 亚洲乱码一区二区免费版| 爱豆传媒免费全集在线观看| 一级毛片 在线播放| 亚洲综合色惰| 国产亚洲91精品色在线| 视频中文字幕在线观看| 免费观看无遮挡的男女| 国产亚洲精品久久久com| 在线免费观看不下载黄p国产| 久久精品久久久久久久性| 亚洲国产精品国产精品| 观看免费一级毛片| 丝袜喷水一区| 在线观看人妻少妇| 国产在视频线在精品| 最近最新中文字幕大全电影3| 国产淫语在线视频| 午夜免费男女啪啪视频观看| 日日摸夜夜添夜夜添av毛片| 国产成人freesex在线| 精品熟女少妇av免费看| 国产午夜精品久久久久久一区二区三区| 国内揄拍国产精品人妻在线| 97超视频在线观看视频| 国产精品三级大全| 成人漫画全彩无遮挡| 欧美xxxx黑人xx丫x性爽| 亚洲高清免费不卡视频| 99久国产av精品国产电影| 成年人午夜在线观看视频 | 国产成人a区在线观看| 美女被艹到高潮喷水动态| 欧美高清成人免费视频www| 看免费成人av毛片| 高清午夜精品一区二区三区| 中文天堂在线官网| 男女边摸边吃奶| 深夜a级毛片| 草草在线视频免费看| 少妇的逼好多水| 国产精品女同一区二区软件| 美女内射精品一级片tv| 18禁在线无遮挡免费观看视频| 亚洲天堂国产精品一区在线| 精品久久久久久久人妻蜜臀av| 特级一级黄色大片| 黄色日韩在线| 国产精品.久久久| 日本免费在线观看一区| 国产黄色免费在线视频| 在线观看一区二区三区| 国产男人的电影天堂91| 成人漫画全彩无遮挡| 在线 av 中文字幕| 日韩欧美精品免费久久| 国产美女午夜福利| 欧美丝袜亚洲另类| 狠狠精品人妻久久久久久综合| 免费黄网站久久成人精品| 免费黄频网站在线观看国产| 中文字幕av在线有码专区| 少妇熟女aⅴ在线视频| 看十八女毛片水多多多| 边亲边吃奶的免费视频| 国产精品综合久久久久久久免费| 五月伊人婷婷丁香| 在线免费观看的www视频| 亚洲高清免费不卡视频| 亚洲最大成人手机在线| 极品教师在线视频| 亚洲精品一区蜜桃| 免费看光身美女| 国产老妇女一区| 我的老师免费观看完整版| 日韩欧美精品v在线| 亚洲av男天堂| 国产成人精品婷婷| 色播亚洲综合网| 在线免费观看不下载黄p国产| 日韩欧美国产在线观看| 久久久欧美国产精品| 亚洲成色77777| 男的添女的下面高潮视频| 久久久久久伊人网av| 国产高清有码在线观看视频| 亚洲综合精品二区| 99热全是精品| 亚洲美女搞黄在线观看| 国产高清三级在线| 成人综合一区亚洲| 亚洲经典国产精华液单| 国产高潮美女av| 亚洲国产精品专区欧美| 淫秽高清视频在线观看| 久久6这里有精品| 插阴视频在线观看视频| 亚洲精品第二区| 99久久人妻综合| 亚洲人成网站在线观看播放| 白带黄色成豆腐渣| 国产高潮美女av| 欧美bdsm另类| 91久久精品国产一区二区三区| 亚洲aⅴ乱码一区二区在线播放| 国产人妻一区二区三区在| 午夜精品国产一区二区电影 | 18禁在线无遮挡免费观看视频| 国产中年淑女户外野战色| 99久久人妻综合| 欧美人与善性xxx| 夫妻性生交免费视频一级片| 91aial.com中文字幕在线观看| 国产成人精品婷婷| 亚洲婷婷狠狠爱综合网| 麻豆av噜噜一区二区三区| 国产成人精品一,二区| 国产精品一区二区三区四区免费观看| 国产女主播在线喷水免费视频网站 | 久久久久精品久久久久真实原创| 国产成人福利小说| av在线天堂中文字幕| 两个人视频免费观看高清| 亚洲人与动物交配视频| 亚洲四区av| 国产69精品久久久久777片| 麻豆成人午夜福利视频| 观看免费一级毛片| 一个人看的www免费观看视频| 少妇被粗大猛烈的视频| 免费大片18禁| 午夜爱爱视频在线播放| 午夜精品一区二区三区免费看| .国产精品久久| 日本色播在线视频| 舔av片在线| 久久国内精品自在自线图片| 又黄又爽又刺激的免费视频.| 亚洲av二区三区四区| 亚洲av.av天堂| 一夜夜www| 亚洲精品中文字幕在线视频 | 国产色婷婷99| 网址你懂的国产日韩在线| av免费观看日本| 国产亚洲午夜精品一区二区久久 | 能在线免费观看的黄片| 中文字幕人妻熟人妻熟丝袜美| 免费无遮挡裸体视频| 亚洲成人中文字幕在线播放| 精品久久久久久久久av| 国产精品久久久久久精品电影小说 | 蜜桃久久精品国产亚洲av| 日本熟妇午夜| 黄片wwwwww| 日本av手机在线免费观看| 欧美丝袜亚洲另类| a级一级毛片免费在线观看| 日日啪夜夜撸| 又粗又硬又长又爽又黄的视频| 插阴视频在线观看视频| 久久精品国产鲁丝片午夜精品| 精品少妇黑人巨大在线播放| 久久99蜜桃精品久久| 淫秽高清视频在线观看| 少妇人妻精品综合一区二区| 女的被弄到高潮叫床怎么办| 久久亚洲国产成人精品v| 欧美xxⅹ黑人| 亚洲av福利一区| 国语对白做爰xxxⅹ性视频网站| 最近的中文字幕免费完整| 男女那种视频在线观看| 伦理电影大哥的女人| 日产精品乱码卡一卡2卡三| 国产亚洲91精品色在线| 少妇猛男粗大的猛烈进出视频 | 18禁动态无遮挡网站| 欧美人与善性xxx| 国产黄频视频在线观看| 欧美性猛交╳xxx乱大交人| 午夜爱爱视频在线播放| 欧美97在线视频| 国产高潮美女av| 九色成人免费人妻av| 日韩亚洲欧美综合| 成人一区二区视频在线观看| av.在线天堂| 国产精品久久久久久精品电影| 观看美女的网站| 国产乱来视频区| 久久这里有精品视频免费| av网站免费在线观看视频 | eeuss影院久久| 亚洲经典国产精华液单| 内地一区二区视频在线| 亚洲最大成人av| 在线a可以看的网站| 国产视频内射| 国产亚洲一区二区精品| 免费黄频网站在线观看国产| 精品国产一区二区三区久久久樱花 | 欧美另类一区| 日本欧美国产在线视频| 成人一区二区视频在线观看| 亚洲精品aⅴ在线观看| 一边亲一边摸免费视频| 亚洲国产精品成人久久小说| 午夜精品一区二区三区免费看| 亚洲欧美一区二区三区黑人 | 身体一侧抽搐| 亚州av有码| 禁无遮挡网站| 十八禁国产超污无遮挡网站| 国产一级毛片七仙女欲春2| 床上黄色一级片| 亚洲av免费高清在线观看| 美女被艹到高潮喷水动态| 久久久久久九九精品二区国产| 欧美另类一区| 少妇人妻精品综合一区二区| 欧美日本视频| 免费高清在线观看视频在线观看| 深爱激情五月婷婷| 国产精品一二三区在线看| 欧美激情久久久久久爽电影| 久久精品久久久久久噜噜老黄| av卡一久久| 久久精品国产鲁丝片午夜精品| 九九在线视频观看精品| 亚洲av.av天堂| 一本一本综合久久| 黄色欧美视频在线观看| 91午夜精品亚洲一区二区三区| 国产高清国产精品国产三级 | 国产一区二区在线观看日韩| 毛片一级片免费看久久久久| 中文字幕亚洲精品专区| 精品酒店卫生间| 色综合亚洲欧美另类图片| 丝袜美腿在线中文| 精品人妻一区二区三区麻豆| 欧美区成人在线视频| 美女黄网站色视频| 国产亚洲av嫩草精品影院| 成人欧美大片| 国产黄色免费在线视频| 午夜亚洲福利在线播放| 女人被狂操c到高潮| 国产男女超爽视频在线观看| 色吧在线观看| 99久久精品国产国产毛片| 联通29元200g的流量卡| 人妻夜夜爽99麻豆av| 国产成人精品婷婷| 免费在线观看成人毛片| 亚洲性久久影院| 日韩伦理黄色片| 大片免费播放器 马上看| videos熟女内射| 久久久久久九九精品二区国产| 舔av片在线| 亚洲精品国产av蜜桃| 色网站视频免费| 国产精品三级大全| 日韩一区二区三区影片| 久久精品国产亚洲av涩爱| 久久精品夜夜夜夜夜久久蜜豆| 七月丁香在线播放| 日韩强制内射视频| 国产一级毛片在线| 欧美日韩综合久久久久久| 国产男人的电影天堂91| 人妻系列 视频| 久久久久国产网址| 黑人高潮一二区| 午夜爱爱视频在线播放| 久久久久免费精品人妻一区二区| 91aial.com中文字幕在线观看| 日韩av在线大香蕉| 国产大屁股一区二区在线视频| 中文天堂在线官网| 综合色丁香网| 亚洲av电影不卡..在线观看| 色综合色国产| 日本wwww免费看| 久久久久精品性色| 日日摸夜夜添夜夜添av毛片| 国内揄拍国产精品人妻在线| 亚洲精品国产av成人精品| 尤物成人国产欧美一区二区三区| 亚洲国产高清在线一区二区三| 国产精品爽爽va在线观看网站| 晚上一个人看的免费电影| 亚洲国产高清在线一区二区三| 一级a做视频免费观看| 亚洲精品中文字幕在线视频 | 精品国内亚洲2022精品成人| 午夜福利高清视频| 一个人免费在线观看电影| 男人狂女人下面高潮的视频| 最近中文字幕2019免费版| 久久久久久伊人网av| 成人性生交大片免费视频hd| 亚洲在久久综合| 国产精品av视频在线免费观看| 能在线免费观看的黄片| 内射极品少妇av片p| 精品午夜福利在线看| 韩国av在线不卡| 亚洲无线观看免费| 免费观看a级毛片全部| 国产在视频线精品| 黄色欧美视频在线观看| 日韩欧美精品v在线| 一个人观看的视频www高清免费观看| 国产色爽女视频免费观看| 亚洲精品视频女| 成人午夜高清在线视频| 久久草成人影院| 国产精品一区二区性色av| 我要看日韩黄色一级片| 高清毛片免费看| 高清午夜精品一区二区三区| 国产视频首页在线观看| 如何舔出高潮| 欧美三级亚洲精品| 在线观看美女被高潮喷水网站| 国产精品蜜桃在线观看| 三级国产精品欧美在线观看| 免费大片18禁| 尤物成人国产欧美一区二区三区| 麻豆成人av视频| 亚洲av免费在线观看| 人妻系列 视频| 久久久久久久久久人人人人人人| 99九九线精品视频在线观看视频| 国产综合懂色| 国产乱人视频| 一二三四中文在线观看免费高清| 成人综合一区亚洲| av免费观看日本| 久久久久网色| h日本视频在线播放| 夫妻午夜视频| 91精品国产九色| 大话2 男鬼变身卡| 国产黄色小视频在线观看| 日韩欧美国产在线观看| 久久精品人妻少妇| 一区二区三区乱码不卡18| 日韩成人伦理影院| 国产成人免费观看mmmm| 国产视频内射| 成人亚洲欧美一区二区av| 久久综合国产亚洲精品| 亚洲久久久久久中文字幕| 中文天堂在线官网| 色播亚洲综合网| 亚洲精华国产精华液的使用体验| 免费电影在线观看免费观看| 美女高潮的动态| 精品一区二区三区人妻视频| 国产精品精品国产色婷婷| 99热这里只有精品一区| 麻豆久久精品国产亚洲av| 一级毛片久久久久久久久女| 亚洲三级黄色毛片| 久久久成人免费电影| 久久精品国产亚洲网站| 免费观看在线日韩| 久久精品久久精品一区二区三区| 亚洲在线自拍视频| 亚洲精品成人久久久久久| 亚洲国产精品成人久久小说| 白带黄色成豆腐渣| 日韩强制内射视频| 精品国内亚洲2022精品成人| 日韩不卡一区二区三区视频在线| 久久久久久久午夜电影| 亚洲婷婷狠狠爱综合网| 少妇的逼水好多| 日韩不卡一区二区三区视频在线| 秋霞在线观看毛片| 日本黄大片高清| 久久精品久久精品一区二区三区| 秋霞伦理黄片| 高清日韩中文字幕在线| 少妇熟女欧美另类| 肉色欧美久久久久久久蜜桃 | 亚洲国产欧美人成| 日日啪夜夜爽| 汤姆久久久久久久影院中文字幕 | 久久精品国产鲁丝片午夜精品| 国产精品福利在线免费观看| 人妻少妇偷人精品九色| 国产精品.久久久| 久久精品综合一区二区三区| 真实男女啪啪啪动态图| 久久精品久久久久久噜噜老黄| 99热6这里只有精品| 国产熟女欧美一区二区| 久久久久久国产a免费观看| 伊人久久国产一区二区| 免费高清在线观看视频在线观看| 禁无遮挡网站| 国产一区二区三区综合在线观看 | 免费在线观看成人毛片| 麻豆久久精品国产亚洲av| 91精品伊人久久大香线蕉| 国产伦精品一区二区三区四那| av国产免费在线观看| 午夜免费激情av| 成人高潮视频无遮挡免费网站| 亚洲国产av新网站| 搡老乐熟女国产| 蜜桃久久精品国产亚洲av| 在线天堂最新版资源| 国产高清有码在线观看视频| 亚洲精品一二三| 天堂√8在线中文| 中文字幕人妻熟人妻熟丝袜美| 少妇人妻一区二区三区视频| 大香蕉97超碰在线| 一级毛片黄色毛片免费观看视频| 久久韩国三级中文字幕| 搞女人的毛片| 激情五月婷婷亚洲| 欧美极品一区二区三区四区| 国产 一区精品| 精品人妻熟女av久视频| 麻豆久久精品国产亚洲av| 午夜免费激情av| av福利片在线观看| 久久久精品94久久精品| 日韩欧美一区视频在线观看 | 一级毛片 在线播放| 精品99又大又爽又粗少妇毛片| 亚洲精品日韩av片在线观看| 成人性生交大片免费视频hd| 精品熟女少妇av免费看| 国产精品久久视频播放| 人妻制服诱惑在线中文字幕| 欧美日韩视频高清一区二区三区二| 国产精品人妻久久久久久| 激情五月婷婷亚洲| 少妇的逼水好多| 中国国产av一级| 国产精品不卡视频一区二区| 99视频精品全部免费 在线| 欧美成人一区二区免费高清观看| av福利片在线观看| 日韩av在线大香蕉| 免费看光身美女| 国产精品一区二区三区四区久久| 欧美日韩亚洲高清精品| 身体一侧抽搐| 99久国产av精品国产电影| 国产又色又爽无遮挡免| or卡值多少钱| 大又大粗又爽又黄少妇毛片口| 亚洲av成人精品一二三区| 久久久久久久国产电影| 夜夜看夜夜爽夜夜摸| 日本一二三区视频观看| 欧美最新免费一区二区三区| 亚洲婷婷狠狠爱综合网| 老司机影院毛片| 日韩欧美精品免费久久| 国产熟女欧美一区二区| 毛片女人毛片| 国产精品综合久久久久久久免费| 久久国内精品自在自线图片| 亚洲第一区二区三区不卡| 日本欧美国产在线视频| 你懂的网址亚洲精品在线观看| 一区二区三区免费毛片| 九九爱精品视频在线观看| 亚洲怡红院男人天堂| 高清在线视频一区二区三区| 国产成人免费观看mmmm| 精品亚洲乱码少妇综合久久| 精品一区二区免费观看| 国产色爽女视频免费观看| 亚洲av福利一区| 亚洲三级黄色毛片| 日韩强制内射视频| 少妇的逼水好多| 久久综合国产亚洲精品| 精品少妇黑人巨大在线播放| 日日干狠狠操夜夜爽| 国产一区二区三区综合在线观看 | 免费黄频网站在线观看国产| 最近最新中文字幕免费大全7| 99热这里只有是精品50| 国国产精品蜜臀av免费| 久久久a久久爽久久v久久| 国内揄拍国产精品人妻在线| 高清在线视频一区二区三区| 国产国拍精品亚洲av在线观看| 久久久欧美国产精品| 欧美一区二区亚洲| av在线亚洲专区| 久久99热这里只频精品6学生| 亚洲成人av在线免费| 亚洲精品日本国产第一区| 国产亚洲5aaaaa淫片| 亚洲四区av| 日韩在线高清观看一区二区三区| 国产综合精华液| 国产在视频线在精品| 草草在线视频免费看| 中文字幕久久专区| 色哟哟·www| 欧美高清成人免费视频www| 偷拍熟女少妇极品色| 青春草视频在线免费观看| 久久久精品欧美日韩精品| 成年女人在线观看亚洲视频 | 久久人人爽人人爽人人片va| 91精品一卡2卡3卡4卡| 最近中文字幕高清免费大全6| 亚洲成色77777| 亚洲欧美一区二区三区国产| 国产亚洲5aaaaa淫片| 九九爱精品视频在线观看| 国产成人91sexporn| 极品少妇高潮喷水抽搐| 国产精品日韩av在线免费观看| 97超碰精品成人国产| 舔av片在线| 99久久精品一区二区三区| 美女xxoo啪啪120秒动态图| 伦理电影大哥的女人| 婷婷色综合www|