A highly wave dissipation offshore breakwater

The aim of this paper is to develop an offshore breakwater, for which coefficients of both the wave reflection and transmission have low values. The breakwater is suggested to compose of n layers of porous materials with different porosities. A complex eigen function method is used in the theoretical analysis. Continuities of both mass flux and fluid pressure are assumed at interfaces between every two adjoining porous materials and at the interface between end materials and water region. Following a series of mathematical processes, the coefficients of the wave transmission and reflection along with the wave energy loss are calculated. The porosity of materials is varied in computations; and results are compared among structures composing of different layers of porous materials. A single layer offshore breakwater is shown to reduce simultaneously the coefficients of transmission and reflection only when the structure is very wide in the direction of wave propagation, and the structure material has a high porosity. A multilayer breakwater, however, can function well in reducing both coefficients at a much narrower width; structure having more layers can be more effective at narrower width. Finally, several experiments are conducted; theoretical computations and experimental results agree well.     

Modified moving particle method for modeling wave interaction with multi layered porous structures

Modified Moving Particle method in Porous media (MMPP) is introduced in this study for simulating a flow interaction with porous structures. By making use of the sub-particle scale (SPS) turbulence model, a unified set of equations are introduced for the entire computational domain and a proper boundary treatment is suggested at the interfaces between fluid and the porous media. Similar to the Incompressible Smoothed Particle Hydrodynamic (ISPH) method, a robust two-step semi-implicit scheme is utilized to satisfy the incompressibility criterion. By means of the introduced model, different flow regimes through multi-layered porous structures with arbitrary shapes can be simulated and there is no need to implement calibration factors.The developed MMPP model is then validated via simulating the experiments of Liu et al. (1999) i.e. linear and turbulent flows through porous dams and the experiments of Sakakiyama and Liu (2001) i.e. wave overtopping on a caisson breakwater protected by multi layered porous materials. Good agreements between numerical and laboratory data present the ability of the introduced model in simulating various flow regimes through multi-layered porous structures. It is concluded that the turbulent flow is an important issue particularly at the interface between the free fluid and porous media and consequently, the accuracy of the previous Lagrangian models that were based on neglecting the turbulence effect can be improved significantly by means of the present model. In addition, to satisfy the continuity criteria in the SPH models, it is necessary to modify density of particles in accordance with their porosity.     

Short-crested waves interaction with a concentric porous cylinder system with partially porous outer cylinder

In this paper, based on the linear wave theory, the interaction of short-crested waves with a concentric dual cylindrical system with a partially porous outer cylinder is studied by using the scaled boundary finite element method (SBFEM), which is a novel semi-analytical method with the advantages of combining the finite element method (FEM) with the boundary element method (BEM). The whole solution domain is divided into one unbounded sub-domain and one bounded sub-domain by the exterior cylinder. By weakening the governing differential equation in the circumferential direction, the SBFEM equations for both domains can be solved analytically in the radial direction. Only the boundary on the circumference of the exterior porous cylinder is discretized with curved surface finite elements. Meanwhile, by introducing a variable porous-effect parameter G, non-homogeneous materials caused by the complex configuration of the exterior cylinder are modeled without additional efforts. Comparisons clearly demonstrate the excellent accuracy and computational efficiency associated with the present SBFEM. The effects of the wide range wave parameters and the structure configuration are examined. This parametric study will help determine the various hydrodynamic effects of the concentric porous cylindrical structure.     

波浪与外壁透空双方形沉箱相互作用

在线性波浪理论下,利用复合边界元素法(composite BEM)数值解析在等水深、规则波浪入射二种外壁透空双方形沉箱的无因次波力及Kd绕射分布图,并分别与其他研究者所作双圆筒内、外圆柱无因次波力及绕射分布图的计算结果进行比较,都说明本模式的合理性和可行性。在考虑不同透水参数下,分别计算波浪作用在内、外结构物的波力及沉箱四周绕射系数大小分布图。计算结果显示:波浪作用于外壁透空全透水双方形沉箱,外方柱无因次最大波力值会随着透水参数增加而降低;内方柱无因次最大波力值则随着透水参数增加而增加。本研究结果可供设计外海透水方形沉箱结构交互作用参考。     

A coupled numerical model of wave interaction with porous medium

