首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 125 毫秒
1.
In this paper a series of numerical simulations are performed to investigate the vortex shedding mechanism for a solitary wave propagating over a submerged breakwater by use of Reynolds averaged Navier-Stokes (RANS) model combined with a k-ε model. Flows of different Reynolds numbers up to Re=1.4×105 corresponding to varying incident wave heights are considered in which the characteristic fluid velocity is represented by the maximum horizontal velocity above the submerged breakwater. For the verification of the accuracy of the numerical model, the incident waves and the velocity field in the vicinity of the breakwater are compared with experimental data. The result shows that the model is capable of describing vortex shedding for a solitary wave propagating over a rectangular submerged breakwater. Key features of vortex generation, evolution and dissipation are investigated. It is found that the vortex shedding and their evolution due to separated boundary layer over the breakwater are strongly related to the Reynolds number. A considerable number of vortices and complicated vortex pattern are observed as the Reynolds number increases.  相似文献   

2.
- A one-layer time-invariant eddy viscosity model is specified to develop a mathematical model for describing the essential features of the turbulent wave boundary layer over a rough bed. The functional form of the eddy viscosity is evaluated based on computational results from a two-equation turbulence model in which the eddy viscosity varies with time and space. The present eddy viscosity model simplifies much of the mathematical complexity in many existing models. Predictions from the present model have been compared with a wide range of experimental data. It is found that the eddy viscosity model adopted in the present study is physically reasonable.  相似文献   

3.
Effects of Stokes production on summer ocean shelf dynamics   总被引:1,自引:0,他引:1  
A two-dimensional numerical model,which is configured on the basis of Princeton ocean model(POM),is used to study the effect of Stokes production(SP) of the turbulent kinetic energy on a density profile and Ekman transport in an idealized shelf region in summer.The energy input from SP is parameterized and included into the Mellor-Yamada turbulence closure submodel.Results reveal that the intensity of wind-driven upwelling fronts near the sea surface is weakened by the SP-associated turbulent kinetic energy input.The vertical eddy viscosity coefficient in the surface boundary layer is enhanced greatly owing to the impact of SP,which decreases the alongshore velocity and changes the distribution of upwelling.In addition,the SP-induced mixing easily suppresses the strong stratification and significantly increases the depth of the upper mixed layer(ML) under strong winds.  相似文献   

4.
In this paper a series of numerical simulations are performed to investigate the vortex shedding mechanism for a solitary wave propagating over a submerged breakwater by use of Reynolds averaged Navier-Stokes (RANS) model combined with a k-ε model. Flows of different Reynolds numbers up to Re=1.4×105 corresponding to varying incident wave heights are considered in which the characteristic fluid velocity is represented by the maximum horizontal velocity above the submerged breakwater. For the verification of...  相似文献   

5.
The effects of various types of waves on vertical plane turbulent jets are studied numerically in the paper. A three-dimensional numerical model in σ -coordinate is developed to study these problems by use of large eddy simulation method. Turbulence is modeled by a dynamic coherent eddy model. Numerical results including the distribution of velocity, the decay law of the mean velocity along axis, the turbulent Reynolds stresses and the volume flux per unit width without wave, in the first-order Stokes waves, in the second-order Stokes waves, in the fifth-order Stokes waves, in the solitary waves and in random waves are compared and analyzed. A focus on coherent structures, probability density functions and correlation functions of jets is also investigated. The numerical results are of great theoretical importance for understanding jet turbulent behaviors in different types of waves.  相似文献   

6.
A three-dimensional numerical model is established to simulate the turbulent oscillatory boundary layer over a fixed and rough bed composed by randomly arrayed solid spheres based on the lattice Boltzmann method and the large eddy simulation model.The equivalent roughness height,the location of the theoretical bed and the time variation of the friction velocity are investigated using the log-fit method.The time series of turbulent intensity and Reynolds stress are also investigated.The equivalent roughness height of cases with Reynolds numbers of 1×10~4–6×10~4 is approximately 2.81 d(grain size).The time variation of the friction velocity in an oscillatory cycle exhibits sinusoidal-like behavior.The friction factor depends on the relative roughness in the rough turbulent regime,and the pattern of solid particles arrayed as the rough bed in the numerical simulations has no obvious effect on the friction factor.  相似文献   

