The mass-front velocities of granular flows results from the joint action of particle size gradations and the underlying surfaces.However,because of the complexity of friction during flow movement,details such as the ...The mass-front velocities of granular flows results from the joint action of particle size gradations and the underlying surfaces.However,because of the complexity of friction during flow movement,details such as the slope-toe impedance effects and momentum-transfer mechanisms have not been completely explained by theoretical analyses,numerical simulations,or field investigations.To study the mass-front velocity of dry granular flows influenced by the angle of the slope to the runout plane and particle size gradations we conducted model experiments that recorded the motion of rapid and long-runout rockslides or avalanches.Flume tests were conducted using slope angles of 25°,35°,45°,and 55° and three particle size gradations.The resulting mass-front motions consisted of three stages:acceleration,velocity maintenance,and deceleration.The existing methods of velocity prediction could not explain the slowing effect of the slope toe or the momentum-transfer steady velocity stage.When the slope angle increased from 25° to 55°,the mass-front velocities dropped significantly to between 44.4% and59.6% of the peak velocities and energy lossesincreased from 69.1% to 83.7% of the initial,respectively.The velocity maintenance stages occurred after the slope-toe and mass-front velocity fluctuations.During this stage,travel distances increased as the angles increased,but the average velocity was greatest at 45°.At a slope angle of 45°,as the median particle size increased,energy loss around the slope toe decreased,the efficiency of momentum transfer increased,and the distance of the velocity maintenance stage increased.We presented an improved average velocity formula for granular flow and a geometrical model of the energy along the flow line.展开更多
This study was conducted to investigate the flow field characteristics of right-angled flow passage with various cavities in the typical hydraulic manifold block.A low-speed visualization test rig was developed and th...This study was conducted to investigate the flow field characteristics of right-angled flow passage with various cavities in the typical hydraulic manifold block.A low-speed visualization test rig was developed and the flow field of the right-angled flow passage with different cavity structures was measured using 2D-PIV technique.Numerical model was established to simulate the three-dimensional flow field.Seven eddy viscosity turbulence models were investigated in predicting the flow field by comparing against the particle image relocimetry(PIV)measurement results.By defining the weight error function K,the S-A model was selected as the appropriate turbulence model.Then,a three-factor,three-level response surface numerical test was conducted to investigate the influence of flow passage connection type,cavity diameter and cavity length-diameter ratio on pressure loss.The results show that the Box-Benhnken Design(BBD)model can predict the total pressure loss accurately.The optimal factor level appeared in flow passage connection type II,14.64 mm diameter and 67.53%cavity length-diameter ratio.The total pressure loss decreased by 11.15%relative to the worst factor level,and total pressure loss can be reduced by 64.75%when using an arc transition right-angled flow passage,which indicates a new direction for the optimization design of flow passage in hydraulic manifold blocks.展开更多
The rheological properties of most liquid in nature are between liquids and solids, including both elastic changes and viscosity changes, that is socalled "viscoelastic". Dynamic oscillatory test was used to quantit...The rheological properties of most liquid in nature are between liquids and solids, including both elastic changes and viscosity changes, that is socalled "viscoelastic". Dynamic oscillatory test was used to quantitatively study the distinct viscoelastic behaviors of debris flow slurry in the shear stress conditions for the first time in this study. The debris flow slurry samples were from Jiangjiagou Ravine, Yunnan Province, China. The experimental results were found that at the low and middle stages of shearing, when the angular velocity 09〈72.46 s-1, the loss modulus (G was greater than the storage modulus (G3, i.e. G"〉G'. At the late stage of shearing, when the angular velocity co-72.46 s-x, the storage modulus was greater than or equal to the loss G = G, tan -〈 1 (where phase-shift modulus, i.e. ' 〉 " angle 5=G",/G3, and the debris flow slurry was in a gel state. Therefore, the progress of this experimental study further reveals the mechanism of hyperconcentrated debris flows with a high velocity on low-gradient ravines.展开更多
Long-term trends of yearly and seasonal averages of tropospheric ozone over the whole country and some important regions of China during 1979-2005 are analyzed,based on tropospheric ozone residue(TOR) data retrieved f...Long-term trends of yearly and seasonal averages of tropospheric ozone over the whole country and some important regions of China during 1979-2005 are analyzed,based on tropospheric ozone residue(TOR) data retrieved from satellite measurements.The relationship between the TOR and Southern Oscillation Index(SOI) is studied.The results show that,over the whole country,there is a slight increase of TOR in summer and a decrease in other seasons,while the overall trend for the whole period is insignificant.There are decreasing trends of TOR over the Pearl River Delta and the Sichuan Basin.Significant increasing trends of TOR are found over the North China Plain(NCP) for all seasons except for winter,with a maximum rate of 1.10 DU per decade for summer.There are significant correlations between TOR and SOI for some regions in China but not for the NCP,suggesting that the observed increasing trend of TOR over the NCP may not be linked with changes in atmospheric circulations.展开更多
The international standard IEC 61400-12-1 Wind turbines--Part 12-1: Power performance measurements of electricity producing wind turbines" aims to provide a uniform methodology that will ensure consistency, accuracy...The international standard IEC 61400-12-1 Wind turbines--Part 12-1: Power performance measurements of electricity producing wind turbines" aims to provide a uniform methodology that will ensure consistency, accuracy and reproducibility in the measurement and analysis of power performance by wind turbines. Annex G of this standard provides a methodology for the appropriate arrangement of instruments on the meteorological mast to ensure accurate measurement. For cup anemometers it provides recommendations about their location relative to the mast so that the effect of mast and boom interference on their output may be minimised. These recommendations are given for both tubular masts and lattice masts. This paper compares the flow distortion predicted by the IEC standard and the results of a 3D CFD (computational fluid dynamics) simulation of a triangular lattice mast. Based on the results of wind tunnel and CFD simulation it was found that the flow distortion surrounding the lattice mast was overpredicted by the method suggested in appendix G oflEC61400-12-1. Using the CFD data it was possible to determine, for a range of flow directions and mast heights, the distance from the mast that anemometers would need to be in order to be outside the flow distortion field.展开更多
Vortex flow is applied to a cyclone dust collector, a vortex combustion chamber, and a vortex diode for vortex control. In order to apply the vortex flow to the industries, it is necessary to keep the stable flow cond...Vortex flow is applied to a cyclone dust collector, a vortex combustion chamber, and a vortex diode for vortex control. In order to apply the vortex flow to the industries, it is necessary to keep the stable flow condition and to estimate the response time of the transient flow process and also the intensity of the vortex flow. For control vortex flow, two types of vortex chamber with two inlet pipes were designed. One of them is to promote the vortex flow named as Co-Rotating Flow System and another one is to hinder the vortex flow named as Counter-Rotating Flow System. The pressure drops and the velocity distributions were measured for these vortex chambers. The estimation of the tangential velocity by the application of the angular momentum flux is compared with the measured velocity by a cylindrical Pitot-tube. The characteristics of the total pressure drop could be explained by introducing the circulation.展开更多
This study is aimed at using blade 3-D optimization to control corner flows in the high through-flow fan/booster of a high bypass ratio commercial turbofan engine. Two kinds of blade 3-D optimization, end-bending and ...This study is aimed at using blade 3-D optimization to control corner flows in the high through-flow fan/booster of a high bypass ratio commercial turbofan engine. Two kinds of blade 3-D optimization, end-bending and bow, are focused on. On account of the respective operation mode and environment, the approach to 3-D aerodynamic modeling of rotor blades is different from stator vanes. Based on the understanding of the mechanism of the corner flow and the consideration of intensity problem for rotors, this paper uses a variety of blade 3-D optimization approaches, such as loading distribution optimization, perturbation of departure angles and stacking-axis manipulation, which are suitable for rotors and stators respectively. The obtained 3-D blades and vanes can improve the corner flow features by end-bending and bow effects. The results of this study show that flows in corners of the fan/booster, such as the fan hub region, the tip and hub of the vanes of the booster, are very complex and dominated by 3-D effects. The secondary flows there are found to have a strong detrimental effect on the compressor performance. The effects of both end-bending and bow can improve the flow separation in corners, but the specific ways they work and application scope are somewhat different. Redesigning the blades via blade 3-D optimization to control the corner flow has effectively reduced the loss generation and improved the stall margin by a large amount.展开更多
This paper mainly deals with the type II singularities of the mean curvature flow from a symplectic surface or from an almost calibrated Lagrangian surface in a K¨ahler surface.The relation between the maximum of...This paper mainly deals with the type II singularities of the mean curvature flow from a symplectic surface or from an almost calibrated Lagrangian surface in a K¨ahler surface.The relation between the maximum of the Kahler angle and the maximum of |H|2 on the limit flow is studied.The authors also show the nonexistence of type II blow-up flow of a symplectic mean curvature flow which is normal flat or of an almost calibrated Lagrangian mean curvature flow which is flat.展开更多
The present study experimentally investigated the near-field flow mixing characteristics of two turbulent jets issuing from equilateral triangular and circular orifice plates into effectively unbounded surroundings,re...The present study experimentally investigated the near-field flow mixing characteristics of two turbulent jets issuing from equilateral triangular and circular orifice plates into effectively unbounded surroundings,respectively.Planar particle image velocimetry(PIV) was applied to measure the velocity field at the same Reynolds number of Re=50,000,where Re = UeDe /with Ue being the exit bulk velocity and the kinematic viscosity of fluid,D e the equivalent diameters.The instantaneous velocity,mean velocity,Reynolds stresses were obtained.From the mean velocity field,the centreline velocity decay rate and half-velocity width were derived.Comparing the mixing characteristics of the two jets,it is found that the triangular jet has a faster mixing rate than the circular counterpart.The triangular jet entrainments with the ambient fluid at a higher rate in the near field.This is evidenced by a shorter unmixed core,faster Reynolds stress and centreline turbulence intensity growth.The primary coherent structures in the near field are found to break down more rapidly in the triangular jet as compared to the circular jet.Over the entire measurement region,the triangular jet maintained a higher rate of decay and spread.Moreover,all components of Reynolds stress of the triangular jet appear to reach their peaks earlier,and then decay more rapidly than those of the circular jet.In addition,the axis-switching phenomenon is observed in the triangular jet.展开更多
Many researches on critical nozzles have been performed to accurately measure the mass flow rate of gas flow,and to standardize the performance as a flow meter.Recently,much interest is being paid on the measurement o...Many researches on critical nozzles have been performed to accurately measure the mass flow rate of gas flow,and to standardize the performance as a flow meter.Recently,much interest is being paid on the measurement of very small mass flow rate in industry fields such as MEMS applications.However,the design and performance data of the critical nozzles obtained so far have been applied mainly to the critical nozzles with comparatively large diameters,and the works available on miniature critical nozzles are lacking.In the present study,a computational fluid dynamics method has been applied to investigate the influence of the diffuser angle on discharge coefficient of the miniature critical nozzles.In computations,the throat diameter of critical nozzle is varied from 0.2 mm to 5.0 mm and the diffuser angle is changed from 2 deg to 8 deg.The computational results are validated with some experimental data available.The results show that the present computational results predict appropriately the discharge coefficient of the gas flows through miniature critical nozzles.It is known that the discharge coefficient is considerably influenced by the diffuser angle,as the throat diameter of nozzle becomes small below a certain value.This implies that the miniature critical nozzles should be carefully designed.展开更多
In this paper,firstly,by solving the Riemann problem of the zero-pressure flow in gas dynamics with a flux approximation,we construct parameterized delta-shock and constant density solutions,then we show that,as the f...In this paper,firstly,by solving the Riemann problem of the zero-pressure flow in gas dynamics with a flux approximation,we construct parameterized delta-shock and constant density solutions,then we show that,as the flux perturbation vanishes,they converge to the delta-shock and vacuum state solutions of the zero-pressure flow,respectively.Secondly,we solve the Riemann problem of the Euler equations of isentropic gas dynamics with a double parameter flux approximation including pressure.Furthermore,we rigorously prove that,as the two-parameter flux perturbation vanishes,any Riemann solution containing two shock waves tends to a delta-shock solution to the zero-pressure flow;any Riemann solution containing two rarefaction waves tends to a two-contact-discontinuity solution to the zero-pressure flow and the nonvacuum intermediate state in between tends to a vacuum state.Finally,numerical results are given to present the formation processes of delta shock waves and vacuum states.展开更多
The natural convective heat transfer performance and thermo-fluidic characteristics of honeycombs with/without chimney extensions are numerically investigated.The present numerical simulations are validated by the pur...The natural convective heat transfer performance and thermo-fluidic characteristics of honeycombs with/without chimney extensions are numerically investigated.The present numerical simulations are validated by the purposely-designed experimental measurements on honeycombs with/without chimney.Good agreement between numerical simulation and experimental measurement is obtained.The influences of inclination angle and geometric parameters such as cell shape,streamwise and spanwise length are also numerically quantified.With the increment in inclination angle,the overall heat transfer rate decreases for the honeycombs with/without chimney.For honeycombs with the same void volume fraction but different cell shapes,there is little difference on the overall heat transfer rate.To enhance the natural convective heat transfer of honeycombs,these techniques including increasing the length of honeycomb in the streamwise/spanwise direction,increasing the thermal conductivity of hon-eycomb structure or adding a chimney extension may be helpful.展开更多
基金supported by the National Natural Science Foundation of China (Grant Nos.41272297,41401195)the Applied Basic Research Fund of the Science and Technology Department of Sichuan Province (2014JY0121)the Key Research Fund of the Education Department of Sichuan Province (14ZA0095)
文摘The mass-front velocities of granular flows results from the joint action of particle size gradations and the underlying surfaces.However,because of the complexity of friction during flow movement,details such as the slope-toe impedance effects and momentum-transfer mechanisms have not been completely explained by theoretical analyses,numerical simulations,or field investigations.To study the mass-front velocity of dry granular flows influenced by the angle of the slope to the runout plane and particle size gradations we conducted model experiments that recorded the motion of rapid and long-runout rockslides or avalanches.Flume tests were conducted using slope angles of 25°,35°,45°,and 55° and three particle size gradations.The resulting mass-front motions consisted of three stages:acceleration,velocity maintenance,and deceleration.The existing methods of velocity prediction could not explain the slowing effect of the slope toe or the momentum-transfer steady velocity stage.When the slope angle increased from 25° to 55°,the mass-front velocities dropped significantly to between 44.4% and59.6% of the peak velocities and energy lossesincreased from 69.1% to 83.7% of the initial,respectively.The velocity maintenance stages occurred after the slope-toe and mass-front velocity fluctuations.During this stage,travel distances increased as the angles increased,but the average velocity was greatest at 45°.At a slope angle of 45°,as the median particle size increased,energy loss around the slope toe decreased,the efficiency of momentum transfer increased,and the distance of the velocity maintenance stage increased.We presented an improved average velocity formula for granular flow and a geometrical model of the energy along the flow line.
基金Projects(51705446,51890881) supported by the National Natural Science Foundation of China
文摘This study was conducted to investigate the flow field characteristics of right-angled flow passage with various cavities in the typical hydraulic manifold block.A low-speed visualization test rig was developed and the flow field of the right-angled flow passage with different cavity structures was measured using 2D-PIV technique.Numerical model was established to simulate the three-dimensional flow field.Seven eddy viscosity turbulence models were investigated in predicting the flow field by comparing against the particle image relocimetry(PIV)measurement results.By defining the weight error function K,the S-A model was selected as the appropriate turbulence model.Then,a three-factor,three-level response surface numerical test was conducted to investigate the influence of flow passage connection type,cavity diameter and cavity length-diameter ratio on pressure loss.The results show that the Box-Benhnken Design(BBD)model can predict the total pressure loss accurately.The optimal factor level appeared in flow passage connection type II,14.64 mm diameter and 67.53%cavity length-diameter ratio.The total pressure loss decreased by 11.15%relative to the worst factor level,and total pressure loss can be reduced by 64.75%when using an arc transition right-angled flow passage,which indicates a new direction for the optimization design of flow passage in hydraulic manifold blocks.
基金supported by the Youth Talent Team Program of Institute of Mountain Hazards and Environment,CASthe National Natural Science Foundation of China (Grant No.406710260)
文摘The rheological properties of most liquid in nature are between liquids and solids, including both elastic changes and viscosity changes, that is socalled "viscoelastic". Dynamic oscillatory test was used to quantitatively study the distinct viscoelastic behaviors of debris flow slurry in the shear stress conditions for the first time in this study. The debris flow slurry samples were from Jiangjiagou Ravine, Yunnan Province, China. The experimental results were found that at the low and middle stages of shearing, when the angular velocity 09〈72.46 s-1, the loss modulus (G was greater than the storage modulus (G3, i.e. G"〉G'. At the late stage of shearing, when the angular velocity co-72.46 s-x, the storage modulus was greater than or equal to the loss G = G, tan -〈 1 (where phase-shift modulus, i.e. ' 〉 " angle 5=G",/G3, and the debris flow slurry was in a gel state. Therefore, the progress of this experimental study further reveals the mechanism of hyperconcentrated debris flows with a high velocity on low-gradient ravines.
基金supported by the Chinese National Science Foundation(40775074)China Meteorological Administration(GYHY(QX)200706005).
文摘Long-term trends of yearly and seasonal averages of tropospheric ozone over the whole country and some important regions of China during 1979-2005 are analyzed,based on tropospheric ozone residue(TOR) data retrieved from satellite measurements.The relationship between the TOR and Southern Oscillation Index(SOI) is studied.The results show that,over the whole country,there is a slight increase of TOR in summer and a decrease in other seasons,while the overall trend for the whole period is insignificant.There are decreasing trends of TOR over the Pearl River Delta and the Sichuan Basin.Significant increasing trends of TOR are found over the North China Plain(NCP) for all seasons except for winter,with a maximum rate of 1.10 DU per decade for summer.There are significant correlations between TOR and SOI for some regions in China but not for the NCP,suggesting that the observed increasing trend of TOR over the NCP may not be linked with changes in atmospheric circulations.
文摘The international standard IEC 61400-12-1 Wind turbines--Part 12-1: Power performance measurements of electricity producing wind turbines" aims to provide a uniform methodology that will ensure consistency, accuracy and reproducibility in the measurement and analysis of power performance by wind turbines. Annex G of this standard provides a methodology for the appropriate arrangement of instruments on the meteorological mast to ensure accurate measurement. For cup anemometers it provides recommendations about their location relative to the mast so that the effect of mast and boom interference on their output may be minimised. These recommendations are given for both tubular masts and lattice masts. This paper compares the flow distortion predicted by the IEC standard and the results of a 3D CFD (computational fluid dynamics) simulation of a triangular lattice mast. Based on the results of wind tunnel and CFD simulation it was found that the flow distortion surrounding the lattice mast was overpredicted by the method suggested in appendix G oflEC61400-12-1. Using the CFD data it was possible to determine, for a range of flow directions and mast heights, the distance from the mast that anemometers would need to be in order to be outside the flow distortion field.
文摘Vortex flow is applied to a cyclone dust collector, a vortex combustion chamber, and a vortex diode for vortex control. In order to apply the vortex flow to the industries, it is necessary to keep the stable flow condition and to estimate the response time of the transient flow process and also the intensity of the vortex flow. For control vortex flow, two types of vortex chamber with two inlet pipes were designed. One of them is to promote the vortex flow named as Co-Rotating Flow System and another one is to hinder the vortex flow named as Counter-Rotating Flow System. The pressure drops and the velocity distributions were measured for these vortex chambers. The estimation of the tangential velocity by the application of the angular momentum flux is compared with the measured velocity by a cylindrical Pitot-tube. The characteristics of the total pressure drop could be explained by introducing the circulation.
基金supported by National Natural Science Foundation of China (51006005,50736007)"Fan-Zhou" Youth Foundation(20100401)
文摘This study is aimed at using blade 3-D optimization to control corner flows in the high through-flow fan/booster of a high bypass ratio commercial turbofan engine. Two kinds of blade 3-D optimization, end-bending and bow, are focused on. On account of the respective operation mode and environment, the approach to 3-D aerodynamic modeling of rotor blades is different from stator vanes. Based on the understanding of the mechanism of the corner flow and the consideration of intensity problem for rotors, this paper uses a variety of blade 3-D optimization approaches, such as loading distribution optimization, perturbation of departure angles and stacking-axis manipulation, which are suitable for rotors and stators respectively. The obtained 3-D blades and vanes can improve the corner flow features by end-bending and bow effects. The results of this study show that flows in corners of the fan/booster, such as the fan hub region, the tip and hub of the vanes of the booster, are very complex and dominated by 3-D effects. The secondary flows there are found to have a strong detrimental effect on the compressor performance. The effects of both end-bending and bow can improve the flow separation in corners, but the specific ways they work and application scope are somewhat different. Redesigning the blades via blade 3-D optimization to control the corner flow has effectively reduced the loss generation and improved the stall margin by a large amount.
基金Project supported by the National Natural Science Foundation of China (Nos. 10901088, 11001268)
文摘This paper mainly deals with the type II singularities of the mean curvature flow from a symplectic surface or from an almost calibrated Lagrangian surface in a K¨ahler surface.The relation between the maximum of the Kahler angle and the maximum of |H|2 on the limit flow is studied.The authors also show the nonexistence of type II blow-up flow of a symplectic mean curvature flow which is normal flat or of an almost calibrated Lagrangian mean curvature flow which is flat.
基金the support of the Fundamental Research Funds for the Central Universities (Grant No. 3132013029)the National Natural Science Foundation of China (Grant Nos. 10921202 and11072005)
文摘The present study experimentally investigated the near-field flow mixing characteristics of two turbulent jets issuing from equilateral triangular and circular orifice plates into effectively unbounded surroundings,respectively.Planar particle image velocimetry(PIV) was applied to measure the velocity field at the same Reynolds number of Re=50,000,where Re = UeDe /with Ue being the exit bulk velocity and the kinematic viscosity of fluid,D e the equivalent diameters.The instantaneous velocity,mean velocity,Reynolds stresses were obtained.From the mean velocity field,the centreline velocity decay rate and half-velocity width were derived.Comparing the mixing characteristics of the two jets,it is found that the triangular jet has a faster mixing rate than the circular counterpart.The triangular jet entrainments with the ambient fluid at a higher rate in the near field.This is evidenced by a shorter unmixed core,faster Reynolds stress and centreline turbulence intensity growth.The primary coherent structures in the near field are found to break down more rapidly in the triangular jet as compared to the circular jet.Over the entire measurement region,the triangular jet maintained a higher rate of decay and spread.Moreover,all components of Reynolds stress of the triangular jet appear to reach their peaks earlier,and then decay more rapidly than those of the circular jet.In addition,the axis-switching phenomenon is observed in the triangular jet.
文摘Many researches on critical nozzles have been performed to accurately measure the mass flow rate of gas flow,and to standardize the performance as a flow meter.Recently,much interest is being paid on the measurement of very small mass flow rate in industry fields such as MEMS applications.However,the design and performance data of the critical nozzles obtained so far have been applied mainly to the critical nozzles with comparatively large diameters,and the works available on miniature critical nozzles are lacking.In the present study,a computational fluid dynamics method has been applied to investigate the influence of the diffuser angle on discharge coefficient of the miniature critical nozzles.In computations,the throat diameter of critical nozzle is varied from 0.2 mm to 5.0 mm and the diffuser angle is changed from 2 deg to 8 deg.The computational results are validated with some experimental data available.The results show that the present computational results predict appropriately the discharge coefficient of the gas flows through miniature critical nozzles.It is known that the discharge coefficient is considerably influenced by the diffuser angle,as the throat diameter of nozzle becomes small below a certain value.This implies that the miniature critical nozzles should be carefully designed.
基金supported by National Natural Science Foundation of China(Grant No.11361073)
文摘In this paper,firstly,by solving the Riemann problem of the zero-pressure flow in gas dynamics with a flux approximation,we construct parameterized delta-shock and constant density solutions,then we show that,as the flux perturbation vanishes,they converge to the delta-shock and vacuum state solutions of the zero-pressure flow,respectively.Secondly,we solve the Riemann problem of the Euler equations of isentropic gas dynamics with a double parameter flux approximation including pressure.Furthermore,we rigorously prove that,as the two-parameter flux perturbation vanishes,any Riemann solution containing two shock waves tends to a delta-shock solution to the zero-pressure flow;any Riemann solution containing two rarefaction waves tends to a two-contact-discontinuity solution to the zero-pressure flow and the nonvacuum intermediate state in between tends to a vacuum state.Finally,numerical results are given to present the formation processes of delta shock waves and vacuum states.
基金supported by the National 111 Project of China(Grant No.B06024)the National Basic Research Program of China("973"Project)(Grant No.2011CB610305)the National Natural Science Foundation of China(Grant No.51206128)
文摘The natural convective heat transfer performance and thermo-fluidic characteristics of honeycombs with/without chimney extensions are numerically investigated.The present numerical simulations are validated by the purposely-designed experimental measurements on honeycombs with/without chimney.Good agreement between numerical simulation and experimental measurement is obtained.The influences of inclination angle and geometric parameters such as cell shape,streamwise and spanwise length are also numerically quantified.With the increment in inclination angle,the overall heat transfer rate decreases for the honeycombs with/without chimney.For honeycombs with the same void volume fraction but different cell shapes,there is little difference on the overall heat transfer rate.To enhance the natural convective heat transfer of honeycombs,these techniques including increasing the length of honeycomb in the streamwise/spanwise direction,increasing the thermal conductivity of hon-eycomb structure or adding a chimney extension may be helpful.