A new coupling model of wave interaction with porous medium is established in which the wave field solver is based on the two dimensional Reynolds Averaged Navier-Stokes (RANS) equations with a kε closure. Incident waves, which could be linear waves, cnoidal waves or solitary waves, are produced by a piston-type wave maker in the computational domain and the free surface is traced through the Piecewise Linear Interface Construction-Volume of Fluid (PLIC-VOF) method. Nonlinear Forchheimer equations are adopted to calculate the flow field within the porous media. By introducing a velocity–pressure correction equation, the wave field and the porous flow field are highly and efficiently coupled. The two fields are solved simultaneously and no boundary condition is needed at the interface of the internal porous flow and the external wave. The newly developed numerical model is used to simulate wave interaction with porous seabed and the numerical results agree well with the experimental data. The additional numerical tests are also conducted to study the effects of seabed thickness, porosity and permeability coefficient on wave damping and the pore water pressure responses.     

Boundary element analysis of liquid sloshing characteristics in axisymmetric tanks with various porous baffles

This paper aims to investigate the effects of the porous baffles on the suppression of sloshing for the tanks with axisymmetric geometries under lateral excitation. Based on the assumptions of inviscid, irrotational, incompressible liquid and small amplitude sloshing, an axisymmetric boundary element method (BEM) for 3D Laplace equation is derived by using the Green's theorem together with the weighted residual method. And a zoning method is employed to model fluid domain in the tanks with complex porous baffles. Meanwhile, the porous baffles are treated motioning together with the tanks, and the velocity across the porous baffle is assumed to be linearly proportional to the pressure gradient between each side of the porous baffle. And the mechanism of suppressing the sloshing response is mainly the energy dissipation of the fluid passing through the porous baffle. Moreover, the linear free surface boundary conditions are also used to solve the governing equations. Compared with other numerical methods, the most prominent advantage of the BEM in solving axisymmetric potential problem is that only the boundaries of half the cross-section instead of the entire problem domain should be discretized, which can cut down large amount of memory and time costs. The present method is verified by comparing the numerical results with the existing literatures, and excellent agreements are obtained. Meanwhile, the proposed models are applied to investigate the effects of the porous baffles on sloshing response in circular cylindrical, annular cylindrical and conical tanks. The effects of the porous baffle length, porous-effect parameter, installation angle and baffle height on the sloshing force, natural frequency and surface elevation are studied. Additionally, some typical sloshing pressure distributions, velocity potential contours and velocity fields are plotted. The results show that swirls at the tips of the baffles can be observed in many cases, and the top-mounted porous baffle makes more significant suppression effects on sloshing response than that of bottom-mounted porous baffle, while increasing the number of ring porous baffles can achieve better restraint effects on sloshing response. And increasing the baffle length of the horizontal wall-mounted ring porous baffle can significantly decrease the sloshing frequencies, as well as the first non-dimensional natural frequency decreases with decrease in porous-effect parameter of the coaxial porous baffle. In addition, remarkable effects on sloshing can be obtained when reasonable designed by selecting the optimal porous-effect parameter, installation angle and baffle height. And this paper can be a useful guide for the seismic design and analysis of many actual liquid storage tanks (such as the Advanced Passive PWR, large water cooling tower, etc.).     

波浪与多孔介质结构相互作用SPH模拟

应用经CSPM法和黎曼解修正后的光滑粒子流体动力学(SPH)方法,建立了主动吸收无反射数值波浪水槽,研究波浪作用下多孔介质结构的水动力特性。流体运动控制方程采用N-S方程,多孔渗水结构内流体的运动控制方程考虑渗流力的影响。数值计算结果给出了水槽内不同位置测点的波面历时曲线和越浪量随时间变化曲线,并同试验结果和Philip Liu的数值计算结果进行了比较。并对一个波浪周期内斜坡堤多孔介质结构内外的速度场和压力分布进行了讨论分析。计算分析表明,数值计算波面较Philip Liu的计算结果与试验结果吻合更好。说明应用SPH方法建立的二维数学模型能够较好地模拟破碎波在多孔渗水斜坡上的爬坡和越浪。     

波浪在局部可渗透水平海床上传播的解析解

