Numerical study on cnoidal wave run-up around a vertical circular cylinder

A finite element model of Boussinesq-type equations was set up, and a direct numerical method is proposed so that the full reflection boundary condition is exactly satisfied at a curved wall surface. The accuracy of the model was verified in tests. The present model was used to further examine cnoidal wave propagation and run-up around the cylinder. The results showed that the Ursell number is a nonlinear parameter that indicates the normalized profile of cnoidal waves and has a significant effect on the wave run-up. Cnoidal waves with the same Ursell number have the same normalized profile, but a difference in the relative wave height can still cause differences in the wave run-up between these waves. The maximum dimensionless run-up was predicted under various conditions. Cnoidal waves hold entirely distinct properties from Stokes waves under the influence of the water depth, and the nonlinearity of cnoidal waves enhances rather than weakens with increasing wavelength. Thus, the variations in the maximum run-up with the wavelength for cnoidal waves are completely different from those for Stokes waves, and there are even significant differences in the variation between different cnoidal waves.     

Numerical simulation of nonlinear wave radiation by a moving vertical cylinder

In this paper the fully nonlinear potential model based on a finite element method is used to investigate the nonlinear wave motion around a moving circular cylinder. The results for the cylinder in transient motion are compared with the experimental data and a much better agreement than the linear theory is found. Further simulation for a circular cylinder in sinusoidal motion is made. It is found that when the ratio of the cylinder diameter D to the wavelength L is relatively small at a fixed motion amplitude the nonlinear components of the runup on the cylinder surface at the second- and third-harmonic frequencies become more important and this is confirmed by the experimental data. Results for the hydrodynamic force are also provided for a cylinder oscillating in a channel. It is noticed that when the frequency of the cylinder motion in a channel is between the first and the second natural frequencies of the symmetric mode, the time history has components not only at the frequency of the cylinder motion but also at the first natural frequency. The latter remains significant over the period that the simulation is made. This has important implications to model testing. If measurement is to be made at such a frequency it may take long time for the motion to become periodic at the frequency of the cylinder motion.     

透空式防波堤周围的非线性波浪传播的数值模拟

对非线形波浪在透空式防波堤周围的波浪变形进行了数值模拟,在Boussinesq波浪方程中加入与透空建筑物有关的新的耗散项,从而界定了透空建筑物引起的部分反射和透射,波浪折射衍射的传播过程通过控制方程求解。波浪控制方程通过有限差分方法求解。模型应用于模拟波浪经过具有部分反射的群桩式透空结构,结果表明透空式防波堤可以有效地衰减波浪,是重力式结构的一种替代形式。     

Vortex shedding suppression for flow over a circular cylinder near a plane boundary

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.     

片状分布植被对孤立波爬高影响的数值模拟研究

近岸植被对波浪爬坡具有一定的衰减作用。在自然界中,由于植被的死亡、再生或人为破坏等原因,近岸植被通常呈片状分布,且其内部分布也是不均匀的。本文以完全非线性Boussinesq方程为基础,引入植被作用项,建立了模拟近岸植被区波浪传播的数值模型,验证了模型可靠性,进而采用该模型模拟分析了片状分布植被对孤立波爬高的影响。数值模拟结果表明,片状分布植被能有效减小孤立波爬高;对于均匀分布的片状植被,高密度片状植被对孤立波爬高的消减效果优于低密度片状植被;对于相同密度、不同分布形式的片状植被,均匀分布的片状植被对孤立波的消减效果优于不均匀分布的片状植被;对于不均匀分布的片状植被,前密后疏的片状植被对孤立波的消减效果优于前疏后密的片状植被。     

Numerical investigation of tsunami-like wave hydrodynamic characteristics and its comparison with solitary wave

Solitary waves have been commonly used as an initial condition in the experimental and numerical modelling of tsunamis for decades. However, the main component of a tsunami waves acts at completely different spatial and temporal scales than solitary waves. Thus, use of solitary waves as approximation of a tsunami wave may not yield realistic model results, especially in the coastal region where the shoaling effect restrains the development of the tsunami wave. Alternatively, N-shaped waves may be used to give a more realistic approximation of the tsunami wave profile. Based on the superposition of the sech²(*) waves, the observed tsunami wave profile could be approximated with the N-shaped wave method, and this paper presents numerical simulation results based on the tsunami-like wave generated based on the observed tsunami wave profile measured in the Tohoku tsunami. This tsunami-like wave was numerically generated with an internal wave source method based on the two-phase incompressible flow model with a Volume of Fluid (VOF) method to capture the free surface, and a finite volume scheme was used to solve all the governing equations. The model is first validated for the case of a solitary wave propagating within a straight channel, by comparing its analytical solutions to model results. Further, model comparisons between the solitary and tsunami-like wave are then made for (a) the simulation of wave run-up on shore and (b) wave transport over breakwater. Comparisons show that use of these largely different waveform shapes as inputs produces significant differences in overall wave evolution, hydrodynamic load characteristics as well as velocity and vortex fields. Further, it was found that the solitary wave uses underestimated the total energy and hence underestimated the run-up distance.     

