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.     

Wave-induced response of seabed: Various formulations and their applicability

In this study, a set of generalized analytical solutions are developed for the wave-induced response of a saturated porous seabed under plane strain condition. When considering the water waves originating in deep water and travelling towards the shore, their velocities, lengths and heights vary. Depending on the characteristics of the wave and the properties of the seabed, different formulations (fully dynamic, partly dynamic, quasi-static) for the wave-induced response of the seabed are possible. The solutions for the response with these formulations are established in terms of non-dimensional parameters. The results are presented in terms of pore pressure, shear stress and vertical effective stress distributions within the seabed. For typical values of wave period and seabed permeability, the regions of applicability of the three formulations are identified and plotted in parametric spaces. With given wave and seabed characteristics, these regions provide quick identification of the appropriate formulation for an adequate evaluation of the wave-induced seabed response.     

不同破碎波对沙质海床作用的实验研究

破碎波对近海海岸地形以及海岸建筑物影响强烈,通过物理模型实验对孤立波、规则波作用下破碎带的床面形态以及孔隙水压力进行分析。破碎波冲击海床,破碎处床面上形成沙坝和沙坑,与规则波相比,孤立波破碎时对床面的冲刷更加剧烈,床面形成的沙坝和沙坑尺度更大,且土体内孔隙水压力幅值也较大。同时研究了波面变化对孔隙水压力的影响,发现波面变化历时曲线与孔隙水压力历时曲线相似,与孔隙水压力梯度历时曲线更为相似,说明波面变化更能反映海床内部孔隙水压力梯度的变化。通过探讨波浪与海床之间相互耦合作用,发现破碎带地形变化使得波浪出现不同破碎类型,分析得出卷破波比崩破波作用下孔隙水压力幅值大。     

Numerical study on the interaction among a nonlinear wave,composite breakwater and sandy seabed

Most previous investigations related to composite breakwaters have focused on the wave forces acting on the structure itself from a hydrodynamic aspect. The foundational aspects of a composite breakwater under wave-induced cyclic loading are also important in studying the stability of a composite breakwater. In this study, numerical simulations were performed to investigate the wave-induced pore water pressure and flow changes inside the rubble mound of the composite breakwater and seabed foundation. The validity and applicability of the numerical model were demonstrated by comparing numerical results with existing experimental data. Moreover, the present model clearly has shown that the instantaneous directions of pore water flow motion inside the seabed induced by surface waves are in good agreement with the general wave-induced pore water flow inside the seabed. The model is further used to discuss the stability of a composite breakwater, i.e., the interaction among nonlinear waves, composite breakwater and seabed. Numerical results suggest that the stability of a composite breakwater is affected by not only downward shear flow generating on the seaward slope face of the rubble mound but, also, a high and dense pore water pressure gradient inside the rubble mound and seabed foundation.     

Stability and liquefaction analysis of porous seabed subjected to cnoidal wave

Cnoidal wave theory is appropriate to periodic wave progressing in water whose depth is less than 1/10 wavelength. However, the cnoidal wave theory has not been widely applied in practical engineering because the formula for wave profile involves Jacobian elliptic function. In this paper, a cnoidal wave-seabed system is modeled and discussed in detail. The seabed is treated as porous medium and characterized by Biot's partly dynamic equations (u–p model). A simple and useful calculating technique for Jacobian elliptic function is presented. Upon specification of water depth, wave height and wave period, Taylor's expression and precise integration method are used to estimate Jacobian elliptic function and cnoidal wave pressure. Based on the numerical results, the effects of cnoidal wave and seabed characteristics, such as water depth, wave height, wave period, permeability, elastic modulus, and degree of saturation, on the cnoidal wave-induced excess pore pressure and liquefaction phenomenon are studied.     

Wave-induced seabed instability in front of a breakwater

