共查询到20条相似文献，搜索用时 62 毫秒
1.
When fluid flow passes a cylinder, the drag crisis phenomenon occurs between the sub-critical and the super-critical Reynolds numbers. The focus of the present studies was on the numerical prediction of the drag crisis based on CFD methods. In this work, block structured meshes with refined grids near the cylinder surface and in the downstream were employed. Both 2D and 3D simulations were performed using various turbulence models, including the SST k − ω model, the k − ϵ model, the SST with LCTM, the DES model, and the LES model. In the convergence studies, the effects of the grid size, the time step, the first grid size and the aspect ratio (for 3D simulations) on the solutions were examined. The errors due to spatial and time discretizations were quantified according to a V&V procedure. Validation studies were carried out for various Reynolds numbers between Re = 6.31 × 104 and 7.57 × 105. The averaged drag force, the RMS of lift force and the Strouhal number were compared with experimental data. The studies indicated that standard 2D and 3D RANS methods were inadequate to capture the drag crisis phenomenon. The LES method however has the potential to address the problem. 相似文献
2.
Large eddy simulations of the flow around a circular cylinder at high Reynolds numbers are reported. Five Reynolds numbers were chosen, such that the drag crisis was captured. A total of 18 cases were computed to investigate the effect of gridding strategy, turbulence modelling, numerical schemes and domain width on the results. It was found that unstructured grids provide better resolution of key flow features, when a ‘reasonable’ grid size is to be maintained.When using coarse grids for large eddy simulation, the effect of turbulence models and numerical schemes becomes more pronounced. The dynamic mixed Smagorinsky model was found to be superior to the Smagorinsky model, since the model coefficient is allowed to dynamically adjust based on the local flow and grid size. A blended upwind-central convection scheme was also found to provide the best accuracy, since a fully central scheme exhibits artificial wiggles, due to dispersion errors, which pollute the solution.Mean drag, fluctuating lift Strouhal number and base pressure are compared to experiments and empirical estimates for Reynolds numbers ranging from 6.31 × 104 to 5.06 × 105. In terms of the drag coefficient, the drag crisis is well captured by the present simulations, although the other integral quantities (rms lift and Strouhal number) show larger discrepancies. For the lowest Reynolds number, the drag is seen to be more sensitive to the domain width than the spanwise grid spacing, while at the higher Reynolds numbers the grid resolution plays a more important role, due to the larger extent of the turbulent boundary layer. 相似文献
3.
The aim of this paper is to evaluate the accuracy, stability and efficiency of the overset grid approach coupled with the RANS (Reynolds Averaged Navier-Stokes) model via the benchmark computations of flows around a stationary smooth circular cylinder. Two dimensional numerical results are presented within a wide range of Reynolds numbers (6.31 × 104 ∼ 7.57 × 105) including the critical flow regime. All the simulations are carried out using the RANS solver pimpleFoam provided by OpenFOAM, an open source CFD (Computational Fluid Dynamics) toolkit. Firstly, a grid convergence study is performed. The results of the time-averaged drag and lift force coefficients, root-mean square value of lift force coefficient and Strouhal number (St number) are then compared with the experimental data. The velocity, vorticity fields and pressure distribution are also given. One main conclusion is that the numerical solutions in regard to a fixed cylinderare not deteriorated due to the implementation of the overset grid. Furthermore, it can be an appealing approach to facilitate simulations of Vortex Induced Vibrations (VIV), which involves grid deformation. The present study is a good start to implement the overset grid to solve VIV problems in the future. 相似文献
4.
The response of an oscillating circular cylinder at the wake of an upstream fixed circular cylinder was classified by different researchers as galloping, wake induced galloping or wake induced vibration. Furthermore it is already known that a sharp edge square cylinder would undergo galloping if it is subjected to uniform flow. In this study the influence of the wake of a fixed circular cylinder on the response of a downstream square cylinder at different spacing ratios (S/D = 4, 8, 11) is experimentally investigated. The subject appears not to have received previous attention. The lateral displacements, lift forces and the pressure data from gauges mounted in the wake of the oscillating cylinder are recorded and analyzed. The single degree of freedom vibrating system has a low mass-damping parameter and the Reynolds number ranges from 7.7 × 102 to 3.7 × 104.In contrast to that for two circular cylinders in tandem arrangement, the freely mounted downstream square cylinder displays a VIV type of response at all spacing ratios tested. There is no sign of galloping or wake induced galloping with the square cylinder. With increase at the spacing ratio the cross-flow oscillations decrease. It is shown that the vortices arriving from the upstream fixed circular cylinder play a major role on the shedding mechanism behind the downstream square cylinder and cause the square cylinder to shed vortices with frequencies above Strouhal frequency of the fixed square cylinder (St = 0.13). The VIV type of oscillations in the downstream square cylinder is most probably caused by the vortices newly generated behind the square cylinder. 相似文献
5.
A practical method to account for the influence of sinkage and trim on the drag of a freely floating (free to sink and trim) common monohull ship at a Froude number F ≤ 0.45 is considered. The sinkage and the trim are estimated via two alternative simple methods, considered previously. The drag is also estimated in a simple way, based on the classical Froude decomposition into viscous and wave components. Specifically, well-known semiempirical expressions for the friction drag, the viscous pressure drag and the drag due to hull roughness are used, and the wave drag is evaluated via a practical linear potential flow method. This simple approach can be used for ship models as well as full-scale ships with smooth or rough hull surfaces, and is well suited for early ship design and optimization. The method considered here to determine the sinkage and the trim, and their influence on the drag, yields theoretical predictions of the drag of the Wigley, S60 and DTMB5415 hulls that are much closer to experimental measurements than the corresponding predictions for the hull surfaces of the ships in equilibrium position at rest. These numerical results suggest that sinkage and trim effects, significant at Froude numbers 0.25 < F, on the drag of a typical freely floating monohull ship can be realistically accounted for in a practical manner that only requires simple potential flow computations without iterative computations for a sequence of hull positions. 相似文献
6.
Experiments employing a low-mass-damping cylinder have been conducted to determine the vortex-induced vibration (VIV) response of four suppressors of the flexible-shroud family. The VIV suppressors were inspired in the concept of the Ventilated Trousers (VT), a flexible shroud composed of a flexible net fitted with three-dimensional bobbins. Reynolds number varied between 5 × 103 and 25 × 103, while reduced velocity varied from 2 to 26. The VIV dynamic response showed that the VT suppressed the peak amplitude of vibration down to 40% of that of a bare cylinder. Other flexible shrouds also achieved suppression, but not as efficiently. Drag was reduced during the VIV synchronization range, but remained above the value for a bare static cylinder thereafter. Spectral analysis of displacement and lift revealed that, depending on the geometry and distribution of the bobbins, the flexible shroud can develop an unstable behavior, capturing energy from the wake and sustaining vibrations for higher reduced velocities. PIV measurements of the wake revealed that the entrainment flow through the mesh is necessary to extend the vortex-formation length of the wake; this mechanism only occurs for the VT mesh. 相似文献
7.
Flow past a circular cylinder with multiple small control rods is studied by numerical simulation for ReD ranging from 1161.3 to 6387.1. The Reynolds-Averaged-Navier–Stokes (RANS) equations and shear stress transport (SST) k − ω turbulence model are used to calculate the vortex field, while a fourth-order Runge–Kutta method is employed for evaluating the structure dynamics of the cylinder group. Comparisons with experimental results demonstrate the validation of this method. This study is concerned with the vortex induced vibration (VIV) suppression efficacy of small control rods placed around a main cylinder. The effects of control rod number, diameter ratio, spacing ratio and Reynolds number on the hydrodynamics and vibration responses of the main cylinder are investigated. The reduced percents of in-line and cross-flow amplitudes and the increased percents of the whole cross-sectional area of cylinders and the drag coefficient are used to give a comprehensive evaluation. Results of simulation indicate that placing small rods with appropriate number at appropriate locations can achieve good suppression effectiveness at a wide range of Reynolds number. The numerical result for the case with nine control rods, diameter ratio of 0.15 and spacing ratio of 0.6 shows the best suppression effect among the cases investigated in this study. 相似文献
8.
Experimental studies were carried out to investigate the response features of an inclined flexible bare cylinder as well as a straked cylinder in a towing tank, with the main purpose of further improving the understanding of the effect of yaw angle on vortex-induced vibration (VIV) suppression. Four yaw angles (a = 0°, 15°, 30°, 45°), which is defined as the angle between the cylinder axis and the plane orthogonal to the oncoming fluid flow, were tested. The cylinder model was towed along the tank to generate a uniform fluid flow. The towing velocity was in the range of 0.05–1.0 m/s with an interval of 0.05 m/s. The corresponding Reynolds number ranged from 800 to 16000. The strakes selected for the experiments had a pitch of 17.5D and a height of 0.25D, which is generally considered as the most effective configuration for VIV suppression of a flexible cylinder in water. The experimental results indicate that VIV suppression effectiveness of the inclined flexible straked cylinder is closely related to the yaw angle. The displacement amplitudes are significantly suppressed in both cross-flow (CF) and in-line (IL) directions at a = 0°. However, with increasing yaw angle, the suppression efficiencies of the CF and IL displacement amplitudes gradually decrease. In addition, the CF dominant frequencies of the straked cylinder obviously deviate from those of the bare cylinder at a = 0° and 15°. This deviation is substantially alleviated with increasing yaw angle. The IL dominant frequencies show less dependency on the yaw angle. Similar trends are also observed on the dominant modes of vibration and the mean drag coefficients. 相似文献
9.
为研究四柱体布置情况下倒角半径变化对柱体绕流水动力特性的影响,使用Fluent软件,采用大涡模拟方法研究了在雷诺数Re=3 900下6种不同倒角半径的柱体在方形四柱体布置时的三维流场。在模型分析验证有效后,分析了柱体后方瞬时流场、水动力参数、时均流场的变化情况。分析结果表明:随着倒角半径的增大,上游柱体的平均阻力系数逐渐减小,下游柱体的平均阻力系数除了在R~+=0.1处增幅很大以外,其余均随倒角半径变大而平稳变大;各柱体的升力系数均方根变化趋势基本相同;R~+=0.1、0.5时,上下游两柱体的升力系数曲线相位相反,而在R~+=0.2、0.3和0.4时,上下游两柱体的升力系数曲线相位相同。 相似文献
10.
Vegetation canopies control mean and turbulent flow structure as well as surface wave processes in coastal regions. A non-hydrostatic RANS model based on NHWAVE (Ma et al., 2012) is developed to study turbulent mixing, surface wave attenuation and nearshore circulation induced by vegetation. A nonlinear k − ϵ model accounting for vegetation-induced turbulence production is implemented to study turbulent flow within the vegetation field. The model is calibrated and validated using experimental data from vegetated open channel flow, as well as nonbreaking and breaking random wave propagation in vegetation fields. It is found that the drag-related coefficients in the k − ϵ model Cfk and Cfϵ can greatly affect turbulent flow structure, but seldom change the wave attenuation rate. The bulk drag coefficient CD is the major parameter controlling surface wave damping by vegetation canopies. Using the empirical formula of Mendez and Losada (2004), the present model provides accurate predictions of vegetation-induced wave energy dissipation. Wave propagation through a finite patch of vegetation in the surf zone is investigated as well. It is found that the presence of a finite patch of vegetation may generate strong pressure-driven nearshore currents, with an onshore mean flow in the unvegetated zone and an offshore return flow in the vegetated zone. 相似文献
11.
The purpose of this research work is to study the effect of specific surface s, the fluid–solid contact surface per volume unit, on the wave energy dissipation by porous structures consisting in dense arrays of emergent vertical cylinders. Experiments have been carried out in a 10 m long wave flume. Three cylinder diameters D are considered in order to study the effects of the specific surface while keeping the porosity constant. In a first series, the length of the porous zone is kept constant for the three cylinder diameters tested. The measurements, which include various wave steepness conditions, demonstrate the role of specific surface s on both wave attenuation and interference processes. The larger the specific surface is, the stronger the wave damping is. Damping is found to be almost proportional to 1/D when laminar, turbulent and inertial effects are of same order. Results are compared to numerical calculations based on either a constant rate of wave damping within the porous medium per unit wavelength or a quadratic damping developed using a force expression based on the work of [26]. This latter model, calibrated with drag and inertia coefficients, shows a good agreement with measurements. In a second series, both porous length and water depth are kept proportional to the cylinder diameter for the three diameters. Scale effects are then discussed and underline the importance of the flow regime within the porous medium. 相似文献
12.
Vengatesan Venugopal Kamlesh S. Varyani Nigel D.P. Barltrop 《Ocean Engineering》2006,33(11-12):1669-1704
The results of wave force measurements carried out on a section of horizontally submerged rectangular cylinders, which are used as pontoons in many offshore structures, are reported in this paper. Two rectangular cylinders with aspect (depth–breadth) ratios equal to and and a square section (aspect ratio=1.0) cylinder are chosen for this study. Experiments are carried out in a wave tank at a water depth of 2.2 m at low Keulegan–Carpenter (KC) numbers to measure the horizontal and vertical wave forces acting on a 100 mm section, located at mid-length of the cylinders. For each cylinder, tests are carried out for two relative depths of submergence of 2.68 and 4.68. Measured wave forces in regular and irregular waves are then used to derive drag (CD) and inertia coefficients (CM). The analysis show that at very low KC numbers the inertia coefficients for all cylinders approached the potential flow values for both horizontal and vertical forces. The drag coefficients at low KC numbers exhibited large values and they decreased sharply with increase in KC number. For the square cylinder, where relatively a large KC number is obtained compared to other cylinders, inertia coefficients reached minimum values in the range of KC of about 3–4 and increased thereafter. In this range, CM values are about 50% or so, smaller than the same at KC close to zero. The results of the experiments reveal that aspect ratio has large influence on hydrodynamic coefficients. 相似文献
13.
《Coastal Engineering》2006,53(9):781-792
This paper presents results of a series of detailed measurements of geometric and migrating characteristics of ripples superimposed upon sandwaves under the action of combined waves and currents. Velocity measurements within the fluid, surface wave characteristics and 3D mapping of the bottom were recorded with an Acoustic Doppler Velocimeter (ADV), an acoustic water level sensor and a 32 composite element array of sub-aquatic acoustic sensors, respectively. Bottom records were statistically analyzed to obtain height, length and migration rates of ripples. Experiments examined ripple heights and wavelengths for the mobility factors (as defined in Eq. (4)) and the Reynolds wave number within the ranges 10 < ψ < 88 and 16 × 103 < Rew < 5 × 105, respectively. Measured values were compared with laboratory and field data together with semi-empirical and analytical formulae from the literature. Good correlation was obtained when plotting measured ripple length and length in dimensionless form as a function of the Reynolds wave number Rew. Under a given hydraulic condition, it was observed that ripples with different geometric characteristics may coexist at different locations over the sandwave. Ripple steepness is presented as a function of the Shields parameter although characterized with rather large scatter. Finally, average ripple migration speed is presented as a function of the Shields parameter and the mobility number. 相似文献
14.
15.
We describe experiments with multi-directional focused waves interacted with a vertical circular cylinder in a 3D wave basin. The focus of this study is on the run-up of multi-directional focused waves, wave forces, and wave pressures on the cylinder. Part I, the study on wave run-up, has already been presented by Li et al. (2012). In this paper, the analysis of the wave force on the vertical cylinder is presented.In this experiment, a cylinder with 0.25 m in diameter was adopted and different wave parameters, such as focused wave amplitude, peak frequency, frequency bandwidth and directional spreading index, are considered. The model scale kpa (kp is the wave number corresponding to peak frequency, a is the radium of the cylinder) varies from 0.32 to 0.65. The maximum forces of multi-directional focused wave on cylinder were measured and investigated. The results showed that the wave parameters have a significant influence on the wave force, and that the spatial profile of the surface of multi-directional focused wave can also affect its force on the cylinder, which is different from two-dimensional wave. In addition, the ‘secondary loading cycle’ phenomenon was also observed and discussed. In our experiments, the ‘secondary loading cycles’ occur when kA > 0.36 for all cases. While in some referred small scale experiments, the secondary load cycles are observed even for kA = 0.2, when the waves are longer enough. To larger model scale, the pronounced secondary load cycle occurs with larger wave steepness waves. 相似文献
16.
A renewable energy harvester using the piezoelectric effect is developed for the ocean tidal and wind flow. The harvester is made of connected driving blades to an octo-generator, which has a rotator with n blades and a stator attached by eight mass-spring-piston-cylinder-piezoelectricity devices. The resonance and force magnification are utilized to increase the power output of the harvester. A corresponding mathematical model is developed to calculate the root mean square of the generated electric power. The simulation results indicate that the generated power is largely enhanced when the near-resonant condition is established. The power increases with increases in the magnetic flux density, the large-to-small diameter ratio of the cylinder, the size of magnetic bar face, and decreases in the gap between two magnetic faces and the size of the piezoelectric bar face. A generated power of 5 kW is realized by the harvester working under an ocean tidal speed, V = 1.75 m/s, and its geometric and material properties of driving length L = 7.5 m, spring constant kv = 65000 N/m, gap between the two magnets s = 0.0015 m, large to small diameter ratio of the cylinder z = 6, and magnetic flux density Br = 1.45 T. 相似文献
17.
The paper examines the dependency between total sediment transport, q, and grain size, D (i.e. q ∝ Dp) under dam break generated swash flows. Experiments were performed in a dam break flume over a sloping mobile sand bed with median grain sizes ranging from 0.22 mm to 2.65 mm. The total sediment transport was measured by truncating the flume bed and collecting the sediment transported over the edge. The experiments were designed to exclude pre-generated turbulence and pre-suspended sediment so as to focus solely on the swash flow. The magnitude and nature of the grain size dependency (i.e. p value) were inferred for different flow parameters; the initial dam depth, do, the integrated depth averaged velocity cubed, ∫ u3dt, and against the predicted transport potential, qp, using the Meyer-Peter Muller (MPM) transport model and variations of that model. The data show that negative dependencies (p < 0) are obtained for do and qp, whilst positive dependencies (p > 0) are obtained for ∫ u3dt. This indicates that a given do and qp transport less sediment as grain size increases, whereas transport increases with grain size for a given ∫ u3dt. The p value is found to be narrowly ranged, 0.5 ≤ p ≤ − 0.5. On average, the incorporation of a pressure gradient term via the piezometric head into the MPM formulation reduces qp by 4% (fine sand) to 18% (coarse sand). The measured total transport for fine and coarse sands is best predicted using MPM and MPM + dp*/dx respectively. However, the inferred optimum transport coefficient in the MPM formulation is about 30, much higher than the standard coefficient in a steady flow and this is not due to the presence of the pre-suspended sediment. The optimum transport coefficient indicates some sensitivity to grain size, suggesting that some transport processes remain unaccounted for in the model. 相似文献
18.
A horizontal, circular cylinder fitted with one bilge keel is forced to rotate harmonically around its axis. The bilge keel load and hull pressure distribution are investigated. A fully submerged condition (infinite fluid), and three partly-submerged conditions are considered. A two-dimensional numerical study is performed, and the results are validated against recently published experimental data by van’t Veer et al. [30]. In addition, comparisons for mass and drag coefficients are also made with experimental data for plate in infinite fluid (Keulegan and Carpenter [8]), and wall-mounted plate (Sarpkaya and O’Keefe [9]) in oscillatory flow.A Navier–Stokes solver based on the Finite Volume Method is adopted for solving laminar flow of incompressible water. The free-surface condition is linearized by neglecting the nonlinear free-surface terms and the influence of viscous stresses in the free surface zone, while the body-boundary condition is exact. This simplified modeling of the problem required the mesh to be fine only around the bilge keels, leading to a total number of cells around N ≃ 1 ×104, which reduced computational cost significantly.The influence of draft and amplitude of oscillations on the bilge keel force and hull pressure distribution are considered. The bilge keel force is presented in terms of non-dimensional drag and mass coefficients including higher harmonic components. The numerical results are also compared with the industry standard empirical method for calculation of roll damping proposed by Ikeda et al. [4]. In general, a good agreement between the results of the present numerical method and the experimental data is obtained and the differences with those predicted by the empirical method are addressed. 相似文献
19.
Laminar flow past a circular cylinder with 3 small control rods is investigated by numerical simulation. This study is concerned with the suppression efficacy of vortex induced vibration by small control rods located around a main cylinder. The effects of the attack angle and rod-to-cylinder gap ratio on the hydrodynamics and vibration responses of the main cylinder are investigated. The attack angle of α = 45° is performed as the critical angle for VIV suppression of 3 control rods. The 3 control rods have no effect on VIV suppression when the attack angle is less than the critical angle. The 3 control rods have an excellent VIV suppression efficacy when the attack angle is larger than the critical angle. The transverse vibration frequency of the cylinder with 3 control rods is less than that for an isolated cylinder for all the configurations. The numerical results for the configurations of α = 45° & 60°, G/D = 0.6–1.2 show excellent suppression efficient among the cases investigated in this study. The best suppression efficient is found at α = 45°, G/D = 0.9 for 3 control rods. 2 rods in behind of the main cylinder perform more efficient than that of 1 rod in front for VIV suppression as the gap ratio of G/D less than 1.0. 相似文献
20.
《Estuarine, Coastal and Shelf Science》2006,66(1-2):1-5
Salt intrusion in estuaries is important for ecological reasons as well as water extraction purposes. The distance salt intrudes upstream depends on a number of factors, including river discharge, tidal and wind mixing and gravitational circulation. In this paper, an analytical solution is presented for the salt intrusion in a well mixed, funnel-shaped estuary whose cross sectional area decreases exponentially (with decay coefficient β) with distance, x, inland, and in which longitudinal mixing is constant along the length of the estuary. The solution predicts that a graph of the logarithm of salinity against exp (βx) should be a straight line, with slope proportional to the mixing coefficient Kx. The solution is tested against observations from 15 surveys over a four-year period in the Incomati estuary. Good straight line fits, as predicted, are observed on all surveys, with a mean R2 = 0.97. The average value of Kx for all surveys is 38 m2 s−1. The solution is used to make predictions about the minimum river flow required to prevent salt intruding to an extent where it causes a detrimental effect on water extraction. The minimum recommended river flow required to prevent this is 35 m3 s−1. In recent years, flow has fallen below this level for several months each year. 相似文献