本文建立了波浪在局部可渗透水平海床上传播的解析解,并研究了波浪在局部可渗透海床上的透射、反射问题。研究中将计算域划分为4个区域,中间区域为流域,海底可渗透,其下区域为多孔介质海床,左右两个区域也为流域,但海底不可渗透。应用线性波浪理论,建立了各流域包含非传播模态的速度势表达式,给出了海床内部的压强表达式,并利用交界面上匹配条件,求解了表达式中的待定系数。基于该解析模型,探讨了海床渗透系数、相对水深、渗透海床长度对波浪传播变形的影响。结果表明,波高沿程衰减,强度随渗透系数、渗透海床长度的增加以及相对水深的减小而变大;局部可渗透海床会引起波浪的反射和透射,随着海床长度的增加,反射系数振荡变化,并最终趋于常数,透射系数指数衰减,并最终趋于0。     

Moving particle method for modeling wave interaction with porous structures

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.     

Wave-induced mean flows in vertical rubble mound structures

Waves impinging on rubble mound breakwaters and seawalls induce a mean flow within the breakwater, analogous to the so-called undertow within the surf zone. Here, using a plane wave approximation (kh<1.5), a second-order problem is solved for an idealized breakwater with a rectangular cross-section to show the origin and the nature of the mean flow within the porous structure. The mean flow is expressed in terms of a mean stream function analytically derived, obtained based on the mass flux balance between the incident, reflected and transmitted waves. Furthermore, the evolution of other second-order magnitudes such as mean water level and mass flux is analyzed under different incident wave conditions, structure geometry and porous material characteristics. Results show that the evolution of the different mean quantities is controlled mainly by reflection and consequently depends highly on structure geometry and porous material characteristics. Furthermore, it is shown that the return flow is stronger with increasing mass flux decay. Some qualitative experiments to show the described mechanism are also presented.     

Pore pressure loads on gravity structures

It is common practice to compute wave-induced loads on the immersed surface of gravity structures exposed to the wave motion and disregard the pore-water pressure variation on the foundation surface. However, when the soil is porous, wave-induced pressures propagate within the soil under the structure and result in a rather significant contribution to overall loads. This paper describes a practical method for numerical modeling of the pore pressure under a gravity platform foundation for compressible water and a rigid, but porous soil. The porous soil may be bounded by an impermeable horizontal layer at some arbitrary depth.

The paper outlines the basic boundary element procedure for pore pressure analysis and presents numerical results for a typical gravity structure as well as results for comparison with an existing analytical solution for a vertical circular cylinder.     

A new approximate analytic solution for water wave scattering by a submerged horizontal porous disk

This study gives a new approximate analytic solution for water wave scattering by a submerged horizontal porous disk in the context of the linear potential theory. The solution is based on the domain decomposition method. The velocity potentials are determined by two different approaches. One approach is to adopt decompositions for velocity potentials, and the other is to expand the vertical derivative of the velocity potential on the porous disk along the radial direction. Hence the velocity potentials are determined by the matched eigenfunction expansions. Differing from previous solutions with respect to the porous disk, the present solution needs no complex dispersion relations. Thus the new solution is easier for numerical implementation. According to numerical examples, the convergence of the present solution is satisfactory. In addition, the present predictions of the wave surface elevation and the vertical wave force on the disk agree very well with previous results by different approaches. The present solution can also be extended to other structures involving disks, such as a fish cage, a porous disk with finite thickness, and a submerged elastic disk.     

Wave motion over a breakwater system of a horizontal plate and a vertical porous wall

A two-dimensional analytical solution is presented to study the reflection and transmission of linear water waves propagating past a submerged horizontal plate and through a vertical porous wall. The velocity potential in each fluid domain is formulated using three sets of orthogonal eigenfunctions and the unknown coefficients are determined from the matching conditions. Wave elevations and hydrodynamic forces acting on the porous wall are computed. Reflection and transmission coefficients are presented to examine the performance of the breakwater system. The present analytical solutions are found in fairly good agreement with the available laboratory data. The results indicate that the plate length, the porous-effect, the gap between plate and porous wall, and the submerged depth of the plate all show a significant influence on the reflected and transmitted wave fields. It is also interesting to note that the submerged plate plays an important role in reducing the transmitted wave height, especially for long incident waves.     

Wave reflection from vertical breakwater with porous structure

This paper presents numerical solutions for the wave reflection from submerged porous structures in front of the impermeable vertical breakwater. A new time-dependent mild-slope equation involves the parameters of the porous medium including the porosity, the friction factor and the inertia coefficient, etc. is derived for solving the boundary value problem. A comprehensive comparison between the present model and the existing analytical solution for the case of simple rectangular geometries of the submerged structure is performed first. Then, more complicated cases such as the inclined and trapezoidal submerged porous structures in front of the vertical breakwater with sloping bottom are considered. This study also examines the effects of the permeable properties and the geometric configurations of the porous structure to the wave reflection. It is found that the submerged porous structure with trapezoidal shape has more efficiency to reduce the wave reflection than that of triangular shape. The numerical results show that the minimum wave reflection is occurred when the breakwater is located at the intermediate water depth.     

