首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 140 毫秒
1.
A numerical model describing the propagation and run-up process of nearshore tsunamis in the vicinity of shorelines is developed based on an approximate Riemann solver. The governing equations of the model are the nonlinear shallow-water equations. The governing equations are discretized explicitly by using a finite volume method. The nonlinear terms in the momentum equations are solved with the Harten-Lax-van Leer-Contact (HLLC) approximate Riemann solver. The developed model is first applied to prediction of water motions in a parabolic basin, and propagation and subsequent run-up process of nearshore tsunamis around a circular island. Computed results are then compared with available analytical solutions and laboratory measurements. Very reasonable agreements are observed.  相似文献   

2.
The interaction of a solitary wave with an array of surface-piercing vertical circular cylinders is investigated numerically. The wave motion is modeled by a set of generalized Boussinesq equations. The governing equations are discretized using a finite element method. The numerical model is validated against the experimental data of solitary wave reflection from a vertical wall and solitary wave scattering by a vertical circular cylinder respectively. The predicted wave surface elevation and the wave forces on the cylinder agree well with the experimental data. The numerical model is then employed to study solitary wave scattering by arrays of two circular cylinders and four circular cylinders respectively. The effect of wave direction on the wave forces and the wave runup on the cylinders is quantified.  相似文献   

3.
四叉树网格下的椭圆型缓坡方程数值模型研究   总被引:1,自引:1,他引:0  
唐军  李巧生  沈永明 《海洋学报》2013,35(5):162-168
波浪是近岸海域关键的水动力因素之一。考虑到近岸地形复杂、波浪演化显著的特点,建立了四叉树网格体系下的椭圆型缓坡方程数值模型,采用有限体积法对模型进行数值离散,应用GPBiCG(m, n)算法求解离散后的控制方程。模型中根据波浪波长布局计算网格,生成多层次四叉树网格,对复杂计算域有较好的适应性,并且在离散和方程求解中无需引入形函数、不产生复杂的交叉项,节约了存储空间和计算时间。将模型成功应用于物理模型实验及Acapulco海湾的波浪场数值模拟,结果表明该模型能够准确、高效地模拟近岸波浪场,可为近岸波浪场的模拟提供一定的理论和技术支持。  相似文献   

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

5.
《Ocean Modelling》1999,1(2-4):71-80
Ocean general circulation models (OGCMs) which represent the governing equations on a finite difference grid require shorter time steps with increasing resolution. Thus, until now, in the absence of filtering, the time step length has been determined by the smallest grid spacing within the model domain. Here we present a method for reducing the time step length (and increasing the number of time steps taken) at selected points in the grid, so as to minimise the computational cost of integrating the OGCM, whilst achieving numerical stability throughout the model domain without filtering. This variable time stepping method can be used to overcome numerical constraints associated with the convergence of longitude–latitude grids at the poles, and also to allow efficient integration of model domains with variable resolution. Examples of the computational saving are given.  相似文献   

6.
The benchmark simulations of wave run-up on a fixed single truncated circular cylinder and four circular cylinders are presented in this paper. Our in-house CFD solver naoe-FOAM-SJTU is adopted which is an unsteady two-phase CFD code based on the open source package OpenFOAM. The Navier-Stokes equations are employed as the governing equations, and the volume of fluid (VOF) method is applied for capturing the free surface. Monochromatic incident waves with the specified wave period and wave height are simulated and wave run-up heights around the cylinder are computed and recorded with numerical virtual wave probes. The relationship between the wave run-up heights and the incident wave parameters are analyzed. The numerical results indicate that the presented naoe-FOAM-SJTU solver can provide accurate predictions for the wave run-up on one fixed cylinder and four cylinders, which has been proved by the comparison of simulated results with experimental data.  相似文献   

7.
赵明  滕斌 《中国海洋工程》2004,18(2):267-280
The incompressible viscous uniform and shear flow past a circular cylinder is studied. The two-dimensional Navier-Stokes equations are solved by a finite element method. The governing equations are discretized by a weighted residual method in space. The stable three-step scheme is applied to the momentum equations in the time integration. The numerical model is firstly applied to the computation of the lid-driven cavity flow for its validation. The computed results agree well with the measured data and other numerical results. Then, it is used to simulate the viscous uniform and shear flow over a circular cylinder for Reynolds numbers from lO0 to lO00. The transient time interval before the vortex shedding occurs is shortened considerably by introduction of artificial perturbation. The computed Strouhal number, drag and lift coefficients agree well with the experimental data. The computation shows that the finite element model can be successfully applied to the viscous flow problem.  相似文献   