Steady streaming around a circular cylinder in an oscillatory flow

Steady streaming induced by an oscillatory flow around a circular cylinder is investigated using a numerical method. Two-dimensional Reynolds-averaged Navier-Stokes equations are solved using a finite element method with a k-ω turbulent model closure. The range of the Keulegan-Carpenter (KC) number investigated is between 2 and 40, which is substantially higher than those reported in literature related to steady streaming to date. A constant value of Stokes number (β) of 196 is chosen in this study. The steady streaming structures and velocity distribution are analysed in detail. It is found that the characteristics of steady streaming are strongly related to the vortex shedding flow regimes.     

基于局部冲刷研究的大型圆柱周围波浪特性分析

波浪环境下,大型圆柱的存在改变了周围原有波浪条件,直接导致局部冲刷。利用MacCamy和Fuchs的绕射理论对相对柱径较大(0.15相似文献   

Numerical simulation of solitary wave propagation based on MPS method

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.     

基于MPS无网格方法的孤立波传播的数值模拟

在MPS无网格方法中,引进预定候选粒子集概念用以生成邻接粒子集矩阵,使该部分的机时耗费缩短为引进前的1/11;采用Bi-CGSTAB方法求解压力泊松方程,显著地提高了求解速度.模拟了孤立波在数值波浪水槽中的传播及其与直墙作用时的爬升、回落过程,结果表明模拟波面结果与解析值及实测结果基本相符,针对不同波高的孤立波计算得到的墙前最大爬升值与实测结果也是一致的.     

孤立波与带窄缝双箱相互作用模拟研究

针对孤立波与带窄缝双箱的作用问题,应用时域高阶边界元方法建立了二维数值水槽。其中,自由水面满足完全非线性运动学和动力学边界条件,对瞬时自由表面流体质点采用混合欧拉-拉格朗日法追踪,采用四阶龙格库塔法对下一时刻的自由水面的速度势和波面升高进行更新。采用加速度势法求解物体湿表面的瞬时波浪力。采用推板方法生成孤立波。通过模拟孤立波在直墙上的爬高以及施加在直墙上的波浪力,并与已发表的实验和数值结果对比,验证本数值模型的准确性。通过数值模拟计算研究了窄缝宽度、方箱尺寸对波浪在箱体迎浪侧爬高,窄缝内波面升高,箱体背浪侧透射波高及箱体受波浪荷载的影响。同时研究了有一定时间间隔的双孤立波与带窄缝双箱系统作用问题。     

完全非线性波浪破碎模型沿岸流数值模拟

探讨一种基于完全非线性Boussinesq方程的波浪破碎模型在沿岸流计算中的应用问题。针对控制方程中的完全非线性项对沿岸流成长过程的影响进行了深入讨论。数学模型计算结果表明,完全非线性项有使平均流局部化的作用;通过数模实验还发现,垂向高阶涡度项可以有效抑制破波区外回流;运用Visser的实验室沿岸流实测资料从沿岸流速度、波高和平均水位几方面对所提模型进行了验证,并给出了紊动参数的计算结果。     

Large-eddy simulation of the turbulent near wake behind a circular cylinder: Reynolds number effect

The purpose of the present work is <!–<query>The highlights are in an incorrect format. Hence they have been deleted. Please refer the online instructions: http://www.elsevier.com/highlights and provide 3-5 bullet points.</query>–>to study the effect of the Reynolds number on the near-wake structure and separating shear layers behind a circular cylinder. Three-dimensional unsteady large-eddy simulation is carried out and two different subgrid scale models are applied in order to evaluate the turbulent wake reasonably. The Reynolds number based on the free-stream velocity and the cylinder diameter is ranging from Re = 5500–41,300, corresponding to the full development of the shear-layer instability in the intermediate subcritical flow regime. For a complete validation of this numerical study, hydrodynamic bulk coefficients are computed and compared to experimental measurements and numerical studies in the literature. Special focus is made on the variations of both the large-scale near-wake structure and the small-scale shear-layer instability with increasing Reynolds numbers. The present numerical study clearly shows the broadband nature of the shear-layer instability as well as the dependence of the shear-layer frequency especially on the high Reynolds numbers.     

A new form of generalized Boussinesq equations for varying water depth

A new set of equations of motion for wave propagation in water with varying depth is derived in this study. The equations expressed by the velocity potentials and the wave surface elevations include first-order non-linearity of waves and have the same dispersion characteristic to the extended Boussinesq equations. Compared to the extended Boussinesq equations, the equations have only two unknown scalars and do not contain spatial derivatives with an order higher than 2. The wave equations are solved by a finite element method. Fourth-order predictor–corrector method is applied in the time integration and a damping layer is applied at the open boundary for absorbing the outgoing waves. The model is applied to several examples of wave propagation in variable water depth. The computational results are compared with experimental data and other numerical results available in literature. The comparison demonstrates that the new form of the equations is capable of calculating wave transformation from relative deep water to shallow water.     

Breaking solitary wave evolution over a porous underwater step

Solitary wave evolution over a shelf including porous damping is investigated using Volume-Averaged Reynolds Averaged Navier–Stokes equations. Porous media induced damping is determined based on empirical formulations for relevant parameters, and numerical results are compared with experimental information available in the literature. The aim of this work is to investigate the effect of wave damping on soliton disintegration and evolution along the step for both breaking and non-breaking solitary waves. The influence of several parameters such as geometrical configuration (step height and still water level), porous media properties (porosity and nominal diameter) or solitary wave characteristics (wave height) is analyzed. Numerical simulations show the porous bed induced wave damping is able to modify wave evolution along the step. Step height is observed as a relevant parameter to influence wave evolution. Depth ratio upstream and downstream of the edge appears to be the more relevant parameter in the transmission and reflection coefficients than porosity or the ratio of wave height–water depth. Porous step also modifies the fission and the solitary wave disintegration process although the number of solitons is observed to be the same in both porous and impermeable steps. In the absence of breaking, porous bed triggers a faster fission of the incident wave into a second and a third soliton, and the leading and the second soliton reduces their amplitude while propagating. This decrement is observed to increase with porosity. Moreover, the second soliton is released before on an impermeable step. Breaking process is observed to dominate over the wave dissipation at the porous bottom. Fission is first produced on a porous bed revealing a clear influence of the bottom characteristics on the soliton generation. The amplitude of the second and third solitons is very similar in both impermeable and porous steps but they evolved differently due to the effect of bed damping.     

Large-eddy simulation of the flow past a circular cylinder at sub- to super-critical Reynolds numbers

Large-eddy simulation of turbulent flow past a circular cylinder at sub- to super-critical Reynolds numbers is performed using a high-fidelity orthogonal curvilinear grid solver. Verification studies investigate the effects of grid resolution, aspect ratio and convection scheme. Monotonic convergence is achieved in grid convergence studies. Validation studies use all available experimental benchmark data. Although the grids are relatively large and fine enough for sufficiently resolved turbulence near the cylinder, the grid uncertainties are large indicating the need for even finer grids. Large aspect ratio is required for sub-critical Reynolds number cases, whereas small aspect ratio is sufficient for critical and super-critical Reynolds number cases. All the experimental trends were predicted with reasonable accuracy, in consideration the large facility bias, age of most of the data, and differences between experimental and computational setup in particular free stream turbulence and roughness. The largest errors were for under prediction of turbulence separation.     

Wave run-up on a coaxial perforated circular cylinder

This paper describes a plane regular wave interaction with a combined cylinder which consists of a solid inner column and a coaxial perforated outer cylinder.The outer perforated surface is a thin porous cylinder with an annular gap between it and the inner cylinder.The non-linear boundary condition at the perforated wall is a prime focus in the study;energy dissipation at the perforated wall occurs through the resistance to the fluid across the perforated wall.Explicit analytical formulae are presented to calculate the wave run-up on the outer and inner surfaces of the perforated cylinder and the surface of the inner column.The theoretical results of the wave run-up are compared with previous experimental data.Numerical results have also been obtained:when the ratio of the annular gap between the two cylinders to incident wavelength(b-a)/L≤0.1,the wave run-up on the inner surface of the perforated cylinder and the surface of inner column can partially or completely exceed the incident wave height.     

Numerical simulation of irregular wave overtopping against a smooth sea dike

Based on the filtered Navier-Stokes equations and Smagorinsky turbulence model,a numerical wave flume is developed to investigate the overtopping process of irregular waves over smooth sea dikes.Simulations of fully nonlinear standing wave and regular wave’s run-up on a sea dike are carried out to validate the implementation of the numerical wave flume with wave generation and absorbing modules.To model stationary ergodic stochastic processes,several cases with different random seeds are computed for each specified irregular wave spectrum.It turns out that the statistical mean overtopping discharge shows good agreement with empirical formulas,other numerical results and experimental data.     

多向不规则波绕射的数值模拟

采用有限元求解Ｂｏｕｓｓｉｎｅｓｑ方程数值计算模型,对单突堤和双突堤的波浪绕射进行了数值模拟,并与物模试验进行对比验证。模拟结果表明入射波浪的方向分布对波浪的绕射具有明显的影响,使防波堤后的波浪增大,开阔区域的浪浪减小,波浪的方向分布对港内波浪的分布具有一定的均化作用。