共查询到20条相似文献，搜索用时 58 毫秒
1.
The apparent bed roughness, the roughness value experienced by a mean flow outside the wave-boundary layer, is deduced from the physical bed roughness and the wave–current interaction mechanism. Both the physical bed roughness and the wave–current interaction are described by a (combination of) model(s). Modelling of the apparent bed roughness leads to realistic results, however, the final results are rather sensitive to the particular choice of these models. Four bed form models and two wave–current interaction models were implemented in a 1-DV flow model to calculate near-bed velocities. A comparison between measured and predicted velocities shows that reasonable results can be obtained in this way. A constant bed roughness of 0.1 m, however, leads to even better results at this site during all conditions. This can be explained by the reversed influence of the form roughness and the wave–current interaction on the apparent bed roughness value for varying wave conditions. 相似文献
2.
This paper discusses numerical results from three-dimensional large eddy simulations of an oscillating cylinder under prescribed movements in uniform flow. Six cases,namely pure in-line,pure cross-flow and two groups of 'Figure of Eight' oscillation patterns are under investigation at Reynolds number Re=24000. The 'Figure of Eight' pattern in each group is with identical shape but opposite orbital directions. The numerical results on hydrodynamic forces,higher order force components,and vortex shedding mode... 相似文献
3.
The flow patterns in the near wake of a cylinder (either circular or square in shape, D=25 mm) placed in the proximity of a fully developed turbulent boundary layer (thickness δ=0.4D) are investigated experimentally using particle image velocimetry (PIV). The effects of changing the gap height (S) between the cylinder bottom and the wall surface, over the gap ratio range S/D=0.1–1.0, have been investigated. The results show that both the ensemble-averaged and instantaneous flow fields are strongly dependent on S/D. The flow patterns for the two types of cylinders share many similarities with respect to the change in S/D, such as the reduced recirculation length and increased velocity fluctuation in the near wake with increasing S/D, as well as the trend of suppression of vortex shedding at small S/D and onset of vortex shedding at large S/D. However, developments of the shear layers, in terms of wake width, flow curvature, etc., are considerably different for these two types of cylinders. In general, the wake development and momentum exchange for the square cylinder are slower those for the circular cylinder at the same gap ratio. Correspondingly, it is shown that the periodic vortex shedding is delayed and weakened in the case of square cylinder, as compared to that of the circular cylinder at the same S/D. 相似文献
4.
Andrs E. Tejada-Mart&#x; nez Chester E. Grosch Ann E. Gargett Jeff A. Polton Jerome A. Smith J.A. MacKinnon 《Ocean Modelling》2009,30(2-3):115-142
A three-dimensional numerical model for large-eddy simulation (LES) of oceanic turbulent processes is described. The numerical formulation comprises a spectral discretization in the horizontal directions and a high-order compact finite-difference discretization in the vertical direction. Time-stepping is accomplished via a second-order accurate fractional-step scheme. LES subgrid-scale (SGS) closure is given by a traditional Smagorinsky eddy-viscosity parametrization for which the model coefficient is derived following similarity theory in the near-surface region. Alternatively, LES closure is given by the dynamic Smagorinsky parametrization for which the model coefficient is computed dynamically as a function of the flow. Validation studies are presented demonstrating the temporal and spatial accuracy of the formulation for laminar flows with analytical solutions. Further validation studies are described involving direct numerical simulation (DNS) and LES of turbulent channel flow and LES of decaying isotropic turbulence. Sample flow problems include surface Ekman layers and wind-driven shallow water flows both with and without Langmuir circulation (LC), generated by wave effects parameterized via the well-known Craik–Leibovich (C–L) vortex force. In the case of the surface Ekman layers, the inner layer (where viscous effects are important) is not resolved and instead is parameterized with the Smagorinsky models previously described. The validity of the dynamic Smagorinsky model (DSM) for parameterizing the surface inner layer is assessed and a modification to the surface stress boundary condition based on log-layer behavior is introduced improving the performance of the DSM. Furthermore, in Ekman layers with wave effects, the implicit LES grid filter leads to LC subgrid-scales requiring ad hoc modeling via an explicit spatial filtering of the C–L force in place of a suitable SGS parameterization. 相似文献
5.
The flow about a circular cylinder placed centrally inside a channel is studied numerically with an unstructured collocated grid finite volume method based on the primitive variable formulation. The distance between the channel walls is allowed to vary to change the blockage ratio. Simulations are carried out over a range of Reynolds numbers that are consistent with the two-dimensional assumption. The study confirms that transition to vortex shedding regime is delayed when the channel walls are close to the cylinder because of the interaction between the vortices from the channel wall and cylinder wake. In the unsteady vortex shedding regime, the wake pattern is opposite to the classic Karman street in respect of the positions of the shed vortices. The cylinder drag coefficient and Strouhal number are considerably increased at smaller gaps while the root-mean-squared lift coefficient is significantly decreased. Several important flow parameters are correlated with the input parameters, namely Reynolds number and blockage ratio. 相似文献
6.
In this study, the Navier-Stokes equations and the pressure Poisson equation for two-dimensional time-dependent viscous flows are solved with a finite difference method in a curvilinear coordinate system. With this numerical procedure, the vortex shedding flow past a circular cylinder near a wall is investigated. The flow is calculated for a broad range of gap ratios for different Reynolds numbers ranging from 80 to 1000. Based on the numerical solutions, the vortex shedding is observed using various methods, and the mechanism for the vortex shedding suppression at small gap ratios is analyzed. The critical gap ratio at which the vortex shedding is suppressed is identified at different Reynolds numbers. 相似文献
7.
Numerical study on the characteristics of flow field and wave propagation near submerged breakwater on slope 总被引：2，自引：1，他引：1
In this study,characteristics of flow field and wave propagation near submerged breakwater on a sloping bed are investigated with numerical model. The governing equations of the vertical twodimensional model are Reynolds Averaged Navier Stokes equations. The Reynolds stress terms are closed by a nonlinear k ε turbulence transportation model. The free surface is traced through the PILC-VOF method. The proposed numerical model is verified with experimental results. The numerical result shows that the wave profile may become more asymmetrical when wave propa-gates over breakwater. When wave crest propagates over breakwater,the anticlockwise vortex may generate. On the contrary,when wave hollow propagates over breakwater,the clockwise vortex may generate. Meanwhile,the influenced zone of vortex created by wave crest is larger than that created by wave hollow. All the maximum values of the turbulent kinetic energy,turbulent dissi-pation and eddy viscosity occur on the top of breakwater. Both the turbulent dissipation and eddy viscosity increase as the turbulent kinetic energy increases. Wave energy may rapidly decrease near the breakwater because turbulent dissipation increases and energy in lower harmonics is transferred into higher harmonics. 相似文献
8.
基于计算流体力学(CFD)开源代码OpenFOAM开展了不同雷诺数(Re=100、1500和3900)和倾斜角度(-60°≤α≤60°)工况下倾斜圆柱绕流流场的三维数值模拟,研究倾斜圆柱绕流的三维瞬时及时均尾流流场、流线拓扑、升阻力系数与旋涡脱落频率随雷诺数及倾斜角度变化的规律,探讨在顺流向及逆流向情况下独立性原则对倾斜圆柱绕流的适用性。研究结果表明:随着圆柱倾角的增大,倾斜圆柱尾流产生较为明显的轴向流,尾流旋涡脱落受到明显抑制,细碎旋涡逐渐消失,尾流宽度随之减小;随着雷诺数的增大,圆柱尾流涡管发生明显的变形,展向掺混使得大量细碎旋涡产生,呈现出明显的三维特性。在不同雷诺数下,阻力系数均值、升力系数均方根及无量纲涡脱频率在一定倾角范围内符合独立性原则。 相似文献
9.
The unsteady flow past a circular cylinder which starts translating and transversely oscillating from rest in a viscous fluid is investigated at a Reynolds numbers of R=103 and at a Strouhal number of π/4 and for the maximum oscillatory to translational velocity ratios between 0.1 and 1.0. This study is based on numerical solutions of the two-dimensional unsteady Navier–Stokes equations. The object of the study is to examine the effect of increase of velocity ratio on the near-wake structure as well as the hydrodynamic forces acting on the cylinder. For all velocity rates a periodic structure of vortex evaluation and shedding develops which is repeated exactly as time advances. Vortex dynamics close behind the body are affected by changing acceleration of the cylinder and a changeover from one mode to a different mode of vortex formation is observed with increase in velocity ratio. A comparison of the present results with the impulsively started translating case has been included to illustrate the effect of velocity ratio on drag at small values of velocity ratio. 相似文献
10.
涡旋核心线是海洋中尺度涡旋结构的重要组成要素,涡旋核心线提取和可视化对于切入中尺度涡三维结构研究、开展海洋物质能量垂直运输分析具有重要意义。本文基于客观参考框架和准则,提出了使用客观化的流场参数得到区域,并提取其山谷线作为涡旋核心线的方法,实现了对海洋三维结构中尺度涡旋核心线的提取和可视化。首先,引入了最优局部参考系,使速度、速度梯度等测度转换为在运动的参考系下保持不变的客观量,提升了在海洋科学实践中的可靠性和实用性。其次,针对含有垂向速度的海洋三维流场数据,计算其空间雅可比矩阵,展示了涡旋核心区域的三维结构,实现了海洋涡旋研究从二维到三维的提升。最后,分别在多个半径大小播撒流线种子点,分析不同旋转方向的涡旋,对已提取的涡旋核心线实验结果进行验证,证明了客观海洋三维涡旋核心线提取方法的有效性及可行性。 相似文献
11.
The slightly compressible flow formulation is applied to the free-surface, three-dimensional turbulent flow around a Wigley hull. Two turbulence models (large eddy simulation and Baldwin–Lomax) are used and compared. The simulation conditions are the ones for which experimental and numerical results exist. The computational grid is built using an algebraic grid generator with the model fixed in space. The codes use the interface-capturing technique for computing the free-surface displacements and the Beam and Warming scheme for marching in time the numerical model. The results compare well with the experimental data available. 相似文献
12.
Based on the second-order solutions obtained for the three-dimensional weakly nonlinear random waves propagating over a steady uniform current in finite water depth, the joint statistical distribution of the velocity and acceleration of the fluid particle in the current direction is derived using the characteristic function expansion method. From the joint distribution and the Morison equation, the theoretical distributions of drag forces, inertia forces and total random forces caused by waves propagating over a steady uniform current are determined. The distribution of inertia forces is Gaussian as that derived using the linear wave model, whereas the distributions of drag forces and total random forces deviate slightly from those derived utilizing the linear wave model. The distributions presented can be determined by the wave number spectrum of ocean waves, current speed and the second order wave–wave and wave–current interactions. As an illustrative example, for fully developed deep ocean waves, the parameters appeared in the distributions near still water level are calculated for various wind speeds and current speeds by using Donelan–Pierson–Banner spectrum and the effects of the current and the nonlinearity of ocean waves on the distribution are studied. 相似文献
13.
The problem of unsteady, laminar flow past a circular cylinder which starts translating and oscillating impulsively from rest in a viscous fluid is numerically investigated at a Reynolds number of R = 103. The flow is incompressible and two-dimensional, and the cylinder oscillations are harmonic. The transverse oscillations are only allowed when the maximum oscillatory-to-translational velocity ratio is 0.5. The investigation is based on an implicit finite difference scheme for integrating the unsteady Navier-Stokes equations together with the mass-conservation equation in their vorticity stream function formulation. A non-inertial coordinate transformation is used so that the grid mesh remains fixed relative to the accelerating cylinder. Present calculations are performed within the range of sufficiently large oscillation amplitude to induce separation. The time variation of the in-line and transverse force coefficients are presented. The study also focuses on the laminar asymmetric flow structure in the near-wake region. In this flow regime, it is found that there is alternate shedding of vortices from either side of the cylinder over an oscillation cycle (as predicted experimentally); this is the classical mode of vortex shedding leading to formation of the Kármán street. 相似文献
14.
In this study the basic characteristics of the dynamic response and vortex shedding from an elastically mounted circular cylinder in laminar flow is numerically investigated. The Reynolds number ranges from 80 to 160, a regime that is fully laminar. The governing equations of fluid flow are cast in terms of vorticity. The two-dimensional vorticity transport equation is solved using a vortex method. Effects of important parameters on the system response and vortex shedding are investigated; these include: mass ratio, damping ratio, Reynolds number and reduced velocity. The numerical results show that a decrease in either the mass ratio or damping ratio of the system can lead to an increase in both the oscillation amplitude and the reduced velocity range over which lock-in occurs. The results also suggest that the mass-damping parameter may characterize the system response adequately, although the effect of changing mass ratio appears to be a little more profound compared to damping ratio. Vorticity contour plots suggest that the vortex shedding occurs in the 2S mode, although a wake structure similar to the C(2S) mode appears at distances 15–20 diameters downstream in the lock-in region. The simulation results are in good agreement with previously published data. 相似文献
15.
为研究四柱体布置情况下倒角半径变化对柱体绕流水动力特性的影响,使用Fluent软件,采用大涡模拟方法研究了在雷诺数Re=3 900下6种不同倒角半径的柱体在方形四柱体布置时的三维流场。在模型分析验证有效后,分析了柱体后方瞬时流场、水动力参数、时均流场的变化情况。分析结果表明:随着倒角半径的增大,上游柱体的平均阻力系数逐渐减小,下游柱体的平均阻力系数除了在R~+=0.1处增幅很大以外,其余均随倒角半径变大而平稳变大;各柱体的升力系数均方根变化趋势基本相同;R~+=0.1、0.5时,上下游两柱体的升力系数曲线相位相反,而在R~+=0.2、0.3和0.4时,上下游两柱体的升力系数曲线相位相同。 相似文献
16.
Alan M. Davies 《Coastal Engineering》1987,11(5-6)
A transformation method is presented by which current profiles (of tidal or wind-induced origin) can be extracted at any horizontal position and moment in time from a vertically integrated, two-dimensional, hydrodynamic numerical model. An arbitrary vertical variation of eddy viscosity can be included in the method, which can incorporate a no-slip bottom boundary condition. The technique assumes that the sea is homogeneous.The method is used to improve the representation of bottom stress within the two-dimensional model, whereby the bottom stress is no longer related simply to the depth-mean current as in the “conventional” two-dimensional, vertically integrated model.Idealized calculations for a range of eddy viscosity profiles, show that elevations, current profiles, and time series of current extracted from this “enhanced” two-dimensional numerical model are in good agreement with currents obtained from a full three-dimensional model. 相似文献
17.
A coupled wave–tide–surge model has been developed in this study in order to investigate the effect of the interactions among tides, storm surges, and wind waves. The coupled model is based on the synchronous dynamic coupling of a third-generation wave model, WAM cycle 4, and the two-dimensional tide–surge model. The surface stress, which is generated by interactions between wind and wave, is calculated by using the WAM model directly based on an analytical approximation of the results using the quasi-linear theory of wave generation. The changes in bottom friction are created by the interactions between waves and currents and calculated by using simplified bottom boundary layer model. In consequence, the combined wave–current-induced bottom velocity and effective bottom drag coefficient were increased in the shallow waters during the strong storm conditions. 相似文献
18.
《Coastal Engineering》2006,53(5-6):441-462
The structure of large-scale turbulence under a broken solitary wave on a 1 in 50 plane slope was studied. Three-component velocity measurements were taken at different heights above a smooth bed in the middle surf zone using an acoustic Doppler velocimeter. The measured data showed that turbulent velocity components were well correlated in the middle part of the water column. The velocity correlations could be produced by an oblique vortex similar to the obliquely descending eddy observed previously by other investigators. The vertical distributions of the relative values of the components of the Reynolds stress tensor showed that the structure of turbulence evolved continuously between the free surface and the bottom. The evolution was related to transition from two-dimensional to three-dimensional flow structures and the effect of the solid bottom on flow structures. Time histories of measured turbulent kinetic energy and turbulence stresses showed episodic turbulent events near the free surface but more sporadic turbulence in the lower layer. Large or intense turbulent events were found to have short duration and time lag relative to the wave crest point. These events also maintained good correlations between the turbulence velocity components close to the bottom.Instantaneous turbulent velocity fields were measured near the bottom at the same cross-shore location by using a stereoscopic particle image velocimetry system. These measurements showed that the near-bed flow field was characterized by large-scale, coherent flow structures that were the sources of most of the turbulent kinetic energy and turbulence stresses. The types of organized flow structures observed included vortices and downbursts of turbulence descending directly from above, lateral spreading of turbulent fluid along the bed, and formation of vortices in shear layers between fluid streams. A common feature of the organized flow structures near the bed was the large turbulence velocities in the longitudinal and transverse directions, which reflected the influence of a solid bottom on the breaking-wave-generated turbulence arriving at the bed. 相似文献
19.
Direct numerical simulation was conducted to investigate the flow past a slotted cylinder at low Reynolds number (Re) of 100. The slotting of cylinder affects the boundary layer separation, vortex formation position, recirculation region length and wake width, which are determined by the type of slit. The streamwise slit (SS1), T-shaped slit (SS3) and Y-shaped slit (SS4) act as passive jets, while the transverse slit (SS2) achieves an alternate self-organized boundary layer suction and blowing. The flow rate in slits fluctuates over time due to the alternate vortex shedding and fluctuating pressure distribution around the cylinder surface. One fluctuation cycle of flow rate is caused by a pair of vortices shedding for SS2, SS3 and SS4, while it is created by each vortex shedding for SS1. The wall shear stress and flow impact on the slit wall partly contribute to the hydrodynamic forces acting on the slotted cylinder. Taking into account the internal wall of slit, the transverse slit plays the best role in suppressing the fluid forces with drag reduction of 1.7% and lift reduction of 17%. 相似文献
20.
两层粘性流体中圆柱体受迫振荡数值模拟 总被引：1，自引：1，他引：1
研究两层粘性流体中无限长水平圆柱体的受迫振荡问题。在湍流模式下,采用VOF方法追踪两层流体的内界面,基于动网格技术模拟圆柱体的运动边界,对均匀流中横向振荡圆柱体的绕流场进行了数值模拟。计算受迫振荡圆柱体的升力系数、阻力系数随时间的演化曲线和圆柱体的尾涡分布,以及圆柱体的受迫振荡激发两层流体内界面的扰动,并与均匀流体的情况进行了比较分析。研究表明,流体的两层分层效应对受迫振荡圆柱体的升阻力系数和尾涡分布特性都有显著影响,在水下输油气管道涡激振动特性的工程评估中,应考虑流体的密度分层效应。 相似文献