8.
近壁圆柱绕流水动力特性数值模拟与实验研究   总被引:3,自引:0,他引:3  
通过数值模拟和物理模型实验,对距壁面一定高度的圆柱绕流水动力特性进行了研究。数值模拟采用有限体积法对标准k-ε模式方程进行离散,采用SIMPLE算法进行求解,模拟绕流流场。在物理模型实验中,将PVC圆管制作的实验模型安放在水槽内,在圆管的跨中沿表面周向均匀布置水下压力传感器,用于测量绕流圆柱体表面动水压力分布。通过改变Re数和间隙比来分析它们对近壁圆柱绕流水动力特性的影响。基于数值流动显示技术,给出了近壁绕流流场的尾流流态分析。通过数值结果与实验结果的对比,对近壁绕流圆柱体的升力系数及其表面动水压力分布进行了研究,对比结果显示了较好的一致性。  相似文献   

9.
岛屿岛礁海域海浪能谱模型研究进展   总被引:2,自引:0,他引:2  
毛科峰  陈希  王亮 《海洋学报》2014,36(5):161-169
波浪能谱模型在岛屿岛礁海域的波浪预报研究和海洋工程中应用广泛,但存在模式计算格点无法充分体现岛屿岛礁的复杂地形特征和很难刻画波浪受到岛屿岛礁影响发生变形物理过程等两个关键问题。多重网格嵌套方案、岛屿次网格地形效应计算方案以及非结构网格、无网格、动态自适应四叉树网格等技术在体现岛屿岛礁复杂地形方面取得了较好的效果;将相位解析模型与波浪能谱模型优势互补是提高能谱模型对岛屿近岸波浪变形物理过程计算能力的一个有效方法。开展球坐标系下波作用密度谱方程的自适应四叉树网格求解方法研究,借鉴相位解析模型最新成果完善能谱模式的绕射、反射、底摩擦等物理过程,是提高岛屿岛礁海域海浪精细预报技术水平的前沿性、探索性研究方向。  相似文献   

10.
北部湾是一个半封闭的超浅海。本文在数值试验的基础上,提出了一个研究该海域台风风暴潮的数值模型。数值模式为二维深度平均流模型,采用嵌套细网格技术,细网格分辨率为沿经纬方向0.1°,细网格边界值由粗网格提供。台风风场计算采用Jelesnianski模型风场.模式方程组的数值解由交替方向隐式(ADI)方法积分得到。本文对该海域最常见的两种台风移行所引起的风暴潮进行了数值模拟。与几个潮汐观测站的增水记录比较,计算结果基本上反映了台风引起的水位变化,对研究和模拟该海域台风风暴潮是适用的,可用于该海域风暴潮数值预报试验。  相似文献   

11.
At the circular Babi Island in the Flores tsunami (1992) and pear shaped island in the Okushiri event (1993), unexpectedly large tsunami run‐up heights in the lee of conic islands were observed. The flume and basin physical model studies were conducted in the Coastal Hydraulic Laboratory, Engineering Research and Development Center, U.S. Army Corps of Engineers to provide a better understanding of the physical phenomena and verify numerical models used in predicting tsunami wave run‐up on beaches, islands, and vertical walls. Reasonably accurate comparison of run‐up height of solitary waves on a circular island has been obtained between laboratory experimental results and two‐dimensional computation model results. In this study we apply three‐dimensional RANS model to simulate wave run‐up on conical island. In the run‐up computation we obtain that 3D calculations are in very good comparison with laboratory and 2D numerical results. A close examination of the three‐dimensional velocity distribution around conical island to compare with depth‐integrated model is performed. It is shown that the velocity distribution along the vertical coordinate is not uniform: and velocity field is weaker in the bottom layer and higher on the sea surface. The maximum difference (about 40%) appears at the time when solitary wave reached the circular island.  相似文献   

12.
适于模拟不规则水域波浪的缓坡方程两种数值模型比较   总被引:1,自引:1,他引:0  
本文分析比较了适于不规则水域波浪模拟的椭圆型缓坡方程两种数值模型。两种数值模型均采用有限体积法离散,分别基于四叉树网格和非结构化三角形网格建立。首先结合近岸缓坡地形上波浪传播的经典物理模型实验对两种数值模型分别进行了验证,并结合计算结果对比分析了两种模型的计算精度和效率。计算结果表明,两种数值模型均可有效地模拟近岸波浪的传播变形;相对非结构化三角形网格下的模型,基于四叉树网格建立的数值模型在数值离散和求解过程中无需引入形函数、不产生复杂的交叉项,离散简单,易于程序实现,且节约计算存储空间,计算效率高。  相似文献   