The wave-induced soil response in a porous seabed has become an important factor for the stability of offshore facilities, because many marine structures may have failed due to seabed instability and concomitant subsidence. An analytical solution is presented for the wave-induced soil response under the action of a three-dimensional wave system. Based on this general solution, the mechanism of seabed instability is then investigated. The general solutions for pore pressure and effective stresses are readily reducible to two dimensions for progressive waves, and are compared to theoretical and experimental work available. Some dominant factors affecting the wave-induced seabed instability are discussed; including permeability, seabed thickness and degree of saturation.     

波浪作用下孔隙海床-管线动力相互作用分析

波浪作用下海床中的孔隙水压力与有效应力是影响海底管线稳定性的主要因素。然而,在目前的海床响应分析中一般将管线假定为刚性,并不能合理地考虑海床与管线之间的相互作用效应,同时也没有考虑土体和管线加速度对海床动力响应的惯性影响,从而无法确定由此所引起的管线内应力。为此考虑管线的柔性,分别采用饱和孔隙介质的Biot动力固结理论和弹性动力学理论列出了海床与管线的控制方程,进而采用摩擦接触理论考虑海床与管线之间的相互作用效应,基于有限元方法建立了海床-管线相互作用的计算模型及其数值算法。通过变动参数对比计算讨论了管线几何尺寸、海床土性参数对波浪所引起的管线周围海床孔隙水压力和管线内应力的影响。     

Transient soil response in a porous seabed with variable permeability

The soil permeability of many natural marine sediments decreases with depth because of consolidation under overburden pressure. This is accompanied by a decrease in porosity and void ratio that also affect the permeability. Conventional theories for wave-induced soil response have assumed a homogeneous porous seabed. This paper presents a new approach for the wave-induced response in a soil matrix, with variable permeability as a function of burial depth. The soil matrix considered is unsaturated and anisotropic, and is subject to a three-dimensional wave system. The pore pressure and effective stresses induced by such a system are obtained from a set of equations incorporating a variable permeability. Verification is available through reduction to the simple case of uniform permeability. The results indicate that the effect of variable soil permeability on pore pressure and vertical effective stress may be significant, especially in a gravelled seabed and for unsaturated sandy soils.     

Parametric study of the wave-induced residual liquefaction around an embedded pipeline

Wave-induced liquefaction in a porous seabed around submarine pipeline may cause catastrophic consequences such as large horizontal displacements of pipelines on the seabed, sinking or floatation of buried pipelines. Most previous studies in relation to the wave and seabed interactions with embedded pipeline dealt with the wave-induced instaneous seabed response and possible resulting momentary liquefaction (where the soil is liquefied instantaneously during the passage of a wave trough), using theory of poro-elasticity. Studies for the interactions between a buried pipeline and a soil undergoing build-up of pore pressure and residual liquefaction have been comparatively rare. In this paper, this complicated process was investigated by using a new developed integrated numerical model with RANS (Reynolds averaged Navier–Stokes) equations used for governing the incompressible flow in the wave field and Biot consolidation equations used for linking the solid–pore fluid interactions in a porous seabed with embedded pipeline. Regarding the wave-induced residual soil response, a two-dimensional poro-elastoplastic solution with the new definition of the source term was developed, where the pre-consolidation analysis of seabed foundation under gravitational forces including the body forces of a pipeline was incorporated. The proposed numerical model was verified with laboratory experiment to demonstrate its accuracy and effectiveness. The numerical results indicate that residual liquefaction is more likely to occur in the vicinity of the pipeline compared to that in the far-field. The inclusion of body forces of a pipeline in the pre-consolidation analysis of seabed foundation significantly affects the potential for residual liquefaction in the vicinity of the pipeline, especially for a shallow-embedded case. Parametric studies reveal that the gradients of maximum liquefaction depth with various wave and soil characteristics become steeper as pipeline burial depth decreases.     

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

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

Evaluation of storm wave-induced silty seabed instability and geo-hazards: A case study in the Yellow River delta

