Propagation of a solitary wave over rigid porous beds

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.     

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.     

The numerical study of wave-induced pore water pressure response in highly permeable seabed

The coupling numerical model of wave interaction with porous medium is used to study waveinduced pore water pressure in high permeability seabed.In the model,the wave field solver is based on the two dimensional Reynolds-averaged Navier-Stokes(RANS) equations with a k-ε closure,and Forchheimer equations are adopted for flow within the porous media.By introducing a Velocity-Pressure Correction equation for the wave flow and porous flow,a highly efficient coupling between the two flows is implemented.The numerical tests are conducted to study the effects of seabed thickness,porosity,particle size and intrinsic permeability coefficient on regular wave and solitary wave-induced pore water pressure response.The results indicate that,as compared with regular wave-induced,solitary wave-induced pore water pressure has larger values and stronger action on seabed with different parameters.The results also clearly show the flow characteristics of pore water flow within seabed and water wave flow on seabed.The maximum pore water flow velocities within seabed under solitary wave action are higher than those under regular wave action.     

Verification of a VOF-based two-phase flow model for wave breaking and wave–structure interactions

The objective of the present study is to develop a volume of fluid (VOF)-based two-phase flow model and to discuss the applicability of the model to the simulation of wave–structure interactions. First, an overview of the development of VOF-type models for applications in the field of coastal engineering is presented. The numerical VOF-based two-phase flow model has been developed and applied to the simulations of wave interactions with a submerged breakwater as well as of wave breaking on a slope. Numerical results are then compared with laboratory experimental data in order to verify the applicability of the numerical model to the simulations of complex interactions of waves and permeable coastal structures, including the effects of wave breaking. It is concluded that the two-phase flow model with the aid of the advanced VOF technique can provide with acceptably accurate numerical results on the route to practical purposes.     

波浪与抛石潜堤相互作用的MPS法数值模型研究

为了研究波浪与抛石潜堤相互作用过程中大自由表面变形和堤内渗流等强非线性紊流运动问题,利用改进的MPS法,建立了模拟波浪与抛石潜堤相互作用的MPS法数值计算模型。模型将抛石潜堤假定为均质多孔介质,采用Drew的二相流运动方程描述多孔介质内外的流体运动;通过在动量方程中增加非线性阻力项,并引入亚粒子尺度紊流模型,模拟波浪与可渗结构物相互作用过程中的紊流运动。选取“U”型管中多孔介质内渗流过程和孤立波与可渗潜堤相互作用两个典型的渗流问题,通过将数值计算结果与理论解和实测值的对比分析,对所提出的MPS法紊流渗流模型的模拟精度进行验证。结果表明:基于改进的MPS法构建的垂向二维紊流渗流模型可以很好地再现“U”型管中多孔介质内渗流以及波浪作用下可渗潜堤内外的复杂流场,显著缓解流-固界面处的压力震荡与粒子分布不均匀问题,实现了较高的模拟精度。     

Numerical modeling of nonlinear water waves over heterogeneous porous beds

The transformation of nonlinear water waves over porous beds is studied by applying a numerical model based on Chen's [2006. Fully nonlinear Boussinesq-type equations for waves and currents over porous beds. Journal of Engineering Mechanics, 132:2, 220–230] Boussinesq-type equations for highly nonlinear waves on permeable beds. The numerical model uses a high-order time-marching solution and fourth-order finite-difference schemes for discretization of first-order spatial derivatives to obtain a computational accuracy consistent with the model equations. By forcing the wave celerity and spatial porous-damping rate of the linearized model to match the exact linear theory for horizontal porous bed over a prescribed range of relative depths, the values of the model parameters are optimally determined. Numerical simulations of the damped wave propagation over finite-thickness porous layer demonstrate the accuracy of both the numerical model and governing equations, which have been shown by prior theoretical analyses to be accurate for both nominal and thick porous layers. These simulations also elucidate on the significance of the higher-order porous-damping terms and the influence of the hydraulic parameters. Application of the model to the simulation of the wave field around a laboratory-scale submerged porous mound provides a measure of its capability, as well as useful insight into the scaling of the porous-resistance coefficients. For application to heterogeneous porous beds, the assumption of weak spatial variation of the porous resistance is examined using truncated forms of the governing equations. The results indicate that the complete set of Boussinesq-type equations is applicable to porous beds of nonhomogeneous makeup.     

