首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 328 毫秒
1.
Incompressible SPH flow model for wave interactions with porous media   总被引:1,自引:0,他引:1  
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.  相似文献   

2.
内孤立波中半潜平台动力响应特性   总被引:1,自引:1,他引:0  
基于内孤立波mKdV理论,采用时域有限位移运动方程,结合改进的Morison公式,研究了两层流体中内孤立波与带分段式悬链系泊约束半潜平台的相互作用问题。针对东沙群岛南部海域附近某实测内孤立波特征参数,计算分析了在该内孤立波作用下SEDCO-700型半潜平台的内孤立波载荷、运动响应及其系泊张力的变化特性。研究表明,内孤立波不仅会对半潜平台产生突发性冲击载荷作用,使其产生大幅度水平漂移运动,并导致其系泊张力显著增大,因此在半潜平台等深海平台的设计与应用中,内孤立波的影响是不可忽视的。  相似文献   

3.
Propagation of a solitary wave over rigid porous beds   总被引:1,自引:0,他引:1  
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.  相似文献   

4.
为了研究泻湖型海湾内经常出现的湾中岛的形成机理,应用水平二维潮流、泥沙输移和地形演变耦合模型,对具有典型的沙坝泻湖地貌形态的水东湾的湾中岛的形成和演化进行了数值模拟,成功地模拟出了涨潮三角洲(大洲岛)的形成过程,计算结果与实际地貌形态总体符合。模型中分别考虑了全沙输沙和推移质输沙两种输沙情况。结果表明,水东湾湾中岛是由涨潮流引起的泻湖内泥沙不断淤积而形成的,是一种涨潮三角洲的地貌形态;全沙输沙模式比推移质输沙模式更适合这一地形演化过程的模拟。模拟结果也再现了湾中岛在25 a期间的地形变化过程,这一研究结果为通过计算数值模拟来研究一般海湾的地貌形态的形成机理和演化过程提供了实际算例。  相似文献   

5.
This study aims to develop a robust, accurate and computationally efficient hydrodynamic and sediment transport model for dam break flows. The two dimensional shallow water equations are resolved based on the finite volume method with an unstructured quadtree mesh. The sediment transport and bed evolution modules are coupled with hydrodynamic module to predict simultaneously the hydrodynamics, sediment concentrations and morphological changes. The interface flux is computed by the HLL approximate Riemann solver with second order accuracy. The effects of pressure and gravity are included in source term in this model, which can simplify the computation and eliminate numerical imbalance between source and flux terms. For dam break flows occurring in complicated geometries, the quadtree rectangular mesh is used to refine the interesting area and important part. The model is first verified against results from laboratory experiments, existing numerical models and real life case. It is then used to simulate dam break flows over a mobile bed to investigate the bed evolution. The results are compared with experimental data and field data with good agreement. The method is simple, efficient, and conservative. It shows promise for handling hydrodynamic simulation and sediment transport for a wide range of dam break flows.  相似文献   

6.
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 sech2(*) 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.  相似文献   

7.
LU  Yongjun 《中国海洋工程》2002,16(1):107-122
A 2-D mathematical model of tidal current and sediment has been developed for the Oujiang Estuary and the WenzhouBay. This model accomodates complicated features including multiple islands, existence of turbidity, and significant differ-ence in size distribution of bed material. The governing equations for non-uniform suspended load and bed load transport arepresented in a boundary-fitted orthogonal curvilinear coordinate system. The numerical solution procedures along with theirinitial conditions, boundary conditions, and movable boundary technique are presented. Strategies for computation of thecritical condition of deposition or erosion, sediment transport capacity, non-uniform bed load discharge, etc. are suggested.The model verification computation shows that, the tidal levels computed from the model are in good agreement with the fielddata at the 18 tidal gauge stations. The computed velocities and flow directions also agree well with the values measuredalong the totally 52 synchronously observed verticals distributed over 8 cross sections. The computed tidal water throughputsthrough the Huangda‘ao cross section are close to the measured data. And the computed values of bed deformation fromYangfushan to the estuary outfall and in the outer-sea area are in good agreement with the data observed from 1986 to 1992.The changes of tidal volumes through the estuary, velocities in different channels and the bed form due to the influence of thereclamation project on the Wenzhou shoal are predicted by means of this model.  相似文献   

