Transition of flows past a pair of side-by-side circular cylinders are investigated by numerical simulations and the bifurcation analysis of the numerical results. Various flow patterns behind the cylinder-pair have b...Transition of flows past a pair of side-by-side circular cylinders are investigated by numerical simulations and the bifurcation analysis of the numerical results. Various flow patterns behind the cylinder-pair have been identified by the gap ratio (G) and Reynolds number (Re). This study focus on transition of in-phase and anti-phase vortex shedding synchronized forms. A nested Cartesian-grid formulation, in combination with an effective immersed boundary method and a two-step fractional-step procedure, has been adopted to simulate the flows. Numerical results reveal that the in-phase and anti-phase vortex shedding flows at Re = 100 can co-exist at 2.08 ≤G≤ 2.58. Hysteresis loop with increasing/decreasing G at constant Reynolds number Re = 100 is reported.展开更多
A robust Reynolds-Averaged Navier-Stokes(RANS)based solver is established to predict the complex unsteady aerodynamic characteristics of the Active Flap Control(AFC)rotor.The complex motion with multiple degrees of fr...A robust Reynolds-Averaged Navier-Stokes(RANS)based solver is established to predict the complex unsteady aerodynamic characteristics of the Active Flap Control(AFC)rotor.The complex motion with multiple degrees of freedom of the Trailing Edge Flap(TEF)is analyzed by employing an inverse nested overset grid method.Simulation of non-rotational and rotational modes of blade motion are carried out to investigate the formation and development of TEF shedding vortex with high-frequency deflection of TEF.Moreover,the mechanism of TEF deflection interference with blade tip vortex and overall rotor aerodynamics is also explored.In nonrotational mode,two bundles of vortices form at the gap ends of TEF and the main blade and merge into a single TEF vortex.Dynamic deflection of the TEF significantly interferes with the blade tip vortex.The position of the blade tip vortex consistently changes,and its frequency is directly related to the frequency of TEF deflection.In rotational mode,the tip vortex forms a helical structure.The end vortices at the gap sides co-swirl and subsequently merge into the concentrated beam of tip vortices,causing fluctuations in the vorticity and axial position of the tip vortex under the rotor.This research concludes with the investigation on suppression of Blade Vortex Interaction(BVI),showing an increase in miss distance and reduction in the vorticity of tip vortex through TEF phase control at a particular control frequency.Through this mechanism,a designed TEF deflection law increases the miss distance by 34.7%and reduces vorticity by 11.9%at the target position,demonstrating the effectiveness of AFC in mitigating BVI.展开更多
Vortex-acoustic coupling is one of the most important potential sources of combustion instability in solid rocket motors (SRMs). Based on the Von Karman Institute for Fluid Dynamics (VKI) experimental motor, the i...Vortex-acoustic coupling is one of the most important potential sources of combustion instability in solid rocket motors (SRMs). Based on the Von Karman Institute for Fluid Dynamics (VKI) experimental motor, the influence of the thermal inhibitor position and temperature on vortex-shedding-driven pressure oscillations is numerically studied via the large eddy simulation (LES) method. The simulation results demonstrate that vortex shedding is a periodic process and its accurate frequency can be numerically obtained. Acoustic modes could be easily excited by vortex shedding. The vortex shedding frequency and second acoustic frequency dominate the pressure oscillation characteristics in the chamber. Thermal inhibitor position and gas temperature have little effect on vortex shedding frequency, but have great impact on pressure oscillation amplitude. Pressure amplitude is much higher when the thermal inhibitor locates at the acoustic velocity anti-nodes. The farther the thermal inhibitor is to the nozzle head, the more vortex energy would be dissipated by the turbulence. Therefore, the vortex shedding amplitude at the second acoustic velocity antinode near 3/4L (L is chamber length) is larger than those of others. Besides, the natural acoustic frequencies increase with the gas temperature. As the vortex shedding frequency departs from the natural acoustic frequency, the vortex-acoustic feedback loop is decoupled. Consequently, both the vortex shedding and acoustic amplitudes decrease rapidly.展开更多
The flow past a circular-section cylinder with a conic shroud perforated with four holes at the peak was simulated numerically at Re = 100, considering two factors, viz.the angle of attack and the diameter of the hole...The flow past a circular-section cylinder with a conic shroud perforated with four holes at the peak was simulated numerically at Re = 100, considering two factors, viz.the angle of attack and the diameter of the holes. The effects of the perforated conic shroud on the vortex shedding pattern in the near wake was mainly investigated, as well as the time history of the drag and lift forces. In the investigated parameter space, three flow regimes were generally identified, corresponding to weak, moderate, and strong disturbance effects.In regime I, the wake can mainly be described by alternately shedding Kármán or Kármán-like vortices. In regime II, the spanwise vortices are obviously disturbed along the span due to the appearance of additional vorticity components and their interactions with the spanwise vortices, but still shed in synchronization along the spanwise direction. In regime III,the typical Kármán vortices partially or totally disappear,and some new vortex shedding patterns appear, such as-type, obliquely shedding, and crossed spanwise vortices with opposite sign. Corresponding to these complex vortex shedding patterns in the near wake, the fluid forces no longer oscillate regularly at a single vortex shedding frequency, but rather with a lower modulation frequency and multiple amplitudes. An overview of these flow regimes is presented.展开更多
Vortex street flowmeter has been used in steady flow measurement for about three decades The benefits of this type of flowmeter include high accuracy,good linearty,wide measuring range,and excellent reliability Ho...Vortex street flowmeter has been used in steady flow measurement for about three decades The benefits of this type of flowmeter include high accuracy,good linearty,wide measuring range,and excellent reliability However,in unsteady flow measurement,the pressure disturbance as well as the noise from the system or surrounding can reduce the signal to noise ratio of the flowmeter seriously Aimed to use vortex street flowmeters in unsteady flow measurement,the characteristics of the vortex shedding induced hydrodynamic vibration around the prism bluff body in a vortex street flowmeter are investigated numerically and by expriments The results show that the hydrodynamic vibrations with 180° phase shift occur at the axisymmetric points of the channel around the bluff body The most intense vibration occurs at the points on the lateral faces close to the base of the prism The results provide therefore a useful reference for developing an anti interference vortex flowmeter using the different ial sensing technique.展开更多
In this paper we have made a numerical study on the control of vortex shedding and drag reduction of a cylinder by attaching thin splitter plates. The wake structure of the cylinder of square cross-section with attach...In this paper we have made a numerical study on the control of vortex shedding and drag reduction of a cylinder by attaching thin splitter plates. The wake structure of the cylinder of square cross-section with attached splitter plates is analyzed for a range of Reynolds number, based on the incident stream and height of the cylinder, in the laminar range. The Navier-Stokes equations governing the flow are solved by the control volume method over a staggered grid arrangement. We have used the semi-implicit method for pressure-linked equation (SIMPLE) algorithm for computation. Our results show that the presence of a splitter plate upstream of the cylinder reduces the drag, but it has a small impact on the vortex shedding frequency when the plate length is beyond 1.5 time the height of the cylinder. The presence of a downstream splitter plate dampens the vortex shedding frequency. The entrainment of fluid into the inner side of the separated shear layers is obstructed by the downstream splitter plate. Our results suggest that by attaching in-line splitter plates both upstream and downstream of the cylinder, the vortex shedding can be suppressed, as well as a reduction in drag be obtained. We made a parametric study to determine the optimal length of these splitter plates so as to achieve low drag and low vortex shedding frequency.展开更多
An attempt has been made to explore whether the power relation can be obtained from theoretical considerations. The classical laminar and turbulent boundary layer concepts have been employed to determine appropriate v...An attempt has been made to explore whether the power relation can be obtained from theoretical considerations. The classical laminar and turbulent boundary layer concepts have been employed to determine appropriate values of the scaling lengths associated with vortex shedding and shear layer frequencies to predict the power law relationship with Reynolds number. The predicted results are in good agreement with experimental results. The findings will provide a greater insight into the overall phenomenon involved.展开更多
A comprehensive hot wire investigation of the flow around a circular cylinder is carried out in an 18" × 18" wind tunnel to look into the dominant frequencies at the stagnation, separation and separated shear l...A comprehensive hot wire investigation of the flow around a circular cylinder is carried out in an 18" × 18" wind tunnel to look into the dominant frequencies at the stagnation, separation and separated shear layers in the transition Reynolds number range. The majority of the experiments are carried out at Reynolds number of 4.5×104, with additional transition frequency tests at Reynolds numbers of 2.9×104, 3.3×104 and 9.7×104 respectively. The results are analysed in terms of power spectral density. While the frequency associated with stagnation is found to be essentially due to vortex shedding, frequency doubling of vortex shedding is also evident in the separated shear layers. Two peaks associated with transition frequencies are detected and their possible implications are presented.展开更多
The present study aims at the investigation of the effects of turbulence-chemistry interaction on combus- tion instabilities using a probability density function (PDF) method. The instantaneous quantities in the flo...The present study aims at the investigation of the effects of turbulence-chemistry interaction on combus- tion instabilities using a probability density function (PDF) method. The instantaneous quantities in the flow field were decomposed into the Favre-averaged variables and the stochastic fluctuations, which were calculated by unsteady Reynolds averaged Navier-Stokes (U-RANS) equations and the PDF model, respectively. A joint fluctuating velocity- frequency-composition PDF was used. The governing equa- tions are solved by a consistent hybrid finite volume/Monte- Carlo algorithm on triangular unstructured meshes. A non- reacting flow behind a triangular-shaped bluff body flame stabilizer in a rectilinear combustor was simulated by the present method. The results demonstrate the capability of the present method to capture the large-scale coherent struc- tures. The triple decomposition was performed, by divid- ing the coherent Favre-averaged velocity into time-averaged value and periodical coherent part, to analyze the coherent and incoherent contributions to Reynolds stresses. A sim- ple modification to the coefficients in the turbulent frequency model will help to improve the simulation results. Unsteady flow fields were depicted by streamlines and vorticity con- tours. Moreover, the association between turbulence produc- tion and vorticity saddle points is illustrated.展开更多
The suppression method of vortex shedding from a circular cylinder has been studied experimentally in the Reynolds number range from 300 to 1600. The test is performed in a water channel. The model cylinder is 1 cm in...The suppression method of vortex shedding from a circular cylinder has been studied experimentally in the Reynolds number range from 300 to 1600. The test is performed in a water channel. The model cylinder is 1 cm in diameter and 38 cm in length. A row of small rods of 0.18 cm in diameter and 1.5 cm in length are perpendicularly connected to the surface of the model cylinder and distributed along the meridian, The distance between the neighboring rods and the angle of attack of the rods can be changed so that the suppression effect on vortex shedding can be adjusted. The results show that vortex shedding can be suppressed effectively if the distance between the neighboring rods is smaller than 3 times and the cylinder diameter and the angle of attack is in the range of 30degreesless than or equal tobeta<90&DEG;.展开更多
The wake of a two-dimensional rotary oscillating plate of chord to thickness ratio B/H=3.0 was investigated by smoke wire visualization and hot-wire measurement in a wind tunnel at Re=1500 and oscillation amplitudes a...The wake of a two-dimensional rotary oscillating plate of chord to thickness ratio B/H=3.0 was investigated by smoke wire visualization and hot-wire measurement in a wind tunnel at Re=1500 and oscillation amplitudes and reduced frequencies in the ranges ofα=2.5°-50°and f_(e)H/V_(∞)=0.02-0.25,respectively,where fe is oscillation frequency,V_(∞) is on-coming wind speed.There are six modes of synchronized vortex shedding observed,including three new modes,i.e.,3P—three vortex pairs shed in one period,2S_(2)—two pairs of vortex twins shed in one period,and T—three vortices shed in one period of plate oscillation.The zone of every mode was found out and drawn on the amplitude-reduced frequency plane.Numerical simulation was carried out,and the simulated vorticity field,velocity profiles,and wake width were verified by experiment at typical oscillation conditions.At higher oscillation frequencies,the linear local unsteady term of inertial force dominates,and the simulated moment is an accurate linear function of rotary displacement and angular velocity.At lower frequencies,where 2T,3P,and 2S_(2) modes of vortex shedding occur,however,the nonlinear convective term of inertial force dominates,and the simulated fluctuating moment deviates from Scanlan’s expression.Statistical results of numerical data show that the root-mean-square(r.m.s)value of fluctuating moment is an exponential function of single parameter η=αf_(e)H/V_(∞).展开更多
Owing to the influence of the viscosity of the flow field,the strength of the shedding vortex decreases gradually in the process of backward propagation.Large-scale vortexes constantly break up,forming smaller vortexe...Owing to the influence of the viscosity of the flow field,the strength of the shedding vortex decreases gradually in the process of backward propagation.Large-scale vortexes constantly break up,forming smaller vortexes.In engineering,when numerical simulation of vortex evolution process is carried out,a large grid is needed to be arranged in the area of outflow field far from the boundary layer in order to ensure the calculation efficiency.As a result,small scale vortexes at the far end of the flow field cannot be captured by the sparse grid in this region,resulting in the dissipation or even disappearance of vortexes.In this paper,the effect of grid scale is quantified and compared with the viscous effect through theoretical derivation.The theoretical relationship between the mesh viscosity and the original viscosity of the flow field is established,and the viscosity term in the turbulence model is modified.This method proves to be able to effectively improve the intensity of small-scale shedding vortexes at the far end of the flow field under the condition of sparse grid.The error between the simulation results and the results obtained by using fine mesh is greatly reduced,the calculation time is shortened,and the high-precision and efficient simulation of the flow field is realized.展开更多
The effect of rigid bed proximity on flow parameters and hydrodynamic loads in offshore pipelines exposed to turbulent flow is investigated numerically. The Galerkin finite volume method is employed to solve the unste...The effect of rigid bed proximity on flow parameters and hydrodynamic loads in offshore pipelines exposed to turbulent flow is investigated numerically. The Galerkin finite volume method is employed to solve the unsteady incompressible 2D Navier–Stokes equations. The large eddy simulation turbulence model is solved using the artificial compressibility method and dual time-stepping approach. The proposed algorithm is developed for a wide range of turbulent flows with Reynolds numbers of 9500 to 1.5×10^4.Evaluation of the developed numerical model shows that the proposed technique is capable of properly predicting hydrodynamic forces and simulating the flow pattern. The obtained results show that the lift and drag coefficients are strongly affected by the gap ratio. The mean drag coefficient slightly increases as the gap ratio increases, although the mean lift coefficient rapidly decreases. The vortex shedding suppression happen at the gap ratio of less than 0.2.展开更多
Flows around rectangular cylinders with a series of width-to-height ratios are calculated by means of the Improved Finite Analytic Method (IFAM) and the formation, development and shedding of vortices from the cylinde...Flows around rectangular cylinders with a series of width-to-height ratios are calculated by means of the Improved Finite Analytic Method (IFAM) and the formation, development and shedding of vortices from the cylinders are simulated successfully. According to these results of time-dependent processes the physical phenomena in the flows are investigated in detail, and the discontinuity of Strouhal number is explained in the case of the width-to-height ratio equal to 3.0. The numer- ical solutions here show good agreement with the experimental results. In addition, based on several hundreds of the calculated flow patterns a moving picture is made by the computer image processing technology and recorded on a video tape, and then the vivid pictures of the physical process of vortex-shedding can be replayed later and analysed in detail.展开更多
Vortex-shedding flow induced by the vertical oscillation of a cylinder with bottom-attached disks of different diameter ratio Dd/Dc and thickness ratio td/Dc is studied by a 3D (three-dimensional) numerical model de...Vortex-shedding flow induced by the vertical oscillation of a cylinder with bottom-attached disks of different diameter ratio Dd/Dc and thickness ratio td/Dc is studied by a 3D (three-dimensional) numerical model developed in this paper, and compared with the results obtained through 2D (two-dimensional) numerical model. The high-order upwind scheme is applied to stabilize the computation, and convergence is accelerated by the multi-grid method. Qualitative and quantitative analyses of the differences between the 2D and 3D simulation results reveal the 3D effect on the flow field characteristics and hydrodynamic coefficients of the vertically oscillating cylinder with a bottom-attached disk. The 3D effect on the fluid field is mainly reflected in the significance of three vortex-shedding patterns: ωx has a greater effect on the flow fields around the sharp edges relative to the vortices generated in the 2D simulation. In the slice along the axial orientation, the vortex effect of ωy along the radial axis is smaller than that of ωx along the circumferential direction, indicating the radial effect on the velocity more pronounced than the circumferential effect around the sharp edges of the disk. The rotational interaction ωz of the fluid in the horizontal plane during the heave motion is insignificant. Based on the 2D and 3D simulation results, the turning points that separate the increasing regimes of the added mass coefficient and damping ratio are identified. The dependence of the turning point on the diameter ratio Dd/Dc and thickness ratio td/Dc are discussed in detail.展开更多
A study to investigate the effects of taper on vortex shedding coherence on High Mast Lighting Towers (HMLT) with models of eight-, twelve-, and sixteen-sided polygonal cross-section was performed in Purdue’s Boeing ...A study to investigate the effects of taper on vortex shedding coherence on High Mast Lighting Towers (HMLT) with models of eight-, twelve-, and sixteen-sided polygonal cross-section was performed in Purdue’s Boeing Low-Speed Wind Tunnel. Partial tower models were mounted on springs to recreate a flutter phenomenon seen on high mast lighting towers and data was taken using a stationary configuration within the wind tunnel. The model was later oscillated at specified frequencies and amplitudes and the resulting wake and surface pressures were recorded and compared to the stationary cases. The researchers aim to study the characteristics of a “lock-in” phenomenon, that is, a region of pole height where there is a vortex cell with a single shedding frequency, instead of different shedding frequencies for different diameters as Strouhal theory dictates. Results show the existence of vortex cell shedding for clamped models. Using a motor and a forcing cam to recreate the elastic movement of the HMLT in ambient conditions has yielded a specific range of diameters to determine the size of the locked in vortex cells. According to standard Department of Transportation manufacturing standards for tapered HMLT, the lock in distance for small excitations (0.254 cm) would be approximately 305 cm in tower height.展开更多
采用基于剪切应力输运模型的改进延迟分离涡模拟(improvement of delay detached-eddy simulation,IDDES)方法对两串列圆锥柱体高雷诺数流动干扰问题进行了数值计算,分析了不同间距比对流场结构和受扰圆锥柱体气动力特性的影响。结果表...采用基于剪切应力输运模型的改进延迟分离涡模拟(improvement of delay detached-eddy simulation,IDDES)方法对两串列圆锥柱体高雷诺数流动干扰问题进行了数值计算,分析了不同间距比对流场结构和受扰圆锥柱体气动力特性的影响。结果表明:间距比对两串列圆锥柱体的流场特性影响明显;当间距比为2时,下游圆锥柱体受到上游圆锥柱体的遮挡影响,其平均阻力系数较低;当间距比为4时,上游圆锥柱体的尾涡在下游圆锥柱体表面发生再附,下游(受扰)圆锥柱体的脉动阻力系数比小间距比时有所增加;随着间距比的继续增加,当间距比为8时,上、下游圆锥柱体形成各自的旋涡脱落,下游圆锥柱体的脉动升力系数明显高于间距比为2~4时的情况。展开更多
文摘Transition of flows past a pair of side-by-side circular cylinders are investigated by numerical simulations and the bifurcation analysis of the numerical results. Various flow patterns behind the cylinder-pair have been identified by the gap ratio (G) and Reynolds number (Re). This study focus on transition of in-phase and anti-phase vortex shedding synchronized forms. A nested Cartesian-grid formulation, in combination with an effective immersed boundary method and a two-step fractional-step procedure, has been adopted to simulate the flows. Numerical results reveal that the in-phase and anti-phase vortex shedding flows at Re = 100 can co-exist at 2.08 ≤G≤ 2.58. Hysteresis loop with increasing/decreasing G at constant Reynolds number Re = 100 is reported.
基金supported by the National Natural Science Foundation of China(No.11972190)。
文摘A robust Reynolds-Averaged Navier-Stokes(RANS)based solver is established to predict the complex unsteady aerodynamic characteristics of the Active Flap Control(AFC)rotor.The complex motion with multiple degrees of freedom of the Trailing Edge Flap(TEF)is analyzed by employing an inverse nested overset grid method.Simulation of non-rotational and rotational modes of blade motion are carried out to investigate the formation and development of TEF shedding vortex with high-frequency deflection of TEF.Moreover,the mechanism of TEF deflection interference with blade tip vortex and overall rotor aerodynamics is also explored.In nonrotational mode,two bundles of vortices form at the gap ends of TEF and the main blade and merge into a single TEF vortex.Dynamic deflection of the TEF significantly interferes with the blade tip vortex.The position of the blade tip vortex consistently changes,and its frequency is directly related to the frequency of TEF deflection.In rotational mode,the tip vortex forms a helical structure.The end vortices at the gap sides co-swirl and subsequently merge into the concentrated beam of tip vortices,causing fluctuations in the vorticity and axial position of the tip vortex under the rotor.This research concludes with the investigation on suppression of Blade Vortex Interaction(BVI),showing an increase in miss distance and reduction in the vorticity of tip vortex through TEF phase control at a particular control frequency.Through this mechanism,a designed TEF deflection law increases the miss distance by 34.7%and reduces vorticity by 11.9%at the target position,demonstrating the effectiveness of AFC in mitigating BVI.
基金the National Natural Science Foundation of China (Grant No.51076015)
文摘Vortex-acoustic coupling is one of the most important potential sources of combustion instability in solid rocket motors (SRMs). Based on the Von Karman Institute for Fluid Dynamics (VKI) experimental motor, the influence of the thermal inhibitor position and temperature on vortex-shedding-driven pressure oscillations is numerically studied via the large eddy simulation (LES) method. The simulation results demonstrate that vortex shedding is a periodic process and its accurate frequency can be numerically obtained. Acoustic modes could be easily excited by vortex shedding. The vortex shedding frequency and second acoustic frequency dominate the pressure oscillation characteristics in the chamber. Thermal inhibitor position and gas temperature have little effect on vortex shedding frequency, but have great impact on pressure oscillation amplitude. Pressure amplitude is much higher when the thermal inhibitor locates at the acoustic velocity anti-nodes. The farther the thermal inhibitor is to the nozzle head, the more vortex energy would be dissipated by the turbulence. Therefore, the vortex shedding amplitude at the second acoustic velocity antinode near 3/4L (L is chamber length) is larger than those of others. Besides, the natural acoustic frequencies increase with the gas temperature. As the vortex shedding frequency departs from the natural acoustic frequency, the vortex-acoustic feedback loop is decoupled. Consequently, both the vortex shedding and acoustic amplitudes decrease rapidly.
基金supported by the National Key Scientific Instrument and Equipment Development Program of China (Grant 2011YQ120048)
文摘The flow past a circular-section cylinder with a conic shroud perforated with four holes at the peak was simulated numerically at Re = 100, considering two factors, viz.the angle of attack and the diameter of the holes. The effects of the perforated conic shroud on the vortex shedding pattern in the near wake was mainly investigated, as well as the time history of the drag and lift forces. In the investigated parameter space, three flow regimes were generally identified, corresponding to weak, moderate, and strong disturbance effects.In regime I, the wake can mainly be described by alternately shedding Kármán or Kármán-like vortices. In regime II, the spanwise vortices are obviously disturbed along the span due to the appearance of additional vorticity components and their interactions with the spanwise vortices, but still shed in synchronization along the spanwise direction. In regime III,the typical Kármán vortices partially or totally disappear,and some new vortex shedding patterns appear, such as-type, obliquely shedding, and crossed spanwise vortices with opposite sign. Corresponding to these complex vortex shedding patterns in the near wake, the fluid forces no longer oscillate regularly at a single vortex shedding frequency, but rather with a lower modulation frequency and multiple amplitudes. An overview of these flow regimes is presented.
基金This project is supported by National Natural Science Foundation of China(No.59835160).
文摘Vortex street flowmeter has been used in steady flow measurement for about three decades The benefits of this type of flowmeter include high accuracy,good linearty,wide measuring range,and excellent reliability However,in unsteady flow measurement,the pressure disturbance as well as the noise from the system or surrounding can reduce the signal to noise ratio of the flowmeter seriously Aimed to use vortex street flowmeters in unsteady flow measurement,the characteristics of the vortex shedding induced hydrodynamic vibration around the prism bluff body in a vortex street flowmeter are investigated numerically and by expriments The results show that the hydrodynamic vibrations with 180° phase shift occur at the axisymmetric points of the channel around the bluff body The most intense vibration occurs at the points on the lateral faces close to the base of the prism The results provide therefore a useful reference for developing an anti interference vortex flowmeter using the different ial sensing technique.
文摘In this paper we have made a numerical study on the control of vortex shedding and drag reduction of a cylinder by attaching thin splitter plates. The wake structure of the cylinder of square cross-section with attached splitter plates is analyzed for a range of Reynolds number, based on the incident stream and height of the cylinder, in the laminar range. The Navier-Stokes equations governing the flow are solved by the control volume method over a staggered grid arrangement. We have used the semi-implicit method for pressure-linked equation (SIMPLE) algorithm for computation. Our results show that the presence of a splitter plate upstream of the cylinder reduces the drag, but it has a small impact on the vortex shedding frequency when the plate length is beyond 1.5 time the height of the cylinder. The presence of a downstream splitter plate dampens the vortex shedding frequency. The entrainment of fluid into the inner side of the separated shear layers is obstructed by the downstream splitter plate. Our results suggest that by attaching in-line splitter plates both upstream and downstream of the cylinder, the vortex shedding can be suppressed, as well as a reduction in drag be obtained. We made a parametric study to determine the optimal length of these splitter plates so as to achieve low drag and low vortex shedding frequency.
文摘An attempt has been made to explore whether the power relation can be obtained from theoretical considerations. The classical laminar and turbulent boundary layer concepts have been employed to determine appropriate values of the scaling lengths associated with vortex shedding and shear layer frequencies to predict the power law relationship with Reynolds number. The predicted results are in good agreement with experimental results. The findings will provide a greater insight into the overall phenomenon involved.
文摘A comprehensive hot wire investigation of the flow around a circular cylinder is carried out in an 18" × 18" wind tunnel to look into the dominant frequencies at the stagnation, separation and separated shear layers in the transition Reynolds number range. The majority of the experiments are carried out at Reynolds number of 4.5×104, with additional transition frequency tests at Reynolds numbers of 2.9×104, 3.3×104 and 9.7×104 respectively. The results are analysed in terms of power spectral density. While the frequency associated with stagnation is found to be essentially due to vortex shedding, frequency doubling of vortex shedding is also evident in the separated shear layers. Two peaks associated with transition frequencies are detected and their possible implications are presented.
基金supported by the National Natural Science Foundation of China (50936005)the Fundamental Research Fundsfor the Central Universities (WK2090130008,WK2090130011)
文摘The present study aims at the investigation of the effects of turbulence-chemistry interaction on combus- tion instabilities using a probability density function (PDF) method. The instantaneous quantities in the flow field were decomposed into the Favre-averaged variables and the stochastic fluctuations, which were calculated by unsteady Reynolds averaged Navier-Stokes (U-RANS) equations and the PDF model, respectively. A joint fluctuating velocity- frequency-composition PDF was used. The governing equa- tions are solved by a consistent hybrid finite volume/Monte- Carlo algorithm on triangular unstructured meshes. A non- reacting flow behind a triangular-shaped bluff body flame stabilizer in a rectilinear combustor was simulated by the present method. The results demonstrate the capability of the present method to capture the large-scale coherent struc- tures. The triple decomposition was performed, by divid- ing the coherent Favre-averaged velocity into time-averaged value and periodical coherent part, to analyze the coherent and incoherent contributions to Reynolds stresses. A sim- ple modification to the coefficients in the turbulent frequency model will help to improve the simulation results. Unsteady flow fields were depicted by streamlines and vorticity con- tours. Moreover, the association between turbulence produc- tion and vorticity saddle points is illustrated.
基金This study is financially supported by the National Natural Science Foundaion of China(Grant No.10172087)
文摘The suppression method of vortex shedding from a circular cylinder has been studied experimentally in the Reynolds number range from 300 to 1600. The test is performed in a water channel. The model cylinder is 1 cm in diameter and 38 cm in length. A row of small rods of 0.18 cm in diameter and 1.5 cm in length are perpendicularly connected to the surface of the model cylinder and distributed along the meridian, The distance between the neighboring rods and the angle of attack of the rods can be changed so that the suppression effect on vortex shedding can be adjusted. The results show that vortex shedding can be suppressed effectively if the distance between the neighboring rods is smaller than 3 times and the cylinder diameter and the angle of attack is in the range of 30degreesless than or equal tobeta<90&DEG;.
基金supported by the National Natural Science Foundation of China(Grant No.11572305).
文摘The wake of a two-dimensional rotary oscillating plate of chord to thickness ratio B/H=3.0 was investigated by smoke wire visualization and hot-wire measurement in a wind tunnel at Re=1500 and oscillation amplitudes and reduced frequencies in the ranges ofα=2.5°-50°and f_(e)H/V_(∞)=0.02-0.25,respectively,where fe is oscillation frequency,V_(∞) is on-coming wind speed.There are six modes of synchronized vortex shedding observed,including three new modes,i.e.,3P—three vortex pairs shed in one period,2S_(2)—two pairs of vortex twins shed in one period,and T—three vortices shed in one period of plate oscillation.The zone of every mode was found out and drawn on the amplitude-reduced frequency plane.Numerical simulation was carried out,and the simulated vorticity field,velocity profiles,and wake width were verified by experiment at typical oscillation conditions.At higher oscillation frequencies,the linear local unsteady term of inertial force dominates,and the simulated moment is an accurate linear function of rotary displacement and angular velocity.At lower frequencies,where 2T,3P,and 2S_(2) modes of vortex shedding occur,however,the nonlinear convective term of inertial force dominates,and the simulated fluctuating moment deviates from Scanlan’s expression.Statistical results of numerical data show that the root-mean-square(r.m.s)value of fluctuating moment is an exponential function of single parameter η=αf_(e)H/V_(∞).
基金Project supported by the National Key Project,China(Grant No.GJXM92579)the National Natural Science Foundation of China(Grant No.12072232)。
文摘Owing to the influence of the viscosity of the flow field,the strength of the shedding vortex decreases gradually in the process of backward propagation.Large-scale vortexes constantly break up,forming smaller vortexes.In engineering,when numerical simulation of vortex evolution process is carried out,a large grid is needed to be arranged in the area of outflow field far from the boundary layer in order to ensure the calculation efficiency.As a result,small scale vortexes at the far end of the flow field cannot be captured by the sparse grid in this region,resulting in the dissipation or even disappearance of vortexes.In this paper,the effect of grid scale is quantified and compared with the viscous effect through theoretical derivation.The theoretical relationship between the mesh viscosity and the original viscosity of the flow field is established,and the viscosity term in the turbulence model is modified.This method proves to be able to effectively improve the intensity of small-scale shedding vortexes at the far end of the flow field under the condition of sparse grid.The error between the simulation results and the results obtained by using fine mesh is greatly reduced,the calculation time is shortened,and the high-precision and efficient simulation of the flow field is realized.
基金Supported by the Technology Innovation Program(Grant number:10053121)funded by the Ministry of Trade,Industry&Energy(MI,Korea)by the Energy Efficiency&Resource of Korea Institute of Energy Technology Evaluation and Planning(KETEP)grant funded by the Ministry of Knowledge Economy of Korea(Grant number:2014301002-1870)
文摘The effect of rigid bed proximity on flow parameters and hydrodynamic loads in offshore pipelines exposed to turbulent flow is investigated numerically. The Galerkin finite volume method is employed to solve the unsteady incompressible 2D Navier–Stokes equations. The large eddy simulation turbulence model is solved using the artificial compressibility method and dual time-stepping approach. The proposed algorithm is developed for a wide range of turbulent flows with Reynolds numbers of 9500 to 1.5×10^4.Evaluation of the developed numerical model shows that the proposed technique is capable of properly predicting hydrodynamic forces and simulating the flow pattern. The obtained results show that the lift and drag coefficients are strongly affected by the gap ratio. The mean drag coefficient slightly increases as the gap ratio increases, although the mean lift coefficient rapidly decreases. The vortex shedding suppression happen at the gap ratio of less than 0.2.
基金The project supported by the National Natural Science Foundation of China
文摘Flows around rectangular cylinders with a series of width-to-height ratios are calculated by means of the Improved Finite Analytic Method (IFAM) and the formation, development and shedding of vortices from the cylinders are simulated successfully. According to these results of time-dependent processes the physical phenomena in the flows are investigated in detail, and the discontinuity of Strouhal number is explained in the case of the width-to-height ratio equal to 3.0. The numer- ical solutions here show good agreement with the experimental results. In addition, based on several hundreds of the calculated flow patterns a moving picture is made by the computer image processing technology and recorded on a video tape, and then the vivid pictures of the physical process of vortex-shedding can be replayed later and analysed in detail.
基金financially supported by the National Natural Science Foundation of China(Grant Nos.51239007 and 51179077)the Sino-UK Higher Education Research Partnership for Ph.D.Studies
文摘Vortex-shedding flow induced by the vertical oscillation of a cylinder with bottom-attached disks of different diameter ratio Dd/Dc and thickness ratio td/Dc is studied by a 3D (three-dimensional) numerical model developed in this paper, and compared with the results obtained through 2D (two-dimensional) numerical model. The high-order upwind scheme is applied to stabilize the computation, and convergence is accelerated by the multi-grid method. Qualitative and quantitative analyses of the differences between the 2D and 3D simulation results reveal the 3D effect on the flow field characteristics and hydrodynamic coefficients of the vertically oscillating cylinder with a bottom-attached disk. The 3D effect on the fluid field is mainly reflected in the significance of three vortex-shedding patterns: ωx has a greater effect on the flow fields around the sharp edges relative to the vortices generated in the 2D simulation. In the slice along the axial orientation, the vortex effect of ωy along the radial axis is smaller than that of ωx along the circumferential direction, indicating the radial effect on the velocity more pronounced than the circumferential effect around the sharp edges of the disk. The rotational interaction ωz of the fluid in the horizontal plane during the heave motion is insignificant. Based on the 2D and 3D simulation results, the turning points that separate the increasing regimes of the added mass coefficient and damping ratio are identified. The dependence of the turning point on the diameter ratio Dd/Dc and thickness ratio td/Dc are discussed in detail.
文摘A study to investigate the effects of taper on vortex shedding coherence on High Mast Lighting Towers (HMLT) with models of eight-, twelve-, and sixteen-sided polygonal cross-section was performed in Purdue’s Boeing Low-Speed Wind Tunnel. Partial tower models were mounted on springs to recreate a flutter phenomenon seen on high mast lighting towers and data was taken using a stationary configuration within the wind tunnel. The model was later oscillated at specified frequencies and amplitudes and the resulting wake and surface pressures were recorded and compared to the stationary cases. The researchers aim to study the characteristics of a “lock-in” phenomenon, that is, a region of pole height where there is a vortex cell with a single shedding frequency, instead of different shedding frequencies for different diameters as Strouhal theory dictates. Results show the existence of vortex cell shedding for clamped models. Using a motor and a forcing cam to recreate the elastic movement of the HMLT in ambient conditions has yielded a specific range of diameters to determine the size of the locked in vortex cells. According to standard Department of Transportation manufacturing standards for tapered HMLT, the lock in distance for small excitations (0.254 cm) would be approximately 305 cm in tower height.
文摘采用基于剪切应力输运模型的改进延迟分离涡模拟(improvement of delay detached-eddy simulation,IDDES)方法对两串列圆锥柱体高雷诺数流动干扰问题进行了数值计算,分析了不同间距比对流场结构和受扰圆锥柱体气动力特性的影响。结果表明:间距比对两串列圆锥柱体的流场特性影响明显;当间距比为2时,下游圆锥柱体受到上游圆锥柱体的遮挡影响,其平均阻力系数较低;当间距比为4时,上游圆锥柱体的尾涡在下游圆锥柱体表面发生再附,下游(受扰)圆锥柱体的脉动阻力系数比小间距比时有所增加;随着间距比的继续增加,当间距比为8时,上、下游圆锥柱体形成各自的旋涡脱落,下游圆锥柱体的脉动升力系数明显高于间距比为2~4时的情况。