A coupled CFD-DEM simulation of upward seepage flow in coarse sands

Abstract

The suction anchor becomes more popular for offshore oil and gas industry in deeper water. For suction anchor–soil interaction, the prediction of hydraulic conductivity of porous materials is a long-standing problem in offshore engineering. To investigate the hydraulic characteristics, an upward seepage flow through saturated sands is considered in this study. A numerical approach, which is able to describe the fluid–particle interaction at particle scale, has been employed to analyse fluid flow in sands. This approach is constructed by adopting a coupled discrete element method and computational fluid dynamic approach (CFD-DEM numerical model). The coupled CFD-DEM approach is first benchmarked by a classic geomechanics problem where analytical solutions are available, and then employed to investigate the characteristics of upward seepage flow in coarse sand columns. Through numerical modelling, the predicted relation between hydraulic gradient and flow velocity is obtained and it is compared with the classical analytical correlation. The effect of several bulk and micromechanical parameters including packing porosity, particle size combination and inter-particle rolling resistance on the flow characteristics is numerically examined. The results show that the particle polydispersity and packing porosity have significant effect on the hydraulic conductivity in the seepage flow. The introduction of inter-particle rolling resistance can change initial packing structure of particle assembly in some extent rather than the hydraulic conductivity from the particle shape effect perspective. A further development of numerical model, in which the effect of non-spherical particles on the seepage flow, will be carried out later.     

2-D numerical analysis of near-field flow at low-crested permeable breakwaters