13.
在二阶 Boussinesq 方程基础上,通过引入含水深导数项对该方程进行了理论上的改进,使得该方程在应用于无限沙坝 Bragg反射问题时与理论解析解在更大范围内符合.基于该改进的高阶 Boussinesq 方程,在非交错网格下建立了混合 4 阶的Adams-Bashforth- Moulton 格式的数学模型.将数值模型应用到有限个连续沙坝上波浪传播变形问题的数值模拟中,通过两点法给出数值波浪反射系数,将这些反射系数与已有的实验数据进行对比,对比表明改进后的模型计算出的反射系数与实验结果吻合更好,这验证了本文理论改进的有效性.  相似文献   

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

15.
Numerical modeling of nonlinear water waves over heterogeneous porous beds   总被引:1,自引:0,他引:1  
Eric C. Cruz  Qin Chen   《Ocean Engineering》2007,34(8-9):1303-1321
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.  相似文献   

16.
Although the finite difference method is computationally efficient, it is acknowledged to be inferior when dealing with flow-over on structures with a complex geometry because of its rectilinear grid system. Therefore, we developed a numerical procedure that can cope with flow over structures with complex shapes while, at the same time, retaining the simplicity and efficiency of a rectilinear grid system. We used the immersed boundary method, which involves application of immersed boundary forces at solid boundaries rather than conventional boundary conditions, to investigate wave interactions with coastal structures in a three-dimensional numerical wave tank by solving the Navier–Stokes equations for two-phase flows. We simulated the run-up of a solitary wave around a circular island. Maximum run-up heights were computed around the island and compared with available laboratory measurements and previous numerical results. The three-dimensional features of the run-up process were analyzed in detail and compared with those of depth-integrated equations models.  相似文献   

17.
《Coastal Engineering》2001,42(3):219-239
This paper describes an adaptive quadtree-based 2DH wave–current interaction model for evaluating nearly horizontal wave-induced currents in the surf-zone. The model accounts for wave breaking, shoaling, refraction, diffraction, wave–current interaction, set-up and set-down, mixing processes (turbulent diffusion), bottom frictional effects, and movement of the land–water interface at the shoreline. The wave period- and depth-averaged governing equations, which conserve mass, momentum, energy and wave action, are discretised explicitly by means of an Adams–Bashforth second-order finite difference technique on adaptive hierarchical staggered quadtree grids. Grid adaptation is achieved through seeding points distributed according to flow criteria (e.g. local current gradients). The model is verified for nearshore circulation at a sinusoidal beach and nearshore currents at a multi-cusped beach. Reasonable agreement is obtained with experimental data from da Silva Lima [da Silva Lima, S.S.L., 1981. Wave-induced Nearshore Currents. PhD Thesis, Department of Civil Engineering, University of Liverpool] and Borthwick et al. [Borthwick, A.G.L., Foote, Y.L.M., Ridehalgh, A., 1997. Nearshore measurements at a cusped beach in the UK Coastal Research Facility, Coastal Dynamics '97, Plymouth, 953–962]. The modelling approach presented herein should be useful in simulating nearshore processes in complicated natural coastal domains. Of particular value is the local grid enrichment capability, which permits refined modelling of important localised flow behaviour such as rip currents and surf-zone circulation systems.  相似文献   

18.
A numerical scheme for solving the class of extended Boussinesq equations is presented. Unlike previous schemes, where the governing equations are integrated through time using a fourth-order method, a second-order Godunov-type scheme is used thus saving storage and computational resources. The spatial derivatives are discretised using a combination of finite-volume and finite-difference methods. A fourth-order MUSCL reconstruction technique is used to compute the values at the cell interfaces for use in the local Riemann problems, whilst the bed source and dispersion terms are discretised using centred finite-differences of up to fourth-order accuracy. Numerical results show that the class of extended Boussinesq equations can be accurately solved without the need for a fourth-order time discretisation, thus improving the computational speed of Boussinesq-type numerical models. The numerical scheme has been applied to model a number of standard test cases for the extended Boussinesq equations and comparisons made to physical wave flume experiments.  相似文献   

19.
黄海、渤海盐度的垂直结构具有典型的自模性,而其水平分布又受平流、水平扩散效应及径流等因素的影响。本文根据黄海、渤海实测资料拟合了盐度垂直剖面的自模函数,并结合描述表、底层盐度及上均匀层厚度这3个特征量水平分布的方程,给出盐度三维结构的准三维模式。在模式中,综合考虑了海面风和热输入的强迫作用以及流场的平流、侧向混合及底层混合的影响因素,同时还考虑了径流、蒸发及降水的作用,较客观地反映了盐度的三维分布及其变化的物理过程。试报结果分析表明,模式的功能较好,结果令人满意。  相似文献   

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

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

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