Models based on the theoretical framework of soil mechanics are presented to evaluate storm wave-induced silty seabed instability and geo-hazards through a case study in the Yellow River delta. First, the transient and residual mechanisms of wave-induced pore pressure are analyzed. Three typical models (i.e., elastic model, pore pressure development mode and elasto-plastic model) are proposed to calculate wave-induced stresses in the seabed. Next, mechanisms and calculation methods of wave-induced seabed instability modes such as scour, liquefaction, seepage instability and shear slide are proposed. Typical results of storm wave-induced excess pore pressure and seabed instability are given and relevant discussions are made. At last, the formation mechanism of geo-hazards in the Yellow River delta is analyzed based on the proposed mechanism and calculated results. Results and analysis indicate that both transient and residual mechanisms are important to storm wave-induced response of silty seabed and hence the elasto-plastic model is more appropriate. Complete liquefaction does not happen, while other types of instability occur mostly within 2–6 m under the seabed surface. Wave-induced scour, seepage instability and shear slide are all possible instability modes under the 1-year storm waves, and scour is predominant for the 50-year storm waves. The formation mechanism of geo-hazards such as shallow slide and storm wave reactivation, pockmarks, silt flow and gully, disturbed stratum and hard crust in the Yellow River are well explained based on the proposed mechanisms and calculated results of storm wave-induced silty seabed instability.     

波浪作用下斜坡上单桩周围海床孔隙水压力响应特性试验研究

单桩基础周围斜坡海床中的波致孔隙水压力响应与纯斜坡海床存在较大差异。为了解不同波高、波周期条件下,单桩基础周围波浪传播变形及其对斜坡海床孔压振荡响应的影响,在波浪水槽末端铺设了长6 m、坡度1∶16的斜坡砂床进行试验。通过改变桩身位置和波浪参数,测量斜坡段各处波面形态,采集单桩周围孔隙水压力,分析了桩身位置及波浪参数对斜坡海床孔压响应的影响。结果表明：相同入射波条件下,随距坡脚水平距离增加,波高、近底流速和桩周孔隙水压力幅值都随之增大;桩周孔隙水压力幅值分布规律为：桩前孔压幅值明显大于桩侧与桩后孔压幅值。当Keulegan-Carpenter数大于6时,随着波高和波周期增大,马蹄涡产生的负压区使得桩侧海床孔隙水压力与纯斜坡海床孔隙水压力差值迅速增加。     

Pore pressure response of seabed in standing waves and its mechanism

Experiments were carried out in a newly-developed rectangular flume to investigate the pore pressure response of silty seabed in standing waves, and to emphatically discuss the physical mechanism of wave-induced pore pressure. We first analyzed the features of wave-induced stress, laid out the experimental setup, and particularly monitored the temporal and spatial variation of pore pressure. Then, we summarized the physical essence of wave-induced pore pressure, discussed the applicability of the existing models, and proposed several key points about developing new models. Results indicate that regular standing waves were obtained in the rectangular flume. The oscillating pore pressure (OPP) and the residual pore pressure (RPP) were observed simultaneously in the seabed under wave nodes or antinodes, and were found with coupling effect. The OPP and RPP were related to the elastic and plastic volumetric strains, respectively. Both elastic and plastic volumetric strains may be caused by the wave-induced cyclic spherical stress or deviatoric stress. The elastic model could only simulate the OPP and the pore pressure development mode could only simulate the RPP. The yielding under isotropic compression and the shear expansion of geotechnical materials, as well as the wave-induced stress paths, should be considered in the elastoplastic model to simulate the wave-induced pore pressure accurately.     

Numerical simulation of the nonlinear wave-induced dynamic response of anisotropic poro-elastoplastic seabed

Abstract