This paper describes the capability of a numerical model named COrnell BReaking waves And Structures (COBRAS) [Lin, P., Liu, P.L.-F., 1998. A numerical study of breaking waves in the surf zone. Journal of Fluid Mechanics 359, 239–264; Liu, P.L.-F., Lin, P., Chang, K.A., Sakakiyama, T., 1999. Numerical modeling of wave interaction with porous structures. Journal of Waterway, Port, Coastal and Ocean Engineering 125, 322–330, Liu, P.L.-F., Lin, P., Hsu, T., Chang, K., Losada, I.J., Vidal, C., Sakakiyama, T., 2000. A Reynolds averaged Navier–Stokes equation model for nonlinear water wave and structure interactions. Proc. Coastal Structures '99, 169–174] based on the Reynolds Averaged Navier–Stokes (RANS) equations to simulate the most relevant hydrodynamic near-field processes that take place in the interaction between waves and low-crested breakwaters. The model considers wave reflection, transmission, overtopping and breaking due to transient nonlinear waves including turbulence in the fluid domain and in the permeable regions for any kind of geometry and number of layers. Small-scale laboratory tests were conducted in order to validate the model, with different wave conditions and breakwater configurations. In the present study, regular waves of different heights and periods impinging on a wide-crested structure are considered. Three different water depths were tested in order to examine the influence of the structure freeboard. The experimental set-up includes a flow recirculation system aimed at preventing water piling-up at the lee of the breakwater due to overtopping. The applicability and validity of the model are examined by comparing the results of the numerical computations with experimental data. The model is proved to simulate with a high degree of agreement all the studied magnitudes, free surface displacement, pressure inside the porous structure and velocity field. The results obtained show that this model represents a substantial improvement in the numerical modelling of low-crested structures (LCS) since it includes many processes neglected previously by existing models. The information provided by the model can be useful to analyse structure functionality, structure stability, scour and many other hydrodynamic processes of interest.     

Numerical Simulation on Hydraulic Performances of Quarter Circular Breakwater

Quarter circular breakwater （QCB） is a new-type breakwater developed from senti-circular breakwater （SCB）. The superstructure of QCB is composed of a quarter circular front wall, a horizontal base slab and a vertical rear wall. The width of QCB＇ s base slab is about half that of SCB, which makes QCB suitable to be used on relatively finn soil foundation. The numerical wave flume based on the Reynolds averaged Navier-Stokes equations for impressible viscosity fluid is adopted in this paper to simulate the hydraulic performances of QCB. Since the geometry of both breakwaters is similar and SCB has been studied in depth, the hydraulic performances of QCB are given in comparison with those of SCB.     

Incompressible SPH flow model for wave interactions with porous media

The paper presents an Incompressible Smoothed Particle Hydrodynamics (ISPH) method to simulate wave interactions with a porous medium. The SPH method is a mesh free particle modeling approach that is capable of tracking the large deformation of free surfaces in an easy and accurate manner. The ISPH method employs a strict incompressible hydrodynamic formulation to solve the fluid pressure and the numerical solution is obtained by using a two-step semi-implicit scheme. The ISPH flow model solves the unsteady 2D Navier–Stokes (NS) equations for the flows outside the porous media and the NS type model equations for the flows inside the porous media. The presence of porous media is considered by including additional friction forces into the equations. The developed ISPH model is first validated by the solitary and regular waves damping over a porous bed and the solitary wave interacting with a submerged porous breakwater. The convergence of the method and the sensitivity of relevant model parameters are discussed. Then the model is applied to the breaking wave interacting with a breakwater covered with a layer of porous materials. The computational results demonstrate that the ISPH flow model could provide a promising simulation tool in coastal hydrodynamic applications.     

Numerical simulation of wave damping over porous seabeds

This paper presents a study of wave damping over porous seabeds by using a two-dimensional numerical model. In this model, the flow outside of porous media is described by the Reynolds Averaged Navier–Stokes equations. The spatially averaged Navier–Stokes equations, in which the presence of porous media is considered by including additional inertia and nonlinear friction forces, is derived and implemented for the porous flow. Unlike the earlier models, the present model explicitly represents the flow resistance dependency on Reynolds number in order to cover wider ranges of porous flows. The numerical model is validated against available theories and experimental data. The comparison between the numerical results and the theoretical results indicates that the omission or linearization of the nonlinear resistance terms in porous flow models, which is the common practice in most of analytical models, can lead to significant errors in estimating wave damping rate. The present numerical model is used to simulate nonlinear wave interaction with porous seabeds and it is found that the numerical results compare well with the experimental data for different wave nonlinearity. The additional numerical tests are also conducted to study the effects of wavelength, seabed thickness and Reynolds number on wave damping.     

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.     

Numerical study of air chamber for oscillating water column wave energy convertor

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.     

Fully Coupled Fluid-Structure Interaction Model Based on Distributed Lagrange Multiplier/Fictitious Domain Method

This paper,with a finite element method,studies the interaction of a coupled incompressible fluid-rigid structure system with a free surface subjected to external wave excitations.With this fully coupled model,the rigid structure is taken as "fictitious" fluid with zero strain rate.Both fluid and structure are described by velocity and pressure.The whole domain,including fluid region and structure region,is modeled by the incompressible Navier-Stokes equations which are discretized with fixed Eulerian mesh.However,to keep the structure's rigid body shape and behavior,a rigid body constraint is enforced on the "fictitious" fluid domain by use of the Distributed Lagrange Multiplier/Fictitious Domain(DLM/FD) method which is originally introduced to solve particulate flow problems by Glowinski et al.For the verification of the model presented herein,a 2D numerical wave tank is established to simulate small amplitude wave propagations,and then numerical results are compared with analytical solutions.Finally,a 2D example of fluid-structure interaction under wave dynamic forces provides convincing evidences for the method excellent solution quality and fidelity.     

Water waves within a porous medium on an undulating bed

Water wave refraction–diffraction within a porous medium on an undulating seabed is considered based on linear wave theory. Using the model of wave-induced flow within a porous medium and Galerkin eigenfunction expansions, refraction–diffraction equations for surface waves are derived. With these equations, the wave reflection from a porous structure on a sloping beach is investigated and numerical results of reflection coefficients are obtained. A comparison between the present results with those in the literature is made for a special case and the agreement is satisfactory. This structure can be viewed as an idealized model of rubble-mound seawalls along coastlines.     

Propagation of water waves over permeable rippled beds

Unsteady two-dimensional Navier-Stokes equations and Navier-Stokes type model equations for porous flow were solved numerically to simulate the propagation of water waves over a permeable rippled bed. A boundary-fitted coordinate system was adopted to make the computational meshes consistent with the rippled bed. The accuracy of the numerical scheme was confirmed by comparing the numerical results concerning the spatial distribution of wave amplitudes over impermeable and permeable rippled beds with the analytical solutions. For periodic incident waves, the flow field over the wavy wall is discussed in terms of the steady Eulerian streaming velocity. The trajectories of the fluid particles that are initially located close to the ripples were also determined. One of the main results herein is that under the action of periodic water waves, fluid particles on an impermeable rippled bed initially moved back and forth around the ripple crest, with increasing vertical distance from the rippled wall. After one or two wave periods, they are then lifted towards the next ripple crest. All of the marked particles on a permeable rippled bed were shifted onshore with a much larger displacement than those on an impermeable bed. Finally, the flow fields and the particle motions close to impermeable and permeable beds induced by a solitary wave are elucidated.     

波浪与大孔隙多孔介质相互作用的耦合数学模型

建立了波浪与大孔隙多孔介质相互作用的耦合数学模型,波浪域的控制方程为雷诺时均方程和k-ε紊流模型。对于计算域的入射波采用推板式造波,它可以是线性波、椭圆余弦波和孤立波。采用PLIC-VOF法追踪波浪自由表面。对于多孔介质内的孔隙流场采用非线性Forchheimer方程,两区域共享连续方程,最后导出的波浪域与孔隙流域的压力修正方程具有完全相同的形式,利用这个方程能够同时而不是分别求解波浪场和孔隙流场,避免了在内部边界上给定匹配条件,实现了波浪场与孔隙流场的同步耦合。波浪与粗颗粒海床、平底床面上抛石潜堤及斜坡上抛石潜堤相互作用的验证计算结果表明该模型可用于研究波浪与大孔隙多孔介质相互作用的问题。     

Modelling of multi-bodies in close proximity under water waves—Fluid forces on floating bodies

This work presents two-dimensional numerical results of the dependence of wave forces of multiple floating bodies in close proximity on the incident wave frequency, gap width, body draft, body breadth and body number based on both viscous fluid and potential flow models. The numerical models were validated by the available experimental data of fluid oscillation in narrow gaps. Numerical investigations show that the large amplitude responses of horizontal and vertical wave forces appear around the fluid resonant frequencies. The convectional potential flow model is observed to un-physically overestimate the magnitudes of wave forces as the fluid resonance takes place. By introducing artificial damping term with appropriate damping coefficients μ∈[0.4, 0.5], the potential flow model may work as well as the viscous fluid model, which agree with the damping coefficients used in our previous work for the predication of wave height under gap resonance. In addition, the numerical results of viscous fluid model suggest that the horizontal wave force is highly dependent on the water level difference between the opposite sides of an individual body and the overall horizontal wave force on the floating system is generally smaller than the summation of wave force on each body.     

Numerical simulation of random wave slamming on structures in the splash zone

The numerical investigation of random wave slamming on superstructures of marine structures in the splash zone is presented in this paper. The impact pressures on the underside of the structure are computed based on the improved volume of fluid method (VOF). The governing equations are Reynolds time-averaged equations and the two equation k– model. The third order upwind difference scheme is applied to the convection term to reduce the effect of numerical viscosity. The numerical wave flume with random wave-maker suitable for VOF is established. Appropriate moving contact-line boundary conditions are introduced to the model wave in contact with and separated from the underside of structure. Parametric studies have been carried out for different incident waves, structure dimensions and structure clearance. The numerical results are verified by the experimental results.     

Simulation of the nonlinear dynamic interactions between waves, a submerged breakwater and the seabed

This study employed direct numerical simulation to simulate the fully nonlinear interaction between the water waves, the submerged breakwater, and the seabed under differing wave conditions. In the numerical simulation, the laminar flow condition in the seabed was applied to evaluate the more exact fluid resistance acting on the porous media. Varying incident wave conditions were applied to the flow field resulting from the wave–structure–seabed interaction, and the variation in the pore water pressure beneath the submerged breakwater was investigated along the cross-section of the submerged breakwater. Structural safety and scouring were also considered on the basis of the numerical results for the flow field around the structure and the variation of the pore water pressure.