7.
Large Eddy Simulation for Wave Breaking in the Surf Zone   总被引:1,自引:0,他引:1  
In this paper, (he large eddy simulation method is used combined with the marker and cell method to study the wave propagation or shoaling and breaking process. As wave propagates into shallow water, the shoaling leads lo the increase of wave height, and then at a certain position, the wave will be breaking. The breaking wave is a powerful agent for generating turbulence, which plays an important role in most of the fluid dynamic processes throughout the surf zone, such as transformation of wave energy, generation of near-shore current and diffusion of materials. So a proper numerical model for describing the turbulence effect is needed. In this paper, a revised Smagorinsky subgrid-scale mode! is used to describe the turbulence effect. The present study reveals that the coefficient of the Smagorinsky model for wave propagation or breaking simulation may be taken as a varying function of the water depth and distance away from the wave breaking point. The large eddy simulation model presented in this pape  相似文献   

8.
The effects of various types of waves on vertical plane turbulent jets are studied numerically in the paper. A three-dimensional numerical model in σ-coordinate is developed to study these problems by use of large eddy simulation method. Turbulence is modeled by a dynamic coherent eddy model. Numerical results including the distribution of velocity,the decay law of the mean velocity along axis,the turbulent Reynolds stresses and the volume flux per unit width without wave,in the first-order Stokes waves,in ...  相似文献   

9.
刘臻 《中国海洋工程》2011,25(1):169-178
Oscillating Water Column (OWC) wave energy converting system is one of the most widely used facilities all over the world.The air chamber is utilized to convert the wave energy into the pneumatic energy.The numerical wave tank based on the two-phase VOF model is established in the present study to investigate the operating performance of OWC air chamber.The RANS equations,standard k-ε turbulence model and dynamic mesh technology are employed in the numerical model.The effects of incident wave conditions and shape parameters on the wave energy converting efficiency are studied and the capability of the present numerical wave tank on the corresponding engineering application is validated.  相似文献   

10.
This paper aims to propose an improved numerical model for wave breaking in the nearshore region based on the fully nonlinear form of Boussinesq equations. The model uses the k equation turbulence scheme to determine the eddy viscosity in the Boussinesq equations. To calculate the turbulence production term in the equation, a new formula is derived based on the concept of surface roller. By use of this formula, the turbulence production in the one-equation turbulence scheme is directly related to the difference between the water portide velocity and the wave celerity. The model is verified by Hansen and Svendsen‘s experimental data (1979) in terms of wave height and setup and sctdown. The comparison between the model and experimental results of wave height and setup and setdown shows satisfactory agreement. The modeled turbulence energy decreases as waves attenuate in the surf zone. The modeled production term peaks at the breaking point and decreases as waves propagate shoreward. It is also suggested that both convection and diffusion play their important roles in the transport of turbulence energy immediately after wave breaking. When waves approach to the shoreline, the production and dissipation of turbulence energy are almost balanced. By use of the slot technique for the simulation of the movable shoreline boundary, wave ranup in the swash zone is well simulated by the present model.  相似文献   

11.
《Coastal Engineering》2002,46(1):25-50
A mathematical model based on the Volume-Averaged/Reynolds Averaged Navier-Stokes (VARANS) equations is developed to describe surface wave motions in the vicinity of a coastal structure, which could be either a rigid solid structure or a permeable structure or a combination of both. In the VARANS equations, the volume-averaged Reynolds stress is modeled by adopting the nonlinear eddy viscosity assumption. The model equations for the volume-averaged turbulent kinetic energy and its dissipation rate are derived by taking the volume-average of the standard kϵ equations. Because of the volume-averaging process, the effects of the small-scale turbulence in porous media are introduced. The performance of the model is checked by comparing numerical solutions with the experimental data related to a composite breakwater reported by Sakakiyama and Liu [Coast. Eng. 121 (2001) 117].  相似文献   

12.
Numerical simulations are performed to investigate the influence of variable front slopes on flow evolution and waveform inversion of a depression ISW (internal solitary wave) over an idealized shelf with variable front slopes. A finite volume based on Cartesian grid method is adopted to solve the Reynolds averaged Navier-Stokes equations using a k-ε model for the turbulent closure. Numerical results exhibit the variations of several pertinent properties of the flow field, in the case with or without waveform inversion on the horizontal plateau of an obstacle. The clockwise vortex is stronger than the counterclockwise one, almost throughout the wave-obstacle interaction. Analysis of the turbulent energy budget reveals that the turbulent production term in the governing equations dominates the wave evolution during a wave-obstacle interaction; otherwise the buoyancy production term and the dissipation term due to viscosity within turbulent eddies play a major role in energy dissipation. In addition, the front slope affects mainly the process and reflection of the wave evolution but has less influence than other physical parameters. Moreover, total wave energy of the leading crest is smaller than that of the leading trough even in the cases with waveform inversion on the plateau.  相似文献   

13.
14.
为模拟潜堤上破碎波浪传播时产生能量的耗散这一特性,在改进的具有四阶色散的Boussinesq水波方程中中入二阶紊动粘性项,建立了考虑波浪破碎的水波数学模型.在非交错网格下建立了有限差分数值模型,并利用三阶Adams-Bash forth格式预报、四阶Adams-Mouton格式校正对数值模型进行求解.通过数值试验,模拟...  相似文献   

15.
We study the interactions between a non-breaking solitary wave and a submerged permeable breakwater experimentally and numerically. The particle image velocimetry (PIV) technique is employed to measure instantaneous free surface displacements and velocity fields in the vicinity of a porous dike. The porous medium, consisting of uniform glass spheres, is mounted on the seafloor. Due to the limited size of each field of view (FOV) for high spatial resolution purposes, four FOVs are set in order to form a continuous flow field around the structure. Quantitative mean properties are obtained by ensemble averaging 30 repeated instantaneous measurements. The Reynolds decomposition method is then adopted to separate the velocity fluctuations for each trial to estimate the turbulent kinetic energy. In addition, a highly accurate two-dimensional model with the volume of fluid interface tracking technique is used to simulate an idealized volume-averaged porous medium. The model is based on the Volume-Averaged Reynolds Averaged Navier–Stokes equations coupled with the non-linear kε turbulence closure solver. Comparisons are performed between measurements and numerical results for the time histories of the free surface elevation recorded by wave gauges and the spatial distributions of free surface displacement with the corresponding velocity and turbulent kinetic energy around the permeable object imaged by the PIV system. Fairly good agreements are obtained. It is found that the measured and modeled turbulent intensities on the weather side are much larger than those on the lee side of the object, and that the magnitude of the turbulent intensity increases with increasing wave height of a solitary wave at a constant water depth. The verified numerical model is then used to estimate the energy reflection, transmission and dissipation using the energy integral method by varying the aspect ratio and the grain size of the permeable obstacle.  相似文献   

16.
A smoothed particle hydrodynamic (SPH) model is developed to simulate wave interaction with porous structures. The mean flow outside the porous structures is obtained by solving Reynolds Averaged Navier–Stokes (RANS) equations and the turbulence field is calculated by a large eddy simulation (LES) model. The porous flow is described by the spatially averaged Navier–Stokes type equations with the resistance effect of the porous media being represented by an empirical frictional source term. The interface boundaries between the porous flow and the outside flow are modeled by means of specifying a transition zone along the interface. The model is validated against other available numerical results and experimental data for wave damping over porous seabed with different levels of permeability. The validated model is then employed to investigate wave breaking over a submerged porous breakwater and good agreements between the SPH model results and the experimental data are obtained in terms of free surface displacement. In addition the predicted velocity, vorticity and pressure fields near the porous breakwater and in the breaking wave zone are also analyzed.  相似文献   

17.
基于非静压数值计算模型,本文系统研究了聚焦波作用下透水潜堤的消波特性,通过设置合理的计算工况,详细分析了波高、堤顶水深、谱峰周期、孔隙率以及堤顶宽度5种因素对透水潜堤消波特性的影响。与此同时,本文将透水潜堤的计算结果同不透水潜堤的计算结果进行了对比分析。计算结果表明:透水潜堤对聚焦波的消减作用要强于不透水潜堤,从而说明,透水潜堤能更有效地降低畸形波对海岸基础设施的影响;波高和堤顶水深是影响潜堤消波特性的重要因素,随入射波高增加、堤顶水深减小,透水潜堤对波浪的消减作用逐渐增强。透水潜堤对长周期波浪的消波效果较差。在本文考虑的孔隙率范围内,孔隙率越大,透水潜堤消波效果越好;当孔隙率为0.4,堤顶宽度为0.612 5 m时,透水潜堤可消减54%的入射波能,比不透水潜堤对入射波能的消减增加36.1%。本文研究结果可为进一步认识透水潜堤的消波特性和海岸防护工程设计提供相应的参考。  相似文献   

18.
Pengzhi Lin  C. W. Li   《Ocean Engineering》2003,30(7):855-876
A three-dimensional numerical model is developed in this study to investigate the problem of wave–current–body interaction. The model solves the spatially averaged Navier–Stokes equations. Turbulence effects are modeled by a subgrid-scale (SGS) model using the concept of large eddy simulation (LES). The model is employed to study the wave–current interaction with a square cylinder that is mounted on the bottom and vertically pierces the free surface. The force analysis demonstrates that the presence of waves can reduce both the strength and frequency of vortex shedding induced by a uniform current due to the nonlinear wave–current interaction. The free surface elevation, strain rates of the mean flow, and eddy viscosity are found to closely correlate with the mechanism of vortex shedding. It is also shown that when the vortex shedding is neglected in the calculation such as by the potential flow approach, one may significantly underestimate the magnitude of in-line force. The energy spectral analysis reveals that there exist initiating, growing, and decaying regions for shedding vortices around the cylinder. In the vortex initiating region, both coherent and turbulent structures are nearly two-dimensional that become three-dimensional in the vortex growing region. The kinetic energy of both coherent and turbulent motions is dissipated in the vortex decaying region, within which the mean flow gradually returns back to two-dimensional.  相似文献   

19.
《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.  相似文献   

20.
Breaking wave induced nearsurface turbulence has important consequences for many physical and biochemical processes including water column and nutrients mixing,heat and gases exchange across air-sea interface.The energy loss from wave breaking and the bubble plume penetration depth are estimated.As a consequence,the vertical distribution of the turbulent kinetic energy(TKE),the TKE dissipation rate and the eddy viscosity induced by wave breaking are also provided.It is indicated that model results are found to be consistent with the observational evidence that most TKE generated by wave breaking is lost within a depth of a few meters near the sea surface.High turbulence level with intensities of eddy viscosity induced by breaking is nearly four orders larger than υwl(=κuwz),the value predicted for the wall layer scaling close to the surface,where uw is the friction velocity in water,κ with 0.4 is the von Kármán constant,and z is the water depth,and the strength of the eddy viscosity depends both on wind speed and sea state,and decays rapidly through the depth.This leads to the conclusion that the breaking wave induced vertical mixing is mainly limited to the near surface layer,well above the classical values expected from the similarity theory.Deeper down,however,the effects of wave breaking on the vertical mixing become less important.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号