共查询到17条相似文献，搜索用时 156 毫秒
1.
应用经CSPM法和黎曼解修正后的光滑粒子流体动力学(SPH)方法,建立了主动吸收无反射数值波浪水槽,研究波浪作用下多孔介质结构的水动力特性。流体运动控制方程采用N-S方程,多孔渗水结构内流体的运动控制方程考虑渗流力的影响。数值计算结果给出了水槽内不同位置测点的波面历时曲线和越浪量随时间变化曲线,并同试验结果和Philip Liu的数值计算结果进行了比较。并对一个波浪周期内斜坡堤多孔介质结构内外的速度场和压力分布进行了讨论分析。计算分析表明,数值计算波面较Philip Liu的计算结果与试验结果吻合更好。说明应用SPH方法建立的二维数学模型能够较好地模拟破碎波在多孔渗水斜坡上的爬坡和越浪。 相似文献
2.
近岸植被对波浪爬坡具有一定的衰减作用。在自然界中，由于植被的死亡、再生或人为破坏等原因，近岸植被通常呈片状分布，且其内部分布也是不均匀的。本文以完全非线性Boussinesq方程为基础，引入植被作用项，建立了模拟近岸植被区波浪传播的数值模型，验证了模型可靠性，进而采用该模型模拟分析了片状分布植被对孤立波爬高的影响。数值模拟结果表明，片状分布植被能有效减小孤立波爬高；对于均匀分布的片状植被，高密度片状植被对孤立波爬高的消减效果优于低密度片状植被；对于相同密度、不同分布形式的片状植被，均匀分布的片状植被对孤立波的消减效果优于不均匀分布的片状植被；对于不均匀分布的片状植被，前密后疏的片状植被对孤立波的消减效果优于前疏后密的片状植被。 相似文献
3.
波浪对透空式结构物冲击作用的光滑粒子流体动力学数值模拟 总被引：3，自引：0，他引：3
基于光滑粒子流体动力学(SPH)方法的二维数值波浪水槽模拟了规则波对透空式结构物的冲击作用。通过黎曼解和CSPM相结合的方法对连续方程和动量方程进行了修正。对造波边界采用虚粒子法模拟。提出了一种耦合计算方法来消除造波边界附近的压力波动现象,在结构物边界处设置了适合的冲击边界条件。应用修正的SPH模型模拟了规则波对浪溅区结构物的冲击作用,给出了结构物附近流场和压力场的变化特征,并应用物理模型试验结果对数值模型进行了验证。 相似文献
4.
5.
6.
7.
本文采用圆柱体阵列来模拟珊瑚礁面的大糙率,通过波浪水槽实验研究礁面糙率对孤立波传播变形及岸滩爬高的影响。结果表明,粗糙礁面的存在显著削弱了礁坪上孤立波的首峰和礁后岸滩反射造成的次峰,同时降低了波浪在珊瑚礁面的传播速度;垂直于岸线方向沿礁相对波高随着入射波增大而减小,随着礁坪水深的增大而增大,粗糙礁面上波高沿礁的衰减更为显著;礁前斜坡的无量纲反射波高随无量纲入射波高的变化与礁坪水深有关,当入射波高足够大时其趋于常值,粗糙礁面略微增大了礁前斜坡的反射;无量纲透射波高和岸滩爬高随着无量纲入射波高的增大而减小,特别是礁坪水深较大时更为显著,粗糙礁面时的无量纲岸滩爬高相对于光滑礁面平均减小46%;通过回归分析得出了同时适合于光滑和粗糙礁面的预测孤立波岸滩爬高的经验关系式。 相似文献
8.
为准确模拟孤立波在岸礁地形上的传播和爬坡,采用基于完全非线性Boussinesq方程开发的Funwave-TVD模型,探究模型的可行性,并利用验证后的模型进一步研究岸礁各地形因素对孤立波爬高的影响。研究结果表明:模型能准确模拟孤立波在岸礁陡变地形上的传播及变形,摩擦系数对礁前陡坡及礁坪上的波浪传播模拟影响不大,但对爬坡预测的敏感性较强;模型空间步长可适当增大,提高计算效率;随着礁坪宽度的增大以及礁后斜坡的变缓,孤立波爬坡高度下降明显,而礁前陡坡坡度变化对孤立波爬坡高度影响不大。 相似文献
9.
孤立波与带窄缝双箱相互作用模拟研究 总被引：1，自引：1，他引：0
针对孤立波与带窄缝双箱的作用问题,应用时域高阶边界元方法建立了二维数值水槽。其中,自由水面满足完全非线性运动学和动力学边界条件,对瞬时自由表面流体质点采用混合欧拉-拉格朗日法追踪,采用四阶龙格库塔法对下一时刻的自由水面的速度势和波面升高进行更新。采用加速度势法求解物体湿表面的瞬时波浪力。采用推板方法生成孤立波。通过模拟孤立波在直墙上的爬高以及施加在直墙上的波浪力,并与已发表的实验和数值结果对比,验证本数值模型的准确性。通过数值模拟计算研究了窄缝宽度、方箱尺寸对波浪在箱体迎浪侧爬高,窄缝内波面升高,箱体背浪侧透射波高及箱体受波浪荷载的影响。同时研究了有一定时间间隔的双孤立波与带窄缝双箱系统作用问题。 相似文献
10.
波浪的方向分布对波浪的传播及其与工程结构物的作用都具有明显影响，目前现有的研究大多是基于单向波浪进行的。为了研究方向分布对群墩结构上的爬高影响，基于规则波浪与群墩作用的理论解，结合多向不规则波浪的造波方法，建立了多向不规则波浪与群墩作用的计算模型，同时进行了物理模型试验对模型的有效性进行了验证。系统地对群墩周围及表面上的波浪爬高进行了计算分析，结果表明，方向分布对波浪爬高具有较大的影响，且不同位置处的影响并不相同，在实际的工程设计中如果按照单向波浪计算，可能低估或者高估群墩周围的爬高。 相似文献
11.
The concept of candidate particle set is introduced in the MPS gridless numerical method to generate neighboring particle set matrix, which can reduce the CPU time to 1/11 of that before introduction. The Bi-CGSTAB (bi-conjugate gradient stabilized) algorithm is applied to solving the Poisson pressure equation, by which the solving speed is significantly accelerated. The process of solitary waves propagating over a numerical flume and interacting with a vertical wall is simulated. The simulated results of water surface elevation are in good agreement with the analytical solution as well as the measured data. The predicted maximum values of the run-up of solitary waves with various relative incident wave heights agree well with the measured results. 相似文献
12.
Although the finite difference method is computationally efficient, it is acknowledged to be inferior when dealing with flow-over on structures with a complex geometry because of its rectilinear grid system. Therefore, we developed a numerical procedure that can cope with flow over structures with complex shapes while, at the same time, retaining the simplicity and efficiency of a rectilinear grid system. We used the immersed boundary method, which involves application of immersed boundary forces at solid boundaries rather than conventional boundary conditions, to investigate wave interactions with coastal structures in a three-dimensional numerical wave tank by solving the Navier–Stokes equations for two-phase flows. We simulated the run-up of a solitary wave around a circular island. Maximum run-up heights were computed around the island and compared with available laboratory measurements and previous numerical results. The three-dimensional features of the run-up process were analyzed in detail and compared with those of depth-integrated equations models. 相似文献
13.
The boundary integral equation method (BIEM) is developed as a tool for studying two-dimensional, nonlinear water wave problems, including the phenomena of wave generation, propagation and run-up. The wave motions are described by a potential flow theory. Nonlinear free-surface boundary conditions are incorporated in the numerical formulation. Examples are given for either a solitary wave or two successive solitary waves. Special treatment is developed to trace the run-up and run-down along a shoreline. The accuracy of the present scheme is verified by comparing numerical results with experimental data of maximum run-up. 相似文献
14.
Smoothed Particle Hydrodynamics method (SPH) has a good adaptability for simulating of free surface flow problems. However, there are some shortcomings of SPH which are still in open discussion. This paper presents a corrected solid boundary handling method for weakly compressible SPH. This improved method is very helpful for numerical stability and pressure distribution. Compared with other solid boundary handling methods, this corrected method is simpler for virtual ghost particle interpolation and the ghost particle evaluation relationship is clearer. Several numerical tests are given, like dam breaking, solitary wave impact and sloshing tank waves. The results show that the corrected solid boundary processing method can recover the spurious oscillations of pressure distribution when simulating the problems with complex geometry boundary. 相似文献
15.
This paper presents a numerical model for simulating wave interaction with porous structures. Incompressible smoothed particle hydrodynamics in porous media (ISPHP) method is introduced in this study as a mesh free particle approach that is capable of efficiently tracking the large deformation of free surfaces in a Lagrangian coordinate system. The developed model solves two porous and pure fluid flows simultaneously by means of one equation that is equivalent to the unsteady 2D Navier–Stokes (NS) equations for the flows outside the porous media and the extended Forchheimer equation for the flows inside the porous media. Interface boundary between pure fluid and porous media is effectively modeled by the SPH integration technique. A two-step semi-implicit scheme is also used to solve the fluid pressure satisfying the fluid incompressibility criterion.The developed ISPHP model is then validated via different experimental and numerical data. Fluid flow pattern through porous dam with different porosities is studied and regular wave attenuation over porous seabed is investigated. As a practical case, wave running up and overtopping on a caisson breakwater protected by a porous armor layer are modeled. The results show good agreements between numerical and laboratory data in terms of free surface displacement, overtopping rate and pressure distribution. Based on this study, ISPHP model is an efficient method for simulating the coastal applications with porous structures. 相似文献
16.
Benchmark computations of wave run-up on single cylinder and four cylinders by naoe-FOAM-SJTU solver
The benchmark simulations of wave run-up on a fixed single truncated circular cylinder and four circular cylinders are presented in this paper. Our in-house CFD solver naoe-FOAM-SJTU is adopted which is an unsteady two-phase CFD code based on the open source package OpenFOAM. The Navier-Stokes equations are employed as the governing equations, and the volume of fluid (VOF) method is applied for capturing the free surface. Monochromatic incident waves with the specified wave period and wave height are simulated and wave run-up heights around the cylinder are computed and recorded with numerical virtual wave probes. The relationship between the wave run-up heights and the incident wave parameters are analyzed. The numerical results indicate that the presented naoe-FOAM-SJTU solver can provide accurate predictions for the wave run-up on one fixed cylinder and four cylinders, which has been proved by the comparison of simulated results with experimental data. 相似文献
17.
The form of Boussinesq equation derived by Nwogu （1993） using velocity at an arbitrary distance and surface elevation as variables is used to simulate wave surface elevation changes. In the numerical experiment, water depth was divided into five layers with six layer interfaces to simulate velocity at each layer interface. Besides, a physical experiment was carried out to validate numerical model and study solitary wave propagation.“Water column collapsing”method （WCCM） was used to generate solitary wave. A series of wave gauges around an impervious breakwater were set-up in the flume to measure the solitary wave shoaling, run-up, and breaking processes. The results show that the measured data and simulated data are in good agreement. Moreover, simulated and measured surface elevations were analyzed by the wavelet transform method. It shows that different wave frequencies stratified in the wavelet amplitude spectrum. Finally, horizontal and vertical velocities of each layer interface were analyzed in the process of solitary wave propagation through submerged breakwater. 相似文献