In this paper, a 2D poro-elastoplastic model for wave-induced dynamic response in an anisotropic seabed is derived analytically. The seabed is treated as a porous medium and characterized by Biot’s consolidation equations. The soil plasticity and wave non-linearity are included in the model and both the pore fluid and the soil skeleton are assumed to be compressible. The nonlinear ocean waves are respectively considered as progressive and standing waves. The previous experimental data is used to validate the proposed model. Numerical results demonstrate that the influence of nonlinear wave components should not be ignored without committing substantial error. A significant difference between progressive and standing waves is also observed for the development of residual pore pressure, as well as the distribution of liquefied zone. A detailed parametric investigation reveals that the nonlinear wave-induced seabed response is also affected significantly by cross-anisotropic soil parameters.     

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.     

An analytical solution for wave-induced seabed response in a multi-layered poro-elastic seabed

In this study, a new analytical solution for the wave-induced seabed response in a multi-layered poro-elastic seabed is developed. The seabed is treated as a multi-layered porous medium and characterized by Biot’s theory. The displacements of the solid skeleton and the pore pressure are expressed in terms of two scalar potentials and one vector. Then, the Biot’s dynamic equation can be solved using Fourier transformation and reducing to Helmholtz equations. To obtain the general solutions for the multi-layered poro-elastic seabed in the frequency-wave-number domain, the transmission and reflection matrices (TRM) method is used to form the equivalent stiffness. Using the boundary conditions and continuous conditions, the frequency-wave-number domain solutions are obtained. Finally, the time-space domain solutions for the multi-layered poro-elastic seabed are obtained by means of the inverse Fourier transformation with respect to the horizontal coordinate. Based on the new solution, a parametric study is carried out to examine the effects of soil characteristics (number of layers, permeability and shear modulus) and wave characteristics (water depth and wave steepness) on seabed responses. The results indicate that the seabed response is affected significantly by permeability, shear modulus and relative water depth.     

Effects of cross-correlated multiple spatially random soil properties on wave-induced oscillatory seabed response

The evaluation of seabed response under wave loading is important for prediction of stability of foundations of offshore structures. In this study, a stochastic finite element model which integrates the Karhunen-Loève expansion random field simulation and finite element modeling of wave-induced seabed response is established. The wave-induced oscillatory response in a spatially random heterogeneous porous seabed considering cross-correlated multiple soil properties is investigated. The effects of multiple spatial random soil properties, correlation length and the trend function (the relation of the mean value versus depth) on oscillatory pore water pressure and momentary liquefaction are discussed. The stochastic analyses show that the uncertainty bounds of oscillatory pore water pressure are wider for the case with multiple spatially random soil properties compared with those with the single random soil property. The mean pore water pressure of the stochastic analysis is greater than the one obtained by the deterministic analysis. Therefore, the average momentary liquefaction zone in the stochastic analysis is shallower than the deterministic one. The median of momentary liquefaction depth generally decreases with the increase of vertical correlation length. When the slope of the trend function increases, the uncertainty of pore water pressure is greatly reduced at deeper depth of the seabed. Without considering the trend of soil properties, the wave-induced momentary liquefaction potential may be underestimated.     

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.     

Experimental investigation on the wave-induced pore pressure around shallowly embedded pipelines

A series of regular wave experiments have been done in a large-scale wave flume to investigate the wave-induced pore pressure around the submarine shallowly embedded pipelines.The model pipelines are buried in three kinds of soils,including gravel,sand and silt with different burial depth.The input waves change with height and period.The results show that the amplitudes of wave-induced pore pressure increase as the wave period increase,and decay from the surface to the bottom of seabed.Higher pore pressures are recorded at the pipeline top and the lower pore pressures at the bottom,especially in the sand seabed.The normalized pressure around pipeline decreases as the relative water depth,burial depth or scattering parameters increase.For the silt seabed,the wavelet transform has been successfully used to analyze the signals of wave-induced pore pressure,and the oscillatory and residual pore pressure can be extracted by wavelet analysis.Higher oscillatory pressures are recorded at the bottom and the lower pressures at the top of the pipeline.However,higher residual pressures are recorded at the top and the lower pressures at the bottom of the pipeline.