Wave interaction with cylindrical porous piles

The interaction of waves with arrays of porous circular cylinders is studied theoretically and, under the assumption of potential flow and linear wave theory, an analytical solution is derived. The solution is valid for either submerged or emerged structures. The extension to the cases of unidirectional and multidirectional waves is obtained by means of a transfer function. For specific conditions the model gives the same solution as those previously presented by other authors. Numerical results are presented which exemplify diverse wave and mechanical parameters on the wave transformation due to the presence of a system of circular cylinders.     

On waves propagating over a submerged poro-elastic structure

In this paper, the problem of incident waves propagating over a submerged poro-elastic structure is studied theoretically. A linear wave theory is used to describe the wave motion. The submerged poro-elastic structure is modeled based on Biot's theory, in which the fluid motion is described using the potential wave theory of Sollitt and Cross (1972). In the present approach, the problem domain is divided into four subregions. Using general solutions for each region and matching dynamic and kinematic conditions for neighboring regions, analytic solutions are derived for the wave fields and poro-elastic structure. The present analytic solutions compare very well with simplified cases of impermeable, rigid structures, and with those of porous structures. Using the present analytic solution, the effects of a poro-elastic submerged structure on waves are studied. The results show that softer poro-elastic structures can induce higher reflection and lower transmission from incident waves. For low permeability conditions, the elasticity of the structure can induce resonance, while higher permeability can depress the resonant effects.     

Wave motion over a submerged breakwater with an upper horizontal porous plate and a lower horizontal solid plate

This study examines the hydrodynamic performance of a modified two-layer horizontal-plate breakwater. The breakwater consists of an upper submerged horizontal porous plate and a lower submerged horizontal solid plate. By means of the matched eigenfunction expansion method, a linear analytical solution is developed for the interaction of water waves with the structure. Then the reflection coefficient, the transmission coefficient, the energy-loss coefficient and the wave forces acting on the plates are calculated. The numerical results obtained for limiting cases are exactly the same as previous predictions for a single submerged horizontal solid plate and a single submerged horizontal porous plate. Numerical results show that with a suitable geometrical porosity of the upper plate, the uplift wave forces on both plates can be controlled at a low level. Numerical results also show that the transmission coefficient will be always small if the dimensionless plate length (plate length versus incident wavelength) exceeds a certain moderate value. This is rather significant for practical engineering, as the incident wavelength varies over a wide range in practice. Moreover, it is found that the hydrodynamic performance of the present structure may be further enhanced if the lower plate is also perforated.     

Wave reflection and transmission by porous breakwaters: A new analytical solution

This study gives a new analytical solution for wave reflection and transmission by a surface-piercing porous breakwater. Velocity potential decompositions in the breakwater are used to obtain the solution. Different from traditional solutions, the new solution needs no complex dispersion relations (complex wave numbers) for wave motion through porous media. Thus, difficult procedures in traditional solutions for solving complex dispersion relations and handling non-self-adjoint eigenvalue problem are avoided. The calculated results of the new solution are in very good agreement with those of traditional solutions and multi-domain boundary element method solutions.     

深埋海底隧洞内水压力下饱和土的动力特性

将衬砌分别视为多孔柔性材料和弹性介质,在频率域内研究内水压力作用下饱和分数导数黏弹性土-深埋圆形隧洞多孔柔性或弹性衬砌系统的耦合简谐振动。土体的宏观力学特性采用多孔介质理论来模拟。通过引入与土体体积分数相关的应力系数,利用衬砌和土体界面处位移连续,分别得到饱和黏弹性土和衬砌的位移、应力和孔隙水压力等的解析表达式。在此基础上,分析了多孔柔性衬砌和弹性衬砌结构的差异,并考察了应力系数、渗透系数、分数阶导数本构参数等对系统动力响应的影响,结果表明:多孔柔性衬砌材料条件下系统的动力响应明显大于弹性衬砌材料条件下系统的动力响应;随着分数导数阶数和材料参数比的增加,系统共振效应明显减弱;衬砌边界透水和不透水只是反映边界渗透性的两种极限状态。