8.
A 2DH numerical, model which is capable of computing nearshore circulation and morphodynamics, including dune erosion, breaching and overwash, is used to simulate overwash caused by Hurricane Ivan (2004) on a barrier island. The model is forced using parametric wave and surge time series based on field data and large-scale numerical model results. The model predicted beach face and dune erosion reasonably well as well as the development of washover fans. Furthermore, the model demonstrated considerable quantitative skill (upwards of 66% of variance explained, maximum bias − 0.21 m) in hindcasting the post-storm shape and elevation of the subaerial barrier island when a sheet flow sediment transport limiter was applied. The prediction skill ranged between 0.66 and 0.77 in a series of sensitivity tests in which several hydraulic forcing parameters were varied. The sensitivity studies showed that the variations in the incident wave height and wave period affected the entire simulated island morphology while variations in the surge level gradient between the ocean and back barrier bay affected the amount of deposition on the back barrier and in the back barrier bay. The model sensitivity to the sheet flow sediment transport limiter, which served as a proxy for unknown factors controlling the resistance to erosion, was significantly greater than the sensitivity to the hydraulic forcing parameters. If no limiter was applied the simulated morphological response of the barrier island was an order of magnitude greater than the measured morphological response.  相似文献   

9.
The paper focuses on the numerical simulation of erosion of plane sloping beaches by irregular wave attack in three wave flumes of different scales. One of the prime objectives of the tests was to provide a consistent data set for the improvement of numerical beach profile models. A practical application of this research with wave attack on plane sloping beaches is the erosion of the plane beaches after nourishment. Three models (CROSMOR, UNIBEST-TC and DELFT3D) have been used to simulate the flume experimental results focusing on the wave height distribution and the morphological development (erosion and deposition) along the beach profiles. Overall, the model predictions for wave heights show consistent results. Generally, the computed wave heights (Hrms and H1/3) are within 10% to 15% of the measured values for all tests (under-prediction of the largest wave heights close to the shore). The three models can simulate the beach erosion of the wave flume tests (erosive tests) reasonably well using default values of the sand transport parameters. The model performance for the accretive tests is less good than that for the erosive tests. A practical field application of this research is the erosion of nourished beaches, as these beaches generally have rather plane beach slopes immediately after nourishment. Various graphs are given to estimate the beach erosion of nourished beaches.  相似文献   

10.
Wave-induced currents may drive nearshore mixing and transport processes, including coastal pollutant dispersion, littoral drift, and long-term morphological changes through beach erosion and accretion. In this study, a numerical model is newly developed to simulate wave climate and localized currents in complicated coastal environments. The model developed is based on a quadtree grid system. The two-dimensional hydrodynamic governing equations are solved by using an explicit Adams-Bashforth finite difference scheme. Effects of wave breaking, shoaling, refraction, diffraction, wave-current interaction, set-up and set-down, turbulent mixing, bed friction, and shoreline movement are incorporated in the model. Results are presented for set-up at a beach in a flume due to normally incident waves, and longshore currents generated by oblique waves on a plane beach.  相似文献   

11.
- In this paper, the solitary wave deformation along a gentle slope and the impact pressure, on the wall are investigated experimentally and the results are compared with numerical results obtained based on the volume of fluid (VOF) method. The topography used in the experiment consists of three segments. The left segment is a 1:4 slope, the middle segment 1 :SO slope and the right segment a horizontal bed. Both the wave heights and breaking points obtained from numerical simulation and experiments are in good agreement. Numerical results give reasonable pressure distributions of breaking waves on the wall.  相似文献   

