共查询到20条相似文献，搜索用时 265 毫秒
1.
A three-dimensional (3D) large-eddy-simulation model with macroscopic model equations of porous flow is proposed to investigate solitary waves interacting with permeable breakwaters. The major objective of this paper is twofold. First, we seek to evaluate the present model through the comparison with available simulated and measured data in the literature. The second aim, given the 3D nature of flow past a permeable breakwater, the variations of permeable breakwater modeled on both macroscopic and microscopic scales are examined. First validation is carried out with experiments on solitary wave propagation in a 3D wave basin and then runup on a vertical permeable breakwater with a gap in the lateral direction. A satisfactory agreement on the free surface elevation time series is obtained between model and measured results. Second, we replicate the experiments on a solitary wave interaction with a submerged permeable breakwater in a two-dimensional narrow wave flume. The porous medium is composed of spheres with a uniform size and arranged in a non-staggered regular pattern such that the porous medium can thus be modeled on macroscopic and microscopic scales. The numerical calculations indicate that the results obtained with macroscopic and microscopic modeling both fit the measurements fairly well in terms of the free surface elevations and velocity fields. Specifically, the microscopic modeling better simulates detailed phenomena such as flow injection from the porous medium and the initial stage of the formation of the main vortex in the leeward face of the obstacle. After the solitary wave completely propagates over the permeable object, the discrepancies between macroscopic and microscopic model results are insignificant. More accurate 3D results are used to determine the trajectories of fluid particles around the porous object to help understand the possible sediment movements in suspensions. 相似文献
2.
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. 相似文献
3.
A Numerical Investigation of the Reduction of Solitary Wave Runup by A Row of Vertical Slotted Piles
To improve the current understanding of the reduction of tsunami-like solitary wave runup by the pile breakwater on a sloping beach, we developed a 3D numerical wave tank based on the CFD tool OpenFOAM in this study. The Navier Stokes equations were applied to solve the two-phase incompressible flow, combined with an LES model to solve the turbulence and a VOF method to capture the free surface. The adopted model was firstly validated with existing empirical formulas for solitary wave runup on the slope without the pile structure. It is then validated using our new laboratory observations of the free surface elevation, the velocity and the pressure around a row of vertical slotted piles subjected to solitary waves, as well as the wave runup on the slope behind the piles. Subsequently, a set of numerical simulations were implemented to analyze the wave reflection, the wave transmission, and the shoreline runup with various offshore wave heights, offshore water depths, adjacent pile spaces and beach slopes. Finally, an improved empirical equation accounting for the maximum wave runup on the slope was proposed by taking the presence of the pile breakwater into consideration. 相似文献
4.
内孤立波具有振幅尺度大、能量集中的特点,其引起流场和密度场的迅速变化可能对海洋工程结构物以及水下潜体造成严重威胁。因此研究不同造波条件下生成的内孤立波运动的流场特征具有重要的学术意义和实际应用价值。采用直接数值模拟方法和给定的初始密度场密度跃迁函数,对重力塌陷激发内孤立波的运动过程进行研究,探讨了不同造波条件下,激发产生的内孤立波波型、涡度、振幅和水平速度等流场特征。结果表明:(1)直接模拟数值方法能够模拟内孤立波传播过程中的密度界面波型反转现象;(2)从定性和定量的角度,证实了不稳定内孤立波传播过程中存在能量的向后传递;(3)对于相同的台阶深度(水闸两侧初始密度界面的高度差),初始涡流保持相同,但是随着上下层水深比的减小,其强度下降显著;(4)台阶深度对初始涡流的垂直结构的影响要大于上下层水深比,且台阶深度对内孤立波的振幅、水平速度的影响显著。 相似文献
5.
Based on a parallel SPH-LES model, a three dimensional numerical wave basin is developed to study wave interaction with coastal structures. The OpenMP programming technology combined with an existing MPI program contained in the parallel version of SPHysics code has been implemented to enable the simulation of hundred millions of particles running on a computer cluster. As part of the numerical basin development work an active absorbing wave maker and a sponge layer are introduced. The dynamic boundary conditions are also corrected to reduce the spurious effects. Wave generation and propagation in the numerical wave basin is first tested and confirmed with analytical results. Then, the model is applied to simulate wave interactions with a vertical breakwater and a vertical cylinder in order to further assess the capability of the numerical wave basin. The predicted free surface elevation near the vertical breakwater is compared with the experimental data while the horizontal forces and overturning moments acting on the vertical cylinder are verified with the analytical results. In all these cases the model results show excellent agreement with the validation data. 相似文献
6.
Comparative study on the accuracy of solitary wave generations in an ISPH-based numerical wave flume
In the present study, six solitary wave generations by different mathematical approximations are investigated using a piston type wave maker at dimensionless amplitudes ranging from 0.1 to 0.6 and two water depths. Incompressible smoothed particle hydrodynamics is used to simulate solitary wave propagation along the fixed depth channel. The present numerical results are compared with analytical results and experimental data in terms of free surface displacements, fluid particle velocity, phase speed, paddle motion, etc. The present mesh-free numerical results of wave profile variations over time proved that “Rayleigh” has the lowest relative wave height variation. However, its solitary wave has notable phase lead, while “Third order” and “Ninth order” have the least wave lags. Furthermore, the record of present numerical free surface elevation at different distances and the loss of amplitude of the main pulse showed that regarding both of them, “Ninth order” has supremacy over five others. Considering the numerical velocity components of generated solitary wave, “Third order” and “Ninth order” trace analytical results more accurately than other four ones, whereas “Rayleigh” is the most accurate one in predicting the maximum runup. Finally, the paddle motion, its velocity, and displacement, as well as phase speed and outskirts decay coefficient are also compared and discussed intensely. 相似文献
7.
孤立波与带窄缝双箱相互作用模拟研究 总被引：1，自引：1，他引：0
针对孤立波与带窄缝双箱的作用问题,应用时域高阶边界元方法建立了二维数值水槽。其中,自由水面满足完全非线性运动学和动力学边界条件,对瞬时自由表面流体质点采用混合欧拉-拉格朗日法追踪,采用四阶龙格库塔法对下一时刻的自由水面的速度势和波面升高进行更新。采用加速度势法求解物体湿表面的瞬时波浪力。采用推板方法生成孤立波。通过模拟孤立波在直墙上的爬高以及施加在直墙上的波浪力,并与已发表的实验和数值结果对比,验证本数值模型的准确性。通过数值模拟计算研究了窄缝宽度、方箱尺寸对波浪在箱体迎浪侧爬高,窄缝内波面升高,箱体背浪侧透射波高及箱体受波浪荷载的影响。同时研究了有一定时间间隔的双孤立波与带窄缝双箱系统作用问题。 相似文献
8.
In this study, solitary waves passing over a submerged breakwater are investigated both experimentally and numerically. A total of 9 experimental conditions are carried out, including different incident wave heights and water depths. Numerical simulations are performed using a high-order finite-difference model solving Navier-Stokes (N-S) equations. The predicted water wave elevation, velocity and pressure show good agreement with experimental data, verifying the accuracy and capacity of the numerical model. Furthermore, parametric studies are conducted by numerical modelling to examine the effects of the geometrical features of submerged dike on hydrodynamic characteristics around the breakwater.
相似文献9.
堤前远破波运动与冲淤形态 总被引：1，自引：0，他引：1
利用数值波浪水槽,对远破波作用下堤前波浪运动进行了数值模拟。分析了堤前速度场的变化,指出堤前波浪运动在一个波周期内可分为向堤波和离堤波。波浪在距堤约L/2处破碎后向堤冲击并沿堤面爬升,称为向堤波;波浪从堤面最高处回落与向堤波相遇并破碎,称为离堤波。堤前远破波的前半个周期的向堤波速度场大于后半个周期的离堤波速度场,且堤前波浪运动在一个波周期内的前半个周期和后半个周期完全不同。床面某点(L/4,3L/4)的速度在一个波浪周期内随时间呈向堤和离堤变化,向堤速度最大值大于离堤速度的最大值。利用堤前速度场的研究成果,进一步分析了远破波作用下堤前海床冲淤形成过程,完善了冲淤机理。 相似文献
10.
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. 相似文献
11.
Vegetation damping effects on propagating water waves have been investigated by many researchers. This paper investigates the effects of damping due to vegetation on solitary water wave run-up via numerical simulation. The numerical model is based on an implementation of Morison's formulation for vegetation induced inertia and drag stresses in the nonlinear shallow water equations. The numerical model is solved via a finite volume method on a Cartesian cut cell mesh. The accuracy of the numerical scheme and the effects of the vegetation terms in the present model are validated by comparison with experiment results. The model is then applied to simulate a solitary wave propagating on a plane slope with vegetation. The sensitivity of solitary wave run-up to plant height, diameter and stem density is investigated by comparison of the numerical results for different patterns of vegetation. The numerical results show that vegetation can effectively reduce solitary wave propagation velocity and that solitary wave run-up is decreased with increase of plant height in water and also diameter and stem density. 相似文献
12.
A numerical assessment study of tsunami attack on the rubble mound breakwater of Haydarpasa Port, located at the southern entrance of the Istanbul Bosphorus Strait in the Sea of Marmara, Turkey, is carried out in this study using a Volume-Averaged Reynolds-Averaged Navier-Stokes solver, IHFOAM, developed in OpenFOAM® environment. The numerical model is calibrated with and validated against the data from solitary wave and tsunami overflow experiments representing tsunami attack. Furthermore, attack of a potential tsunami near Haydarpasa Port is simulated to investigate effects of a more realistic tsunami that cannot be generated in a wave flume with the present state of the art technology. Discussions on practical engineering applications of this type of numerical modeling studies are given focusing on pressure distributions around the crown-wall of the rubble mound breakwater, and the forces acting on the single stone located behind the crown-wall at the rear side of the breakwater. Numerical modeling of stability/failure mechanism of the overall cross-section is studied throughout the paper.The present study shows that hydrodynamics along the wave flume and over the breakwater can be simulated properly for both solitary wave and tsunami overflow experiments. Stability of the overall cross-section can only be simulated qualitatively for solitary wave cases; on the other hand, the effect of the time elapsed during tsunami overflow cannot be reflected in the simulations using the present numerical tool. However, the stability of the overall cross-section under tsunami overflow is assessed by evaluating forces acting on the rear side armor unit supporting the crown-wall of the rubble mound breakwater as a practical engineering application in the present paper. Furthermore, two non-dimensional parameters are derived to discuss the stability of this armor unit; and thus, the stability condition of the overall cross-section. Approximate threshold values for these non-dimensional parameters are presented comparing experimental and numerical results as a starting point for engineers in practice. Finally, investigations on the solitary wave and tsunami overflow experiments/simulations are extended to the potential tsunami simulation in the scope of both representation of a realistic tsunami in a wave flume and stability of the rubble mound breakwater. 相似文献
13.
Propagation of a solitary wave over rigid porous beds 总被引：1，自引：0，他引：1
The unsteady two-dimensional Navier–Stokes equations and Navier–Stokes type model equations for porous flows were solved numerically to simulate the propagation of a solitary wave over porous beds. The free surface boundary conditions and the interfacial boundary conditions between the water region and the porous bed are in complete form. The incoming waves were generated using a piston type wavemaker set up in the computational domain. Accuracy of the numerical model was verified by comparing the numerical results with the theoretical solutions. The main characteristics of the flow fields in both the water region and the porous bed were discussed by specifying the velocity fields. Behaviors of boundary layer flows in both fluid and porous bed regions were also revealed. Effects of different parameters on the wave height attenuation were studied and discussed. The results of this numerical model indicate that for the investigated incident wave as the ratio of the porous bed depth to the fluid depth exceeds 10, any further increase of the porous bed depth has no effect on wave height attenuation. 相似文献
14.
为研究内孤立波的地形和背景流共振机制,用地形和背景流共振机制计算了3个潜标观测的内孤立波(不同模态、不同波长)的流速和传播速度,并与观测到的内孤立波进行比较。潜标观测的第一模态内孤立波(波长分别为6.4和3.3km)都是下凹型内孤立波,2个内孤立波的传播速度约为1.4m/s、最大振幅约为48m,水平流向结构都是上层西北向、下层东南向,波长3.3km 的内孤立波波峰前后有更明显的下降流和上升流。用共振机制计算出的第一模态和第二模态纬向流速的垂向结构与观测相同,最大纬向流速出现的深度与观测一致,分别相差5和12m。用共振机制计算出的内孤立波传播速度与用 KdV 方程计算的传播速度相当,共振机制计算波速为0.66~1.21m/s,KdV 方程计算波速为0.79~1.40m/s。 相似文献
15.
利用时域高阶边界元方法建立了模拟极限波浪运动的完全非线性数值模型,其中自由水面满足完全非线性自由水面条件.采用半混合欧拉-拉格朗日方法追踪流体瞬时水面,运用四阶Runge-Kutta方法更新下一时间步的波面和速度势,同时应用镜像格林函数消除水槽两个侧面和底面上的积分.研究中利用波浪聚焦的方法产生极限波浪,并且在水槽中开展了物理模型实验,将测点试验数据与数值结果进行了对比,两者吻合得很好.对极限波浪运动的非线性和流域内速度分布进行了研究. 相似文献
16.
17.
《Coastal Engineering》2001,44(2):117-139
In this paper, laboratory data for free surface displacements and velocity fields in front of a caisson breakwater covered with wave-dissipating blocks, together with wave pressures acting on the caisson, are presented and discussed. The core of the breakwater is made of a concrete caisson with a vertical front wall. The caisson is protected by a thick layer of tetrapods and is supported by a rubble mound. The breakwater is placed on the 1/25 impermeable slope. Two types of incident waves are used in the experiments: nonbreaking waves and spilling-type breaking waves. In the breaking wave case, the incident wave breaks offshore before it reaches the breakwater. The velocity data are obtained by using both the Laser Doppler Velocimeter (LDV) and the Electromagnetic Current Meter (EMCM). The raw data are analyzed using a numerical-filtering scheme so that turbulent fluctuations are separated from the phase-dependent wave motions. The vertical profiles of the time-averaged (over a wave period) turbulent velocity components at several vertical cross-sections in front of the breakwater are then analyzed. The spatial variations of the time-averaged turbulence velocity suggest that turbulence is generated inside the protective armor layer and transported into the flow region in front of the breakwater. The wave pressures on the vertical face and on the bottom of the caisson are also reported. 相似文献
18.
缓坡方程是描述近岸波浪运动较好的数学模型之一。在发展的自适应有限元求解缓坡方程的基础上,采用迭代求解的方法,确定波浪相对于边界的入射方向,从而对边界条件进行改进,建立了求解缓坡方程的数值计算模型。典型算例表明,考虑波浪相对于边界的入射角度后,模型可以更好地模拟吸收波浪边界,同时对多向波对双突堤的绕射进行了模拟研究,与试验结果比较表明,所建立的数值计算模型能够适用于多向不规则波传播过程的模拟研究。 相似文献
19.
《Coastal Engineering》2006,53(2-3):181-190
Two-dimensional depth-averaged Boussinesq-type equations were presented with the consideration of slowly varying bathymetry and effects of bottom viscous boundary layer. These Boussinesq-type equations were written in terms of the horizontal velocity components evaluated at an arbitrary elevation in the water depth and the free surface displacement. The leading order effects of the bottom boundary layer were represented by a convolution integral in the depth-integrated continuity equation. To test the validity of the theory, a set of laboratory experiments was performed to measure the viscous damping and shoaling of a solitary wave propagating in a wave tank. The time histories of the free surface profiles were measured at several locations along the centerline of the flume. To compare these laboratory data with theoretical results, the two-dimensional Boussinesq-type equations were integrated across the wave tank, resulting in a set of one-dimensional equations, while the side-wall boundary layers were properly considered. The agreement between the experimental data and numerical results was very good. The bottom shear stress formula was also given and its impact on the sediment transport rate was discussed. 相似文献
20.
植被斜坡岸滩海啸波消减数值模拟研究 总被引：1，自引：0，他引：1
An explicit one-dimensional model based on the shallow water equations(SWEs) was established in this work to simulate tsunami wave propagation on a vegetated beach. This model adopted the finite-volume method(FVM)for maintaining the mass balance of these equations. The resistance force caused by vegetation was taken into account as a source term in the momentum equation. The Harten–Lax–van Leer(HLL) approximate Riemann solver was applied to evaluate the interface fluxes for tracing the wet/dry transition boundary. This proposed model was used to simulate solitary wave run-up and long-periodic wave propagation on a sloping beach. The calibration process suitably compared the calculated results with the measured data. The tsunami waves were also simulated to discuss the water depth, tsunami force, as well as the current speed in absence of and in presence of forest domain. The results indicated that forest growth at the beach reduced wave energy loss caused by tsunamis. A series of sensitivity analyses were conducted with respect to variable parameters(such as vegetation densities, wave heights, wave periods, bed resistance, and beach slopes) to identify important influences on mitigating tsunami damage on coastal forest beach. 相似文献