12.
《Coastal Engineering》2006,53(4):363-379
Low crested structures for beach defence purposes are frequently configured in the form of cells. This paper investigates the hydrodynamic and morphological response of such defence schemes by means of laboratory experiments and numerical simulations.Experiments were carried out at the large mobile-bed basin of the “Laboratorio di Idraulica Costiera” of the “Politecnico di Bari”, IT. Tests peculiarity is the simultaneous simulation of three different types of lateral confinement: a channel, representative of an indefinitely long structure, and two areas with narrow and wide gaps, representative of a common scheme and of an almost isolated structure respectively. Submerged, zero freeboard and emerged conditions were tested. A typical storm of the Northern Adriatic sea was reproduced by 6 wave attacks.Water level and currents were measured in front and behind the structures, at gaps and roundheads. Results, focusing on rip current intensities and piling-up, are presented and discussed. An extensive plan view of wave intensities and flow patterns inside the basin is provided through numerical simulations with MIKE 21 numerical suite.The sand bed was surveyed along 30 profiles covering the narrow gap area and the isolated structure roundheads. For the submerged conditions tests, 15 of these profiles and the shoreline position were monitored after each step of the simulated storm.Average and maximum values of local erosion at the structure toe are presented. Time evolution of erosion at gaps is related to the main hydrodynamic variables through Bijker sediment transport formula. Evaluation of volume changes in the protected area shows that sand is in average trapped behind the barriers, being transported offshore from gaps and returning inshore over the crest. The observed and reconstructed transport mechanism is a combination of bed load and suspended transport, the latter being slightly more important.Bed level changes are well predicted by morphodynamic simulations performed with MIKE 21 CAMS.  相似文献   

13.
A 2-D mathematical model of tidal current and sediment has been developed for the Oujiang Estuary and the Wenzhou Bay. This model accomodates complicated features including multiple islands, existence of turbidity, and significant differ-ence in size distribution of bed material. The governing equations for non-uniform suspended load and bed load transport are presented in a boundary-fitted orthogonal curvilinear coordinate system. The numerical solution procedures along with their initial conditions, boundary conditions, and movable boundary technique are presented. Strategies for computation of the critical condition of deposition or erosion, sediment transport capacity, non-uniform bed load discharge, etc. are suggested. The model verification computation shows that, the tidal levels computed from the model are in good agreement with the field data at the 18 tidal gauge stations. The computed velocities and flow directions also agree well with the values measured along the totally 52 synchronously ob  相似文献   

14.
In this study an Euler-Euler two-phase model was developed to investigate the tunnel erosion beneath a submarine pipeline exposed to unidirectional flow. Both of the fluid and sediment phases were described via the Navier-Stokes equations, i.e. the model was implemented using time-averaged continuity and momentum equations for the fluid and sediment phases and a modified kε turbulence closure for the fluid phase. The fluid and sediment phases were coupled by considering the drag and lift interaction forces. The model was employed to simulate the tunnel erosion around the pipeline laid on an erodible bed. Comparison between the numerical result and experimental measurement confirms that the numerical model successfully predicts the bed profile and velocity field during the tunnel erosion. It is evident that the sediments are transported as the sheet-flow mode in the tunnel erosion stage. Also the transport rate under the pipe increases rapidly at the early stage and then reduces gradually at the end of the tunnel erosion beneath pipelines.  相似文献   

15.
《Coastal Engineering》1999,36(3):171-195
A morphological stability analysis is carried out for a long straight coast with a longshore bar. The situation with oblique wave incidence and a wave-driven longshore current is considered. The flow and sediment transport are described by a numerical modelling system. The models comprise: (i) a wave model with depth refraction, shoaling and wave breaking, (ii) a depth integrated model for wave driven currents and (iii) a sediment transport model for the bed load transport and the suspended load transport in combined waves and current. The direction of the sediment transport is taken to be parallel to the depth integrated mean current velocity, neglecting the effects of a bed slope and secondary currents. An instability is found to develop around the bar crest. The instability is periodic in the alongshore direction, and tends to form rip channels and to steepen the offshore face of the bar between the rip channels. The alongshore wave length of the most unstable perturbation is determined for different combinations of the wave conditions and the geometry of the profile.  相似文献   

16.
The three-dimensional numerical model with σ-coordinate transformation in the vertical direction is applied to the simulation of surface water waves and wave-induced laminar boundary layers. Unlike most of the previous investigations that solved the simplified one-dimensional boundary layer equation of motion and neglected the interaction between boundary layer and outside flow, the present model solves the full Navier–Stokes equations (NSE) in the entire domain from bottom to free surface. A non-uniform mesh system is used in the vertical direction to resolve the thin boundary layer. Linear wave, Stokes wave, cnoidal wave and solitary wave are considered. The numerical results are compared to analytical solutions and available experimental data. The numerical results agree favorably to all of the experimental data. It is found that the analytical solutions are accurate for both linear wave and Stokes wave but inadequate for cnoidal wave or solitary wave. The possible reason is that the existing analytical solutions for cnoidal and solitary waves adopt the first-order approximation for free stream velocity and thus overestimate the near bottom velocity. Besides velocity, the present model also provides accurate results for wave-induced bed shear stress.  相似文献   

17.
内孤立波作用下Spar平台的运动响应   总被引:9,自引:0,他引:9       下载免费PDF全文
建立了内孤立波作用下Spar平台运动响应的时域数值模型,并利用这一模型计算了内孤立波作用下Spar平台的运动响应。其中根据K-dv方程解出内孤立波的波面及水质点速度表达式,应用Morison公式求出内孤立波对Spar平台的水平作用力;由伯努力方程求解压强,通过在平台底面积分求出内孤立波对Spar平台的垂向作用力;在悬链线理论的基础上,分析锚链对平台的作用力,采用四阶Runge-Kutta方法求解运动方程。数值结果表明,当内孤立波不断接近海洋平台时,Spar平台受到的作用力随之增大,海洋平台会产生远大于表面波作用下的水平位移,从而证实了内孤立波是影响海洋平台稳定的重要因素,在海洋平台的设计和评估中内孤立波的作用是不可忽略的。  相似文献   

18.
波浪溢流现象使得海堤受到了越浪和溢流的联合作用,复杂的水动力过程会引起海堤后坡产生严重的侵蚀破坏。基于FLUENT软件建立了二维数值波浪溢流水槽模型,该模型运用UDF速度边界造波法分析在不同超高条件下海堤后坡流量和水流厚度的水力学特性。结果表明数值模拟结果与前人物理模型试验结果吻合,该模型可以真实地模拟出海堤波浪溢流现象。在此基础上进一步研究了波浪溢流中越浪和溢流在不同相对超高条件下的主导性作用,而后建立了十分准确的波浪溢流海堤后坡稳定水流厚度计算公式。  相似文献   

19.
人工鱼礁在波浪作用下水动力特性的研究对于人工鱼礁的工程安全与设计具有非常重要的意义。基于有限体积法,采用推板造波形式,通过利用VOF方法求解波面的方式建立了模拟人工鱼礁与波浪相互作用的三维数值波浪水槽。基于该数值模型研究了波浪作用下三角型镂空人工鱼礁的受力情况,并与物理模型实验结果进行比较,结果显示模拟得到的人工鱼礁受力和波浪形态均与实验结果吻合良好。根据数值模拟结果,拟合出波浪作用下人工鱼礁的水动力系数。单体三角型镂空人工鱼礁的速度力系数Cd随着Kc数、Re数的增加呈现减小的趋势;惯性力系数Cm则随着Kc数、Re数的增加呈现波动趋势。数值模拟结果显示在一个波浪周期内,人工鱼礁周围产生了较强的上升流和回流,其内部产生了明显的涡旋结构。研究结果为人工鱼礁的设计优化提供了理论依据。  相似文献   

20.
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.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号